Phased diploid genome assemblies and pan-genomes provide insights into the genetic history of apple domestication

Genomic and Bioinformatic Resources for Perennial Fruit Species

In the post-genomic era, data management and development of bioinformatic tools are critical for the adequate exploitation of genomics data. In this review, we address the actual situation for the subset of crops represented by the perennial fruit species. The agronomical singularity of these species compared to plant and crop model species provides significant challenges on the implementation of good practices generally not addressed in other species. Studies are usually performed over several years in non-controlled environments, usage of rootstock is common, and breeders heavily rely on vegetative propagation. A reference genome is now available for all the major species as well as many members of the economically important genera for breeding purposes. Development of pangenome for these species is beginning to gain momentum which will require a substantial effort in term of bioinformatic tool development. The available tools for genome annotation and functional analysis will also be presented.

The apple FERONIA receptor-like kinase MdMRLK2 negatively regulates Valsa canker resistance by suppressing defence responses and hypersensitive reaction

Valsa canker, caused by the fungus Valsa mali, is one of the most destructive diseases of apple trees in China and other East Asian countries. The plant receptor-like kinase FERONIA is involved in plant cell growth, development, and immunity. However, little is known about the function of FERONIA in apple defence against V. mali. In this study, we found that MdMRLK2 was highly induced by V. mali in twigs of V. mali-susceptible Malus mellana but not in those of the resistant species Malus yunnaensis. 35S:MdMRLK2 apple plants showed compromised resistance relative to wild-type (WT) plants. Further analyses indicated that 35S:MdMRLK2 apple plants had enhanced abscisic acid (ABA) levels and reduced salicylic acid (SA) levels relative to the WT on V. mali infection. MdMRLK2 overexpression also suppressed polyphenol accumulation and inhibited the activities of phenylalanine ammonia-lyase (PAL), β-1,3-glucanase (GLU), and chitinase (CHT) during V. mali infection. Moreover, MdMRLK2 interacted with MdHIR1, a hypersensitive-induced response protein, and suppressed the MdHIR1-mediated hypersensitive reaction (HR), probably by impairing MdHIR1 self-interaction. Collectively, these findings demonstrate that overexpression of MdMRLK2 compromises Valsa canker resistance, probably by (a) altering ABA and SA levels, (b) suppressing polyphenol accumulation, (c) inhibiting PAL, GLU, and CHT activities, and (d) blocking MdHIR1-mediated HR by disrupting MdHIR1 self-interaction.

Cultivated hawthorn (Crataegus pinnatifida var. major) genome sheds light on the evolution of Maleae (apple tribe)

Cultivated hawthorn (Crataegus pinnatifida var. major) is an important medicinal and edible plant with a long history of use for health protection in China. Herein, we provide a de novo chromosome-level genome sequence of the hawthorn cultivar "Qiu Jinxing." We assembled an 823.41 Mb genome encoding 40 571 genes and further anchored the 779.24 Mb sequence into 17 pseudo-chromosomes, which account for 94.64% of the assembled genome. Phylogenomic analyses revealed that cultivated hawthorn diverged from other species within the Maleae (apple tribe) at approximately 35.4 Mya. Notably, genes involved in the flavonoid and triterpenoid biosynthetic pathways have been significantly amplified in the hawthorn genome. In addition, our results indicated that the Maleae share a unique ancient tetraploidization event; however, no recent independent whole-genome duplication event was specifically detected in hawthorn. The amplification of non-specific long terminal repeat retrotransposons contributed the most to the expansion of the hawthorn genome. Furthermore, we identified two paleo-sub-genomes in extant species of Maleae and found that these two sub-genomes showed different rearrangement mechanisms. We also reconstructed the ancestral chromosomes of Rosaceae and discussed two possible paleo-polyploid origin patterns (autopolyploidization or allopolyploidization) of Maleae. Overall, our study provides an improved context for understanding the evolution of Maleae species, and this new high-quality reference genome provides a useful resource for the horticultural improvement of hawthorn.

Allele-aware chromosome-level genome assembly of Artemisia annua reveals the correlation between ADS expansion and artemisinin yield

Artemisia annua is the major natural source of artemisinin, an anti-malarial medicine commonly used worldwide. Here, we present chromosome-level haploid maps for two A. annua strains with different artemisinin contents to explore the relationships between genomic organization and artemisinin production. High-fidelity sequencing, optical mapping, and chromatin conformation capture sequencing were used to assemble the heterogeneous and repetitive genome and resolve the haplotypes of A. annua. Approximately 50,000 genes were annotated for each haplotype genome, and a triplication event that occurred approximately 58.12 million years ago was examined for the first time in this species. A total of 3,903,467-5,193,414 variants (SNPs, indels, and structural variants) were identified in the 1.5-Gb genome during pairwise comparison between haplotypes, consistent with the high heterozygosity of this species. Genomic analyses revealed a correlation between artemisinin concents and the copy number of amorpha-4,11-diene synthase genes. This correlation was further confirmed by resequencing of 36 A. annua samples with varied artemisinin contents. Circular consensus sequencing of transcripts facilitated the detection of paralog expression. Collectively, our study provides chromosome-level allele-aware genome assemblies for two A. annua strains and new insights into the biosynthesis of artemisinin and its regulation, which will contribute to conquering malaria worldwide.

The roles of WRKY transcription factors in Malus spp. and Pyrus spp

The WRKY transcription factor gene family is known to be involved in plant defense against pathogens and in tolerance to different environmental stresses at different stages of development. The response mechanisms through which these genes act can be influenced by different phytohormones as well as by many trans- and cis-acting elements, making this network an important topic for analysis, but still something complex to fully understand. According to available reports, these genes can also perform important roles in pome species (Malus spp. and Pyrus spp.) metabolism, especially in adaptation of these plants to stressful conditions. Here, we present a quick review of what is known about WRKY genes in Malus and Pyrus genomes offering a simple way to understand what is already known about this topic. We also add information connecting the evolution of these transcription factors with others that can also be found in pomes.

Assembly of complete diploid-phased chromosomes from draft genome sequences

