Molecular footprints of selection effects and whole genome duplication (WGD) events in three blueberry species: detected by transcriptome dataset

Genome-wide identification of MYB gene family and exploration of selenium metabolism-related candidates in paper mulberry (Broussonetia papyrifera)

KEY MESSAGE

Genome-wide identified 144 MYB family members in B. papyrifera. Integrated correlation analysis and target gene-binding motif prediction indicate that BpMYB135 is vital in regulating selenium metabolism. Selenium is an essential micronutrient for maintaining the health of humans and animals. Broussonetia papyrifera, a forage tree with high nutritional value, exhibits a remarkable ability to accumulate selenium. Although previous studies have preliminarily unfolded the molecular mechanisms underlying selenium accumulation, the roles of transcription factors in regulating selenium uptake and transformation remain poorly understood. This study used various strategies including bioinformatic, physiological, and molecular experiments to explore candidates regarding Se metabolism. Briefly, 144 MYB transcription factor family members were identified and classified into four types (R1, R2R3, R1R2R3, and R4), with phylogenetic analysis further dividing them into 58 subfamilies. The promoters of those BpMYBs contain numerous cis-acting elements associated with plant growth, development, and stress response. qRT-PCR assay confirmed 8 of 15 BpMYBs exhibit a remarkable correlation with selenium content at the threshold absolute value of 0.5. Additionally, foliar application of exogenous abscisic acid (ABA), methyl jasmonate (MeJA), and salicylic acid (SA) reveals different response patterns of BpMYBs. The subcellular localization assay simultaneously verifies that the candidate BpMYB135 functions within the nucleus. Overall, this funding highlights the potential regulatory mechanisms of selenium metabolism in B. papyrifera, providing a foundation for improving its forage value through genetic modification.

Identification of R2R3-MYB family in blueberry and its potential involvement of anthocyanin biosynthesis in fruits

BACKGROUND

Blueberries (Vaccinium corymbosum) are regarded as "superfoods" attributed to large amounts of anthocyanins, a group of flavonoid metabolites, which provide pigmentation in plant and beneficial effects for human health. MYB transcription factor is one of vital components in the regulation of plant secondary metabolism, which occupies a dominant position in the regulatory network of anthocyanin biosynthesis. However, the role of MYB family in blueberry responding to anthocyanin biosynthesis remains elusive.

RESULTS

In this study, we conducted a comprehensive analysis of VcMYBs in blueberry based on the genome data, including phylogenetic relationship, conserved motifs, identification of differentially expressed MYB genes during fruit development and their expression profiling, etc. A total of 437 unique MYB sequences with two SANT domains were identified in blueberry, which were divided into 3 phylogenetic trees. Noticeably, there are many trigenic and tetragenic VcMYBs pairs with more than 95% identity to each other. Meanwhile, the transcript accumulations of VcMYBs were surveyed underlying blueberry fruit development, and they showed diverse expression patterns, suggesting various functional roles in fruit ripening. More importantly, distinct transcript profiles between skin and pulp of ripe fruit were observed for several VcMYBs, such as VcMYB437, implying the potential roles in anthocyanin biosynthesis.

CONCLUSIONS

Totally, 437 VcMYBs were identified and characterized. Subsequently, their transcriptional patterns were explored during fruit development and fruit tissues (skin and pulp) closely related to anthocyanin biosynthesis. These genome-wide data and findings will contribute to demonstrating the functional roles of VcMYBs and their regulatory mechanisms for anthocyanins production and accumulation in blueberry in the future study.

Genome-Wide Identification and Analysis of the MADS-Box Transcription Factor Genes in Blueberry (Vaccinium spp.) and Their Expression Pattern during Fruit Ripening

MADS-box is a class of transcriptional regulators that are ubiquitous in plants and plays important roles in the process of plant growth and development. Identification and analysis of blueberry MADS-box genes can lay a foundation for their function investigations. In the present study, 249 putative MADS-box genes were identified in the blueberry genome. Those MADS-box genes were distributed on 47 out of 48 chromosomes. The phylogenetic and evolutionary analyses showed that blueberry MADS-box genes were divided into 131 type I members and 118 type II members. The type I genes contained an average of 1.89 exons and the type II genes contained an average of 7.83 exons. Motif analysis identified 15 conserved motifs, of which 4 were related to the MADS domain and 3 were related to the K-box domain. A variety of cis-acting elements were found in the promoter region of the blueberry MADS-box gene, indicating that the MADS-box gene responded to various hormones and environmental alterations. A total of 243 collinear gene pairs were identified, most of which had a Ka/Ks value of less than 1. Nine genes belonging to SEP, AP3/PI, and AGL6 subfamilies were screened based on transcriptomic data. The expression patterns of those nine genes were also verified using quantitative PCR, suggesting that VcMADS6, VcMADS35, VcMADS44, VcMADS58, VcMADS125, VcMADS188, and VcMADS212 had potential functions in blueberry fruit ripening. The results of this study provide references for an in-depth understanding of the biological function of the blueberry MADS-box genes and the mechanism of blueberry fruit ripening.

