Using RNA-Seq to assemble a rose transcriptome with more than 13,000 full-length expressed genes and to develop the WagRhSNP 68k Axiom SNP array for rose (Rosa L.)

Genomic prediction and genome-wide association study using combined genotypic data from different genotyping systems: application to apple fruit quality traits

With advances in next-generation sequencing technologies, various marker genotyping systems have been developed for genomics-based approaches such as genomic selection (GS) and genome-wide association study (GWAS). As new genotyping platforms are developed, data from different genotyping platforms must be combined. However, the potential use of combined data for GS and GWAS has not yet been clarified. In this study, the accuracy of genomic prediction (GP) and the detection power of GWAS increased for most fruit quality traits of apples when using combined data from different genotyping systems, Illumina Infinium single-nucleotide polymorphism array and genotyping by random amplicon sequencing-direct (GRAS-Di) systems. In addition, the GP model, which considered the inbreeding effect, further improved the accuracy of the seven fruit traits. Runs of homozygosity (ROH) islands overlapped with the significantly associated regions detected by the GWAS for several fruit traits. Breeders may have exploited these regions to select promising apples by breeders, increasing homozygosity. These results suggest that combining genotypic data from different genotyping platforms benefits the GS and GWAS of fruit quality traits in apples. Information on inbreeding could be beneficial for improving the accuracy of GS for fruit traits of apples; however, further analysis is required to elucidate the relationship between the fruit traits and inbreeding depression (e.g. decreased vigor).

The role of invasive plant species in drought resilience in agriculture: the case of sweet briar (Rosa rubiginosa L.)

Sweet briar (Rosa rubiginosa) belongs to the group of wild roses. Under natural conditions it grows throughout Europe, and was introduced also into the southern hemisphere, where it has efficiently adapted to dry lands. This review focuses on the high adaptation potential of sweet briar to soil drought in the context of global climatic changes, especially considering steppe formation and desertification of agricultural, orchard, and horticultural areas. We provide a comprehensive overview of current knowledge on sweet briar traits associated with drought tolerance and particularly water use efficiency, sugar accumulation, accumulation of CO2 in intercellular spaces, stomatal conductance, gibberellin level, effective electron transport between photosystem II and photosystem I, and protein content. We discuss the genetics and potential applications in plant breeding and suggest future directions of study concerning invasive populations of R. rubiginosa. Finally, we point out that sweet briar can provide new genes for breeding in the context of depleting gene pools of the crop plants.

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.

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.

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.

Detection of Reproducible Major Effect QTL for Petal Traits in Garden Roses

The detection of QTL by association genetics depends on the genetic architecture of the trait under study, the size and structure of the investigated population and the availability of phenotypic and marker data of sufficient quality and quantity. In roses, we previously demonstrated that major QTL could already be detected in small association panels. In this study, we analyzed petal number, petal size and fragrance in a small panel of 95 mostly tetraploid garden rose genotypes. After genotyping the panel with the 68 K Axiom WagRhSNP chip we detected major QTL for all three traits. Each trait was significantly influenced by several genomic regions. Some of the QTL span genomic regions that comprise several candidate genes. Selected markers from some of these regions were converted into KASP markers and were validated in independent populations of up to 282 garden rose genotypes. These markers demonstrate the robustness of the detected effects independent of the set of genotypes analyzed. Furthermore, the markers can serve as tools for marker-assisted breeding in garden roses. Over an extended timeframe, they may be used as a starting point for the isolation of the genes underlying the QTL.

Recent Large-Scale Genotyping and Phenotyping of Plant Genetic Resources of Vegetatively Propagated Crops

Several recent national and international projects have focused on large-scale genotyping of plant genetic resources in vegetatively propagated crops like fruit and berries, potatoes and woody ornamentals. The primary goal is usually to identify true-to-type plant material, detect possible synonyms, and investigate genetic diversity and relatedness among accessions. A secondary goal may be to create sustainable databases that can be utilized in research and breeding for several years ahead. Commonly applied DNA markers (like microsatellite DNA and SNPs) and next-generation sequencing each have their pros and cons for these purposes. Methods for large-scale phenotyping have lagged behind, which is unfortunate since many commercially important traits (yield, growth habit, storability, and disease resistance) are difficult to score. Nevertheless, the analysis of gene action and development of robust DNA markers depends on environmentally controlled screening of very large sets of plant material. Although more time-consuming, co-operative projects with broad-scale data collection are likely to produce more reliable results. In this review, we will describe some of the approaches taken in genotyping and/or phenotyping projects concerning a wide variety of vegetatively propagated crops.

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.

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.

