Identification of Pyrus single nucleotide polymorphisms (SNPs) and evaluation for genetic mapping in European pear and interspecific Pyrus hybrids

A multiplexed plant-animal SNP array for selective breeding and species conservation applications

Reliable and high-throughput genotyping platforms are of immense importance for identifying and dissecting genomic regions controlling important phenotypes, supporting selection processes in breeding programs, and managing wild populations and germplasm collections. Amongst available genotyping tools, single nucleotide polymorphism arrays have been shown to be comparatively easy to use and generate highly accurate genotypic data. Single-species arrays are the most commonly used type so far; however, some multi-species arrays have been developed for closely related species that share single nucleotide polymorphism markers, exploiting inter-species cross-amplification. In this study, the suitability of a multiplexed plant-animal single nucleotide polymorphism array, including both closely and distantly related species, was explored. The performance of the single nucleotide polymorphism array across species for diverse applications, ranging from intra-species diversity assessments to parentage analysis, was assessed. Moreover, the value of genotyping pooled DNA of distantly related species on the single nucleotide polymorphism array as a technique to further reduce costs was evaluated. Single nucleotide polymorphism performance was generally high, and species-specific single nucleotide polymorphisms proved suitable for diverse applications. The multi-species single nucleotide polymorphism array approach reported here could be transferred to other species to achieve cost savings resulting from the increased throughput when several projects use the same array, and the pooling technique adds another highly promising advancement to additionally decrease genotyping costs by half.

SNP-based high-density linkage map construction and QTL mapping of black spot disease resistance in Chinese sand pear

Black spot disease (PBS) caused by Alternaria alternata is an economic disease of pear (Pyrus pyrifolia Nakai). Developing cultivars with durable PBS resistance traits is an important research objective for improving pear germplasm. The Deshengxiang is a popular pear variety in China and resistant to PBS. This study aimed to detect quantitative trait loci (QTL) associated with PBS resistance trait in pear and determine closely linked molecular markers by specific locus amplified fragment sequencing (SLAF-seq). F1 population resulting from a cross between "Deshengxiang" (female) and "Guiguan," a susceptible (male) variety, was developed and evaluated in 2016 and 2017. SLAF technology was used to discover SNPs in the F1 individuals and subsequently a high-density genetic linkage map for PBS resistance was constructed which contained 17,604 SNP markers. Based on the linkage map, the markers were distributed into 17 linkage groups, spanning 1548.48 cM, with a mean marker distance of 0.09 cM, representing the densest genetic map of the genus Pyrus. QTL analysis of PBS resistance identified a locus strongly related to PBS resistance at 77.68 ~ 112.99 cM on linkage group 15, which was further narrowed down to 93.79 ~ 112.99 cM. Two markers, Marker94293 and Marker94206, located at 97.47 and 102.93 cM, were closely associated with PBS resistance, with a Δ (SNP index) value of 0.46. Co-localization of QTL interval, bioinformatics analysis, and functional annotation revealed PBS putative candidate genes. Overall, the high-density pear linkage map is a suitable reference for mapping PBS resistance trait, QTL, and genes identified in this study contribute information that could be useful for PBS improvement in pear.

A high-density integrated map for grapevine based on three mapping populations genotyped by the Vitis18K SNP chip

We present a high-density integrated map for grapevine, allowing refinement and improved understanding of the grapevine genome, while demonstrating the applicability of the Vitis18K SNP chip for linkage mapping. The improvement of grapevine through biotechnology requires identification of the molecular bases of target traits by studying marker-trait associations. The Vitis18K SNP chip provides a useful genotyping tool for genome-wide marker analysis. Most linkage maps are based on single mapping populations, but an integrated map can increase marker density and show order conservation. Here we present an integrated map based on three mapping populations. The parents consist of the well-known wine cultivars 'Cabernet Sauvignon', 'Corvina' and 'Rhine Riesling', the lesser-known wine variety 'Deckrot', and a table grape selection, G1-7720. Three high-density population maps with an average inter-locus gap ranging from 0.74 to 0.99 cM were developed. These maps show high correlations (0.9965-0.9971) with the reference assembly, containing only 93 markers with large order discrepancies compared to expected physical positions, of which a third is consistent across multiple populations. Moreover, the genetic data aid the further refinement of the grapevine genome assembly, by anchoring 104 yet unanchored scaffolds. From these population maps, an integrated map was constructed which includes 6697 molecular markers and reduces the inter-locus gap distance to 0.60 cM, resulting in the densest integrated map for grapevine thus far. A small number of discrepancies, mainly of short distance, involve 88 markers that remain conflictual across maps. The integrated map shows similar collinearity to the reference assembly (0.9974) as the single maps. This high-density map increases our understanding of the grapevine genome and provides a useful tool for its further characterization and the dissection of complex traits.

Molecular characterization of intergeneric hybrids between Malus and Pyrus

Apple (Malus) and pear (Pyrus) are economically important fruit crops well known for their unique textures, flavours, and nutritional qualities. Both genera are characterised by a distinct pattern of secondary metabolites, which directly affect not only resistance to certain diseases, but also have significant impacts on the flavour and nutritional value of the fruit. The identical chromosome numbers, similar genome size, and their recent divergence date, together with DNA markers have shown that apple and pear genomes are highly co-linear. This study utilized comparative genomic approaches, including simple sequence repeats, high resolution single nucleotide polymorphism melting analysis, and single nucleotide polymorphism chip analysis to identify genetic differences among hybrids of Malus and Pyrus, and F2 offspring. This research has demonstrated and validated that these three marker types, along with metabolomics analysis are very powerful tools to detect and confirm hybridity of progeny derived from crosses between apple and pear in both cross directions. Furthermore, this work analysed the genus-specific metabolite patterns and the resistance to fire blight (Erwinia amylovora) in progeny. The findings of this work will enhance and accelerate the breeding of novel tree fruit crops that benefit producers and consumers, by enabling marker assisted selection of desired traits introgressed between pear and apple.

