Gene expression atlas of fruit ripening and transcriptome assembly from RNA-seq data in octoploid strawberry (Fragaria × ananassa)

YABBY transcription factor family in the octoploid Fragaria × ananassa and five diploid Fragaria species

YABBY genes encode specific TFs of seed plants involved in development and formation of leaves, flowers, and fruit. In the present work, genome-wide and expression analyses of the YABBY gene family were performed in six species of the Fragaria genus: Fragaria × ananassa, F. daltoniana, F. nilgerrensis, F. pentaphylla, F. viridis, and F. vesca. The chromosomal location, synteny pattern, gene structure, and phylogenetic analyses were carried out. By combining RNA-seq data and RT-qPCR analysis we explored specific expression of YABBYs in F. × ananassa and F. vesca. We also analysed the promoter regions of FaYABBYs and performed MeJA application to F. × ananassa fruit to observe effects on gene expression. We identified and characterized 25 YABBY genes in F. × ananassa and six in each of the other five species, which belong to FIL/YAB3 (YABBY1), YAB2 (YABBY2), YAB5 (YABBY5), CRC, and INO clades previously described. Division of the YABBY1 clade into YABBY1.1 and YABBY1.2 subclades is reported. We observed differential expression according to tissue, where some FaYABBYs are expressed mainly in leaves and flowers and to a minor extent during fruit development of F. × ananassa. Specifically, the FaINO genes contain jasmonate-responsive cis-acting elements in their promoters which may be functional since FaINOs are upregulated in F. × ananassa fruit under MeJA treatment. This study suggests that YABBY TFs play an important role in the development- and environment-associated responses of the Fragaria genus.

Genome-Wide Analysis of Polygalacturonase Gene Family Reveals Its Role in Strawberry Softening

Fruit softening is a prominent attribute governing both longevity on shelves and commercial worth. Polygalacturonase (PG) plays a major role in strawberry fruit softening. However, the PG gene family in strawberry has not been comprehensively analyzed. In this study, 75 FaPG genes were identified in the octoploid strawberry genome, which were classified into three groups according to phylogenetic analysis. Subcellular localization prediction indicated that FaPGs are mostly localized to the plasma membrane, cytoplasm, and chloroplasts. Moreover, the expression of FaPGs during strawberry development and ripening of 'Benihoppe' and its softer mutant was estimated. The results showed that among all 75 FaPGs, most genes exhibited low expression across developmental stages, while two group c members (FxaC_21g15770 and FxaC_20g05360) and one group b member, FxaC_19g05040, displayed relatively higher and gradual increases in their expression trends during strawberry ripening and softening. FxaC_21g15770 was selected for subsequent silencing to validate its role in strawberry softening due to the fact that it exhibited the highest and most changed expression level across different developmental stages in 'Benihoppe' and its mutant. Silencing FxaC_21g15770 could significantly improve strawberry fruit firmness without affecting fruit color, soluble solids, cellulose, and hemicellulose. Conversely, silencing FxaC_21g15770 could significantly suppress the expression of other genes related to pectin degradation such as FaPG-like, FaPL, FaPME, FaCX, FaCel, FaGlu, FaXET, and FaEG. These findings provide basic information on the FaPG gene family for further functional research and indicate that FxaC_21g15770 plays a vital role in strawberry fruit softening.

Advances in genomics and genome editing for improving strawberry (Fragaria ×ananassa)

The cultivated strawberry, Fragaria ×ananassa, is a recently domesticated fruit species of economic interest worldwide. As such, there is significant interest in continuous varietal improvement. Genomics-assisted improvement, including the use of DNA markers and genomic selection have facilitated significant improvements of numerous key traits during strawberry breeding. CRISPR/Cas-mediated genome editing allows targeted mutations and precision nucleotide substitutions in the target genome, revolutionizing functional genomics and crop improvement. Genome editing is beginning to gain traction in the more challenging polyploid crops, including allo-octoploid strawberry. The release of high-quality reference genomes and comprehensive subgenome-specific genotyping and gene expression profiling data in octoploid strawberry will lead to a surge in trait discovery and modification by using CRISPR/Cas. Genome editing has already been successfully applied for modification of several strawberry genes, including anthocyanin content, fruit firmness and tolerance to post-harvest disease. However, reports on many other important breeding characteristics associated with fruit quality and production are still lacking, indicating a need for streamlined genome editing approaches and tools in Fragaria ×ananassa. In this review, we present an overview of the latest advancements in knowledge and breeding efforts involving CRISPR/Cas genome editing for the enhancement of strawberry varieties. Furthermore, we explore potential applications of this technology for improving other Rosaceous plant species.

Non-climacteric fruit development and ripening regulation: 'the phytohormones show'

Fruit ripening involves numerous physiological, structural, and metabolic changes that result in the formation of edible fruits. This process is controlled at different molecular levels, with essential roles for phytohormones, transcription factors, and epigenetic modifications. Fleshy fruits are classified as either climacteric or non-climacteric species. Climacteric fruits are characterized by a burst in respiration and ethylene production at the onset of ripening, while regulation of non-climacteric fruit ripening has been commonly attributed to abscisic acid (ABA). However, there is controversy as to whether mechanisms regulating fruit ripening are shared between non-climacteric species, and to what extent other hormones contribute alongside ABA. In this review, we summarize classic and recent studies on the accumulation profile and role of ABA and other important hormones in the regulation of non-climacteric fruit development and ripening, as well as their crosstalk, paying special attention to the two main non-climacteric plant models, strawberry and grape. We highlight both the common and different roles of these regulators in these two crops, and discuss the importance of the transcriptional and environmental regulation of fruit ripening, as well as the need to optimize genetic transformation methodologies to facilitate gene functional analyses.

Molecular bases of strawberry fruit quality traits: Advances, challenges, and opportunities

The strawberry is one of the world's most popular fruits, providing humans with vitamins, fibers, and antioxidants. Cultivated strawberry (Fragaria × ananassa) is an allo-octoploid and highly heterozygous, making it a challenge for breeding, quantitative trait locus (QTL) mapping, and gene discovery. Some wild strawberry relatives, such as Fragaria vesca, have diploid genomes and are becoming laboratory models for the cultivated strawberry. Recent advances in genome sequencing and CRISPR-mediated genome editing have greatly improved the understanding of various aspects of strawberry growth and development in both cultivated and wild strawberries. This review focuses on fruit quality traits that are most relevant to the consumers, including fruit aroma, sweetness, color, firmness, and shape. Recently available phased-haplotype genomes, single nucleotide polymorphism (SNP) arrays, extensive fruit transcriptomes, and other big data have made it possible to locate key genomic regions or pinpoint specific genes that underlie volatile synthesis, anthocyanin accumulation for fruit color, and sweetness intensity or perception. These new advances will greatly facilitate marker-assisted breeding, the introgression of missing genes into modern varieties, and precise genome editing of selected genes and pathways. Strawberries are poised to benefit from these recent advances, providing consumers with fruit that is tastier, longer-lasting, healthier, and more beautiful.

Storage of halved strawberry fruits affects aroma, phytochemical content and gene expression, and is affected by pre-harvest factors

Introduction

Strawberry fruit are highly valued for their aroma which develops during ripening. However, they have a short shelf-life. Low temperature storage is routinely used to extend shelf-life for transport and storage in the supply chain, however cold storage can also affect fruit aroma. Some fruit continue to ripen during chilled storage; however, strawberries are a non-climacteric fruit and hence ripening postharvest is limited. Although most strawberry fruit is sold whole, halved fruit is also used in ready to eat fresh fruit salads which are of increasing consumer demand and pose additional challenges to fresh fruit storage.

Methods

To better understand the effects of cold storage, volatilomic and transcriptomic analyses were applied to halved Fragaria x ananassa cv. Elsanta fruit stored at 4 or 8°C for up to 12 days over two growing seasons.

Results and discussion

The volatile organic compound (VOC) profile differed between 4 or 8°C on most days of storage. Major differences were detected between the two different years of harvest indicating that aroma change at harvest and during storage is highly dependent on environmental factors during growth. The major component of the aroma profile in both years was esters. Over 3000 genes changed in expression over 5 days of storage at 8°C in transcriptome analysis. Overall, phenylpropanoid metabolism, which may also affect VOCs, and starch metabolism were the most significantly affected pathways. Genes involved in autophagy were also differentially expressed. Expression of genes from 43 different transcription factor (TF) families changed in expression: mostly they were down-regulated but NAC and WRKY family genes were mainly up-regulated. Given the high ester representation amongst VOCs, the down-regulation of an alcohol acyl transferase (AAT) during storage is significant. A total of 113 differentially expressed genes were co-regulated with the AAT gene, including seven TFs. These may be potential AAT regulators.

A molecular phenology scale of grape berry development

