Hi-C Maps Research Articles

A fast and adaptive detection framework for genome-wide chromatin loop mapping from Hi-C data.

Chromatin loop identification plays an important role in molecular biology and 3D genomics research, as it constitutes a fundamental process in transcription and gene regulation. Such precise chromatin structures can be identified across genome-wide interaction matrices via Hi-C data analysis, which is essential for unraveling the intricacies of transcriptional regulation. Given the increasing number of genome-wide contact maps, derived from both in situ Hi-C and single-cell Hi-C experiments, there is a pressing need for efficient and resilient algorithms capable of processing data from diverse experiments rapidly and adaptively. Here, we propose YOLOOP, a novel detection-based framework that is different from the conventional paradigm. YOLOOP stands out for its speed, surpassing the performance of previous state-of-the-art (SOTA) chromatin loop detection methods. It achieves a 30-fold acceleration compared with classification-based methods, up to 20-fold acceleration compared with the SOTA kernel-based framework, and a fivefold acceleration compared with statistical algorithms. Furthermore, the proposed framework is capable of generalizing across various cell types, multiresolution Hi-C maps, and diverse experimental protocols. Compared with the existing paradigms, YOLOOP shows up to a 10% increase in recall and a 15% increase in F1-score, particularly noteworthy in the GM12878 cell line. YOLOOP also offers fast adaptability with straightforward fine-tuning, making it readily applicable to extremely sparse single-cell Hi-C contact maps. It maintains its exceptional speed, completing genome-wide detection at a 10 kb resolution for a single-cell contact map within 1 min and for a 900-cell-superimposed contact map within 3 min, enabling fast analysis of large-scale single-cell data.

Genome sequence of the sugarcane aphid, Melanaphis sacchari (Hemiptera: Aphididae).

The sugarcane aphid, Melanaphis sacchari (Zehntner, 1897), is an agricultural pest that causes damage to plants in the Poaceae (the grasses) family, such as sorghum and sugarcane. Here, we used Nanopore long reads and Hi-C interaction map to generate a chromosome-level assembly with a total length of 356.1 Mb, of which 85.5% (304.6 Mb) is contained within the three autosomes and the X chromosome. Repetitive sequences accounted for 16.29% of the chromosomes and a total of 12,530 protein-coding genes were annotated, achieving 95.8% benchmarking universal single-copy orthologs (BUSCO) gene completeness. This offers a substantial improvement compared to previous low-quality genomic resources. Phylogenomic analysis by comparing M. sacchari with twenty-four published aphid genomes representing three aphid tribes reveals that M. sacchari belongs to the tribe Aphidini and maintained a conserved chromosome structure with other Aphidini species. The high-quality genomic resources reported in this study will be useful for understanding the evolution of aphid genomes and studying pest management of M. sacchari.

De novo Whole-Genome Assembly of the 10-Gigabase Fokienia Hodginsii Genome to Reveal Differential Epigenetic Events Between Callus and Xylem.

Fokienia hodginsii (F. hodginsii), belonging to the genus Fokienia of the Cupressaceae. F. hodginsii has significant application value due to its wood properties and great research value in evolutionary studies as a gymnosperm. However, the genome of F. hodginsii remains unknown due to the large size of gymnosperms genome. Pacific Bioscience sequencing, Hi-C mapping, whole-genome Bisulfite Sequencing (BS-Seq), long-read isoform sequencing (Iso-Seq), direct RNA sequencing (DRS), quantitative proteomics, and metabonomics analysis are employed to facilitate genome assembly, gene annotation, and investigation into epigenetic mechanisms. In this study, the 10G F. hodginsii genome is assembled into 11 chromosomes. Furthermore, 50521 protein-coding genes are annotated and determined that 65% of F. hodginsii genome comprises repetitive sequences. It is discovered that transposable element (TE)-including introns is associated with higher expression. The DNA methylome of F. hodginsii reveals that xylem has a higher DNA methylation level compared to callus. Moreover, DRS reveals the significant alterations in RNA full-length ratio, which potentially associated with poly(A) length (PAL) and alternative polyadenylation (APA). Finally, the morphology measurement and metabonomics analysis revealed the difference of 14 cultivars. In summary, the genomes and epigenetics datasets provide a molecular basis for callus formation in the gymnosperm family.

