Sergio Lanteri

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

BackgroundThe eggplant (Solanum melongena L.) genome is relatively unexplored, especially compared to those of the other major Solanaceae crops tomato and potato. In particular, no SNP markers are publicly available; on the other hand, over 1,000 SSR markers were developed and publicly available. We have combined the recently developed Restriction-site Associated DNA (RAD) approach with Illumina DNA sequencing for rapid and mass discovery of both SNP and SSR markers for eggplant.ResultsRAD tags were generated from the genomic DNA of a pair of eggplant mapping parents, and sequenced to produce ~17.5 Mb of sequences arrangeable into ~78,000 contigs. The resulting non-redundant genomic sequence dataset consisted of ~45,000 sequences, of which ~29% were putative coding sequences and ~70% were in common between the mapping parents. The shared sequences allowed the discovery of ~10,000 SNPs and nearly 1,000 indels, equivalent to a SNP frequency of 0.8 per Kb and an indel frequency of 0.07 per Kb. Over 2,000 of the SNPs are likely to be mappable via the Illumina GoldenGate assay. A subset of 384 SNPs was used to successfully fingerprint a panel of eggplant germplasm, producing a set of informative diversity data. The RAD sequences also included nearly 2,000 putative SSRs, and primer pairs were designed to amplify 1,155 loci.ConclusionThe high throughput sequencing of the RAD tags allowed the discovery of a large number of DNA markers, which will prove useful for extending our current knowledge of the genome organization of eggplant, for assisting in marker-aided selection and for carrying out comparative genomic analyses within the Solanaceae family.

Gene-based microsatellite development for mapping and phylogeny studies in eggplant.

BackgroundEggplant (Solanum melongena L.) is a member of the Solanaceae family. In spite of its widespread cultivation and nutritional and economic importance, its genome has not as yet been extensively investigated. Few analyses have been carried out to determine the genetic diversity of eggplant at the DNA level, and linkage relationships have not been well characterised. As for the other Solanaceae crop species (potato, tomato and pepper), the level of intra-specific polymorphism appears to be rather limited, and so it is important that an effort is made to develop more informative DNA markers to make progress in understanding the genetics of eggplant and to advance its breeding. The aim of the present work was to develop a set of functional microsatellite (SSR) markers, via an in silico analysis of publicly available DNA sequence.ResultsFrom >3,300 genic DNA sequences, 50 SSR-containing candidates suitable for primer design were recovered. Of these, 39 were functional, and were then applied to a panel of 44 accessions, of which 38 were cultivated eggplant varieties, and six were from related Solanum species. The usefulness of the SSR assays for diversity analysis and taxonomic discrimination was demonstrated by constructing a phylogeny based on SSR polymorphisms, and by the demonstration that most were also functional when tested with template from tomato, pepper and potato. As a results of BLASTN analyses, several eggplant SSRs were found to have homologous counterparts in the phylogenetically related species, which carry microsatellite motifs in the same position.ConclusionThe set of eggplant EST-SSR markers was informative for phylogenetic analysis and genetic mapping. Since EST-SSRs lie within expressed sequence, they have the potential to serve as perfect markers for genes determining variation in phenotype. Their high level of transferability to other Solanaceae species can be used to provide anchoring points for the integration of genetic maps across species.

Amplified fragment length polymorphism for genetic diversity assessment in globe artichoke

Globe artichoke (Cynara cardunculus L. var. scolymus L.) is a diploid (2n=2x=34), predominantly cross-pollinated plant native to the Mediterranean basin, and Italy contains the richest primary cultivated ‘gene pool’. Commercial production is mainly based on perennial cultivation of vegetatively propagated clones that are highly heterozygous and segregate widely when progeny-tested. Analysis of the artichoke genome by means of molecular markers has been limited to a few studies; here we report on the genetic relatedness among 118 artichoke accessions, including clones belonging to the same varietal type, two accessions of cultivated cardoon (C. cardunculus L. var. altilis DC.) and four accessions of wild cardoon [C. cardunculus L. var. sylvestris (Lamk) Fiori] as measured by amplified fragment length polymorphism (AFLP). Eight primer combinations yielded a total of 667 bands, of which 519 were polymorphic. Genetic similarities among accessions were calculated according to Jaccard’s Similarity Index and used to construct a dendrogram based on the unweighted pair group method using arithmetic averages. Our results demonstrate that AFLP markers can be useful in evaluating Cynara cardunculus genetic diversity and in classifying accessions to phylogenetic groups based on their genetic similarity values. Genetic variation among artichoke clones belonging to the same varietal type was in some cases higher than that found among accessions differently named and coming from different areas. The lowest Jaccard’s Similarity Index found within a varietal type can be considered as a threshold for the identification of accessions which share an analogous genetic background. This will enable the selection of representatives in order to develop and manage a germplasm ‘core collection’ as well as the identification of suitable material for future artichoke breeding efforts.

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

