Recent researches on long noncoding RNAs (lncRNAs) have expanded our horizon of gene regulation and the cellular complexity. However, the number, characteristics and expression patterns of lncRNAs remain poorly characterized and how these lncRNAs biogenesis are regulated in response to drought stress in cotton are still largely unclear. In the study, using a reproducibility-based RNA-sequencing and bioinformatics strategy to analyze the lncRNAs of 9 samples under three different environment stresses (control, drought stress and re-watering, three replications), we totally identified 10,820 lncRNAs of high-confidence through five strict steps filtration, of which 9,989 were lincRNAs, 153 were inronic lncRNAs, 678 were anti-sense lncRNAs. Coding function analysis showed 6,470 lncRNAs may have the ability to code proteins. Small RNAs precursor analysis revealed that 196 lncRNAs may be the precursors to small RNAs, most of which (35.7%, 70) were miRNAs. Expression patterns analysis showed that most of lncRNAs were expressed at a low level and most inronic lncRNAs (75.95%) had a consistent expression pattern with their adjacent protein-coding genes. Further analysis of transcriptome data uncovered that lncRNAs XLOC_063105 and XLOC_115463 probably function in regulating two adjacent coding genes CotAD_37096 and CotAD_12502, respectively. Investigations of the content of plant hormones and proteomics analysis under drought stress also complemented the prediction. We analyzed the characteristics and the expression patterns of lncRNAs under drought stress and re-watering treatment, and found lncRNAs may be likely to involve in regulating plant hormones pathway in response to drought stress.

BackgroundDNA methylation, with a cryptic role in genome stability, gene transcription and expression, is involved in the drought response process in plants, but the complex regulatory mechanism is still largely unknown.ResultsHere, we performed whole-genome bisulfite sequencing (WGBS) and identified long non-coding RNAs on cotton leaves under drought stress and re-watering treatments. We obtained 31,223 and 30,997 differentially methylated regions (representing 2.48% of the genome) after drought stress and re-watering treatments, respectively. Our data also showed that three sequence contexts, including mCpG, mCHG, mCHH, all presented a hyper-methylation pattern under drought stress and were nearly restored to normal levels after the re-watering treatment. Among all the methylation variations, asymmetric CHH methylation was the most consistent with external environments, suggesting that methylation/demethylation in a CHH context may constitute a novel epigenetic modification in response to drought stress. Combined with the targets of long non-coding RNAs, we found that long non-coding RNAs may mediate variations in methylation patterns by splicing into microRNAs. Furthermore, the many hormone-related genes with methylation variations suggested that plant hormones might be a potential mechanism in the drought response.ConclusionsFuture crop-improvement strategies may benefit by taking into account not only the DNA genetic variations in cotton varieties but also the epigenetic modifications of the genome.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-017-3681-y) contains supplementary material, which is available to authorized users.

BackgroundTrehalose (a-D-glucopyranosyl a-D-glucopyranoside) is a nonreducing disaccharide and is widely distributed in bacteria, fungi, algae, plants and invertebrates. In the study, the identification of trehalose-6-phosphate synthase (TPS) genes stress-related in cotton, and the genetic structure analysis and molecular evolution analysis of TPSs were conducted with bioinformatics methods, which could lay a foundation for further research of TPS functions in cotton.ResultsThe genome information of Gossypium raimondii (group D), G. arboreum L. (group A), and G. hirsutum L. (group AD) was used in the study. Fifty-three TPSs were identified comprising 15 genes in group D, 14 in group A, and 24 in group AD. Bioinformatics methods were used to analyze the genetic structure and molecular evolution of TPSs. Real-time PCR analysis was performed to investigate the expression patterns of gene family members. All TPS family members in cotton can be divided into two subfamilies: Class I and Class II. The similarity of the TPS sequence is high within the same species and close within their family relatives. The genetic structures of two TPS subfamily members are different, with more introns and a more complicated gene structure in Class I. There is a TPS domain(Glyco transf_20) at the N-terminal in all TPS family members and a TPP domain(Trehalose_PPase) at the C-terminal in all except GrTPS6, GhTPS4, and GhTPS9. All Class II members contain a UDP-forming domain. The responses to environmental stresses showed that stresses could induce the expression of TPSs but the expression patterns vary with different stresses.ConclusionsThe distribution of TPSs varies with different species but is relatively uniform on chromosomes. Genetic structure varies with different gene members, and expression levels vary with different stresses and exhibit tissue specificity. The upregulated genes in upland cotton TM-1 is significantly more than that in G. raimondii and G. arboreum L. Shixiya 1.Electronic supplementary materialThe online version of this article (doi:10.1186/s12863-016-0360-y) contains supplementary material, which is available to authorized users.