Applications of Genomic Tools in Plant Breeding: Crop Biofortification

Crop breeding has mainly been focused on increasing productivity, either directly or by decreasing the losses caused by biotic and abiotic stresses (that is, incorporating resistance to diseases and enhancing tolerance to adverse conditions, respectively). Quite the opposite, little attention has been paid to improve the nutritional value of crops. It has not been until recently that crop biofortification has become an objective within breeding programs, through either conventional methods or genetic engineering. There are many steps along this long path, from the initial evaluation of germplasm for the content of nutrients and health-promoting compounds to the development of biofortified varieties, with the available and future genomic tools assisting scientists and breeders in reaching their objectives as well as speeding up the process. This review offers a compendium of the genomic technologies used to explore and create biodiversity, to associate the traits of interest to the genome, and to transfer the genomic regions responsible for the desirable characteristics into potential new varieties. Finally, a glimpse of future perspectives and challenges in this emerging area is offered by taking the present scenario and the slow progress of the regulatory framework as the starting point.

Construction of transgenic detection system of Brassica napus L. based on single nucleotide polymorphism chip

Brassica napus L. is a vital oil crop in China. As auxiliary tools for rapeseed breeding, transgenic technologies play a considerable role in heterosis, variety improvement, and pest resistance. Research on transgenic detection technologies is of great significance for the introduction, supervision, and development of transgenic rapeseed in China. However, the transgenic detection methods currently in use are complex and time-consuming, with low output. A single nucleotide polymorphism (SNP) chip can effectively overcome such limitations. In the present study, we collected 40 transgenic elements and designed 291 probes. The probe sequences were submitted to Illumina Company, and the Infinium chip technology was used to prepare SNP chips. In the present Brassica napus transgenic detection experiment, 84 high-quality probes of 17 transgenic elements were preliminarily screened, and genotyping effect was optimised for the probe signal value. Ultimately, a transgenic detection system for B. napus was developed. The developed system has the advantages of simple operation, minimal technical errors, and stable detection outcomes. A transgenic detection sensitivity test revealed that the probe designed could accurately detect 1% of transgenic samples and had high detection sensitivity. In addition, in repeatability tests, the CaMV35S promoter coefficient of variation was approximately 3.58%. Therefore, the SNP chip had suitable repeatability in transgene detection. The SNP chip developed could be used to construct transgenic detection systems for B. napus.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03062-6.

Dissecting Genetic Resistance to Fire Blight in Three Pear Populations

Fire blight, caused by the bacterial pathogen Erwinia amylovora, is a persistent problem for pear (Pyrus spp.) growers in the United States. Growing resistant cultivars is one of the best options for managing fire blight. The cultivars Potomac and Old Home and the selection NJA2R59T69 display resistance to fire blight. As such, three mapping populations (El Dorado × Potomac, Old Home × Bartlett, and NJA2R59T69 × Bartlett) were developed to identify genomic regions associated with resistance to fire blight. Progeny were phenotyped during 2017 and 2018 by inoculating multiple actively growing shoots of field-grown seedling trees with E. amylovora isolate E153n via the cut-leaf method. Genotyping was conducted using the recently developed Axiom Pear 70 K Genotyping Array and chromosomal linkage groups were created for each population. An integrated two-way pseudo-testcross approach was used to map quantitative trait loci (QTLs). Resistance QTLs were identified on chromosome 2 for each population. The QTLs identified in the El Dorado × Potomac and Old Home × Bartlett populations are in the same region as QTLs that were previously identified in Harrow Sweet and Moonglow. The QTL in NJA2R59T69 mapped proximally to the previously identified QTLs and originated from an unknown Asian or occidental source. Future research will focus on further characterizing the resistance regions and developing tools for DNA-informed breeding.

A de novo genome assembly of the dwarfing pear rootstock Zhongai 1

'Zhongai 1' [(Pyrus ussuriensis × communis) × spp.] is an excellent pear dwarfing rootstock common in China. It is dwarf itself and has high dwarfing efficiency on most of main Pyrus cultivated species when used as inter-stock. Here we describe the draft genome sequences of 'Zhongai 1' which was assembled using PacBio long reads, Illumina short reads and Hi-C technology. We estimated the genome size is approximately 511.33 Mb by K-mer analysis and obtained a final genome of 510.59 Mb with a contig N50 size of 1.28 Mb. Next, 506.31 Mb (99.16%) of contigs were clustered into 17 chromosomes with a scaffold N50 size of 23.45 Mb. We further predicted 309.86 Mb (60.68%) of repetitive sequences and 43,120 protein-coding genes. The assembled genome will be a valuable resource and reference for future pear breeding, genetic improvement, and comparative genomics among related species. Moreover, it will help identify genes involved in dwarfism, early flowering, stress tolerance, and commercially desirable fruit characteristics.

Development of an integrated 200K SNP genotyping array and application for genetic mapping, genome assembly improvement and genome wide association studies in pear (Pyrus)

