Structural variations (SVs) in repetitive sequences could only be detected within a broad region due to imprecise breakpoints, leading to classification errors and inaccurate trait analysis. Through manual inspection at 4532 variant regions identified by integrating 14 detection pipelines between two tomato genomes, we generated an SV benchmark at base-pair resolution. Evaluation of all pipelines yielded F1-scores below 53.77% with this benchmark, underscoring the urgent need for advanced detection algorithms in plant genomics. Analyzing the alignment features of the repetitive sequences in each region, we summarized four patterns of SV breakpoints and revealed that deviations in breakpoint identification were primarily due to copy misalignment. According to the similarities among copies, we identified 1635 bona fide SVs with precise breakpoints, including substitutions (223), which should be taken as a fundamental SV type, alongside insertions (780), deletions (619), and inversions (13), all showing preferences for SV occurrence within AT-repeat regions of regulatory loci. This precise resolution of complex SVs will foster genome analysis and crop improvement.

NAC [NO APICAL MERISTEM (NAM), ARABIDOPSIS TRANSCRIPTION ACTIVATOR FACTOR 1/2 (ATAF1/2), and CUP-SHAPED COTYLEDON (CUC2)] transcription factors are key regulators of plant growth, development, and stress responses but were also crucial players during land plant adaptation and crop domestication. Using representative members of green algae, bryophytes, lycophytes, gymnosperms, and angiosperms, we expanded the evolutionary history of NAC transcription factors to unveil the relationships among members of this gene family. We found a massive increase in the number of NAC transcription factors from green algae to lycophytes and an even larger increase in flowering plants. Many of the NAC clades arose later during evolution since we found eudicot- and monocot-specific clades. Cis-elements analysis in NAC promoters showed the presence of abiotic and biotic stress as well as hormonal response elements, which indicate the ancestral function of NAC transcription factor genes in response to environmental stimuli and in plant development. At the transcriptional level, the expression of NAC transcription factors was low or absent in male reproduction, particularly mature pollen, across the plant kingdom. We also identified NAC genes with conserved expression patterns in response to heat stress in Marchantia polymorpha and Oryza sativa. Our study provides further evidence that transcriptional mechanisms associated with stress responses and development emerged early during plant land adaptation and are still conserved in flowering plants and domesticated crops.

Arbuscular mycorrhizal (AM) symbiosis can prime plant defenses, leading to mycorrhiza-induced resistance (MIR) against different attackers, including insect herbivores. Still, our knowledge of the complex molecular regulation leading to MIR is very limited. Here, we showed that the AM fungus Funneliformis mosseae protects tomato plants against two different chewing herbivores, Spodoptera exigua and Manduca sexta. We explored the underlying molecular mechanism through genome-wide transcriptional profiling, bioinformatics network analyses, and functional bioassays. Herbivore-triggered jasmonate (JA)-regulated defenses were primed in leaves of mycorrhizal plants. Likewise, ethylene (ET) biosynthesis and signaling were also higher in leaves of mycorrhizal plants both before and after herbivory. We hypothesized that fine-tuned ET signaling is required for the primed defense response leading to MIR. ET is a complex regulator of plant responses to stress and is generally considered a negative regulator of plant defenses against herbivory. However, ET-deficient or insensitive lines did not show AM-primed JA biosynthesis or defense response, and were unable to develop MIR against any of the herbivores. Thus, we demonstrate that hormone crosstalk is central to the priming of plant immunity by beneficial microbes, with ET fine-tuning being essential for the primed JA biosynthesis and boosted defenses leading to MIR in tomato.

Comparative genomic studies of closely related taxa are important for our understanding of the causes of divergence on a changing Earth. This being said, the genomic resources available for marine intertidal molluscs are limited and currently, there are few publicly available high-quality annotated genomes for intertidal species and for molluscs in general. Here we report transcriptome assemblies for six species of Patellogastropoda and genome assemblies and annotations for three of these species (Scurria scurra, Scurria viridula and Scurria zebrina). Comparative analysis using these genomic resources suggest that and recently diverging lineages (10-20 Mya) have experienced similar amounts of contractions and expansions but across different gene families. Furthermore, differences among recently diverged species are reflected in variation in the amount of coding and noncoding material in genomes, such as amount of repetitive elements and lengths of transcripts and introns and exons. Additionally, functional ontologies of species-specific and duplicated genes together with demographic inference support the finding that recent divergence among members of the genus Scurria aligns with their unique ecological characteristics. Overall, the resources presented here will be valuable for future studies of adaptation in molluscs and in intertidal habitats as a whole.

