Interchromosomal Interactions Research Articles

Epigenetic control and inheritance of rDNA arrays.

Ribosomal RNA (rRNA) genes exist in multiple copies arranged in tandem arrays known as ribosomal DNA (rDNA). The total number of gene copies is variable, and the mechanisms buffering this copy number variation remain unresolved. We surveyed the number, distribution, and activity of rDNA arrays at the level of individual chromosomes across multiple human and primate genomes. Each individual possessed a unique fingerprint of copy number distribution and activity of rDNA arrays. In some cases, entire rDNA arrays were transcriptionally silent. Silent rDNA arrays showed reduced association with the nucleolus and decreased interchromosomal interactions, indicating that the nucleolar organizer function of rDNA depends on transcriptional activity. Methyl-sequencing of flow-sorted chromosomes, combined with long read sequencing, showed epigenetic modification of rDNA promoter and coding region by DNA methylation. Silent arrays were in a closed chromatin state, as indicated by the accessibility profiles derived from Fiber-seq. Removing DNA methylation restored the transcriptional activity of silent arrays. Array activity status remained stable through the iPS cell re-programming. Family trio analysis demonstrated that the inactive rDNA haplotype can be traced to one of the parental genomes, suggesting that the epigenetic state of rDNA arrays may be heritable. We propose that the dosage of rRNA genes is epigenetically regulated by DNA methylation, and these methylation patterns specify nucleolar organizer function and can propagate transgenerationally.

Ploidy variation induces butterfly effect on chromatin topology in wheat.

Polyploidy is a prominent driver of plant diversification, accompanied with dramatic chromosomal rearrangement and epigenetic changes that affect gene expression. How chromatin interactions within and between subgenomes adapt to ploidy transition remains poorly understood. We generate open chromatin interaction maps for natural hexaploid wheat (AABBDD), extracted tetraploid wheat (AABB), diploid wheat progenitor Aegilops tauschii (DD) and resynthesized hexaploid wheat (RHW, AABBDD). Thousands of intra- and interchromosomal loops are de novo established or disappeared in AB subgenomes after separation of D subgenome, in which 37-95% of novel loops are lost again in RHW after merger of D genome. Interestingly, more than half of novel loops are formed by cascade reactions that are triggered by disruption of chromatin interaction between AB and D subgenomes. The interaction repressed genes in RHW relative to DD are expression suppressed, resulting in more balanced expression of the three homoeologs in RHW. The interaction levels of cascade anchors are decreased step-by-step. Leading single nucleotide polymorphisms of yield- and plant architecture-related quantitative trait locus are significantly enriched in cascade anchors. The expression of 116 genes interacted with these anchors are significantly correlated with the corresponding traits. Our findings reveal trans-regulation of intrachromosomal loops by interchromosomal interactions during genome merger and separation in polyploid species.

Multi-omics analysis reveals the dynamic interplay between Vero host chromatin structure and function during vaccinia virus infection

The genome folds into complex configurations and structures thought to profoundly impact its function. The intricacies of this dynamic structure-function relationship are not well understood particularly in the context of viral infection. To unravel this interplay, here we provide a comprehensive investigation of simultaneous host chromatin structural (via Hi-C and ATAC-seq) and functional changes (via RNA-seq) in response to vaccinia virus infection. Over time, infection significantly impacts global and local chromatin structure by increasing long-range intra-chromosomal interactions and B compartmentalization and by decreasing chromatin accessibility and inter-chromosomal interactions. Local accessibility changes are independent of broad-scale chromatin compartment exchange (~12% of the genome), underscoring potential independent mechanisms for global and local chromatin reorganization. While infection structurally condenses the host genome, there is nearly equal bidirectional differential gene expression. Despite global weakening of intra-TAD interactions, functional changes including downregulated immunity genes are associated with alterations in local accessibility and loop domain restructuring. Therefore, chromatin accessibility and local structure profiling provide impactful predictions for host responses and may improve development of efficacious anti-viral counter measures including the optimization of vaccine design.

Preferential Co-Expression and Colocalization of rDNA-Contacting Genes with LincRNAs Suggest Their Involvement in Shaping Inter-Chromosomal Interactions with Nucleoli.

