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

Development of a next generation SNP genotyping array for wheat

High-throughput genotyping arrays have provided a cost-effective, reliable and interoperable system for genotyping hexaploid wheat and its relatives. Existing, highly cited arrays including our 35K Wheat Breeder's array and the Illumina 90K array were designed based on a limited amount of varietal sequence diversity and with imperfect knowledge of SNP positions. Recent progress in wheat sequencing has given us access to a vast pool of SNP diversity, whilst technological improvements have allowed us to fit significantly more probes onto a 384-well format Axiom array than previously possible. Here we describe a novel Axiom genotyping array, the 'Triticum aestivum Next Generation' array (TaNG), largely derived from whole genome skim sequencing of 204 elite wheat lines and 111 wheat landraces taken from the Watkins 'Core Collection'. We used a novel haplotype optimization approach to select SNPs with the highest combined varietal discrimination and a design iteration step to test and replace SNPs which failed to convert to reliable markers. The final design with 43 372 SNPs contains a combination of haplotype-optimized novel SNPs and legacy cross-platform markers. We show that this design has an improved distribution of SNPs compared to previous arrays and can be used to generate genetic maps with a significantly higher number of distinct bins than our previous array. We also demonstrate the improved performance of TaNGv1.1 for Genome-wide association studies (GWAS) and its utility for Copy Number Variation (CNV) analysis. The array is commercially available with supporting marker annotations and initial genotyping results freely available.

Diversity analyses in two ornamental and large-genome Ranunculaceae species based on a low-cost Klenow NGS-based protocol

Persian buttercup (Ranunculus asiaticus L.) and poppy anemone (Anemone coronaria L.) are ornamental, outcrossing, perennial species belonging to the Ranunculaceae family, characterized by large and highly repetitive genomes. We applied K-seq protocol in both species to generate high-throughput sequencing data and produce a large number of genetic polymorphisms. The technique entails the application of Klenow polymerase-based PCR using short primers designed by analyzing k-mer sets in the genome sequence. To date the genome sequence of both species has not been released, thus we designed primer sets based on the reference the genome sequence of the related species Aquilegia oxysepala var. kansuensis (Brühl). A whole of 11,542 SNPs were selected for assessing genetic diversity of eighteen commercial varieties of R. asiaticus, while 1,752 SNPs for assessing genetic diversity in six cultivars of A. coronaria. UPGMA dendrograms were constructed and in R. asiaticus integrated in with PCA analysis. This study reports the first molecular fingerprinting within Persian buttercup, while the results obtained in poppy anemone were compared with a previously published SSR-based fingerprinting, proving K-seq to be an efficient protocol for the genotyping of complex genetic backgrounds.

Development of a 135K SNP genotyping array for Actinidia arguta and its applications for genetic mapping and QTL analysis in kiwifruit

Kiwifruit (Actinidia spp) is a woody, perennial and deciduous vine. In this genus, there are multiple ploidy levels but the main cultivated cultivars are polyploid. Despite the availability of many genomic resources in kiwifruit, SNP genotyping is still a challenge given these different levels of polyploidy. Recent advances in SNP array technologies have offered a high-throughput genotyping platform for genome-wide DNA polymorphisms. In this study, we developed a high-density SNP genotyping array to facilitate genetic studies and breeding applications in kiwifruit. SNP discovery was performed by genome-wide DNA sequencing of 40 kiwifruit genotypes. The identified SNPs were stringently filtered for sequence quality, predicted conversion performance and distribution over the available Actinidia chinensis genome. A total of 134 729 unique SNPs were put on the array. The array was evaluated by genotyping 400 kiwifruit individuals. We performed a multidimensional scaling analysis to assess the diversity of kiwifruit germplasm, showing that the array was effective to distinguish kiwifruit accessions. Using a tetraploid F1 population, we constructed an integrated linkage map covering 3060.9 cM across 29 linkage groups and performed QTL analysis for the sex locus that has been identified on Linkage Group 3 (LG3) in Actinidia arguta. Finally, our dataset presented evidence of tetrasomic inheritance with partial preferential pairing in A. arguta. In conclusion, we developed and evaluated a 135K SNP genotyping array for kiwifruit. It has the advantage of a comprehensive design that can be an effective tool in genetic studies and breeding applications in this high-value crop.

First genetic maps development and QTL mining in Ranunculus asiaticus L. through ddRADseq