Recent environmental changes in the Mediterranean region, attributable to anthropogenic climate change, present a substantial challenge to the adaptive evaluation of crops and the development of novel improvement strategies. In this study, we established a hydroponic tomato cultivation protocol under in vitro conditions to analyze the transcriptomic profile of seedlings exposed to salinity stress. The study also examined the impact of Atlanticell Micomix, a biological stimulant derived from a mixture of mycorrhizal microorganisms and rhizobacteria, on plant growth and development under standard conditions and in response to moderate salinity. Our transcriptomic analysis indicated a differential effect of biostimulant inoculation compared to the effect induced by salinity stress, involving genes such as GOX3 or DIR1, which are associated with the plant's defense response to adverse conditions. In addition, the presence of a cross-regulatory module between jasmonic acid and auxin, involving potential orthologs of IAA29 and JAZ, was proposed. The application of the biostimulant demonstrated a potential priming effect on the tomato seedlings, which might be useful in reversing the transcriptomic effects caused by salt stress. A comprehensive analysis of the pathways differentially affected by the treatments facilitates further investigation into the mechanisms underlying these effects.

Potato (Solanum tuberosum) is the third-most important food crop in the world. Although the potato genome has been fully sequenced, functional genomics research of potato lags behind that of other major food crops, largely due to the lack of a model experimental potato line. Here, we present a diploid potato line, 'Jan,' which possesses all essential characteristics for facile functional genomics studies. Jan exhibits a high level of homozygosity after seven generations of self-pollination. Jan is vigorous, highly fertile and produces tubers with outstanding traits. Additionally, it demonstrates high regeneration rates and excellent transformation efficiencies. We generated a chromosome-scale genome assembly for Jan, annotated its genes and identified syntelogs relative to the potato reference genome assembly DMv6.1 to facilitate functional genomics. To miniaturize plant architecture, we developed two 'mini-Jan' lines with compact and dwarf plant stature through CRISPR/Cas9-mediated mutagenesis targeting the Dwarf and Erecta genes involved in growth. One mini-Jan mutant, mini-JanE, is fully fertile and will permit higher-throughput studies in limited growth chamber and greenhouse space. Thus, Jan and mini-Jan offer a robust model system that can be leveraged for gene editing and functional genomics research in potato.

Pan-genomics and genome-editing technologies are revolutionizing breeding of global crops1,2. A transformative opportunity lies in exchanging genotype-to-phenotype knowledge between major crops (that is, those cultivated globally) and indigenous crops (that is, those locally cultivated within a circumscribed area)3-5 to enhance our food system. However, species-specific genetic variants and their interactions with desirable natural or engineered mutations pose barriers to achieving predictable phenotypic effects, even between related crops6,7. Here, by establishing a pan-genome of the crop-rich genus Solanum8 and integrating functional genomics and pan-genetics, we show that gene duplication and subsequent paralogue diversification are major obstacles to genotype-to-phenotype predictability. Despite broad conservation of gene macrosynteny among chromosome-scale references for 22 species, including 13 indigenous crops, thousands of gene duplications, particularly within key domestication gene families, exhibited dynamic trajectories in sequence, expression and function. By augmenting our pan-genome with African eggplant cultivars9 and applying quantitative genetics and genome editing, we dissected an intricate history of paralogue evolution affecting fruit size. The loss of a redundant paralogue of the classical fruit size regulator CLAVATA3 (CLV3)10,11 was compensated by a lineage-specific tandem duplication. Subsequent pseudogenization of the derived copy, followed by a large cultivar-specific deletion, created a single fused CLV3 allele that modulates fruit organ number alongside an enzymatic gene controlling the same trait. Our findings demonstrate that paralogue diversifications over short timescales are underexplored contingencies in trait evolvability. Exposing and navigating these contingencies is crucial for translating genotype-to-phenotype relationships across species.

Late blight caused by Phytophthora infestans is a major threat to global potato and tomato production. Sustainable management of late blight requires the development of resistant crop cultivars. This process can be facilitated by high-throughput identification of functional genes involved in late blight pathogenesis. In this study, we generated a high-quality transcriptomic time-course dataset focusing on the initial 24 h of contact between P. infestans and 3 solanaceous plant species, tobacco (Nicotiana benthamiana), tomato (Solanum lycopersicum), and potato (Solanum tuberosum). Our results demonstrate species-specific transcriptional regulation in early stages of the infection. Transient silencing of putative RIBOSE-5-PHOSPHATE ISOMERASE and HMG-CoA REDUCTASE genes in N. benthamiana lowered plant resistance against P. infestans. Furthermore, heterologous expression of a putative tomato Golgi-localized nucleosugar transporter-encoding gene exacerbated P. infestans infection of N. benthamiana. In comparison, bioassays using transgenic tomato lines showed that the quantitative disease resistance genes were required but insufficient for late blight resistance; genetic knock-out of the susceptibility gene enhanced resistance. The same RNA-seq dataset was exploited to examine the transcriptional landscape of P. infestans and revealed host-specific gene expression patterns in the pathogen. This temporal transcriptomic diversity, in combination with genomic distribution features, identified the P. infestans IPI-B family GLYCINE-RICH PROTEINs as putative virulence factors that promoted disease severity or induced plant tissue necrosis when transiently expressed in N. benthamiana. These functional genes underline the effectiveness of functional gene-mining through a time-course dual RNA-seq approach and provide insight into the molecular interactions between solanaceous plants and P. infestans.