An Axiom SNP genotyping array for Douglas-fir

BACKGROUND

In forest trees, genetic markers have been used to understand the genetic architecture of natural populations, identify quantitative trait loci, infer gene function, and enhance tree breeding. Recently, new, efficient technologies for genotyping thousands to millions of single nucleotide polymorphisms (SNPs) have finally made large-scale use of genetic markers widely available. These methods will be exceedingly valuable for improving tree breeding and understanding the ecological genetics of Douglas-fir, one of the most economically and ecologically important trees in the world.

RESULTS

We designed SNP assays for 55,766 potential SNPs that were discovered from previous transcriptome sequencing projects. We tested the array on ~ 2300 related and unrelated coastal Douglas-fir trees (Pseudotsuga menziesii var. menziesii) from Oregon and Washington, and 13 trees of interior Douglas-fir (P. menziesii var. glauca). As many as ~ 28 K SNPs were reliably genotyped and polymorphic, depending on the selected SNP call rate. To increase the number of SNPs and improve genome coverage, we developed protocols to 'rescue' SNPs that did not pass the default Affymetrix quality control criteria (e.g., 97% SNP call rate). Lowering the SNP call rate threshold from 97 to 60% increased the number of successful SNPs from 20,669 to 28,094. We used a subset of 395 unrelated trees to calculate SNP population genetic statistics for coastal Douglas-fir. Over a range of call rate thresholds (97 to 60%), the median call rate for SNPs in Hardy-Weinberg equilibrium ranged from 99.2 to 99.7%, and the median minor allele frequency ranged from 0.198 to 0.233. The successful SNPs also worked well on interior Douglas-fir.

CONCLUSIONS

Based on the original transcriptome assemblies and comparisons to version 1.0 of the Douglas-fir reference genome, we conclude that these SNPs can be used to genotype about 10 K to 15 K loci. The Axiom genotyping array will serve as an excellent foundation for studying the population genomics of Douglas-fir and for implementing genomic selection. We are currently using the array to construct a linkage map and test genomic selection in a three-generation breeding program for coastal Douglas-fir.

Characterization of Indian bred rose cultivars using morphological and molecular markers for conservation and sustainable management

Rose (Rosa × hybrid L.) is one of the most important commercial ornamental crops cultivated worldwide for its beauty, fragrance and nutraceutical values. Characterization of rose germplasm provides precise information about the extent of diversity present among the cultivars. It also helps in cultivar identification, intellectual property right protection, variety improvement and genetic diversity conservation. In the present study, 109 Indian bred rose cultivars were characterized using 59 morphological and 48 SSR markers. Out of 48 SSRs used, 31 markers exhibited polymorphism and 96 alleles were identified with an average of 3.9 alleles per locus. Nei's expected heterozygosity value of each locus ranged from 0.08 (with SSR ABRII/RPU32) to 0.78 (SSR Rh58). The similarity coefficient values ranged from 0.42 to 0.90 which indicated presence of moderated diversity among Indian cultivars. The neighbor-joining tree based on morphological data grouped the cultivars into two major clusters and several minor clusters based on their morphological resemblance. However, UPGMA dendrogram constructed using matching coefficient values grouped the cultivars into eight different clusters. Interpopulation analysis revealed higher genetic similarities between Hybrid Tea and Floribunda cultivars. An analysis for presence of population sub-structure grouped the Indian cultivars into eight different genetic groups. Analysis of molecular variance revealed apportioning of 97.59% of the variation to within subgroup diversity and 3.07% to between the cultivar groups. We have demonstrated here successful utilization of robust SSR to distinguish cultivars and assess genetic diversity among Indian bred rose cultivars. The information provided here is useful for cultivar identification and protection, cultivar improvement and genetic diversity conservation.

Expected Genotype Quality and Diploidized Marker Data from Genotyping-by-Sequencing of Urochloa spp. Tetraploids

