Haplotype-resolved sweet potato genome traces back its hexaploidization history

Genome-Wide Identification and Expression Analysis of the Xyloglucan Endotransglucosylase/Hydrolase Gene Family in Sweet Potato [Ipomoea batatas (L.) Lam]

The xyloglucan endotransglucosylase/hydrolase (XET/XEH, also named XTH) family is a multigene family, the function of which plays a significant role in cell-wall rebuilding and stress tolerance in plants. However, the specific traits of the XTH gene family members and their expression pattern in different tissues and under stress have not been carried out in sweet potato. Thirty-six XTH genes were identified in I. batatas, all of which had conserved structures (Glyco_hydro_16). Based on Neighbor-Joining phylogenetic analysis the IbXTHs can be divided into three subfamilies-the I/II, IIIA, and IIIB subfamilies, which were unevenly distributed on 13 chromosomes, with the exception of Chr9 and Chr15. Multiple cis-acting regions related to growth and development, as well as stress responses, may be found in the IbXTH gene promoters. The segmental duplication occurrences greatly aided the evolution of IbXTHs. The results of a collinearity analysis showed that the XTH genes of sweet potato shared evolutionary history with three additional species, including A. thaliana, G. max, and O. sativa. Additionally, based on the transcriptome sequencing data, the results revealed that the IbXTHs have different expression patterns in leaves, stems, the root body (RB), the distal end (DE), the root stock (RS), the proximal end (PE), the initiative storage root (ISR), and the fibrous root (FR), and many of them are well expressed in the roots. Differentially expressed gene (DEG) analysis of FRs after hormone treatment of the roots indicated that IbXTH28 and IbXTH30 are up-regulated under salicylic acid (SA) treatment but down-regulated under methyl jasmonate (MeJA) treatment. Attentionally, there were only two genes showing down-regulation under the cold and drought treatment. Collectively, all of the findings suggested that genes from the XTH family are crucial for root specificity. This study could provide a theoretical basis for further research on the molecular function of sweet potato XTH genes.

A systematical genome-wide analysis and screening of WRKY transcription factor family engaged in abiotic stress response in sweetpotato

BACKGROUND

WRKY transcription factors play pivotal roles in regulating plant multiple abiotic stress tolerance, however, a genome-wide systematical analysis of WRKY genes in sweetpotato is still missing.

RESULTS

Herein, 84 putative IbWRKYs with WRKY element sequence variants were identified in sweetpotato reference genomes. Fragment duplications, rather than tandem duplications, were shown to play prominent roles in IbWRKY gene expansion. The collinearity analysis between IbWRKYs and the related orthologs from other plants further depicted evolutionary insights into IbWRKYs. Phylogenetic relationships displayed that IbWRKYs were divided into three main groups (I, II and III), with the support of the characteristics of exon-intron structures and conserved protein motifs. The IbWRKY genes, mainly from the group Ib, displayed remarkable and diverse expression profiles under multiple abiotic stress (NaCl, PEG6000, cold and heat) and hormone (ABA, ACC, JA and SA) treatments, which were determined by RNA-seq and qRT-PCR assays, suggesting their potential roles in mediating particular stress responses. Moreover, IbWRKY58L could interact with IbWRKY82 as revealed by yeast two-hybrid based on the protein interaction network screening. And abiotic stress-remarkably induced IbWRKY21L and IbWRKY51 were shown to be localized in the nucleus and had no transactivation activities.

CONCLUSION

These results provide valuable insights into sweetpotato IbWRKYs and will lay a foundation for further exploring functions and possible regulatory mechanisms of IbWRKYs in abiotic stress tolerance.

Genome-Wide Identification and Expression Analysis of SWEET Family Genes in Sweet Potato and Its Two Diploid Relatives

Sugar Will Eventually be Exported Transporter (SWEET) proteins are key transporters in sugar transportation. They are involved in the regulation of plant growth and development, hormone crosstalk, and biotic and abiotic stress responses. However, SWEET family genes have not been explored in the sweet potato. In this study, we identified 27, 27, and 25 SWEETs in cultivated hexaploid sweet potato (Ipomoea batatas, 2n = 6x = 90) and its two diploid relatives, Ipomoea trifida (2n = 2x = 30) and Ipomoea triloba (2n = 2x = 30), respectively. These SWEETs were divided into four subgroups according to their phylogenetic relationships with Arabidopsis. The protein physiological properties, chromosome localization, phylogenetic relationships, gene structures, promoter cis-elements, protein interaction networks, and expression patterns of these 79 SWEETs were systematically investigated. The results suggested that homologous SWEETs are differentiated in sweet potato and its two diploid relatives and play various vital roles in plant growth, tuberous root development, carotenoid accumulation, hormone crosstalk, and abiotic stress response. This work provides a comprehensive comparison and furthers our understanding of the SWEET genes in the sweet potato and its two diploid relatives, thereby supplying a theoretical foundation for their functional study and further facilitating the molecular breeding of sweet potato.

Natural allelic variation confers high resistance to sweet potato weevils in sweet potato

Sweet potato (Ipomoea batatas L.) is a major root crop worldwide. Sweet potato weevils (SPWs) pose one of the most significant challenges to sweet potato production in tropical and subtropical regions, causing deleterious economic and environmental effects. Characterizing the mechanisms underlying natural resistance to SPWs is therefore crucial; however, the genetic basis of host SPW resistance (SPWR) remains unclear. Here we obtained two sweet potato germplasm with high SPWR and, by map-based cloning, revealed two major SPW-resistant genes-SPWR1 and SPWR2-that are important regulators of natural defence against SPWs. The SPW-induced WRKY transcriptional factor SPWR1 directly activates the expression of SPWR2, and SPWR2, the conserved dehydroquinate synthase, promotes the accumulation of quinate derivative metabolites that confer SPWR in sweet potato. Generally, our results provide new insights into the molecular mechanism underlying sweet potato-SPW interactions and will aid future efforts to achieve eco-friendly SPW management.

Cut-dip-budding delivery system enables genetic modifications in plants without tissue culture

Of the more than 370 000 species of higher plants in nature, fewer than 0.1% can be genetically modified due to limitations of the current gene delivery systems. Even for those that can be genetically modified, the modification involves a tedious and costly tissue culture process. Here, we describe an extremely simple cut-dip-budding (CDB) delivery system, which uses Agrobacterium rhizogene to inoculate explants, generating transformed roots that produce transformed buds due to root suckering. We have successfully used CDB to achieve the heritable transformation of plant species in multiple plant families, including two herbaceous plants (Taraxacum kok-saghyz and Coronilla varia), a tuberous root plant (sweet potato), and three woody plant species (Ailanthus altissima, Aralia elata, and Clerodendrum chinense). These plants have previously been difficult or impossible to transform, but the CDB method enabled efficient transformation or gene editing in them using a very simple explant dipping protocol, under non-sterile conditions and without the need for tissue culture. Our work suggests that large numbers of plants could be amenable to genetic modifications using the CDB method.

Integrated analysis of carotenoid metabolites and transcriptome identifies key genes controlling carotenoid compositions and content in sweetpotato tuberous roots (Ipomoea batatas L.)

Sweetpotato (Ipomoea batatas L.) with different depths of yellow color contains different compositions of carotenoids, which are beneficial for human health. In this study, we performed an integrated analysis of metabolomic and transcriptomic to identify key genes playing a major role in carotenoid coloration in sweetpotato tuberous roots. Herein, 14 carotenoids were identified in five sweetpotatoes. Orange-red and orange cultivars were dominated by β-carotene (385.33 μg/g and 85.07 μg/g), yellow cultivar had a high β-cryptoxanthin (11.23 μg/g), light-yellow cultivar was rich in zeaxanthin (5.12 μg/g), whereas lutein (3.34 μg/g) was the main carotenoid in white cultivar. Furthermore, 27 differentially expressed genes involved in carotenoid metabolism were identified based on comparative transcriptome. Weighted gene co-expression network analysis identified 15 transcription factors highly associated with carotenoid content in sweetpotatoes. These results provide valuable information for revealing the regulatory mechanism of carotenoid metabolism in different-colored sweetpotato tuberous roots.

Resequencing of sweetpotato germplasm resources reveals key loci associated with multiple agronomic traits