Comparative transcriptome analyses of three Gentiana species provides signals for the molecular footprints of selection effects and the phylogenetic relationships

Species in Gentiana section Cruciata are important alpine plants with a center of diversity and speciation in Qinghai-Tibet Plateau (QTP), and some of these species are sympatrically distributed in northeastern QTP. Studies on genome features and natural selection signatures of sympatric species in section Crucata have been impeded by a lack of genomic resources. Here, we showed transcript characterizations and molecular footprints of selection effects on G. straminea, G. dahurica and G. officinalis based on the comparative transcriptome. A total of 62.97 Gb clean reads were obtained with unigene numbers per species ranging from 141,819 to 236,408 after assembly. We found that these three species had similar distribution of functional categories in different databases, and key enzyme-encoding genes involved in the iridoids biosynthesis were also obtained. The selective pressure analyses indicated that most paired orthologs between these three species were subject to negative selection, and only a low proportion of the orthologs that underwent positive selection were detected. We found that some positive selected genes were involved in "catalytic activity", "metabolic process", "response to stimulus" and "response to stress". Besides, large numbers of SSR primer pairs with transferabilities were successfully designed based on the available transcriptome datasets of three Gentiana species. The phylogenetic relationships reconstructed based on 352 single-copy nuclear genes provided a rough phylogenetic framework for this genus and confirmed the monophyly of section Cruciata. Our study not only provides insights for the natural selection effects on sympatric Gentiana species, but also enhances future genetic breeding or evolutionary studies on Qinjiao species.

Identification of ARF family in blueberry and its potential involvement of fruit development and pH stress response

Background

Auxin responsive factor (ARF) family is one of core components in auxin signalling pathway, which governs diverse developmental processes and stress responses. Blueberry is an economically important berry-bearing crop and prefers to acidic soil. However, the understandings of ARF family has not yet been reported in blueberry.

Results

In the present study, 60 ARF genes (VcARF) were identified in blueberry, and they showed diverse gene structures and motif compositions among the groups and similar within each group in the phylogenetic tree. Noticeably, 9 digenic, 5 trigenic and 6 tetragenic VcARF pairs exhibited more than 95% identity to each other. Computational analysis indicated that 23 VcARFs harbored the miRNA responsive element (MRE) of miR160 or miR167 like other plant ARF genes. Interestingly, the MRE of miR156d/h-3p was observed in the 5’UTR of 3 VcARFs, suggesting a potentially novel post-transcriptional control. Furthermore, the transcript accumulations of VcARFs were investigated during fruit development, and three categories of transcript profiles were observed, implying different functional roles. Meanwhile, the expressions of VcARFs to different pH conditions (pH4.5 and pH6.5) were surveyed in pH-sensitive and tolerant blueberry species, and a number of VcARFs showed different transcript accumulations. More importantly, distinct transcriptional response to pH stress (pH6.5) were observed for several VcARFs (such as VcARF6s and VcARF19-3/19–4) between pH-sensitive and tolerant species, suggesting their potential roles in adaption to pH stress.

Conclusions

Sixty VcARF genes were identified and characterized, and their transcript profiles were surveyed during fruit development and in response to pH stress. These findings will contribute to future research for eliciting the functional roles of VcARFs and regulatory mechanisms, especially fruit development and adaption to pH stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08556-y.

There and back again; historical perspective and future directions for Vaccinium breeding and research studies

The genus Vaccinium L. (Ericaceae) contains a wide diversity of culturally and economically important berry crop species. Consumer demand and scientific research in blueberry (Vaccinium spp.) and cranberry (Vaccinium macrocarpon) have increased worldwide over the crops' relatively short domestication history (~100 years). Other species, including bilberry (Vaccinium myrtillus), lingonberry (Vaccinium vitis-idaea), and ohelo berry (Vaccinium reticulatum) are largely still harvested from the wild but with crop improvement efforts underway. Here, we present a review article on these Vaccinium berry crops on topics that span taxonomy to genetics and genomics to breeding. We highlight the accomplishments made thus far for each of these crops, along their journey from the wild, and propose research areas and questions that will require investments by the community over the coming decades to guide future crop improvement efforts. New tools and resources are needed to underpin the development of superior cultivars that are not only more resilient to various environmental stresses and higher yielding, but also produce fruit that continue to meet a variety of consumer preferences, including fruit quality and health related traits.

