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.

A chromosome-length genome assembly and annotation of blackberry (Rubus argutus, cv. "Hillquist").

Blackberries (Rubus spp.) are the fourth most economically important berry crop worldwide. Genome assemblies and annotations have been developed for Rubus species in subgenus Idaeobatus, including black raspberry (R. occidentalis), red raspberry (R. idaeus), and R. chingii, but very few genomic resources exist for blackberries and their relatives in subgenus Rubus. Here we present a chromosome-length assembly and annotation of the diploid blackberry germplasm accession "Hillquist" (R. argutus). "Hillquist" is the only known source of primocane-fruiting (annual-fruiting) in tetraploid fresh-market blackberry breeding programs and is represented in the pedigree of many important cultivars worldwide. The "Hillquist" assembly, generated using Pacific Biosciences long reads scaffolded with high-throughput chromosome conformation capture sequencing, consisted of 298 Mb, of which 270 Mb (90%) was placed on 7 chromosome-length scaffolds with an average length of 38.6 Mb. Approximately 52.8% of the genome was composed of repetitive elements. The genome sequence was highly collinear with a novel maternal haplotype-resolved linkage map of the tetraploid blackberry selection A-2551TN and genome assemblies of R. chingii and red raspberry. A total of 38,503 protein-coding genes were predicted, of which 72% were functionally annotated. Eighteen flowering gene homologs within a previously mapped locus aligning to an 11.2 Mb region on chromosome Ra02 were identified as potential candidate genes for primocane-fruiting. The utility of the "Hillquist" genome has been demonstrated here by the development of the first genotyping-by-sequencing-based linkage map of tetraploid blackberry and the identification of possible candidate genes for primocane-fruiting. This chromosome-length assembly will facilitate future studies in Rubus biology, genetics, and genomics and strengthen applied breeding programs.

A chromosome-level genome sequence assembly of the red raspberry (Rubus idaeus L.).

Rubus idaeus L. (red raspberry), is a perennial woody plant species of the Rosaceae family that is widely cultivated in the temperate regions of world and is thus an economically important soft fruit species. It is prized for its flavour and aroma, as well as a high content of healthful compounds such as vitamins and antioxidants. Breeding programs exist globally for red raspberry, but variety development is a long and challenging process. Genomic and molecular tools for red raspberry are valuable resources for breeding. Here, a chromosome-length genome sequence assembly and related gene predictions for the red raspberry cultivar 'Anitra' are presented, comprising PacBio long read sequencing scaffolded using Hi-C sequence data. The assembled genome sequence totalled 291.7 Mbp, with 247.5 Mbp (84.8%) incorporated into seven sequencing scaffolds with an average length of 35.4 Mbp. A total of 39,448 protein-coding genes were predicted, 75% of which were functionally annotated. The seven chromosome scaffolds were anchored to a previously published genetic linkage map with a high degree of synteny and comparisons to genomes of closely related species within the Rosoideae revealed chromosome-scale rearrangements that have occurred over relatively short evolutionary periods. A chromosome-level genomic sequence of R. idaeus will be a valuable resource for the knowledge of its genome structure and function in red raspberry and will be a useful and important resource for researchers and plant breeders.

De novo sequencing and analysis of the transcriptome of two highbush blueberry (Vaccinium corymbosum L.) cultivars 'Bluecrop' and 'Legacy' at harvest and following post-harvest storage.

