Partial preferential chromosome pairing is genotype dependent in tetraploid rose

Rate of double reduction and genetic variability in yield, quality, and senescence related traits in tetraploid potato (Solanum tuberosum L.)

The amount of genetic variability is the foundation for genetic change in any plant breeding program, and the amount of double reduction can influence genetic gain and the amount of future genetic diversity in polyploid species. Our study investigates these factors using variance components analysis on a dataset comprising 13,131 potato breeding lines and phenotypic data from Scandinavian environments spanning 17 years (2003 to 2021). Pedigree information was used in quantitative genetic models to estimate additive genetic variance and the relative importance of additive and non-additive genetic variance. We used two models, a baseline model (M1) without effects due to specific combining ability (SCA) and M2 (including SCA due to interaction between parental genomes). Two cross-validation (CV) schemes [5-Fold and leave-one-breeding-cycle-out (LBCO)] were used to evaluate the prediction ability (PA) of each model. We estimated the rate of double reduction phenomenon (DRP) by determining the rate best fitting the data using a marginal likelihood approach. Our findings showed a wide range of variation in different traits, with very large proportion of additive genetic variance in dry matter content (DMC), but intermediate additive genetic variance for relative yield (RY), germination (GR), and withering (WNG). All traits showed modest non-additive genetic variance. Furthermore, genotype x environment interaction played a significant role in trait variability but is still much smaller than the additive genetic variance. After using different DRP rates, we found that a model with a 0.05 DRP rate provided the best fit to the data. Heritability estimates indicated a strong genetic basis for DMC, while other traits showed more moderate heritability, which shows contributions from both additive and interaction factors. Model comparison by 5-Fold CV and LBCO and the log likelihood ratio test (LRT) highlighted the importance of considering SCA when capturing trait variability. In 5-Fold CV, PA ranged from 0.296 to 0.812 in M1 and 0.300 to 0.813 in M2. Under LBCO CV, PA ranged from 0.180 to 0.726 in M1 and 0.180 to 0.728 in M2. However, an increase in PA in Model 2, which incorporates SCA, compared to Model 1, can be attributed to the inclusion of SCA effects. Furthermore, the LRT results indicated a highly significant difference between the models. CV and LRT suggest the need for genetic models that account for both additive and SCA effects. Our analysis also showed that genotype x environment interactions should be accounted for in order to maximize the accuracy of predicted breeding values of tetraploid potato clones. The rate of double reductions was small and insignificant.

Tests for segregation distortion in tetraploid F1 populations

KEY MESSAGE

In tetraploid F1 populations, traditional segregation distortion tests often inaccurately flag SNPs due to ignoring polyploid meiosis processes and genotype uncertainty. We develop tests that account for these factors. Genotype data from tetraploid F1 populations are often collected in breeding programs for mapping and genomic selection purposes. A common quality control procedure in these groups is to compare empirical genotype frequencies against those predicted by Mendelian segregation, where SNPs detected to have segregation distortion are discarded. However, current tests for segregation distortion are insufficient in that they do not account for double reduction and preferential pairing, two meiotic processes in polyploids that naturally change gamete frequencies, leading these tests to detect segregation distortion too often. Current tests also do not account for genotype uncertainty, again leading these tests to detect segregation distortion too often. Here, we incorporate double reduction, preferential pairing, and genotype uncertainty in likelihood ratio and Bayesian tests for segregation distortion. Our methods are implemented in a user-friendly R package, menbayes. We demonstrate the superiority of our methods to those currently used in the literature on both simulations and real data.

Haplotype-resolved genome assembly and resequencing provide insights into the origin and breeding of modern rose

Modern rose (Rosa hybrida) is a recently formed interspecific hybrid and has become one of the most important and widely cultivated ornamentals. Here we report the haplotype-resolved chromosome-scale genome assembly of the tetraploid R. hybrida 'Samantha' ('JACmantha') and a genome variation map of 233 Rosa accessions involving various wild species, and old and modern cultivars. Homologous chromosomes of 'Samantha' exhibit frequent homoeologous exchanges. Population genomic and genomic composition analyses reveal the contributions of wild Rosa species to modern roses and highlight that R. odorata and its derived cultivars are important contributors to modern roses, much like R. chinensis 'Old Blush'. Furthermore, selective sweeps during modern rose breeding associated with major agronomic traits, including continuous and recurrent flowering, double flower, flower senescence and disease resistance, are identified. This study provides insights into the genetic basis of modern rose origin and breeding history, and offers unprecedented genomic resources for rose improvement.

Genetic linkage mapping and quantitative trait locus (QTL) analysis of sweet basil (Ocimum basilicum L.) to identify genomic regions associated with cold tolerance and major volatiles

Chilling sensitivity is one of the greatest challenges affecting the marketability and profitability of sweet basil (Ocimum basilicum L.) in the US and worldwide. Currently, there are no sweet basils commercially available with significant chilling tolerance and traditional aroma profiles. This study was conducted to identify quantitative trait loci (QTLs) responsible for chilling tolerance and aroma compounds in a biparental mapping population, including the Rutgers advanced breeding line that served as a chilling tolerant parent, 'CB15', the chilling sensitive parent, 'Rutgers Obsession DMR' and 200 F2 individuals. Chilling tolerance was assessed by percent necrosis using machine learning and aroma profiling was evaluated using gas chromatography (GC) mass spectrometry (MS). Single nucleotide polymorphism (SNP) markers were generated from genomic sequences derived from double digestion restriction-site associated DNA sequencing (ddRADseq) and converted to genotype data using a reference genome alignment. A genetic linkage map was constructed and five statistically significant QTLs were identified in response to chilling temperatures with possible interactions between QTLs. The QTL on LG24 (qCH24) demonstrated the largest effect for chilling response and was significant in all three replicates. No QTLs were identified for linalool, as the population did not segregate sufficiently to detect this trait. Two significant QTLs were identified for estragole (also known as methyl chavicol) with only qEST1 on LG1 being significant in the multiple-QTL model (MQM). QEUC26 was identified as a significant QTL for eucalyptol (also known as 1,8-cineole) on LG26. These QTLs may represent key mechanisms for chilling tolerance and aroma in basil, providing critical knowledge for future investigation of these phenotypic traits and molecular breeding.