Sweetpotato is an important crop that exhibits hexaploidy and high heterozygosity, which limits gene mining for important agronomic traits. Here, 314 sweetpotato germplasm resources were deeply resequenced, and 4 599 509 SNPs and 846 654 InDels were generated, among which 196 124 SNPs were nonsynonymous and 9690 InDels were frameshifted. Based on the Indels, genome-wide marker primers were designed, and 3219 of 40 366 primer pairs were selected to construct the core InDel marker set. The molecular ID of 104 sweetpotato samples verified the availability of these primers. The sweetpotato population structures were then assessed through multiple approaches using SNPs, and diverse approaches demonstrated that population stratification was not obvious for most Chinese germplasm resources. As many as 20 important agronomic traits were evaluated, and a genome-wide association study was conducted on these traits. A total of 19 high-confidence loci were detected in both models. These loci included several candidate genes, such as IbMYB1, IbZEP1, and IbYABBY1, which might be involved in anthocyanin metabolism, carotenoid metabolism, and leaf morphogenesis, respectively. Among them, IbZEP1 and IbYABBY1 were first reported in sweetpotato. The variants in the promoter and the expression levels of IbZEP1 were significantly correlated with flesh color (orange or not orange) in sweetpotato. The expression levels of IbYABBY1 were also correlated with leaf shape. These results will assist in genetic and breeding studies in sweetpotato.

AtGCS promoter-driven clustered regularly interspaced short palindromic repeats/Cas9 highly efficiently generates homozygous/biallelic mutations in the transformed roots by Agrobacterium rhizogenes-mediated transformation

Agrobacterium rhizogenes-mediated (ARM) transformation is an efficient and powerful tool to generate transgenic roots to study root-related biology. For loss-of-function studies, transgenic-root-induced indel mutations by CRISPR/Cas9 only with homozygous/biallelic mutagenesis can exhibit mutant phenotype(s) (excluding recessive traits). However, a low frequency of homozygous mutants was produced by a constitutive promoter to drive Cas9 expression. Here, we identified a highly efficient Arabidopsis thaliana gamma-glutamylcysteine synthetase promoter, termed AtGCSpro, with strong activity in the region where the root meristem will initiate and in the whole roots in broad eudicots species. AtGCSpro achieved higher homozygous/biallelic mutation efficiency than the most widely used CaMV 35S promoter in driving Cas9 expression in soybean, Lotus japonicus, and tomato roots. Using the pAtGCSpro-Cas9 system, the average homozygous/biallelic mutation frequency is 1.7-fold and 8.3-fold higher than the p2 × 35Spro-Cas9 system for single and two target site(s) in the genome, respectively. Our results demonstrate the advantage of the pAtGCSpro-Cas9 system used in ARM transformation, especially its great potential in diploids with multiple-copy genes targeted mutations and polyploid plants with multiplex genome editing. AtGCSpro is conservatively active in various eudicots species, suggesting that AtGCSpro might be applied in a wide range of dicots species.

Novel Insights into the Nature of Intraspecific Genome Size Diversity in Cannabis sativa L

Cannabis sativa has been used for millennia in traditional medicine for ritual purposes and for the production of food and fibres, thus, providing important and versatile services to humans. The species, which currently has a worldwide distribution, strikes out for displaying a huge morphological and chemical diversity. Differences in Cannabis genome size have also been found, suggesting it could be a useful character to differentiate between accessions. We used flow cytometry to investigate the extent of genome size diversity across 483 individuals belonging to 84 accessions, with a wide range of wild/feral, landrace, and cultivated accessions. We also carried out sex determination using the MADC2 marker and investigated the potential of flow cytometry as a method for early sex determination. All individuals were diploid, with genome sizes ranging from 1.810 up to 2.152 pg/2C (1.189-fold variation), apart from a triploid, with 2.884 pg/2C. Our results suggest that the geographical expansion of Cannabis and its domestication had little impact on its overall genome size. We found significant differences between the genome size of male and female individuals. Unfortunately, differences were, however, too small to be discriminated using flow cytometry through the direct processing of combined male and female individuals.

Plant metabolic gene clusters in the multi-omics era

Secondary metabolism in plants gives rise to a vast array of small-molecule natural products. The discovery of operon-like gene clusters in plants has provided a new perspective on the evolution of specialized metabolism and the opportunity to rapidly advance the metabolic engineering of natural product production. Here, we review historical aspects of the study of plant metabolic gene clusters as well as general strategies for identifying plant metabolic gene clusters in the multi-omics era. We also emphasize the exploration of their natural variation and evolution, as well as new strategies for the prospecting of plant metabolic gene clusters and a deeper understanding of how their structure influences their function.

Genome-Wide Identification of the A20/AN1 Zinc Finger Protein Family Genes in Ipomoea batatas and Its Two Relatives and Function Analysis of IbSAP16 in Salinity Tolerance

Stress-associated protein (SAP) genes—encoding A20/AN1 zinc-finger domain-containing proteins—play pivotal roles in regulating stress responses, growth, and development in plants. They are considered suitable candidates to improve abiotic stress tolerance in plants. However, the SAP gene family in sweetpotato (Ipomoea batatas) and its relatives is yet to be investigated. In this study, 20 SAPs in sweetpotato, and 23 and 26 SAPs in its wild diploid relatives Ipomoea triloba and Ipomoea trifida were identified. The chromosome locations, gene structures, protein physiological properties, conserved domains, and phylogenetic relationships of these SAPs were analyzed systematically. Binding motif analysis of IbSAPs indicated that hormone and stress responsive cis-acting elements were distributed in their promoters. RT-qPCR or RNA-seq data revealed that the expression patterns of IbSAP, ItbSAP, and ItfSAP genes varied in different organs and responded to salinity, drought, or ABA (abscisic acid) treatments differently. Moreover, we found that IbSAP16 driven by the 35 S promoter conferred salinity tolerance in transgenic Arabidopsis. These results provided a genome-wide characterization of SAP genes in sweetpotato and its two relatives and suggested that IbSAP16 is involved in salinity stress responses. Our research laid the groundwork for studying SAP-mediated stress response mechanisms in sweetpotato.

Two independent allohexaploidizations and genomic fractionation in Solanales

Solanales, an order of flowering plants, contains the most economically important vegetables among all plant orders. To date, many Solanales genomes have been sequenced. However, the evolutionary processes of polyploidization events in Solanales and the impact of polyploidy on species diversity remain poorly understood. We compared two representative Solanales genomes (Solanum lycopersicum L. and Ipomoea triloba L.) and the Vitis vinifera L. genome and confirmed two independent polyploidization events. Solanaceae common hexaploidization (SCH) and Convolvulaceae common hexaploidization (CCH) occurred ∼43-49 and ∼40-46 million years ago (Mya), respectively. Moreover, we identified homologous genes related to polyploidization and speciation and constructed multiple genomic alignments with V. vinifera genome, providing a genomic homology framework for future Solanales research. Notably, the three polyploidization-produced subgenomes in both S. lycopersicum and I. triloba showed significant genomic fractionation bias, suggesting the allohexaploid nature of the SCH and CCH events. However, we found that the higher genomic fractionation bias of polyploidization-produced subgenomes in Solanaceae was likely responsible for their more abundant species diversity than that in Convolvulaceae. Furthermore, through genomic fractionation and chromosomal structural variation comparisons, we revealed the allohexaploid natures of SCH and CCH, both of which were formed by two-step duplications. In addition, we found that the second step of two paleohexaploidization events promoted the expansion and diversity of β-amylase (BMY) genes in Solanales. These current efforts provide a solid foundation for future genomic and functional exploration of Solanales.

TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies

