High-quality reference genome and annotation aids understanding of berry development for evergreen blueberry (Vaccinium darrowii)

Genome-wide identification of the WRKY gene family in blueberry (Vaccinium spp.) and expression analysis under abiotic stress

Introduction

The WRKY transcription factor (TF) family is one of the largest TF families in plants and is widely involved in responses to both biotic and abiotic stresses.

Methods

To clarify the function of the WRKY family in blueberries, this study identified the WRKY genes within the blueberry genome and systematically analyzed gene characteristics, phylogenetic evolution, promoter cis-elements, expression patterns, and subcellular localization of the encoded products.

Results

In this study, 57 VcWRKY genes were identified, and all encoding products had a complete WRKY heptapeptide structure and zinc-finger motif. The VcWRKY genes were divided into three subgroups (I-III) by phylogenetic analysis. Group II was divided into five subgroups: IIa, IIb, IIc, IId, and IIe. 57 VcWRKY genes were distributed unevenly across 32 chromosomes. The amino acids ranged from 172 to 841, and molecular weights varied from 19.75 to 92.28 kD. Intra-group syntenic analysis identified 12 pairs of duplicate segments. Furthermore, 34 cis-element recognition sites were identified in the promoter regions of VcWRKY genes, primarily comprising phytohormone-responsive and light-responsive elements. Comparative syntenic maps were generated to investigate the evolutionary relationships of VcWRKY genes, revealing the closest homology to dicotyledonous WRKY gene families. VcWRKY genes were predominantly expressed in the fruit flesh and roots of blueberries. Gene expression analysis showed that the responses of VcWRKY genes to stress treatments were more strongly in leaves than in roots. Notably, VcWRKY13 and VcWRKY25 exhibited significant upregulation under salt stress, alkali stress, and saline-alkali stress, and VcWRKY1 and VcWRKY13 showed notable induction under drought stress. Subcellular localization analysis confirmed that VcWRKY13 and VcWRKY25 function within the nucleus.

Conclusion

These findings establish a foundation for further investigation into the functions and regulatory mechanisms of VcWRKY genes and provide guidance for selecting stress-tolerant genes in the development of blueberry cultivars.

Genome report: First reference genome of Vaccinium floribundum Kunth, an emblematic Andean species

Vaccinium floribundum Kunth, known as "mortiño," is an endemic shrub species of the Andean region adapted to harsh conditions in high-altitude ecosystems. It plays an important ecological role as a pioneer species in the aftermath of deforestation and human-induced fires within paramo ecosystems, emphasizing its conservation value. While previous studies have offered insights into the genetic diversity of mortiño, comprehensive genomic studies are still missing to fully understand the unique adaptations of this species and its population status, highlighting the importance of generating a reference genome for this plant. ONT and Illumina sequencing were used to establish a reference genome for this species. Three different de novo genome assemblies were generated and compared for quality, continuity and completeness. The Flye assembly was selected as the best and refined by filtering out short ONT reads, screening for contaminants and genome scaffolding. The final assembly has a genome size of 529 Mb, containing 1,317 contigs and 97% complete BUSCOs, indicating a high level of integrity of the genome. Additionally, the LTR Assembly Index of 12.93 further categorizes this assembly as a reference genome. The genome of V. floribundum reported in this study is the first reference genome generated for this species, providing a valuable tool for further studies. This high-quality genome, based on the quality and completeness parameters obtained, will not only help uncover the genetic mechanisms responsible for its unique traits and adaptations to high-altitude ecosystems but will also contribute to conservation strategies for a species endemic to the Andes.

Integrated Transcriptomic and Proteomic Analysis Reveals Molecular Mechanisms of the Cold Stress Response during the Overwintering Period in Blueberries (Vaccinium spp.)