De novo genome assembly is essential for genomic research. High-quality genomes assembled into phased pseudomolecules are challenging to produce and often contain assembly errors because of repeats, heterozygosity, or the chosen assembly strategy. Although algorithms that produce partially phased assemblies exist, haploid draft assemblies that may lack biological information remain favored because they are easier to generate and use. We developed HaploSync, a suite of tools that produces fully phased, chromosome-scale diploid genome assemblies, and performs extensive quality control to limit assembly artifacts. HaploSync scaffolds sequences from a draft diploid assembly into phased pseudomolecules guided by a genetic map and/or the genome of a closely related species. HaploSync generates a report that visualizes the relationships between current and legacy sequences, for both haplotypes, and displays their gene and marker content. This quality control helps the user identify misassemblies and guides Haplosync’s correction of scaffolding errors. Finally, HaploSync fills assembly gaps with unplaced sequences and resolves collapsed homozygous regions. In a series of plant, fungal, and animal kingdom case studies, we demonstrate that HaploSync efficiently increases the assembly contiguity of phased chromosomes, improves completeness by filling gaps, corrects scaffolding, and correctly phases highly heterozygous, complex regions.

Research Progress on Genetic Basis of Fruit Quality Traits in Apple (Malus × domestica)

Identifying the genetic variation characteristics of phenotypic traits is important for fruit tree breeding. During the long-term evolution of fruit trees, gene recombination and natural mutation have resulted in a high degree of heterozygosity. Apple (Malus × domestica Borkh.) shows strong ecological adaptability and is widely cultivated, and is among the most economically important fruit crops worldwide. However, the high level of heterozygosity and large genome of apple, in combination with its perennial life history and long juvenile phase, complicate investigation of the genetic basis of fruit quality traits. With continuing augmentation in the apple genomic resources available, in recent years important progress has been achieved in research on the genetic variation of fruit quality traits. This review focuses on summarizing recent genetic studies on apple fruit quality traits, including appearance, flavor, nutritional, ripening, and storage qualities. In addition, we discuss the mapping of quantitative trait loci, screening of molecular markers, and mining of major genes associated with fruit quality traits. The overall aim of this review is to provide valuable insights into the mechanisms of genetic variation and molecular breeding of important fruit quality traits in apple.

Genome-Wide Identification, Characterization, and Comparative Analysis of NLR Resistance Genes in Coffea spp

The largest family of disease resistance genes in plants are nucleotide-binding site leucine-rich repeat genes (NLRs). The products of these genes are responsible for recognizing avirulence proteins (Avr) of phytopathogens and triggering specific defense responses. Identifying NLRs in plant genomes with standard gene annotation software is challenging due to their multidomain nature, sequence diversity, and clustered genomic distribution. We present the results of a genome-wide scan and comparative analysis of NLR loci in three coffee species (Coffea canephora, Coffea eugenioides and their interspecific hybrid Coffea arabica). A total of 1311 non-redundant NLR loci were identified in C. arabica, 927 in C. canephora, and 1079 in C. eugenioides, of which 809, 562, and 695 are complete loci, respectively. The NLR-Annotator tool used in this study showed extremely high sensitivities and specificities (over 99%) and increased the detection of putative NLRs in the reference coffee genomes. The NLRs loci in coffee are distributed among all chromosomes and are organized mostly in clusters. The C. arabica genome presented a smaller number of NLR loci when compared to the sum of the parental genomes (C. canephora, and C. eugenioides). There are orthologous NLRs (orthogroups) shared between coffee, tomato, potato, and reference NLRs and those that are shared only among coffee species, which provides clues about the functionality and evolutionary history of these orthogroups. Phylogenetic analysis demonstrated orthologous NLRs shared between C. arabica and the parental genomes and those that were possibly lost. The NLR family members in coffee are subdivided into two main groups: TIR-NLR (TNL) and non-TNL. The non-TNLs seem to represent a repertoire of resistance genes that are important in coffee. These results will support functional studies and contribute to a more precise use of these genes for breeding disease-resistant coffee cultivars.

Plant pan-genomics: recent advances, new challenges, and roads ahead

Pan-genomics can encompass most of the genetic diversity of a species or population and has proved to be a powerful tool for studying genomic evolution and the origin and domestication of species, and for providing information for plant improvement. Plant genomics has greatly progressed because of improvements in sequencing technologies and the rapid reduction of sequencing costs. Nevertheless, pan-genomics still presents many challenges, including computationally intensive assembly methods, high costs with large numbers of samples, ineffective integration of big data, and difficulty in applying it to downstream multi-omics analysis and breeding research. In this review, we summarize the definition and recent achievements of plant pan-genomics, computational technologies used for pan-genome construction, and the applications of pan-genomes in plant genomics and molecular breeding. We also discuss challenges and perspectives for future pan-genomics studies and provide a detailed pipeline for sample selection, genome assembly and annotation, structural variation identification, and construction and application of graph-based pan-genomes. The aim is to provide important guidance for plant pan-genome research and a better understanding of the genetic basis of genome evolution, crop domestication, and phenotypic diversity for future studies.

Hard versus soft selective sweeps during domestication and improvement in soybean

Genome scans for selection can provide an efficient way to dissect the genetic basis of domestication traits and understand mechanisms of adaptation during crop evolution. Selection involving soft sweeps (simultaneous selection for multiple alleles) is probably common in plant genomes but is under-studied, and few if any studies have systematically scanned for soft sweeps in the context of crop domestication. Using genome resequencing data from 302 wild and domesticated soybean accessions, we conducted selection scans using five widely employed statistics to identify selection candidates under classical (hard) and soft sweeps. Across the genome, inferred hard sweeps are predominant in domesticated soybean landraces and improved varieties, whereas soft sweeps are more prevalent in a representative subpopulation of the wild ancestor. Six domestication-related genes, representing both hard and soft sweeps and different stages of domestication, were used as positive controls to assess the detectability of domestication-associated sweeps. Performance of various test statistics suggests that differentiation-based (FST ) methods are robust for detecting complete hard sweeps, and that LD-based strategies perform well for identifying recent/ongoing sweeps; however, none of the test statistics detected a known soft sweep we previously documented at the domestication gene Dt1. Genome scans yielded a set of 66 candidate loci that were identified by both differentiation-based and LD-based (iHH) methods; notably, this shared set overlaps with many previously identified QTLs for soybean domestication/improvement traits. Collectively, our results will help to advance genetic characterizations of soybean domestication traits and shed light on selection modes involved in adaptation in domesticated plant species.

De novo phasing resolves haplotype sequences in complex plant genomes