Pear (Pyrus; 2n = 34), the third most important temperate fruit crop, has great nutritional and economic value. Despite the availability of many genomic resources in pear, it is challenging to genotype novel germplasm resources and breeding progeny in a timely and cost-effective manner. Genotyping arrays can provide fast, efficient and high-throughput genetic characterization of diverse germplasm, genetic mapping and breeding populations. We present here 200K AXIOM^® PyrSNP, a large-scale single nucleotide polymorphism (SNP) genotyping array to facilitate genotyping of Pyrus species. A diverse panel of 113 re-sequenced pear genotypes was used to discover SNPs to promote increased adoption of the array. A set of 188 diverse accessions and an F₁ population of 98 individuals from 'Cuiguan' × 'Starkrimson' was genotyped with the array to assess its effectiveness. A large majority of SNPs (166 335 or 83%) are of high quality. The high density and uniform distribution of the array SNPs facilitated prediction of centromeric regions on 17 pear chromosomes, and significantly improved the genome assembly from 75.5% to 81.4% based on genetic mapping. Identification of a gene associated with flowering time and candidate genes linked to size of fruit core via genome wide association studies showed the usefulness of the array in pear genetic research. The newly developed high-density SNP array presents an important tool for rapid and high-throughput genotyping in pear for genetic map construction, QTL identification and genomic selection.

Development of a highly efficient Axiom™ 70 K SNP array for Pyrus and evaluation for high-density mapping and germplasm characterization

Background

Both a source of diversity and the development of genomic tools, such as reference genomes and molecular markers, are equally important to enable faster progress in plant breeding. Pear (Pyrus spp.) lags far behind other fruit and nut crops in terms of employment of available genetic resources for new cultivar development. To address this gap, we designed a high-density, high-efficiency and robust single nucleotide polymorphism (SNP) array for pear, with the main objectives of conducting genetic diversity and genome-wide association studies.

Results

By applying a two-step design process, which consisted of the construction of a first ‘draft’ array for the screening of a small subset of samples, we were able to identify the most robust and informative SNPs to include in the Applied Biosystems™ Axiom™ Pear 70 K Genotyping Array, currently the densest SNP array for pear. Preliminary evaluation of this 70 K array in 1416 diverse pear accessions from the USDA National Clonal Germplasm Repository (NCGR) in Corvallis, OR identified 66,616 SNPs (93% of all the tiled SNPs) as high quality and polymorphic (PolyHighResolution). We further used the Axiom Pear 70 K Genotyping Array to construct high-density linkage maps in a bi-parental population, and to make a direct comparison with available genotyping-by-sequencing (GBS) data, which suggested that the SNP array is a more robust method of screening for SNPs than restriction enzyme reduced representation sequence-based genotyping.

Conclusions

The Axiom Pear 70 K Genotyping Array, with its high efficiency in a widely diverse panel of Pyrus species and cultivars, represents a valuable resource for a multitude of molecular studies in pear. The characterization of the USDA-NCGR collection with this array will provide important information for pear geneticists and breeders, as well as for the optimization of conservation strategies for Pyrus.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5712-3) contains supplementary material, which is available to authorized users.

Genetic Mapping of the Incompatibility Locus in Olive and Development of a Linked Sequence-Tagged Site Marker

The genetic control of self-incompatibility (SI) has been recently disclosed in olive. Inter-varietal crossing confirmed the presence of only two incompatibility groups (G1 and G2), suggesting a simple Mendelian inheritance of the trait. A double digest restriction associated DNA (ddRAD) sequencing of a biparental population segregating for incompatibility groups has been performed and high-density linkage maps were constructed in order to map the SI locus and identify gene candidates and linked markers. The progeny consisted of a full-sib family of 229 individuals derived from the cross 'Leccino' (G1) × 'Dolce Agogia' (G2) varieties, segregating 1:1 (G1:G2), in accordance with a diallelic self-incompatibility (DSI) model. A total of 16,743 single nucleotide polymorphisms was identified, 7,006 in the female parent 'Leccino' and 9,737 in the male parent 'Dolce Agogia.' Each parental map consisted of 23 linkage groups and showed an unusual large size (5,680 cM in 'Leccino' and 3,538 cM in 'Dolce Agogia'). Recombination was decreased across all linkage groups in pollen mother cells of 'Dolce Agogia,' the parent with higher heterozygosity, compared to megaspore mother cells of 'Leccino,' in a context of a species that showed exceptionally high recombination rates. A subset of 109 adult plants was assigned to either incompatibility group by a stigma test and the diallelic self-incompatibility (DSI) locus was mapped to an interval of 5.4 cM on linkage group 18. This region spanned a size of approximately 300 Kb in the olive genome assembly. We developed a sequence-tagged site marker in the DSI locus and identified five haplotypes in 57 cultivars with known incompatibility group assignment. A combination of two single-nucleotide polymorphisms (SNPs) was sufficient to predict G1 or G2 phenotypes in olive cultivars, enabling early marker-assisted selection of compatible genotypes and allowing for a rapid screening of inter-compatibility among cultivars in order to guarantee effective fertilization and increase olive production. The construction of high-density linkage maps has led to the development of the first functional marker in olive and provided positional candidate genes in the SI locus.

Chromosome level high-density integrated genetic maps improve the Pyrus bretschneideri 'DangshanSuli' v1.0 genome

Background

