Hi-C 2.0: An optimized Hi-C procedure for high-resolution genome-wide mapping of chromosome conformation

Polymer physics models reveal structural folding features of single-molecule gene chromatin conformations

Here, we employ polymer physics models of chromatin to investigate the 3D folding of a 2Mb wide genomic region encompassing the human LTN1 gene, a crucial DNA locus involved in key cellular functions. Through extensive Molecular Dynamics simulations, we reconstruct in-silico the ensemble of single-molecule LTN1 3D structures, which we benchmark against recent in-situ Hi-C 2.0 data. The model-derived single molecules are then used to predict structural folding features at the single-cell level, providing testable predictions for super-resolution microscopy experiments.

Connecting Chromatin Structures to Gene Regulation Using Dynamic Polymer Simulations

The transfer of regulatory information between distal loci on chromatin is thought to involve physical proximity, but key biophysical features of these contacts remain unclear. For instance, it is unknown how close and for how long two loci need to be in order to productively interact. The main challenge is that it is currently impossible to measure chromatin dynamics with high spatiotemporal resolution at scale. Polymer simulations provide an accessible and rigorous way to test biophysical models of chromatin regulation, yet there is a lack of simple and general methods for extracting the values of model parameters. Here we adapt the Nelder-Mead simplex optimization algorithm to select the best polymer model matching a given Hi-C dataset, using the MYC locus as an example. The model's biophysical parameters predict a compartmental rearrangement of the MYC locus in leukemia, which we validate with single-cell measurements. Leveraging trajectories predicted by the model, we find that loci with similar Hi-C contact frequencies can exhibit widely different contact dynamics. Interestingly, the frequency of productive interactions between loci exhibits a non-linear relationship with their Hi-C contact frequency when we enforce a specific capture radius and contact duration. These observations are consistent with recent experimental observations and suggest that the dynamic ensemble of chromatin configurations, rather than average contact matrices, is required to fully predict productive long-range chromatin interactions.

Regulating epithelial-mesenchymal plasticity from 3D genome organization

Epithelial-mesenchymal transition (EMT) is a dynamic process enabling polarized epithelial cells to acquire mesenchymal features implicated in development and carcinoma progression. As our understanding evolves, it is clear the reversible execution of EMT arises from complex epigenomic regulation involving histone modifications and 3-dimensional (3D) genome structural changes, leading to a cascade of transcriptional events. This review summarizes current knowledge on chromatin organization in EMT, with a focus on hierarchical structures of the 3D genome and chromatin accessibility changes.

Genome wide nucleosome landscape shapes 3D chromatin organization

The hierarchical chromatin organization begins with formation of nucleosomes, which fold into chromatin domains punctuated by boundaries and ultimately chromosomes. In a hierarchal organization, lower levels shape higher levels. However, the dependence of higher-order 3D chromatin organization on the nucleosome-level organization has not been studied in cells. We investigated the relationship between nucleosome-level organization and higher-order chromatin organization by perturbing nucleosomes across the genome by deleting Imitation SWItch (ISWI) and Chromodomain Helicase DNA-binding (CHD1) chromatin remodeling factors in budding yeast. We find that changes in nucleosome-level properties are accompanied by changes in 3D chromatin organization. Short-range chromatin contacts up to a few kilo-base pairs decrease, chromatin domains weaken, and boundary strength decreases. Boundary strength scales with accessibility and moderately with width of nucleosome-depleted region. Change in nucleosome positioning seems to alter the stiffness of chromatin, which can affect formation of chromatin contacts. Our results suggest a biomechanical "bottom-up" mechanism by which nucleosome distribution across genome shapes 3D chromatin organization.

TAD border deletion at the Kit locus causes tissue-specific ectopic activation of a neighboring gene

Topologically associated domains (TADs) restrict promoter-enhancer interactions, thereby maintaining the spatiotemporal pattern of gene activity. However, rearrangements of the TADs boundaries do not always lead to significant changes in the activity pattern. Here, we investigated the consequences of the TAD boundaries deletion on the expression of developmentally important genes encoding tyrosine kinase receptors: Kit, Kdr, Pdgfra. We used genome editing in mice to delete the TADs boundaries at the Kit locus and characterized chromatin folding and gene expression in pure cultures of fibroblasts, mast cells, and melanocytes. We found that although Kit is highly active in both mast cells and melanocytes, deletion of the TAD boundary between the Kit and Kdr genes results in ectopic activation only in melanocytes. Thus, the epigenetic landscape, namely the mutual arrangement of enhancers and actively transcribing genes, is important for predicting the consequences of the TAD boundaries removal. We also found that mice without a TAD border between the Kit and Kdr genes have a phenotypic manifestation of the mutation - a lighter coloration. Thus, the data obtained shed light on the principles of interaction between the 3D chromatin organization and epigenetic marks in the regulation of gene activity.

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.

Rules of engagement for condensins and cohesins guide mitotic chromosome formation

During mitosis, interphase chromatin is rapidly converted into rod-shaped mitotic chromosomes. Using Hi-C, imaging, proteomics and polymer modeling, we determine how the activity and interplay between loop-extruding SMC motors accomplishes this dramatic transition. Our work reveals rules of engagement for SMC complexes that are critical for allowing cells to refold interphase chromatin into mitotic chromosomes. We find that condensin disassembles interphase chromatin loop organization by evicting or displacing extrusive cohesin. In contrast, condensin bypasses cohesive cohesins, thereby maintaining sister chromatid cohesion while separating the sisters. Studies of mitotic chromosomes formed by cohesin, condensin II and condensin I alone or in combination allow us to develop new models of mitotic chromosome conformation. In these models, loops are consecutive and not overlapping, implying that condensins do not freely pass one another but stall upon encountering each other. The dynamics of Hi-C interactions and chromosome morphology reveal that during prophase loops are extruded in vivo at ~1-3 kb/sec by condensins as they form a disordered discontinuous helical scaffold within individual chromatids.

Mitotic chromosomes are self-entangled and disentangle through a topoisomerase-II-dependent two-stage exit from mitosis

The topological state of chromosomes determines their mechanical properties, dynamics, and function. Recent work indicated that interphase chromosomes are largely free of entanglements. Here, we use Hi-C, polymer simulations, and multi-contact 3C and find that, by contrast, mitotic chromosomes are self-entangled. We explore how a mitotic self-entangled state is converted into an unentangled interphase state during mitotic exit. Most mitotic entanglements are removed during anaphase/telophase, with remaining ones removed during early G1, in a topoisomerase-II-dependent process. Polymer models suggest a two-stage disentanglement pathway: first, decondensation of mitotic chromosomes with remaining condensin loops produces entropic forces that bias topoisomerase II activity toward decatenation. At the second stage, the loops are released, and the formation of new entanglements is prevented by lower topoisomerase II activity, allowing the establishment of unentangled and territorial G1 chromosomes. When mitotic entanglements are not removed in experiments and models, a normal interphase state cannot be acquired.