Introduced concept of expected genotype quality (EGQ) and software to calculate it Provided read depth guidelines for GBS in tetraploids Developed software to generate diploidized genotype calls from VCF files Demonstrated value of aligning GBS reads to a mock reference genome for SNP discovery Recommend filtering based on GQ and read depth to prevent genotype bias Although genotyping-by-sequencing (GBS) is a well-established marker technology in diploids, the development of best practices for tetraploid species is a topic of current research. We determined the theoretical relationship between read depth and the phred-scaled probability of genotype misclassification conditioned on the true genotype, which we call expected genotype quality (EGQ). If the GBS method has 0.5% allelic error, then 17 reads are needed to classify simplex tetraploids as heterozygous with 95% accuracy (EGQ = 13) vs. 61 reads to determine allele dosage. We developed an R script to convert tetraploid GBS data in variant call format (VCF) into diploidized genotype calls and applied it to 267 interspecific hybrids of the tetraploid forage grass Urochloa. When reads were aligned to a mock reference genome created from GBS data of the Urochloa brizantha (Hochst. ex A. Rich.) R. D. Webster cultivar Marandu, 25,678 biallelic single nucleotide polymorphism (SNPs) were discovered, compared with ∼3000 SNPs when aligning to the closest true reference genomes, Setaria viridis (L.) P. Beauv. and S. italica (L.) P. Beauv. Cross-validation revealed that missing genotypes were imputed by the random forest method with a median accuracy of 0.85 regardless of heterozygote frequency. Using the Urochloa spp. hybrids, we illustrated how filtering samples based only on genotype quality (GQ) creates genotype bias; a depth threshold based on EGQ is also needed regardless of whether genotypes are called using a diploidized or allele dosage model.

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.

Hybrid-Transcriptome Sequencing and Associated Metabolite Analysis Reveal Putative Genes Involved in Flower Color Difference in Rose Mutants

Gene mutation is a common phenomenon in nature that often leads to phenotype differences, such as the variations in flower color that frequently occur in roses. With the aim of revealing the genomic information and inner mechanisms, the differences in the levels of both transcription and secondary metabolism between a pair of natural rose mutants were investigated by using hybrid RNA-sequencing and metabolite analysis. Metabolite analysis showed that glycosylated derivatives of pelargonidin, e.g., pelargonidin 3,5 diglucoside and pelargonidin 3-glucoside, which were not detected in white flowers (Rosa 'Whilte Mrago Koster'), constituted the major pigments in pink flowers. Conversely, the flavonol contents of petal, such as kaempferol-3-glucoside, quercetin 3-glucoside, and rutin, were higher in white flowers. Hybrid RNA-sequencing obtained a total of 107,280 full-length transcripts in rose petal which were annotated in major databases. Differentially expressed gene (DEG) analysis showed that the expression of genes involved in the flavonoid biosynthesis pathway was significantly different, e.g., CHS, FLS, DFR, LDOX, which was verified by qRT-PCR during flowering. Additionally, two MYB transcription factors were found and named RmMYBAN2 and RmMYBPA1, and their expression patterns during flowering were also analyzed. These findings indicate that these genes may be involved in the flower color difference in the rose mutants, and competition between anthocyanin and flavonol biosynthesis is a primary cause of flower color variation, with its regulation reflected by transcriptional and secondary metabolite levels.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

CYP79D73 Participates in Biosynthesis of Floral Scent Compound 2-Phenylethanol in Plumeria rubra

