Chromosome-scale genome sequence assemblies of the 'Autumn Bliss' and 'Malling Jewel' cultivars of the highly heterozygous red raspberry (Rubus idaeus L.) derived from long-read Oxford Nanopore sequence data.

Red raspberry (Rubus idaeus L.) is an economically valuable soft-fruit species with a relatively small (~300 Mb) but highly heterozygous diploid (2n = 2x = 14) genome. Chromosome-scale genome sequences are a vital tool in unravelling the genetic complexity controlling traits of interest in crop plants such as red raspberry, as well as for functional genomics, evolutionary studies, and pan-genomics diversity studies. In this study, we developed genome sequences of a primocane fruiting variety ('Autumn Bliss') and a floricane variety ('Malling Jewel'). The use of long-read Oxford Nanopore Technologies sequencing data yielded long read lengths that permitted well resolved genome sequences for the two cultivars to be assembled. The de novo assemblies of 'Malling Jewel' and 'Autumn Bliss' contained 79 and 136 contigs respectively, and 263.0 Mb of the 'Autumn Bliss' and 265.5 Mb of the 'Malling Jewel' assembly could be anchored unambiguously to a previously published red raspberry genome sequence of the cultivar 'Anitra'. Single copy ortholog analysis (BUSCO) revealed high levels of completeness in both genomes sequenced, with 97.4% of sequences identified in 'Autumn Bliss' and 97.7% in 'Malling Jewel'. The density of repetitive sequence contained in the 'Autumn Bliss' and 'Malling Jewel' assemblies was significantly higher than in the previously published assembly and centromeric and telomeric regions were identified in both assemblies. A total of 42,823 protein coding regions were identified in the 'Autumn Bliss' assembly, whilst 43,027 were identified in the 'Malling Jewel' assembly. These chromosome-scale genome sequences represent an excellent genomics resource for red raspberry, particularly around the highly repetitive centromeric and telomeric regions of the genome that are less complete in the previously published 'Anitra' genome sequence.

SSR-Based Analysis of Genetic Diversity and Structure of Sweet Cherry (Prunus avium L.) from 19 Countries in Europe.

Sweet cherry (Prunus avium L.) is a temperate fruit species whose production might be highly impacted by climate change in the near future. Diversity of plant material could be an option to mitigate these climate risks by enabling producers to have new cultivars well adapted to new environmental conditions. In this study, subsets of sweet cherry collections of 19 European countries were genotyped using 14 SSR. The objectives of this study were (i) to assess genetic diversity parameters, (ii) to estimate the levels of population structure, and (iii) to identify germplasm redundancies. A total of 314 accessions, including landraces, early selections, and modern cultivars, were monitored, and 220 unique SSR genotypes were identified. All 14 loci were confirmed to be polymorphic, and a total of 137 alleles were detected with a mean of 9.8 alleles per locus. The average number of alleles (N = 9.8), PIC value (0.658), observed heterozygosity (Ho = 0.71), and expected heterozygosity (He = 0.70) were higher in this study compared to values reported so far. Four ancestral populations were detected using STRUCTURE software and confirmed by Principal Coordinate Analysis (PCoA), and two of them (K1 and K4) could be attributed to the geographical origin of the accessions. A N-J tree grouped the 220 sweet cherry accessions within three main clusters and six subgroups. Accessions belonging to the four STRUCTURE populations roughly clustered together. Clustering confirmed known genealogical data for several accessions. The large genetic diversity of the collection was demonstrated, in particular within the landrace pool, justifying the efforts made over decades for their conservation. New sources of diversity will allow producers to face challenges, such as climate change and the need to develop more sustainable production systems.

Towards a Joint International Database: Alignment of SSR Marker Data for European Collections of Cherry Germplasm.

The objective of our study was the alignment of microsatellite or simple sequence repeat (SSR) marker data across germplasm collections of cherry within Europe. Through the European Cooperative program for Plant Genetic Resources ECPGR, a number of European germplasm collections had previously been analysed using standard sets of SSR loci. However, until now these datasets remained unaligned. We used a combination of standard reference genotypes and ad-hoc selections to compile a central dataset representing as many alleles as possible from national datasets produced in France, Great Britain, Germany, Italy, Sweden and Switzerland. Through the comparison of alleles called in data from replicated samples we were able to create a series of alignment factors, supported across 448 different allele calls, that allowed us to align a dataset of 2241 SSR profiles from six countries. The proportion of allele comparisons that were either in agreement with the alignment factor or confounded by null alleles ranged from 67% to 100% and this was further improved by the inclusion of a series of allele-specific adjustments. The aligned dataset allowed us to identify groups of previously unknown matching accessions and to identify and resolve a number of errors in the prior datasets. The combined and aligned dataset represents a significant step forward in the co-ordinated management of field collections of cherry in Europe.