The chromosome-level genome assembly of the giant dobsonfly Acanthacorydalis orientalis (McLachlan, 1899)

Acanthacorydalis orientalis (McLachlan, 1899) (Megaloptera: Corydalidae) is an important freshwater-benthic invertebrate species that serves as an indicator for water-quality biomonitoring and is valuable for conservation from East Asia. Here, a high-quality reference genome for A. orientalis was constructed using Oxford Nanopore sequencing and High throughput Chromosome Conformation Capture (Hi-C) technology. The final genome size is 547.98 Mb, with the N50 values of contig and scaffold being 7.77 Mb and 50.53 Mb, respectively. The longest contig and scaffold are 20.57 Mb and 62.26 Mb in length, respectively. There are 99.75% contigs anchored onto 13 pseudo-chromosomes. Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis showed that the completeness of the genome assembly is 99.01%. There are 10,977 protein-coding genes identified, of which 84.00% are functionally annotated. The genome contains 44.86% repeat sequences. This high-quality genome provides substantial data for future studies on population genetics, aquatic adaptation, and evolution of Megaloptera and other related insect groups.

Towards Development of the 4C-Based Method Detecting Interactions of Plasmid DNA with Host Genome

Chromosome conformation capture techniques have revolutionized our understanding of chromatin architecture and dynamics at the genome-wide scale. In recent years, these methods have been applied to a diverse array of species, revealing fundamental principles of chromosomal organization. However, structural organization of the extrachromosomal entities, like viral genomes or plasmids, and their interactions with the host genome, remain relatively underexplored. In this work, we introduce an enhanced 4C-protocol tailored for probing plasmid DNA interactions. We design specific plasmid vector and optimize protocol to allow high detection rate of contacts between the plasmid and host DNA.

Modification of the Hi-C Technology for Molecular Genetic Analysis of Formalin-Fixed Paraffin-Embedded Sections of Tumor Tissues

Molecular genetic analysis of tumor tissues is the most important step towards understanding the mechanisms of cancer development; it is also necessary for the choice of targeted therapy. The Hi-C (high-throughput chromatin conformation capture) technology can be used to detect various types of genomic variants, including balanced chromosomal rearrangements, such as inversions and translocations. We propose a modification of the Hi-C method for the analysis of chromatin contacts in formalin-fixed paraffin-embedded (FFPE) sections of tumor tissues. The developed protocol allows to generate high-quality Hi-C data and detect all types of chromosomal rearrangements. We have analyzed various databases to compile a comprehensive list of translocations that hold clinical importance for the targeted therapy selection. The practical value of molecular genetic testing is its ability to influence the treatment strategies and to provide prognostic insights. Detecting specific chromosomal rearrangements can guide the choice of the targeted therapies, which is a critical aspect of personalized medicine in oncology.

High resolution maps of chromatin reorganization through mouse meiosis reveal novel features of the 3D meiotic structure

When germ cells transition from the mitotic cycle into meiotic prophase I (MPI), chromosomes condense into an array of chromatin loops that are required to promote homolog pairing and genetic recombination. To identify the changes in chromosomal conformation, we isolated nuclei on a trajectory from spermatogonia to the end of MPI. At each stage along this trajectory, we built genomic interaction maps with the highest temporal and spatial resolution to date. The changes in chromatin folding coincided with a concurrent decline in mitotic cohesion and a rise in meiotic cohesin complexes. We found that the stereotypical large-scale A and B compartmentalization was lost during meiotic prophase I alongside the loss of topological associating domains (TADs). Still, local subcompartments were detected and maintained throughout meiosis. The enhanced Micro-C resolution revealed that, despite the loss of TADs, higher frequency contact sites between two loci were detectable during meiotic prophase I coinciding with CTCF bound sites. The pattern of interactions around these CTCF sites with their neighboring loci showed that CTCF sites were often anchoring the meiotic loops. Additionally, the localization of CTCF to the meiotic axes indicated that these anchors were at the base of loops. Strikingly, even in the face of the dramatic reconfiguration of interphase chromatin into a condensed loop-array, the interactions between regulatory elements remained well preserved. This establishes a potential mechanism for how the meiotic chromatin maintains active transcription within a highly structured genome. In summary, the high temporal and spatial resolution of these data revealed previously unappreciated aspects of mammalian meiotic chromatin organization.

Machine and deep learning methods for predicting 3D genome organization

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

Matrin3 mediates differentiation through stabilizing chromatin loop-domain interactions and YY1 mediated enhancer-promoter interactions

Although emerging evidence indicates that alterations in proteins within nuclear compartments elicit changes in chromosomal architecture and differentiation, the underlying mechanisms are not well understood. Here we investigate the direct role of the abundant nuclear complex protein Matrin3 (Matr3) in chromatin architecture and development in the context of myogenesis. Using an acute targeted protein degradation platform (dTAG-Matr3), we reveal the dynamics of development-related chromatin reorganization. High-throughput chromosome conformation capture (Hi-C) experiments revealed substantial chromatin loop rearrangements soon after Matr3 depletion. Notably, YY1 binding was detected, accompanied by the emergence of novel YY1-mediated enhancer-promoter loops, which occurred concurrently with changes in histone modifications and chromatin-level binding patterns. Changes in chromatin occupancy by Matr3 also correlated with these alterations. Overall, our results suggest that Matr3 mediates differentiation through stabilizing chromatin accessibility and chromatin loop-domain interactions, and highlight a conserved and direct role for Matr3 in maintenance of chromosomal architecture.

High-quality Momordica balsamina genome elucidates its potential use in improving stress resilience and therapeutic properties of bitter gourd

Introduction

Momordica balsamina is the closest wild species that can be crossed with an important fruit vegetable crop, Momordica charantia, has immense medicinal value, and placed under II subclass of primary gene pool of bitter gourd. M. balsamina is tolerant to major biotic and abiotic stresses. Genome characterization of Momordica balsamina as a wild relative of bitter gourd will contribute to the knowledge of the gene pool available for improvement in bitter gourd. There is potential to transfer gene/s related to biotic resistance and medicinal importance from M. balsamina to M. charantia to produce high-quality, better yielding and stress tolerant bitter gourd genotypes.

Methods

The present study provides the first and high-quality chromosome-level genome assembly of M. balsamina with size 384.90 Mb and N50 30.96 Mb using sequence data from 10x Genomics, Nanopore, and Hi-C platforms.

Results

