Genomic consequences of apple improvement

Genomic prediction and genome-wide association study using combined genotypic data from different genotyping systems: application to apple fruit quality traits

With advances in next-generation sequencing technologies, various marker genotyping systems have been developed for genomics-based approaches such as genomic selection (GS) and genome-wide association study (GWAS). As new genotyping platforms are developed, data from different genotyping platforms must be combined. However, the potential use of combined data for GS and GWAS has not yet been clarified. In this study, the accuracy of genomic prediction (GP) and the detection power of GWAS increased for most fruit quality traits of apples when using combined data from different genotyping systems, Illumina Infinium single-nucleotide polymorphism array and genotyping by random amplicon sequencing-direct (GRAS-Di) systems. In addition, the GP model, which considered the inbreeding effect, further improved the accuracy of the seven fruit traits. Runs of homozygosity (ROH) islands overlapped with the significantly associated regions detected by the GWAS for several fruit traits. Breeders may have exploited these regions to select promising apples by breeders, increasing homozygosity. These results suggest that combining genotypic data from different genotyping platforms benefits the GS and GWAS of fruit quality traits in apples. Information on inbreeding could be beneficial for improving the accuracy of GS for fruit traits of apples; however, further analysis is required to elucidate the relationship between the fruit traits and inbreeding depression (e.g. decreased vigor).

The genetic basis of apple shape and size unraveled by digital phenotyping

Great diversity of shape, size, and skin color is observed among the fruits of different apple genotypes. These traits are critical for consumers and therefore interesting targets for breeding new apple varieties. However, they are difficult to phenotype and their genetic basis, especially for fruit shape and ground color, is largely unknown. We used the FruitPhenoBox to digitally phenotype 525 genotypes of the apple reference population (apple REFPOP) genotyped for 303,148 single nucleotide polymorphism (SNP) markers. From the apple images, 573 highly heritable features describing fruit shape and size as well as 17 highly heritable features for fruit skin color were extracted to explore genotype-phenotype relationships. Out of these features, seven principal components (PCs) and 16 features with the Pearson's correlation r < 0.75 (selected features) were chosen to carry out genome-wide association studies (GWAS) for fruit shape and size. Four PCs and eight selected features were used in GWAS for fruit skin color. In total, 69 SNPs scattered over all 17 apple chromosomes were significantly associated with round, conical, cylindrical, or symmetric fruit shapes and fruit size. Novel associations with major effect on round or conical fruit shapes and fruit size were identified on chromosomes 1 and 2. Additionally, 16 SNPs associated with PCs and selected features related to red overcolor as well as green and yellow ground color were found on eight chromosomes. The identified associations can be used to advance marker-assisted selection in apple fruit breeding to systematically select for desired fruit appearance.

Crop Landraces and Indigenous Varieties: A Valuable Source of Genes for Plant Breeding

Landraces and indigenous varieties comprise valuable sources of crop species diversity. Their utilization in plant breeding may lead to increased yield and enhanced quality traits, as well as resilience to various abiotic and biotic stresses. Recently, new approaches based on the rapid advancement of genomic technologies such as deciphering of pangenomes, multi-omics tools, marker-assisted selection (MAS), genome-wide association studies (GWAS), and CRISPR/Cas9 gene editing greatly facilitated the exploitation of landraces in modern plant breeding. In this paper, we present a comprehensive overview of the implementation of new genomic technologies and highlight their importance in pinpointing the genetic basis of desirable traits in landraces and indigenous varieties of annual, perennial herbaceous, and woody crop species cultivated in the Mediterranean region. The need for further employment of advanced -omic technologies to unravel the full potential of landraces and indigenous varieties underutilized genetic diversity is also indicated. Ultimately, the large amount of genomic data emerging from the investigation of landraces and indigenous varieties reveals their potential as a source of valuable genes and traits for breeding. The role of landraces and indigenous varieties in mitigating the ongoing risks posed by climate change in agriculture and food security is also highlighted.

Multi-omics analysis reveals the biosynthesis of flavonoids during the browning process of Malus sieversii explants

