Proof of concept of Bayesian integrated QTL analyses across pedigree-related families from breeding programs of an outbreeding species. Results include QTL confidence intervals, individuals' genotype probabilities and genomic breeding values. Bayesian QTL linkage mapping approaches offer the flexibility to study multiple full sib families with known pedigrees simultaneously. Such a joint analysis increases the probability of detecting these quantitative trait loci (QTL) and provide insight of the magnitude of QTL across different genetic backgrounds. Here, we present an improved Bayesian multi-QTL pedigree-based approach on an outcrossing species using progenies with different (complex) genetic relationships. Different modeling assumptions were studied in the QTL analyses, i.e., the a priori expected number of QTL varied and polygenic effects were considered. The inferences include number of QTL, additive QTL effect sizes and supporting credible intervals, posterior probabilities of QTL genotypes for all individuals in the dataset, and QTL-based as well as genome-wide breeding values. All these features have been implemented in the FlexQTL(™) software. We analyzed fruit firmness in a large apple dataset that comprised 1,347 individuals forming 27 full sib families and their known ancestral pedigrees, with genotypes for 87 SSR markers on 17 chromosomes. We report strong or positive evidence for 14 QTL for fruit firmness on eight chromosomes, validating our approach as several of these QTL were reported previously, though dispersed over a series of studies based on single mapping populations. Interpretation of linked QTL was possible via individuals' QTL genotypes. The correlation between the genomic breeding values and phenotypes was on average 90 %, but varied with the number of detected QTL in a family. The detailed posterior knowledge on QTL of potential parents is critical for the efficiency of marker-assisted breeding.

BackgroundThe amount and structure of genetic diversity in dessert apple germplasm conserved at a European level is mostly unknown, since all diversity studies conducted in Europe until now have been performed on regional or national collections. Here, we applied a common set of 16 SSR markers to genotype more than 2,400 accessions across 14 collections representing three broad European geographic regions (North + East, West and South) with the aim to analyze the extent, distribution and structure of variation in the apple genetic resources in Europe.ResultsA Bayesian model-based clustering approach showed that diversity was organized in three groups, although these were only moderately differentiated (FST = 0.031). A nested Bayesian clustering approach allowed identification of subgroups which revealed internal patterns of substructure within the groups, allowing a finer delineation of the variation into eight subgroups (FST = 0.044). The first level of stratification revealed an asymmetric division of the germplasm among the three groups, and a clear association was found with the geographical regions of origin of the cultivars. The substructure revealed clear partitioning of genetic groups among countries, but also interesting associations between subgroups and breeding purposes of recent cultivars or particular usage such as cider production. Additional parentage analyses allowed us to identify both putative parents of more than 40 old and/or local cultivars giving interesting insights in the pedigree of some emblematic cultivars.ConclusionsThe variation found at group and subgroup levels may reflect a combination of historical processes of migration/selection and adaptive factors to diverse agricultural environments that, together with genetic drift, have resulted in extensive genetic variation but limited population structure. The European dessert apple germplasm represents an important source of genetic diversity with a strong historical and patrimonial value. The present work thus constitutes a decisive step in the field of conservation genetics. Moreover, the obtained data can be used for defining a European apple core collection useful for further identification of genomic regions associated with commercially important horticultural traits in apple through genome-wide association studies.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-016-0818-0) contains supplementary material, which is available to authorized users.

Efficient breeding and selection of high-quality apple cultivars requires knowledge and understanding of the underlying genetics. The availability of genetic linkage maps constructed with molecular markers enables the detection and analysis of major genes and quantitative trait loci contributing to the quality traits of a genotype. A segregating population of the cross between the apple varieties 'Fiesta' (syn. 'Red Pippin') and 'Discovery' has been observed over three years at three different sites in Switzerland and data on growth habit, blooming behaviour, juvenile period and fruit quality has been recorded. QTL analyses were performed, based on a genetic linkage map consisting of 804 molecular markers and covering all 17 apple chromosomes. With the maximum likelihood based interval mapping method, the investigated complex traits could be dissected into a number of QTLs affecting the observed characters. Genomic regions participating in the genetic control of stem diameter, plant height increment, leaf size, blooming time, blooming intensity, juvenile phase length, time of fruit maturity, number of fruit, fruit size and weight, fruit flesh firmness, sugar content and fruit acidity were identified and compared with previously mapped QTLs in apple. Although 'Discovery' fruit displayed a higher acid content, both acidity QTLs were attributed to the sweeter parent 'Fiesta'. This indicated homozygosity at the acidity loci in 'Discovery' preventing their detection in the progeny due to the lack of segregation.

Breeding for scab-resistant apple cultivars by pyramiding several resistance genes in the same genetic background is a promising way to control apple scab caused by the fungus Venturia inaequalis. To achieve this goal, DNA markers linked to the genes of interest are required in order to select seedlings with the desired resistance allele combinations. For several apple scab resistance genes, molecular markers are already available; but until now, none existed for the apple scab resistance gene Vbj originating from the crab apple Malus baccata jackii. Using bulk segregant analysis, three RAPD markers linked to Vbj were first identified. These markers were transformed into more reliable sequence-characterised amplified region (SCAR) markers that proved to be codominant. In addition, three SSR markers and one SCAR were identified by comparing homologous linkage groups of existing genetic maps. Discarding plants showing genotype-phenotype incongruence (GPI plants) plants, a linkage map was calculated. Vbj mapped between the markers CH05e03 (SSR) and T6-SCAR, at 0.6 cM from CH05e03 and at 3.9 cM from T6-SCAR. Without the removal of the GPI plants, Vbj was placed 15 cM away from the closest markers. Problems and pitfalls due to GPI plants and the consequences for mapping the resistance gene accurately are discussed. Finally, the usefulness of co-dominant markers for pedigree analysis is also demonstrated.

Fire blight is the most important bacterial disease in apple (Malus×domestica) and pear (Pyrus communis) production. Today, the causal bacterium Erwinia amylovora is present in many apple-and pear-growing areas. We investigated the natural resistance of the wild apple Malus×robusta 5 against E. amylovora, previously mapped to linkage group 3. With a fine-mapping approach on a population of 2,133 individuals followed by phenotyping of the recombinants from the region of interest, we developed flanking markers useful for marker-assisted selection. Open reading frames were predicted on the sequence of a BAC spanning the resistance locus. One open reading frame coded for a protein belonging to the NBS-LRR family. The in silico investigation of the structure of the candidate resistance gene against fire blight of M.×robusta 5, FB_MR5, led us hypothesize the presence of a coiled-coil region followed by an NBS and an LRR-like structure with the consensus 'LxxLx[IL]xxCxxLxxL'. The function of FB_MR5 was predicted in agreement with the decoy/guard model, that FB_MR5 monitors the transcribed RIN4_MR5, a homolog of RIN4 of Arabidopsis thaliana that could interact with the previously described effector AvrRpt2 EA of E. amylovora.