Fruit growth and development consist of a continuous succession of physical, biochemical, and physiological changes driven by a genetic program that dynamically responds to environmental cues. Establishing recognizable stages over the whole fruit lifetime represents a fundamental requirement for research and fruit crop cultivation. This is especially relevant in perennial crops like grapevine (Vitis vinifera L.) to scale the development of its fruit across genotypes and growing conditions. In this work, molecular-based information from several grape berry transcriptomic datasets was exploited to build a molecular phenology scale (MPhS) and to map the ontogenic development of the fruit. The proposed statistical pipeline consisted of an unsupervised learning procedure yielding an innovative combination of semiparametric, smoothing, and dimensionality reduction tools. The transcriptomic distance between fruit samples was precisely quantified by means of the MPhS that also enabled to highlight the complex dynamics of the transcriptional program over berry development through the calculation of the rate of variation of MPhS stages by time. The MPhS allowed the alignment of time-series fruit samples proving to be a complementary method for mapping the progression of grape berry development with higher detail compared to classic time- or phenotype-based approaches.

QTL analysis for ascorbic acid content in strawberry fruit reveals a complex genetic architecture and association with GDP-L-galactose phosphorylase

Strawberry (Fragaria × ananassa) fruits are an excellent source of L-ascorbic acid (AsA), a powerful antioxidant for plants and humans. Identifying the genetic components underlying AsA accumulation is crucial for enhancing strawberry nutritional quality. Here, we unravel the genetic architecture of AsA accumulation using an F₁ population derived from parental lines 'Candonga' and 'Senga Sengana', adapted to distinct Southern and Northern European areas. To account for environmental effects, the F₁ and parental lines were grown and phenotyped in five locations across Europe (France, Germany, Italy, Poland and Spain). Fruit AsA content displayed normal distribution typical of quantitative traits and ranged five-fold, with significant differences among genotypes and environments. AsA content in each country and the average in all of them was used in combination with 6,974 markers for quantitative trait locus (QTL) analysis. Environmentally stable QTLs for AsA content were detected in linkage group (LG) 3A, LG 5A, LG 5B, LG 6B and LG 7C. Candidate genes were identified within stable QTL intervals and expression analysis in lines with contrasting AsA content suggested that GDP-L-Galactose Phosphorylase FaGGP(3A), and the chloroplast-located AsA transporter gene FaPHT4;4(7C) might be the underlying genetic factors for QTLs on LG 3A and 7C, respectively. We show that recessive alleles of FaGGP(3A) inherited from both parental lines increase fruit AsA content. Furthermore, expression of FaGGP(3A) was two-fold higher in lines with high AsA. Markers here identified represent a useful resource for efficient selection of new strawberry cultivars with increased AsA content.

Comprehensive profiling of endogenous phytohormones and expression analysis of 1-aminocyclopropane-1-carboxylic acid synthase gene family during fruit development and ripening in octoploid strawberry (Fragaria× ananassa)

The non-climacteric octoploid strawberry (Fragaria × ananassa Duchesne ex Rozier) was used as a model to study its regulation during fruit ripening. High performance liquid chromatography electrospray tandem-mass spectrometry (HPLC-ESI-MS/MS) was employed to profile 28 different endogenous phytohormones in strawberry. These include auxins, cytokinins (CKs), abscisic acid (ABA), ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC), jasmonates, and phenolic compounds salicylic acid (SA), benzoic acid (BzA) and phenylacetic acid (PAA) together with their various metabolic forms that have remained largely unexplored thus far. ABA, ACC and CK N⁶-(Δ²-isopentenyl)adenine (iP) were found to be associated with ripening while ABA catabolites 9-hydroxy-ABA and phaseic acid mimicked the pattern of climacteric decline at the turning phase of strawberry ripening. The content of other CK forms except iP decreased as fruit ripened, as also that of auxins indole-3-acetic acid (IAA) and oxo-IAA, and of jasmonates. Data presented here also suggest that both the transition and progression of strawberry fruit ripening are associated with N⁶-(Δ²-isopentenyl)adenosine-5'-monophosphate (iPRMP) → N⁶-(Δ²-isopentenyl)adenosine (iPR) → iP as the preferred CK metabolic pathway. In contrast, the ethylene precursor ACC was present at higher levels, with its abundance increasing from the onset of ripening to the red ripe stage. Further investigation of ripening-specific ACC accumulation revealed the presence of a large ACC synthase (ACS) encoding gene family in octoploid strawberry that was previously unknown. Seventeen ACS genes were found differentially expressed in fruit tissues, while six of them showed induced expression during strawberry fruit ripening. These data suggest a possible role(s) of ACC, ABA, and iP in strawberry fruit ripening. These data add new dimension to the existing knowledge of the interplay of different endogenous phytohormones in octoploid strawberry, paving the way for further investigation of their individual role(s) in fruit ripening.

Transcriptomic analysis of seed development in Paysonia auriculata (Brassicaceae) identifies genes involved in hydroxy fatty acid biosynthesis

Paysonia auriculata (Brassicaceae) produces multiple hydroxy fatty acids as major components of the seed oil. We tracked the changes in seed oil composition and gene expression during development, starting 14 days after flowers had been pollinated. Seed oil changes showed initially higher levels of saturated and unsaturated fatty acids (FAs) but little accumulation of hydroxy fatty acids (HFAs). Starting 21 days after pollination (DAP) HFA content sharply increased, and reached almost 30% at 28 DAP. Total seed oil also increased from a low of approximately 2% at 14 DAP to a high of approximately 20% by 42 DAP. We identified almost all of the fatty acid synthesis and modification genes that are known from Arabidopsis, and, in addition, a strong candidate for the hydroxylase gene that mediates the hydroxylation of fatty acids to produce valuable hydroxy fatty acids (HFAs) in this species. The gene expression network revealed is very similar to that of the emerging oil crop, Physaria fendleri, in the sister genus to Paysonia. Phylogenetic analyses indicate the hydroxylase enzyme, FAH12, evolved only once in Paysonia and Physaria, and that the enzyme is closely related to FAD2 enzymes. Phylogenetic analyses of FAD2 and FAH12 in the Brassicaceae and outgroup genera suggest that the branch leading to the hydroxylase clade of Paysonia and Physaria is under relaxed selection, compared with the strong purifying selection found across the FAD2 lineages.

High-quality haplotype-resolved genome assembly of cultivated octoploid strawberry

Cultivated strawberry (Fragaria × ananassa), a perennial herb belonging to the family Rosaceae, is a complex octoploid with high heterozygosity at most loci. However, there is no research on the haplotype of the octoploid strawberry genome. Here we aimed to obtain a high-quality genome of the cultivated strawberry cultivar, "Yanli", using single molecule real-time sequencing and high-throughput chromosome conformation capture technology. The "Yanli" genome was 823 Mb in size, with a long terminal repeat assembly index of 14.99. The genome was phased into two haplotypes, Hap1 (825 Mb with contig N50 of 26.70 Mb) and Hap2 (808 Mb with contig N50 of 27.51 Mb). Using the combination of Hap1 and Hap2, we obtained for the first time a haplotype-resolved genome with 56 chromosomes for the cultivated octoploid strawberry. We identified a ~ 10 Mb inversion and translocation on chromosome 2-1. 104 957 and 102 356 protein-coding genes were annotated in Hap1 and Hap2, respectively. Analysis of the genes related to the anthocyanin biosynthesis pathway revealed the structural diversity and complexity in the expression of the alleles in the octoploid F. × ananassa genome. In summary, we obtained a high-quality haplotype-resolved genome assembly of F. × ananassa, which will provide the foundation for investigating gene function and evolution of the genome of cultivated octoploid strawberry.

A multi-omics framework reveals strawberry flavor genes and their regulatory elements

Flavor is essential to consumer preference of foods and is an increasing focus of plant breeding programs. In fruit crops, identifying genes underlying volatile organic compounds has great promise to accelerate flavor improvement, but polyploidy and heterozygosity in many species have slowed progress. Here we use octoploid cultivated strawberry to demonstrate how genomic heterozygosity, transcriptomic intricacy and fruit metabolomic diversity can be treated as strengths and leveraged to uncover fruit flavor genes and their regulatory elements. Multi-omics datasets were generated including an expression quantitative trait loci map with 196 diverse breeding lines, haplotype-phased genomes of a highly-flavored breeding selection, a genome-wide structural variant map using five haplotypes, and volatile genome-wide association study (GWAS) with > 300 individuals. Overlaying regulatory elements, structural variants and GWAS-linked allele-specific expression of numerous genes to variation in volatile compounds important to flavor. In one example, the functional role of anthranilate synthase alpha subunit 1 in methyl anthranilate biosynthesis was supported via fruit transient gene expression assays. These results demonstrate a framework for flavor gene discovery in fruit crops and a pathway to molecular breeding of cultivars with complex and desirable flavor.

