Whole genome resequencing of four Italian sweet pepper landraces provides insights on sequence variation in genes of agronomic value

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.

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.

A Compendium for Novel Marker-Based Breeding Strategies in Eggplant

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

First genetic maps development and QTL mining in Ranunculus asiaticus L. through ddRADseq

Persian Buttercup (Ranunculus asiaticus L.; 2x=2n=16; estimated genome size: 7.6Gb) is an ornamental and perennial crop native of Asia Minor and Mediterranean basin, marketed both as cut flower or potted plant. Currently new varieties are developed by selecting plants carrying desirable traits in segregating progenies obtained by controlled mating, which are propagated through rhizomes or micro-propagated in vitro. In order to escalate selection efficiency and respond to market requests, more knowledge of buttercup genetics would facilitate the identification of markers associated with loci and genes controlling key ornamental traits, opening the way for molecular assisted breeding programs. Reduced-representation sequencing (RRS) represents a powerful tool for plant genotyping, especially in case of large genomes such as the one of buttercup, and have been applied for the development of high-density genetic maps in several species. We report on the development of the first molecular-genetic maps in R. asiaticus based on of a two-way pseudo-testcross strategy. A double digest restriction-site associated DNA (ddRAD) approach was applied for genotyping two F₁ mapping populations, whose female parents were a genotype of a so called 'ponpon' and of a 'double flower' varieties, while the common male parental ('Cipro') was a genotype producing a simple flower. The ddRAD generated a total of ~2Gb demultiplexed reads, resulting in an average of 8,3M reads per line. The sstacks pipeline was applied for the construction of a mock reference genome based on sequencing data, and SNP markers segregating in only one of the parents were retained for map construction by treating the F₁ population as a backcross. The four parental maps (two of the female parents and two of the common male parent) were aligned with 106 common markers and 8 linkage groups were identified, corresponding to the haploid chromosome number of the species. An average of 586 markers were associated with each parental map, with a marker density ranging from 1 marker/cM to 4.4 markers/cM. The developed maps were used for QTL analysis for flower color, leading to the identification of major QTLs for purple pigmentation. These results contribute to dissect on the genetics of Persian buttercup, enabling the development of new approaches for future varietal development.

High-quality chromosome-scale genomes facilitate effective identification of large structural variations in hot and sweet peppers

Pepper (Capsicum annuum) is an important vegetable crop that has been subjected to intensive breeding, resulting in limited genetic diversity, especially for sweet peppers. Previous studies have reported pepper draft genome assemblies using short read sequencing, but their capture of the extent of large structural variants (SVs), such as presence-absence variants (PAVs), inversions, and copy-number variants (CNVs) in the complex pepper genome falls short. In this study, we sequenced the genomes of representative sweet and hot pepper accessions by long-read and/or linked-read methods and advanced scaffolding technologies. First, we developed a high-quality reference genome for the sweet pepper cultivar 'Dempsey' and then used the reference genome to identify SVs in 11 other pepper accessions and constructed a graph-based pan-genome for pepper. We annotated an average of 42 972 gene families in each pepper accession, defining a set of 19 662 core and 23 115 non-core gene families. The new pepper pan-genome includes informative variants, 222 159 PAVs, 12 322 CNVs, and 16 032 inversions. Pan-genome analysis revealed PAVs associated with important agricultural traits, including potyvirus resistance, fruit color, pungency, and pepper fruit orientation. Comparatively, a large number of genes are affected by PAVs, which is positively correlated with the high frequency of transposable elements (TEs), indicating TEs play a key role in shaping the genomic landscape of peppers. The datasets presented herein provide a powerful new genomic resource for genetic analysis and genome-assisted breeding for pepper improvement.

Monograph of wild and cultivated chili peppers (Capsicum L., Solanaceae)

