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.

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 F1 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 F1 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.

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.

Analysis of a Multi-Environment Trial for Black Raspberry (Rubus occidentalis L.) Quality Traits.

U.S. black raspberry (BR) production is currently limited by narrowly adapted, elite germplasm. An improved understanding of genetic control and the stability of pomological traits will inform the development of improved BR germplasm and cultivars. To this end, the analysis of a multiple-environment trial of a BR mapping population derived from a cross that combines wild ancestors introgressed with commercial cultivars on both sides of its pedigree has provided insights into genetic variation, genotype-by-environment interactions, quantitative trait loci (QTL), and QTL-by-environment interactions (QEI) of fruit quality traits among diverse field environments. The genetic components and stability of four fruit size traits and six fruit biochemistry traits were characterized in this mapping population following their evaluation over three years at four distinct locations representative of current U.S. BR production. This revealed relatively stable genetic control of the four fruit size traits across the tested production environments and less stable genetic control of the fruit biochemistry traits. Of the fifteen total QTL, eleven exhibited significant QEI. Closely overlapping QTL revealed the linkage of several fruit size traits: fruit mass, drupelet count, and seed fraction. These and related findings are expected to guide further genetic characterization of BR fruit quality, management of breeding germplasm, and development of improved BR cultivars for U.S. production.

Discovery of the REN11 Locus From Vitis aestivalis for Stable Resistance to Grapevine Powdery Mildew in a Family Segregating for Several Unstable and Tissue-Specific Quantitative Resistance Loci.

Race-specific resistance loci, whether having qualitative or quantitative effects, present plant-breeding challenges for phenotypic selection and deciding which loci to select or stack with other resistance loci for improved durability. Previously, resistance to grapevine powdery mildew (GPM, caused by Erysiphe necator) was predicted to be conferred by at least three race-specific loci in the mapping family B37-28 × C56-11 segregating for GPM resistance from Vitis aestivalis. In this study, 9 years of vineyard GPM disease severity ratings plus a greenhouse and laboratory assays were genetically mapped, using a rhAmpSeq core genome marker platform with 2,000 local haplotype markers. A new qualitative resistance locus, named REN11, on the chromosome (Chr) 15 was found to be effective in nearly all (11 of 12) vineyard environments on leaves, rachis, berries, and most of the time (7 of 12) stems. REN11 was independently validated in a pseudo-testcross with the grandparent source of resistance, "Tamiami." Five other loci significantly predicted GPM severity on leaves in only one or two environments, which could indicate race-specific resistance or their roles in different timepoints in epidemic progress. Loci on Chr 8 and 9 reproducibly predicted disease severity on stems but not on other tissues and had additive effects with REN11 on the stems. The rhAmpSeq local haplotype sequences published in this study for REN11 and Chr 8 and 9 stem quantitative trait locus (QTL) can be used directly for marker-assisted selection or converted to SNP assays. In screening for REN11 in a diversity panel of 20,651 vines representing the diversity of Vitis, this rhAmpSeq haplotype had a false positive rate of 0.034% or less. The effects of the other foliar resistance loci detected in this study seem too unstable for genetic improvement regardless of quantitative effect size, whether due to race specificity or other environmental variables.

Integrating genomics and multiplatform metabolomics enables metabolite quantitative trait loci detection in breeding-relevant apple germplasm.

Apple (Malus × domestica) has commercial and nutritional value, but breeding constraints of tree crops limit varietal improvement. Marker-assisted selection minimises these drawbacks, but breeders lack applications for targeting fruit phytochemicals. To understand genotype-phytochemical associations in apples, we have developed a high-throughput integration strategy for genomic and multiplatform metabolomics data. Here, 124 apple genotypes, including members of three pedigree-connected breeding families alongside diverse cultivars and wild selections, were genotyped and phenotyped. Metabolite genome-wide association studies (mGWAS) were conducted with c. 10 000 single nucleotide polymorphisms and phenotypic data acquired via LC-MS and 1 H NMR untargeted metabolomics. Putative metabolite quantitative trait loci (mQTL) were then validated via pedigree-based analyses (PBA). Using our developed method, 519, 726 and 177 putative mQTL were detected in LC-MS positive and negative ionisation modes, and NMR, respectively. mQTL were indicated on each chromosome, with hotspots on linkage groups 16 and 17. A chlorogenic acid mQTL was discovered on chromosome 17 via mGWAS and validated with a two-step PBA, enabling discovery of novel candidate gene-metabolite relationships. Complementary data from three metabolomics approaches and dual genomics analyses increased confidence in validity of compound annotation and mQTL detection. Our platform demonstrates the utility of multiomic integration to advance data-driven, phytochemical-based plant breeding.

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.