Insights into transcription factors controlling strawberry fruit development and ripening

Fruit ripening is a highly regulated and complex process involving a series of physiological and biochemical changes aiming to maximize fruit organoleptic traits to attract herbivores, maximizing therefore seed dispersal. Furthermore, this process is of key importance for fruit quality and therefore consumer acceptance. In fleshy fruits, ripening involves an alteration in color, in the content of sugars, organic acids and secondary metabolites, such as volatile compounds, which influence flavor and aroma, and the remodeling of cell walls, resulting in the softening of the fruit. The mechanisms underlying these processes rely on the action of phytohormones, transcription factors and epigenetic modifications. Strawberry fruit is considered a model of non-climacteric species, as its ripening is mainly controlled by abscisic acid. Besides the role of phytohormones in the regulation of strawberry fruit ripening, a number of transcription factors have been identified as important regulators of these processes to date. In this review, we present a comprehensive overview of the current knowledge on the role of transcription factors in the regulation of strawberry fruit ripening, as well as in compiling candidate regulators that might play an important role but that have not been functionally studied to date.

RNAseq, transcriptome analysis and identification of DEGs involved in development and ripening of Fragaria chiloensis fruit

Fragaria chiloensis (Chilean strawberry) is a native species that produces fruit with an exotic pinkish color and a fruity aroma. It has a non-climacteric pattern of fruit ripening, and it is the mother of the commercial Fragaria x ananassa. The ripening of F. chiloensis fruit seems stimulated by ABA, and a complete set of genes participate in its softening, color, and aroma development. In addition, a set of transcription factors regulate the entire process, but few of them have been described. Over the last two decades, RNA-seq was used to identify genes at three fruit development/ripening stages, named C2 (unripe, large green) to C4 (full ripe), in whole fruit and fruit without achenes. A total of 204,754 contigs were assembled considering all samples, obtaining an N50 of 1.125 bp. Differentially expressed genes (DEGs) between two samples were identified, obtaining a total of 77,181 DEGs. Transcripts for genes involved in ABA biosynthesis present high and differential expression during the C2, C3, and C4 stages. Besides, contigs corresponding to ABA receptors, which interact with a regulatory network, are also differentially expressed. Genes associated with cell wall remodeling and those involved in flavonoid synthesis were also differentially expressed. An interaction network was built considering differentially expressed genes for the phenylpropanoid and flavonoid molecular pathways and having FcMYB1 as a transcription factor regulator. Identifying key genes could give an option to control the ripening of this non-climacteric fruit.

Genome-Wide Analysis of Ascorbic Acid Metabolism Related Genes in Fragaria × ananassa and Its Expression Pattern Analysis in Strawberry Fruits

Ascorbic acid (AsA) is an important antioxidant for scavenging reactive oxygen species and it is essential for human health. Strawberry (Fragaria × ananassa) fruits are rich in AsA. In recent years, strawberry has been regarded as a model for non-climacteric fruit ripening. However, in contrast to climacteric fruits, such as tomato, the regulatory mechanism of AsA accumulation in strawberry fruits remains largely unknown. In this study, we first identified 125 AsA metabolism-related genes from the cultivated strawberry "Camarosa" genome. The expression pattern analysis using an available RNA-seq data showed that the AsA biosynthetic-related genes in the D-mannose/L-galactose pathway were downregulated remarkably during fruit ripening which was opposite to the increasing AsA content in fruits. The D-galacturonate reductase gene (GalUR) in the D-Galacturonic acid pathway was extremely upregulated in strawberry receptacles during fruit ripening. The FaGalUR gene above belongs to the aldo-keto reductases (AKR) superfamily and has been proposed to participate in AsA biosynthesis in strawberry fruits. To explore whether there are other genes in the AKR superfamily involved in regulating AsA accumulation during strawberry fruit ripening, we further implemented a genome-wide analysis of the AKR superfamily using the octoploid strawberry genome. A total of 80 FaAKR genes were identified from the genome and divided into 20 subgroups based on phylogenetic analysis. These FaAKR genes were unevenly distributed on 23 chromosomes. Among them, nine genes showed increased expression in receptacles as the fruit ripened, and notably, FaAKR23 was the most dramatically upregulated FaAKR gene in receptacles. Compared with fruits at green stage, its expression level increased by 142-fold at red stage. The qRT-PCR results supported that the expression of FaAKR23 was increased significantly during fruit ripening. In particular, the FaAKR23 was the only FaAKR gene that was significantly upregulated by abscisic acid (ABA) and suppressed by nordihydroguaiaretic acid (NDGA, an ABA biosynthesis blocker), indicating FaAKR23 might play important roles in ABA-mediated strawberry fruit ripening. In a word, our study provides useful information on the AsA metabolism during strawberry fruit ripening and will help understand the mechanism of AsA accumulation in strawberry fruits.

A Comprehensive Study of the WRKY Transcription Factor Family in Strawberry

WRKY transcription factors play critical roles in plant growth and development or stress responses. Using up-to-date genomic data, a total of 64 and 257 WRKY genes have been identified in the diploid woodland strawberry, Fragaria vesca, and the more complex allo-octoploid commercial strawberry, Fragaria × ananassa cv. Camarosa, respectively. The completeness of the new genomes and annotations has enabled us to perform a more detailed evolutionary and functional study of the strawberry WRKY family members, particularly in the case of the cultivated hybrid, in which homoeologous and paralogous FaWRKY genes have been characterized. Analysis of the available expression profiles has revealed that many strawberry WRKY genes show preferential or tissue-specific expression. Furthermore, significant differential expression of several FaWRKY genes has been clearly detected in fruit receptacles and achenes during the ripening process and pathogen challenged, supporting a precise functional role of these strawberry genes in such processes. Further, an extensive analysis of predicted development, stress and hormone-responsive cis-acting elements in the strawberry WRKY family is shown. Our results provide a deeper and more comprehensive knowledge of the WRKY gene family in strawberry.

Transcriptional Basis for Haustorium Formation and Host Establishment in Hemiparasitic Psittacanthus schiedeanus Mistletoes

The mistletoe Psittacanthus schiedeanus, a keystone species in interaction networks between plants, pollinators, and seed dispersers, infects a wide range of native and non-native tree species of commercial interest. Here, using RNA-seq methodology we assembled the whole circularized quadripartite structure of P. schiedeanus chloroplast genome and described changes in the gene expression of the nuclear genomes across time of experimentally inoculated seeds. Of the 140,467 assembled and annotated uniGenes, 2,000 were identified as differentially expressed (DEGs) and were classified in six distinct clusters according to their expression profiles. DEGs were also classified in enriched functional categories related to synthesis, signaling, homoeostasis, and response to auxin and jasmonic acid. Since many orthologs are involved in lateral or adventitious root formation in other plant species, we propose that in P. schiedeanus (and perhaps in other rootless mistletoe species), these genes participate in haustorium formation by complex regulatory networks here described. Lastly, and according to the structural similarities of P. schiedeanus enzymes with those that are involved in host cell wall degradation in fungi, we suggest that a similar enzymatic arsenal is secreted extracellularly and used by mistletoes species to easily parasitize and break through tissues of the host.

Transcriptomic Complexity in Strawberry Fruit Development and Maturation Revealed by Nanopore Sequencing

The use of alternative transcription start or termination sites (aTSS or aTTS) as well as alternative splicing (AS) produce diverse transcript isoforms, playing indispensable roles in the plant development and environmental adaptations. Despite the advances in the finding of the genome-wide alternatively spliced genes in strawberry, it remains unexplored how AS responds to the developmental cues and what relevance do these outcomes have to the gene function. In this study, we have systematically investigated the transcriptome complexity using long-read Oxford Nanopore Technologies along the four successive developmental stages. The full-length cDNA sequencing results unraveled thousands of previously unexplored transcript isoforms raised from aTSS, aTTS, and AS. The relative contributions of these three processes to the complexity of strawberry fruit transcripts were compared. The aTSS and aTTS were more abundant than the AS. Differentially expressed transcripts unraveled the key transitional role of the white fruit stage. Isoform switches of transcripts from 757 genes were observed. They were associated with protein-coding potential change and domain gain or loss as the main consequences. Those genes with switched isoforms take part in the key processes of maturation in the late stages. A case study using yeast two hybrid analysis supported the functional divergence of the two isoforms of the B-box protein 22. Our results provided a new comprehensive overview of the dynamic transcriptomic landscape during strawberry fruit development and maturation.

Genome-Wide Analysis of MYB10 Transcription Factor in Fragaria and Identification of QTLs Associated with Fruit Color in Octoploid Strawberry