Plumeria (Plumeria rubra), well known for its brightly colored and fragrant flowers, emits a number of floral volatile organic compounds (VOCs). Plumeria flowers emit a total of 43 VOCs including nine phenylpropanoids/benzenoids, such as 2-phenylethanol (2PE), benzaldehyde, 2-phenylacetaldehyde (PAld), (E/Z)-phenylacetaldoxime (PAOx), benzyl nitrile (BN), and 2-phenylnitroethane (PN). To identify genes and pathways involved in the production of the major compound 2PE, we analyzed the plumeria floral transcriptome and found a highly expressed, flower-specific gene encoding a cytochrome P450 family 79D protein (PrCYP79D73), which catalyzed the formation of (E/Z)-PAOx. Feeding experiments with deuterated phenylalanine or deuterated (E/Z)-PAOx showed that (E/Z)-PAOx is an intermediate in the biosynthesis of 2PE, as are two nitrogen-containing volatiles, BN and PN, in plumeria flowers. Crude enzyme extracts from plumeria flowers converted l-phenylalanine to (E/Z)-PAOx, PAld, 2PE, BN, and PN. The biosynthesis of these compounds increased with addition of PrCYP79D73-enriched microsomes but was blocked by pretreatment with 4-phenylimidazole, an inhibitor of cytochrome P450 enzymes. Moreover, overexpression of PrCYP79D73 in Nicotiana benthamiana resulted in the emission of (E/Z)-PAOx as well as PAld, 2PE, BN, and PN, all of which were also found among plumeria floral VOCs. Taken together, our results demonstrate that PrCYP79D73 is a crucial player in the biosynthesis of the major floral VOC 2PE and other nitrogen-containing volatiles. These volatiles may be required for plant defense as well as to attract pollinators for the successful reproduction of plumeria.

15 years of GDR: New data and functionality in the Genome Database for Rosaceae

The Genome Database for Rosaceae (GDR, https://www.rosaceae.org) is an integrated web-based community database resource providing access to publicly available genomics, genetics and breeding data and data-mining tools to facilitate basic, translational and applied research in Rosaceae. The volume of data in GDR has increased greatly over the last 5 years. The GDR now houses multiple versions of whole genome assembly and annotation data from 14 species, made available by recent advances in sequencing technology. Annotated and searchable reference transcriptomes, RefTrans, combining peer-reviewed published RNA-Seq as well as EST datasets, are newly available for major crop species. Significantly more quantitative trait loci, genetic maps and markers are available in MapViewer, a new visualization tool that better integrates with other pages in GDR. Pathways can be accessed through the new GDR Cyc Pathways databases, and synteny among the newest genome assemblies from eight species can be viewed through the new synteny browser, SynView. Collated single-nucleotide polymorphism diversity data and phenotypic data from publicly available breeding datasets are integrated with other relevant data. Also, the new Breeding Information Management System allows breeders to upload, manage and analyze their private breeding data within the secure GDR server with an option to release data publicly.

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.

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.

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.

Comparative RNA-Seq analysis reveals a critical role for brassinosteroids in rose (Rosa hybrida) petal defense against Botrytis cinerea infection

Background

One of the most popular ornamental plants worldwide, roses (Rosa sp.), are very susceptible to Botrytis gray mold disease. The necrotrophic infection of rose petals by B. cinerea causes the collapse and death of these tissues in both the growth and post-harvest stages, resulting in serious economic losses. To understand the molecular basis of rose resistance against B. cinerea, we profiled the petal transcriptome using RNA-Seq technology.

Results

We identified differentially transcribed genes (DTGs) in petals during B. cinerea infection at 30 h post inoculation (hpi) and/or 48 hpi. Gene ontology term enrichment and pathway analyses revealed that metabolic, secondary metabolite biosynthesis, plant-pathogen interaction, and plant hormone signal transduction pathways were involved. The expression of 370 cell-surface immune receptors was upregulated during infection. In addition, 188 genes encoding transcription factors were upregulated, particularly in the ERF, WRKY, bHLH, MYB, and NAC families, implying their involvement in resistance against B. cinerea. We further identified 325 upregulated DTGs in the hormone signal transduction pathways. Among them, the brassinosteroid (BR)-related genes were the most significantly enriched. To confirm the role of BR in Botrytis resistance, exogenous BR was applied to rose flowers before the inoculation of B. cinerea, which enhanced the defense response in these petals.

Conclusions

Our global transcriptome profiling provides insights into the complex gene regulatory networks mediating the rose petal response to B. cinerea. We further demonstrated the role of the phytohormone BR in the resistance of petals to necrotrophic fungal pathogens.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0668-x) contains supplementary material, which is available to authorized users.

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.

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.

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.

Revolution in Genotyping Platforms for Crop Improvement