Chromosomal level reference genomes provide a crucial foundation for genomics research such as genome-wide association studies (GWAS) and whole genome selection. The chromosomal-level sequences of both the European (Pyrus communis) and Chinese (P. bretschneideri) pear genomes have not been published in public databases so far.

Results

To anchor the scaffolds of P. bretschneideri ‘DangshanSuli’ (DS) v1.0 genome into pseudo-chromosomes, two genetic maps (MH and YM maps) were constructed using half sibling populations of Chinese pear crosses, ‘Mantianhong’ (MTH) × ‘Hongxiangsu’ (HXS) and ‘Yuluxiang’ (YLX) × MTH, from 345 and 162 seedlings, respectively, which were prepared for SNP discovery using genotyping-by-sequencing (GBS) technology. The MH and YM maps, each with 17 linkage groups (LGs), were constructed from 2606 and 2489 SNP markers and spanned 1847 and 1668 cM, respectively, with average marker intervals of 0.7. The two maps were further merged with a previously published genetic map (BD) based on the cross ‘Bayuehong’ (BYH) × ‘Dangshansuli’ (DS) to build a new integrated MH-YM-BD map. By using 7757 markers located on the integrated MH-YM-BD map, 898 scaffolds (400.57 Mb) of the DS v1.0 assembly were successfully anchored into 17 pseudo-chromosomes, accounting for 78.8% of the assembled genome size. About 88.31% of them (793 scaffolds) were directionally anchored with two or more markers on the pseudo-chromosomes. Furthermore, the errors in each pseudo-chromosome (especially 1, 5, 7 and 11) were manually corrected and pseudo-chromosomes 1, 5 and 7 were extended by adding 19, 12 and 14 scaffolds respectively in the newly constructed DS v1.1 genome. Synteny analyses revealed that the DS v1.1 genome had high collinearity with the apple genome, and the homologous fragments between pseudo-chromosomes were similar to those found in previous studies. Moreover, the red-skin trait of Asian pear was mapped to an identical locus as identified previously.

Conclusions

The accuracy of DS v1.1 genome was improved by using larger mapping populations and merged genetic map. With more than 400 MB anchored to 17 pseudo-chromosomes, the new DS v1.1 genome provides a critical tool that is essential for studies of pear genetics, genomics and molecular breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5224-6) contains supplementary material, which is available to authorized users.

Genome-wide genetic variation and comparison of fruit-associated traits between kumquat (Citrus japonica) and Clementine mandarin (Citrus clementina)

The comprehensive genetic variation of two citrus species were analyzed at genome and transcriptome level. A total of 1090 differentially expressed genes were found during fruit development by RNA-sequencing. Fruit size (fruit equatorial diameter) and weight (fresh weight) are the two most important components determining yield and consumer acceptability for many horticultural crops. However, little is known about the genetic control of these traits. Here, we performed whole-genome resequencing to reveal the comprehensive genetic variation of the fruit development between kumquat (Citrus japonica) and Clementine mandarin (Citrus clementina). In total, 5,865,235 single-nucleotide polymorphisms (SNPs) and 414,447 insertions/deletions (InDels) were identified in the two citrus species. Based on integrative analysis of genome and transcriptome of fruit, 640,801 SNPs and 20,733 InDels were identified. The features, genomic distribution, functional effect, and other characteristics of these genetic variations were explored. RNA-sequencing identified 1090 differentially expressed genes (DEGs) during fruit development of kumquat and Clementine mandarin. Gene Ontology revealed that these genes were involved in various molecular functional and biological processes. In addition, the genetic variation of 939 DEGs and 74 multiple fruit development pathway genes from previous reports were also identified. A global survey identified 24,237 specific alternative splicing events in the two citrus species and showed that intron retention is the most prevalent pattern of alternative splicing. These genome variation data provide a foundation for further exploration of citrus diversity and gene-phenotype relationships and for future research on molecular breeding to improve kumquat, Clementine mandarin and related species.

Crop Breeding Chips and Genotyping Platforms: Progress, Challenges, and Perspectives

There is a rapidly rising trend in the development and application of molecular marker assays for gene mapping and discovery in field crops and trees. Thus far, more than 50 SNP arrays and 15 different types of genotyping-by-sequencing (GBS) platforms have been developed in over 25 crop species and perennial trees. However, much less effort has been made on developing ultra-high-throughput and cost-effective genotyping platforms for applied breeding programs. In this review, we discuss the scientific bottlenecks in existing SNP arrays and GBS technologies and the strategies to develop targeted platforms for crop molecular breeding. We propose that future practical breeding platforms should adopt automated genotyping technologies, either array or sequencing based, target functional polymorphisms underpinning economic traits, and provide desirable prediction accuracy for quantitative traits, with universal applications under wide genetic backgrounds in crops. The development of such platforms faces serious challenges at both the technological level due to cost ineffectiveness, and the knowledge level due to large genotype-phenotype gaps in crop plants. It is expected that such genotyping platforms will be achieved in the next ten years in major crops in consideration of (a) rapid development in gene discovery of important traits, (b) deepened understanding of quantitative traits through new analytical models and population designs, (c) integration of multi-layer -omics data leading to identification of genes and pathways responsible for important breeding traits, and (d) improvement in cost effectiveness of large-scale genotyping. Crop breeding chips and genotyping platforms will provide unprecedented opportunities to accelerate the development of cultivars with desired yield potential, quality, and enhanced adaptation to mitigate the effects of climate change.

