Genome architecture and tetrasomic inheritance of autotetraploid potato

Bridging the gap: unravelling plant centromeres in the telomere-to-telomere era

Centromeres are specific regions of the chromosomes that play a pivotal role in the segregation of chromosomes, by facilitating the loading of the kinetochore, which forms the link between the chromosomes to the spindle fibres during cell division. In plants and animals, these regions often form megabase-scale loci of tandemly repeated DNA sequences, which have presented a challenge to genomic studies even in model species. The functional designation of centromeres is determined epigenetically by the incorporation of a centromere-specific variant of histone H3. Recent developments in long-read sequencing technology have allowed the assembly of these regions for the first time and have prompted a reassessment of fidelity of centromere function and the evolutionary dynamics of these regions.

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

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

Genome assembly in the telomere-to-telomere era

Genome sequences largely determine the biology and encode the history of an organism, and de novo assembly - the process of reconstructing the genome sequence of an organism from sequencing reads - has been a central problem in bioinformatics for four decades. Until recently, genomes were typically assembled into fragments of a few megabases at best, but now technological advances in long-read sequencing enable the near-complete assembly of each chromosome - also known as telomere-to-telomere assembly - for many organisms. Here, we review recent progress on assembly algorithms and protocols, with a focus on how to derive near-telomere-to-telomere assemblies. We also discuss the additional developments that will be required to resolve remaining assembly gaps and to assemble non-diploid genomes.

Potato: from functional genomics to genetic improvement

Potato is the most widely grown non-grain crop and ranks as the third most significant global food crop following rice and wheat. Despite its long history of cultivation over vast areas, slow breeding progress and environmental stress have led to a scarcity of high-yielding potato varieties. Enhancing the quality and yield of potato tubers remains the ultimate objective of potato breeding. However, conventional breeding has faced challenges due to tetrasomic inheritance, high genomic heterozygosity, and inbreeding depression. Recent advancements in molecular biology and functional genomic studies of potato have provided valuable insights into the regulatory network of physiological processes and facilitated trait improvement. In this review, we present a summary of identified factors and genes governing potato growth and development, along with progress in potato genomics and the adoption of new breeding technologies for improvement. Additionally, we explore the opportunities and challenges in potato improvement, offering insights into future avenues for potato research.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

Chromosome-level scaffolding of haplotype-resolved assemblies using Hi-C data without reference genomes

Scaffolding is crucial for constructing most chromosome-level genomes. The high-throughput chromatin conformation capture (Hi-C) technology has become the primary scaffolding strategy due to its convenience and cost-effectiveness. As sequencing technologies and assembly algorithms advance, constructing haplotype-resolved genomes is increasingly preferred because haplotypes can provide additional genetic information on allelic and non-allelic variations. ALLHiC is a widely used allele-aware scaffolding tool designed for this purpose. However, its dependence on chromosome-level reference genomes and a higher chromosome misassignment rate still impede the unravelling of haplotype-resolved genomes. Here we present HapHiC, a reference-independent allele-aware scaffolding tool with superior performance on chromosome assignment as well as contig ordering and orientation. In addition, we provide new insights into the challenges in allele-aware scaffolding by conducting comprehensive analyses on various adverse factors. Finally, with the help of HapHiC, we constructed the haplotype-resolved allotriploid genome for Miscanthus × giganteus, an important lignocellulosic bioenergy crop.

Haplotype-resolved genome assembly for tetraploid Chinese cherry (Prunus pseudocerasus) offers insights into fruit firmness

