Transcript Profile of Flowering Regulatory Genes in VcFT-Overexpressing Blueberry Plants

In order to identify genetic components in flowering pathways of highbush blueberry (Vaccinium corymbosum L.), a transcriptome reference composed of 254,396 transcripts and 179,853 gene contigs was developed by assembly of 72.7 million reads using Trinity. Using this transcriptome reference and a query of flowering pathway genes of herbaceous plants, we identified potential flowering pathway genes/transcripts of blueberry. Transcriptome analysis of flowering pathway genes was then conducted on leaf tissue samples of transgenic blueberry cv. Aurora (‘VcFT-Aurora’), which overexpresses a blueberry FLOWERING LOCUS T-like gene (VcFT). Sixty-one blueberry transcripts of 40 genes showed high similarities to 33 known flowering-related genes of herbaceous plants, of which 17 down-regulated and 16 up-regulated genes were identified in ‘VcFT-Aurora’. All down-regulated genes encoded transcription factors/enzymes upstream in the signaling pathway containing VcFT. A blueberry CONSTANS-LIKE 5-like (VcCOL5) gene was down-regulated and associated with five other differentially expressed (DE) genes in the photoperiod-mediated flowering pathway. Three down-regulated genes, i.e., a MADS-AFFECTING FLOWERING 2-like gene (VcMAF2), a MADS-AFFECTING FLOWERING 5-like gene (VcMAF5), and a VERNALIZATION1-like gene (VcVRN1), may function as integrators in place of FLOWERING LOCUS C (FLC) in the vernalization pathway. Because no CONSTAN1-like or FLOWERING LOCUS C-like genes were found in blueberry, VcCOL5 and VcMAF2/VcMAF5 or VRN1 might be the major integrator(s) in the photoperiod- and vernalization-mediated flowering pathway, respectively. The major down-stream genes of VcFT, i.e., SUPPRESSOR of Overexpression of Constans 1-like (VcSOC1), LEAFY-like (VcLFY), APETALA1-like (VcAP1), CAULIFLOWER 1-like (VcCAL1), and FRUITFULL-like (VcFUL) genes were present and showed high similarity to their orthologues in herbaceous plants. Moreover, overexpression of VcFT promoted expression of all of these VcFT downstream genes. These results suggest that VcFT’s down-stream genes appear conserved in blueberry.