Integrated high-density consensus genetic map of Pyrus and anchoring of the 'Bartlett' v1.0 (Pyrus communis) genome

Genetic maps are essential tools for pear genetics and genomics research. In this study, we first constructed an integrated simple sequence repeat (SSR) and single nucleotide polymorphism (SNP)-based consensus genetic map for pear based on common SSR markers between nine published maps. A total of 5,085 markers, including 1,232 SSRs and 3,853 SNPs, were localized on a consensus map spanning 3,266.0 cM in total, with an average marker interval of 0.64 cM, which represents the highest density consensus map of pear to date. Using three sets of high-density SNP-based genetic maps with European pear genetic backgrounds, we anchored a total of 291.5 Mb of the ‘Bartlett’ v1.0 (Pyrus communis L.) genome scaffolds into 17 pseudo-chromosomes. This accounted for 50.5% of the genome assembly, which was a great improvement on the 29.7% achieved originally. Intra-genome and inter-genome synteny analyses of the new ‘Bartlett’ v1.1 genome assembly with the Asian pear ‘Dangshansuli’ (Pyrus bretschneideri Rehd.) and apple (Malus × domestica Borkh.) genomes uncovered four new segmental duplication regions. The integrated high-density SSR and SNP-based consensus genetic map provided new insights into the genetic structure patterns of pear and assisted in the genome assembly of ‘Bartlett’ through further exploration of different pear genetic maps.

Genome mapping of postzygotic hybrid necrosis in an interspecific pear population

Deleterious epistatic interactions in plant inter- and intraspecific hybrids can cause a phenomenon known as hybrid necrosis, characterized by a typical seedling phenotype whose main distinguishing features are dwarfism, tissue necrosis and in some cases lethality. Identification of the chromosome regions associated with this type of incompatibility is important not only to increase our understanding of the evolutionary diversification that led to speciation but also for breeding purposes. Development of molecular markers linked to the lethal genes will allow breeders to avoid incompatible inbred combinations that could affect the expression of important agronomic tratis co-segregating with these genes. Although hybrid necrosis has been reported in several plant taxa, including Rosaceae species, this phenomenon has not been described previously in pear. In the interspecific pear population resulting from a cross between PEAR3 (Pyrus bretschneideri × Pyrus communis) and 'Moonglow' (P. communis), we observed two types of hybrid necrosis, expressed at different stages of plant development. Using a combination of previously mapped and newly developed genetic markers, we identified three chromosome regions associated with these two types of lethality, which were genetically independent. One type resulted from a negative epistatic interaction between a locus on linkage group 5 (LG5) of PEAR3 and a locus on LG1 of 'Moonglow', while the second type was due to a gene that maps to LG2 of PEAR3 and which either acts alone or more probably interacts with another gene of unknown location inherited from 'Moonglow'.

Genome-Wide Sequence Variation Identification and Floral-Associated Trait Comparisons Based on the Re-sequencing of the 'Nagafu No. 2' and 'Qinguan' Varieties of Apple (Malus domestica Borkh.)

Apple (Malus domestica Borkh.) is a commercially important fruit worldwide. Detailed information on genomic DNA polymorphisms, which are important for understanding phenotypic traits, is lacking for the apple. We re-sequenced two elite apple varieties, 'Nagafu No. 2' and 'Qinguan,' which have different characteristics. We identified many genomic variations, including 2,771,129 single nucleotide polymorphisms (SNPs), 82,663 structural variations (SVs), and 1,572,803 insertion/deletions (INDELs) in 'Nagafu No. 2' and 2,262,888 SNPs, 63,764 SVs, and 1,294,060 INDELs in 'Qinguan.' The 'SNP,' 'INDEL,' and 'SV' distributions were non-random, with variation-rich or -poor regions throughout the genomes. In 'Nagafu No. 2' and 'Qinguan' there were 171,520 and 147,090 non-synonymous SNPs spanning 23,111 and 21,400 genes, respectively; 3,963 and 3,196 SVs in 3,431 and 2,815 genes, respectively; and 1,834 and 1,451 INDELs in 1,681 and 1,345 genes, respectively. Genetic linkage maps of 190 flowering genes associated with multiple flowering pathways in 'Nagafu No. 2,' 'Qinguan,' and 'Golden Delicious,' identified complex regulatory mechanisms involved in floral induction, flower bud formation, and flowering characteristics, which might reflect the genetic variation of the flowering genes. Expression profiling of key flowering genes in buds and leaves suggested that the photoperiod and autonomous flowering pathways are major contributors to the different floral-associated traits between 'Nagafu No. 2' and 'Qinguan.' The genome variation data provided a foundation for the further exploration of apple diversity and gene-phenotype relationships, and for future research on molecular breeding to improve apple and related species.

Genomics of pear and other Rosaceae fruit trees

The family Rosaceae includes many economically important fruit trees, such as pear, apple, peach, cherry, quince, apricot, plum, raspberry, and loquat. Over the past few years, whole-genome sequences have been released for Chinese pear, European pear, apple, peach, Japanese apricot, and strawberry. These sequences help us to conduct functional and comparative genomics studies and to develop new cultivars with desirable traits by marker-assisted selection in breeding programs. These genomics resources also allow identification of evolutionary relationships in Rosaceae, development of genome-wide SNP and SSR markers, and construction of reference genetic linkage maps, which are available through the Genome Database for the Rosaceae website. Here, we review the recent advances in genomics studies and their practical applications for Rosaceae fruit trees, particularly pear, apple, peach, and cherry.