Machine learning unravels inherent structural patterns in Escherichia coli Hi-C matrices and predicts chromosome dynamics.

High dimensional nature of the chromosomal conformation contact map ('Hi-C Map'), even for microscopically small bacterial cell, poses challenges for extracting meaningful information related to its complex organization. Here we first demonstrate that an artificial deep neural network-based machine-learnt (ML) low-dimensional representation of a recently reported Hi-C interaction map of archetypal bacteria Escherichia coli can decode crucial underlying structural pattern. The ML-derived representation of Hi-C map can automatically detect a set of spatially distinct domains across E. coli genome, sharing reminiscences of six putative macro-domains previously posited via recombination assay. Subsequently, a ML-generated model assimilates the intricate relationship between large array of Hi-C-derived chromosomal contact probabilities and respective diffusive dynamics of each individual chromosomal gene and identifies an optimal number of functionally important chromosomal contact-pairs that are majorly responsible for heterogenous, coordinate-dependent sub-diffusive motions of chromosomal loci. Finally, the ML models, trained on wild-type E. coli show-cased its predictive capabilities on mutant bacterial strains, shedding light on the structural and dynamic nuances of ΔMatP30MM and ΔMukBEF22MM chromosomes. Overall our results illuminate the power of ML techniques in unraveling the complex relationship between structure and dynamics of bacterial chromosomal loci, promising meaningful connections between ML-derived insights and biological phenomena.

Identifying topologically associating domains using differential kernels.

Chromatin is a polymer complex of DNA and proteins that regulates gene expression. The three-dimensional (3D) structure and organization of chromatin controls DNA transcription and replication. High-throughput chromatin conformation capture techniques generate Hi-C maps that can provide insight into the 3D structure of chromatin. Hi-C maps can be represented as a symmetric matrix [Formula: see text], where each element represents the average contact probability or number of contacts between chromatin loci i and j. Previous studies have detected topologically associating domains (TADs), or self-interacting regions in [Formula: see text] within which the contact probability is greater than that outside the region. Many algorithms have been developed to identify TADs within Hi-C maps. However, most TAD identification algorithms are unable to identify nested or overlapping TADs and for a given Hi-C map there is significant variation in the location and number of TADs identified by different methods. We develop a novel method to identify TADs, KerTAD, using a kernel-based technique from computer vision and image processing that is able to accurately identify nested and overlapping TADs. We benchmark this method against state-of-the-art TAD identification methods on both synthetic and experimental data sets. We find that the new method consistently has higher true positive rates (TPR) and lower false discovery rates (FDR) than all tested methods for both synthetic and manually annotated experimental Hi-C maps. The TPR for KerTAD is also largely insensitive to increasing noise and sparsity, in contrast to the other methods. We also find that KerTAD is consistent in the number and size of TADs identified across replicate experimental Hi-C maps for several organisms. Thus, KerTAD will improve automated TAD identification and enable researchers to better correlate changes in TADs to biological phenomena, such as enhancer-promoter interactions and disease states.

MDB-19. THE 3D CHROMATIN LANDSCAPE OF MEDULLOBLASTOMA

Abstract BACKGROUND Medulloblastoma (MB) is a highly malignant childhood cerebellar tumor comprised of four molecularly and clinically distinct subgroups. Sporadic yet subgroup-specific genetic and epigenetic landscapes suggest that higher-order chromatin topology may contribute to distinct subgroup identities and disease biology. Herein, we investigated the 3D chromatin landscape of MB by performing genome-wide chromosome conformation capture (Hi-C) across molecular subgroups using a cohort of n=24 patient-derived xenografts and primary tumor specimens, aiming to elucidate regulatory mechanisms facilitated by genome organization. METHODS Integrative analyses were adapted to incorporate a wide variety of sample-paired datasets. Chromatin contact matrices generated from Hi-C allowed for the mapping of active/inactive compartments (A/B compartments), topologically associated domains (TADs), and distal pairwise chromatin contacts (loops). Topological insights of chromatin architecture were complemented with whole-genome sequencing, RNA-seq, ATAC-seq, DNA methylation, CTCF and histone modification (H3K27ac/H3K27me3) ChIP-seq datasets to dissect mechanisms of gene regulation. RESULTS Unsupervised analysis revealed that chromatin compartmentalization is subgroup stratified. Gene loci contained in subgroup-unique active compartments correspondingly exhibited subgroup-specific chromatin contact patterns, adopted open chromatin conformations, and remained transcriptionally active. As an example, LIN28B resides within a Group 3-specific ‘A’ compartment where its local topology displays higher complexity compared to that of other subgroups. This consequentially facilitates a Group 3-unique loop connecting the LIN28B promoter with a distal enhancer, resulting in subgroup-specific transcriptional upregulation. Targeted displacement of CTCF that anchors and stabilizes chromatin looping at the LIN28B locus was sufficient to alter its transcriptional activity. Using similar lines of investigation, mechanistic insights into the topological determinants of known enhancer-hijacking events, including GFI1, GFI1B and PRDM6, were resolved. CONCLUSIONS Pioneering this integrative approach centered on comprehensive Hi-C maps, our study provides a high-resolution blueprint of the 3D chromatin landscape defining MB subgroups, informing the mechanistic basis of recurrent driver gene activation, and disclosing putative subgroup-specific vulnerabilities.