Chinese cherry (Prunus pseudocerasus) holds considerable importance as one of the primary stone fruit crops in China. However, artificially improving its traits and genetic analysis are challenging due to lack of high-quality genomic resources, which mainly result from difficulties associated with resolving its tetraploid and highly heterozygous genome. Herein, we assembled a chromosome-level, haplotype-resolved genome of the cultivar 'Zhuji Duanbing', comprising 993.69 Mb assembled into 32 pseudochromosomes using PacBio HiFi, Oxford Nanopore, and Hi-C. Intra-haplotype comparative analyses revealed extensive intra-genomic sequence and expression consistency. Phylogenetic and comparative genomic analyses demonstrated that P. pseudocerasus was a stable autotetraploid species, closely related to wild P. pusilliflora, with the two diverging ~18.34 million years ago. Similar to other Prunus species, P. pseudocerasus underwent a common whole-genome duplication event that occurred ~139.96 million years ago. Because of its low fruit firmness, P. pseudocerasus is unsuitable for long-distance transportation, thereby restricting its rapid development throughout China. At the ripe fruit stage, P. pseudocerasus cv. 'Zhuji Duanbing' was significantly less firm than P. avium cv. 'Heizhenzhu'. The difference in firmness is attributed to the degree of alteration in pectin, cellulose, and hemicellulose contents. In addition, comparative transcriptomic analyses identified GalAK-like and Stv1, two genes involved in pectin biosynthesis, which potentially caused the difference in firmness between 'Zhuji Duanbing' and 'Heizhenzhu'. Transient transformations of PpsGalAK-like and PpsStv1 increase protopectin content and thereby enhance fruit firmness. Our study lays a solid foundation for functional genomic studies and the enhancement of important horticultural traits in Chinese cherries.

Haplotype-resolved genomes of octoploid species in Phyllanthaceae family reveal a critical role for polyploidization and hybridization in speciation

The Phyllanthaceae family comprises a diverse range of plants with medicinal, edible, and ornamental value, extensively cultivated worldwide. Polyploid species commonly occur in Phyllanthaceae. Due to the rather complex genomes and evolutionary histories, their speciation process has been still lacking in research. In this study, we generated chromosome-scale haplotype-resolved genomes of two octoploid species (Phyllanthus emblica and Sauropus spatulifolius) in Phyllanthaceae family. Combined with our previously reported one tetraploid (Sauropus androgynus) and one diploid species (Phyllanthus cochinchinensis) from the same family, we explored their speciation history. The three polyploid species were all identified as allopolyploids with subgenome A/B. Each of their two distinct subgenome groups from various species was uncovered to independently share a common diploid ancestor (Ancestor-AA and Ancestor-BB). Via different evolutionary routes, comprising various scenarios of bifurcating divergence, allopolyploidization (hybrid polyploidization), and autopolyploidization, they finally evolved to the current tetraploid S. androgynus, and octoploid S. spatulifolius and P. emblica, respectively. We further discuss the variations in copy number of alleles and the potential impacts within the two octoploids. In addition, we also investigated the fluctuation of metabolites with medical values and identified the key factor in its biosynthesis process in octoploids species. Our study reconstructed the evolutionary history of these Phyllanthaceae species, highlighting the critical roles of polyploidization and hybridization in their speciation processes. The high-quality genomes of the two octoploid species provide valuable genomic resources for further research of evolution and functional genomics.

Origin and evolution of a new tetraploid mangrove species in an intertidal zone

Polyploidy is a major factor in the evolution of plants, yet we know little about the origin and evolution of polyploidy in intertidal species. This study aimed to identify the evolutionary transitions in three true-mangrove species of the genus Acanthus distributed in the Indo-West Pacific region. For this purpose, we took an integrative approach that combined data on morphology, cytology, climatic niche, phylogeny, and biogeography of 493 samples from 42 geographic sites. Our results show that the Acanthus ilicifolius lineage distributed east of the Thai-Malay Peninsula possesses a tetraploid karyotype, which is morphologically distinct from that of the lineage on the west side. The haplotype networks and phylogenetic trees for the chloroplast genome and eight nuclear genes reveal that the tetraploid species has two sub-genomes, one each from A. ilicifolius and A . ebracteatus, the paternal and maternal parents, respectively. Population structure analysis also supports the hybrid speciation history of the new tetraploid species. The two sub-genomes of the tetraploid species diverged from their diploid progenitors during the Pleistocene. Environmental niche models revealed that the tetraploid species not only occupied the near-entire niche space of the diploids, but also expanded into novel environments. Our findings suggest that A. ilicifolius species distributed on the east side of the Thai-Malay Peninsula should be regarded as a new species, A. tetraploideus, which originated from hybridization between A. ilicifolius and A. ebracteatus, followed by chromosome doubling. This is the first report of a true-mangrove allopolyploid species that can reproduce sexually and clonally reproduction, which explains the long-term adaptive potential of the species.