The genus Fragaria encompass fruits with diverse colors influenced by the distribution and accumulation of anthocyanin. Particularly, the fruit colors of strawberries with different ploidy levels are determined by expression and natural variations in the vital structural and regulatory genes involved in the anthocyanin pathway. Among the regulatory genes, MYB10 transcription factor is crucial for the expression of structural genes in the anthocyanin pathway. In the present study, we performed a genome wide investigation of MYB10 in the diploid and octoploid Fragaria species. Further, we identified seven quantitative trait loci (QTLs) associated with fruit color in octoploid strawberry. In addition, we predicted 20 candidate genes primarily influencing the fruit color based on the QTL results and transcriptome analysis of fruit skin and flesh tissues of light pink, red, and dark red strawberries. Moreover, the computational and transcriptome analysis of MYB10 in octoploid strawberry suggests that the difference in fruit colors could be predominantly influenced by the expression of MYB10 from the F. iinumae subgenome. The outcomes of the present endeavor will provide a platform for the understanding and tailoring of anthocyanin pathway in strawberry for the production of fruits with aesthetic colors.

An Evolutionary Analysis of B-Box Transcription Factors in Strawberry Reveals the Role of FaBBx28c1 in the Regulation of Flowering Time

Flowering connects vegetative and generative developmental phases and plays a significant role in strawberry production. The mechanisms that regulate strawberry flowering time are unclear. B-box transcription factors (BBXs) play important roles in the flowering time regulation of plants. Nevertheless, BBXs in octoploid cultivated strawberry (Fragaria ananassa) and their functions in flowering time regulation have not been identified. Here, we identified 51 FaBBXs from cultivated strawberry and 16 FvBBXs from diploid wild strawberry (Fragaria vesca), which were classified into five groups according to phylogenetic analysis. Further evolutionary analysis showed that whole-genome duplication or segmental duplication is a crucial factor that leads to the expansion of the BBX gene family in two strawberry species. Moreover, some loss and acquisition events of FaBBX genes were identified in the genome of cultivated strawberry that could have affected traits of agronomic interest, such as fruit quality. The promoters of FaBBX genes showed an enrichment in light-responsive, cis-regulatory elements, with 16 of these genes showing changes in their transcriptional activity in response to blue light treatment. On the other hand, FaBBX28c1, whose transcriptional activity is reduced in response to blue light, showed a delay in flowering time in Arabidopsis transgenic lines, suggesting its role in the regulation of flowering time in cultivated strawberry. Our results provide new evolutionary insight into the BBX gene family in cultivated strawberry and clues regarding their function in flowering time regulation in plants.

Systems Biology Applied to the Study of Papaya Fruit Ripening: The Influence of Ethylene on Pulp Softening

Papaya is a fleshy fruit that undergoes fast ethylene-induced modifications. The fruit becomes edible, but the fast pulp softening is the main factor that limits the post-harvest period. Papaya fast pulp softening occurs due to cell wall disassembling coordinated by ethylene triggering that massively expresses pectinases. In this work, RNA-seq analysis of ethylene-treated and non-treated papayas enabled a wide transcriptome overview that indicated the role of ethylene during ripening at the gene expression level. Several families of transcription factors (AP2/ERF, NAC, and MADS-box) were differentially expressed. ACO, ACS, and SAM-Mtase genes were upregulated, indicating a high rate of ethylene biosynthesis after ethylene treatment. The correlation among gene expression and physiological data demonstrated ethylene treatment can indeed simulate ripening, and regulation of changes in fruit color, aroma, and flavor could be attributed to the coordinated expression of several related genes. Especially about pulp firmness, the identification of 157 expressed genes related to cell wall metabolism demonstrated that pulp softening is accomplished by a coordinated action of several different cell wall-related enzymes. The mechanism is different from other commercially important fruits, such as strawberry, tomato, kiwifruit, and apple. The observed behavior of this new transcriptomic data confirms ethylene triggering is the main event that elicits fast pulp softening in papayas.

Autophagy Is Required for Strawberry Fruit Ripening

Autophagy is a catabolic and recycling pathway that maintains cellular homeostasis under normal growth and stress conditions. Two major types of autophagy, microautophagy and macroautophagy, have been described in plants. During macroautophagy, cellular content is engulfed by a double-membrane vesicle called autophagosome. This vesicle fuses its outer membrane with the tonoplast and releases the content into the vacuole for degradation. During certain developmental processes, autophagy is enhanced by induction of several autophagy-related genes (ATG genes). Autophagy in crop development has been studied in relation to leaf senescence, seed and reproductive development, and vascular formation. However, its role in fruit ripening has only been partially addressed. Strawberry is an important berry crop, representative of non-climacteric fruit. We have analyzed the occurrence of autophagy in developing and ripening fruits of the cultivated strawberry. Our data show that most ATG genes are conserved in the genome of the cultivated strawberry Fragaria x ananassa and they are differentially expressed along the ripening of the fruit receptacle. ATG8-lipidation analysis proves the presence of two autophagic waves during ripening. In addition, we have confirmed the presence of autophagy at the cellular level by the identification of autophagy-related structures at different stages of the strawberry ripening. Finally, we show that blocking autophagy either biochemically or genetically dramatically affects strawberry growth and ripening. Our data support that autophagy is an active and essential process with different implications during strawberry fruit ripening.

A Detailed Spatial Expression Analysis of Wing Phenotypes Reveals Novel Patterns of Odorant Binding Proteins in the Soybean Aphid, Aphis glycines

The wide range of insect niches has led to a rapid expansion of chemosensory gene families as well as their relatively independent evolution and a high variation. Previous studies have revealed some functions for odorant-binding proteins (OBPs) in processes beyond olfaction, such as gustation and reproduction. In this study, a comparative transcriptomic analysis strategy was applied for the soybean aphid, Aphis glycines, focusing on various functional tissues and organs of winged aphids, including the antenna, head, leg, wing, thorax, cauda, and cornicle. Detailed spatial OBP expression patterns in winged and wingless parthenogenetic aphids were detected by RT-qPCR. Twelve OBPs were identified, and three new OBPs in A. glycines are first reported. All OBPs showed comparatively higher expression in sensory organs and tissues, such as the antenna, head, or leg. Additionally, we found some novel expression patterns for aphid OBPs (Beckendorf et al., 2008). Five OBPs exhibited high-expression levels in the cauda and four in the cornicle (Biasio et al., 2015). Three genes (OBP2/3/15) were highly expressed in the wing (Calvello et al., 2003). Two (OBP3/15) were significantly more highly expressed in the wingless thorax than in the winged thorax with the wings removed, and these transcripts were significantly enriched in the removed wings. More details regarding OBP spatial expression were revealed under our strategy. These findings supported the existence of carrier transport functions other than for foreign chemicals and therefore broader ligand ranges of aphid OBPs. It is important for understanding how insect OBPs function in chemical perception as well as their other potential physiological functions.

The NAC transcription factor FaRIF controls fruit ripening in strawberry

In contrast to climacteric fruits such as tomato, the knowledge on key regulatory genes controlling the ripening of strawberry, a nonclimacteric fruit, is still limited. NAC transcription factors (TFs) mediate different developmental processes in plants. Here, we identified and characterized Ripening Inducing Factor (FaRIF), a NAC TF that is highly expressed and induced in strawberry receptacles during ripening. Functional analyses based on stable transgenic lines aimed at silencing FaRIF by RNA interference, either from a constitutive promoter or the ripe receptacle-specific EXP2 promoter, as well as overexpression lines showed that FaRIF controls critical ripening-related processes such as fruit softening and pigment and sugar accumulation. Physiological, metabolome, and transcriptome analyses of receptacles of FaRIF-silenced and overexpression lines point to FaRIF as a key regulator of strawberry fruit ripening from early developmental stages, controlling abscisic acid biosynthesis and signaling, cell-wall degradation, and modification, the phenylpropanoid pathway, volatiles production, and the balance of the aerobic/anaerobic metabolism. FaRIF is therefore a target to be modified/edited to control the quality of strawberry fruits.

Evolution of the MLO gene families in octoploid strawberry (Fragaria ×ananassa) and progenitor diploid species identified potential genes for strawberry powdery mildew resistance