Learning Micro-C from Hi-C with diffusion models.

In the last few years, Micro-C has shown itself as an improved alternative to Hi-C. It replaced the restriction enzymes in Hi-C assays with micrococcal nuclease (MNase), resulting in capturing nucleosome resolution chromatin interactions. The signal-to-noise improvement of Micro-C allows it to detect more chromatin loops than high-resolution Hi-C. However, compared with massive Hi-C datasets available in the literature, there are only a limited number of Micro-C datasets. To take full advantage of these Hi-C datasets, we present HiC2MicroC, a computational method learning and then predicting Micro-C from Hi-C based on the denoising diffusion probabilistic models (DDPM). We trained our DDPM and other regression models in human foreskin fibroblast (HFFc6) cell line and evaluated these methods in six different cell types at 5-kb and 1-kb resolution. Our evaluations demonstrate that both HiC2MicroC and regression methods can markedly improve Hi-C towards Micro-C, and our DDPM-based HiC2MicroC outperforms regression in various terms. First, HiC2MicroC successfully recovers most of the Micro-C loops even those not detected in Hi-C maps. Second, a majority of the HiC2MicroC-recovered loops anchor CTCF binding sites in a convergent orientation. Third, HiC2MicroC loops share genomic and epigenetic properties with Micro-C loops, including linking promoters and enhancers, and their anchors are enriched for structural proteins (CTCF and cohesin) and histone modifications. Lastly, we find our recovered loops are also consistent with the loops identified from promoter capture Micro-C (PCMicro-C) and Chromatin Interaction Analysis by Paired-End Tag Sequencing (ChIA-PET). Overall, HiC2MicroC is an effective tool for further studying Hi-C data with Micro-C as a template. HiC2MicroC is publicly available at https://github.com/zwang-bioinformatics/HiC2MicroC/.

OA30 Primary osteoarthritis chondrocyte map of chromatin conformation reveals novel effector genes

Abstract Background/Aims Osteoarthritis is a complex disease with a huge public health burden. Genome-wide association studies (GWAS) have identified hundreds of osteoarthritis-associated sequence variants, but the effector genes underpinning these signals remain largely elusive. Understanding chromosome organisation in 3D space is essential for identifying long-range contacts between distant genomic features (e.g., between genes and regulatory elements), in a tissue-specific manner. Here, we generate the first whole genome chromosome conformation analysis (Hi-C) map of primary osteoarthritis chondrocytes and identify novel effector genes for the disease. Methods Primary chondrocytes collected from eight knee osteoarthritis patients underwent Hi-C analysis to link chromosomal structure to genomic sequence. The identified loops were then combined with osteoarthritis GWAS results to identify variants involved in gene regulation via enhancer-promoter interactions. Results We identified 36 genetic variants residing within chromatin loop anchors that are associated with 13 osteoarthritis GWAS signals. Five of these variants reside directly in enhancer regions of 3 enhancer-promoter loops newly described in chondrocytes, pointing to high-confidence effector genes: PAPPA, NSD2 and NCOA6. PAPPA is directly associated with the turnover of IGF-1 transport proteins, and IGF-1 is an important factor in the repair of damaged chondrocytes. Conclusion We have constructed the first Hi-C map of primary human chondrocytes and have made it available as a resource for the scientific community. By integrating 3D-genomics with large-scale genetic association data, we identify high-confidence osteoarthritis effector genes which enhance our understanding of disease and can serve as putative high-value novel targets. Disclosure N. Bittner: None. C. Shi: None. D. Zhao: None. J. Ding: None. L. Southam: None. D. Swift: None. P. Kreitmaier: None. M. Tutino: None. J. Hankinson: None. M. Wilkinson: None. G. Orozco: None. E. Zeggini: None.

