Search

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 F1 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.

Significance This study provides a deep population genomic analysis of 10,000 Capsicum accessions held in genebanks and representing a frame of the global diversity of the genus. By combining single nucleotide polymorphisms (SNPs) based data and passport information, we investigated the genomic diversity and population structure of wild and domesticated peppers, tracing back to routes of evolution and providing a model of Capsicum annuum distribution, which reflects human trade and historical/cultural influences. Our results highlight west–east routes of expansion, shedding light on the links between South and Mesoamerica, Africa, and East/South Asia, the latter two constituting important diversification centers of pepper diversity. Finally, we outline a roadmap for genebank management and future direction for better exploitation of germplasm resources., Genebanks collect and preserve vast collections of plants and detailed passport information, with the aim of preserving genetic diversity for conservation and breeding. Genetic characterization of such collections has the potential to elucidate the genetic histories of important crops, use marker–trait associations to identify loci controlling traits of interest, search for loci undergoing selection, and contribute to genebank management by identifying taxonomic misassignments and duplicates. We conducted a genomic scan with genotyping by sequencing (GBS) derived single nucleotide polymorphisms (SNPs) of 10,038 pepper (Capsicum spp.) accessions from worldwide genebanks and investigated the recent history of this iconic staple. Genomic data detected up to 1,618 duplicate accessions within and between genebanks and showed that taxonomic ambiguity and misclassification often involve interspecific hybrids that are difficult to classify morphologically. We deeply interrogated the genetic diversity of the commonly consumed Capsicum annuum to investigate its history, finding that the kinds of peppers collected in broad regions across the globe overlap considerably. The method ReMIXTURE—using genetic data to quantify the similarity between the complement of peppers from a focal region and those from other regions—was developed to supplement traditional population genetic analyses. The results reflect a vision of pepper as a highly desirable and tradable cultural commodity, spreading rapidly throughout the globe along major maritime and terrestrial trade routes. Marker associations and possible selective sweeps affecting traits such as pungency were observed, and these traits were shown to be distributed nonuniformly across the globe, suggesting that human preferences exerted a primary influence over domesticated pepper genetic structure.

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.

In a search for natural compounds with anti-HIV-1 activity, we studied the effect of the ethanolic extract obtained from leaves, bark, and peels of Punica granatum L. for the inhibition of the HIV-1 reverse transcriptase (RT)-associated ribonuclease H (RNase H) and integrase (IN) LEDGF-dependent activities. The chemical analyses led to the detection of compounds belonging mainly to the phenolic and flavonoid chemical classes. Ellagic acid, flavones, and triterpenoid molecules were identified in leaves. The bark and peels were characterized by the presence of hydrolyzable tannins, such as punicalins and punicalagins, together with ellagic acid. Among the isolated compounds, the hydrolyzable tannins and ellagic acid showed a very high inhibition (IC50 values ranging from 0.12 to 1.4 µM and 0.065 to 0.09 µM of the RNase H and IN activities, respectively). Of the flavonoids, luteolin and apigenin were found to be able to inhibit RNase H and IN functions (IC50 values in the 3.7–22 μM range), whereas luteolin 7-O-glucoside showed selective activity for HIV-1 IN. In contrast, betulinic acid, ursolic acid, and oleanolic acid were selective for the HIV-1 RNase H activity. Our results strongly support the potential of non-edible P. granatum organs as a valuable source of anti-HIV-1 compounds.

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.

[EN] 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-1 and T-2 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., Research was financially supported by the project CRISPR/Cas9-mediated gene knock-out in eggplant financed by Compagnia San Paolo.

