Hypermethylation and small RNA expression are associated with increased age in almond (Prunus dulcis [Mill.] D.A. Webb) accessions

The focus of this study is to profile changes in DNA methylation and small RNA expression occurring with increased age in almond breeding germplasm to identify possible biomarkers of age that can be used to assess the potential of individuals to develop aging-related disorders. To profile DNA methylation in almond germplasm, 70 methylomes were generated from almond individuals representing three age cohorts (11, 7, and 2 years old) using an enzymatic methyl-seq approach followed by analysis to call differentially methylated regions (DMRs) within these cohorts. Small RNA (sRNA) expression was profiled in three breeding selections, each from two age cohorts (1 and 6 years old), using sRNA-Seq followed by differential expression analysis. Weighted chromosome-level methylation analysis reveals hypermethylation in 11-year-old almond breeding selections when compared to 2-year-old selections in the CG and CHH contexts. Seventeen consensus DMRs were identified in all age contrasts. sRNA expression differed significantly between the two age cohorts tested, with significantly decreased expression in sRNAs in the 6-year-old selections compared to the 1-year-old. Almond shows a pattern of hypermethylation and decreased sRNA expression with increased age. Identified DMRs and differentially expressed sRNAs could function as putative biomarkers of age following validation in additional age groups.

A review of plant epigenetics through the lens of almond

While genomes were originally seen as static entities that stably held and organized genetic information, recent advances in sequencing have uncovered the dynamic nature of the genome. New conceptualizations of the genome include complex relationships between the environment and gene expression that must be maintained, regulated, and sometimes even transmitted over generations. The discovery of epigenetic mechanisms has allowed researchers to understand how traits like phenology, plasticity, and fitness can be altered without changing the underlying deoxyribonucleic acid sequence. While many discoveries were first made in animal systems, plants provide a particularly complex set of epigenetic mechanisms due to unique aspects of their biology and interactions with human selective breeding and cultivation. In the plant kingdom, annual plants have received the most attention; however, perennial plants endure and respond to their environment and human management in distinct ways. Perennials include crops such as almond, for which epigenetic effects have long been linked to phenomena and even considered relevant for breeding. Recent discoveries have elucidated epigenetic phenomena that influence traits such as dormancy and self-compatibility, as well as disorders like noninfectious bud failure, which are known to be triggered by the environment and influenced by inherent aspects of the plant. Thus, epigenetics represents fertile ground to further understand almond biology and production and optimize its breeding. Here, we provide our current understanding of epigenetic regulation in plants and use almond as an example of how advances in epigenetics research can be used to understand biological fitness and agricultural performance in crop plants.

Whole-genome sequence and methylome profiling of the almond [Prunus dulcis (Mill.) D.A. Webb] cultivar 'Nonpareil'

Almond [Prunus dulcis (Mill.) D.A. Webb] is an economically important, specialty nut crop grown almost exclusively in the United States. Breeding and improvement efforts worldwide have led to the development of key, productive cultivars, including 'Nonpareil,' which is the most widely grown almond cultivar. Thus far, genomic resources for this species have been limited, and a whole-genome assembly for 'Nonpareil' is not currently available despite its economic importance and use in almond breeding worldwide. We generated a 571X coverage genome sequence using Illumina, PacBio, and optical mapping technologies. Gene prediction revealed 49,321 putative genes using MinION Oxford nanopore and Illumina RNA sequencing, and genome annotation found that 68% of predicted models are associated with at least one biological function. Furthermore, epigenetic signatures of almond, namely DNA cytosine methylation, have been implicated in a variety of phenotypes including self-compatibility, bud dormancy, and development of noninfectious bud failure. In addition to the genome sequence and annotation, this report also provides the complete methylome of several almond tissues, including leaf, flower, endocarp, mesocarp, exocarp, and seed coat. Comparisons between methylation profiles in these tissues revealed differences in genome-wide weighted % methylation and chromosome-level methylation enrichment.

Identification of Putative Markers of Non-infectious Bud Failure in Almond [Prunus dulcis (Mill.) D.A. Webb] Through Genome Wide DNA Methylation Profiling and Gene Expression Analysis in an Almond × Peach Hybrid Population

