The Role of Italy in the Use of Advanced Plant Genomic Techniques on Fruit Trees: State of the Art and Future Perspectives.

Climate change is deeply impacting the food chain production, lowering quality and yield. In this context, the international scientific community has dedicated many efforts to enhancing resilience and sustainability in agriculture. Italy is among the main European producers of several fruit trees; therefore, national research centers and universities undertook several initiatives to maintain the specificity of the 'Made in Italy' label. Despite their importance, fruit crops are suffering from difficulties associated with the conventional breeding approaches, especially in terms of financial commitment, land resources availability, and long generation times. The 'new genomic techniques' (NGTs), renamed in Italy as 'technologies for assisted evolution' (TEAs), reduce the time required to obtain genetically improved cultivars while precisely targeting specific DNA sequences. This review aims to illustrate the role of the Italian scientific community in the use of NGTs, with a specific focus on Citrus, grapevine, apple, pear, chestnut, strawberry, peach, and kiwifruit. For each crop, the key genes and traits on which the scientific community is working, as well as the technological improvements and advancements on the regeneration of local varieties, are presented. Lastly, a focus is placed on the legal aspects in the European and in Italian contexts.

Quantitative Trait Loci Mapping and Identification of Candidate Genes Linked to Fruit Acidity in Apricot (Prunus armeniaca L.).

Apricot breeding programs could be strongly improved by the availability of molecular markers linked to the main fruit quality traits. Fruit acidity is one of the key factors in consumer acceptance, but despite its importance, the molecular bases of this trait are still poorly understood. In order to increase the genetic knowledge on the fruit acidity, an F1 apricot population ('Lito' × 'BO81604311') has been phenotyped for titratable acidity and juice pH for the three following years. In addition, the contents of the main organic acids of the juice (malate, citrate, and quinate) were also evaluated. A Gaussian distribution was observed for most of the traits in this progeny, confirming their quantitative inheritance. An available simple sequence repeat (SSR)-based molecular map, implemented with new markers in specific genomic regions, was used to perform a quantitative trait loci (QTL) analysis. The molecular map was also anchored to the recently published apricot genome sequence of 'Stella.' Several major QTLs linked to fruit acidity-related traits have been identified both in the 'Lito' (no. 21) and 'BO81604311' (no. 13), distributed in five linkage groups (LG 4, 5, 6, 7, and 8). Some of these QTLs show good stability between years and their linked markers were used to identify candidate genes in specific QTLs genomic regions.

Resistance to Sharka in Apricot: Comparison of Phase-Reconstructed Resistant and Susceptible Haplotypes of 'Lito' Chromosome 1 and Analysis of Candidate Genes.

Sharka, a common disease among most stone fruit crops, is caused by the Plum Pox Virus (PPV). Resistant genotypes have been found in apricot (Prunus armeniaca L.), one of which-the cultivar 'Lito' heterozygous for the resistance-has been used to map a major quantitative trait locus (QTL) on linkage group 1, following a pseudo-test-cross mating design with 231 individuals. In addition, 19 SNP markers were selected from among the hundreds previously developed, which allowed the region to be limited to 236 kb on chromosome 1. A 'Lito' bacterial artificial chromosome (BAC) library was produced, screened with markers of the region, and positive BAC clones were sequenced. Resistant (R) and susceptible (S) haplotypes were assembled independently. To refine the assembly, the whole genome of 'Lito' was sequenced to high coverage (98×) using PacBio technology, enabling the development of a detailed assembly of the region that was able to predict and annotate the genes in the QTL region. The selected cultivar 'Lito' allowed not only to discriminate structural variants between the two haplotypic regions but also to distinguish specific allele expression, contributing towards mining the PPVres locus. In light of these findings, genes previously indicated (i.e., MATHd genes) to have a possible role in PPV resistance were further analyzed, and new candidates were discussed. Although the results are not conclusive, the accurate and independent assembly of R and S haplotypes of 'Lito' is a valuable resource to predict and test alternative transcription and regulation mechanisms underpinning PPV resistance.

Characterization of 25 full-length S-RNase alleles, including flanking regions, from a pool of resequenced apple cultivars.