Persian Buttercup (Ranunculus asiaticus L.; 2x=2n=16; estimated genome size: 7.6Gb) is an ornamental and perennial crop native of Asia Minor and Mediterranean basin, marketed both as cut flower or potted plant. Currently new varieties are developed by selecting plants carrying desirable traits in segregating progenies obtained by controlled mating, which are propagated through rhizomes or micro-propagated in vitro. In order to escalate selection efficiency and respond to market requests, more knowledge of buttercup genetics would facilitate the identification of markers associated with loci and genes controlling key ornamental traits, opening the way for molecular assisted breeding programs. Reduced-representation sequencing (RRS) represents a powerful tool for plant genotyping, especially in case of large genomes such as the one of buttercup, and have been applied for the development of high-density genetic maps in several species. We report on the development of the first molecular-genetic maps in R. asiaticus based on of a two-way pseudo-testcross strategy. A double digest restriction-site associated DNA (ddRAD) approach was applied for genotyping two F₁ mapping populations, whose female parents were a genotype of a so called 'ponpon' and of a 'double flower' varieties, while the common male parental ('Cipro') was a genotype producing a simple flower. The ddRAD generated a total of ~2Gb demultiplexed reads, resulting in an average of 8,3M reads per line. The sstacks pipeline was applied for the construction of a mock reference genome based on sequencing data, and SNP markers segregating in only one of the parents were retained for map construction by treating the F₁ population as a backcross. The four parental maps (two of the female parents and two of the common male parent) were aligned with 106 common markers and 8 linkage groups were identified, corresponding to the haploid chromosome number of the species. An average of 586 markers were associated with each parental map, with a marker density ranging from 1 marker/cM to 4.4 markers/cM. The developed maps were used for QTL analysis for flower color, leading to the identification of major QTLs for purple pigmentation. These results contribute to dissect on the genetics of Persian buttercup, enabling the development of new approaches for future varietal development.

Morphological studies of rose prickles provide new insights

Prickles are common structures in plants that play a key role in defense against herbivores. In the Rosa genus, prickles are widely present with great diversity in terms of form and density. For cut rose production, prickles represent an important issue, as they can damage the flower and injure workers. Our objectives were to precisely describe the types of prickles that exist in roses, their tissues of origin and their development. We performed a detailed histological analysis of prickle initiation and development in a rose F1 population. Based on the prickle investigation of 110 roses, we proposed the first categorization of prickles in the Rosa genus. They are mainly divided into two categories, nonglandular prickles (NGPs) and glandular prickles (GPs), and subcategories were defined based on the presence/absence of hairs and branches. We demonstrated that NGPs and GPs both originate from multiple cells of the ground meristem beneath the protoderm. For GPs, the gland cells originate from the protoderm of the GP at the early developmental stage. Our findings clearly demonstrate that prickles are not modified trichomes (which originate from the protoderm). These conclusions are different from the current mainstream hypothesis. These results provide a foundation for further studies on prickle initiation and development in plants.

Analysis of allelic variants of RhMLO genes in rose and functional studies on susceptibility to powdery mildew related to clade V homologs

Rose has 19 MLO genes. Of these, RhMLO1 and RhMLO2 were shown to be required for powdery mildew infection, which suggests their potential as susceptibility targets towards disease resistance. Powdery mildew, caused by Podosphaera pannosa, is one of the most serious and widespread fungal diseases for roses, especially in greenhouse-grown cut roses. It has been shown that certain MLO genes are involved in powdery mildew susceptibility and that loss of function in these genes in various crops leads to broad-spectrum, long-lasting resistance against this fungal disease. For this reason, these MLO genes are called susceptibility genes. We carried out a genome-wide identification of the MLO gene family in the Rosa chinensis genome, and screened for allelic variants among 22 accessions from seven different Rosa species using re-sequencing and transcriptome data. We identified 19 MLO genes in rose, of which four are candidate genes for functional homologs in clade V, which is the clade containing all dicot MLO susceptibility genes. We detected a total of 198 different allelic variants in the set of Rosa species and accessions, corresponding to 5-15 different alleles for each of the genes. Some diploid Rosa species shared alleles with tetraploid rose cultivars, consistent with the notion that diploid species have contributed to the formation of tetraploid roses. Among the four RhMLO genes in clade V, we demonstrated using expression study, virus-induced gene silencing as well as transient RNAi silencing that two of them, RhMLO1 and RhMLO2, are required for infection by P. pannosa and suggest their potential as susceptibility targets for powdery mildew resistance breeding in rose.

Identification and QTL Analysis of Flavonoids and Carotenoids in Tetraploid Roses Based on an Ultra-High-Density Genetic Map

Roses are highly valuable within the flower industry. The metabolites of anthocyanins, flavonols, and carotenoids in rose petals are not only responsible for the various visible petal colors but also important bioactive compounds that are important for human health. In this study, we performed a QTL analysis on pigment contents to locate major loci that determine the flower color traits. An F1 population of tetraploid roses segregating for flower color was used to construct an ultra-high-density genetic linkage map using whole-genome resequencing technology to detect genome-wide SNPs. Previously developed SSR and SNP markers were also utilized to increase the marker density. Thus, a total of 9,259 markers were mapped onto seven linkage groups (LGs). The final length of the integrated map was 1285.11 cM, with an average distance of 0.14 cM between adjacent markers. The contents of anthocyanins, flavonols and carotenoids of the population were assayed to enable QTL analysis. Across the 33 components, 46 QTLs were detected, explaining 11.85-47.72% of the phenotypic variation. The mapped QTLs were physically clustered and primarily distributed on four linkage groups, namely LG2, LG4, LG6, and LG7. These results improve the basis for flower color marker-assisted breeding of tetraploid roses and guide the development of rose products.