Malus sieversii is a precious apple germplasm resource. Browning of explants is one of the most important factors limiting the survival rate of plant tissue culture. In order to explore the molecular mechanism of the browning degree of different strains of Malus sieversii, we compared the dynamic changes of Malus sieversii and Malus robusta Rehd. during the whole browning process using a multi-group method. A total of 44 048 differentially expressed genes (DEGs) were identified by transcriptome analysis on the DNBSEQ-T7 sequencing platform. KEGG enrichment analysis showed that the DEGs were significantly enriched in the flavonoid biosynthesis pathway. In addition, metabonomic analysis showed that (-)-epicatechin, astragalin, chrysin, irigenin, isoquercitrin, naringenin, neobavaisoflavone and prunin exhibited different degrees of free radical scavenging ability in the tissue culture browning process, and their accumulation in different varieties led to differences in the browning degree among varieties. Comprehensive transcriptome and metabonomics analysis of the data related to flavonoid biosynthesis showed that PAL, 4CL, F3H, CYP73A, CHS, CHI, ANS, DFR and PGT1 were the key genes for flavonoid accumulation during browning. In addition, WGCNA analysis revealed a strong correlation between the known flavonoid structure genes and the selected transcriptional genes. Protein interaction predictions demonstrated that 19 transcription factors (7 MYBs and 12 bHLHs) and 8 flavonoid structural genes had targeted relationships. The results show that the interspecific differential expression of flavonoid genes is the key influencing factor of the difference in browning degree between Malus sieversii and Malus robusta Rehd., providing a theoretical basis for further study on the regulation of flavonoid biosynthesis.

Genomic analysis of fruit size and shape traits in apple: unveiling candidate genes through GWAS analysis

Genomic tools facilitate the efficient selection of improved genetic materials within a breeding program. Here, we focus on two apple fruit quality traits: shape and size. We utilized data from 11 fruit morphology parameters gathered across three years of harvest from 355 genotypes of the apple REFPOP collection, which serves as a representative sample of the genetic variability present in European-cultivated apples. The data were then employed for genome-wide association studies (GWAS) using the FarmCPU and the BLINK models. The analysis identified 59 SNPs associated with fruit size and shape traits (35 with FarmCPU and 45 with BLINK) responsible for 71 QTNs. These QTNs were distributed across all chromosomes except for chromosomes 10 and 15. Thirty-four QTNs, identified by 27 SNPs, were related for size traits, and 37 QTNs, identified by 26 SNPs, were related to shape attributes. The definition of the haploblocks containing the most relevant SNPs served to propose candidate genes, among them the genes of the ovate family protein MdOFP17 and MdOFP4 that were in a 9.7kb haploblock on Chromosome 11. RNA-seq data revealed low or null expression of these genes in the oblong cultivar "Skovfoged" and higher expression in the flat "Grand'mere." The Gene Ontology enrichment analysis support a role of OFPs and hormones in shape regulation. In conclusion, this comprehensive GWAS analysis of the apple REFPOP collection has revealed promising genetic markers and candidate genes associated with apple fruit shape and size attributes, providing valuable insights that could enhance the efficiency of future breeding programs.

An efficient screening system of disease-resistant genes from wild apple, Malus sieversii in response to Valsa mali pathogenic fungus

For molecular breeding of future apples, wild apple (Malus sieversii), the primary progenitor of domesticated apples, provides abundant genetic diversity and disease-resistance traits. Valsa canker (caused by the fungal pathogen Valsa mali) poses a major threat to wild apple population as well as to cultivated apple production in China. In the present study, we developed an efficient system for screening disease-resistant genes of M. sieversii in response to V. mali. An optimal agrobacterium-mediated transient transformation of M. sieversii was first used to manipulate in situ the expression of candidate genes. After that, the pathogen V. mali was inoculated on transformed leaves and stems, and 3 additional methods for slower disease courses were developed for V. mali inoculation. To identify the resistant genes, a series of experiments were performed including morphological (incidence, lesion area/length, fungal biomass), physiological (H2O2 content, malondialdehyde content), and molecular (Real-time quantitative Polymerase Chain Reaction) approaches. Using the optimized system, we identified two transcription factors with high resistance to V. mali, MsbHLH41 and MsEIL3. Furthermore, 35 and 45 downstream genes of MsbHLH41 and MsEIL3 were identified by screening the V. mali response gene database in M. sieversii, respectively. Overall, these results indicate that the disease-resistant gene screening system has a wide range of applications for identifying resistant genes and exploring their immune regulatory networks.

Leveraging Single-Cell Populations to Uncover the Genetic Basis of Complex Traits

The ease and throughput of single-cell genomics have steadily improved, and its current trajectory suggests that surveying single-cell populations will become routine. We discuss the merger of quantitative genetics with single-cell genomics and emphasize how this synergizes with advantages intrinsic to plants. Single-cell population genomics provides increased detection resolution when mapping variants that control molecular traits, including gene expression or chromatin accessibility. Additionally, single-cell population genomics reveals the cell types in which variants act and, when combined with organism-level phenotype measurements, unveils which cellular contexts impact higher-order traits. Emerging technologies, notably multiomics, can facilitate the measurement of both genetic changes and genomic traits in single cells, enabling single-cell genetic experiments. The implementation of single-cell genetics will advance the investigation of the genetic architecture of complex molecular traits and provide new experimental paradigms to study eukaryotic genetics.