A total of 6,32,098 transposons elements; 2,15,379 simple sequence repeats; 5,67,483 transcription factor binding sites; 3,376 noncoding RNA genes; and 41,652 protein-coding genes were identified, and 4,347 disease resistance, 67 heat stress-related, 05 carotenoid-related, 15 salt stress-related, 229 cucurbitacin-related, 19 terpenes-related, 37 antioxidant activity, and 06 sex determination-related genes were characterized.

Conclusion

Genome sequencing of M. balsamina will facilitate interspecific introgression of desirable traits. This information is cataloged in the form of webgenomic resource available at http://webtom.cabgrid.res.in/mbger/. Our finding of comparative genome analysis will be useful to get insights into the patterns and processes associated with genome evolution and to uncover functional regions of cucurbit genomes.

Mapping nucleosome-resolution chromatin organization and enhancer-promoter loops in plants using Micro-C-XL

Although chromatin organizations in plants have been dissected at the scales of compartments and topologically associating domain (TAD)-like domains, there remains a gap in resolving fine-scale structures. Here, we use Micro-C-XL, a high-throughput chromosome conformation capture (Hi-C)-based technology that involves micrococcal nuclease (instead of restriction enzymes) and long cross-linkers, to dissect single nucleosome-resolution chromatin organization in Arabidopsis. Insulation analysis reveals more than 14,000 boundaries, which mostly include chromatin accessibility, epigenetic modifications, and transcription factors. Micro-C-XL reveals associations between RNA Pols and local chromatin organizations, suggesting that gene transcription substantially contributes to the establishment of local chromatin domains. By perturbing Pol II both genetically and chemically at the gene level, we confirm its function in regulating chromatin organization. Visible loops and stripes are assigned to super-enhancers and their targeted genes, thus providing direct insights for the identification and mechanistic analysis of distal CREs and their working modes in plants. We further investigate possible factors regulating these chromatin loops. Subsequently, we expand Micro-C-XL to soybean and rice. In summary, we use Micro-C-XL for analyses of plants, which reveal fine-scale chromatin organization and enhancer-promoter loops and provide insights regarding three-dimensional genomes in plants.

A distant global control region is essential for normal expression of anterior HOXA genes during mouse and human craniofacial development

Craniofacial abnormalities account for approximately one third of birth defects. The regulatory programs that build the face require precisely controlled spatiotemporal gene expression, achieved through tissue-specific enhancers. Clusters of coactivated enhancers and their target genes, known as superenhancers, are important in determining cell identity but have been largely unexplored in development. In this study we identified superenhancer regions unique to human embryonic craniofacial tissue. To demonstrate the importance of such regions in craniofacial development and disease, we focused on an ~600 kb noncoding region located between NPVF and NFE2L3. We identified long range interactions with this region in both human and mouse embryonic craniofacial tissue with the anterior portion of the HOXA gene cluster. Mice lacking this superenhancer exhibit perinatal lethality, and present with highly penetrant skull defects and orofacial clefts phenocopying Hoxa2-/- mice. Moreover, we identified two cases of de novo copy number changes of the superenhancer in humans both with severe craniofacial abnormalities. This evidence suggests we have identified a critical noncoding locus control region that specifically regulates anterior HOXA genes and copy number changes are pathogenic in human patients.

Combinatorial single-cell profiling of major chromatin types with MAbID

Gene expression programs result from the collective activity of numerous regulatory factors. Studying their cooperative mode of action is imperative to understand gene regulation, but simultaneously measuring these factors within one sample has been challenging. Here we introduce Multiplexing Antibodies by barcode Identification (MAbID), a method for combinatorial genomic profiling of histone modifications and chromatin-binding proteins. MAbID employs antibody-DNA conjugates to integrate barcodes at the genomic location of the epitope, enabling combined incubation of multiple antibodies to reveal the distributions of many epigenetic markers simultaneously. We used MAbID to profile major chromatin types and multiplexed measurements without loss of individual data quality. Moreover, we obtained joint measurements of six epitopes in single cells of mouse bone marrow and during mouse in vitro differentiation, capturing associated changes in multifactorial chromatin states. Thus, MAbID holds the potential to gain unique insights into the interplay between gene regulatory mechanisms, especially for low-input samples and in single cells.

Hi-C2B: Optimized Detection of Chromosomal Contacts Within Synchronized Meiotic S. cerevisiae Cells

Hi-C, a genome-wide chromosome conformation capture assay, is a powerful tool used to study three-dimensional genome organization by converting physical pairwise interactions into counts of pairwise interactions. To study the many temporally regulated facets of meiotic recombination in S. cerevisiae, the Hi-C assay must be robust such that fine- and wide-scale comparisons between genetic datasets can be made. Here we describe an updated protocol for Hi-C (Hi-C2B) that generates reproducible libraries of interaction data with low noise and for a relatively low cost.

Chromosome-scale genomes of commercially important mahoganies, Swietenia macrophylla and Khaya senegalensis

Mahogany species (family Meliaceae) are highly valued for their aesthetic and durable wood. Despite their economic and ecological importance, genomic resources for mahogany species are limited, hindering genetic improvement and conservation efforts. Here we perform chromosome-scale genome assemblies of two commercially important mahogany species: Swietenia macrophylla and Khaya senegalensis. By combining 10X sequencing and Hi-C data, we assemble high-quality genomes of 274.49 Mb (S. macrophylla) and 406.50 Mb (K. senegalensis), with scaffold N50 lengths of 8.51 Mb and 7.85 Mb, respectively. A total of 99.38% and 98.05% of the assembled sequences are anchored to 28 pseudo-chromosomes in S. macrophylla and K. senegalensis, respectively. We predict 34,129 and 31,908 protein-coding genes in S. macrophylla and K. senegalensis, respectively, of which 97.44% and 98.49% are functionally annotated. The chromosome-scale genome assemblies of these mahogany species could serve as a vital genetic resource, especially in understanding the properties of non-model woody plants. These high-quality genomes could support the development of molecular markers for breeding programs, conservation efforts, and the sustainable management of these valuable forest resources.

Chromosome-level genome assembly of the deep-sea snail Phymorhynchus buccinoides provides insights into the adaptation to the cold seep habitat

BACKGROUND

The deep-sea snail Phymorhynchus buccinoides belongs to the genus Phymorhynchus (Neogastropoda: Raphitomidae), and it is a dominant specie in the cold seep habitat. As the environment of the cold seep is characterized by darkness, hypoxia and high concentrations of toxic substances such as hydrogen sulfide (H2S), exploration of the diverse fauna living around cold seeps will help to uncover the adaptive mechanisms to this unique habitat. In the present study, a chromosome-level genome of P. buccinoides was constructed and a series of genomic and transcriptomic analyses were conducted to explore its molecular adaptation mechanisms to the cold seep environments.