Genome phasing is a recently developed assembly method that separates heterozygous eukaryotic genomic regions and builds haplotype-resolved assemblies. Because differences between haplotypes are ignored in most published de novo genomes, assemblies are available as consensus genomes consisting of haplotype mixtures, thus increasing the need for genome phasing. Here, we review the operating principles and characteristics of several freely available and widely used phasing tools (TrioCanu, FALCON-Phase, and ALLHiC). An examination of downstream analyses using haplotype-resolved genome assemblies in plants indicated significant differences among haplotypes regarding chromosomal rearrangements, sequence insertions, and expression of specific alleles that contribute to the acquisition of the biological characteristics of plant species. Finally, we suggest directions to solve addressing limitations of current genome-phasing methods. This review provides insights into the current progress, limitations, and future directions of de novo genome phasing, which will enable researchers to easily access and utilize genome-phasing in studies involving highly heterozygous complex plant genomes.

A Solanum lycopersicoides reference genome facilitates insights into tomato specialized metabolism and immunity

Wild relatives of tomato are a valuable source of natural variation in tomato breeding, as many can be hybridized to the cultivated species (Solanum lycopersicum). Several, including Solanum lycopersicoides, have been crossed to S. lycopersicum for the development of ordered introgression lines (ILs), facilitating breeding for desirable traits. Despite the utility of these wild relatives and their associated ILs, few finished genome sequences have been produced to aid genetic and genomic studies. Here we report a chromosome-scale genome assembly for S. lycopersicoides LA2951, which contains 37 938 predicted protein-coding genes. With the aid of this genome assembly, we have precisely delimited the boundaries of the S. lycopersicoides introgressions in a set of S. lycopersicum cv. VF36 × LA2951 ILs. We demonstrate the usefulness of the LA2951 genome by identifying several quantitative trait loci for phenolics and carotenoids, including underlying candidate genes, and by investigating the genome organization and immunity-associated function of the clustered Pto gene family. In addition, syntenic analysis of R2R3MYB genes sheds light on the identity of the Aubergine locus underlying anthocyanin production. The genome sequence and IL map provide valuable resources for studying fruit nutrient/quality traits, pathogen resistance, and environmental stress tolerance. We present a new genome resource for the wild species S. lycopersicoides, which we use to shed light on the Aubergine locus responsible for anthocyanin production. We also provide IL boundary mappings, which facilitated identifying novel carotenoid quantitative trait loci of which one was likely driven by an uncharacterized lycopene β-cyclase whose function we demonstrate.

Root Breeding in the Post-Genomics Era: From Concept to Practice in Apple

The development of rootstocks with a high-quality dwarf-type root system is a popular research topic in the apple industry. However, the precise breeding of rootstocks is still challenging, mainly because the root system is buried deep underground, roots have a complex life cycle, and research on root architecture has progressed slowly. This paper describes ideas for the precise breeding and domestication of wild apple resources and the application of key genes. The primary goal of this research is to combine the existing rootstock resources with molecular breeding and summarize the methods of precision breeding. Here, we reviewed the existing rootstock germplasm, high-quality genome, and genetic resources available to explain how wild resources might be used in modern breeding. In particular, we proposed the 'from genotype to phenotype' theory and summarized the difficulties in future breeding processes. Lastly, the genetics governing root diversity and associated regulatory mechanisms were elaborated on to optimize the precise breeding of rootstocks.

SunUp and Sunset genomes revealed impact of particle bombardment mediated transformation and domestication history in papaya

Transgenic papaya is widely publicized for controlling papaya ringspot virus. However, the impact of particle bombardment on the genome remains unknown. The transgenic SunUp and its progenitor Sunset genomes were assembled into 351.5 and 350.3 Mb in nine chromosomes, respectively. We identified a 1.64 Mb insertion containing three transgenic insertions in SunUp chromosome 5, consisting of 52 nuclear-plastid, 21 nuclear-mitochondrial and 1 nuclear genomic fragments. A 591.9 kb fragment in chromosome 5 was translocated into the 1.64 Mb insertion. We assembled a gapless 9.8 Mb hermaphrodite-specific region of the Y^h chromosome and its 6.0 Mb X counterpart. Resequencing 86 genomes revealed three distinct groups, validating their geographic origin and breeding history. We identified 147 selective sweeps and defined the essential role of zeta-carotene desaturase in carotenoid accumulation during domestication. Our findings elucidated the impact of particle bombardment and improved our understanding of sex chromosomes and domestication to expedite papaya improvement.

Expression Characteristics in Roots, Phloem, Leaves, Flowers and Fruits of Apple circRNA

Circular RNAs (circRNAs) are covalently closed non-coding RNAs that play pivotal roles in various biological processes. However, circRNAs' roles in different tissues of apple are currently unknown. A total of 6495 unique circRNAs were identified from roots, phloem, leaves, flowers and fruits; 65.99% of them were intergenic circRNAs. Similar to other plants, tissue-specific expression was also observed for apple circRNAs; only 175 (2.69%) circRNAs were prevalently expressed in all five different tissues, while 1256, 1064, 912, 904 and 1080 circRNAs were expressed only in roots, phloem, leaves, flowers and fruit, respectively. The hosting-genes of circRNAs showed significant differences enriched in COG, GO terms or KEGG pathways in five tissues, suggesting the special functions of circRNAs in different tissues. Potential binding interactions between circRNAs and miRNAs were investigated using TargetFinder; 2989 interactions between 647 circRNAs and 192 miRNA were predicated in the present study. It also predicted that Chr00:18744403|18744580-mdm-miR160 might play an important role in the formation of flowers or in regulating the coloration of flowers, Chr10:6857496|6858910-mdm-miR168 might be involved in response to drought stress in roots, and Chr03:1226434|1277176 may absorb mdm-miR482a-3p and play a major role in disease resistance. Two circRNAs were experimentally analyzed by qRT-PCR with divergent primers, the expression levels were consistent with RNA-seq, which indicates that the RNA-seq datasets were reliable.

Phenotypic divergence between the cultivated apple (Malus domestica) and its primary wild progenitor (Malus sieversii)

An understanding of the relationship between the cultivated apple (Malus domestica) and its primary wild progenitor species (M. sieversii) not only provides an understanding of how apples have been improved in the past, but may be useful for apple improvement in the future. We measured 10 phenotypes in over 1000 unique apple accessions belonging to M. domestica and M. sieversii from Canada's Apple Biodiversity Collection. Using principal components analysis (PCA), we determined that M. domestica and M. sieversii differ significantly in phenotypic space and are nearly completely distinguishable as two separate groups. We found that M. domestica had a shorter juvenile phase than M. sieversii and that cultivated trees produced flowers and ripe fruit later than their wild progenitors. Cultivated apples were also 3.6 times heavier, 43% less acidic, and had 68% less phenolic content than wild apples. Using historical records, we found that apple breeding over the past 200 years has resulted in a trend towards apples that have higher soluble solids, are less bitter, and soften less during storage. Our results quantify the significant changes in phenotype that have taken place since apple domestication, and provide evidence that apple breeding has led to continued phenotypic divergence of the cultivated apple from its wild progenitor species.

