Abstract
The genomes and annotated genes of allotetraploid cotton Gossypium hirsutum have been extensively studied in recent years. However, the expression, regulation, and evolution of intergenic genes (ITGs) have not been completely deciphered. In this study, we identified a novel set of actively expressed ITGs in G. hirsutum cotton, through transcriptome profiling based on deep sequencing data, as well as chromatin immunoprecipitation, followed by sequencing (ChIP-seq) of histone modifications and how the ITGs evolved. Totals of 17,567 and 8249 ITGs were identified in G. hirsutum and Gossypium arboreum, respectively. The expression of ITGs in G. hirsutum was significantly higher than that in G. arboreum. Moreover, longer exons were observed in G. hirsutum ITGs. Notably, 42.3% of the ITGs from G. hirsutum were generated by the long terminal repeat (LTR) insertions, while their proportion in genic genes was 19.9%. The H3K27ac and H3K4me3 modification proportions and intensities of ITGs were equivalent to genic genes. The H3K4me1 modifications were lower in ITGs. Additionally, evolution analyses revealed that the ITGs from G. hirsutum were mainly produced around 6.6 and 1.6 million years ago (Mya), later than the pegged time for genic genes, which is 7.0 Mya. The characterization of ITGs helps to elucidate the evolution of cotton genomes and shed more light on their biological functions in the transcriptional regulation of eukaryotic genes, along with the roles of histone modifications in speciation and diversification.