BackgroundThe globe artichoke (Cynara cardunculus L. var. scolymus) genome is relatively poorly explored, especially compared to those of the other major Asteraceae crops sunflower and lettuce. No SNP markers are in the public domain. We have combined the recently developed restriction-site associated DNA (RAD) approach with the Illumina DNA sequencing platform to effect the rapid and mass discovery of SNP markers for C. cardunculus.ResultsRAD tags were sequenced from the genomic DNA of three C. cardunculus mapping population parents, generating 9.7 million reads, corresponding to ~1 Gbp of sequence. An assembly based on paired ends produced ~6.0 Mbp of genomic sequence, separated into ~19,000 contigs (mean length 312 bp), of which ~21% were fragments of putative coding sequence. The shared sequences allowed for the discovery of ~34,000 SNPs and nearly 800 indels, equivalent to a SNP frequency of 5.6 per 1,000 nt, and an indel frequency of 0.2 per 1,000 nt. A sample of heterozygous SNP loci was mapped by CAPS assays and this exercise provided validation of our mining criteria. The repetitive fraction of the genome had a high representation of retrotransposon sequence, followed by simple repeats, AT-low complexity regions and mobile DNA elements. The genomic k-mers distribution and CpG rate of C. cardunculus, compared with data derived from three whole genome-sequenced dicots species, provided a further evidence of the random representation of the C. cardunculus genome generated by RAD sampling.ConclusionThe RAD tag sequencing approach is a cost-effective and rapid method to develop SNP markers in a highly heterozygous species. Our approach permitted to generate a large and robust SNP datasets by the adoption of optimized filtering criteria.

A RAD Tag Derived Marker Based Eggplant Linkage Map and the Location of QTLs Determining Anthocyanin Pigmentation

Both inter- and intra-specific maps have been developed in eggplant (Solanum melongena L.). The former benefit from an enhanced frequency of marker polymorphism, but their relevance to marker-assisted crop breeding is limited. Combining the restriction-site associated DNA strategy with high throughput sequencing has facilitated the discovery of a large number of functional single nucleotide polymorphism (SNP) markers discriminating between the two eggplant mapping population parental lines ‘305E40’ and ‘67/3’. A set of 347 de novo SNPs, together with 84 anchoring markers, were applied to the F2 mapping population bred from the cross ‘305E40’ x ‘67/3’ to construct a linkage map. In all, 415 of the 431 markers were assembled into twelve major and one minor linkage group, spanning 1,390 cM, and the inclusion of established markers allowed each linkage group to be assigned to one of the 12 eggplant chromosomes. The map was then used to discover the genetic basis of seven traits associated with anthocyanin content. Each of the traits proved to be controlled by between one and six quantitative trait loci (QTL), of which at least one was a major QTL. Exploitation of syntenic relationships between the eggplant and tomato genomes facilitated the identification of potential candidate genes for the eggplant QTLs related to anthocyanin accumulation. The intra-specific linkage map should have utility for elucidating the genetic basis of other phenotypic traits in eggplant.

Stress-Induced Biosynthesis of Dicaffeoylquinic Acids in Globe Artichoke

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

The isolation and mapping of a novel hydroxycinnamoyltransferase in the globe artichoke chlorogenic acid pathway

BackgroundThe leaves of globe artichoke and cultivated cardoon (Cynara cardunculus L.) have significant pharmaceutical properties, which mainly result from their high content of polyphenolic compounds such as monocaffeoylquinic and dicaffeoylquinic acid (DCQ), and a range of flavonoid compounds.ResultsHydroxycinnamoyl-CoA:quinate hydroxycinnamoyltransferase (HQT) encoding genes have been isolated from both globe artichoke and cultivated cardoon (GenBank accessions DQ915589 and DQ915590, respectively) using CODEHOP and PCR-RACE. A phylogenetic analysis revealed that their sequences belong to one of the major acyltransferase groups (anthranilate N-hydroxycinnamoyl/benzoyltransferase). The heterologous expression of globe artichoke HQT in E. coli showed that this enzyme can catalyze the esterification of quinic acid with caffeoyl-CoA or p-coumaroyl-CoA to generate, respectively, chlorogenic acid (CGA) and p-coumaroyl quinate. Real time PCR experiments demonstrated an increase in the expression level of HQT in UV-C treated leaves, and established a correlation between the synthesis of phenolic acids and protection against damage due to abiotic stress. The HQT gene, together with a gene encoding hydroxycinnamoyl-CoA:shikimate/quinate hydroxycinnamoyltransferase (HCT) previously isolated from globe artichoke, have been incorporated within the developing globe artichoke linkage maps.ConclusionA novel acyltransferase involved in the biosynthesis of CGA in globe artichoke has been isolated, characterized and mapped. This is a good basis for our effort to understand the genetic basis of phenylpropanoid (PP) biosynthesis in C. cardunculus.

