Genome reconstruction in Cynara cardunculus taxa gains access to chromosome-scale DNA variation

The Dark Side of the pollen: BSA-seq identified genomic regions linked to male sterility in globe artichoke

Globe artichoke (Cynara cardunculus var. scolymus; 2n = 2x = 34) is a food crop consumed for its immature flower heads. Traditionally, globe artichoke varietal types are vegetatively propagated. However, seed propagation makes it possible to treat the crop as annual, increasing field uniformity and reducing farmers costs, as well as pathogens diffusion. Despite globe artichoke's significant agricultural value and the critical role of heterosis in the development of superior varieties, the production of hybrids remains challenging without a reliable system for large-scale industrial seed production. Male sterility (MS) presents a promising avenue for overcoming these challenges by simplifying the hybridization process and enabling cost-effective seed production. However, within the Cynara genus, genic male sterility has been linked to three recessive loci in globe artichoke, with no definitive genetic mechanism elucidated to date. A 250 offsprings F2 population, derived from a cross between a MS globe artichoke and a male fertile (MF) cultivated cardoon (C. cardunculus var. altilis) and fitting a monogenic segregation model (3:1), was analyzed through BSA-seq, aiming at the identification of genomic regions/genes affecting male sterility. Four QTL regions were identified on chromosomes 4, 12, and 14. By analyzing the sequence around the highest pick on chromosome 14, a cytochrome P450 (CYP703A2) was identified, carrying a deleterious substitution (R/Q) fixed in the male sterile parent. A single dCAPS marker was developed around this SNP, allowing the discrimination between MS and MF genotypes within the population, suitable for applications in plant breeding programs. A 3D model of the protein was generated by homology modeling, revealing that the mutated amino acid is part of a highly conserved motif crucial for protein folding.

Bioactive Compounds from Artichoke and Application Potential

Cynara cardunculus L. var. scolymus, known as the artichoke, originated in the Mediterranean region and is now cultivated in several countries. The artichoke has leaves, a stem, and a head, also called a floral capitulum, covered with green and pointed bracts. It is rich in polyphenols, flavonoids, anthocyanins, phenolic compounds, inulin, coumarins, terpenes, dietary fibre, enzymes, polysaccharides, minerals and vitamins, and therefore has a wide range of uses, including in the food industry, medicine and biofuels. Several studies have shown that artichokes have properties such as antioxidant, anti-inflammatory, antimicrobial, anticancer, hypocholesterolaemic, anti-HIV, cardioprotective, hepatoprotective and lipid-lowering effects. The aim of this study is to provide a literature review on the phytochemical composition, bioactivity and applications, focusing on the methods of extraction, purification and concentration of enzymes present in artichoke.

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.

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.

Towards a Cardoon (Cynara cardunculus var. altilis)-Based Biorefinery: A Case Study of Improved Cell Cultures via Genetic Modulation of the Phenylpropanoid Pathway

Cultivated cardoon (Cynara cardunculus var. altilis L.) is a promising candidate species for the development of plant cell cultures suitable for large-scale biomass production and recovery of nutraceuticals. We set up a protocol for Agrobacterium tumefaciens-mediated transformation, which can be used for the improvement of cardoon cell cultures in a frame of biorefinery. As high lignin content determines lower saccharification yields for the biomass, we opted for a biotechnological approach, with the purpose of reducing lignin content; we generated transgenic lines overexpressing the Arabidopsis thaliana MYB4 transcription factor, a known repressor of lignin/flavonoid biosynthesis. Here, we report a comprehensive characterization, including metabolic and transcriptomic analyses of AtMYB4 overexpression cardoon lines, in comparison to wild type, underlining favorable traits for their use in biorefinery. Among these, the improved accessibility of the lignocellulosic biomass to degrading enzymes due to depletion of lignin content, the unexpected increased growth rates, and the valuable nutraceutical profiles, in particular for hydroxycinnamic/caffeoylquinic and fatty acids profiles.

A holistic approach to soil contamination and sustainable phytoremediation with energy crops in the Aegean Region of Turkey

