A SNP-Enabled Assessment of Genetic Diversity, Evolutionary Relationships and the Identification of Candidate Genes in Chrysanthemum

Transcriptome analysis of flower colour reveals the correlation between SNP and differential expression genes in Phalaenopsis

BACKGROUND

Phalaenopsis is an important ornamental plant that has great economic value in the world flower market as one of the most popular flower resources.

OBJECTIVE

To investigate the flower colour formation of Phalaenopsis at the transcription level, the genes involved in flower color formation were identified from RNA-seq in this study.

METHODS

In this study, white and purple petals of Phalaenopsis were collected and analyzed to obtained (1) differential expression genes (DEGs) between white and purple flower color and (2) the association between single nucleotide polymorphisms (SNP) mutations and DEGs at the transcriptome level.

RESULTS

The results indicated that a total of 1,175 DEGs were identified, and 718 and 457 of them were up- and down-regulated genes, respectively. Gene Ontology and pathway enrichment showed that the biosynthesis of the secondary metabolites pathway was key to color formation, and the expression of 12 crucial genes (C4H, CCoAOMT, F3'H, UA3'5'GT, PAL, 4CL, CCR, CAD, CALDH, bglx, SGTase, and E1.11.17) that are involved in the regulation of flower color in Phalaenopsis.

CONCLUSION

This study reported the association between the SNP mutations and DEGs for color formation at RNA level, and provides a new insight to further investigate the gene expression and its relationship with genetic variants from RNA-seq data in other species.

Analyses of a chromosome-scale genome assembly reveal the origin and evolution of cultivated chrysanthemum

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a globally important ornamental plant with great economic, cultural, and symbolic value. However, research on chrysanthemum is challenging due to its complex genetic background. Here, we report a near-complete assembly and annotation for C. morifolium comprising 27 pseudochromosomes (8.15 Gb; scaffold N50 of 303.69 Mb). Comparative and evolutionary analyses reveal a whole-genome triplication (WGT) event shared by Chrysanthemum species approximately 6 million years ago (Mya) and the possible lineage-specific polyploidization of C. morifolium approximately 3 Mya. Multilevel evidence suggests that C. morifolium is likely a segmental allopolyploid. Furthermore, a combination of genomics and transcriptomics approaches demonstrate the C. morifolium genome can be used to identify genes underlying key ornamental traits. Phylogenetic analysis of CmCCD4a traces the flower colour breeding history of cultivated chrysanthemum. Genomic resources generated from this study could help to accelerate chrysanthemum genetic improvement.

Phenotypic differences among and within extant populations of Chrysanthemum arcticum L. and C. a. subsp. arcticum

BACKGROUND

Chrysanthemum arcticum, arctic daisy and its two subspecies (Chrysanthemum arcticum subsp. arcticum, Chrysanthemum arcticum subsp. polaré) are the only chrysanthemum species native to North America. A study on species' variation in morphological and diagnostic traits is important to link morphological traits with previously described single nucleotide polymorphism (SNP) markers, particularly when the genomes are sequenced. The purpose of this study was to establish phenotypic differences and soil conditions among wild C. arcticum and C. a. subsp. arcticum populations, when grown in a uniform environment for two years, for potential linkages with our SNP library. Sixteen quantitative morphological traits and five qualitative morphological traits were investigated for 255 individuals from nine C. arcticum populations and 326 individuals from 21 C. a. subsp. arcticum populations.

RESULTS

In long-day controlled environment, C. arcticum flowering rate was 0% in Year 1, increased to 2.7% in Year 2, while C. a. subsp. arcticum flowering rate was 98.5% in Year 2. Two distinct clusters, distributed by taxonomic classification, were detected by Principal component analysis (PCoA) for 551 individuals from C. arcticum and C. a. subsp. arcticum. Pearson's correlation coefficient analysis indicated a positive and significant correlation between plant height, flower fresh and dry weights. Flower fresh weights were correlated with Δflower weight, while inflorescence length had showed a negative correlation with leaf number. Soil samples had high Na levels along with heavy metals. Thus, the species are salt-tolerant.

CONCLUSION

A high level of salt tolerance (Na) is tolerated by these maritime species which is a unique trait in Chrysanthemum. A new diagnostic trait of inflorescence length was discovered to distinguish among C. arcticum and C. a. subsp. arcticum. Significant flowering differences occurred among the species C. arcticum and C. a. subsp. arcticum under same photoperiodic environment, including flowering rates and visible bud date. This study on the species' variation in morphological and diagnostic traits is of importance to link morphological traits with single nucleotide polymorphism (SNP) markers.

Towards the Improvement of Ornamental Attributes in Chrysanthemum: Recent Progress in Biotechnological Advances

Incessant development and introduction of novel cultivars with improved floral attributes are vital in the dynamic ornamental industry. Chrysanthemum (Chrysanthemum morifolium) is a highly favored ornamental plant, ranking second globally in the cut flower trade, after rose. Development of new chrysanthemum cultivars with improved and innovative modifications in ornamental attributes, including floral color, shape, plant architecture, flowering time, enhanced shelf life, and biotic and abiotic stress tolerance, is a major goal in chrysanthemum breeding. Despite being an economically important ornamental plant, the application of conventional and molecular breeding approaches to various key traits of chrysanthemum is hindered owing to its genomic complexity, heterozygosity, and limited gene pool availability. Although classical breeding of chrysanthemum has resulted in the development of several hundreds of cultivars with various morphological variations, the genetic and transcriptional control of various important ornamental traits remains unclear. The coveted blue colored flowers of chrysanthemums cannot be achieved through conventional breeding and mutation breeding due to technical limitations. However, blue-hued flower has been developed by genetic engineering, and transgenic molecular breeding has been successfully employed, leading to substantial progress in improving various traits. The recent availability of whole-genome sequences of chrysanthemum offers a platform to extensively employ MAS to identify a large number of markers for QTL mapping, and GWAS to dissect the genetic control of complex traits. The combination of NGS, multi-omic platforms, and genome editing technologies has provided a tremendous scope to decipher the molecular and regulatory mechanisms. However, the application and integration of these technologies remain inadequate for chrysanthemum. This review, therefore, details the significance of floral attributes, describes the efforts of recent advancements, and highlights the possibilities for future application towards the improvement of crucial ornamental traits in the globally popular chrysanthemum plant.

Descriptive Genomic Analysis and Sequence Genotyping of the Two Papaya Species (Vasconcellea pubescens and Vasconcellea chilensis) Using GBS Tools

A genotyping by sequencing (GBS) approach was used to analyze the organization of genetic diversity in V. pubescens and V. chilensis. GBS identified 4675 and 4451 SNPs/INDELs in two papaya species. The cultivated orchards of V. pubescens exhibited scarce genetic diversity and low but significant genetic differentiation. The neutrality test yielded a negative and significant result, suggesting that V. pubescens suffered a selective sweep or a rapid expansion after a bottleneck during domestication. In contrast, V. chilensis exhibited a high level of genetic diversity. The genetic differentiation among the populations was slight, but it was possible to distinguish the two genetic groups. The neutrality test indicated no evidence that natural selection and genetic drift affect the natural population of V. chilensis. Using the Carica papaya genome as a reference, we identified critical SNPs/INDELs associated with putative genes. Most of the identified genes are related to stress responses (salt and nematode) and vegetative and reproductive development. These results will be helpful for future breeding and conservation programs of the Caricaceae family.

Genome-wide association study, haplotype analysis, and genomic prediction reveal the genetic basis of yield-related traits in soybean (Glycine max L.)

Identifying the genetic components underlying yield-related traits in soybean is crucial for improving its production and productivity. Here, 211 soybean genotypes were evaluated across six environments for four yield-related traits, including seed yield per plant (SYP), number of pods per plant number of seeds per plant and 100-seed weight (HSW). Genome-wide association study (GWAS) and genomic prediction (GP) analyses were performed using 12,617 single nucleotide polymorphism markers from NJAU 355K SoySNP Array. A total of 57 SNPs were significantly associated with four traits across six environments and a combined environment using five Genome-wide association study models. Out of these, six significant SNPs were consistently identified in more than three environments using multiple GWAS models. The genomic regions (±670 kb) flanking these six consistent SNPs were considered stable QTL regions. Gene annotation and in silico expression analysis revealed 15 putative genes underlying the stable QTLs that might regulate soybean yield. Haplotype analysis using six significant SNPs revealed various allelic combinations regulating diverse phenotypes for the studied traits. Furthermore, the GP analysis revealed that accurate breeding values for the studied soybean traits is attainable at an earlier generation. Our study paved the way for increasing soybean yield performance within a short breeding cycle.

The genus Chrysanthemum: Phylogeny, biodiversity, phytometabolites, and chemodiversity

The ecologically and economically important genus Chrysanthemum contains around 40 species and many hybrids and cultivars. The dried capitulum of Chrysanthemum morifolium (CM) Ramat. Tzvel, i.e., Flos Chrysanthemi, is frequently used in traditional Chinese medicine (TCM) and folk medicine for at least 2,200 years. It has also been a popular tea beverage for about 2,000 years since Han Dynasty in China. However, the origin of different cultivars of CM and the phylogenetic relationship between Chrysanthemum and related Asteraceae genera are still elusive, and there is a lack of comprehensive review about the association between biodiversity and chemodiversity of Chrysanthemum. This article aims to provide a synthetic summary of the phylogeny, biodiversity, phytometabolites and chemodiversity of Chrysanthemum and related taxonomic groups, focusing on CM and its wild relatives. Based on extensive literature review and in light of the medicinal value of chrysanthemum, we give some suggestions for its relationship with some genera/species and future applications. Mining chemodiversity from biodiversity of Chrysanthemum containing subtribe Artemisiinae, as well as mining therapeutic efficacy and other utilities from chemodiversity/biodiversity, is closely related with sustainable conservation and utilization of Artemisiinae resources. There were eight main cultivars of Flos Chrysanthemi, i.e., Hangju, Boju, Gongju, Chuju, Huaiju, Jiju, Chuanju and Qiju, which differ in geographical origins and processing methods. Different CM cultivars originated from various hybridizations between multiple wild species. They mainly contained volatile oils, triterpenes, flavonoids, phenolic acids, polysaccharides, amino acids and other phytometabolites, which have the activities of antimicrobial, anti-viral, antioxidant, anti-aging, anticancer, anti-inflammatory, and closely related taxonomic groups could also be useful as food, medicine and tea. Despite some progresses, the genetic/chemical relationships among varieties, species and relevant genera have yet to be clarified; therefore, the roles of pharmacophylogeny and omics technology are highlighted.

Genetic Diversity and Genome-Wide Association Study of Architectural Traits of Spray Cut Chrysanthemum Varieties

The architecture of spray cut chrysanthemum is crucial for the quality and quantity of cut flower production. However, the mechanism underlying plant architecture still needs to be clarified. In this study, we measured nine architecture-related traits of 195 spray cut chrysanthemum varieties during a two-year period. The results showed that the number of upper primary branches, number of lateral flower buds and primary branch length widely varied. Additionally, plant height had a significant positive correlation with number of leaf nodes and total number of lateral buds. Number of upper primary branches had a significant negative correlation with primary branch diameter, primary branch angle and primary branch length. Plant height, total number of lateral buds, number of upper primary branches, stem diameter, primary branch diameter and primary branch length were vulnerable to environmental impacts. All varieties could be divided into five categories according to cluster analysis, and the typical plant architecture of the varieties was summarized. Finally, a genome-wide association study (GWAS) was performed to find potential functional genes.

Uncovering phylogenetic relationships and genetic diversity of water dropwort using phenotypic traits and SNP markers

The water dropworts Oenanthe linearis Wall. ex DC. and O. javanica (Blume) DC. are aquatic perennial herbs that have been used in China as vegetables and traditional medicines. However, their phylogenetic relationships and genetic diversity are poorly understood. Here, we presented the phenotypic traits and genome-wide DNA marker-based analysis of 158 water dropwort accessions representing both species. The analysis revealed that Oenanthe linearis was readily segregated into linear-leaf and deep-cleft leaf water dropworts according to their leaf shapes at flowering. Oenanthe javanica was classified by clustering analysis into two clusters based mainly on the morphological characteristics of their ultimate segments (leaflets). A set of 11 493 high-quality single-nucleotide polymorphisms was identified and used to construct a phylogenetic tree. There was strong discrimination between O. linearis and O. javanica, which was consistent with their phenotype diversification. The population structure and phylogenetic tree analyses suggested that the O. linearis accessions formed two major groups, corresponding to the linear-leaf and deep-cleft leaf types. The most obvious phenotypic differences between them were fully expressed at the reproductive growth stage. A single-nucleotide polymorphism-based analysis revealed that the O. javanica accessions could be categorized into groups I andII. However, this finding did not entirely align with the clusters revealed by morphological classification. Landraces were clustered into one group along with the remaining wild accessions. Hence, water dropwort domestication was short in duration. The level of genetic diversity for O. linearis (π = 0.1902) was slightly lower than that which was estimated for O. javanica (π = 0.2174). There was a low level of genetic differentiation between O. linearis and O. javanica (Fst = 0.0471). The mean genetic diversity among accessions ranged from 0.1818 for the linear-leaf types to 0.2318 for the groupII accessions. The phenotypic traits and the single-nucleotide polymorphism markers identified here lay empirical foundation for future genomic studies on water dropwort.

Anthocyanins in Floral Colors: Biosynthesis and Regulation in Chrysanthemum Flowers

Chrysanthemum (Chrysanthemum morifolium) is an economically important ornamental crop across the globe. As floral color is the major factor determining customer selection, manipulation of floral color has been a major objective for breeders. Anthocyanins are one of the main pigments contributing to a broad variety of colors in the ray florets of chrysanthemum. Manipulating petal pigments has resulted in the development of a vast range of floral colors. Although the candidate genes involved in anthocyanin biosynthesis have been well studied, the genetic and transcriptional control of floral color remains unclear. Despite advances in multi-omics technology, these methods remain in their infancy in chrysanthemum, owing to its large complex genome and hexaploidy. Hence, there is a need to further elucidate and better understand the genetic and molecular regulatory mechanisms in chrysanthemum, which can provide a basis for future advances in breeding for novel and diverse floral colors in this commercially beneficial crop. Therefore, this review describes the significance of anthocyanins in chrysanthemum flowers, and the mechanism of anthocyanin biosynthesis under genetic and environmental factors, providing insight into the development of novel colored ray florets. Genetic and molecular regulatory mechanisms that control anthocyanin biosynthesis and the various breeding efforts to modify floral color in chrysanthemum are detailed.

Identification of Ear Morphology Genes in Maize (Zea mays L.) Using Selective Sweeps and Association Mapping

The performance of maize hybrids largely depend on two parental inbred lines. Improving inbred lines using artificial selection is a key task in breeding programs. However, it is important to elucidate the effects of this selection on inbred lines. Altogether, 208 inbred lines from two maize heterosis groups, named Shaan A and Shaan B, were sequenced by the genotype-by-sequencing to detect genomic changes under selection pressures. In addition, we completed genome-wide association analysis in 121 inbred lines to identify candidate genes for ear morphology related traits. In a genome-wide selection scan, the inbred lines from Shaan A and Shaan B groups showed obvious population divergences and different selective signals distributed in 337 regions harboring 772 genes. Meanwhile, functional enrichment analysis showed those selected genes are mainly involved in regulating cell development. Interestingly, some ear morphology related traits showed significant differentiation between the inbred lines from the two heterosis groups. The genome-wide association analysis of ear morphology related traits showed that four associated genes were co-localized in the selected regions with high linkage disequilibrium. Our spatiotemporal pattern and gene interaction network results for the four genes further contribute to our understanding of the mechanisms behind ear and fruit length development. This study provides a novel insight into digging a candidate gene for complex traits using breeding materials. Our findings in relation to ear morphology will help accelerate future maize improvement.

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.

Identification of favorable SNP alleles and candidate genes responsible for inflorescence-related traits via GWAS in chrysanthemum

81 SNPs were identified for three inflorescence-related traits, in which 15 were highly favorable. Two dCAPS markers were developed for future MAS breeding, and six candidate genes were predicted. Chrysanthemum is a leading ornamental species worldwide and demonstrates a wealth of morphological variation. Knowledge about the genetic basis of its phenotypic variation for key horticultural traits can contribute to its effective management and genetic improvement. In this study, we conducted a genome-wide association study (GWAS) based on two years of phenotype data and a set of 92,617 single nucleotide polymorphisms (SNPs) using a panel of 107 diverse cut chrysanthemums to dissect the genetic control of three inflorescence-related traits. A total of 81 SNPs were significantly associated with the three inflorescence-related traits (capitulum diameter, number of ray florets and flowering time) in at least one environment, with an individual allele explaining 22.72-38.67% of the phenotypic variation. Fifteen highly favorable alleles were identified for the three target traits by computing the phenotypic effect values for the stable associations detected in 2 year-long trials at each locus. Dosage pyramiding effects of the highly favorable SNP alleles and significant linear correlations between highly favorable allele numbers and corresponding phenotypic performance were observed. Two highly favorable SNP alleles correlating to flowering time and capitulum diameter were converted to derived cleaved amplified polymorphic sequence (dCAPS) markers to facilitate future breeding. Finally, six putative candidate genes were identified that contribute to flowering time and capitulum diameter. These results serve as a foundation for analyzing the genetic mechanisms underlying important horticultural traits and provide valuable insights into molecular marker-assisted selection (MAS) in chrysanthemum breeding programs.

Genome-wide association study identifies favorable SNP alleles and candidate genes for waterlogging tolerance in chrysanthemums

Chrysanthemums are sensitive to waterlogging stress, and the development of screening methods for tolerant germplasms or genes and the breeding of tolerant new varieties are of great importance in chrysanthemum breeding. To understand the genetic basis of waterlogging tolerance (WT) in chrysanthemums, we performed a genome-wide association study (GWAS) using 92,811 single nucleotide polymorphisms (SNPs) in a panel of 88 chrysanthemum accessions, including 64 spray cut and 24 disbud chrysanthemums. The results showed that the average MFVW (membership function value of waterlogging) of the disbud type (0.65) was significantly higher than that of the spray type (0.55) at P < 0.05, and the MFVW of the Asian accessions (0.65) was significantly higher than that of the European accessions (0.48) at P < 0.01. The GWAS performed using the general linear model (GLM) and mixed linear model (MLM) identified 137 and 14 SNP loci related to WT, respectively, and 11 associations were commonly predicted. By calculating the phenotypic effect values for 11 common SNP loci, six highly favorable SNP alleles that explained 12.85-21.85% of the phenotypic variations were identified. Furthermore, the dosage-pyramiding effects of the favorable alleles and the significant linear correlations between the numbers of highly favorable alleles and phenotypic values were identified (r ² = 0.45; P < 0.01). A major SNP locus (Marker6619-75) was converted into a derived cleaved amplified polymorphic sequence (dCAPS) marker that cosegregated with WT with an average efficiency of 78.9%. Finally, four putative candidate genes in the WT were identified via quantitative real-time PCR (qRT-PCR). The results presented in this study provide insights for further research on WT mechanisms and the application of molecular marker-assisted selection (MAS) in chrysanthemum WT breeding programs.

Genome-wide association study identifies favorable SNP alleles and candidate genes for waterlogging tolerance in chrysanthemums

Chrysanthemums are sensitive to waterlogging stress, and the development of screening methods for tolerant germplasms or genes and the breeding of tolerant new varieties are of great importance in chrysanthemum breeding. To understand the genetic basis of waterlogging tolerance (WT) in chrysanthemums, we performed a genome-wide association study (GWAS) using 92,811 single nucleotide polymorphisms (SNPs) in a panel of 88 chrysanthemum accessions, including 64 spray cut and 24 disbud chrysanthemums. The results showed that the average MFVW (membership function value of waterlogging) of the disbud type (0.65) was significantly higher than that of the spray type (0.55) at P < 0.05, and the MFVW of the Asian accessions (0.65) was significantly higher than that of the European accessions (0.48) at P < 0.01. The GWAS performed using the general linear model (GLM) and mixed linear model (MLM) identified 137 and 14 SNP loci related to WT, respectively, and 11 associations were commonly predicted. By calculating the phenotypic effect values for 11 common SNP loci, six highly favorable SNP alleles that explained 12.85-21.85% of the phenotypic variations were identified. Furthermore, the dosage-pyramiding effects of the favorable alleles and the significant linear correlations between the numbers of highly favorable alleles and phenotypic values were identified (r ² = 0.45; P < 0.01). A major SNP locus (Marker6619-75) was converted into a derived cleaved amplified polymorphic sequence (dCAPS) marker that cosegregated with WT with an average efficiency of 78.9%. Finally, four putative candidate genes in the WT were identified via quantitative real-time PCR (qRT-PCR). The results presented in this study provide insights for further research on WT mechanisms and the application of molecular marker-assisted selection (MAS) in chrysanthemum WT breeding programs.

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

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

Investigation of Differences in Fertility among Progenies from Self-Pollinated Chrysanthemum

Most chrysanthemum cultivars are self-incompatible, so it is very difficult to create pure lines that are important in chrysanthemum breeding and theoretical studies. In our previous study, we obtained a self-compatible chrysanthemum cultivar and its self-pollinated seed set was 56.50%. It was interesting that the seed set of its ten progenies ranged from 0% to 37.23%. Examination of the factors causing the differences in the seed set will lead to an improved understanding of chrysanthemum self-incompatibility, and provide valuable information for creating pure lines. Pollen morphology, pollen germination percentage, pistil receptivity and embryo development were investigated using the in vitro culture method, the paraffin section technique, scanning electron microscopy and transmission electron microscopy. Moreover, RNA sequencing and bioinformatics were applied to analyzing the transcriptomic profiles of mature stigmas and anthers. It was found that the self-pollinated seed set of "Q10-33-1①","Q10-33-1③","Q10-33-1④" and "Q10-33-1⑩" were 37.23%, 26.77%, 7.97% and 0%, respectively. The differences in fertility among four progenies were mainly attributable to differences in pollen germination percentage and pistil receptivity. Failure of the seed set in "Q10-33-1⑩" was possibly due to self-incompatibility. In the transcriptomic files, 22 potential stigma S genes and 8 potential pollen S genes were found out.