Comparative Analysis of Transposable Elements and the Identification of Candidate Centromeric Elements in the Prunus Subgenus Cerasus and Its Relatives

The subgenus Cerasus and its relatives include many crucial economic drupe fruits and ornamental plants. Repetitive elements make up a large part of complex genomes, and some of them play an important role in gene regulation that can affect phenotypic variation. However, the variation in their genomes remains poorly understood. This work conducted a comprehensive repetitive sequence identification across the draft genomes of eight taxa of the genus Prunus, including four of the Prunus subgenus Cerasus (Prunus pseudocerasus, P. avium, P. yedoensis and P. × yedoensis) as well as congeneric species (Prunus salicina, P. armeniaca, P. dulcis and P. persica). Annotation results showed high proportions of transposable elements in their genomes, ranging from 52.28% (P. armeniaca) to 61.86% (P. pseudocerasus). The most notable differences in the contents of long terminal repeat retrotransposons (LTR-RTs) and tandem repeats (TRs) were confirmed with de novo identification based on the structure of each genome, which significantly contributed to their genome size variation, especially in P. avium and P.salicina. Sequence comparisons showed many similar LTR-RTs closely related to their phylogenetic relationships, and a highly similar monomer unit of the TR sequence was conserved among species. Additionally, the predicted centromere-associated sequence was located in centromeric regions with FISH in the 12 taxa of Prunus. It presented significantly different signal intensities, even within the diverse interindividual phenotypes for Prunus tomentosa. This study provides insight into the LTR-RT and TR variation within Prunus and increases our knowledge about its role in genome evolution.

Twenty years of plant genome sequencing: achievements and challenges

Publication of the complete genome sequence of Arabidopsis thaliana, the first plant reference genome, in December 2000 heralded the beginning of the plant genome era. Over the past 20 years reference genomes have been generated for hundreds of plant species, spanning non-vascular to flowering plants. Releasing these plant genomes has dramatically advanced studies in all disciplines of plant biology. Importantly, multiple reference-level genomes have been generated for the major crops and their progenitors, enabling the creation of pan-genomes and exploration of domestication history and natural variations that can be adopted by modern crop breeding. We summarize the progress of plant genome sequencing and the challenges of sequencing more complex plant genomes and generating pan-genomes.

The haplotype-resolved chromosome pairs of a heterozygous diploid African cassava cultivar reveal novel pan-genome and allele-specific transcriptome features

BACKGROUND

Cassava (Manihot esculenta) is an important clonally propagated food crop in tropical and subtropical regions worldwide. Genetic gain by molecular breeding has been limited, partially because cassava is a highly heterozygous crop with a repetitive and difficult-to-assemble genome.

FINDINGS

Here we demonstrate that Pacific Biosciences high-fidelity (HiFi) sequencing reads, in combination with the assembler hifiasm, produced genome assemblies at near complete haplotype resolution with higher continuity and accuracy compared to conventional long sequencing reads. We present 2 chromosome-scale haploid genomes phased with Hi-C technology for the diploid African cassava variety TME204. With consensus accuracy >QV46, contig N50 >18 Mb, BUSCO completeness of 99%, and 35k phased gene loci, it is the most accurate, continuous, complete, and haplotype-resolved cassava genome assembly so far. Ab initio gene prediction with RNA-seq data and Iso-Seq transcripts identified abundant novel gene loci, with enriched functionality related to chromatin organization, meristem development, and cell responses. During tissue development, differentially expressed transcripts of different haplotype origins were enriched for different functionality. In each tissue, 20-30% of transcripts showed allele-specific expression (ASE) differences. ASE bias was often tissue specific and inconsistent across different tissues. Direction-shifting was observed in <2% of the ASE transcripts. Despite high gene synteny, the HiFi genome assembly revealed extensive chromosome rearrangements and abundant intra-genomic and inter-genomic divergent sequences, with large structural variations mostly related to LTR retrotransposons. We use the reference-quality assemblies to build a cassava pan-genome and demonstrate its importance in representing the genetic diversity of cassava for downstream reference-guided omics analysis and breeding.

CONCLUSIONS

The phased and annotated chromosome pairs allow a systematic view of the heterozygous diploid genome organization in cassava with improved accuracy, completeness, and haplotype resolution. They will be a valuable resource for cassava breeding and research. Our study may also provide insights into developing cost-effective and efficient strategies for resolving complex genomes with high resolution, accuracy, and continuity.

Complete chloroplast genome studies of different apple varieties indicated the origin of modern cultivated apples from Malus sieversii and Malus sylvestris

Background

Apple is one of the most important temperate deciduous fruit trees worldwide, with a wide range of cultivation. In this study, we assessed the variations and phylogenetic relationships between the complete chloroplast genomes of wild and cultivated apples (Malus spp.).

Method

We obtained the complete chloroplast genomes of 24 apple varieties using next-generation sequencing technology and compared them with genomes of (downloaded from NCBI) the wild species.

Result

The chloroplast genome of Malus is highly conserved, with a genome length of 160,067-160,290 bp, and all have a double-stranded circular tetrad structure. The gene content and sequences of genomes of wild species and cultivated apple were almost the same, but several mutation hotspot regions (psbI-atpA, psbM-psbD, and ndhC-atpE) were detected in these genomes. These regions can provide valuable information for solving specific molecular markers in taxonomic research. Phylogenetic analysis revealed that Malus formed a new clade and four cultivated varieties clustered into a branch with M. sylvestris and M. sieversii, which indicated that M. sylvestris and M. sieversii were the ancestor species of the cultivated apple.

Allele-aware chromosome-scale assembly of the allopolyploid genome of hexaploid Ma bamboo (Dendrocalamus latiflorus Munro)