The complex hexaploid oil-Camellia genome traces back its phylogenomic history and multi-omics analysis of Camellia oil biosynthesis

Oil-Camellia (Camellia oleifera), belonging to the Theaceae family Camellia, is an important woody edible oil tree species. The Camellia oil in its mature seed kernels, mainly consists of more than 90% unsaturated fatty acids, tea polyphenols, flavonoids, squalene and other active substances, which is one of the best quality edible vegetable oils in the world. However, genetic research and molecular breeding on oil-Camellia are challenging due to its complex genetic background. Here, we successfully report a chromosome-scale genome assembly for a hexaploid oil-Camellia cultivar Changlin40. This assembly contains 8.80 Gb genomic sequences with scaffold N50 of 180.0 Mb and 45 pseudochromosomes comprising 15 homologous groups with three members each, which contain 135 868 genes with an average length of 3936 bp. Referring to the diploid genome, intragenomic and intergenomic comparisons of synteny indicate homologous chromosomal similarity and changes. Moreover, comparative and evolutionary analyses reveal three rounds of whole-genome duplication (WGD) events, as well as the possible diversification of hexaploid Changlin40 with diploid occurred approximately 9.06 million years ago (MYA). Furthermore, through the combination of genomics, transcriptomics and metabolomics approaches, a complex regulatory network was constructed and allows to identify potential key structural genes (SAD, FAD2 and FAD3) and transcription factors (AP2 and C2H2) that regulate the metabolism of Camellia oil, especially for unsaturated fatty acids biosynthesis. Overall, the genomic resource generated from this study has great potential to accelerate the research for the molecular biology and genetic improvement of hexaploid oil-Camellia, as well as to understand polyploid genome evolution.

Scalable telomere-to-telomere assembly for diploid and polyploid genomes with double graph

Despite advances in long-read sequencing technologies, constructing a near telomere-to-telomere assembly is still computationally demanding. Here we present hifiasm (UL), an efficient de novo assembly algorithm combining multiple sequencing technologies to scale up population-wide near telomere-to-telomere assemblies. Applied to 22 human and two plant genomes, our algorithm produces better diploid assemblies at a cost of an order of magnitude lower than existing methods, and it also works with polyploid genomes.

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

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

Genome-wide identification and expression analysis of the SET domain-containing gene family in potato (Solanum tuberosum L.)

Genes containing the SET domain can catalyse histone lysine methylation, which in turn has the potential to cause changes to chromatin structure and regulation of the transcription of genes involved in diverse physiological and developmental processes. However, the functions of SET domain-containing (StSET) genes in potato still need to be studied. The objectives of our study can be summarized as in silico analysis to (i) identify StSET genes in the potato genome, (ii) systematically analyse gene structure, chromosomal distribution, gene duplication events, promoter sequences, and protein domains, (iii) perform phylogenetic analyses, (iv) compare the SET domain-containing genes of potato with other plant species with respect to protein domains and orthologous relationships, (v) analyse tissue-specific expression, and (vi) study the expression of StSET genes in response to drought and heat stresses. In this study, we identified 57 StSET genes in the potato genome, and the genes were physically mapped onto eleven chromosomes. The phylogenetic analysis grouped these StSET genes into six clades. We found that tandem duplication through sub-functionalisation has contributed only marginally to the expansion of the StSET gene family. The protein domain TDBD (PFAM ID: PF16135) was detected in StSET genes of potato while it was absent in all other previously studied species. This study described three pollen-specific StSET genes in the potato genome. Expression analysis of four StSET genes under heat and drought in three potato clones revealed that these genes might have non-overlapping roles under different abiotic stress conditions and durations. The present study provides a comprehensive analysis of StSET genes in potatoes, and it serves as a basis for further functional characterisation of StSET genes towards understanding their underpinning biological mechanisms in conferring stress tolerance.

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

KEY MESSAGE

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

Genome-wide identification of GATA transcription factors in tetraploid potato and expression analysis in differently colored potato flesh