Almond [Prunus dulcis (Mill.) D.A. Webb] is an economically important nut crop susceptible to the genetic disorder, Non-infectious Bud Failure (NBF). Despite the severity of exhibition in several prominent almond cultivars, no causal mechanism has been identified underlying NBF development. The disorder is hypothesized to be associated with differential DNA methylation patterns based on patterns of inheritance (i.e., via sexual reproduction and clonal propagation) and previous work profiling methylation in affected trees. Peach (Prunus persica L. Batsch) is a closely related species that readily hybridizes with almond; however, peach is not known to exhibit NBF. A cross between an NBF-exhibiting 'Carmel' cultivar and early flowering peach ('40A17') produced an F₁ where ∼50% of progeny showed signs of NBF, including canopy die-back, erratic branching patterns (known as "crazy-top"), and rough bark. In this study, whole-genome DNA methylation profiles were generated for three F₁ progenies exhibiting NBF and three progenies considered NBF-free. Subsequent alignment to both the almond and peach reference genomes showed an increase in genome-wide methylation levels in NBF hybrids in CG and CHG contexts compared to no-NBF hybrids when aligned to the almond genome but no difference in methylation levels when aligned to the peach genome. Significantly differentially methylated regions (DMRs) were identified by comparing methylation levels across the genome between NBF- and no-NBF hybrids in each methylation context. In total, 115,635 DMRs were identified based on alignment to the almond reference genome, and 126,800 DMRs were identified based on alignment to the peach reference genome. Nearby genes were identified as associated with the 39 most significant DMRs occurring either in the almond or peach alignments alone or occurring in both the almond and peach alignments. These DMR-associated genes include several uncharacterized proteins and transposable elements. Quantitative PCR was also performed to analyze the gene expression patterns of these identified gene targets to determine patterns of differential expression associated with differential DNA methylation. These DMR-associated genes, particularly those showing corresponding patterns of differential gene expression, represent key targets for almond breeding for future cultivars and mitigating the effects of NBF-exhibition in currently affected cultivars.

Resistance to Aspergillus flavus and Aspergillus parasiticus in Almond Advanced Selections and Cultivars and Its Interaction with the Aflatoxin Biocontrol Strategy

Aflatoxin contamination of almond kernels, caused by Aspergillus flavus and A. parasiticus, is a severe concern for growers because of its high toxicity. In California, the global leader of almond production, aflatoxin can be managed by applying the biological control strain AF36 of A. flavus and selecting resistant cultivars. Here, we classified the almond genotypes by K-Means cluster analysis into three groups (susceptible [S], moderately susceptible [MS], or resistant [R]) based on aflatoxin content of inoculated kernels. The protective effects of the shell and seedcoat in preventing aflatoxin contamination were also examined. The presence of intact shells reduced aflatoxin contamination >100-fold_. The seedcoat provided a layer of protection but not complete protection. In kernel inoculation assays, none of the studied almond genotypes showed a total resistance to the pathogen. However, nine traditional cultivars and four advanced selections were classified as R. Because these advanced selections contained germplasm derived from peach, we compared the kernel resistance of three peach cultivars to that shown by kernels of an R (Sonora) and an S (Carmel) almond cultivar and five pistachio cultivars. Overall, peach kernels were significantly more resistant to the pathogen than almond kernels, which were more resistant than pistachio kernels. Finally, we studied the combined effect of the cultivar resistance and the biocontrol strain AF36 in limiting aflatoxin contamination. For this, we coinoculated almond kernels of R Sonora and S Carmel with AF36 72 h before or 48 h after inoculating with an aflatoxin-producing strain of A. flavus. The percentage of aflatoxin reduction by AF36 strain was greater in kernels of Carmel (98%) than in those of Sonora (83%). Cultivar resistance also affected the kernel colonization by the biological control strain. AF36 strain limited aflatoxin contamination in almond kernels even when applied 48 h after the aflatoxin-producing strain. Our results show that biocontrol combined with the use of cultivars with resistance to aflatoxin contamination can result in a more robust protection strategy than the use of either practice in isolation.

Advances in Rootstock Breeding of Nut Trees: Objectives and Strategies

The production and consumption of nuts are increasing in the world due to strong economic returns and the nutritional value of their products. With the increasing role and importance given to nuts (i.e., walnuts, hazelnut, pistachio, pecan, almond) in a balanced and healthy diet and their benefits to human health, breeding of the nuts species has also been stepped up. Most recent fruit breeding programs have focused on scion genetic improvement. However, the use of locally adapted grafted rootstocks also enhanced the productivity and quality of tree fruit crops. Grafting is an ancient horticultural practice used in nut crops to manipulate scion phenotype and productivity and overcome biotic and abiotic stresses. There are complex rootstock breeding objectives and physiological and molecular aspects of rootstock-scion interactions in nut crops. In this review, we provide an overview of these, considering the mechanisms involved in nutrient and water uptake, regulation of phytohormones, and rootstock influences on the scion molecular processes, including long-distance gene silencing and trans-grafting. Understanding the mechanisms resulting from rootstock × scion × environmental interactions will contribute to developing new rootstocks with resilience in the face of climate change, but also of the multitude of diseases and pests.

Relative Telomere Length and Telomerase Reverse Transcriptase (TERT) Expression Are Associated with Age in Almond (Prunus dulcis [Mill.] D.A.Webb)

While all organisms age, our understanding of how aging occurs varies among species. The aging process in perennial plants is not well-defined, yet can have implications on production and yield of valuable fruit and nut crops. Almond exhibits an age-related disorder known as non-infectious bud failure (BF) that affects vegetative bud development, indirectly affecting kernel yield. This species and disorder present an opportunity to address aging in a commercially relevant and vegetatively propagated perennial crop. The hypothesis tested in this study was that relative telomere length and/or telomerase reverse transcriptase (TERT) expression can serve as biomarkers of aging in almond. Relative telomere lengths and expression of TERT, a subunit of the enzyme telomerase, were measured via qPCR methods using bud and leaf samples collected from distinct age cohorts over a two-year period. Results from this work show a marginal but significant association between both relative telomere length and TERT expression, and age, suggesting that as almonds age, telomeres shorten and TERT expression decreases. This work provides information on potential biomarkers of perennial plant aging, contributing to our knowledge of this process. In addition, these results provide opportunities to address BF in almond breeding and nursery propagation.