KEY MESSAGE: Data obtained from Illumina resequencing of 63 apple cultivars were used to obtain full-length S-RNase sequences using a strategy based on both alignment and de novo assembly of reads. The reproductive biology of apple is regulated by the S-RNase-based gametophytic self-incompatibility system, that is genetically controlled by the single, multi-genic and multi-allelic S locus. Resequencing of apple cultivars provided a huge amount of genetic data, that can be aligned to the reference genome in order to characterize variation to a genome-wide level. However, this approach is not immediately adaptable to the S-locus, due to some peculiar features such as the high degree of polymorphism, lack of colinearity between haplotypes and extensive presence of repetitive elements. In this study we describe a dedicated procedure aimed at characterizing S-RNase alleles from resequenced cultivars. The S-genotype of 63 apple accessions is reported; the full length coding sequence was determined for the 25 S-RNase alleles present in the 63 resequenced cultivars; these included 10 previously incomplete sequences (S 5 , S 6a , S 6b , S 8 , S 11 , S 23 , S 39 , S 46 , S 50 and S 58 ). Moreover, sequence divergence clearly suggests that alleles S 6a and S 6b , proposed to be neutral variants of the same alleles, should be instead considered different specificities. The promoter sequences have also been analyzed, highlighting regions of homology conserved among all the alleles.

Genome-Wide Association Mapping of Flowering and Ripening Periods in Apple.

Deciphering the genetic control of flowering and ripening periods in apple is essential for breeding cultivars adapted to their growing environments. We implemented a large Genome-Wide Association Study (GWAS) at the European level using an association panel of 1,168 different apple genotypes distributed over six locations and phenotyped for these phenological traits. The panel was genotyped at a high-density of SNPs using the Axiom®Apple 480 K SNP array. We ran GWAS with a multi-locus mixed model (MLMM), which handles the putatively confounding effect of significant SNPs elsewhere on the genome. Genomic regions were further investigated to reveal candidate genes responsible for the phenotypic variation. At the whole population level, GWAS retained two SNPs as cofactors on chromosome 9 for flowering period, and six for ripening period (four on chromosome 3, one on chromosome 10 and one on chromosome 16) which, together accounted for 8.9 and 17.2% of the phenotypic variance, respectively. For both traits, SNPs in weak linkage disequilibrium were detected nearby, thus suggesting the existence of allelic heterogeneity. The geographic origins and relationships of apple cultivars accounted for large parts of the phenotypic variation. Variation in genotypic frequency of the SNPs associated with the two traits was connected to the geographic origin of the genotypes (grouped as North+East, West and South Europe), and indicated differential selection in different growing environments. Genes encoding transcription factors containing either NAC or MADS domains were identified as major candidates within the small confidence intervals computed for the associated genomic regions. A strong microsynteny between apple and peach was revealed in all the four confidence interval regions. This study shows how association genetics can unravel the genetic control of important horticultural traits in apple, as well as reduce the confidence intervals of the associated regions identified by linkage mapping approaches. Our findings can be used for the improvement of apple through marker-assisted breeding strategies that take advantage of the accumulating additive effects of the identified SNPs.

The Peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity.

BACKGROUND: The availability of the peach genome sequence has fostered relevant research in peach and related Prunus species enabling the identification of genes underlying important horticultural traits as well as the development of advanced tools for genetic and genomic analyses. The first release of the peach genome (Peach v1.0) represented a high-quality WGS (Whole Genome Shotgun) chromosome-scale assembly with high contiguity (contig L50 214.2 kb), large portions of mapped sequences (96%) and high base accuracy (99.96%). The aim of this work was to improve the quality of the first assembly by increasing the portion of mapped and oriented sequences, correcting misassemblies and improving the contiguity and base accuracy using high-throughput linkage mapping and deep resequencing approaches. RESULTS: Four linkage maps with 3,576 molecular markers were used to improve the portion of mapped and oriented sequences (from 96.0% and 85.6% of Peach v1.0 to 99.2% and 98.2% of v2.0, respectively) and enabled a more detailed identification of discernible misassemblies (10.4 Mb in total). The deep resequencing approach fixed 859 homozygous SNPs (Single Nucleotide Polymorphisms) and 1347 homozygous indels. Moreover, the assembled NGS contigs enabled the closing of 212 gaps with an improvement in the contig L50 of 19.2%. CONCLUSIONS: The improved high quality peach genome assembly (Peach v2.0) represents a valuable tool for the analysis of the genetic diversity, domestication, and as a vehicle for genetic improvement of peach and related Prunus species. Moreover, the important phylogenetic position of peach and the absence of recent whole genome duplication (WGD) events make peach a pivotal species for comparative genomics studies aiming at elucidating plant speciation and diversification processes.

