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.

Boosting grapevine breeding for climate-smart viticulture: from genetic resources to predictive genomics

The multifaceted nature of climate change is increasing the urgency to select resilient grapevine varieties, or generate new, fitter cultivars, to withstand a multitude of new challenging conditions. The attainment of this goal is hindered by the limiting pace of traditional breeding approaches, which require decades to result in new selections. On the other hand, marker-assisted breeding has proved useful when it comes to traits governed by one or few genes with great effects on the phenotype, but its efficacy is still restricted for complex traits controlled by many loci. On these premises, innovative strategies are emerging which could help guide selection, taking advantage of the genetic diversity within the Vitis genus in its entirety. Multiple germplasm collections are also available as a source of genetic material for the introgression of alleles of interest via adapted and pioneering transformation protocols, which present themselves as promising tools for future applications on a notably recalcitrant species such as grapevine. Genome editing intersects both these strategies, not only by being an alternative to obtain focused changes in a relatively rapid way, but also by supporting a fine-tuning of new genotypes developed with other methods. A review on the state of the art concerning the available genetic resources and the possibilities of use of innovative techniques in aid of selection is presented here to support the production of climate-smart grapevine genotypes.

Structural and expression analysis of polyphenol oxidases potentially involved in globe artichoke (C. cardunculus var. scolymus L.) tissue browning

Globe artichoke capitula are susceptible to browning due to oxidation of phenols caused by the activity of polyphenol oxidases (PPOs), this reduces their suitability for fresh or processed uses. A genome-wide analysis of the globe artichoke PPO gene family was performed. Bioinformatics analyses identified eleven PPOs and their genomic and amino acidic features were annotated. Cis-acting element analysis identified a gene regulatory and functional profile associated to plant growth and development as well as stress response. For some PPOs, phylogenetic analyses revealed a structural and functional conservation with different Asteraceae PPOs, while the allelic variants of the eleven PPOs investigated across four globe artichoke varietal types identified several SNP/Indel variants, some of which having impact on gene translation. By RTqPCR were assessed the expression patterns of PPOs in plant tissues and in vitro calli characterized by different morphologies. Heterogeneous PPO expression profiles were observed and three of them (PPO6, 7 and 11) showed a significant increase of transcripts in capitula tissues after cutting. Analogously, the same three PPOs were significantly up-regulated in calli showing a brown phenotype due to oxidation of phenols. Our results lay the foundations for a future application of gene editing aimed at disabling the three PPOs putatively involved in capitula browning.

Genomic Analysis Highlights Putative Defective Susceptibility Genes in Tomato Germplasm

Tomato (Solanum lycopersicum L.) is one of the most widely grown vegetables in the world and is impacted by many diseases which cause yield reduction or even crop failure. Breeding for disease resistance is thus a key objective in tomato improvement. Since disease arises from a compatible interaction between a plant and a pathogen, a mutation which alters a plant susceptibility (S) gene facilitating compatibility may induce broad-spectrum and durable plant resistance. Here, we report on a genome-wide analysis of a set of 360 tomato genotypes, with the goal of identifying defective S-gene alleles as a potential source for the breeding of resistance. A set of 125 gene homologs of 10 S-genes (PMR 4, PMR5, PMR6, MLO, BIK1, DMR1, DMR6, DND1, CPR5, and SR1) were analyzed. Their genomic sequences were examined and SNPs/indels were annotated using the SNPeff pipeline. A total of 54,000 SNPs/indels were identified, among which 1300 were estimated to have a moderate impact (non-synonymous variants), while 120 were estimated to have a high impact (e.g., missense/nonsense/frameshift variants). The latter were then analyzed for their effect on gene functionality. A total of 103 genotypes showed one high-impact mutation in at least one of the scouted genes, while in 10 genotypes, more than 4 high-impact mutations in as many genes were detected. A set of 10 SNPs were validated through Sanger sequencing. Three genotypes carrying high-impact homozygous SNPs in S-genes were infected with Oidium neolycopersici, and two highlighted a significantly reduced susceptibility to the fungus. The existing mutations fall within the scope of a history of safe use and can be useful to guide risk assessment in evaluating the effect of new genomic techniques.

Diversity analyses in two ornamental and large-genome Ranunculaceae species based on a low-cost Klenow NGS-based protocol

Persian buttercup (Ranunculus asiaticus L.) and poppy anemone (Anemone coronaria L.) are ornamental, outcrossing, perennial species belonging to the Ranunculaceae family, characterized by large and highly repetitive genomes. We applied K-seq protocol in both species to generate high-throughput sequencing data and produce a large number of genetic polymorphisms. The technique entails the application of Klenow polymerase-based PCR using short primers designed by analyzing k-mer sets in the genome sequence. To date the genome sequence of both species has not been released, thus we designed primer sets based on the reference the genome sequence of the related species Aquilegia oxysepala var. kansuensis (Brühl). A whole of 11,542 SNPs were selected for assessing genetic diversity of eighteen commercial varieties of R. asiaticus, while 1,752 SNPs for assessing genetic diversity in six cultivars of A. coronaria. UPGMA dendrograms were constructed and in R. asiaticus integrated in with PCA analysis. This study reports the first molecular fingerprinting within Persian buttercup, while the results obtained in poppy anemone were compared with a previously published SSR-based fingerprinting, proving K-seq to be an efficient protocol for the genotyping of complex genetic backgrounds.

CRISPR/Cas9-Based Knock-Out of the PMR4 Gene Reduces Susceptibility to Late Blight in Two Tomato Cultivars

Phytophthora infestans, the causal agent of late blight (LB) in tomato (Solanum lycopersicum L.), is a devastating disease and a serious concern for plant productivity. The presence of susceptibility (S) genes in plants facilitates pathogen proliferation; thus, disabling these genes may help provide a broad-spectrum and durable type of tolerance/resistance. Previous studies on Arabidopsis and tomato have highlighted that knock-out mutants of the PMR4 susceptibility gene are tolerant to powdery mildew. Moreover, PMR4 knock-down in potato has been shown to confer tolerance to LB. To verify the same effect in tomato in the present study, a CRISPR-Cas9 vector containing four single guide RNAs (sgRNAs: sgRNA1, sgRNA6, sgRNA7, and sgRNA8), targeting as many SlPMR4 regions, was introduced via Agrobacterium-tumefaciens-mediated transformation into two widely grown Italian tomato cultivars: 'San Marzano' (SM) and 'Oxheart' (OX). Thirty-five plants (twenty-six SM and nine OX) were selected and screened to identify the CRISPR/Cas9-induced mutations. The different sgRNAs caused mutation frequencies ranging from 22.1 to 100% and alternatively precise insertions (sgRNA6) or deletions (sgRNA7, sgRNA1, and sgRNA8). Notably, sgRNA7 induced in seven SM genotypes a -7 bp deletion in the homozygous status, whereas sgRNA8 led to the production of fifteen SM genotypes with a biallelic mutation (-7 bp and -2 bp). Selected edited lines were inoculated with P. infestans, and four of them, fully knocked out at the PMR4 locus, showed reduced disease symptoms (reduction in susceptibility from 55 to 80%) compared to control plants. The four SM lines were sequenced using Illumina whole-genome sequencing for deeper characterization without exhibiting any evidence of mutations in the candidate off-target regions. Our results showed, for the first time, a reduced susceptibility to Phytophtora infestans in pmr4 tomato mutants confirming the role of KO PMR4 in providing broad-spectrum protection against pathogens.

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.