Insights into glucosinolate accumulation and metabolic pathways in Isatis indigotica Fort

BACKGROUND

Glucosinolates (GSLs) play important roles in defending against exogenous damage and regulating physiological activities in plants. However, GSL accumulation patterns and molecular regulation mechanisms are largely unknown in Isatis indigotica Fort.

RESULTS

Ten GSLs were identified in I. indigotica, and the dominant GSLs were epiprogoitrin (EPI) and indole-3-methyl GSL (I3M), followed by progoitrin (PRO) and gluconapin (GNA). The total GSL content was highest (over 20 μmol/g) in reproductive organs, lowest (less than 1.0 μmol/g) in mature organs, and medium in fresh leaves (2.6 μmol/g) and stems (1.5 μmol/g). In the seed germination process, the total GSL content decreased from 27.2 μmol/g (of seeds) to 2.7 μmol/g (on the 120th day) and then increased to 4.0 μmol/g (180th day). However, the content of indole GSL increased rapidly in the first week after germination and fluctuated between 1.13 μmol/g (28th day) and 2.82 μmol/g (150th day). Under the different elicitor treatments, the total GSL content increased significantly, ranging from 2.9-fold (mechanical damage, 3 h) to 10.7-fold (MeJA, 6 h). Moreover, 132 genes were involved in GSL metabolic pathways. Among them, no homologs of AtCYP79F2 and AtMAM3 were identified, leading to a distinctive GSL profile in I. indigotica. Furthermore, most genes involved in the GSL metabolic pathway were derived from tandem duplication, followed by dispersed duplication and segmental duplication. Purifying selection was observed, although some genes underwent relaxed selection. In addition, three tandem-arrayed GSL-OH genes showed different expression patterns, suggesting possible subfunctionalization during evolution.

CONCLUSIONS

Ten different GSLs with their accumulation patterns and 132 genes involved in the GSL metabolic pathway were explored, which laid a foundation for the study of GSL metabolism and regulatory mechanisms in I. indigotica.

Comprehensive Analysis of the SBP Family in Blueberry and Their Regulatory Mechanism Controlling Chlorophyll Accumulation

SQUAMOSA Promoter Binding Protein (SBP) family genes act as central players to regulate plant growth and development with functional redundancy and specificity. Addressing the diversity of the SBP family in crops is of great significance to precisely utilize them to improve agronomic traits. Blueberry is an important economic berry crop. However, the SBP family has not been described in blueberry. In the present study, twenty VcSBP genes were identified through data mining against blueberry transcriptome databases. These VcSBPs could be clustered into eight groups, and the gene structures and motif compositions are divergent among the groups and similar within each group. The VcSBPs were differentially expressed in various tissues. Intriguingly, 10 VcSBPs were highly expressed at green fruit stages and dramatically decreased at the onset of fruit ripening, implying that they are important regulators during early fruit development. Computational analysis showed that 10 VcSBPs were targeted by miR156, and four of them were further verified by degradome sequencing. Moreover, their functional diversity was studied in Arabidopsis. Noticeably, three VcSBPs significantly increased chlorophyll accumulation, and qRT-PCR analysis indicated that VcSBP13a in Arabidopsis enhanced the expression of chlorophyll biosynthetic genes such as AtDVR, AtPORA, AtPORB, AtPORC, and AtCAO. Finally, the targets of VcSBPs were computationally identified in blueberry, and the Y1H assay showed that VcSBP13a could physically bind to the promoter region of the chlorophyll-associated gene VcLHCB1. Our findings provided an overall framework for individually understanding the characteristics and functions of the SBP family in blueberry.

Comparing de novo transcriptome assembly tools in di- and autotetraploid non-model plant species

Background

Polyploidy is very common in plants and can be seen as one of the key drivers in the domestication of crops and the establishment of important agronomic traits. It can be the main source of genomic repatterning and introduces gene duplications, affecting gene expression and alternative splicing. Since fully sequenced genomes are not yet available for many plant species including crops, de novo transcriptome assembly is the basis to understand molecular and functional mechanisms. However, in complex polyploid plants, de novo transcriptome assembly is challenging, leading to increased rates of fused or redundant transcripts. Since assemblers were developed mainly for diploid organisms, they may not well suited for polyploids. Also, comparative evaluations of these tools on higher polyploid plants are extremely rare. Thus, our aim was to fill this gap and to provide a basic guideline for choosing the optimal de novo assembly strategy focusing on autotetraploids, as the scientific interest in this type of polyploidy is steadily increasing.

Results

