Strawberry is an economically grown horticulture crop required for fruit consumption. The ripening of its fruit is a complex biological process regulated by various hormones. Abscisic acid (ABA) is a critical phytohormone involved in fruit ripening. However, little is known about the long non-coding RNAs (LncRNAs), especially transposon-derived LncRNA (TE-lncRNA), response to hormones during fruit ripening in octoploid strawberry. In the study, the transcriptome data of developing strawberry fruits treated with ABA and its inhibitor Nordihydroguaiaretic acid (NGDA) were analyzed to identify responsive LncRNAs and coding genes. A total of 14,552 LncRNAs were identified, including 8,617 transposon-derived LncRNAs (TE-LncRNAs), 412 LncRNAs (282 TE-LncRNAs), and 382 ABA-sensitive LncRNAs (231 TE-LncRNAs). Additionally, a weighted co-expression network analysis constructed 27 modules containing coding RNAs and LncRNAs. Seven modules, including “MEdarkorange” and “MElightyellow” were significantly correlated with ABA/NDGA treatments, resulting in 247 hub genes, including 21 transcription factors and 22 LncRNAs (15 TE-LncRNAs). Gene ontology enrichment analysis further revealed that ABA/NDGA-responsive modules, including LncRNAs, were associated with various metabolic pathways involved in strawberry fruit development and ripening, including lipid metabolism, organic acid metabolism, and phenylpropanoid metabolism. The current study identifies many high-confidence LncRNAs in strawberry, with a percentage of them being ABA pathway-specific and 22 hub-responsive LncRNAs, providing new insight into strawberry or other Rosaceae crop fruit ripening.

Quinoa is a gluten-free food crop that contains all the essential amino acids and vitamins. The selection of proper housekeeping and tissue-specific genes is the crucial prerequisite for gene expression analysis using the common approach, real-time quantitative PCR (RT-qPCR). In this study, we identified 40 novel candidate housekeeping genes by the minimum transcript per million (TPM), coefficient of variation (CV) and maximum fold change (MFC) methods and 19 candidate tissue-specific genes by the co-expression network method based on an RNA-seq dataset that included 53 stem, leaf, flower and seed samples, as well as additional shoot and root samples under different stresses. The expression stability of 12 housekeeping and tissue-specific genes, as well as that of another two traditionally used housekeeping genes, was further evaluated using qPCR and ranked using NormFinder, BestKeeper and the comparative delta-Ct method. The results demonstrated that MIF, RGGA, VATE and UBA2B were ranked as the top four most stable candidate housekeeping genes. qPCR analysis also revealed three leaf-specific genes and five root-specific genes, but no stem-specific gene was identified. Gene Ontology (GO) enrichment analysis identified that housekeeping genes were mainly enriched in the small molecule metabolic process, organonitrogen compound metabolic process, NAD binding and ligase activity. In addition, tissue-specific genes are closely associated with the major functions of a specific tissue. Specifically, GO terms “photosynthesis” and “thylakoid” were most significantly overrepresented in candidate leaf-specific genes. The novel housekeeping and tissue-specific genes in our study will enable better normalization and quantification of transcript levels in quinoa.

Several studies have mined short-read RNA sequencing datasets to identify lncRNAs, and others have focused on the function of individual lncRNA in abiotic stress response. However, our understanding of the complement, function and origin of long-non-coding RNA (lncRNAs) response to abiotic stress, especially transposon derived lncRNA (TE-lncRNA), is still in its infancy. To discover and study lncRNAs in maize ( Zea mays ssp. mays), we utilized a dataset of 127 RNA sequencing samples that included PacBio fl-cDNA and total RNA-Seq datasets. Overall, we identified 23,309 candidate lncRNAs, 60% of which were identified in polyadenylated (polyA+) samples. The majority (65%) of the 23,309 lncRNAs had sequence similarity to transposable elements (TEs). Most had similarity to long-terminal-repeat retrotransposons from the Copia and Gypsy superfamilies, representing the high proportion of these elements in the genome, but class II, DNA transposons were enriched for lncRNAs relative to their genomic representation by 2-fold. By assessing the fraction of lncRNAs that respond to abiotic stresses like heat, cold, salt and drought, we identified 1,077 differentially expressed lncRNA transcripts. Their expression was correlated (r ² =0.48) with their nearest gene, suggesting that lncRNAs are subject to some of the cis regulatory features as neighboring genes. By inferring co-expression networks across our large dataset, we found that 39 lncRNAs act as major hubs in co-expression networks, of which 18 appeared to be derived from TEs. These results suggest that lncRNAs, especially TE-lncRNAs, may play key regulatory roles in moderating abiotic responses.