The Population Structure of a Globe Artichoke Worldwide Collection, as Revealed by Molecular and Phenotypic Analyzes

The knowledge of the organization of the domesticated gene pool of crop species is an essential requirement to understand crop evolution, to rationalize conservation programs, and to support practical decisions in plant breeding. Here, we integrate simple sequence repeat (SSR) analysis and phenotypic characterization to investigate a globe artichoke collection that comprises most of the varieties cultivated worldwide. We show that the cultivated gene pool of globe artichoke includes five distinct genetic groups associated with the major phenotypic typologies: Catanesi (which based on our analysis corresponds to Violetti di Provenza), Spinosi, Violetti di Toscana, Romaneschi, and Macau. We observed that 17 and 11% of the molecular and phenotypic variance, respectively, is between these groups, while within groups, strong linkage disequilibrium and heterozygote excess are evident. The divergence between groups for quantitative traits correlates with the average broad-sense heritability within the groups. The phenotypic divergence between groups for both qualitative and quantitative traits is strongly and positively correlated with SSR divergence (F_ST) between groups. All this implies a low population size and strong bottleneck effects, and indicates a long history of clonal propagation and selection during the evolution of the domesticated gene pool of globe artichoke. Moreover, the comparison between molecular and phenotypic population structures suggests that harvest time, plant architecture (i.e., plant height, stem length), leaf spininess, head morphology (i.e., head shape, bract shape, spininess) together with the number of heads per plant were the main targets of selection during the evolution of the cultivated germplasm. We emphasize our findings in light of the potential exploitation of this collection for association mapping studies.

In-Depth Characterization of greenflesh Tomato Mutants Obtained by CRISPR/Cas9 Editing: A Case Study With Implications for Breeding and Regulation

Gene editing has already proved itself as an invaluable tool for the generation of mutants for crop breeding, yet its ultimate impact on agriculture will depend on how crops generated by gene editing technologies are regulated, and on our ability to characterize the impact of mutations on plant phenotype. A starting operational strategy for evaluating gene editing-based approaches to plant breeding might consist of assessing the effect of the induced mutations in a crop- and locus-specific manner: this involves the analysis of editing efficiency in different cultivars of a crop, the assessment of potential off-target mutations, and a phenotypic evaluation of edited lines carrying different mutated alleles. Here, we targeted the GREENFLESH (GF) locus in two tomato cultivars ('MoneyMaker' and 'San Marzano') and evaluated the efficiency, specificity and mutation patterns associated with CRISPR/Cas9 activity for this gene. The GF locus encodes a Mg-dechelatase responsible for initiating chlorophyll degradation; in gf mutants, ripe fruits accumulate both carotenoids and chlorophylls. Phenotypic evaluations were conducted on two transgene-free T₂ 'MoneyMaker' gf lines with different mutant alleles (a small insertion of 1 nucleotide and a larger deletion of 123 bp). Both lines, in addition to reduced chlorophyll degradation, showed a notable increase in carotenoid and tocopherol levels during fruit ripening. Infection of gf leaves and fruits with Botrytis cinerea resulted in a significant reduction of infected area and pathogen proliferation compared to the wild type (WT). Our data indicates that the CRISPR/Cas9-mediated mutation of the GF locus in tomato is efficient, specific and reproducible and that the resulting phenotype is robust and consistent with previously characterized greenflesh mutants obtained with different breeding techniques, while also shedding light on novel traits such as vitamin E overaccumulation and pathogen resistance. This makes GF an appealing target for breeding tomato cultivars with improved features for cultivation, as well as consumer appreciation and health.

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.

Oleosin Cor a 15 is a novel allergen for Italian hazelnut allergic children

BACKGROUND

Hazelnut allergy, which is characterized by symptoms that range from mild to severe, is one of the most common allergies in children throughout Europe, and an accurate diagnosis of this allergy is therefore essential. However, lipophilic allergens, such as oleosins, are generally underrepresented in diagnostic tests. We therefore sought to characterize the IgE reactivity of raw and roasted hazelnut oleosins, using the sera of hazelnut-allergic pediatric patients.

METHODS

Raw and roasted hazelnut oil body-associated proteins were analyzed by means of 1D and 2D electrophoresis and MS. Oleosin IgE reactivity was assessed by immunoblotting with the sera of 27 children who have confirmed hazelnut allergies and from 10 tolerant subjects. A molecular characterization of the oleosins was performed by interrogating the C. avellana cv. Jefferson and cv. TGL genomes, and through expression and purification of the recombinant new allergen.

RESULTS

A proteomic and genomic investigation allowed two new oleosins to be identified, in addition to Cor a 12 and Cor a 13, in hazelnut oil bodies. One of the new oleosins was registered as a new allergen, according to the WHO/IUIS Allergen Nomenclature Subcommittee criteria, and termed Cor a 15. Cor a 15 was the most frequently immunorecognized oleosin in our cohort. Oleosins resulted to be the only immunorecognized allergens in a subgroup of allergic patients who showed low ImmunoCAP assay IgE values and positive OFC and PbP. Hazelnut roasting resulted in an increase in oleosin immunoreactivity.

CONCLUSION

A novel hazelnut oleosin, named Cor a 15, has been discovered. Cor a 15 could play a role in eliciting an allergic reaction in a subgroup of pediatric patients that exclusively immunorecognize oleosins. The high prevalence of hazelnut oleosin sensitization here reported further confirms the need to include oleosins in routine diagnostic procedures.

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.

Mapping the Genetic Regions Responsible for Key Phenology-Related Traits in the European Hazelnut

An increasing interest in the cultivation of (European) hazelnut (Corylus avellana) is driving a demand to breed cultivars adapted to non-conventional environments, particularly in the context of incipient climate change. Given that plant phenology is so strongly determined by genotype, a rational approach to support these breeding efforts will be to identify quantitative trait loci (QTLs) and the genes underlying the basis for adaptation. The present study was designed to map QTLs for phenology-related traits, such as the timing of both male and female flowering, dichogamy, and the period required for nuts to reach maturity. The analysis took advantage of an existing linkage map developed from a population of F₁ progeny bred from the cross "Tonda Gentile delle Langhe" × "Merveille de Bollwiller," consisting in 11 LG. A total of 42 QTL-harboring regions were identified. Overall, 71 QTLs were detected, 49 on the TGdL map and 22 on the MB map; among these, 21 were classified as major; 13 were detected in at least two of the seasons (stable-major QTL). In detail, 20 QTLs were identified as contributing to the time of male flowering, 15 to time of female flowering, 25 to dichogamy, and 11 to time of nut maturity. LG02 was found to harbor 16 QTLs, while 15 QTLs mapped to LG10 and 14 to LG03. Many of the QTLs were clustered with one another. The major cluster was located on TGdL_02 and consisted of mainly major QTLs governing all the analyzed traits. A search of the key genomic regions revealed 22 candidate genes underlying the set of traits being investigated. Many of them have been described in the literature as involved in processes related to flowering, control of dormancy, budburst, the switch from vegetative to reproductive growth, or the morphogenesis of flowers and seeds.