Capsicum L. (tribe Capsiceae, Solanaceae) is an American genus distributed ranging from the southern United States of America to central Argentina and Brazil. The genus includes chili peppers, bell peppers, ajíes, habaneros, jalapeños, ulupicas and pimientos, well known for their economic importance around the globe. Within the Solanaceae, the genus can be recognised by its shrubby habit, actinomorphic flowers, distinctive truncate calyx with or without appendages, anthers opening by longitudinal slits, nectaries at the base of the ovary and the variously coloured and usually pungent fruits. The highest diversity of this genus is located along the northern and central Andes. Although Capsicum has been extensively studied and great advances have been made in the understanding of its taxonomy and the relationships amongst species, there is no monographic treatment of the genus as a whole. Based on morphological and molecular evidence studied from field and herbarium specimens, we present here a comprehensive taxonomic treatment for the genus, including updated information about morphology, anatomy, karyology, phylogeny and distribution. We recognise 43 species and five varieties, including C.mirum Barboza, sp. nov. from São Paulo State, Brazil and a new combination C.muticum (Sendtn.) Barboza, comb. nov.; five of these taxa are cultivated worldwide (C.annuumL.var.annuum, C.baccatumL.var.pendulum (Willd.) Eshbaugh, C.baccatumL.var.umbilicatum (Vell.) Hunz. & Barboza, C.chinense Jacq. and C.frutescens L.). Nomenclatural revision of the 265 names attributed to chili peppers resulted in 89 new lectotypifications and five new neotypifications. Identification keys and detailed descriptions, maps and illustrations for all taxa are provided.

Whole-genome resequencing reveals genomic footprints of Italian sweet and hot pepper heirlooms giving insight into genes underlying key agronomic and qualitative traits

BACKGROUND

Pepper is a major crop species of the Solanaceae family, largely appreciated for its high nutritional and healthy contribution to human diets. In the Mediterranean basin, the favorable pedoclimatic conditions enhanced the selection of several diversified landraces cultivated pepper (Capsicum annuum), for whom Italy can be considered a main pole of diversification. Hence, a survey of traditional C. annuum genetic resources is essential for deep understanding of such diversity and for applications in genomics assisted breeding. Here, we report whole-genome resequencing analyses of two sweet and two pungent genotypes highly diffused in South Italy and representative of the variability for shape, colour and nutritional properties.

RESULTS

The four genomes were reconstructed at a chromosomal scale using a reference-guided approach, based on a dataset of 2.6 billion paired-end reads, corresponding to 20× genome coverage and a mapping rate above 99% for a final genomes size of approximately 3 Gb. After five iterations of variant calling, a total of 29,258,818 single nucleotide polymorphisms (SNPs) and 1,879,112 InDels, were identified. Substantial differences were observed among the four genomes based on geographical origin, with chromosomes 9 and 11 showing more polymorphisms in the accessions with higher fruit weight and absence of pungency. Among the identified variants, a small private indel (T - > TA) shared between sweet and big fruits accessions induces a frameshift with the generation of a new stop codon in a gene annotated as extensin, whereas two private SNPs within hot types were identified in 1-aminocyclopropane-1-carboxylate oxidase (ACO), a key gene involved in fruit ripening. The estimation of repetitive elements highlights a preponderant presence of Long Terminal Repeats (LTRs), the majority of which belonged to Gypsy superfamily. By comparing the four genomes with publicly available references including 'CM334' and Zunla-1 highlight the presence of 49,475 shared gene families.

CONCLUSIONS

The new genomic sequences aim to enrich the whole genome information of pepper local varieties, providing a valuable tool for precision gene mapping, marker discovery, comparative studies. Such knowledge widens the frontiers to understand the selection history of Italian pepper landraces toward the recognition of specificity local agri-food products marks.

Genome-Wide Survey and Development of the First Microsatellite Markers Database (AnCorDB) in Anemone coronaria L

