Transcriptome-Based Identification of Candidate Flowering-Associated Genes of Blueberry in a Plant Factory with Artificial Lighting (PFAL) under Short-Day-Length Conditions.

In blueberry (Vaccinium corymbosum L.), a perennial shrub, flower bud initiation is mediated by a short-day (SD) photoperiod and buds bloom once the chilling requirement is satisfied. A plant factory with artificial lighting (PFAL) is a planting system that can provide a stable and highly efficient growing environment for blueberry production. However, the characteristics of bud differentiation of blueberry plants inside PFAL systems are poorly understood. To better understand flower bud initiation and the flowering mechanism of blueberry in PFAL systems, the anatomical structure of apical buds under SD conditions in a PFAL system was observed using the southern highbush cultivar 'Misty' and a transcriptomic analysis was performed to identify the candidate flowering genes. The results indicated that the apical bud of 'Misty' differentiated gradually along with SD time course and swelled obviously when chilling was introduced. A total of 39.28 Gb clean data were generated, and about 20.31-24.11 Mb high-quality clean reads were assembled, yielding a total of 17370 differentially expressed genes (DEGs), of which 9637 were up-regulated and 7733 were down-regulated. Based on the functional annotation, 26 DEGs were identified including 20 flowering-related and 6 low-temperature DEGs, out of which the expressive level of four flowering-related DEGs (VcFT2, VcFPA, VcFMADS1, and VcCOP1) and two low-temperature-induced DEGs (VcTIL-1 and VcLTI 65-like) were confirmed by qRT-PCR with a good consistency with the pattern of transcriptome. Functional analysis indicated that VcFT2 was highly conserved with nuclear and cytoplasmic subcellular localization and was expressed mainly in blueberry leaf tissue. In Arabidopsis, ectopic overexpression of VcFT2 results in an early flowering phenotype, indicating that VcFT2 is a vital regulator of the SD-mediated flowering pathway in blueberry. These results contribute to the investigation of photoperiod-mediated flowering mechanisms of blueberry in PFAL systems.

Integrated transcriptome and metabolome analysis reveals the anthocyanin biosynthesis mechanisms in blueberry (Vaccinium corymbosum L.) leaves under different light qualities.

Introduction: Blueberry (Vaccinium corymbosum L.) is a popular fruit with an abundance of anthocyanins in its leaves and fruits. Light is one of the pivotal environmental elements that affects plant growth and development, but the regulatory mechanism between light quality and anthocyanin formation is poorly understood. Methods: An integrated transcriptome and metabolome analysis was performed to investigate the effects of white (control), blue (B), red (R), and red/blue (60R/40B) light on blueberry growth and reveal the potential pathway controlling anthocyanin biosynthesis in blueberry leaves. Results: The anthocyanin content was significantly improved by the blue and red/blue light when compared with white light, whereas there was a significant reduction in the photosynthesis under the blue light, showing an inverse trend to that of anthocyanin accumulation. Transcriptomic analysis resulted in the assembly of 134,709 unigenes. Of these, 22 were differentially expressed genes (DEGs) that participate in the anthocyanin biosynthesis pathway, with the majority being significantly up-regulated under the blue light. Most of the photosynthesis-related genes that were down-regulated were expressed during anthocyanin accumulation. Targeted metabolome profiling identified 44 metabolites associated with anthocyanin biosynthesis. The contents of most of these metabolites were higher under blue light than the other light conditions, which was consistent with the transcriptome results. The integrated transcriptome and metabolome analysis suggested that, under blue light, leucoanthocyanidin dioxygenase (LDOX), O-methyltransferase (OMT), and UDP-glucose flavonoid glucosyltransferase (UFGT) were the most significantly expressed, and they promoted the synthesis of cyanidin (Cy), malvidin (Mv), and pelargonidin (Pg) anthocyanidins, respectively. The expression levels of dihydroflavonol 4-reductase (DFR) and OMT, as well as the accumulation of delphinidin (Dp), peonidin (Pn), and petunidin (Pt), were significantly increased by the red/blue light. Discussion: The blue and red/blue lights promoted anthocyanin biosynthesis via inducing the expression of key structural genes and accumulation of metabolites involved in anthocyanin synthesis pathway. Moreover, there was a possible feedback regulating correlation between anthocyanin biosynthesis and photosynthesis under different light qualities in blueberry leaves. This study would provide a theoretical basis for elucidating the underlying regulatory mechanism of anthocyanin biosynthesis of V. corymbosum.