Genetic Mapping in Autohexaploid Sweet Potato with Low-Coverage NGS-Based Genotyping Data

Next-generation sequencing (NGS)-based genotyping methods can generate numerous genetic markers in a single experiment and have contributed to plant genetic mapping. However, for high precision genetic analysis, the complicated genetic segregation mode in polyploid organisms requires high-coverage NGS data and elaborate analytical algorithms. In the present study, we propose a simple strategy for the genetic mapping of polyploids using low-coverage NGS data. The validity of the strategy was investigated using simulated data. Previous studies indicated that accurate allele dosage estimation from low-coverage NGS data (read depth < 40) is difficult. Therefore, we used allele dosage probabilities calculated from read counts in association analyses to detect loci associated with phenotypic variations. The allele dosage probabilities showed significant detection power, although higher allele dosage estimation accuracy resulted in higher detection power. On the contrary, differences in the segregation patterns between the marker and causal genes resulted in a drastic decrease in detection power even if the marker and casual genes were in complete linkage and the allele dosage estimation was accurate. These results indicated that the use of a larger number of markers is advantageous, even if the accuracy of allele dosage estimation is low. Finally, we applied the strategy for the genetic mapping of autohexaploid sweet potato (Ipomoea batatas) populations to detect loci associated with agronomic traits. Our strategy could constitute a cost-effective approach for preliminary experiments done performed to large-scale studies.

Molecular markers from the chloroplast genome of rose provide a complementary tool for variety discrimination and profiling

The rose is one of the most important ornamental woody plants because of its extensive use and high economic value. Herein, we sequenced a complete chloroplast genome of the miniature rose variety Rosa 'Margo Koster' and performed comparative analyses with sequences previously published for other species in the Rosaceae family. The chloroplast genome of Rosa 'Margo Koster', with a size of 157,395 bp, has a circular quadripartite structure typical of angiosperm chloroplast genomes and contains a total of 81 protein-coding genes, 30 tRNA genes and 4 rRNA genes. Conjunction regions in the chloroplast genome of Rosa 'Margo Koster' were verified and manually corrected by Sanger sequencing. Comparative genome analysis showed that the IR contraction and expansion events resulted in rps19 and ycf1 pseudogenes. The phylogenetic analysis within the Rosa genus showed that Rosa 'Margo Koster' is closer to Rosa odorata than to other Rosa species. Additionally, we identified and screened highly divergent sequences and cpSSRs and compared their power to discriminate rose varieties by Sanger sequencing and capillary electrophoresis. The results showed that 15 cpSSRs are polymorphic, but their discriminating power is only moderate among a set of rose varieties. However, more than 150 single nucleotide variations (SNVs) were discovered in the flanking region of cpSSRs, and the results indicated that these SNVs have a higher divergence and stronger power for profiling rose varieties. These findings suggest that nucleotide mutations in the chloroplast genome may be an effective and powerful tool for rose variety discrimination and DNA profiling. These molecular markers in the chloroplast genome sequence of Rosa spp. will facilitate population and phylogenetic studies and other related studies of this species.

Molecular mapping and identification of quantitative trait loci for domestication traits in the field cress (Lepidium campestre L.) genome

Lepidium campestre (L.) or field cress is a multifaceted oilseed plant, which is not yet domesticated. Moreover, the molecular and genetic mechanisms underlying the domestication traits of field cress remain largely elusive. The overarching goal of this study is to identify quantitative trait loci (QTL) that are fundamental for domestication of field cress. Mapping and dissecting quantitative trait variation may provide important insights into genomic trajectories underlying field cress domestication. We used 7624 single nucleotide polymorphism (SNP) markers for QTL mapping in 428 F₂ interspecific hybrid individuals, while field phenotyping was conducted in F_2:3 segregating families. We applied multiple QTL mapping algorithms to detect and estimate the QTL effects for seven important domestication traits of field cress. Verification of pod shattering across sites revealed that the non-shattering lines declined drastically whereas the shattering lines increased sharply, possibly due to inbreeding followed by selection events. In total, 1461 of the 7624 SNP loci were mapped to eight linkage groups (LGs), spanning 571.9 cM map length. We identified 27 QTL across all LGs of field cress genome, which captured medium to high heritability, implying that genomics-assisted selection could deliver domesticated lines in field cress breeding. The use of high throughput genotyping can accelerate the process of domestication in novel crop species. This is the first QTL mapping analysis in the field cress genome that may lay a foundational framework for positional or functional QTL cloning, introgression as well as genomics-assisted breeding in field cress domestication.

Field cress genome mapping: Integrating linkage and comparative maps with cytogenetic analysis for rDNA carrying chromosomes