Dendrocalamus latiflorus Munro is a woody clumping bamboo with rapid shoot growth. Both genetic transformation and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) gene editing techniques are available for D. latiflorus, enabling reverse genetic approaches. Thus, D. latiflorus has the potential to be a model bamboo species. However, the genome sequence of D. latiflorus has remained unreported due to its polyploidy and large genome size. Here, we sequenced the D. latiflorus genome and assembled it into three allele-aware subgenomes (AABBCC), representing the largest genome of a major bamboo species. We assembled 70 allelic chromosomes (2, 737 Mb) for hexaploid D. latiflorus using both single-molecule sequencing from the Pacific Biosciences (PacBio) Sequel platform and chromosome conformation capture sequencing (Hi-C). Repetitive sequences comprised 52.65% of the D. latiflorus genome. We annotated 135 231 protein-coding genes in the genome based on transcriptomes from eight different tissues. Transcriptome sequencing using RNA-Seq and PacBio single-molecule real-time long-read isoform sequencing revealed highly differential alternative splicing (AS) between non-abortive and abortive shoots, suggesting that AS regulates the abortion rate of bamboo shoots. This high-quality hexaploid genome and comprehensive strand-specific transcriptome datasets for this Poaceae family member will pave the way for bamboo research using D. latiflorus as a model species.

Pangenomics in crop improvement-from coding structural variations to finding regulatory variants with pangenome graphs

Since the first reported crop pangenome in 2014, advances in high-throughput and cost-effective DNA sequencing technologies facilitated multiple such studies including the pangenomes of oilseed rape (Brassica napus L.), soybean [Glycine max (L.) Merr.], rice (Oryza sativa L.), wheat (Triticum aestivum L.), and barley (Hordeum vulgare L.). Compared with single-reference genomes, pangenomes provide a more accurate representation of the genetic variation present in a species. By combining the genomic data of multiple accessions, pangenomes allow for the detection and annotation of complex DNA polymorphisms such as structural variations (SVs), one of the major determinants of genetic diversity within a species. In this review we summarize the current literature on crop pangenomics, focusing on their application to find candidate SVs involved in traits of agronomic interest. We then highlight the potential of pangenomes in the discovery and functional characterization of noncoding regulatory sequences and their variations. We conclude with a summary and outlook on innovative data structures representing the complete content of plant pangenomes including annotations of coding and noncoding elements and outcomes of transcriptomic and epigenomic experiments.

Pangenomes as a Resource to Accelerate Breeding of Under-Utilised Crop Species

Pangenomes are a rich resource to examine the genomic variation observed within a species or genera, supporting population genetics studies, with applications for the improvement of crop traits. Major crop species such as maize (Zea mays), rice (Oryza sativa), Brassica (Brassica spp.), and soybean (Glycine max) have had pangenomes constructed and released, and this has led to the discovery of valuable genes associated with disease resistance and yield components. However, pangenome data are not available for many less prominent crop species that are currently under-utilised. Despite many under-utilised species being important food sources in regional populations, the scarcity of genomic data for these species hinders their improvement. Here, we assess several under-utilised crops and review the pangenome approaches that could be used to build resources for their improvement. Many of these under-utilised crops are cultivated in arid or semi-arid environments, suggesting that novel genes related to drought tolerance may be identified and used for introgression into related major crop species. In addition, we discuss how previously collected data could be used to enrich pangenome functional analysis in genome-wide association studies (GWAS) based on studies in major crops. Considering the technological advances in genome sequencing, pangenome references for under-utilised species are becoming more obtainable, offering the opportunity to identify novel genes related to agro-morphological traits in these species.

Genetic architecture and genomic predictive ability of apple quantitative traits across environments

Implementation of genomic tools is desirable to increase the efficiency of apple breeding. Recently, the multi-environment apple reference population (apple REFPOP) proved useful for rediscovering loci, estimating genomic predictive ability, and studying genotype by environment interactions (G × E). So far, only two phenological traits were investigated using the apple REFPOP, although the population may be valuable when dissecting genetic architecture and reporting predictive abilities for additional key traits in apple breeding. Here we show contrasting genetic architecture and genomic predictive abilities for 30 quantitative traits across up to six European locations using the apple REFPOP. A total of 59 stable and 277 location-specific associations were found using GWAS, 69.2% of which are novel when compared with 41 reviewed publications. Average genomic predictive abilities of 0.18-0.88 were estimated using main-effect univariate, main-effect multivariate, multi-environment univariate, and multi-environment multivariate models. The G × E accounted for up to 24% of the phenotypic variability. This most comprehensive genomic study in apple in terms of trait-environment combinations provided knowledge of trait biology and prediction models that can be readily applied for marker-assisted or genomic selection, thus facilitating increased breeding efficiency.

Expanding Gene-Editing Potential in Crop Improvement with Pangenomes

Pangenomes aim to represent the complete repertoire of the genome diversity present within a species or cohort of species, capturing the genomic structural variance between individuals. This genomic information coupled with phenotypic data can be applied to identify genes and alleles involved with abiotic stress tolerance, disease resistance, and other desirable traits. The characterisation of novel structural variants from pangenomes can support genome editing approaches such as Clustered Regularly Interspaced Short Palindromic Repeats and CRISPR associated protein Cas (CRISPR-Cas), providing functional information on gene sequences and new target sites in variant-specific genes with increased efficiency. This review discusses the application of pangenomes in genome editing and crop improvement, focusing on the potential of pangenomes to accurately identify target genes for CRISPR-Cas editing of plant genomes while avoiding adverse off-target effects. We consider the limitations of applying CRISPR-Cas editing with pangenome references and potential solutions to overcome these limitations.

Comparative chloroplast genome analyses of cultivated spinach and two wild progenitors shed light on the phylogenetic relationships and variation

Spinacia is a genus of important leafy vegetable crops worldwide and includes cultivated Spinacia oleracea and two wild progenitors, Spinacia turkestanica and Spinacia tetrandra. However, the chloroplast genomes of the two wild progenitors remain unpublished, limiting our knowledge of chloroplast genome evolution among these three Spinacia species. Here, we reported the complete chloroplast genomes of S. oleracea, S. turkestanica, and S. tetrandra obtained via Illumina sequencing. The three chloroplast genomes exhibited a typical quadripartite structure and were 150,739, 150,747, and 150,680 bp in size, respectively. Only three variants were identified between S. oleracea and S. turkestanica, whereas 690 variants were obtained between S. oleracea and S. tetrandra, strongly demonstrating the close relationship between S. turkestanica and S. oleracea. This was further supported by phylogenetic analysis. We reported a comprehensive variant dataset including 503 SNPs and 83 Indels using 85 Spinacia accessions containing 61 S. oleracea, 16 S. turkestanica, and eight S. tetrandra accessions. Thirteen S. oleracea accessions were derived through introgression from S. turkestanica that acts as the maternal parent. Together, these results provide a valuable resource for spinach breeding programs and improve our understanding of the phylogenetic relationships within Amaranthaceae.

Genome sequence and transcriptome of Sorbus pohuashanensis provides insights into population evolution and leaf sunburn response