Different developmental genes shape frequent dynamic inter-chromosomal contacts with rDNA units in human and Drosophila cells. In the course of differentiation, changes in these contacts occur, coupled with changes in the expression of hundreds of rDNA-contacting genes. The data suggest a possible role of nucleoli in the global regulation of gene expression. However, the mechanism behind the specificity of these inter-chromosomal contacts, which are rebuilt in every cell cycle, is not yet known. Here, we describe the strong association of rDNA-contacting genes with numerous long intergenic non-coding RNAs (lincRNAs) in HEK293T cells and in initial and differentiated K562 cells. We observed that up to 600 different lincRNAs were preferentially co-expressed with multiple overlapping sets of rDNA-contacting developmental genes, and there was a strong correlation between the genomic positions of rDNA-contacting genes and lincRNA mappings. These two findings suggest that lincRNAs might guide the corresponding developmental genes toward rDNA clusters. We conclude that the inter-chromosomal interactions of rDNA-contacting genes with nucleoli might be guided by lincRNAs, which might physically link particular genomic regions with rDNA clusters.

DNA methylation-based high-resolution mapping of long-distance chromosomal interactions in nucleosome-depleted regions.

3C-based methods have significantly advanced our understanding of 3D genome organization. However, it remains a formidable task to precisely capture long-range chromosomal interactions between individual loci, such as those between promoters and distal enhancers. Here, we present Methyltransferase Targeting-based chromosome Architecture Capture (MTAC), a method that maps the contacts between a target site (viewpoint) and the rest of the genome in budding yeast with high resolution and sensitivity. MTAC detects hundreds of intra- and inter-chromosomal interactions within nucleosome-depleted regions (NDRs) that cannot be captured by 4C, Hi-C, or Micro-C. By applying MTAC to various viewpoints, we find that (1) most long-distance chromosomal interactions detected by MTAC reflect tethering by the nuclear pore complexes (NPCs), (2) genes co-regulated by methionine assemble into inter-chromosomal clusters near NPCs upon activation, (3) mediated by condensin, the mating locus forms a highly specific interaction with the recombination enhancer (RE) in a mating-type specific manner, and (4) correlation of MTAC signals among NDRs reveal spatial mixing and segregation of the genome. Overall, these results demonstrate MTAC as a powerful tool to resolve fine-scale long-distance chromosomal interactions and provide insights into the 3D genome organization.

Multiscale 3D genome organization underlies duck fatty liver with no adipose inflammation or serious injury

Non-alcoholic fatty liver disease (NAFLD) is the most common form of chronic liver disease. Little is known about how gene expression and chromatin structure are regulated in NAFLD due to lack of suitable model. Ducks naturally develop fatty liver similar to serious human non-alcoholic fatty liver (NAFL) without adipose inflammation and liver fibrosis, thus serves as a good model for investigating molecular mechanisms of adipose metabolism and anti-inflammation. Here, we constructed a NAFLD model without adipose inflammation and liver fibrosis in ducks. By performing dynamic pathological and transcriptomic analyses, we identified critical genes involving in regulation of the NF-κB and MHCII signaling, which usually lead to adipose inflammation and liver fibrosis. We further generated dynamic three-dimensional chromatin maps during liver fatty formation and recovery. This showed that ducks enlarged hepatocyte cell nuclei to reduce inter-chromosomal interaction, decompress chromatin structure, and alter strength of intra-TAD and loop interactions during fatty liver formation. These changes partially contributed to the tight control the NF-κB and the MHCII signaling. Our analysis uncovers duck chromatin reorganization might be advantageous to maintain liver regenerative capacity and reduce adipose inflammation. These findings shed light on new strategies for NAFLD control.

Potential roles of inter-chromosomal interactions in cell fate determination.

Mammalian genomic DNA is packed in a small nucleus, and its folding and organization in the nucleus are critical for gene regulation and cell fate determination. In interphase, chromosomes are compartmentalized into certain nuclear spaces and territories that are considered incompatible with each other. The regulation of gene expression is influenced by the epigenetic characteristics of topologically associated domains and A/B compartments within chromosomes (intrachromosomal). Previously, interactions among chromosomes detected via chromosome conformation capture-based methods were considered noise or artificial errors. However, recent studies based on newly developed ligation-independent methods have shown that inter-chromosomal interactions play important roles in gene regulation. This review summarizes the recent understanding of spatial genomic organization in mammalian interphase nuclei and discusses the potential mechanisms that determine cell identity. In addition, this review highlights the potential role of inter-chromosomal interactions in early mouse development.

The effect of trisomic chromosomes on spatial genome organization and global transcription in embryonic stem cells.