Pedigree analysis of 220 almond genotypes reveals two world mainstream breeding lines based on only three different cultivars

Loss of genetic variability is an increasing challenge in tree breeding programs due to the repeated use of a reduced number of founder genotypes. However, in almond, little is known about the genetic variability in current breeding stocks, although several cases of inbreeding depression have been reported. To gain insights into the genetic structure in modern breeding programs worldwide, marker-verified pedigree data of 220 almond cultivars and breeding selections were analyzed. Inbreeding coefficients, pairwise relatedness, and genetic contribution were calculated for these genotypes. The results reveal two mainstream breeding lines based on three cultivars: "Tuono", "Cristomorto", and "Nonpareil". Descendants from "Tuono" or "Cristomorto" number 76 (sharing 34 descendants), while "Nonpareil" has 71 descendants. The mean inbreeding coefficient of the analyzed genotypes was 0.041, with 14 genotypes presenting a high inbreeding coefficient, over 0.250. Breeding programs from France, the USA, and Spain showed inbreeding coefficients of 0.075, 0.070, and 0.037, respectively. According to their genetic contribution, modern cultivars from Israel, France, the USA, Spain, and Australia trace back to a maximum of six main founding genotypes. Among the group of 65 genotypes carrying the S_f allele for self-compatibility, the mean relatedness coefficient was 0.125, with "Tuono" as the main founding genotype (24.7% of total genetic contribution). The results broaden our understanding about the tendencies followed in almond breeding over the last 50 years and will have a large impact into breeding decision-making process worldwide. Increasing current genetic variability is required in almond breeding programs to assure genetic gain and continuing breeding progress.

Redomesticating Almond to Meet Emerging Food Safety Needs

Almond is a desirable and high-quality food source where the presence of nut allergens and a vulnerability to aflatoxin and Salmonella contamination represent threats to consumer safety. In 2019, over 1 billion kg. of almonds, representing over 80% of the world total, were produced in California from a relatively few varieties with a very narrow genetic base. To address emerging needs mandated by cultural and climate changes, new germplasm has been introduced combining peach as well as wild peach and wild almond species. Advanced breeding selections incorporating exotic germplasm into a genetic background compatible with commercial production in California have demonstrated sizable reductions in level of kernel immunoreactivity as well as opportunities for improved control of aflatoxin and Salmonella. Breeding strategies employed include direct selection for reduced kernel immunoreactivity from an introgression enriched germplasm, the integration and pyramiding of resistance to multiple components of the aflatoxin disease-insect complex, and introduction of novel nut and tree traits to facilitate mechanized catch-frame field harvesting to avoid contamination with soil-borne pathogens such as Salmonella and Escherichia coli, as well as agrochemical residues.

An roGFP2-Based Bacterial Bioreporter for Redox Sensing of Plant Surfaces

The reduction-oxidation (redox) environment of the phytobiome (i.e., the plant-microbe interface) can strongly influence the outcome of the interaction between microbial pathogens, commensals, and their host. We describe a noninvasive method using a bacterial bioreporter that responds to reactive oxygen species and redox-active chemicals to compare microenvironments perceived by microbes during their initial encounter of the plant surface. A redox-sensitive variant of green fluorescent protein (roGFP2), responsive to changes in intracellular levels of reduced and oxidized glutathione, was expressed under the constitutive SP6 and fruR promoters in the epiphytic bacterium Pantoea eucalypti 299R (Pe299R/roGFP2). Analyses of Pe299R/roGFP2 cells by ratiometric fluorometry showed concentration-dependent responses to several redox active chemicals, including hydrogen peroxide (H₂O₂), dithiothreitol (DTT), and menadione. Changes in intracellular redox were detected within 5 min of addition of the chemical to Pe299R/roGFP2 cells, with approximate detection limits of 25 and 6 μM for oxidation by H₂O₂ and menadione, respectively, and 10 μM for reduction by DTT. Caffeic acid, chlorogenic acid, and ascorbic acid mitigated the H₂O₂-induced oxidation of the roGFP2 bioreporter. Aqueous washes of peach and rose flower petals from young blossoms created a lower redox state in the roGFP2 bioreporter than washes from fully mature blossoms. The bioreporter also detected differences in surface washes from peach fruit at different stages of maturity and between wounded and nonwounded sites. The Pe299R/roGFP2 reporter rapidly assesses differences in redox microenvironments and provides a noninvasive tool that may complement traditional redox-sensitive chromophores and chemical analyses of cell extracts.

Comparison of Laboratory-Developed and Commercial Monoclonal Antibody-Based Sandwich Enzyme-Linked Immunosorbent Assays for Almond (Prunus dulcis) Detection and Quantification

