A high-density, multi-parental SNP genetic map on apple validates a new mapping approach for outcrossing species.

Quantitative trait loci (QTL) mapping approaches rely on the correct ordering of molecular markers along the chromosomes, which can be obtained from genetic linkage maps or a reference genome sequence. For apple (Malus domestica Borkh), the genome sequence v1 and v2 could not meet this need; therefore, a novel approach was devised to develop a dense genetic linkage map, providing the most reliable marker-loci order for the highest possible number of markers. The approach was based on four strategies: (i) the use of multiple full-sib families, (ii) the reduction of missing information through the use of HaploBlocks and alternative calling procedures for single-nucleotide polymorphism (SNP) markers, (iii) the construction of a single backcross-type data set including all families, and (iv) a two-step map generation procedure based on the sequential inclusion of markers. The map comprises 15 417 SNP markers, clustered in 3 K HaploBlock markers spanning 1 267 cM, with an average distance between adjacent markers of 0.37 cM and a maximum distance of 3.29 cM. Moreover, chromosome 5 was oriented according to its homoeologous chromosome 10. This map was useful to improve the apple genome sequence, design the Axiom Apple 480 K SNP array and perform multifamily-based QTL studies. Its collinearity with the genome sequences v1 and v3 are reported. To our knowledge, this is the shortest published SNP map in apple, while including the largest number of markers, families and individuals. This result validates our methodology, proving its value for the construction of integrated linkage maps for any outbreeding species.

High-density SNP-based genetic maps for the parents of an outcrossed and a selfed tetraploid garden rose cross, inferred from admixed progeny using the 68k rose SNP array.

Dense genetic maps create a base for QTL analysis of important traits and future implementation of marker-assisted breeding. In tetraploid rose, the existing linkage maps include <300 markers to cover 28 linkage groups (4 homologous sets of 7 chromosomes). Here we used the 68k WagRhSNP Axiom single-nucleotide polymorphism (SNP) array for rose, in combination with SNP dosage calling at the tetraploid level, to genotype offspring from the garden rose cultivar 'Red New Dawn'. The offspring proved to be not from a single bi-parental cross. In rose breeding, crosses with unintended parents occur regularly. We developed a strategy to separate progeny into putative populations, even while one of the parents was unknown, using principle component analysis on pairwise genetic distances based on sets of selected SNP markers that were homozygous, and therefore uninformative for one parent. One of the inferred populations was consistent with self-fertilization of 'Red New Dawn'. Subsequently, linkage maps were generated for a bi-parental and a self-pollinated population with 'Red New Dawn' as the common maternal parent. The densest map, for the selfed parent, had 1929 SNP markers on 25 linkage groups, covering 1765.5 cM at an average marker distance of 0.9 cM. Synteny with the strawberry (Fragaria vesca) genome was extensive. Rose ICM1 corresponded to F. vesca pseudochromosome 7 (Fv7), ICM4 to Fv4, ICM5 to Fv3, ICM6 to Fv2 and ICM7 to Fv5. Rose ICM2 corresponded to parts of F. vesca pseudochromosomes 1 and 6, whereas ICM3 is syntenic to the remainder of Fv6.

PediHaplotyper: software for consistent assignment of marker haplotypes in pedigrees.

In the study of large outbred pedigrees with many founders, individual bi-allelic markers, such as SNP markers, carry little information. After phasing the marker genotypes, multi-allelic loci consisting of groups of closely linked markers can be identified, which are called "haploblocks". Here, we describe PediHaplotyper, an R package capable of assigning consistent alleles to such haploblocks, allowing for missing and incorrect SNP data. These haploblock genotypes are much easier to interpret by the human investigator than the original SNP data and also allow more efficient QTL analyses that require less memory and computation time.

FaRXf1: a locus conferring resistance to angular leaf spot caused by Xanthomonas fragariae in octoploid strawberry.