The objective of this current review article is to evaluate the current knowledge of the contaminated soil in the study area based on reports and the results of previous experimental studies in the literature and to discuss the feasibility of phytoremediation with biofuel production using energy crops. The results indicated that the soil contamination was related mainly to the thermal power plant and mining activities in Kütahya, high industrial activity in İzmir, heavy metal and radioactive pollution in Manisa and Muğla. Moreover, the sources of the contamination are geothermal resources and transportation in Aydın and Denizli, respectively. However, soil pollution in Afyonkarahisar and Uşak provinces has not been discussed due to a lack of detailed reports and data in the literature. Besides, energy crops such as Zea mays, Ricinus communis, and Gossypium hirsitum were identified as appropriate candidates for İzmir, Denizli, Manisa, and Aydın due to being resistant to the arid climate. In Muğla province, Eucalyptus grandis and Eucalyptus bicostata can be cultivated because of having adaptation to moderate climatic conditions. Ricinus communis and Helianthus annuus were determined to be very suitable energy crops for the phytoremediation of many heavy metals in Kütahya. The review promotes the development of economic, environmental, and social benefits to regain the contaminated areas through phytoremediation. The findings of the study are important for creating sustainable solutions for remediation of polluted soils in Turkey, as well as for shedding light on the process of establishing appropriate policies to make soils contaminated suitable for agriculture.

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.

Active Packaging-Releasing System with Foeniculum vulgare Essential Oil for the Quality Preservation of Ready-to-Cook (RTC) Globe Artichoke Slices

Two globe artichoke genotypes, “Spinoso sardo” and “Opera F1”, have been processed as ready-to-cook (RTC) slices and refrigerated at 4 °C for 12 days (i) to evaluate the suitability to be processed as RTC slices; (ii) to evaluate the effect of a Foeniculum vulgare essential oil (EO) emitter, within an active package system, to delay quality decay, thus extending shelf life; (iii) to estimate the impact of EO emitter on the sensory profile of the RTC slices after cooking. Results revealed that both globe artichoke genotypes possess a good attitude to be processed as RTC product. “Opera F1” showed the best performances for color parameters, texture and chemical indexes, while “Spinoso sardo” showed lower mass loss (ML) over the storage time. The addition of EO emitter slowed down the consumption of O₂, better preserved texture when compared to the control and more effectively control polyphenol oxidase (PPO) activity and antioxidants’ retention during the cold storage. Microbial counts in control globe artichoke RTC slices were significantly higher than those packed with EO emitter, confirming the inhibiting role played by EO of F. vulgare. In addition, the EO emitter did not influence negatively the sensory profile of RTC globe artichoke slices after microwave cooking.

Genome-wide identification and analysis of DNA methyltransferase and demethylase gene families in Dendrobium officinale reveal their potential functions in polysaccharide accumulation

BACKGROUND

DNA methylation is a conserved and important epigenetic modification involved in the regulation of numerous biological processes, including plant development, secondary metabolism, and response to stresses. However, no information is available regarding the identification of cytosine-5 DNA methyltransferase (C5-MTase) and DNA demethylase (dMTase) genes in the orchid Dendrobium officinale.

RESULTS

In this study, we performed a genome-wide analysis of DoC5-MTase and DodMTase gene families in D. officinale. Integrated analysis of conserved motifs, gene structures and phylogenetic analysis showed that eight DoC5-MTases were divided into four subfamilies (DoCMT, DoDNMT, DoDRM, DoMET) while three DodMTases were divided into two subfamilies (DoDML3, DoROS1). Multiple cis-acting elements, especially stress-responsive and hormone-responsive ones, were found in the promoter region of DoC5-MTase and DodMTase genes. Furthermore, we investigated the expression profiles of DoC5-MTase and DodMTase in 10 different tissues, as well as their transcript abundance under abiotic stresses (cold and drought) and at the seedling stage, in protocorm-like bodies, shoots, and plantlets. Interestingly, most DoC5-MTases were downregulated whereas DodMTases were upregulated by cold stress. At the seedling stage, DoC5-MTase expression decreased as growth proceeded, but DodMTase expression increased.

CONCLUSIONS

These results provide a basis for elucidating the role of DoC5-MTase and DodMTase in secondary metabolite production and responses to abiotic stresses in D. officinale.

"Mind the Gap": Hi-C Technology Boosts Contiguity of the Globe Artichoke Genome in Low-Recombination Regions