Sorbus pohuashanensis is a potential horticulture and medicinal plant, but its genomic and genetic backgrounds remain unknown. Here, we sequenced and assembled the S. pohuashanensis (Hance) Hedl. reference genome using PacBio long reads. Based on the new reference genome, we resequenced a core collection of 22 Sorbus spp. samples, which were divided into 2 groups (G1 and G2) based on phylogenetic and PCA analyses. These phylogenetic clusters were highly consistent with their classification based on leaf shape. Natural hybridization between the G1 and G2 groups was evidenced by a sample (R21) with a highly heterozygous genotype. Nucleotide diversity (π) analysis showed that G1 had a higher diversity than G2 and that G2 originated from G1. During the evolution process, the gene families involved in photosynthesis pathways expanded and the gene families involved in energy consumption contracted. RNA-seq data suggested that flavonoid biosynthesis and heat-shock protein (HSP)-heat-shock factor (HSF) pathways play important roles in protection against sunburn. This study provides new insights into the evolution of Sorbus spp. genomes. In addition, the genomic resources, and the identified genetic variations, especially those related to stress resistance, will help future efforts to produce and breed Sorbus spp.

Plant genome sequence assembly in the era of long reads: Progress, challenges and future directions

Third-generation long-read sequencing is transforming plant genomics. Oxford Nanopore Technologies and Pacific Biosciences are offering competing long-read sequencing technologies and enable plant scientists to investigate even large and complex plant genomes. Sequencing projects can be conducted by single research groups and sequences of smaller plant genomes can be completed within days. This also resulted in an increased investigation of genomes from multiple species in large scale to address fundamental questions associated with the origin and evolution of land plants. Increased accessibility of sequencing devices and user-friendly software allows more researchers to get involved in genomics. Current challenges are accurately resolving diploid or polyploid genome sequences and better accounting for the intra-specific diversity by switching from the use of single reference genome sequences to a pangenome graph.

Combining novel technologies with interdisciplinary basic research to enhance horticultural crops

Horticultural crops mainly include fruits, vegetables, ornamental trees and flowers, and tea trees (Melaleuca alternifolia). They produce a variety of nutrients for the daily human diet in addition to the nutrition provided by staple crops, and some of them additionally possess ornamental and medicinal features. As such, horticultural crops make unique and important contributions to both food security and a colorful lifestyle. Under the current climate change scenario, the growing population and limited arable land means that agriculture, and especially horticulture, has been facing unprecedented challenges to meet the diverse demands of human daily life. Breeding horticultural crops with high quality and adaptability and establishing an effective system that combines cultivation, post-harvest handling, and sales becomes increasingly imperative for horticultural production. This review discusses characteristic and recent research highlights in horticultural crops, focusing on the breeding of quality traits and the mechanisms that underpin them. It additionally addresses challenges and potential solutions in horticultural production and post-harvest practices. Finally, we provide a prospective as to how emerging technologies can be implemented alongside interdisciplinary basic research to enhance our understanding and exploitation of horticultural crops.

Genomic Approaches for Improvement of Tropical Fruits: Fruit Quality, Shelf Life and Nutrient Content

The breeding of tropical fruit trees for improving fruit traits is complicated, due to the long juvenile phase, generation cycle, parthenocarpy, polyploidy, polyembryony, heterozygosity and biotic and abiotic factors, as well as a lack of good genomic resources. Many molecular techniques have recently evolved to assist and hasten conventional breeding efforts. Molecular markers linked to fruit development and fruit quality traits such as fruit shape, size, texture, aroma, peel and pulp colour were identified in tropical fruit crops, facilitating Marker-assisted breeding (MAB). An increase in the availability of genome sequences of tropical fruits further aided in the discovery of SNP variants/Indels, QTLs and genes that can ascertain the genetic determinants of fruit characters. Through multi-omics approaches such as genomics, transcriptomics, metabolomics and proteomics, the identification and quantification of transcripts, including non-coding RNAs, involved in sugar metabolism, fruit development and ripening, shelf life, and the biotic and abiotic stress that impacts fruit quality were made possible. Utilizing genomic assisted breeding methods such as genome wide association (GWAS), genomic selection (GS) and genetic modifications using CRISPR/Cas9 and transgenics has paved the way to studying gene function and developing cultivars with desirable fruit traits by overcoming long breeding cycles. Such comprehensive multi-omics approaches related to fruit characters in tropical fruits and their applications in breeding strategies and crop improvement are reviewed, discussed and presented here.

High-quality reference genome and annotation aids understanding of berry development for evergreen blueberry (Vaccinium darrowii)

Vaccinium darrowii Camp (2n = 2x = 24) is a native North American blueberry species and an important source of traits such as low chill requirement in commercial southern highbush blueberry breeding (Vaccinium corymbosum, 2n = 4x = 48). We present a chromosomal-scale genome of V. darrowii generated by the combination of PacBio sequencing and high throughput chromatin conformation capture (Hi-C) scaffolding technologies, yielding a total length of 1.06 Gigabases (Gb). Over 97.8% of the genome sequences are scaffolded into 24 chromosomes representing the two haplotypes. The primary haplotype assembly of V. darrowii contains 34,809 protein-coding genes. Comparison to a V. corymbosum haplotype assembly reveals high collinearity between the two genomes with small intrachromosomal rearrangements in eight chromosome pairs. With small RNA sequencing, the annotation was further expanded to include more than 200,000 small RNA loci and 638 microRNAs expressed in berry tissues. Transcriptome analysis across fruit development stages indicates that genes involved in photosynthesis are downregulated, while genes involved in flavonoid and anthocyanin biosynthesis are significantly increased at the late stage of berry ripening. A high-quality reference genome and accompanying annotation of V. darrowii is a significant new resource for assessing the evergreen blueberry contribution to the breeding of southern highbush blueberries.

Current status of structural variation studies in plants

Structural variations (SVs) including gene presence/absence variations and copy number variations are a common feature of genomes in plants and, together with single nucleotide polymorphisms and epigenetic differences, are responsible for the heritable phenotypic diversity observed within and between species. Understanding the contribution of SVs to plant phenotypic variation is important for plant breeders to assist in producing improved varieties. The low resolution of early genetic technologies and inefficient methods have previously limited our understanding of SVs in plants. However, with the rapid expansion in genomic technologies, it is possible to assess SVs with an ever-greater resolution and accuracy. Here, we review the current status of SV studies in plants, examine the roles that SVs play in phenotypic traits, compare current technologies and assess future challenges for SV studies.

The application of pangenomics and machine learning in genomic selection in plants

