Genomic analyses provide insights into sex differentiation of tetraploid strawberry (Fragaria moupinensis).

The strawberry genus, Fragaria, exhibits a wide range of sexual systems and natural ploidy variation. Nearly, all polyploid strawberry species exhibit separate sexes (dioecy). Research has identified the sex-determining sequences as roughly conserved but with repeatedly changed genomic locations across octoploid strawberries. However, it remains unclear whether tetraploid wild strawberries evolved dioecy independently or shared a common origin with octoploid strawberries. In this study, we investigated the sex determinants of F. moupinensis, a dioecious plant with heterogametic females (ZW). Utilizing a combination of haplotype-resolved genome sequencing of the female F. moupinensis, k-mer-based and coverage-based genome-wide association studies (GWAS), and transcriptomic analysis, we discovered a non-recombining, approximately 33.6 kb W-specific region on chromosome 2a. Within this region, only one candidate sex-determining gene (FmoAFT) was identified. Furthermore, an extensive resequencing of the entire Fragaria genus indicated that the W-specific region displays conservative female specificity across all tetraploid species. This observation suggests that dioecy evolved independently in tetraploid and octoploid strawberries. Moreover, employing virus-induced gene silencing (VIGS), we knocked down the expression of the FmoAFT homologue transcript in cultivated strawberries, revealing its potential role in promoting female functions during early carpel development. We also applied DNA affinity purification sequencing (DAP-seq) and yeast one-hybrid assays to identify potential direct targets of FmoAFT. These insights shed new light on the genetic basis and evolutionary history of sex determination in strawberries, thereby facilitating the formulation of strategies to manipulate sex determination in breeding programs.

Segmented Flow Strategies for Integrating Liquid Chromatography-Mass Spectrometry with Nuclear Magnetic Resonance for Lipidomics.

Building an accurate lipid inventory relies on coordinated information from orthogonal analytical capabilities. Integrating the familiar workflow of liquid chromatography (LC), high-resolution mass spectrometry (HRMS), and tandem mass spectrometry (MS/MS) with proton nuclear magnetic resonance spectroscopy (1H NMR) would be ideal for building that inventory. For absolute lipid structural elucidation, LC-HRMS/MS can provide lower-level structural information with superior sensitivity, while 1H NMR can provide invaluable higher-order structural information for the disambiguation of isomers with absolute chemical specificity. Digitization of the LC eluent followed by splitting the microfractions into two flow paths in a defined ratio for HRMS/MS and NMR would be the ideal strategy to permit correlation of the MS and NMR data as a function of chromatographic retention time. Here, we report an active segmentation platform to transform analytical flow rate LC eluent into parallel microliter segmented flow queues for high confidence correlation of the MS, MS/MS, and NMR data. The practical details in implementing this strategy to achieve an integrated LC-MS-NMR platform are presented, including the development of an active segmentation technology using a four-port two-way valve to transform the LC eluent into parallel segmented flows for online MS analysis followed by offline segment-specific 1H NMR and optimization of the detector response toward segmented flow. To demonstrate the practicality of this novel platform, it was tested using lipid mixture samples.

Effects of coastal aquaculture on sediment organic matter: Assessed with multiple spectral and isotopic indices.

Sediment organic matter (SOM) is important in the biogeochemical cycling of carbon, nutrients, and pollutants in the coastal environment, which is increasingly disturbed by aquaculture that is particularly intense in China. However, the identification of aquaculture signals in SOM is rather challenging in the complex coastal environment that receives materials from a variety of sources. This was studied in a typical culture area of shellfish and algae in SE China from July 2019 to October 2020, using a combination of elemental (OC, TN, N/C), isotopic (δ13C and δ15N), spectral (absorption spectroscopy and fluorescence EEMs-PARAFAC), and statistical analysis (principal component analysis, PCA). All indices of SOM quantity and several spectral indices for the SOM composition correlated significantly with grain size, with lower SOM quantity and higher autochthonous contribution in coarse sediments. The strong correlations between elemental and spectral indices suggested that optical analysis could provide valuable indices for assessing the quantity of bulk organic matter. The comparison of SOM indices between different zones and between different months showed an overall limited influence of shellfish and laver culture. This indicated the sustainability of these types of aquaculture that require no manual addition of feeds and thus are generally clean. The further applications of end-member mixing analysis using the IsoSource program and PCA were more sensitive, which identified the removal of SOM by shellfish in the growing season and the contribution from shellfish residuals after the harvest and the cultured laver at some locations. Overall, our results have implications for a better understanding of the biogeochemical processes and ecosystem sustainability in the coastal environment under intense aquaculture activities.