In China, the Liaodong Peninsula is an important growing area for blueberries because of the high organic matter content in the soil, the abundance of light, and the large temperature difference between day and night. However, the low temperature and relative humidity of the air during the winter and early spring in the Liaodong Peninsula are the main reasons for the damage to blueberry plants. Here, we documented the transcriptome and proteome dynamics in response to cold stress in three blueberry cultivars ('Northland', 'Bluecrop', and 'Berkeley'). Functional enrichment analysis indicated that many differentially expressed genes (DEGs) and differentially abundant proteins (DAPs) were mainly involved in the pathways of protein processing in the endoplasmic reticulum, the glutathione metabolism pathway, and ribosomes. We identified 12,747 transcription factors (TFs) distributed in 20 families. Based on our findings, we speculated that cold tolerance development was caused by the expression of calcium-related genes (CDPKs and CMLs), glutathione proteins, and TFs (NAC, WRKY, and ERF). Our investigation found that three cultivars experienced cold damage when exposed to temperatures between -9 °C and -15 °C in the field. Therefore, the cold resistance of blueberries during overwintering should not only resist the influence of low temperatures but also complex environmental factors such as strong winds and low relative humidity in the air. The order of cold resistance strength in the three blueberry cultivars was 'Berkeley', 'Bluecrop', and 'Northland'. These results provide a comprehensive profile of the response to cold stress, which has the potential to be used as a selection marker for programs to improve cold tolerance in blueberries.

Genetic and epigenetic signatures for improved breeding of cultivated blueberry

Blueberry belongs to the Vaccinium genus and is a highly popular fruit crop with significant economic importance. It was not until the early twentieth century that they began to be domesticated through extensive interspecific hybridization. Here, we collected 220 Vaccinium accessions from various geographical locations, including 154 from the United States, 14 from China, eight from Australia, and 29 from Europe and other countries, comprising 164 Vaccinium corymbosum, 15 Vaccinium ashei, 10 lowbush blueberries, seven half-high blueberries, and others. We present the whole-genome variation map of 220 accessions and reconstructed the hundred-year molecular history of interspecific hybridization of blueberry. We focused on the two major blueberry subgroups, the northern highbush blueberry (NHB) and southern highbush blueberry (SHB) and identified candidate genes that contribute to their distinct traits in climate adaptability and fruit quality. Our analysis unveiled the role of gene introgression from Vaccinium darrowii and V. ashei into SHB in driving the differentiation between SHB and NHB, potentially facilitating SHB's adaptation to subtropical environments. Assisted by genome-wide association studies, our analysis suggested VcTBL44 as a pivotal gene regulator governing fruit firmness in SHB. Additionally, we conducted whole-genome bisulfite sequencing on nine NHB and 12 SHB cultivars, and characterized regions that are differentially methylated between the two subgroups. In particular, we discovered that the β-alanine metabolic pathway genes were enriched for DNA methylation changes. Our study provides high-quality genetic and epigenetic variation maps for blueberry, which offer valuable insights and resources for future blueberry breeding.

Genome-Wide Identification Analysis of GST Gene Family in Wild Blueberry Vaccinium duclouxii and Their Impact on Anthocyanin Accumulation

Vaccinium duclouxii, a wild blueberry species native to the mountainous regions of southwestern China, is notable for its exceptionally high anthocyanin content, surpassing that of many cultivated varieties and offering significant research potential. Glutathione S-transferases (GSTs) are versatile enzymes crucial for anthocyanin transport in plants. Yet, the GST gene family had not been previously identified in V. duclouxii. This study utilized a genome-wide approach to identify and characterize the GST gene family in V. duclouxii, revealing 88 GST genes grouped into seven distinct subfamilies. This number is significantly higher than that found in closely related species, with these genes distributed across 12 chromosomes and exhibiting gene clustering. A total of 46 members are classified as tandem duplicates. The gene structure of VdGST is relatively conserved among related species, showing closer phylogenetic relations to V. bracteatum and evidence of purifying selection. Transcriptomic analysis and qRT-PCR indicated that VdGSTU22 and VdGSTU38 were highly expressed in flowers, VdGSTU29 in leaves, and VdGSTF11 showed significant expression in ripe and fully mature fruits, paralleling trends seen with anthocyanin accumulation. Subcellular localization identified VdGSTF11 primarily in the plasma membrane, suggesting a potential role in anthocyanin accumulation in V. duclouxii fruits. This study provides a foundational basis for further molecular-level functional analysis of the transport and accumulation of anthocyanins in V. duclouxii, enhancing our understanding of the molecular mechanisms underlying anthocyanin metabolism in this valuable species.