Inheritance pattern of tetraploids pummelo, mandarin, and their interspecific hybrid sour orange is highly influenced by their phylogenomic structure

Citrus polyploidy is associated with a wide range of morphological, genetic, and physiological changes that are often advantageous for breeding. Citrus triploid hybrids are very interesting as new seedless varieties. However, tetraploid rootstocks promote adaptation to different abiotic stresses and promote resilience. Triploid and tetraploid hybrids can be obtained through sexual hybridizations using tetraploid parents (2x × 4x, 4x × 2x, or 4x × 4x), but more knowledge is needed about the inheritance pattern of tetraploid parents to optimize the efficiency of triploid varieties and tetraploid rootstock breeding strategies. In this work, we have analyzed the inheritance pattern of three tetraploid genotypes: 'Chandler' pummelo (Citrus maxima) and 'Cleopatra' mandarin (Citrus reticulata), which represent two clear examples of autotetraploid plants constituted by the genome of a single species, and the 'Sevillano' sour orange, which is an allotetraploid interspecific hybrid between C. maxima and C. reticulata. Polymorphic simple sequence repeat (SSR) and single-nucleotide polymorphism (SNP) markers were used to estimate parental heterozygosity restitution, and allele frequencies for centromeric loci were used to calculate the preferential pairing rate related to the proportion of disomic and tetrasomic segregation. The tetraploid pummelo and mandarin displayed tetrasomic segregation. Sour orange evidenced a clear intermediate inheritance for five of the nine chromosomes (1, 2, 5, 7, and 8), a slight tendency toward tetrasomic inheritance on chromosome 3, and intermediate inheritance with a tendency toward disomy for chromosomes 4, 6, and 9. These results indicate that the interspecific versus intraspecific phylogenomic origin affects preferential pairing and, therefore, the inheritance patterns. Despite its high level of heterozygosity, the important preferential chromosome pairing observed in sour orange results in a limited diversity of the genotypic variability of its diploid gametes, and consequently, a large part of the genetic value of the original diploid sour orange is transferred to the tetraploid progenies.

QTL discovery for resistance to black spot and cercospora leaf spot, and defoliation in two interconnected F1 bi-parental tetraploid garden rose populations

Garden roses are an economically important horticultural crop worldwide, and two major fungal pathogens, black spot (Diplocarpon rosae F.A. Wolf) and cercospora leaf spot of rose (Rosisphaerella rosicola Pass.), affect both the health and ornamental value of the plant. Most studies on black spot disease resistance have focused on diploid germplasm, and little work has been performed on cercospora leaf spot resistance. With the use of newly developed software tools for autopolyploid genetics, two interconnected tetraploid garden rose F1 populations (phenotyped over the course of 3 years) were used for quantitative trait locus (QTL) analysis of black spot and cercospora leaf spot resistance as well as plant defoliation. QTLs for black spot resistance were mapped to linkage groups (LGs) 1-6. QTLs for cercospora resistance and susceptibility were found in LGs 1, 4, and 5 and for defoliation in LGs 1, 3, and 5. The major locus on LG 5 for black spot resistance coincides with the previously discovered Rdr4 locus inherited from Rosa L. 'Radbrite' (Brite Eyes™), the common parent used in these mapping populations. This work is the first report of any QTL for cercospora resistance/susceptibility in tetraploid rose germplasm and the first report of defoliation QTL in roses. A major QTL for cercospora susceptibility coincides with the black spot resistance QTL on LG 5 (Rdr4). A major cercospora resistance QTL was found on LG 1. These populations provide a genetic resource that will further the knowledge base of rose genetics as more traits are studied. Studying more traits from these populations will allow for the stacking of various QTLs for desirable traits.

Autopolyploid inheritance and a heterozygous reciprocal translocation shape chromosome genetic behavior in tetraploid blueberry (Vaccinium corymbosum)

Understanding chromosome recombination behavior in polyploidy species is key to advancing genetic discoveries. In blueberry, a tetraploid species, the line of evidences about its genetic behavior still remain poorly understood, owing to the inter-specific, and inter-ploidy admixture of its genome and lack of in depth genome-wide inheritance and comparative structural studies. Here we describe a new high-quality, phased, chromosome-scale genome of a diploid blueberry, clone W85. The genome was integrated with cytogenetics and high-density, genetic maps representing six tetraploid blueberry cultivars, harboring different levels of wild genome admixture, to uncover recombination behavior and structural genome divergence across tetraploid and wild diploid species. Analysis of chromosome inheritance and pairing demonstrated that tetraploid blueberry behaves as an autotetraploid with tetrasomic inheritance. Comparative analysis demonstrated the presence of a reciprocal, heterozygous, translocation spanning one homolog of chr-6 and one of chr-10 in the cultivar Draper. The translocation affects pairing and recombination of chromosomes 6 and 10. Besides the translocation detected in Draper, no other structural genomic divergences were detected across tetraploid cultivars and highly inter-crossable wild diploid species. These findings and resources will facilitate new genetic and comparative genomic studies in Vaccinium and the development of genomic assisted selection strategy for this crop.

Rose Rosette Disease Resistance Loci Detected in Two Interconnected Tetraploid Garden Rose Populations

Rose rosette disease (RRD), caused by the Rose rosette emaravirus (RRV), is a major threat to the garden rose industry in the United States. There has been limited work on the genetics of host plant resistance to RRV. Two interconnected tetraploid garden rose F1 biparental mapping populations were created to develop high-quality tetraploid rose linkage maps that allowed the discovery of RRD resistance quantitative trait loci (QTLs) on linkage groups (LGs) 5, 6, and 7. These QTLs individually accounted for around 18-40% of the phenotypic variance. The locus with the greatest effect on partial resistance was found in LG 5. Most individuals with the LG 5 QTL were in the simplex configuration; however, two individuals were duplex (likely due to double reduction). Identification of resistant individuals and regions of interest can help the development of diagnostic markers for marker-assisted selection in a breeding program.

Genome architecture and tetrasomic inheritance of autotetraploid potato