QTL analysis of plant development and fruit traits in pepper and performance of selective phenotyping

A QTL analysis was performed to determine the genetic basis of 13 horticultural traits conditioning yield in pepper (Capsicum annuum). The mapping population was a large population of 297 recombinant inbred lines (RIL) originating from a cross between the large-fruited bell pepper cultivar ‘Yolo Wonder’ and the small-fruited chilli pepper ‘Criollo de Morelos 334’. A total of 76 QTLs were detected for 13 fruit and plant traits, grouped in 28 chromosome regions. These QTLs explained together between 7% (internode growth time) and 91% (fruit diameter) of the phenotypic variation. The QTL analysis was also performed on two subsets of 141 and 93 RILs sampled using the MapPop software. The smaller populations allowed for the detection of a reduced set of QTLs and reduced the overall percentage of trait variation explained by QTLs. The frequency of false positives as well as the individual effect of QTLs increased in reduced population sets as a result of reduced sampling. The results from the QTL analysis permitted an overall glance over the genetic architecture of traits considered by breeders for selection. Colinearities between clusters of QTLs controlling fruit traits and/or plant development in distinct pepper species and in related solanaceous crop species (tomato and eggplant) suggests that shared mechanisms control the shape and growth of different organs throughout these species.

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

BackgroundCynara cardunculus L. is an edible plant of pharmaceutical interest, in particular with respect to the polyphenolic content of its leaves. It includes three taxa: globe artichoke, cultivated cardoon, and wild cardoon. The dominating phenolics are the di-caffeoylquinic acids (such as cynarin), which are largely restricted to Cynara species, along with their precursor, chlorogenic acid (CGA). The scope of this study is to better understand CGA synthesis in this plant.ResultsA gene sequence encoding a hydroxycinnamoyltransferase (HCT) involved in the synthesis of CGA, was identified. Isolation of the gene sequence was achieved by using a PCR strategy with degenerated primers targeted to conserved regions of orthologous HCT sequences available. We have isolated a 717 bp cDNA which shares 84% aminoacid identity and 92% similarity with a tobacco gene responsible for the biosynthesis of CGA from p-coumaroyl-CoA and quinic acid. In silico studies revealed the globe artichoke HCT sequence clustering with one of the main acyltransferase groups (i.e. anthranilate N-hydroxycinnamoyl/benzoyltransferase). Heterologous expression of the full length HCT (GenBank accession DQ104740) cDNA in E. coli demonstrated that the recombinant enzyme efficiently synthesizes both chlorogenic acid and p-coumaroyl quinate from quinic acid and caffeoyl-CoA or p-coumaroyl-CoA, respectively, confirming its identity as a hydroxycinnamoyl-CoA: quinate HCT. Variable levels of HCT expression were shown among wild and cultivated forms of C. cardunculus subspecies. The level of expression was correlated with CGA content.ConclusionThe data support the predicted involvement of the Cynara cardunculus HCT in the biosynthesis of CGA before and/or after the hydroxylation step of hydroxycinnamoyl esters.

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

BackgroundThe globe artichoke (Cynara cardunculus var. scolymus L.) is a significant crop in the Mediterranean basin. Despite its commercial importance and its both dietary and pharmaceutical value, knowledge of its genetics and genomics remains scant. Microsatellite markers have become a key tool in genetic and genomic analysis, and we have exploited recently acquired EST (expressed sequence tag) sequence data (Composite Genome Project - CGP) to develop an extensive set of microsatellite markers.ResultsA unigene assembly was created from over 36,000 globe artichoke EST sequences, containing 6,621 contigs and 12,434 singletons. Over 12,000 of these unigenes were functionally assigned on the basis of homology with Arabidopsis thaliana reference proteins. A total of 4,219 perfect repeats, located within 3,308 unigenes was identified and the gene ontology (GO) analysis highlighted some GO terms enrichments among different classes of microsatellites with respect to their position. Sufficient flanking sequence was available to enable the design of primers to amplify 2,311 of these microsatellites, and a set of 300 was tested against a DNA panel derived from 28 C. cardunculus genotypes. Consistent amplification and polymorphism was obtained from 236 of these assays. Their polymorphic information content (PIC) ranged from 0.04 to 0.90 (mean 0.66). Between 176 and 198 of the assays were informative in at least one of the three available mapping populations.ConclusionEST-based microsatellites have provided a large set of de novo genetic markers, which show significant amounts of polymorphism both between and within the three taxa of C. cardunculus. They are thus well suited as assays for phylogenetic analysis, the construction of genetic maps, marker-assisted breeding, transcript mapping and other genomic applications in the species.