A commercially available monoclonal antibody (mAb)-based direct sandwich enzyme-linked immunosorbent assay (ELISA) kit (BioFront Technologies, Tallahassee, Fla., U.S.A.) was compared with an in-house developed mAb 4C10-based ELISA for almond detection. The assays were comparable in sensitivity (limit of detection < 1 ppm full fat almond, limit of quantification < 5 ppm full fat almond), specificity (no cross-reactivity with 156 tested foods at a concentration of 100000 ppm whole sample), and reproducibility (intra- and interassay variability < 15% CV). The target antigens were stable and detectable in whole almond seeds subjected to autoclaving, blanching, frying, microwaving, and dry roasting. The almond recovery ranges for spiked food matrices were 84.3% to 124.6% for 4C10 ELISA and 81.2% to 127.4% for MonoTrace ELISA. The almond recovery ranges for commercial and laboratory prepared foods with declared/known almond amount were 30.9% to 161.2% for 4C10 ELISA and 38.1% to 207.6% for MonoTrace ELISA. Neither assay registered any false-positive or negative results among the tested commercial and laboratory prepared samples.

PRACTICAL APPLICATION

Ability to detect and quantify trace amounts of almonds is important for improving safety of almond sensitive consumers. Two monoclonal antibody-based ELISAs were compared for almond detection. The information is useful to food industry, regulatory agencies, scientific community, and almond consumers.

Application of a Bayesian ordinal animal model for the estimation of breeding values for the resistance to Monilinia fruticola (G.Winter) Honey in progenies of peach [Prunus persica (L.) Batsch]

Fruit brown rot caused by Monilinia spp. is the most important fungal disease of stone fruits worldwide. Several phenotyping protocols to accurately characterize and evaluate brown rot infection have been proposed; however, the outcomes from those studies have not led to consistent advances in resistance breeding programs. Breeding for disease resistance is one of the most challenging objectives for crop improvement because disease expression is tetrahedral: it is simultaneously influenced by agent, host, environment, and human management. The present study presents a strategy based on Bayesian inference to analyze a peach breeding progeny for resistance to brown rot, evaluated using a polytomous ordinal scale. A pedigree containing two sources of resistance, one from peach and the other from almond, several commercial cultivars, and two segregating populations were analyzed to estimate the narrow-sense heritability (h² ) and breeding values (EBVs) for brown rot resistance in progenies. Results show promise for genetic improvement of disease resistance and other traits characterized by strong environmental interactions.

Effects of Peach Cultivar on Enzymatic Browning Following Cell Damage from High-Pressure Processing

Peach cultivars contribute to unique product characteristics and may affect the degree of browning after high-pressure processing (HPP). Nine peach cultivars were subjected to HPP at 0, 100, and 400 MPa for 10 min. Proton nuclear magnetic resonance (¹H NMR) relaxometry, light microscopy, color, polyphenol oxidase (PPO) activity, and total phenols were evaluated. The development of enzymatic browning during refrigerated storage occurred because of damage during HPP that triggered loss of cell integrity, allowing substrates to interact with enzymes. Increasing pressure levels resulted in greater damage, as determined by shifts in transverse relaxation time (T₂) and by light micrographs. Discoloration was triggered by membrane decompartmentalization but limited by PPO activity, which was found to correlate to cultivar harvest time (early, mid, and late season). Outcomes from the microstructure, ¹H NMR ,and PPO activity evaluation were an effective means of determining membrane decompartmentalization and allowed for prediction of browning scenarios.

The Impact of Maturity Stage on Cell Membrane Integrity and Enzymatic Browning Reactions in High Pressure Processed Peaches (Prunus persica)

Fruit maturity is an important factor associated with final product quality, and it may have an effect on the level of browning in peaches that are high pressure processed (HPP). Peaches from three different maturities, as determined by firmness (M1 = 50-55 N, M2 = 35-40 N, and M3 = 15-20 N), were subjected to pressure levels at 0.1, 200, and 400 MPa for 10 min. The damage from HPP treatment results in loss of fruit integrity and the development of browning during storage. Increasing pressure levels of HPP treatment resulted in greater damage, particularly in the more mature peaches, as determined by shifts in transverse relaxation time (T2) of the vacuolar component and by light microscopy. The discoloration of peach slices of different maturities processed at the same pressure was comparable, indicating that the effect of pressure level is greater than that of maturity in the development of browning.

Genome-wide view of genetic diversity reveals paths of selection and cultivar differentiation in peach domestication

Domestication and cultivar differentiation are requisite processes for establishing cultivated crops. These processes inherently involve substantial changes in population structure, including those from artificial selection of key genes. In this study, accessions of peach (Prunus persica) and its wild relatives were analysed genome-wide to identify changes in genetic structures and gene selections associated with their differentiation. Analysis of genome-wide informative single-nucleotide polymorphism loci revealed distinct changes in genetic structures and delineations among domesticated peach and its wild relatives and among peach landraces and modern fruit (F) and modern ornamental (O-A) cultivars. Indications of distinct changes in linkage disequilibrium extension/decay and of strong population bottlenecks or inbreeding were identified. Site frequency spectrum- and extended haplotype homozygosity-based evaluation of genome-wide genetic diversities supported selective sweeps distinguishing the domesticated peach from its wild relatives and each F/O-A cluster from the landrace clusters. The regions with strong selective sweeps harboured promising candidates for genes subjected to selection. Further sequence-based evaluation further defined the candidates and revealed their characteristics. All results suggest opportunities for identifying critical genes associated with each differentiation by analysing genome-wide genetic diversity in currently established populations. This approach obviates the special development of genetic populations, which is particularly difficult for long-lived tree crops.