In the past decade, the application of high-throughput sequencing to crop genotyping has given rise to novel platforms capable of genotyping tens of thousands of genome-wide DNA markers. Coupled with the decreasing costs of sequencing, this rapid increase in markers allows accelerated and highly accurate genotyping of entire crop populations and diversity sets using single nucleotide polymorphisms (SNPs). These revolutionary advances accelerate crop improvement by facilitating a more precise connection of phenotype to genotype through association studies, linkage mapping and diversity analysis. The platforms driving the advances in genotyping are array technologies and genotyping by sequencing (GBS) methods, which include both low-coverage whole genome resequencing (skim sequencing) and reduced representation sequencing (RRS) approaches. Here, we outline and compare these genotyping platforms and provide a perspective on the promising future of crop genotyping. While SNP arrays provide high quality, simple handling, and unchallenging analysis, the lower cost of RRS and the greater data volume produced by skim sequencing suggest that use of GBS will become more prevalent in crop genomics as sequencing costs decrease and data analysis becomes more streamlined. Graphical Abstract.

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.

MIKCC-type MADS-box genes in Rosa chinensis: the remarkable expansion of ABCDE model genes and their roles in floral organogenesis

MIKC^C-type MADS-box (MIKC^C) genes encode transcription factors that have crucial roles in controlling floral organogenesis and flowering time in plants. Although this gene family has been well characterized in many plant species, its evolutionary and comprehensive functional analysis in rose is lacking. In this study, 58 non-redundant MIKC^C uni-transcripts were extensively identified from rose transcriptomes. Phylogenetic analysis placed these genes into 12 clades with their Arabidopsis and strawberry counterparts, and revealed that ABCDE model (including AP1/FUL, AP3/PI, AG, and SEP clades), and SOC1 and AGL6 clade genes have remarkably expanded in Rosa chinensis, whereas genes from the FLC and AGL17 clades were undetectable. Sequence alignments suggest that the AP3/PI clade may contribute to more specific functions in rose due to a high variation of amino acid residues within its MADS-box domains. A comparative analysis of gene expression in specific floral organ differentiation stages and floral organs between R. chinensis cv. Old Blush and the closely related mutant genotype R. chinensis cv. Viridiflora (floral organs mutated into leaf-like structures) further revealed the roles of ABCDE model genes during floral organogenesis in rose. Analysis of co-expression networks provided an overview of the regulatory mechanisms of rose MIKC^C genes and shed light on both the prominent roles of AP3/PI clade genes in floral organogenesis and the roles of RcAGL19, RcAGL24, and RcSOC1 in regulating floral transition in rose. Our analyses provide an overall insight of MIKC^C genes in rose and their potential roles in floral organogenesis.

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.

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

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

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.

Automated tetraploid genotype calling by hierarchical clustering

New software to make tetraploid genotype calls from SNP array data was developed, which uses hierarchical clustering and multiple F1 populations to calibrate the relationship between signal intensity and allele dosage. SNP arrays are transforming breeding and genetics research for autotetraploids. To fully utilize these arrays, the relationship between signal intensity and allele dosage must be calibrated for each marker. We developed an improved computational method to automate this process, which is provided as the R package ClusterCall. In the training phase of the algorithm, hierarchical clustering within an F1 population is used to group samples with similar intensity values, and allele dosages are assigned to clusters based on expected segregation ratios. In the prediction phase, multiple F1 populations and the prediction set are clustered together, and the genotype for each cluster is the mode of the training set samples. A concordance metric, defined as the proportion of training set samples equal to the mode, can be used to eliminate unreliable markers and compare different algorithms. Across three potato families genotyped with an 8K SNP array, ClusterCall scored 5729 markers with at least 0.95 concordance (94.6% of its total), compared to 5325 with the software fitTetra (82.5% of its total). The three families were used to predict genotypes for 5218 SNPs in the SolCAP diversity panel, compared with 3521 SNPs in a previous study in which genotypes were called manually. One of the additional markers produced a significant association for vine maturity near a well-known causal locus on chromosome 5. In conclusion, when multiple F1 populations are available, ClusterCall is an efficient method for accurate, autotetraploid genotype calling that enables the use of SNP data for research and plant breeding.

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.

Gene expression and physiological responses associated to stomatal functioning in Rosa×hybrida grown at high relative air humidity

