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.

The integrated genomics of crop domestication and breeding.

As a major event in human civilization, wild plants were successfully domesticated to be crops, largely owing to continuing artificial selection. Here, we summarize new discoveries made during the past decade in crop domestication and breeding. The construction of crop genome maps and the functional characterization of numerous trait genes provide foundational information. Approaches to read, interpret, and write complex genetic information are being leveraged in many plants for highly efficient de novo or re-domestication. Understanding the underlying mechanisms of crop microevolution and applying the knowledge to agricultural productions will give possible solutions for future challenges in food security.

Genome design of hybrid potato.

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

Rewiring of the Fruit Metabolome in Tomato Breeding.

Humans heavily rely on dozens of domesticated plant species that have been further improved through intensive breeding. To evaluate how breeding changed the tomato fruit metabolome, we have generated and analyzed a dataset encompassing genomes, transcriptomes, and metabolomes from hundreds of tomato genotypes. The combined results illustrate how breeding globally altered fruit metabolite content. Selection for alleles of genes associated with larger fruits altered metabolite profiles as a consequence of linkage with nearby genes. Selection of five major loci reduced the accumulation of anti-nutritional steroidal glycoalkaloids in ripened fruits, rendering the fruit more edible. Breeding for pink tomatoes modified the content of over 100 metabolites. The introgression of resistance genes from wild relatives in cultivars also resulted in major and unexpected metabolic changes. The study reveals a multi-omics view of the metabolic breeding history of tomato, as well as provides insights into metabolome-assisted breeding and plant biology.

Genome evolution and diversity of wild and cultivated potatoes.

Potato (Solanum tuberosum L.) is the world's most important non-cereal food crop, and the vast majority of commercially grown cultivars are highly heterozygous tetraploids. Advances in diploid hybrid breeding based on true seeds have the potential to revolutionize future potato breeding and production1-4. So far, relatively few studies have examined the genome evolution and diversity of wild and cultivated landrace potatoes, which limits the application of their diversity in potato breeding. Here we assemble 44 high-quality diploid potato genomes from 24 wild and 20 cultivated accessions that are representative of Solanum section Petota, the tuber-bearing clade, as well as 2 genomes from the neighbouring section, Etuberosum. Extensive discordance of phylogenomic relationships suggests the complexity of potato evolution. We find that the potato genome substantially expanded its repertoire of disease-resistance genes when compared with closely related seed-propagated solanaceous crops, indicative of the effect of tuber-based propagation strategies on the evolution of the potato genome. We discover a transcription factor that determines tuber identity and interacts with the mobile tuberization inductive signal SP6A. We also identify 561,433 high-confidence structural variants and construct a map of large inversions, which provides insights for improving inbred lines and precluding potential linkage drag, as exemplified by a 5.8-Mb inversion that is associated with carotenoid content in tubers. This study will accelerate hybrid potato breeding and enrich our understanding of the evolution and biology of potato as a global staple food crop.

Graph pangenome captures missing heritability and empowers tomato breeding.

Missing heritability in genome-wide association studies defines a major problem in genetic analyses of complex biological traits1,2. The solution to this problem is to identify all causal genetic variants and to measure their individual contributions3,4. Here we report a graph pangenome of tomato constructed by precisely cataloguing more than 19 million variants from 838 genomes, including 32 new reference-level genome assemblies. This graph pangenome was used for genome-wide association study analyses and heritability estimation of 20,323 gene-expression and metabolite traits. The average estimated trait heritability is 0.41 compared with 0.33 when using the single linear reference genome. This 24% increase in estimated heritability is largely due to resolving incomplete linkage disequilibrium through the inclusion of additional causal structural variants identified using the graph pangenome. Moreover, by resolving allelic and locus heterogeneity, structural variants improve the power to identify genetic factors underlying agronomically important traits leading to, for example, the identification of two new genes potentially contributing to soluble solid content. The newly identified structural variants will facilitate genetic improvement of tomato through both marker-assisted selection and genomic selection. Our study advances the understanding of the heritability of complex traits and demonstrates the power of the graph pangenome in crop breeding.

Building pangenome graphs.