RESULTS

The assembled genome size of the P. buccinoides was approximately 2.1 Gb, which is larger than most of the reported snail genomes, possibly due to the high proportion of repetitive elements. About 92.0% of the assembled base pairs of contigs were anchored to 34 pseudo-chromosomes with a scaffold N50 size of 60.0 Mb. Compared with relative specie in the shallow water, the glutamate regulative and related genes were expanded in P. buccinoides, which contributes to the acclimation to hypoxia and coldness. Besides, the relatively high mRNA expression levels of the olfactory/chemosensory genes in osphradium indicate that P. buccinoides might have evolved a highly developed and sensitive olfactory organ for its orientation and predation. Moreover, the genome and transcriptome analyses demonstrate that P. buccinoides has evolved a sulfite-tolerance mechanism by performing H2S detoxification. Many genes involved in H2S detoxification were highly expressed in ctenidium and hepatopancreas, suggesting that these tissues might be critical for H2S detoxification and sulfite tolerance.

CONCLUSIONS

In summary, our report of this chromosome-level deep-sea snail genome provides a comprehensive genomic basis for the understanding of the adaptation strategy of P. buccinoides to the extreme environment at the deep-sea cold seeps.

Choreography of lamina-associated domains: structure meets dynamics

Lamina-associated domains are large regions of heterochromatin positioned at the nuclear periphery. These domains have been implicated in gene repression, especially in the context of development. In mammals, LAD organization is dependent on nuclear lamins, inner nuclear membrane proteins, and chromatin state. In addition, chromatin readers and modifier proteins have been implicated in this organization, potentially serving as molecular tethers that interact with both nuclear envelope proteins and chromatin. More recent studies have focused on teasing apart the rules that govern dynamic LAD organization and how LAD organization, in turn, relates to gene regulation and overall 3D genome organization. This review highlights recent studies in mammalian cells uncovering factors that instruct the choreography of LAD organization, re-organization, and dynamics at the nuclear lamina, including LAD dynamics in interphase and through mitotic exit, when LAD organization is re-established, as well as intra-LAD subdomain variations.

Genome and transcriptome of Ips nitidus provide insights into high-altitude hypoxia adaptation and symbiosis

Ips nitidus is a well-known conifer pest that has contributed significantly to spruce forest disturbance in the Qinghai-Tibet Plateau and seriously threatens the ecological balance of these areas. We report a chromosome-level genome of I. nitidus determined by PacBio and Hi-C technology. Phylogenetic inference showed that it diverged from the common ancestor of I. typographus ∼2.27 mya. Gene family expansion in I. nitidus was characterized by DNA damage repair and energy metabolism, which may facilitate adaptation to high-altitude hypoxia. Interestingly, differential gene expression analysis revealed upregulated genes associated with high-altitude hypoxia adaptation and downregulated genes associated with detoxification after feeding and tunneling in fungal symbiont Ophiostoma bicolor-colonized substrates. Our findings provide evidence of the potential adaptability of I. nitidus to conifer host, high-altitude hypoxia and insight into how fungal symbiont assist in this process. This study enhances our understanding of insect adaptation, symbiosis, and pest management.

Three-Dimensional Gene Regulation Network in Glioblastoma Ferroptosis

Ferroptosis is an iron-dependent form of cell death, which is reported to be associated with glioma progression and drug sensitivity. Targeting ferroptosis is a potential therapeutic approach for glioma. However, the molecular mechanism of glioma cell ferroptosis is not clear. In this study, we profile the change of 3D chromatin structure in glioblastoma ferroptosis by using HiChIP and study the 3D gene regulation network in glioblastoma ferroptosis. A combination of an analysis of HiChIP and RNA-seq data suggests that change of chromatin loops mediated by 3D chromatin structure regulates gene expressions in glioblastoma ferroptosis. Genes that are regulated by 3D chromatin structures include genes that were reported to function in ferroptosis, like HDM2 and TXNRD1. We propose a new regulatory mechanism governing glioblastoma cell ferroptosis by 3D chromatin structure.

Chromosome-level genome assembly of an important wolfberry fruit fly (Neoceratitis asiatica Becker)

Tephritidae pests are significant agricultural pests with a notable impact on the economy, with a wide range of species and most having broad host ranges and strong reproductive abilities. However, the wolfberry fruit fly, Neoceratitis asiatica (Becker), is a Tephritidae fly which only harms wolfberry. Here, we assembled and annotated N. asiatica genome at the chromosome level and compared it with the genomic and transcriptomic information from other Tephritidae flies. The assembled genome of N. asiatica had a size of 563.8 Mb and achieved a completeness level of 99.1%, 18,387 genes were annotated totally. All contigs were assembled into 7 linkage groups with an N50 of 93.166 Mb assisted by the Hi-C technique. The high-quality genome developed here will provide a significant resource for exploring the genetic basis of the adaptive and reproductive differences among various Tephritidae pests, and provides an important theoretical basis for the prevention and control of Tephritidae pests.

The dynamic role of cohesin in maintaining human genome architecture

Recent advances in genomic and imaging techniques have revealed the complex manner of organizing billions of base pairs of DNA necessary for maintaining their functionality and ensuring the proper expression of genetic information. The SMC proteins and cohesin complex primarily contribute to forming higher-order chromatin structures, such as chromosomal territories, compartments, topologically associating domains (TADs) and chromatin loops anchored by CCCTC-binding factor (CTCF) protein or other genome organizers. Cohesin plays a fundamental role in chromatin organization, gene expression and regulation. This review aims to describe the current understanding of the dynamic nature of the cohesin-DNA complex and its dependence on cohesin for genome maintenance. We discuss the current 3C technique and numerous bioinformatics pipelines used to comprehend structural genomics and epigenetics focusing on the analysis of Cohesin-centred interactions. We also incorporate our present comprehension of Loop Extrusion (LE) and insights from stochastic modelling.

Pan-genome analysis sheds light on structural variation-based dissection of agronomic traits in melon crops

Sweetness and appearance of fresh fruits are key palatable and preference attributes for consumers and are often controlled by multiple genes. However, fine-mapping the key loci or genes of interest by single genome-based genetic analysis is challenging. Herein, we present the chromosome-level genome assembly of 1 landrace melon accession (Cucumis melo ssp. agrestis) with wild morphologic features and thus construct a melon pan-genome atlas via integrating sequenced melon genome datasets. Our comparative genomic analysis reveals a total of 3.4 million genetic variations, of which the presence/absence variations (PAVs) are mainly involved in regulating the function of genes for sucrose metabolism during melon domestication and improvement. We further resolved several loci that are accountable for sucrose contents, flesh color, rind stripe, and suture using a structural variation (SV)-based genome-wide association study. Furthermore, via bulked segregation analysis (BSA)-seq and map-based cloning, we uncovered that a single gene, (CmPIRL6), determines the edible or inedible characteristics of melon fruit exocarp. These findings provide important melon pan-genome information and provide a powerful toolkit for future pan-genome-informed cultivar breeding of melon.