Fragaria mitogenomes evolve rapidly in structure but slowly in sequence and incur frequent multinucleotide mutations mediated by microinversions.

Plant mitochondrial DNA has been described as evolving rapidly in structure but slowly in sequence. However, many of the noncoding portions of plant mitogenomes are not homologous among species, raising questions about the rate and spectrum of mutations in noncoding regions. Recent studies have suggested that the lack of homology in noncoding regions could be due to increased sequence divergence. We compared 30 kb of coding and 200 kb of noncoding DNA from 13 sequenced Fragaria mitogenomes, followed by analysis of the rate of sequence divergence, microinversion events and structural variations. Substitution rates in synonymous sites and nongenic sites are nearly identical, suggesting that the genome-wide point mutation rate is generally consistent. A surprisingly high number of large multinucleotide substitutions were detected in Fragaria mitogenomes, which may have resulted from microinversion events and could affect phylogenetic signal and local rate estimates. Fragaria mitogenomes preferentially accumulate deletions relative to insertions and substantial genomic arrangements, whereas mutation rates could positively associate with these sequence and structural changes among species. Together, these observations suggest that plant mitogenomes exhibit low point mutations genome-wide but exceptionally high structural variations, and our results favour a gain-and-loss model for the rapid loss of homology among plant mitogenomes.

Comparative Transcriptome Analysis Uncovers the Regulatory Roles of MicroRNAs Involved in Petal Color Change of Pink-Flowered Strawberry.

The pink-flowered strawberry is popular in China due to its high ornamental value. In the present study, sRNAome, transcriptome, and degradome sequencing were performed to understand the functions of microRNAs (miRNAs) and their target genes during flower development in pink-flowered strawberry. Nine small RNA libraries and a mixed degradome library from flower petals at different developmental stages were constructed and sequenced. A total of 739 known miRNAs and 964 novel miRNAs were identified via small RNA sequencing, and 639 miRNAs were identified to cleave 2,816 target genes based on the degradome data. Additionally, 317 differentially expressed miRNAs among the various stages of flower development were identified, which regulated 2,134 differentially expressed target genes. These target genes were significantly enriched in the transcriptional regulation, phenylpropanoid biosynthesis, and plant hormone signal transduction pathways. Furthermore, integrated microRNAomic and transcriptomic analyses suggested that 98 miRNAs targeted several transcription factors, including MYBs (26), bHLHs (12), NACs (14), and SPLs (19), related to anthocyanin accumulation. In addition, 27 differentially expressed miRNAs might affect anthocyanin biosynthesis by regulating 23 targets involved in the hormone signal transduction pathway. The quantitative real-time PCR (qRT-PCR) analysis confirmed the expression changes of 21 miRNA-target pairs. Furthermore, the transient expression of candidate miRNAs was performed in the pink-flowered strawberry cultivar "Fenyun" at the bud stage. Introduction of FamiR156a, FamiR396e, and FamiR858_R-2 in the "Fenyun" increased flower color intensity, while transient expression of FamiR828a decreased flower color intensity. Overall, the present study uncovers the regulatory functions of microRNAs, including anthocyanin biosynthesis, hormone signaling, and regulation factors during flower development and coloration in pink-flowered strawberry. This work expands the knowledge of miRNAs affecting coloration in strawberry and provides rich resources for future functional studies.

Transcriptome and de novo analysis of Rosa xanthina f. spontanea in response to cold stress.

BACKGROUND: Rose is one of most popular ornamental plants worldwide and is of high economic value and great cultural importance. However, cold damage restricts its planting application in cold areas. To elucidate the metabolic response of rose under low temperature stress, we conducted transcriptome and de novo analysis of Rosa xanthina f. spontanea. RESULTS: A total of 124,106 unigenes from 9 libraries were generated by de novo assembly, with N50 length was 1470 bp, under 4 °C and - 20 °C stress (23 °C was used as a control). Functional annotation and prediction analyses identified 55,084 unigenes, and 67.72% of these unigenes had significant similarity (BLAST, E ≤ 10- 5) to those in the public databases. A total of 3031 genes were upregulated and 3891 were downregulated at 4 °C compared with 23 °C, and 867 genes were upregulated and 1763 were downregulated at - 20 °C compared with 23 °C. A total of 468 common DEGs were detected under cold stress, and the matched DEGs were involved in three functional categories: biological process (58.45%), cellular component (11.27%) and molecular function (30.28%). Based on KEGG functional annotations, four pathways were significantly enriched: metabolic pathway, response to plant pathogen interaction (32 genes); starch and sucrose metabolism (21 genes); circadian rhythm plant (8 genes); and photosynthesis antenna proteins (7 genes). CONCLUSIONS: Our study is the first to report the response to cold stress at the transcriptome level in R. xanthina f. spontanea. The results can help to elucidate the molecular mechanism of cold resistance in rose and provide new insights and candidate genes for genetically enhancing cold stress tolerance.