KEY MESSAGE: Angular leaf spot is a devastating bacterial disease of strawberry. Resistance from two wild accessions is highly heritable and controlled by a major locus on linkage group 6D. Angular leaf spot caused by Xanthomonas fragariae is the only major bacterial disease of cultivated strawberry (Fragaria ×ananassa). While this disease may cause reductions of up to 8 % of marketable yield in Florida winter annual production, no resistant cultivars have been commercialized. Wild accessions US4808 and US4809 were previously identified as resistant to the four genetic clades of X. fragariae, and introgression of the trait into commercial quality perennial-type germplasm was initiated. Previous reports indicated high heritability for the trait but proposed both single-locus and multi-locus inheritance models. The objective of this study was to determine the mode of inheritance of resistance, to identify causal loci, and to begin introgression of resistance into Florida-adapted germplasm. Resistance was observed in two years of field trials with inoculated plants that assayed four full-sib families descended from US4808 to US4809. Resistance segregated 1:1 in all families indicating control by a dominant allele at a single locus. Using a selective genotyping approach with the IStraw90 Axiom(®) SNP array and pedigree-based QTL detection, a single major-effect QTL was identified in two full-sib families, one descended from each resistant accession. High-resolution melt curve analysis validated the presence of the QTL in separate populations. The QTL was delimited to the 33.1-33.6 Mbp (F. vesca vesca v1.1 reference) and 34.8-35.3 Mbp (F. vesca bracteata v2.0 reference) regions of linkage group 6D for both resistance sources and was designated FaRXf1. Characterization of this locus will facilitate marker-assisted selection toward the development of resistant cultivars.

Using RNA-Seq to assemble a rose transcriptome with more than 13,000 full-length expressed genes and to develop the WagRhSNP 68k Axiom SNP array for rose (Rosa L.).

In order to develop a versatile and large SNP array for rose, we set out to mine ESTs from diverse sets of rose germplasm. For this RNA-Seq libraries containing about 700 million reads were generated from tetraploid cut and garden roses using Illumina paired-end sequencing, and from diploid Rosa multiflora using 454 sequencing. Separate de novo assemblies were performed in order to identify single nucleotide polymorphisms (SNPs) within and between rose varieties. SNPs among tetraploid roses were selected for constructing a genotyping array that can be employed for genetic mapping and marker-trait association discovery in breeding programs based on tetraploid germplasm, both from cut roses and from garden roses. In total 68,893 SNPs were included on the WagRhSNP Axiom array. Next, an orthology-guided assembly was performed for the construction of a non-redundant rose transcriptome database. A total of 21,740 transcripts had significant hits with orthologous genes in the strawberry (Fragaria vesca L.) genome. Of these 13,390 appeared to contain the full-length coding regions. This newly established transcriptome resource adds considerably to the currently available sequence resources for the Rosaceae family in general and the genus Rosa in particular.

Accuracy and responses of genomic selection on key traits in apple breeding.

The application of genomic selection in fruit tree crops is expected to enhance breeding efficiency by increasing prediction accuracy, increasing selection intensity and decreasing generation interval. The objectives of this study were to assess the accuracy of prediction and selection response in commercial apple breeding programmes for key traits. The training population comprised 977 individuals derived from 20 pedigreed full-sib families. Historic phenotypic data were available on 10 traits related to productivity and fruit external appearance and genotypic data for 7829 SNPs obtained with an Illumina 20K SNP array. From these data, a genome-wide prediction model was built and subsequently used to calculate genomic breeding values of five application full-sib families. The application families had genotypes at 364 SNPs from a dedicated 512 SNP array, and these genotypic data were extended to the high-density level by imputation. These five families were phenotyped for 1 year and their phenotypes were compared to the predicted breeding values. Accuracy of genomic prediction across the 10 traits reached a maximum value of 0.5 and had a median value of 0.19. The accuracies were strongly affected by the phenotypic distribution and heritability of traits. In the largest family, significant selection response was observed for traits with high heritability and symmetric phenotypic distribution. Traits that showed non-significant response often had reduced and skewed phenotypic variation or low heritability. Among the five application families the accuracies were uncorrelated to the degree of relatedness to the training population. The results underline the potential of genomic prediction to accelerate breeding progress in outbred fruit tree crops that still need to overcome long generation intervals and extensive phenotyping costs.

A qRT-PCR assay for the expression of all Mal d 1 isoallergen genes.