Discovery of non-climacteric and suppressed climacteric bud sport mutations originating from a climacteric Japanese plum cultivar (Prunus salicina Lindl.)

Japanese plums are classified as climacteric; however, some economically important cultivars selected in California produce very little ethylene and require long ripening both "on" and "off" the tree to reach eating-ripe firmness. To unravel the ripening behavior of different Japanese plum cultivars, ripening was examined in the absence (air) or in the presence of ethylene or propylene (an ethylene analog) following a treatment or not with 1-methylcyclopropene (1-MCP, an ethylene action inhibitor). Detailed physiological studies revealed for the first time three distinct ripening types in plum fruit: climacteric, suppressed-climacteric, and non-climacteric. Responding to exogenous ethylene or propylene, the slow-softening supressed-climacteric cultivars produced detectable amounts of ethylene, in contrast to the novel non-climacteric cultivar that produced no ethylene and softened extremely slowly. Genetic analysis using microsatellite markers produced identical DNA profiles for the climacteric cultivars "Santa Rosa" and "July Santa Rosa," the suppressed-climacteric cultivars "Late Santa Rosa," "Casselman," and "Roysum" and the novel non-climacteric "Sweet Miriam," as expected since historic records present most of these cultivars as bud-sport mutations derived initially from "Santa Rosa." This present study provides a novel fruit system to address the molecular basis of ripening and to develop markers that assist breeders in providing high-quality stone fruit cultivars that can remain "on-tree," increasing fruit flavor, saving harvesting costs, and potentially reducing the need for low-temperature storage during postharvest handling.

Methylation of the S f locus in almond is associated with S-RNase loss of function

Self-compatibility in almond (Prunus dulcis) is attributed to the presence of the S f haplotype, allelic to and dominant over the series of S-alleles controlling self-incompatibility. Some forms of the S f haplotype, however, are phenotypically self-incompatible even though their nucleotide sequences are identical. DNA from leaves and styles from genetically diverse almond samples was cloned and sequenced and then analyzed for changes affecting S f -RNase variants. Epigenetic changes in several cytosine residues were detected in a fragment of 4,700 bp of the 5' upstream region of all self-compatible samples of the S f -RNases, differentiating them from all self-incompatible samples of S f -RNases analyzed. This is the first report of DNA methylation in a Rosaceae species and appears to be strongly associated with inactivation of the S f allele. Results facilitate an understanding of the evolution of self-compatibility/self-incompatibility in almond and other Prunus species, and suggest novel approaches for future crop improvement.

Application of genomic and quantitative genetic tools to identify candidate resistance genes for brown rot resistance in peach

The availability of a complete peach genome assembly and three different peach genome sequences created by our group provide new opportunities for application of genomic data and can improve the power of the classical Quantitative Trait Loci (QTL) approaches to identify candidate genes for peach disease resistance. Brown rot caused by Monilinia spp., is the most important fungal disease of stone fruits worldwide. Improved levels of peach fruit rot resistance have been identified in some cultivars and advanced selections developed in the UC Davis and USDA breeding programs. Whole genome sequencing of the Pop-DF parents lead to discovery of high-quality SNP markers for QTL genome scanning in this experimental population. Pop-DF created by crossing a brown rot moderately resistant cultivar 'Dr. Davis' and a brown rot resistant introgression line, 'F8,1-42', derived from an initial almond × peach interspecific hybrid, was evaluated for brown rot resistance in fruit of harvest maturity over three seasons. Using the SNP linkage map of Pop-DF and phenotypic data collected with inoculated fruit, a genome scan for QTL identified several SNP markers associated with brown rot resistance. Two of these QTLs were placed on linkage group 1, covering a large (physical) region on chromosome 1. The genome scan for QTL and SNP effects predicted several candidate genes associated with disease resistance responses in other host-pathogen systems. Two potential candidate genes, ppa011763m and ppa026453m, may be the genes primarily responsible for M. fructicola recognition in peach, activating both PAMP-triggered immunity (PTI) and effector-triggered immunity (ETI) responses. Our results provide a foundation for further genetic dissection, marker assisted breeding for brown rot resistance, and development of peach cultivars resistant to brown rot.

Heterogeneity in the entire genome for three genotypes of peach [Prunus persica (L.) Batsch] as distinguished from sequence analysis of genomic variants

Background

Peach [Prunus persica (L.) Batsch] is an economically important fruit crop that has become a genetic-genomic model for all Prunus species in the family Rosaceae. A doubled haploid reference genome sequence length of 227.3 Mb, a narrow genetic base contrasted by a wide phenotypic variability, the generation of cultivars through hybridization with subsequent clonal propagation, and the current accessibility of many founder genotypes, as well as the pedigree of modern commercial cultivars make peach a model for the study of inter-cultivar genomic heterogeneity and its shaping by artificial selection.

