Cristian Del Fabbro

The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla

The analysis of the first plant genomes provided unexpected evidence for genome duplication events in species that had previously been considered as true diploids on the basis of their genetics. These polyploidization events may have had important consequences in plant evolution, in particular for species radiation and adaptation and for the modulation of functional capacities. Here we report a high-quality draft of the genome sequence of grapevine (Vitis vinifera) obtained from a highly homozygous genotype. The draft sequence of the grapevine genome is the fourth one produced so far for flowering plants, the second for a woody species and the first for a fruit crop (cultivated for both fruit and beverage). Grapevine was selected because of its important place in the cultural heritage of humanity beginning during the Neolithic period. Several large expansions of gene families with roles in aromatic features are observed. The grapevine genome has not undergone recent genome duplication, thus enabling the discovery of ancestral traits and features of the genetic organization of flowering plants. This analysis reveals the contribution of three ancestral genomes to the grapevine haploid content. This ancestral arrangement is common to many dicotyledonous plants but is absent from the genome of rice, which is a monocotyledon. Furthermore, we explain the chronology of previously described whole-genome duplication events in the evolution of flowering plants.

The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution

Rosaceae is the most important fruit-producing clade, and its key commercially relevant genera (Fragaria, Rosa, Rubus and Prunus) show broadly diverse growth habits, fruit types and compact diploid genomes. Peach, a diploid Prunus species, is one of the best genetically characterized deciduous trees. Here we describe the high-quality genome sequence of peach obtained from a completely homozygous genotype. We obtained a complete chromosome-scale assembly using Sanger whole-genome shotgun methods. We predicted 27,852 protein-coding genes, as well as noncoding RNAs. We investigated the path of peach domestication through whole-genome resequencing of 14 Prunus accessions. The analyses suggest major genetic bottlenecks that have substantially shaped peach genome diversity. Furthermore, comparative analyses showed that peach has not undergone recent whole-genome duplication, and even though the ancestral triplicated blocks in peach are fragmentary compared to those in grape, all seven paleosets of paralogs from the putative paleoancestor are detectable.

An Extensive Evaluation of Read Trimming Effects on Illumina NGS Data Analysis

Next Generation Sequencing is having an extremely strong impact in biological and medical research and diagnostics, with applications ranging from gene expression quantification to genotyping and genome reconstruction. Sequencing data is often provided as raw reads which are processed prior to analysis 1 of the most used preprocessing procedures is read trimming, which aims at removing low quality portions while preserving the longest high quality part of a NGS read. In the current work, we evaluate nine different trimming algorithms in four datasets and three common NGS-based applications (RNA-Seq, SNP calling and genome assembly). Trimming is shown to increase the quality and reliability of the analysis, with concurrent gains in terms of execution time and computational resources needed.

High throughput approaches reveal splicing of primary microRNA transcripts and tissue specific expression of mature microRNAs in Vitis vinifera

BackgroundMicroRNAs are short (~21 base) single stranded RNAs that, in plants, are generally coded by specific genes and cleaved specifically from hairpin precursors. MicroRNAs are critical for the regulation of multiple developmental, stress related and other physiological processes in plants. The recent annotation of the genome of the grapevine (Vitis vinifera L.) allowed the identification of many putative conserved microRNA precursors, grouped into multiple gene families.ResultsHere we use oligonucleotide arrays to provide the first indication that many of these microRNAs show differential expression patterns between tissues and during the maturation of fruit in the grapevine. Furthermore we demonstrate that whole transcriptome sequencing and deep-sequencing of small RNA fractions can be used both to identify which microRNA precursors are expressed in different tissues and to estimate genomic coordinates and patterns of splicing and alternative splicing for many primary miRNA transcripts.ConclusionOur results show that many microRNAs are differentially expressed in different tissues and during fruit maturation in the grapevine. Furthermore, the demonstration that whole transcriptome sequencing can be used to identify candidate splicing events and approximate primary microRNA transcript coordinates represents a significant step towards the large-scale elucidation of mechanisms regulating the expression of microRNAs at the transcriptional and post-transcriptional levels.

rNA: a fast and accurate short reads numerical aligner

SUMMARY The advent of high-throughput sequencers (HTS) introduced the need of new tools in order to analyse the large amount of data that those machines are able to produce. The mandatory first step for a wide range of analyses is the alignment of the sequences against a reference genome. We present a major update to our rNA (randomized Numerical Aligner) tool. The main feature of rNA is the fact that it achieves an accuracy greater than the majority of other tools in a feasible amount of time. rNA executables and source codes are freely downloadable at http://iga-rna.sourceforge.net/. CONTACT vezzi@appliedgenomics.org; delfabbro@appliedgenomics.org SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