Fruit firmness and in particular the individual components of texture and moisture loss, are considered the key quality traits when describing blueberry fruit quality, and whilst these traits are genetically regulated, the mechanisms governing their control are not clearly understood. In this investigation, RNAseq was performed on fruits of two blueberry cultivars with very different storage properties, 'Bluecrop' and 'Legacy', at harvest, three weeks storage in a non-modified environment at 4 °C and after three weeks storage at 4 °C followed by three days at 21 °C, with the aim of understanding the transcriptional changes that occur during storage in cultivars with very different post-harvest fruit quality. De novo assemblies of the transcriptomes of the two cultivars were performed separately and a total of 39,335 and 41,896 unigenes for 'Bluecrop' and 'Legacy' respectively were resolved. Differential gene expression analyses were grouped into four cluster profiles based on changes in transcript abundance between harvest and 24 days post-harvest. A total of 290 unigenes were up-regulated in 'Legacy' only, 685 were up-regulated in 'Bluecrop', 252 were up-regulated in both cultivars and 948 were down-regulated in both cultivars between harvest and 24 days post-harvest. Unigenes showing significant differential expression between harvest and following post-harvest cold-storage were grouped into classes of biological processes including stress responses, cell wall metabolism, wax metabolism, calcium metabolism, cellular components, and biological processes. In total 21 differentially expressed unigenes with a putative role in regulating the response to post-harvest cold-storage in the two cultivars were identified from the de novo transcriptome assemblies performed. The results presented provide a stable foundation from which to perform further analyses with which to functionally validate the candidate genes identified, and to begin to understand the genetic mechanisms controlling changes in firmness in blueberry fruits post-harvest.

Genetic mapping and identification of a QTL determining tolerance to freezing stress in Fragaria vesca L.

Extreme cold and frost cause significant stress to plants which can potentially be lethal. Low temperature freezing stress can cause significant and irreversible damage to plant cells and can induce physiological and metabolic changes that impact on growth and development. Low temperatures cause physiological responses including winter dormancy and autumn cold hardening in strawberry (Fragaria) species, and some diploid F. vesca accessions have been shown to have adapted to low-temperature stresses. To study the genetics of freezing tolerance, a F. vesca mapping population of 143 seedlings segregating for differential responses to freezing stress was raised. The progeny was mapped using 'Genotyping-by-Sequencing' and a linkage map of 2,918 markers at 851 loci was resolved. The mapping population was phenotyped for freezing tolerance response under controlled and replicated laboratory conditions and subsequent quantitative trait loci analysis using interval mapping revealed a single significant quantitative trait locus on Fvb2 in the physical interval 10.6 Mb and 15.73 Mb on the F. vesca v4.0 genome sequence. This physical interval contained 896 predicted genes, several of which had putative roles associated with tolerance to abiotic stresses including freezing. Differential expression analysis of the 896 QTL-associated gene predictions in the leaves and crowns from 'Alta' and 'NCGR1363' parental genotypes revealed genotype-specific changes in transcript accumulation in response to low temperature treatment as well as expression differences between genotypes prior to treatment for many of the genes. The putative roles, and significant interparental differential expression levels of several of the genes reported here identified them as good candidates for the control of the effects of freezing tolerance at the QTL identified in this investigation and the possible role of these candidate genes in response to freezing stress is discussed.

Major-effect candidate genes identified in cultivated strawberry (Fragaria × ananassa Duch.) for ellagic acid deoxyhexoside and pelargonidin-3-O-malonylglucoside biosynthesis, key polyphenolic compounds.

Strawberries are rich in polyphenols which impart health benefits when metabolized by the gut microbiome, including anti-inflammatory, neuroprotective, and antiproliferative effects. In addition, polyphenolic anthocyanins contribute to the attractive color of strawberry fruits. However, the genetic basis of polyphenol biosynthesis has not been extensively studied in strawberry. In this investigation, ripe fruits from three cultivated strawberry populations were characterized for polyphenol content using HPLC-DAD-MSn and genotyped using the iStraw35k array. GWAS and QTL analyses identified genetic loci controlling polyphenol biosynthesis. QTL were identified on four chromosomes for pelargonidin-3-O-malonylglucoside, pelargonidin-3-O-acetylglucoside, cinnamoyl glucose, and ellagic acid deoxyhexoside biosynthesis. Presence/absence of ellagic acid deoxyhexoside and pelargonidin-3-O-malonylglucoside was found to be under the control of major gene loci on LG1X2 and LG6b, respectively, on the F. × ananassa linkage maps. Interrogation of gene predictions in the F. vesca reference genome sequence identified a single candidate gene for ellagic acid deoxyhexoside biosynthesis, while seven malonyltransferase genes were identified as candidates for pelargonidin-3-O-malonylglucoside biosynthesis. Homologous malonyltransferase genes were identified in the F. × ananassa 'Camarosa' genome sequence but the candidate for ellagic acid deoxyhexoside biosynthesis was absent from the 'Camarosa' sequence. This study demonstrated that polyphenol biosynthesis in strawberry is, in some cases, under simple genetic control, supporting previous observations of the presence or absence of these compounds in strawberry fruits. It has also shed light on the mechanisms controlling polyphenol biosynthesis and enhanced the knowledge of these biosynthesis pathways in strawberry. The above findings will facilitate breeding for strawberries enriched in compounds with beneficial health effects.