Field cress (Lepidium campestre L.), despite its potential as a sustainable alternative oilseed plant, has been underutilized, and no prior attempts to characterize the genome at the genetic or molecular cytogenetic level have been conducted. Genetic maps are the foundation for anchoring and orienting annotated genome assemblies and positional cloning of candidate genes. Our principal goal was to construct a genetic map using integrated approaches of genetic, comparative and cytogenetic map analyses. In total, 503 F2 interspecific hybrid individuals were genotyped using 7,624 single nucleotide polymorphism markers. Comparative analysis demonstrated that ~57% of the sequenced loci in L. campestre were congruent with Arabidopsis thaliana (L.) genome and suggested a novel karyotype, which predates the ancestral crucifer karyotype. Aceto-orcein chromosome staining and fluorescence in situ hybridization (FISH) analyses confirmed that L. campestre, L. heterophyllum Benth. and their hybrids had a chromosome number of 2n = 2x = 16. Flow cytometric analysis revealed that both species possess 2C roughly 0.4 picogram DNA. Integrating linkage and comparative maps with cytogenetic map analyses assigned two linkage groups to their particular chromosomes. Future work could incorporate FISH utilizing A. thaliana mapped BAC clones to allow the chromosomes of field cress to be identified reliably.

polymapR-linkage analysis and genetic map construction from F1 populations of outcrossing polyploids

Motivation

Polyploid species carry more than two copies of each chromosome, a condition found in many of the world’s most important crops. Genetic mapping in polyploids is more complex than in diploid species, resulting in a lack of available software tools. These are needed if we are to realize all the opportunities offered by modern genotyping platforms for genetic research and breeding in polyploid crops.

Results

polymapR is an R package for genetic linkage analysis and integrated genetic map construction from bi-parental populations of outcrossing autopolyploids. It can currently analyse triploid, tetraploid and hexaploid marker datasets and is applicable to various crops including potato, leek, alfalfa, blueberry, chrysanthemum, sweet potato or kiwifruit. It can detect, estimate and correct for preferential chromosome pairing, and has been tested on high-density marker datasets from potato, rose and chrysanthemum, generating high-density integrated linkage maps in all of these crops.

Availability and implementation

polymapR is freely available under the general public license from the Comprehensive R Archive Network (CRAN) at http://cran.r-project.org/package=polymapR.

Supplementary information

Supplementary data are available at Bioinformatics online.

Current achievements and future prospects in the genetic breeding of chrysanthemum: a review

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a leading flower with applied value worldwide. Developing new chrysanthemum cultivars with novel characteristics such as new flower colors and shapes, plant architectures, flowering times, postharvest quality, and biotic and abiotic stress tolerance in a time- and cost-efficient manner is the ultimate goal for breeders. Various breeding strategies have been employed to improve the aforementioned traits, ranging from conventional techniques, including crossbreeding and mutation breeding, to a series of molecular breeding methods, including transgenic technology, genome editing, and marker-assisted selection (MAS). In addition, the recent extensive advances in high-throughput technologies, especially genomics, transcriptomics, proteomics, metabolomics, and microbiomics, which are collectively referred to as omics platforms, have led to the collection of substantial amounts of data. Integration of these omics data with phenotypic information will enable the identification of genes/pathways responsible for important traits. Several attempts have been made to use emerging molecular and omics methods with the aim of accelerating the breeding of chrysanthemum. However, applying the findings of such studies to practical chrysanthemum breeding remains a considerable challenge, primarily due to the high heterozygosity and polyploidy of the species. This review summarizes the recent achievements in conventional and modern molecular breeding methods and emerging omics technologies and discusses their future applications for improving the agronomic and horticultural characteristics of chrysanthemum.

The development of a high-density genetic map significantly improves the quality of reference genome assemblies for rose

Roses are important woody plants featuring a set of important traits that cannot be investigated in traditional model plants. Here, we used the restriction-site associated DNA sequencing (RAD-seq) technology to develop a high-density linkage map of the backcross progeny (BC1F1) between Rosa chinensis 'Old Blush' (OB) and R. wichuraiana 'Basyes' Thornless' (BT). We obtained 643.63 million pair-end reads and identified 139,834 polymorphic tags that were distributed uniformly in the rose genome. 2,213 reliable markers were assigned to seven linkage groups (LGs). The length of the genetic map was 1,027.425 cM in total with a mean distance of 0.96 cM per marker locus. This new linkage map allowed anchoring an extra of 1.21/23.14 Mb (12.18/44.52%) of the unassembled OB scaffolds to the seven reference pseudo-chromosomes, thus significantly improved the quality of assembly of OB reference genome. We demonstrate that, while this new linkage map shares high collinearity level with strawberry genome, it also features two chromosomal rearrangements, indicating its usefulness as a resource for understanding the evolutionary scenario among Rosaceae genomes. Together with the newly released genome sequences for OB, this linkage map will facilitate the identification of genetic components underpinning key agricultural and biological traits, hence should greatly advance the studies and breeding efforts of rose.