Aneuploidy frequently occurs in cancer and developmental diseases such as Down syndrome, with its functional consequences implicated in dosage effects on gene expression and global perturbation of stress response and cell proliferation pathways. However, how aneuploidy affects spatial genome organization remains less understood. In this study, we addressed this question by utilizing the previously established isogenic wild-type (WT) and trisomic mouse embryonic stem cells (mESCs). We employed a combination of Hi-C, RNA-seq, chromosome painting and nascent RNA imaging technologies to compare the spatial genome structures and gene transcription among these cells. We found that trisomy has little effect on spatial genome organization at the level of A/B compartment or topologically associating domain (TAD). Inter-chromosomal interactions are associated with chromosome regions with high gene density, active histone modifications and high transcription levels, which are confirmed by imaging. Imaging also revealed contracted chromosome volume and weakened transcriptional activity for trisomic chromosomes, suggesting potential implications for the transcriptional output of these chromosomes. Our data resources and findings may contribute to a better understanding of the consequences of aneuploidy from the angle of spatial genome organization.

The zinc-finger protein Z4 cooperates with condensin II to regulate somatic chromosome pairing and 3D chromatin organization.

Chromosome pairing constitutes an important level of genome organization, yet the mechanisms that regulate pairing in somatic cells and the impact on 3D chromatin organization are still poorly understood. Here, we address these questions in Drosophila, an organism with robust somatic pairing. In Drosophila, pairing preferentially occurs at loci consisting of numerous architectural protein binding sites (APBSs), suggesting a role of architectural proteins (APs) in pairing regulation. Amongst these, the anti-pairing function of the condensin II subunit CAP-H2 is well established. However, the factors that regulate CAP-H2 localization and action at APBSs remain largely unknown. Here, we identify two factors that control CAP-H2 occupancy at APBSs and, therefore, regulate pairing. We show that Z4, interacts with CAP-H2 and is required for its localization at APBSs. We also show that hyperosmotic cellular stress induces fast and reversible unpairing in a Z4/CAP-H2 dependent manner. Moreover, by combining the opposite effects of Z4 depletion and osmostress, we show that pairing correlates with the strength of intrachromosomal 3D interactions, such as active (A) compartment interactions, intragenic gene-loops, and polycomb (Pc)-mediated chromatin loops. Altogether, our results reveal new players in CAP-H2-mediated pairing regulation and the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions.

Spatially coordinated heterochromatinization of long synaptic genes in fragile X syndrome

Short tandem repeat (STR) instability causes transcriptional silencing in several repeat expansion disorders. In fragile X syndrome (FXS), mutation-length expansion of a CGG STR represses FMR1 via local DNA methylation. Here, we find megabase-scale H3K9me3 domains on autosomes and encompassing FMR1 on the X chromosome in FXS patient-derived iPSCs, iPSC-derived neural progenitors, EBV-transformed lymphoblasts, and brain tissue with mutation-length CGG expansion. H3K9me3 domains connect via inter-chromosomal interactions and demarcate severe misfolding of TADs and loops. They harbor long synaptic genes replicating at the end of S phase, replication-stress-induced double-strand breaks, and STRs prone to stepwise somatic instability. CRISPR engineering of the mutation-length CGG to premutation length reverses H3K9me3 on the X chromosome and multiple autosomes, refolds TADs, and restores gene expression. H3K9me3 domains can also arise in normal-length iPSCs created with perturbations linked to genome instability, suggesting their relevance beyond FXS. Our results reveal Mb-scale heterochromatinization and trans interactions among loci susceptible to instability.

Interchromosomal Colocalization with Parental Genes Is Linked to the Function and Evolution of Mammalian Retrocopies.