Tomato is a major global crop. However, its production is limited by Botrytis cinerea. Due to the toxicity of postharvest pesticide application, alternative control methods such as priming are being investigated. Plants were treated with β-aminobutyric acid (BABA) at two developmental stages and resistance against B. cinerea was tested in fruit tissue and in progenies. DNA methylation and RNA sequencing were conducted to characterise the (epi)genetic changes associated with long-lasting resistance. Grafting experiments were done to assess the systemic nature of this signal, which was further characterised by small RNA (sRNA) sequencing of scions. Only BABA-treated seedlings displayed induced resistance (IR). DNA methylation analysis revealed seedling-specific changes, which occurred in the context of lower basal methylation. BABA-IR was found to be transmissible from primed rootstock to grafted unprimed scions. In these scions, we identified a subset of mobile 24 nt sRNAs associated with genes showing primed expression during infection in fruit. Our results demonstrate the functional association of a systemic signal with long-lasting IR and priming. Through integrated omics approaches, we have identified markers of long-lasting priming in tomato fruit which could also serve as targets for durable resistance in other crops.

Understanding the evolution of pathogen resistance genes (nucleotide-binding site-leucine-rich repeats, NLRs) within a species requires a comprehensive examination of factors that affect gene loss and gain. We present a new reference genome of Solanum chilense, which leads to an increased number and more accurate annotation of NLRs. Using a target capture approach, we quantify the presence-absence variation (PAV) of NLR loci across 20 populations from different habitats. We build a rigorous pipeline to validate the identification of PAV of NLRs and then show that PAV is larger within populations than between populations, suggesting that maintenance of NLR diversity is linked to population dynamics. The amount of PAV appears not to be correlated with the NLR presence in gene clusters in the genome, but rather with the past demographic history of the species, with loss of NLRs in diverging (smaller) populations at the distribution edges. Finally, using a redundancy analysis, we find limited evidence of PAV being linked to environmental gradients. Our results suggest that random processes (genetic drift and demography) and weak positive selection for local adaptation shape the evolution of NLRs at the single nucleotide polymorphism and PAV levels in an outcrossing plant with high nucleotide diversity.

Controlled environment agriculture (CEA) represents one of the fastest-growing sectors of horticulture. Production in controlled environments ranges from highly controlled indoor environments with 100% artificial lighting (vertical farms or plant factories) to high-tech greenhouses with or without supplemental lighting, to simpler greenhouses and high tunnels. Although food production occurs in the soil inside high tunnels, most CEA operations use various hydroponic systems to meet crop irrigation and fertility needs. The expansion of CEA offers promise as a tool for increasing food production in and near urban systems as these systems do not rely on arable agricultural land. In addition, CEA offers resilience to climate instability by growing inside protective structures. Products harvested from CEA systems tend to be of high quality, both internal and external, and are sought after by consumers. Currently, CEA producers rely on cultivars bred for production in open-field agriculture. Because of high energy and other production costs in CEA, only a limited number of food crops have proven themselves to be profitable to produce. One factor contributing to this situation may be a lack of optimized cultivars. Indoor growing operations offer opportunities for breeding cultivars that are ideal for these systems. To facilitate breeding these specialized cultivars, a wide range of tools are available for plant breeders to help speed this process and increase its efficiency. This review aims to cover breeding opportunities and needs for a wide range of horticultural crops either already being produced in CEA systems or with potential for CEA production. It also reviews many of the tools available to breeders including genomics-informed breeding, marker-assisted selection, precision breeding, high-throughput phenotyping, and potential sources of germplasm suitable for CEA breeding. The availability of published genomes and trait-linked molecular markers should enable rapid progress in the breeding of CEA-specific food crops that will help drive the growth of this industry.