polymapR-linkage analysis and genetic map construction from F1 populations of outcrossing polyploids

Motivation

Polyploid species carry more than two copies of each chromosome, a condition found in many of the world’s most important crops. Genetic mapping in polyploids is more complex than in diploid species, resulting in a lack of available software tools. These are needed if we are to realise all the opportunities offered by modern genotyping platforms for genetic research and breeding in polyploid crops.

Results

polymapR is an R package for genetic linkage analysis and integrated genetic map construction from bi-parental populations of outcrossing autopolyploids. It can currently analyse triploid, tetraploid and hexaploid marker datasets and is applicable to various crops including potato, leek, alfalfa, blueberry, chrysanthemum, sweet potato or kiwifruit. It can detect, estimate and correct for preferential chromosome pairing, and has been tested on high-density marker datasets from potato, rose and chrysanthemum, generating high-density integrated linkage maps in all of these crops.

Availability and Implementation

polymapR is freely available under the general public license from the Comprehensive R Archive Network (CRAN) athttp://cran.r-project.org/packages=polymapR.

Contact

Chris Maliepaard chris.maliepaard@wur.nl or Roeland E. Voorrips roeland.voorrips@wur.nl

In the name of the rose: a roadmap for rose research in the genome era

The recent completion of the rose genome sequence is not the end of a process, but rather a starting point that opens up a whole set of new and exciting activities. Next to a high-quality genome sequence other genomic tools have also become available for rose, including transcriptomics data, a high-density single-nucleotide polymorphism array and software to perform linkage and quantitative trait locus mapping in polyploids. Rose cultivars are highly heterogeneous and diverse. This vast diversity in cultivated roses can be explained through the genetic potential of the genus, introgressions from wild species into commercial tetraploid germplasm and the inimitable efforts of historical breeders. We can now investigate how this diversity can best be exploited and refined in future breeding work, given the rich molecular toolbox now available to the rose breeding community. This paper presents possible lines of research now that rose has entered the genomics era, and attempts to partially answer the question that arises after the completion of any draft genome sequence: 'Now that we have "the" genome, what's next?'. Having access to a genome sequence will allow both (fundamental) scientific and (applied) breeding-orientated questions to be addressed. We outline possible approaches for a number of these questions.

Comparative transcriptomics identifies patterns of selection in roses

BACKGROUND

Roses are important plants for human beings with pivotal economical and biological traits like continuous flowering, flower architecture, color and scent. Due to frequent hybridization and high genome heterozygosity, classification of roses and their relatives remains a big challenge.

RESULTS

Here, to identify potential markers for phylogenetic reconstruction and to reveal the patterns of natural selection in roses, we generated sets of high quality and comprehensive reference transcriptomes for Rosa chinensis 'Old Blush' (OB) and R. wichuriana 'Basye's Thornless' (BT), two species exhibiting contrasted traits of high economical importance. The assembled reference transcriptomes showed transcripts N50 above 2000 bp. Two roses shared about 10,073 transcripts (N50 = 2282 bp), in which a set of 5959 transcripts was conserved within genera of Rosa. Further comparison with species in Rosaceae identified 4447 transcripts being common (Rosaceae-common) in Rosa, Malus, Prunus, Rubus, and Fragaria, while a pool of 164 transcripts being specific for roses (Rosa-specific). Among the Rosaceae-common transcripts, 409 transcripts showed a signature of positive selection and a clustered expression in different tissues. Interestingly, nine of these rapidly evolving genes were related to DNA damage repair and responses to environmental stimulus, a potential associated with genome confliction post hybridization. Coincident with this fast evolution pattern in rose genes, 24 F-box and four TMV resistant proteins were significantly enriched in the Rosa-specific genes.

CONCLUSIONS

We expect that these Rosaceae-common and Rosa-specific transcripts should facilitate the phylogenetic analysis of Rosaceae plants as well as investigations of Rosa-specific biology. The data reported here could provide fundamental genomic tools and knowledge critical for understanding the biology and domestication of roses and for roses breeding.

polymapR-linkage analysis and genetic map construction from F1 populations of outcrossing polyploids

MOTIVATION

Polyploid species carry more than two copies of each chromosome, a condition found in many of the world's most important crops. Genetic mapping in polyploids is more complex than in diploid species, resulting in a lack of available software tools. These are needed if we are to realize all the opportunities offered by modern genotyping platforms for genetic research and breeding in polyploid crops.

RESULTS

polymapR is an R package for genetic linkage analysis and integrated genetic map construction from bi-parental populations of outcrossing autopolyploids. It can currently analyse triploid, tetraploid and hexaploid marker datasets and is applicable to various crops including potato, leek, alfalfa, blueberry, chrysanthemum, sweet potato or kiwifruit. It can detect, estimate and correct for preferential chromosome pairing, and has been tested on high-density marker datasets from potato, rose and chrysanthemum, generating high-density integrated linkage maps in all of these crops.