ERNE-BS5: aligning BS-treated sequences by multiple hits on a 5-letters alphabet

Cytosine methylation is a DNA modification that has great impact on the regulation of gene expression and important implications for the biology and health of several living beings, including humans. Bisulfite conversion followed by next generation sequencing (BS-seq) of DNA is the gold standard technique used to detect DNA methylation at single-base resolution on a genome scale through the identification of 5-methylcytosine (5-mC). However, by converting unmethylated cytosines into thymines, BS-seq poses computational challenges to read alignment and aggravates the issue of multiple hits due to the ambiguity raised by the reduced sequence complexity. Here we present ERNE-BS5 (Extended Randomized Numerical alignEr - BiSulfite 5), an aligning program developed to efficiently map BS-treated reads against large genomes (e.g., human). To achieve this goal we have implemented three different ideas: (i) we use a 5-letters alphabet for storing methylation information, (ii) we use a weighted context-aware Hamming distance to identify a T coming from an unmethylated C context, and (iii) we use an iterative process to position multiple-hit reads starting from a preliminary map built using single-hit alignments. The map is corrected and extended at each cycle using the alignments added in the previous iteration. ERNE-BS5 is based on a new improved version of the rNA [20] aligning software with a more efficient core. ERNE (Extended Randomized Numerical alignEr) is a short string alignment package whose goal is to provide an all-inclusive set of tools to handle short reads. ERNE comprises: ERNE-MAP, ERNE-DMAP, ERNE-FILTER, ERNE-VISUAL, and, from now on, ERNE-BS5. ERNE is free software and distributed with an Open Source License (GPL V3) and can be downloaded at: http://erne.sourceforge.net

GAM: Genomic Assemblies Merger: A Graph Based Method to Integrate Different Assemblies

Many software tools are currently available to solve the hard goal of assembling millions of fragments produced in sequencing projects. Such a variety includes packages for long and short reads, generated by classical and next-generation sequencing technologies. Often the result produced by different tools can diverge---sometime significantly---for many reasons: the underlying algorithm, the data structures employed, the heuristics implemented, default parameters, etc. On the ground of the above considerations, we were motivated in developing a methodology which may both guide in a comparison of different assemblers output and improve the overall quality of the genome assembly sequences,by merging the sequences produced by different assembly programs.

Grapevine field experiments reveal the contribution of genotype, the influence of environment and the effect of their interaction (GxE) on berry transcriptome

Changes in the performance of genotypes in different environments are defined as genotype × environment (G×E) interactions. In grapevine (Vitis vinifera), complex interactions between different genotypes and climate, soil and farming practices yield unique berry qualities. However, the molecular basis of this phenomenon remains unclear. To dissect the basis of grapevine G×E interactions we characterized berry transcriptome plasticity, the genome methylation landscape and within-genotype allelic diversity in two genotypes cultivated in three different environments over two vintages. We identified, through a novel data-mining pipeline, genes with expression profiles that were: unaffected by genotype or environment, genotype-dependent but unaffected by the environment, environmentally-dependent regardless of genotype, and G×E-related. The G×E-related genes showed different degrees of within-cultivar allelic diversity in the two genotypes and were enriched for stress responses, signal transduction and secondary metabolism categories. Our study unraveled the mutual relationships between genotypic and environmental variables during G×E interaction in a woody perennial species, providing a reference model to explore how cultivated fruit crops respond to diverse environments. Also, the pivotal role of vineyard location in determining the performance of different varieties, by enhancing berry quality traits, was unraveled.

mrNA: The MPI Randomized Numerical Aligner

The advent of Next Generation Sequencers (NGS)has driven the necessity to design new and more sophisticated tools in order to cope with the huge amount of data produced by these novel technologies. String alignment against a genome reference is the first and most important phase in every(re)-sequencing project. Recently, distributed tools able to align large amounts of sequences using clusters or clouds of computers, have been put forward. The aim of this work is to propose a new tool named mrNA (the MPI version of the original rNA program) able to align NGS data using a cluster of computers. mrNA was designed to tackle the main computational bottleneck of all classical parallel implementation of aligners: references longer than 4 Gbp. mrNA, together with rNA, are open source programs downloadable athttp://iga-rna.sourceforge.net/.