Steroidal glycoalkaloids (SGAs) are specialized metabolites produced by hundreds of Solanum species including food crops, such as tomato, potato and eggplant. Unlike true alkaloids, nitrogen is introduced at a late stage of SGA biosynthesis through an unknown transamination reaction. Here, we reveal the mechanism by which GLYCOALKALOID METABOLISM12 (GAME12) directs the biosynthesis of nitrogen-containing steroidal alkaloid aglycone in Solanum. We report that GAME12, a neofunctionalized γ-aminobutyric acid (GABA) transaminase, undergoes changes in both active site specificity and subcellular localization to switch from its renown and generic activity in core metabolism to function in a specialized metabolic pathway. Moreover, overexpression of GAME12 alone in engineered S. nigrum leaves is sufficient for de novo production of nitrogen-containing SGAs. Our results highlight how hijacking a core metabolism GABA shunt enzyme is crucial in numerous Solanum species for incorporating a nitrogen to a steroidal-specialized metabolite backbone and form defensive alkaloids.

Potato (Solanum tuberosum) is the third most important food crop in the world. Although the potato genome has been fully sequenced, functional genomics research of potato lags relative to other major food crops due primarily to the lack of a model experimental potato line. Here, we present a diploid potato line, 'Jan', which possesses all essential characteristics for facile functional genomics studies. Jan has a high level of homozygosity after seven generations of self-pollination. Jan is vigorous and highly fertile with outstanding tuber traits, high regeneration rates, and excellent transformation efficiencies. We generated a chromosome-scale genome assembly for Jan, annotated genes, and identified syntelogs relative to the potato reference genome assembly DMv6.1 to facilitate functional genomics. To miniaturize plant architecture, we developed two "mini-Jan" lines with compact and dwarf plant stature using CRISPR/Cas9-mediated mutagenesis targeting the Dwarf and Erecta genes related to growth. Mini-Jan mutants are fully fertile and will permit higher-throughput studies in limited growth chamber and greenhouse space. Thus, Jan and mini-Jan provide an outstanding model system that can be leveraged for gene editing and functional genomics research in potato.

Plant-specialized metabolism represents an inexhaustible source of active molecules, some of which have been used in human health for decades. Among these, monoterpene indole alkaloids (MIAs) include a wide range of valuable compounds with anticancer, antihypertensive, or neuroactive properties. This is particularly the case for the pachysiphine derivatives which show interesting antitumor and anti-Alzheimer activities but accumulate at very low levels in several Tabernaemontana species. Unfortunately, genome data in Tabernaemontanaceae are lacking and knowledge on the biogenesis of pachysiphine-related MIAs in planta remains scarce, limiting the prospects for the biotechnological supply of many pachysiphine-derived biopharmaceuticals. Here, we report a raw version of the toad tree (Tabernaemontana elegans) genome sequence. These new genomic resources led to the identification and characterization of a couple of genes encoding cytochrome P450 with pachysiphine synthase activity. Our phylogenomic and docking analyses highlight the different evolutionary processes that have been recruited to epoxidize the pachysiphine precursor tabersonine at a specific position and in a dedicated orientation, thus enriching our understanding of the diversification and speciation of the MIA metabolism in plants. These gene discoveries also allowed us to engineer the synthesis of MIAs in yeast through the combinatorial association of metabolic enzymes resulting in the tailor-made synthesis of non-natural MIAs. Overall, this work represents a step forward for the future supply of pachysiphine-derived drugs by microbial cell factories.

The recently described canal-1 tomato mutant, which has a variegated leaf phenotype, has been shown to affect canalization of yield. The corresponding protein is orthologous to AtSCO2 -SNOWY COTYLEDON 2, which has suggested roles in thylakoid biogenesis. Here we characterize the canal-1 mutant through a multi-omics approach, by comparing mutant to wild-type tissues. While white canal-1 leaves are devoid of chlorophyll, green leaves of the mutant appear wild-type-like, despite an impaired protein function. Transcriptomic data suggest that green mutant leaves compensate for this impaired protein function by upregulation of transcription of photosystem assembly and photosystem component genes, thereby allowing adequate photosystem establishment, which is reflected in their wild-type-like proteome. White canal-1 leaves, however, likely fail to reach a certain threshold enabling this overcompensation, and plastids get trapped in an undeveloped state, while additionally suffering from high light stress, indicated by the overexpression of ELIP homolog genes. The metabolic profile of white and to a lesser degree also green tissues revealed upregulation of amino acid levels, that was at least partially mediated by transcriptional and proteomic upregulation. These combined changes are indicative of a stress response and suggest that white tissues behave as carbon sinks. In summary, our work demonstrates the relevance of the SCO2 protein in both photosystem assembly and as a consequence in the canalization of yield.