Powdery mildew (PM) caused by Podosphaera aphanis is a major fungal disease of cultivated strawberry. Mildew Resistance Locus O (MLO) is a gene family described for having conserved seven-transmembrane domains. Induced loss-of-function in specific MLO genes can confer durable and broad resistance against PM pathogens. However, the genomic structure and potential role of MLO genes for PM resistance have not been characterized yet in the octoploid cultivated strawberry. In the present study, MLO gene families were characterized in four diploid progenitor species (Fragaria vesca, F. iinumae, F. viridis, and F. nipponica) and octoploid cultivated (Fragaria ×ananassa) strawberry, and potential sources of MLO-mediated susceptibility were identified. Twenty MLO sequences were identified in F. vesca and 68 identified in F. ×ananassa. Phylogenetic analysis divided diploid and octoploid strawberry MLO genes into eight different clades, in which three FveMLO (MLO10, MLO17, and MLO20) and their twelve orthologs of FaMLO were grouped together with functionally characterized MLO genes conferring PM susceptibility. Copy number variations revealed differences in MLO composition among homoeologous chromosomes, supporting the distinct origin of each subgenome during the evolution of octoploid strawberry. Dissecting genomic sequence and structural variations in candidate FaMLO genes revealed their potential role associated with genetic controls and functionality in strawberry against PM pathogen. Furthermore, the gene expression profiling and RNAi silencing of putative FaMLO genes in response to the pathogen indicate the function in PM resistance. These results are a critical first step in understanding the function of strawberry MLO genes and will facilitate further genetic studies of PM resistance in cultivated strawberry.

Single Nucleotide Polymorphisms in Genes Encoding Toll-Like Receptors 7 and 8 and Their Association with Proviral Load of SRLVs in Goats of Polish Carpathian Breed

Toll-like receptors (TLRs) 7 and 8 are important in single-stranded viral RNA recognition, so genetic variation of these genes may play a role in SRLVs infection and disease progression. Present study aimed to identify SNPs in genes encoding TLR7 and TLR8 in goats of Carpathian breed and analyze their association with the SRLVs provirus concentration as index of disease progression. A total of 14 SNPs were detected, 6 SNPs in the TLR7 gene locus and 8 SNPs in the TLR8 gene. Nine of the 14 identified polymorphisms, 4 in the TLR7 gene and 5 in TLR8 gene, were significantly associated with the SRLVs proviral concentration. These SNPs were located in 3'UTR, 5'UTR and intron sequences as well as in the coding sequences, but they led to silent changes. Homozygous genotypes of three TLR7 SNPs (synonymous variant 1:50703293, 3'UTR variant 1:50701297 and 5'UTR variant 1:50718645) were observed in goats with lower provirus copy number as well as in seronegative animals. The results obtained in this study suggest that SNPs of TLR7/TLR8 genes may induce differential innate immune response towards SRLVs affecting proviral concentration and thereby disease pathogenesis and progression. These findings support a role for genetic variations of TLR7 and TLR8 in SRLVs infection and warrants further studies on the effect of TLR7/TLR8 polymorphisms on SRLVs infection in different populations.

Organ-Specific Transcriptome Analysis Identifies Candidate Genes Involved in the Stem Specialization of Bermudagrass (Cynodon dactylon L.)

As an important warm-season turfgrass and forage grass species with wide applications, bermudagrass (Cynodon dactylon L.) simultaneously has shoot, stolon and rhizome, three types of stems with different physiological functions. To better understand how the three types of stems differentiate and specialize, we generated an organ-specific transcriptome dataset of bermudagrass encompassing 114,169 unigenes, among which 100,878 and 65,901 could be assigned to the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the Gene Ontology (GO) terms, respectively. Using the dataset, we comprehensively analyzed the gene expression of different organs, especially the shoot, stolon and rhizome. The results indicated that six organs of bermudagrass all contained more than 52,000 significantly expressed unigenes, however, only 3,028 unigenes were enrich-expressed in different organs. Paired comparison analyses further indicated that 11,762 unigenes were differentially expressed in the three types of stems. Gene enrichment analysis revealed that 39 KEGG pathways were enriched with the differentially expressed unigenes (DEGs). Specifically, 401 DEGs were involved in plant hormone signal transduction, whereas 1,978 DEGs were transcription factors involved in gene expression regulation. Furthermore, in agreement with the starch content and starch synthase assay results, DEGs encoding starch synthesis-related enzymes all showed the highest expression level in the rhizome. These results not only provided new insights into the specialization of stems in bermudagrass but also made solid foundation for future gene functional studies in this important grass species and other stoloniferous/rhizomatous plants.

Genetic Analysis of Methyl Anthranilate, Mesifurane, Linalool, and Other Flavor Compounds in Cultivated Strawberry (Fragaria × ananassa)

The cultivated strawberry (Fragaria × ananassa) is an economically important fruit crop that is intensively bred for improved sensory qualities. The diversity of fruit flavors and aromas in strawberry results mainly from the interactions of sugars, acids, and volatile organic compounds (VOCs) that are derived from diverse biochemical pathways influenced by the expression of many genes. This study integrates multiomic analyses to identify QTL and candidate genes for multiple aroma compounds in a complex strawberry breeding population. Novel fruit volatile QTL was discovered for methyl anthranilate, methyl 2-hexenoate, methyl 2-methylbutyrate, mesifurane, and a shared QTL on Chr 3 was found for nine monoterpene and sesquiterpene compounds, including linalool, 3-carene, β-phellandrene, α-limonene, linalool oxide, nerolidol, α-caryophellene, α-farnesene, and β-farnesene. Fruit transcriptomes from a subset of 64 individuals were used to support candidate gene identification. For methyl esters including the grape-like methyl anthranilate, a novel ANTHANILIC ACID METHYL TRANSFERASE-like gene was identified. Two mesifurane QTL correspond with the known biosynthesis gene O-METHYL TRANSFERASE 1 and a novel FURANEOL GLUCOSYLTRANSFERASE. The shared terpene QTL contains multiple fruit-expressed terpenoid pathway-related genes including NEROLIDOL SYNTHASE 1 (FanNES1). The abundance of linalool and other monoterpenes is partially governed by a co-segregating expression-QTL (eQTL) for FanNES1 transcript variation, and there is additional evidence for quantitative effects from other terpenoid-pathway genes in this narrow genomic region. These QTLs present new opportunities in breeding for improved flavor in commercial strawberry.

Genetic Analysis of Methyl Anthranilate, Mesifurane, Linalool, and Other Flavor Compounds in Cultivated Strawberry (Fragaria × ananassa)

The cultivated strawberry (Fragaria × ananassa) is an economically important fruit crop that is intensively bred for improved sensory qualities. The diversity of fruit flavors and aromas in strawberry results mainly from the interactions of sugars, acids, and volatile organic compounds (VOCs) that are derived from diverse biochemical pathways influenced by the expression of many genes. This study integrates multiomic analyses to identify QTL and candidate genes for multiple aroma compounds in a complex strawberry breeding population. Novel fruit volatile QTL was discovered for methyl anthranilate, methyl 2-hexenoate, methyl 2-methylbutyrate, mesifurane, and a shared QTL on Chr 3 was found for nine monoterpene and sesquiterpene compounds, including linalool, 3-carene, β-phellandrene, α-limonene, linalool oxide, nerolidol, α-caryophellene, α-farnesene, and β-farnesene. Fruit transcriptomes from a subset of 64 individuals were used to support candidate gene identification. For methyl esters including the grape-like methyl anthranilate, a novel ANTHANILIC ACID METHYL TRANSFERASE-like gene was identified. Two mesifurane QTL correspond with the known biosynthesis gene O-METHYL TRANSFERASE 1 and a novel FURANEOL GLUCOSYLTRANSFERASE. The shared terpene QTL contains multiple fruit-expressed terpenoid pathway-related genes including NEROLIDOL SYNTHASE 1 (FanNES1). The abundance of linalool and other monoterpenes is partially governed by a co-segregating expression-QTL (eQTL) for FanNES1 transcript variation, and there is additional evidence for quantitative effects from other terpenoid-pathway genes in this narrow genomic region. These QTLs present new opportunities in breeding for improved flavor in commercial strawberry.

A comprehensive RNA-Seq-based gene expression atlas of the summer squash (Cucurbita pepo) provides insights into fruit morphology and ripening mechanisms

BACKGROUND

Summer squash (Cucurbita pepo: Cucurbitaceae) are a popular horticultural crop for which there is insufficient genomic and transcriptomic information. Gene expression atlases are crucial for the identification of genes expressed in different tissues at various plant developmental stages. Here, we present the first comprehensive gene expression atlas for a summer squash cultivar, including transcripts obtained from seeds, shoots, leaf stem, young and developed leaves, male and female flowers, fruits of seven developmental stages, as well as primary and lateral roots.

RESULTS

In total, 27,868 genes and 2352 novel transcripts were annotated from these 16 tissues, with over 18,000 genes common to all tissue groups. Of these, 3812 were identified as housekeeping genes, half of which assigned to known gene ontologies. Flowers, seeds, and young fruits had the largest number of specific genes, whilst intermediate-age fruits the fewest. There also were genes that were differentially expressed in the various tissues, the male flower being the tissue with the most differentially expressed genes in pair-wise comparisons with the remaining tissues, and the leaf stem the least. The largest expression change during fruit development was early on, from female flower to fruit two days after pollination. A weighted correlation network analysis performed on the global gene expression dataset assigned 25,413 genes to 24 coexpression groups, and some of these groups exhibited strong tissue specificity.

