Riccardo Velasco

The genome of woodland strawberry ( Fragaria vesca )

The woodland strawberry, Fragaria vesca (2n = 2x = 14), is a versatile experimental plant system. This diminutive herbaceous perennial has a small genome (240 Mb), is amenable to genetic transformation and shares substantial sequence identity with the cultivated strawberry (Fragaria × ananassa) and other economically important rosaceous plants. Here we report the draft F. vesca genome, which was sequenced to ×39 coverage using second-generation technology, assembled de novo and then anchored to the genetic linkage map into seven pseudochromosomes. This diploid strawberry sequence lacks the large genome duplications seen in other rosids. Gene prediction modeling identified 34,809 genes, with most being supported by transcriptome mapping. Genes critical to valuable horticultural traits including flavor, nutritional value and flowering time were identified. Macrosyntenic relationships between Fragaria and Prunus predict a hypothetical ancestral Rosaceae genome that had nine chromosomes. New phylogenetic analysis of 154 protein-coding genes suggests that assignment of Populus to Malvidae, rather than Fabidae, is warranted.

Genome-wide transcriptional analysis of grapevine berry ripening reveals a set of genes similarly modulated during three seasons and the occurrence of an oxidative burst at vèraison

BackgroundGrapevine (Vitis species) is among the most important fruit crops in terms of cultivated area and economic impact. Despite this relevance, little is known about the transcriptional changes and the regulatory circuits underlying the biochemical and physical changes occurring during berry development.ResultsFruit ripening in the non-climacteric crop species Vitis vinifera L. has been investigated at the transcriptional level by the use of the Affymetrix Vitis GeneChip® which contains approximately 14,500 unigenes. Gene expression data obtained from berries sampled before and after véraison in three growing years, were analyzed to identify genes specifically involved in fruit ripening and to investigate seasonal influences on the process. From these analyses a core set of 1477 genes was found which was similarly modulated in all seasons. We were able to separate ripening specific isoforms within gene families and to identify ripening related genes which appeared strongly regulated also by the seasonal weather conditions. Transcripts annotation by Gene Ontology vocabulary revealed five overrepresented functional categories of which cell wall organization and biogenesis, carbohydrate and secondary metabolisms and stress response were specifically induced during the ripening phase, while photosynthesis was strongly repressed. About 19% of the core gene set was characterized by genes involved in regulatory processes, such as transcription factors and transcripts related to hormonal metabolism and signal transduction. Auxin, ethylene and light emerged as the main stimuli influencing berry development. In addition, an oxidative burst, previously not detected in grapevine, characterized by rapid accumulation of H2O2 starting from véraison and by the modulation of many ROS scavenging enzymes, was observed.ConclusionThe time-course gene expression analysis of grapevine berry development has identified the occurrence of two well distinct phases along the process. The pre-véraison phase represents a reprogramming stage of the cellular metabolism, characterized by the expression of numerous genes involved in hormonal signalling and transcriptional regulation. The post-véraison phase is characterized by the onset of a ripening-specialized metabolism responsible for the phenotypic traits of the ripe berry. Between the two phases, at véraison, an oxidative burst and the concurrent modulation of the anti-oxidative enzymatic network was observed. The large number of regulatory genes we have identified represents a powerful new resource for dissecting the mechanisms of fruit ripening control in non-climacteric plants.

Genome-wide SNP detection, validation, and development of an 8K SNP array for apple.

As high-throughput genetic marker screening systems are essential for a range of genetics studies and plant breeding applications, the International RosBREED SNP Consortium (IRSC) has utilized the Illumina Infinium® II system to develop a medium- to high-throughput SNP screening tool for genome-wide evaluation of allelic variation in apple (Malus×domestica) breeding germplasm. For genome-wide SNP discovery, 27 apple cultivars were chosen to represent worldwide breeding germplasm and re-sequenced at low coverage with the Illumina Genome Analyzer II. Following alignment of these sequences to the whole genome sequence of ‘Golden Delicious’, SNPs were identified using SoapSNP. A total of 2,113,120 SNPs were detected, corresponding to one SNP to every 288 bp of the genome. The Illumina GoldenGate® assay was then used to validate a subset of 144 SNPs with a range of characteristics, using a set of 160 apple accessions. This validation assay enabled fine-tuning of the final subset of SNPs for the Illumina Infinium® II system. The set of stringent filtering criteria developed allowed choice of a set of SNPs that not only exhibited an even distribution across the apple genome and a range of minor allele frequencies to ensure utility across germplasm, but also were located in putative exonic regions to maximize genotyping success rate. A total of 7867 apple SNPs was established for the IRSC apple 8K SNP array v1, of which 5554 were polymorphic after evaluation in segregating families and a germplasm collection. This publicly available genomics resource will provide an unprecedented resolution of SNP haplotypes, which will enable marker-locus-trait association discovery, description of the genetic architecture of quantitative traits, investigation of genetic variation (neutral and functional), and genomic selection in apple.

Mitochondrial DNA of Vitis vinifera and the Issue of Rampant Horizontal Gene Transfer