Range-wide differential adaptation and genomic offset in critically endangered Asian rosewoods

In the billion-dollar global illegal wildlife trade, rosewoods have been the world's most trafficked wild product since 2005. Dalbergia cochinchinensis and Dalbergia oliveri are the most sought-after rosewoods in the Greater Mekong Subregion. They are exposed to significant genetic risks and the lack of knowledge on their adaptability limits the effectiveness of conservation efforts. Here, we present genome assemblies and range-wide genomic scans of adaptive variation, together with predictions of genomic offset to climate change. Adaptive genomic variation was differentially associated with temperature and precipitation-related variables between the species, although their natural ranges overlap. The findings are consistent with differences in pioneering ability and in drought tolerance. We predict their genomic offsets will increase over time and with increasing carbon emission pathway but at a faster pace in D. cochinchinensis than in D. oliveri. These results and the distinct gene-environment association in the eastern coastal edge of Vietnam suggest species-specific conservation actions: germplasm representation across the range in D. cochinchinensis and focused on hotspots of genomic offset in D. oliveri. We translated our genomic models into a seed source matching application, seedeR, to rapidly inform restoration efforts. Our ecological genomic research uncovering contrasting selection forces acting in sympatric rosewoods is of relevance to conserving tropical trees globally and combating risks from climate change.

Fundamental insights into the correlation between chromosome configuration and transcription

Eukaryotic chromosomes exhibit a hierarchical organization that spans a spectrum of length scales, ranging from sub-regions known as loops, which typically comprise hundreds of base pairs, to much larger chromosome territories that can encompass a few mega base pairs. Chromosome conformation capture experiments that involve high-throughput sequencing methods combined with microscopy techniques have enabled a new understanding of inter- and intra-chromosomal interactions with unprecedented details. This information also provides mechanistic insights on the relationship between genome architecture and gene expression. In this article, we review the recent findings on three-dimensional interactions among chromosomes at the compartment, topologically associating domain, and loop levels and the impact of these interactions on the transcription process. We also discuss current understanding of various biophysical processes involved in multi-layer structural organization of chromosomes. Then, we discuss the relationships between gene expression and genome structure from perturbative genome-wide association studies. Furthermore, for a better understanding of how chromosome architecture and function are linked, we emphasize the role of epigenetic modifications in the regulation of gene expression. Such an understanding of the relationship between genome architecture and gene expression can provide a new perspective on the range of potential future discoveries and therapeutic research.

Female naïve human pluripotent stem cells carry X chromosomes with Xa-like and Xi-like folding conformations

Three-dimensional (3D) genomics shows immense promise for studying X chromosome inactivation (XCI) by interrogating changes to the X chromosomes' 3D states. Here, we sought to characterize the 3D state of the X chromosome in naïve and primed human pluripotent stem cells (hPSCs). Using chromatin tracing, we analyzed X chromosome folding conformations in these cells with megabase genomic resolution. X chromosomes in female naïve hPSCs exhibit folding conformations similar to the active X chromosome (Xa) and the inactive X chromosome (Xi) in somatic cells. However, naïve X chromosomes do not exhibit the chromatin compaction typically associated with these somatic X chromosome states. In H7 naïve human embryonic stem cells, XIST accumulation observed on damaged X chromosomes demonstrates the potential for naïve hPSCs to activate XCI-related mechanisms. Overall, our findings provide insight into the X chromosome status of naïve hPSCs with a single-chromosome resolution and are critical in understanding the unique epigenetic regulation in early embryonic cells.

Chromosome-scale genomes of commercial timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis)

Wood is the most important natural and endlessly renewable source of energy. Despite the ecological and economic importance of wood, many aspects of its formation have not yet been investigated. We performed chromosome-scale genome assemblies of three timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis) which exhibit different wood properties such as wood density, hardness, growth rate, and fiber cell wall thickness. The combination of 10X, stLFR, Hi-Fi sequencing and HiC data led us to assemble high-quality genomes evident by scaffold N50 length of 55.97 Mb (O. pyramidale), 22.37 Mb (M. ferrea) and 14.55 Mb (T. grandis) with >97% BUSCO completeness of the assemblies. A total of 35774, 24027, and 44813 protein-coding genes were identified in M. ferrea, T. grandis and O. pyramidale, respectively. The data generated in this study is anticipated to serve as a valuable genetic resource and will promote comparative genomic analyses, and it is of practical importance in gaining a further understanding of the wood properties in non-model woody species.

An upgraded method of high-throughput chromosome conformation capture (Hi-C 3.0) in cotton (Gossypium spp.)

High-throughput chromosome conformation capture (Hi-C) technology has been applied to explore the chromatin interactions and shed light on the biological functions of three-dimensional genomic features. However, it remains challenging to guarantee the high quality of Hi-C library in plants and hence the reliable capture of chromatin structures, especially loops, due to insufficient fragmentation and low efficiency of proximity ligations. To overcome these deficiencies, we optimized the parameters of the Hi-C protocol, principally the cross-linking agents and endonuclease fragmentation strategy. The double cross-linkers (FA+DSG) and double restriction enzymes (DpnII+DdeI) were utilized. Thus, a systematic in situ Hi-C protocol was designed using plant tissues embedded with comprehensive quality controls to monitor the library construction. This upgraded method, termed Hi-C 3.0, was applied to cotton leaves for trial. In comparison with the conventional Hi-C 2.0, Hi-C 3.0 can obtain more than 50% valid contacts at a given sequencing depth to improve the signal-to-noise ratio. Hi-C 3.0 can furthermore enhance the capturing of loops almost as twice as that of Hi-C 2.0. In addition, Hi-C 3.0 showed higher efficiency of compartment detection and identified compartmentalization more accurately. In general, Hi-C 3.0 contributes to the advancement of the Hi-C method in plants by promoting its capability on decoding the chromatin organization.

High-quality chromosome-level scaffolds of the plant bug Pachypeltis micranthus provide insights into the availability of Mikania micrantha control

BACKGROUND

The plant bug, Pachypeltis micranthus Mu et Liu (Hemiptera: Miridae), is an effective potential biological control agent for Mikania micrantha H.B.K. (Asteraceae; one of the most notorious invasive weeds worldwide). However, limited knowledge about this species hindered its practical application and research. Accordingly, sequencing the genome of this mirid bug holds great significance in controlling M. micrantha.

RESULTS