Retrocopies are gene duplicates arising from reverse transcription of mature mRNA transcripts and their insertion back into the genome. While long being regarded as processed pseudogenes, more and more functional retrocopies have been discovered. How the stripped-down retrocopies recover expression capability and become functional paralogs continually intrigues evolutionary biologists. Here, we investigated the function and evolution of retrocopies in the context of 3D genome organization. By mapping retrocopy-parent pairs onto sequencing-based and imaging-based chromatin contact maps in human and mouse cell lines and onto Hi-C interaction maps in 5 other mammals, we found that retrocopies and their parental genes show a higher-than-expected interchromosomal colocalization frequency. The spatial interactions between retrocopies and parental genes occur frequently at loci in active subcompartments and near nuclear speckles. Accordingly, colocalized retrocopies are more actively transcribed and translated and are more evolutionarily conserved than noncolocalized ones. The active transcription of colocalized retrocopies may result from their permissive epigenetic environment and shared regulatory elements with parental genes. Population genetic analysis of retroposed gene copy number variants in human populations revealed that retrocopy insertions are not entirely random in regard to interchromosomal interactions and that colocalized retroposed gene copy number variants are more likely to reach high frequencies, suggesting that both insertion bias and natural selection contribute to the colocalization of retrocopy-parent pairs. Further dissection implies that reduced selection efficacy, rather than positive selection, contributes to the elevated allele frequency of colocalized retroposed gene copy number variants. Overall, our results hint a role of interchromosomal colocalization in the "resurrection" of initially neutral retrocopies.

Histone deacetylation and cytosine methylation compartmentalize heterochromatic regions in the genome organization of Neurospora crassa.

Chromosomes must correctly fold in eukaryotic nuclei for proper genome function. Eukaryotic organisms hierarchically organize their genomes, including in the fungus Neurospora crassa, where chromatin fiber loops compact into Topologically Associated Domain-like structures formed by heterochromatic region aggregation. However, insufficient data exist on how histone posttranslational modifications (PTMs), including acetylation, affect genome organization. In Neurospora, the HCHC complex [composed of the proteins HDA-1, CDP-2 (Chromodomain Protein-2), Heterochromatin Protein-1, and CHAP (CDP-2 and HDA-1 Associated Protein)] deacetylates heterochromatic nucleosomes, as loss of individual HCHC members increases centromeric acetylation, and alters the methylation of cytosines in DNA. Here, we assess whether the HCHC complex affects genome organization by performing Hi-C in strains deleted of the cdp-2 or chap genes. CDP-2 loss increases intra- and interchromosomal heterochromatic region interactions, while loss of CHAP decreases heterochromatic region compaction. Individual HCHC mutants exhibit different patterns of histone PTMs genome-wide, as CDP-2 deletion increases heterochromatic H4K16 acetylation, yet smaller heterochromatic regions lose H3K9 trimethylation and gain interheterochromatic region interactions; CHAP loss produces minimal acetylation changes but increases heterochromatic H3K9me3 enrichment. Loss of both CDP-2 and the DIM-2 DNA methyltransferase causes extensive genome disorder as heterochromatic-euchromatic contacts increase despite additional H3K9me3 enrichment. Our results highlight how the increased cytosine methylation in HCHC mutants ensures genome compartmentalization when heterochromatic regions become hyperacetylated without HDAC activity.

Three-dimensional chromatin landscapes in hepatocellular carcinoma associated with hepatitis B virus.

Three-dimensional (3D) chromatin architecture frequently altered in cancer. However, its changes during the pathogenesis of hepatocellular carcinoma (HCC) remained elusive. Hi-C and RNA-seq were applied to study the 3D chromatin landscapes and gene expression of HCC and ANHT. Hi-C Pro was used to generate genome-wide raw interaction matrices, which were normalized via iterative correction (ICE). Moreover, the chromosomes were divided into different compartments according to the first principal component (E1). Furthermore, topologically associated domains (TADs) were visualized via WashU Epigenome Browser. Furthermore, differential expression analysis of ANHT and HCC was performed using the DESeq2 R package. Additionally, dysregulated genes associated with 3D genome architecture altered were confirmed using TCGA, qRT-PCR, immunohistochemistry (IHC), etc. RESULTS: First, the intrachromosomal interactions of chr1, chr2, chr5, and chr11 were significantly different, and the interchromosomal interactions of chr4-chr10, chr13-chr21, chr15-chr22, and chr16-chr19 are remarkably different between ANHT and HCC, which resulted in the up-regulation of TP53I3 and ZNF738 and the down-regulation of APOC3 and APOA5 in HCC. Second, 49 compartment regions on 18 chromosomes have significantly switched (A-B or B-A) during HCC tumorigenesis, contributing to up-regulation of RAP2A. Finally, a tumor-specific TAD boundary located on chr5: 6271000-6478000 and enhancer hijacking were identified in HCC tissues, potentially associated with the elevated expression of MED10, whose expression were associated with poor prognosis of HCC patients. This study demonstrates the crucial role of chromosomal structure variation in HCC oncogenesis and potential novel biomarkers of HCC, laying a foundation for cancer precision medicine development.