Development of High-Density Genetic Linkage Maps and Identification of Loci for Chestnut Gall Wasp Resistance in Castanea spp

Castanea sativa is an important multipurpose species in Europe for nut and timber production as well as for its role in the landscape and in the forest ecosystem. This species has low tolerance to chestnut gall wasp (Dryocosmus kuriphilus Yasumatsu), which is a pest that was accidentally introduced into Europe in early 2000 and devastated forest and orchard trees. Resistance to the gall wasp was found in the hybrid cultivar ‘Bouche de Bétizac’ (C. sativa × C. crenata) and studied by developing genetic linkage maps using a population derived from a cross between ‘Bouche de Bétizac’ and the susceptible cultivar ‘Madonna’ (C. sativa). The high-density genetic maps were constructed using double-digest restriction site-associated DNA-seq and simple sequence repeat markers. The map of ‘Bouche de Bétizac’ consisted of 1459 loci and spanned 809.6 cM; the map of ‘Madonna’ consisted of 1089 loci and spanned 753.3 cM. In both maps, 12 linkage groups were identified. A single major QTL was recognized on the ‘Bouche de Bétizac’ map, explaining up to 67–69% of the phenotypic variance of the resistance trait (Rdk1). The Rdk1 quantitative trait loci (QTL) region included 11 scaffolds and two candidate genes putatively involved in the resistance response were identified. This study will contribute to C. sativa breeding programs and to the study of Rdk1 genes.

Identification of a caleosin associated with hazelnut (Corylus avellana L.) oil bodies

Caleosins are involved in several cellular and biological processes that are closely associated with the synthesis, degradation and stability of oil bodies (OB). Because of the importance and the multiple roles of these OB-associated proteins, in silico identification of sequences corresponding to putative caleosins in the hazelnut genome has been performed, and the association with seed OB was verified using a proteomic approach. Five full-length sequences (CavCLO-H1, CavCLO-H2, CavCLO-H3, CavCLO-L1, CavCLO-L2), belonging to the two groups of caleosins (H and L), have been identified in the hazelnut genome. The number of identified caleosins is in agreement with that previously observed in other plant species, confirming that caleosins comprise small gene families in plants. A proteomic approach allowed us to verify only the presence of CavCLO-H1 in hazelnut OB, suggesting that several members inside this family could have different roles during plant growth and development. In silico analysis also suggests that CavCLO-H1 may act as a peroxygenase.

Different Phenotypes, Similar Genomes: Three Newly Sequenced Fusarium fujikuroi Strains Induce Different Symptoms in Rice Depending on Temperature

Bakanae, caused by the hemibiotrophic fungus Fusarium fujikuroi, is one of the most important diseases of rice and is attributed to up to 75% of losses, depending on the strain and environmental conditions. Some strains cause elongation and thin leaves, whereas others induce stunting and chlorotic seedlings. Differences in symptoms are attributed to genetic differences in the strains. F. fujikuroi strains Augusto2, CSV1, and I1.3 were sequenced with Illumina MiSeq, and pathogenicity trials were conducted on rice cultivar Galileo, which is susceptible to bakanae. By performing gene prediction, single nucleotide polymorphism (SNP) calling, and structural variant analysis with a reference genome, we show how an extremely limited number of polymorphisms in genes not commonly associated with bakanae disease can cause strong differences in phenotype. CSV1 and Augusto2 were particularly close, with only 21,887 SNPs between them, but they differed in virulence, reaction to temperature, induced symptoms, colony morphology and color, growth speed, fumonisin, and gibberellin production. Genes potentially involved in the shift in phenotype were identified. Furthermore, we show how temperature variation may result in different symptoms even in rice plants inoculated with the same F. fujikuroi strain. Moreover, all of the F. fujikuroi strains became more virulent at higher temperatures. Significant differences were likewise observed in gibberellic acid production and in the expression of both fungal and plant gibberellin biosynthetic genes.

MCSeEd (Methylation Context Sensitive Enzyme ddRAD): A New Method to Analyze DNA Methylation

Methylation context sensitive enzyme ddRAD (MCSeEd) is a NGS-based method for genome-wide investigations of DNA methylation at different contexts requiring only low to moderate sequencing depth. It is particularly useful for identifying methylation changes in experimental systems challenged by biotic or abiotic stresses or at different developmental stages.

Simultaneous CRISPR/Cas9 Editing of Three PPO Genes Reduces Fruit Flesh Browning in Solanum melongena L

Polyphenol oxidases (PPOs) catalyze the oxidization of polyphenols, which in turn causes the browning of the eggplant berry flesh after cutting. This has a negative impact on fruit quality for both industrial transformation and fresh consumption. Ten PPO genes (named SmelPPO1-10) were identified in eggplant thanks to the recent availability of a high-quality genome sequence. A CRISPR/Cas9-based mutagenesis approach was applied to knock-out three target PPO genes (SmelPPO4, SmelPPO5, and SmelPPO6), which showed high transcript levels in the fruit after cutting. An optimized transformation protocol for eggplant cotyledons was used to obtain plants in which Cas9 is directed to a conserved region shared by the three PPO genes. The successful editing of the SmelPPO4, SmelPPO5, and SmelPPO6 loci of in vitro regenerated plantlets was confirmed by Illumina deep sequencing of amplicons of the target sites. Besides, deep sequencing of amplicons of the potential off-target loci identified in silico proved the absence of detectable non-specific mutations. The induced mutations were stably inherited in the T₁ and T₂ progeny and were associated with a reduced PPO activity and browning of the berry flesh after cutting. Our results provide the first example of the use of the CRISPR/Cas9 system in eggplant for biotechnological applications and open the way to the development of eggplant genotypes with low flesh browning which maintain a high polyphenol content in the berries.

Identification of DNA methyltransferases and demethylases in Solanum melongena L., and their transcription dynamics during fruit development and after salt and drought stresses

DNA methylation through the activity of cytosine-5-methyltransferases (C5-MTases) and DNA demethylases plays important roles in genome protection as well as in regulating gene expression during plant development and plant response to environmental stresses. In this study, we report on a genome-wide identification of six C5-MTases (SmelMET1, SmelCMT2, SmelCMT3a, SmelCMT3b, SmelDRM2, SmelDRM3) and five demethylases (SmelDemethylase_1, SmelDemethylase_2, SmelDemethylase_3, SmelDemethylase_4, SmelDemethylase_5) in eggplant. Gene structural characteristics, chromosomal localization and phylogenetic analyses are also described. The transcript profiling of both C5-MTases and demethylases was assessed at three stages of fruit development in three eggplant commercial F1 hybrids: i.e. 'Clara', 'Nite Lady' and 'Bella Roma', representative of the eggplant berry phenotypic variation. The trend of activation of C5-MTases and demethylase genes varied in function of the stage of fruit development and was genotype dependent. The transcription pattern of C5MTAses and demethylases was also assessed in leaves of the F1 hybrid 'Nite Lady' subjected to salt and drought stresses. A marked up-regulation and down-regulation of some C5-MTases and demethylases was detected, while others did not vary in their expression profile. Our results suggest a role for both C5-MTases and demethylases during fruit development, as well as in response to abiotic stresses in eggplant, and provide a starting framework for supporting future epigenetic studies in the species.