Inflorescence Development and Floral Organogenesis in Taraxacum kok-saghyz.

Rubber dandelion (Taraxacum kok-saghyz Rodin; TK) has received attention for its natural rubber content as a strategic biomaterial, and a promising, sustainable, and renewable alternative to synthetic rubber from fossil carbon sources. Extensive research on the domestication and rubber content of TK has demonstrated TK's potential in industrial applications as a relevant natural rubber and latex-producing alternative crop. However, many aspects of its biology have been neglected in published studies. For example, floral development is still poorly characterized. TK inflorescences were studied by scanning electron microscopy. Nine stages of early inflorescence development are proposed, and floral micromorphology is detailed. Individual flower primordia development starts at the periphery and proceeds centripetally in the newly-formed inflorescence meristem. Floral organogenesis begins in the outermost flowers of the capitulum, with corolla ring and androecium formation. Following, pappus primordium-forming a ring around the base of the corolla tube-and gynoecium are observed. The transition from vegetative to inflorescence meristem was observed 21 days after germination. This description of inflorescence and flower development in TK sheds light on the complex process of flowering, pollination, and reproduction. This study will be useful for genetics, breeding, systematics, and development of agronomical practices for this new rubber-producing crop.

Computational Analysis of AmpSeq Data for Targeted, High-Throughput Genotyping of Amplicons.

Amplicon sequencing (AmpSeq) is a practical, intuitive strategy with a semi-automated computational pipeline for analysis of highly multiplexed PCR-derived sequences. This genotyping platform is particularly cost-effective when multiplexing 96 or more samples with a few amplicons up to thousands of amplicons. Amplicons can target from a single nucleotide to the upper limit of the sequencing platform. The flexibility of AmpSeq's wet lab methods make it a tool of broad interest for diverse species, and AmpSeq excels in flexibility, high-throughput, low-cost, accuracy, and semi-automated analysis. Here we provide an open science framework procedure to output data out of an AmpSeq project, with an emphasis on the bioinformatics pipeline to generate SNPs, haplotypes and presence/absence variants in a set of diverse genotypes. Open-access tutorial datasets with actual data and a containerization open source software instance are provided to enable training in each of these genotyping applications. The pipelines presented here should be applicable to the analysis of various target-enriched (e.g., amplicon or sequence capture) Illumina sequence data.

Genotyping-by-sequencing of Brassica oleracea vegetables reveals unique phylogenetic patterns, population structure and domestication footprints.

Brassica oleracea forms a diverse and economically significant crop group. Improvement efforts are often hindered by limited knowledge of diversity contained within available germplasm. Here, we employ genotyping-by-sequencing to investigate a diverse panel of 85 landrace and improved B. oleracea broccoli, cauliflower, and Chinese kale entries. Ultimately, 21,680 high-quality SNPs were used to reveal a complex and admixed population structure and clarify phylogenetic relationships among B. oleracea groups. Each broccoli landrace contained, on average, 8.4 times as many unique alleles as an improved broccoli and landraces collectively represented 81% of all broccoli-specific alleles. Commercial broccoli hybrids were largely represented by a single subpopulation identified within a complex population structure. Greater allelic diversity in landrace broccoli and 96.1% of SNPs differentiating improved cauliflower from landrace cauliflower were common to the larger pool of broccoli germplasm, supporting a parallel or later development of cauliflower due to introgression events from broccoli. Chinese kale was readily distinguished by principal coordinate analysis. Genotyping was accomplished with and without reliance upon a reference genome producing 141,317 and 20,815 filtered SNPs, respectively, supporting robust SNP discovery methods in neglected or unimproved crop groups that lack a reference genome. This work clarifies the population structure, phylogeny, and domestication footprints of landrace and improved B. oleracea broccoli using many genotyping-by-sequencing markers. Additionally, a large pool of genetic diversity contained in broccoli landraces is described which may enhance future breeding efforts.