Genome-wide identification and comprehensive analysis reveal potential roles of long non-coding RNAs in fruit development of southern highbush blueberry (Vaccinium corymbosum L.).

Introduction: Blueberries have a high antioxidant content and are produced as healthy food worldwide. Long non-coding RNAs (lncRNAs) are a type of regulatory RNAs that play a variety of roles in plants. Nonetheless, information on lncRNAs and their functions during blueberry fruit development is scarce in public databases. Methods: In the present study, we performed genome-wide identification of lncRNAs in a southern highbush blueberry using strand-specific RNA sequencing (ssRNA-Seq). Differentially expressed lncRNAs (DE-lncRNAs) and their potential target genes were analyzed at four stages of fruit development. Cis-regulatory DE-lncRNAs were predicted using co-localization analysis. Results: These findings included a total of 25,036 lncRNAs from 17,801 loci. Blueberry lncRNAs had shorter transcript lengths, smaller open reading frame (ORF) sizes, fewer exons, and fewer isoforms than protein-coding RNAs, as well as lower expression levels and higher stage-specificity during fruit development. A total of 105 DE-lncRNAs were identified among the comparison group of PAD vs. CUP, 443 DE-lncRNAs were detected when comparing CUP with PINK fruits, and 285 DE-lncRNAs were revealed when comparing PINK and BLUE fruits. According to Kyoto Encyclopedia of Genes and Genomes annotation, target genes of DE-lncRNAs were primarily enriched in the "Autophagy-other", "DNA replication", "Endocytosis", 'photosynthesis' and 'chlorophyll metabolism' pathways, suggesting that lncRNAs may pay potential roles in fruit expansion and ripening. Moreover, several lncRNAs have been proposed as cis-regulators of the key genes involved in flavonoid biosynthesis. MSTRG.107242.6, and its putative target gene, BTB/POZ and TAZ domain-containing protein, might play critical roles in anthocyanin accumulation in blueberries. Discussion: These findings highlight the regulatory function of lncRNAs and aid in elucidating the molecular mechanism underlying blueberry fruit growth.

Transcriptomic Analysis Reveals Potential Gene Regulatory Networks Under Cold Stress of Loquat (Eriobotrya japonica Lindl.).

Loquat (Eriobotrya japonica Lindl. ) is one of the most economically important evergreen fruit crops in China, while it often suffered the injury of cold stress in winter and earlier spring, and the annual yield loss of loquat fruits caused by cold or freezing stress was immeasurable. However, knowledge about the physiological response and molecular mechanism under cold stress is still limited. To investigate the potential regulation mechanism pre- and post-cold stress in loquat and the changes in physiological indicators, a comparative transcriptome analysis was performed against a cold-resistant cv. "Huoju" and a cold-sensitive cv. "Ninghaibai". The results of physiological indicators related to cold resistance indicated that rachis was most sensitive to cold stress and was considered as the representative organ to directly evaluate cold resistance of loquat based on subordinate function analysis. Here, we compared the transcriptome profiles of rachis pre- and under cold stress in "Huoju" and "Ninghaibai". A total of 4,347 and 3,513 differentially expressed genes (DEGs) were detected in "Ninghaibai" and "Huoju", among which 223 and 166 were newly identified genes, respectively, most of them were functionally enriched in plant hormone signal transduction (Huoju: 142; Ninghaibai: 200), and there were higher plant hormone content and related DEG expression levels in "Huoju" than that of "Ninghaibai". Moreover, a total of 3,309 differentially expressed transcription factors (DETFs) were identified, and some DEGs and DETFs were screened to be subjected to co-expression network analysis based on the gene expression profile data. Some candidate DEGs, including UDP-glycosyltransferase (UGT), glycosyltransferase (GT), sugar phosphate/phosphate translocator (SPT), sugar transport protein (STP), proline-rich receptor-like protein kinase (PERK), and peroxidise (POD), were significantly affected by cold stress, and the expression level of these genes obtained from real-time quantitative RT-PCR was consistent with the pattern of transcriptome profile, which suggested that these genes might play the vital roles in cold resistance of loquat. Our results provide an invaluable resource for the identification of specific genes and TFs and help to clarify gene transcription during the cold stress response of loquat.