Pangenome graphs can represent all variation between multiple reference genomes, but current approaches to build them exclude complex sequences or are based upon a single reference. In response, we developed the PanGenome Graph Builder, a pipeline for constructing pangenome graphs without bias or exclusion. The PanGenome Graph Builder uses all-to-all alignments to build a variation graph in which we can identify variation, measure conservation, detect recombination events and infer phylogenetic relationships.

Engineering a universal and efficient platform for terpenoid synthesis in yeast.

Engineering microbes for the production of valuable natural products is often hindered by the regulation of native competing metabolic networks in host. This is particularly evident in the case of terpenoid synthesis in yeast, where the canonical terpenoid precursors are tightly coupled to the biosynthesis of sterols essential for yeast viability. One way to circumvent this limitation is by engineering product pathways less connected to the host native metabolism. Here, we introduce a two-step isopentenol utilization pathway (IUP) in Saccharomyces cerevisiae to augment the native mevalonate pathway by providing a shortcut to the synthesis of the common terpenoid precursors, isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP). As such, the IUP was capable of elevating the IPP/DMAPP pool by 147-fold compared with the native pathway. We further demonstrate that cofeeding isoprenol and prenol enhances geranyl diphosphate (GPP) content for monoterpene biosynthesis. More importantly, we established a synthetic three-step route for efficient synthesis of di-and tetraterpene precursor geranylgeranyl diphosphate (GGPP), circumventing the competition with farnesyl diphosphate (FPP) for sterol biosynthesis and elevating the GGPP level by 374-fold. We combine these IUP-supported precursor-forming platforms with downstream terpene synthases to harness their potential and improve the production of industrially relevant terpenoids by several fold. Our exploration provides a universal and effective platform for supporting terpenoid synthesis in yeast.

Biosynthesis of bridged tricyclic sesquiterpenes in Inula lineariifolia.

Presilphiperfolane-type sesquiterpenes represent a unique group of atypical sesquiterpenoids characterized by their distinctive tricyclic structure. They have significant potential as lead compounds for pharmaceutical and agrochemical development. Herein, we utilized a transcriptomic approach to identify a terpene synthase (TPS) gene responsible for the biosynthesis of rare presilphiperfolane-type sesquiterpenes in Inula lineariifolia, designated as IlTPS1. Through phylogenetic analysis, we have identified the evolutionary conservation of key motifs, including RR(x)8W, DDxxD, and NSE/DTE in IlTPS1, which are shared with other tricyclic sesquiterpene synthases in the TPS-a subfamily of Asteraceae plants. Subsequent biochemical characterization of recombinant IlTPS1 revealed it to be a multiproduct enzyme responsible for the synthesis of various tricyclic sesquiterpene alcohols from farnesyl diphosphate (FPP), resulting in production of seven distinct sesquiterpenes. Mass spectrometry and nuclear magnetic resonance (NMR) spectrometry identified presilphiperfolan-8ß-ol and presilphiperfol-7-ene as predominant products. Furthermore, biological activity assays revealed that the products from IlTPS1 exhibited a potent antifungal activity against Nigrospora oryzae. Our study represents a significant advancement as it not only functionally identifies the first step enzyme in presilphiperfolane biosynthesis but also establishes the heterologous bioproduction of these unique sesquiterpenes.

Association Analysis Provides Insights into Plant Mitonuclear Interactions.

Cytonuclear interaction refers to the complex and ongoing process of coevolution between nuclear and organelle genomes, which are responsible for cellular respiration, photosynthesis, lipid metabolism, etc. and play a significant role in adaptation and speciation. There have been a large number of studies to detect signatures of cytonuclear interactions. However, identification of the specific nuclear and organelle genetic polymorphisms that are involved in these interactions within a species remains relatively rare. The recent surge in whole genome sequencing has provided us an opportunity to explore cytonuclear interaction from a population perspective. In this study, we analyzed a total of 3,439 genomes from 7 species to identify signals of cytonuclear interactions by association (linkage disequilibrium) analysis of variants in both the mitochondrial and nuclear genomes across flowering plants. We also investigated examples of nuclear loci identified based on these association signals using subcellular localization assays, gene editing, and transcriptome sequencing. Our study provides a novel perspective on the investigation of cytonuclear coevolution, thereby enriching our understanding of plant fitness and offspring sterility.