Chromatin architecture of two different strains of Artemisia annua reveals the alterations in interaction and gene expression.

The hierarchical architecture of chromatins affects the gene expression level of glandular secreting trichomes and the artemisinin biosynthetic pathway-related genes, consequently bringing on huge differences in the content of artemisinin and its derivatives of A. annua. The plant of traditional Chinese medicine "Qinghao" is called Artemisia annua L. in Chinese Pharmacopoeia. High content and the total amount of artemisinin is the main goal of A. annua breeding, nevertheless, the change of chromatin organization during the artemisinin synthesis process has not been discovered yet. This study intended to find the roles of chromatin structure in the production of artemisinin through bioinformatics and experimental validation. Chromosome conformation capture analysis was used to scrutinize the interactions among chromosomes and categorize various scales of chromatin during artemisinin synthesis in A. annua. To confirm the effect of the changes in chromatin structure, Hi-C and RNA-sequencing were performed on two different strains to find the correlation between chromatin structure and gene expression levels on artemisinin synthesis progress and regulation. Our results revealed that the frequency of intra-chromosomal interactions was higher in the inter-chromosomal interactions between the root and leaves on a high artemisinin production strain (HAP) compared to a low artemisinin production strain (LAP). We found that compartmental transition was connected with interactions among different chromatins. Interestingly, glandular secreting trichomes (GSTs) and the artemisinin biosynthetic pathway (ABP) related genes were enriched in the areas which have the compartmental transition, reflecting the regulation of artemisinin synthesis. Topologically associated domain boundaries were associated with various distributions of genes and expression levels. Genes associated with ABP and GST in the adjacent loop were highly expressed, suggesting that epigenetic regulation plays an important role during artemisinin synthesis and glandular secreting trichomes production process. Chromatin structure could show an important status in the mechanisms of artemisinin synthesis process in A. annua.

Intra- and interchromosomal contact mapping reveals the Igh locus has extensive conformational heterogeneity and interacts with B-lineage genes

To produce a diverse antibody repertoire, immunoglobulin heavy-chain (Igh) loci undergo large-scale alterations in structure to facilitate juxtaposition and recombination of spatially separated variable (VH), diversity (DH), and joining (JH) genes. These chromosomal alterations are poorly understood. Uncovering their patterns shows how chromosome dynamics underpins antibody diversity. Using tiled Capture Hi-C, we produce a comprehensive map of chromatin interactions throughout the 2.8-Mb Igh locus in progenitor B cells. We find that the Igh locus folds into semi-rigid subdomains and undergoes flexible looping of the VH genes to its 3' end, reconciling two views of locus organization. Deconvolution of single Igh locus conformations using polymer simulations identifies thousands of different structures. This heterogeneity may underpin the diversity of V(D)J recombination events. All three immunoglobulin loci also participate in a highly specific, developmentally regulated network of interchromosomal interactions with genes encoding B cell-lineage factors. This suggests a model of interchromosomal coordination of B cell development.

Noncoding SNP at rs1663689 represses ADGRG6 via interchromosomal interaction and reduces lung cancer progression.

A previous genome-wide association study (GWAS) revealed an association of the noncoding SNP rs1663689 with susceptibility to lung cancer in the Chinese population. However, the underlying mechanism is unknown. In this study, using allele-specific 4C-seq in heterozygous lung cancer cells combined with epigenetic information from CRISPR/Cas9-edited cell lines, we show that the rs1663689 C/C variant represses the expression of ADGRG6, a gene located on a separate chromosome, through an interchromosomal interaction of the rs1663689 bearing region with the ADGRG6 promoter. This reduces downstream cAMP-PKA signaling and subsequently tumor growth both in vitro and in xenograft models. Using patient-derived organoids, we show that rs1663689 T/T-but not C/C-bearing lung tumors are sensitive to the PKA inhibitor H89, potentially informing therapeutic strategies. Our study identifies a genetic variant-mediated interchromosomal interaction underlying ADGRG6 regulation and suggests that targeting the cAMP-PKA signaling pathway may be beneficial in lung cancer patients bearing the homozygous risk genotype at rs1663689.

Probabilistic establishment of speckle-associated inter-chromosomal interactions.