Genome-wide analysis of LysM gene family members and their expression in response to Colletotrichum fructicola infection in Octoploid strawberry(Fragaria × ananassa).

The cultivated octoploid strawberry (Fragaria × ananassa) is an economically important fruit that is planted worldwide. The lysin motif (LysM) protein family is composed of the major class of plant pattern recognition receptors, which play important roles in sensing pathogen-associated molecular patterns (PAMPs), and subsequently triggers downstream plant immunity. In the present study, a comprehensive, genome-wide analysis of F. × ananassa LysM (FaLysM) genes was performed to investigate gene structures, phylogenic relationships, chromosome location, collinear relationships, transcription factor binding sites, and protein model analysis. We aimed to identify the LysM genes involved in the defense against plant pathogens. A total of 14 FaLysM genes were identified in the F. × ananassa genome and divided into 2 subgroups (LYP and LYK) on the basis of the phylogenetic analysis. The Ka/Ks ratio for the duplicated pair of most FaLysM genes was less than 1, which indicates that the selection pressure was mostly subject to the purifying selection during evolution. The protein model analysis revealed that FaLysM2-10 contain conserved mode of chitin binding, which suggest the potential role of FaLysM2-10 in pathogen perception and plant immunity. The RNA-Seq results showed the differential regulation of 14 FaLysM genes in response to Colletotrichum fructicola infection, implying the complex interaction between C. fructicola and strawberry. Knockout of candidate effector gene CfLysM2, which was previously proved to be highly expressed during C. fructicola infection, resulted in the up-regulation of six FaLysM genes (FaLysM1, FaLysM2, FaLysM3, FaLysM7, FaLysM8, and FaLysM12), indicating the competitive relations between CfLysM2 and FaLysM genes. Overall, this study provides fundamental information on the roles of LysM proteins in octoploid strawberry and its interaction with C. fructicola, laying useful information for further investigation on the C. fructicola-strawberry interaction and strawberry resistance breeding.

Transcriptomic profiling and discovery of key genes involved in adventitious root formation from green cuttings of highbush blueberry (Vaccinium corymbosum L.).

BACKGROUND: Propagation of cuttings is frequently used in various plant species, including blueberry, which shows special root characteristics that may hinder adventitious root (AR) formation. AR formation is influenced by various factors, and auxin is considered to play a central role; however, little is known of the related regulatory mechanisms. In this study, a comparative transcriptome analysis of green cuttings treated with or without indole-butyric acid (IBA) was performed via RNA_seq to identify candidate genes associated with IBA-induced AR formation. RESULTS: Rooting phenotypes, especially the rooting rate, were significantly promoted by exogenous auxin in the IBA application. Blueberry AR formation was an auxin-induced process, during which adventitious root primordium initiation (rpi) began at 14 days after cutting (DAC), root primordium (rp) was developed at 21 DAC, mature AR was observed at 28 DAC and finally outgrowth from the stem occurred at 35 DAC. Higher IAA levels and lower ABA and zeatin contents might facilitate AR formation and development. A time series transcriptome analysis identified 14,970 differentially expressed genes (DEGs) during AR formation, of which there were 7467 upregulated and 7503 downregulated genes. Of these, approximately 35 candidate DEGs involved in the auxin-induced pathway and AR formation were further identified, including 10 auxin respective genes (ARFs and SAURs), 13 transcription factors (LOB domain-containing protein (LBDs)), 6 auxin transporters (AUX22, LAX3/5 and PIN-like 6 (PIL6s)) and 6 rooting-associated genes (root meristem growth factor 9 (RGF9), lateral root primordium 1 (LRP1s), and dormancy-associated protein homologue 3 (DRMH3)). All these identified DEGs were highly upregulated in certain stages during AR formation, indicating their potential roles in blueberry AR formation. CONCLUSIONS: The transcriptome profiling results indicated candidate genes or major regulatory factors that influence adventitious root formation in blueberry and provided a comprehensive understanding of the rooting mechanism underlying the auxin-induced AR formation from blueberry green cuttings.

Chromosome-level genome assembly and annotation of the loquat (Eriobotrya japonica) genome.

