Analysis of evolution and genetic diversity of sweetpotato and its related different polyploidy wild species I. trifida using RAD-seq

Genetic and morphological evidence support the specific status of the endemic Ericaandevalensis (Ericales, Ericaceae)

Assessing the taxonomic status of closely related taxa is crucial in plant systematics and can have important implications for conservation and human plant use. Ericaandevalensis Cabezudo & Rivera is a metallophyte endemic species from highly metal-polluted soils of SW Iberian Peninsula, an area with a mining history going back more than 5,000 years. Ericaandevalensis is closely related to Ericamackayana Bab., a northern Iberian species also present in western Ireland. The status of E.andevalensis as a species or subspecies subordinated to E.mackayana is subject to debate. Here, we assessed the genetic and phenotypic relationship between both species, including the population structure of E.andevalensis. We used high throughput sequencing to determine genome-wide Single Nucleotide Polymorphisms (SNPs), and morphometric analyses from 35 reproductive and vegetative traits. The morphological analysis showed at least eight characters that can discriminate the two species, from which ovary hairiness and the size of leaf glandular hairs were the most informative. Genetic analyses showed that each species formed a monophyletic cluster with full support, separated by an interspecific genetic distance >4-fold higher than intra-specific distance. Population genetic analyses of E.andevalensis shows that populations are highly structured, with the Portuguese one as the most isolated and less variable. These results support the recognition of E.andevalensis as a distinct species with a highly constrained ecological requirements and a narrow geographic distribution, but with a limited gene flow between populations. We discuss the implications of these outcomes in conservation policies and potential uses of E.andevalensis such as decontamination of polluted soils.

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.

Random mutagenesis in vegetatively propagated crops: opportunities, challenges and genome editing prospects

In order to meet the growing human food and nutrition demand a perpetual process of crop improvement is idealized. It has seen changing trends and varying concepts throughout human history; from simple selection to complex gene-editing. Among these techniques, random mutagenesis has been shown to be a promising technology to achieve desirable genetic gain with less time and minimal efforts. Over the decade, several hundred varieties have been released through random mutagenesis, but the production is falling behind the demand. Several food crops like banana, potato, cassava, sweet potato, apple, citrus, and others are vegetatively propagated. Since such crops are not propagated through seed, genetic improvement through classical breeding is impractical for them. Besides, in the case of polyploids, accomplishment of allelic homozygosity requires a considerable land area, extensive fieldwork with huge manpower, and hefty funding for an extended period of time. Apart from induction, mapping of induced genes to facilitate the knowledge of biological processes has been performed only in a few selected facultative vegetative crops like banana and cassava which can form a segregating population. During the last few decades, there has been a shift in the techniques used for crop improvement. With the introduction of the robust technologies like meganucleases, zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and clustered regularly interspaced short palindromic repeats (CRISPR) more and more crops are being subjected to gene editing. However, more work needs to be done in case of vegetatively propagated crops.

Plastid genome sequencing, identification of nuclear SNP markers, and quality assessment of medicinal rhizomatous herb Polygonatum odoratum (Asparagaceae) cultivars

Polygonatum odoratum (Mill.) Druce (Asparagaceae, Asparagales) is a widely cultivated medicinal herb in China. However, this useful herb is understudied despite being known as a medicinal resource with top grade medical and edible properties since long. In this study, P. odoratum and four cultivars were investigated. The variations in morphological characteristics and vegetative phases of each cultivar were observed. For genetic aspect, the plastid genome of P. odoratum varies in length from 154,569 bp to 155,491 bp, containing a large single-copy region of 83,486-84,459 bp, a small single-copy region of 18,292-18,471 bp, and two inverted repeats of 26,302-26,370 bp. A total of 131 genes were predicted, including 85 protein-coding, 38 tRNA, and eight rRNA genes. Genome comparisons revealed a slight variation in the sequence across the five accessions, but two highly variable regions (trnC-petN and rpl32-trnL) were detected when comparing the four different cultivars. For the RAD-seq markers, a total of 33.64 Gb of clean data, with an average value of 1.08 Gb per sample, were analyzed for the presence of single nucleotide polymorphisms (SNPs). Well-resolved phylogenies of the P. odoratum cultivars are constructed; the nonmonophyletic relationship in the plastome-based phylogenetic trees, yet monophyletic form in the RAD-based linkage map suggested possibility of hybrid cultivar for P. odoratum "Dazhu" (GDDZ), which was further supported by morphological observations. Quality assessment based on the standards of the Chinese Pharmacopoeia on Polygonati Odorati Rhizoma (POR) on the four cultivars used in this study recorded that PORs from P. odoratum "Zhongzhu" (GDZZ) met the minimum criteria for the acceptance as raw material for medicinal drug production. This study has provided insights on the morphological variations, genetic background, and medicinal qualities of P. odoratum cultivars that could be explored for future genetic improvement as well as breeding programs of P. odoratum for POR production.