A haplotype resolved chromosome-scale assembly of North American wild apple Malus fusca and comparative genomics of the fire blight Mfu10 locus

The Pacific crabapple (Malus fusca) is a wild relative of the commercial apple (Malus × domestica). With a range extending from Alaska to Northern California, M. fusca is extremely hardy and disease resistant. The species represents an untapped genetic resource for the development of new apple cultivars with enhanced stress resistance. However, gene discovery and utilization of M. fusca have been hampered by the lack of genomic resources. Here, we present a high‐quality, haplotype‐resolved, chromosome‐scale genome assembly and annotation for M. fusca. The genome was assembled using high‐fidelity long‐reads and scaffolded using genetic maps and high‐throughput chromatin conformation capture sequencing, resulting in one of the most contiguous apple genomes to date. We annotated the genome using public transcriptomic data from the same species taken from diverse plant structures and developmental stages. Using this assembly, we explored haplotypic structural variation within the genome of M. fusca, identifying thousands of large variants. We further showed high sequence co‐linearity with other domesticated and wild Malus species. Finally, we resolve a known quantitative trait locus associated with resistance to fire blight (Erwinia amylovora). Insights gained from the assembly of a reference‐quality genome of this hardy wild apple relative will be invaluable as a tool to facilitate DNA‐informed introgression breeding.

The East Asian wild apples, Malus baccata (L.) Borkh and Malus hupehensis (Pamp.) Rehder., are additional contributors to the genomes of cultivated European and Chinese varieties

The domestication process in long-lived plant perennials differs dramatically from that of annuals, with a huge amount of genetic exchange between crop and wild populations. Though apple is a major fruit crop grown worldwide, the contribution of wild apple species to the genetic makeup of the cultivated apple genome remains a topic of intense study. We used population genomics approaches to investigate the contributions of several wild apple species to European and Chinese rootstock and dessert genomes, with a focus on the extent of wild-crop gene flow. Population genetic structure inferences revealed that the East Asian wild apples, Malus baccata (L.) Borkh and M. hupehensis (Pamp.), form a single panmictic group, and that the European dessert and rootstock apples form a specific gene pool whereas the Chinese dessert and rootstock apples were a mixture of three wild gene pools, suggesting different evolutionary histories of European and Chinese apple varieties. Coalescent-based inferences and gene flow estimates indicated that M. baccata - M. hupehensis contributed to the genome of both European and Chinese cultivated apples through wild-to-crop introgressions, and not as an initial contributor as previously supposed. We also confirmed the contribution through wild-to-crop introgressions of Malus sylvestris Mill. to the cultivated apple genome. Apple tree domestication is therefore one example in woody perennials that involved gene flow from several wild species from multiple geographical areas. This study provides an example of a complex protracted process of domestication in long-lived plant perennials, and is a starting point for apple breeding programmes.

Genetic diversity, population structure, and genome-wide association analysis of ginkgo cultivars

Ginkgo biloba is an economically valuable tree worldwide. The species has nearly become extinct during the Quaternary, which has likely resulted in reduction of its genetic variability. The genetic variability is now conserved in few natural populations in China and a number of cultivars that are, however, derived from a few ancient trees, helping the species survive in China through medieval times. Despite the recent interest in ginkgo, however, detailed knowledge of its genetic diversity, conserved in cultivated trees and cultivars, has remained poor. This limits efficient conservation of its diversity as well as efficient use of the existing germplasm resources. Here we performed genotyping-by-sequencing (GBS) on 102 cultivated germplasms of ginkgo collected to explore their genetic structure, kinship, and inbreeding prediction. For the first time in ginkgo, a genome-wide association analysis study (GWAS) was used to attempt gene mapping of seed traits. The results showed that most of the germplasms did not show any obvious genetic relationship. The size of the ginkgo germplasm population expanded significantly around 1500 years ago during the Sui and Tang dynasties. Classification of seed cultivars based on a phylogenetic perspective does not support the current classification criteria based on phenotype. Twenty-four candidate genes were localized after performing GWAS on the seed traits. Overall, this study reveals the genetic basis of ginkgo seed traits and provides insights into its cultivation history. These findings will facilitate the conservation and utilization of the domesticated germplasms of this living fossil plant.

Ecological and evolutionary drivers of phenotypic and genetic variation in the European crabapple [Malus sylvestris (L.) Mill.], a wild relative of the cultivated apple

BACKGROUND AND AIMS

Studying the relationship between phenotypic and genetic variation in populations distributed across environmental gradients can help us to understand the ecological and evolutionary processes involved in population divergence. We investigated the patterns of genetic and phenotypic diversity in the European crabapple, Malus sylvestris, a wild relative of the cultivated apple (Malus domestica) that occurs naturally across Europe in areas subjected to different climatic conditions, to test for divergence among populations.

METHODS