Analysis of the first Taraxacum kok-saghyz transcriptome reveals potential rubber yield related SNPs.

Taraxacum kok-saghyz (TK) is a potential alternative crop for natural rubber (NR) production, due to its high molecular weight rubber, short breeding cycle, and diverse environmental adaptation. However, improvements in rubber yield and agronomically relevant traits are still required before it can become a commercially-viable crop. An RNA-Seq based transcriptome was developed from a pool of roots from genotypes with high and low rubber yield. A total of 55,532 transcripts with lengths over 200 bp were de novo assembled. As many as 472 transcripts were significantly homologous to 49 out of 50 known plant putative rubber biosynthesis related genes. 158 transcripts were significantly differentially expressed between high rubber and low rubber genotypes. 21,036 SNPs were different in high and low rubber TK genotypes. Among these, 50 SNPs were found within 39 transcripts highly homologous to 49 publically-searched rubber biosynthesis related genes. 117 SNPs were located within 36 of the differentially expressed gene sequences. This comprehensive TK transcriptomic reference, and large set of SNPs including putative exonic markers associated with rubber related gene homologues and differentially expressed genes, provides a solid foundation for further genetic dissection of rubber related traits, comparative genomics and marker-assisted selection for the breeding of TK.

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 (h2 ) 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.

Genome-wide DNA-(de)methylation is associated with Noninfectious Bud-failure exhibition in Almond (Prunus dulcis [Mill.] D.A.Webb).

Noninfectious bud-failure (BF) remains a major threat to almond production in California, particularly with the recent rapid expansion of acreage and as more intensive cultural practices and modern cultivars are adopted. BF has been shown to be inherited in both vegetative and sexual progeny, with exhibition related to the age and propagation history of scion clonal sources. These characteristics suggest an epigenetic influence, such as the loss of juvenility mediated by DNA-(de)methylation. Various degrees of BF have been reported among cultivars as well as within sources of clonal propagation of the same cultivar. Genome-wide methylation profiles for different clones within almond genotypes were developed to examine their association with BF levels and association with the chronological time from initial propagation. The degree of BF exhibition was found to be associated with DNA-(de)methylation and clonal age, which suggests that epigenetic changes associated with ageing may be involved in the differential exhibition of BF within and among almond clones. Research is needed to investigate the potential of DNA-(de)methylation status as a predictor for BF as well as for effective strategies to improve clonal selection against age related deterioration. This is the first report of an epigenetic-related disorder threatening a major tree crop.

Genetic dissection of powdery mildew resistance in interspecific half-sib grapevine families using SNP-based maps.

Quantitative trait locus (QTL) identification in perennial fruit crops is impeded largely by their lengthy generation time, resulting in costly and labor-intensive maintenance of breeding programs. In a grapevine (genus Vitis) breeding program, although experimental families are typically unreplicated, the genetic backgrounds may contain similar progenitors previously selected due to their contribution of favorable alleles. In this study, we investigated the utility of joint QTL identification provided by analyzing half-sib families. The genetic control of powdery mildew was studied using two half-sib F1 families, namely GE0711/1009 (MN1264 × MN1214; N = 147) and GE1025 (MN1264 × MN1246; N = 125) with multiple species in their ancestry. Maternal genetic maps consisting of 1077 and 1641 single nucleotide polymorphism (SNP) markers, respectively, were constructed using a pseudo-testcross strategy. Ratings of field resistance to powdery mildew were obtained based on whole-plant evaluation of disease severity. This 2-year analysis uncovered two QTLs that were validated on a consensus map in these half-sib families with improved precision relative to the parental maps. Examination of haplotype combinations based on the two QTL regions identified strong association of haplotypes inherited from 'Seyval blanc', through MN1264, with powdery mildew resistance. This investigation also encompassed the use of microsatellite markers to establish a correlation between 206-bp (UDV-015b) and 357-bp (VViv67) fragment sizes with resistance-carrying haplotypes. Our work is one of the first reports in grapevine demonstrating the use of SNP-based maps and haplotypes for QTL identification and tagging of powdery mildew resistance in half-sib families.

A next-generation marker genotyping platform (AmpSeq) in heterozygous crops: a case study for marker-assisted selection in grapevine.

