Gene order and location in chromosomes of species are non-random. Neighboring gene pairs tend to display some similarities, such as co-expression and co-modification. However, the contribution of linear proximity, spatial proximity, and evolutionary proximity to these similarities remain unclear, together with whether the presence of several types of proximity can strengthens the similarities. In this study, we investigated the properties of three kinds of colocalized gene pairs: intrachromosomal colocalized gene pairs, always-neighboring gene pairs, and evolutionary neighboring gene pairs. Our analysis showed that (1) Different types of colocalized genes differentially contribute to co-expression, co-modifications and conservation across species; (2) Intrachromosomal colocalization can strengthen co-expression and co-modification of neighboring gene pairs and their conservation across species; (3) The combination of the three kinds of colocalization can lead to the strongest co-modification and is most strongly conserved across species. (4) Colocalized gene pairs are indicative of phylogenetic relationships and whole genome duplications (WGDs). These results provide valuable clues for future efforts to understand the characteristics of colocalized gene pairs and how the neighborhood affects their interactions.

Background: Theoretically, paralogous genes generated through whole genome duplications should share identical expression levels due to their identical sequences and chromatin environments. However, functional divergences and expression differences have arisen due to selective pressures throughout evolution. A comprehensive investigation of the expression patterns of paralogous gene pairs in response to various stresses and a study of correlations between the expression levels and sequence divergences of the paralogs are needed. Results: In this study, we analyzed the expression patterns of paralogous genes under different types of stress and investigated the correlations between the expression levels and sequence divergences of the paralogs. We analyzed the differential expression patterns of the paralogs under four different types of stress (drought, cold, infection, and herbivory) and classified them into three main types according to their expression patterns. We then further analyzed the differential expression patterns under various degrees of stress and constructed corresponding co-expression networks of differentially expressed paralogs and transcription factors. Finally, we investigated the correlations between the expression levels and sequence divergences of the paralogs and identified positive correlations between expression level and sequence divergence. With regard to sequence divergence, we identified correlations between selective pressures and phylogenetic relationships. Conclusions: These results shed light on differential expression patterns of paralogs in response to environmental stresses and are helpful for understanding the relationships between expression levels and sequences divergences.

Tea plants are continuously confronted with a wide range of biotic and abiotic stressors in the field, which can occur concurrently or sequentially. To elucidate the molecular mechanisms in responses to such individual and combined stresses, we used RNAseq to compare the temporal changes in the transcriptome of Camellia sinensis to Ectropis oblique Prout alone or in combination with exposure to drought and heat. Compared with the individual stress, tea plants exhibit significant differences in transcriptome profiles under the combined stresses. Additionally, many unique genes exhibited significant differences in expression in individual and combined stress conditions. Our research showed novel insights into the molecular mechanisms of E. oblique Prout resistance in tea plants and provided a valuable resource for developing tea varieties with broad spectrum stress tolerance.

Background: RNA-binding proteins (RBPs) play important roles in cellular homeostasis by regulating the expression of thousands of transcripts, which have been reported to be involved in human tumorigenesis. Despite previous reports of the dysregulation of RBPs in cancers, the degree of dysregulation of RBPs in cancers and the intrinsic relevance between dysregulated RBPs and clinical TNM information remains unknown. Furthermore, the co-expressed networks of dysregulated RBPs with transcriptional factors and lncRNAs also require further investigation. Results: Here, we firstly analyzed the deviations of expression levels of 1,542 RBPs from 20 cancer types and found that (1) RBPs are dysregulated in almost all 20 cancer types, especially in BLCA, COAD, READ, STAD, LUAD, LUSC and GBM with proportion of deviation larger than 300% compared with non-RBPs in normal tissues. (2) Up- and down-regulated RBPs also show opposed patterns of differential expression in cancers and normal tissues. In addition, down-regulated RBPs show a greater degree of dysregulated expression than up-regulated RBPs do. Secondly, we analyzed the intrinsic relevance between dysregulated RBPs and clinical TNM information and found that (3) Clinical TNM information for two cancer types—CHOL and KICH—is shown to be closely related to patterns of differentially expressed RBPs (DE RBPs) by co-expression cluster analysis. Thirdly, we identified ten key RBPs (seven down-regulated and three up-regulated) in CHOL and seven key RBPs (five down-regulated and two up-regulated) in KICH by analyzing co-expression correlation networks. Fourthly, we constructed the co-expression networks of key RBPs between 1,570 TFs and 4,147 lncRNAs for CHOL and KICH, respectively. Conclusions: These results may provide an insight into the understanding of the functions of RBPs in human carcinogenesis. Furthermore, key RBPs and the co-expressed networks offer useful information for potential prognostic biomarkers and therapeutic targets for patients with cancers at the N and M stages in two cancer types CHOL and KICH.