Mexican wild diploid potato species are reproductively isolated from A-genome species, including cultivated potatoes; thus, their genomic relationships remain unknown. Solanum stoloniferum Schlechtd. et Bché. (2n = 4x = 48, AABB) is a Mexican allotetraploid species frequently used in potato breeding. We constructed a chromosome-scale assembly of the S. stoloniferum genome using PacBio long-read sequencing and Hi-C scaffolding technologies. The final assembly consisted of 1742 Mb, among which 745 Mb and 713 Mb were anchored to the 12 A-genome and 12 B-genome chromosomes, respectively. Using the RNA-seq datasets, we detected 20 994 and 19 450 genes in the A and B genomes, respectively. Among these genes, 5138 and 3594 were specific to the A and B genomes, respectively, and 15 856 were homoeologous, of which 18.6-25.4% were biasedly expressed. Structural variations such as large pericentromeric inversions were frequently found between the A- and B-genome chromosomes. A comparison of the gene sequences from 38 diverse genomes of the related Solanum species revealed that the S. stoloniferum B genome and Mexican diploid species, with the exception of S. verrucosum, were monophyletically distinct from the S. stoloniferum A genome and the other A-genome species, indicating that the Mexican diploid species share the B genome. The content and divergence of transposable elements (TEs) revealed recent bursts and transpositions of TEs after polyploidization. Thus, the S. stoloniferum genome has undergone dynamic structural differentiation and TE mobilization and reorganization to stabilize the genomic imbalance. This study provides new insights into polyploid evolution and the efficient use of allotetraploid species in potato breeding.

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.

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.

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.

Late blight (LB), caused by oomycete Phytophthora infestans, is one of the most destructive diseases of the cultivated tomato, Solanum lycopersicum. Since new and aggressive clonal lineages of P. infestans, many of which overcoming formerly effective fungicides or host resistance genes, have continued to emerge, it is crucial to identify, characterize, and utilize new sources of host resistance in tomato breeding. A recent screening of tomato germplasm identified Solanum pimpinellifolium accession PI 224710 with very strong resistance to several current P. infestans clonal lineages. The present study aimed to identify and characterize QTLs associated with LB resistance in PI 224710. Disease screening of a large F2 population (n = 1721), derived from a cross between PI 224710 and LB-susceptible tomato breeding line Fla. 8059, followed by F3 progeny testing, resulted in the identification of 43 highly-resistant and 27 highly-susceptible F2 individuals. A selective genotyping approach, using 469 non-identical SNP markers, resulted in the construction of a genetic linkage map and identification of three LB-resistance QTLs on chromosomes 6, 9 and 10 of PI 224710. A comparison of the QTLs genomic locations with the tomato physical map resulted in the identification of several candidate genes, which might be underpinning the LB-resistance QTLs in PI 224710. The identified markers associated with the LB-resistance QTLs can be utilized in breeding programs to transfer resistance from PI 224710 into tomato breeding lines and hybrid cultivars via marker-assisted breeding; they also can be used to develop near-isogenic lines for fine mapping of the QTLs.

The online version contains supplementary material available at 10.1007/s11032-024-01498-1.

In this study, we explore the interplay between the plant hormones gibberellins (GA), brassinosteroids (BR), and Indole-3-Acetic Acid (IAA) in their collective impact on plant shade avoidance elongation under varying light conditions. We focus particularly on low Red:Far-red (R:FR) light conditions achieved by supplementing the background light with FR. We characterized the tomato internode response to low R:FR and, with RNA-seq analysis, we were able to identify some of the potential regulatory hormonal pathways. Through a series of exogenous pharmacological modulations of GA, IAA, and BR, we demonstrate that GA and BR are sufficient but also necessary for inducing stem elongation under low R:FR light conditions. Intriguingly, while IAA alone shows limited effects, its combination with GA yields significant elongation, suggesting a nuanced hormonal balance. Furthermore, we unveil the complex interplay of these hormones under light with low R:FR, where the suppression of one hormone's effect can be compensated by the others. This study provides insights into the hormonal mechanisms governing plant adaptation to light, highlighting the intricate and adaptable nature of plant growth responses. Our findings have far-reaching implications for agricultural practices, offering potential strategies for optimizing plant growth and productivity in various lighting environments.