Growth rates and traits related to carbon uptake in seedlings collected across Europe were measured in controlled conditions and associated with the genetic status of the seedlings, which was assessed using 13 microsatellite loci and the Bayesian clustering method. Isolation-by-distance, isolation-by-climate and isolation-by-adaptation patterns, which can explain genetic and phenotypic differentiation among M. sylvestris populations, were also tested.

KEY RESULTS

A total of 11.6 % of seedlings were introgressed by M. domestica, indicating that crop-wild gene flow is ongoing in Europe. The remaining seedlings (88.4 %) belonged to seven M. sylvestris populations. Significant phenotypic trait variation among M. sylvestris populations was observed. We did not observe significant isolation by adaptation; however, the significant association between genetic variation and the climate during the Last Glacial Maximum suggests that there has been local adaptation of M. sylvestris to past climates.

CONCLUSIONS

This study provides insight into the phenotypic and genetic differentiation among populations of a wild relative of the cultivated apple. This might help us to make better use of its diversity and provide options for mitigating the impact of climate change on the cultivated apple through breeding.

Evaluation of Malus Germplasm Identifies Genetic Sources of Powdery Mildew and Frogeye Leaf Spot Resistance for Apple Breeding

Apple is an important fruit crop of temperate regions. The narrow genetic base of commercially cultivated apples has resulted in its vulnerability to a large number of fungal, bacterial, and viral pathogens. Apple breeders are always seeking new sources of resistance within the cross-compatible Malus species that can be deployed into elite genetic backgrounds. We have evaluated resistance to two major fungal diseases of apples: powdery mildew and frogeye leaf spot, using a germplasm collection of 174 Malus accessions to identify novel sources of genetic resistance. In 2020 and 2021, we evaluated these accessions for the incidence and severity of powdery mildew and frogeye leaf spot diseases at Cornell AgriTech, Geneva, New York, in a partially managed orchard. The severity and incidence of powdery mildew and frogeye leaf spot, as well as weather parameters were recorded in June, July, and August. Total incidence of powdery mildew and frogeye leaf spot infections increased from 33 to 38%, and 56 to 97% in 2020 and 2021, respectively. Our analysis showed that relative humidity and precipitation correlate with powdery mildew and frogeye leaf spot susceptibility. The predictor variables with highest impact to the variability of powdery mildew were accessions and relative humidity in May. A total of 65 Malus accessions were found to be resistant to powdery mildew, and only one accession showed moderate resistance to frogeye leaf spot. Several of these accessions belong to Malus hybrid species and domesticated apples and can therefore be potential sources of novel resistance alleles for apple breeding.

Alternative Splicing Analysis Revealed the Role of Alpha-Linolenic Acid and Carotenoids in Fruit Development of Osmanthus fragrans

Alternative splicing refers to the process of producing different splicing isoforms from the same pre-mRNA through different alternative splicing events, which almost participates in all stages of plant growth and development. In order to understand its role in the fruit development of Osmanthus fragrans, transcriptome sequencing and alternative splicing analysis was carried out on three stages of O. fragrans fruit (O. fragrans "Zi Yingui"). The results showed that the proportion of skipping exon events was the highest in all three periods, followed by a retained intron, and the proportion of mutually exclusive exon events was the lowest and most of the alternative splicing events occurred in the first two periods. The results of enrichment analysis of differentially expressed genes and differentially expressed isoforms showed that alpha-Linolenic acid metabolism, flavonoid biosynthesis, carotenoid biosynthesis, photosynthesis, and photosynthetic-antenna protein pathways were significantly enriched, which may play an important role in the fruit development of O. fragrans. The results of this study lay the foundation for further study of the development and maturation of O. fragrans fruit and further ideas for controlling fruit color and improving fruit quality and appearance.

Quantitative Traits of Interest in Apple Breeding and Their Implications for Selection

Apple breeding is a laborious and long-lasting process that requires qualified resources, land, time, and funds. In this study, more than 5000 F₁ apple hybrids from direct and testcrosses were analyzed. The results revealed how the phenotypic expression of the main quantitative traits of interest assessed in five half-sib families was controlled by the additive genetic effects and by non-additive effects of dominance and epistasis. The statistical number of hybrids required to ensure efficient selection increased exponentially with the number of desirable traits. The minimum number of progenies required to obtain a hybrid with associated quantitative traits of agronomic interest was highly variable. For two independent traits essential in selection (fruit size and quality), but incorporated together in the same hybrid, the statistical number was between about 30 and 300. If three more cumulative traits were added (a large number of fruits per tree, resistance/tolerance to apple scab, and powdery mildew attack), the limits increased to between 1500 and 18,000. The study highlighted the need for new apple varieties due to the narrowing of the genetic diversity of the cultivated species and how the choice of parents used in hybridizations (as well as the objectives pursued in the selection) can increase the efficiency of apple breeding.