A modified Agrobacterium-mediated transformation for two oomycete pathogens.

Oomycetes are a group of filamentous microorganisms that include some of the biggest threats to food security and natural ecosystems. However, much of the molecular basis of the pathogenesis and the development in these organisms remains to be learned, largely due to shortage of efficient genetic manipulation methods. In this study, we developed modified transformation methods for two important oomycete species, Phytophthora infestans and Plasmopara viticola, that bring destructive damage in agricultural production. As part of the study, we established an improved Agrobacterium-mediated transformation (AMT) method by prokaryotic expression in Agrobacterium tumefaciens of AtVIP1 (VirE2-interacting protein 1), an Arabidopsis bZIP gene required for AMT but absent in oomycetes genomes. Using the new method, we achieved an increment in transformation efficiency in two P. infestans strains. We further obtained a positive GFP transformant of P. viticola using the modified AMT method. By combining this method with the CRISPR/Cas12a genome editing system, we successfully performed targeted mutagenesis and generated loss-of-function mutations in two P. infestans genes. We edited a MADS-box transcription factor-encoding gene and found that a homozygous mutation in MADS-box results in poor sporulation and significantly reduced virulence. Meanwhile, a single-copy avirulence effector-encoding gene Avr8 in P. infestans was targeted and the edited transformants were virulent on potato carrying the cognate resistance gene R8, suggesting that loss of Avr8 led to successful evasion of the host immune response by the pathogen. In summary, this study reports on a modified genetic transformation and genome editing system, providing a potential tool for accelerating molecular genetic studies not only in oomycetes, but also other microorganisms.

The genomic route to tomato breeding: Past, present, and future.

Over the past 10,000 years, tomato species have undergone both unintentional and intentional selection to enhance their favorable traits for human consumption and manufacturing. These selection processes have significantly influenced the genomes of tomato species and have played a critical role in improving tomato varieties. In this review, we summarize recent advances in tomato genome sequencing, explore the impact of human-driven selection, and recapitulate key genes associated with important agronomic traits in tomato breeding. We provide several examples of genomics-guided tomato breeding to highlight the potential of genome resources in facilitating tomato improvement. Furthermore, we elaborate the progress and strategies of tomato breeding through genome design and present how such efforts can help future enhancements of tomato to align with the demands of sustainability and evolving human societies.

Genomes of autotetraploid wild and cultivated Ziziphus mauritiana reveal polyploid evolution and crop domestication.

Indian jujube (Ziziphus mauritiana) holds a prominent position in the global fruit and pharmaceutical markets. Here, we report the assemblies of haplotype-resolved, telomere-to-telomere genomes of autotetraploid wild and cultivated Indian jujube plants using a two-stage assembly strategy. The generation of these genomes permitted in-depth investigations into the divergence and evolutionary history of this important fruit crop. Using a graph-based pan-genome constructed from eight monoploid genomes, we identified structural variation (SV)-FST hotspots and SV hotspots. Gap-free genomes provide a means to obtain a global view of centromere structures. We identified presence-absence variation-related genes in four monoploid genomes (cI, cIII, wI, and wIII) and resequencing populations. We also present the population structure and domestication trajectory of the Indian jujube based on the resequencing of 73 wild and cultivated accessions. Metabolomic and transcriptomic analyses of mature fruits of wild and cultivated accessions unveiled the genetic basis underlying loss of fruit astringency during domestication of Indian jujube. This study reveals mechanisms underlying the divergence, evolution, and domestication of the autotetraploid Indian jujube and provides rich and reliable genetic resources for future research.

CuGenDBv2: an updated database for cucurbit genomics.