Soil microbiome data of two apple orchards in the UK.

The microbial communities in two apple orchards were characterised using amplicon-based metabarcoding. Samples were taken from tree station locations along a linear transect and from adjacent grass aisles, at both orchards. Comparison was made between the communities occurring at tree station locations and the grass aisles, and between orchards. Further discussion of these datasets is given in https://doi.org/10.1016/j.apsoil.2018.05.015 (Deakin et al., 2018).

Spatial structuring of soil microbial communities in commercial apple orchards.

Characterising spatial microbial community structure is important to understand and explain the consequences of continuous plantation of one crop species on the performance of subsequent crops, especially where this leads to reduced growth vigour and crop yield. We investigated the spatial structure, specifically distance-decay of similarity, of soil bacterial and fungal communities in two long-established orchards with contrasting agronomic characteristics. A spatially explicit sampling strategy was used to collect soil from under recently grubbed rows of apple trees and under the grassed aisles. Amplicon-based metabarcoding technology was used to characterise the soil microbial communities. The results suggested that (1) most of the differences in soil microbial community structure were due to large-scale differences (i.e. between orchards), (2) within-orchard, small-scale (1-5 m) spatial variability was also present, but spatial relationships in microbial community structure differed between orchards and were not predictable, and (3) vegetation type (i.e. trees or grass and their associated management) can significantly alter the structure of soil microbial communities, affecting a large proportion of microbial groups. The discontinuous nature of soil microbial community structure in the tree stations and neighbouring grass aisles within an orchard illustrate the importance of vegetation type and allied weed and nutrient management on soil microbial community structure.

Analysis of the genetic diversity and structure across a wide range of germplasm reveals prominent gene flow in apple at the European level.

BACKGROUND: The 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. RESULTS: A 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. CONCLUSIONS: The 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.

A new three-locus model for rootstock-induced dwarfing in apple revealed by genetic mapping of root bark percentage.

Rootstock-induced dwarfing of apple scions revolutionized global apple production during the twentieth century, leading to the development of modern intensive orchards. A high root bark percentage (the percentage of the whole root area constituted by root cortex) has previously been associated with rootstock-induced dwarfing in apple. In this study, the root bark percentage was measured in a full-sib family of ungrafted apple rootstocks and found to be under the control of three loci. Two quantitative trait loci (QTLs) for root bark percentage were found to co-localize to the same genomic regions on chromosome 5 and chromosome 11 previously identified as controlling dwarfing, Dw1 and Dw2, respectively. A third QTL was identified on chromosome 13 in a region that has not been previously associated with dwarfing. The development of closely linked sequence-tagged site markers improved the resolution of allelic classes, thereby allowing the detection of dominance and epistatic interactions between loci, with high root bark percentage only occurring in specific allelic combinations. In addition, we report a significant negative correlation between root bark percentage and stem diameter (an indicator of tree vigour), measured on a clonally propagated grafted subset of the mapping population. The demonstrated link between root bark percentage and rootstock-induced dwarfing of the scion leads us to propose a three-locus model that is able to explain levels of dwarfing from the dwarf 'M.27' to the semi-invigorating rootstock 'M.116'. Moreover, we suggest that the QTL on chromosome 13 (Rb3) might be analogous to a third dwarfing QTL, Dw3, which has not previously been identified.

Contributions of roots and rootstocks to sustainable, intensified crop production.