The genome sequence and transcriptome of Potentilla micrantha and their comparison to Fragaria vesca (the woodland strawberry).

Background: The genus Potentilla is closely related to that of Fragaria, the economically important strawberry genus. Potentilla micrantha is a species that does not develop berries but shares numerous morphological and ecological characteristics with Fragaria vesca. These similarities make P. micrantha an attractive choice for comparative genomics studies with F. vesca. Findings: In this study, the P. micrantha genome was sequenced and annotated, and RNA-Seq data from the different developmental stages of flowering and fruiting were used to develop a set of gene predictions. A 327 Mbp sequence and annotation of the genome of P. micrantha, spanning 2674 sequence contigs, with an N50 size of 335,712, estimated to cover 80% of the total genome size of the species was developed. The genus Potentilla has a characteristically larger genome size than Fragaria, but the recovered sequence scaffolds were remarkably collinear at the micro-syntenic level with the genome of F. vesca, its closest sequenced relative. A total of 33,602 genes were predicted, and 95.1% of bench-marking universal single-copy orthologous genes were complete within the presented sequence. Thus, we argue that the majority of the gene-rich regions of the genome have been sequenced. Conclusions: Comparisons of RNA-Seq data from the stages of floral and fruit development revealed genes differentially expressed between P. micrantha and F. vesca.The data presented are a valuable resource for future studies of berry development in Fragaria and the Rosaceae and they also shed light on the evolution of genome size and organization in this family.

Genotyping-by-sequencing enables linkage mapping in three octoploid cultivated strawberry families.

Genotyping-by-sequencing (GBS) was used to survey genome-wide single-nucleotide polymorphisms (SNPs) in three biparental strawberry (Fragaria × ananassa) populations with the goal of evaluating this technique in a species with a complex octoploid genome. GBS sequence data were aligned to the F. vesca 'Fvb' reference genome in order to call SNPs. Numbers of polymorphic SNPs per population ranged from 1,163 to 3,190. Linkage maps consisting of 30-65 linkage groups were produced from the SNP sets derived from each parent. The linkage groups covered 99% of the Fvb reference genome, with three to seven linkage groups from a given parent aligned to any particular chromosome. A phylogenetic analysis performed using the POLiMAPS pipeline revealed linkage groups that were most similar to ancestral species F. vesca for each chromosome. Linkage groups that were most similar to a second ancestral species, F. iinumae, were only resolved for Fvb 4. The quantity of missing data and heterogeneity in genome coverage inherent in GBS complicated the analysis, but POLiMAPS resolved F. × ananassa chromosomal regions derived from diploid ancestor F. vesca.

Additive QTLs on three chromosomes control flowering time in woodland strawberry (Fragaria vesca L.).

Flowering time is an important trait that affects survival, reproduction and yield in both wild and cultivated plants. Therefore, many studies have focused on the identification of flowering time quantitative trait locus (QTLs) in different crops, and molecular control of this trait has been extensively investigated in model species. Here we report the mapping of QTLs for flowering time and vegetative traits in a large woodland strawberry mapping population that was phenotyped both under field conditions and in a greenhouse after flower induction in the field. The greenhouse experiment revealed additive QTLs in three linkage groups (LG), two on both LG4 and LG7, and one on LG6 that explain about half of the flowering time variance in the population. Three of the QTLs were newly identified in this study, and one co-localized with the previously characterized FvTFL1 gene. An additional strong QTL corresponding to previously mapped PFRU was detected in both field and greenhouse experiments indicating that gene(s) in this locus can control the timing of flowering in different environments in addition to the duration of flowering and axillary bud differentiation to runners and branch crowns. Several putative flowering time genes were identified in these QTL regions that await functional validation. Our results indicate that a few major QTLs may control flowering time and axillary bud differentiation in strawberries. We suggest that the identification of causal genes in the diploid strawberry may enable fine tuning of flowering time and vegetative growth in the closely related octoploid cultivated strawberry.