We present a comparison of two common (SOAPdenovo-Trans, Trinity) and one recently published transcriptome assembler (TransLiG) on diploid and autotetraploid species of the genera Acer and Vaccinium using Arabidopsis thaliana as a reference. The number of assembled transcripts was up to 11 and 14 times higher with an increased number of short transcripts for Acer and Vaccinium, respectively, compared to A. thaliana. In diploid samples, Trinity and TransLiG performed similarly good while in autotetraploids, TransLiG assembled most complete transcriptomes with an average of 1916 assembled BUSCOs vs. 1705 BUSCOs for Trinity. Of all three assemblers, SOAPdenovo-Trans performed worst (1133 complete BUSCOs).

Conclusion

All three assembly tools produced complete assemblies when dealing with the model organism A. thaliana, independently of its ploidy level, but their performances differed extremely when it comes to non-model autotetraploids, where specifically TransLiG and Trinity produced a high number of redundant transcripts. The recently published assembler TransLiG has not been tested yet on any plant organism but showed highest completeness and full-length transcriptomes, especially in autotetraploids. Including such species during the development and testing of new assembly tools is highly appreciated and recommended as many important crops are polyploid.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04078-8.

Chromosome-Level Genome Assembly of the American Cranberry (Vaccinium macrocarpon Ait.) and Its Wild Relative Vaccinium microcarpum

The American cranberry (Vaccinium macrocarpon Ait.) is an iconic North American fruit crop of great cultural and economic importance. Cranberry can be considered a fruit crop model due to its unique fruit nutrient composition, overlapping generations, recent domestication, both sexual and asexual reproduction modes, and the existence of cross-compatible wild species. Development of cranberry molecular resources started very recently; however, further genetic studies are now being limited by the lack of a high-quality genome assembly. Here, we report the first chromosome-scale genome assembly of cranberry, cultivar Stevens, and a draft genome of its close wild relative species Vaccinium microcarpum. More than 92% of the estimated cranberry genome size (492 Mb) was assembled into 12 chromosomes, which enabled gene model prediction and chromosome-level comparative genomics. Our analysis revealed two polyploidization events, the ancient γ-triplication, and a more recent whole genome duplication shared with other members of the Ericaeae, Theaceae and Actinidiaceae families approximately 61 Mya. Furthermore, comparative genomics within the Vaccinium genus suggested cranberry-V. microcarpum divergence occurred 4.5 Mya, following their divergence from blueberry 10.4 Mya, which agrees with morphological differences between these species and previously identified duplication events. Finally, we identified a cluster of subgroup-6 R2R3 MYB transcription factors within a genomic region spanning a large QTL for anthocyanin variation in cranberry fruit. Phylogenetic analysis suggested these genes likely act as anthocyanin biosynthesis regulators in cranberry. Undoubtedly, these new cranberry genomic resources will facilitate the dissection of the genetic mechanisms governing agronomic traits and further breeding efforts at the molecular level.

High-Density Linkage Map Construction and QTL Identification in a Diploid Blueberry Mapping Population

Genotyping by sequencing approaches have been widely applied in major crops and are now being used in horticultural crops like berries and fruit trees. As the original and largest producer of cultivated blueberry, the United States maintains the most diverse blueberry germplasm resources comprised of many species of different ploidy levels. We previously constructed an interspecific mapping population of diploid blueberry by crossing the parent F₁#10 (Vaccinium darrowii Fla4B × diploid V. corymbosum W85-20) with the parent W85-23 (diploid V. corymbosum). Employing the Capture-Seq technology developed by RAPiD Genomics, with an emphasis on probes designed in predicted gene regions, 117 F₁ progeny, the two parents, and two grandparents of this population were sequenced, yielding 131.7 Gbp clean sequenced reads. A total of 160,535 single nucleotide polymorphisms (SNPs), referenced to 4,522 blueberry genome sequence scaffolds, were identified and subjected to a parent-dependent sliding window approach to further genotype the population. Recombination breakpoints were determined and marker bins were deduced to construct a high density linkage map. Twelve blueberry linkage groups (LGs) consisting of 17,486 SNP markers were obtained, spanning a total genetic distance of 1,539.4 cM. Among 18 horticultural traits phenotyped in this population, quantitative trait loci (QTLs) that were significant over at least 2 years were identified for chilling requirement, cold hardiness, and fruit quality traits of color, scar size, and firmness. Interestingly, in 1 year, a QTL associated with timing of early bloom, full bloom, petal fall, and early green fruit was identified in the same region harboring the major QTL for chilling requirement. In summary, we report here the first high density bin map of a diploid blueberry mapping population and the identification of several horticulturally important QTLs.