Molecular mechanism of abscisic acid signaling response factor VcbZIP55 to promote anthocyanin biosynthesis in blueberry (Vaccinium corymbosum)

A high content of anthocyanin in blueberry (Vaccinium corymbosum) is an important indicator to evaluate fruit quality. Abscisic acid (ABA) can promote anthocyanin biosynthesis, but since the molecular mechanism is unclear, clarifying the mechanism will improve for blueberry breeding and cultivation regulation. VcbZIP55 regulating anthocyanin synthesis in blueberry were screened and mined using the published Isoform-sequencing, RNA-Seq and qRT-PCR at different fruit developmental stages. Blueberry genetic transformation and transgenic experiments confirmed that VcbZIP55 could promote anthocyanin biosynthesis in blueberry adventitious buds, tobacco leaves, blueberry leaves and blueberry fruit. VcbZIP55 responded to ABA signals and its expression was upregulated in blueberry fruit. In addition, using VcbZIP55 for Yeast one hybrid assay (Y1H) and transient expression in tobacco leaves demonstrated an interaction between VcbZIP55 and a G-Box motif on the VcMYB1 promoter to activate the expression of VcMYB1. This study will lay the theoretical foundation for the molecular mechanisms of phytohormone regulation responsible for anthocyanin synthesis and provide theoretical support for blueberry quality improvement.

Unveiling the evolutionary history of lingonberry (Vaccinium vitis-idaea L.) through genome sequencing and assembly of European and North American subspecies

Lingonberry (Vaccinium vitis-idaea L.) produces tiny red berries that are tart and nutty in flavor. It grows widely in the circumpolar region, including Scandinavia, northern parts of Eurasia, Alaska, and Canada. Although cultivation is currently limited, the plant has a long history of cultural use among indigenous communities. Given its potential as a food source, genomic resources for lingonberry are significantly lacking. To advance genomic knowledge, the genomes for 2 subspecies of lingonberry (V. vitis-idaea ssp. minus and ssp. vitis-idaea var. 'Red Candy') were sequenced and de novo assembled into contig-level assemblies. The assemblies were scaffolded using the bilberry genome (Vaccinium myrtillus) to generate a chromosome-anchored reference genome consisting of 12 chromosomes each with a total length of 548.07 Mb [contig N50 = 1.17 Mb, BUSCO (C%) = 96.5%] for ssp. vitis-idaea and 518.70 Mb [contig N50 = 1.40 Mb, BUSCO (C%) = 96.9%] for ssp. minus. RNA-seq-based gene annotation identified 27,243 and 25,718 genes on the respective assembly, and transposable element detection methods found that 45.82 and 44.58% of the genome were repeats. Phylogenetic analysis confirmed that lingonberry was most closely related to bilberry and was more closely related to blueberries than cranberries. Estimates of past effective population size suggested a continuous decline over the past 1-3 MYA, possibly due to the impacts of repeated glacial cycles during the Pleistocene leading to frequent population fragmentation. The genomic resource created in this study can be used to identify industry-relevant genes (e.g. anthocyanin production), infer phylogeny, and call sequence-level variants (e.g. SNPs) in future research.

Two haplotype-resolved genomes reveal important flower traits in bigleaf hydrangea (Hydrangea macrophylla) and insights into Asterid evolution

The Hydrangea genus belongs to the Hydrangeaceae family, in the Cornales order of flowering plants, which early diverged among the Asterids, and includes several species that are commonly used ornamental plants. Of them, Hydrangea macrophylla is one of the most valuable species in the nursery trade, yet few genomic resources are available for this crop or closely related Asterid species. Two high-quality haplotype-resolved reference genomes of hydrangea cultivars 'Veitchii' and 'Endless Summer' [highest quality at 2.22 gigabase pairs (Gb), 396 contigs, N50 22.8 megabase pairs (Mb)] were assembled and scaffolded into the expected 18 pseudochromosomes. Utilizing the newly developed high-quality reference genomes along with high-quality genomes of other related flowering plants, nuclear data were found to support a single divergence point in the Asterids clade where both the Cornales and Ericales diverged from the euasterids. Genetic mapping with an F1 hybrid population demonstrated the power of linkage mapping combined with the new genomic resources to identify the gene for inflorescence shape, CYP78A5 located on chromosome 4, and a novel gene, BAM3 located on chromosome 17, for causing double flower. Resources developed in this study will not only help to accelerate hydrangea genetic improvement but also contribute to understanding the largest group of flowering plants, the Asterids.