Genomic selection approaches have increased the speed of plant breeding, leading to growing crop yields over the last decade. However, climate change is impacting current and future yields, resulting in the need to further accelerate breeding efforts to cope with these changing conditions. Here we present approaches to accelerate plant breeding and incorporate nonadditive effects in genomic selection by applying state-of-the-art machine learning approaches. These approaches are made more powerful by the inclusion of pangenomes, which represent the entire genome content of a species. Understanding the strengths and limitations of machine learning methods, compared with more traditional genomic selection efforts, is paramount to the successful application of these methods in crop breeding. We describe examples of genomic selection and pangenome-based approaches in crop breeding, discuss machine learning-specific challenges, and highlight the potential for the application of machine learning in genomic selection. We believe that careful implementation of machine learning approaches will support crop improvement to help counter the adverse outcomes of climate change on crop production.

Insights into the effect of human civilization on Malus evolution and domestication

The evolutionary history of the Malus genus has not been well studied. In the current study, we presented genetic evidence on the origin of the Malus genus based on genome sequencing of 297 Malus accessions, revealing the genetic relationship between wild species and cultivated apples. Our results demonstrated that North American and East Asian wild species are closer to the outgroup (pear) than Central Asian species, and hybrid species including natural (separated before the Pleistocene, about 2.5 Mya) and artificial hybrids (including ornamental trees and rootstocks) are between East and Central Asian wild species. Introgressions from M. sylvestris in cultivated apples appeared to be more extensive than those from M. sieversii, whose genetic background flowed westward across Eurasia and eastward to wild species including M. prunifolia, M. × asiatica, M. × micromalus, and M. × robust. Our results suggested that the loss of ancestral gene flow from M. sieversii in cultivated apples accompanied the movement of European traders around the world since the Age of Discovery. Natural SNP variations showed that cultivated apples had higher nucleotide diversity than wild species and more unique SNPs than other apple groups. An apple ERECTA-like gene that underwent selection during domestication on 15^th chromosome was identified as a likely major determinant of fruit length and diameter, and an NB-ARC domain-containing gene was found to strongly affect anthocyanin accumulation using a genome-wide association approach. Our results provide new insights into the origin and domestication of apples and will be useful in new breeding programmes and efforts to increase fruit crop productivity.

The long noncoding RNA MdLNC499 bridges MdWRKY1 and MdERF109 function to regulate early-stage light-induced anthocyanin accumulation in apple fruit

Anthocyanin pigments contribute to plant coloration and are valuable sources of antioxidants in the human diet as components of fruits and vegetables. Their production is known to be induced by light in apple fruit (Malus domestica); however, the underlying molecular mechanism responsible for early-stage light-induced anthocyanin biosynthesis remains unclear. Here, we identified an ethylene response factor (ERF) protein, ERF109, involved in light-induced anthocyanin biosynthesis and found that it promotes coloration by directly binding to anthocyanin-related gene promoters. Promoter::β-glucuronidase reporter analysis and Hi-C sequencing showed that a long noncoding RNA, MdLNC499, located nearby MdERF109, induces the expression of MdERF109. A W-box cis-element in the MdLNC499 promoter was found to be regulated by a transcription factor, MdWRKY1. Transient expression in apple fruit and stable transformation of apple calli allowed us to reconstruct a MdWRKY1-MdLNC499-MdERF109 transcriptional cascade in which MdWRKY1 is activated by light to increase the transcription of MdLNC499, which in turn induces MdERF109. The MdERF109 protein induces the expression of anthocyanin-related genes and the accumulation of anthocyanins in the early stages of apple coloration. Our results provide a platform for better understanding the various regulatory mechanisms involved in light-induced apple fruit coloration.

Advances in genomics and genome editing for breeding next generation of fruit and nut crops

Fruit tree crops are an essential part of the food production systems and are key to achieve food and nutrition security. Genetic improvement of fruit trees by conventional breeding has been slow due to the long juvenile phase. Advancements in genomics and molecular biology have paved the way for devising novel genetic improvement tools like genome editing, which can accelerate the breeding of these perennial crops to a great extent. In this article, advancements in genomics of fruit trees covering genome sequencing, transcriptome sequencing, genome editing technologies (GET), CRISPR-Cas system based genome editing, potential applications of CRISPR-Cas9 in fruit tree crops improvement, the factors influencing the CRISPR-Cas editing efficiency and the challenges for CRISPR-Cas9 applications in fruit tree crops improvement are reviewed. Besides, base editing, a recently emerging more precise editing system, and the future perspectives of genome editing in the improvement of fruit and nut crops are covered.

Unraveling a genetic roadmap for improved taste in the domesticated apple

Although taste is an important aspect of fruit quality, an understanding of its genetic control remains elusive in apple and other fruit crops. In this study, we conducted genomic sequence analysis of 497 Malus accessions and revealed erosion of genetic diversity caused by apple breeding and possible independent domestication events of dessert and cider apples. Signatures of selection for fruit acidity and size, but not for fruit sugar content, were detected during the processes of both domestication and improvement. Furthermore, we found that single mutations in major genes affecting fruit taste, including Ma1, MdTDT, and MdSOT2, dramatically decrease malate, citrate, and sorbitol accumulation, respectively, and correspond to important domestication events. Interestingly, Ma1 was identified to have pleiotropic effects on both organic acid content and sugar:acid ratio, suggesting that it plays a vital role in determining fruit taste. Fruit taste is unlikely to have been negatively affected by linkage drag associated with selection for larger fruit that resulted from the pyramiding of multiple genes with minor effects on fruit size. Collectively, our study provides new insights into the genetic basis of fruit quality and its evolutionary roadmap during apple domestication, pinpointing several candidate genes for genetic manipulation of fruit taste in apple.

Domestication of Crop Metabolomes: Desired and Unintended Consequences

The majority of the crops and vegetables of today were domesticated from their wild progenitors within the past 12 000 years. Considerable research effort has been expended on characterizing the genes undergoing positive and negative selection during the processes of crop domestication and improvement. Many studies have also documented how the contents of a handful of metabolites have been altered during human selection, but we are only beginning to unravel the true extent of the metabolic consequences of breeding. We highlight how crop metabolomes have been wittingly or unwittingly shaped by the processes of domestication, and highlight how we can identify new targets for metabolite engineering for the purpose of de novo domestication of crop wild relatives.

Polyploidy underlies co-option and diversification of biosynthetic triterpene pathways in the apple tribe