The GATA gene family belongs to a kind of transcriptional regulatory protein featuring a zinc finger motif, which is essential for plant growth and development. However, the identification of the GATA gene family in tetraploid potato is still not performed. In the present research, a total of 88 GATA genes in the tetraploid potato C88.v1 genome were identified by bioinformatics methods. These StGATA genes had an uneven distribution on 44 chromosomes, and the corresponding StGATA proteins were divided into four subfamilies (I-IV) based on phylogenetic analysis. The cis-elements of StGATA genes were identified, including multiple cis-elements related to light-responsive and hormone-responsive. The collinearity analysis indicates that segmental duplication is a key driving force for the expansion of GATA gene family in tetraploid potato, and that the GATA gene families of tetraploid potato and Arabidopsis share a closer evolutionary relationship than rice. The transcript profiling analysis showed that all 88 StGATA genes had tissue-specific expression, indicating that the StGATA gene family members participate in the development of multiple potato tissues. The RNA-seq analysis was also performed on the tuber flesh of two potato varieties with different color, and 18 differentially expressed GATA transcription factor genes were screened, of which eight genes were validated through qRT-PCR. In this study, we identified and characterized StGATA transcription factors in tetraploid potato for the first time, and screened differentially expressed genes in potato flesh with different color. It provides a theoretical basis for further understanding the StGATA gene family and its function in anthocyanin biosynthesis.

The structure, function, and evolution of plant centromeres

Centromeres are essential regions of eukaryotic chromosomes responsible for the formation of kinetochore complexes, which connect to spindle microtubules during cell division. Notably, although centromeres maintain a conserved function in chromosome segregation, the underlying DNA sequences are diverse both within and between species and are predominantly repetitive in nature. The repeat content of centromeres includes high-copy tandem repeats (satellites), and/or specific families of transposons. The functional region of the centromere is defined by loading of a specific histone 3 variant (CENH3), which nucleates the kinetochore and shows dynamic regulation. In many plants, the centromeres are composed of satellite repeat arrays that are densely DNA methylated and invaded by centrophilic retrotransposons. In some cases, the retrotransposons become the sites of CENH3 loading. We review the structure of plant centromeres, including monocentric, holocentric, and metapolycentric architectures, which vary in the number and distribution of kinetochore attachment sites along chromosomes. We discuss how variation in CENH3 loading can drive genome elimination during early cell divisions of plant embryogenesis. We review how epigenetic state may influence centromere identity and discuss evolutionary models that seek to explain the paradoxically rapid change of centromere sequences observed across species, including the potential roles of recombination. We outline putative modes of selection that could act within the centromeres, as well as the role of repeats in driving cycles of centromere evolution. Although our primary focus is on plant genomes, we draw comparisons with animal and fungal centromeres to derive a eukaryote-wide perspective of centromere structure and function.

Comprehensive assessment of 11 de novo HiFi assemblers on complex eukaryotic genomes and metagenomes

Pacific Biosciences (PacBio) HiFi sequencing technology generates long reads (>10 kbp) with very high accuracy (<0.01% sequencing error). Although several de novo assembly tools are available for HiFi reads, there are no comprehensive studies on the evaluation of these assemblers. We evaluated the performance of 11 de novo HiFi assemblers on (1) real data for three eukaryotic genomes; (2) 34 synthetic data sets with different ploidy, sequencing coverage levels, heterozygosity rates, and sequencing error rates; (3) one real metagenomic data set; and (4) five synthetic metagenomic data sets with different composition abundance and heterozygosity rates. The 11 assemblers were evaluated using quality assessment tool (QUAST) and benchmarking universal single-copy ortholog (BUSCO). We also used several additional criteria, namely, completion rate, single-copy completion rate, duplicated completion rate, average proportion of largest category, average distance difference, quality value, run-time, and memory utilization. Results show that hifiasm and hifiasm-meta should be the first choice for assembling eukaryotic genomes and metagenomes with HiFi data. We performed a comprehensive benchmarking study of commonly used assemblers on complex eukaryotic genomes and metagenomes. Our study will help the research community to choose the most appropriate assembler for their data and identify possible improvements in assembly algorithms.

Haplotype-resolved chromosome-level genome of hexaploid Jerusalem artichoke provides insights into its origin, evolution, and inulin metabolism