Transcriptome characterization and expression profiling in chestnut cultivars resistant or susceptible to the gall wasp Dryocosmus kuriphilus

The oriental gall wasp Dryocosmus kuriphilus represents a limiting pest for the European Chestnut (Castanea sativa, Fagaceae) as it creates severe yield losses. The European Chestnut is a deciduous tree, having major social, economic and environmental importance in Southern Europe, covering an area of 2.53 million hectares, including 75,000 ha devoted to fruit production. Cultivars show different susceptibility and very few are resistant to gall wasp. To deeply investigate the plant response and understand which factors can lead the plant to develop or not the gall, the study of transcriptome is basic (fundamental). To date, little transcriptomic information are available for C. sativa species. Hence, we present a de novo assembly of the chestnut transcriptome of the resistant Euro-Japanese hybrid 'Bouche de Bétizac' (BB) and the susceptible cultivar 'Madonna' (M), collecting RNA from buds at different stages of budburst. The two transcriptomes were assembled into 34,081 (BB) and 30,605 (M) unigenes, respectively. The former was used as a reference sequence for further characterization analyses, highlighting the presence of 1444 putative resistance gene analogs (RGAs) and about 1135 unigenes, as putative MiRNA targets. A global quantitative transcriptome profiling comparing the resistant and the susceptible cultivars, in the presence or not of the gall wasp, revealed some GO enrichments as "response to stimulus" (GO:0050896), and "developmental processes" (e.g., post-embryonic development, GO:0009791). Many up-regulated genes appeared to be transcription factors (e.g., RAV1, AP2/ERF, WRKY33) or protein regulators (e.g., RAPTOR1B) and storage proteins (e.g., LEA D29) involved in "post-embryonic development". Our analysis was able to provide a large amount of information, including 7k simple sequence repeat (SSR) and 335k single-nucleotide polymorphism (SNP)/INDEL markers, and generated the first reference unigene catalog for the European Chestnut. The transcriptome data for C. sativa will contribute to understand the genetic basis of the resistance to gall wasp and will provide useful information for next molecular genetic studies of this species and its relatives.

A chromosome-anchored eggplant genome sequence reveals key events in Solanaceae evolution

With approximately 450 species, spiny Solanum species constitute the largest monophyletic group in the Solanaceae family, but a high-quality genome assembly from this group is presently missing. We obtained a chromosome-anchored genome assembly of eggplant (Solanum melongena), containing 34,916 genes, confirming that the diploid gene number in the Solanaceae is around 35,000. Comparative genomic studies with tomato (S. lycopersicum), potato (S. tuberosum) and pepper (Capsicum annuum) highlighted the rapid evolution of miRNA:mRNA regulatory pairs and R-type defense genes in the Solanaceae, and provided a genomic basis for the lack of steroidal glycoalkaloid compounds in the Capsicum genus. Using parsimony methods, we reconstructed the putative chromosomal complements of the key founders of the main Solanaceae clades and the rearrangements that led to the karyotypes of extant species and their ancestors. From 10% to 15% of the genes present in the four genomes were syntenic paralogs (ohnologs) generated by the pre-γ, γ and T paleopolyploidy events, and were enriched in transcription factors. Our data suggest that the basic gene network controlling fruit ripening is conserved in different Solanaceae clades, and that climacteric fruit ripening involves a differential regulation of relatively few components of this network, including CNR and ethylene biosynthetic genes.

Single Primer Enrichment Technology (SPET) for High-Throughput Genotyping in Tomato and Eggplant Germplasm

Single primer enrichment technology (SPET) is a new, robust, and customizable solution for targeted genotyping. Unlike genotyping by sequencing (GBS), and like DNA chips, SPET is a targeted genotyping technology, relying on the sequencing of a region flanking a primer. Its reliance on single primers, rather than on primer pairs, greatly simplifies panel design, and allows higher levels of multiplexing than PCR-based genotyping. Thanks to the sequencing of the regions surrounding the target SNP, SPET allows the discovery of thousands of closely linked, novel SNPs. In order to assess the potential of SPET for high-throughput genotyping in plants, a panel comprising 5k target SNPs, designed both on coding regions and introns/UTRs, was developed for tomato and eggplant. Genotyping of two panels composed of 400 tomato and 422 eggplant accessions, comprising both domesticated material and wild relatives, generated a total of 12,002 and 30,731 high confidence SNPs, respectively, which comprised both target and novel SNPs in an approximate ratio of 1:1.6, and 1:5.5 in tomato and eggplant, respectively. The vast majority of the markers was transferrable to related species that diverged up to 3.4 million years ago (Solanum pennellii for tomato and S. macrocarpon for eggplant). Maximum Likelihood phylogenetic trees and PCA outputs obtained from the whole dataset highlighted genetic relationships among accessions and species which were congruent with what was previously reported in literature. Better discrimination among domesticated accessions was achieved by using the target SNPs, while better discrimination among wild species was achieved using the whole SNP dataset. Our results reveal that SPET genotyping is a robust, high-throughput technology for genetic fingerprinting, with a high degree of cross-transferability between crops and their cultivated and wild relatives, and allows identification of duplicates and mislabeled accessions in genebanks.

Analysis of DNA Methylation Patterns Associated with In Vitro Propagated Globe Artichoke Plants Using an EpiRADseq-Based Approach

Globe artichoke represents one of the main horticultural species of the Mediterranean basin, and 'Spinoso sardo' is the most widespread and economically relevant varietal type in Sardinia, Italy. In the last decades, in vitro culture of meristematic apices has increased the frequency of aberrant plants in open-field production. These off-type phenotypes showed highly pinnate-parted leaves and late inflorescence budding, and emerged from some branches of the true-to-type 'Spinoso sardo' plants. This phenomenon cannot be foreseen and is reversible through generations, suggesting the occurrence of epigenetic alterations. Here, we report an exploratory study on DNA methylation patterns in off-type/true-to-type globe artichoke plants, using a modified EpiRADseq technology, which allowed the identification of 2,897 differentially methylated loci (DML): 1,998 in CG, 458 in CHH, and 441 in CHG methylation contexts of which 720, 88, and 152, respectively, were in coding regions. Most of them appeared involved in primary metabolic processes, mostly linked to photosynthesis, regulation of flower development, and regulation of reproductive processes, coherently with the observed phenotype. Differences in the methylation status of some candidate genes were integrated with transcriptional analysis to test whether these two regulation levels might interplay in the emergence and spread of the 'Spinoso sardo' non-conventional phenotype.

High density SNP mapping and QTL analysis for time of leaf budburst in Corylus avellana L