Extensive long-range polycomb interactions and weak compartmentalization are hallmarks of human neuronal 3D genome.

Chromatin architecture regulates gene expression and shapes cellular identity, particularly in neuronal cells. Specifically, polycomb group (PcG) proteins enable establishment and maintenance of neuronal cell type by reorganizing chromatin into repressive domains that limit the expression of fate-determining genes and sustain distinct gene expression patterns in neurons. Here, we map the 3D genome architecture in neuronal and non-neuronal cells isolated from the Wernicke's area of four human brains and comprehensively analyze neuron-specific aspects of chromatin organization. We find that genome segregation into active and inactive compartments is greatly reduced in neurons compared to other brain cells. Furthermore, neuronal Hi-C maps reveal strong long-range interactions, forming a specific network of PcG-mediated contacts in neurons that is nearly absent in other brain cells. These interacting loci contain developmental transcription factors with repressed expression in neurons and other mature brain cells. But only in neurons, they are rich in bivalent promoters occupied by H3K4me3 histone modification together with H3K27me3, which points to a possible functional role of PcG contacts in neurons. Importantly, other layers of chromatin organization also exhibit a distinct structure in neurons, characterized by an increase in short-range interactions and a decrease in long-range ones.

Abstract 3551: Continued analysis of extensive data towards Genome in a Bottle benchmarks for a new tumor normal pair

Abstract Here, we detail in-progress genome-scale measurements from a variety of technologies of the first tumor normal benchmark from the Genome in a Bottle (GIAB) consortium. We created the first broadly-consented tumor cell line from a pancreatic ductal adenocarcinoma with matched normal pancreatic and duodenal tissue. Data is being collected from a large homogeneous batch of the tumor cell line and paired normal tissues. As we receive data we make it publicly available on the NCBI hosted GIAB FTP site, https://ftp.ncbi.nlm.nih.gov/ReferenceSamples/giab/data_somatic/. When complete, this dataset will include WGS measurements from Illumina, ONT, Hi-C, Bionano, single cell WGS, PacBio HiFi and Onso, and Element. Analysis and development of a variant calling benchmark using this data is ongoing in collaboration with an open working group of the GIAB consortium. Initial Hi-C and optical mapping data from the tumor cell line indicates substantial aneuploidy from translocations that cause large deletions. We found roughly 17 large inversions and translocations and 16 chromosomes with extensive loss of heterozygosity due to missing >30% of one copy, in addition to a few smaller duplications. Low coverage single cell sequencing that provides ploidy estimates across chromosomes showed that most of the large deletions appear in all or nearly all cells with some variation in a few cells. These observations are consistent with bulk WGS analyses from two batches of cells. Additionally, many of the observed large deletions correspond to deletions seen in the population of TCGA chromosomal Loss of Heterogeneity samples. Examining somatic SNVs in non-repetitive regions, we find that close to 60% occur in almost all cells in diploid regions, 30% occur in almost all cells in haploid regions, and 5% in only some cells in diploid and haploid regions, respectively. We take these results to indicate the cell line is relatively homogeneous and stable with most CNVs being deletions. As such, we plan to explore using long reads, ultralong reads, and Hi-C data to generate a near-complete genome assembly of the dominant tumor clone as well as a complete diploid assembly of the normal. With this personalized genome assembly, we will explore aligning tumor and normal reads to each haplotype of the normal to characterize somatic variants, including variants in minor clones. Building on the methods GIAB and T2T have developed to polish diploid assemblies of GIAB’s normal genomes and mosaic variant characterization, we will develop benchmarks for somatic variants against the normal assembly as well as standard references. Citation Format: Justin Wagner, Jennifer McDaniel, Gail L. Rosen, Nathan D. Olson, Vaidehi Patel, Chunlin Xiao, Andrew Liss, Justin Zook. Continued analysis of extensive data towards Genome in a Bottle benchmarks for a new tumor normal pair [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 3551.

