New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication

The genetic architecture of the pepper metabolome and the biosynthesis of its signature capsianoside metabolites

Capsicum (pepper) is among the most economically important species worldwide, and its fruits accumulate specialized metabolites with essential roles in plant environmental interaction and human health benefits as well as in conferring their unique taste. However, the genetics underlying differences in metabolite presence/absence and/or accumulation remain largely unknown. In this study, we carried out a genome-wide association study as well as generating and characterizing a novel backcross inbred line mapping population to determine the genetic architecture of the pepper metabolome. This genetic analysis provided over 1,000 metabolic quantitative trait loci (mQTL) for over 250 annotated metabolites. We identified 92 candidate genes involved in various mQTLs. Among the identified loci, we described and validated a gene cluster of eleven UDP-glycosyltransferases (UGTs) involved in monomeric capsianoside biosynthesis. We additionally constructed the gene-by-gene-based biosynthetic pathway of pepper capsianoside biosynthesis, including both core and decorative reactions. Given that one of these decorative pathways, namely the glycosylation of acyclic diterpenoid glycosides, contributes to plant resistance, these data provide new insights and breeding resources for pepper. They additionally provide a blueprint for the better understanding of the biosynthesis of species-specific natural compounds in general.

NRC Immune receptor networks show diversified hierarchical genetic architecture across plant lineages

Plants' complex immune systems include nucleotide-binding domain and leucine-rich repeat-containing (NLR) proteins, which help recognize invading pathogens. In solanaceous plants, the NRC (NLR required for cell death) family includes helper NLRs that form a complex genetic network with multiple sensor NLRs to provide resistance against pathogens. However, the evolution and function of NRC networks outside solanaceous plants are currently unclear. Here, we conducted phylogenomic and macroevolutionary analyses comparing NLRs identified from different asterid lineages and found that NRC networks expanded significantly in most lamiids but not in Ericales and campanulids. Using transient expression assays in Nicotiana benthamiana, we showed that NRC networks are simple in Ericales and campanulids, but have high complexity in lamiids. Phylogenetic analyses grouped the NRC helper NLRs into three NRC0 subclades that are conserved, and several family-specific NRC subclades of lamiids that show signatures of diversifying selection. Functional analyses revealed that members of the NRC0 subclades are partially interchangeable, whereas family-specific NRC members in lamiids lack interchangeability. Our findings highlight the distinctive evolutionary patterns of the NRC networks in asterids and provide potential insights into transferring disease resistance across plant lineages.

Full-length transcriptome sequencing of pepper fruit during development and construction of a transcript variation database