Vegetables are an indispensable part of the daily diet of humans. Therefore, it is vital to systematically study the genomic data of vegetables and build a platform for data sharing and analysis. In this study, a comprehensive platform for vegetables with a user-friendly Web interface-The Vegetable Information Resource (TVIR, http://tvir.bio2db.com)-was built based on the genomes of 59 vegetables. TVIR database contains numerous important functional genes, including 5215 auxin genes, 2437 anthocyanin genes, 15 002 flowering genes, 79 830 resistance genes, and 2639 glucosinolate genes of 59 vegetables. In addition, 2597 N6-methyladenosine (m6A) genes were identified, including 513 writers, 1058 erasers, and 1026 readers. A total of 2 101 501 specific clustered regularly interspaced short palindromic repeat (CRISPR) guide sequences and 17 377 miRNAs were detected and deposited in TVIR database. Information on gene synteny, duplication, and orthologs is also provided for 59 vegetable species. TVIR database contains 2 346 850 gene annotations by the Swiss-Prot, TrEMBL, Gene Ontology (GO), Pfam, and Non-redundant (Nr) databases. Synteny, Primer Design, Blast, and JBrowse tools are provided to facilitate users in conducting comparative genomic analyses. This is the first large-scale collection of vegetable genomic data and bioinformatic analysis. All genome and gene sequences, annotations, and bioinformatic results can be easily downloaded from TVIR. Furthermore, transcriptome data of 98 vegetables have been collected and collated, and can be searched by species, tissues, or different growth stages. TVIR is expected to become a key hub for vegetable research globally. The database will be updated with newly assembled vegetable genomes and comparative genomic studies in the future.

Transcriptome analysis of sweet potato responses to potassium deficiency

Background

As one of three essential nutrients, potassium is regarded as a main limiting factor for growth and development in plant. Sweet potato (Ipomoea batatas L.) is one of seven major food crops grown worldwide, and is both a nutrient-rich food and a bioenergy crop. It is a typical ‘K-favoring’ crop, and the level of potassium ion (K⁺) supplementation directly influences its production. However, little is known about the transcriptional changes in sweet potato genes under low-K⁺ conditions. Here, we analyzed the transcriptomic profiles of sweet potato roots in response to K⁺ deficiency to determine the effect of low-K⁺ stress on this economically important crop.

Results

The roots of sweet potato seedlings with or without K⁺ treatment were harvested and used for transcriptome analyses. The results showed 559 differently expressed genes (DEGs) in low and high K⁺ groups. Among the DEGs, 336 were upregulated and 223 were downregulated. These DEGs were involved in transcriptional regulation, calcium binding, redox-signaling, biosynthesis, transport, and metabolic process. Further analysis revealed previously unknow genes involved in low-K⁺ stress, which could be investigated further to improve low K⁺ tolerance in plants. Confirmation of RNA-sequencing results using qRT-PCR displayed a high level of consistency between the two experiments. Analysis showed that many auxin-, ethylene- and jasmonic acid-related genes respond to K⁺ deficiency, suggesting that these hormones have important roles in K⁺ nutrient signaling in sweet potato.

Conclusions

According to the transcriptome data of sweet potato, various DEGs showed transcriptional changes in response to low-K⁺ stress. However, the expression level of some kinases, transporters, transcription factors (TFs), hormone-related genes, and plant defense-related genes changed significantly, suggesting that they have important roles during K⁺ deficiency. Thus, this study identifies potential genes for genetic improvement of responses to low-K⁺ stress and provides valuable insight into the molecular mechanisms regulating low K⁺ tolerance in sweet potato. Further research is required to clarify the function of these DEGs under low-K⁺ stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08870-5.

Exploring and exploiting genetics and genomics for sweetpotato improvement: Status and perspectives

Sweetpotato (Ipomoea batatas (L.) Lam.) is one of the most important root crops cultivated worldwide. Because of its adaptability, high yield potential, and nutritional value, sweetpotato has become an important food crop, particularly in developing countries. To ensure adequate crop yields to meet increasing demand, it is essential to enhance the tolerance of sweetpotato to environmental stresses and other yield-limiting factors. The highly heterozygous hexaploid genome of I. batatas complicates genetic studies and limits improvement of sweetpotato through traditional breeding. However, application of next-generation sequencing and high-throughput genotyping and phenotyping technologies to sweetpotato genetics and genomics research has provided new tools and resources for crop improvement. In this review, we discuss the genomics resources that are available for sweetpotato, including the current reference genome, databases, and available bioinformatics tools. We systematically review the current state of knowledge on the polyploid genetics of sweetpotato, including studies of its origin and germplasm diversity and the associated mapping of important agricultural traits. We then outline the conventional and molecular breeding approaches that have been applied to sweetpotato. Finally, we discuss future goals for genetic studies of sweetpotato and crop improvement via breeding in combination with state-of-the-art multi-omics approaches such as genomic selection and gene editing. These approaches will advance and accelerate genetic improvement of this important root crop and facilitate its sustainable global production.

A chromosome-level reference genome of a Convolvulaceae species Ipomoea cairica

Ipomoea cairica is a perennial creeper that has been widely introduced as a garden ornamental across tropical, subtropical, and temperate regions. Because it grows extremely fast and spreads easily, it has been listed as an invasive species in many countries. Here, we constructed the chromosome-level reference genome of Ipomoea cairica by Pacific Biosciences HiFi and Hi-C sequencing, with the assembly size of 733.0 Mb, the contig N50 of 43.8 Mb, the scaffold N50 of 45.7 Mb, and the Benchmarking Universal Single-Copy Orthologs complete rate of 98.0%. Hi-C scaffolding assigned 97.9% of the contigs to 15 pseudo-chromosomes. Telomeric repeat analysis reveals that 7 of the 15 pseudo-chromosomes are gapless and telomere to telomere. The transposable element content of Ipomoea cairica is 73.4%, obviously higher than that of other Ipomoea species. A total of 38,115 protein-coding genes were predicted, with the Benchmarking Universal Single-Copy Orthologs complete rate of 98.5%, comparable to that of the genome assembly, and 92.6% of genes were functional annotated. In addition, we identified 3,039 tRNA genes and 2,403 rRNA genes in the assembled genome. Phylogenetic analysis showed that Ipomoea cairica formed a clade with Ipomoea aquatica, and they diverged from each other 8.1 million years ago. Through comparative genome analysis, we reconfirmed that a whole genome triplication event occurred specific to Convolvulaceae family and in the ancestor of the genus Ipomoea and Cuscuta. This high-quality reference genome of Ipomoea cairica will greatly facilitate the studies on the molecular mechanisms of its rapid growth and invasiveness.

IbMYB308, a Sweet Potato R2R3-MYB Gene, Improves Salt Stress Tolerance in Transgenic Tobacco

The MYB (v-myb avian myeloblastosis viral oncogene homolog) transcription factor family plays an important role in plant growth, development, and response to biotic and abiotic stresses. However, the gene functions of MYB transcription factors in sweet potato (Ipomoea batatas (L.) Lam) have not been elucidated. In this study, an MYB transcription factor gene, IbMYB308, was identified and isolated from sweet potato. Multiple sequence alignment showed that IbMYB308 is a typical R2R3-MYB transcription factor. Further, quantitative real-time PCR (qRT-PCR) analysis revealed that IbMYB308 was expressed in root, stem, and, especially, leaf tissues. Moreover, it showed that IbMYB308 had a tissue-specific profile. The experiment also showed that the expression of IbMYB308 was induced by different abiotic stresses (20% PEG-6000, 200 mM NaCl, and 20% H₂O₂). After a 200 mM NaCl treatment, the expression of several stress-related genes (SOD, POD, APX, and P5CS) was upregulation in transgenic plants, and the CAT activity, POD activity, proline content, and protein content in transgenic tobacco had increased, while MDA content had decreased. In conclusion, this study demonstrated that IbMYB308 could improve salt stress tolerance in transgenic tobacco. These findings lay a foundation for future studies on the R2R3-MYB gene family of sweet potato and suggest that IbMYB308 could potentially be used as an important positive factor in transgenic plant breeding to improve salt stress tolerance in sweet potato plants.

Genome-Wide Characterization of Nitrogenase Reductase (nifH) Genes in the Sweet Potato [Ipomoea batatas (L.) Lam] and Its Wild Ancestors

The sweet potato (Ipomoea batatas (L.) Lam.) is an important and widely grown crop, and the nitrogenase reductase (nifH) gene is the most widely sequenced marker gene used to identify nitrogen-fixing bacteria and archaea. There have been many examples of the isolation of the diazotrophic endophytes in sweet potatoes, and there has been no report on whether sweet potatoes and their wild ancestors harbored nifH genes. In this study, a comprehensive analysis of nifH genes has been conducted on these species by using bioinformatics and molecular biology methods. A total of 20, 19 and 17 nifH genes were identified for the first time in sweet potatoes, I. trifida and I. triloba, respectively. Based on a phylogenetic analysis, all of the nifH genes, except for g10233.t1, itf14g14040.t1 and itb14g15470.t1, were clustered into five independent clades: I, II, III, IV and V. The nifH genes clustered in the same phylogenetic branch showed a more similar distribution of conserved motifs and exons–introns than those of the other ones. All of the identified genes were further mapped on the 15 chromosomes of the sweet potato, I. trifida and I. triloba. No segmental duplication was detected in each genome of three Ipomoea species, and 0, 8 and 7 tandemly duplicated gene pairs were detected in the genome of the sweet potato, I. trifida and I. triloba, respectively. Synteny analysis between the three Ipomoea species revealed that there were 7, 7 and 8 syntenic gene pairs of nifH genes detected between the sweet potato and I. trifida, between the sweet potato and I. triloba and between I. trifida and I. triloba, respectively. All of the duplicated and syntenic nifH genes were subjected to purifying selection inside duplicated genomic elements during speciation, except for the tandemly duplicated gene pair itf11g07340.t2_itf11g07340.t3, which was subjected to positive selection. Different expression profiles were detected in the sweet potato, I. trifida and I. triloba. According to the above results, four nifH genes of the sweet potato (g950, g16683, g27094 and g33987) were selected for quantitative real-time polymerase chain reaction (qRT-PCR) analysis in two sweet potato cultivars (Eshu 15 and Long 9) under nitrogen deficiency (N0) and normal (N1) conditions. All of them were upregulated in the N1 treatment and were consistent with the analysis of the RNA-seq data. We hope that these results will provide new insights into the nifH genes in the sweet potato and its wild ancestors and will contribute to the molecular breeding of sweet potatoes in the future.

Comparative transcriptome analysis of purple-fleshed sweet potato and its yellow-fleshed mutant provides insight into the transcription factors involved in anthocyanin biosynthesis in tuberous root

In various plant species, many transcription factors (TFs), such as MYB, bHLH, and WD40, have been identified as regulators of anthocyanin biosynthesis in underground organs. However, the regulatory elements of anthocyanin biosynthesis in the tuberous roots of sweet potato have not been elucidated yet. Here, we selected the purple-fleshed sweet potato cultivar "Zhezi1" (ZZ ^P ) and its spontaneous yellow-fleshed mutant "Xinli" (XL ^Y ) to investigate the regulatory mechanism of the anthocyanin biosynthesis in the tuberous roots of sweet potato. By analyzing the IbMYB1 genotype in ZZ ^P and XL ^Y , we found that the IbMYB1-2, a MYB TF involved in anthocyanin biosynthesis, was missing in the XL ^Y genome, which might lead to an extreme decrease in anthocyanins in XL ^Y . A comparative transcriptome analysis of ZZ ^P and XL ^Y was conducted to find the TFs involved in anthocyanin biosynthesis in ZZ ^P and XL ^Y . The anthocyanin structural genes were significantly enriched among the differentially expressed genes. Moreover, one MYB activator (IbMYB1), one bHLH (IbbHLH2), three WRKY activator candidates (IbWRKY21, IbWRKY24, and IbWRKY44), and two MYB repressors (IbMYB27 and IbMYBx-ZZ) were highly expressed in ZZ ^P accompanied with anthocyanin structural genes. We also tested the expression of these TFs in six purple- and two orange-fleshed sweet potato cultivars. Interestingly, most of these TFs were significantly positively correlated with anthocyanin contents in these cultivars. The function of the anthocyanin biosynthesis repression of IbMYB27 and IbMYBx-ZZ was verified through transient co-transformation with IbMYB1 into tobacco leaves. Further functional verification of the above TFs was conducted by Y2H, BiFC, and dual-luciferase assays. These tests showed that the MYB-bHLH-WD40/MYB-bHLH-WD40-WRKY complex activated the promoter of anthocyanin structural gene IbDFR and promoters for IbWRKY44, IbMYB27, and IbMYBx-ZZ, indicating reinforcement and feedback regulation to maintain the level of anthocyanin accumulation in the tuberous roots of purple-fleshed sweet potato. These results may provide new insights into the regulatory mechanism of anthocyanin biosynthesis and accumulation in underground organs of sweet potatoes.

Transcriptomic analysis of tuberous root in two sweet potato varieties reveals the important genes and regulatory pathways in tuberous root development

Background

Tuberous root formation and development is a complex process in sweet potato, which is regulated by multiple genes and environmental factors. However, the regulatory mechanism of tuberous root development is unclear.

Results

In this study, the transcriptome of fibrous roots (R0) and tuberous roots in three developmental stages (Rl, R2, R3) were analyzed in two sweet potato varieties, GJS-8 and XGH. A total of 22,914 and 24,446 differentially expressed genes (DEGs) were identified in GJS-8 and XGH respectively, 15,920 differential genes were shared by GJS-8 and XGH. KEGG pathway enrichment analysis showed that the DEGs shared by GJS-8 and XGH were mainly involved in “plant hormone signal transduction” “starch and sucrose metabolism” and “MAPK signal transduction”. Trihelix transcription factor (Tai6.25300) was found to be closely related to tuberous root enlargement by the comprehensive analysis of these DEGs and weighted gene co-expression network analysis (WGCNA).

Conclusion

A hypothetical model of genetic regulatory network for tuberous root development of sweet potato is proposed, which emphasizes that some specific signal transduction pathways like “plant hormone signal transduction” “Ca²⁺signal” “MAPK signal transduction” and metabolic processes including “starch and sucrose metabolism” and “cell cycle and cell wall metabolism” are related to tuberous root development in sweet potato. These results provide new insights into the molecular mechanism of tuberous root development in sweet potato.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08670-x.

Genome-Wide Identification and Expression Analysis of Expansin Gene Family in the Storage Root Development of Diploid Wild Sweetpotato Ipomoea trifida

Expansins play important roles in root growth and development, but investigation of the expansin gene family has not yet been reported in Ipomoea trifida, and little is known regarding storage root (SR) development. In this work, we identified a total of 37 expansins (ItrEXPs) in our previously reported SR-forming I. trifida strain Y22 genome, which included 23 ItrEXPAs, 4 ItrEXPBs, 2 ItrEXLAs and 8 ItrEXLBs. The phylogenetic relationship, genome localization, subcellular localization, gene and protein structure, promoter cis-regulating elements, and protein interaction network were systematically analyzed to reveal the possible roles of ItrEXPs in the SR development of I. trifida. The gene expression profiling in Y22 SR development revealed that ItrEXPAs and ItrEXLBs were down-regulated, and ItrEXPBs were up-regulated while ItrEXLAs were not obviously changed during the critical period of SR expansion, and might be beneficial to SR development. Combining the tissue-specific expression in young SR transverse sections of Y22 and sweetpotato tissue, we deduced that ItrEXLB05, ItrEXLB07 and ItrEXLB08 might be the key genes for initial SR formation and enlargement, and ItrEXLA02 might be the key gene for root growth and development. This work provides new insights into the functions of the expansin gene family members in I. trifida, especially for EXLA and EXLB subfamilies genes in SR development.

De novo assembly of the complete mitochondrial genome of sweet potato (Ipomoea batatas [L.] Lam) revealed the existence of homologous conformations generated by the repeat-mediated recombination

Sweet potato (Ipomoea batatas [L.] Lam) is an important food crop, an excellent fodder crop, and a new type of industrial raw material crop. The lack of genomic resources could affect the process of industrialization of sweet potato. Few detailed reports have been completed on the mitochondrial genome of sweet potato. In this research, we sequenced and assembled the mitochondrial genome of sweet potato and investigated its substructure. The mitochondrial genome of sweet potato is 270,304 bp with 23 unique core genes and 12 variable genes. We detected 279 pairs of repeat sequences and found that three pairs of direct repeats could mediate the homologous recombination into four independent circular molecules. We identified 70 SSRs in the whole mitochondrial genome of sweet potato. The longest dispersed repeat in mitochondrial genome was a palindromic repeat with a length of 915 bp. The homologous fragments between the chloroplast and mitochondrial genome account for 7.35% of the mitochondrial genome. We also predicted 597 RNA editing sites and found that the rps3 gene was edited 54 times, which occurred most frequently. This study further demonstrates the existence of multiple conformations in sweet potato mitochondrial genomes and provides a theoretical basis for the evolution of higher plants and cytoplasmic male sterility breeding.

Beyond a reference genome: pangenomes and population genomics of underutilized and orphan crops for future food and nutrition security

Underutilized crops are, by definition, under-researched compared to staple crops yet come with traits that may be especially important given climate change and the need to feed a globally increasing population. These crops are often stress-tolerant, and this combined with unique and beneficial nutritional profiles. Whilst progress is being made by generating reference genome sequences, in this Tansley Review, we show how this is only the very first step. We advocate that going 'beyond a reference genome' should be a priority, as it is only at this stage one can identify the specific genes and the adaptive alleles that underpin the valuable traits. We sum up how population genomic and pangenomic approaches have led to the identification of stress- and disease-tolerant alleles in staple crops and compare this to the small number of examples from underutilized crops. We also demonstrate how previously underutilized crops have benefitted from genomic advances and that many breeding targets in underutilized crops are often well studied in staple crops. This cross-crop population-level resequencing could lead to an understanding of the genetic basis of adaptive traits in underutilized crops. This level of investment may be crucial for fully understanding the value of these crops before they are lost.

A Genome for Bidens hawaiensis: A Member of a Hexaploid Hawaiian Plant Adaptive Radiation

The plant genus Bidens (Asteraceae or Compositae; Coreopsidae) is a species-rich and circumglobally distributed taxon. The 19 hexaploid species endemic to the Hawaiian Islands are considered an iconic example of adaptive radiation, of which many are imperiled and of high conservation concern. Until now, no genomic resources were available for this genus, which may serve as a model system for understanding the evolutionary genomics of explosive plant diversification. Here, we present a high-quality reference genome for the Hawai'i Island endemic species B. hawaiensis A. Gray reconstructed from long-read, high-fidelity sequences generated on a Pacific Biosciences Sequel II System. The haplotype-aware, draft genome assembly consisted of ~6.67 Giga bases (Gb), close to the holoploid genome size estimate of 7.56 Gb (±0.44 SD) determined by flow cytometry. After removal of alternate haplotigs and contaminant filtering, the consensus haploid reference genome was comprised of 15 904 contigs containing ~3.48 Gb, with a contig N50 value of 422 594. The high interspersed repeat content of the genome, approximately 74%, along with hexaploid status, contributed to assembly fragmentation. Both the haplotype-aware and consensus haploid assemblies recovered >96% of Benchmarking Universal Single-Copy Orthologs. Yet, the removal of alternate haplotigs did not substantially reduce the proportion of duplicated benchmarking genes (~79% vs. ~68%). This reference genome will support future work on the speciation process during adaptive radiation, including resolving evolutionary relationships, determining the genomic basis of trait evolution, and supporting ongoing conservation efforts.

Genome-wide identification and characterization of NBS-encoding genes in the sweet potato wild ancestor Ipomoea trifida (H.B.K.)

The most predominant type of resistance (R) genes contain nucleotide-binding sites and leucine-rich repeat (NBS-LRR) domains, characterization of which is helpful for plant resistance improvement. However, the NBS genes of Ipomoea trifida (H.B.K.) remain insufficient to date. In this study, a genome-wide analysis of the NBS-encoding gene in I. trifida (H.B.K.) was carried out. A total of 442 NBS encoding genes were identified, amounting to 1.37% of the total genes of I. trifida (H.B.K.). Based on the analysis of the domains, the identified ItfNBS genes were further classified into seven groups: CNL, NL, CN, N, TNL, TN, and RNL. Phylogenetic analysis showed that the I. trifida NBS genes clustered into three independent clades: RNL, TNL, and CNL. Chromosome location analysis revealed that the distribution of ItfNBS genes in chromosomes was uneven, with a number ranging from 3 to 45. Multiple stress-related regulatory elements were detected in the promoters of the NBS-encoding genes, and their expression profiles were obtained. The qRT-PCR analysis revealed that IbNBS10, IbNBS20, IbNBS258, and IbNBS88 responded to stem nematode infection. These results provide critical proof for further characterization and analysis of NBS-encoding genes with important functions.

Genome-wide survey and expression analysis of GRAS transcription factor family in sweetpotato provides insights into their potential roles in stress response

BACKGROUND

The plant-specific GRAS transcription factors play pivotal roles in various adverse environmental conditions. Numerous GRAS genes have been explored and characterized in different plants, however, comprehensive survey on GRASs in sweetpotato is lagging.

RESULTS

In this study, 72 putative sweetpotato IbGRAS genes with uneven distribution were isolated on 15 chromosomes and classified into 12 subfamilies supported by gene structures and motif compositions. Moreover, both tandem duplication and segmental duplication events played critical roles in the expansion of sweetpotato GRAS genes, and the collinearity between IbGRAS genes and the related orthologs from nine other plants further depicted evolutionary insights into GRAS gene family. RNA-seq analysis under salt stress and qRT-PCR detection of 12 selected IbGRAS genes demonstrated their significant and varying inductions under multiple abiotic stresses (salt, drought, heat and cold) and hormone treatments (ABA, ACC and JA). Consistently, the promoter regions of IbGRAS genes harbored a series of stress- and hormone-associated cis-acting elements. Among them, IbGRAS71, the potential candidate for breeding tolerant plants, was characterized as having transactivation activity in yeasts, while IbGRAS-2/-4/-9 did not. Moreover, a complex interaction relationship between IbGRASs was observed through the interaction network analysis and yeast two-hybrid assays.

CONCLUSIONS

Our results laid a foundation for further functional identifications of IbGRAS genes, and multiple members may serve as potential regulators for molecular breeding of tolerant sweetpotato.

Discovery and characterization of sweetpotato's closest tetraploid relative

The origin of sweetpotato, a hexaploid species, is poorly understood, partly because the identity of its tetraploid progenitor remains unknown. In this study, we identify, describe and characterize a new species of Ipomoea that is sweetpotato's closest tetraploid relative known to date and probably a direct descendant of its tetraploid progenitor. We integrate morphological, phylogenetic, and genomic analyses of herbarium and germplasm accessions of the hexaploid sweetpotato, its closest known diploid relative Ipomoea trifida, and various tetraploid plants closely related to them from across the American continent. We identify wild autotetraploid plants from Ecuador that are morphologically distinct from Ipomoea batatas and I. trifida, but monophyletic and sister to I. batatas in phylogenetic analysis of nuclear data. We describe this new species as Ipomoea aequatoriensis T. Wells & P. Muñoz sp. nov., distinguish it from hybrid tetraploid material collected in Mexico; and show that it likely played a direct role in the origin of sweetpotato's hexaploid genome. This discovery transforms our understanding of sweetpotato's origin.

Sweet Potato (Ipomoea batatas L.) Phenotypes: From Agroindustry to Health Effects

Sweet potato (SP; Ipomoea batatas (L.) Lam) is an edible tuber native to America and the sixth most important food crop worldwide. China leads its production in a global market of USD 45 trillion. SP domesticated varieties differ in specific phenotypic/genotypic traits, yet all of them are rich in sugars, slow digestible/resistant starch, vitamins, minerals, bioactive proteins and lipids, carotenoids, polyphenols, ascorbic acid, alkaloids, coumarins, and saponins, in a genotype-dependent manner. Individually or synergistically, SP’s phytochemicals help to prevent many illnesses, including certain types of cancers and cardiovascular disorders. These and other topics, including the production and market diversification of raw SP and its products, and SP’s starch as a functional ingredient, are briefly discussed in this review.

Genome-Wide Identification and Characterization of CDPK Family Reveal Their Involvements in Growth and Development and Abiotic Stress in Sweet Potato and Its Two Diploid Relatives

Calcium-dependent protein kinase (CDPKs) is one of the calcium-sensing proteins in plants. They are likely to play important roles in growth and development and abiotic stress responses. However, these functions have not been explored in sweet potato. In this study, we identified 39 CDPKs in cultivated hexaploid sweet potato (Ipomoea batatas, 2n = 6x = 90), 35 CDPKs in diploid relative Ipomoea trifida (2n = 2x = 30), and 35 CDPKs in Ipomoea triloba (2n = 2x = 30) via genome structure analysis and phylogenetic characterization, respectively. The protein physiological property, chromosome localization, phylogenetic relationship, gene structure, promoter cis-acting regulatory elements, and protein interaction network were systematically investigated to explore the possible roles of homologous CDPKs in the growth and development and abiotic stress responses of sweet potato. The expression profiles of the identified CDPKs in different tissues and treatments revealed tissue specificity and various expression patterns in sweet potato and its two diploid relatives, supporting the difference in the evolutionary trajectories of hexaploid sweet potato. These results are a critical first step in understanding the functions of sweet potato CDPK genes and provide more candidate genes for improving yield and abiotic stress tolerance in cultivated sweet potato.

Chromosome painting reveals genomic structure of three polyploid species of Ipomoea

Cultivated sweetpotato (Ipomoea batatas (L.) Lam.) from the family Convolvulaceae is a hexaploid species with 2n = 6x = 90, and has been controversial regarding its nature as an autopolyploid arising within a species or allopolyploid forming between species. Here, we developed oligonucleotide-based painting probes for two chromosomes of I.nil, a model diploid Ipomoea species. Using these probes, we revealed pairing behavior of homoeologous chromosomes in I. batatas and its two possible polyploid ancestral species, tetraploid I. tabascana (2n = 4x = 60) and hexaploid I. trifida (2n = 6x = 90). Chromosome painting analysis revealed a high percentage of quadrivalent formation in zygotene-pachytene cells of I. tabascana, which supported that I. tabascana was an autotetraploid likely derived by doubling of structurally similar and homologous genomes rather than a hybrid between I. batatas and I. trifida (2x). The high frequency of hexavalent/bivalent and tetravalent pairing was observed in I. trifida (6x) and I. batatas. However, the percentage of hexavalent pairing in I. trifida (6x) was far higher than that in I. batatas. Thus, present results tended to support that I. trifida (6x) was an autohexaploid, while I. batatas was more likely to be a segmental allohexaploid.

Metabolomic and Transcriptomic Analyses of the Flavonoid Biosynthetic Pathway for the Accumulation of Anthocyanins and Other Flavonoids in Sweetpotato Root Skin and Leaf Vein Base

Sweetpotato [Ipomoea batatas (L.) Lam.] is a major tuberous root crop that is rich in flavonoids. Here, we discovered a spontaneous mutation in the color of the leaf vein base (LVB) and root skin (RS) in the Zheshu 81 cultivar. The flavonoid and anthocyanin metabolites and molecular mechanism were analyzed using metabolome and transcriptome data. Compared to the wild type, 13 differentially accumulated metabolites (DAMs) in the LVB and 59 DAMs in the RS were all significantly downregulated. Moreover, all the anthocyanin metabolites decreased significantly. The differentially expressed genes (DEGs) encoding the key enzymes in the later enzymatic reaction of anthocyanin and flavonoid were significantly downregulated in the mutant. The expression trends of the transcription factor MYB were evidently related to the anthocyanin content. These results offer insights into the coloration in the LVB and RS and a theoretical basis for determining the regulation of flavonoid and anthocyanin synthesis in sweetpotato.

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.

exTREEmaTIME: a method for incorporating uncertainty into divergence time estimates

We present a method of divergence time estimation (exTREEmaTIME) that aims to effectively account for uncertainty in divergence time estimates. The method requires a minimal set of assumptions, and, based on these assumptions, estimates the oldest possible divergence times and youngest possible divergence times that are consistent with the assumptions. We use a series of simulations and empirical analyses to illustrate that exTREEmaTIME is effective at representing uncertainty. We then describe how exTREEmaTIME can act as a basis to determine the implications of the more stringent assumptions that are incorporated into other methods of divergence time estimation that produce more precise estimates. This is critically important given that many of the assumptions that are incorporated into these methods are highly complex, difficult to justify biologically, and as such can lead to estimates that are highly inaccurate. This article has an associated First Person interview with the first author of the paper.

Salinity responses and tolerance mechanisms in underground vegetable crops: an integrative review

The present review gives an insight into the salinity stress tolerance responses and mechanisms of underground vegetable crops. Phytoprotectants, agronomic practices, biofertilizers, and modern biotechnological approaches are crucial for salinity stress management. Underground vegetables are the source of healthy carbohydrates, resistant starch, antioxidants, vitamins, mineral, and nutrients which benefit human health. Soil salinity is a serious threat to agriculture that severely affects the growth, development, and productivity of underground vegetable crops. Salt stress induces several morphological, anatomical, physiological, and biochemical changes in crop plants which include reduction in plant height, leaf area, and biomass. Also, salinity stress impedes the growth of the underground organs, which ultimately reduces crop yield. Moreover, salt stress is detrimental to photosynthesis, membrane integrity, nutrient balance, and leaf water content. Salt tolerance mechanisms involve a complex interplay of several genes, transcription factors, and proteins that are involved in the salinity tolerance mechanism in underground crops. Besides, a coordinated interaction between several phytoprotectants, phytohormones, antioxidants, and microbes is needed. So far, a comprehensive review of salinity tolerance responses and mechanisms in underground vegetables is not available. This review aims to provide a comprehensive view of salt stress effects on underground vegetable crops at different levels of biological organization and discuss the underlying salt tolerance mechanisms. Also, the role of multi-omics in dissecting gene and protein regulatory networks involved in salt tolerance mechanisms is highlighted, which can potentially help in breeding salt-tolerant underground vegetable crops.

The IbBBX24-IbTOE3-IbPRX17 module enhances abiotic stress tolerance by scavenging reactive oxygen species in sweet potato

Soil salinity and drought limit sweet potato yield. Scavenging of reactive oxygen species (ROS) by peroxidases (PRXs) is essential during plant stress responses, but how PRX expression is regulated under abiotic stress is not well understood. Here, we report that the B-box (BBX) family transcription factor IbBBX24 activates the expression of the class III peroxidase gene IbPRX17 by binding to its promoter. Overexpression of IbBBX24 and IbPRX17 significantly improved the tolerance of sweet potato to salt and drought stresses, whereas reducing IbBBX24 expression increased their susceptibility. Under abiotic stress, IbBBX24- and IbPRX17-overexpression lines showed higher peroxidase activity and lower H₂ O₂ accumulation compared with the wild-type. RNA sequencing analysis revealed that IbBBX24 modulates the expression of genes encoding ROS scavenging enzymes, including PRXs. Moreover, interaction between IbBBX24 and the APETALA2 (AP2) protein IbTOE3 enhances the ability of IbBBX24 to activate IbPRX17 transcription. Overexpression of IbTOE3 improved the tolerance of tobacco plants to salt and drought stresses by scavenging ROS. Together, our findings elucidate the mechanism underlying the IbBBX24-IbTOE3-IbPRX17 module in response to abiotic stress in sweet potato and identify candidate genes for developing elite crop varieties with enhanced abiotic stress tolerance.

flopp: Extremely Fast Long-Read Polyploid Haplotype Phasing by Uniform Tree Partitioning

Resolving haplotypes in polyploid genomes using phase information from sequencing reads is an important and challenging problem. We introduce two new mathematical formulations of polyploid haplotype phasing: (1) the min-sum max tree partition problem, which is a more flexible graphical metric compared with the standard minimum error correction (MEC) model in the polyploid setting, and (2) the uniform probabilistic error minimization model, which is a probabilistic analogue of the MEC model. We incorporate both formulations into a long-read based polyploid haplotype phasing method called flopp. We show that flopp compares favorably with state-of-the-art algorithms-up to 30 times faster with 2 times fewer switch errors on 6 × ploidy simulated data. Further, we show using real nanopore data that flopp can quickly reveal reasonable haplotype structures from the autotetraploid Solanum tuberosum (potato).

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

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

Genome-wide comparative analysis of the nucleotide-binding site-encoding genes in four Ipomoea species

The nucleotide-binding site (NBS)-encoding gene is a major type of resistance (R) gene, and its diverse evolutionary patterns were analyzed in different angiosperm lineages. Until now, no comparative studies have been done on the NBS encoding genes in Ipomoea species. In this study, various numbers of NBS-encoding genes were identified across the whole genome of sweet potato (Ipomoea batatas) (#889), Ipomoea trifida (#554), Ipomoea triloba (#571), and Ipomoea nil (#757). Gene analysis showed that the CN-type and N-type were more common than the other types of NBS-encoding genes. The phylogenetic analysis revealed that the NBS-encoding genes formed three monophyletic clades: CNL, TNL, and RNL, which were distinguished by amino acid motifs. The distribution of the NBS-encoding genes among the chromosomes was non-random and uneven; 83.13, 76.71, 90.37, and 86.39% of the genes occurred in clusters in sweet potato, I. trifida, I. triloba, and I. nil, respectively. The duplication pattern analysis reveals the presence of higher segmentally duplicated genes in sweet potatoes than tandemly duplicated ones. The opposite trend was found for the other three species. A total of 201 NBS-encoding orthologous genes were found to form synteny gene pairs between any two of the four Ipomea species, suggesting that each of the synteny gene pairs was derived from a common ancestor. The gene expression patterns were acquired by analyzing using the published datasets. To explore the candidate resistant genes in sweet potato, transcriptome analysis has been carried out using two resistant (JK20 and JK274) and susceptible cultivars (Tengfei and Santiandao) of sweet potato for stem nematodes and Ceratocystis fimbriata pathogen, respectively. A total of 11 differentially expressed genes (DEGs) were found in Tengfei and JK20 for stem nematodes and 19 DEGs in Santiandao and JK274 for C. fimbriata. Moreover, six DEGs were further selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis. The results may provide new insights into the evolution of NBS-encoding genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

Enhancing crop diversity for food security in the face of climate uncertainty

Global agriculture is dominated by a handful of species that currently supply a huge proportion of our food and feed. It additionally faces the massive challenge of providing food for 10 billion people by 2050, despite increasing environmental deterioration. One way to better plan production in the face of current and continuing climate change is to better understand how our domestication of these crops included their adaptation to environments that were highly distinct from those of their centre of origin. There are many prominent examples of this, including the development of temperate Zea mays (maize) and the alteration of day-length requirements in Solanum tuberosum (potato). Despite the pre-eminence of some 15 crops, more than 50 000 species are edible, with 7000 of these considered semi-cultivated. Opportunities afforded by next-generation sequencing technologies alongside other methods, including metabolomics and high-throughput phenotyping, are starting to contribute to a better characterization of a handful of these species. Moreover, the first examples of de novo domestication have appeared, whereby key target genes are modified in a wild species in order to confer predictable traits of agronomic value. Here, we review the scale of the challenge, drawing extensively on the characterization of past agriculture to suggest informed strategies upon which the breeding of future climate-resilient crops can be based.

Isolation and functional characterization of a glucose-6-phosphate/phosphate translocator (IbG6PPT1) from sweet potato (Ipomoea batatas (L.) Lam.)

Sweet potato (Ipomoea batatas (L.) Lam.) is a good source of carbohydrates, an excellent raw material for starch-based industries, and a strong candidate for biofuel production due to its high starch content. However, the molecular basis of starch biosynthesis and accumulation in sweet potato is still insufficiently understood. Glucose-6-phosphate/phosphate translocators (GPTs) mediate the import of glucose-6-phosphate (Glc6P) into plastids for starch synthesis. Here, we report the isolation of a GPT-encoding gene, IbG6PPT1, from sweet potato and the identification of two additional IbG6PPT1 gene copies in the sweet potato genome. IbG6PPT1 encodes a chloroplast membrane-localized GPT belonging to the GPT1 group and highly expressed in storage root of sweet potato. Heterologous expression of IbG6PPT1 resulted in increased starch content in the leaves, root tips, and seeds and soluble sugar in seeds of Arabidopsis thaliana, but a reduction in soluble sugar in the leaves. These findings suggested that IbG6PPT1 might play a critical role in the distribution of carbon sources in source and sink and the accumulation of carbohydrates in storage tissues and would be a good candidate gene for controlling critical starch properties in sweet potato.

Genome-Wide Identification and Analysis of NAC Transcription Factor Family in Two Diploid Wild Relatives of Cultivated Sweet Potato Uncovers Potential NAC Genes Related to Drought Tolerance

NAC (NAM, ATAF1/2, and CUC2) proteins play a pivotal role in modulating plant development and offer protection against biotic and abiotic stresses. Until now, no systematic knowledge of NAC family genes is available for the food security crop, sweet potato. Here, a comprehensive genome-wide survey of NAC domain-containing proteins identified 130 ItbNAC and 144 ItfNAC genes with full length sequences in the genomes of two diploid wild relatives of cultivated sweet potato, Ipomoea triloba and Ipomoea trifida, respectively. These genes were physically mapped onto 15 I. triloba and 16 I. trifida chromosomes, respectively. Phylogenetic analysis divided all 274 NAC proteins into 20 subgroups together with NAC transcription factors (TFs) from Arabidopsis. There were 9 and 15 tandem duplication events in the I. triloba and I. trifida genomes, respectively, indicating an important role of tandem duplication in sweet potato gene expansion and evolution. Moreover, synteny analysis suggested that most NAC genes in the two diploid sweet potato species had a similar origin and evolutionary process. Gene expression patterns based on RNA-Seq data in different tissues and in response to various hormone, biotic or abiotic treatments revealed their possible involvement in organ development and response to various biotic/abiotic stresses. The expression of 36 NAC TFs, which were upregulated in the five tissues and in response to mannitol treatment, was also determined by real-time quantitative polymerase chain reaction (RT-qPCR) in hexaploid cultivated sweet potato exposed to drought stress. Those results largely corroborated the expression profile of mannitol treatment uncovered by the RNA-Seq data. Some significantly up-regulated genes related to drought stress, such as ItbNAC110, ItbNAC114, ItfNAC15, ItfNAC28, and especially ItfNAC62, which had a conservative spatial conformation with a closely related paralogous gene, ANAC019, may be potential candidate genes for a sweet potato drought tolerance breeding program. This analysis provides comprehensive and systematic information about NAC family genes in two diploid wild relatives of cultivated sweet potato, and will provide a blueprint for their functional characterization and exploitation to improve the tolerance of sweet potato to abiotic stresses.

Genome-Wide Identification and Expression Analysis of the MADS-Box Gene Family in Sweet Potato [Ipomoea batatas (L.) Lam]

MADS-box gene, one of the largest transcription factor families in plants, is a class of transcription factors widely present in eukaryotes. It plays an important role in plant growth and development and participates in the growth and development of flowers and fruits. Sweet potato is the seventh most important food crop in the world. Its tuberous roots, stems, and leaves contain a large number of proteins, lipids, carotenoids, anthocyanins, conjugated phenolic acids, and minerals, which have high edible, forage, and medicinal value, and is also an important energy crop. At present, MADS-box genes in sweet potato have rarely been reported, and there has been no study on the genome-wide identification and classification of MADS-box genes in Ipomoea batatas. This study provided the first comprehensive analysis of sweet potato MADS-box genes. We identified 95 MADS-box genes, analyzed the structure and protein of sweet potato MADS-box genes, and categorized them based on phylogenetic analysis with Arabidopsis MADS-box proteins. Chromosomal localization indicated an unequal number of MADS-box genes in all 14 chromosomes except LG3, with more than 10 MADS-box genes located on chromosomes LG7, LG11, and LG15. The MADS domain and core motifs of the sweet potato MADS-box genes were identified by motif analysis. We identified 19 MADS-box genes with collinear relationships and analyzed duplication events. Cis-acting elements, such as light-responsive, auxin-responsive, drought-inducible, and MeJA-responsive elements, were found in the promoter region of the MADS-box genes in sweet potato, which further indicates the basis of MADS-box gene regulation in response to environmental changes and hormones. RNA-seq suggested that sweet potato MADS-box genes exhibit tissue-specific expression patterns, with 34 genes highly expressed in sweet potato flowers and fruits, and 19 genes highly expressed in the tuberous root, pencil root, or fibrous root. qRT-PCR again validated the expression levels of the 10 genes and found that IbMADS1, IbMADS18, IbMADS19, IbMADS79, and IbMADS90 were highly expressed in the tuberous root or fibrous root, and IbMADS18, IbMADS31, and IbMADS83 were highly expressed in the fruit. In this study, the molecular basis of MADS-box genes of sweet potato was analyzed from various angles. The effects of MADS-box genes on the growth and development of sweet potato were investigated, which may provide a certain theoretical basis for molecular breeding of sweet potato.

Genome-Wide Identification and Analysis of NAC Transcription Factor Family in Two Diploid Wild Relatives of Cultivated Sweet Potato Uncovers Potential NAC Genes Related to Drought Tolerance

NAC (NAM, ATAF1/2, and CUC2) proteins play a pivotal role in modulating plant development and offer protection against biotic and abiotic stresses. Until now, no systematic knowledge of NAC family genes is available for the food security crop, sweet potato. Here, a comprehensive genome-wide survey of NAC domain-containing proteins identified 130 ItbNAC and 144 ItfNAC genes with full length sequences in the genomes of two diploid wild relatives of cultivated sweet potato, Ipomoea triloba and Ipomoea trifida, respectively. These genes were physically mapped onto 15 I. triloba and 16 I. trifida chromosomes, respectively. Phylogenetic analysis divided all 274 NAC proteins into 20 subgroups together with NAC transcription factors (TFs) from Arabidopsis. There were 9 and 15 tandem duplication events in the I. triloba and I. trifida genomes, respectively, indicating an important role of tandem duplication in sweet potato gene expansion and evolution. Moreover, synteny analysis suggested that most NAC genes in the two diploid sweet potato species had a similar origin and evolutionary process. Gene expression patterns based on RNA-Seq data in different tissues and in response to various hormone, biotic or abiotic treatments revealed their possible involvement in organ development and response to various biotic/abiotic stresses. The expression of 36 NAC TFs, which were upregulated in the five tissues and in response to mannitol treatment, was also determined by real-time quantitative polymerase chain reaction (RT-qPCR) in hexaploid cultivated sweet potato exposed to drought stress. Those results largely corroborated the expression profile of mannitol treatment uncovered by the RNA-Seq data. Some significantly up-regulated genes related to drought stress, such as ItbNAC110, ItbNAC114, ItfNAC15, ItfNAC28, and especially ItfNAC62, which had a conservative spatial conformation with a closely related paralogous gene, ANAC019, may be potential candidate genes for a sweet potato drought tolerance breeding program. This analysis provides comprehensive and systematic information about NAC family genes in two diploid wild relatives of cultivated sweet potato, and will provide a blueprint for their functional characterization and exploitation to improve the tolerance of sweet potato to abiotic stresses.

Dynamic network biomarker analysis discovers IbNAC083 in the initiation and regulation of sweet potato root tuberization

The initiation and development of storage roots (SRs) are intricately regulated by a transcriptional regulatory network. One key challenge is to accurately pinpoint the tipping point during the transition from pre-swelling to SRs and to identify the core regulators governing such a critical transition. To solve this problem, we performed a dynamic network biomarker (DNB) analysis of transcriptomic dynamics during root development in Ipomoea batatas (sweet potato). First, our analysis identified stage-specific expression patterns for a significant proportion (>9%) of the sweet potato genes and unraveled the chronology of events that happen at the early and later stages of root development. Then, the results showed that different root developmental stages can be depicted by co-expressed modules of sweet potato genes. Moreover, we identified the key components and transcriptional regulatory network that determine root development. Furthermore, through DNB analysis an early stage, with a root diameter of 3.5 mm, was identified as the critical period of SR swelling initiation, which is consistent with morphological and metabolic changes. In particular, we identified a NAM/ATAF/CUC (NAC) domain transcription factor, IbNAC083, as a core regulator of this initiation in the DNB-associated network. Further analyses and experiments showed that IbNAC083, along with its associated differentially expressed genes, induced dysfunction of metabolism processes, including the biosynthesis of lignin, flavonol and starch, thus leading to the transition to swelling roots.

Genetic Identity, Diversity, and Population Structure of CIP's Sweetpotato (I. batatas) Germplasm Collection

The in trust sweetpotato collection housed by the International Center of Potato (CIP) is one of the largest assemblages of plant material representing the genetic resources of this important staple crop. The collection currently contains almost 6,000 accessions of Ipomoea batatas (cultivated sweetpotato) and over 1,000 accessions of sweetpotato crop wild relatives (CWRs). In this study, the entire cultivated collection (5,979 accessions) was genotyped with a panel of 20 simple sequence repeat (SSR) markers to assess genetic identity, diversity, and population structure. Genotyping and phenotyping of in vitro plantlets and mother plants were conducted simultaneously on 2,711 accessions (45% of the total collection) to identify and correct possible genetic identity errors which could have occurred at any time over the thirty plus years of maintenance in the in vitro collection. Within this group, 533 accessions (19.6%) had errors in identity. Field evaluations of morphological descriptors were carried out to confirm the marker data. A phylogenetic tree was constructed to reveal the intraspecific relationships in the population which uncovered high levels of redundancy in material from Peru and Latin America. These genotypic data were supported by morphological data. Population structure analysis demonstrated support for four ancestral populations with many of the accessions having lower levels of gene flow from the other populations. This was especially true of germplasm derived from Peru, Ecuador, and Africa. The set of 20 SSR markers was subsequently utilized to examine a subset of 189 accessions from the USDA sweetpotato germplasm collection and to identify and reconcile potential errors in the identification of clones shared between these collections. Marker analysis demonstrated that the USDA subset of material had 65 unique accessions that were not found in the larger CIP collection. As far as the authors are aware, this is the first report of genotyping an entire sweetpotato germplasm collection in its entirety.

A single amino acid mutant in the EAR motif of IbMYB44.2 reduced the inhibition of anthocyanin accumulation in the purple-fleshed sweetpotato

Purple-fleshed sweetpotato (Ipomoea batatas（L.）Lam.) is rich in anthocyanins. R2R3-type MYB transcription factors（TFs）with EAR motifs inhibiting anthocyanin biosynthesis have been reported, and there is still a lack of information on how mutations in the EAR motifs of MYBs affect anthocyanin accumulation. In this study, we obtained three IbMYB44 TFs by bioinformatics. Among these TFs, IbMYB44.1, IbMYB44.3 with a complete EAR motif and IbMYB44.2 with a single amino acid mutant in the EAR motif caused an amino acid substitution from leucine to valine. RT-qPCR analysis showed that IbMYB44s was expressed at lower levels in the purple-fleshed sweetpotato than in nonpurple-fleshed sweetpotato (P < 0.01). Transient expression assays showed that the inhibitory effect of IbMYB44.1/3 was stronger than IbMYB44.2 in tobacco leaves and red-skinned pears. RT-qPCR analysis further proved that IbMYB44.1/3 significantly inhibited the expression of anthocyanin biosynthesis-related genes compared with IbMYB44.2 in tobacco leaves and red-skinned pears. A dual luciferase reporter assay showed that IbMYB44s cannot directly activate the IbANS promoter, and the result was also verified by yeast one-hybrid (Y1H) experiments. Moreover, we identified the interaction of IbMYB340 with IbMYB44.1, IbMYB44.2 and IbMYB44.3 via yeast two-hybrid (Y2H) assays. Thus, IbMYB44.1/3 could interact with IbMYB340 to negatively regulate anthocyanin biosynthesis. This study enriched the regulatory network of anthocyanins and also provided a theoretical basis for a single amino acid mutant from leucine to valine in the EAR motif of IbMYB44.2 affecting anthocyanin biosynthesis in the purple-fleshed sweetpotato.

Meiosis in Polyploids and Implications for Genetic Mapping: A Review

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

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.

Genome-Wide Identification and Expression Analysis of JAZ Family Involved in Hormone and Abiotic Stress in Sweet Potato and Its Two Diploid Relatives

Jasmonate ZIM-domain (JAZ) proteins are key repressors of a jasmonic acid signaling pathway. They play essential roles in the regulation of plant growth and development, as well as environmental stress responses. However, this gene family has not been explored in sweet potato. In this study, we identified 14, 15, and 14 JAZs in cultivated hexaploid sweet potato (Ipomoea batatas, 2n = 6x = 90), and its two diploid relatives Ipomoea trifida (2n = 2x = 30) and Ipomoea triloba (2n = 2x = 30), respectively. These JAZs were divided into five subgroups according to their phylogenetic relationships with Arabidopsis. The protein physiological properties, chromosome localization, phylogenetic relationship, gene structure, promoter cis-elements, protein interaction network, and expression pattern of these 43 JAZs were systematically investigated. The results suggested that there was a differentiation between homologous JAZs, and each JAZ gene played different vital roles in growth and development, hormone crosstalk, and abiotic stress response between sweet potato and its two diploid relatives. Our work provided comprehensive comparison and understanding of the JAZ genes in sweet potato and its two diploid relatives, supplied a theoretical foundation for their functional study, and further facilitated the molecular breeding of sweet potato.