The genome sequence of the apple, Malus domestica (Suckow) Borkh., 1803

We present genome assemblies from four Malus domestica cultivars (the apple; Streptophyta; Magnoliopsida; Rosales; Rosaceae). The genome sequences are 643–653 megabases in span. The greater part of each assembly length (99.24–99.74%) is scaffolded into 17 chromosomal pseudomolecules. The mitochondrial and plastid genomes were also assembled and are 400 kilobases and 167 kilobases in length respectively.

A single QTL harboring multiple genetic variations leads to complicated phenotypic segregation in apple flesh firmness and crispness

Within a QTL, the genetic recombination and interactions among five and two functional variations at MdbHLH25 and MdWDR5A caused much complicated phenotype segregation in apple FFR and FCR. The storability of climacteric fruit like apple is a quantitative trait. We previously identified 62 quantitative trait loci (QTLs) associating flesh firmness retainability (FFR) and flesh crispness retainability (FCR), but only a few functional genetic variations were identified and validated. The genetic variation network controlling fruit storability is far to be understood and diagnostic markers are needed for molecular breeding. We previously identified overlapped QTLs F16.1/H16.2 for FFR and FCR using an F1 population derived from 'Zisai Pearl' × 'Red Fuji'. In this study, five and two single-nucleotide polymorphisms (SNPs) were identified on the candidate genes MdbHLH25 and MdWDR5A within the QTL region. The SNP1 A allele at MdbHLH25 promoter reduced the expression and SNP2 T allele and/or SNP4/5 GT alleles at the exons attenuated the function of MdbHLH25 by downregulating the expression of the target genes MdACS1, which in turn led to a reduction in ethylene production and maintenance of higher flesh crispness. The SNPs did not alter the protein-protein interaction between MdbHLH25 and MdWDR5A. The joint effect of SNP genotype combinations by the SNPs on MdbHLH25 (SNP1, SNP2, and SNP4) and MdWDR5A (SNPi and SNPii) led to a much broad spectrum of phenotypic segregation in FFR and FCR. Together, the dissection of these genetic variations contributes to understanding the complicated effects of a QTL and provides good potential for marker development in molecular breeding.

DNA profiling with the 20K apple SNP array reveals Malus domestica hybridization and admixture in M. sieversii, M. orientalis, and M. sylvestris genebank accessions

The USDA-ARS National Plant Germplasm System (NPGS) apple collection in Geneva, NY, USA maintains accessions of the primary Malus domestica (Suckow) Borkh. progenitor species M. sieversii (Ledeb.) M. Roem., M. orientalis Uglitzk., and M. sylvestris (L.) Mill. Many of these accessions originated from seeds that were collected from wild populations in the species' centers of diversity. Some of these accessions have fruit phenotypes that suggest recent M. domestica hybridization, which if true would represent crop contamination of wild species populations and mislabeled species status of NPGS accessions. Pedigree connections and admixture between M. domestica and its progenitor species can be readily identified with apple SNP array data, despite such arrays not being designed for these purposes. To investigate species purity, most (463 accessions) of the NPGS accessions labeled as these three progenitor species were genotyped using the 20K apple SNP array. DNA profiles obtained were compared with a dataset of more than 5000 unique M. domestica apple cultivars. Only 212 accessions (151 M. sieversii, 26 M. orientalis, and 35 M. sylvestris) were identified as "pure" species representatives because their DNA profiles did not exhibit genotypic signatures of recent hybridization with M. domestica. Twenty-one accessions (17 M. sieversii, 1 M. orientalis, and 3 M. sylvestris) previously labeled as wild species were instead fully M. domestica. Previously unrealized hybridization and admixture between wild species and M. domestica was identified in 230 accessions (215 M. sieversii, 9 M. orientalis, and 6 M. sylvestris). Among these species-mislabeled accessions, 'Alexander', 'Gold Reinette', 'Charlamoff', 'Rosmarina Bianca', and 'King of the Pippins' were the most frequently detected M. domestica parents or grandparents. These results have implications for collection management, including germplasm distribution, and might affect conclusions of previous research focused on these three progenitor species in the NPGS apple collection. Specifically, accessions received from the NPGS for breeding and genomics, genetics, and evolutionary biology research might not be truly representative of their previously assigned species.

Genetic structure and geneflow of Malus across the Korean Peninsula using genotyping-by-sequencing