AVAILABILITY AND IMPLEMENTATION

polymapR is freely available under the general public license from the Comprehensive R Archive Network (CRAN) at http://cran.r-project.org/package=polymapR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome structure of Rosa multiflora, a wild ancestor of cultivated roses

The draft genome sequence of a wild rose (Rosa multiflora Thunb.) was determined using Illumina MiSeq and HiSeq platforms. The total length of the scaffolds was 739,637,845 bp, consisting of 83,189 scaffolds, which was close to the 711 Mbp length estimated by k-mer analysis. N50 length of the scaffolds was 90,830 bp, and extent of the longest was 1,133,259 bp. The average GC content of the scaffolds was 38.9%. After gene prediction, 67,380 candidates exhibiting sequence homology to known genes and domains were extracted, which included complete and partial gene structures. This large number of genes for a diploid plant may reflect heterogeneity of the genome originating from self-incompatibility in R. multiflora. According to CEGMA analysis, 91.9% and 98.0% of the core eukaryotic genes were completely and partially conserved in the scaffolds, respectively. Genes presumably involved in flower color, scent and flowering are assigned. The results of this study will serve as a valuable resource for fundamental and applied research in the rose, including breeding and phylogenetic study of cultivated roses.

Tools for Genetic Studies in Experimental Populations of Polyploids

Polyploid organisms carry more than two copies of each chromosome, a condition rarely tolerated in animals but which occurs relatively frequently in the plant kingdom. One of the principal challenges faced by polyploid organisms is to evolve stable meiotic mechanisms to faithfully transmit genetic information to the next generation upon which the study of inheritance is based. In this review we look at the tools available to the research community to better understand polyploid inheritance, many of which have only recently been developed. Most of these tools are intended for experimental populations (rather than natural populations), facilitating genomics-assisted crop improvement and plant breeding. This is hardly surprising given that a large proportion of domesticated plant species are polyploid. We focus on three main areas: (1) polyploid genotyping; (2) genetic and physical mapping; and (3) quantitative trait analysis and genomic selection. We also briefly review some miscellaneous topics such as the mode of inheritance and the availability of polyploid simulation software. The current polyploid analytic toolbox includes software for assigning marker genotypes (and in particular, estimating the dosage of marker alleles in the heterozygous condition), establishing chromosome-scale linkage phase among marker alleles, constructing (short-range) haplotypes, generating linkage maps, performing genome-wide association studies (GWAS) and quantitative trait locus (QTL) analyses, and simulating polyploid populations. These tools can also help elucidate the mode of inheritance (disomic, polysomic or a mixture of both as in segmental allopolyploids) or reveal whether double reduction and multivalent chromosomal pairing occur. An increasing number of polyploids (or associated diploids) are being sequenced, leading to publicly available reference genome assemblies. Much work remains in order to keep pace with developments in genomic technologies. However, such technologies also offer the promise of understanding polyploid genomes at a level which hitherto has remained elusive.

Dynamic and epistatic QTL mapping reveals the complex genetic architecture of waterlogging tolerance in chrysanthemum

37 unconditional QTLs, 51 conditional QTLs and considerable epistatic QTLs were detected for waterlogging tolerance, and six favourable combinations were selected accelerating the possible application of MAS in chrysanthemum breeding. Chrysanthemum is seriously impacted by soil waterlogging. To determine the genetic characteristics of waterlogging tolerance (WAT) in chrysanthemum, a population of 162 F₁ lines was used to construct a genetic map to identify the dynamic and epistatic quantitative trait loci (QTLs) for four WAT traits: wilting index (WI), dead leaf ratio (DLR), chlorosis score (Score) and membership function value of waterlogging (MFVW). The h _B² for the WAT traits ranged from 0.49 to 0.64, and transgressive segregation was observed in both directions. A total of 37 unconditional consensus QTLs with 5.81-18.21% phenotypic variation explanation (PVE) and 51 conditional consensus QTLs with 5.90-24.56% PVE were detected. Interestingly, three unconditional consensus QTLs were consistently identified across different stages, whereas no conditional consensus QTLs were consistently expressed. In addition, considerable epistatic QTLs, all with PVE values ranging from 0.01 to 8.87%, were detected by a joint analysis of WAT phenotypes. These results illustrated that the QTLs (genes) controlling WAT were environmentally dependent and selectively expressed at different times and indicated that both additive and epistatic effects underlie the inheritance of WAT in chrysanthemum. The findings of the current study provide insights into the complex genetic architecture of WAT, and the identification of favourable alleles represents an important step towards the application of molecular marker-assisted selection (MAS) and QTL pyramiding in chrysanthemum WAT breeding programmes.

Genotyping-by-sequencing application on diploid rose and a resulting high-density SNP-based consensus map