Potato (Solanum tuberosum) is the most consumed non-cereal food crop. Most commercial potato cultivars are autotetraploids with highly heterozygous genomes, severely hampering genetic analyses and improvement. By leveraging the state-of-the-art sequencing technologies and polyploid graph binning, we achieved a chromosome-scale, haplotype-resolved genome assembly of a cultivated potato, Cooperation-88 (C88). Intra-haplotype comparative analyses revealed extensive sequence and expression differences in this tetraploid genome. We identified haplotype-specific pericentromeres on chromosomes, suggesting a distinct evolutionary trajectory of potato homologous centromeres. Furthermore, we detected double reduction events that are unevenly distributed on haplotypes in 1021 of 1034 selfing progeny, a feature of autopolyploid inheritance. By distinguishing maternal and paternal haplotype sets in C88, we simulated the origin of heterosis in cultivated tetraploid with a survey of 3110 tetra-allelic loci with deleterious mutations, which were masked in the heterozygous condition by two parents. This study provides insights into the genomic architecture of autopolyploids and will guide their breeding.

Identification of QTLs for Reduced Susceptibility to Rose Rosette Disease in Diploid Roses

Resistance to rose rosette disease (RRD), a fatal disease of roses (Rosa spp.), is a high priority for rose breeding. As RRD resistance is time-consuming to phenotype, the identification of genetic markers for resistance could expedite breeding efforts. However, little is known about the genetics of RRD resistance. Therefore, we performed a quantitative trait locus (QTL) analysis on a set of inter-related diploid rose populations phenotyped for RRD resistance and identified four QTLs. Two QTLs were found in multiple years. The most consistent QTL is qRRV_TX2WSE_ch5, which explains approximately 20% and 40% of the phenotypic variation in virus quantity and severity of RRD symptoms, respectively. The second, a QTL on chromosome 1, qRRD_TX2WSE_ch1, accounts for approximately 16% of the phenotypic variation for severity. Finally, a third QTL on chromosome 3 was identified only in the multiyear analysis, and a fourth on chromosome 6 was identified in data from one year only. In addition, haplotypes associated with significant changes in virus quantity and severity were identified for qRRV_TX2WSE_ch5 and qRRD_TX2WSE_ch1. This research represents the first report of genetic determinants of resistance to RRD. In addition, marker trait associations discovered here will enable better parental selection when breeding for RRD resistance and pave the way for marker-assisted selection for RRD resistance.

Molecular Authentication, Phytochemical Evaluation and Asexual Propagation of Wild-Growing Rosa canina L. (Rosaceae) Genotypes of Northern Greece for Sustainable Exploitation

Dogroses belong to a taxonomically difficult genus and family and represent important phytogenetic resources associated with high ornamental, pharmaceutical-cosmetic and nutritional values, thus suggesting a potentially high exploitation merit. Triggered by these prospects, wild-growing Rosa canina populations of Greece were selected for investigation and evaluation of their potential for integrated domestication. We collected ripe rosehips from Greek native wild-growing populations (samples from seven genotypes) for phytochemical analysis (total phenolics, total flavonoids, antioxidant activity and vitamin C content), leaf samples for DNA analysis using the ITS2 sequence (nine genotypes) and fresh soft-wood stem cuttings for propagation trials (seven genotypes). After evaluation of these materials, this study reports for the first-time distinct DNA-fingerprinted genotypes from Greece with interesting phytochemical profiles mainly in terms of Vitamic C content (up to 500.22 ± 0.15 mg of ascorbic acid equivalents/100 g of sample) as well as effective asexual propagation protocols for prioritized R. canina genotypes via cuttings. The latter highlights the importance of the levels of external hormone application (2000 ppm of indole-3-butyric acid), the effect of season (highly-effective spring trials) and genotype-specific differences in rooting capacities of the studied genotypes. All inclusive, this study offers new artificially selected material of Greek native R. canina with a consolidated identity and interesting phytochemical profile. These materials are currently under ex-situ conservation for further evaluation and characterization in pilot field studies, thus facilitating its sustainable exploitation for applications in the agro-alimentary, medicinal-cosmetic, and ornamental sectors.

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.

Detecting quantitative trait loci and exploring chromosomal pairing in autopolyploids using polyqtlR

Motivation

The investigation of quantitative trait loci (QTL) is an essential component in our understanding of how organisms vary phenotypically. However, many important crop species are polyploid (carrying more than two copies of each chromosome), requiring specialized tools for such analyses. Moreover, deciphering meiotic processes at higher ploidy levels is not straightforward, but is necessary to understand the reproductive dynamics of these species, or uncover potential barriers to their genetic improvement.

Results

Here, we present polyqtlR, a novel software tool to facilitate such analyses in (auto)polyploid crops. It performs QTL interval mapping in F₁ populations of outcrossing polyploids of any ploidy level using identity-by-descent probabilities. The allelic composition of discovered QTL can be explored, enabling favourable alleles to be identified and tracked in the population. Visualization tools within the package facilitate this process, and options to include genetic co-factors and experimental factors are included. Detailed information on polyploid meiosis including prediction of multivalent pairing structures, detection of preferential chromosomal pairing and location of double reduction events can be performed.

Availabilityand implementation

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

Supplementary information

Supplementary data are available at Bioinformatics online.

Molecular characteristics of S-RNase alleles as the determinant of self-incompatibility in the style of Fragaria viridis

Strawberry (Fragaria spp.) is a member of the Rosoideae subfamily in the family Rosaceae. The self-incompatibility (SI) of some diploid species is a key agronomic trait that acts as a basic pollination barrier; however, the genetic mechanism underlying SI control in strawberry remains unclear. Two candidate S-RNases (S_a- and S_b-RNase) identified in the transcriptome of the styles of the self-incompatible Fragaria viridis 42 were confirmed to be SI determinants at the S locus following genotype identification and intraspecific hybridization using selfing progenies. Whole-genome collinearity and RNase T2 family analysis revealed that only an S locus exists in Fragaria; however, none of the compatible species contained S-RNase. Although the results of interspecific hybridization experiments showed that F. viridis (SI) styles could accept pollen from F. mandshurica (self-compatible), the reciprocal cross was incompatible. S_a and S_b-RNase contain large introns, and their noncoding sequences (promotors and introns) can be transcribed into long noncoding RNAs (lncRNAs). Overall, the genus Fragaria exhibits S-RNase-based gametophytic SI, and S-RNase loss occurs at the S locus of compatible germplasms. In addition, a type of SI-independent unilateral incompatibility exists between compatible and incompatible Fragaria species. Furthermore, the large introns and neighboring lncRNAs in S-RNase in Fragaria could offer clues about S-RNase expression strategies.

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.