Inter-chromosomal interactions play a crucial role in genome organization, yet the organizational principles remain elusive. Here, we introduce a novel computational method to systematically characterize inter-chromosomal interactions using in situ Hi-C results from various cell types. Our method successfully identifies two apparently hub-like inter-chromosomal contacts associated with nuclear speckles and nucleoli, respectively. Interestingly, we discover that nuclear speckle-associated inter-chromosomal interactions are highly cell-type invariant with a marked enrichment of cell-type common super-enhancers (CSEs). Validation using DNA Oligopaint fluorescence in situ hybridization (FISH) shows a strong but probabilistic interaction behavior between nuclear speckles and CSE-harboring genomic regions. Strikingly, we find that the likelihood of speckle-CSE associations can accurately predict two experimentally measured inter-chromosomal contacts from Hi-C and Oligopaint DNA FISH. Our probabilistic establishment model well describes the hub-like structure observed at the population level as a cumulative effect of summing individual stochastic chromatin-speckle interactions. Lastly, we observe that CSEs are highly co-occupied by MAZ binding and MAZ depletion leads to significant disorganization of speckle-associated inter-chromosomal contacts. Taken together, our results propose a simple organizational principle of inter-chromosomal interactions mediated by MAZ-occupied CSEs.

EccDB: a comprehensive repository for eccDNA-mediated chromatin contacts in multi-species.

We developed the eccDB database to integrate available resources for extrachromosomal circular DNA (eccDNA) data. eccDB is a comprehensive repository for storing, browsing, searching, and analyzing eccDNAs from multispecies. The database provides regulatory and epigenetic information on eccDNAs, with a focus on analyzing intrachromosomal and interchromosomal interactions to predict their transcriptional regulatory functions. Moreover, eccDB identifies eccDNAs from unknown DNA sequences and analyzes the functional and evolutionary relationships of eccDNAs among different species. Overall, eccDB offers web-based analytical tools and a comprehensive resource for biologists and clinicians to decipher the molecular regulatory mechanisms of eccDNAs. eccDB is freely available at http://www.xiejjlab.bio/eccDB.

Long range inter-chromosomal interaction of Oct4 distal enhancer loci regulates ESCs pluripotency

Nuclear architecture underlies the transcriptional programs within the cell to establish cell identity. As previously demonstrated, long-range chromatin interactions of the Oct4 distal enhancer (DE) are correlated with active transcription in naïve state embryonic stem cells. Here, we identify and characterize extreme long-range interactions of the Oct4 DE through a novel CRISPR labeling technique we developed and chromosome conformation capture to identify lethal giant larvae 2 (Llgl2) and growth factor receptor-bound protein 7 (Grb7) as putative functional interacting target genes in different chromosomes. We show that the Oct4 DE directly regulates expression of Llgl2 and Grb7 in addition to Oct4. Expression of Llgl2 and Grb7 closely correlates with the pluripotent state, where knock down of either result in loss of pluripotency, and overexpression enhances somatic cell reprogramming. We demonstrated that biologically important interactions of the Oct4 DE can occur at extreme distances that are necessary for the maintenance of the pluripotent state.

Investigating human chromosome organization by whole-genome simulations.

Both theoreticians and experimentalists have been dedicating significant effort to studying chromosome organization and its relation to gene expression and regulation. Many experimental procedures were designed, such as Hi-C or ChIA-PET, to capture structural information of chromosomes. Several theoretical models were proposed to understand the principles underlying chromosome compartmentalization. MiChroM is a model that describes the compartmentalization behavior of interphase chromosomes. Compared to experimental data, the agreement of the in silico contact maps suggested that phase separation of A/B compartments is a key feature determining the structure of individual chromosomes. In this work, we use the Open-MiChroM software to simulate the 46 chromosomes of a GM12878 human cell at 50 kb resolution, seeking to understand whether the nuclear organization of the whole genome would rely on the same principles observed for the organization of a single chromosome. The ensemble of structures generated by our simulations captures chromosome organization within territories and phase separation of chromatin compartments. Inside a chromosome territory, active chromatin localizes preferentially in the periphery, while inactive chromatin lies on the inside. The intra-chromosomal contact maps and the shape metrics of each territory are consistent with experimental data. The resulting genome architecture resembles the inverted nucleus configuration, with active chromatin towards the periphery. In the interface between chromosomes, we also observe phase separation of compartments. However, the inter-chromosomal maps present higher contact probabilities for heterochromatin loci when compared to the experimental data. These results suggest that the interactions of chromatin with other nuclear features like the lamina greatly influences inter-chromosomal interactions while also affecting the overall positioning of chromosomes in the nucleus.