Small secreted peptides (SSPs) are essential for defense mechanisms in plant-microbe interactions, acting as danger-associated molecular patterns (DAMPs). Despite the first discovery of SSPs over three decades ago, only a limited number of SSP families, particularly within Solanaceae plants, have been identified due to inefficient approaches. This study employed comparative genomics screens with Solanaceae proteomes (tomato, tobacco, and pepper) to discover a novel SSP family, SolP. Bioinformatics analysis suggests that SolP may serve as an endogenous signal initiating the plant PTI response. Interestingly, SolP family members from tomato, tobacco, and pepper share an identical sequence (VTSNALALVNRFAD), named SlSolP12 (also referred to as NtSolP15 or CaSolP1). Biochemical and phenotypic analyses revealed that synthetic SlSolP12 peptide triggers multiple defense responses: ROS burst, MAPK activation, callose deposition, stomatal closure, and expression of immune defense genes. Furthermore, SlSolP12 enhances systemic resistance against Botrytis cinerea infection in tomato plants and interferes with classical peptides, flg22 and Systemin, which modulate the immune response. Remarkably, SolP12 activates ROS in diverse plant species, such as Arabidopsis thaliana, soybean, and rice, showing a broad spectrum of biological activities. This study provides valuable approaches for identifying endogenous SSPs and highlights SlSolP12 as a novel DAMP that could serve as a useful target for crop protection.

Fruit weight is an important agronomic trait in pepper production and is closely related to yield. At present, many quantitative trait loci (QTL) related to fruit weight have been found in pepper; however, the genes affecting fruit weight remain unknown. We analyzed the fruit weight-related quantitative traits in an intraspecific Capsicum annuum cross between the cultivated species blocky-type pepper, cv. Qiemen, and the bird pepper accession, "129-1" (Capsicum annuum var. glatriusculum), which was the wild progenitor of C. annuum. Using the QTL-seq combined with the linkage-based QTL mapping approach, QTL detection was performed; and two major effects of QTL related to fruit weight, qFW2.1 and qFW3.1, were identified on chromosomes 2 and 3. The qFW2.1 maximum explained 12.28% of the phenotypic variance observed in two F2 generations, with the maximum LOD value of 11.02, respectively; meanwhile, the qFW3.1 maximum explained 15.50% of the observed phenotypic variance in the two F2 generations, with the maximum LOD value of 11.36, respectively. qFW2.1 was narrowed down to the 1.22 Mb region using homozygous recombinant screening from BC2S2 and BC2S3 populations, while qFW3.1 was narrowed down to the 4.61Mb region. According to the transcriptome results, a total of 47 and 86 differentially expressed genes (DEGs) in the candidate regions of qFW2.1 and qFW3.1 were identified. Further, 19 genes were selected for a qRT-PCR analysis based on sequence difference combined with the gene annotation. Finally, Capana02g002938 and Capana02g003021 are the most likely candidate genes for qFW2.1, and Capana03g000903 may be a candidate gene for qFW3.1. Taken together, our results identified and fine-mapped two major QTL for fruit weight in pepper that will facilitate marker-assistant breeding for the manipulation of yield in pepper.

Rising temperature extremes during critical reproductive periods threaten the yield of major grain and fruit crops. Flowering plant reproduction depends on development of sufficient numbers of pollen grains and on their ability to generate a cellular extension, the pollen tube, which elongates through the pistil to deliver sperm cells to female gametes for double fertilization. These critical phases of the life cycle are sensitive to temperature and limit productivity under high temperature (HT). Previous studies have investigated the effects of HT on pollen development, but little is known about how HT applied during the pollen tube growth phase affects fertility. Here, we used tomato as a model fruit crop to determine how HT affects the pollen tube growth phase, taking advantage of cultivars noted for fruit production in exceptionally hot growing seasons. We found that exposure to HT solely during the pollen tube growth phase limits fruit biomass and seed set more significantly in thermosensitive cultivars than in thermotolerant cultivars. Importantly, we found that pollen tubes from the thermotolerant Tamaulipas cultivar have enhanced growth in vivo and in vitro under HT. Analysis of the pollen tube transcriptome's response to HT allowed us to develop hypotheses for the molecular basis of cellular thermotolerance in the pollen tube and we define two response modes (enhanced induction of stress responses, and higher basal levels of growth pathways repressed by heat stress) associated with reproductive thermotolerance. Importantly, we define key components of the pollen tube stress response identifying enhanced ROS homeostasis and pollen tube callose synthesis and deposition as important components of reproductive thermotolerance in Tamaulipas. Our work identifies the pollen tube growth phase as a viable target to enhance reproductive thermotolerance and delineates key pathways that are altered in crop varieties capable of fruiting under HT conditions.

An enduring question in evolutionary biology concerns the degree to which episodes of convergent trait evolution depend on the same genetic programs, particularly over long timescales. In this work, we genetically dissected repeated origins and losses of prickles-sharp epidermal projections-that convergently evolved in numerous plant lineages. Mutations in a cytokinin hormone biosynthetic gene caused at least 16 independent losses of prickles in eggplants and wild relatives in the genus Solanum. Homologs underlie prickle formation across angiosperms that collectively diverged more than 150 million years ago, including rice and roses. By developing new Solanum genetic systems, we leveraged this discovery to eliminate prickles in a wild species and an indigenously foraged berry. Our findings implicate a shared hormone activation genetic program underlying evolutionarily widespread and recurrent instances of plant morphological innovation.