This study was to understand the genetic structure and diversity of the Korean Malus species. We used genotyping-by-sequencing (GBS) technology to analyze samples of 112 individuals belonging to 18 populations of wild Malus spp. Using GBS, we identified thousands of single nucleotide polymorphisms in the species analyzed. M. baccata and M. toringo, two dominant mainland species of the Korean Peninsula, were distinguishable based on their genetic structure. However, M. toringo collected from Jeju Island exhibited a different genetic profile than that from the mainland. We identified M. cf. micromalus as a hybrid resulting from the Jeju Island M. toringo (pollen donor) and the mainland M. baccata, (pollen recipient). Putative M. mandshurica distributed on the Korean Peninsula showed a high structural and genetic similarity with M. baccata, indicating that it might be an ecotype. Overall, this study contributes to the understanding of the population history and genetic structure of Malus in the Korean Peninsula.

A phased, chromosome-scale genome of 'Honeycrisp' apple (Malus domestica)

The apple cultivar 'Honeycrisp' has superior fruit quality traits, cold hardiness, and disease resistance, making it a popular breeding parent. However, it suffers from several physiological disorders, production, and postharvest issues. Despite several available apple genome sequences, understanding of the genetic mechanisms underlying cultivar-specific traits remains lacking. Here, we present a highly contiguous, fully phased, chromosome-level genome of 'Honeycrisp' apples, using PacBio HiFi, Omni-C, and Illumina sequencing platforms, with two assembled haplomes of 674 Mbp and 660 Mbp, and contig N50 values of 32.8 Mbp and 31.6 Mbp, respectively. Overall, 47,563 and 48,655 protein-coding genes were annotated from each haplome, capturing 96.8-97.4% complete BUSCOs in the eudicot database. Gene family analysis reveals most 'Honeycrisp' genes are assigned into orthogroups shared with other genomes, with 121 'Honeycrisp'-specific orthogroups. This resource is valuable for understanding the genetic basis of important traits in apples and related Rosaceae species to enhance breeding efforts.

Exploring epigenetic variation for breeding climate resilient apple crops

Climate change with warmer winter and spring temperatures poses major challenges to apple fruit production. Long-term observations confirm the trend toward earlier flowering, which leads to an increased risk of frost damage. New breeding strategies are needed to generate cultivars that are able to stay largely unaffected by warmer temperatures. Recently, epigenetic variation has been proposed as a new resource for breeding purposes and seems suitable in principle for apple breeding. However, to serve as a new resource for apple breeding, it is necessary to clarify whether epigenetic variation can be induced by the environment, whether it can create phenotypic variation, and whether this variation is stable across generations. In this brief review, we summarize the impact of climate change on the timing of apple phenology, highlight how epigenetic variation can potentially support novel breeding strategies, and point out important features of epigenetic variation that are required for its application in breeding programs.

Multi-omics analyses reveal MdMYB10 hypermethylation being responsible for a bud sport of apple fruit color

Apple bud sports offer a rich resource for clonal selection of numerous elite cultivars. The accumulation of somatic mutations as plants develop may potentially impact the emergence of bud sports. Previous studies focused on somatic mutation in the essential genes associated with bud sports. However, the rate and function of genome-wide somatic mutations that accumulate when a bud sport arises remain unclear. In this study, we identified a branch from a 10-year-old tree of the apple cultivar 'Oregon Spur II' as a bud sport. The mutant branch showed reduced red coloration on fruit skin. Using this plant material, we assembled a high-quality haplotype reference genome consisting of 649.61 Mb sequences with a contig N50 value of 2.04 Mb. We then estimated the somatic mutation rate of the apple tree to be 4.56 × 10 -8 per base per year, and further identified 253 somatic single-nucleotide polymorphisms (SNPs), including five non-synonymous SNPs, between the original type and mutant samples. Transcriptome analyses showed that 69 differentially expressed genes between the original type and mutant fruit skin were highly correlated with anthocyanin content. DNA methylation in the promoter of five anthocyanin-associated genes was increased in the mutant compared with the original type as determined using DNA methylation profiling. Among the genetic and epigenetic factors that directly and indirectly influence anthocyanin content in the mutant apple fruit skin, the hypermethylated promoter of MdMYB10 is important. This study indicated that numerous somatic mutations accumulated at the emergence of a bud sport from a genome-wide perspective, some of which contribute to the low coloration of the bud sport.

Integer programming for selecting set of informative markers in paternity inference

Background

Parentage information is fundamental to various life sciences. Recent advances in sequencing technologies have made it possible to accurately infer parentage even in non-model species. The optimization of sets of genome-wide markers is valuable for cost-effective applications but requires extremely large amounts of computation, which presses for the development of new efficient algorithms.

Results

Here, for a closed half-sib population, we generalized the process of marker loci selection as a binary integer programming problem. The proposed systematic formulation considered marker localization and the family structure of the potential parental population, resulting in an accurate assignment with a small set of markers. We also proposed an efficient heuristic approach, which effectively improved the number of markers, localization, and tolerance to missing data of the set. Applying this method to the actual genotypes of apple (Malus × domestica) germplasm, we identified a set of 34 SNP markers that distinguished 300 potential parents crossed to a particular cultivar with a greater than 99% accuracy.