Here, 712.72 Mb high-quality chromosome-level scaffolds of P. micranthus were generated, of which 707.51 Mb (99.27%) of assembled sequences were anchored onto 15 chromosome-level scaffolds with contig N50 of 16.84 Mb. The P. micranthus genome had the highest GC content (42.43%) and the second highest proportion of repetitive sequences (375.82 Mb, 52.73%) than the three other mirid bugs (i.e., Apolygus lucorum, Cyrtorhinus lividipennis, and Nesidiocoris tenuis). Phylogenetic analysis showed that P. micranthus clustered with other mirid bugs and diverged from the common ancestor approximately 200 million years ago. Gene family expansion and/or contraction were analyzed, and significantly expanded gene families associated with P. micranthus feeding and adaptation to M. micrantha were manually identified. Compared with the whole body, transcriptome analysis of the salivary gland revealed that most of the upregulated genes were significantly associated with metabolism pathways and peptidase activity, particularly among cysteine peptidase, serine peptidase, and polygalacturonase; this could be one of the reasons for precisely and highly efficient feeding by the oligophagous bug P. micranthus on M. micrantha.

CONCLUSION

Collectively, this work provides a crucial chromosome-level scaffolds resource to study the evolutionary adaptation between mirid bug and their host. It is also helpful in searching for novel environment-friendly biological strategies to control M. micrantha.

Novel chromosomes and genomes provide new insights into evolution and adaptation of the whole genome duplicated yeast-like fungus TN3-1 isolated from natural honey

Aureobasidium melanogenum TN3-1 strain and A. melanogenum P16 strain were isolated from the natural honey and the mangrove ecosystem, respectively. The former can produce much higher pullulan from high concentration of glucose than the latter. In order to know what happened to their genomes, the PacBio sequencing and Hi-C technologies were used to create the first high-quality chromosome-level reference genome assembly of A. melanogenum TN3-1 (51.61 Mb) and A. melanogenum P16 (25.82 Mb) with the contig N50 of 2.19 Mb and 2.26 Mb, respectively. Based on the Hi-C results, a total of 93.33% contigs in the TN3-1 strain and 92.31% contigs in the P16 strain were anchored onto 24 and 12 haploid chromosomes, respectively. The genomes of the TN3-1 strain had two subgenomes A and B. Synteny analysis showed that the genomic contents of the two subgenomes were asymmetric with many structural variations. Intriguingly, the TN3-1 strain was revealed as a recent hybrid/fusion between the ancestor of A. melanogenum CBS105.22/CBS110374 and the ancestor of another unidentified strain of A. melanogenum similar to P16 strain. We estimated that the two ancient progenitors diverged around 18.38 Mya and merged around 10.66-9.98 Mya. It was found that in the TN3-1 strain, telomeres of each chromosome contained high level of long interspersed nuclear elements (LINEs), but had low level of the telomerase encoding gene. Meanwhile, there were high level of transposable elements (TEs) inserted in the chromosomes of the TN3-1 strain. In addition, the positively selected genes of the TN3-1 strain were mainly enriched in the metabolic processes related to harsh environmental adaptability. Most of the stress-related genes were found to be related to the adjacent LTRs, and the glucose derepression was caused by the mutation of the Glc7-2 in the Snf-Mig1 system. All of these could contribute to its genetic instability, genome evolution, high stress resistance, and high pullulan production from glucose.

Protein arginine methyltransferase 5 (Prmt5) localizes to chromatin loop anchors and modulates expression of genes at TAD boundaries during early adipogenesis

Protein arginine methyltransferase 5 (Prmt5) is an essential regulator of embryonic development and adult progenitor cell functions. Prmt5 expression is mis-regulated in many cancers, and the development of Prmt5 inhibitors as cancer therapeutics is an active area of research. Prmt5 functions via effects on gene expression, splicing, DNA repair, and other critical cellular processes. We examined whether Prmt5 functions broadly as a genome-wide regulator of gene transcription and higher-order chromatin interactions during the initial stages of adipogenesis using ChIP-Seq, RNA-seq, and Hi-C using 3T3-L1 cells, a frequently utilized model for adipogenesis. We observed robust genome-wide Prmt5 chromatin-binding at the onset of differentiation. Prmt5 localized to transcriptionally active genomic regions, acting as both a positive and a negative regulator. A subset of Prmt5 binding sites co-localized with mediators of chromatin organization at chromatin loop anchors. Prmt5 knockdown decreased insulation strength at the boundaries of topologically associating domains (TADs) adjacent to sites with Prmt5 and CTCF co-localization. Genes overlapping such weakened TAD boundaries showed transcriptional dysregulation. This study identifies Prmt5 as a broad regulator of gene expression, including regulation of early adipogenic factors, and reveals an unappreciated requirement for Prmt5 in maintaining strong insulation at TAD boundaries and overall chromatin organization.

Evaluating the mouse neural precursor line, SN4741, as a suitable proxy for midbrain dopaminergic neurons

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques - karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq - to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo, mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes.

High-Quality Chromosome-Level Genome Assembly of the Corsac Fox (Vulpes corsac) Reveals Adaptation to Semiarid and Harsh Environments

The Corsac fox (Vulpes corsac) is a species of fox distributed in the arid prairie regions of Central and Northern Asia, with distinct adaptations to dry environments. Here, we applied Oxford-Nanopore sequencing and a chromosome structure capture technique to assemble the first Corsac fox genome, which was then assembled into chromosome fragments. The genome assembly has a total length of 2.2 Gb with a contig N50 of 41.62 Mb and a scaffold N50 of 132.2 Mb over 18 pseudo-chromosomal scaffolds. The genome contained approximately 32.67% of repeat sequences. A total of 20,511 protein-coding genes were predicted, of which 88.9% were functionally annotated. Phylogenetic analyses indicated a close relation to the Red fox (Vulpes vulpes) with an estimated divergence time of ~3.7 million years ago (MYA). We performed separate enrichment analyses of species-unique genes, the expanded and contracted gene families, and positively selected genes. The results suggest an enrichment of pathways related to protein synthesis and response and an evolutionary mechanism by which cells respond to protein denaturation in response to heat stress. The enrichment of pathways related to lipid and glucose metabolism, potentially preventing stress from dehydration, and positive selection of genes related to vision, as well as stress responses in harsh environments, may reveal adaptive evolutionary mechanisms in the Corsac fox under harsh drought conditions. Additional detection of positive selection for genes associated with gustatory receptors may reveal a unique desert diet strategy for the species. This high-quality genome provides a valuable resource for studying mammalian drought adaptation and evolution in the genus Vulpes.

High-resolution Hi-C maps highlight multiscale chromatin architecture reorganization during cold stress in Brachypodium distachyon

BACKGROUND