Anemone coronaria L. (2n = 2x = 16) is a perennial, allogamous, highly heterozygous plant marketed as a cut flower or in gardens. Due to its large genome size, limited efforts have been made in order to develop species-specific molecular markers. We obtained the first draft genome of the species by Illumina sequencing an androgenetic haploid plant of the commercial line “MISTRAL^® Magenta”. The genome assembly was obtained by applying the MEGAHIT pipeline and consisted of 2 × 10⁶ scaffolds. The SciRoKo SSR (Simple Sequence Repeats)-search module identified 401.822 perfect and 188.987 imperfect microsatellites motifs. Following, we developed a user-friendly “Anemone coronaria Microsatellite DataBase” (AnCorDB), which incorporates the Primer3 script, making it possible to design couples of primers for downstream application of the identified SSR markers. Eight genotypes belonging to eight cultivars were used to validate 62 SSRs and a subset of markers was applied for fingerprinting each cultivar, as well as to assess their intra-cultivar variability. The newly developed microsatellite markers will find application in Breeding Rights disputes, developing genetic maps, marker assisted breeding (MAS) strategies, as well as phylogenetic studies.

Transcriptome-Based Identification and Functional Characterization of NAC Transcription Factors Responsive to Drought Stress in Capsicum annuum L

Capsicum annuum L. is one of the most cultivated Solanaceae species, and in the open field, water limitation leading to drought stress affects its fruit quality, fruit setting, fruit size and ultimately yield. We identified stage-specific and a common core set of differentially expressed genes, following RNA-seq transcriptome analyses of a breeding line subjected to acute drought stress followed by recovery (rewatering), at three stages of plant development. Among them, two NAC transcription factor (TF) genes, i.e., CaNAC072 and CaNAC104, were always upregulated after drought stress and downregulated after recovery. The two TF proteins were observed to be localized in the nucleus following their transient expression in Nicotiana benthamiana leaves. The expression of the two NACs was also induced by NaCl, polyethylene glycol (PEG) and abscisic acid (ABA) treatments, suggesting that CaNAC072 is an early, while CaNAC104 is a late abiotic stress-responsive gene. Virus-induced gene silencing (VIGS) of CaNAC104 did not affect the pepper plantlet’s tolerance to drought stress, while VIGS of CaNAC072 increased drought tolerance. Heterologous expression of CaNAC072 in Arabidopsis thaliana as well as in plants mutated for its homolog ANAC072 did not increase drought stress tolerance. This highlights a different role of the two NAC homologs in the two species. Here, we discuss the complex role of NACs as transcriptional switches in the response to drought stress in bell pepper.

Identification of Novel miRNAs and Their Target Genes in the Response to Abscisic Acid in Arabidopsis

miRNAs are involved in various biological processes, including adaptive responses to abiotic stress. To understand the role of miRNAs in the response to ABA, ABA-responsive miRNAs were identified by small RNA sequencing in wild-type Arabidopsis, as well as in abi1td, mkkk17, and mkkk18 mutants. We identified 10 novel miRNAs in WT after ABA treatment, while in abi1td, mkkk17, and mkkk18 mutants, three, seven, and nine known miRNAs, respectively, were differentially expressed after ABA treatment. One novel miRNA (miRn-8) was differentially expressed in the mkkk17 mutant. Potential target genes of the miRNA panel were identified using psRNATarget. Sequencing results were validated by quantitative RT-PCR of several known and novel miRNAs in all genotypes. Of the predicted targets of novel miRNAs, seven target genes of six novel miRNAs were further validated by 5' RLM-RACE. Gene ontology analyses showed the potential target genes of ABA-responsive known and novel miRNAs to be involved in diverse cellular processes in plants, including development and stomatal movement. These outcomes suggest that a number of the identified miRNAs have crucial roles in plant responses to environmental stress, as well as in plant development, and might have common regulatory roles in the core ABA signaling pathway.

Improved genome assembly and pan-genome provide key insights into eggplant domestication and breeding