BACKGROUND: A considerable number of individuals suffer from oral allergy syndrome (OAS) to apple, resulting in the avoidance of apple consumption. Apple cultivars differ greatly in their allergenic properties, but knowledge of the causes for such differences is incomplete. Mal d 1 is considered the major apple allergen. For Mal d 1, a wide range of isoallergens and variants exist, and they are encoded by a large gene family. To identify the specific proteins/genes that are potentially involved in the allergy, we developed a PCR assay to monitor the expression of each individual Mal d 1 gene. Gene-specific primer pairs were designed for the exploitation of sequence differences among Mal d 1 genes. The specificity of these primers was validated using both in silico and in vitro techniques. Subsequently, this assay was applied to the peel and flesh of fruits from the two cultivars 'Florina' and 'Gala'. RESULTS: We successfully developed gene-specific primer pairs for each of the 31 Mal d 1 genes and incorporated them into a qRT-PCR assay. The results from the application of the assay showed that 11 genes were not expressed in fruit. In addition, differential expression was observed among the Mal d 1 genes that were expressed in the fruit. Moreover, the expression levels were tissue and cultivar dependent. CONCLUSION: The assay developed in this study facilitated the first characterisation of the expression levels of all known Mal d 1 genes in a gene-specific manner. Using this assay on different fruit tissues and cultivars, we obtained knowledge concerning gene relevance in allergenicity. This study provides new perspectives for research on both plant breeding and immunotherapy.

Pedimap: software for the visualization of genetic and phenotypic data in pedigrees.

Pedimap is a user-friendly software tool for visualizing phenotypic and genotypic data for related individuals linked in pedigrees. Genetic data can include marker scores, Identity-by-Descent probabilities, and marker linkage map positions, allowing the visualization of haplotypes through lineages. The pedigrees can accommodate all types of inheritance, including selfing, cloning, and repeated backcrossing, and all ploidy levels are supported. Visual association of the genetic data with phenotypic data simplifies the exploration of large data sets, thereby improving breeding decision making. Data are imported from text files; in addition data exchange with other software packages (FlexQTL(TM) and GenomeStudio(TM)) is possible. Instructions for use and an executable version compatible with the Windows platform are available for free from http://www.plantbreeding.wur.nl/UK/software_pedimap.html.

Diversity arrays technology (DArT) markers in apple for genetic linkage maps.

Diversity Arrays Technology (DArT) provides a high-throughput whole-genome genotyping platform for the detection and scoring of hundreds of polymorphic loci without any need for prior sequence information. The work presented here details the development and performance of a DArT genotyping array for apple. This is the first paper on DArT in horticultural trees. Genetic mapping of DArT markers in two mapping populations and their integration with other marker types showed that DArT is a powerful high-throughput method for obtaining accurate and reproducible marker data, despite the low cost per data point. This method appears to be suitable for aligning the genetic maps of different segregating populations. The standard complexity reduction method, based on the methylation-sensitive PstI restriction enzyme, resulted in a high frequency of markers, although there was 52-54% redundancy due to the repeated sampling of highly similar sequences. Sequencing of the marker clones showed that they are significantly enriched for low-copy, genic regions. The genome coverage using the standard method was 55-76%. For improved genome coverage, an alternative complexity reduction method was examined, which resulted in less redundancy and additional segregating markers. The DArT markers proved to be of high quality and were very suitable for genetic mapping at low cost for the apple, providing moderate genome coverage. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-011-9579-5) contains supplementary material, which is available to authorized users.

QTL dynamics for fruit firmness and softening around an ethylene-dependent polygalacturonase gene in apple (Malus x domestica Borkh.).

Apple fruit are well known for their storage life, although a wide range of flesh softening occurs among cultivars. Loss of firmness is genetically coordinated by the action of several cell wall enzymes, including polygalacturonase (PG) which depolymerizes cell wall pectin. By the analysis of 'Fuji' (Fj) and 'Mondial Gala' (MG), two apple cultivars characterized by a distinctive ripening behaviour, the involvement of Md-PG1 in the fruit softening process was confirmed to be ethylene dependent by its transcript being down-regulated by 1-methylcyclopropene treatment in MG and in the low ethylene-producing cultivar Fj. Comparing the PG sequence of MG and Fj, a single nucleotide polymorphism (SNP) was discovered. Segregation of the Md-PG1(SNP) marker within a full-sib population, obtained by crossing Fj and MG, positioned Md-PG1 in the linkage group 10 of MG, co-located with a quantitative trait locus (QTL) identified for fruit firmness in post-harvest ripening. Fruit firmness and softening analysed in different stages, from harvest to post-storage, determined a shift of the QTL from the top of this linkage group to the bottom, where Md-ACO1, a gene involved in ethylene biosynthesis in apple, is mapped. This PG-ethylene-related gene has beeen positioned in the apple genome on chromosome 10, which contains several QTLs controlling fruit firmness and softening, and the interplay among the allelotypes of the linked loci should be considered in the design of a marker-assisted selection breeding scheme for apple texture.