The adaptation of plants to cold stress involves changes in gene expression profiles that are associated with epigenetic regulation. Although the three-dimensional (3D) genome architecture is considered an important epigenetic regulator, the role of 3D genome organization in the cold stress response remains unclear.

RESULTS

In this study, we developed high-resolution 3D genomic maps using control and cold-treated leaf tissue of the model plant Brachypodium distachyon using Hi-C to determine how cold stress affects the 3D genome architecture. We generated ~ 1.5 kb resolution chromatin interaction maps and showed that cold stress disrupts different levels of chromosome organization, including A/B compartment transition, a reduction in chromatin compartmentalization and the size of topologically associating domains (TADs), and loss of long-range chromatin loops. Integrating RNA-seq information, we identified cold-response genes and revealed that transcription was largely unaffected by the A/B compartment transition. The cold-response genes were predominantly localized in compartment A. In contrast, transcriptional changes are required for TAD reorganization. We demonstrated that dynamic TAD events were associated with H3K27me3 and H3K27ac state alterations. Moreover, a loss of chromatin looping, rather than a gain of looping, coincides with alterations in gene expression, indicating that chromatin loop disruption may play a more important role than loop formation in the cold-stress response.

CONCLUSIONS

Our study highlights the multiscale 3D genome reprogramming that occurs during cold stress and expands our knowledge of the mechanisms underlying transcriptional regulation in response to cold stress in plants.

Mitotic chromosomes scale to nuclear-cytoplasmic ratio and cell size in Xenopus

During the rapid and reductive cleavage divisions of early embryogenesis, subcellular structures such as the nucleus and mitotic spindle scale to decreasing cell size. Mitotic chromosomes also decrease in size during development, presumably to scale coordinately with mitotic spindles, but the underlying mechanisms are unclear. Here we combine in vivo and in vitro approaches using eggs and embryos from the frog Xenopus laevis to show that mitotic chromosome scaling is mechanistically distinct from other forms of subcellular scaling. We found that mitotic chromosomes scale continuously with cell, spindle, and nuclear size in vivo. However, unlike for spindles and nuclei, mitotic chromosome size cannot be reset by cytoplasmic factors from earlier developmental stages. In vitro, increasing nuclear-cytoplasmic (N/C) ratio is sufficient to recapitulate mitotic chromosome scaling, but not nuclear or spindle scaling, through differential loading of maternal factors during interphase. An additional pathway involving importin α scales mitotic chromosomes to cell surface area/volume ratio (SA/V) during metaphase. Finally, single-chromosome immunofluorescence and Hi-C data suggest that mitotic chromosomes shrink during embryogenesis through decreased recruitment of condensin I, resulting in major rearrangements of DNA loop architecture to accommodate the same amount of DNA on a shorter chromosome axis. Together, our findings demonstrate how mitotic chromosome size is set by spatially and temporally distinct developmental cues in the early embryo.

SATB1 regulates 3D genome architecture in T cells by constraining chromatin interactions surrounding CTCF-binding sites

Special AT-rich sequence binding protein 1 (SATB1) has long been proposed to act as a global chromatin loop organizer in T cells. However, the exact functions of SATB1 in spatial genome organization remain elusive. Here we show that the depletion of SATB1 in human and murine T cells leads to transcriptional dysregulation for genes involved in T cell activation, as well as alterations of 3D genome architecture at multiple levels, including compartments, topologically associating domains, and loops. Importantly, SATB1 extensively colocalizes with CTCF throughout the genome. Depletion of SATB1 leads to increased chromatin contacts among and across the SATB1/CTCF co-occupied sites, thereby affecting the transcription of critical regulators of T cell activation. The loss of SATB1 does not affect CTCF occupancy but significantly reduces the retention of CTCF in the nuclear matrix. Collectively, our data show that SATB1 contributes to 3D genome organization by constraining chromatin topology surrounding CTCF-binding sites.

Chromosome-level genome assembly of Dastarcus helophoroides provides insights into CYP450 genes expression upon insecticide exposure

BACKGROUND

Dastarcus helophoroides is an important natural enemy of cerambycids, and is wildly used in biological control of pests. Nevertheless, the absence of complete genomic information limits the investigation of the underlying molecular mechanisms. Here, a chromosome-level of Dastarcus helophoroides genome is assembled using a combination strategy of Illumina, PacBio, 10x™ Genomics, and Hi-C.

RESULTS

The final assembly is 609.09 Mb with contig N50, scaffold N50 and GC content of 5.46 Mb, 42.56 Mb and 31.50%, respectively, and 95.25% of the contigs anchor into 13 chromosomes. In total 14 890 protein-coding genes and 65.37% repeat sequences are predicted in the assembly genome. The phylogenetic analysis of single-copy gene families shared among 20 insect species indicates that Dastarcus helophoroides is placed as the sister species to clade (Nitidulidae+Curculionoidea+Chrysomeloidea) + Tenebrionoidea, and diverges from the related species ~242.9 Mya. In total 36 expanded gene families are identified in Dastarcus helophoroides genome, and are functionally related to drug metabolism and metabolism of xenobiotics by cytochrome P450. Some members of CYP4 Clade and CYP6 Clade are up-regulated in Dastarcus helophoroides adults upon insecticide exposure, of which expressions of DhCYP4Q, DhCYP6A14X1 and DhCYP4C1 are significantly up-regulated. The silencing of the three genes leads to adults more sensitive to insecticide and increased knocked-down rate, which may indicate their critical roles in stress resistance and detoxication.

CONCLUSION

Our study systematically integrated the chromosome-level genome, transcriptome and gene expression of Dastarcus helophoroides, which will provide valuable resources for understanding mechanisms of pesticide metabolism, growth and development, and utilization of the natural enemy in integrated control. © 2022 Society of Chemical Industry.

The little skate genome and the evolutionary emergence of wing-like fins

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

3D genomics and its applications in precision medicine

Three-dimensional (3D) genomics is an emerging discipline that studies the three-dimensional structure of chromatin and the three-dimensional and functions of genomes. It mainly studies the three-dimensional conformation and functional regulation of intranuclear genomes, such as DNA replication, DNA recombination, genome folding, gene expression regulation, transcription factor regulation mechanism, and the maintenance of three-dimensional conformation of genomes. Self-chromosomal conformation capture (3C) technology has been developed, and 3D genomics and related fields have developed rapidly. In addition, chromatin interaction analysis techniques developed by 3C technologies, such as paired-end tag sequencing (ChIA-PET) and whole-genome chromosome conformation capture (Hi-C), enable scientists to further study the relationship between chromatin conformation and gene regulation in different species. Thus, the spatial conformation of plant, animal, and microbial genomes, transcriptional regulation mechanisms, interaction patterns of chromosomes, and the formation mechanism of spatiotemporal specificity of genomes are revealed. With the help of new experimental technologies, the identification of key genes and signal pathways related to life activities and diseases is sustaining the rapid development of life science, agriculture, and medicine. In this paper, the concept and development of 3D genomics and its application in agricultural science, life science, and medicine are introduced, which provides a theoretical basis for the study of biological life processes.