Genome-wide characterization and expression analysis of HAK K+ transport family in Ipomoea

The potassium transporter high-affinity K⁺ transporter/K⁺ uptake permease/K⁺ transporter (HAK/KUP/KT) family plays a vital role in potassium uptake, and potassium ion (K⁺)-mediated environmental stress. In the present study, we identified 22 IbHAK/KUP/KT (HAK) genes in sweet potato [Ipomoea batata (L.) Lam] and the same number of HAK genes from sweet potato wild relative Ipomoea trifida. Phylogeny analysis indicated that the HAKs can be divided into five clades. Chromosomal distribution and genome synteny analyses revealed two tandem-duplicated gene pairs IbHAK16/17 and IbHAK17/18 on chromosomes 13 and eight segmental-duplicated gene pairs on chromosomes 1, 3, 5, 8, 10, 12, 14 among the IbHAK gene family. Eleven orthologous HAK gene pairs between I. batata and I. trifida were involved in the duplication of genomic blocks based on comparative genomic analysis. The Ka/Ks ratios of these IbHAK genes ranged from 0.02 to 0.55(< 1), further indicated that purifying selection was the primary force driving the evolution of HAKs in Ipomoea. A heat map based on RNA-seq data showed that 13 HAKs in Xushu32 (a K⁺-tolerant sweet potato genotype) and 10 HAKs in Ningzi1 (a K⁺-sensitive sweet potato genotype) in response to K⁺ deficiency stress. Quantitative real-time PCR (qRT-PCR) analysis revealed IbHAK2, -3, -8, -10, -11, -18, -19, and -21 were induced in both Xushu32 and Ningzi1 under low K⁺ stress. Compared with other IbHAK genes, IbHAK8 showed more strongly upregulation after exposure to drought and salt stress. Furthermore, co-expression analysis showed that only IbHAK8 of 22 IbHAK genes involved in network interactions with 30 genes related to abiotic and biotic stresses. Taken together, these results are helpful for further functional studies on IbHAK and molecular breeding of sweet potato.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-020-02552-3.

Phylogenomic Relationships and Evolution of Polyploid Salix Species Revealed by RAD Sequencing Data

Polyploidy is common in the genus Salix. However, little is known about the origin, parentage and genomic composition of polyploid species because of a lack of suitable molecular markers and analysis tools. We established a phylogenomic framework including species of all described sections of Eurasian shrub willows. We analyzed the genomic composition of seven polyploid willow species in comparison to putative diploid parental species to draw conclusions on their origin and the effects of backcrossing and post-origin evolution. We applied recently developed programs like SNAPP, HyDe, and SNiPloid to establish a bioinformatic pipeline for unravelling the complexity of polyploid genomes. RAD sequencing revealed 23,393 loci and 320,010 high quality SNPs for the analysis of relationships of 35 species of Eurasian shrub willows (Salix subg. Chamaetia/Vetrix). Polyploid willow species appear to be predominantly of allopolyploid origin. More ancient allopolyploidization events were observed for two hexaploid and one octoploid species, while our data suggested a more recent allopolyploid origin for the included tetraploids and identified putative parental taxa. SNiPloid analyses disentangled the different genomic signatures resulting from hybrid origin, backcrossing, and secondary post-origin evolution in the polyploid species. Our RAD sequencing data demonstrate that willow genomes are shaped by ancient and recent reticulate evolution, polyploidization, and post-origin divergence of species.