Sustainable intensification is seen as the main route for meeting the world's increasing demands for food and fibre. As demands mount for greater efficiency in the use of resources to achieve this goal, so the focus on roots and rootstocks and their role in acquiring water and nutrients, and overcoming pests and pathogens, is increasing. The purpose of this review is to explore some of the ways in which understanding root systems and their interactions with soils could contribute to the development of more sustainable systems of intensive production. Physical interactions with soil particles limit root growth if soils are dense, but root-soil contact is essential for optimal growth and uptake of water and nutrients. X-ray microtomography demonstrated that maize roots elongated more rapidly with increasing root-soil contact, as long as mechanical impedance was not limiting root elongation, while lupin was less sensitive to changes in root-soil contact. In addition to selecting for root architecture and rhizosphere properties, the growth of many plants in cultivated systems is profoundly affected by selection of an appropriate rootstock. Several mechanisms for scion control by rootstocks have been suggested, but the causal signals are still uncertain and may differ between crop species. Linkage map locations for quantitative trait loci for disease resistance and other traits of interest in rootstock breeding are becoming available. Designing root systems and rootstocks for specific environments is becoming a feasible target.

Saturated linkage map construction in Rubus idaeus using genotyping by sequencing and genome-independent imputation.

BACKGROUND: Rapid development of highly saturated genetic maps aids molecular breeding, which can accelerate gain per breeding cycle in woody perennial plants such as Rubus idaeus (red raspberry). Recently, robust genotyping methods based on high-throughput sequencing were developed, which provide high marker density, but result in some genotype errors and a large number of missing genotype values. Imputation can reduce the number of missing values and can correct genotyping errors, but current methods of imputation require a reference genome and thus are not an option for most species. RESULTS: Genotyping by Sequencing (GBS) was used to produce highly saturated maps for a R. idaeus pseudo-testcross progeny. While low coverage and high variance in sequencing resulted in a large number of missing values for some individuals, a novel method of imputation based on maximum likelihood marker ordering from initial marker segregation overcame the challenge of missing values, and made map construction computationally tractable. The two resulting parental maps contained 4521 and 2391 molecular markers spanning 462.7 and 376.6 cM respectively over seven linkage groups. Detection of precise genomic regions with segregation distortion was possible because of map saturation. Microsatellites (SSRs) linked these results to published maps for cross-validation and map comparison. CONCLUSIONS: GBS together with genome-independent imputation provides a rapid method for genetic map construction in any pseudo-testcross progeny. Our method of imputation estimates the correct genotype call of missing values and corrects genotyping errors that lead to inflated map size and reduced precision in marker placement. Comparison of SSRs to published R. idaeus maps showed that the linkage maps constructed with GBS and our method of imputation were robust, and marker positioning reliable. The high marker density allowed identification of genomic regions with segregation distortion in R. idaeus, which may help to identify deleterious alleles that are the basis of inbreeding depression in the species.

Development of a dense SNP-based linkage map of an apple rootstock progeny using the Malus Infinium whole genome genotyping array.

BACKGROUND: A whole-genome genotyping array has previously been developed for Malus using SNP data from 28 Malus genotypes. This array offers the prospect of high throughput genotyping and linkage map development for any given Malus progeny. To test the applicability of the array for mapping in diverse Malus genotypes, we applied the array to the construction of a SNP-based linkage map of an apple rootstock progeny. RESULTS: Of the 7,867 Malus SNP markers on the array, 1,823 (23.2%) were heterozygous in one of the two parents of the progeny, 1,007 (12.8%) were heterozygous in both parental genotypes, whilst just 2.8% of the 921 Pyrus SNPs were heterozygous. A linkage map spanning 1,282.2 cM was produced comprising 2,272 SNP markers, 306 SSR markers and the S-locus. The length of the M432 linkage map was increased by 52.7 cM with the addition of the SNP markers, whilst marker density increased from 3.8 cM/marker to 0.5 cM/marker. Just three regions in excess of 10 cM remain where no markers were mapped. We compared the positions of the mapped SNP markers on the M432 map with their predicted positions on the 'Golden Delicious' genome sequence. A total of 311 markers (13.7% of all mapped markers) mapped to positions that conflicted with their predicted positions on the 'Golden Delicious' pseudo-chromosomes, indicating the presence of paralogous genomic regions or mis-assignments of genome sequence contigs during the assembly and anchoring of the genome sequence. CONCLUSIONS: We incorporated data for the 2,272 SNP markers onto the map of the M432 progeny and have presented the most complete and saturated map of the full 17 linkage groups of M. pumila to date. The data were generated rapidly in a high-throughput semi-automated pipeline, permitting significant savings in time and cost over linkage map construction using microsatellites. The application of the array will permit linkage maps to be developed for QTL analyses in a cost-effective manner, and the identification of SNPs that have been assigned erroneous positions on the 'Golden Delicious' reference sequence will assist in the continued improvement of the genome sequence assembly for that variety.