Evolutionary history and pan-genome dynamics of strawberry (Fragaria spp.).

Strawberry (Fragaria spp.) has emerged as a model system for various fundamental and applied research in recent years. In total, the genomes of five different species have been sequenced over the past 10 y. Here, we report chromosome-scale reference genomes for five strawberry species, including three newly sequenced species' genomes, and genome resequencing data for 128 additional accessions to estimate the genetic diversity, structure, and demographic history of key Fragaria species. Our analyses obtained fully resolved and strongly supported phylogenies and divergence times for most diploid strawberry species. These analyses also uncovered a new diploid species (Fragaria emeiensis Jia J. Lei). Finally, we constructed a pan-genome for Fragaria and examined the evolutionary dynamics of gene families. Notably, we identified multiple independent single base mutations of the MYB10 gene associated with white pigmented fruit shared by different strawberry species. These reference genomes and datasets, combined with our phylogenetic estimates, should serve as a powerful comparative genomic platform and resource for future studies in strawberry.

Genome-Wide Analysis of the R2R3-MYB Gene Family in Fragaria × ananassa and Its Function Identification During Anthocyanins Biosynthesis in Pink-Flowered Strawberry.

The strawberry (Fragaria × ananassa) is an economically important fruit throughout the world. The large R2R3-MYB gene family participates in a variety of plant functions, including anthocyanin biosynthesis. The present study is the first genome-wide analysis of the MYB gene family in the octoploid strawberry and describes the identification and characterization of the family members using the recently sequenced F. × ananassa genome. Specifically, we aimed to identify the key MYBs involved in petal coloration in the pink-flowered strawberry, which increases its ornamental value. A comprehensive, genome-wide analysis of F. × ananassa R2R3-FaMYBs was performed, investigating gene structures, phylogenic relationships, promoter regions, chromosomal locations, and collinearity. A total of 393 R2R3-FaMYB genes were identified in the F. × ananassa genome and divided into 36 subgroups based on phylogenetic analysis. Most genes with similar functions in the same subgroup exhibited similar exon-intron structures and motif compositions. These R2R3-FaMYBs were unevenly distributed over 28 chromosomes. The expansion of the R2R3-FaMYB gene family in the F. × ananassa genome was found to be caused mainly by segmental duplication. The Ka/Ks analysis indicated that duplicated R2R3-FaMYBs mostly experienced purifying selection and showed limited functional divergence after the duplication events. To elucidate which R2R3-FaMYB genes were associated with anthocyanin biosynthesis in the petals of the pink-flowered strawberry, we compared transcriptional changes in different flower developmental stages using RNA-seq. There were 131 differentially expressed R2R3-FaMYB genes identified in the petals, of which three genes, FaMYB28, FaMYB54, and FaMYB576, appeared likely, based on the phylogenetic analysis, to regulate anthocyanin biosynthesis. The qRT-PCR showed that these three genes were more highly expressed in petals than in other tissues (fruit, leaf, petiole and stolon) and their expressions were higher in red compared to pink and white petals. These results facilitate the clarification on the roles of the R2R3-FaMYB genes in petal coloration in the pink-flowered strawberry. This work provides useful information for further functional analysis on the R2R3-FaMYB gene family in F. × ananassa.

The long noncoding RNA FRILAIR regulates strawberry fruit ripening by functioning as a noncanonical target mimic.

Long noncoding RNAs (lncRNAs) are emerging as important regulators in plant development, but few of them have been functionally characterized in fruit ripening. Here, we have identified 25,613 lncRNAs from strawberry ripening fruits based on RNA-seq data from poly(A)-depleted libraries and rRNA-depleted libraries, most of which exhibited distinct temporal expression patterns. A novel lncRNA, FRILAIR harbours the miR397 binding site that is highly conserved in diverse strawberry species. FRILAIR overexpression promoted fruit maturation in the Falandi strawberry, which was consistent with the finding from knocking down miR397, which can guide the mRNA cleavage of both FRILAIR and LAC11a (encoding a putative laccase-11-like protein). Moreover, LAC11a mRNA levels were increased in both FRILAIR overexpressing and miR397 knockdown fruits, and accelerated fruit maturation was also found in LAC11a overexpressing fruits. Overall, our study demonstrates that FRILAIR can act as a noncanonical target mimic of miR397 to modulate the expression of LAC11a in the strawberry fruit ripening process.