The growing area of European hazelnut (Corylus avellana L.) is increasing, as well as the number of producing countries, and there is a pressing need for new improved cultivars. Hazelnut conventional breeding process is slow, due to the length of juvenile phase and the high heterozygosity level. The development of genetic linkage maps and the identification of molecular markers tightly linked to QTL (quantitative trait loci) of agronomic interest are essential tools for speeding up the selection of seedlings carrying desired traits through marker-assisted selection. The objectives of this study were to enrich a previous linkage map and confirm QTL related to time of leaf budburst, using an F₁ population obtained by crossing Tonda Gentile delle Langhe with Merveille de Bollwiller. Genotyping-by-Sequencing was used to identify a total of 9,999 single nucleotide polymorphism markers. Well saturated linkage maps were constructed for each parent using the double pseudo-testcross mapping strategy. A reciprocal translocation was detected in Tonda Gentile delle Langhe between two non-homologous chromosomes. Applying a bioinformatic approach, we were able to disentangle ‘pseudo-linkage’ between markers, removing markers around the translocation breakpoints and obtain a linear order of the markers for the two chromosomes arms, for each linkage group involved in the translocation. Twenty-nine QTL for time of leaf budburst were identified, including a stably expressed region on LG_02 of the Tonda Gentile delle Langhe map. The stability of these QTL and their coding sequence content indicates promise for the identification of specific chromosomal regions carrying key genes involved in leaf budburst.

Comprehensive Characterization of Simple Sequence Repeats in Eggplant (Solanum melongena L.) Genome and Construction of a Web Resource

We have characterized the simple sequence repeat (SSR) markers of the eggplant (Solanum melongena) using a recent high quality sequence of its whole genome. We found nearly 133,000 perfect SSRs, a density of 125.5 SSRs/Mbp, and also about 178,400 imperfect SSRs. Of the perfect SSRs, 15.6% were complex, with two stretches of repeats separated by an intervening block of <100 nt. Di- and trinucleotide SSRs accounted, respectively, for 43 and 37% of the total. The SSRs were classified according to their number of repeats and overall length, and were assigned to their linkage group. We found 2,449 of the perfect SSRs in 2,086 genes, with an overall density of 18.5 SSRs/Mbp across the gene space; 3,524 imperfect SSRs were present in 2,924 genes at a density of 26.7 SSRs/Mbp. Putative functions were assigned via ontology to genes containing at least one SSR. Using this data we developed an "Eggplant Microsatellite DataBase" (EgMiDB) which permits identification of SSR markers in terms of their location on the genome, type of repeat (perfect vs. imperfect), motif type, sequence, repeat number and genomic/gene context. It also suggests forward and reverse primers. We employed an in silico PCR analysis to validate these SSR markers, using as templates two CDS sets and three assembled transcriptomes obtained from diverse eggplant accessions.

Genomic assessment in Lactobacillus plantarum links the butyrogenic pathway with glutamine metabolism

The butyrogenic capability of Lactobacillus (L.) plantarum is highly dependent on the substrate type and so far not assigned to any specific metabolic pathway. Accordingly, we compared three genomes of L. plantarum that showed a strain-specific capability to produce butyric acid in human cells growth media. Based on the genomic analysis, butyric acid production was attributed to the complementary activities of a medium-chain thioesterase and the fatty acid synthase of type two (FASII). However, the genomic islands of discrepancy observed between butyrogenic L. plantarum strains (S2T10D, S11T3E) and the non-butyrogenic strain O2T60C do not encompass genes of FASII, but several cassettes of genes related to sugar metabolism, bacteriocins, prophages and surface proteins. Interestingly, single amino acid substitutions predicted from SNPs analysis have highlighted deleterious mutations in key genes of glutamine metabolism in L. plantarum O2T60C, which corroborated well with the metabolic deficiency suffered by O2T60C in high-glutamine growth media and its consequent incapability to produce butyrate. In parallel, the increase of glutamine content induced the production of butyric acid by L. plantarum S2T10D. The present study reveals a previously undescribed metabolic route for butyric acid production in L. plantarum, and a potential involvement of the glutamine uptake in its regulation.

The genome-wide identification and transcriptional levels of DNA methyltransferases and demethylases in globe artichoke

Changes to the cytosine methylation status of DNA, driven by the activity of C5 methyltransferases (C5-MTases) and demethylases, exert an important influence over development, transposon movement, gene expression and imprinting. Three groups of C5-MTase enzymes have been identified in plants, namely MET (methyltransferase 1), CMT (chromomethyltransferases) and DRM (domains rearranged methyltransferases). Here the repertoire of genes encoding C5-MTase and demethylase by the globe artichoke (Cynara cardunculus var. scolymus) is described, based on sequence homology, a phylogenetic analysis and a characterization of their functional domains. A total of ten genes encoding C5-MTase (one MET, five CMTs and four DRMs) and five demethylases was identified. An analysis of their predicted product's protein structure suggested an extensive level of conservation has been retained by the C5-MTases. Transcriptional profiling based on quantitative real time PCR revealed a number of differences between the genes encoding maintenance and de novo methyltransferases, sometimes in a tissue- or development-dependent manner, which implied a degree of functional specialization.

Coding SNPs analysis highlights genetic relationships and evolution pattern in eggplant complexes

Brinjal (Solanum melongena), scarlet (S. aethiopicum) and gboma (S. macrocarpon) eggplants are three Old World domesticates. The genomic DNA of a collection of accessions belonging to the three cultivated species, along with a representation of various wild relatives, was characterized for the presence of single nucleotide polymorphisms (SNPs) using a genotype-by-sequencing approach. A total of 210 million useful reads were produced and were successfully aligned to the reference eggplant genome sequence. Out of the 75,399 polymorphic sites identified among the 76 entries in study, 12,859 were associated with coding sequence. A genetic relationships analysis, supported by the output of the FastSTRUCTURE software, identified four major sub-groups as present in the germplasm panel. The first of these clustered S. aethiopicum with its wild ancestor S. anguivi; the second, S. melongena, its wild progenitor S. insanum, and its relatives S. incanum, S. lichtensteinii and S. linneanum; the third, S. macrocarpon and its wild ancestor S. dasyphyllum; and the fourth, the New World species S. sisymbriifolium, S. torvum and S. elaeagnifolium. By applying a hierarchical FastSTRUCTURE analysis on partitioned data, it was also possible to resolve the ambiguous membership of the accessions of S. campylacanthum, S. violaceum, S. lidii, S. vespertilio and S. tomentsum, as well as to genetically differentiate the three species of New World Origin. A principal coordinates analysis performed both on the entire germplasm panel and also separately on the entries belonging to sub-groups revealed a clear separation among species, although not between each of the domesticates and their respective wild ancestors. There was no clear differentiation between either distinct cultivar groups or different geographical provenance. Adopting various approaches to analyze SNP variation provided support for interpretation of results. The genotyping-by-sequencing approach showed to be highly efficient for both quantifying genetic diversity and establishing genetic relationships among and within cultivated eggplants and their wild relatives. The relevance of these results to the evolution of eggplants, as well as to their genetic improvement, is discussed.

Genome-Wide Identification of BAHD Acyltransferases and In vivo Characterization of HQT-like Enzymes Involved in Caffeoylquinic Acid Synthesis in Globe Artichoke