Fine genetic mapping of the Co locus controlling columnar growth habit in apple.

Tree architecture is an important, complex and dynamic trait affected by diverse genetic, ontogenetic and environmental factors. 'Wijcik McIntosh', a columnar (reduced branching) sport of 'McIntosh' and a valuable genetic resource, has been used intensively in apple-breeding programs for genetic improvement of tree architecture. The columnar growth habit is primarily controlled by the dominant allele of gene Co (columnar) on linkage group-10. But the Co locus is not well mapped and the Co gene remains unknown. To precisely map the Co locus and to identify candidate genes of Co, a sequence-based approach using both peach and apple genomes was used to develop new markers linked more tightly to Co. Five new simple sequence repeats markers were developed (C1753-3520, C18470-25831, C6536-31519, C7223-38004 and C7629-22009). The first four markers were obtained from apple genomic sequences on chromosome-10, whereas the last (C7629-22009) was from an unanchored apple contig that contains an apple expressed sequence tag CV082943, which was identified through synteny analysis between the peach and apple genomes. Genetic mapping of these five markers in four F(1) populations of 528 genotypes and 290 diverse columnar selections/cultivars (818 genotypes in total) delimited the Co locus in a genetic interval with 0.37 % recombination between markers C1753-3520 and C7629-22009. Marker C18470-25831 co-segregates with Co in the 818 genotypes studied. The Co region is estimated to be 193 kb and contains 26 predicted gene in the 'Golden Delicious' genome. Among the 26 genes, three are putative LATERAL ORGAN BOUNDARIES (LOB) DOMAIN (LBD) containing transcription factor genes known of essential roles in plant lateral organ development, and are therefore considered as strong candidates of Co, designated MdLBD1, MdLBD2, and MdLBD3. Although more comprehensive studies are required to confirm the function of MdLBD1-3, the present work represents an important step forward to better understand the genetic and molecular control of tree architecture in apple.

A new gene for resistance to Dysaphis pyri in pear and identification of flanking microsatellite markers.

Dysaphis pyri is an important aphid pest of European pear (Pyrus communis) cultivars, none of which are currently reported to be resistant. In this study, we produced a progeny of the European pear Comice crossed with an accession of snow pear (Pyrus nivalis) that segregated for resistance to D. pyri in a Mendelian fashion, indicating the presence of a major gene, Dp-1. Following screening of the parents and seedlings with microsatellite markers, cosegregation analysis indicated that Dp-1 is flanked by NH006b and NH014a on linkage group 17, 2.3 and 3.6 cM away, respectively. Evidence is also presented for the duplication of linkage groups 9 and 17, which is a consequence of the allopolyploid origin of pear.

Mapping of A1 conferring resistance to the aphid Amphorophora idaei and dw (dwarfing habit) in red raspberry (Rubus idaeus L.) using AFLP and microsatellite markers.

BACKGROUND: Raspberry breeding programmes worldwide aim to produce improved cultivars to satisfy market demands and within these programmes there are many targets, including increased fruit quality, yield and season, and improved pest and disease resistance and plant habit. The large raspberry aphid, Amphorophora idaei, transmits four viruses and vector resistance is an objective in raspberry breeding. The development of molecular tools that discriminate between aphid resistance genes from different sources will allow the pyramiding of such genes and the development of raspberry varieties with superior pest resistance. We have raised a red raspberry (Rubus idaeus) F1 progeny from the cross 'Malling Jewel' x 'Malling Orion' (MJ x MO), which segregates for resistance to biotype 1 of the aphid Amphorophora idaei and for a second phenotypic trait, dwarf habit. These traits are controlled by single genes, denoted (A1) and (dw) respectively. RESULTS: The progeny of 94 seedlings was scored for the segregation of 95 AFLP and 22 SSR markers and a linkage map was constructed that covers a total genetic distance of 505 cM over seven linkage groups. The average linkage group length was 72.2 cM and there was an average of 17 markers per linkage group, of which at least two were codominant SSRs, allowing comparisons with previously published maps of raspberry. The two phenotypic traits, A1 and dw, mapped to linkage groups 3 and 6 respectively. CONCLUSION: The mapping of A1 will facilitate the discrimination of resistance genes from different sources and the pyramiding of aphid resistance genes in new raspberry cultivars; the mapping of dw will allow further investigations into the genetics of dwarfing habit in Rubus.