CONCLUSIONS

These findings enrich our understanding about the transcriptomic events associated with summer squash development and ripening. This comprehensive gene expression atlas is expected not only to provide a global view of gene expression patterns in all major tissues in C. pepo but to also serve as a valuable resource for functional genomics and gene discovery in Cucurbitaceae.

Strawberry sweetness and consumer preference are enhanced by specific volatile compounds

Breeding crops for improved flavor is challenging due to the high cost of sensory evaluation and the difficulty of connecting sensory experience to chemical composition. The main goal of this study was to identify the chemical drivers of sweetness and consumer liking for fresh strawberries (Fragaria × ananassa). Fruit of 148 strawberry samples from cultivars and breeding selections were grown and harvested over seven years and were subjected to both sensory and chemical analyses. Each panel consisted of at least 100 consumers, resulting in more than 15,000 sensory data points per descriptor. Three sugars, two acids and 113 volatile compounds were quantified. Consumer liking was highly associated with sweetness intensity, texture liking, and flavor intensity, but not sourness intensity. Partial least square analyses revealed 20 volatile compounds that increased sweetness perception independently of sugars; 18 volatiles that increased liking independently of sugars; and 15 volatile compounds that had positive effects on both. Machine learning-based predictive models including sugars, acids, and volatiles explained at least 25% more variation in sweetness and liking than models accounting for sugars and acids only. Volatile compounds such as γ-dodecalactone; 5-hepten-2-one, 6-methyl; and multiple medium-chain fatty acid esters may serve as targets for breeding or quality control attributes for strawberry products. A genetic association study identified two loci controlling ester production, both on linkage group 6 A. Co-segregating makers in these regions can be used for increasing multiple esters simultaneously. This study demonstrates a paradigm for improvement of fruit sweetness and flavor in which consumers drive the identification of the most important chemical targets, which in turn drives the discovery of genetic targets for marker-assisted breeding.

Down-regulation of Fra a 1.02 in strawberry fruits causes transcriptomic and metabolic changes compatible with an altered defense response

The strawberry Fra a 1 proteins belong to the class 10 Pathogenesis-Related (PR-10) superfamily. In strawberry, a large number of members have been identified, but only a limited number is expressed in the fruits. In this organ, Fra a 1.01 and Fra a 1.02 are the most abundant Fra proteins in the green and red fruits, respectively, however, their function remains unknown. To know the function of Fra a 1.02 we have generated transgenic lines that silence this gene, and performed metabolomics, RNA-Seq, and hormonal assays. Previous studies associated Fra a 1.02 to strawberry fruit color, but the analysis of anthocyanins in the ripe fruits showed no diminution in their content in the silenced lines. Gene ontology (GO) analysis of the genes differentially expressed indicated that oxidation/reduction was the most represented biological process. Redox state was not apparently altered since no changes were found in ascorbic acid and glutathione (GSH) reduced/oxidized ratio, but GSH content was reduced in the silenced fruits. In addition, a number of glutathione-S-transferases (GST) were down-regulated as result of Fra a 1.02-silencing. Another highly represented GO category was transport which included a number of ABC and MATE transporters. Among the regulatory genes differentially expressed WRKY33.1 and WRKY33.2 were down-regulated, which had previously been assigned a role in strawberry plant defense. A reduced expression of the VQ23 gene and a diminished content of the hormones JA, SA, and IAA were also found. These data might indicate that Fra a 1.02 participates in the defense against pathogens in the ripe strawberry fruits.

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.

Reannotation of the cultivated strawberry genome and establishment of a strawberry genome database

Cultivated strawberry (Fragaria × ananassa) is an important fruit crop species whose fruits are enjoyed by many worldwide. An octoploid of hybrid origin, the complex genome of this species was recently sequenced, serving as a key reference genome for cultivated strawberry and related species of the Rosaceae family. The current annotation of the F. ananassa genome mainly relies on ab initio predictions and, to a lesser extent, transcriptome data. Here, we present the structure and functional reannotation of the F. ananassa genome based on one PacBio full-length RNA library and ninety-two Illumina RNA-Seq libraries. This improved annotation of the F. ananassa genome, v1.0.a2, comprises a total of 108,447 gene models, with 97.85% complete BUSCOs. The models of 19,174 genes were modified, 360 new genes were identified, and 11,044 genes were found to have alternatively spliced isoforms. Additionally, we constructed a strawberry genome database (SGD) for strawberry gene homolog searching and annotation downloading. Finally, the transcriptome of the receptacles and achenes of F. ananassa at four developmental stages were reanalyzed and qualified, and the expression profiles of all the genes in this annotation are also provided. Together, this study provides an updated annotation of the F. ananassa genome, which will facilitate genomic analyses across the Rosaceae family and gene functional studies in cultivated strawberry.

Contrasting dynamics in abscisic acid metabolism in different Fragaria spp. during fruit ripening and identification of the enzymes involved

Abscisic acid (ABA) is a key hormone in non-climacteric Fragaria spp, regulating multiple physiological processes throughout fruit ripening. Its concentration increases during ripening, and it promotes fruit (receptacle) development. However, its metabolism in the fruit is largely unknown. We analyzed the concentrations of ABA and its catabolites at different developmental stages of strawberry ripening in diploid and octoploid genotypes and identified two functional ABA-glucosyltransferases (FvUGT71A49 and FvUGT73AC3) and two regiospecific ABA-8'-hydroxylases (FaCYP707A4a and FaCYP707A1/3). ABA-glucose ester content increased during ripening in diploid F. vesca varieties but decreased in octoploid F.×ananassa. Dihydrophaseic acid content increased throughout ripening in all analyzed receptacles, while 7'-hydroxy-ABA and neo-phaseic acid did not show significant changes during ripening. In the studied F. vesca varieties, the receptacle seems to be the main tissue for ABA metabolism, as the concentration of ABA and its metabolites in the receptacle was generally 100 times higher than in achenes. The accumulation patterns of ABA catabolites and transcriptomic data from the literature show that all strawberry fruits produce and metabolize considerable amounts of the plant hormone ABA during ripening, which is therefore a conserved process, but also illustrate the diversity of this metabolic pathway which is species, variety, and tissue dependent.

Genome-wide identification of ubiquitin proteasome subunits as superior reference genes for transcript normalization during receptacle development in strawberry cultivars

Background

Gene transcripts that show invariant abundance during development are ideal as reference genes (RGs) for accurate gene expression analyses, such as RNA blot analysis and reverse transcription–quantitative real time PCR (RT-qPCR) analyses. In a genome-wide analysis, we selected three “Commonly used” housekeeping genes (HKGs), fifteen “Traditional” HKGs, and nine novel genes as candidate RGs based on 80 publicly available transcriptome libraries that include data for receptacle development in eight strawberry cultivars.

Results

The results of the multifaceted assessment consistently revealed that expression of the novel RGs showed greater stability compared with that of the “Commonly used” and “Traditional” HKGs in transcriptome and RT-qPCR analyses. Notably, the majority of stably expressed genes were associated with the ubiquitin proteasome system. Among these, two 26 s proteasome subunits, RPT6A and RPN5A, showed superior expression stability and abundance, and are recommended as the optimal RGs combination for normalization of gene expression during strawberry receptacle development.

Conclusion

These findings provide additional useful and reliable RGs as resources for the accurate study of gene expression during receptacle development in strawberry cultivars.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07393-9.

Chromosome Level Assembly of Homozygous Inbred Line 'Wongyo 3115' Facilitates the Construction of a High-Density Linkage Map and Identification of QTLs Associated With Fruit Firmness in Octoploid Strawberry (Fragaria × ananassa)

Strawberry is an allo-octoploid crop with high genome heterozygosity and complexity, which hinders the sequencing and the assembly of the genome. However, in the present study, we have generated a chromosome level assembly of octoploid strawberry sourced from a highly homozygous inbred line 'Wongyo 3115', using long- and short-read sequencing technologies. The assembly of 'Wongyo 3115' produced 805.6 Mb of the genome with 323 contigs scaffolded into 208 scaffolds with an N50 of 27.3 Mb after further gap filling. The whole genome annotation resulted in 151,892 genes with a gene density of 188.52 (genes/Mb) and validation of a genome, using BUSCO analysis resulted in 94.10% complete BUSCOs. Firmness is one of the vital traits in strawberry, which facilitate the postharvest shelf-life qualities. The molecular and genetic mechanisms that contribute the firmness in strawberry remain unclear. We have constructed a high-density genetic map based on the 'Wongyo 3115' reference genome to identify loci associated with firmness in the present study. For the quantitative trait locus (QTL) identification, the 'BS F₂' populations developed from two inbred lines were genotyped, using an Axiom 35K strawberry chip, and marker positions were analyzed based on the 'Wongyo 3115' genome. Genetic maps were constructed with 1,049 bin markers, spanning the 3,861 cM. Using firmness data of 'BS F₂' obtained from 2 consecutive years, five QTLs were identified on chromosomes 3-3, 5-1, 6-1, and 6-4. Furthermore, we predicted the candidate genes associated with firmness in strawberries by utilizing transcriptome data and QTL information. Overall, we present the chromosome-level assembly and annotation of a homozygous octoploid strawberry inbred line and a linkage map constructed to identify QTLs associated with fruit firmness.