BACKGROUND: The loquat (Eriobotrya japonica) is a species of flowering plant in the family Rosaceae that is widely cultivated in Asian, European, and African countries. It blossoms in the winter and ripens in the early summer. The genome of loquat has to date not been published, which limits the study of molecular biology in this cultivated species. Here, we used the third-generation sequencing technology of Nanopore and Hi-C technology to sequence the genome of Eriobotrya. FINDINGS: We generated 100.10 Gb of long reads using Oxford Nanopore sequencing technologies. Three types of Illumina high-throughput sequencing data, including genome short reads (47.42 Gb), transcriptome short reads (11.06 Gb), and Hi-C short reads (67.25 Gb), were also generated to help construct the loquat genome. All data were assembled into a 760.1-Mb genome assembly. The contigs were mapped to chromosomes by using Hi-C technology based on the contacts between contigs, and then a genome was assembled exhibiting 17 chromosomes and a scaffold N50 length of 39.7 Mb. A total of 45,743 protein-coding genes were annotated in the Eriobotrya genome, and we investigated the phylogenetic relationships between the Eriobotrya and 6 other Rosaceae species. Eriobotrya shows a close relationship with Malus and Pyrus, with the divergence time of Eriobotrya and Malus being 6.76 million years ago. Furthermore, chromosome rearrangement was found in Eriobotrya and Malus. CONCLUSIONS: We constructed the first high-quality chromosome-level Eriobotrya genome using Illumina, Nanopore, and Hi-C technologies. This work provides a valuable reference genome for molecular studies of the loquat and provides new insight into chromosome evolution in this species.

Proteome analysis provides new insight into major proteins involved in gibberellin-induced fruit setting in triploid loquat (Eriobotrya japonica).

BACKGROUND: Parthenocarpy can be induced by gibberellin (GA) treatment in plants. The fruits of the loquat exhibit many seeds. GA treatment can induce the development of seedless fruit and increase fruit quality during production. However, the molecular mechanism of fruit setting under GA treatment is still unclear. OBJECTIVE: Our aim was to explore GA-induced parthenocarpy in triploid loquat by proteome analysis to identify the differentially expressed proteins. METHODS: A proteome analysis was performed using TMT protein labeling and LC-MS/MS in triploid loquat. RESULTS: A total of 7290 protein groups were identified in the two stages of fruit setting. The quantitative results showed that 923 differentially expressed proteins (DEPs) were isolated, which were enriched in five pathways: ribosome, citrate cycle (TCA cycle), pentose phosphate, carbon metabolism, and carbon fixation. Twenty-four DEPs were annotated as putative key regulatory proteins involved in fruit setting, which were related to the auxin response, gibberellin metabolism, ethylene synthesis, and cell division. In addition, thirty-five DEPs were involved in the formation of the cell wall, which might be downstream proteins involved in cell proliferation during fruit setting. CONCLUSION: Our report reveals new insight into the protein dynamics of loquat fruit setting induced by GA treatment via the analysis of proteome profiles and provides a reference for other Rosaceae species.

Comparative Analysis of Transcriptomes to Identify Genes Associated with Fruit Size in the Early Stage of Fruit Development in Pyrus pyrifolia.

Pear (Pyrus L.) is an important commercial fruit in the world. The fruit size is one of the important characters in fruit quality. The previous research reported that the fruit size of pear was mainly caused by the number of cell in about 40 days after blossom (DAB) in nature. However, studies about the mechanisms underlying cell division in young fruit development are very limited in pear. Two pear accessions codenamed 'GH59B' with big fruit and 'GH81S' with small fruit in three stages were sampled and the RNA-seq high-throughput sequencing was used to evaluate changes of gene transcription levels in the early stage of fruit development. The difference of cell size among two samples was little in 40 DAB, implying that the difference of the fruit size was caused by the number of the cell. More than 274,517,982 high quality reads from six libraries of fruit development were sequenced. A total of 797 differentially expressed genes (DEGs) were identified. Three cytokinin dehydrogenase genes and two gibberellin 2-beta-dioxygenase gene were identified in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to zeatin and gibberellin. Their expression was upregulated at 20 DAB in 'GH81S' and at 30 DAB in 'GH59B', suggesting that the small fruit size might be related to the early degradation of cytokinin and gibberellin inducing a short period of cell division. A total of 38 DEGs of transcription factors were found and 23 DEGs including NAC, ERF and bHLH transcription factors were highly related with cytokinin dehydrogenase and gibberellin dioxygenase genes. Altogether, the results of the present study provide information from a comprehensive gene expression analysis and insight into the molecular mechanism underlying the difference of fruit size in Pyrus pyrifolia.