A haplotype-resolved chromosome-scale genome for Quercus rubra L. provides insights into the genetics of adaptive traits for red oak species

Northern red oak (Quercus rubra L.) is an ecologically and economically important forest tree native to North America. We present a chromosome-scale genome of Q. rubra generated by the combination of PacBio sequences and chromatin conformation capture (Hi-C) scaffolding. This is the first reference genome from the red oak clade (section Lobatae). The Q. rubra assembly spans 739 Mb with 95.27% of the genome in 12 chromosomes and 33,333 protein-coding genes. Comparisons to the genomes of Quercus lobata and Quercus mongolica revealed high collinearity, with intrachromosomal structural variants present. Orthologous gene family analysis with other tree species revealed that gene families associated with defense response were expanding and contracting simultaneously across the Q. rubra genome. Quercus rubra had the most CC-NBS-LRR and TIR-NBS-LRR resistance genes out of the 9 species analyzed. Terpene synthase gene family comparisons further reveal tandem gene duplications in TPS-b subfamily, similar to Quercus robur. Phylogenetic analysis also identified 4 subfamilies of the IGT/LAZY gene family in Q. rubra important for plant structure. Single major QTL regions were identified for vegetative bud break and marcescence, which contain candidate genes for further research, including a putative ortholog of the circadian clock constituent cryptochrome (CRY2) and 8 tandemly duplicated genes for serine protease inhibitors, respectively. Genome-environment associations across natural populations identified candidate abiotic stress tolerance genes and predicted performance in a common garden. This high-quality red oak genome represents an essential resource to the oak genomic community, which will expedite comparative genomics and biological studies in Quercus species.

Blueberry and cranberry pangenomes as a resource for future genetic studies and breeding efforts

Domestication of cranberry and blueberry began in the United States in the early 1800s and 1900s, respectively, and in part owing to their flavors and health-promoting benefits are now cultivated and consumed worldwide. The industry continues to face a wide variety of production challenges (e.g. disease pressures), as well as a demand for higher-yielding cultivars with improved fruit quality characteristics. Unfortunately, molecular tools to help guide breeding efforts for these species have been relatively limited compared with those for other high-value crops. Here, we describe the construction and analysis of the first pangenome for both blueberry and cranberry. Our analysis of these pangenomes revealed both crops exhibit great genetic diversity, including the presence-absence variation of 48.4% genes in highbush blueberry and 47.0% genes in cranberry. Auxiliary genes, those not shared by all cultivars, are significantly enriched with molecular functions associated with disease resistance and the biosynthesis of specialized metabolites, including compounds previously associated with improving fruit quality traits. The discovery of thousands of genes, not present in the previous reference genomes for blueberry and cranberry, will serve as the basis of future research and as potential targets for future breeding efforts. The pangenome, as a multiple-sequence alignment, as well as individual annotated genomes, are publicly available for analysis on the Genome Database for Vaccinium-a curated and integrated web-based relational database. Lastly, the core-gene predictions from the pangenomes will serve useful to develop a community genotyping platform to guide future molecular breeding efforts across the family.

The telomere-to-telomere gap-free reference genome of wild blueberry (Vaccinium duclouxii) provides its high soluble sugar and anthocyanin accumulation

Vaccinium duclouxii, endemic to southwestern China, is a berry-producing shrub or small tree belonging to the Ericaceae family, with high nutritive, medicinal, and ornamental value, abundant germplasm resources, and good edible properties. In addition, V. duclouxii exhibits strong tolerance to adverse environmental conditions, making it a promising candidate for research and offering wide-ranging possibilities for utilization. However, the lack of V. duclouxii genome sequence has hampered its development and utilization. Here, a high-quality telomere-to-telomere genome sequence of V. duclouxii was de novo assembled and annotated. All of 12 chromosomes were assembled into gap-free single contigs, providing the highest integrity and quality assembly reported so far for blueberry. The V. duclouxii genome is 573.67 Mb, which encodes 41 953 protein-coding genes. Combining transcriptomics and metabolomics analyses, we have uncovered the molecular mechanisms involved in sugar and acid accumulation and anthocyanin biosynthesis in V. duclouxii. This provides essential molecular information for further research on the quality of V. duclouxii. Moreover, the high-quality telomere-to-telomere assembly of the V. duclouxii genome will provide insights into the genomic evolution of Vaccinium and support advancements in blueberry genetics and molecular breeding.