Comparison of marker-based pairwise relatedness estimators on a pedigreed plant population.

Several estimators have been proposed that use molecular marker data to infer the degree of relatedness for pairs of individuals. The objective of this study was to evaluate the performance of seven estimators when applied to marker data of a set of 33 key individuals from a large complex apple pedigree. The evaluation considered different scenarios of allele frequencies and different numbers of marker loci. The method of moments estimators were Similarity, Queller-Goodknight, Lynch-Ritland and Wang. The maximum likelihood estimators were Thompson, Anderson-Weir and Jacquard. The pedigree-based coancestry coefficients were taken as the point of reference in calculating correlations and root mean square error (RMSE). The marker data comprised 86 multi-allelic SSR markers on 17 linkage groups, covering 11 Morgans. Additionally, we simulated 10 datasets conditional on the real pedigree to support the results on the real dataset. None of the estimators outperformed the others. Knowledge of allele frequencies appeared to be the most influential, i.e., the highest correlations and lowest RMSE were found when frequencies from the founder population were available. When equal allele frequencies were used, all estimators resulted in very similar, but on average lower, correlations. The use of allele frequencies estimated from the set of 33 individuals gave, on average, the poorest results. The maximum likelihood estimators and the Lynch-Ritland estimator were the most sensitive to allele frequencies. The results from the simulation study fully supported the trends in results of the real dataset. This study indicated that high correlations (up to 0.90) and small RMSE (below 0.03), may be obtained when population allelic frequencies are available. In this scenario, the performances of the various estimators were similar, but seemed to favor the maximum likelihood estimators. In the absence of reliable allele frequencies the method of moments estimators were shown to be more robust. The number of marker loci influenced the average performance of the estimators; however, the ranking was not affected. Correlations up to 0.80 were obtained when two markers per chromosome and appropriate allele frequencies were available. Adding more markers to the current dataset may lead to marginal improvements.

Linked vs unlinked markers: multilocus microsatellite haplotype-sharing as a tool to estimate gene flow and introgression.

We have explored the use of multilocus microsatellite haplotypes to study introgression from cultivated (Malus domestica) into wild apple (Malus sylvestris), and to study gene flow among remnant populations of M. sylvestris. A haplotype consisted of alleles at microsatellite loci along one chromosome. As destruction of haplotypes through recombination occurs much faster than loss of alleles due to genetic drift, the lifespan of a multilocus haplotype is much shorter than that of the underlying alleles. When different populations share the same haplotype, this may indicate recent gene flow between populations. Similarly, haplotypes shared between two species would be a strong signal for introgression. As the expected lifespan of a haplotype depends on the strength of the linkage, the length [in centiMorgans (cM)] of the haplotype shared contains information on the number of generations passed. This application of shared haplotypes is distinct from using haplotype-sharing to detect association between markers and a certain trait. We inferred haplotypes for four to eight microsatellite loci on Linkage Group 10 of apple from genotype data using the program phase, and then identified those haplotypes shared between populations and species. Compared with a Bayesian analysis of unlinked microsatellite loci using the program structure, haplotype-sharing detected a partially different set of putative hybrids. Cultivated haplotypes present in M. sylvestris were short (< 1.5 cM), indicating that introgression had taken place many generations ago, except for two Belgian plants that contained a haplotype of 47.1 cM, indicating recent introgression. In the estimation of gene flow, F(ST) based on unlinked loci indicated small (0.032-0.058) but statistically significant differentiation between some populations only. However, various M. sylvestris haplotypes were shared in nearly all pairwise comparisons of populations, and their length indicated recent gene flow. Hence, all Dutch populations should be considered as one conservation unit. The added value of using sharing of multilocus microsatellite haplotypes as a source of population genetic information is discussed.

In vivo assessment with prick-to-prick testing and double-blind, placebo-controlled food challenge of allergenicity of apple cultivars.