Zika virus, an arthropod-borne flavivirus, is an emerging healthcare threat worldwide. Zika virus is responsible for severe neurological effects, such as paralytic Guillain-Barrè syndrome, in adults, and also congenital malformations, especially microcephaly. No specific antiviral drugs and vaccines are currently available, and treatments are palliative, but medicinal plants show great potential as natural sources of anti-Zika phytochemicals. This study deals with the investigation of the composition, cytotoxicity, and anti-Zika activity of Punica granatum leaf ethanolic extract, fractions, and phytoconstituents. P. granatum leaves were collected from different areas in Italy and Greece in different seasons. Crude extracts were analyzed and fractionated, and the pure compounds were isolated. The phytochemical and biomolecular fingerprint of the pomegranate leaves was determined. The antiviral activities of the leaf extract, fractions, and compounds were investigated against the MR766 and HPF2013 Zika virus strains in vitro. Both the extract and its fractions were found to be active against Zika virus infection. Of the compounds isolated, ellagic acid showed particular anti-Zika activities, with EC50 values of 30.86 µM for MR766 and 46.23 µM for HPF2013. The mechanism of action was investigated using specific antiviral assays, and it was demonstrated that ellagic acid was primarily active as it prevented Zika virus infection and was able to significantly reduce Zika virus progeny production. Our data demonstrate the anti-Zika activity of pomegranate leaf extract and ellagic acid for the first time. These findings identify ellagic acid as a possible anti-Zika candidate compound that can be used for preventive and therapeutic interventions.

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 N50 from 126 kb to 44.8 Mb (∼356-fold increase) and N90 from 29 kb to 17.8 Mb (∼685-fold increase). While the L90 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.

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 &lsquo, Bouche de B&eacute, tizac&rsquo, (C. sativa &times, C. crenata) and studied by developing genetic linkage maps using a population derived from a cross between &lsquo, and the susceptible cultivar &lsquo, Madonna&rsquo, (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 &lsquo, consisted of 1459 loci and spanned 809.6 cM, the map of &lsquo, 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 &lsquo, map, explaining up to 67&ndash, 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.

Hazelnut breeding is a slow activity due to the long lifecycle of the plant, the long time between pollination and reliable phenotype observation in the progeny, the presence of self-incompatibility and the high level of heterozygosity. Improved knowledge on the hazelnut genome and the development of marker-assisted selection would greatly facilitate breeding programs. Objective of this study was the construction of high-density genetic maps and the detection of QTL (quantitative trait loci) related to phenological, vegetative and productive traits, by using an F1 progeny obtained by crossing ‘Tonda Gentile delle Langhe’ (female parent, hereafter TGdL) with ‘Merveille de Bollwiller’ (syn. ‘Hall’s Giant’ male parent, MB). Genotyping-by-sequencing (GBS), followed by SNP mining using a Samtools-based pipeline, yielded 7,980 SNP markers. Independent linkage maps were constructed for each parent on the basis of the double pseudo-testcross mapping strategy. The linkage map for TGdL consisted of 1,236 markers covering a total genetic length of 900.4 cM, with a mean inter-marker distance of 0.83 cM. A reciprocal translocation was detected in TGdL between two non-homologous chromosomes. The linkage map for MB was based on segregation at 1,211 markers covering a total genetic length of 899.1 cM, with a mean inter-marker distance of 0.82 cM. Plant and nut traits were recorded during five years and statistically analyzed with molecular data to identify marker/QTL associations for flowering and bud burst time, sucker habit and nut traits, such as nut and seed size, roundness index, percent kernel and ease of pellicle removal.

Globe artichoke (Cynara cardunculus var. scolymus) is a cross-pollinated, highly heterozygous species, which is conventionally propagated vegetatively. A scheme is described here which combines phenotypic with genotypic selection to fast track the development of a seed-propagated variety. The scheme was tested by making three selections, on a phenotypic basis, from a Brazilian seed-propagated variety showing an high phenotypic variation. The genetic relatedness as well as the heterozygosity of the material in study, in respect to standard variety representatives, was initially assessed with a wide set of microsatellite markers. Afterwards, an AFLP-based selection demonstrated to provide a practical and cheap means of conducting marker assisted breeding, which can be easily adopted also in laboratories of small seed companies. The selection approach described here could be readily adopted also to convert current vegetatively propagated landraces into seed-propagated varieties.

[EN] 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., This work has been funded by the European Union's Horizon 2020 Research and Innovation Programme under the grant agreement number 677379 (G2P-SOL project: Linking genetic resources, genomes, and phenotypes of solanaceous crops).

BACKGROUND: Interferon-induced expression of programmed cell death ligands (PD-L1/PD-L2) may sustain tumour immuneevasion. Patients featuring MET amplification, a genetic lesion driving transformation, may benefit from anti-MET treatment. We explored if MET-targeted therapy interferes with Interferon-gamma modulation of PD-L1/PD-L2 in MET-amplified tumours.METHODS: PD-L1/PD-L2 expression and signalling pathways downstream of MET or Interferon-gamma were analysed in MET-amplified tumour cell lines and in patient-derived tumour organoids, in basal condition, upon Interferon-gamma stimulation, and after anti-MET therapy.RESULTS: PD-L1 and PD-L2 were upregulated in MET-amplified tumour cells upon Interferon-gamma treatment. This induction was impaired by JNJ-605, a selective inhibitor of MET kinase activity, and MvDN30, an antibody inducing MET proteolytic cleavage. We found that activation of JAKs/STAT1, signal transducers downstream of the Interferon-gamma receptor, was neutralised by MET inhibitors. Moreover, JAK2 and MET associated in the same signalling complex depending on MET phosphorylation. Results were confirmed in MET-amplified organoids derived from human colorectal tumours, where JNJ-605 treatment revoked Interferon-gamma induced PD-L1 expression.CONCLUSIONS: These data show that in MET-amplified cancers, treatment with MET inhibitors counteracts the induction of PD-1 ligands by Interferon-gamma. Thus, therapeutic use of anti-MET drugs may provide additional clinical benefit over and above the intended inhibition of the target oncogene.

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.

Pepper (Capsicum spp.) belongs to the Solanaceae, which is an economically important family of flowering plants consisting of about 102 genera and over 2500 species. The Solanaceae family includes crops of agronomic importance for which the efforts in genome sequencing are ongoing by almost 10 years (https://www.solgenomics.net/organism/sol100/view). Since the beginning of 2014, various consortia have released the genome sequences of domesticated and wild Capsicum species. The first effort was focused on the whole-genome sequencing of Capsicum annuum CM334 and of Capsicum chinense PI159236, which were widely used as founders of mapping populations and carry important disease resistance traits. Just a couple of months later, the genome sequences of C. annuum Zunla-1 and of the wild species Chiltepin (C. annuum var. glabriusculum) were published. Both studies reported a pepper genome size of ~3–3.5 Gb, rich in repetitive elements (over 80%) with about 35 thousand genes. The improved version of the reference genome CM334 as well as of C. chinense PI159236 together with the sequencing of the domesticated Capsicum baccatum revealed evolutionary relationships and estimated lineage divergence times occurring in Capsicum. Recently, the linked-read sequencing technology has been applied for the sequencing of a C. annuum accession that was an F1 cross hybrid of CM334 and a non-pungent pepper breeding line. Furthermore, genome resequencing studies have been performed with the aim to analyze loci of interest related to biotic/abiotic stresses and to qualitative features. In this chapter, we provide an overview of the genome sequencing and annotation strategies and describe the main results disclosed by all the whole and targeted genome sequencing projects in Capsicum.

The genome sequence of globe artichoke (Cynara cardunculus L. var. scolymus, 2n = 2x = 34) is now available for use. A survey of C. cardunculus genetic resources is essential for understanding the evolution of the species, carrying out genetic studies and for application of breeding strategies. We report on the resequencing analyses (~35×) of four globe artichoke genotypes, representative of the core varietal types, as well as a genotype of the related taxa cultivated cardoon. The genomes were reconstructed at a chromosomal scale and structurally/functionally annotated. Gene prediction indicated a similar number of genes, while distinctive variations in miRNAs and resistance gene analogues (RGAs) were detected. Overall, 23,5 M SNP/indel were discovered (range 6,34 M –14,50 M). The impact of some missense SNPs on the biological functions of genes involved in the biosynthesis of phenylpropanoid and sesquiterpene lactone secondary metabolites was predicted. The identified variants contribute to infer on globe artichoke domestication of the different varietal types, and represent key tools for dissecting the path from sequence variation to phenotype. The new genomic sequences are fully searchable through independent Jbrowse interfaces (www.artichokegenome.unito.it), which allow the analysis of collinearity and the discovery of genomic variants, thus representing a one-stop resource for C. cardunculus genomics.

Caleosins are involved in several cellular and biological processes that are closely associated with the synthesis, degradation, and stability of oil bodies (OBs). 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 OBs has been 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 to verify only the presence of CavCLO-H1 in hazelnut OBs, 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 peroxygenase.

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 Betizac’ (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.

The cultivation of nut species, such as hazelnut and chestnut, still largely relies on traditional cultivars that were selected across centuries. Up to date, only a few breeding programs are active, yielding a limited number of new accessions, due to the particularly long time required for the selection process and the high heterozygosity of these species. There is an urgent need for improved breeding techniques based on marker-assisted selection (MAS) and biotechnological tools, to face new threats, such as the spread of alien species and the climatic changes. In this research we present the results of 10 years of activity addressed to the development of genetic maps in hazelnut (Corylus avellana L., progeny 'Tonda Gentile delle Langhe' × 'Merveille de Bollwiller' syn. 'Hall's Giant') and chestnut (progeny Castanea sativa × (C. sativa × C. crenata)) aimed at the detection of markers putatively useful for MAS. Saturated maps for the parent cultivars of the two progenies, accounting about 180 individuals each, were constructed using both SSR markers and SNP markers obtained by genotyping by sequencing, using the reference genomes as template. In hazelnut, phenological traits and fruit characters of commercial and technological relevance were mapped as QTLs; in chestnut, a trait of resistance to Dryocosmus kuriphilus Yasumatsu was mapped as a simple Mendelian character, and its genomic region is under deep investigation to find putative genes of resistance. Further breeding strategies are being carried out for hazelnut using in vitro techniques aimed at developing efficient systems of genetic transformation and regeneration for the modification or introduction of genes, such as those related to self-incompatibility and resistance to major pathogens.

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.

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.

Methods for investigating DNA methylation nowadays either require a reference genome and high coverage, or investigate only CG methylation. Moreover, no large-scale analysis can be performed for N6-methyladenosine (6mA). Here we describe the methylation content sensitive enzyme double-digest restriction-site-associated DNA (ddRAD) technique (MCSeEd), a reduced-representation, reference-free, cost-effective approach for characterizing whole genome methylation patterns across different methylation contexts (e.g., CG, CHG, CHH, 6mA). MCSeEd can also detect genetic variations among hundreds of samples. MCSeEd is based on parallel restrictions carried out by combinations of methylation insensitive and sensitive endonucleases, followed by next-generation sequencing. Moreover, we present a robust bioinformatic pipeline (available at https://bitbucket.org/capemaster/mcseed/src/master/) for differential methylation analysis combined with single nucleotide polymorphism calling without or with a reference genome.

Globe artichoke represents one of the main horticultural species of the Mediterranean basin, and &lsquo, Spinoso sardo&rsquo, 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 &lsquo, 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 &lsquo, non-conventional phenotype.

Background: DNA copy number alterations (CNAs) and gene expression changes have amply been encountered in colorectal cancers (CRCs), but the extent at which CNAs affect gene expression, as well as their relevance for tumor development, are still poorly defined. Here we aimed at assessing the clinical relevance of these parameters in a 10 year follow-up study. Methods: Tumors and normal adjacent colon mucosa, obtained at primary surgery from 21 CRC patients, were subjected to (i) high-resolution array CGH (a-CGH) for the detection of CNAs and (ii) microarray-based transcriptome profiling for the detection of gene expression (GE) changes. Correlations between these genomic and transcriptomic changes and their associations with clinical and histopathological parameters were assessed with the aim to identify molecular signatures associated with disease-free survival of the CRC patients during a 10 year follow-up. Results: DNA copy number gains were frequently detected in chromosomes 7, 8q, 13, 19, 20q and X, whereas DNA copy number losses were frequently detected in chromosomes 1p, 4, 8p, 15, 17p, 18, 19 and 22q. None of these alterations were observed in all samples. In addition, we found that 2,498 genes were up- and that 1,094 genes were down-regulated in the tumor samples compared to their corresponding normal mucosa (p < 0.01). The expression of 65 genes was found to be significantly associated with prognosis (p < 0.01). Specifically, we found that up-regulation of the IL17RA, IGF2BP2 and ABCC2 genes, and of genes acting in the mTOR and cytokine receptor pathways, were strongly associated with a poor survival. Subsequent integrated analyses revealed that increased expression levels of the MMP9, BMP7, UBE2C, I-CAM, NOTCH3, NOTCH1, PTGES2, HMGB1 and ERBB3 genes were associated with copy number gains, whereas decreased expression levels of the MUC1, E2F2, HRAS and SIRT3 genes were associated with copy number losses. Pathways related to cell cycle progression, eicosanoid metabolism, and TGF-? and apoptosis signaling, were found to be most significantly affected. Conclusions: Our results suggest that CNAs in CRC tumor tissues are associated with concomitant changes in the expression of cancer-related genes. In other genes epigenetic mechanism may be at work. Up-regulation of the IL17RA, IGF2BP2 and ABCC2 genes, and of genes acting in the mTOR and cytokine receptor pathways, appear to be associated with a poor survival. These alterations may, in addition to Dukes' staging, be employed as new prognostic biomarkers for the prediction of clinical outcome in CRC patients.

The first high-quality genome assembly of the globe artichoke has been produced within the Compositae Genome Project and the resequencing analyses of four globe artichoke genotypes, representative of the core varietal types, as well as a genotype of the related taxa cultivated cardoon was, later on, carried out. The Web site “www.artichokegenome.unito.it” hosts all the available genomic sequences, together with their structural/functional annotations and project information are presented to users via the open-source tool JBrowse, allowing the analysis of collinearity and the discovery of genomic variants, thus representing a one-stop resource for Cynara cardunculus genomics. Pseudomolecules as well as unmapped scaffolds were used for the bulk mining of SSR markers and for the construction of the first globe artichoke microsatellite marker database. A database, called “Cynara cardunculus MicroSatellite DataBase” (CyMSatDB) was developed to provide a searchable interface to the SSR data. CyMSatDB facilitates the retrieval of SSR markers, as well as suggested forward and reverse primers, on the basis of genomic location, genomic versus genic context, perfect versus imperfect repeat, motif type, motif sequence, and repeat number.

Genetic variation within globe artichoke landraces contributes to the variable yield performance detected in cultivation. It was built up over many generations of vegetative propagation as a consequence of the limited selection applied by farmers as well as the accumulation of mutations. This chapter reports the results obtained by applying clonal selection programmes within local landraces, with the goal to identify genotypes retaining desirable traits and able to ensure a higher and stable yield. The selected high-yielding genotypes while contributing to increase productivity may also implement the conservation ‘in situ’ of valuable germplasm at risk of extinction. With the goal to move to breeding programme based on crossing strategies, it is pivotal the development of genetic maps with high levels of genome coverage. Genetic maps represent the first step for localizing genes or quantitative trait loci (QTL) that are associated with economically important trait and open the way for marker-assisted selection, comparative mapping between different species, a framework for anchoring physical maps and the basis for map-based cloning of genes. In this chapter, the history of the development of genetic maps in Cynara cardunculus is retraced, starting from the first maps based on a limited number of mostly dominant molecular markers, up to the most recently developed maps allowing the identification of genomic region influencing key agronomic traits.

A correct botanical identification and analytical quality control of volatile key-markers responsible for aroma and biological activities is necessary to monitor volatile compounds transferred from a plant to the related herbal tea and human intake to guarantee their safe use. This is mainly true for markers limited by regulations or by a recommended maximum amount of consumption per day. GC–MS is the elective technique to analyze volatiles, provided that for aqueous samples (herbal teas) an appropriate sample preparation procedure, and/or a water-compatible GC stationary phases are applied. Solid Phase Micro Extraction (SPME) on-line coupled to GC–MS in a fully automatic approach is here applied to sample and quantify key markers in plant material (headspace) and in the corresponding herbal tea (direct immersion). In parallel, a new generation of GC columns coated with ionic liquid based stationary phases compatible with aqueous samples (Watercol™) was applied to test direct injection of aqueous samples (DAI-GC-FID). The latter approach fully bypasses sample preparation thus speeding up quality control. This study deals with the quantitation of menthol, α- and β-thujone, estragole, and anethole contained in several plant species commonly used for herbal teas (i.e. peppermint, sage, wormwood, fennel, aniseed) and regulated by International Organizations. The two methods gave comparable results and are characterized by high repeatability, linearity and accuracy, although, as expected, their sensitivity was different because DAI-GC-FID implies injection of the sample as such without analyte concentration as for DI-SPME-GC–MS. For instance, LOD and LOQ of estragole were 0.03 and 0.1 mg L−1 with DI-SPME-GC–MS and 0.1 and 0.8 mg L−1 with DAI-GC-FID. The two methods are fully complementary and their adoption depends on the amount of marker(s) to be quantified.

Understanding the population genetic structure of crop gene pool is a fundamental prerequisite in plant evolutionary biology and plant breeding. Based on this consideration, in this chapter we pursued two main aims. The first was to review the status of the research pointed at deciphering the population structure of the gene pools of globe artichoke. The second was to show how the use of a well-characterized germplasm collection together with deep genotyping could serve the identification of QTLs for traits of agronomic and economic relevance in globe artichoke.

Methods for investigating DNA methylation nowadays either require a reference genome and high coverage, or investigate only CG methylation. Moreover, no large-scale analysis can be performed for N6-methyladenosine (6 mA) at an affordable price. Here we describe the methylation content sensitive enzyme double-digest restriction-site-associated DNA (ddRAD) technique (MCSeEd), a reduced-representation, reference-free, cost-effective approach for characterizing whole genome methylation patterns across different methylation contexts (e.g., CG, CHG, CHH, 6 mA). MCSeEd can also detect genetic variations among hundreds of samples. MCSeEd is based on parallel restrictions carried out by combinations of methylation insensitive and sensitive endonucleases, followed by next-generation sequencing. Moreover, we present a robust bioinformatic pipeline (available at https://bitbucket.org/capemaster/mcseed/src/master/) for differential methylation analysis combined with single nucleotide polymorphism calling without or with a reference genome.

The bracts of the globe artichoke inflorescence vary in their pigmentation and the extent of fleshy thorn development. Here, a genetic analysis of these two traits is presented, based on a pre-existing and well-characterized mapping population derived from a cross between a globe artichoke variety and a cultivated cardoon. While both traits appeared to be simply inherited when the number of trait classes was limited to three (pigmentation) or two (thorniness), extending the classification to a larger number of states generated continuous distributions appropriate for a quantitative trait locus (QTL) analysis. The seven QTL identified by this means, mapping to five different linkage groups, were all expressed across two growing seasons. The largest effect QTL mapped onto homologues linkage group in each parental maps, and explained up to 69 and 73 % of the phenotypic variance respectively for bracts pigmentation and thorniness.

We have identified a new t(1;21)(p32;q22) chromosomal translocation in a MDS/AML patient that results in expression of an aberrant C-terminally truncated RUNX1 protein lacking several regulatory domains. As similar truncated RUNX1 proteins are generated by genetic aberrations including chromosomal translocations and point mutations, we used the t(1;21)(p32;q22) chromosomal translocation as a model to explore whether C-terminally truncated RUNX1 proteins trigger effects similar to those induced by well-characterized leukemogenic RUNX1 fusion genes. In vitro analysis of transduced human hematopoietic/progenitor stem cells showed that truncated RUNX1 proteins increase proliferation and self-renewal and disrupt the differentiation program by interfering with RUNX1b. These effects are similar to but milder than those induced by the RUNX1/ETO fusion protein. GSEA analysis confirmed similar altered gene expression patterns in the truncated RUNX1 and RUNX1/ETO models, with both models showing alterations in genes involved in self-renewal and leukemogenesis, including homeobox genes, primitive erythroid genes and leukemogenic transcription factors. We propose that C-terminally truncated RUNX1 proteins can contribute to leukemogenesis in a similar way to RUNX1 fusion genes but through a milder phenotype.

Cynara cardunculus, a member of the Asteraceae family, comprises the three taxa var. scolymus (globe artichoke), var. altilis (cultivated cardoon) and the ancestral var. sylvestris (wild cardoon). The substantial quantities of lignocellulosic biomass produced by these plants (up to 30.0 t/ha year−1 dry matter by the cultivated cardoon) can be used either as a source of bioenergy and/or as raw material for paper pulp production. Here, genotyping-by-sequencing to an F1 population derived from a cross between a globe artichoke (C3) and a cultivated cardoon (ALT) genotypes has been used to perform a genome-wide linkage analysis, leading to the elaboration of a pair of highly dense genetic maps, each derived from one of the two highly heterozygous parental genotypes. In both maps, the number of linkage groups (17) matched the species’ haploid chromosome number. The F1 population was phenotyped over two seasons with respect to plant height, stem number, capitulum number, leaf and stem fresh weight, and the dry weight of the whole plant, the leaves, the stems, the capitula and the achenes. The phenotypic data were combined with the linkage maps to identify 81 quantitative trait loci, of which 50 were placed on the C3 map and 31 on the ALT map. The loci were scattered over 13 linkage groups, and were clustered within 27 genomic regions, 22 of which harboured two or more QTL. Ten of these regions were specific to the C3 map and six to the ALT map, while the other 11 were represented on both maps. The 27 regions harboured in all 1960 genes, 83% of which could be functionally annotated. An enrichment for certain gene ontology terms was noted for the gene content of the genomic regions harbouring loci influencing seed yield and the number/weight of stems.

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 F1 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.

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

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., Scientific Reports, 7 (1), ISSN:2045-2322

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.

Acute myeloid leukemia (AML) and myelodysplastic syndromes (MDS) represent a group of heterogeneous diseases with a considerable diversity in molecular pathogenesis and clinical outcomes. Recurrent...

Globe artichoke (Cynara cardunculus var. scolymus), whose immature inflorescences (capitula) are consumed as a vegetable all over the world, contributes significantly to the agricultural economy of the Mediterranean basin. Here, we describe a QTL (quantitative trait loci) analysis aimed at elucidating the mode of inheritance of seven main and first-order capitulum traits. Mapping was carried out in an F1 population obtained by crossing a globe artichoke with a cultivated cardoon (C. cardunculus var. altilis). A total of 100 QTL associated with the seven capitulum traits were mapped to 23 chromosomal regions, scattered over 12 of the 17 linkage groups. Among these, 73 were expressed in both growing seasons, while the others were only detected in one season. Up to nine QTL per trait were identified, and major QTL, responsible for some 20 % of the phenotypic variation, were detected for capitulum length, diameter, shape index and fresh weight. The QTL for correlated traits frequently co-localized, most likely due to pleiotropy. This study represents the first report on yield traits QTL in globe artichoke. The QTL identified, along with linked markers, particularly those located in four hot-spot QTL regions are of practical interest for crop improvement based on marker-assisted breeding.

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a very rare disease that currently lacks genomic and genetic biomarkers to assist in its clinical management. We performed whole-exome sequencing (WES) of three BPDCN cases. Based on these data, we designed a resequencing approach to identify mutations in 38 selected genes in 25 BPDCN samples. WES revealed 37-99 deleterious gene mutations per exome with no common affected genes between patients, but with clear overlap in terms of molecular and disease pathways (hematological and dermatological disease). We identified for the first time deleterious mutations in IKZF3, HOXB9, UBE2G2 and ZEB2 in human leukemia. Target sequencing identified 29 recurring genes, ranging in prevalence from 36% for previously known genes, such as TET2, to 12-16% for newly identified genes, such as IKZF3 or ZEB2. Half of the tumors had mutations affecting either the DNA methylation or chromatin remodeling pathways. The clinical analysis revealed that patients with mutations in DNA methylation pathway had a significantly reduced overall survival (P=0.047). We provide the first mutational profiling of BPDCN. The data support the current WHO classification of the disease as a myeloid disorder and provide a biological rationale for the incorporation of epigenetic therapies for its treatment.

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 (VIGS) 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.

DEK is an oncoprotein involved in a variety of cellular functions, such as DNA repair, replication, and transcriptional control. DEK is preferentially expressed in actively proliferating and malignant cells, including melanoma cell lines in which DEK was previously demonstrated to play a critical role in proliferation and chemoresistance. Still, the impact of this protein in melanoma progression remains unclear. Thus, we performed a comprehensive analysis of DEK expression in different melanocytic tumors. The immunostaining results of 303 tumors demonstrated negligible DEK expression in benign lesions. Conversely, malignant lesions, particularly in metastatic cases, were largely positive for DEK expression, which was partially associated with genomic amplification. Importantly, DEK overexpression was correlated with histological features of aggressiveness in primary tumors and poor prognosis in melanoma patients. In conclusion, our study provides new insight into the involvement of DEK in melanoma progression, as well as proof of concept for its potential application as a marker and therapeutic target of melanoma.

Globe artichoke (Cynara cardunculus var. scolymus) is an out-crossing, perennial, multi-use crop species that is grown worldwide and belongs to the Compositae, one of the most successful Angiosperm families. We describe the first genome sequence of globe artichoke. The assembly, comprising of 13,588 scaffolds covering 725 of the 1,084 Mb genome, was generated using ~133-fold Illumina sequencing data and encodes 26,889 predicted genes. Re-sequencing (30×) of globe artichoke and cultivated cardoon (C. cardunculus var. altilis) parental genotypes and low-coverage (0.5 to 1×) genotyping-by-sequencing of 163 F1 individuals resulted in 73% of the assembled genome being anchored in 2,178 genetic bins ordered along 17 chromosomal pseudomolecules. This was achieved using a novel pipeline, SOILoCo (Scaffold Ordering by Imputation with Low Coverage), to detect heterozygous regions and assign parental haplotypes with low sequencing read depth and of unknown phase. SOILoCo provides a powerful tool for de novo genome analysis of outcrossing species. Our data will enable genome-scale analyses of evolutionary processes among crops, weeds, and wild species within and beyond the Compositae, and will facilitate the identification of economically important genes from related species.

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.