The Cucurbitaceae (cucurbit) family consists of about 1,000 species in 95 genera, including many economically important and popular fruit and vegetable crops. During the past several years, reference genomes have been generated for >20 cucurbit species, and variome and transcriptome profiling data have been rapidly accumulated for cucurbits. To efficiently mine, analyze and disseminate these large-scale datasets, we have developed an updated version of Cucurbit Genomics Database. The updated database, CuGenDBv2 (http://cucurbitgenomics.org/v2), currently hosts 34 reference genomes from 27 cucurbit species/subspecies belonging to 10 different genera. Protein-coding genes from these genomes have been comprehensively annotated by comparing their protein sequences to various public protein and domain databases. A novel 'Genotype' module has been implemented to facilitate mining and analysis of the functionally annotated variome data including SNPs and small indels from large-scale genome sequencing projects. An updated 'Expression' module has been developed to provide a comprehensive gene expression atlas for cucurbits. Furthermore, synteny blocks between any two and within each of the 34 genomes, representing a total of 595 pair-wise genome comparisons, have been identified and can be explored and visualized in the database.

Loss of salt tolerance during tomato domestication conferred by variation in a Na+ /K+ transporter.

Domestication has resulted in reduced salt tolerance in tomato. To identify the genetic components causing this deficiency, we performed a genome-wide association study (GWAS) for root Na+ /K+ ratio in a population consisting of 369 tomato accessions with large natural variations. The most significant variations associated with root Na+ /K+ ratio were identified within the gene SlHAK20 encoding a member of the clade IV HAK/KUP/KT transporters. We further found that SlHAK20 transports Na+ and K+ and regulates Na+ and K+ homeostasis under salt stress conditions. A variation in the coding sequence of SlHAK20 was found to be the causative variant associated with Na+ /K+ ratio and confer salt tolerance in tomato. Knockout mutations in tomato SlHAK20 and the rice homologous genes resulted in hypersensitivity to salt stress. Together, our study uncovered a previously unknown molecular mechanism of salt tolerance responsible for the deficiency in salt tolerance in cultivated tomato varieties. Our findings provide critical information for molecular breeding to improve salt tolerance in tomato and other crops.

Genomic analyses provide insights into peach local adaptation and responses to climate change.

The environment has constantly shaped plant genomes, but the genetic bases underlying how plants adapt to environmental influences remain largely unknown. We constructed a high-density genomic variation map of 263 geographically representative peach landraces and wild relatives. A combination of whole-genome selection scans and genome-wide environmental association studies (GWEAS) was performed to reveal the genomic bases of peach adaptation to diverse climates. A total of 2092 selective sweeps that underlie local adaptation to both mild and extreme climates were identified, including 339 sweeps conferring genomic pattern of adaptation to high altitudes. Using genome-wide environmental association studies (GWEAS), a total of 2755 genomic loci strongly associated with 51 specific environmental variables were detected. The molecular mechanism underlying adaptive evolution of high drought, strong UVB, cold hardiness, sugar content, flesh color, and bloom date were revealed. Finally, based on 30 yr of observation, a candidate gene associated with bloom date advance, representing peach responses to global warming, was identified. Collectively, our study provides insights into molecular bases of how environments have shaped peach genomes by natural selection and adds candidate genes for future studies on evolutionary genetics, adaptation to climate changes, and breeding.

Comprehensive dissection of meiotic DNA double-strand breaks and crossovers in cucumber.

Meiotic recombination drives genetic diversity and crop genome optimization. In plant breeding, parents with favorable traits are crossed to create elite varieties. Different hybridizations produce diverse types of segment reshuffling between homologous chromosomes. However, little is known about the factors that cause hybrid-specific changes in crossovers (COs). Here, we constructed 2 F2 populations from crosses between a semiwild and 2 domesticated cucumber (Cucumis sativus) accessions and examined CO events. COs mainly occurred around genes and differed unevenly along chromosomes between the 2 hybrids. Fine-scale CO distributions were suppressed in regions of heterozygous structural variations (SVs) and were accelerated by high sequence polymorphism. C. sativus RADiation sensitive 51A (CsRAD51A) binding, histone H3 lysine 4 trimethylation (H3K4me3) modification, chromatin accessibility, and hypomethylation were positively associated with global CO landscapes and in local DNA double-strand break (DSB) hotspots and genes. The frequency and suppression of COs could be roughly predicted based on multiomic information. Differences in CO events between hybrids could be partially traced to distinct genetic and epigenetic features and were significantly associated with specific DSB hotspots and heterozygous SVs. Our findings identify the genomic and epigenetic features that contribute to CO formation and hybrid-specific divergence in cucumber and provide theoretical support for selecting parental combinations and manipulating recombination events at target genomic regions during plant breeding.

Gain-of-function of the 1-aminocyclopropane-1-carboxylate synthase gene ACS1G induces female flower development in cucumber gynoecy.

Unisexual flowers provide a useful system for studying plant sex determination. In cucumber (Cucumis sativus L.), three major Mendelian loci control unisexual flower development, Female (F), androecious [a; 1-aminocyclopropane-1-carboxylate {ACC} synthase 11, acs11], and Monoecious (M; ACS2), referred to here as the Female, Androecious, Monoecious (FAM) model, in combination with two genes, gynoecious (g, the WIP family C2H2 zinc finger transcription factor gene WIP1) and the ethylene biosynthetic gene ACC oxidase 2 (ACO2). The F locus, conferring gynoecy and the potential for increasing fruit yield, is defined by a 30.2-kb tandem duplication containing three genes. However, the gene that determines the Female phenotype, and its mechanism, remains unknown. Here, we created a set of mutants and revealed that ACS1G is responsible for gynoecy conferred by the F locus. The duplication resulted in ACS1G acquiring a new promoter and expression pattern; in plants carrying the F locus duplication, ACS1G is expressed early in floral bud development, where it functions with ACO2 to generate an ethylene burst. The resulting ethylene represses WIP1 and activates ACS2 to initiate gynoecy. This early ACS1G expression bypasses the need for ACS11 to produce ethylene, thereby establishing a dominant pathway for female floral development. Based on these findings, we propose a model for how these ethylene biosynthesis genes cooperate to control unisexual flower development in cucumber.

A natural promoter variation of SlBBX31 confers enhanced cold tolerance during tomato domestication.

Cold stress affects crop growth and productivity worldwide. Understanding the genetic basis of cold tolerance in germplasms is critical for crop improvement. Plants can coordinate environmental stimuli of light and temperature to regulate cold tolerance. However, it remains unknown which gene in germplasms could have such function. Here, we utilized genome-wide association study (GWAS) to investigate the cold tolerance of wild and cultivated tomato accessions and discovered that increased cold tolerance is accompanied with tomato domestication. We further identified a 27-bp InDel in the promoter of the CONSTANS-like transcription factor (TF) SlBBX31 is significantly linked with cold tolerance. Coincidentally, a key regulator of light signalling, SlHY5, can directly bind to the SlBBX31 promoter to activate SlBBX31 transcription while the 27-bp InDel can prevent S1HY5 from transactivating SlBBX31. Parallel to these findings, we observed that the loss of function of SlBBX31 results in impaired tomato cold tolerance. SlBBX31 can also modulate the cold-induced expression of several ERF TFs including CBF2 and DREBs. Therefore, our study has uncovered that SlBBX31 is possibly selected during tomato domestication for cold tolerance regulation, providing valuable insights for the development of hardy tomato varieties.

Revisiting ancient polyploidy in leptosporangiate ferns.

Ferns, and particularly homosporous ferns, have long been assumed to have experienced recurrent whole-genome duplication (WGD) events because of their substantially large genome sizes, surprisingly high chromosome numbers, and high degrees of polyploidy among many extant members. As the number of sequenced fern genomes is limited, recent studies have employed transcriptome data to find evidence for WGDs in ferns. However, they have reached conflicting results concerning the occurrence of ancient polyploidy, for instance, in the lineage of leptosporangiate ferns. Because identifying WGDs in a phylogenetic context is the foremost step in studying the contribution of ancient polyploidy to evolution, we here revisited earlier identified WGDs in leptosporangiate ferns, mainly the core leptosporangiate ferns, by building KS -age distributions and applying substitution rate corrections and by conducting statistical gene tree-species tree reconciliation analyses. Our integrative analyses not only identified four ancient WGDs in the sampled core leptosporangiate ferns but also identified false positives and false negatives for WGDs that recent studies have reported earlier. In conclusion, we underscore the significance of substitution rate corrections and uncertainties in gene tree-species tree reconciliations in calling WGD events and advance an exemplar workflow to overcome such often-overlooked issues.