Genotyping by RAD Sequencing Analysis Assessed the Genetic Distinctiveness of Experimental Lines and Narrowed Down the Genomic Region Responsible for Leaf Shape in Endive (Cichorium endivia L.)

The characterization of genetic diversity in elite breeding stocks is crucial for the registration and protection of new varieties. Moreover, experimental population structure analysis and information about the genetic distinctiveness of commercial materials are essential for crop breeding programs. The purpose of our research was to assess the genetic relationships of 32 endive (Cichorium endivia L.) breeding lines, 18 from var. latifolium (escarole) and 14 from var. crispum (curly), using heterologous Cichorium intybus-derived simple sequence repeats (SSR) markers and single-nucleotide polymorphisms (SNP) markers. We found that 14 out of 29 SSR markers were successfully amplified, but only 8 of them were related to polymorphic loci. To overcome the limitation of the low number of informative SSR marker loci, an alternative SNP-based approach was employed. The 4621 SNPs produced by a restriction site-associated DNA marker sequencing approach were able to fully discriminate the 32 endive accessions; most importantly, as many as 50 marker loci were found to distinguish the curly group from the escarole group. Interestingly, 24 of the marker loci mapped within a peripheral segment of chromosome 8 of lettuce (Lactuca sativa L.), spanning a chromosomal region of 49.6 Mb. Following Sanger sequencing-based validation, three genes were determined to carry nonsynonymous SNPs, and one of them matched a putative ortholog of AtELP1, subunit 1 of the Elongator complex. Considering that several previously characterized Elongator complex subunit mutants exhibited elongated and/or curly leaf phenotypes, this gene should be taken into consideration for a better understanding of the underlying mechanism controlling leaf shape in endive.

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

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

The Molecular Determination of Hybridity and Homozygosity Estimates in Breeding Populations of Lettuce (Lactuca sativa L.)

The development of new varieties of horticultural crops benefits from the integration of conventional and molecular marker-assisted breeding schemes in order to combine phenotyping and genotyping information. In this study, a selected panel of 16 microsatellite markers were used in different steps of a breeding programme of lettuce (Lactuca sativa L., 2 n = 18). Molecular markers were first used to genotype 71 putative parental lines and to plan 89 controlled crosses designed to maximise recombination potentials. The resulting 871 progeny plants were then molecularly screened, and their marker allele profiles were compared with the profiles expected based on the parental lines. The average cross-pollination success rate was 68 ± 33%, so 602 F1 hybrids were completely identified. Unexpected genotypes were detected in 5% of cases, consistent with this species' spontaneous out-pollination rate. Finally, in a later step of the breeding programme, 47 different F3 progenies, selected by phenotyping for a number of morphological descriptors, were characterised in terms of their observed homozygosity and within-population genetic uniformity and stability. Ten of these populations had a median homozygosity above 90% and a median genetic similarity above 95% and are, therefore, particularly suitable for pre-commercial trials. In conclusion, this study shows the synergistic effects and advantages of conventional and molecular methods of selection applied in different steps of a breeding programme aimed at developing new varieties of lettuce.

The wild sweetpotato (Ipomoea trifida) genome provides insights into storage root development

BACKGROUND

Sweetpotato (Ipomoea batatas (L.) Lam.) is the seventh most important crop in the world and is mainly cultivated for its underground storage root (SR). The genetic studies of this species have been hindered by a lack of high-quality reference sequence due to its complex genome structure. Diploid Ipomoea trifida is the closest relative and putative progenitor of sweetpotato, which is considered a model species for sweetpotato, including genetic, cytological, and physiological analyses.

RESULTS