Results

The quantitative genomic differences among the three genotypes studied as genomic variants, included small variants (SNPs and InDels) and structural variants (SV) (duplications, inversions and translocations). The heirloom cultivar 'Georgia Belle’ and an almond by peach introgression breeding line 'F8,1-42’ are more heterogeneous than is the modern cultivar 'Dr. Davis’ when compared to the peach reference genome ('Lovell’). A pair-wise comparison of consensus genome sequences with 'Lovell’ showed that 'F8,1-42’ and 'Georgia Belle’ were more divergent than were 'Dr. Davis’ and 'Lovell’.

Conclusions

A novel application of emerging bioinformatics tools to the analysis of ongoing genome sequencing project outputs has led to the identification of a range of genomic variants. Results can be used to delineate the genomic and phenotypic differences among peach genotypes. For crops such as fruit trees, the availability of old cultivars, breeding selections and their pedigrees, make them suitable models for the study of genome shaping by artificial selection. The findings from the study of such genomic variants can then elucidate the control of pomological traits and the characterization of metabolic pathways, thus facilitating the development of protocols for the improvement of Prunus crops.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-14-750) contains supplementary material, which is available to authorized users.

Whole genome sequencing of peach (Prunus persica L.) for SNP identification and selection

Background

The application of next generation sequencing technologies and bioinformatic scripts to identify high frequency SNPs distributed throughout the peach genome is described. Three peach genomes were sequenced using Roche 454 and Illumina/Solexa technologies to obtain long contigs for alignment to the draft 'Lovell' peach sequence as well as sufficient depth of coverage for 'in silico' SNP discovery.

Description

The sequences were aligned to the 'Lovell' peach genome released April 01, 2010 by the International Peach Genome Initiative (IPGI). 'Dr. Davis', 'F8, 1-42' and 'Georgia Belle' were sequenced to add SNPs segregating in two breeding populations, Pop DF ('Dr. Davis' × 'F8, 1-42') and Pop DG ('Dr. Davis' × 'Georgia Belle'). Roche 454 sequencing produced 980,000 total reads with 236 Mb sequence for 'Dr. Davis' and 735,000 total reads with 172 Mb sequence for 'F8, 1-42'. 84 bp × 84 bp paired end Illumina/Solexa sequences yielded 25.5, 21.4, 25.5 million sequences for 'Dr. Davis', 'F8, 1-42' and 'Georgia Belle', respectively. BWA/SAMtools were used for alignment of raw reads and SNP detection, with custom PERL scripts for SNP filtering. Velvet's Columbus module was used for sequence assembly. Comparison of aligned and overlapping sequences from both Roche 454 and Illumina/Solexa resulted in the selection of 6654 high quality SNPs for 'Dr. Davis' vs. 'F8, 1-42' and 'Georgia Belle', distributed on eight major peach genome scaffolds as defined from the 'Lovell' assembly.

Conclusion

The eight scaffolds contained about 215-225 Mb of peach genomic sequences with one SNP/~ 40,000 bases. All sequences from Roche 454 and Illumina/Solexa have been submitted to NCBI for public use in the Short Read Archive database. SNPs have been deposited in the NCBI SNP database.

A fruit quality gene map of Prunus

BACKGROUND

Prunus fruit development, growth, ripening, and senescence includes major biochemical and sensory changes in texture, color, and flavor. The genetic dissection of these complex processes has important applications in crop improvement, to facilitate maximizing and maintaining stone fruit quality from production and processing through to marketing and consumption. Here we present an integrated fruit quality gene map of Prunus containing 133 genes putatively involved in the determination of fruit texture, pigmentation, flavor, and chilling injury resistance.

RESULTS

A genetic linkage map of 211 markers was constructed for an intraspecific peach (Prunus persica) progeny population, Pop-DG, derived from a canning peach cultivar 'Dr. Davis' and a fresh market cultivar 'Georgia Belle'. The Pop-DG map covered 818 cM of the peach genome and included three morphological markers, 11 ripening candidate genes, 13 cold-responsive genes, 21 novel EST-SSRs from the ChillPeach database, 58 previously reported SSRs, 40 RAFs, 23 SRAPs, 14 IMAs, and 28 accessory markers from candidate gene amplification. The Pop-DG map was co-linear with the Prunus reference T x E map, with 39 SSR markers in common to align the maps. A further 158 markers were bin-mapped to the reference map: 59 ripening candidate genes, 50 cold-responsive genes, and 50 novel EST-SSRs from ChillPeach, with deduced locations in Pop-DG via comparative mapping. Several candidate genes and EST-SSRs co-located with previously reported major trait loci and quantitative trait loci for chilling injury symptoms in Pop-DG.

CONCLUSION

The candidate gene approach combined with bin-mapping and availability of a community-recognized reference genetic map provides an efficient means of locating genes of interest in a target genome. We highlight the co-localization of fruit quality candidate genes with previously reported fruit quality QTLs. The fruit quality gene map developed here is a valuable tool for dissecting the genetic architecture of fruit quality traits in Prunus crops.