Early blight and brown leaf spot are often cited as the most problematic pathogens of tomato in many agricultural regions. Their causal agents are Alternaria spp., a genus of Ascomycota containing numerous necrotrophic pathogens. Breeding programs have yielded quantitatively resistant commercial cultivars, but fungicide application remains necessary to mitigate the yield losses. A major hindrance to resistance breeding is the complexity of the genetic determinants of resistance and susceptibility. In the absence of sufficiently resistant germplasm, we sequenced the transcriptomes of Heinz 1706 tomatoes treated with strongly virulent and weakly virulent isolates of Alternaria spp. 3 h post infection. We expanded existing functional gene annotations in tomato and using network statistics, we analyzed the transcriptional modules associated with defense and susceptibility.

The induced responses are very distinct. The weakly virulent isolate induced a defense response of calcium-signaling, hormone responses, and transcription factors. These defense-associated processes were found in a single transcriptional module alongside secondary metabolite biosynthesis genes, and other defense responses. Co-expression and gene regulatory networks independently predicted several D clade ethylene response factors to be early regulators of the defense transcriptional module, as well as other transcription factors both known and novel in pathogen defense, including several JA-associated genes. In contrast, the strongly virulent isolate elicited a much weaker response, and a separate transcriptional module bereft of hormone signaling.

Our findings have predicted major defense regulators and several targets for downstream functional analyses. Combined with our improved gene functional annotation, they suggest that defense is achieved through induction of Alternaria-specific immune pathways, and susceptibility is mediated by modulating hormone responses. The implication of multiple specific clade D ethylene response factors and upregulation of JA-associated genes suggests that host defense in this pathosystem involves ethylene response factors to modulate jasmonic acid signaling.

Viroids, one of the smallest known infectious agents, induce symptoms of varying severity, ranging from latent to severe, based on the combination of viroid isolates and host plant species. Because viroids are transmissible between plant species, asymptomatic viroid-infected plants may serve as latent sources of infection for other species that could exhibit severe symptoms, occasionally leading to agricultural and economic losses. Therefore, predicting the symptoms induced by viroids in host plants without biological experiments could remarkably enhance control measures against viroid damage. Here, we developed an algorithm using unsupervised machine learning to predict the severity of disease symptoms caused by viroids (e.g., potato spindle tuber viroid; PSTVd) in host plants (e.g., tomato). This algorithm, mimicking the RNA silencing mechanism thought to be linked to viroid pathogenicity, requires only the genome sequences of the viroids and host plants. It involves three steps: alignment of synthetic short sequences of the viroids to the host plant genome, calculation of the alignment coverage, and clustering of the viroids based on coverage using UMAP and DBSCAN. Validation through inoculation experiments confirmed the effectiveness of the algorithm in predicting the severity of disease symptoms induced by viroids. As the algorithm only requires the genome sequence data, it may be applied to any viroid and plant combination. These findings underscore a correlation between viroid pathogenicity and the genome sequences of viroid isolates and host plants, potentially aiding in the prevention of viroid outbreaks and the breeding of viroid-resistant crops.

Vegetables, as indispensable non-staple foods in people's daily diet, provide a variety of essential vitamins, minerals, and other nutrients, as well as special phytochemicals, which are recognized as functional components for human nutritional balance or medicinal purposes [...].

The domain of unknown function 668 (DUF668) is a gene family that may play a key role in plant growth and development as well as in responding to adversity coercion stresses. However, the DUF668 gene family has not yet been well identified and characterized in tomato. In this study, a total of nine putative SlDUF668 genes were identified in tomato, distributed on six chromosomes. Phylogenetic analyses revealed that SlDUF668 proteins were classified into two major groups. Members within the same group largely displayed analogous gene structure and conserved motif compositions. Several cis-elements were exhibited in the upstream sequences of the SlDUF668 genes, including elements implicated in plant growth and development processes, abiotic stress and hormone responses. Further, the study assessed the expression patterns of the SlDUF668 gene family in various tomato tissues, five plant hormones treatments, three abiotic stresses using qRT-PCR. The SlDUF668 genes expressed ubiquitously in various tissues, and five genes (SlDUF668-04, SlDUF668-06, SlDUF668-07, SlDUF668-08 and SlDUF668-09) showed tissue specificity. And SlDUF668 genes responded to abiotic stresses such as salt, drought and cold to varying degrees. Overall, our study provided a base for the tomato DUF668 gene family and laid a foundation for further understanding the functional characteristics of DUF668 genes in tomato plants.