Roses, which have been cultivated for at least 5000 years, are one of the most important ornamental crops in the world. Because of the interspecific nature and high heterozygosity in commercial roses, the genetic resources available for rose are limited. To effectively identify markers associated with QTL controlling important traits, such as disease resistance, abundant markers along the genome and careful phenotyping are required. Utilizing genotyping by sequencing technology and the strawberry genome (Fragaria vesca v2.0.a1) as a reference, we generated thousands of informative single nucleotide polymorphism (SNP) markers. These SNPs along with known bridge simple sequence repeat (SSR) markers allowed us to create the first high-density integrated consensus map for diploid roses. Individual maps were first created for populations J06-20-14-3×"Little Chief" (J14-3×LC), J06-20-14-3×"Vineyard Song" (J14-3×VS) and "Old Blush"×"Red Fairy" (OB×RF) and these maps were linked with 824 SNPs and 13 SSR bridge markers. The anchor SSR markers were used to determine the numbering of the rose linkage groups. The diploid consensus map has seven linkage groups (LGs), a total length of 892.2 cM, and an average distance of 0.25 cM between 3527 markers. By combining three individual populations, the marker density and the reliability of the marker order in the consensus map was improved over a single population map. Extensive synteny between the strawberry and diploid rose genomes was observed. This consensus map will serve as the tool for the discovery of marker-trait associations in rose breeding using pedigree-based analysis. The high level of conservation observed between the strawberry and rose genomes will help further comparative studies within the Rosaceae family and may aid in the identification of candidate genes within QTL regions.

Development and Applications of a High Throughput Genotyping Tool for Polyploid Crops: Single Nucleotide Polymorphism (SNP) Array

Polypoid species play significant roles in agriculture and food production. Many crop species are polyploid, such as potato, wheat, strawberry, and sugarcane. Genotyping has been a daunting task for genetic studies of polyploid crops, which lags far behind the diploid crop species. Single nucleotide polymorphism (SNP) array is considered to be one of, high-throughput, relatively cost-efficient and automated genotyping approaches. However, there are significant challenges for SNP identification in complex, polyploid genomes, which has seriously slowed SNP discovery and array development in polyploid species. Ploidy is a significant factor impacting SNP qualities and validation rates of SNP markers in SNP arrays, which has been proven to be a very important tool for genetic studies and molecular breeding. In this review, we (1) discussed the pros and cons of SNP array in general for high throughput genotyping, (2) presented the challenges of and solutions to SNP calling in polyploid species, (3) summarized the SNP selection criteria and considerations of SNP array design for polyploid species, (4) illustrated SNP array applications in several different polyploid crop species, then (5) discussed challenges, available software, and their accuracy comparisons for genotype calling based on SNP array data in polyploids, and finally (6) provided a series of SNP array design and genotype calling recommendations. This review presents a complete overview of SNP array development and applications in polypoid crops, which will benefit the research in molecular breeding and genetics of crops with complex genomes.

Recent progress in whole genome sequencing, high-density linkage maps, and genomic databases of ornamental plants

Genome information is useful for functional analysis of genes, comparative genomic analysis, breeding of new varieties by marker-assisted selection, and map-based gene isolation. Genome-related research in ornamentals plants has been relatively slow to develop because of their heterozygosity or polyploidy. Advances in analytical instruments, such as next-generation sequencers and information processing technologies have revolutionized biology, and have been applied in a large number and variety of species, including ornamental plants. Recently, high-quality whole genome sequences have been reported in plant genetics and physiology studies of model ornamentals, such as those in genus Petunia and Japanese morning glory (Ipomoea nil). In this review, whole genome sequencing and construction of high-density genetic linkage maps based on SNP markers of ornamentals will be discussed. The databases that store this information for ornamentals are also described.

Genetic dissection of adventitious shoot regeneration in roses by employing genome-wide association studies

We analysed the capacity to regenerate adventitious shoots in 96 rose genotypes and found 88 SNP markers associated with QTLs, some of which are derived from candidate genes for shoot regeneration. In an association panel of 96 rose genotypes previously analysed for petal colour, we conducted a genome-wide association study on the capacity of leaf petioles for direct shoot regeneration. Shoot regeneration rate and shoot ratio (number of shoots/total number of explants) were used as phenotypic descriptors for regeneration capacity. Two independent experiments were carried out with six replicates of ten explants each. We found significant variation between the genotypes ranging from 0.88 to 88.33% for the regeneration rate and from 0.008 to 1.2 for the shoot ratio, which exceeded the rates reported so far. Furthermore, we found 88 SNP markers associated with either the shoot regeneration rate or the shoot ratio. In this association analysis, we found 12 SNP markers from ESTs (expressed sequence tags) matching known candidate genes that are involved in shoot morphogenesis. The best markers explained more than 51% of the variance in the shoot regeneration rate and more than 0.65 of the variance in the shoot regeneration ratio between the homozygote marker classes. The genes underlying some of the best markers such as a GT-transcription factor or an LRR receptor-like protein kinase are novel candidate genes putatively involved in the observed phenotypic differences. The associated markers were mapped to the closely related genome of Fragaria vesca and revealed many distinct clusters, which also comprised the known candidate genes that functioned in the organogenesis of plant shoots. However, the validation of candidate genes and their functional relationship to shoot regeneration require further analysis in independent rose populations and functional analyses.