The mitochondrial genome of grape (Vitis vinifera), the largest organelle genome sequenced so far, is presented. The genome is 773,279 nt long and has the highest coding capacity among known angiosperm mitochondrial DNAs (mtDNAs). The proportion of promiscuous DNA of plastid origin in the genome is also the largest ever reported for an angiosperm mtDNA, both in absolute and relative terms. In all, 42.4% of chloroplast genome of Vitis has been incorporated into its mitochondrial genome. In order to test if horizontal gene transfer (HGT) has also contributed to the gene content of the grape mtDNA, we built phylogenetic trees with the coding sequences of mitochondrial genes of grape and their homologs from plant mitochondrial genomes. Many incongruent gene tree topologies were obtained. However, the extent of incongruence between these gene trees is not significantly greater than that observed among optimal trees for chloroplast genes, the common ancestry of which has never been in doubt. In both cases, we attribute this incongruence to artifacts of tree reconstruction, insufficient numbers of characters, and gene paralogy. This finding leads us to question the recent phylogenetic interpretation of Bergthorsson et al. (2003, 2004) and Richardson and Palmer (2007) that rampant HGT into the mtDNA of Amborella best explains phylogenetic incongruence between mitochondrial gene trees for angiosperms. The only evidence for HGT into the Vitis mtDNA found involves fragments of two coding sequences stemming from two closteroviruses that cause the leaf roll disease of this plant. We also report that analysis of sequences shared by both chloroplast and mitochondrial genomes provides evidence for a previously unknown gene transfer route from the mitochondrion to the chloroplast.

A Dense Single-Nucleotide Polymorphism-Based Genetic Linkage Map of Grapevine (Vitis vinifera L.) Anchoring Pinot Noir Bacterial Artificial Chromosome Contigs

The construction of a dense genetic map for Vitis vinifera and its anchoring to a BAC-based physical map is described: it includes 994 loci mapped onto 19 linkage groups, corresponding to the basic chromosome number of Vitis. Spanning 1245 cM with an average distance of 1.3 cM between adjacent markers, the map was generated from the segregation of 483 single-nucleotide polymorphism (SNP)-based genetic markers, 132 simple sequence repeats (SSRs), and 379 AFLP markers in a mapping population of 94 F1 individuals derived from a V. vinifera cross of the cultivars Syrah and Pinot Noir. Of these markers, 623 were anchored to 367 contigs that are included in a physical map produced from the same clone of Pinot Noir and covering 352 Mbp. On the basis of contigs containing two or more genetically mapped markers, region-dependent estimations of physical and recombinational distances are presented. The markers used in this study include 118 SSRs common to an integrated map derived from five segregating populations of V. vinifera. The positions of these SSR markers in the two maps are conserved across all Vitis linkage groups. The addition of SNP-based markers introduces polymorphisms that are easy to database, are useful for evolutionary studies, and significantly increase the density of the map. The map provides the most comprehensive view of the Vitis genome reported to date and will be relevant for future studies on structural and functional genomics and genetic improvement.

Ripening and Genotype Control Stilbene Accumulation in Healthy Grapes

In grapes, stilbene synthesis occurs in the skin, and it is induced by biotic and abiotic stresses. To date, experimental evidence of a constitutive production of resveratrols in healthy grape is scarce and not conclusive. The aim of the present work was to investigate stilbene biosynthesis in healthy grapes both at biochemical and molecular levels. By measuring the concentration of resveratrols in ripe berries of 78 Vitis vinifera varieties for 3 years, we could identify significant differences among genotypes, providing the first tentative varietal classification based on resveratrol content. Furthermore, an increasing stilbene accumulation from veraison to ripening phase was also observed. Using real-time reverse transcription polymerase chain reaction and a berry-specific cDNA array, gene expression analysis was carried out on two distinct pools of berries belonging to the high and low resveratrol producers and on three berry developmental stages. The stilbene synthase, phenylalanine ammonia-lyase, and 4-coumarate-CoA ligase expression profiles showed an increasing concentration of these transcripts from véraison to maturity and a higher accumulation in the grape of high resveratrol producers. Macroarray data analysis revealed that high resveratrol levels are also accompanied by the up-regulation of genes involved in plant defense and the concomitant underexpression of genes related to the ripening process and to indole alkaloid synthesis.

Saturated linkage map construction in Rubus idaeus using genotyping by sequencing and genome-independent imputation.