Globe artichoke (Cynara cardunculus L. var. scolymus) is a rich source of compounds promoting human health (phytonutrients), among them caffeoylquinic acids (CQAs), mainly represented by chlorogenic acid (CGA), and dicaffeoylquinic acids (diCQAs). The enzymes involved in their biosynthesis belong to the large family of BAHD acyltransferases. Following a survey of the globe artichoke genome, we identified 69 BAHD proteins carrying the catalytic site (HXXXD). Their phylogenetic analysis together with another 43 proteins, from 21 species, representative of the BAHD family, highlighted their grouping in seven major clades. Nine globe artichoke acyltransferases clustered in a sub-group of Clade V, with 3 belonging to hydroxycinnamoyl-CoA:quinate hydroxycinnamoyl transferase (HQT) and 2 to hydroxycinnamoyl-CoA:shikimate/quinate hydroxycinnamoyl transferase (HCT) like proteins. We focused our attention on the former, HQT1, HQT2, and HQT3, as they are known to play a key role in CGA biosynthesis. The expression of genes coding for the three HQTs and correlation of expression with the CQA content is reported for different globe artichoke tissues. For the first time in the globe artichoke, we developed and applied the virus-induced gene silencing approach with the goal of assessing in vivo the effect of HQT1 silencing, which resulted in a marked reduction of both CGA and diCQAs. On the other hand, when the role of the three HQTs was assessed in leaves of Nicotiana benthamiana through their transient overexpression, significant increases in mono- and diCQAs content were observed. Using transient GFP fusion proteins expressed in N. benthamiana leaves we also established the sub-cellular localization of these three enzymes.

Large-scale transcriptome characterization and mass discovery of SNPs in globe artichoke and its related taxa

Cynara cardunculus (2n = 2× = 34) is a member of the Asteraceae family that contributes significantly to the agricultural economy of the Mediterranean basin. The species includes two cultivated varieties, globe artichoke and cardoon, which are grown mainly for food. Cynara cardunculus is an orphan crop species whose genome/transcriptome has been relatively unexplored, especially in comparison to other Asteraceae crops. Hence, there is a significant need to improve its genomic resources through the identification of novel genes and sequence-based markers, to design new breeding schemes aimed at increasing quality and crop productivity. We report the outcome of cDNA sequencing and assembly for eleven accessions of C. cardunculus. Sequencing of three mapping parental genotypes using Roche 454-Titanium technology generated 1.7 × 10⁶ reads, which were assembled into 38,726 reference transcripts covering 32 Mbp. Putative enzyme-encoding genes were annotated using the KEGG-database. Transcription factors and candidate resistance genes were surveyed as well. Paired-end sequencing was done for cDNA libraries of eight other representative C. cardunculus accessions on an Illumina Genome Analyzer IIx, generating 46 × 10⁶ reads. Alignment of the IGA and 454 reads to reference transcripts led to the identification of 195,400 SNPs with a Bayesian probability exceeding 95%; a validation rate of 90% was obtained by Sanger-sequencing of a subset of contigs. These results demonstrate that the integration of data from different NGS platforms enables large-scale transcriptome characterization, along with massive SNP discovery. This information will contribute to the dissection of key agricultural traits in C. cardunculus and facilitate the implementation of marker-assisted selection programs.

2-D DIGE analysis of UV-C radiation-responsive proteins in globe artichoke leaves

Plants respond to ultraviolet stress inducing a self-defence through the regulation of specific gene family members. The UV acclimation is the result of biochemical and physiological processes, such as enhancement of the antioxidant enzymatic system and accumulation of UV-absorbing phenolic compounds (e.g. flavonoids). Globe artichoke is an attractive species for studying the protein network involved in UV stress response, being characterized by remarkable levels of inducible antioxidants. Proteomic tools can assist the evaluation of the expression patterns of UV-responsive proteins and we applied the difference in-gel electrophoresis (DIGE) technology for monitoring the globe artichoke proteome variation at four time points following an acute UV-C exposure. A total of 145 UV-C-modulated proteins were observed and 119 were identified by LC-MS/MS using a ∼144,000 customized Compositae protein database, which included about 19,000 globe artichoke unigenes. Proteins were Gene Ontology (GO) categorized, visualized on their pathways and their behaviour was discussed. A predicted protein interaction network was produced and highly connected hub-like proteins were highlighted. Most of the proteins differentially modulated were chloroplast located, involved in photosynthesis, sugar metabolisms, protein folding and abiotic stress. The identification of UV-C-responsive proteins may contribute to shed light on the molecular mechanisms underlying plant responses to UV stress.

RAD tag sequencing as a source of SNP markers in Cynara cardunculus L

Background

The globe artichoke (Cynara cardunculus L. var. scolymus) genome is relatively poorly explored, especially compared to those of the other major Asteraceae crops sunflower and lettuce. No SNP markers are in the public domain. We have combined the recently developed restriction-site associated DNA (RAD) approach with the Illumina DNA sequencing platform to effect the rapid and mass discovery of SNP markers for C. cardunculus.

Results

RAD tags were sequenced from the genomic DNA of three C. cardunculus mapping population parents, generating 9.7 million reads, corresponding to ~1 Gbp of sequence. An assembly based on paired ends produced ~6.0 Mbp of genomic sequence, separated into ~19,000 contigs (mean length 312 bp), of which ~21% were fragments of putative coding sequence. The shared sequences allowed for the discovery of ~34,000 SNPs and nearly 800 indels, equivalent to a SNP frequency of 5.6 per 1,000 nt, and an indel frequency of 0.2 per 1,000 nt. A sample of heterozygous SNP loci was mapped by CAPS assays and this exercise provided validation of our mining criteria. The repetitive fraction of the genome had a high representation of retrotransposon sequence, followed by simple repeats, AT-low complexity regions and mobile DNA elements. The genomic k-mers distribution and CpG rate of C. cardunculus, compared with data derived from three whole genome-sequenced dicots species, provided a further evidence of the random representation of the C. cardunculus genome generated by RAD sampling.

Conclusion

The RAD tag sequencing approach is a cost-effective and rapid method to develop SNP markers in a highly heterozygous species. Our approach permitted to generate a large and robust SNP datasets by the adoption of optimized filtering criteria.

Strain dependent expression of stress response and virulence genes of Listeria monocytogenes in meat juices as determined by microarray

A subgenomic array, encompassing 54 probes targeting genes responsible for virulence, adhesion and stress response in Listeria monocytogenes, was used in order to study their expression in food systems. RNA extracted from L. monocytogenes inoculated in BHI and in situ (i.e. in minced meat and fermented sausage juices) and incubated at 4°C, was hybridized on the array and the results obtained were compared in order to understand the effect that the food juice has on the expression. Three different strains of L. monocytogenes were tested, in order to determine the effect of the strain provenience. As determined by cluster analysis, each strain behaved in a different way when inoculated in food juices. The goal was to respond to acidic and osmotic stresses encountered in the food, particularly in the fermented sausage juice. No differences in the expression profile between the three strains were observed, when they were inoculated in BHI. On the other hand, in the meat and sausage juices, the iap, gadC and gadE genes, together with different internalin encoding genes, were significantly differentially expressed in the three strains.

Identification of SNP and SSR markers in eggplant using RAD tag sequencing

Background

The eggplant (Solanum melongena L.) genome is relatively unexplored, especially compared to those of the other major Solanaceae crops tomato and potato. In particular, no SNP markers are publicly available; on the other hand, over 1,000 SSR markers were developed and publicly available. We have combined the recently developed Restriction-site Associated DNA (RAD) approach with Illumina DNA sequencing for rapid and mass discovery of both SNP and SSR markers for eggplant.

