The tomato pan-genome uncovers new genes and a rare allele regulating fruit flavor

Comprehensive regulatory networks for tomato organ development based on the genome and RNAome of MicroTom tomato

MicroTom has a short growth cycle and high transformation efficiency, and is a prospective model plant for studying organ development, metabolism, and plant-microbe interactions. Here, with a newly assembled reference genome for this tomato cultivar and abundant RNA-seq data derived from tissues of different organs/developmental stages/treatments, we constructed multiple gene co-expression networks, which will provide valuable clues for the identification of important genes involved in diverse regulatory pathways during plant growth, e.g. arbuscular mycorrhizal symbiosis and fruit development. Additionally, non-coding RNAs, including miRNAs, lncRNAs, and circRNAs were also identified, together with their potential targets. Interacting networks between different types of non-coding RNAs (miRNA-lncRNA), and non-coding RNAs and genes (miRNA-mRNA and lncRNA-mRNA) were constructed as well. Our results and data will provide valuable information for the study of organ differentiation and development of this important fruit. Lastly, we established a database (http://eplant.njau.edu.cn/microTomBase/) with genomic and transcriptomic data, as well as details of gene co-expression and interacting networks on MicroTom, and this database should be of great value to those who want to adopt MicroTom as a model plant for research.

Solanum americanum genome-assisted discovery of immune receptors that detect potato late blight pathogen effectors

Potato (Solanum tuberosum) and tomato (Solanum lycopersicon) crops suffer severe losses to late blight caused by the oomycete pathogen Phytophthora infestans. Solanum americanum, a relative of potato and tomato, is globally distributed and most accessions are highly blight resistant. We generated high-quality reference genomes of four S. americanum accessions, resequenced 52 accessions, and defined a pan-NLRome of S. americanum immune receptor genes. We further screened for variation in recognition of 315P. infestans RXLR effectors in 52 S. americanum accessions. Using these genomic and phenotypic data, we cloned three NLR-encoding genes, Rpi-amr4, R02860 and R04373, that recognize cognate P. infestans RXLR effectors PITG_22825 (AVRamr4), PITG_02860 and PITG_04373. These genomic resources and methodologies will support efforts to engineer potatoes with durable late blight resistance and can be applied to diseases of other crops.

Natural variations in the Sl-AKR9 aldo/keto reductase gene impact fruit flavor volatile and sugar contents

The unique flavors of different fruits depend upon complex blends of soluble sugars, organic acids, and volatile organic compounds. 2-Phenylethanol and phenylacetaldehyde are major contributors to flavor in many foods, including tomato. In the tomato fruit, glucose, and fructose are the chemicals that most positively contribute to human flavor preferences. We identified a gene encoding a tomato aldo/keto reductase, Sl-AKR9, that is associated with phenylacetaldehyde and 2-phenylethanol contents in fruits. Two distinct haplotypes were identified; one encodes a chloroplast-targeted protein while the other encodes a transit peptide-less protein that accumulates in the cytoplasm. Sl-AKR9 effectively catalyzes reduction of phenylacetaldehyde to 2-phenylethanol. The enzyme can also metabolize sugar-derived reactive carbonyls, including glyceraldehyde and methylglyoxal. CRISPR-Cas9-induced loss-of-function mutations in Sl-AKR9 significantly increased phenylacetaldehyde and lowered 2-phenylethanol content in ripe fruit. Reduced fruit weight and increased soluble solids, glucose, and fructose contents were observed in the loss-of-function fruits. These results reveal a previously unidentified mechanism affecting two flavor-associated phenylalanine-derived volatile organic compounds, sugar content, and fruit weight. Modern varieties of tomato almost universally contain the haplotype associated with larger fruit, lower sugar content, and lower phenylacetaldehyde and 2-phenylethanol, likely leading to flavor deterioration in modern varieties.

Population analysis reveals the roles of DNA methylation in tomato domestication and metabolic diversity

DNA methylation is an important epigenetic marker, yet its diversity and consequences in tomato breeding at the population level are largely unknown. We performed whole-genome bisulfite sequencing (WGBS), RNA sequencing, and metabolic profiling on a population comprising wild tomatoes, landraces, and cultivars. A total of 8,375 differentially methylated regions (DMRs) were identified, with methylation levels progressively decreasing from domestication to improvement. We found that over 20% of DMRs overlapped with selective sweeps. Moreover, more than 80% of DMRs in tomato were not significantly associated with single-nucleotide polymorphisms (SNPs), and DMRs had strong linkages with adjacent SNPs. We additionally profiled 339 metabolites from 364 diverse accessions and further performed a metabolic association study based on SNPs and DMRs. We detected 971 and 711 large-effect loci via SNP and DMR markers, respectively. Combined with multi-omics, we identified 13 candidate genes and updated the polyphenol biosynthetic pathway. Our results showed that DNA methylation variants could complement SNP profiling of metabolite diversity. Our study thus provides a DNA methylome map across diverse accessions and suggests that DNA methylation variation can be the genetic basis of metabolic diversity in plants.

Genomic basis of selective breeding from the closest wild relative of large-fruited tomato

The long and intricate domestication history of the tomato (Solanum lycopersicum) includes selection sweeps that have not been fully explored, and these sweeps show significant evolutionary trajectories of domestication traits. Using three distinct selection strategies, we represented comprehensive selected sweeps from 53 Solanum pimpinellifolium (PIM) and 166 S. lycopersicum (BIG) accessions, which are defined as pseudo-domestication in this study. We identified 390 potential selection sweeps, some of which had a significant impact on fruit-related traits and were crucial to the pseudo-domestication process. During tomato pseudo-domestication, we discovered a minor-effect allele of the SlLEA gene related to fruit weight (FW), as well as the major haplotypes of fw2.2/cell number regulator (CNR), fw3.2/SlKLUH, and fw11.3/cell size regulator (CSR) in cultivars. Furthermore, 18 loci were found to be significantly associated with FW and six fruit-related agronomic traits in genome-wide association studies. By examining population differentiation, we identified the causative variation underlying the divergence of fruit flavonoids across the large-fruited tomatoes and validated BRI1-EMS-SUPPRESSOR 1.2 (SlBES1.2), a gene that may affect flavonoid content by modulating the MYB12 expression profile. Our results provide new research routes for the genetic basis of fruit traits and excellent genomic resources for tomato genomics-assisted breeding.

The Progression in Developing Genomic Resources for Crop Improvement

Sequencing technologies have rapidly evolved over the past two decades, and new technologies are being continually developed and commercialized. The emerging sequencing technologies target generating more data with fewer inputs and at lower costs. This has also translated to an increase in the number and type of corresponding applications in genomics besides enhanced computational capacities (both hardware and software). Alongside the evolving DNA sequencing landscape, bioinformatics research teams have also evolved to accommodate the increasingly demanding techniques used to combine and interpret data, leading to many researchers moving from the lab to the computer. The rich history of DNA sequencing has paved the way for new insights and the development of new analysis methods. Understanding and learning from past technologies can help with the progress of future applications. This review focuses on the evolution of sequencing technologies, their significant enabling role in generating plant genome assemblies and downstream applications, and the parallel development of bioinformatics tools and skills, filling the gap in data analysis techniques.

An Integrative Transcriptomics and Proteomics Approach to Identify Putative Genes Underlying Fruit Ripening in Tomato near Isogenic Lines with Long Shelf Life

The elucidation of the ripening pathways of climacteric fruits helps to reduce postharvest losses and improve fruit quality. Here, we report an integrative study on tomato ripening for two near-isogenic lines (NIL115 and NIL080) with Solanum pimpinellifolium LA0722 introgressions. A comprehensive analysis using phenotyping, molecular, transcript, and protein data were performed. Both NILs show improved fruit firmness and NIL115 also has longer shelf life compared to the cultivated parent. NIL115 differentially expressed a transcript from the APETALA2 ethylene response transcription factor family (AP2/ERF) with a potential role in fruit ripening. E4, another ERF, showed an upregulated expression in NIL115 as well as in the wild parent, and it was located physically close to a wild introgression. Other proteins whose expression levels changed significantly during ripening were identified, including an ethylene biosynthetic enzyme (ACO3) and a pectate lyase (PL) in NIL115, and an alpha-1,4 glucan phosphorylase (Pho1a) in NIL080. In this study, we provide insights into the effects of several genes underlying tomato ripening with potential impact on fruit shelf life. Data integration contributed to unraveling ripening-related genes, providing opportunities for assisted breeding.

Long-read, chromosome-scale assembly of Vitis rotundifolia cv. Carlos and its unique resistance to Xylella fastidiosa subsp. fastidiosa

BACKGROUND

Muscadine grape (Vitis rotundifolia) is resistant to many of the pathogens that negatively impact the production of common grape (V. vinifera), including the bacterial pathogen Xylella fastidiosa subsp. fastidiosa (Xfsf), which causes Pierce's Disease (PD). Previous studies in common grape have indicated Xfsf delays host immune response with a complex O-chain antigen produced by the wzy gene. Muscadine cultivars range from tolerant to completely resistant to Xfsf, but the mechanism is unknown.

RESULTS

We assembled and annotated a new, long-read genome assembly for 'Carlos', a cultivar of muscadine that exhibits tolerance, to build upon the existing genetic resources available for muscadine. We used these resources to construct an initial pan-genome for three cultivars of muscadine and one cultivar of common grape. This pan-genome contains a total of 34,970 synteny-constrained entries containing genes of similar structure. Comparison of resistance gene content between the 'Carlos' and common grape genomes indicates an expansion of resistance (R) genes in 'Carlos.' We further identified genes involved in Xfsf response by transcriptome sequencing 'Carlos' plants inoculated with Xfsf. We observed 234 differentially expressed genes with functions related to lipid catabolism, oxidation-reduction signaling, and abscisic acid (ABA) signaling as well as seven R genes. Leveraging public data from previous experiments of common grape inoculated with Xfsf, we determined that most differentially expressed genes in the muscadine response were not found in common grape, and three of the R genes identified as differentially expressed in muscadine do not have an ortholog in the common grape genome.

CONCLUSIONS

Our results support the utility of a pan-genome approach to identify candidate genes for traits of interest, particularly disease resistance to Xfsf, within and between muscadine and common grape.

k-mer-Based Genome-Wide Association Studies in Plants: Advances, Challenges, and Perspectives

Genome-wide association studies (GWAS) have allowed the discovery of marker-trait associations in crops over recent decades. However, their power is hampered by a number of limitations, with the key one among them being an overreliance on single-nucleotide polymorphisms (SNPs) as molecular markers. Indeed, SNPs represent only one type of genetic variation and are usually derived from alignment to a single genome assembly that may be poorly representative of the population under study. To overcome this, k-mer-based GWAS approaches have recently been developed. k-mer-based GWAS provide a universal way to assess variation due to SNPs, insertions/deletions, and structural variations without having to specifically detect and genotype these variants. In addition, k-mer-based analyses can be used in species that lack a reference genome. However, the use of k-mers for GWAS presents challenges such as data size and complexity, lack of standard tools, and potential detection of false associations. Nevertheless, efforts are being made to overcome these challenges and a general analysis workflow has started to emerge. We identify the priorities for k-mer-based GWAS in years to come, notably in the development of user-friendly programs for their analysis and approaches for linking significant k-mers to sequence variation.

The Effect of Methodological Considerations on the Construction of Gene-Based Plant Pan-genomes

Pan-genomics is an emerging approach for studying the genetic diversity within plant populations. In contrast to common resequencing studies that compare whole genome sequencing data with a single reference genome, the construction of a pan-genome (PG) involves the direct comparison of multiple genomes to one another, thereby enabling the detection of genomic sequences and genes not present in the reference, as well as the analysis of gene content diversity. Although multiple studies describing PGs of various plant species have been published in recent years, a better understanding regarding the effect of the computational procedures used for PG construction could guide researchers in making more informed methodological decisions. Here, we examine the effect of several key methodological factors on the obtained gene pool and on gene presence-absence detections by constructing and comparing multiple PGs of Arabidopsis thaliana and cultivated soybean, as well as conducting a meta-analysis on published PGs. These factors include the construction method, the sequencing depth, and the extent of input data used for gene annotation. We observe substantial differences between PGs constructed using three common procedures (de novo assembly and annotation, map-to-pan, and iterative assembly) and that results are dependent on the extent of the input data. Specifically, we report low agreement between the gene content inferred using different procedures and input data. Our results should increase the awareness of the community to the consequences of methodological decisions made during the process of PG construction and emphasize the need for further investigation of commonly applied methodologies.

Exploiting genetic and genomic resources to enhance productivity and abiotic stress adaptation of underutilized pulses

Underutilized pulses and their wild relatives are typically stress tolerant and their seeds are packed with protein, fibers, minerals, vitamins, and phytochemicals. The consumption of such nutritionally dense legumes together with cereal-based food may promote global food and nutritional security. However, such species are deficient in a few or several desirable domestication traits thereby reducing their agronomic value, requiring further genetic enhancement for developing productive, nutritionally dense, and climate resilient cultivars. This review article considers 13 underutilized pulses and focuses on their germplasm holdings, diversity, crop-wild-crop gene flow, genome sequencing, syntenic relationships, the potential for breeding and transgenic manipulation, and the genetics of agronomic and stress tolerance traits. Recent progress has shown the potential for crop improvement and food security, for example, the genetic basis of stem determinacy and fragrance in moth bean and rice bean, multiple abiotic stress tolerant traits in horse gram and tepary bean, bruchid resistance in lima bean, low neurotoxin in grass pea, and photoperiod induced flowering and anthocyanin accumulation in adzuki bean have been investigated. Advances in introgression breeding to develop elite genetic stocks of grass pea with low β-ODAP (neurotoxin compound), resistance to Mungbean yellow mosaic India virus in black gram using rice bean, and abiotic stress adaptation in common bean, using genes from tepary bean have been carried out. This highlights their potential in wider breeding programs to introduce such traits in locally adapted cultivars. The potential of de-domestication or feralization in the evolution of new variants in these crops are also highlighted.

Pan-genome of Citrullus genus highlights the extent of presence/absence variation during domestication and selection

The rich genetic diversity in Citrullus lanatus and the other six species in the Citrullus genus provides important sources in watermelon breeding. Here, we present the Citrullus genus pan-genome based on the 400 Citrullus genus resequencing data, showing that 477 Mb contigs and 6249 protein-coding genes were absent in the Citrullus lanatus reference genome. In the Citrullus genus pan-genome, there are a total of 8795 (30.5%) genes that exhibit presence/absence variations (PAVs). Presence/absence variation (PAV) analysis showed that a lot of gene PAV were selected during the domestication and improvement, such as 53 favorable genes and 40 unfavorable genes were identified during the C. mucosospermus to C. lanatus landrace domestication. We also identified 661 resistance gene analogs (RGAs) in the Citrullus genus pan-genome, which contains 90 RGAs (89 variable and 1 core gene) located on the pangenome additional contigs. By gene PAV-based GWAS, 8 gene presence/absence variations were found associated with flesh color. Finally, based on the results of gene PAV selection analysis between watermelon populations with different fruit colors, we identified four non-reference candidate genes associated with carotenoid accumulation, which had a significantly higher frequency in the white flesh. These results will provide an important source for watermelon breeding.

Insights into the Transcriptomics of Crop Wild Relatives to Unravel the Salinity Stress Adaptive Mechanisms

The narrow genomic diversity of modern cultivars is a major bottleneck for enhancing the crop's salinity stress tolerance. The close relatives of modern cultivated plants, crop wild relatives (CWRs), can be a promising and sustainable resource to broaden the diversity of crops. Advances in transcriptomic technologies have revealed the untapped genetic diversity of CWRs that represents a practical gene pool for improving the plant's adaptability to salt stress. Thus, the present study emphasizes the transcriptomics of CWRs for salinity stress tolerance. In this review, the impacts of salt stress on the plant's physiological processes and development are overviewed, and the transcription factors (TFs) regulation of salinity stress tolerance is investigated. In addition to the molecular regulation, a brief discussion on the phytomorphological adaptation of plants under saline environments is provided. The study further highlights the availability and use of transcriptomic resources of CWR and their contribution to pangenome construction. Moreover, the utilization of CWRs' genetic resources in the molecular breeding of crops for salinity stress tolerance is explored. Several studies have shown that cytoplasmic components such as calcium and kinases, and ion transporter genes such as Salt Overly Sensitive 1 (SOS1) and High-affinity Potassium Transporters (HKTs) are involved in the signaling of salt stress, and in mediating the distribution of excess Na⁺ ions within the plant cells. Recent comparative analyses of transcriptomic profiling through RNA sequencing (RNA-Seq) between the crops and their wild relatives have unraveled several TFs, stress-responsive genes, and regulatory proteins for generating salinity stress tolerance. This review specifies that the use of CWRs transcriptomics in combination with modern breeding experimental approaches such as genomic editing, de novo domestication, and speed breeding can accelerate the CWRs utilization in the breeding programs for enhancing the crop's adaptability to saline conditions. The transcriptomic approaches optimize the crop genomes with the accumulation of favorable alleles that will be indispensable for designing salt-resilient crops.

The dissection of tomato flavor: biochemistry, genetics, and omics

Flavor and quality are the major drivers of fruit consumption in the US. However, the poor flavor of modern commercial tomato varieties is a major cause of consumer dissatisfaction. Studies in flavor research have informed the role of volatile organic compounds in improving overall liking and sweetness of tomatoes. These studies have utilized and applied the tools of molecular biology, genetics, biochemistry, omics, machine learning, and gene editing to elucidate the compounds and biochemical pathways essential for good tasting fruit. Here, we discuss the progress in identifying the biosynthetic pathways and chemical modifications of important tomato volatile compounds. We also summarize the advances in developing highly flavorful tomato varieties and future steps toward developing a "perfect tomato".

Identification and validation of new MADS-box homologous genes in 3010 rice pan-genome

KEY MESSAGE

Identification and validation of ten new MADS-box homologous genes in 3010 rice pan-genome for rice breeding. The functional genome is significant for rice breeding. MADS-box genes encode transcription factors that are indispensable for rice growth and development. The reported 15,362 novel genes in the rice pan-genome (RPAN) of Asian cultivated rice accessions provided a useful gene reservoir for the identification of more MADS-box candidates to overcome the limitation for the usage of only 75 MADS-box genes identified in Nipponbare for rice breeding. Here, we report the identification and validation of ten MADS-box homologous genes in RPAN. Origin and identity analysis indicated that they are originated from different wild rice accessions and structure of motif analysis revealed high variations in their amino acid sequences. Phylogenetic results with 277 MADS-box genes in 41 species showed that all these ten MADS-box homologous genes belong to type I (SRF-like, M-type). Gene expression analysis confirmed the existence of these ten MADS-box genes in IRIS_313-10,394, all of them were expressed in flower tissues, and six of them were highly expressed during seed development. Altogether, we identified and validated experimentally, for the first time, ten novel MADS-box genes in RPAN, which provides new genetic sources for rice improvement.

A comprehensive metabolic map reveals major quality regulations in red-flesh kiwifruit (Actinidia chinensis)

Kiwifruit (Actinidia chinensis) is one of the popular fruits world-wide, and its quality is mainly determined by key metabolites (sugars, flavonoids, and vitamins). Previous works on kiwifruit are mostly done via a single omics approach or involve only limited metabolites. Consequently, the dynamic metabolomes during kiwifruit development and ripening and the underlying regulatory mechanisms are poorly understood. In this study, using high-resolution metabolomic and transcriptomic analyses, we investigated kiwifruit metabolic landscapes at 11 different developmental and ripening stages and revealed a parallel classification of 515 metabolites and their co-expressed genes into 10 distinct metabolic vs gene modules (MM vs GM). Through integrative bioinformatics coupled with functional genomic assays, we constructed a global map and uncovered essential transcriptomic and transcriptional regulatory networks for all major metabolic changes that occurred throughout the kiwifruit growth cycle. Apart from known MM vs GM for metabolites such as soluble sugars, we identified novel transcription factors that regulate the accumulation of procyanidins, vitamin C, and other important metabolites. Our findings thus shed light on the kiwifruit metabolic regulatory network and provide a valuable resource for the designed improvement of kiwifruit quality.

Limitations and advantages of using metabolite-based genome-wide association studies: focus on fruit quality traits

In the last decades, linkage mapping has help in the location of metabolite quantitative trait loci (QTL) in many species; however, this approach shows some limitations. Recently, thanks to the most recent advanced in high-throughput genotyping technologies like next-generation sequencing, metabolite genome-wide association study (mGWAS) has been proposed a powerful tool to identify the genetic variants in polygenic agrinomic traits. Fruit flavor is a complex interaction of aroma volatiles and taste being sugar and acid ratio key parameter for flavor acceptance. Here, we review recent progress of mGWAS in pinpoint gene polymorphisms related to flavor-related metabolites in fruits. Despite clear successes in discovering novel genes or regions associated with metabolite accumulation affecting sensory attributes in fruits, GWAS incurs in several limitations summarized in this review. In addition, in our own work, we performed mGWAS on 194 Citrus grandis accessions to investigate the genetic control of individual primary and lipid metabolites in ripe fruit. We have identified a total of 667 associations for 14 primary metabolites including amino acids, sugars, and organic acids, and 768 associations corresponding to 47 lipids. Furthermore, candidate genes related to important metabolites related to fruit quality such as sugars, organic acids and lipids were discovered.

Prioritizing Colombian plant genetic resources for investment in research using indicators about the geographic origin, vulnerability status, economic benefits, and food security importance

Germplasm banks are the most significant repository for plant genetic resources for food and agriculture (PGRFA) worldwide. Despite their strategic importance, national germplasm banks of tropical megadiverse developing countries such as Colombia have extremely limited funds. Therefore, making strategic decisions about research investment is essential. Here, we designed a data-driven approach to build an index that sorts Colombian PGRFA into three groups: high, medium, or low priority, based on four pillars of information from open-access databases and aligned with the sustainable goals of no poverty and zero hunger: Geographic origin, vulnerability status, economic benefits, and food security importance. We analyzed 345 PGRFA using the index, separating them into two groups, 275 already conserved in the Colombian germplasm bank (BGVCOL group) and 70 not currently conserved in the BGVCOL (NCB group). We used fuzzy logic to classify each PGRFA by each pillar and integrate it to obtain a priority index. Missing data for native crops were frequent in the BGVCOL group. Therefore we adopted an imputation strategy to fill the gaps and calculated the uncertainty. After applying the index, PGRFA with higher priority were 24 (8.72%) from the BGVCOL (i.e., 15 potatoes, three tomatoes, two tree tomatoes, pineapple, cocoa, papaya, and yacon) and one from NCB (i.e., coffee). We concluded that this methodology successfully prioritized PGRFA in Colombia and shows the big holes of knowledge for future research and alternatives to improve this index. The versatility of this methodology could be helpful in other genebanks with budget limitations for research investment.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10531-023-02599-7.

A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?

As large-scale genomic studies have progressed, it has been revealed that a single reference genome pattern cannot represent genetic diversity at the species level. While domestic animals tend to have complex routes of origin and migration, suggesting a possible omission of some population-specific sequences in the current reference genome. Conversely, the pangenome is a collection of all DNA sequences of a species that contains sequences shared by all individuals (core genome) and is also able to display sequence information unique to each individual (variable genome). The progress of pangenome research in humans, plants and domestic animals has proved that the missing genetic components and the identification of large structural variants (SVs) can be explored through pangenomic studies. Many individual specific sequences have been shown to be related to biological adaptability, phenotype and important economic traits. The maturity of technologies and methods such as third-generation sequencing, Telomere-to-telomere genomes, graphic genomes, and reference-free assembly will further promote the development of pangenome. In the future, pangenome combined with long-read data and multi-omics will help to resolve large SVs and their relationship with the main economic traits of interest in domesticated animals, providing better insights into animal domestication, evolution and breeding. In this review, we mainly discuss how pangenome analysis reveals genetic variations in domestic animals (sheep, cattle, pigs, chickens) and their impacts on phenotypes and how this can contribute to the understanding of species diversity. Additionally, we also go through potential issues and the future perspectives of pangenome research in livestock and poultry.

Super-pangenome analyses highlight genomic diversity and structural variation across wild and cultivated tomato species

Effective utilization of wild relatives is key to overcoming challenges in genetic improvement of cultivated tomato, which has a narrow genetic basis; however, current efforts to decipher high-quality genomes for tomato wild species are insufficient. Here, we report chromosome-scale tomato genomes from nine wild species and two cultivated accessions, representative of Solanum section Lycopersicon, the tomato clade. Together with two previously released genomes, we elucidate the phylogeny of Lycopersicon and construct a section-wide gene repertoire. We reveal the landscape of structural variants and provide entry to the genomic diversity among tomato wild relatives, enabling the discovery of a wild tomato gene with the potential to increase yields of modern cultivated tomatoes. Construction of a graph-based genome enables structural-variant-based genome-wide association studies, identifying numerous signals associated with tomato flavor-related traits and fruit metabolites. The tomato super-pangenome resources will expedite biological studies and breeding of this globally important crop.

Deregulation of ζ-carotene desaturase in Arabidopsis and tomato exposes a unique carotenoid-derived redundant regulation of floral meristem identity and function

A level of redundancy and interplay among the transcriptional regulators of floral development safeguards a plant's reproductive success and ensures crop production. In the present study, an additional layer of complexity in the regulation of floral meristem (FM) identity and flower development is elucidated linking carotenoid biosynthesis and metabolism to the regulation of determinate flowering. The accumulation and subsequent cleavage of a diverse array of ζ-carotenes in the chloroplast biogenesis 5 (clb5) mutant of Arabidopsis results in the reprogramming of meristematic gene regulatory networks establishing FM identity mirroring that of the FM identity master regulator, APETALA1 (AP1). The immediate transition to floral development in clb5 requires long photoperiods in a GIGANTEA-independent manner, whereas AP1 is essential for the floral organ development of clb5. The elucidation of this link between carotenoid metabolism and floral development translates to tomato exposing a regulation of FM identity redundant to and initiated by AP1 and proposed to be dependent on the E class floral initiation and organ identity regulator, SEPALLATA3 (SEP3).

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Genetics and breeding of phenolic content in tomato, eggplant and pepper fruits

Phenolic acids and flavonoids are large groups of secondary metabolites ubiquitous in the plant kingdom. They are currently in the spotlight due to the numerous health benefits associated with their consumption, as well as for their vital roles in plant biological processes and in plant-environment interaction. Tomato, eggplant and pepper are in the top ten most consumed vegetables in the world, and their fruit accumulation profiles have been extensively characterized, showing substantial differences. A broad array of genetic and genomic tools has helped to identify QTLs and candidate genes associated with the fruit biosynthesis of phenolic acids and flavonoids. The aim of this review was to synthesize the available information making it easily available for researchers and breeders. The phenylpropanoid pathway is tightly regulated by structural genes, which are conserved across species, along with a complex network of regulatory elements like transcription factors, especially of MYB family, and cellular transporters. Moreover, phenolic compounds accumulate in tissue-specific and developmental-dependent ways, as different paths of the metabolic pathway are activated/deactivated along with fruit development. We retrieved 104 annotated putative orthologues encoding for key enzymes of the phenylpropanoid pathway in tomato (37), eggplant (29) and pepper (38) and compiled 267 QTLs (217 for tomato, 16 for eggplant and 34 for pepper) linked to fruit phenolic acids, flavonoids and total phenolics content. Combining molecular tools and genetic variability, through both conventional and genetic engineering strategies, is a feasible approach to improve phenolics content in tomato, eggplant and pepper. Finally, although the phenylpropanoid biosynthetic pathway has been well-studied in the Solanaceae, more research is needed on the identification of the candidate genes behind many QTLs, as well as their interactions with other QTLs and genes.

Tomato MED25 regulates fruit ripening by interacting with EIN3-like transcription factors

Fruit ripening relies on the precise spatiotemporal control of RNA polymerase II (Pol II)-dependent gene transcription, and the evolutionarily conserved Mediator (MED) coactivator complex plays an essential role in this process. In tomato (Solanum lycopersicum), a model climacteric fruit, ripening is tightly coordinated by ethylene and several key transcription factors. However, the mechanism underlying the transmission of context-specific regulatory signals from these ripening-related transcription factors to the Pol II transcription machinery remains unknown. Here, we report the mechanistic function of MED25, a subunit of the plant Mediator transcriptional coactivator complex, in controlling the ethylene-mediated transcriptional program during fruit ripening. Multiple lines of evidence indicate that MED25 physically interacts with the master transcription factors of the ETHYLENE-INSENSITIVE 3 (EIN3)/EIN3-LIKE (EIL) family, thereby playing an essential role in pre-initiation complex formation during ethylene-induced gene transcription. We also show that MED25 forms a transcriptional module with EIL1 to regulate the expression of ripening-related regulatory as well as structural genes through promoter binding. Furthermore, the EIL1-MED25 module orchestrates both positive and negative feedback transcriptional circuits, along with its downstream regulators, to fine-tune ethylene homeostasis during fruit ripening.

From tetraploid to diploid, a pangenomic approach to identify genes lost during synthetic diploidization of Eragrostis curvula

INTRODUCTION

In Eragrostis curvula, commonly known as weeping lovegrass, a synthetic diploidization event of the facultative apomictic tetraploid Tanganyika INTA cv. originated from the sexual diploid Victoria cv. Apomixis is an asexual reproduction by seeds in which the progeny is genetically identical to the maternal plant.

METHODS

To assess the genomic changes related to ploidy and to the reproductive mode occurring during diploidization, a mapping approach was followed to obtain the first E. curvula pangenome assembly. In this way, gDNA of Tanganyika INTA was extracted and sequenced in 2x250 Illumina pair-end reads and mapped against the Victoria genome assembly. The unmapped reads were used for variant calling, while the mapped reads were assembled using Masurca software.

RESULTS

The length of the assembly was 28,982,419 bp distributed in 18,032 contigs, and the variable genes annotated in these contigs rendered 3,952 gene models. Functional annotation of the genes showed that the reproductive pathway was differentially enriched. PCR amplification in gDNA and cDNA of Tanganyika INTA and Victoria was conducted to validate the presence/absence variation in five genes related to reproduction and ploidy. The polyploid nature of the Tanganyika INTA genome was also evaluated through the variant calling analysis showing the single nucleotide polymorphism (SNP) coverage and allele frequency distribution with a segmental allotetraploid pairing behavior.

DISCUSSION

The results presented here suggest that the genes were lost in Tanganyika INTA during the diploidization process that was conducted to suppress the apomictic pathway, affecting severely the fertility of Victoria cv.

A Compendium for Novel Marker-Based Breeding Strategies in Eggplant

The worldwide production of eggplant is estimated at about 58 Mt, with China, India and Egypt being the major producing countries. Breeding efforts in the species have mainly focused on increasing productivity, abiotic and biotic tolerance/resistance, shelf-life, the content of health-promoting metabolites in the fruit rather than decreasing the content of anti-nutritional compounds in the fruit. From the literature, we collected information on mapping quantitative trait loci (QTLs) affecting eggplant's traits following a biparental or multi-parent approach as well as genome-wide association (GWA) studies. The positions of QTLs were lifted according to the eggplant reference line (v4.1) and more than 700 QTLs were identified, here organized into 180 quantitative genomic regions (QGRs). Our findings thus provide a tool to: (i) determine the best donor genotypes for specific traits; (ii) narrow down QTL regions affecting a trait by combining information from different populations; (iii) pinpoint potential candidate genes.

Citrus genomic resources unravel putative genetic determinants of Huanglongbing pathogenicity

Citrus HLB caused by Candidatus Liberibacter asiaticus is a pathogen-triggered immune disease. Here, we identified putative genetic determinants of HLB pathogenicity by integrating citrus genomic resources to characterize the pan-genome of accessions that differ in their response to HLB. Genome-wide association mapping and analysis of allele-specific expression between susceptible, tolerant, and resistant accessions further refined candidates underlying the response to HLB. We first developed a phased diploid assembly of Citrus sinensis 'Newhall' genome and produced resequencing data for 91 citrus accessions that differ in their response to HLB. These data were combined with previous resequencing data from 356 accessions for genome-wide association mapping of the HLB response. Genes determinants for HLB pathogenicity were associated with host immune response, ROS production, and antioxidants. Overall, this study has provided a significant resource of citrus genomic data and identified candidate genes to be further explored to understand the genetic determinants of HLB pathogenicity.

Heterologous mogrosides biosynthesis in cucumber and tomato by genetic manipulation

Mogrosides are widely used as high-value natural zero-calorie sweeteners that exhibit an array of biological activities and allow for vegetable flavour breeding by modern molecular biotechnology. In this study, we developed an In-fusion based gene stacking strategy for transgene stacking and a multi-gene vector harbouring 6 mogrosides biosynthesis genes and transformed it into Cucumis sativus and Lycopersicon esculentum. Here we show that transgenic cucumber can produce mogroside V and siamenoside I at 587 ng/g FW and 113 ng/g FW, respectively, and cultivated transgenic tomato with mogroside III. This study provides a strategy for vegetable flavour improvement, paving the way for heterologous biosynthesis of mogrosides.

Resequencing of durian genomes reveals large genetic variations among different cultivars

Durian (Durio zibethinus), which yields the fruit known as the "King of Fruits," is an important economic crop in Southeast Asia. Several durian cultivars have been developed in this region. In this study, we resequenced the genomes of three popular durian cultivars in Thailand, including Kradumthong (KD), Monthong (MT), and Puangmanee (PM) to investigate genetic diversities of cultivated durians. KD, MT, and PM genome assemblies were 832.7, 762.6, and 821.6 Mb, and their annotations covered 95.7, 92.4, and 92.7% of the embryophyta core proteins, respectively. We constructed the draft durian pangenome and analyzed comparative genomes with related species in Malvales. Long terminal repeat (LTR) sequences and protein families in durian genomes had slower evolution rates than that in cotton genomes. However, protein families with transcriptional regulation function and protein phosphorylation function involved in abiotic and biotic stress responses appeared to evolve faster in durians. The analyses of phylogenetic relationships, copy number variations (CNVs), and presence/absence variations (PAVs) suggested that the genome evolution of Thai durians was different from that of the Malaysian durian, Musang King (MK). Among the three newly sequenced genomes, the PAV and CNV profiles of disease resistance genes and the expressions of methylesterase inhibitor domain containing genes involved in flowering and fruit maturation in MT were different from those in KD and PM. These genome assemblies and their analyses provide valuable resources to gain a better understanding of the genetic diversity of cultivated durians, which may be useful for the future development of new durian cultivars.

Innovative Advances in Plant Genotyping

Over the past decade, advances in plant genotyping have been critical in enabling the identification of genetic diversity, in understanding evolution, and in dissecting important traits in both crops and native plants. The widespread popularity of single-nucleotide polymorphisms (SNPs) has prompted significant improvements to SNP-based genotyping, including SNP arrays, genotyping by sequencing, and whole-genome resequencing. More recent approaches, including genotyping structural variants, utilizing pangenomes to capture species-wide genetic diversity and exploiting machine learning to analyze genotypic data sets, are pushing the boundaries of what plant genotyping can offer. In this chapter, we highlight these innovations and discuss how they will accelerate and advance future genotyping efforts.

The SlDOG1 Affect Biosynthesis of Steroidal Glycoalkaloids by Regulating GAME Expression in Tomato

Steroidal alkaloids (SAs) and steroidal glycoalkaloids (SGAs) are common constituents of plant species belonging to the Solanaceae family. However, the molecular mechanism regulating the formation of SAs and SGAs remains unknown. Here, genome-wide association mapping was used to elucidate SA and SGA regulation in tomatoes: a SlGAME5-like glycosyltransferase (Solyc10g085240) and the transcription factor SlDOG1 (Solyc10g085210) were significantly associated with steroidal alkaloid composition. In this study, it was found that rSlGAME5-like can catalyze a variety of substrates for glycosidation and can catalyze SA and flavonol pathways to form O-glucoside and O-galactoside in vitro. The overexpression of SlGAME5-like promoted α-tomatine, hydroxytomatine, and flavonol glycoside accumulation in tomatoes. Furthermore, assessments of natural variation combined with functional analyses identified SlDOG1 as a major determinant of tomato SGA content, which also promoted SA and SGA accumulation via the regulation of GAME gene expression. This study provides new insights into the regulatory mechanisms underlying SGA production in tomatoes.

Biosynthesis of Chlorophyll and Other Isoprenoids in the Plastid of Red Grape Berry Skins

Despite current knowledge showing that fruits like tomato and grape berries accumulate different components of the light reactions and Calvin cycle, the role of green tissues in fruits is not yet fully understood. In mature tomato fruits, chlorophylls are degraded and replaced by carotenoids through the conversion of chloroplasts in chromoplasts, while in red grape berries, chloroplasts persist at maturity and chlorophylls are masked by anthocyanins. To study isoprenoid and lipid metabolism in grape skin chloroplasts, metabolites of enriched organelle fractions were analyzed by high-performance liquid chromatography-high-resolution mass spectrometry (HPLC-HRMS) and the expression of key genes was evaluated by real-time polymerase chain reaction (PCR) in berry skins and leaves. Overall, the results indicated that chloroplasts of the grape berry skins, as with leaf chloroplasts, share conserved mechanisms of synthesis (and degradation) of important components of the photosynthetic machinery. Some of these components, such as chlorophylls and their precursors, and catabolites, carotenoids, quinones, and lipids have important roles in grape and wine sensory characteristics.

Targeting ripening regulators to develop fruit with high quality and extended shelf life

Fruit quality directly impacts fruit marketability and consumer acceptance. Breeders have focused on fruit quality traits to extend shelf life, primarily through fruit texture, but, in some cases, have neglected other qualities such as flavor and nutrition. In recent years, integrative biotechnology and consumer-minded approaches have surfaced, aiding in the development of flavorful, long-lasting fruit. Here, we discussed how specific transcription factors and hormones involved in fruit ripening can be targeted to generate high-quality fruit through traditional breeding and bioengineering. We highlight regulators that can be used to generate novel-colored fruit or biofortify fresh produce with health-promoting nutrients, such as vitamin C. Overall, we argue that addressing grower and industry needs must be balanced with consumer-based traits.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Era of gapless plant genomes: innovations in sequencing and mapping technologies revolutionize genomics and breeding

Whole-genome sequencing and assembly have revolutionized plant genetics and molecular biology over the last two decades. However, significant shortcomings in first- and second-generation technology resulted in imperfect reference genomes: numerous and large gaps of low quality or undeterminable sequence in areas of highly repetitive DNA along with limited chromosomal phasing restricted the ability of researchers to characterize regulatory noncoding elements and genic regions that underwent recent duplication events. Recently, advances in long-read sequencing have resulted in the first gapless, telomere-to-telomere (T2T) assemblies of plant genomes. This leap forward has the potential to increase the speed and confidence of genomics and molecular experimentation while reducing costs for the research community.

Going broad and deep: sequencing-driven insights into plant physiology, evolution, and crop domestication

Deep sequencing is a term that has become embedded in the plant genomic literature in recent years and with good reason. A torrent of (largely) high-quality genomic and transcriptomic data has been collected and most of this has been publicly released. Indeed, almost 1000 plant genomes have been reported (www.plabipd.de) and the 2000 Plant Transcriptomes Project has long been completed. The EarthBioGenome project will dwarf even these milestones. That said, massive progress in understanding plant physiology, evolution, and crop domestication has been made by sequencing broadly (across a species) as well as deeply (within a single individual). We will outline the current state of the art in genome and transcriptome sequencing before we briefly review the most visible of these broad approaches, namely genome-wide association and transcriptome-wide association studies, as well as the compilation of pangenomes. This will include both (i) the most commonly used methods reliant on single nucleotide polymorphisms and short InDels and (ii) more recent examples which consider structural variants. We will subsequently present case studies exemplifying how their application has brought insight into either plant physiology or evolution and crop domestication. Finally, we will provide conclusions and an outlook as to the perspective for the extension of such approaches to different species, tissues, and biological processes.

A pangenome analysis pipeline provides insights into functional gene identification in rice

BACKGROUND

A pangenome aims to capture the complete genetic diversity within a species and reduce bias in genetic analysis inherent in using a single reference genome. However, the current linear format of most plant pangenomes limits the presentation of position information for novel sequences. Graph pangenomes have been developed to overcome this limitation. However, bioinformatics analysis tools for graph format genomes are lacking.

RESULTS

To overcome this problem, we develop a novel strategy for pangenome construction and a downstream pangenome analysis pipeline (PSVCP) that captures genetic variants' position information while maintaining a linearized layout. Using PSVCP, we construct a high-quality rice pangenome using 12 representative rice genomes and analyze an international rice panel with 413 diverse accessions using the pangenome as the reference. We show that PSVCP successfully identifies causal structural variations for rice grain weight and plant height. Our results provide insights into rice population structure and genomic diversity. We characterize a new locus (qPH8-1) associated with plant height on chromosome 8 undetected by the SNP-based genome-wide association study (GWAS).

CONCLUSIONS

Our results demonstrate that the pangenome constructed by our pipeline combined with a presence and absence variation-based GWAS can provide additional power for genomic and genetic analysis. The pangenome constructed in this study and the associated genome sequence and genetic variants data provide valuable genomic resources for rice genomics research and improvement in future.

Open avenues for carotenoid biofortification of plant tissues

Plant carotenoids are plastidial isoprenoids that function as photoprotectants, pigments, and precursors of apocarotenoids such as the hormones abscisic acid and strigolactones. Humans do not produce carotenoids but need to obtain them from their diet as precursors of retinoids, including vitamin A. Carotenoids also provide numerous other health benefits. Multiple attempts to improve the carotenoid profile of different crops have been carried out by manipulating carotenoid biosynthesis, degradation, and/or storage. Here, we will focus on open questions and emerging subjects related to the use of biotechnology for carotenoid biofortification. After impressive achievements, new efforts should be directed to extend the use of genome-editing technologies to overcome regulatory constraints and improve consumer acceptance of the carotenoid-enriched products. Another challenge is to prevent off-target effects like those resulting from altered hormone levels and metabolic homeostasis. Research on biofortification of green tissues should also look for new ways to deal with the negative impact that altered carotenoid contents may have on photosynthesis. Once a carotenoid-enriched product has been obtained, additional effort should be devoted to confirming that carotenoid intake from the engineered food is also improved. This work involves ensuring post-harvest stability and assessing bioaccessibility of the biofortified product to confirm that release of carotenoids from the food matrix has not been negatively affected. Successfully addressing these challenges will ensure new milestones in carotenoid biotechnology and biofortification.

CRISPR/Cas9-Mediated Enrichment Coupled to Nanopore Sequencing Provides a Valuable Tool for the Precise Reconstruction of Large Genomic Target Regions

Complete and accurate identification of genetic variants associated with specific phenotypes can be challenging when there is a high level of genomic divergence between individuals in a study and the corresponding reference genome. We have applied the Cas9-mediated enrichment coupled to nanopore sequencing to perform a targeted de novo assembly and accurately reconstruct a genomic region of interest. This approach was used to reconstruct a 250-kbp target region on chromosome 5 of the common bean genome (Phaseolus vulgaris) associated with the shattering phenotype. Comparing a non-shattering cultivar (Midas) with the reference genome revealed many single-nucleotide variants and structural variants in this region. We cut five 50-kbp tiled sub-regions of Midas genomic DNA using Cas9, followed by sequencing on a MinION device and de novo assembly, generating a single contig spanning the whole 250-kbp region. This assembly increased the number of Illumina reads mapping to genes in the region, improving their genotypability for downstream analysis. The Cas9 tiling approach for target enrichment and sequencing is a valuable alternative to whole-genome sequencing for the assembly of ultra-long regions of interest, improving the accuracy of downstream genotype-phenotype association analysis.

Establishing Physalis as a Solanaceae model system enables genetic reevaluation of the inflated calyx syndrome

The highly diverse Solanaceae family contains several widely studied models and crop species. Fully exploring, appreciating, and exploiting this diversity requires additional model systems. Particularly promising are orphan fruit crops in the genus Physalis, which occupy a key evolutionary position in the Solanaceae and capture understudied variation in traits such as inflorescence complexity, fruit ripening and metabolites, disease and insect resistance, self-compatibility, and most notable, the striking inflated calyx syndrome (ICS), an evolutionary novelty found across angiosperms where sepals grow exceptionally large to encapsulate fruits in a protective husk. We recently developed transformation and genome editing in Physalis grisea (groundcherry). However, to systematically explore and unlock the potential of this and related Physalis as genetic systems, high-quality genome assemblies are needed. Here, we present chromosome-scale references for P. grisea and its close relative Physalis pruinosa and use these resources to study natural and engineered variations in floral traits. We first rapidly identified a natural structural variant in a bHLH gene that causes petal color variation. Further, and against expectations, we found that CRISPR-Cas9-targeted mutagenesis of 11 MADS-box genes, including purported essential regulators of ICS, had no effect on inflation. In a forward genetics screen, we identified huskless, which lacks ICS due to mutation of an AP2-like gene that causes sepals and petals to merge into a single whorl of mixed identity. These resources and findings elevate Physalis to a new Solanaceae model system and establish a paradigm in the search for factors driving ICS.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

Involvement of IDA-HAE Module in Natural Development of Tomato Flower Abscission

The unwanted detachment of organs such as flowers, leaves, and fruits from the main body of a plant (abscission) has significant effects on agricultural practice. Both timely and precise regulation of organ abscission from a plant is crucial as it influences the agricultural yield. The tomato (Solanum lycopersicum) has become a model system for research on organ abscission. Here, we characterized four tomato natural abscission variants named jointless (j), functionally impaired jointless (fij), functionally impaired jointless like (fij like), and normal joint (NJ), based on their cellular features within the flower abscission zones (AZ). Using eight INFLORESCENCE DEFICIENT IN ABSCISSION (SlIDA) genes and eight HAESA genes (SlHAE) identified in the genome sequence of tomato, we analyzed the pattern of gene expression during flower abscission. The AZ-specific expression for three tomato abscission polygalacturonases (SlTAPGs) in the development of flower AZ, and the progression of abscission validated our natural abscission system. Compared to that of j, fij, and fij like variants, the AZ-specific expression for SlIDA, SlIDL2, SlIDL3, SlIDL4, and SlIDL5 in the NJ largely corelated and increased with the process of abscission. Of eight SlHAE genes examined, the expression for SlHSL6 and SlHSL7 were found to be AZ-specific and increased as abscission progressed in the NJ variant. Unlike the result of gene expression obtained from natural abscission system, an in silico analysis of transcriptional binding sites uncovered that SlIDA genes (SlIDA, SlIDL6, and SlIDL7) are predominantly under the control of environmental stress, while most of the SlHSL genes are affiliated with the broader context in developmental processes and stress responses. Our result presents the potential bimodal transcriptional regulation of the tomato IDA-HAE module associated with flower abscission in tomatoes.

Graph-based pan-genomes: increased opportunities in plant genomics

Due to the development of sequencing technology and the great reduction in sequencing costs, an increasing number of plant genomes have been assembled, and numerous genomes have revealed large amounts of variations. However, a single reference genome does not allow the exploration of species diversity, and therefore the concept of pan-genome was developed. A pan-genome is a collection of all sequences available for a species, including a large number of consensus sequences, large structural variations, and small variations including single nucleotide polymorphisms and insertions/deletions. A simple linear pan-genome does not allow these structural variations to be intuitively characterized, so graph-based pan-genomes have been developed. These pan-genomes store sequence and structural variation information in the form of nodes and paths to store and display species variation information in a more intuitive manner. The key role of graph-based pan-genomes is to expand the coordinate system of the linear reference genome to accommodate more regions of genetic diversity. Here, we review the origin and development of graph-based pan-genomes, explore their application in plant research, and further highlight the application of graph-based pan-genomes for future plant breeding.

Rice Pangenome Genotyping Array: an efficient genotyping solution for pangenome-based accelerated genetic improvement in rice

The advent of the pangenome era has unraveled previously unknown genetic variation existing within diverse crop plants, including rice. This untapped genetic variation is believed to account for a major portion of phenotypic variation existing in crop plants. However, the use of conventional single reference-guided genotyping often fails to capture a large portion of this genetic variation leading to a reference bias. This makes it difficult to identify and utilize novel population/cultivar-specific genes for crop improvement. Thus, we developed a Rice Pangenome Genotyping Array (RPGA) harboring probes assaying 80K single-nucleotide polymorphisms (SNPs) and presence-absence variants spanning the entire 3K rice pangenome. This array provides a simple, user-friendly and cost-effective (60-80 USD per sample) solution for rapid pangenome-based genotyping in rice. The genome-wide association study (GWAS) conducted using RPGA-SNP genotyping data of a rice diversity panel detected a total of 42 loci, including previously known as well as novel genomic loci regulating grain size/weight traits in rice. Eight of these identified trait-associated loci (dispensable loci) could not be detected with conventional single reference genome-based GWAS. A WD repeat-containing PROTEIN 12 gene underlying one of such dispensable locus on chromosome 7 (qLWR7) along with other non-dispensable loci were subsequently detected using high-resolution quantitative trait loci mapping confirming authenticity of RPGA-led GWAS. This demonstrates the potential of RPGA-based genotyping to overcome reference bias. The application of RPGA-based genotyping for population structure analysis, hybridity testing, ultra-high-density genetic map construction and chromosome-level genome assembly, and marker-assisted selection was also demonstrated. A web application (http://www.rpgaweb.com) was further developed to provide an easy to use platform for the imputation of RPGA-based genotyping data using 3K rice reference panel and subsequent GWAS.

Natural genetic variation in the HAIRS ABSENT (H) gene increases type-VI glandular trichomes in both wild and domesticated tomatoes

Glandular trichomes produce and exude secondary metabolites, conferring insect resistance in many crop species. Whereas some of its wild relatives are insect-resistant, tomato (Solanum lycopersicum) is not. Identifying the genetic changes that altered trichome development and biochemistry during tomato domestication would contribute to breeding for insect resistance. A mutation in the HAIRS ABSENT (H) gene, which encodes a C2H2 zinc finger protein (ZFP8), leads to reduced trichome density. Several geographic accessions of S. pimpinellifolium, the wild ancestor of domesticated tomato, have glabrous organs that resemble the phenotype caused by h. Here, we investigated allelic diversity for H in tomato and S. pimpinellifolium accessions and their associated trichome phenotypes. We also evaluated how the developmental stage can affect trichome development in glabrous and non-glabrous plants. We found that glabrous accessions of S. pimpinellifolium have different ZFP8 nucleotide sequence changes, associated with altered trichome development and density. We also found that while the glabrous appearance of h mutants is caused by a lower density of long trichomes, the density of type-VI glandular trichomes is increased, particularly in the adult stages of plant development. These insights on the genetic control of trichome development may contribute to breeding for insect resistance in tomatoes and other crops.

Rewiring of the seed metabolome during Tartary buckwheat domestication

Crop domestication usually leads to the narrowing genetic diversity. However, human selection mainly focuses on visible traits, such as yield and plant morphology, with most metabolic changes being invisible to the naked eye. Buckwheat accumulates abundant bioactive substances, making it a dual-purpose crop with excellent nutritional and medical value. Therefore, examining the wiring of these invisible metabolites during domestication is of major importance. The comprehensive profiling of 200 Tartary buckwheat accessions exhibits 540 metabolites modified as a consequence of human selection. Metabolic genome-wide association study illustrates 384 mGWAS signals for 336 metabolites are under selection. Further analysis showed that an R2R3-MYB transcription factor FtMYB43 positively regulates the synthesis of procyanidin. Glycoside hydrolase gene FtSAGH1 is characterized as responsible for the release of active salicylic acid, the precursor of aspirin and indispensably in plant defence. UDP-glucosyltransferase gene FtUGT74L2 is characterized as involved in the glycosylation of emodin, a major medicinal component specific in Polygonaceae. The lower expression of FtSAGH1 and FtUGT74L2 were associated with the reduction of salicylic acid and soluble EmG owing to domestication. This first large-scale metabolome profiling in Tartary buckwheat will facilitate genetic improvement of medicinal properties and disease resistance in Tartary buckwheat.

Two independent loss-of-function mutations in anthocyanidin synthase homeologous genes are responsible for the all-green phenotype of sweet basil

Sweet basil, Ocimum basilicum L., is an important culinary herb grown worldwide. Although basil is green, many landraces, breeding lines, and exotic cultivars have purple stems and flowers. This anthocyanin pigmentation is unacceptable in traditional Italian basil used for Pesto sauce production. In the current study, we aimed to resolve the genetics that underlines the different colors. We used the recently published sweet basil genome to map quantitative trait loci (QTL) for flower and stem color in a bi-parental F₂ population. It was found that the pigmentation is governed by a single QTL, harboring an anthocyanidin synthase (ANS) gene (EC 1.14.20.4). Further analysis revealed that the basil genome harbors two homeologous ANS genes, each carrying a loss-of-function mutation. ObANS1 carries a single base pair insertion resulting in a frameshift, and ObANS2 carries a missense mutation within the active site. In the purple-flower parent, ANS1 is functional, and ANS2 carries a nonsense mutation. The functionality of the ObANS1 active allele was validated by complementation assay in an Arabidopsis ANS mutant. Moreover, we have restored the functionality of the missense-mutated ObANS2 using site-directed activation. We found that the non-functional alleles were expressed to similar levels as the functional allele, suggesting polyploids invest futile effort in expressing non-functional genes, offsetting their advantageous redundancy. This work demonstrated the usefulness of the genomics and genetics of basil to understand the basic mechanism of metabolic traits and raise fundamental questions in polyploid plant biology.