High relative air humidity (RH≥85%) during growth disturbs stomatal functioning, resulting in excessive water loss in conditions of high evaporative demand. We investigated the expression of nine abscisic acid (ABA)-related genes (involved in ABA biosynthesis, oxidation and conjugation) and two non-ABA related genes (involved in the water stress response) aiming to better understand the mechanisms underlying contrasting stomatal functioning in plants grown at high RH. Four rose genotypes with contrasting sensitivity to high RH (one sensitive, one tolerant and two intermediate) were grown at moderate (62±3%) or high (89±4%) RH. The sensitive genotype grown at high RH showed a significantly higher stomatal conductance (g_s) and water loss in response to closing stimuli as compared to the other genotypes. Moreover, high RH reduced the leaf ABA concentration and its metabolites to a greater extent in the sensitive genotype as compared to the tolerant one. The large majority of the studied genes had a relevant role on stomatal functioning (NCED1, UGT75B2, BG2, OST1, ABF3 and Rh-APX) while two others showed a minor contribution (CYP707A3 and BG1) and AAO3, CYP707A1 and DREB1B did not contribute to the tolerance trait. These results show that multiple genes form a highly complex regulatory network acting together towards the genotypic tolerance to high RH.

Current status and biotechnological advances in genetic engineering of ornamental plants

Cut flower markets are developing in many countries as the international demand for cut flowers is rapidly growing. Developing new varieties with modified characteristics is an important aim in floriculture. Production of transgenic ornamental plants can shorten the time required in the conventional breeding of a cultivar. Biotechnology tools in combination with conventional breeding methods have been used by cut flower breeders to change flower color, plant architecture, post-harvest traits, and disease resistance. In this review, we describe advances in genetic engineering that have led to the development of new cut flower varieties.

Development and validation of the Axiom(®) Apple480K SNP genotyping array

Cultivated apple (Malus × domestica Borkh.) is one of the most important fruit crops in temperate regions, and has great economic and cultural value. The apple genome is highly heterozygous and has undergone a recent duplication which, combined with a rapid linkage disequilibrium decay, makes it difficult to perform genome-wide association (GWA) studies. Single nucleotide polymorphism arrays offer highly multiplexed assays at a relatively low cost per data point and can be a valid tool for the identification of the markers associated with traits of interest. Here, we describe the development and validation of a 487K SNP Affymetrix Axiom(®) genotyping array for apple and discuss its potential applications. The array has been built from the high-depth resequencing of 63 different cultivars covering most of the genetic diversity in cultivated apple. The SNPs were chosen by applying a focal points approach to enrich genic regions, but also to reach a uniform coverage of non-genic regions. A total of 1324 apple accessions, including the 92 progenies of two mapping populations, have been genotyped with the Axiom(®) Apple480K to assess the effectiveness of the array. A large majority of SNPs (359 994 or 74%) fell in the stringent class of poly high resolution polymorphisms. We also devised a filtering procedure to identify a subset of 275K very robust markers that can be safely used for germplasm surveys in apple. The Axiom(®) Apple480K has now been commercially released both for public and proprietary use and will likely be a reference tool for GWA studies in apple.

Transcriptome Analysis of Gerbera hybrida Including in silico Confirmation of Defense Genes Found

For the ornamental crop Gerbera hybrida, breeding at the moment is done using conventional methods. As this has drawbacks in breeding speed and efficiency, especially for complex traits like disease resistance, we set out to develop genomic resources. The leaf and flower bud transcriptomes of four parents, used to generate two gerbera populations, were sequenced using Illumina paired-end sequencing. In total, 36,770 contigs with an average length of 1397 bp were generated and these have been the starting point for SNP identification and annotation. The consensus contig sequences were used to map reads of individual parents, to identify genotype specific SNPs, and to assess the presence of common SNPs between genotypes. Comparison with the non-redundant protein database (nr) showed that 29,146 contigs gave BLAST hits. Of sequences with blast results, 73.3% obtained a clear gene ontology (GO) annotation. EST contigs coding for enzymes were found in Kyoto Encyclopedia of Genes and Genomes maps (KEGG). Through, these annotated data and KEGG molecular interaction network, transcripts associated with the phenylpropanoid metabolism, other secondary metabolite biosynthesis pathways, phytohormone biosynthesis and signal transduction were analyzed in more detail. Identifying genes involved in these processes could provide genetic and genomic resources for studying the mechanism of disease resistance in gerbera.