Genetic control of pear rootstock-induced dwarfing and precocity is linked to a chromosomal region syntenic to the apple Dw1 loci

BACKGROUND

The vigour and precocity of trees highly influences their efficiency in commercial production. In apple, dwarfing rootstocks allow high-density plantings while their precocious flowering enables earlier fruit production. Currently, there is a lack of pear (Pyrus communis L.) rootstocks that are equivalent to the high yielding apple rootstock 'M9'. For the efficient breeding of new Pyrus rootstocks it is crucial to understand the genetic determinants of vigour control and precocity. In this study we used quantitative trait loci (QTLs) analysis to identify genetic loci associated with the desired traits, using a segregating population of 405 F1 P. communis seedlings from a cross between 'Old Home' and 'Louise Bonne de Jersey' (OHxLBJ). The seedlings were grafted as rootstocks with 'Doyenne du Comice' scions and comprehensively phenotyped over four growing seasons for traits related to tree architecture and flowering, in order to describe the growth of the scions.

RESULTS

A high density single nucleotide polymorphism (SNP)-based genetic map comprising 597 polymorphic pear and 113 apple markers enabled the detection of QTLs influencing expression of scion vigour and precocity located on linkage groups (LG)5 and LG6 of 'Old Home'. The LG5 QTL maps to a position that is syntenic to the apple 'Malling 9' ('M9') Dw1 locus at the upper end of LG5. An allele of a simple sequence repeat (SSR) associated with apple Dw1 segregated with dwarfing and precocity in pear and was identified in other pear germplasm accessions. The orthology of the vigour-controlling LG5 QTL between apple and pear raises the possibility that the dwarfing locus Dw1 arose before the divergence of apple and pear, and might therefore be present in other Rosaceae species.

CONCLUSION

We report the first QTLs associated with vigour control and flowering traits in pear rootstocks. Orthologous loci were found to control scion growth and precocity in apple and pear rootstocks. The application of our results may assist in the breeding process of a pear rootstock that confers both vigour control and precocity to the grafted scion cultivar.

Evaluating the Influence of the Microsatellite Marker Set on the Genetic Structure Inferred in Pyrus communis L

Fingerprinting information can be used to elucidate in a robust manner the genetic structure of germplasm collections, allowing a more rational and fine assessment of genetic resources. Bayesian model-based approaches are nowadays majorly preferred to infer genetic structure, but it is still largely unresolved how marker sets should be built in order to obtain a robust inference. The objective was to evaluate, in Pyrus germplasm collections, the influence of the SSR marker set size on the genetic structure inferred, also evaluating the influence of the criterion used to select those markers. Inferences were performed considering an increasing number of SSR markers that ranged from just two up to 25, incorporated one at a time into the analysis. The influence of the number of SSR markers used was evaluated comparing the number of populations and the strength of the signal detected, and also the similarity of the genotype assignments to populations between analyses. In order to test if those results were influenced by the criterion used to select the SSRs, several choosing scenarios based on the discrimination power or the fixation index values of the SSRs were tested. Our results indicate that population structure could be inferred accurately once a certain SSR number threshold was reached, which depended on the underlying structure within the genotypes, but the method used to select the markers included on each set appeared not to be very relevant. The minimum number of SSRs required to provide robust structure inferences and adequate measurements of the differentiation, even when low differentiation levels exist within populations, was proved similar to that of the complete list of recommended markers for fingerprinting. When a SSR set size similar to the minimum marker sets recommended for fingerprinting it is used, only major divisions or moderate (FST>0.05) differentiation of the germplasm are detected.

A genetic linkage map of black raspberry (Rubus occidentalis) and the mapping of Ag(4) conferring resistance to the aphid Amphorophora agathonica

We have constructed a densely populated, saturated genetic linkage map of black raspberry and successfully placed a locus for aphid resistance. Black raspberry (Rubus occidentalis L.) is a high-value crop in the Pacific Northwest of North America with an international marketplace. Few genetic resources are readily available and little improvement has been achieved through breeding efforts to address production challenges involved in growing this crop. Contributing to its lack of improvement is low genetic diversity in elite cultivars and an untapped reservoir of genetic diversity from wild germplasm. In the Pacific Northwest, where most production is centered, the current standard commercial cultivar is highly susceptible to the aphid Amphorophora agathonica Hottes, which is a vector for the Raspberry mosaic virus complex. Infection with the virus complex leads to a rapid decline in plant health resulting in field replacement after only 3-4 growing seasons. Sources of aphid resistance have been identified in wild germplasm and are used to develop mapping populations to study the inheritance of these valuable traits. We have constructed a genetic linkage map using single-nucleotide polymorphism and transferable (primarily simple sequence repeat) markers for F1 population ORUS 4305 consisting of 115 progeny that segregate for aphid resistance. Our linkage map of seven linkage groups representing the seven haploid chromosomes of black raspberry consists of 274 markers on the maternal map and 292 markers on the paternal map including a morphological locus for aphid resistance. This is the first linkage map of black raspberry and will aid in developing markers for marker-assisted breeding, comparative mapping with other Rubus species, and enhancing the black raspberry genome assembly.

Development and preliminary evaluation of a 90 K Axiom® SNP array for the allo-octoploid cultivated strawberry Fragaria × ananassa

BACKGROUND