Chili pepper is an important spice and a model plant for fruit development studies. Large-scale omics information on chili pepper plant development continues to be gathered for understanding development as well as capsaicin biosynthesis. In this study, a full-spectrum transcriptome data of eight chili pepper tissues at five growth stages using the Oxford Nanopore long-read sequencing approach was generated. Of the 485 351 transcripts, 35 336 were recorded as reference transcripts (genes), while 450 015 were novel including coding, lnc, and other non-coding RNAs. These novel transcripts belonged to unknown/intergenic (347703), those retained introns (26336), and had multi-exons with at least one junction match (20333). In terms of alternative splicing, retained intron had the highest proportion (14795). The number of tissue-specific expressed transcripts ranged from 22 925 (stem) to 40 289 (flower). The expression changes during fruit and placenta development are discussed in detail. Integration of gene expression and capsaicin content quantification throughout the placental development clarifies that capsaicin biosynthesis in pepper is mainly derived from valine, leucin, and isoleucine degradation as well as citrate cycle and/or pyrimidine metabolism pathways. Most importantly, a user-friendly Pepper Full-Length Transcriptome Variation Database (PFTVD 1.0) (http://pepper-database.cn/) has been developed. PFTVD 1.0 provides transcriptomics and genomics information and allows users to analyse the data using various tools implemented. This work highlights the potential of long-read sequencing to discover novel genes and transcripts and their diversity in plant developmental biology.

Comparative Genomics Screens Identify a Novel Small Secretory Peptide, SlSolP12, which Activates Both Local and Systemic Immune Response in Tomatoes and Exhibits Broad-Spectrum Activity

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.

A chromosome-scale assembly of Brassica carinata (BBCC) accession HC20 containing resistance to multiple pathogens and an early generation assessment of introgressions into B. juncea (AABB)

Brassica carinata (BBCC) commonly referred to as Ethiopian mustard is a natural allotetraploid containing the genomes of Brassica nigra (BB) and Brassica oleracea (CC). It is an oilseed crop endemic to the northeastern regions of Africa. Although it is under limited cultivation, B. carinata is valuable as it is resistant/highly tolerant to most of the pathogens affecting widely cultivated Brassica species of the U's triangle. We report a chromosome-scale genome assembly of B. carinata accession HC20 using long-read Oxford Nanopore sequencing and Bionano optical maps. The assembly has a scaffold N50 of ~39.8 Mb and covers ~1.11 Gb of the genome. We compared the long-read genome assemblies of the U's triangle species and found extensive gene collinearity between the diploids and allopolyploids with no evidence of major gene losses. Therefore, B. juncea (AABB), B. napus (AACC), and B. carinata can be regarded as strict allopolyploids. We cataloged the nucleotide-binding and leucine-rich repeat immune receptor (NLR) repertoire of B. carinata and, identified 465 NLRs, and compared these with the NLRs in the other Brassica species. We investigated the extent and nature of early-generation genomic interactions between the constituent genomes of B. carinata and B. juncea in interspecific crosses between the two species. Besides the expected recombination between the constituent B genomes, extensive homoeologous exchanges were observed between the A and C genomes. Interspecific crosses, therefore, can be used for transferring disease resistance from B. carinata to B. juncea and broadening the genetic base of the two allotetraploid species.

Involvement of the GH38 Family Exoglycosidase α-Mannosidase in Strawberry Fruit Ripening

Exoglycosidase enzymes hydrolyze the N-glycosylations of cell wall enzymes, releasing N-glycans that act as signal molecules and promote fruit ripening. Vesicular exoglycosidase α-mannosidase enzymes of the GH38 family (EC 3.2.1.24; α-man) hydrolyze N-glycans in non-reduced termini. Strawberry fruit (Fragaria × ananassa) is characterized by rapid softening as a result of cell wall modifications during the fruit ripening process. Enzymes acting on cell wall polysaccharides explain the changes in fruit firmness, but α-man has not yet been described in F. × ananassa, meaning that the indirect effects of N-glycan removal on its fruit ripening process are unknown. The present study identified 10 GH38 α-man sequences in the F. × ananassa genome with characteristic conserved domains and key residues. A phylogenetic tree built with the neighbor-joining method and three groups of α-man established, of which group I was classified into three subgroups and group III contained only Poaceae spp. sequences. The real-time qPCR results demonstrated that FaMAN genes decreased during fruit ripening, a trend mirrored by the total enzyme activity from the white to ripe stages. The analysis of the promoter regions of these FaMAN genes was enriched with ripening and phytohormone response elements, and contained cis-regulatory elements related to stress responses to low temperature, drought, defense, and salt stress. This study discusses the relevance of α-man in fruit ripening and how it can be a useful target to prolong fruit shelf life.

Two telomere-to-telomere gapless genomes reveal insights into Capsicum evolution and capsaicinoid biosynthesis

Chili pepper (Capsicum) is known for its unique fruit pungency due to the presence of capsaicinoids. The evolutionary history of capsaicinoid biosynthesis and the mechanism of their tissue specificity remain obscure due to the lack of high-quality Capsicum genomes. Here, we report two telomere-to-telomere (T2T) gap-free genomes of C. annuum and its wild nonpungent relative C. rhomboideum to investigate the evolution of fruit pungency in chili peppers. We precisely delineate Capsicum centromeres, which lack high-copy tandem repeats but are extensively invaded by CRM retrotransposons. Through phylogenomic analyses, we estimate the evolutionary timing of capsaicinoid biosynthesis. We reveal disrupted coding and regulatory regions of key biosynthesis genes in nonpungent species. We also find conserved placenta-specific accessible chromatin regions, which likely allow for tissue-specific biosynthetic gene coregulation and capsaicinoid accumulation. These T2T genomic resources will accelerate chili pepper genetic improvement and help to understand Capsicum genome evolution.

Genetic tapestry of Capsicum fruit colors: a comparative analysis of four cultivated species

KEY MESSAGE

Genome-wide association study of color spaces across the four cultivated Capsicum spp. revealed a shared set of genes influencing fruit color, suggesting mechanisms and pathways across Capsicum species are conserved during the speciation. Notably, Cytochrome P450 of the carotenoid pathway, MYB transcription factor, and pentatricopeptide repeat-containing protein are the major genes responsible for fruit color variation across the Capsicum species. Peppers (Capsicum spp.) rank among the most widely consumed spices globally. Fruit color, serving as a determinant for use in food colorants and cosmeceuticals and an indicator of nutritional contents, significantly influences market quality and price. Cultivated Capsicum species display extensive phenotypic diversity, especially in fruit coloration. Our study leveraged the genetic variance within four Capsicum species (Capsicum baccatum, Capsicum chinense, Capsicum frutescens, and Capsicum annuum) to elucidate the genetic mechanisms driving color variation in peppers and related Solanaceae species. We analyzed color metrics and chromatic attributes (Red, Green, Blue, L*, a*, b*, Luminosity, Hue, and Chroma) on samples cultivated over six years (2015-2021). We resolved genomic regions associated with fruit color diversity through the sets of SNPs obtained from Genotyping by Sequencing (GBS) and genome-wide association study (GWAS) with a Multi-Locus Mixed Linear Model (MLMM). Significant SNPs with FDR correction were identified, within the Cytochrome P450, MYB-related genes, Pentatricopeptide repeat proteins, and ABC transporter family were the most common among the four species, indicating comparative evolution of fruit colors. We further validated the role of a pentatricopeptide repeat-containing protein (Chr01:31,205,460) and a cytochrome P450 enzyme (Chr08:45,351,919) via competitive allele-specific PCR (KASP) genotyping. Our findings advance the understanding of the genetic underpinnings of Capsicum fruit coloration, with developed KASP assays holding potential for applications in crop breeding and aligning with consumer preferences. This study provides a cornerstone for future research into exploiting Capsicum's diverse fruit color variation.

Elucidating the callus-to-shoot-forming mechanism in Capsicum annuum 'Dempsey' through comparative transcriptome analyses

BACKGROUND

The formation of shoots plays a pivotal role in plant organogenesis and productivity. Despite its significance, the underlying molecular mechanism of de novo regeneration has not been extensively elucidated in Capsicum annuum 'Dempsey', a bell pepper cultivar. To address this, we performed a comparative transcriptome analysis focusing on the differential expression in C. annuum 'Dempsey' shoot, callus, and leaf tissue. We further investigated phytohormone-related biological processes and their interacting genes in the C. annuum 'Dempsey' transcriptome based on comparative transcriptomic analysis across five species.

RESULTS

We provided a comprehensive view of the gene networks regulating shoot formation on the callus, revealing a strong involvement of hypoxia responses and oxidative stress. Our comparative transcriptome analysis revealed a significant conservation in the increase of gene expression patterns related to auxin and defense mechanisms in both callus and shoot tissues. Consequently, hypoxia response and defense mechanism emerged as critical regulators in callus and shoot formation in C. annuum 'Dempsey'. Current transcriptome data also indicated a substantial decline in gene expression linked to photosynthesis within regenerative tissues, implying a deactivation of the regulatory system governing photosynthesis in C. annuum 'Dempsey'.

CONCLUSION

Coupled with defense mechanisms, we thus considered spatial redistribution of auxin to play a critical role in the shoot morphogenesis via primordia outgrowth. Our findings shed light on shoot formation mechanisms in C. annuum 'Dempsey' explants, important information for regeneration programs, and have broader implications for precise molecular breeding in recalcitrant crops.

Mapping and Omics Integration: Towards Precise Rice Disease Resistance Breeding

Rice (Oryza sativa), as a staple crop feeding a significant portion of the global population, particularly in Asian countries, faces constant threats from various diseases jeopardizing global food security. A precise understanding of disease resistance mechanisms is crucial for developing resilient rice varieties. Traditional genetic mapping methods, such as QTL mapping, provide valuable insights into the genetic basis of diseases. However, the complex nature of rice diseases demands a holistic approach to gain an accurate knowledge of it. Omics technologies, including genomics, transcriptomics, proteomics, and metabolomics, enable a comprehensive analysis of biological molecules, uncovering intricate molecular interactions within the rice plant. The integration of various mapping techniques using multi-omics data has revolutionized our understanding of rice disease resistance. By overlaying genetic maps with high-throughput omics datasets, researchers can pinpoint specific genes, proteins, or metabolites associated with disease resistance. This integration enhances the precision of disease-related biomarkers with a better understanding of their functional roles in disease resistance. The improvement of rice breeding for disease resistance through this integration represents a significant stride in agricultural science because a better understanding of the molecular intricacies and interactions underlying disease resistance architecture leads to a more precise and efficient development of resilient and productive rice varieties. In this review, we explore how the integration of mapping and omics data can result in a transformative impact on rice breeding for enhancing disease 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.

The landscape of abiotic and biotic stress-responsive splice variants with deep RNA-seq datasets in hot pepper

Alternative splicing (AS) is a widely observed phenomenon in eukaryotes that plays a critical role in development and stress responses. In plants, the large number of RNA-seq datasets in response to different environmental stressors can provide clues for identification of condition-specific and/or common AS variants for preferred agronomic traits. We report RNA-seq datasets (350.7 Gb) from Capsicum annuum inoculated with one of three bacteria, one virus, or one oomycete and obtained additional existing transcriptome datasets. In this study, we investigated the landscape of AS in response to environmental stressors, signaling molecules, and tissues from 425 total samples comprising 841.49 Gb. In addition, we identified genes that undergo AS under specific and shared stress conditions to obtain potential genes that may be involved in enhancing tolerance to stressors. We uncovered 1,642,007 AS events and identified 4,354 differential alternative splicing genes related to environmental stressors, tissues, and signaling molecules. This information and approach provide useful data for basic-research focused on enhancing tolerance to environmental stressors in hot pepper or establishing breeding programs.

Genetic characterization of a locus responsible for low pungency using EMS-induced mutants in Capsicum annuum L

KEY MESSAGE

The pepper mutants ('221-2-1a' and '1559-1-2h') with very low pungency were genetically characterized. The Pun4 locus, responsible for the reduced pungency of the mutant fruits, was localized to a 208 Mb region on chromosome 6. DEMF06G16460, encoding 3-ketoacyl-CoA synthase, was proposed as a strong candidate gene based on the genetic analyses of bulked segregants, DEG, and expression analyses. Capsaicinoids are unique alkaloids present in pepper (Capsicum spp.), synthesized through the condensation of by-products from the phenylpropanoid and branched-chain fatty acid pathways, and accumulating in the placenta. In this study, we characterized two allelic ethyl methanesulfonate-induced mutant lines with extremely low pungency ('221-2-1a' and '1559-1-2h'). These mutants, derived from the pungent Korean landrace 'Yuwolcho,' exhibited lower capsaicinoid content than Yuwolcho but still contained a small amount of capsaicinoid with functional capsaicinoid biosynthetic genes. Genetic crosses between the mutants and Yuwolcho or pungent lines indicated that a single recessive mutation was responsible for the low-pungency phenotype of mutant 221-2-1a; we named the causal locus Pungency 4 (Pun4). To identify Pun4, we combined genome-wide polymorphism analysis and transcriptome analysis with bulked-segregant analysis. We narrowed down the location of Pun4 to a 208-Mb region on chromosome 6 containing five candidate genes, of which DEMF06G16460, encoding a 3-ketoacyl-CoA synthase associated with branched-chain fatty acid biosynthesis, is the most likely candidate for Pun4. The expression of capsaicinoid biosynthetic genes in placental tissues in Yuwolcho and the mutant was consistent with the branched-chain fatty acid pathway playing a pivotal role in the lower pungency observed in the mutant. We also obtained a list of differentially expressed genes in placental tissues between the mutant and Yuwolcho, from which we selected candidate genes using gene co-expression analysis. In summary, we characterized the capsaicinoid biosynthesis-related locus Pun4 through integrated of genetic, genomic, and transcriptome analyses. These findings will contribute to our understanding of capsaicinoid biosynthesis in pepper.

Graft incompatibility between pepper and tomato can be attributed to genetic incompatibility between diverged immune systems

Graft compatibility is the capacity of two plants to form cohesive vascular connections. Tomato and pepper are incompatible graft partners; however, the underlying cause of graft rejection between these two species remains unknown.We diagnosed graft incompatibility between tomato and diverse pepper varieties based on weakened biophysical stability, decreased growth, and persistent cell death using trypan blue and TUNEL assays. Transcriptomic analysis of cell death in the junction was performed using RNA-sequencing, and molecular signatures for incompatible graft response were characterized based on meta-transcriptomic comparisons with other biotic processes.We show that tomato is broadly incompatible with diverse pepper cultivars. These incompatible graft partners activate prolonged transcriptional changes that are highly enriched for defense processes. Amongst these processes was broad NLR upregulation and hypersensitive response. Using transcriptomic datasets for a variety of biotic stress treatments, we identified a significant overlap in the genetic profile of incompatible grafting and plant parasitism. In addition, we found over 1000 genes that are uniquely upregulated in incompatible grafts.Based on NLR overactivity, DNA damage, and prolonged cell death we have determined that tomato and pepper graft incompatibility is likely caused by a form of genetic incompatibility, which triggers a hyperimmune-response.

Identification and expression analyses of B3 genes reveal lineage-specific evolution and potential roles of REM genes in pepper

BACKGROUND

The B3 gene family, one of the largest plant-specific transcription factors, plays important roles in plant growth, seed development, and hormones. However, the B3 gene family, especially the REM subfamily, has not been systematically and functionally studied.

RESULTS

In this study, we performed genome-wide re-annotation of B3 genes in five Solanaceae plants, Arabidopsis thaliana, and Oryza sativa, and finally predicted 1,039 B3 genes, including 231 (22.2%) newly annotated genes. We found a striking abundance of REM genes in pepper species (Capsicum annuum, Capsicum baccatum, and Capsicum chinense). Comparative motif analysis revealed that REM and other subfamilies (ABI3/VP1, ARF, RAV, and HSI) consist of different amino acids. We verified that the large number of REM genes in pepper were included in the specific subgroup (G8) through the phylogenetic analysis. Chromosome location and evolutionary analyses suggested that the G8 subgroup genes evolved mainly via a pepper-specific recent tandem duplication on chromosomes 1 and 3 after speciation between pepper and other Solanaceae. RNA-seq analyses suggested the potential functions of REM genes under salt, heat, cold, and mannitol stress conditions in pepper (C. annuum).

CONCLUSIONS

Our study provides evolutionary and functional insights into the REM gene family in pepper.

An evolutionary view of vanillylamine synthase pAMT, a key enzyme of capsaicinoid biosynthesis pathway in chili pepper

Pungent capsaicinoid is synthesized only in chili pepper (Capsicum spp.). The production of vanillylamine from vanillin is a unique reaction in the capsaicinoid biosynthesis pathway. Although putative aminotransferase (pAMT) has been isolated as the vanillylamine synthase gene, it is unclear how Capsicum acquired pAMT. Here, we present a phylogenetic overview of pAMT and its homologs. The Capsicum genome contained 5 homologs, including pAMT, CaGABA-T1, CaGABA-T3, and two pseudogenes. Phylogenetic analysis indicated that pAMT is a member of the Solanaceae cytoplasmic GABA-Ts. Comparative genome analysis found that multiple copies of GABA-T exist in a specific Solanaceae genomic region, and the cytoplasmic GABA-Ts other than pAMT are located in the region. The cytoplasmic GABA-T was phylogenetically close to pseudo-GABA-T harboring a plastid transit peptide (pseudo-GABA-T3). This suggested that Solanaceae cytoplasmic GABA-Ts occurred via duplication of a chloroplastic GABA-T ancestor and subsequent loss of the plastid transit signal. The cytoplasmic GABA-T may have been translocated from the specific Solanaceae genomic region during Capsicum divergence, resulting in the current pAMT locus. A recombinant protein assay demonstrated that pAMT had higher vanillylamine synthase activity than those of other plant GABA-Ts. pAMT was expressed exclusively in the placental septum of mature green fruit, whereas tomato orthologs SlGABA-T2/4 exhibit a ubiquitous expression pattern in plants. These findings suggested that both the increased catalytic efficiency and transcriptional changes in pAMT may have contributed to establish vanillylamine synthesis in the capsaicinoid biosynthesis pathway. This study provides insights into the establishment of pungency in the evolution of chili peppers.

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.

Chromosome level genome assembly of endangered medicinal plant Anisodus tanguticus

Anisodus tanguticus is a medicinal herb that belongs to the Anisodus genus of the Solanaceae family. This endangered herb is mainly distributed in Qinghai-Tibet Plateau. In this study, we combined the Illumina short-read, Nanopore long-read and high-throughput chromosome conformation capture (Hi-C) sequencing technologies to de novo assemble the A. tanguticus genome. A high-quality chromosomal-level genome assembly was obtained with a genome size of 1.26 Gb and a contig N50 of 25.07 Mb. Of the draft genome sequences, 97.47% were anchored to 24 pseudochromosomes with a scaffold N50 of 51.28 Mb. In addition, 842.14 Mb of transposable elements occupying 66.70% of the genome assembly were identified and 44,252 protein-coding genes were predicted. The genome assembly of A. tanguticus will provide genetic repertoire to understand the adaptation strategy of Anisodus species in the plateau, which will further promote the conservation of endangered A. tanguticus resources.

Multi-environment association study highlights candidate genes for robust agronomic quantitative trait loci in a novel worldwide Capsicum core collection

Investigating crop diversity through genome-wide association studies (GWAS) on core collections helps in deciphering the genetic determinants of complex quantitative traits. Using the G2P-SOL project world collection of 10 038 wild and cultivated Capsicum accessions from 10 major genebanks, we assembled a core collection of 423 accessions representing the known genetic diversity. Since complex traits are often highly dependent upon environmental variables and genotype-by-environment (G × E) interactions, multi-environment GWAS with a 10 195-marker genotypic matrix were conducted on a highly diverse subset of 350 Capsicum annuum accessions, extensively phenotyped in up to six independent trials from five climatically differing countries. Environment-specific and multi-environment quantitative trait loci (QTLs) were detected for 23 diverse agronomic traits. We identified 97 candidate genes potentially implicated in 53 of the most robust and high-confidence QTLs for fruit flavor, color, size, and shape traits, and for plant productivity, vigor, and earliness traits. Investigating the genetic architecture of agronomic traits in this way will assist the development of genetic markers and pave the way for marker-assisted selection. The G2P-SOL pepper core collection will be available upon request as a unique and universal resource for further exploitation in future gene discovery and marker-assisted breeding efforts by the pepper community.

Characterization of the FLA Gene Family in Tomato (Solanum lycopersicum L.) and the Expression Analysis of SlFLAs in Response to Hormone and Abiotic Stresses

Fasciclin-like arabinogalactan proteins (FLAs), a subclass of arabinogalactan proteins (AGPs), participate in mediating plant growth, development, and response to abiotic stress. However, the characterization and function of FLAs in tomato are currently unknown. In this study, members of the tomato FLA family are characterized and analyzed in relation to their response to phytohormonal and abiotic stresses. The results show that a total of 24 FLA members were characterized in tomato. The structural domain analysis showed that these members have a high protein similarity. The expression profiles of different tissues indicated that the genes of most members of the tomato FLA gene family are highly expressed in roots, but to a lower extent in fruits. qRT-PCR analysis revealed that all 24 tomato FLA genes are responsive to ABA and MeJA. SlFLAs showed a positive response to salt and cold stress. SlFLA1, SlFLA12, and SlFLA14 are significantly induced under darkness. SlFLA1 and SlFLA3 are significantly induced under drought stress. This study provides a basis for a further understanding of the role of tomato FLA homologous genes in plant response to abiotic stress and lays the foundation for further research on the function of FLAs in tomato.

Expression of Genes Involved in Anthracnose Resistance in Chili (Capsicum baccatum) 'PBC80'-Derived Recombinant Inbred Lines

Chili anthracnose has long been a threat to chili production worldwide. Capsicum baccatum 'PBC80' has been identified as a source of resistance to anthracnose. Recently, a QTL for ripe fruit resistance from 'PBC80'-derived RILs was located on chromosome 4 (123 Mb) and contained over 80 defense-related genes. To identify the genes most related to anthracnose resistance, a fine map of the QTL region was developed using single-marker analysis. Nine genes were selected from the new QTL (1.12 Mb) to study their expression after being challenged with Colletotrichum scovillei 'MJ5' in two different RIL genotypes (Resistance/Resistance or R/R and Susceptible/Susceptible or S/S) at 0, 6 and 12 h. Of the nine genes, LYM2, CQW23_09597, CLF, NFXL1, and PR-14 were significantly up-regulated, compared to the control, in the R/R genotype. ERF was up-regulated in both chili genotypes. However, the expression was relatively and constantly low in the S/S genotype. Most up-regulated genes reached the highest peak (2.3-4.5 fold) at 6 h, except for ERF, which had the highest peak at 12 h (6.4 fold). The earliest and highest expressed gene was a pathogen receptor, LYM2.

The Capsicum terpenoid biosynthetic module is affected by spider-mite herbivory

In response to herbivory, Capsicum annuum leaves adapt their specialized metabolome that may protect the plant against herbivore feeding either directly or indirectly through volatile metabolites acting as cues for natural enemies of the herbivore. The volatile blend of spider-mite infested leaves differs from non-challenged leaves predominantly by a higher contribution of mono- and sesquiterpenes. In addition to these terpenoids released into the headspace, the terpenoid composition of the leaves alters upon herbivory. All this suggests an important role for terpenoids and their biosynthetic machinery in the defence against herbivory. Here, we show that the C. annuum genome contains a terpene synthase (TPS) gene family of 103 putative members of which structural analysis revealed that 27 encode functional enzymes. Transcriptome analysis showed that several TPS loci were differentially expressed upon herbivory in leaves of two C. annuum genotypes, that differ in susceptibility towards spider mites. The relative expression of upstream biosynthetic genes from the mevalonate and the methylerythritol phosphate pathway also altered upon herbivory, revealing a shift in the metabolic flux through the terpene biosynthetic module. The expression of multiple genes potentially acting downstream of the TPSs, including cytochrome P450 monooxygenases, UDP-glucosyl transferases, and transcription factors strongly correlated with the herbivory-induced TPS genes. A selection of herbivory-induced TPS genes was functionally characterized through heterologous expression and the products that these enzymes catalysed matched with the volatile and non-volatile terpenoids induced in response to herbivory.

Expression of alcohol acyltransferase is a potential determinant of fruit volatile ester variations in Capsicum

KEY MESSAGE

The transcript level of alcohol acyltransferase 1 (AAT1) may be the main factor influencing the variations in volatile esters that characterizing the fruity/exotic aroma of pepper fruit. Volatile esters are key components for characterizing the fruity/exotic aroma of pepper (Capsicum spp.) fruit. In general, the volatile ester content in the fruit is the consequence of a delicate balance between their synthesis by alcohol acyltransferases (AATs) and degradation by carboxylesterases (CXEs). However, the precise role of these families of enzymes with regard to volatile ester content remains unexplored in Capsicum. In this study, we found that the volatile ester content was relatively low in C. annuum and much higher in C. chinense, particularly in pungent varieties. Additionally, fruits collected from multiple non-pungent C. chinense varieties, which harbor loss-of-function mutations in capsaicinoid biosynthetic genes, acyltransferase (Pun1), putative aminotransferase (pAMT), or putative ketoacyl-ACP reductase (CaKR1) were analyzed. The volatile ester contents of non-pungent C. chinense varieties (pamt/pamt) were equivalent to those of pungent varieties, but their levels were significantly lower in non-pungent NMCA30036 (pun12/pun12) and C. chinense (Cakr1/Cakr1) varieties. Multiple AAT-like sequences were identified from the pepper genome sequences, whereas only one CXE-like sequence was identified. Among these, AAT1, AAT2, and CXE1 were isolated from fruits of C. chinense and C. annuum. Gene expression analysis revealed that the AAT1 transcript level is a potential determinant of fruit volatile ester variations in Capsicum. Furthermore, enzymatic assays demonstrated that AAT1 is responsible for the biosynthesis of volatile esters in pepper fruit. Identification of a key gene for aroma biosynthesis in pepper fruit will provide a theoretical basis for the development of molecular tools for flavor improvement.

Genome assembly and analysis of Lactuca virosa: implications for lettuce breeding

Lettuce (Lactuca sativa L.) is a leafy vegetable crop with ongoing breeding efforts related to quality, resilience, and innovative production systems. To breed resilient and resistant lettuce in the future, valuable genetic variation found in close relatives could be further exploited. Lactuca virosa (2x = 2n = 18), a wild relative assigned to the tertiary lettuce gene pool, has a much larger genome (3.7 Gbp) than Lactuca sativa (2.5 Gbp). It has been used in interspecific crosses and is a donor to modern crisphead lettuce cultivars. Here, we present a de novo reference assembly of L. virosa with high continuity and complete gene space. This assembly facilitated comparisons to the genome of L. sativa and to that of the wild species L. saligna, a representative of the secondary lettuce gene pool. To assess the diversity in gene content, we classified the genes of the 3 Lactuca species as core, accessory, and unique. In addition, we identified 3 interspecific chromosomal inversions compared to L. sativa, which each may cause recombination suppression and thus hamper future introgression breeding. Using 3-way comparisons in both reference-based and reference-free manners, we show that the proliferation of long-terminal repeat elements has driven the genome expansion of L. virosa. Further, we performed a genome-wide comparison of immune genes, nucleotide-binding leucine-rich repeat, and receptor-like kinases among Lactuca spp. and indicated the evolutionary patterns and mechanisms behind their expansions. These genome analyses greatly facilitate the understanding of genetic variation in L. virosa, which is beneficial for the breeding of improved lettuce varieties.

Construction of a High-Density Paulownia Genetic Map and QTL Mapping of Important Phenotypic Traits Based on Genome Assembly and Whole-Genome Resequencing

Quantitative trait locus (QTL) mapping based on a genetic map is a very effective method of marker-assisted selection in breeding, and whole-genome resequencing is one of the useful methods to obtain high-density genetic maps. In this study, the hybrid assembly of Illumina, PacBio, and chromatin interaction mapping data was used to construct high-quality chromosomal genome sequences of Paulownia fortunei, with a size of 476.82 Mb, a heterozygosity of 0.52%, and a contig and scaffold N50s of 7.81 Mb and 21.81 Mb, respectively. Twenty scaffolds with a total length of 437.72 Mb were assembled into 20 pseudochromosomes. Repeat sequences with a total length of 243.96 Mb accounted for 51.16% of the entire genome. In all, 26,903 protein-coding gene loci were identified, and 26,008 (96.67%) genes had conserved functional motifs. Further comparative genomics analysis preliminarily showed that the split of P. fortunei with Tectona grandis likely occurred 38.8 (33.3-45.1) million years ago. Whole-genome resequencing was used to construct a merged genetic map of 20 linkage groups, with 2993 bin markers (3,312,780 SNPs), a total length of 1675.14 cm, and an average marker interval of 0.56 cm. In total, 73 QTLs for important phenotypic traits were identified (19 major QTLs with phenotypic variation explained ≥ 10%), including 10 for the diameter at breast height, 7 for the main trunk height, and 56 for branch-related traits. These results not only enrich P. fortunei genomic data but also form a solid foundation for fine QTL mapping and key marker/gene mining of Paulownia, which is of great significance for the directed genetic improvement of these species.

Exploring horticultural traits and disease resistance in Capsicum baccatum through segmental introgression lines

KEY MESSAGE

Segmental introgression and advanced backcross lines were developed and validated as important tools for improving agronomically important traits in pepper, offering improved sensitivity in detecting quantitative trait loci for breeding. Segmental introgression lines (SILs) and advanced backcross lines (ABs) can accelerate genetics and genomics research and breeding in crop plants. This study presents the development of a complete collection of SILs and ABs in pepper using Capsicum annuum cv. 'CM334' as the recipient parent and Capsicum baccatum 'PBC81', which displays various agronomically important traits including powdery mildew and anthracnose resistance, as donor parent. Using embryo rescue to overcome abortion in interspecific crosses, and marker-assisted selection with genotyping-in-thousands by sequencing (GT-seq) to develop SILs and ABs containing different segments of the C. baccatum genome, we obtained 63 SILs and 44 ABs, covering 94.8% of the C. baccatum genome. We characterized them for traits including powdery mildew resistance, anthracnose resistance, anthocyanin accumulation, trichome density, plant architecture, and fruit morphology. We validated previously known loci for these traits and discovered new sources of variation and quantitative trait loci (QTLs). A total of 15 QTLs were identified, including four for anthracnose resistance with three novel loci, seven for plant architecture, and four for fruit morphology. This is the first complete collection of pepper SILs and ABs validated for agronomic traits and will enhance QTL detection and serve as valuable breeding resources. Further, these SILs and ABs will be useful for comparative genomics and to better understand the genetic mechanisms underlying important agronomic traits in pepper, ultimately leading to improved crop productivity and sustainability.

Global co-expression network for key factor selection on environmental stress RNA-seq dataset in Capsicum annuum

Environmental stresses significantly affect plant growth, development, and productivity. Therefore, a deeper understanding of the underlying stress responses at the molecular level is needed. In this study, to identify critical genetic factors associated with environmental stress responses, the entire 737.3 Gb clean RNA-seq dataset across abiotic, biotic stress, and phytohormone conditions in Capsicum annuum was used to perform individual differentially expressed gene analysis and to construct gene co-expression networks for each stress condition. Subsequently, gene networks were reconstructed around transcription factors to identify critical factors involved in the stress responses, including the NLR gene family, previously implicated in resistance. The abiotic and biotic stress networks comprise 233 and 597 hubs respectively, with 10 and 89 NLRs. Each gene within the NLR groups in the network exhibited substantial expression to particular stresses. The integrated analysis strategy of the transcriptome network revealed potential key genes for complex environmental conditions. Together, this could provide important clues to uncover novel key factors using high-throughput transcriptome data in other species as well as plants.

Accurate, automated taxonomic assignment of genebank accessions: a new method demonstrated using high-throughput marker data from 10,000 Capsicum spp. accessions

KEY MESSAGE

We demonstrate how an algorithm that uses cheap genetic marker data can ensure the taxonomic assignments of genebank samples are complete, intuitive, and consistent-which enhances their value. To maximise the benefit of genebank resources, accurate and complete taxonomic assignments are imperative. The rise of genebank genomics allows genetic methods to be used to ensure this, but these need to be largely automated since the number of samples dealt with is too great for efficient manual recategorisation, however no clearly optimal method has yet arisen. A recent landmark genebank genomic study sequenced over 10,000 genebank accessions of peppers (Capsicum spp.), a species of great commercial, cultural, and scientific importance, which suffers from much taxonomic ambiguity. Similar datasets will, in coming decades, be produced for hundreds of plant taxa, affording a perfect opportunity to develop automated taxonomic correction methods in advance of the incipient genebank genomics explosion, alongside providing insights into pepper taxonomy in general. We present a marker-based taxonomic assignment approach that combines ideas from several standard classification algorithms, resulting in a highly flexible and customisable classifier suitable to impose intuitive assignments, even in highly reticulated species groups with complex population structures and evolutionary histories. Our classifier performs favourably compared with key alternative methods. Possible sensible alterations to pepper taxonomy based on the results are proposed for discussion by the relevant communities.

The Past Is Present: Coevolution of Viruses and Host Resistance Within Geographic Centers of Plant Diversity

Understanding the coevolutionary history of plants, pathogens, and disease resistance is vital for plant pathology. Here, I review Francis O. Holmes's work with tobacco mosaic virus (TMV) framed by the foundational work of Nikolai Vavilov on the geographic centers of origin of plants and crop wild relatives (CWRs) and T. Harper Goodspeed's taxonomy of the genus Nicotiana. Holmes developed a hypothesis that the origin of host resistance to viruses was due to coevolution of both at a geographic center. In the 1950s, Holmes proved that genetic resistance to TMV, especially dominant R-genes, was centered in South America for Nicotiana and other solanaceous plants, including Capsicum, potato, and tomato. One seeming exception was eggplant (Solanum melongena). Not until the acceptance of plate tectonics in the 1960s and recent advances in evolutionary taxonomy did it become evident that northeast Africa was the home of eggplant CWRs, far from Holmes's geographic center for TMV-R-gene coevolution. Unbeknownst to most plant pathologists, Holmes's ideas predated those of H.H. Flor, including experimental proof of the gene-for-gene interaction, identification of R-genes, and deployment of dominant host genes to protect crop plants from virus-associated yield losses.

A comprehensive review on genomic resources in medicinally and industrially important major spices for future breeding programs: Status, utility and challenges

In the global market, spices possess a high-value but low-volume commodities of commerce. The food industry depends largely on spices for taste, flavor, and therapeutic properties in replacement of cheap synthetic ones. The estimated growth rate for spices demand in the world is ∼3.19%. Since spices grow in limited geographical regions, India is one of the leading producer of spices, contributing 25-30 percent of total world trade. Hitherto, there has been no comprehensive review of the genomic resources of industrially important major medicinal spices to overcome major impediments in varietal improvement and management. This review focuses on currently available genomic resources of 24 commercially significant spices, namely, Ajwain, Allspice, Asafoetida, Black pepper, Cardamom large, Cardamom small, Celery, Chillies, Cinnamon, Clove, Coriander, Cumin, Curry leaf, Dill seed, Fennel, Fenugreek, Garlic, Ginger, Mint, Nutmeg, Saffron, Tamarind, Turmeric and Vanilla. The advent of low-cost sequencing machines has contributed immensely to the voluminous data generation of these spices, cracking the complex genomic architecture, marker discovery, and understanding comparative and functional genomics. This review of spice genomics resources concludes the perspective and way forward to provide footprints by uncovering genome assemblies, sequencing and re-sequencing projects, transcriptome-based studies, non-coding RNA-mediated regulation, organelles-based resources, developed molecular markers, web resources, databases and AI-directed resources in candidate spices for enhanced breeding potential in them. Further, their integration with molecular breeding could be of immense use in formulating a strategy to protect and expand the production of the spices due to increased global demand.

Extrachromosomal circular DNA and structural variants highlight genome instability in Arabidopsis epigenetic mutants

Abundant extrachromosomal circular DNA (eccDNA) is associated with transposable element (TE) activity. However, how the eccDNA compartment is controlled by epigenetic regulations and what is its impact on the genome is understudied. Here, using long reads, we sequence both the eccDNA compartment and the genome of Arabidopsis thaliana mutant plants affected in DNA methylation and post-transcriptional gene silencing. We detect a high load of TE-derived eccDNA with truncated and chimeric forms. On the genomic side, on top of truncated and full length TE neo-insertions, we detect complex structural variations (SVs) notably at a disease resistance cluster being a natural hotspot of SV. Finally, we serendipitously identify large tandem duplications in hypomethylated plants, suggesting that SVs could have been overlooked in epigenetic mutants. We propose that a high eccDNA load may alter DNA repair pathways leading to genome instability and the accumulation of SVs, at least in plants.

Tissues and mechanisms associated with Verticillium wilt resistance in tomato using bi-grafted near-isogenic lines

Host resistance is the primary means to control Verticillium dahliae, a soil-borne pathogen causing major losses on a broad range of plants, including tomato. The tissues and mechanisms responsible for resistance remain obscure. In the field, resistant tomato used as rootstocks does not confer resistance. Here, we created bi-grafted plants with near-isogenic lines (NILs) exhibiting (Ve1) or lacking (ve1) resistance to V. dahliae race 1. Ten days after inoculation, scion and rootstock tissues were subjected to differential gene expression and co-expression network analyses. Symptoms only developed in susceptible scions regardless of the rootstock. Infection caused more dramatic alteration of tomato gene expression in susceptible compared with resistant tissues, including pathogen receptor, signaling pathway, pathogenesis-related protein, and cell wall modification genes. Differences were observed between scions and rootstocks, primarily related to physiological processes in these tissues. Gene expression in scions was influenced by the rootstock genotype. A few genes were associated with the Ve1 genotype, which was independent of infection or tissue type. Several were physically clustered, some near the Ve1 locus on chromosome 9. Transcripts mapped to V. dahliae were dominated by secreted candidate effector proteins. These findings advance knowledge of molecular mechanisms underlying the tomato-V. dahliae interaction.

A physical map of traits of agronomic importance based on potato and tomato genome sequences

Potato, tomato, pepper, and eggplant are worldwide important crop and vegetable species of the Solanaceae family. Molecular linkage maps of these plants have been constructed and used to map qualitative and quantitative traits of agronomic importance. This research has been undertaken with the vision to identify the molecular basis of agronomic characters on the one hand, and on the other hand, to assist the selection of improved varieties in breeding programs by providing DNA-based markers that are diagnostic for specific agronomic characters. Since 2011, whole genome sequences of tomato and potato became available in public databases. They were used to combine the results of several hundred mapping and map-based cloning studies of phenotypic characters between 1988 and 2022 in physical maps of the twelve tomato and potato chromosomes. The traits evaluated were qualitative and quantitative resistance to pathogenic oomycetes, fungi, bacteria, viruses, nematodes, and insects. Furthermore, quantitative trait loci for yield and sugar content of tomato fruits and potato tubers and maturity or earliness were physically mapped. Cloned genes for pathogen resistance, a few genes underlying quantitative trait loci for yield, sugar content, and maturity, and several hundred candidate genes for these traits were included in the physical maps. The comparison between the physical chromosome maps revealed, in addition to known intrachromosomal inversions, several additional inversions and translocations between the otherwise highly collinear tomato and potato genomes. The integration of the positional information from independent mapping studies revealed the colocalization of qualitative and quantitative loci for resistance to different types of pathogens, called resistance hotspots, suggesting a similar molecular basis. Synteny between potato and tomato with respect to genomic positions of quantitative trait loci was frequently observed, indicating eventual similarity between the underlying genes.

Genome-Wide Identification of WRKY Gene Family and Functional Characterization of CcWRKY25 in Capsicum chinense

Pepper is renowned worldwide for its distinctive spicy flavor. While the gene expression characteristics of the capsaicinoid biosynthesis pathway have been extensively studied, there are already a few reports regarding transcriptional regulation in capsaicin biosynthesis. In this study, 73 WRKYs were identified in the genome of Capsicum chinense, and their physicochemical traits, DNA, and protein sequence characteristics were found to be complex. Combining RNA-seq and qRT-PCR data, the WRKY transcription factor CA06g13580, which was associated with the accumulation tendency of capsaicinoid, was screened and named CcWRKY25. CcWRKY25 was highly expressed in the placenta of spicy peppers. The heterologous expression of CcWRKY25 in Arabidopsis promoted the expression of genes PAL, 4CL1, 4CL2, 4CL3, CCR, and CCoAOMT and led to the accumulation of lignin and flavonoids. Furthermore, the expression of the capsaicinoid biosynthesis pathway genes (CBGs) pAMT, AT3, and KAS was significantly reduced in CcWRKY25-silenced pepper plants, resulting in a decrease in the amount of capsaicin. However, there was no noticeable difference in lignin accumulation. The findings suggested that CcWRKY25 could be involved in regulating capsaicinoid synthesis by promoting the expression of genes upstream of the phenylpropanoid pathway and inhibiting CBGs' expression. Moreover, the results highlighted the role of CcWRKY25 in controlling the pungency of pepper and suggested that the competitive relationship between lignin and capsaicin could also regulate the spiciness of the pepper.

Duplicate Genes Contribute to Variability in Abiotic Stress Resistance in Allopolyploid Wheat

Gene duplication is a universal biological phenomenon that drives genomic variation and diversity, plays a crucial role in plant evolution, and contributes to innovations in genetic engineering and crop development. Duplicated genes participate in the emergence of novel functionality, such as adaptability to new or more severe abiotic stress resistance. Future crop research will benefit from advanced, mechanistic understanding of the effects of gene duplication, especially in the development and deployment of high-performance, stress-resistant, elite wheat lines. In this review, we summarize the current knowledge of gene duplication in wheat, including the principle of gene duplication and its effects on gene function, the diversity of duplicated genes, and how they have functionally diverged. Then, we discuss how duplicated genes contribute to abiotic stress response and the mechanisms of duplication. Finally, we have a future prospects section that discusses the direction of future efforts in the short term regarding the elucidation of replication and retention mechanisms of repetitive genes related to abiotic stress response in wheat, excellent gene function research, and practical applications.

Investigation of Resistance Genes in Genus Vigna Reveals Highly Variable NLRome in Parallel Domesticated Member Species

Vigna is a unique genus that consist of multiple crop species that are domesticated in parallel fashion between 7-10 thousand years ago. Here we studied the evolution of nucleotide-binding site leucine-rich repeat receptor (NLR) genes across five crop species of genus Vigna. In total identified 286, 350, 234, 250, 108 and 161 NLR genes were from Phaseolous vulgaris, Vigna. unguiculata, Vigna mungo, Vigna radiata, Vigna angularis and Vigna umbellata respectively. Comprehensive phylogenetic and clusterization analysis reveals the presence of seven subgroups of Coiled coil like NLRs (CC-NLR) genes and four distinct lineages of Toll interleukin receptor like NLRs (TIR-NLR). Subgroup CC_G10-NLR shows large scale diversification among Vigna species suggesting genus specific distinct duplication pattern in Vigna species. Mainly birth of new NLR gene families and higher rate of terminal duplication is the major determinants for expansion of NLRome in genus Vigna. Recent expansion of NLRome in V. anguiculata and V. radiata was also observed which might suggest that domestication have supported their duplication of lineage specific NLR genes. In short, large scale difference in the architecture of NLRome were observed in diploid plant species. Our findings allowed us to hypothesized that independent parallel domestication is the major drivers of highly divergent evolution of NLRome in genus Vigna.

Mapping of the bs5 and bs6 non-race-specific recessive resistances against bacterial spot of pepper

Bacterial spot caused by Xanthomonas euvesicatoria is a major disease of pepper (Capsicum annuum L.) in warm and humid production environments. Use of genetically resistant cultivars is an effective approach to manage bacterial spot. Two recessive resistance genes, bs5 and bs6, confer non-race-specific resistance against bacterial spot. The objective of our study was to map these two loci in the pepper genome. We used a genotyping-by-sequencing approach to initially map the position of the two resistances. Segregating populations for bs5 and bs6 were developed by crossing susceptible Early CalWonder (ECW) with near-isogenic lines ECW50R (bs5 introgression) or ECW60R (bs6 introgression). Following fine-mapping, bs5 was delimited to a ~535 Kbp interval on chromosome 3, and bs6 to a ~666 Kbp interval in chromosome 6. We identified 14 and 8 candidate resistance genes for bs5 and bs6, respectively, based on predicted protein coding polymorphisms between ECW and the corresponding resistant parent. This research enhances marker-assisted selection of bs5 and bs6 in breeding programs and is a crucial step towards elucidating the molecular mechanisms underlying the resistances.

Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding

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

Nuclear phylogeny and insights into whole-genome duplications and reproductive development of Solanaceae plants

Solanaceae, the nightshade family, have ∼2700 species, including the important crops potato and tomato, ornamentals, and medicinal plants. Several sequenced Solanaceae genomes show evidence for whole-genome duplication (WGD), providing an excellent opportunity to investigate WGD and its impacts. Here, we generated 93 transcriptomes/genomes and combined them with 87 public datasets, for a total of 180 Solanaceae species representing all four subfamilies and 14 of 15 tribes. Nearly 1700 nuclear genes from these transcriptomic/genomic datasets were used to reconstruct a highly resolved Solanaceae phylogenetic tree with six major clades. The Solanaceae tree supports four previously recognized subfamilies (Goetzeioideae, Cestroideae, Nicotianoideae, and Solanoideae) and the designation of three other subfamilies (Schizanthoideae, Schwenckioideae, and Petunioideae), with the placement of several previously unassigned genera. We placed a Solanaceae-specific whole-genome triplication (WGT1) at ∼81 million years ago (mya), before the divergence of Schizanthoideae from other Solanaceae subfamilies at ∼73 mya. In addition, we detected two gene duplication bursts (GDBs) supporting proposed WGD events and four other GDBs. An investigation of the evolutionary histories of homologs of carpel and fruit developmental genes in 14 gene (sub)families revealed that 21 gene clades have retained gene duplicates. These were likely generated by the Solanaceae WGT1 and may have promoted fleshy fruit development. This study presents a well-resolved Solanaceae phylogeny and a new perspective on retained gene duplicates and carpel/fruit development, providing an improved understanding of Solanaceae evolution.

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.

Revealing evolution of tropane alkaloid biosynthesis by analyzing two genomes in the Solanaceae family

Tropane alkaloids (TAs) are widely distributed in the Solanaceae, while some important medicinal tropane alkaloids (mTAs), such as hyoscyamine and scopolamine, are restricted to certain species/tribes in this family. Little is known about the genomic basis and evolution of TAs biosynthesis and specialization in the Solanaceae. Here, we present chromosome-level genomes of two representative mTAs-producing species: Atropa belladonna and Datura stramonium. Our results reveal that the two species employ a conserved biosynthetic pathway to produce mTAs despite being distantly related within the nightshade family. A conserved gene cluster combined with gene duplication underlies the wide distribution of TAs in this family. We also provide evidence that branching genes leading to mTAs likely have evolved in early ancestral Solanaceae species but have been lost in most of the lineages, with A. belladonna and D. stramonium being exceptions. Furthermore, we identify a cytochrome P450 that modifies hyoscyamine into norhyoscyamine. Our results provide a genomic basis for evolutionary insights into the biosynthesis of TAs in the Solanaceae and will be useful for biotechnological production of mTAs via synthetic biology approaches.

Genetic Analyses of Seed Longevity in Capsicum annuum L. in Cold Storage Conditions

Seed longevity is the most important trait in the genebank management system. No seed can remain infinitely viable. There are 1241 accessions of Capsicum annuum L. available at the German Federal ex situ genebank at IPK Gatersleben. C. annuum (Capsicum) is the most economically important species of the genus Capsicum. So far, there is no report that has addressed the genetic basis of seed longevity in Capsicum. Here, we convened a total of 1152 Capsicum accessions that were deposited in Gatersleben over forty years (from 1976 to 2017) and assessed their longevity by analyzing the standard germination percentage after 5-40 years of storage at -15/-18 °C. These data were used to determine the genetic causes of seed longevity, along with 23,462 single nucleotide polymorphism (SNP) markers covering all of the 12 Capsicum chromosomes. Using the association-mapping approach, we identified a total of 224 marker trait associations (MTAs) (34, 25, 31, 35, 39, 7, 21 and 32 MTAs after 5-, 10-, 15-, 20-, 25-, 30-, 35- and 40-year storage intervals) on all the Capsicum chromosomes. Several candidate genes were identified using the blast analysis of SNPs, and these candidate genes are discussed.

Diversity and expression analysis of ZIP transporters and associated metabolites under zinc and iron stress in Capsicum

The members of ZRT, IRT-like protein (ZIP) family are involved in the uptake and transportation of several metal ions. Here, we report a comprehensive identification of ZIP transporter genes from Capsicum annuum, C. chinense, and C. baccatum, and their expression analysis under Zn and Fe stress. Changes in root morphology and differential accumulation of several metabolites from sugars, amino acids, carboxylic acids, and fatty acids in root and leaf tissues of plants in the absence of Zn and Fe were observed. Further, metabolites such as L-aspartic acid, 2-ketoglutaric acids, β-L-fucopyranose, quininic acid, chlorogenic acid, and aucubin were significantly upregulated in root and leaf tissues under Zn/Fe deprived conditions. qRT-PCR analysis of 17 CaZIPs in different tissues revealed tissue-specific expression of CaZIP1-2, CaZIP4-8, CaZIP13, and CaZIP16-17 under normal conditions. However, the absence of Zn and Fe significantly induced the expression of CaZIP4-5, CaZIP7-9, and CaZIP14 genes in root and leaf tissues. Additionally, in the absence of Fe, upregulation of CaZIP4-5 and CaZIP8 and increased uptake of mineral elements Cu, Zn, Mg, P, and S were observed in roots, suggesting their potential role in metal-ion uptake in Capsicum. The identified genes provide the basis for future studies of mineral uptake and their biofortification to increase the nutritional values in Capsicum.

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.

A 163-bp insertion in the Capana10g000198 encoding a MYB transcription factor causes male sterility in pepper (Capsicum annuum L.)

Male sterility provides an efficient approach for commercial exploitation of heterosis. Despite more than 20 genic male sterile (GMS) mutants documented in pepper (Capsicum annuum L.), only two causal genes have been successfully identified. Here, a novel spontaneous recessive GMS mutant, designated msc-3, is identified and characterized at both phenotypic and histological levels. Pollen abortion of msc-3 mutant may be due to the delayed tapetum degradation, leading to the non-degeneration of tetrads callosic wall. Then, a modified MutMap method and molecular marker linkage analysis were employed to fine mapping the msc-3 locus, which was delimited to the ~139.91-kb region harboring 10 annotated genes. Gene expression and structure variation analyses indicate the Capana10g000198, encoding a R2R3-MYB transcription factor, is the best candidate gene for the msc-3 locus. Expression profiling analysis shows the Capana10g000198 is an anther-specific gene, and a 163-bp insertion in the Capana10g000198 is highly correlated with the male sterile (MS) phenotype. Additionally, downregulation of Capana10g000198 in male fertile plants through virus-induced gene silencing resulted in male sterility. Finally, possible regulatory relationships of the msc-3 gene with the other two reported pepper GMS genes, msc-1 and msc-2, have been studied, and comparative transcriptome analysis reveals the expression of 16 GMS homologs are significantly downregulated in the MS anthers. Overall, our results reveal that Capana10g000198 is the causal gene underlying the msc-3 locus, providing important theoretical clues and basis for further in-depth study on the regulatory mechanisms of pollen development in pepper.

Comparative and expression analyses of AP2/ERF genes reveal copy number expansion and potential functions of ERF genes in Solanaceae

BACKGROUND

The AP2/ERF gene family is a superfamily of transcription factors that are important in the response of plants to abiotic stress and development. However, comprehensive research of the AP2/ERF genes in the Solanaceae family is lacking.

RESULTS

Here, we updated the annotation of AP2/ERF genes in the genomes of eight Solanaceae species, as well as Arabidopsis thaliana and Oryza sativa. We identified 2,195 AP2/ERF genes, of which 368 (17%) were newly identified. Based on phylogenetic analyses, we observed expansion of the copy number of these genes, especially those belonging to specific Ethylene-Responsive Factor (ERF) subgroups of the Solanaceae. From the results of chromosomal location and synteny analyses, we identified that the AP2/ERF genes of the pepper (Capsicum annuum), the tomato (Solanum lycopersicum), and the potato (Solanum tuberosum) belonging to ERF subgroups form a tandem array and most of them are species-specific without orthologs in other species, which has led to differentiation of AP2/ERF gene repertory among Solanaceae. We suggest that these genes mainly emerged through recent gene duplication after the divergence of these species. Transcriptome analyses showed that the genes have a putative function in the response of the pepper and tomato to abiotic stress, especially those in ERF subgroups.

CONCLUSIONS

Our findings will provide comprehensive information on AP2/ERF genes and insights into the structural, evolutionary, and functional understanding of the role of these genes in the Solanaceae.

Chromosome-scale genome assembly of a Japanese chili pepper landrace, Capsicum annuum 'Takanotsume'

Here, we report the genome sequence of a popular Japanese chili pepper landrace, Capsicum annuum 'Takanotsume'. We used long-read sequencing and optical mapping, together with the genetic mapping technique, to obtain the chromosome-scale genome assembly of 'Takanotsume'. The assembly consists of 12 pseudomolecules, which corresponds to the basic chromosome number of C. annuum, and is 3,058.5 Mb in size, spanning 97.0% of the estimated genome size. A total of 34,324 high-confidence genes were predicted in the genome, and 83.4% of the genome assembly was occupied by repetitive sequences. Comparative genomics of linked-read sequencing-derived de novo genome assemblies of two Capsicum chinense lines and whole-genome resequencing analysis of Capsicum species revealed not only nucleotide sequence variations but also genome structure variations (i.e. chromosomal rearrangements and transposon-insertion polymorphisms) between 'Takanotsume' and its relatives. Overall, the genome sequence data generated in this study will accelerate the pan-genomics and breeding of Capsicum, and facilitate the dissection of genetic mechanisms underlying the agronomically important traits of 'Takanotsume'.

GlPS1 overexpression accumulates coumarin secondary metabolites in transgenic Arabidopsis

The dried root of Glehnia littoralis is a traditional Chinese herbal medicine mainly used to treat lung diseases and plays an important role in fighting coronavirus disease 2019 pneumonia in China. This study focused on the key enzyme gene GlPS1 for furanocoumarin synthesis in G. littoralis. In the 35S:GlPS1 transgenic Arabidopsis study, the Arabidopsis thaliana-overexpressing GlPS1 gene was more salt-tolerant than Arabidopsis in the blank group. Metabolomics analysis showed 30 differential metabolites in Arabidopsis, which overexpressed the GlPS1 gene. Twelve coumarin compounds were significantly upregulated, and six of these coumarin compounds were not detected in the blank group. Among these differential coumarin metabolites, isopimpinellin and aesculetin have been annotated by the Kyoto Encyclopedia of Genes and Genomes and isopimpinellin was not detected in the blank group. Through structural comparison, imperatorin was formed by dehydration and condensation of zanthotoxol and a molecule of isoprenol, and the difference between them was only one isoprene. Results showed that the GlPS1 gene positively regulated the synthesis of coumarin metabolites in A. thaliana and at the same time improved the salt tolerance of A. thaliana.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11240-022-02427-w.

Strive or thrive: Trends in Phytophthora capsici gene expression in partially resistant pepper

Partial resistance in plants generally exerts a low selective pressure on pathogens, and thus ensuring their durability in agrosystems. However, little is known about the effect of partial resistance on the molecular mechanisms of pathogenicity, a knowledge that could advance plant breeding for sustainable plant health. Here we investigate the gene expression of Phytophthora capsici during infection of pepper (Capsicum annuum L.), where only partial genetic resistance is reported, using Illumina RNA-seq. Comparison of transcriptomes of P. capsici infecting susceptible and partially resistant peppers identified a small number of genes that redirected its own resources into lipid biosynthesis to subsist on partially resistant plants. The adapted and non-adapted isolates of P. capsici differed in expression of genes involved in nucleic acid synthesis and transporters. Transient ectopic expression of the RxLR effector genes CUST_2407 and CUST_16519 in pepper lines differing in resistance levels revealed specific host-isolate interactions that either triggered local necrotic lesions (hypersensitive response or HR) or elicited leave abscission (extreme resistance or ER), preventing the spread of the pathogen to healthy tissue. Although these effectors did not unequivocally explain the quantitative host resistance, our findings highlight the importance of plant genes limiting nutrient resources to select pepper cultivars with sustainable resistance to P. capsici.