Jerusalem artichoke (Helianthus tuberosus) is a global multifunctional crop. It has wide applications in the food, health, feed, and biofuel industries and in ecological protection; it also serves as a germplasm pool for breeding of the global oil crop common sunflower (Helianthus annuus). However, biological studies of Jerusalem artichoke have been hindered by a lack of genome sequences, and its high polyploidy and large genome size have posed challenges to genome assembly. Here, we report a 21-Gb chromosome-level assembly of the hexaploid Jerusalem artichoke genome, which comprises 17 homologous groups, each with 6 pseudochromosomes. We found multiple large-scale chromosome rearrangements between Jerusalem artichoke and common sunflower, and our results show that the hexaploid genome of Jerusalem artichoke was formed by a hybridization event between a tetraploid and a diploid Helianthus species, followed by chromosome doubling of the hybrid, which occurred approximately 2 million years ago. Moreover, we identified more copies of actively expressed genes involved in inulin metabolism and showed that these genes may still be undergoing loss of function or sub- or neofunctionalization. These genomic resources will promote further biological studies, breeding improvement, and industrial utilization of Helianthus crops.

iFLAS: positive-unlabeled learning facilitates full-length transcriptome-based identification and functional exploration of alternatively spliced isoforms in maize

The advent of full-length transcriptome sequencing technologies has accelerated the discovery of novel splicing isoforms. However, existing alternative splicing (AS) tools are either tailored for short-read RNA-Seq data or designed for human and animal studies. The disparities in AS patterns between plants and animals still pose a challenge to the reliable identification and functional exploration of novel isoforms in plants. Here, we developed integrated full-length alternative splicing analysis (iFLAS), a plant-optimized AS toolkit that introduced a semi-supervised machine learning method known as positive-unlabeled (PU) learning to accurately identify novel isoforms. iFLAS also enables the investigation of AS functions from various perspectives, such as differential AS, poly(A) tail length, and allele-specific AS (ASAS) analyses. By applying iFLAS to three full-length transcriptome sequencing datasets, we systematically identified and functionally characterized maize (Zea mays) AS patterns. We found intron retention not only introduces premature termination codons, resulting in lower expression levels of isoforms, but may also regulate the length of 3'UTR and poly(A) tail, thereby affecting the functional differentiation of isoforms. Moreover, we observed distinct ASAS patterns in two genes within heterosis offspring, highlighting their potential value in breeding. These results underscore the broad applicability of iFLAS in plant full-length transcriptome-based AS research.

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

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

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

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

Screening for Resistance Resources against Bacterial Wilt in Wild Potato

Potato is an important crop, used not only for food production but also for various industrial applications. With the introduction of the potato as a staple food strategy, the potato industry in China has grown rapidly. However, issues related to bacterial wilt, exacerbated by factors such as seed potato transportation and continuous cropping, have become increasingly severe in the primary potato cultivation regions of China, leading to significant economic losses. The extensive genetic diversity of Ralstonia solanacearum (R. solanacearum), which is the pathogen of bacterial wilt, has led to a lack of highly resistant potato genetic resources. There is a need to identify and cultivate potato varieties with enhanced resistance to reduce the adverse impact of this disease on the industry. We screened 55 accessions of nine different wild potato species against the bacterial wilt pathogen R. solanacearum PO2-1, which was isolated from native potato plants and belongs to phylotype II. Three accessions of two species (ACL24-2, PNT880-3, and PNT204-23) were identified with high resistance phenotypes to the tested strains. We found these accessions also showed high resistance to different phylotype strains. Among them, only PNT880-3 was capable of flowering and possessed viable pollen, and it was diploid. Consistent with the high resistance, decreased growth of R. solanacearum was detected in PNT880-3. All these findings in our study reveal that the wild potato PNT880-3 was a valuable resistance source to bacterial wilt with breeding potential.

The haplotype-resolved genome assembly of autotetraploid rhubarb Rheum officinale provides insights into its genome evolution and massive accumulation of anthraquinones