A high-throughput genotyping platform is needed to enable marker-assisted breeding in the allo-octoploid cultivated strawberry Fragaria × ananassa. Short-read sequences from one diploid and 19 octoploid accessions were aligned to the diploid Fragaria vesca 'Hawaii 4' reference genome to identify single nucleotide polymorphisms (SNPs) and indels for incorporation into a 90 K Affymetrix® Axiom® array. We report the development and preliminary evaluation of this array.

RESULTS

About 36 million sequence variants were identified in a 19 member, octoploid germplasm panel. Strategies and filtering pipelines were developed to identify and incorporate markers of several types: di-allelic SNPs (66.6%), multi-allelic SNPs (1.8%), indels (10.1%), and ploidy-reducing "haploSNPs" (11.7%). The remaining SNPs included those discovered in the diploid progenitor F. iinumae (3.9%), and speculative "codon-based" SNPs (5.9%). In genotyping 306 octoploid accessions, SNPs were assigned to six classes with Affymetrix's "SNPolisher" R package. The highest quality classes, PolyHigh Resolution (PHR), No Minor Homozygote (NMH), and Off-Target Variant (OTV) comprised 25%, 38%, and 1% of array markers, respectively. These markers were suitable for genetic studies as demonstrated in the full-sib family 'Holiday' × 'Korona' with the generation of a genetic linkage map consisting of 6,594 PHR SNPs evenly distributed across 28 chromosomes with an average density of approximately one marker per 0.5 cM, thus exceeding our goal of one marker per cM.

CONCLUSIONS

The Affymetrix IStraw90 Axiom array is the first high-throughput genotyping platform for cultivated strawberry and is commercially available to the worldwide scientific community. The array's high success rate is likely driven by the presence of naturally occurring variation in ploidy level within the nominally octoploid genome, and by effectiveness of the employed array design and ploidy-reducing strategies. This array enables genetic analyses including generation of high-density linkage maps, identification of quantitative trait loci for economically important traits, and genome-wide association studies, thus providing a basis for marker-assisted breeding in this high value crop.

Development and preliminary evaluation of a 90 K Axiom® SNP array for the allo-octoploid cultivated strawberry Fragaria × ananassa

BACKGROUND

A high-throughput genotyping platform is needed to enable marker-assisted breeding in the allo-octoploid cultivated strawberry Fragaria × ananassa. Short-read sequences from one diploid and 19 octoploid accessions were aligned to the diploid Fragaria vesca 'Hawaii 4' reference genome to identify single nucleotide polymorphisms (SNPs) and indels for incorporation into a 90 K Affymetrix® Axiom® array. We report the development and preliminary evaluation of this array.

RESULTS

About 36 million sequence variants were identified in a 19 member, octoploid germplasm panel. Strategies and filtering pipelines were developed to identify and incorporate markers of several types: di-allelic SNPs (66.6%), multi-allelic SNPs (1.8%), indels (10.1%), and ploidy-reducing "haploSNPs" (11.7%). The remaining SNPs included those discovered in the diploid progenitor F. iinumae (3.9%), and speculative "codon-based" SNPs (5.9%). In genotyping 306 octoploid accessions, SNPs were assigned to six classes with Affymetrix's "SNPolisher" R package. The highest quality classes, PolyHigh Resolution (PHR), No Minor Homozygote (NMH), and Off-Target Variant (OTV) comprised 25%, 38%, and 1% of array markers, respectively. These markers were suitable for genetic studies as demonstrated in the full-sib family 'Holiday' × 'Korona' with the generation of a genetic linkage map consisting of 6,594 PHR SNPs evenly distributed across 28 chromosomes with an average density of approximately one marker per 0.5 cM, thus exceeding our goal of one marker per cM.

CONCLUSIONS

The Affymetrix IStraw90 Axiom array is the first high-throughput genotyping platform for cultivated strawberry and is commercially available to the worldwide scientific community. The array's high success rate is likely driven by the presence of naturally occurring variation in ploidy level within the nominally octoploid genome, and by effectiveness of the employed array design and ploidy-reducing strategies. This array enables genetic analyses including generation of high-density linkage maps, identification of quantitative trait loci for economically important traits, and genome-wide association studies, thus providing a basis for marker-assisted breeding in this high value crop.

High-density genetic linkage map construction and identification of fruit-related QTLs in pear using SNP and SSR markers

Pear (Pyrus spp) is an important fruit crop, grown in all temperate regions of the world, with global production ranked after grape and apples among deciduous tree crops. A high-density linkage map is a valuable tool for fine mapping quantitative trait loci (QTL) and map-based gene cloning. In this study, we firstly constructed a high-density linkage map of pear using SNPs integrated with SSRs, developed by the rapid and robust technology of restriction-associated DNA sequencing (RADseq). The linkage map consists of 3143 SNP markers and 98 SSRs, 3241 markers in total, spanning 2243.4 cM, with an average marker distance of 0.70 cM. Anchoring SSRs were able to anchor seventeen linkage groups to their corresponding chromosomes. Based on this high-density integrated pear linkage map and two years of fruit phenotyping, a total of 32 potential QTLs for 11 traits, including length of pedicel (LFP), single fruit weight (SFW), soluble solid content (SSC), transverse diameter (TD), vertical diameter (VD), calyx status (CS), flesh colour (FC), juice content (JC), number of seeds (NS), skin colour (SC), and skin smooth (SS), were identified and positioned on the genetic map. Among them, some important fruit-related traits have for the first time been identified, such as calyx status, length of pedicel, and flesh colour, and reliable localization of QTLs were verified repeatable. This high-density linkage map of pear is a worthy reference for mapping important fruit traits, QTL identification, and comparison and combination of different genetic maps.