Introduction
Blueberry contains high amounts of antioxidants known to be important for human health [1]. The highbush blueberry (2n = 4x = 48) (Vaccinium corymbosum L.), including northern and southern ecotypes, is the major cultivated Vaccinium fruit crop [2,3]. The northern highbush cultivars have better winter hardiness, but their adoption in warm areas is often limited by their requirement for more chill units (CU) (generally > 800) to break dormancy in the spring. In contrast, the southern highbush blueberry cultivars are derivatives of the northern highbush blueberry with additional genes from other southern Vaccinum species [e.g., V. darrowi (2n = 2x = 24), V. ashei (2n = 6x = 72), and V. tenellum (2n = 2x = 24)] [2]. These southern cultivars have better summer heat tolerance and require fewer chill units (150 to 600 CU) to induce flowering, but they are often less cold/freeze-tolerant. Developing new cultivars with different flowering times and chilling requirements, high cold/heat tolerances, and high yields are the top priorities in breeding for sustainable blueberry production, particularly in anticipation of climate changes.
Abundant genetic resources of the genus Vaccinium (i.e. 450 species) contain high genetic diversity but also lead to great complexity making it difficult to improve the cultivated tetraploid highbush cultivars through conventional breeding [2]. Alternatively, new genomic and biotechnological tools will facilitate studies on blueberry genetics and genomics and improve breeding efficiency [2,[4][5][6][7][8][9][10]. For blueberry genomic studies, a draft blueberry genome assembly with 358 million base pairs (Mb) of diploid blueberry accession 'W8520' was generated and used for genetic diversity studies of cultivated blueberry cultivars as well as marker-assisted breeding and comparative genomics of Vaccinium species [2]. Using RNA-Seq data and a draft blueberry genome assembly, candidate genes involved in fruit ripening and biosynthesis of bioactive compounds have been reported [7]. In addition, 22,401 blueberry (V. corymbosum) expressed sequence tags (ESTs) are available in GenBank. Generation and analyses of the transcriptomes of flower and fruit tissues of V. corymbosum 'Northland' yielded 34,464 NCBI Unigene clusters out of 64 million sequencing reads for identifying genes involved in antioxidant biosynthesis [8]. Transcriptomes of the major highbush blueberry cultivar Bluecrop were generated from leaves, flower buds at different stages of cold acclimation, and fruit tissues at different stages of development. About 15,000 contigs were derived from 600,000 reads and are considered to be useful for analyzing differentially expressed (DE) genes in flower bud formation, vernalization, cold acclimation, and fruit development [9].
Reliable and highly effective Agrobacterium tumefaciens-mediated transformation and regeneration methods have been developed for several highbush blueberry cultivars [11,12]. The availability of blueberry transcriptome information and genetic transformation tools enables gene cloning and functional gene analysis. For example, overexpressing the endogenous blueberry C-repeat binding factor (CBF) gene (VcCBF) (GenBank AF234316) resulted in a significant increase of freeze tolerance in leaves and floral buds [10]. These results demonstrate the potential of genomics and genetic transformation for blueberry breeding.
Flowering is controlled and modulated by flowering genes and environmental signals (e.g., light and temperature) and is of importance for sustainable blueberry production. To date, the molecular basis of flowering in Arabidopsis thaliana (L.) and other herbaceous plants has been extensively studied through investigations into the networks of flowering-related genes [13][14][15][16][17][18][19][20][21]. In contrast, few similarly detailed studies have been conducted to reveal genetic control of seasonal flowering of woody perennials [14,22]. Of the reported flowering-related genes of A. thaliana, the major ones include CONSTANS 1 (CO1), FLOWERING LOCUS C (FLC), FLOWERING LOCUS T (FT), SUPPRESSOR of Overexpression of Constans 1 (SOC1), LEAFY (LFY), APE-TALA1 (AP1), and FRUITFULL (FUL) [23]. To study the vernalization-mediated flowering mechanisms in blueberry, we isolated one blueberry FT-like gene (VcFT), three AP1-like genes, three SOC1-like genes and one LFY-like gene (VcLFY) from the highbush blueberry cvs. Bluecrop and Legacy. Overexpressing VcFT caused early and continuous flowering in in vitro shoots and in one-year old 'Aurora' plants [5]. However, our observation on two-and three-year old plants showed that the VcFT-overexpressing plants did not flower normally and the majority of the flower buds did not open under greenhouse conditions without chilling. Analyses of DE genes in VcFT-overexpressing blueberry plants should reveal VcFT-regulated flowering gene networks.
The aim of this study is to identify flowering-related gene networks of highbush blueberry. We developed a transcriptome reference of blueberry, assembled using Trinity [24], and annotated flowering pathway genes by comparative genomics of A. thaliana, rice (Oryza sativa), and some cereals, were used to identify flowering-related genes of blueberry [25]. DE analysis of flowering-related genes was subsequently conducted on the transcriptome from leaf tissues of transgenic 'Aurora' plants (herein 'VcFT-Aurora') overexpressing VcFT. In this way, we expected to identify functional flowering pathway genes of blueberry, which respond to the major flowering pathway integrator VcFT.

Plant materials
The northern highbush blueberry cv. Aurora and the southern highbush cv. Legacy were used. Both cultivars are tetraploid and require at least 800 CU to flower. The reference transcriptome was developed using eight 'Legacy' tissue samples derived from a wild type plant, one plant from each of two transgenic events containing an overexpressed VcCBF (herein 'VcCBF-Legacy') obtained in our previous study [10]. These tissue samples included unvernalized leaves and floral buds and flower tissues with corollas removed. Three plants of wild type 'Aurora' and three plants from one transgenic line of 'VcFT-Aurora', a transgenic line overexpressing VcFT obtained in our previous study [5], were used to analyze differential expression resulting from VcFT overexpression. Differential expression comparisons across tissue types (leaves, floral buds, and flower) utilized three plants of wild type 'Legacy'.
Growth and flowering of non-transgenic 'Aurora' and transgenic 'VcFT-Aurora' plants were investigated using 12 non-transgenic plants and 12 plants from each of the five transgenic events obtained in our previous study [5]. Phenotypic data (i.e., plant height, flowering time, the number of floral buds that flowered, and flower number per flowered bud) were collected. For differential expression analyses, young leaf tissues of 3 two-year old, unvernalized 'Aurora' plants and 1 transgenic 'VcFT-Aurora' event of were sampled for RNA isolation and sequencing.
All plants were from in vitro cultured shoots and grown in the greenhouse (heated for winter) under natural light and ambient temperature conditions, and a regular schedule of irrigation and fertilization using 0.2 g/L fertilizer (Nitrogen: Phosphorus: Potassium = 21: 7: 7) [11]. To get full vernalization, one-year old plants were grown in controlled environment chambers at 4°C with a 12 hour photoperiod for two months; two-and three-year old plants were exposed to the natural environment in winter in a secured courtyard between our greenhouses.

RNA preparation and sequencing
Floral buds were collected in November 2013 before the plants were exposed to a non-heated greenhouse for chilling. Newly emerged young leaf and flower tissues of 'Legacy' were collected in February 2014. Flower tissues consisted of all parts from newly opened flowers except corollas. Young leaf tissues of two-year old 'Aurora' and 'VcFT-Aurora' plants were collected in June 2014 from unvernalized plants. All tissues collected were frozen immediately in liquid nitrogen and stored at -80°C.
Total RNA was isolated from 0.5 g tissues for each sample, using a cetyltrimethylammonium bromide (CTAB) method [26]. The samples were purified using RNeasy Mini Kit and On-Column DNase digestion with the RNase-free DNase Set (Qiagen, Valencia, CA, USA). Integrity of the RNA samples was assessed using the Agilent RNA 6000 Pico Kit (Agilent Technologies, Inc. Germany). All samples had a RNA quality score above 8.0 before submission for sequencing (100-bp paired-end reads) using the Illumina HiSeq2500 platform at the Research Technology Support Facility of Michigan State University (East Lansing, MI, USA).

de novo transcriptome assembly
FastQC program (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) was used to assess the quality of sequencing reads for the per base quality scores range from 26-40. Generation of the transcriptome reference for highbush blueberry from RNA-seq was carried out using the Trinity platform (trinity/20140413p1) [24]. The command for running about 72 million reads (MR) was "Trinity-seqType fq-JM 50G -left reads_R1_CBFVF.fq-right reads_R2_CBFVF.fq-SS_lib_type RF-CPU 24". The memory request was 250G and the wall-time request was 48 hr. Left/Right reads were a combination of the reads from eight samples of non-transgenic and transgenic 'Legacy'. The quality of the assembly was assessed by alignment with existing EST singlets and contigs in the Unigene v1 library obtained from www.vaccinium.org. These sequences were derived from a variety of tissues and cultivars of highbush blueberry. BLASTplus version 2.2.27 was used to run blastn with default parameters except dust was off, percent identity ranged from 80-95%, with an evalue threshold of 1e-05. Unigene sequences were used as query with refTrinity as database. Results were filtered for completeness, calculated as (qend-(qstart-1))/qlen.

Differential expression analysis
RNA-seq reads of three biological replicates for each of the 'Legacy' tissues, plus 'Aurora' and 'VcFT-Aurora' leaves were analyzed. Two technical replicates were sequenced for each biological replicate and combined together for analysis. The paired reads, two sets for each biological replicate, were aligned to the transcriptome reference developed for 'Legacy' and the abundance of each read was estimated using the Trinity command "align_and_estimate_abundance.pl". This generated 12 "genes.counts.matrix" and 12 "transcripts.counts.matrix" files. The Trinity command "abundance_estimates_to_matrix.pl-est_method RSEM" was then used to merge "genes.counts.matrix" files and "transcripts.counts.matrix" files, separately. The resultant two matrices were used for running DE analysis using the Trinity command "run_-DE_analysis.pl-method edgeR". The DE genes or transcripts with false discovery rate (FDR) values below 0.05 were used for further analyses of flowering pathway genes of blueberry.
For RT-PCR analysis of differentially expressed genes/transcripts, reverse-transcription of RNA to cDNA was performed using SuperScript II reverse transcriptase (Invitrogen, Carlsbad, CA, USA). The resulting cDNA was diluted (volume 1:10) in water, and 2 μl/sample, was used for each PCR reaction. Three RNA samples from leaf tissues for each of non-transgenic 'Aurora' and transgenic 'VcFT-Aurora' were used. The primers were included in S1 Table. The reaction conditions for all primer pairs were 94°C for 2 min, 35 cycles of 45 s at 94°C, 60 s at 60°C and 90 s at 72°C, with a final 10 min extension at 72°C. RT-PCR products were separated on 1.0% agarose gel containing ethidium bromide visualized, and photographed under UV light.

Sequence sources of plant flowering pathway genes
Three protein datasets of the flowering pathway genes of Arabidopsis, rice, and cereal crops [i.e., bread wheat (Triticum aestivum), einkorn wheat (T. monococcum), barley (Hordeum vulgare), and maize (Zea mays)] were developed according to the gene IDs reported by [25] (S2 Table). Each dataset was used for orthologue identification using tblastn against the transcriptome reference of blueberry. The top hits (i.e., transcripts and genes) with e-values less than e-20, were retained to do comparative genomics of blueberry flowering genes and to search for the DE profiles of VcFT-overexpressing 'Aurora'.

Annotation of DE genes/transcripts
Trinotate was used to annotate the transcriptome reference and DE transcripts. Blueberry flowering pathway transcripts were also used as queries to search the A. thaliana protein database and GenBank using blastx. All Trinity and Trinotate analyses were performed using the resources at the High Performance Computing Center of Michigan State University.

Phenotypes of VcFT overexpressing plants
'VcFT-Aurora' constitutively expresses VcFT driven by the cauliflower mosaic virus (CaMV) 35S promoter [5]. Overexpressing VcFT in 'VcFT-Aurora' resulted in dwarf plants, early floral bud formation, and continuous flowering (Fig 1A-1C). For example, the height of one-or twoyear old 'VcFT-Aurora' plants was about half that of non-transgenic 'Aurora' and some newly formed buds flowered during the growing season ( Fig 1A). 'Aurora' plants did not develop floral buds until they were three-years old. Fully chilled plants flowered with a bloom period of about one week and each bud had 5-10 flowers, while unvernalized plants did not flower. In contrast, 'VcFT-Aurora' shoots started to flower while they were in in vitro culture [5]. Fully chilled 'VcFT-Aurora' plants flowered normally with 5-10 flowers/bud (Fig 1D). The plants then formed new shoots and reverted to continuous floral bud formation and flowering. For unvernalized 'VcFT-Aurora' plants, 20-50% of floral buds flowered with 2-3 flowers/bud ( Fig  1A). These results suggest that over-expressing VcFT alone promotes early flowering, but is not sufficient to completely substitute the need for vernalization for blueberries to flower normally under non-chilling conditions.

Development of a transcriptome reference for blueberry
Approximately 8-10 million pair-end reads (MR) were generated for each sample. Using Trinity [24], 72.7 MR obtained from eight samples of leaf, floral bud, and flower tissues were assembled. The total Trinity 'genes' and 'transcripts' were 179,853 and 254,196, respectively ( Table 1). All these transcript contigs, a total of 192 million bases, was designated as refTrinity and were used as the transcriptome reference for identifying candidate genes and differential expression analyses. The quality of the assembly was assessed by aligning to known EST and gene contig sequences obtained from the database at www.vaccinium.org. refTrinity aligned to 61% of 2,955 ESTs and 70% of 748 contig sequences with at least 90% identity and 90% completeness (S3 Table). Considering that the query sequences originated from varied cultivars, tissues, and treatments, the results suggest that Trinity is an efficient platform in generating transcriptome references for tetraploid blueberry plants.

Identification of flowering pathway genes of blueberry
Abundant sequence data (i.e., DNA, cDNA, EST, and protein) of flowering pathway genes from many plant species are deposited in GenBank. To identify orthologues of flowering pathway genes of blueberry, we used protein sequences derived from the previous reports and a cutoff e-value of above e-20 to identify flowering pathway genes in blueberry (S2 Table) [25].
When 156 sequences of 81 flowering-related genes of Arabidopsis (TAIR 10) were used to blast against refTrinity, 127 sequences of 57 genes got hits. Of the query of 160 sequences of 103 flowering pathway genes of rice (RGAP 7), only OsF gene (LOC_Os02g08150.1) did not have a match. With the query of 49 flowering-related genes of cereal crops, 48 genes had matching sequences. Overall, we pooled 1,223 transcript contigs from 720 gene contigs that showed high similarities (e-values e-20) to the query sequences, which allowed us to generate a reference of potential flowering pathway genes of blueberry (herein Trinity_floral_ref) (S4 Table).
Arabidopsis CO1 and rice HEADING DATE (HD1) are the major integrator genes in their respective photoperiod-regulated flowering pathways [25]. Using the CO1 as a query to search online Vaccinium sequences, we did not get any sequence with a high similarity. In Trinity_flor-al_ref, we found three CONSTANS-LIKE (COL)-like genes (i.e., VcCOL2, VcCOL4, and VcCOL5), but no CO1-like (VcCO1) or HD1-like genes.

Tissue-specific expression of different flowering pathway genes
Comparative analyses of three pools of Trinity transcripts (i.e., leaf vs. floral bud, flower vs. floral bud, and leaf vs. flower) of non-transgenic blueberry plants showed that 1,193 out of 1,223 transcript contigs of the Trinity_floral_ref were found in the combined contigs of leaf, bud, and flower tissues of non-transgenic 'Legacy'. Furthermore, 1,122 (out of 1,193) contigs were shared among three tissue types, of which 397 contigs showed no differential expression between the tissues. The remaining 725 transcripts (69 genes) showed differential expression in the comparisons of leaf vs. bud, flower vs. bud, or leaf vs. flower (Table 2). For example, of the photoperiod pathway genes, VcCOL5 is a DE gene with high expression in leaf; but VcCOL2 and VcCOL4 did not show differential expression. Of the genes of the vernalization pathway, homologues of MAF2, VRN1, and TaVRN1 showed the highest expression in leaf tissues and VcVRN5 had the highest expression in buds; but VcFRI and VcMAF5 were not differentially expressed in the comparisons ( Table 2). The results indicate that expression of some flowering-  related transcripts/genes in blueberry is tissue specific. In addition, RNA-seq transcriptome analysis provides a very efficient approach to confirm tissue specificity of gene expression.
In addition, 1,130 transcript contigs (92.4% of the 1,223) of blueberry flowering-related genes were identified. The absence of some gene/transcript contigs is attributed to either tissue-specific or genotype-dependent expression. Differential expression analyses of transcriptome of 'VcFT-Aurora' indicated that overexpressed VcFT had a broad impact on global gene expression. We found expression levels of 3,023 Trinity gene contigs (2.4% of 124,560) were significantly altered in 'VcFT-Aurora' plants (p-value<0.001 and FDR<0.05), including 1,498 up-regulated and 1,525 down-regulated gene contigs (Fig 2). At the transcript level, 4,844 Trinity transcripts (2.6% of 187,779) were significantly altered, of which 2,270 were up-regulated and 2,574 were down-regulated (Fig 2). Sixtyone transcripts of 40 gene contigs showed high similarities to 33 known flowering-related genes, of which 17 were down-regulated and 16 were up-regulated (Table 3). All down-regulated flowering-related genes were the upstream genes of VcFT. The 35S promoter drove high expression of the VcFT, which was associated with up-regulated TWIN SISTER of FT (TSF)like (VcTSF) as well as the downstream genes VcAP1, VcFUL, SEPALLATA-like (VcSPL), VcLFY, and VcSOC1. The up-regulation of the selected VcAP1, VcLFY, and VcFUL was confirmed in our RT-PCR analysis (Fig 3). Interestingly, nine VcFT upstream genes were also upregulated (Table 3).
In 'VcFT-Aurora', VcCOL5 is a major down-regulated flowering pathway gene associated with five other DE genes in the photoperiod-regulated flowering pathway (Table 3). VcMAF2 and VcMAF5 or VcVRN1 may function as integrators of the vernalization pathway, in which six other genes [i.e., VcFRI, VcHUA2, SUPPRESSOR OF FRIGIDA4-like (VcSUF4), 20S PRO-TEASOME SUBUNIT PAF2-like (VcPAF2), and ACTIN-RELATED PROTEIN 6 (ARP6)-like gene (VcARP6)] were involved (Table 3). No VcCO1 or VcFLC were found in blueberry. These results suggest VcCOL5 and VcMAF2/VcMAF5 or VcVRN1 in blueberry may be the major integrator(s) in the photoperiod and vernalization pathways, respectively. Comparison of transcript expression profiles between 'VcFT-Aurora' and 'Aurora' to identify differentially expressed genes (a, b) and isoforms (c, d) in leaf tissues. a, c, MA plot for differential expression analysis generated by EdgeR: for each gene or isoform, the log 2 (fold change) (log 2 (Aurora/VcFT-Aurora)) between the two samples is plotted (A, y axis) against the gene's log 2 (average expression) in the two samples (M, x axis). b, d, Volcano plot reporting false discovery rate (−logFDR, y axis) as a function of log 2 (fold change) between the samples (logFC, x axis). Genes or isoforms that are identified as significantly differentially expressed at most 0.1% FDR are colored in red.

Discussion
The genetic nature of many fruit crops, e.g., polyploidy, high heterozygosity, long juvenility, and clonal propagation, makes traditional breeding a lengthy and expensive process. Although a few genetic markers have been developed to facilitate the selection process [9] blueberry breeding is mainly guided by breeders' experiences. Manipulation of flowering pathway genes is anticipated to change the regime of blueberry flowering and to increase fruit production [4,5]. To this end, we here have used comparative genomics and biotechnology tools to investigate flowering pathway genes and flowering mechanisms of blueberry.

Transcriptome reference of blueberry
In a previous study of anthocyanin biosynthesis in blueberry, over 64 million RNA-sequencing reads were generated from blueberry skin and pulp, 34,464 unigenes were obtained through de novo assembly [8]. The first whole genome sequence of blueberry with 15,129 scaffolds for 358 Mb was reported on a diploid V. corymbosum accession 'W8520' [2,7]. Using this draft genome reference and RNA-seq data, great efforts have been made to identify genes involved in blueberry ripening and over 400 biosynthetic pathways (e.g., ethylene and anthocyanin) [7]. Currently, there is an online blueberry genomics database (BBDG) that contains a collection of EST sequences (http://bioinformatics.towson.edu/BBGD/)), and some published RNA-seq data are available in the GenBank [8]. To date, no blueberry genome reference has been made available to the public.
We demonstrated that the Trinity platform [24] is amenable to use and is very efficient in generating transcriptome references for tetraploid blueberry. With 72.7 MR (8-10 MR/sample), we believe our transcriptome reference (refTrinity) covers the majority of transcripts in leaves, flowers, and flower buds. Although this transcriptome reference does not include rootand fruit-specific transcripts, it provides a good coverage of the majority of blueberry flowering pathway genes for studying gene expression in blueberry plants.

Conserved flowering pathway genes in blueberry
In this study, it is not surprising that the orthologues of the majority of the flowering pathway genes found in Arabidopsis, rice, and cereals crops were identified in blueberry, although some of these orthologues may have distinctive functions to meet the specific needs of blueberry floral transitions or flowering. Meanwhile, there could be other genes (e.g., VcCO1 and VcFLC) that might have eluded our comparative genomic analysis due to their high specificity or uniqueness to blueberry.

Flowering gene network of blueberry
Based on profiles of DE flowering-related genes in 'VcFT-Aurora', we propose a flowering gene network in blueberry (Fig 4). In this network, blueberry has a universal FT-mediated flowering pathway, where the photoperiod pathway (e.g., VcCOL2, and VcCOL5), and the vernalization and autonomous pathways (e.g., VcFRI, VcMAF2, VcMAF5, or VcVRN1) function through VcFT and its downstream integrators VcSOC1 and VcLFY, although no orthologues of CO1 and FLC were identified. In this universal pathway, over-expressing VcFT promotes floral bud formation and flowering. VcFT also appears to promote expression of several SPL (SQUA-MOSA PROMOTER BINDING LIKE) homologs in the aging pathway. However, VcFT overexpression is not sufficient to completely reverse the influence of environmental stimuli for normal plant flowering; for example, 'VcFT-Aurora' plants did not flower normally when grown under non-chilling conditions. It is possible that there are additional genes/pathways independent of VcFT expression that regulate vernalization-mediated blueberry flowering.
In conclusion, we developed a transcriptome reference of highbush blueberry using RNA sequencing data and Trinity. By searching this transcriptome using known flowering pathway gene sequences of Arabidopsis, rice and cereal, we developed a reference of blueberry flowering-related genes. Comparative genomics of flowering pathway genes in blueberry and herbaceous plants revealed a conserved VcFT-mediated flowering through its down-stream VcSOC1, VcAP1, VcLFY, and VcFUL genes. Interestingly, neither VcCO1 nor VcFLC are present in our transcriptome reference. The profile of DE flowering pathway genes in leaf tissues of VcFToverexpressing plants suggests VcCOL5 and VcMAF2/VcMAF5 may be the major integrator(s) in the photoperiod and vernalization pathway, respectively.
Supporting Information S1  Table. Analysis of refTrinity assembly quality and completeness by alignment with existing unigene sequences. Data from blastn of BLASTplus with e-value threshold of 1e-5 using indicated query and refTrinity as database. (DOCX) S4