Rheum officinale, a member of the Polygonaceae family, is an important medicinal plant that is widely used in traditional Chinese medicine. Here, we report a 7.68-Gb chromosome-scale assembly of R. officinale with a contig N50 of 3.47 Mb, which was clustered into 44 chromosomes across four homologous groups. Comparative genomics analysis revealed that transposable elements have made a significant contribution to its genome evolution, gene copy number variation, and gene regulation and expression, particularly of genes involved in metabolite biosynthesis, stress resistance, and root development. We placed the recent autotetraploidization of R. officinale at ∼0.58 mya and analyzed the genomic features of its homologous chromosomes. Although no dominant monoploid genomes were observed at the overall expression level, numerous allele-differentially-expressed genes were identified, mainly with different transposable element insertions in their regulatory regions, suggesting that they functionally diverged after polyploidization. Combining genomics, transcriptomics, and metabolomics, we explored the contributions of gene family amplification and tetraploidization to the abundant anthraquinone production of R. officinale, as well as gene expression patterns and differences in anthraquinone content among tissues. Our report offers unprecedented genomic resources for fundamental research on the autopolyploid herb R. officinale and guidance for polyploid breeding of herbs.

A chromosome-level genome assembly for Onobrychis viciifolia reveals gene copy number gain underlying enhanced proanthocyanidin biosynthesis

Sainfoin (Onobrychis viciifolia), which belongs to subfamily Papilionoideae of Leguminosae, is a vital perennial forage known as "holy hay" due to its high contents of crude proteins and proanthocyanidins (PAs, also called condensed tannins) that have various pharmacological properties in animal feed, such as alleviating rumen tympanic disease in ruminants. In this study, we select an autotetraploid common sainfoin (2n = 4x = 28) and report its high-quality chromosome-level genome assembly with 28 pseudochromosomes and four haplotypes (~1950.14 Mb, contig N50 = 10.91 Mb). The copy numbers of genes involved in PA biosynthesis in sainfoin are significantly greater than those in four selected Fabales species, namely, autotetraploid Medicago sativa and three other diploid species, Lotus japonicus, Medicago truncatula, and Glycine max. Furthermore, gene expansion is confirmed to be the key contributor to the increased expression of these genes and subsequent PA enhancement in sainfoin. Transcriptomic analyses reveal that the expression of genes involved in the PA biosynthesis pathway is significantly increased in the lines with high PA content compared to the lines with medium and low PA content. The sainfoin genome assembly will improve our understanding of leguminous genome evolution and biosynthesis of secondary metabolites in sainfoin.

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

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

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

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

Dynamic viral integration patterns actively participate in the progression of BK polyomavirus-associated diseases after renal transplantation

The classical lytic infection theory along with large T antigen-mediated oncogenesis cannot explain the BK polyomavirus (BKPyV)-associated tumor secondary to BKPyV-associated nephropathy (BKVAN), viremia/DNAemia, and viruria after renal transplantation. This study performed virome capture sequencing and pathological examination on regularly collected urine sediment and peripheral blood samples, and BKVAN and tumor biopsy tissues of 20 patients with BKPyV-associated diseases of different stages. In the early noncancerous stages, well-amplified integration sites were visualized by in situ polymerase chain reaction, simultaneously with BKPyV inclusion bodies and capsid protein expression. The integration intensity, the proportion of microhomology-mediated end-joining integration, and host PARP-1 and POLQ gene expression levels increased with disease progression. Furthermore, multiomics analysis was performed on BKPyV-associated urothelial carcinoma tissues, identifying tandem-like structures of BKPyV integration using long-read genome sequencing. The carcinogenicity of BKPyV integration was proven to disturb host gene expression and increase viral oncoprotein expression. Fallible DNA double-strand break repair pathways were significantly activated in the parenchyma of BKPyV-associated tumors. Olaparib showed an antitumor activity dose-response effect in the tumor organoids without BRCA1/2 genes mutation. In conclusion, the dynamic viral integration patterns actively participate in the progression of BKPyV-associated diseases and thus could be a potential target for disease monitoring and intervention.

RNA-seq provides insights into potato deubiquitinase responses to drought stress in seedling stage