Results

RAD tags were generated from the genomic DNA of a pair of eggplant mapping parents, and sequenced to produce ~17.5 Mb of sequences arrangeable into ~78,000 contigs. The resulting non-redundant genomic sequence dataset consisted of ~45,000 sequences, of which ~29% were putative coding sequences and ~70% were in common between the mapping parents. The shared sequences allowed the discovery of ~10,000 SNPs and nearly 1,000 indels, equivalent to a SNP frequency of 0.8 per Kb and an indel frequency of 0.07 per Kb. Over 2,000 of the SNPs are likely to be mappable via the Illumina GoldenGate assay. A subset of 384 SNPs was used to successfully fingerprint a panel of eggplant germplasm, producing a set of informative diversity data. The RAD sequences also included nearly 2,000 putative SSRs, and primer pairs were designed to amplify 1,155 loci.

Conclusion

The high throughput sequencing of the RAD tags allowed the discovery of a large number of DNA markers, which will prove useful for extending our current knowledge of the genome organization of eggplant, for assisting in marker-aided selection and for carrying out comparative genomic analyses within the Solanaceae family.

Construction of a reference molecular linkage map of globe artichoke (Cynara cardunculus var. scolymus)

The genome organization of globe artichoke (Cynara cardunculus var. scolymus), unlike other species belonging to Asteraceae (=Compositae) family (i.e. sunflower, lettuce and chicory), remains largely unexplored. The species is highly heterozygous and suffers marked inbreeding depression when forced to self-fertilize. Thus a two-way pseudo-testcross represents the optimal strategy for linkage analysis. Here, we report linkage maps based on the progeny of a cross between globe artichoke (C. cardunculus var. scolymus) and cultivated cardoon (C. cardunculus var. altilis). The population was genotyped using a variety of PCR-based marker platforms, resulting in the identification of 708 testcross markers suitable for map construction. The male map consisted of 177 loci arranged in 17 major linkage groups, spanning 1,015.5 cM, while female map was built with 326 loci arranged into 20 major linkage groups, spanning 1,486.8 cM. The presence of 84 loci shared between these maps and those previously developed from a cross within globe artichoke allowed for map alignment and the definition of 17 homologous linkage groups, corresponding to the haploid number of the species. This will provide a favourable property for QTL scanning; furthermore, as 25 mapped markers (8%) correspond to coding regions, it has an additional value as functional map and might represent an important genetic tool for candidate gene studies in globe artichoke.

Ontology and diversity of transcript-associated microsatellites mined from a globe artichoke EST database

Background

The globe artichoke (Cynara cardunculus var. scolymus L.) is a significant crop in the Mediterranean basin. Despite its commercial importance and its both dietary and pharmaceutical value, knowledge of its genetics and genomics remains scant. Microsatellite markers have become a key tool in genetic and genomic analysis, and we have exploited recently acquired EST (expressed sequence tag) sequence data (Composite Genome Project - CGP) to develop an extensive set of microsatellite markers.

Results

A unigene assembly was created from over 36,000 globe artichoke EST sequences, containing 6,621 contigs and 12,434 singletons. Over 12,000 of these unigenes were functionally assigned on the basis of homology with Arabidopsis thaliana reference proteins. A total of 4,219 perfect repeats, located within 3,308 unigenes was identified and the gene ontology (GO) analysis highlighted some GO term's enrichments among different classes of microsatellites with respect to their position. Sufficient flanking sequence was available to enable the design of primers to amplify 2,311 of these microsatellites, and a set of 300 was tested against a DNA panel derived from 28 C. cardunculus genotypes. Consistent amplification and polymorphism was obtained from 236 of these assays. Their polymorphic information content (PIC) ranged from 0.04 to 0.90 (mean 0.66). Between 176 and 198 of the assays were informative in at least one of the three available mapping populations.

Conclusion

EST-based microsatellites have provided a large set of de novo genetic markers, which show significant amounts of polymorphism both between and within the three taxa of C. cardunculus. They are thus well suited as assays for phylogenetic analysis, the construction of genetic maps, marker-assisted breeding, transcript mapping and other genomic applications in the species.

Isolation and mapping of a C3'H gene (CYP98A49) from globe artichoke, and its expression upon UV-C stress

Globe artichoke represents a natural source of phenolic compounds with dicaffeoylquinic acids along with their biosynthetic precursor chlorogenic acid (5-caffeoylquinic acid) as the predominant molecules. We report the isolation and characterization of a full-length cDNA and promoter of a globe artichoke p-coumaroyl ester 3'-hydroxylase (CYP98A49), which is involved in both chlorogenic acid and lignin biosynthesis. Phylogenetic analyses demonstrated that this gene belongs to the CYP98 family. CYP98A49 was also heterologously expressed in yeast, in order to perform an enzymatic assay with p-coumaroylshikimate and p-coumaroylquinate as substrates. Real Time quantitative PCR analysis revealed that CYP98A49 expression is induced upon exposure to UV-C radiation. A single nucleotide polymorphism in the CYP98A49 gene sequence of two globe artichoke varieties used for genetic mapping allowed the localization of this gene to linkage group 10 within the previously developed maps.

Genetic mapping and annotation of genomic microsatellites isolated from globe artichoke

Cynara cardunculus includes three taxa, the globe artichoke (subsp. scolymus L. Hegi), the cultivated cardoon (var. altilis) and their progenitor, the wild cardoon (var. sylvestris). Globe artichoke is an important component of the Mediterranean rural economy, but its improvement through breeding has been rather limited and its genome organization remains largely unexplored. Here, we report the isolation of 61 new microsatellite loci which amplified a total of 208 alleles in a panel of 22 C. cardunculus genotypes. Of these, 51 were informative for linkage analysis and 39 were used to increase marker density in the available globe artichoke genetic maps. Sequence analysis of the 22 loci associated with genes showed that 9 are located within coding sequence, with the repetitive domain probably being involved in DNA binding or in protein-protein interactions. The expression of the genes associated with 9 of the 22 microsatellite loci was demonstrated by RT-PCR.

Stress-induced biosynthesis of dicaffeoylquinic acids in globe artichoke

Leaf extracts from globe artichoke ( Cynara cardunculus L. var. scolymus) have been widely used in medicine as hepatoprotectant and choleretic agents. Globe artichoke leaves represent a natural source of phenolic acids with dicaffeoylquinic acids, such as cynarin (1,3-dicaffeoylquinic acid), along with its biosynthetic precursor chlorogenic acid (5-caffeoylquinic acid) as the most abundant molecules. This paper reports the development of an experimental system to induce caffeoylquinic acids. This system may serve to study the regulation of the biosynthesis of (poly)phenolic compounds in globe artichoke and the genetic basis of this metabolic regulation. By means of HPLC-PDA and accurate mass LC-QTOF MS and MS/MS analyses, the major phenolic compounds in globe artichoke leaves were identified: four isomers of dicaffeoylquinic acid, three isomers of caffeoylquinic acid, and the flavone luteolin 7-glucoside. Next, plant material was identified in which the concentration of phenolic compounds was comparable in the absence of particular treatments, with the aim to use this material to test the effect of stress application on the regulation of biosynthesis of caffeoylquinic acids. Using this material, the effect of UV-C, methyl jasmonate, and salicylic acid treatments on (poly)phenolic compounds was tested in different globe artichoke genotypes. UV-C exposure consistently increased the levels of dicaffeoylquinic acids in all genotypes, whereas the effect on compounds from the same biosynthetic pathway, for example, chlorogenic acid and luteolin-7-glucoside, was much less pronounced and was not statistically significant. No effect of methyl jasmonate or salicylic acid was found. Time-response experiments indicated that the level of dicaffeoylquinic acids reached a maximum at 24 h after UV radiation. On the basis of these results a role of dicaffeoylquinic acids in UV protection in globe artichoke is hypothesized.