Blueberry and cranberry pangenomes as a resource for future genetic studies and breeding efforts

Domestication of cranberry and blueberry began in the United States in the early 1800s and 1900s, respectively, and in part owing to their flavors and health-promoting benefits are now cultivated and consumed worldwide. The industry continues to face a wide variety of production challenges (e.g. disease pressures) as well as a demand for higher-yielding cultivars with improved fruit quality characteristics. Unfortunately, molecular tools to help guide breeding efforts for these species have been relatively limited compared with those for other high-value crops. Here, we describe the construction and analysis of the first pangenome for both blueberry and cranberry. Our analysis of these pangenomes revealed both crops exhibit great genetic diversity, including the presence-absence variation of 48.4% genes in highbush blueberry and 47.0% genes in cranberry. Auxiliary genes, those not shared by all cultivars, are significantly enriched with molecular functions associated with disease resistance and the biosynthesis of specialized metabolites, including compounds previously associated with improving fruit quality traits. The discovery of thousands of genes, not present in the previous reference genomes for blueberry and cranberry, will serve as the basis of future research and as potential targets for future breeding efforts. The pangenome, as a multiple-sequence alignment, as well as individual annotated genomes, are publicly available for analysis on the Genome Database for Vaccinium - a curated and integrated web-based relational database. Lastly, the core-gene predictions from the pangenomes will serve useful to develop a community genotyping platform to guide future molecular breeding efforts across the family.

Trait Mapping of Phenolic Acids in an Interspecific (Vaccinium corymbosum var. caesariense × V. darrowii) Diploid Blueberry Population

Blueberries (Vaccinium sect. Cyanococcus) are a dietary source of phenolic acids, including chlorogenic acid (CGA) and related compounds such as acetylated caffeoylquinic acid (ACQA) and caffeoylarbutin (CA). These compounds are known to be potent antioxidants with potential health benefits. While the chemistry of these compounds has been extensively studied, the genetic analysis has lagged behind. Understanding the genetic basis for traits with potential health implications may be of great use in plant breeding. By characterizing genetic variation related to fruit chemistry, breeders can make more efficient use of plant diversity to develop new cultivars with higher concentrations of these potentially beneficial compounds. Using a large interspecific F1 population, developed from a cross between the temperate V. corymbosum var. ceasariense and the subtropical V. darrowii, with 1025 individuals genotyped using genotype-by-sequencing methods, of which 289 were phenotyped for phenolic acid content, with data collected across 2019 and 2020, we have identified loci associated with phenolic acid content. Loci for the compounds clustered on the proximal arm of Vc02, suggesting that a single gene or several closely associated genes are responsible for the biosynthesis of all four tested compounds. Within this region are multiple gene models similar to hydroxycinnamoyl CoA shikimate/quinate hydroxycinnamoyltransferase (HCT) and UDP glucose:cinnamate glucosyl transferase (UGCT), genes known to be involved in the CGA biosynthesis pathway. Additional loci on Vc07 and Vc12 were associated with caffeoylarbutin content, suggesting a more complicated biosynthesis of that compound.

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

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

Molecular characterization of SPL gene family during flower morphogenesis and regulation in blueberry