Ubiquitination is a specific protein degradation and reversible post-translational modification process that can be reversed by deubiquitinase (DUBs). DUBs can hydrolyze and release ubiquitin in the substrate protein so that the substrate can avoid degradation or change its activity, and it has an impact on plant growth and development, cell cycle, abiotic stress response, and other biological processes. Transcript sequences of potato varieties "DM1-3", "Atlantic" and "Cooperation-88" downloaded from Potato Genome Resources were used for genome-wide identification of the DUB gene family using Hidden Markov Models and verified in the NCBI CD-Search tool. The characteristics of DUB genes from different potato varieties were analyzed including subcellular localization, gene structural motifs, phylogenetic tree, and sequence homology. Polyethylene glycol 6000 (PEG6000) induced drought stress transcriptome analysis was performed on the "Atlantic", and differentially expressed genes were screened, with emphasis on the characterization of deubiquitinase. DUB genes have a complex gene structure, often with a large number of exons and alternative splicing. Their promoters contain abundant abiotic stress-responsive elements, such as 425 MYC, 325 ABRE, and 320 MYB. There are also a large number of orthologous genes in the DUBs of the three potato varieties, and these genes are often clustered in similar regions on the genome. We performed transcriptome sequencing of the potato under PEG-induced drought stress and analyzed it for the first time using the Atlantic as a reference genome. We identified a total of 6067 down-regulated differentially expressed genes (DEGs) and 4950 up-regulated DEGs under PEG-induced drought stress. We screened the expression of DUBs and observed that 120 DUBs were up-regulated where most of them functioned in the nucleus, and the interacting proteins of DUBs were also localized in the nucleus. We have comprehensively identified and analyzed potato DUBs, and the accurately aligned transcriptome data which will further deepen the understanding of DUBs involved in the regulation of osmotic stress.

Recent Advances in Assembly of Complex Plant Genomes

Over the past 20 years, tremendous advances in sequencing technologies and computational algorithms have spurred plant genomic research into a thriving era with hundreds of genomes decoded already, ranging from those of nonvascular plants to those of flowering plants. However, complex plant genome assembly is still challenging and remains difficult to fully resolve with conventional sequencing and assembly methods due to high heterozygosity, highly repetitive sequences, or high ploidy characteristics of complex genomes. Herein, we summarize the challenges of and advances in complex plant genome assembly, including feasible experimental strategies, upgrades to sequencing technology, existing assembly methods, and different phasing algorithms. Moreover, we list actual cases of complex genome projects for readers to refer to and draw upon to solve future problems related to complex genomes. Finally, we expect that the accurate, gapless, telomere-to-telomere, and fully phased assembly of complex plant genomes could soon become routine.

The Expression of the StNRAMP2 Gene Determined the Accumulation of Cadmium in Different Tissues of Potato

Cadmium (Cd) is a toxic metal that threatens human health when enriched in crops. NRAMPs are a family of natural macrophage proteins reported to play a key role in Cd transport in plants. In order to explore the gene regulation mechanism of potato under Cd stress and the role of NRAMPs family in it, this study analyzed the gene expression differences of two different Cd accumulation levels in potato after 7 days of 50 mg/kg Cd stress and screened out the key genes that may play a major role in the differential accumulation of Cd in different varieties. Additionally, StNRAMP2 was selected for verification. Further verification showed that the StNRAMP2 gene plays an important role in the accumulation of Cd in potato. Interestingly, silencing StNRAMP2 increased Cd accumulation in tubers but significantly decreased Cd accumulation in other sites, suggesting a critical role of StNRAMP2 in Cd uptake and transport in potatoes. To further confirm this conclusion, we performed heterologous expression experiments in which overexpression of StNRAMP2 gene in tomato resulted in a threefold increase in Cd content, which further confirmed the important role of StNRAMP2 in the process of Cd accumulation compared with wild-type plants. In addition, we found that the addition of Cd to the soil increased the activity of the plant antioxidant enzyme system, and silencing StNRAMP2 partially reversed this effect. This suggests that the StNRAMP2 gene plays an important role in plant stress tolerance, and future studies could further explore the role of this gene in other environmental stresses. In conclusion, the results of this study improve the understanding of the mechanism of Cd accumulation in potato and provide experimental basis for remediation of Cd pollution.

Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding

Hybrid potato breeding will transform the crop from a clonally propagated tetraploid to a seed-reproducing diploid. Historical accumulation of deleterious mutations in potato genomes has hindered the development of elite inbred lines and hybrids. Utilizing a whole-genome phylogeny of 92 Solanaceae and its sister clade species, we employ an evolutionary strategy to identify deleterious mutations. The deep phylogeny reveals the genome-wide landscape of highly constrained sites, comprising ∼2.4% of the genome. Based on a diploid potato diversity panel, we infer 367,499 deleterious variants, of which 50% occur at non-coding and 15% at synonymous sites. Counterintuitively, diploid lines with relatively high homozygous deleterious burden can be better starting material for inbred-line development, despite showing less vigorous growth. Inclusion of inferred deleterious mutations increases genomic-prediction accuracy for yield by 24.7%. Our study generates insights into the genome-wide incidence and properties of deleterious mutations and their far-reaching consequences for breeding.

Insights into the genetic determination of tuber shape and eye depth in potato natural population based on autotetraploid potato genome

Potato is one of the world's most important food crops, with a time-consuming breeding process. In this study, we performed a genome-wide association (GWAS) analysis of the two important traits of potato tuber shape and eye depth, using the tetraploid potato genome (2n=4x=48) as a reference. A total of 370 potatoes were divided into three subgroups based on the principal component analysis and evolutionary tree analysis. The genetic diversity within subgroups is low (5.18×10^-5, 4.36×10^-5 and 4.24×10^-5). Genome-wide linkage disequilibrium (LD) analysis showed that their LD is about 60 Kb. GWAS analysis identified that 146 significant single nucleotide polymorphism (SNP) loci at Chr01A1:34.44-35.25 Mb and Chr02A1:28.35-28.54 Mb regions are significantly associated with potato tuber shape, and that three candidate genes that might be related to potato tuber traits, PLATZ transcription factor, UTP-glucose-1-phosphate uridylyltransferase and FAR1 DNA-binding domain, are in the association region of Chr02A1. GWAS analysis identified 53 significant SNP loci at Chr05A2: 49.644-50.146 Mb and Chr06A2: 25.866-26.384 Mb regions with robust associations with potato tuber eye depth. Hydrolase and methyltransferases are present in the association region of Chr05A2, and three CYPs are present in the association region of Chr06A2. Our findings suggested that these genes are closely associated with potato tuber shape and eye depth. Our study identified molecular markers and candidate genes for improving tetraploid potato tuber shape and eye depth and provided ideas and insights for tetraploid potato breeding.

Insights into the Genetic Determination of the Autotetraploid Potato Plant Height

Plant height is an important characteristic, the modification of which can improve the ability of stress adaptation as well as the yield. In this study, genome-wide association analysis was performed for plant height traits in 370 potato cultivars using the tetraploid potato genome as a reference. A total of 92 significant single nucleotide polymorphism (SNP) loci for plant height were obtained, which were particularly significant in haplotypes A3 and A4 on chromosome 1 and A1, A2, and A4 on chromosome 5. Thirty-five candidate genes were identified that were mainly involved in the gibberellin and brassinolide signal transduction pathways, including the FAR1 gene, methyltransferase, ethylene response factor, and ubiquitin protein ligase. Among them, PIF3 and GID1a were only present on chromosome 1, with PIF3 in all four haplotypes and GID1a in haplotype A3. This could lead to more effective genetic loci for molecular marker-assisted selection breeding as well as more precise localization and cloning of genes for plant height traits in potatoes.

Genomic and genetic advances of oiltea-camellia (Camellia oleifera)

Oiltea-camellia (C. oleifera) is a widely cultivated woody oil crop in Southern China and Southeast Asia. The genome of oiltea-camellia was very complex and not well explored. Recently, genomes of three oiltea-camellia species were sequenced and assembled, multi-omic studies of oiltea-camellia were carried out and provided a better understanding of this important woody oil crop. In this review, we summarized the recent assembly of the reference genomes of oiltea-camellia, genes related to economic traits (flowering, photosynthesis, yield and oil component), disease resistance (anthracnose) and environmental stress tolerances (drought, cold, heat and nutrient deficiency). We also discussed future directions of integrating multiple omics for evaluating genetic resources and mining key genes of important traits, and the application of new molecular breeding and gene editing technologies to accelerate the breeding process of oiltea-camellia.