Genome editing of polyploid crops: prospects, achievements and bottlenecks

Plant breeding aims to develop improved crop varieties. Many crops have a polyploid and often highly heterozygous genome, which may make breeding of polyploid crops a real challenge. The efficiency of traditional breeding based on crossing and selection has been improved by using marker-assisted selection (MAS), and MAS is also being applied in polyploid crops, which helps e.g. for introgression breeding. However, methods such as random mutation breeding are difficult to apply in polyploid crops because there are multiple homoeologous copies (alleles) of each gene. Genome editing technology has revolutionized mutagenesis as it enables precisely selecting targets. The genome editing tool CRISPR/Cas is especially valuable for targeted mutagenesis in polyploids, as all alleles and/or copies of a gene can be targeted at once. Even multiple genes, each with multiple alleles, may be targeted simultaneously. In addition to targeted mutagenesis, targeted replacement of undesirable alleles by desired ones may become a promising application of genome editing for the improvement of polyploid crops, in the near future. Several examples of the application of genome editing for targeted mutagenesis are described here for a range of polyploid crops, and achievements and bottlenecks are highlighted.

A tetrasomic inheritance model and likelihood-based method for mapping quantitative trait loci in autotetraploid species

Dissecting the genetic architecture of quantitative traits in autotetraploid species is a methodologically challenging task, but a pivotally important goal for breeding globally important food crops, including potato and blueberry, and ornamental species such as rose. Mapping quantitative trait loci (QTLs) is now a routine practice in diploid species but is far less advanced in autotetraploids, largely due to a lack of analytical methods that account for the complexities of tetrasomic inheritance. We present a novel likelihood-based method for QTL mapping in outbred segregating populations of autotetraploid species. The method accounts properly for sophisticated features of gene segregation and recombination in an autotetraploid meiosis. It may model and analyse molecular marker data with or without allele dosage information, such as that from microarray or sequencing experiments. The method developed outperforms existing bivalent-based methods, which may fail to model and analyse the full spectrum of experimental data, in the statistical power of QTL detection, and accuracy of QTL location, as demonstrated by an intensive simulation study and analysis of data sets collected from a segregating population of potato (Solanum tuberosum). The study enables QTL mapping analysis to be conducted in autotetraploid species under a rigorous tetrasomic inheritance model.

QTL Mapping for Resistance to Cankers Induced by Pseudomonas syringae pv. actinidiae (Psa) in a Tetraploid Actinidia chinensis Kiwifruit Population

Polyploidy is a key driver of significant evolutionary changes in plant species. The genus Actinidia (kiwifruit) exhibits multiple ploidy levels, which contribute to novel fruit traits, high yields and resistance to the canker-causing dieback disease incited by Pseudomonas syringae pv. actinidiae (Psa) biovar 3. However, the genetic mechanism for resistance to Psa observed in polyploid kiwifruit is not yet known. In this study we performed detailed genetic analysis of a tetraploid Actinidia chinensis var. chinensis population derived from a cross between a female parent that exhibits weak tolerance to Psa and a highly Psa-resistant male parent. We used the capture-sequencing approach across the whole kiwifruit genome and generated the first ultra-dense maps in a tetraploid kiwifruit population. We located quantitative trait loci (QTLs) for Psa resistance on these maps. Our approach to QTL mapping is based on the use of identity-by-descent trait mapping, which allowed us to relate the contribution of specific alleles from their respective homologues in the male and female parent, to the control of Psa resistance in the progeny. We identified genes in the diploid reference genome whose function is suggested to be involved in plant defense, which underly the QTLs, including receptor-like kinases. Our study is the first to cast light on the genetics of a polyploid kiwifruit and suggest a plausible mechanism for Psa resistance in this species.

RosBREED: bridging the chasm between discovery and application to enable DNA-informed breeding in rosaceous crops

The Rosaceae crop family (including almond, apple, apricot, blackberry, peach, pear, plum, raspberry, rose, strawberry, sweet cherry, and sour cherry) provides vital contributions to human well-being and is economically significant across the U.S. In 2003, industry stakeholder initiatives prioritized the utilization of genomics, genetics, and breeding to develop new cultivars exhibiting both disease resistance and superior horticultural quality. However, rosaceous crop breeders lacked certain knowledge and tools to fully implement DNA-informed breeding-a "chasm" existed between existing genomics and genetic information and the application of this knowledge in breeding. The RosBREED project ("Ros" signifying a Rosaceae genomics, genetics, and breeding community initiative, and "BREED", indicating the core focus on breeding programs), addressed this challenge through a comprehensive and coordinated 10-year effort funded by the USDA-NIFA Specialty Crop Research Initiative. RosBREED was designed to enable the routine application of modern genomics and genetics technologies in U.S. rosaceous crop breeding programs, thereby enhancing their efficiency and effectiveness in delivering cultivars with producer-required disease resistances and market-essential horticultural quality. This review presents a synopsis of the approach, deliverables, and impacts of RosBREED, highlighting synergistic global collaborations and future needs. Enabling technologies and tools developed are described, including genome-wide scanning platforms and DNA diagnostic tests. Examples of DNA-informed breeding use by project participants are presented for all breeding stages, including pre-breeding for disease resistance, parental and seedling selection, and elite selection advancement. The chasm is now bridged, accelerating rosaceous crop genetic improvement.

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.

Preferential Disomic Segregation and C. micrantha/C. medica Interspecific Recombination in Tetraploid 'Giant Key' Lime; Outlook for Triploid Lime Breeding