Use of AFLP for differentiation of Metschnikowia pulcherrima strains for postharvest disease biological control

Metschnikowia pulcherrima occurs naturally on fruits, buds and floral parts of apple trees. Some strains are effective as biocontrol agents against postharvest decay of apples and other fruits. The usefulness of the amplified fragment length polymorphism (AFLP) technique was evaluated for the genetic analysis of 26 strains of M. pulcherrima, isolated from different sources in different geographical regions. With six AFLP primer pairs, 729 polymorphic bands were scored. The technique showed a high discriminatory power. Genetic relationships between strains were also estimated using AFLP. All the isolates from the carposphere of apple, previously tested as biocontrol agents, were grouped in a single cluster with a high bootstrap value (97), indicating robustness and reproducibility. AFLP patterns could clearly distinguish the different strains and research is in progress to use some putative specific bands for single tag sequence (STS) conversion to develop isolate-specific markers.

Isolation and functional characterization of a cDNA coding a hydroxycinnamoyltransferase involved in phenylpropanoid biosynthesis in Cynara cardunculus L

BACKGROUND

Cynara cardunculus L. is an edible plant of pharmaceutical interest, in particular with respect to the polyphenolic content of its leaves. It includes three taxa: globe artichoke, cultivated cardoon, and wild cardoon. The dominating phenolics are the di-caffeoylquinic acids (such as cynarin), which are largely restricted to Cynara species, along with their precursor, chlorogenic acid (CGA). The scope of this study is to better understand CGA synthesis in this plant.

RESULTS

A gene sequence encoding a hydroxycinnamoyltransferase (HCT) involved in the synthesis of CGA, was identified. Isolation of the gene sequence was achieved by using a PCR strategy with degenerated primers targeted to conserved regions of orthologous HCT sequences available. We have isolated a 717 bp cDNA which shares 84% aminoacid identity and 92% similarity with a tobacco gene responsible for the biosynthesis of CGA from p-coumaroyl-CoA and quinic acid. In silico studies revealed the globe artichoke HCT sequence clustering with one of the main acyltransferase groups (i.e. anthranilate N-hydroxycinnamoyl/benzoyltransferase). Heterologous expression of the full length HCT (GenBank accession DQ104740) cDNA in E. coli demonstrated that the recombinant enzyme efficiently synthesizes both chlorogenic acid and p-coumaroyl quinate from quinic acid and caffeoyl-CoA or p-coumaroyl-CoA, respectively, confirming its identity as a hydroxycinnamoyl-CoA: quinate HCT. Variable levels of HCT expression were shown among wild and cultivated forms of C. cardunculus subspecies. The level of expression was correlated with CGA content.

CONCLUSION

The data support the predicted involvement of the Cynara cardunculus HCT in the biosynthesis of CGA before and/or after the hydroxylation step of hydroxycinnamoyl esters.

Retrotransposon-based S-SAP as a platform for the analysis of genetic variation and linkage in globe artichoke

A high copy number of retrotransposon sequences are present and widely dispersed in plant genomes. Their activity generates a considerable degree of sequence polymorphism. Here, we report the cloning of CYRE-5, a long-terminal repeat carrying retrotransposon-like sequence in Cynara cardunculus L., and its exploitation to develop a DNA fingerprinting assay across 22 accessions, including both cultivated (globe artichoke and cultivated cardoon) and wild (wild cardoon) types. The effectiveness of the sequence-specific amplified polymorphism (S-SAP) platform is compared with that of amplified fragment length polymorphism (AFLP). A genetic linkage analysis, based on a hybrid population between 2 globe artichoke varietal types, resulted in the inclusion of 29 S-SAP loci in the core genetic map, confirming their dispersed distribution across the globe artichoke genome.

A first linkage map of globe artichoke (Cynara cardunculus var. scolymus L.) based on AFLP, S-SAP, M-AFLP and microsatellite markers

We present the first genetic maps of globe artichoke (Cynara cardunculus var. scolymus L. 2n=2x=34), constructed with a two-way pseudo-testcross strategy. A F1 mapping population of 94 individuals was generated between a late-maturing, non-spiny type and an early-maturing spiny type. The 30 AFLP, 13 M-AFLP and 9 S-SAP primer combinations chosen identified, respectively, 352, 38 and 41 polymorphic markers. Of 32 microsatellite primer pairs tested, 12 identified heterozygous loci in one or other parent, and 7 were fully informative as they segregated in both parents. The female parent map comprised 204 loci, spread over 18 linkage groups and spanned 1330.5 cM with a mean marker density of 6.5 cM. The equivalent figures for the male parent map were 180 loci, 17 linkage groups, 1239.4 and 6.9 cM. About 3% of the AFLP and AFLP-derived markers displayed segregation distortion with a P value below 0.01, and were not used for map construction. All the SSR loci were included in the linkage analysis, although one locus did show some segregation distortion. The presence of 78 markers in common to both maps allowed the alignment of 16 linkage groups. The maps generated provide a firm basis for the mapping of agriculturally relevant traits, which will then open the way for the application of a marker-assisted selection breeding strategy in this species.

Development and characterization of microsatellite markers in Cynara cardunculus L

Cynara cardunculus L. is a species native to the Mediterranean basin that comprises 2 crops, globe artichoke (var. scolymus L.) and cultivated cardoon (var. altilis DC), as well as wild cardoon (var. sylvestris (Lamk) Fiori). Globe artichoke represents an important component of the South European agricultural economy but is also cultivated in North Africa, the Near East, South America, the United States, and China. Breeding activities and molecular marker studies have been, to date, extremely limited. Better knowledge of the genome of the species might be gained by developing a range of molecular markers. Here, we report on the development of 14 microsatellites (simple sequence repeats (SSRs)) through a novel approach that we have defined as the microsatellite amplified library (MAL). The approach represents a combination of amplified fragment length polymorphism and a primer extension based enriched library, is rapid, and requires no hybridization enrichment steps. The technique provided a ~40-fold increase in the efficiency of SSR identification compared with conventional library procedures. The developed SSRs were applied for genotyping 36 accessions of C. cardunculus, including a core of 27 varietal types of globe artichoke, 3 accessions of cultivated cardoon, and 6 Sicilian accessions of wild cardoon. Principal coordinates analysis made it possible to differentiate both cultivated and wild forms from each other.Key words: globe artichoke, wild and cultivated cardoon, molecular markers, AFLP, MAL (microsatellite amplified library).