A high-continuity and annotated tomato reference genome

SoIR: a comprehensive Solanaceae information resource for comparative and functional genomic study

The Solanaceae family, which includes economically important crops such as tomatoes, potatoes and peppers, has experienced a rapid expansion in genomic data due to advancements in sequencing technologies. However, existing databases are limited by incomplete species representation, a lack of comprehensive comparative genomic tools and the absence of systematic pan-genomic analyses. To address these gaps, we developed the Solanaceae Information Resource (SoIR, https://soir.bio2db.com), a comprehensive genomics database for the Solanaceae family. SoIR integrates genomic data from 81 species and transcriptomic data from 41 species, encompassing a total of 3 908 408 gene annotations derived from Gene Ontology, nonredundant protein, Pfam, Swiss-Prot and TrEMBL databases. The resource also includes 3 437 115 CRISPR guide sequences, 212 395 transcription factors and 19 086 genes associated with methylation modification. In addition to species-specific analyses, SoIR provides extensive bioinformatics tools for investigating gene family evolution, phylogenetic relationships and karyotype reconstruction across 25 fully sequenced genomes. With advanced tools such as Blast, Synteny and Sequence Alignment, the platform provides users with interactive and intuitive visualizations for conducting cross-species comparative genomics. As the first comprehensive pan-genomic resource for the entire Solanaceae family, SoIR facilitates in-depth cross-species analysis, supporting global research initiatives in plant evolution, functional genomics and crop improvement.

A stepwise guide for pangenome development in crop plants: an alfalfa (Medicago sativa) case study

BACKGROUND

The concept of pangenomics and the importance of structural variants is gaining recognition within the plant genomics community. Due to advancements in sequencing and computational technology, it has become feasible to sequence the entire genome of numerous individuals of a single species at a reasonable cost. Pangenomes have been constructed for many major diploid crops, including rice, maize, soybean, sorghum, pearl millet, peas, sunflower, grapes, and mustards. However, pangenomes for polyploid species are relatively scarce and are available in only few crops including wheat, cotton, rapeseed, and potatoes.

MAIN BODY

In this review, we explore the various methods used in crop pangenome development, discussing the challenges and implications of these techniques based on insights from published pangenome studies. We offer a systematic guide and discuss the tools available for constructing a pangenome and conducting downstream analyses. Alfalfa, a highly heterozygous, cross pollinated and autotetraploid forage crop species, is used as an example to discuss the concerns and challenges offered by polyploid crop species. We conducted a comparative analysis using linear and graph-based methods by constructing an alfalfa graph pangenome using three publicly available genome assemblies. To illustrate the intricacies captured by pangenome graphs for a complex crop genome, we used five different gene sequences and aligned them against the three graph-based pangenomes. The comparison of the three graph pangenome methods reveals notable variations in the genomic variation captured by each pipeline.

CONCLUSION

Pangenome resources are proving invaluable by offering insights into core and dispensable genes, novel gene discovery, and genome-wide patterns of variation. Developing user-friendly online portals for linear pangenome visualization has made these resources accessible to the broader scientific and breeding community. However, challenges remain with graph-based pangenomes including compatibility with other tools, extraction of sequence for regions of interest, and visualization of genetic variation captured in pangenome graphs. These issues necessitate further refinement of tools and pipelines to effectively address the complexities of polyploid, highly heterozygous, and cross-pollinated species.

Solanum pimpinellifolium exhibits complex genetic resistance to Pseudomonas syringae pv. tomato

Pseudomonas syringae pv. tomato (Pst) is the causal agent of bacterial speck disease in tomatoes. The Pto/Prf gene cluster from Solanum pimpinellifolium was introgressed into several modern tomato cultivars and provided protection against Pst race 0 strains for many decades. However, virulent Pst race 1 strains that evade Pto-mediated immunity now predominate in tomato-growing regions worldwide. Here we report the identification of resistance to a Pst race 1 strain (Pst19) in the wild tomato accession S. pimpinellifolium LA1589 (hereafter LA1589), using our rapid high-throughput seedling screen. LA1589 supports less bacterial growth than cultivars, and does not exhibit a hypersensitive response to Pst19. We tested an existing set of 87 Inbred Backcross Lines (IBLs) derived from a cross between susceptible Solanum lycopersicum E-6203 and Solanum pimpinellifolium LA1589 for resistance to Pst19. Using single-marker analysis, we identified three genomic regions associated with resistance. Bacterial growth assays on IBLs confirmed that these regions contribute to resistance in planta. We also mapped candidate genes associated with resistance in a cross between the Solanum lycopersicum var. lycopersicum cultivar Heinz BG-1706 and S. pimpinellifolium LA1589. By comparing candidates from the two mapping approaches, we were able to identify 3 QTL and 5 candidate genes in LA1589 for a role in resistance to Pst19. This work will assist in molecular marker-assisted breeding to protect tomato from bacterial speck disease.

Insights from the first chromosome-level genome assembly of the alpine gentian Gentiana straminea Maxim

Gentiana straminea Maxim. is a perennial herb and mainly distributed in the Qinghai-Tibetan Plateau. To adapt to the extreme environment, it has developed particular morphological, physiological, and genetic structures. Also, rich in iridoids, it is one of the original plants of traditional Chinese herb 'Qinjiao'. Herein, we present its first chromosome-level genome sequence assembly and compare it with the genomes of other Gentiana species to facilitate the analysis of genomic characteristics. The assembled genome size of G. straminea was 1.25 Gb, with a contig N50 of 7.5 Mb. A total of 96.08% of the genome sequences was anchored on 13 pseudochromosomes, with a scaffold N50 of 92.70 Mb. A total of 54,310 protein-coding genes were predicted, 80.25% of which were functionally annotated. Comparative genomic analyses indicated that G. straminea experienced two whole-genome duplication events after the γ whole-genome triplication with other eudicots, and it diverged from other Gentiana species at ~3.2 Mya. A total of 142 enzyme-coding genes related to iridoid biosynthesis were identified in its genome. Additionally, we identified differences in the number and expression patterns of iridoid biosynthetic pathway genes in G. straminea compared with two other Gentiana species by integrating whole-genome sequence and transcriptomic analyses.

Whole-genome landscape of histone H3K4me3 modification during sperm cell lineage development in tomato

BACKGROUND

During male gametogenesis of flowering plants, sperm cell lineage (microspores, generative cells, and sperm cells) differentiated from somatic cells and acquired different cell fates. Trimethylation of histone H3 on lysine 4 (H3K4me3) epigenetically contributes to this process, however, it remained unclear how H3K4me3 influences the gene expression in each cell type. Here, we conducted chromatin immunoprecipitation sequencing (ChIP-seq) to obtain a genome-wide landscape of H3K4me3 during sperm cell lineage development in tomato (Solanum lycopersicum).

RESULTS

We show that H3K4me3 peaks were mainly enriched in the promoter regions, and intergenic H3K4me3 peaks expanded as sperm cell lineage differentiated from somatic cells. H3K4me3 was generally positively associated with transcript abundance and served as a better indicator of gene expression in somatic and vegetative cells, compared to sperm cell lineage. H3K4me3 was mutually exclusive with DNA methylation at 3' proximal of the transcription start sites. The microspore maintained the H3K4me3 features of somatic cells, while generative cells and sperm cells shared an almost identical H3K4me3 pattern which differed from that of the vegetative cell. After microspore division, significant loss of H3K4me3 in genes related to brassinosteroid and cytokinin signaling was observed in generative cells and vegetative cells, respectively.

CONCLUSIONS

Our results suggest the asymmetric division of the microspore significantly reshapes the genome-wide distribution of H3K4me3. Selective loss of H3K4me3 in genes related to hormone signaling may contribute to functional differentiation of sperm cell lineage. This work provides new resource data for the epigenetic studies of gametogenesis in plants.

Chromosome-level genome assembly of Solanum pimpinellifolium

Solanum pimpinellifolium, the closest wild relative of the domesticated tomato, has high potential for use in breeding programs aimed at developing multi-pathogen resistance and quality improvement. We generated a chromosome-level genome assembly of S. pimpinellifolium LA1589, with a size of 833 Mb and a contig N50 of 31 Mb. We anchored 98.80% of the contigs into 12 pseudo-chromosomes, and identified 74.47% of the sequences as repetitive sequences. The genome evaluation revealed BUSCO and LAI score of 98.3% and 14.49, respectively, indicating high quality of this assembly. A total of 41,449 protein-coding genes were predicted in the genome, of which 89.17% were functionally annotated. This high-quality genome assembly serves as a valuable resource for accelerating the biological discovery and molecular breeding of this important horticultural crop.

Genome-wide identification of the TIFY gene family in tobacco and expression analysis in response to Ralstonia solanacearum infection

The TIFY gene family plays an essential role in plant development and abiotic and biotic stress responses. In this study, genome-wide identification of TIFY members in tobacco and their expression pattern analysis in response to Ralstonia solanacearum infection were performed. A total of 33 TIFY genes were identified, including the TIFY, PPD, ZIM&ZML and JAZ subfamilies. Promoter analysis results indicated that a quantity of light-response, drought-response, SA-response and JA-response cis-elements exist in promoter regions. The TIFY gene family exhibited expansion and possessed gene redundancy resulting from tobacco ploidy change. In addition, most NtTIFYs equivalently expressed in roots, stems and leaves, while NtTIFY1, NtTIFY4, NtTIFY18 and NtTIFY30 preferentially expressed in roots. The JAZ III clade showed significant expression changes after inoculation with R. solanacearum, and the expression of NtTIFY7 in resistant varieties, compared with susceptible varieties, was more stably induced. Furthermore, NtTIFY7-silenced plants, compared with the control plants, were more susceptible to bacterial wilt. These results lay a foundation for exploring the evolutionary history of TIFY gene family and revealing gene function of NtTIFYs in tobacco bacterial wilt resistance.

Long-read sequencing of extrachromosomal circular DNA and genome assembly of a Solanum lycopersicum breeding line revealed active LTR retrotransposons originating from S. Peruvianum L. introgressions

Transposable elements (TEs) are a major force in the evolution of plant genomes. Differences in the transposition activities and landscapes of TEs can vary substantially, even in closely related species. Interspecific hybridization, a widely employed technique in tomato breeding, results in the creation of novel combinations of TEs from distinct species. The implications of this process for TE transposition activity have not been studied in modern cultivars. In this study, we used nanopore sequencing of extrachromosomal circular DNA (eccDNA) and identified two highly active Ty1/Copia LTR retrotransposon families of tomato (Solanum lycopersicum), called Salsa and Ketchup. Elements of these families produce thousands of eccDNAs under controlled conditions and epigenetic stress. EccDNA sequence analysis revealed that the major parts of eccDNA produced by Ketchup and Salsa exhibited low similarity to the S. lycopersicum genomic sequence. To trace the origin of these TEs, whole-genome nanopore sequencing and de novo genome assembly were performed. We found that these TEs occurred in a tomato breeding line via interspecific introgression from S. peruvianum. Our findings collectively show that interspecific introgressions can contribute to both genetic and phenotypic diversity not only by introducing novel genetic variants, but also by importing active transposable elements from other species.

Identification of two mutant JASON-RELATED genes associated with unreduced pollen production in potato

KEY MESSAGE

Multiple QTLs control unreduced pollen production in potato. Two major-effect QTLs co-locate with mutant alleles of genes with homology to AtJAS, a known regulator of meiotic spindle orientation. In diploid potato the production of unreduced gametes with a diploid (2n) rather than a haploid (n) number of chromosomes has been widely reported. Besides their evolutionary important role in sexual polyploidisation, unreduced gametes also have a practical value for potato breeding as a bridge between diploid and tetraploid germplasm. Although early articles argued for a monogenic recessive inheritance, the genetic basis of unreduced pollen production in potato has remained elusive. Here, three diploid full-sib populations were genotyped with an amplicon sequencing approach and phenotyped for unreduced pollen production across two growing seasons. We identified two minor-effect and three major-effect QTLs regulating this trait. The two QTLs with the largest effect displayed a recessive inheritance and an additive interaction. Both QTLs co-localised with genes encoding for putative AtJAS homologs, a key regulator of meiosis II spindle orientation in Arabidopsis thaliana. The function of these candidate genes is consistent with the cytological phenotype of mis-oriented metaphase II plates observed in the parental clones. The alleles associated with elevated levels of unreduced pollen showed deleterious mutation events: an exonic transposon insert causing a premature stop, and an amino acid change within a highly conserved domain. Taken together, our findings shed light on the natural variation underlying unreduced pollen production in potato and will facilitate interploidy breeding by enabling marker-assisted selection for this trait.

Revitalizing agriculture: next-generation genotyping and -omics technologies enabling molecular prediction of resilient traits in the Solanaceae family

This review highlights -omics research in Solanaceae family, with a particular focus on resilient traits. Extensive research has enriched our understanding of Solanaceae genomics and genetics, with historical varietal development mainly focusing on disease resistance and cultivar improvement but shifting the emphasis towards unveiling resilience mechanisms in genebank-preserved germplasm is nowadays crucial. Collecting such information, might help researchers and breeders developing new experimental design, providing an overview of the state of the art of the most advanced approaches for the identification of the genetic elements laying behind resilience. Building this starting point, we aim at providing a useful tool for tackling the global agricultural resilience goals in these crops.

Bidirectional Comparisons Revealed Functional Patterns in Interaction between Salmonella enterica and Plants

Plants may harbor the human pathogen Salmonella enterica. Interactions between S. enterica and different plant species have been studied in individual reports. However, disparities arising from the distinct experimental conditions may render a meaningful comparison very difficult. This study explored interaction patterns between different S. enterica strains including serovars Typhimurium 14028s and LT2 and serovar Senftenberg, and different plants (Arabidopsis, lettuce, and tomato) in one approach. Better persistence of S. enterica serovar Typhimurium strains was observed in all tested plants, whereas the resulting symptoms varied depending on plant species. Genes encoding pathogenesis-related proteins were upregulated in plants inoculated with Salmonella. Furthermore, transcriptome of tomato indicated dynamic responses to Salmonella, with strong and specific responses already 24 h after inoculation. By comparing with publicly accessible Arabidopsis and lettuce transcriptome results generated in a similar manner, constants and variables were displayed. Plants responded to Salmonella with metabolic and physiological adjustments, albeit with variability in reprogrammed orthologues. At the same time, Salmonella adapted to plant leaf-mimicking media with changes in biosynthesis of cellular components and adjusted metabolism. This study provides insights into the Salmonella-plant interaction, allowing for a direct comparison of responses and adaptations in both organisms.

Recent Advances in Studying the Regulation of Fruit Ripening in Tomato Using Genetic Engineering Approaches

Tomato (Solanum lycopersicum L.) is one of the most commercially essential vegetable crops cultivated worldwide. In addition to the nutritional value, tomato is an excellent model for studying climacteric fruits' ripening processes. Despite this, the available natural pool of genes that allows expanding phenotypic diversity is limited, and the difficulties of crossing using classical selection methods when stacking traits increase proportionally with each additional feature. Modern methods of the genetic engineering of tomatoes have extensive potential applications, such as enhancing the expression of existing gene(s), integrating artificial and heterologous gene(s), pointing changes in target gene sequences while keeping allelic combinations characteristic of successful commercial varieties, and many others. However, it is necessary to understand the fundamental principles of the gene molecular regulation involved in tomato fruit ripening for its successful use in creating new varieties. Although the candidate genes mediate ripening have been identified, a complete picture of their relationship has yet to be formed. This review summarizes the latest (2017-2023) achievements related to studying the ripening processes of tomato fruits. This work attempts to systematize the results of various research articles and display the interaction pattern of genes regulating the process of tomato fruit ripening.

Divergent Retention of Sucrose Metabolism Genes after Whole Genome Triplication in the Tomato (Solanum lycopersicum)

Sucrose, the primary carbon transport mode and vital carbohydrate for higher plants, significantly impacts plant growth, development, yield, and quality formation. Its metabolism involves three key steps: synthesis, transport, and degradation. Two genome triplication events have occurred in Solanaceae, which have resulted in massive gene loss. In this study, a total of 48 and 65 genes from seven sucrose metabolism gene families in Vitis vinifera and Solanum lycopersicum were identified, respectively. The number of members comprising the different gene families varied widely. And there were significant variations in the pattern of gene duplication and loss in the tomato following two WGD events. Tandem duplication is a major factor in the expansion of the SWEET and Acid INV gene families. All the genes are irregularly distributed on the chromosomes, with the majority of the genes showing collinearity with the grape, particularly the CIN family. And the seven gene families were subjected to a purifying selection. The expression patterns of the different gene families exhibited notable variations. This study presents basic information about the sucrose metabolism genes in the tomato and grape, and paves the way for further investigations into the impact of SCT events on the phylogeny, gene retention duplication, and function of sucrose metabolism gene families in the tomato or Solanaceae, and the adaptive evolution of the tomato.

Leveraging Single-Cell Populations to Uncover the Genetic Basis of Complex Traits

The ease and throughput of single-cell genomics have steadily improved, and its current trajectory suggests that surveying single-cell populations will become routine. We discuss the merger of quantitative genetics with single-cell genomics and emphasize how this synergizes with advantages intrinsic to plants. Single-cell population genomics provides increased detection resolution when mapping variants that control molecular traits, including gene expression or chromatin accessibility. Additionally, single-cell population genomics reveals the cell types in which variants act and, when combined with organism-level phenotype measurements, unveils which cellular contexts impact higher-order traits. Emerging technologies, notably multiomics, can facilitate the measurement of both genetic changes and genomic traits in single cells, enabling single-cell genetic experiments. The implementation of single-cell genetics will advance the investigation of the genetic architecture of complex molecular traits and provide new experimental paradigms to study eukaryotic genetics.

Evolution of the spinach sex-linked region within a rarely recombining pericentromeric region

Sex chromosomes have evolved independently in many different plant lineages. Here, we describe reference genomes for spinach (Spinacia oleracea) X and Y haplotypes by sequencing homozygous XX females and YY males. The long arm of 185-Mb chromosome 4 carries a 13-Mb X-linked region (XLR) and 24.1-Mb Y-linked region (YLR), of which 10 Mb is Y specific. We describe evidence that this reflects insertions of autosomal sequences creating a "Y duplication region" or "YDR" whose presence probably directly reduces genetic recombination in the immediately flanking regions, although both the X and Y sex-linked regions are within a large pericentromeric region of chromosome 4 that recombines rarely in meiosis of both sexes. Sequence divergence estimates using synonymous sites indicate that YDR genes started diverging from their likely autosomal progenitors about 3 MYA, around the time when the flanking YLR stopped recombining with the XLR. These flanking regions have a higher density of repetitive sequences in the YY than the XX assembly and include slightly more pseudogenes compared with the XLR, and the YLR has lost about 11% of the ancestral genes, suggesting some degeneration. Insertion of a male-determining factor would have caused Y linkage across the entire pericentromeric region, creating physically small, highly recombining, terminal pseudoautosomal regions. These findings provide a broader understanding of the origin of sex chromosomes in spinach.

Development of a knowledge graph framework to ease and empower translational approaches in plant research: a use-case on grain legumes

While the continuing decline in genotyping and sequencing costs has largely benefited plant research, some key species for meeting the challenges of agriculture remain mostly understudied. As a result, heterogeneous datasets for different traits are available for a significant number of these species. As gene structures and functions are to some extent conserved through evolution, comparative genomics can be used to transfer available knowledge from one species to another. However, such a translational research approach is complex due to the multiplicity of data sources and the non-harmonized description of the data. Here, we provide two pipelines, referred to as structural and functional pipelines, to create a framework for a NoSQL graph-database (Neo4j) to integrate and query heterogeneous data from multiple species. We call this framework Orthology-driven knowledge base framework for translational research (Ortho_KB). The structural pipeline builds bridges across species based on orthology. The functional pipeline integrates biological information, including QTL, and RNA-sequencing datasets, and uses the backbone from the structural pipeline to connect orthologs in the database. Queries can be written using the Neo4j Cypher language and can, for instance, lead to identify genes controlling a common trait across species. To explore the possibilities offered by such a framework, we populated Ortho_KB to obtain OrthoLegKB, an instance dedicated to legumes. The proposed model was evaluated by studying the conservation of a flowering-promoting gene. Through a series of queries, we have demonstrated that our knowledge graph base provides an intuitive and powerful platform to support research and development programmes.

Next generation sequencing technologies to explore the diversity of germplasm resources: achievements and trends in tomato

Tomato is one of the major vegetable crops grown worldwide and a model species for genetic and biological research. Progress in genomic technologies made possible the development of forefront methods for high-scale sequencing, providing comprehensive insight into the genetic architecture of germplasm resources. This review revisits next-generation sequencing strategies and applications to investigate the diversity of tomato, describing the common platforms used for SNP genotyping of large collections, de novo sequencing, and whole genome resequencing. Significant findings in evolutionary history are outlined, thus discussing how genomics has provided new hints about the processes behind domestication. Finally, achievement and perspectives on pan-genome construction and graphical pan-genome development toward precise mining of the natural variation to be exploited for breeding purposes are presented.

Genomic analyses of rice bean landraces reveal adaptation and yield related loci to accelerate breeding

Rice bean (Vigna umbellata) is an underexploited domesticated legume crop consumed for dietary protein in Asia, yet little is known about the genetic diversity of this species. Here, we present a high-quality reference genome for a rice bean landrace (FF25) built using PacBio long-read data and a Hi-C chromatin interaction map, and assess the phylogenetic position and speciation time of rice bean within the Vigna genus. We sequence 440 landraces (two core collections), and GWAS based on data for growth sites at three widely divergent latitudes reveal loci associated with flowering and yield. Loci harboring orthologs of FUL (FRUITFULL), FT (FLOWERING LOCUS T), and PRR3 (PSEUDO-RESPONSE REGULATOR 3) contribute to the adaptation of rice bean from its low latitude center of origin towards higher latitudes, and the landraces which pyramid early-flowering alleles for these loci display maximally short flowering times. We also demonstrate that copy-number-variation for VumCYP78A6 can regulate seed-yield traits. Intriguingly, 32 landraces collected from a mountainous region in South-Central China harbor a recently acquired InDel in TFL1 (TERMINAL FLOWER1) affecting stem determinacy; these materials also have exceptionally high values for multiple human-desired traits and could therefore substantially advance breeding efforts to improve rice bean.

Rice bean is an underexploited legume crop that has many desirable properties against bio and abiotic stresses. Here, the authors report the genome assembly of this species, conduct population genetics studies and reveal the genetic variations associated with adaptation and yield traits.

Improved and Highly Efficient Agrobacterium rhizogenes-Mediated Genetic Transformation Protocol: Efficient Tools for Functional Analysis of Root-Specific Resistance Genes for Solanum lycopersicum cv. Micro-Tom

Gene function analysis, molecular breeding, and the introduction of new traits in crop plants all require the development of a high-performance genetic transformation system. In numerous crops, including tomatoes, Agrobacterium-mediated genetic transformation is the preferred method. As one of our ongoing research efforts, we are in the process of mapping a broad-spectrum nematode resistance gene (Me1) in pepper. We work to transform tomato plants with candidate genes to confer resistance to nematodes in Solanaceae members. The transformation technology development is designed to produce a reproducible, rapid, and highly effective Agrobacterium-mediated genetic transformation system of Micro-Tom. In our system, a transformation efficiency of over 90% was achieved. The entire procedure, starting from the germination of seeds to the establishment of transformed plants in soil, was completed in 53 days. We confirmed the presence of the NeoR/KanR and DsRed genes in the transformed roots by polymerase chain reaction. The hairy root plants were infected with nematodes, and after 3 months, the presence of DsRed and NeoR/KanR genes was detected in the transformant roots to confirm the long-term effectiveness of the method. The presented study may facilitate root-related research and exploration of root–pathogen interactions.