Fighting Sharka in Peach: Current Limitations and Future Perspectives.

Sharka, caused by Plum Pox Virus (PPV), is by far the most important infectious disease of peach [P. persica (L.) Batsch] and other Prunus species. The progressive spread of the virus in many important growing areas throughout Europe poses serious issues to the economic sustainability of stone fruit crops, peach in particular. The adoption of internationally agreed-upon rules for diagnostic tests, strain-specific monitoring schemes and spatial-temporal modeling of virus spread, are all essential for a more effective sharka containment. The EU regulations on nursery activity should be modified based on the zone delimitation of PPV presence, limiting open-field production of propagation materials only to virus-free areas. Increasing the efficiency of preventive measures should be augmented by the short-term development of resistant cultivars. Putative sources of resistance/tolerance have been recently identified in peach germplasm, although the majority of novel resistant sources to PPV-M have been found in almond. However, the complexity of introgression from related-species imposes the search for alternative strategies. The use of genetic engineering, particularly RNA interference (RNAi)-based approaches, appears as one of the most promising perspectives to introduce a durable resistance to PPV in peach germplasm, notwithstanding the well-known difficulties of in vitro plant regeneration in this species. In this regard, rootstock transformation to induce RNAi-mediated systemic resistance would avoid the transformation of numerous commercial cultivars, and may alleviate consumer resistance to the use of GM plants.

Analysis of the genetic diversity and structure across a wide range of germplasm reveals prominent gene flow in apple at the European level.

BACKGROUND: The amount and structure of genetic diversity in dessert apple germplasm conserved at a European level is mostly unknown, since all diversity studies conducted in Europe until now have been performed on regional or national collections. Here, we applied a common set of 16 SSR markers to genotype more than 2,400 accessions across 14 collections representing three broad European geographic regions (North + East, West and South) with the aim to analyze the extent, distribution and structure of variation in the apple genetic resources in Europe. RESULTS: A Bayesian model-based clustering approach showed that diversity was organized in three groups, although these were only moderately differentiated (FST = 0.031). A nested Bayesian clustering approach allowed identification of subgroups which revealed internal patterns of substructure within the groups, allowing a finer delineation of the variation into eight subgroups (FST = 0.044). The first level of stratification revealed an asymmetric division of the germplasm among the three groups, and a clear association was found with the geographical regions of origin of the cultivars. The substructure revealed clear partitioning of genetic groups among countries, but also interesting associations between subgroups and breeding purposes of recent cultivars or particular usage such as cider production. Additional parentage analyses allowed us to identify both putative parents of more than 40 old and/or local cultivars giving interesting insights in the pedigree of some emblematic cultivars. CONCLUSIONS: The variation found at group and subgroup levels may reflect a combination of historical processes of migration/selection and adaptive factors to diverse agricultural environments that, together with genetic drift, have resulted in extensive genetic variation but limited population structure. The European dessert apple germplasm represents an important source of genetic diversity with a strong historical and patrimonial value. The present work thus constitutes a decisive step in the field of conservation genetics. Moreover, the obtained data can be used for defining a European apple core collection useful for further identification of genomic regions associated with commercially important horticultural traits in apple through genome-wide association studies.

Genetic diversity analysis of the cucurbit powdery mildew fungus Podosphaera xanthii suggests a clonal population structure.

The sexual stage of Podosphaera xanthii is rarely found worldwide. However, chasmothecia are frequently recorded in northern Italy, suggesting the presence of an actively mating population. With the aim of investigating the genetic structure of the Italian population with respect to populations from other countries, genetic diversity analysis was performed both on 92 isolates from European and American countries by multilocus sequence typing (MLST) and on 59 isolates by amplified fragment length polymorphism (AFLP) methods. Mating type frequencies were tested for random mating and two-locus linkage disequilibrium (LD) analysis was performed. Results showed very low levels of genetic diversity: MLST showed no variations in eight housekeeping gene fragments and, accordingly, UPGMA dendrogram from AFLP data showed a high similarity (0.91-1.00 simple matching similarity coefficient) between isolates. Moreover, the random mating test showed no deviations from mating-type 1:1 ratio in the Italian population but deviations were observed in populations from Europe and American countries while two-locus LD analysis showed the presence of significant LD. The results suggest that the populations of P. xanthii are likely to be predominantly clonal, and asexual reproduction, producing a huge amount of conidia, appears to be the predominant type of reproduction of the species.