Conclusions

We present a novel approach for selecting informative markers based on binary integer programming. Since the data generated by high-throughput sequencing technology far exceeds the requirement for parentage assignment, a combination of the systematic marker selection with targeted SNP genotyping, such as KASP, allows flexibly enlarging the analysis up to a scale that has been unrealistic in various species. The method developed in this study can be directly applied to unsolved large-scale problems in breeding, reproduction, and ecological research, and is expected to lead to novel knowledge in various biological fields. The implementation is available at https://github.com/SoNishiyama/IP-SIMPAT.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04801-z.

Genomic insights into positive selection during barley domestication

BACKGROUND

Cultivated barley (Hordeum vulgare) is widely used in animal feed, beverages, and foods and has become a model crop for molecular evolutionary studies. Few studies have examined the evolutionary fates of different types of genes in barley during the domestication process.

RESULTS

The rates of nonsynonymous substitution (Ka) to synonymous substitution (Ks) were calculated by comparing orthologous genes in different barley groups (wild vs. landrace and landrace vs. improved cultivar). The rates of evolution, properties, expression patterns, and diversity of positively selected genes (PSGs) and negatively selected genes (NSGs) were compared. PSGs evolved more rapidly, possessed fewer exons, and had lower GC content than NSGs; they were also shorter and had shorter intron, exon, and first exon lengths. Expression levels were lower, the tissue specificity of expression was higher, and codon usage bias was weaker for PSGs than for NSGs. Nucleotide diversity analysis revealed that PSGs have undergone a more severe genetic bottleneck than NSGs. Several candidate PSGs were involved in plant growth and development, which might make them as excellent targets for the molecular breeding of barley.

CONCLUSIONS

Our comprehensive analysis of the evolutionary, structural, and functional divergence between PSGs and NSGs in barley provides new insight into the evolutionary trajectory of barley during domestication. Our findings also aid future functional studies of PSGs in barley.

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.

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.

An Apple a Day Keeps the Doctor Away: Potential Role of miRNA 146 on Macrophages Treated with Exosomes Derived from Apples

The constant dialogue between the plant world and the animal world (including man among them) has been known since the time of Adam and Eve, where an apple was the origin of the evils of the world. Apart from Snow White—who might have something to object to when it comes to the use of apples—fruits, plants, and natural extracts have been known for millennia as remedies for human health-related ailments. In the light of such evidence, the aim of the present work was to investigate from a biological point of view the potential role of apple exosomes in inflammatory processes on human cells. To this end we isolated and characterized apple exosomes and treated human cells such as macrophages and NCTC L929 as cancer cells in order to evaluate the tumorigenic and anti-inflammatory effect of apple exomes. Microscopic and molecular biology analyses were conducted to characterize exosomes and to assess cell proliferation, death, and miRNA line, as well as gene expression and the uptake of exosomes by cells. The results confirm the absolute biological safety of exosomes and their anti-inflammatory effect, mediated mainly by miRNA146 production by M2 macrophages.

A new method to reconstruct the direction of parent-offspring duo relationships using SNP array data and its demonstration on ancient and modern cultivars in the outcrossing species malus × domestica

Unordered parent-offspring (PO) relationships are an outstanding issue in pedigree reconstruction studies. Resolution of the order of these relationships would expand the results, conclusions, and usefulness of such studies; however, no such PO order resolution (POR) tests currently exist. This study describes two such tests, demonstrated using SNP array data in the outcrossing species apple (Malus × domestica) on a PO relationship of known order ("Keepsake" as a parent of "Honeycrisp") and two PO relationships previously ordered only via provenance information. The first test, POR-1, tests whether some of the extended haplotypes deduced from homozygous SNP calls from one individual in an unordered PO duo are composed of recombinant haplotypes from accurately phased SNP genotypes from the second individual. If so, the first individual would be the offspring of the second individual, otherwise the opposite relationship would be present. The second test, POR-2, does not require phased SNP genotypes and uses similar logic as the POR-1 test, albeit in a different approach. The POR-1 and POR-2 tests determined the correct relationship between "Keepsake" and "Honeycrisp". The POR-2 test confirmed "Reinette Franche" as a parent of "Nonpareil" and "Brabant Bellefleur" as a parent of "Court Pendu Plat". The latter finding conflicted with the recorded provenance information, demonstrating the need for these tests. The successful demonstration of these tests suggests they can add insights to future pedigree reconstruction studies, though caveats, like extreme inbreeding or selfing, would need to be considered where relevant.

Unraveling Genetic Diversity Amongst European Hazelnut (Corylus avellana L.) Varieties in Turkey