Development of ChillPeach genomic tools and identification of cold-responsive genes in peach fruit

The ChillPeach database was developed to facilitate identification of genes controlling chilling injury (CI), a global-scale post-harvest physiological disorder in peach. It contained 7,862 high-quality ESTs (comprising 4,468 unigenes) obtained from mesocarp tissues of two full-sib progeny contrasting for CI, about 48 and 13% of which are unique to Prunus and Arabidopsis, respectively. All ESTs are in the Gateway vector to facilitate functional assessment of the genes. The data set contained several putative SNPs and 184 unigenes with high quality SSRs, of which 42% were novel to Prunus. Microarray slides containing 4,261 ChillPeach unigenes were printed and used in a pilot experiment to identify differentially expressed genes in cold-treated compared to control mesocarp tissues, and in vegetative compared to mesocarp tissues. Quantitative RT-PCR (qRT-PCR) confirmed microarray results for all 13 genes tested. The microarray and qRT-PCR analyses indicated that ChillPeach is rich in putative fruit-specific and novel cold-induced genes. A website ( http://bioinfo.ibmcp.upv.es/genomics/ChillPeachDB ) was created holding detailed information on the ChillPeach database.

Structural and transcriptional analysis of the self-incompatibility locus of almond: identification of a pollen-expressed F-box gene with haplotype-specific polymorphism

Gametophytic self-incompatibility in Rosaceae, Solanaceae, and Scrophulariaceae is controlled by the S locus, which consists of an S-RNase gene and an unidentified "pollen S" gene. An approximately 70-kb segment of the S locus of the rosaceous species almond, the S haplotype-specific region containing the S-RNase gene, was sequenced completely. This region was found to contain two pollen-expressed F-box genes that are likely candidates for pollen S genes. One of them, named SFB (S haplotype-specific F-box protein), was expressed specifically in pollen and showed a high level of S haplotype-specific sequence polymorphism, comparable to that of the S-RNases. The other is unlikely to determine the S specificity of pollen because it showed little allelic sequence polymorphism and was expressed also in pistil. Three other S haplotypes were cloned, and the pollen-expressed genes were physically mapped. In all four cases, SFBs were linked physically to the S-RNase genes and were located at the S haplotype-specific region, where recombination is believed to be suppressed, suggesting that the two genes are inherited as a unit. These features are consistent with the hypothesis that SFB is the pollen S gene. This hypothesis predicts the involvement of the ubiquitin/26S proteasome proteolytic pathway in the RNase-based gametophytic self-incompatibility system.

Characterization of the S-locus region of almond (Prunus dulcis): analysis of a somaclonal mutant and a cosmid contig for an S haplotype

Almond has a self-incompatibility system that is controlled by an S locus consisting of the S-RNase gene and an unidentified "pollen S gene." An almond cultivar "Jeffries," a somaclonal mutant of "Nonpareil" (S(c)S(d)), has a dysfunctional S(c) haplotype both in pistil and pollen. Immunoblot and genomic Southern blot analyses detected no S(c) haplotype-specific signal in Jeffries. Southern blot showed that Jeffries has an extra copy of the S(d) haplotype. These results indicate that at least two mutations had occurred to generate Jeffries: (1) deletion of the S(c) haplotype and (2) duplication of the S(d) haplotype. To analyze the extent of the deletion in Jeffries and gain insight into the physical limit of the S locus region, approximately 200 kbp of a cosmid contig for the S(c) haplotype was constructed. Genomic Southern blot analyses showed that the deletion in Jeffries extends beyond the region covered by the contig. Most cosmid end probes, except those near the S(c)-RNase gene, cross-hybridized with DNA fragments from different S haplotypes. This suggests that regions away from the S(c)-RNase gene can recombine between different S haplotypes, implying that the cosmid contig extends to the borders of the S locus.

Electrophoretic and immunological analyses of almond (Prunusdulcis l.) genotypes and hybrids

Aqueous extracts from sixty almond samples representing various genotypes and interspecies hybrids of almond, including almond-peach, were analyzed for protein and peptide content using electrophoresis, Western immunoblotting, and enzyme-linked immunosorbent assay (ELISA). Nondenaturing nondissociating polyacrylamide gel electrophoresis (NDND-PAGE) of the aqueous extracts indicated that a single major storage protein (almond major protein -- AMP or amandin) dominated the total soluble protein composition. Denaturing SDS--PAGE analyses of the aqueous extracts revealed that the AMP was mainly composed of two sets of polypeptides with estimated molecular masses in the ranges of 38--41 kDa and 20--22 kDa, regardless of the source; however, distinct variations in the intensity and electrophoretic mobility of some bands were noted between samples. In addition to AMP, several minor polypeptides were also present in all the genotypes, and variations were seen in these as well. Regardless of the genotype, AMP was recognized in Western blots by rabbit polyclonal anti-AMP antibodies, mouse monoclonal anti-AMP antibodies (mAbs), and serum IgE from patients displaying strong serum anti-almond IgE reactivity. As with protein staining results, antibody reactivity also revealed common patterns but displayed some variation between samples. An anti-AMP inhibition ELISA was used to quantify and compare aqueous extracts for various samples. All samples (n = 60) reacted in this assay with a mean +/- standard deviation (sigma n) = 0.82 +/- 0.18 when compared to reference aqueous extract from Nonpareil designated as 1.0.

Cloning and characterization of cDNAs encoding S-RNases from almond (Prunus dulcis): primary structural features and sequence diversity of the S-RNases in Rosaceae

cDNAs encoding three S-RNases of almond (Prunus dulcis), which belongs to the family Rosaceae, were cloned and sequenced. The comparison of amino acid sequences between the S-RNases of almond and those of other rosaceous species showed that the amino acid sequences of the rosaceous S-RNases are highly divergent, and intra-subfamilial similarities are higher than inter-subfamilial similarities. Twelve amino acid sequences of the rosaceous S-RNases were aligned to characterize their primary structural features. In spite of their high level of diversification, the rosaceous S-RNases were found to have five conserved regions, C1, C2, C3, C5, and RC4 which is Rosaceae-specific conserved region. Many variable sites fall into one region, named RHV. RHV is located at a similar position to that of the hypervariable region a (HVa) of the solanaceous S-RNases, and is assumed to be involved in recognizing S-specificity of pollen. On the other hand, the region corresponding to another solanaceous hypervariable region (HVb) was not variable in the rosaceous S-RNases. In the phylogenetic tree of the T2/S type RNase, the rosaceous S-RNase fall into two subfamily-specific groups (Amygdaloideae and Maloideae). The results of sequence comparisons and phylogenetic analysis imply that the present S-RNases of Rosaceae have diverged again relatively recently, after the divergence of subfamilies.

Identification of stylar RNases associated with gametophytic self-incompatibility in almond (Prunus dulcis)

Stylar proteins of 13 almond (Prunus dulcis) cultivars with known S-genotypes were surveyed by IEF and 2D-PAGE combined with immunoblot and N-terminal amino acid sequence analyses to identify S-RNases associated with gametophytic self-incompatibility (SI) in this plant species. RNase activities corresponding to Sa and Sb, two of the four S-alleles tested, were identified by IEF and RNase activity staining. The Sa-RNase band reacted with the anti-S4-serum prepared from Japanese pear (Pyrus serotina); no reaction with the antiserum was observed with the Sb-RNase band. When the Sa-RNase band was excised from an IEF gel stained for RNase activity, subjected to SDS-PAGE, and detected by immunoblotting, it appeared that this band consisted of a single protein that reacted with the anti-S4-serum with M(r) of about 28 kDa. With 2D-PAGE and silver staining of the stylar extracts, all four S-proteins could be successfully distinguished from each other in the highly basic zone of the gel. Although Sb-, Sc-, and Sd-proteins had roughly the same M(r) of about 30 kDa, the Sc-protein seemed to be slightly smaller than the Sb-protein and slightly larger than the Sd-protein. In 2D-PAGE profiles as well, the Sa-protein had M(r) of about 28 kDa, apparently smaller than the other three proteins. A bud sport, in which one of the two S-alleles of the original cultivar is impaired, was visualized as a loss of Sc-protein, which is consistent with the previous pollination study. All four S-proteins reacted with the anti-S4-serum, probably because of the differing conformations of these S-proteins in the IEF and 2D-PAGE gels. The Sa-protein in 2D-PAGE appeared to be identical to Sa-RNase in IEF; both had the same M(r) and were reactive with the anti-S4-serum. N-terminal amino acid sequence analysis of the four S-proteins revealed that they were highly homologous to each other and similar to the S-RNases of Malus, Pyrus, Scrophulariaceae, and Solanaceae. Taken together, RNases in the style are strongly suggested to be associated with the gametophytic SI of almond. This is the first report identifying and characterizing S-RNase in almond.

Transient gene expression in maize, rice, and wheat cells using an airgun apparatus

An airgun apparatus has been constructed for transient gene expression studies of monocots. This device utilizes compressed air from a commercial airgun to propel macroprojectile and DNA-coated tungsten particles. The beta-glucuronidase (GUS) reporter gene was used to monitor transient expression in three distinct cell types of maize (Zea mays), rice (Oryza sativa), and wheat (Triticum aestivum). The highest level of GUS activity in cultured maize cells was observed when distance between stopping plate and target cells was adjusted to 4.3 centimeters. Efficiency of transformation was estimated to be 4.4 x 10(-3). In a partial vacuum of 700 millimeters Hg, velocity of macroprojectile was measured at 520 meters per second with a 6% reduction in velocity at atmospheric pressure. A polyethylene film placed in the breech before firing contributed to a 12% increase in muzzle velocity. A 700 millimeters Hg level of vacuum was necessary for maximum number of transfornants. GUS expression was also detected in wheat leaf base tissue of microdissected shoot apices. High levels of transient gene expression were also observed in hard, compact embryogenic callus of rice. These results show that the airgun apparatus is a convenient, safe, and low-cost device for rapid transient gene expression studies in cereals.