BACKGROUND: Apple cultivars have been reported to differ in allergenicity on the basis of in vitro and skin prick tests with apple extracts. OBJECTIVES: We sought to evaluate the efficacy of the prick-to-prick method in assessing differences in allergenicity of apple cultivars and to confirm differences by means of double-blind, placebo-controlled food challenge (DBPCFC). METHODS: Intra-assay and intracultivar variation of prick-to-prick test results were determined in 6 Dutch and 8 Spanish patients with apple allergy by using 5 apples of the cultivars Golden Delicious, Fuji, and Ecolette in duplicate. In addition, 21 cultivars were screened for allergenicity in 15 Dutch patients with birch pollen and apple allergy. Two selected cultivars (Golden Delicious and Santana) were tested with DBPCFCs. The influence of storage conditions on allergenicity was assessed in 5 cultivars. RESULTS: Intra-assay variation of skin prick testing was 3.9%, and intracultivar variation was 4.1%. A ranking of 21 cultivars was made on the basis of prick-to-prick tests in 9 patients. Apple cultivars were classified as of low, intermediate, and high allergenicity, with a significant difference between low and high allergenicity (P < .001). A significant difference in allergenicity determined between Golden Delicious and Santana cultivars (P < .05) was confirmed by means of DBPCFC. With 5 cultivars, controlled atmosphere (2.5% oxygen/1% carbon dioxide) was shown to reduce allergenicity (P < .001) by 15% compared with storage at 2 degrees C. CONCLUSIONS: Prick-to-prick testing with fresh apples is a reproducible method of assessing allergenicity. Apples can be classified as of low or high allergenicity for the majority of patients. This was confirmed by using DBPCFCs. Selection of cultivars and control of storage conditions are both viable strategies for reduction of symptoms in patients with apple allergy.

The Vh8 locus of a new gene-for-gene interaction between Venturia inaequalis and the wild apple Malus sieversii is closely linked to the Vh2 locus in Malus pumila R12740-7A.

The wild apple (Malus sieversii) is a large-fruited species from Central Asia, which is used as a source of scab resistance in cultivar breeding. Phytopathological tests with races of Venturia inaequalis were performed to differentiate scab-resistance genes in Malus as well as an avirulence gene in the pathogen. A novel gene-for-gene interaction between V. inaequalis and Malus was identified. The locus of the scab-resistance gene Vh8 is linked with, or possibly allelic to, that of the Vh2 gene in Malus pumila Russian apple R12740-7A, at the lower end of linkage group 2 of Malus. Race 8 isolate NZ188B.2 is compatible with Vh8, suggesting the loss or modification of the complementary AvrVh8 gene, while isolate 1639 overcomes both Vh2 and Vh8, but is incompatible with at least one other gene not detected by any of the other race isolates tested. Our research is the first to differentiate scab-resistance genes in a putative gene cluster in apple with the aid of races of V. inaequalis.

Silencing the major apple allergen Mal d 1 by using the RNA interference approach.

BACKGROUND: Apple allergy is dominated by IgE antibodies against Mal d 1 in areas where birch pollen is endemic. Apples with significantly decreased levels of Mal d 1 would allow most patients in these areas to eat apples without allergic reactions. OBJECTIVE: The aim of this study was to inhibit the expression of Mal d 1 in apple plants by RNA interference. METHODS: In vitro -grown apple plantlets were transformed with a construct coding for an intron-spliced hairpin RNA containing a Mal d 1-specific inverted repeat sequence separated by a Mal d 1-specific intron sequence. The presence of the construct in transformants was checked by PCR. Expression of Mal d 1 in leaves was monitored by prick-to-prick skin testing in 3 patients allergic to apples and by immunoblotting with a Mal d 1-reactive mAb and with IgE antibodies against Mal d 1. RESULTS: After transformation, plantlets were selected on the basis of having a normal phenotype and growth rate. With PCR, in 6 of 9 selected plantlets, the presence of the gene-silencing construct was demonstrated. By skin prick test it was shown that a wild-type plantlet had significantly ( P < .05) higher allergenicity than 5 of the transformants. Reduction of expression of Mal d 1 was confirmed by immunoblotting. In wild-type and unsuccessful transformants, a strong band was detected with Mal d 1-reactive mAb 5H8 at the expected apparent M r of 17 kDa. This band was virtually absent in the transformants that carried the gene-silencing construct. With human IgE antibodies, the same observations were made. CONCLUSIONS: Mal d 1 expression was successfully reduced by RNA interference. This translated into significantly reduced in vivo allergenicity. These observations support the feasibility of the production by gene silencing of apples hypoallergenic for Mal d 1.