European hazelnut (Corylus avellana) is a diploid (2n = 22), monecious and wind-pollinated species, extensively cultivated for its nuts. Turkey is the world-leading producer of hazelnut, supplying 70-80% of the world's export capacity. Hazelnut is mostly grown in the Black Sea Region, and maintained largely through clonal propagation. Understanding the genetic variation between hazelnut varieties, and defining variety-specific and disease resistance-associated alleles, would facilitate hazelnut breeding in Turkey. Widely grown varieties 'Karafındık' (2), 'Sarıfındık' (5), and 'Yomra' (2) were collected from Akçakoca in the west, while 'Tombul' (8), 'Çakıldak' (3), 'Mincane' (2), 'Allahverdi' (2), 'Sivri' (4), and 'Palaz' (5) were collected from Ordu and Giresun provinces in the east (numbers in parentheses indicate sample sizes for each variety). Powdery mildew resistant and susceptible hazelnut genotypes were collected from the field gene bank and heavily infected orchards in Giresun. Every individual was subjected to double digest restriction enzyme-associated DNA sequencing (ddRAD-seq) and a RADtag library was created. RADtags were aligned to the 'Tombul' reference genome, and Stacks software used to identify polymorphisms. 101 private and six common alleles from nine hazelnut varieties, four private from resistants and only one from susceptible were identified for diagnosis of either a certain hazelnut variety or powdery mildew resistance. Phylogenetic analysis and population structure calculations indicated that 'Mincane', 'Sarıfındık', 'Tombul', 'Çakıldak', and 'Palaz' were genetically close to each other; however, individuals within every varietal group were found in different sub-populations. Our findings indicated that years of clonal propagation of some preferred varieties across the Black Sea Region has resulted in admixed sub-populations and great genetic diversity within each variety. This impedes the development of a true breeding variety. For example, 'Tombul' is the most favored Turkish variety because of its high quality nuts, but an elite 'Tombul' line does not yet exist. This situation continues due to the lack of a breed protection program for commercially valuable hazelnut varieties. This study provides molecular markers suitable for establishing such a program.

Apple Ripening Is Controlled by a NAC Transcription Factor

Softening is a hallmark of ripening in fleshy fruits, and has both desirable and undesirable implications for texture and postharvest stability. Accordingly, the timing and extent of pre-harvest ripening and associated textural changes following harvest are key targets for improving fruit quality through breeding. Previously, we identified a large effect locus associated with harvest date and firmness in apple (Malus domestica) using genome-wide association studies (GWAS). Here, we present additional evidence that polymorphisms in or around a transcription factor gene, NAC18.1, may cause variation in these traits. First, we confirmed our previous findings with new phenotype and genotype data from ∼800 apple accessions. In this population, we compared a genetic marker within NAC18.1 to markers targeting three other firmness-related genes currently used by breeders (ACS1, ACO1, and PG1), and found that the NAC18.1 marker was the strongest predictor of both firmness at harvest and firmness after 3 months of cold storage. By sequencing NAC18.1 across 18 accessions, we revealed two predominant haplotypes containing the single nucleotide polymorphism (SNP) previously identified using GWAS, as well as dozens of additional SNPs and indels in both the coding and promoter sequences. NAC18.1 encodes a protein that is orthogolous to the NON-RIPENING (NOR) transcription factor, a regulator of ripening in tomato (Solanum lycopersicum). We introduced both NAC18.1 transgene haplotypes into the tomato nor mutant and showed that both haplotypes complement the nor ripening deficiency. Taken together, these results indicate that polymorphisms in NAC18.1 may underlie substantial variation in apple firmness through modulation of a conserved ripening program.

Recent Large-Scale Genotyping and Phenotyping of Plant Genetic Resources of Vegetatively Propagated Crops

Several recent national and international projects have focused on large-scale genotyping of plant genetic resources in vegetatively propagated crops like fruit and berries, potatoes and woody ornamentals. The primary goal is usually to identify true-to-type plant material, detect possible synonyms, and investigate genetic diversity and relatedness among accessions. A secondary goal may be to create sustainable databases that can be utilized in research and breeding for several years ahead. Commonly applied DNA markers (like microsatellite DNA and SNPs) and next-generation sequencing each have their pros and cons for these purposes. Methods for large-scale phenotyping have lagged behind, which is unfortunate since many commercially important traits (yield, growth habit, storability, and disease resistance) are difficult to score. Nevertheless, the analysis of gene action and development of robust DNA markers depends on environmentally controlled screening of very large sets of plant material. Although more time-consuming, co-operative projects with broad-scale data collection are likely to produce more reliable results. In this review, we will describe some of the approaches taken in genotyping and/or phenotyping projects concerning a wide variety of vegetatively propagated crops.