The triploid 'Tahiti' lime (C. x latifolia (Yu. Tanaka) Tanaka) naturally originated from a merger between a haploid ovule of lemon (C. x limon (L.) Burm) and a diploid pollen from a 'Mexican' lime (C. x aurantiifolia (Christm.) Swing). The very limited natural inter-varietal diversity and gametic sterility of C. latifolia requires a phylogenomic based reconstruction breeding strategy to insure its diversification. We developed a strategy based on interploid hybridization between diploid lemon and the doubled diploid 'Giant Key' lime. This lime is a doubled diploid of 'Mexican' lime, itself a natural interspecific F1 hybrid between C. medica L. and C. micrantha Wester. For an optimized breeding program, we analyzed the meiotic behavior of the allotetraploid lime, the genetic structure of its diploid gametes, the interspecific recombination between C. medica and C. micrantha, and constructed its genetic map. A population of 272 triploid hybrids was generated using 'Giant Key' lime as pollinator. One hundred fifty-eight SNPs diagnostic of C. micrantha, regularly distributed throughout the citrus genome were successfully developed and applied. The genetic structure of the diploid gametes was examined based on C. micrantha doses along the genome. The diploid gametes transmitted in average 91.17% of the parental interspecific C. medica/C. micrantha heterozygosity. Three chromosomes (2, 8, and 9) showed disomic segregation with high preferential pairing values, while the remaining chromosomes showed an intermediate inheritance with a preferential disomic trend. A total of 131 SNPs were assigned to nine linkage groups to construct the genetic map. It spanned 272.8 cM with a low average recombination rate (0.99 cM Mb-1) and high synteny and colinearity with the reference clementine genome. Our results confirmed that an efficient reconstruction breeding strategy for 'Tahiti' lime is possible, based on interploid hybridization using a doubled diploid of C. aurantiifolia. The tetraploid parent should be selected for favorable agronomic traits and its genetic value should be efficiently inherited by the progeny thanks to transmission of the high level of parental heterozygosity. However, it would require developing numerous progeny to overcome the linkage drag caused by the limited interspecific recombination associated with the predominant disomic inheritance.

Varietal variation and chromosome behaviour during meiosis in Solanum tuberosum

Naturally occurring autopolyploid species, such as the autotetraploid potato Solanum tuberosum, face a variety of challenges during meiosis. These include proper pairing, recombination and correct segregation of multiple homologous chromosomes, which can form complex multivalent configurations at metaphase I, and in turn alter allelic segregation ratios through double reduction. Here, we present a reference map of meiotic stages in diploid and tetraploid S. tuberosum using fluorescence in situ hybridisation (FISH) to differentiate individual meiotic chromosomes 1 and 2. A diploid-like behaviour at metaphase I involving bivalent configurations was predominant in all three tetraploid varieties. The crossover frequency per bivalent was significantly reduced in the tetraploids compared with a diploid variety, which likely indicates meiotic adaptation to the autotetraploid state. Nevertheless, bivalents were accompanied by a substantial frequency of multivalents, which varied by variety and by chromosome (7-48%). We identified possible sites of synaptic partner switching, leading to multivalent formation, and found potential defects in the polymerisation and/or maintenance of the synaptonemal complex in tetraploids. These findings demonstrate the rise of S. tuberosum as a model for autotetraploid meiotic recombination research and highlight constraints on meiotic chromosome configurations and chiasma frequencies as an important feature of an evolved autotetraploid meiosis.

Mapping the black spot resistance locus Rdr3 in the shrub rose 'George Vancouver' allows for the development of improved diagnostic markers for DNA-informed breeding

Rdr3 is a novel resistance gene of black spot in roses that maps to a chromosome 6 homolog. A new DNA test was developed and can be used to pyramid black spot resistance in roses. Diplocarpon rosae, the cause of rose black spot, is one of the most devastating foliar pathogens of cultivated roses (Rosa spp.). The primary method of disease control is fungicides, and they are viewed unfavorably by home gardeners due to potential environmental and health impacts. Planting rose cultivars with genetic resistance to black spot can reduce many of the fungicide applications needed in an integrated pest management system. To date, four resistance genes have been identified in roses (Rdr1, Rdr2, Rdr3, and Rdr4). Rdr3 was never mapped and is thought to be unique from Rdr1 and Rdr2. It is unknown whether it is an allele of Rdr4. To assess the novelty of Rdr3, a mapping population was created by crossing the Rdr3 containing cultivar George Vancouver with the susceptible cultivar Morden Blush. The mapping population was genotyped with the WagRhSNP 68 K Axiom array and mapped using the 'polymapR' package. Rdr3 was mapped to a chromosome 6 homolog confirming it is different from Rdr1 and Rdr2, found on chromosome 1, and from Rdr4, found on chromosome 5. The mapping information was used in conjunction with the Rosa chinensis genome assembly to develop new tightly linked SSRs for marker-assisted breeding. Three markers were able to predict the presence of Rdr3 in a 63-cultivar validation set. Additionally, 12 cultivars appear to have resistance genes other than Rdr3. The improved diagnostic markers will be a great asset to the rose-breeding community toward developing new black spot-resistant cultivars.

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.

Probabilistic reconstruction of genealogies for polyploid plant species

A probabilistic reconstruction of genealogies in a polyploid population (from 2x to 4x) is investigated, by considering genetic data analyzed as the probability of allele presence in a given genotype. Based on the likelihood of all possible crossbreeding patterns, our model enables us to infer and to quantify the whole potential genealogies in the population. We explain in particular how to deal with the uncertain allelic multiplicity that may occur with polyploids. Then we build an ad hoc penalized likelihood to compare genealogies and to decide whether a particular individual brings sufficient information to be included in the taken genealogy. This decision criterion enables us in a next part to suggest a greedy algorithm in order to explore missing links and to rebuild some connections in the genealogies, retrospectively. As a by-product, we also give a way to infer the individuals that may have been favored by breeders over the years. In the last part we highlight the results given by our model and our algorithm, firstly on a simulated population and then on a real population of rose bushes. Most of the methodology relies on the maximum likelihood principle and on graph theory.

pSBVB: A Versatile Simulation Tool To Evaluate Genomic Selection in Polyploid Species