Eggplant (Solanum melongena L.) is an important horticultural crop and one of the most widely grown vegetables from the Solanaceae family. It was domesticated from a wild, prickly progenitor carrying small, round, non-anthocyanic fruits. We obtained a novel, highly contiguous genome assembly of the eggplant '67/3' reference line, by Hi-C retrofitting of a previously released short read- and optical mapping-based assembly. The sizes of the 12 chromosomes and the fraction of anchored genes in the improved assembly were comparable to those of a chromosome-level assembly. We resequenced 23 accessions of S. melongena representative of the worldwide phenotypic, geographic, and genetic diversity of the species, and one each from the closely related species Solanum insanum and Solanum incanum. The eggplant pan-genome contained approximately 51.5 additional megabases and 816 additional genes compared with the reference genome, while the pan-plastome showed little genetic variation. We identified 53 selective sweeps related to fruit color, prickliness, and fruit shape in the nuclear genome, highlighting selection leading to the emergence of present-day S. melongena cultivars from its wild ancestors. Candidate genes underlying the selective sweeps included a MYBL1 repressor and CHALCONE ISOMERASE (for fruit color), homologs of Arabidopsis GLABRA1 and GLABROUS INFLORESCENCE STEMS2 (for prickliness), and orthologs of tomato FW2.2, OVATE, LOCULE NUMBER/WUSCHEL, SUPPRESSOR OF OVATE, and CELL SIZE REGULATOR (for fruit size/shape), further suggesting that selection for the latter trait relied on a common set of orthologous genes in tomato and eggplant.

Whole-genome assembly of Corylus avellana cv'Tonda Gentile delle Langhe' using linked-reads (10X Genomics)

The European hazelnut (Corylus avellana L.; 2n = 2x = 22) is a worldwide economically important tree nut that is cross-pollinated due to sporophytic incompatibility. Therefore, any individual plant is highly heterozygous. Cultivars are clonally propagated using mound layering, rooted suckers, and micropropagation. In recent years, the interest in this crop has increased, due to a growing demand related to the recognized health benefits of nut consumption. C. avellana cv “Tonda Gentile delle Langhe” (“TGdL”) is well-known for its high kernel quality, and the premium price paid for this cultivar is an economic benefit for producers in northern Italy. Assembly of a high-quality genome is a difficult task in many plant species because of the high level of heterozygosity. We assembled a chromosome-level genome sequence of “TGdL” with a two-step approach. First, 10X Genomics Chromium Technology was used to create a high-quality sequence, which was then assembled into scaffolds with cv “Tombul” genome as the reference. Eleven pseudomolecules were obtained, corresponding to 11 chromosomes. A total of 11,046 scaffolds remained unplaced, representing 11% of the genome (46,504,161 bp). Gene prediction, performed with Maker-P software, identified 27,791 genes (AED ≤0.4 and 92% of BUSCO completeness), whose function was analyzed with BlastP and InterProScan software. To characterize “TGdL” specific genetic mechanisms, Orthofinder was used to detect orthologs between hazelnut and closely related species. The “TGdL” genome sequence is expected to be a powerful tool to understand hazelnut genetics and allow detection of markers/genes for important traits to be used in targeted breeding programs.

Cultivated Tomato (Solanum lycopersicum L.) Suffered a Severe Cytoplasmic Bottleneck during Domestication: Implications from Chloroplast Genomes

In various crops, genetic bottlenecks occurring through domestication can limit crop resilience to biotic and abiotic stresses. In the present study, we investigated nucleotide diversity in tomato chloroplast genome through sequencing seven plastomes of cultivated accessions from the Campania region (Southern Italy) and two wild species among the closest (Solanum pimpinellifolium) and most distantly related (S. neorickii) species to cultivated tomatoes. Comparative analyses among the chloroplast genomes sequenced in this work and those available in GenBank allowed evaluating the variability of plastomes and defining phylogenetic relationships. A dramatic reduction in genetic diversity was detected in cultivated tomatoes, nonetheless, a few de novo mutations, which still differentiated the cultivated tomatoes from the closest wild relative S. pimpinellifolium, were detected and are potentially utilizable as diagnostic markers. Phylogenetic analyses confirmed that S. pimpinellifolium is the closest ancestor of all cultivated tomatoes. Local accessions all clustered together and were strictly related with other cultivated tomatoes (S. lycopersicum group). Noteworthy, S. lycopersicum var. cerasiforme resulted in a mixture of both cultivated and wild tomato genotypes since one of the two analyzed accessions clustered with cultivated tomato, whereas the other with S. pimpinellifolium. Overall, our results revealed a very reduced cytoplasmic variability in cultivated tomatoes and suggest the occurrence of a cytoplasmic bottleneck during their domestication.