A ddRAD Based Linkage Map of the Cultivated Strawberry, Fragaria xananassa.

The cultivated strawberry (Fragaria ×ananassa Duch.) is an allo-octoploid considered difficult to disentangle genetically due to its four relatively similar sub-genomic chromosome sets. This has been alleviated by the recent release of the strawberry IStraw90 whole genome genotyping array. However, array resolution relies on the genotypes used in the array construction and may be of limited general use. SNP detection based on reduced genomic sequencing approaches has the potential of providing better coverage in cases where the studied genotypes are only distantly related from the SNP array's construction foundation. Here we have used double digest restriction-associated DNA sequencing (ddRAD) to identify SNPs in a 145 seedling F1 hybrid population raised from the cross between the cultivars Sonata (♀) and Babette (♂). A linkage map containing 907 markers which spanned 1,581.5 cM across 31 linkage groups representing the 28 chromosomes of the species. Comparing the physical span of the SNP markers with the F. vesca genome sequence, the linkage groups resolved covered 79% of the estimated 830 Mb of the F. × ananassa genome. Here, we have developed the first linkage map for F. × ananassa using ddRAD and show that this technique and other related techniques are useful tools for linkage map development and downstream genetic studies in the octoploid strawberry.

Species-specific duplications driving the recent expansion of NBS-LRR genes in five Rosaceae species.

BACKGROUND: Disease resistance (R) genes from different Rosaceae species have been identified by map-based cloning for resistance breeding. However, there are few reports describing the pattern of R-gene evolution in Rosaceae species because several Rosaceae genome sequences have only recently become available. RESULTS: Since most disease resistance genes encode NBS-LRR proteins, we performed a systematic genome-wide survey of NBS-LRR genes between five Rosaceae species, namely Fragaria vesca (strawberry), Malus × domestica (apple), Pyrus bretschneideri (pear), Prunus persica (peach) and Prunus mume (mei) which contained 144, 748, 469, 354 and 352 NBS-LRR genes, respectively. A high proportion of multi-genes and similar Ks peaks (Ks = 0.1- 0.2) of gene families in the four woody genomes were detected. A total of 385 species-specific duplicate clades were observed in the phylogenetic tree constructed using all 2067 NBS-LRR genes. High percentages of NBS-LRR genes derived from species-specific duplication were found among the five genomes (61.81% in strawberry, 66.04% in apple, 48.61% in pear, 37.01% in peach and 40.05% in mei). Furthermore, the Ks and Ka/Ks values of TIR-NBS-LRR genes (TNLs) were significantly greater than those of non-TIR-NBS-LRR genes (non-TNLs), and most of the NBS-LRRs had Ka/Ks ratios less than 1, suggesting that they were evolving under a subfunctionalization model driven by purifying selection. CONCLUSIONS: Our results indicate that recent duplications played an important role in the evolution of NBS-LRR genes in the four woody perennial Rosaceae species. Based on the phylogenetic tree produced, it could be inferred that species-specific duplication has mainly contributed to the expansion of NBS-LRR genes in the five Rosaceae species. In addition, the Ks and Ka/Ks ratios suggest that the rapidly evolved TNLs have different evolutionary patterns to adapt to different pathogens compared with non-TNL resistant genes.

Looking forward to genetically edited fruit crops.

The availability of genome sequences for many fruit crops has redefined the boundaries of genetic engineering and genetically modified (GM) crop plants. However commercialization of GM crops is hindered by numerous regulatory and social hurdles. Here, we focus on recently developed genome-editing tools for fruit crop improvement and their importance from the consumer perspective. Challenges and opportunities for the deployment of new genome-editing tools for fruit plants are also discussed.

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.

Evaluation of SNP Data from the Malus Infinium Array Identifies Challenges for Genetic Analysis of Complex Genomes of Polyploid Origin.