Marker-assisted selection (MAS) is often employed in crop breeding programs to accelerate and enhance cultivar development, via selection during the juvenile phase and parental selection prior to crossing. Next-generation sequencing and its derivative technologies have been used for genome-wide molecular marker discovery. To bridge the gap between marker development and MAS implementation, this study developed a novel practical strategy with a semi-automated pipeline that incorporates trait-associated single nucleotide polymorphism marker discovery, low-cost genotyping through amplicon sequencing (AmpSeq) and decision making. The results document the development of a MAS package derived from genotyping-by-sequencing using three traits (flower sex, disease resistance and acylated anthocyanins) in grapevine breeding. The vast majority of sequence reads (⩾99%) were from the targeted regions. Across 380 individuals and up to 31 amplicons sequenced in each lane of MiSeq data, most amplicons (83 to 87%) had <10% missing data, and read depth had a median of 220-244×. Several strengths of the AmpSeq platform that make this approach of broad interest in diverse crop species include accuracy, flexibility, speed, high-throughput, low-cost and easily automated analysis.

Lessons from a Phenotyping Center Revealed by the Genome-Guided Mapping of Powdery Mildew Resistance Loci.

The genomics era brought unprecedented opportunities for genetic analysis of host resistance, but it came with the challenge that accurate and reproducible phenotypes are needed so that genomic results appropriately reflect biology. Phenotyping host resistance by natural infection in the field can produce variable results due to the uncontrolled environment, uneven distribution and genetics of the pathogen, and developmentally regulated resistance among other factors. To address these challenges, we developed highly controlled, standardized methodologies for phenotyping powdery mildew resistance in the context of a phenotyping center, receiving samples of up to 140 grapevine progeny per F1 family. We applied these methodologies to F1 families segregating for REN1- or REN2-mediated resistance and validated that some but not all bioassays identified the REN1 or REN2 locus. A point-intercept method (hyphal transects) to quantify colony density objectively at 8 or 9 days postinoculation proved to be the phenotypic response most reproducibly predicted by these resistance loci. Quantitative trait locus (QTL) mapping with genotyping-by-sequencing maps defined the REN1 and REN2 loci at relatively high resolution. In the reference PN40024 genome under each QTL, nucleotide-binding site-leucine-rich repeat candidate resistance genes were identified-one gene for REN1 and two genes for REN2. The methods described here for centralized resistance phenotyping and high-resolution genetic mapping can inform strategies for breeding resistance to powdery mildews and other pathogens on diverse, highly heterozygous hosts.

Next Generation Mapping of Enological Traits in an F2 Interspecific Grapevine Hybrid Family.

In winegrapes (Vitis spp.), fruit quality traits such as berry color, total soluble solids content (SS), malic acid content (MA), and yeast assimilable nitrogen (YAN) affect fermentation or wine quality, and are important traits in selecting new hybrid winegrape cultivars. Given the high genetic diversity and heterozygosity of Vitis species and their tendency to exhibit inbreeding depression, linkage map construction and quantitative trait locus (QTL) mapping has relied on F1 families with the use of simple sequence repeat (SSR) and other markers. This study presents the construction of a genetic map by single nucleotide polymorphisms identified through genotyping-by-sequencing (GBS) technology in an F2 mapping family of 424 progeny derived from a cross between the wild species V. riparia Michx. and the interspecific hybrid winegrape cultivar, 'Seyval'. The resulting map has 1449 markers spanning 2424 cM in genetic length across 19 linkage groups, covering 95% of the genome with an average distance between markers of 1.67 cM. Compared to an SSR map previously developed for this F2 family, these results represent an improved map covering a greater portion of the genome with higher marker density. The accuracy of the map was validated using the well-studied trait berry color. QTL affecting YAN, MA and SS related traits were detected. A joint MA and SS QTL spans a region with candidate genes involved in the malate metabolism pathway. We present an analytical pipeline for calling intercross GBS markers and a high-density linkage map for a large F2 family of the highly heterozygous Vitis genus. This study serves as a model for further genetic investigations of the molecular basis of additional unique characters of North American hybrid wine cultivars and to enhance the breeding process by marker-assisted selection. The GBS protocols for identifying intercross markers developed in this study can be adapted for other heterozygous species.

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.

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.