Genomic Selection (GS) is the procedure whereby molecular information is used to predict complex phenotypes and it is standard in many animal and plant breeding schemes. However, only a small number of studies have been reported in horticultural crops, and in polyploid species in particular. In this paper, we have developed a versatile forward simulation tool, called polyploid Sequence Based Virtual Breeding (pSBVB), to evaluate GS strategies in polyploids; pSBVB is an efficient gene dropping software that can simulate any number of complex phenotypes, allowing a very flexible modeling of phenotypes suited to polyploids. As input, it takes genotype data from the founder population, which can vary from single nucleotide polymorphisms (SNP) chips up to sequence, a list of causal variants for every trait and their heritabilities, and the pedigree. Recombination rates between homeologous chromosomes can be specified, so that both allo- and autopolyploid species can be considered. The program outputs phenotype and genotype data for all individuals in the pedigree. Optionally, it can produce several genomic relationship matrices that consider exact or approximate genotype values. pSBVB can therefore be used to evaluate GS strategies in polyploid species (say varying SNP density, genetic architecture or population size, among other factors), or to optimize experimental designs for association studies. We illustrate pSBVB with SNP data from tetraploid potato and partial sequence data from octoploid strawberry, and we show that GS is a promising breeding strategy for polyploid species but that the actual advantage critically depends on the underlying genetic architecture. Source code, examples and a complete manual are freely available in GitHub https://github.com/lauzingaretti/pSBVB.

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.

Genotyping Polyploids from Messy Sequencing Data

Detecting and quantifying the differences in individual genomes (i.e., genotyping), plays a fundamental role in most modern bioinformatics pipelines. Many scientists now use reduced representation next-generation sequencing (NGS) approaches for genotyping. Genotyping diploid individuals using NGS is a well-studied field, and similar methods for polyploid individuals are just emerging. However, there are many aspects of NGS data, particularly in polyploids, that remain unexplored by most methods. Our contributions in this paper are fourfold: (i) We draw attention to, and then model, common aspects of NGS data: sequencing error, allelic bias, overdispersion, and outlying observations. (ii) Many datasets feature related individuals, and so we use the structure of Mendelian segregation to build an empirical Bayes approach for genotyping polyploid individuals. (iii) We develop novel models to account for preferential pairing of chromosomes, and harness these for genotyping. (iv) We derive oracle genotyping error rates that may be used for read depth suggestions. We assess the accuracy of our method in simulations, and apply it to a dataset of hexaploid sweet potato (Ipomoea batatas). An R package implementing our method is available at https://cran.r-project.org/package=updog.

Multi-environment QTL analysis of plant and flower morphological traits in tetraploid rose

KEY MESSAGE

Rose morphological traits such as prickles or petal number are influenced by a few key QTL which were detected across different growing environments-necessary for genomics-assisted selection in non-target environments. Rose, one of the world's most-loved and commercially important ornamental plants, is predominantly tetraploid, possessing four rather than two copies of each chromosome. This condition complicates genetic analysis, and so the majority of previous genetic studies in rose have been performed at the diploid level. However, there may be advantages to performing genetic analyses at the tetraploid level, not least because this is the ploidy level of most breeding germplasm. Here, we apply recently developed methods for quantitative trait loci (QTL) detection in a segregating tetraploid rose population (F1 = 151) to unravel the genetic control of a number of key morphological traits. These traits were measured both in the Netherlands and Kenya. Since ornamental plant breeding and selection are increasingly being performed at locations other than the production sites, environment-neutral QTL are required to maximise the effectiveness of breeding programmes. We detected a number of robust, multi-environment QTL for such traits as stem and petiole prickles, petal number and stem length that were localised on the recently developed high-density SNP linkage map for rose. Our work explores the complex genetic architecture of these important morphological traits at the tetraploid level, while helping to advance the methods for marker-trait exploration in polyploid species.

Genome engineering in ornamental plants: Current status and future prospects

Ornamental plants, like roses, carnations, and chrysanthemums, are economically important and are sold all over the world. In addition, numerous cut and garden flowers add colors to homes and gardens. Various strategies of plant breeding have been employed to improve traits of many ornamental plants. These approaches span from conventional techniques, such as crossbreeding and mutation breeding, to genetically modified plants. Recently, genome editing has become available as an efficient means for modifying traits in plant species. Genome editing technology is useful for genetic analysis and is poised to become a common breeding method for ornamental plants. In this review, we summarize the benefits and limitations of conventional breeding techniques and genome editing methods and discuss their future potential to accelerate the rate breeding programs in ornamental plants.

A miR172 target-deficient AP2-like gene correlates with the double flower phenotype in roses

One of the well-known floral abnormalities in flowering plants is the double-flower phenotype, which corresponds to flowers that develop extra petals, sometimes even containing entire flowers within flowers. Because of their highly priced ornamental value, spontaneous double-flower variants have been found and selected for in a wide range of ornamental species. Previously, double flower formation in roses was associated with a restriction of AGAMOUS expression domain toward the centre of the meristem, leading to extra petals. Here, we characterized the genomic region containing the mutation associated with the switch from simple to double flowers in the rose. An APETALA2-like gene (RcAP2L), a member of the Target Of EAT-type (TOE-type) subfamily, lies within this interval. In the double flower rose, two alleles of RcAP2L are present, one of which harbours a transposable element inserted into intron 8. This insertion leads to the creation of a miR172 resistant RcAP2L variant. Analyses of the presence of this variant in a set of simple and double flower roses demonstrate a correlation between the presence of this allele and the double flower phenotype. These data suggest a role of this miR172 resistant RcAP2L variant in regulating RcAGAMOUS expression and double flower formation in Rosa sp.

A high-quality genome sequence of Rosa chinensis to elucidate ornamental traits