High throughput arrays for the simultaneous genotyping of thousands of single-nucleotide polymorphisms (SNPs) have made the rapid genetic characterisation of plant genomes and the development of saturated linkage maps a realistic prospect for many plant species of agronomic importance. However, the correct calling of SNP genotypes in divergent polyploid genomes using array technology can be problematic due to paralogy, and to divergence in probe sequences causing changes in probe binding efficiencies. An Illumina Infinium II whole-genome genotyping array was recently developed for the cultivated apple and used to develop a molecular linkage map for an apple rootstock progeny (M432), but a large proportion of segregating SNPs were not mapped in the progeny, due to unexpected genotype clustering patterns. To investigate the causes of this unexpected clustering we performed BLAST analysis of all probe sequences against the 'Golden Delicious' genome sequence and discovered evidence for paralogous annealing sites and probe sequence divergence for a high proportion of probes contained on the array. Following visual re-evaluation of the genotyping data generated for 8,788 SNPs for the M432 progeny using the array, we manually re-scored genotypes at 818 loci and mapped a further 797 markers to the M432 linkage map. The newly mapped markers included the majority of those that could not be mapped previously, as well as loci that were previously scored as monomorphic, but which segregated due to divergence leading to heterozygosity in probe annealing sites. An evaluation of the 8,788 probes in a diverse collection of Malus germplasm showed that more than half the probes returned genotype clustering patterns that were difficult or impossible to interpret reliably, highlighting implications for the use of the array in genome-wide association studies.

Whole genome comparisons of Fragaria, Prunus and Malus reveal different modes of evolution between Rosaceous subfamilies.

BACKGROUND: Rosaceae include numerous economically important and morphologically diverse species. Comparative mapping between the member species in Rosaceae have indicated some level of synteny. Recently the whole genome of three crop species, peach, apple and strawberry, which belong to different genera of the Rosaceae family, have been sequenced, allowing in-depth comparison of these genomes. RESULTS: Our analysis using the whole genome sequences of peach, apple and strawberry identified 1399 orthologous regions between the three genomes, with a mean length of around 100 kb. Each peach chromosome showed major orthology mostly to one strawberry chromosome, but to more than two apple chromosomes, suggesting that the apple genome went through more chromosomal fissions in addition to the whole genome duplication after the divergence of the three genera. However, the distribution of contiguous ancestral regions, identified using the multiple genome rearrangements and ancestors (MGRA) algorithm, suggested that the Fragaria genome went through a greater number of small scale rearrangements compared to the other genomes since they diverged from a common ancestor. Using the contiguous ancestral regions, we reconstructed a hypothetical ancestral genome for the Rosaceae 7 composed of nine chromosomes and propose the evolutionary steps from the ancestral genome to the extant Fragaria, Prunus and Malus genomes. CONCLUSION: Our analysis shows that different modes of evolution may have played major roles in different subfamilies of Rosaceae. The hypothetical ancestral genome of Rosaceae and the evolutionary steps that lead to three different lineages of Rosaceae will facilitate our understanding of plant genome evolution as well as have a practical impact on knowledge transfer among member species of Rosaceae.

Comparative genetic mapping between octoploid and diploid Fragaria species reveals a high level of colinearity between their genomes and the essentially disomic behavior of the cultivated octoploid strawberry.

Macrosynteny and colinearity between Fragaria (strawberry) species showing extreme levels of ploidy have been studied through comparative genetic mapping between the octoploid cultivated strawberry (F. xananassa) and its diploid relatives. A comprehensive map of the octoploid strawberry, in which almost all linkage groups are ranged into the seven expected homoeologous groups was obtained, thus providing the first reference map for the octoploid Fragaria. High levels of conserved macrosynteny and colinearity were observed between homo(eo)logous linkage groups and between the octoploid homoeologous groups and their corresponding diploid linkage groups. These results reveal that the polyploidization events that took place along the evolution of the Fragaria genus and the more recent juxtaposition of two octoploid strawberry genomes in the cultivated strawberry did not trigger any major chromosomal rearrangements in genomes involved in F. xananassa. They further suggest the existence of a close relationship between the diploid Fragaria genomes. In addition, despite the possible existence of residual levels of polysomic segregation suggested by the observation of large linkage groups in coupling phase only, the prevalence of linkage groups in coupling/repulsion phase clearly demonstrates that the meiotic behavior is mainly disomic in the cultivated strawberry.