Potato (Solanum tuberosum) is the most popular tuber crop and a model organism. A variety of gene models for potato exist, and despite frequent updates, they are not unified. This hinders the comparison of gene models across versions, limits the ability to reuse experimental data without significant re-analysis, and leads to missing or wrongly annotated genes. Here, we unify the recent potato double monoploid v4 and v6 gene models by developing an automated merging protocol, resulting in a Unified poTato genome model (UniTato). We subsequently established an Apollo genome browser (unitato.nib.si) that enables public access to UniTato and further community-based curation. We demonstrate how the UniTato resource can help resolve problems with missing or misplaced genes and can be used to update or consolidate a wider set of gene models or genome information. The automated protocol, genome annotation files, and a comprehensive translation table are provided at github.com/NIB-SI/unitato.

Ehretia macrophylla Wall, known as wild loquat, is an ecologically, economically, and medicinally significant tree species widely grown in China, Japan, Vietnam, and Nepal. In this study, we have successfully generated a haplotype-resolved chromosome-scale genome assembly of E. macrophylla by integrating PacBio HiFi long-reads, Illumina short-reads, and Hi-C data. The genome assembly consists of two haplotypes, with sizes of 1.82 Gb and 1.58 Gb respectively, and contig N50 lengths of 28.11 Mb and 21.57 Mb correspondingly. Additionally, 99.41% of the assembly was successfully anchored into 40 pseudo-chromosomes. We predicted 58,886 protein-coding genes, of which 99.60% were functionally annotated from databases. We furthermore detected 2.65 Gb repeat sequences, 659,290 rRNAs, 4,931 tRNAs and 4,688 other ncRNAs. The high-quality assembly of the genome offers a solid basis for furthering the fields of molecular breeding and functional genomics of E. macrophylla.

Potato and its wild relatives are distributed mainly in the Mexican highlands and central Andes of South America. The South American A-genome species, including cultivated potatoes, are reproductively isolated from Mexican diploid species. Whole-genome sequencing has disclosed genome structure and similarity, mostly in cultivated potatoes and their closely related species. In this study, we generated a chromosome-scale assembly of the genome of a Mexican diploid species, Solanum bulbocastanum Dun., using PacBio long-read sequencing, optical mapping, and Hi-C scaffolding technologies. The final sequence assembly consisted of 737.9 Mb, among which 647.0 Mb were anchored to the 12 chromosomes. Compared with chromosome-scale assemblies of S. lycopersicum (tomato), S. etuberosum (non-tuber-bearing species with E-genome), S. verrucosum, S. chacoense, S. multidissectum, and S. phureja (all four are A-genome species), the S. bulbocastnum genome was the shortest. It contained fewer transposable elements (56.2%) than A-genome species. A cluster analysis was performed based on pairwise ratios of syntenic regions among the seven chromosome-scale assemblies, showing that the A-genome species were first clustered as a distinct group. Then, this group was clustered with S. bulbocastanum. Sequence similarity in 1,624 single-copy orthologous gene groups among 36 Solanum species and clones separated S. bulbocastanum as a specific group, including other Mexican diploid species, from the A-genome species. Therefore, the S. bulbocastanum genome differs in genome structure and gene sequences from the A-genome species. These findings provide important insights into understanding and utilizing the genetic diversity of S. bulbocastanum and the other Mexican diploid species in potato breeding.

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.

Next-generation sequencing (NGS) library construction often involves using restriction enzymes to decrease genome complexity, enabling versatile polymorphism detection in plants. However, plant leaves frequently contain impurities, such as polyphenols, necessitating DNA purification before enzymatic reactions. To overcome this problem, we developed a PCR-based method for expeditious NGS library preparation, offering flexibility in number of detected polymorphisms. By substituting a segment of the simple sequence repeat sequence in the MIG-seq primer set (MIG-seq being a PCR method enabling library construction with low-quality DNA) with degenerate oligonucleotides, we introduced variability in detectable polymorphisms across various crops. This innovation, named degenerate oligonucleotide primer MIG-seq (dpMIG-seq), enabled a streamlined protocol for constructing dpMIG-seq libraries from unpurified DNA, which was implemented stably in several crop species, including fruit trees. Furthermore, dpMIG-seq facilitated efficient lineage selection in wheat and enabled linkage map construction and quantitative trait loci analysis in tomato, rice, and soybean without necessitating DNA concentration adjustments. These findings underscore the potential of the dpMIG-seq protocol for advancing genetic analyses across diverse plant species.