Rose is the world's most important ornamental plant, with economic, cultural and symbolic value. Roses are cultivated worldwide and sold as garden roses, cut flowers and potted plants. Roses are outbred and can have various ploidy levels. Our objectives were to develop a high-quality reference genome sequence for the genus Rosa by sequencing a doubled haploid, combining long and short reads, and anchoring to a high-density genetic map, and to study the genome structure and genetic basis of major ornamental traits. We produced a doubled haploid rose line ('HapOB') from Rosa chinensis 'Old Blush' and generated a rose genome assembly anchored to seven pseudo-chromosomes (512 Mb with N50 of 3.4 Mb and 564 contigs). The length of 512 Mb represents 90.1-96.1% of the estimated haploid genome size of rose. Of the assembly, 95% is contained in only 196 contigs. The anchoring was validated using high-density diploid and tetraploid genetic maps. We delineated hallmark chromosomal features, including the pericentromeric regions, through annotation of transposable element families and positioned centromeric repeats using fluorescent in situ hybridization. The rose genome displays extensive synteny with the Fragaria vesca genome, and we delineated only two major rearrangements. Genetic diversity was analysed using resequencing data of seven diploid and one tetraploid Rosa species selected from various sections of the genus. Combining genetic and genomic approaches, we identified potential genetic regulators of key ornamental traits, including prickle density and the number of flower petals. A rose APETALA2/TOE homologue is proposed to be the major regulator of petal number in rose. This reference sequence is an important resource for studying polyploidization, meiosis and developmental processes, as we demonstrated for flower and prickle development. It will also accelerate breeding through the development of molecular markers linked to traits, the identification of the genes underlying them and the exploitation of synteny across Rosaceae.

A high-quality genome sequence of Rosa chinensis to elucidate ornamental traits

Rose is the world's most important ornamental plant, with economic, cultural and symbolic value. Roses are cultivated worldwide and sold as garden roses, cut flowers and potted plants. Roses are outbred and can have various ploidy levels. Our objectives were to develop a high-quality reference genome sequence for the genus Rosa by sequencing a doubled haploid, combining long and short reads, and anchoring to a high-density genetic map, and to study the genome structure and genetic basis of major ornamental traits. We produced a doubled haploid rose line ('HapOB') from Rosa chinensis 'Old Blush' and generated a rose genome assembly anchored to seven pseudo-chromosomes (512 Mb with N50 of 3.4 Mb and 564 contigs). The length of 512 Mb represents 90.1-96.1% of the estimated haploid genome size of rose. Of the assembly, 95% is contained in only 196 contigs. The anchoring was validated using high-density diploid and tetraploid genetic maps. We delineated hallmark chromosomal features, including the pericentromeric regions, through annotation of transposable element families and positioned centromeric repeats using fluorescent in situ hybridization. The rose genome displays extensive synteny with the Fragaria vesca genome, and we delineated only two major rearrangements. Genetic diversity was analysed using resequencing data of seven diploid and one tetraploid Rosa species selected from various sections of the genus. Combining genetic and genomic approaches, we identified potential genetic regulators of key ornamental traits, including prickle density and the number of flower petals. A rose APETALA2/TOE homologue is proposed to be the major regulator of petal number in rose. This reference sequence is an important resource for studying polyploidization, meiosis and developmental processes, as we demonstrated for flower and prickle development. It will also accelerate breeding through the development of molecular markers linked to traits, the identification of the genes underlying them and the exploitation of synteny across Rosaceae.

Genomic and transcriptomic sequencing of Rosa hybrida provides microsatellite markers for breeding, flower trait improvement and taxonomy studies

BACKGROUND

Rosa hybrida is a valuable ornamental, food and medicinal crop worldwide, but with relatively limited molecular marker resources, especially for flower-specific markers. In this study, we performed genomic and floral transcriptomic sequencing of modern rose. We obtained comprehensive nucleotide information, from which numerous potential simple sequence repeat (SSR) markers were identified but were found to have high rates of amplification failure and PCR product redundancy.

RESULTS

We applied a filtering strategy for BLAST analysis with the assembled genomic sequence and identified 124,591 genomic and 2,292 EST markers with unique annealing sites. These markers had much greater reliability than those obtained before filtering. Additional BLAST analysis against the transcriptomic sequences uncovered 5225 genomic SSRs associated with 4100 transcripts, 2138 of which were associated with functional genes that were annotated against the non-redundant database. More than 90% of these newly developed molecular markers were polymorphic, based on PCR using a subset of SSRs to analyze tetraploid modern rose accessions, diploid Rosa species and one strawberry accession. The relationships among Rosa species determined by cluster analysis (based on these results) were in agreement with modern rose breeding history, whereas strawberry was isolated in a separate cluster, as expected.

CONCLUSIONS

Our results provide valuable molecular-genetic tools for rose flower trait improvement, breeding and taxonomy. Importantly, we describe a reproducible organ-specific strategy for molecular marker development and selection in plants, which can be applied to other crops.

Segmental allopolyploidy in action: Increasing diversity through polyploid hybridization and homoeologous recombination

PREMISE OF THE STUDY

The genetic bottleneck of polyploid formation can be mitigated by multiple origins, gene flow, and recombination among different lineages. In crop plants with limited origins, efforts to increase genetic diversity have limitations. Here we used lineage recombination to increase genetic diversity in peanut, an allotetraploid likely of single origin, by crossing with a novel allopolyploid genotype and selecting improved lines.

METHODS

Single backcross progeny from cultivated peanut × wild species-derived allotetraploid cross were studied over successive generations. Using genetic assumptions that encompass segmental allotetraploidy, we used single nucleotide polymorphisms and whole-genome sequence data to infer genome structures.

KEY RESULTS

Selected lines, despite a high proportion of wild alleles, are agronomically adapted, productive, and with improved disease resistances. Wild alleles mostly substituted homologous segments of the peanut genome. Regions of dispersed wild alleles, characteristic of gene conversion, also occurred. However, wild chromosome segments sometimes replaced cultivated peanut's homeologous subgenome; A. ipaënsis B sometimes replaced A. hypogaea A subgenome (~0.6%), and A. duranensis replaced A. hypogaea B subgenome segments (~2%). Furthermore, some subgenome regions historically lost in cultivated peanut were "recovered" by wild chromosome segments (effectively reversing the "polyploid ratchet"). These processes resulted in lines with new genome structure variations.

CONCLUSIONS