Whole-genome duplication (WGD) plays important roles in plant evolution and function, yet little is known about how WGD underlies metabolic diversification of natural products that bear significant medicinal properties, especially in nonmodel trees. Here, we reveal how WGD laid the foundation for co-option and differentiation of medicinally important ursane triterpene pathway duplicates, generating distinct chemotypes between species and between developmental stages in the apple tribe. After generating chromosome-level assemblies of a widely cultivated loquat variety and Gillenia trifoliata, we define differentially evolved, duplicated gene pathways and date the WGD in the apple tribe at 13.5 to 27.1 Mya, much more recent than previously thought. We then functionally characterize contrasting metabolic pathways responsible for major triterpene biosynthesis in G. trifoliata and loquat, which pre- and postdate the Maleae WGD, respectively. Our work mechanistically details the metabolic diversity that arose post-WGD and provides insights into the genomic basis of medicinal properties of loquat, which has been used in both traditional and modern medicines.

Research progress of fruit color development in apple (Malus domestica Borkh.)

Apple (Malus domestica Borkh.) is one of the most widely produced and economically important fruits in temperate regions. Fruit color development in apple is a major focus for both breeders and researchers as consumers associate brightly colored red apples with ripeness and a good flavor. In recent years, great progress has been made in the research of apple fruit color development, but its development mechanism has not been systematic dissected from the aspects of genetics, transcription or environmental factors. Here, we summarize research on the coloration of apple fruit, including the development of important genomic databases to identify important genomic regions and genes, genetic and transcriptional factors that regulate pigment accumulation, environmental factors that affect anthocyanin synthesis, and the current breeding progress of red-skinned and red-fleshed apples. We describe key transcription factors, such as MYB, bHLH, and WD40, which are involved in the regulation of anthocyanin synthesis and fruit color development in apple. We also discuss the regulation of apple color by external environmental factors such as light, temperature, and water. The aim of this review is to provide insights into the molecular mechanisms underlying anthocyanin biosynthesis in apple. This information will provide significant guidance for the breeding of high-quality red-skinned and red-fleshed apple varieties.

Back to the Origins: Background and Perspectives of Grapevine Domestication

Domestication is a process of selection driven by humans, transforming wild progenitors into domesticated crops. The grapevine (Vitis vinifera L.), besides being one of the most extensively cultivated fruit trees in the world, is also a fascinating subject for evolutionary studies. The domestication process started in the Near East and the varieties obtained were successively spread and cultivated in different areas. Whether the domestication occurred only once, or whether successive domestication events occurred independently, is a highly debated mystery. Moreover, introgression events, breeding and intense trade in the Mediterranean basin have followed, in the last thousands of years, obfuscating the genetic relationships. Although a succession of studies has been carried out to explore grapevine origin and different evolution models are proposed, an overview of the topic remains pending. We review here the findings obtained in the main phylogenetic and genomic studies proposed in the last two decades, to clarify the fundamental questions regarding where, when and how many times grapevine domestication took place. Finally, we argue that the realization of the pan-genome of grapes could be a useful resource to discover and track the changes which have occurred in the genomes and to improve our understanding about the domestication.

A Reference Genome Assembly of Simmental Cattle, Bos taurus taurus

Genomics research has relied principally on the establishment and curation of a reference genome for the species. However, it is increasingly recognized that a single reference genome cannot fully describe the extent of genetic variation within many widely distributed species. Pangenome representations are based on high-quality genome assemblies of multiple individuals and intended to represent the broadest possible diversity within a species. A Bovine Pangenome Consortium (BPC) has recently been established to begin assembling genomes from more than 600 recognized breeds of cattle, together with other related species to provide information on ancestral alleles and haplotypes. Previously reported de novo genome assemblies for Angus, Brahman, Hereford, and Highland breeds of cattle are part of the initial BPC effort. The present report describes a complete single haplotype assembly at chromosome-scale for a fullblood Simmental cow from an F1 bison-cattle hybrid fetus by trio binning. Simmental cattle, also known as Fleckvieh due to their red and white spots, originated in central Europe in the 1830s as a triple-purpose breed selected for draught, meat, and dairy production. There are over 50 million Simmental cattle in the world, known today for their fast growth and beef yields. This assembly (ARS_Simm1.0) is similar in length to the other bovine assemblies at 2.86 Gb, with a scaffold N50 of 102 Mb (max scaffold 156.8 Mb) and meets or exceeds the continuity of the best Bos taurus reference assemblies to date.

How the pan-genome is changing crop genomics and improvement

Crop genomics has seen dramatic advances in recent years due to improvements in sequencing technology, assembly methods, and computational resources. These advances have led to the development of new tools to facilitate crop improvement. The study of structural variation within species and the characterization of the pan-genome has revealed extensive genome content variation among individuals within a species that is paradigm shifting to crop genomics and improvement. Here, we review advances in crop genomics and how utilization of these tools is shifting in light of pan-genomes that are becoming available for many crop species.

How the pan-genome is changing crop genomics and improvement

Crop genomics has seen dramatic advances in recent years due to improvements in sequencing technology, assembly methods, and computational resources. These advances have led to the development of new tools to facilitate crop improvement. The study of structural variation within species and the characterization of the pan-genome has revealed extensive genome content variation among individuals within a species that is paradigm shifting to crop genomics and improvement. Here, we review advances in crop genomics and how utilization of these tools is shifting in light of pan-genomes that are becoming available for many crop species.

Genomic consequences of apple improvement

The apple (Malus domestica) is one of the world's most commercially important perennial crops and its improvement has been the focus of human effort for thousands of years. Here, we genetically characterise over 1000 apple accessions from the United States Department of Agriculture (USDA) germplasm collection using over 30,000 single-nucleotide polymorphisms (SNPs). We confirm the close genetic relationship between modern apple cultivars and their primary progenitor species, Malus sieversii from Central Asia, and find that cider apples derive more of their ancestry from the European crabapple, Malus sylvestris, than do dessert apples. We determine that most of the USDA collection is a large complex pedigree: over half of the collection is interconnected by a series of first-degree relationships. In addition, 15% of the accessions have a first-degree relationship with one of the top 8 cultivars produced in the USA. With the exception of 'Honeycrisp', the top 8 cultivars are interconnected to each other via pedigree relationships. The cultivars 'Golden Delicious' and 'Red Delicious' were found to have over 60 first-degree relatives, consistent with their repeated use by apple breeders. We detected a signature of intense selection for red skin and provide evidence that breeders also selected for increased firmness. Our results suggest that Americans are eating apples largely from a single family tree and that the apple's future improvement will benefit from increased exploitation of its tremendous natural genetic diversity.