A unique mutation in a MYB gene cosegregates with the nectarine phenotype in peach.

Nectarines play a key role in peach industry; the fuzzless skin has implications for consumer acceptance. The peach/nectarine (G/g) trait was described as monogenic and previously mapped on chromosome 5. Here, the position of the G locus was delimited within a 1.1 cM interval (635 kb) based on linkage analysis of an F2 progeny from the cross 'Contender' (C, peach) x 'Ambra' (A, nectarine). Careful inspection of the genes annotated in the corresponding genomic sequence (Peach v1.0), coupled with variant discovery, led to the identification of MYB gene PpeMYB25 as a candidate for trichome formation on fruit skin. Analysis of genomic re-sequencing data from five peach/nectarine accessions pointed to the insertion of a LTR retroelement in exon 3 of the PpeMYB25 gene as the cause of the recessive glabrous phenotype. A functional marker (indelG) developed on the LTR insertion cosegregated with the trait in the CxA F2 progeny and was validated on a broad panel of genotypes, including all known putative donors of the nectarine trait. This marker was shown to efficiently discriminate between peach and nectarine plants, indicating that a unique mutational event gave rise to the nectarine trait and providing a useful diagnostic tool for early seedling selection in peach breeding programs.

Fine mapping and identification of a candidate gene for a major locus controlling maturity date in peach.

BACKGROUND: Maturity date (MD) is a crucial factor for marketing of fresh fruit, especially those with limited shelf-life such as peach (Prunus persica L. Batsch): selection of several cultivars with differing MD would be advantageous to cover and extend the marketing season. Aims of this work were the fine mapping and identification of candidate genes for the major maturity date locus previously identified on peach linkage group 4. To improve genetic resolution of the target locus two F2 populations derived from the crosses Contender x Ambra (CxA, 306 individuals) and PI91459 (NJ Weeping) x Bounty (WxBy, 103 individuals) were genotyped with the Sequenom and 9K Illumina Peach Chip SNP platforms, respectively. RESULTS: Recombinant individuals from the WxBy F2 population allowed the localisation of maturity date locus to a 220 kb region of the peach genome. Among the 25 annotated genes within this interval, functional classification identified ppa007577m and ppa008301m as the most likely candidates, both encoding transcription factors of the NAC (NAM/ATAF1, 2/CUC2) family. Re-sequencing of the four parents and comparison with the reference genome sequence uncovered a deletion of 232 bp in the upstream region of ppa007577m that is homozygous in NJ Weeping and heterozygous in Ambra, Bounty and the WxBy F1 parent. However, this variation did not segregate in the CxA F2 population being the CxA F1 parent homozygous for the reference allele. The second gene was thus examined as a candidate for maturity date. Re-sequencing of ppa008301m, showed an in-frame insertion of 9 bp in the last exon that co-segregated with the maturity date locus in both CxA and WxBy F2 populations. CONCLUSIONS: Using two different segregating populations, the map position of the maturity date locus was refined from 3.56 Mb to 220 kb. A sequence variant in the NAC gene ppa008301m was shown to co-segregate with the maturity date locus, suggesting this gene as a candidate controlling ripening time in peach. If confirmed on other genetic materials, this variant may be used for marker-assisted breeding of new cultivars with differing maturity date.

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.

Molecular bases and evolutionary dynamics of self-incompatibility in the Pyrinae (Rosaceae).

The molecular bases of the gametophytic self-incompatibility (GSI) system of species of the subtribe Pyrinae (Rosaceae), such as apple and pear, have been widely studied in the last two decades. The characterization of S-locus genes and of the mechanisms underlying pollen acceptance or rejection have been topics of major interest. Besides the single pistil-side S determinant, the S-RNase, multiple related S-locus F-box genes seem to be involved in the determination of pollen S specificity. Here, we collect and review the state of the art of GSI in the Pyrinae. We emphasize recent genomic data that have contributed to unveiling the S-locus structure of the Pyrinae, and discuss their consistency with the models of self-recognition that have been proposed for Prunus and the Solanaceae. Experimental data suggest that the mechanism controlling pollen-pistil recognition specificity of the Pyrinae might fit well with the collaborative 'non-self' recognition system proposed for Petunia (Solanaceae), whereas it presents relevant differences with the mechanism exhibited by the species of the closely related genus Prunus, which uses a single evolutionarily divergent F-box gene as the pollen S determinant. The possible involvement of multiple pollen S genes in the GSI system of Pyrinae, still awaiting experimental confirmation, opens up new perspectives to our understanding of the evolution of S haplotypes, and of the evolution of S-RNase-based GSI within the Rosaceae family. Whereas S-locus genes encode the players determining self-recognition, pollen rejection in the Pyrinae seems to involve a complex cascade of downstream cellular events with significant similarities to programmed cell death.