Development and validation of a 20K single nucleotide polymorphism (SNP) whole genome genotyping array for apple (Malus × domestica Borkh)

High-density SNP arrays for genome-wide assessment of allelic variation have made high resolution genetic characterization of crop germplasm feasible. A medium density array for apple, the IRSC 8K SNP array, has been successfully developed and used for screens of bi-parental populations. However, the number of robust and well-distributed markers contained on this array was not sufficient to perform genome-wide association analyses in wider germplasm sets, or Pedigree-Based Analysis at high precision, because of rapid decay of linkage disequilibrium. We describe the development of an Illumina Infinium array targeting 20K SNPs. The SNPs were predicted from re-sequencing data derived from the genomes of 13 Malus × domestica apple cultivars and one accession belonging to a crab apple species (M. micromalus). A pipeline for SNP selection was devised that avoided the pitfalls associated with the inclusion of paralogous sequence variants, supported the construction of robust multi-allelic SNP haploblocks and selected up to 11 entries within narrow genomic regions of ±5 kb, termed focal points (FPs). Broad genome coverage was attained by placing FPs at 1 cM intervals on a consensus genetic map, complementing them with FPs to enrich the ends of each of the chromosomes, and by bridging physical intervals greater than 400 Kbps. The selection also included ∼3.7K validated SNPs from the IRSC 8K array. The array has already been used in other studies where ∼15.8K SNP markers were mapped with an average of ∼6.8K SNPs per full-sib family. The newly developed array with its high density of polymorphic validated SNPs is expected to be of great utility for Pedigree-Based Analysis and Genomic Selection. It will also be a valuable tool to help dissect the genetic mechanisms controlling important fruit quality traits, and to aid the identification of marker-trait associations suitable for the application of Marker Assisted Selection in apple breeding programs.

Genetic, metabolite and developmental determinism of fruit friction discolouration in pear

BACKGROUND

The unattractive appearance of the surface of pear fruit caused by the postharvest disorder friction discolouration (FD) is responsible for significant consumer dissatisfaction in markets, leading to lower returns to growers. Developing an understanding of the genetic control of FD is essential to enable the full application of genomics-informed breeding for the development of new pear cultivars. Biochemical constituents [phenolic compounds and ascorbic acid (AsA)], polyphenol oxidase (PPO) activity, as well as skin anatomy, have been proposed to play important roles in FD susceptibility in studies on a limited number of cultivars. However, to date there has been no investigation on the biochemical and genetic control of FD, employing segregating populations. In this study, we used 250 seedlings from two segregating populations (POP369 and POP356) derived from interspecific crosses between Asian (Pyrus pyrifolia Nakai and P. bretschneideri Rehd.) and European (P. communis) pears to identify genetic factors associated with susceptibility to FD.

RESULTS

Single nucleotide polymorphism (SNP)-based linkage maps suitable for QTL analysis were developed for the parents of both populations. The maps for population POP369 comprised 174 and 265 SNP markers for the male and female parent, respectively, while POP356 maps comprised 353 and 398 SNP markers for the male and female parent, respectively. Phenotypic data for 22 variables were measured over two successive years (2011 and 2012) for POP369 and one year (2011) only for POP356. A total of 221 QTLs were identified that were linked to 22 phenotyped variables, including QTLs associated with FD for both populations that were stable over the successive years. In addition, clear evidence of the influence of developmental factors (fruit maturity) on FD and other variables was also recorded.

CONCLUSIONS

The QTLs associated with fruit firmness, PPO activity, AsA concentration and concentration of polyphenol compounds as well as FD are the first reported for pear. We conclude that the postharvest disorder FD is controlled by multiple small effect QTLs and that it will be very challenging to apply marker-assisted selection based on these QTLs. However, genomic selection could be employed to select elite genotypes with lower or no susceptibility to FD early in the breeding cycle.

The draft genome sequence of European pear (Pyrus communis L. 'Bartlett')

We present a draft assembly of the genome of European pear (Pyrus communis) 'Bartlett'. Our assembly was developed employing second generation sequencing technology (Roche 454), from single-end, 2 kb, and 7 kb insert paired-end reads using Newbler (version 2.7). It contains 142,083 scaffolds greater than 499 bases (maximum scaffold length of 1.2 Mb) and covers a total of 577.3 Mb, representing most of the expected 600 Mb Pyrus genome. A total of 829,823 putative single nucleotide polymorphisms (SNPs) were detected using re-sequencing of 'Louise Bonne de Jersey' and 'Old Home'. A total of 2,279 genetically mapped SNP markers anchor 171 Mb of the assembled genome. Ab initio gene prediction combined with prediction based on homology searching detected 43,419 putative gene models. Of these, 1219 proteins (556 clusters) are unique to European pear compared to 12 other sequenced plant genomes. Analysis of the expansin gene family provided an example of the quality of the gene prediction and an insight into the relationships among one class of cell wall related genes that control fruit softening in both European pear and apple (Malus × domestica). The 'Bartlett' genome assembly v1.0 (http://www.rosaceae.org/species/pyrus/pyrus_communis/genome_v1.0) is an invaluable tool for identifying the genetic control of key horticultural traits in pear and will enable the wide application of marker-assisted and genomic selection that will enhance the speed and efficiency of pear cultivar development.