Globe artichoke (Cynara cardunculus var. scolymus; 2n2x=34) is cropped largely in the Mediterranean region, being Italy the leading world producer; however, over time, its cultivation has spread to the Americas and China. In 2016, we released the first (v1.0) globe artichoke genome sequence (http://www.artichokegenome.unito.it/). Its assembly was generated using ∼133-fold Illumina sequencing data, covering 725 of the 1,084 Mb genome, of which 526 Mb (73%) were anchored to 17 chromosomal pseudomolecules. Based on v1.0 sequencing data, we generated a new genome assembly (v2.0), obtained from a Hi-C (Dovetail) genomic library, and which improves the scaffold N₅₀ from 126 kb to 44.8 Mb (∼356-fold increase) and N₉₀ from 29 kb to 17.8 Mb (∼685-fold increase). While the L₉₀ of the v1.0 sequence included 6,123 scaffolds, the new v2.0 just 15 super-scaffolds, a number close to the haploid chromosome number of the species. The newly generated super-scaffolds were assigned to pseudomolecules using reciprocal blast procedures. The cumulative size of unplaced scaffolds in v2.0 was reduced of 165 Mb, increasing to 94% the anchored genome sequence. The marked improvement is mainly attributable to the ability of the proximity ligation-based approach to deal with both heterochromatic (e.g.: peri-centromeric) and euchromatic regions during the assembly procedure, which allowed to physically locate low recombination regions. The new high-quality reference genome enhances the taxonomic breadth of the data available for comparative plant genomics and led to a new accurate gene prediction (28,632 genes), thus promoting the map-based cloning of economically important genes.

Hybrid transcriptome sequencing approach improved assembly and gene annotation in Cynara cardunculus (L.)

Background

The investigation of transcriptome profiles using short reads in non-model organisms, which lack of well-annotated genomes, is limited by partial gene reconstruction and isoform detection. In contrast, long-reads sequencing techniques revealed their potential to generate complete transcript assemblies even when a reference genome is lacking. Cynara cardunculus var. altilis (DC) (cultivated cardoon) is a perennial hardy crop adapted to dry environments with many industrial and nutraceutical applications due to the richness of secondary metabolites mostly produced in flower heads. The investigation of this species benefited from the recent release of a draft genome, but the transcriptome profile during the capitula formation still remains unexplored. In the present study we show a transcriptome analysis of vegetative and inflorescence organs of cultivated cardoon through a novel hybrid RNA-seq assembly approach utilizing both long and short RNA-seq reads.

Results

The inclusion of a single Nanopore flow-cell output in a hybrid sequencing approach determined an increase of 15% complete assembled genes and 18% transcript isoforms respect to short reads alone. Among 25,463 assembled unigenes, we identified 578 new genes and updated 13,039 gene models, 11,169 of which were alternatively spliced isoforms. During capitulum development, 3424 genes were differentially expressed and approximately two-thirds were identified as transcription factors including bHLH, MYB, NAC, C2H2 and MADS-box which were highly expressed especially after capitulum opening. We also show the expression dynamics of key genes involved in the production of valuable secondary metabolites of which capitulum is rich such as phenylpropanoids, flavonoids and sesquiterpene lactones. Most of their biosynthetic genes were strongly transcribed in the flower heads with alternative isoforms exhibiting differentially expression levels across the tissues.

Conclusions

This novel hybrid sequencing approach allowed to improve the transcriptome assembly, to update more than half of annotated genes and to identify many novel genes and different alternatively spliced isoforms. This study provides new insights on the flowering cycle in an Asteraceae plant, a valuable resource for plant biology and breeding in Cynara and an effective method for improving gene annotation.

A New Intra-Specific and High-Resolution Genetic Map of Eggplant Based on a RIL Population, and Location of QTLs Related to Plant Anthocyanin Pigmentation and Seed Vigour

Eggplant is the second most important solanaceous berry-producing crop after tomato. Despite mapping studies based on bi-parental progenies and GWAS approaches having been performed, an eggplant intraspecific high-resolution map is still lacking. We developed a RIL population from the intraspecific cross ‘305E40’, (androgenetic introgressed line carrying the locus Rfo-Sa1 conferring Fusarium resistance) x ‘67/3’ (breeding line whose genome sequence was recently released). One hundred and sixty-three RILs were genotyped by a genotype-by-sequencing (GBS) approach, which allowed us to identify 10,361 polymorphic sites. Overall, 267 Gb of sequencing data were generated and ~773 M Illumina paired end (PE) reads were mapped against the reference sequence. A new linkage map was developed, including 7249 SNPs assigned to the 12 chromosomes and spanning 2169.23 cM, with iaci@liberoan average distance of 0.4 cM between adjacent markers. This was used to elucidate the genetic bases of seven traits related to anthocyanin content in different organs recorded in three locations as well as seed vigor. Overall, from 7 to 17 QTLs (at least one major QTL) were identified for each trait. These results demonstrate that our newly developed map supplies valuable information for QTL fine mapping, candidate gene identification, and the development of molecular markers for marker assisted selection (MAS) of favorable alleles.

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

Sweet pepper (Capsicum annuum L.) is a high value crop and one of the most widely grown vegetables belonging to the Solanaceae family. In addition to commercial varieties and F₁ hybrids, a multitude of landraces are grown, whose genetic combination is the result of hundreds of years of random, environmental, and farmer selection. High genetic diversity exists in the landrace gene pool which however has scarcely been studied, thus bounding their cultivation. We re-sequenced four pepper inbred lines, within as many Italian landraces, which representative of as many fruit types: big sized blocky with sunken apex ('Quadrato') and protruding apex or heart shaped ('Cuneo'), elongated ('Corno') and smaller sized sub-spherical ('Tumaticot'). Each genomic sequence was obtained through Illumina platform at coverage ranging from 39 to 44×, and reconstructed at a chromosome scale. About 35.5k genes were predicted in each inbred line, of which 22,017 were shared among them and the reference genome (accession 'CM334'). Distinctive variations in miRNAs, resistance gene analogues (RGAs) and susceptibility genes (S-genes) were detected. A detailed survey of the SNP/Indels occurring in genes affecting fruit size, shape and quality identified the highest frequencies of variation in regulatory regions. Many structural variations were identified as presence/absence variations (PAVs), notably in resistance gene analogues (RGAs) and in the capsanthin/capsorubin synthase (CCS) gene. The large allelic diversity observed in the four inbred lines suggests their potential use as a pre-breeding resource and represents a one-stop resource for C. annuum genomics and a key tool for dissecting the path from sequence variation to phenotype.

Insights into sucrose pathway of chicory stems by integrative transcriptomic and metabolic analyses

The worldwide-cultivated chicory (Cichorium intybus L.) produces food and beneficial compounds, and young pre-flowering inflorescence stems are newly marketed vegetables. These sink-organs undergo growth by metabolizing sugars of leaf origin; the carbohydrate content and sweetness are crucial aspects for consumers' nutrition and acceptance. NMR profiling of 31 hydrosoluble phytochemicals showed that stem contents varied as influenced by genotype, environment and interaction, and that higher sucrose levels were associated with the sweeter of two landraces. Integrative analyses of metabolic and transcriptomic profile variations allowed the dissection of sucrose pathway. Overall, 427 and 23 unigenes respectively fell into the categories of sucrose metabolism and sugar carriers. Among 10 differentially expressed genes, the 11474/sucrose synthase, 53458/fructokinase, 9306 and 17035/hexokinases, and 20171/SWEET-type genes significantly associated to sugar content variation, and deduced proteins were characterised in silico. Correlation analyses encompassing sugar level variation, expressions of the former genes and of computationally assigned transcription factors (10938/NAC, 14712/bHLH, 40133/TALE and 17846/MIKC) revealed a gene network. The latter was minimally affected by the environment and accomplished with markers, representing a resource for biological studies and breeding.

Comparison of the Oxidative Stability and Antioxidant Activity of Extra-Virgin Olive Oil and Oils Extracted from Seeds of Colliguaya integerrima and Cynara cardunculus Under Normal Conditions and After Thermal Treatment

We investigated the potential of two oil extracts from seeds of Colliguaya integerrima (CIO) and Cynara cardunculus (CO) to use as nutritionally edible oils. For this purpose, oil quality was accessed by determining the fatty acid composition, peroxide value, acid value, iodine value, saponification number, phenolic contents, and oxidative stability during thermally induced oxidation of CIO and CO oils and compared to those of extra-virgin olive oil (EVOO). The chemical composition results demonstrated that both oils could be nutritional sources of essential unsaturated fatty acids. Moreover, according to the gravimetric analysis, the main decomposition step occurred in the temperature range of 200-420 °C, showing a similar thermal behavior of EVOO oil. However, CO and EVOO oils showed a higher phenolic content at degradation onset temperature (T0) in contrast with CIO oil. The antioxidant activity of the different studied oils showed a direct correlation with the phenol contents, up to temperatures around 180 °C, where the percentage of free radical scavenging assay for EVOO was higher than CO in contrast with the TPC values. Finally, we analyzed the minor components before and after heating CIO and CO at 180 °C by gas chromatography-mass spectrometry (GC-MS) using library search programs.

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.

RNA-Seq De Novo Assembly and Differential Transcriptome Analysis of Korean Medicinal Herb Cirsium japonicum var. spinossimum

Cirsium japonicum belongs to the Asteraceae or Compositae family and is a medicinal plant in Asia that has a variety of effects, including tumour inhibition, improved immunity with flavones, and antidiabetic and hepatoprotective effects. Silymarin is synthesized by 4-coumaroyl-CoA via both the flavonoid and phenylpropanoid pathways to produce the immediate precursors taxifolin and coniferyl alcohol. Then, the oxidative radicalization of taxifolin and coniferyl alcohol produces silymarin. We identified the expression of genes related to the synthesis of silymarin in C. japonicum in three different tissues, namely, flowers, leaves, and roots, through RNA sequencing. We obtained 51,133 unigenes from transcriptome sequencing by de novo assembly using Trinity v2.1.1, TransDecoder v2.0.1, and CD-HIT v4.6 software. The differentially expressed gene analysis revealed that the expression of genes related to the flavonoid pathway was higher in the flowers, whereas the phenylpropanoid pathway was more highly expressed in the roots. In this study, we established a global transcriptome dataset for C. japonicum. The data shall not only be useful to focus more deeply on the genes related to product medicinal metabolite including flavolignan but also to study the functional genomics for genetic engineering of C. japonicum.

Genotyping-by-sequencing highlights patterns of genetic structure and domestication in artichoke and cardoon

Exploiting the biodiversity of crops and their wild relatives is fundamental for maintaining and increasing food security. The species Cynara cardunculus includes three taxa: the globe artichoke, one of the most important Mediterranean vegetables, the leafy cardoon, and the wild cardoon. In this study, genotyping by sequencing (GBS) was successfully applied to reveal thousands of polymorphisms in a C. cardunculus germplasm collection, including 65 globe artichoke, 9 leafy cardoon, and 21 wild cardoon samples. The collection showed a strong population structure at K = 2, separating the globe artichoke from the leafy and wild cardoon. At higher K values, further substructures were observed, in which the wild cardoon was separated from the leafy cardoon, and the latter included the Spanish wild cardoons, while the wild sample from Portugal was admixed. Moreover, subpopulations within the globe artichoke set were highlighted. Structure analysis restricted to the globe artichoke dataset pointed out genetic differentiation between the ˝Catanesi˝ typology and all the other samples (K = 2). At higher values of K, the separation of the ˝Catanesi˝ group still held true, and green headed landraces from Apulia region, Italy (˝Green Apulian˝) formed a distinct subpopulation. ˝Romaneschi˝ artichoke types fell in a variable group with admixed samples, indicating that they should not be considered as a genetically uniform typology. The results of principal component analysis and Neighbor-Joining hierarchical clustering were consistent with structure results, and in addition provided a measure of genetic relationships among individual genotypes. Both analyses attributed the wild material from Spain and Portugal to the cultivated cardoon group, supporting the idea that this might be indeed a feral form of the leafy cardoon. Different reproductive habit and possibly selective pressure led to a slower LD decay in artichoke compared to cardoon. Genotyping by sequencing has proven a reliable methodology to obtain valuable SNPs and assess population genetics in C. cardunculus.