Conclusive evidence for hexasomic inheritance in chrysanthemum based on analysis of a 183 k SNP array

BACKGROUND

Cultivated chrysanthemum is an outcrossing hexaploid (2n = 6× = 54) with a disputed mode of inheritance. In this paper, we present a single nucleotide polymorphism (SNP) selection pipeline that was used to design an Affymetrix Axiom array with 183 k SNPs from RNA sequencing data (1). With this array, we genotyped four bi-parental populations (with sizes of 405, 53, 76 and 37 offspring plants respectively), and a cultivar panel of 63 genotypes. Further, we present a method for dosage scoring in hexaploids from signal intensities of the array based on mixture models (2) and validation of selection steps in the SNP selection pipeline (3). The resulting genotypic data is used to draw conclusions on the mode of inheritance in chrysanthemum (4), and to make an inference on allelic expression bias (5).

RESULTS

With use of the mixture model approach, we successfully called the dosage of 73,936 out of 183,130 SNPs (40.4%) that segregated in any of the bi-parental populations. To investigate the mode of inheritance, we analysed markers that segregated in the large bi-parental population (n = 405). Analysis of segregation of duplex x nulliplex SNPs resulted in evidence for genome-wide hexasomic inheritance. This evidence was substantiated by the absence of strong linkage between markers in repulsion, which indicated absence of full disomic inheritance. We present the success rate of SNP discovery out of RNA sequencing data as affected by different selection steps, among which SNP coverage over genotypes and use of different types of sequence read mapping software. Genomic dosage highly correlated with relative allele coverage from the RNA sequencing data, indicating that most alleles are expressed according to their genomic dosage.

CONCLUSIONS

The large population, genotyped with a very large number of markers, is a unique framework for extensive genetic analyses in hexaploid chrysanthemum. As starting point, we show conclusive evidence for genome-wide hexasomic inheritance.

Partial preferential chromosome pairing is genotype dependent in tetraploid rose

It has long been recognised that polyploid species do not always neatly fall into the categories of auto- or allopolyploid, leading to the term 'segmental allopolyploid' to describe everything in between. The meiotic behaviour of such intermediate species is not fully understood, nor is there consensus as to how to model their inheritance patterns. In this study we used a tetraploid cut rose (Rosa hybrida) population, genotyped using the 68K WagRhSNP array, to construct an ultra-high-density linkage map of all homologous chromosomes using methods previously developed for autotetraploids. Using the predicted bivalent configurations in this population we quantified differences in pairing behaviour among and along homologous chromosomes, leading us to correct our estimates of recombination frequency to account for this behaviour. This resulted in the re-mapping of 25 695 SNP markers across all homologues of the seven rose chromosomes, tailored to the pairing behaviour of each chromosome in each parent. We confirmed the inferred differences in pairing behaviour among chromosomes by examining repulsion-phase linkage estimates, which also carry information about preferential pairing and recombination. Currently, the closest sequenced relative to rose is Fragaria vesca. Aligning the integrated ultra-dense rose map with the strawberry genome sequence provided a detailed picture of the synteny, confirming overall co-linearity but also revealing new genomic rearrangements. Our results suggest that pairing affinities may vary along chromosome arms, which broadens our current understanding of segmental allopolyploidy.

DNA-informed breeding of rosaceous crops: promises, progress and prospects

Crops of the Rosaceae family provide valuable contributions to rural economies and human health and enjoyment. Sustained solutions to production challenges and market demands can be met with genetically improved new cultivars. Traditional rosaceous crop breeding is expensive and time-consuming and would benefit from improvements in efficiency and accuracy. Use of DNA information is becoming conventional in rosaceous crop breeding, contributing to many decisions and operations, but only after past decades of solved challenges and generation of sufficient resources. Successes in deployment of DNA-based knowledge and tools have arisen when the 'chasm' between genomics discoveries and practical application is bridged systematically. Key steps are establishing breeder desire for use of DNA information, adapting tools to local breeding utility, identifying efficient application schemes, accessing effective services in DNA-based diagnostics and gaining experience in integrating DNA information into breeding operations and decisions. DNA-informed germplasm characterization for revealing identity and relatedness has benefitted many programs and provides a compelling entry point to reaping benefits of genomics research. DNA-informed germplasm evaluation for predicting trait performance has enabled effective reallocation of breeding resources when applied in pioneering programs. DNA-based diagnostics is now expanding from specific loci to genome-wide considerations. Realizing the full potential of this expansion will require improved accuracy of predictions, multi-trait DNA profiling capabilities, streamlined breeding information management systems, strategies that overcome plant-based features that limit breeding progress and widespread training of current and future breeding personnel and allied scientists.