Complete chloroplast genome sequencing of ten wild Fragaria species in China provides evidence for phylogenetic evolution of Fragaria.

Complete chloroplast genomes of ten wild Fragaria species native to China were sequenced. Phylogenetic analysis clustered Fragaria species into two clades: The south clade (F. iinumae, F. chinensis, F. pentaphylla, F. nilgerrensis, F. daltoniana, F. corymbosa, F. moupinensis, F. tibetica, F. nipponica, F. gracilis, and F. nubicola and north clade (F. viridis, F. orientalis, F. moschata, F. mandshurica, F. vesca, F. chiloensis, F. virginiana, and F. × ananassa), while F. iinumae is the oldest extant species. Molecular clock analysis suggested present Fragaria species share a common ancestor 3.57 million years ago (Ma), F. moschata and octoploid species evolve 0.89 and 0.97 Ma, respectively, but F. moschata be not directly involved in current octoploid species formation. Drastic global temperature change since the Palaeocene-Eocene, approx. 55 Ma, especially during uplifting of the Qinghai-Tibet plateau and quaternary glaciation may have driven the formation of Fragaria, separation of two groups and polyploidization.

Rapid Quantitation of 25-Hydroxyvitamin D2 and D3 in Human Serum Using Liquid Chromatography/Drift Tube Ion Mobility-Mass Spectrometry.

Ion mobility was integrated with liquid chromatography/high resolution mass spectrometry (LC/IM-HRMS) to quantify 25-hydroxyvitamin D (25OHD) in human serum. It has previously been shown that 25OHD adopts two gas-phase conformations which are resolved using ion mobility; in contrast, the inactive epimer, 3-epi-25-hydroxyvitamin D (epi25OHD), only adopts one. Interference from epi25OHD was eliminated by filtering the chromatogram to retain the drift time that corresponds to the unique gas-phase conformation of 25OHD. Although ion mobility separates the epimers, some chromatography is required to separate compounds which interfere with ionization or fall at the same nominal m/z. Standards were prepared in 4% albumin solutions and compared against commercial serum quality controls. Standards and quality controls were analyzed and validated using a 2 min LC/IM-MS method. 25-Hydroxyvitamin D3 and D2 were quantified over the range between 2 and 500 ng/mL with bias and precision within 15%. When epi25OHD was spiked into quality control samples, no significant bias was introduced, and analysis of 30 patient samples shows good agreement between this LC/IM-MS and traditional LC/MS/MS methods. This work shows that ion mobility can be incorporated with liquid chromatography and mass spectrometry for rapid quantitation of 25OHD in human serum.

Study on cyanidin metabolism in petals of pink-flowered strawberry based on transcriptome sequencing and metabolite analysis.

BACKGROUND: Pink-flowered strawberry is a promising new ornamental flower derived from intergeneric hybridization (Fragaria × Potentilla) with bright color, a prolonged flowering period and edible fruits. Its flower color ranges from light pink to red. Pigment compounds accumulated in its fruits were the same as in cultivated strawberry fruits, but different from that in its flowers. However, the transcriptional events underlying the anthocyanin biosynthetic pathway have not been fully characterized in petal coloration. To gain insights into the regulatory networks related to anthocyanin biosynthesis and identify the key genes, we performed an integrated analysis of the transcriptome and metabolome in petals of pink-flowered strawberry. RESULTS: The main pigments of red and dark pink petals were anthocyanins, among which cyanidins were the main compound. There were no anthocyanins detected in the white-flowered hybrids. A total of 50,285 non-redundant unigenes were obtained from the transcriptome databases involved in red petals of pink-flowered strawberry cultivar Sijihong at three development stages. Amongst the unigenes found to show significant differential expression, 57 were associated with anthocyanin or other flavonoid biosynthesis, in which they were regulated by 241 differentially expressed members of transcription factor families, such as 40 MYBs, 47 bHLHs, and 41 NACs. Based on a comprehensive analysis relating pigment compounds to gene expression profiles, the mechanism of flower coloration was examined in pink-flowered strawberry. A new hypothesis was proposed to explain the lack of color phenotype of the white-flowered strawberry hybrids based on the transcriptome analysis. The expression patterns of FpDFR and FpANS genes corresponded to the accumulation patterns of cyanidin contents in pink-flowered strawberry hybrids with different shades of pink. Moreover, FpANS, FpBZ1 and FpUGT75C1 genes were the major factors that led to the absence of anthocyanins in the white petals of pink-flowered strawberry hybrids. Meanwhile, the competitive effect of FpFLS and FpDFR genes might further inhibit anthocyanin synthesis. CONCLUSIONS: The data presented herein are important for understanding the molecular mechanisms underlying the petal pigmentation and will be powerful for integrating novel potential target genes to breed valuable pink-flowered strawberry cultivars.

De Novo Assembled Transcriptome Analysis and Identification of Genic SSR Markers in Red-Flowered Strawberry.

Red-flowered strawberry is a new ornamental flower derived from intergeneric hybridization (Fragaria × Potentilla). To date, few molecular markers have been reported for this plant. RNA sequencing provides a relatively fast and low-cost approach for large-scale detection of simple sequence repeats (SSRs). In the present study, we profiled the transcriptome of red-flowered strawberry by Illumina HiSeq 2500 to identify SSRs related to petal color. Based on 2 million clean reads of red and white flowers from red-flowered strawberry hybrids, we assembled 91,835 unigenes with an average length of 717 bp. After functional annotation and prediction, there were 47,058 unigenes; of these, 26,861 had a gene ontology annotation, with 14,264 SSR loci. Mononucleotide SSRs were the predominant repeat type (47.20%, n = 6724), followed by di- (32.50%, n = 4641), tri- (19.10%, n = 2729), tetra- (0.90%, n = 132), hexa- (0.2%, n = 21), and penta- (0.10%, n = 16) nucleotide repeats. The most frequent di-, tri-, and tetra-nucleotide repeats were AG/CT, AAG/CTT, and AAAG/CTTT, respectively. PCR amplification with 105 SSR primer pairs yielded four bands specific to red flowers, namely UgRFsr57622, UgRFsr94149, UgRFsr40142, and UgRFsr54608; corresponding 4 trait-specific markers were found to co-segregate with white and red flower color in hybrid population, demonstrating that the genic SSR marker is useful to discriminate between white and red flowers in strawberry. Markers to discriminate flower color in red-flowered strawberry will be useful for early selection of progeny and for breeding management.

Matrine Mediates Inflammatory Response via Gut Microbiota in TNBS-Induced Murine Colitis.

This study mainly investigated the effect of matrine on TNBS-induced intestinal inflammation in mice. TNBS treatment caused colonic injury and gut inflammation. Matrine (1, 5, and 10 mg/kg) treatment alleviated colonic injury and gut inflammation via reducing bleeding and diarrhea and downregulating cytokines expression (IL-1ß and TNF-α). Meanwhile, serum immunoglobulin G (IgG) was markedly reduced in TNBS treated mice, while 5 and 10 mg/kg matrine alleviated IgG reduction. Fecal microbiota was tested using 16S sequencing and the results showed that TNBS caused gut microbiota dysbiosis, while matrine treatment markedly improved gut microbiota communities (i.e., Bacilli and Mollicutes). Functional analysis showed that cell motility, nucleotide metabolism, and replication and repair were markedly altered in the TNBS group, while matrine treatment significantly affected cell growth and death, membrane transport, nucleotide metabolism, and replication and repair. In conclusion, matrine may serve as a protective mechanism in TNBS-induced colonic inflammation and the beneficial effect may be associated with gut microbiota.

Influence of Experimental Conditions on the Ratio of 25-Hydroxyvitamin D3 Conformers for Validating a Liquid Chromatography/Ion Mobility-Mass Spectrometry Method for Routine Quantitation.

The utility of adding ion mobility (IM) to liquid chromatography/mass spectrometry (LC/MS) for quantitation of 25-hydroxyvitamin D3 (25OHD3) is evaluated. Sodiated 25OHD3 ions adopt both closed and open conformations, whereas the stereoisomer 3-epi-25-hydroxyvitamin D (epi25OHD), when sodiated, adopts only a closed gas-phase conformation. The unique open conformation for sodiated 25OHD3 permits unambiguous quantitation. Nonetheless, the ratio of open versus closed gas-phase conformations for sodiated 25OHD3 can vary with instrumental conditions; conversion from the open to the closed conformer is attributed to radio frequency (rf) heating within the ion accumulation trap. Ion heating becomes significant when space charge in the trap forces ions to larger radii where the rf field is higher. To avoid biasing quantitation, an isotopically labeled internal standard must be used to account for changes in conformer ratio. Thirty-three serum extract samples were evaluated by LC/IM-MS and were found to not be biased by changes in ion conformer ratios, permitting reliable quantitation of 25OHD3 without interference from the epimer.

Characterization of anthocyanins in the hybrid progenies derived from Iris dichotoma and I. domestica by HPLC-DAD-ESI/MS analysis.

Iris dichotoma with different flower colors and I. domestica are beardless wild irises belonging to the family Iridaceae that bloom in the summer and have long flowering periods. In this study, we collected three accessions of I. dichotoma with violet, yellow, and white flowers, respectively, in China, and crossed them with I. domestica individuals. The flower color of the hybrids derived from these crosses was categorized into eight groups: violet, purple, brown, orange, red, pink, yellow, and white. From this population, 45 individuals were selected for analysis, and their fully expanded inner and outer perianths were harvested for extraction of anthocyanins. Using high-performance liquid chromatography-mass spectrometry (HPLC-MS) analysis, 29 anthocyanins were identified by comparing MS and UV-visible spectra and elution order based on published data and guidelines. The 29 anthocyanins were classified into six groups: non-acylated glycosides (3RG, 3RG5G), acetylglycosides (3acRG5G), p-coumaroylglycosides (3pCRG, 3pCRG5G), caffeoylglycosides (3CRG5G), feruloylglycosides (3feRG, 3feRG5G), and acetyl-(p-coumaroyl) glycosides (3ac-pCRG5G). Acylated anthocyanin contents were considerably higher than non-acylated anthocyanin contents in the individuals evaluated, regardless of flower color, except in the yellow-flowered I. dichotoma and its yellow-flowered progeny. We found ten anthocyanins derived from pelargonidin, including pelargonidin 3-O-(caffeoyl)rutinoside-5-O-glucoside (Pg3CRG5G), pelargonidin 3-O-(feruloyl)rutinoside-5-O-glucoside (Pg3feRG5G), and pelargonidin 3-O-(feruloyl)rutinoside (Pg3feRG), that have not yet been reported in other Iris species. Moreover, delphinidin 3-O-(feruloyl) rutinoside-5-O-glucoside (Dp3feRG5G), and delphinidin 3-O-(feruloyl)rutinoside (Dp3feRG) were also characterized for the first time in Iris. Two to five major anthocyanins were detected in the petals of the violet and purple groups, whereas those of the brown group contained three to six anthocyanins. Petals of the orange, red, and pink groups contained two to four major anthocyanins. The total variance explained by the first two principal components (PC) in a PC analysis of the intensities of the compound peaks detected in individuals representing the eight color groups was 44.7%. Individuals were classified into 10 types (A-I to A-X) according to their anthocyanidin composition in the flower petals. Knowledge of the types of anthocyanins determining flower petal color will aid comprehensive understanding of flower color characteristics in hybrid progeny and will contribute to the development of breeding Iris cultivars of diverse flower colors.

Comparative Transcriptomics of Strawberries (Fragaria spp.) Provides Insights into Evolutionary Patterns.

Multiple closely related species with genomic sequences provide an ideal system for studies on comparative and evolutionary genomics, as well as the mechanism of speciation. The whole genome sequences of six strawberry species (Fragaria spp.) have been released, which provide one of the richest genomic resources of any plant genus. In this study, we first generated seven transcriptome sequences of Fragaria species de novo, with a total of 48,557-82,537 unigenes per species. Combined with 13 other species genomes in Rosales, we reconstructed a phylogenetic tree at the genomic level. The phylogenic tree shows that Fragaria closed grouped with Rubus and the Fragaria clade is divided into three subclades. East Asian species appeared in every subclade, suggesting that the genus originated in this area at ∼7.99 Mya. Four species found in mountains of Southwest China originated at ∼3.98 Mya, suggesting that rapid speciation occurred to adapt to changing environments following the uplift of the Qinghai-Tibet Plateau. Moreover, we identified 510 very significantly positively selected genes in the cultivated species F. × ananassa genome. This set of genes was enriched in functions related to specific agronomic traits, such as carbon metabolism and plant hormone signal transduction processes, which are directly related to fruit quality and flavor. These findings illustrate comprehensive evolutionary patterns in Fragaria and the genetic basis of fruit domestication of cultivated strawberry at the genomic/transcriptomic level.