The SPL gene is a plant-specific transcription factor involved in the regulation of plant growth and development, which have been identified in woody plants. The process of floral bud differentiation affects the timing of flowering and fruit set and regulates plant growth, however, the mechanism of regulation of flower development by SPL genes is less studied. In this study, 56 VcSPL genes were identified in the tetraploid blueberry. The VcSPL gene family was classified into six subfamilies, and analysis of cis-elements showed that VcSPL genes were regulated by light, phytohormones (abscisic acid, MeJA), and low temperature. In the evolutionary analysis, segmental replication may play an important role in VcSPL gene amplification. Interestingly, we also studied diploid blueberry (Bilberry), in which 24 SPL genes were identified, and 36 homologous pairs were found, suggesting a high degree of convergence in the syntenic relationship between blueberry (Vaccinium corymbosum L) and bilberry (Vaccinium darrowii). Based on the expression profile, VcSPL genes were expressed at high levels in flowers, shoots, and roots, indicating a diversity of gene functions. Then we selected 20 differentially-expressed SPL genes to further investigate the role of VcSPL in floral induction and initiation. It showed that the genes VcSPL40, VcSPL35, VcSPL45, and VcSPL53 may play a crucial role in the blueberry floral transition phase (from vegetative growth to flower initiation). These results provided important information for understanding and exploring the role of VcSPLs in flower morphogenesis and plant growth.

Full-length fruit transcriptomes of southern highbush (Vaccinium sp.) and rabbiteye (V. virgatum Ait.) blueberry

Background

Blueberries (Vaccinium sp.) are native to North America and breeding efforts to improve blueberry fruit quality are focused on improving traits such as increased firmness, enhanced flavor and greater shelf-life. Such efforts require additional genomic resources, especially in southern highbush and rabbiteye blueberries.

Results

We generated the first full-length fruit transcriptome for the southern highbush and rabbiteye blueberry using the cultivars, Suziblue and Powderblue, respectively. The transcriptome was generated using the Pacific Biosciences single-molecule long-read isoform sequencing platform with cDNA pooled from seven stages during fruit development and postharvest storage. Raw reads were processed through the Isoseq pipeline and full-length transcripts were mapped to the ‘Draper’ genome with unmapped reads collapsed using Cogent. Finally, we identified 16,299 and 15,882 non-redundant transcripts in ‘Suziblue’ and ‘Powderblue’ respectively by combining the reads mapped to Northern Highbush blueberry ‘Draper’ genome and Cogent analysis. In both cultivars, > 80% of sequences were longer than 1,000 nt, with the median transcript length around 1,700 nt. Functionally annotated transcripts using Blast2GO were > 92% in both ‘Suziblue’ and ‘Powderblue’ with overall equal distribution of gene ontology (GO) terms in the two cultivars. Analyses of alternative splicing events indicated that around 40% non-redundant sequences exhibited more than one isoform. Additionally, long non-coding RNAs were predicted to represent 5.6% and 7% of the transcriptomes in ‘Suziblue’ and ‘Powderblue’, respectively. Fruit ripening is regulated by several hormone-related genes and transcription factors. Among transcripts associated with phytohormone metabolism/signaling, the highest number of transcripts were related to abscisic acid (ABA) and auxin metabolism followed by those for brassinosteroid, jasmonic acid and ethylene metabolism. Among transcription factor-associated transcripts, those belonging to ripening-related APETALA2/ethylene-responsive element-binding factor (AP2/ERF), NAC (NAM, ATAF1/2 and CUC2), leucine zipper (HB-zip), basic helix-loop-helix (bHLH), MYB (v-MYB, discovered in avian myeloblastosis virus genome) and MADS-Box gene families, were abundant.

Further we measured three fruit ripening quality traits and indicators [ABA, and anthocyanin concentration, and texture] during fruit development and ripening. ABA concentration increased during the initial stages of fruit ripening and then declined at the Ripe stage, whereas anthocyanin content increased during the final stages of fruit ripening in both cultivars. Fruit firmness declined during ripening in ‘Powderblue’. Genes associated with the above parameters were identified using the full-length transcriptome. Transcript abundance patterns of these genes were consistent with changes in the fruit ripening and quality-related characteristics.

Conclusions

A full-length, well-annotated fruit transcriptome was generated for two blueberry species commonly cultivated in the southeastern United States. The robustness of the transcriptome was verified by the identification and expression analyses of multiple fruit ripening and quality–regulating genes. The full-length transcriptome is a valuable addition to the blueberry genomic resources and will aid in further improving the annotation. It will also provide a useful resource for the investigation of molecular aspects of ripening and postharvest processes.