BackgroundRapid development of highly saturated genetic maps aids molecular breeding, which can accelerate gain per breeding cycle in woody perennial plants such as Rubus idaeus (red raspberry). Recently, robust genotyping methods based on high-throughput sequencing were developed, which provide high marker density, but result in some genotype errors and a large number of missing genotype values. Imputation can reduce the number of missing values and can correct genotyping errors, but current methods of imputation require a reference genome and thus are not an option for most species.ResultsGenotyping by Sequencing (GBS) was used to produce highly saturated maps for a R. idaeus pseudo-testcross progeny. While low coverage and high variance in sequencing resulted in a large number of missing values for some individuals, a novel method of imputation based on maximum likelihood marker ordering from initial marker segregation overcame the challenge of missing values, and made map construction computationally tractable. The two resulting parental maps contained 4521 and 2391 molecular markers spanning 462.7 and 376.6 cM respectively over seven linkage groups. Detection of precise genomic regions with segregation distortion was possible because of map saturation. Microsatellites (SSRs) linked these results to published maps for cross-validation and map comparison.ConclusionsGBS together with genome-independent imputation provides a rapid method for genetic map construction in any pseudo-testcross progeny. Our method of imputation estimates the correct genotype call of missing values and corrects genotyping errors that lead to inflated map size and reduced precision in marker placement. Comparison of SSRs to published R. idaeus maps showed that the linkage maps constructed with GBS and our method of imputation were robust, and marker positioning reliable. The high marker density allowed identification of genomic regions with segregation distortion in R. idaeus, which may help to identify deleterious alleles that are the basis of inbreeding depression in the species.

An Ancient Duplication of Apple MYB Transcription Factors Is Responsible for Novel Red Fruit-Flesh Phenotypes

Anthocyanin accumulation is coordinated in plants by a number of conserved transcription factors. In apple (Malus × domestica), an R2R3 MYB transcription factor has been shown to control fruit flesh and foliage anthocyanin pigmentation (MYB10) and fruit skin color (MYB1). However, the pattern of expression and allelic variation at these loci does not explain all anthocyanin-related apple phenotypes. One such example is an open-pollinated seedling of cv Sangrado that has green foliage and develops red flesh in the fruit cortex late in maturity. We used methods that combine plant breeding, molecular biology, and genomics to identify duplicated MYB transcription factors that could control this phenotype. We then demonstrated that the red-flesh cortex phenotype is associated with enhanced expression of MYB110a, a paralog of MYB10. Functional characterization of MYB110a showed that it was able to up-regulate anthocyanin biosynthesis in tobacco (Nicotiana tabacum). The chromosomal location of MYB110a is consistent with a whole-genome duplication event that occurred during the evolution of apple within the Maloideae family. Both MYB10 and MYB110a have conserved function in some cultivars, but they differ in their expression pattern and response to fruit maturity.

Signaling Pathways Mediating the Induction of Apple Fruitlet Abscission

Apple (Malus × domestica) represents an interesting model tree crop for studying fruit abscission. The physiological fruitlet drop occurring in this species can be easily magnified by using thinning chemicals, such as benzyladenine (BA), to obtain fruits with improved quality and marketability. Despite the economic importance of this process, the molecular determinants of apple fruitlet abscission are still unknown. In this research, BA was used to obtain fruitlet populations with different abscission potentials to be analyzed by means of a newly released 30K oligonucleotide microarray. RNAs were extracted from cortex and seed of apple fruitlets sampled over a 4-d time course, during which BA triggers fruit drop, and used for microarray hybridization. Transcriptomic profiles of persisting and abscising fruitlets were tested for statistical association with abscission potential, allowing us to identify molecular signatures strictly related to fruit destiny. A hypothetical model for apple fruitlet abscission was obtained by putting together available transcriptomic and metabolomic data. According to this model, BA treatment would establish a nutritional stress within the tree that is primarily perceived by the fruitlet cortex whose growth is blocked by resembling the ovary growth inhibition found in other species. In weaker fruits, this stress is soon visible also at the seed level, likely transduced via reactive oxygen species/sugar and hormones signaling cross talk, and followed by a block of embryogenesis and the consequent activation of the abscission zone.

The Draft Genome Sequence of European Pear (Pyrus communis L. 'Bartlett')

We present a draft assembly of the genome of European pear (Pyrus communis) ‘Bartlett’. Our assembly was developed employing second generation sequencing technology (Roche 454), from single-end, 2 kb, and 7 kb insert paired-end reads using Newbler (version 2.7). It contains 142,083 scaffolds greater than 499 bases (maximum scaffold length of 1.2 Mb) and covers a total of 577.3 Mb, representing most of the expected 600 Mb Pyrus genome. A total of 829,823 putative single nucleotide polymorphisms (SNPs) were detected using re-sequencing of ‘Louise Bonne de Jersey’ and ‘Old Home’. A total of 2,279 genetically mapped SNP markers anchor 171 Mb of the assembled genome. Ab initio gene prediction combined with prediction based on homology searching detected 43,419 putative gene models. Of these, 1219 proteins (556 clusters) are unique to European pear compared to 12 other sequenced plant genomes. Analysis of the expansin gene family provided an example of the quality of the gene prediction and an insight into the relationships among one class of cell wall related genes that control fruit softening in both European pear and apple (Malus×domestica). The ‘Bartlett’ genome assembly v1.0 (http://www.rosaceae.org/species/pyrus/pyrus_communis/genome_v1.0) is an invaluable tool for identifying the genetic control of key horticultural traits in pear and will enable the wide application of marker-assisted and genomic selection that will enhance the speed and efficiency of pear cultivar development.