Structural basis for the preservation of a subset of topologically associating domains in interphase chromosomes upon cohesin depletion

Compartment formation in interphase chromosomes is a result of spatial segregation between euchromatin and heterochromatin on a few megabase pairs (Mbp) scale. On the sub-Mbp scales, topologically associating domains (TADs) appear as interacting domains along the diagonal in the ensemble averaged Hi-C contact map. Hi-C experiments showed that most of the TADs vanish upon deleting cohesin, while the compartment structure is maintained, and perhaps even enhanced. However, closer inspection of the data reveals that a non-negligible fraction of TADs is preserved (P-TADs) after cohesin loss. Imaging experiments show that, at the single-cell level, TAD-like structures are present even without cohesin. To provide a structural basis for these findings, we first used polymer simulations to show that certain TADs with epigenetic switches across their boundaries survive after depletion of loops. More importantly, the three-dimensional structures show that many of the P-TADs have sharp physical boundaries. Informed by the simulations, we analyzed the Hi-C maps (with and without cohesin) in mouse liver and human colorectal carcinoma cell lines, which affirmed that epigenetic switches and physical boundaries (calculated using the predicted 3D structures using the data-driven HIPPS method that uses Hi-C as the input) explain the origin of the P-TADs. Single-cell structures display TAD-like features in the absence of cohesin that are remarkably similar to the findings in imaging experiments. Some P-TADs, with physical boundaries, are relevant to the retention of enhancer–promoter/promoter–promoter interactions. Overall, our study shows that preservation of a subset of TADs upon removing cohesin is a robust phenomenon that is valid across multiple cell lines.

Structural basis for the preservation of a subset of topologically associating domains in interphase chromosomes upon cohesin depletion.

Compartment formation in interphase chromosomes is a result of spatial segregation between euchromatin and heterochromatin on a few megabase pairs (Mbp) scale. On the sub-Mbp scales, topologically associating domains (TADs) appear as interacting domains along the diagonal in the ensemble averaged Hi-C contact map. Hi-C experiments showed that most of the TADs vanish upon deleting cohesin, while the compartment structure is maintained, and perhaps even enhanced. However, closer inspection of the data reveals that a non-negligible fraction of TADs is preserved (P-TADs) after cohesin loss. Imaging experiments show that, at the single-cell level, TAD-like structures are present even without cohesin. To provide a structural basis for these findings, we first used polymer simulations to show that certain TADs with epigenetic switches across their boundaries survive after depletion of loops. More importantly, the three-dimensional structures show that many of the P-TADs have sharp physical boundaries. Informed by the simulations, we analyzed the Hi-C maps (with and without cohesin) in mouse liver and human colorectal carcinoma cell lines, which affirmed that epigenetic switches and physical boundaries (calculated using the predicted 3D structures using the data-driven HIPPS method that uses Hi-C as the input) explain the origin of the P-TADs. Single-cell structures display TAD-like features in the absence of cohesin that are remarkably similar to the findings in imaging experiments. Some P-TADs, with physical boundaries, are relevant to the retention of enhancer-promoter/promoter-promoter interactions. Overall, our study shows that preservation of a subset of TADs upon removing cohesin is a robust phenomenon that is valid across multiple cell lines.

Primary osteoarthritis chondrocyte map of chromatin conformation reveals novel candidate effector genes

ObjectivesOsteoarthritis is a complex disease with a huge public health burden. Genome-wide association studies (GWAS) have identified hundreds of osteoarthritis-associated sequence variants, but the effector genes underpinning these signals remain...

Near telomere-to-telomere genome of the model plant Physcomitrium patens.

The model plant Physcomitrium patens has played a pivotal role in enhancing our comprehension of plant evolution and development. However, the current genome harbours numerous regions that remain unfinished and erroneous. To address these issues, we generated an assembly using Oxford Nanopore reads and Hi-C mapping. The assembly incorporates telomeric and centromeric regions, thereby establishing it as a near telomere-to-telomere genome except a region in chromosome 1 that is not fully assembled due to its highly repetitive nature. This near telomere-to-telomere genome resolves the chromosome number at 26 and provides a gap-free genome assembly as well as updated gene models to aid future studies using this model organism.