Elucidating the role of polygalacturonase genes in strawberry fruit softening

To disentangle the role of polygalacturonase (PG) genes in strawberry softening, the two PG genes most expressed in ripe receptacles, FaPG1 and FaPG2, were down-regulated. Transgenic ripe fruits were firmer than those of the wild type when PG genes were silenced individually. Simultaneous silencing of both PG genes by transgene stacking did not result in an additional increase in firmness. Cell walls from ripe fruits were characterized by a carbohydrate microarray. Higher signals of homogalacturonan and rhamnogalacturonan I pectin epitopes in polysaccharide fractions tightly bound to the cell wall were observed in the transgenic genotypes, suggesting a lower pectin solubilization. At the transcriptomic level, the suppression of FaPG1 or FaPG2 alone induced few transcriptomic changes in the ripe receptacle, but the amount of differentially expressed genes increased notably when both genes were silenced. Many genes encoding cell wall-modifying enzymes were down-regulated. The expression of a putative high affinity potassium transporter was induced in all transgenic genotypes, indicating that cell wall weakening and loss of cell turgor could be linked. These results suggest that, besides the disassembly of pectins tightly linked to the cell wall, PGs could play other roles in strawberry softening, such as the release of oligogalacturonides exerting a positive feedback in softening.

Allelic Variation of MYB10 Is the Major Force Controlling Natural Variation in Skin and Flesh Color in Strawberry (Fragaria spp.) Fruit

The fruits of diploid and octoploid strawberry (Fragaria spp) show substantial natural variation in color due to distinct anthocyanin accumulation and distribution patterns. Anthocyanin biosynthesis is controlled by a clade of R2R3 MYB transcription factors, among which MYB10 is the main activator in strawberry fruit. Here, we show that mutations in MYB10 cause most of the variation in anthocyanin accumulation and distribution observed in diploid woodland strawberry (F. vesca) and octoploid cultivated strawberry (F ×ananassa). Using a mapping-by-sequencing approach, we identified a gypsy-transposon in MYB10 that truncates the protein and knocks out anthocyanin biosynthesis in a white-fruited F. vesca ecotype. Two additional loss-of-function mutations in MYB10 were identified among geographically diverse white-fruited F. vesca ecotypes. Genetic and transcriptomic analyses of octoploid Fragaria spp revealed that FaMYB10-2, one of three MYB10 homoeologs identified, regulates anthocyanin biosynthesis in developing fruit. Furthermore, independent mutations in MYB10-2 are the underlying cause of natural variation in fruit skin and flesh color in octoploid strawberry. We identified a CACTA-like transposon (FaEnSpm-2) insertion in the MYB10-2 promoter of red-fleshed accessions that was associated with enhanced expression. Our findings suggest that cis-regulatory elements in FaEnSpm-2 are responsible for enhanced MYB10-2 expression and anthocyanin biosynthesis in strawberry fruit flesh.

Host-Pathogen Interactions between Xanthomonas fragariae and Its Host Fragaria × ananassa Investigated with a Dual RNA-Seq Analysis

Strawberry is economically important and widely grown, but susceptible to a large variety of phytopathogenic organisms. Among them, Xanthomonas fragariae is a quarantine bacterial pathogen threatening strawberry productions by causing angular leaf spots. Using whole transcriptome sequencing, the gene expression of both plant and bacteria in planta was analyzed at two time points, 12 and 29 days post inoculation, in order to compare the pathogen and host response between the stages of early visible and of well-developed symptoms. Among 28,588 known genes in strawberry and 4046 known genes in X. fragariae expressed at both time points, a total of 361 plant and 144 bacterial genes were significantly differentially expressed, respectively. The identified higher expressed genes in the plants were pathogen-associated molecular pattern receptors and pathogenesis-related thaumatin encoding genes, whereas the more expressed early genes were related to chloroplast metabolism as well as photosynthesis related coding genes. Most X. fragariae genes involved in host interaction, recognition, and pathogenesis were lower expressed at late-phase infection. This study gives a first insight into the interaction of X. fragariae with its host. The strawberry plant changed gene expression in order to consistently adapt its metabolism with the progression of infection.

FaMYB9 is involved in the regulation of C6 volatile biosynthesis in strawberry

The large-scale untargeted proteomic and metabolomic studies were conducted in strawberry (Fragaria × ananassa) cv. Akihime fruit at five developmental stages. We found that some C6 volatiles highly contributed to the enrichment of volatiles at the red stage of strawberry fruit. We found that 12 genes involved in LOX pathway for volatile biosynthesis showed multiple patterns in protein abundance during fruit development and ripening, and 9 out of the 12 genes exhibited a significant increase in their relative expression levels at the red stage of fruit. We also found that the MYB9 gene (FaMYB9) expression level was positively correlated with the content of C6 volatiles (R = 0.989) and with the relative expression level and protein abundance of FaLOX5 at different strawberry fruit developmental stages (R = 0.954). The interaction between FaMYB9 and FaLOX5 was detected by yeast two-hybrid, co-immunoprecipitation (Co-IP), bimolecular fluorescence complementation (BiFC), and immunofluorescence (IF) analyses. Transient silencing of FaMYB9 delayed the fruit development and ripening, resulting in a significant decrease in the contents of C6 volatiles, while overexpression of FaMYB9 increased the fruit development and ripening and the contents of C6 volatiles in Akihime fruit. Therefore, FaMYB9 is positively involved in C6 volatile biosynthesis.

Characterizing the involvement of FaMADS9 in the regulation of strawberry fruit receptacle development

FaMADS9 is the strawberry (Fragaria x ananassa) gene that exhibits the highest homology to the tomato (Solanum lycopersicum) RIN gene. Transgenic lines were obtained in which FaMADS9 was silenced. The fruits of these lines did not show differences in basic parameters, such as fruit firmness or colour, but exhibited lower Brix values in three of the four independent lines. The gene ontology MapMan category that was most enriched among the differentially expressed genes in the receptacles at the white stage corresponded to the regulation of transcription, including a high percentage of transcription factors and regulatory proteins associated with auxin action. In contrast, the most enriched categories at the red stage were transport, lipid metabolism and cell wall. Metabolomic analysis of the receptacles of the transformed fruits identified significant changes in the content of maltose, galactonic acid-1,4-lactone, proanthocyanidins and flavonols at the green/white stage, while isomaltose, anthocyanins and cuticular wax metabolism were the most affected at the red stage. Among the regulatory genes that were differentially expressed in the transgenic receptacles were several genes previously linked to flavonoid metabolism, such as MYB10, DIV, ZFN1, ZFN2, GT2, and GT5, or associated with the action of hormones, such as abscisic acid, SHP, ASR, GTE7 and SnRK2.7. The inference of a gene regulatory network, based on a dynamic Bayesian approach, among the genes differentially expressed in the transgenic receptacles at the white and red stages, identified the genes KAN1, DIV, ZFN2 and GTE7 as putative targets of FaMADS9. A MADS9-specific CArG box was identified in the promoters of these genes.

A roadmap for research in octoploid strawberry

The cultivated strawberry (Fragaria × ananassa) is an allo-octoploid species, originating nearly 300 years ago from wild progenitors from the Americas. Since that time the strawberry has become the most widely cultivated fruit crop in the world, universally appealing due to its sensory qualities and health benefits. The recent publication of the first high-quality chromosome-scale octoploid strawberry genome (cv. Camarosa) is enabling rapid advances in genetics, stimulating scientific debate and provoking new research questions. In this forward-looking review we propose avenues of research toward new biological insights and applications to agriculture. Among these are the origins of the genome, characterization of genetic variants, and big data approaches to breeding. Key areas of research in molecular biology will include the control of flowering, fruit development, fruit quality, and plant-pathogen interactions. In order to realize this potential as a global community, investments in genome resources must be continually augmented.

Continuous Exposure to Ethylene Differentially Affects Senescence in Receptacle and Achene Tissues in Strawberry Fruit

