Insight into high-resolution three-dimensional genome organization and its effect on transcription remains largely elusive in plants. Here, using a long-read ChIA-PET approach, we map H3K4me3- and RNA polymerase II (RNAPII)-associated promoter–promoter interactions and H3K9me2-marked heterochromatin interactions at nucleotide/gene resolution in rice. The chromatin architecture is separated into different independent spatial interacting modules with distinct transcriptional potential and covers approximately 82% of the genome. Compared to inactive modules, active modules possess the majority of active loop genes with higher density and contribute to most of the transcriptional activity in rice. In addition, promoter–promoter interacting genes tend to be transcribed cooperatively. In contrast, the heterochromatin-mediated loops form relative stable structure domains in chromatin configuration. Furthermore, we examine the impact of genetic variation on chromatin interactions and transcription and identify a spatial correlation between the genetic regulation of eQTLs and e-traits. Thus, our results reveal hierarchical and modular 3D genome architecture for transcriptional regulation in rice.

Autopolyploidy is widespread in higher plants and important for agricultural yield and quality. However, the effects of genome duplication on the chromatin organization and transcriptional regulation are largely unknown in plants. Using High-throughput Chromosome Conformation Capture (Hi-C), we showed that autotetraploid Arabidopsis presented more inter-chromosomal interactions and fewer short-range chromatin interactions compared with its diploid progenitor. In addition, genome duplication contributed to the switching of some loose and compact structure domains with altered H3K4me3 and H3K27me3 histone modification status. 539 genes were identified with altered transcriptions and chromatin interactions in autotetraploid Arabidopsis. Especially, we found that genome duplication changed chromatin looping and H3K27me3 histone modification in Flowering Locus C. We propose that genome doubling modulates the transcription genome-wide by changed chromatin interactions and at the specific locus by altered chromatin loops and histone modifications.

The complexity of the epigenome landscape and transcriptional regulation is significantly increased during plant polyploidization, which drives genome evolution and contributes to the increased adaptability to diverse environments. However, a comprehensive epigenome map of Brassica napus is still unavailable. In this study, we performed integrative analysis of five histone modifications, RNA polymerase II occupancy, DNA methylation, and transcriptomes in two B. napus lines (2063A and B409), and established global maps of regulatory elements, chromatin states, and their dynamics for the whole genome (including the An and Cn subgenomes) in four tissue types (young leaf, flower bud, silique, and root) of these two lines. Approximately 65.8% of the genome was annotated with different epigenomic signals. Compared with the Cn subgenome, the An subgenome possesses a higher level of active epigenetic marks and lower level of repressive epigenetic marks. Genes from subgenome-unique regions contribute to the major differences between the An and Cn subgenomes. Asymmetric histone modifications between homeologous gene pairs reflect their biased expression patterns. We identified a novel bivalent chromatin state (with H3K4me1 and H3K27me3) in B. napus that is associated with tissue-specific gene expression. Furthermore, we observed that different types of duplicated genes have discrepant patterns of histone modification and DNA methylation levels. Collectively, our findings provide a valuable epigenetic resource for allopolyploid plants.

Single-cell sequencing opens a new era for the investigation of tumor immune microenvironments (TIME). However, at single-cell resolution, a pan-cancer analysis that addresses the identity and diversity of TIMEs is lacking. Here, we first built a pan-cancer single-cell reference of TIMEs with refined subcell types and recognized new cell type–specific transcription factors. We then presented a pan-cancer view of the common features of the TIME and compared the variation of each immune cell type across patients and tumor types in the aspects of abundance, cell states, and cell communications. We found that the abundance and the cell states of dysfunctional T cells were most variable, whereas those of regulatory T cells were relatively stable. A subset of tumor-associated macrophages (TAM), PLTP+C1QC+ TAMs, may regulate the abundance of dysfunctional T cells through cytokine/chemokine signaling. The ligand–receptor communication network of TIMEs was tumor-type specific and dominated by the tumor-enriched immune cells. We additionally developed the single-cell TIME (scTIME) portal (http://scTIME.sklehabc.com) with the scTIME-specific analysis modules and a unified cell annotation. In addition to the immune cell compositions and correlation analysis using refined cell type classifications, the portal also provides cell–cell interaction and cell type–specific gene signature analysis. Our single-cell pan-cancer analysis and scTIME portal will provide more insights into the features of TIMEs, as well as the molecular and cellular mechanisms underlying immunotherapies.