dHPLC efficiency for semi-automated cDNA-AFLP analyses and fragment collection in the apple scab-resistance gene model.

cDNA-AFLP analysis for transcript profiling has been successfully applied to study many plant biological systems, particularly plant-microbe interactions. However, the separation of cDNA-AFLP fragments by gel electrophoresis is reported to be labor-intensive with only limited potential for automation, and the collection of differential bands for gene identification is even more cumbersome. In this work, we present the use of dHPLC (denaturing high performance liquid chromatography) and automated DNA fragment collection using the WAVE(®) System to analyze and recover cDNA-AFLP fragments. The method is successfully applied to the Malus-Venturia inaequalis interaction, making it possible to collect 66 different transcript-derived fragments for apple genes putatively involved in the defense response activated by the HcrVf2 resistance gene. The results, validated by real time quantitative RT-PCR, were consistent with the plant-pathogen interaction under investigation and this further supports the suitability of dHPLC for cDNA-AFLP transcript profiling. Features and advantages of this new approach are discussed, evincing that it is an almost fully automatable and cost-effective solution for processing large numbers of samples and collecting differential genes involved in other biological processes and non-model plants.

A qNMR approach for bitterness phenotyping and QTL identification in an F1 apricot progeny.

In apricot the bitter flavor of seeds is determined by the amount of amygdalin, a cyanogenic glucoside whose cleavage by endogenous enzymes, upon seed crushing, releases toxic hydrogen cyanide. The presence of such a poisonous compound is an obstacle to the use and commercialization of apricot seeds for human or animal nutrition. To investigate the genetic loci involved in the determination of the bitter phenotype a combined genetic and biochemical approach was used, involving a candidate gene analysis and a fine phenotyping via quantitative nuclear magnetic resonance, on an F1 apricot progeny. Seven functional markers were developed and positioned on the genetic maps of the parental lines Lito and BO81604311 and seven putative QTLs for the bitterness level were determined. In conclusion, this analysis has revealed some loci involved in the shaping of the bitterness degree; has proven the complexity of the bitter trait in apricot, reporting an high variance of the QTLs found over the years; has showed the critical importance of the phenotyping step, whose precision and accuracy is a pre-requisite when studying such a multifactorial character.

Mapping of an anthocyanin-regulating MYB transcription factor and its expression in red and green pear, Pyrus communis.

'Max Red Bartlett' is a red bud mutation of the yellow pear (Pyrus communis L.) cultivar 'Williams' (known as 'Bartlett' in North America). Anthocyanins are the most important pigments for red colour in fruits. Synthesis of anthocyanins is mediated by a number of well-characterized enzymes that include chalcone synthase (CHS), flavanone-3-hydroxylase (F3H), dihydroflavonol-4-reductase (DFR), anthocyanidin synthase (ANS), and UDP-glucose:flavonoid-3-O-glucosyltransferase (UFGT). Expression of the genes encoding these five enzymes was examined in pear fruit skin in order to elucidate the molecular mechanism for red coloration. In addition, the gene PcMYB10, encoding an R2R3 MYB transcription factor involved in anthocyanin biosynthetic pathway regulation, was isolated from both 'Williams' and 'Max Red Bartlett'. Analysis of the deduced amino acid sequence suggests that this gene is an ortholog of anthocyanin regulators known in other plant species. Its expression level was significantly higher in 'Max Red Bartlett' (red pear) compared with the original yellow variety 'Williams'. Although the map position of PcMYB10 corresponds to that of MdMYBa and MdMYB10, which control pigmentation of apple fruit skin, PcMYB10 is not directly responsible for red versus yellow colour in the two pear varieties, as the mutation underlying this difference maps to a different region of the pear genome.