Here, we generated the chromosome-scale genome sequence of SR-forming diploid I. trifida var. Y22 with high heterozygosity (2.20%). Although the chromosome-based synteny analysis revealed that the I. trifida shared conserved karyotype with Ipomoea nil after the separation, I. trifida had a much smaller genome than I. nil due to more efficient eliminations of LTR-retrotransposons and lack of species-specific amplification bursts of LTR-RTs. A comparison with four non-SR-forming species showed that the evolution of the beta-amylase gene family may be related to SR formation. We further investigated the relationship of the key gene BMY11 (with identity 47.12% to beta-amylase 1) with this important agronomic trait by both gene expression profiling and quantitative trait locus (QTL) mapping. And combining SR morphology and structure, gene expression profiling and qPCR results, we deduced that the products of the activity of BMY11 in splitting starch granules and be recycled to synthesize larger granules, contributing to starch accumulation and SR swelling. Moreover, we found the expression pattern of BMY11, sporamin proteins and the key genes involved in carbohydrate metabolism and stele lignification were similar to that of sweetpotato during the SR development.

CONCLUSIONS

We constructed the high-quality genome reference of the highly heterozygous I. trifida through a combined approach and this genome enables a better resolution of the genomics feature and genome evolutions of this species. Sweetpotato SR development genes can be identified in I. trifida and these genes perform similar functions and patterns, showed that the diploid I. trifida var. Y22 with typical SR could be considered an ideal model for the studies of sweetpotato SR development.

Root-zone-specific sensitivity of K+-and Ca2+-permeable channels to H2O2 determines ion homeostasis in salinized diploid and hexaploid Ipomoea trifida

Polyploids generally possess superior K+/Na+ homeostasis under saline conditions compared with their diploid progenitors. In this study, we identified the physiological mechanisms involved in the ploidy-related mediation of K+/Na+ homeostasis in the roots of diploid (2x) and hexaploid (6x; autohexaploid) Ipomoea trifida, which is the closest relative of cultivated sweet potato. Results showed that 6x I. trifida retained more K+ and accumulated less Na+ in the root and leaf tissues under salt stress than 2x I. trifida. Compared with its 2x ancestor, 6x I. trifida efficiently prevents K+ efflux from the meristem root zone under salt stress through its plasma membrane (PM) K+-permeable channels, which have low sensitivity to H2O2. Moreover, 6x I. trifida efficiently excludes Na+ from the elongation and mature root zones under salt stress because of the high sensitivity of PM Ca2+-permeable channels to H2O2. Our results suggest the root-zone-specific sensitivity to H2O2 of PM K+- and Ca2+-permeable channels in the co-ordinated control of K+/Na+ homeostasis in salinized 2x and 6x I. trifida. This work provides new insights into the improved maintenance of K+/Na+ homeostasis of polyploids under salt stress.

In silico prediction and segregation analysis of putative virus defense genes based on SSR markers in sweet potato F1 progenies of cultivars 'New Kawogo' and 'Resisto'

In sweet potato, an anti-virus defense mechanism termed reversion has been postulated to lead to virus freedom from once infected plants. The objectives of this study were to identify anti-virus defense genes and evaluate their segregation in progenies. Reference genes from different plant species were used to assemble transcript sequences of each sweet potato defense gene in silico. Sequences were used for evaluate phylogenetic relationships with similar genes from different plant species, mining respective defense genes and thereafter developing simple sequence repeats (SSRs) for segregation analysis. Eight potential defense genes were identified: RNA dependent RNA polymerases 1, 2, 5, and 6; Argonaute 1, and Dicer-like 1, 2, and 4. Identified genes were differentially related to those of other plants and were observed on different chromosomes. The defense genes contained mono-, di-, tri-, tetra, penta-, and hexa-nucleotide repeat motifs. The SSR markers within progenies were segregated in disomic, co-segregation, nullisomic, monosomic, and trisomic modes. These findings indicate the possibility of deriving and utilizing SSRs using published genomic information. Furthermore, and given that the SSR markers were derived from known genes on defined chromosomes, this work will contribute to future molecular breeding and development of resistance gene analogs in this economically important crop.