Genetic diversity was introduced by wild allele introgression, and by introducing new genome structure variations. These results highlight the special possibilities of segmental allotetraploidy and of using lineage recombination to increase genetic diversity in peanut, likely mirroring what occurs in natural segmental allopolyploids with multiple origins.

The Rosa genome provides new insights into the domestication of modern roses

Roses hold high cultural and economic importance as ornamentals and for the perfume industry. We report the rose whole genome sequencing and assembly and resequencing of major genotypes that contributed to rose domestication. We generated a homozygous genotype from a heterozygous diploid modern roses progenitor, Rosa chinensis ‘Old Blush’. Using Single Molecule Real-Time sequencing and a meta-assembly approach we obtained one of the most complete plant genomes to date. Diversity analyses highlighted the mosaic origin of ‘La France’, one of the first hybrids combining the growth vigor of European species and recurrent blooming of Chinese species. Genomic segments of Chinese ancestry revealed new candidate genes for recurrent blooming. Reconstructing regulatory and secondary metabolism pathways allowed us to propose a model of interconnected regulation of scent and flower color. This genome provides a foundation for understanding the mechanisms governing rose traits and will accelerate improvement in roses, Rosaceae and ornamentals.

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.

The fate of ribosomal RNA genes in spontaneous polyploid dogrose hybrids [Rosa L. sect. Caninae (DC.) Ser.] exhibiting non-symmetrical meiosis

Dogroses represent an exceptional system for studying the effects of genome doubling and hybridization: their asymmetrical meiosis enables recombination in bi-parentally inherited chromosomes but prevents it in maternally inherited ones. We employed fluorescent in situ hybridization, genome skimming, amplicon sequencing of genomic and cDNA as well as conventional cloning of nuclear ribosomal DNA in two phylogenetically distinct pentaploid (2n = 5x = 35) species, Rosa canina and Rosa inodora, and their naturally occurring reciprocal hybrids, Rosa dumalis (5x) and Rosa agrestis (5x, 6x). Both progenitor species differed in composition, meiotic behaviour and expression of rDNA loci: R. canina (five 18S and 5-8 5S loci) was dominated by the Canina ribotypes, but R. inodora (four 18S loci and 7-8 5S loci) by the Rubiginosa ribotype. The co-localized 5S/18S loci occurred on either bivalent-forming (R. canina) or univalent-forming (R. inodora) chromosomes. Ribosomal DNA loci were additively inherited; however, the Canina ribotypes were dominantly expressed, even in genotypes with relatively low copy number of these genes. Moreover, we observed rDNA homogenization towards the paternally transmitted Canina ribotype in 6x R. agrestis. The here-observed variation in arrangement and composition of rDNA types between R. canina and R. inodora suggests the involvement of different genomes in bivalent formation. This results supports the hypothesis that the asymmetrical meiosis arose at least twice by independent ancient hybridization events.

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.

Mapping a Novel Black Spot Resistance Locus in the Climbing Rose Brite Eyes™ ('RADbrite')

Rose black spot, caused by Diplocarpon rosae, is one of the most devastating foliar diseases of cultivated roses (Rosa spp.). The globally distributed pathogen has the potential to cause large economic losses in the outdoor cultivation of roses. Fungicides are the primary method to manage the disease, but are often viewed unfavorably by home gardeners due to potential environmental and health impacts. As such, rose cultivars with genetic resistance to black spot are highly desired. The tetraploid climbing rose Brite Eyes^TM ('RADbrite') is known for its resistance to black spot. To better characterize the resistance present in Brite Eyes^TM, phenotyping was conducted on a 94 individual F₁ population developed by crossing Brite Eyes^TM to the susceptible tetraploid rose 'Morden Blush'. Brite Eyes^TM was resistant to all D. rosae races evaluated except for race 12. The progeny were either resistant or susceptible to all races (2, 3, 8, 9, 10, 11, and 13) evaluated. The segregation ratio was 1:1 (χ² = 0.3830, P = 0.5360) suggesting resistance is conferred by a single locus. The roses were genotyped with the WagRhSNP 68K Axiom array and the 'polymapR' package was used to construct a map. A single resistance locus (Rdr4) was identified on the long arm of chromosome 5 homoeolog 4. Three resistance loci have been previously identified (Rdr1, Rdr2, and Rdr3). Both Rdr1 and Rdr2 are located on a chromosome 1 homoeolog. The chromosomal location of Rdr3 is unknown, however, races 3 and 9 are virulent on Rdr3. Rdr4 is either a novel gene or an allele of Rdr3 as it provides resistance to races 3 and 9. Due to its broad resistance, Rdr4 is an excellent gene to introgress into new rose cultivars.

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.

An ultra-dense integrated linkage map for hexaploid chrysanthemum enables multi-allelic QTL analysis

KEY MESSAGE

We constructed the first integrated genetic linkage map in a polysomic hexaploid. This enabled us to estimate inheritance of parental haplotypes in the offspring and detect multi-allelic QTL. Construction and use of linkage maps are challenging in hexaploids with polysomic inheritance. Full map integration requires calculations of recombination frequency between markers with complex segregation types. In addition, detection of QTL in hexaploids requires information on all six alleles at one locus for each individual. We describe a method that we used to construct a fully integrated linkage map for chrysanthemum (Chrysanthemum × morifolium, 2n = 6x = 54). A bi-parental F1 population of 406 individuals was genotyped with an 183,000 SNP genotyping array. The resulting linkage map consisted of 30,312 segregating SNP markers of all possible marker dosage types, representing nine chromosomal linkage groups and 107 out of 108 expected homologues. Synteny with lettuce (Lactuca sativa) showed local colinearity. Overall, it was high enough to number the chrysanthemum chromosomal linkage groups according to those in lettuce. We used the integrated and phased linkage map to reconstruct inheritance of parental haplotypes in the F1 population. Estimated probabilities for the parental haplotypes were used for multi-allelic QTL analyses on four traits with different underlying genetic architectures. This resulted in the identification of major QTL that were affected by multiple alleles having a differential effect on the phenotype. The presented linkage map sets a standard for future genetic mapping analyses in chrysanthemum and closely related species. Moreover, the described methods are a major step forward for linkage mapping and QTL analysis in hexaploids.

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.