Stage-specific dynamic reorganization of genome topology shapes transcriptional neighborhoods in developing human retinal organoids

We have generated a high-resolution Hi-C map of developing human retinal organoids to elucidate spatiotemporal dynamics of genomic architecture and its relationship with gene expression patterns. We demonstrate progressive stage-specific alterations in DNA topology and correlate these changes with transcription of cell-type-restricted gene markers during retinal differentiation. Temporal Hi-C reveals a shift toward A compartment for protein-coding genes and B compartment for non-coding RNAs, displaying high and low expression, respectively. Notably, retina-enriched genes are clustered near lost boundaries of topologically associated domains (TADs), and higher-order assemblages (i.e., TAD cliques) localize in active chromatin regions with binding sites for eye-field transcription factors. These genes gain chromatin contacts at their transcription start site as organoid differentiation proceeds. Our study provides a global view of chromatin architecture dynamics associated with diversification of cell types during retinal development and serves as a foundational resource for in-depth functional investigations of retinal developmental traits.

EPCO-17. 3D GENOME STRUCTURE REGULATES TRANSCRIPTION IN PEDIATRIC HIGH-GRADE GLIOMA

Abstract INTRODUCTION Pediatric high-grade gliomas (pHGGs), including glioblastoma multiforme (GBM) and diffuse intrinsic pontine glioma (DIPG), are highly morbid brain tumors. Up to 80% of DIPGs harbor a somatic missense mutation in genes encoding Histone H3. To investigate whether the H3K27M mutant protein is associated with distinct chromatin structure affecting transcription regulation, we generated the first high-resolution Hi-C and ATAC-Seq maps of pHGG cell lines, and integrated these with tissue and cell genomic data. METHODS We generated sequencing data from patient-derived cell lines (DIPG n=6, GBM n=3, normal n=2) and frozen tissue specimens (DIPG n=1, normal brainstem n=1). Analyses included cell line RNA-Seq, ChIP-Seq (H3K27ac, H3K27me3, H3K27M) and genome-wide chromatin conformation capture (Hi-C), as well as tissue ATAC-Seq. Publicly available pediatric glioma tissue ChIP-Seq data was integrated with cell data. CRISPR knock-down of target enhancer regions was performed. RESULTS We identified tumor-specific enhancers and regulatory networks for known oncogenes in DIPG and GBM. In DIPG, FOX, SOX, STAT and SMAD families were among top H3K27Ac enriched motifs. Significant differences in Topologically Associating Domains (TADs) and DNA looping were observed at OLIG2 and MYCN in H3K27M mutant DIPG, relative to wild-type GBM and normal cells. Pharmacologic treatment targeting H3K27Ac (BET and Bromodomain inhibition) altered these 3D structures. Functional analysis of differentially enriched enhancers in DIPG implicated SOX2, SUZ12, and TRIM24 as top activated upstream regulators. Distinct genomic structural variations leading to enhancer hijacking and gene co-amplification were identified at A2M, JAG2, and FLRT1. CONCLUSION We show genome structural variations enhancer-promoter interactions that impact gene expression in pHGG in the presence and absence of the H3K27M mutation. Our results imply that tridimensional genome alterations may play a critical role in the pHGG epigenetic landscape and thereby contribute to pediatric gliomagenesis. Further studies examining the impact of the alterations is therefore underway.

11. Evaluation of Hi-C versus optical genome mapping for diagnosing constitutional genomic structural variants

11. Evaluation of Hi-C versus optical genome mapping for diagnosing constitutional genomic structural variants

The chromosome-scale reference genome of mirid bugs (Adelphocorissuturalis) genome provides insights into omnivory, insecticide resistance, and survival adaptation