Inheritance and QTL analysis of the determinants of flower color in tetraploid cut roses

The success of cut rose cultivars is a direct result of their aesthetic value. The rose industry thrives on novelty, and the production of novel flower color has been extensively studied. The most popular color is red, and it is, therefore, important for breeders to produce a good red cultivar. The final visible color of the flower is a combination of a number of factors including the type of anthocyanin accumulating, modifications to the anthocyanidin molecule, co-pigmentation and vacuolar pH. Here, we analyze the quantitative variation of the biochemical constituents of flower color in a tetraploid rose population and combine this with marker information in the segregating rose population to map the chromosomal locations of putative QTLs for flower color traits. Within our tetraploid population, we found a number of QTLs that were mapped on ICM 1, 2, 6 and 7. We were able to show the effect of the different QTLs on the final visible color of the flower from salmon to dark red.

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

Dense genetic maps create a base for QTL analysis of important traits and future implementation of marker-assisted breeding. In tetraploid rose, the existing linkage maps include <300 markers to cover 28 linkage groups (4 homologous sets of 7 chromosomes). Here we used the 68k WagRhSNP Axiom single-nucleotide polymorphism (SNP) array for rose, in combination with SNP dosage calling at the tetraploid level, to genotype offspring from the garden rose cultivar 'Red New Dawn'. The offspring proved to be not from a single bi-parental cross. In rose breeding, crosses with unintended parents occur regularly. We developed a strategy to separate progeny into putative populations, even while one of the parents was unknown, using principle component analysis on pairwise genetic distances based on sets of selected SNP markers that were homozygous, and therefore uninformative for one parent. One of the inferred populations was consistent with self-fertilization of 'Red New Dawn'. Subsequently, linkage maps were generated for a bi-parental and a self-pollinated population with 'Red New Dawn' as the common maternal parent. The densest map, for the selfed parent, had 1929 SNP markers on 25 linkage groups, covering 1765.5 cM at an average marker distance of 0.9 cM. Synteny with the strawberry (Fragaria vesca) genome was extensive. Rose ICM1 corresponded to F. vesca pseudochromosome 7 (Fv7), ICM4 to Fv4, ICM5 to Fv3, ICM6 to Fv2 and ICM7 to Fv5. Rose ICM2 corresponded to parts of F. vesca pseudochromosomes 1 and 6, whereas ICM3 is syntenic to the remainder of Fv6.

Genome-Wide Association Analysis of the Anthocyanin and Carotenoid Contents of Rose Petals

Petal color is one of the key characteristics determining the attractiveness and therefore the commercial value of an ornamental crop. Here, we present the first genome-wide association study for the important ornamental crop rose, focusing on the anthocyanin and carotenoid contents in petals of 96 diverse tetraploid garden rose genotypes. Cultivated roses display a vast phenotypic and genetic diversity and are therefore ideal targets for association genetics. For marker analysis, we used a recently designed Axiom SNP chip comprising 68,000 SNPs with additionally 281 SSRs, 400 AFLPs and 246 markers from candidate genes. An analysis of the structure of the rose population revealed three subpopulations with most of the genetic variation between individual genotypes rather than between clusters and with a high average proportion of heterozygous loci. The mapping of markers significantly associated with anthocyanin and carotenoid content to the related Fragaria and Prunus genomes revealed clusters of associated markers indicating five genomic regions associated with the total anthocyanin content and two large clusters associated with the carotenoid content. Among the marker clusters associated with the phenotypes, we found several candidate genes with known functions in either the anthocyanin or the carotenoid biosynthesis pathways. Among others, we identified a glutathione-S-transferase, 4CL, an auxin response factor and F3'H as candidate genes affecting anthocyanin concentration, and CCD4 and Zeaxanthine epoxidase as candidates affecting the concentration of carotenoids. These markers are starting points for future validation experiments in independent populations as well as for functional genomic studies to identify the causal factors for the observed color phenotypes. Furthermore, validated markers may be interesting tools for marker-assisted selection in commercial breeding programmes in that they provide the tools to identify superior parental combinations that combine several associated markers in higher dosages.