Strawberry shelf life is limited, and little is known about the postharvest regulation of senescence in different fruit tissues. Strawberry is classified as a non-climacteric fruit, yet it is known that ethylene affects strawberry ripening. Here the effects of continuous exogenous ethylene (50 µl l^-1) were investigated in cold stored strawberry (5°C). The physiological and biochemical responses of ripe strawberry were evaluated across 6 days, together with hormonal profiles of the whole fruit and individual tissues (achenes and receptacle). Continuous exposure to ethylene induced as a first response an accumulation of abscisic acid (ABA) in the receptacle tissue, followed by an increase in CO₂ production. Ethylene also elicited sucrose hydrolysis and malic acid catabolism, with the major effect seen after 4 days of ethylene exposure. Additionally, accumulation of phenolics (epicatechin and chlorogenic acid) were also observed in ethylene treated strawberry. Achenes did not exhibit a response to ethylene, yet catabolism of both ABA and auxins increased by two thirds during air storage. In contrast, ethylene induced ABA accumulation in the receptacle tissue without ABA catabolism being affected. This hormonal disequilibrium in response to ethylene between the two tissues was maintained during storage, and therefore might be the precursor for the following biochemical variations reported during storage.

The Genetics of Differential Gene Expression Related to Fruit Traits in Strawberry (Fragaria ×ananassa)

Octoploid strawberry (Fragaria ×ananassa) is a major specialty crop under intense annual selection for traits relating to plant vigor and fruit quality. Most functional validation experiments rely on transgenic or transient gene expression assays in the mature receptacle. These findings are not typically translatable to breeding without identifying a natural genetic source of transcript level variation, and developing reliable markers for selection in octoploids. Expression QTL (eQTL) analysis is a genetic/transcriptomic association approach for identifying sequence variants predicting differential expression. This eQTL study analyzed a wide array of mature receptacle-expressed genes, encompassing the majority of total mature receptacle transcript accumulation and almost all strawberry genes described in the literature. These results identified segregating genetic variants associated with the differential expression of hundreds of strawberry genes, many with known interest to breeders. Several of these eQTL pertain to published genes whose expression levels have been demonstrated to influence mature receptacle phenotypes. Many include key genes of the phenylpropanoid pathway, vitamin C, carotenoid, pectin, and receptacle carbohydrate/sugar metabolism. These subgenome-specific genetic markers may allow breeders to select for desired ranges of target gene expression. These results may also guide basic research efforts and facilitate the identification of causal genes underlying trait QTL.

The Genome-Wide Analysis of RALF-Like Genes in Strawberry (Wild and Cultivated) and Five Other Plant Species (Rosaceae)

The rapid alkalinization factor (RALF) gene family is essential for the plant growth and development. However, there is little known about these genes among Rosaceae species. Here, we identify 124 RALF-like genes from seven Rosaceae species, and 39 genes from Arabidopsis, totally 163 genes, divided into four clades according to the phylogenetic analysis, which includes 45 mature RALF genes from Rosaceae species. The YISY motif and RRXL cleavage site are typical features of true RALF genes, but some variants were detected in our study, such as YISP, YIST, NISY, YINY, YIGY, YVGY, FIGY, YIAY, and RRVM. Motif1 is widely distributed among all the clades. According to screening of cis-regulatory elements, GO annotation, expression sequence tags (EST), RNA-seq, and RT-qPCR, we reported that 24 RALF genes coding mature proteins related to tissue development, fungal infection, and hormone response. Purifying selection may play an important role in the evolutionary process of RALF-like genes among Rosaceae species according to the result from ka/ks. The tandem duplication event just occurs in four gene pairs (Fv-RALF9 and Fv-RALF10, Md-RALF7 and Md-RALF8, Pm-RALF2 and Pm-RALF8, and Pp-RALF11 and Pp-RALF14) from four Rosaceae species. Our research provides a wide overview of RALF-like genes in seven Rosaceae species involved in identification, classification, structure, expression, and evolution analysis.

RNA-seq, de novo transcriptome assembly and flavonoid gene analysis in 13 wild and cultivated berry fruit species with high content of phenolics

BACKGROUND

Flavonoids are produced in all flowering plants in a wide range of tissues including in berry fruits. These compounds are of considerable interest for their biological activities, health benefits and potential pharmacological applications. However, transcriptomic and genomic resources for wild and cultivated berry fruit species are often limited, despite their value in underpinning the in-depth study of metabolic pathways, fruit ripening as well as in the identification of genotypes rich in bioactive compounds.

RESULTS

To access the genetic diversity of wild and cultivated berry fruit species that accumulate high levels of phenolic compounds in their fleshy berry(-like) fruits, we selected 13 species from Europe, South America and Asia representing eight genera, seven families and seven orders within three clades of the kingdom Plantae. RNA from either ripe fruits (ten species) or three ripening stages (two species) as well as leaf RNA (one species) were used to construct, assemble and analyse de novo transcriptomes. The transcriptome sequences are deposited in the BacHBerryGEN database (http://jicbio.nbi.ac.uk/berries) and were used, as a proof of concept, via its BLAST portal (http://jicbio.nbi.ac.uk/berries/blast.html) to identify candidate genes involved in the biosynthesis of phenylpropanoid compounds. Genes encoding regulatory proteins of the anthocyanin biosynthetic pathway (MYB and basic helix-loop-helix (bHLH) transcription factors and WD40 repeat proteins) were isolated using the transcriptomic resources of wild blackberry (Rubus genevieri) and cultivated red raspberry (Rubus idaeus cv. Prestige) and were shown to activate anthocyanin synthesis in Nicotiana benthamiana. Expression patterns of candidate flavonoid gene transcripts were also studied across three fruit developmental stages via the BacHBerryEXP gene expression browser (http://www.bachberryexp.com) in R. genevieri and R. idaeus cv. Prestige.

CONCLUSIONS

We report a transcriptome resource that includes data for a wide range of berry(-like) fruit species that has been developed for gene identification and functional analysis to assist in berry fruit improvement. These resources will enable investigations of metabolic processes in berries beyond the phenylpropanoid biosynthetic pathway analysed in this study. The RNA-seq data will be useful for studies of berry fruit development and to select wild plant species useful for plant breeding purposes.

The expanding roles of APETALA2/Ethylene Responsive Factors and their potential applications in crop improvement

Understanding the molecular basis of the gene-regulatory networks underlying agronomic traits or plant responses to abiotic/biotic stresses is very important for crop improvement. In this context, transcription factors, which either singularly or in conjugation directly control the expression of many target genes, are suitable candidates for improving agronomic traits via genetic engineering. In this regard, members of one of the largest class of plant-specific APETALA2/Ethylene Response Factor (AP2/ERF) superfamily, which is implicated in various aspects of development and plant stress adaptation responses, are considered high-value targets for crop improvement. Besides their long-known regulatory roles in mediating plant responses to abiotic stresses such as drought and submergence, the novel roles of AP2/ERFs during fruit ripening or secondary metabolites production have also recently emerged. The astounding functional plasticity of AP2/ERF members is considered to be achieved by their interplay with other regulatory networks and signalling pathways. In this review, we have integrated the recently accumulated evidence from functional genomics studies and described their newly emerged functions in plants. The key structural features of AP2/ERF proteins and the modes of their action are briefly summarized. The importance of AP2/ERFs in plant development and stress responses and a summary of the event of their successful applications in crop improvement programs are also provided. Altogether, we envisage that the synthesized information presented in this review will be useful to design effective strategies for improving agronomic traits in crop plants.

Disease Resistance Genetics and Genomics in Octoploid Strawberry

Octoploid strawberry (Fragaria ×ananassa) is a valuable specialty crop, but profitable production and availability are threatened by many pathogens. Efforts to identify and introgress useful disease resistance genes (R-genes) in breeding programs are complicated by strawberry’s complex octoploid genome. Recently-developed resources in strawberry, including a complete octoploid reference genome and high-resolution octoploid genotyping, enable new analyses in strawberry disease resistance genetics. This study characterizes the complete R-gene collection in the genomes of commercial octoploid strawberry and two diploid ancestral relatives, and introduces several new technological and data resources for strawberry disease resistance research. These include octoploid R-gene transcription profiling, dN/dS analysis, expression quantitative trait loci (eQTL) analysis and RenSeq analysis in cultivars. Octoploid fruit eQTL were identified for 76 putative R-genes. R-genes from the ancestral diploids Fragaria vesca and Fragaria iinumae were compared, revealing differential inheritance and retention of various octoploid R-gene subtypes. The mode and magnitude of natural selection of individual F. ×ananassa R-genes was also determined via dN/dS analysis. R-gene sequencing using enriched libraries (RenSeq) has been used recently for R-gene discovery in many crops, however this technique somewhat relies upon a priori knowledge of desired sequences. An octoploid strawberry capture-probe panel, derived from the results of this study, is validated in a RenSeq experiment and is presented for community use. These results give unprecedented insight into crop disease resistance genetics, and represent an advance toward exploiting variation for strawberry cultivar improvement.