BackgroundAdelphocoris suturalis (Hemiptera: Miridae) is a notorious agricultural pest, which causes serious economic losses to a diverse range of agricultural crops around the world. The poor understanding of its genomic characteristics has seriously hindered the establishment of sustainable and environment-friendly agricultural pest management through biotechnology and biological insecticides.ResultsHere, we report a chromosome-level assembled genome of A. suturalis by integrating Illumina short reads, PacBio, 10x Chromium, and Hi-C mapping technologies. The resulting 1.29 Gb assembly contains twelve chromosomal pseudomolecules with an N50 of 1.4 and 120.6 Mb for the contigs and scaffolds, respectively, and carries 20,010 protein-coding genes. The considerable size of the A. suturalis genome is predominantly attributed to a high amount of retrotransposons, especially long interspersed nuclear elements (LINEs). Transcriptomic and phylogenetic analyses suggest that A. suturalis-specific candidate effectors, and expansion and expression of gene families associated with omnivory, insecticide resistance and reproductive characteristics, such as digestion, detoxification, chemosensory receptors and long-distance migration likely contribute to its strong environmental adaptability and ability to damage crops. Additionally, 19 highly credible effector candidates were identified and transiently overexpressed in Nicotiana benthamiana for functional assays and potential targeting for insect resistance genetic engineering.ConclusionsThe high-quality genome of A. suturalis provides an important genomic landscape for further investigations into the mechanisms of omnivory, insecticide resistance and survival adaptation, and for the development of integrated management strategies.

Genome assembly of Polygala tenuifolia provides insights into its karyotype evolution and triterpenoid saponin biosynthesis.

Polygala tenuifolia is a perennial medicinal plant that has been widely used in traditional Chinese medicine for treating mental diseases. However, the lack of genomic resources limits the insight into its evolutionary and biological characterization. In the present work, we reported the P. tenuifolia genome, the first genome assembly of the Polygalaceae family. We sequenced and assembled this genome by a combination of Illumnina, PacBio HiFi, and Hi-C mapping. The assembly includes 19 pseudochromosomes covering ~92.68% of the assembled genome (~769.62Mb). There are 36 463 protein-coding genes annotated in this genome. Detailed comparative genome analysis revealed that P. tenuifolia experienced two rounds of whole genome duplication that occurred ~39-44 and ~18-20 million years ago, respectively. Accordingly, we systematically reconstructed ancestral chromosomes of P. tenuifolia and inferred its chromosome evolution trajectories from the common ancestor of core eudicots to the present species. Based on the transcriptomics data, enzyme genes and transcription factors involved in the synthesis of triterpenoid saponin in P. tenuifolia were identified. Further analysis demonstrated that whole-genome duplications and tandem duplications play critical roles in the expansion of P450 and UGT gene families, which contributed to the synthesis of triterpenoid saponins. The genome and transcriptome data will not only provide valuable resources for comparative and functional genomic researches on Polygalaceae, but also shed light on the synthesis of triterpenoid saponin.

Mapping the semi-nested community structure of 3D chromosome contact networks.

Mammalian DNA folds into 3D structures that facilitate and regulate genetic processes such as transcription, DNA repair, and epigenetics. Several insights derive from chromosome capture methods, such as Hi-C, which allow researchers to construct contact maps depicting 3D interactions among all DNA segment pairs. These maps show a complex cross-scale organization spanning megabase-pair compartments to short-ranged DNA loops. To better understand the organizing principles, several groups analyzed Hi-C data assuming a Russian-doll-like nested hierarchy where DNA regions of similar sizes merge into larger and larger structures. Apart from being a simple and appealing description, this model explains, e.g., the omnipresent chequerboard pattern seen in Hi-C maps, known as A/B compartments, and foreshadows the co-localization of some functionally similar DNA regions. However, while successful, this model is incompatible with the two competing mechanisms that seem to shape a significant part of the chromosomes' 3D organization: loop extrusion and phase separation. This paper aims to map out the chromosome's actual folding hierarchy from empirical data. To this end, we take advantage of Hi-C experiments and treat the measured DNA-DNA interactions as a weighted network. From such a network, we extract 3D communities using the generalized Louvain algorithm. This algorithm has a resolution parameter that allows us to scan seamlessly through the community size spectrum, from A/B compartments to topologically associated domains (TADs). By constructing a hierarchical tree connecting these communities, we find that chromosomes are more complex than a perfect hierarchy. Analyzing how communities nest relative to a simple folding model, we found that chromosomes exhibit a significant portion of nested and non-nested community pairs alongside considerable randomness. In addition, by examining nesting and chromatin types, we discovered that nested parts are often associated with active chromatin. These results highlight that cross-scale relationships will be essential components in models aiming to reach a deep understanding of the causal mechanisms of chromosome folding.