Using QTL to Identify Genes and Pathways Underlying the Regulation and Production of Milk Components in Cattle

Milk is a complex liquid, and the concentrations of many of its components are under genetic control. Many genes and pathways are known to regulate milk composition, and the purpose of this review is to highlight how the discoveries of quantitative trait loci (QTL) for milk phenotypes can elucidate these pathways. The main body of this review focuses primarily on QTL discovered in cattle (Bos taurus) as a model species for the biology of lactation, and there are occasional references to sheep genetics. The following section describes a range of techniques that can be used to help identify the causative genes underlying QTL when the underlying mechanism involves the regulation of gene expression. As genotype and phenotype databases continue to grow and diversify, new QTL will continue to be discovered, and although proving the causality of underlying genes and variants remains difficult, these new data sets will further enhance our understanding of the biology of lactation.

Influence of human peripheral blood samples preprocessing on the quality of Hi-C libraries

The genome-wide variant of the chromatin conformation capture technique (Hi-C) is a powerful tool for revealing patterns of genome spatial organization, as well as for understanding the effects of their disturbance on disease development. In addition, Hi-C can be used to detect chromosomal rearrangements, including balanced translocations and inversions. The use of the Hi-C method for the detection of chromosomal rearrangements is becoming more widespread. Modern high-throughput methods of genome analysis can effectively reveal point mutations and unbalanced chromosomal rearrangements. However, their sensitivity for determining translocations and inversions remains rather low. The storage of whole blood samples can affect the amount and integrity of genomic DNA, and it can distort the results of subsequent analyses if the storage was not under proper conditions. The Hi-C method is extremely demanding on the input material. The necessary condition for successfully applying Hi-C and obtaining high-quality data is the preservation of the spatial chromatin organization within the nucleus. The purpose of this study was to determine the optimal storage conditions of blood samples for subsequent Hi-C analysis. We selected 10 different conditions for blood storage and sample processing. For each condition, we prepared and sequenced Hi-C libraries. The quality of the obtained data was compared. As a result of the work, we formulated the requirements for the storage and processing of samples to obtain high-quality Hi-C data. We have established the minimum volume of blood sufficient for conducting Hi-C analysis. In addition, we have identified the most suitable methods for isolation of peripheral blood mononuclear cells and their long-term storage. The main requirement we have formulated is not to freeze whole blood.

A Cre-LoxP-based approach for combinatorial chromosome rearrangements in human HAP1 cells

Alterations of human karyotype caused by chromosomal rearrangements are often associated with considerable phenotypic effects. Studying molecular mechanisms underlying these effects requires an efficient and scalable experimental model. Here, we propose a Cre-LoxP-based approach for the generation of combinatorial diversity of chromosomal rearrangements. We demonstrate that using the developed system, both intra- and inter-chromosomal rearrangements can be induced in the human haploid HAP1 cells, although the latter is significantly less effective. The obtained genetically modified HAP1 cell line can be used to dissect genomic effects associated with intra-chromosomal structural variations.

Evaluating the mouse neural precursor line, SN4741, as a suitable proxy for midbrain dopaminergic neurons

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate of disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques - karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq - to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo , mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable a proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes.

Hi-C analysis of genomic contacts revealed karyotype abnormalities in chicken HD3 cell line

BACKGROUND

Karyotype abnormalities are frequent in immortalized continuous cell lines either transformed or derived from primary tumors. Chromosomal rearrangements can cause dramatic changes in gene expression and affect cellular phenotype and behavior during in vitro culture. Structural variations of chromosomes in many continuous mammalian cell lines are well documented, but chromosome aberrations in cell lines from other vertebrate models often remain understudied. The chicken LSCC-HD3 cell line (HD3), generated from erythroid precursors, was used as an avian model for erythroid differentiation and lineage-specific gene expression. However, karyotype abnormalities in the HD3 cell line were not assessed. In the present study, we applied high-throughput chromosome conformation capture to analyze 3D genome organization and to detect chromosome rearrangements in the HD3 cell line.

RESULTS

We obtained Hi-C maps of genomic interactions for the HD3 cell line and compared A/B compartments and topologically associating domains between HD3 and several other cell types. By analysis of contact patterns in the Hi-C maps of HD3 cells, we identified more than 25 interchromosomal translocations of regions ≥ 200 kb on both micro- and macrochromosomes. We classified most of the observed translocations as unbalanced, leading to the formation of heteromorphic chromosomes. In many cases of microchromosome rearrangements, an entire microchromosome together with other macro- and microchromosomes participated in the emergence of a derivative chromosome, resembling "chromosomal fusions'' between acrocentric microchromosomes. Intrachromosomal inversions, deletions and duplications were also detected in HD3 cells. Several of the identified simple and complex chromosomal rearrangements, such as between GGA2 and GGA1qter; GGA5, GGA4p and GGA7p; GGA4q, GGA6 and GGA19; and duplication of the sex chromosome GGAW, were confirmed by FISH.

CONCLUSIONS

In the erythroid progenitor HD3 cell line, in contrast to mature and immature erythrocytes, the genome is organized into distinct topologically associating domains. The HD3 cell line has a severely rearranged karyotype with most of the chromosomes engaged in translocations and can be used in studies of genome structure-function relationships. Hi-C proved to be a reliable tool for simultaneous assessment of the spatial genome organization and chromosomal aberrations in karyotypes of birds with a large number of microchromosomes.

Evaluating the mouse neural precursor line, SN4741, as a suitable proxy for midbrain dopaminergic neurons

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate of disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques - karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq - to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo , mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable a proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes.

Three-dimensional chromatin architecture datasets for aging and Alzheimer's disease

Recently, increasing studies are indicating a close association between dysregulated enhancers and neurodegenerative diseases, such as Alzheimer's disease (AD). However, their contributions were poorly defined for lacking direct links to disease genes. To bridge this gap, we presented the Hi-C datasets of 4 AD patients, 4 dementia-free aged and 3 young subjects, including 30 billion reads. As applications, we utilized them to link the AD risk SNPs and dysregulated epigenetic marks to the target genes. Combining with epigenetic data, we observed more detailed interactions among regulatory regions and found that many known AD risk genes were under long-distance promoter-enhancer interactions. For future AD and aging studies, our datasets provide a reference landscape to better interpret findings of association and epigenetic studies for AD and aging process.