Structural variation in the 3D genome

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 and 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, and TAD boundaries) and analyze their pros and cons. We also point out obstacles to the computational prediction of 3D interactions and suggest future research directions.

Inferring replication timing and proliferation dynamics from single-cell DNA sequencing data

Dysregulated DNA replication is a cause and a consequence of aneuploidy in cancer, yet the interplay between copy number alterations (CNAs), replication timing (RT) and cell cycle dynamics remain understudied in aneuploid tumors. We developed a probabilistic method, PERT, for simultaneous inference of cell-specific replication and copy number states from single-cell whole genome sequencing (scWGS) data. We used PERT to investigate clone-specific RT and proliferation dynamics in  >50,000 cells obtained from aneuploid and clonally heterogeneous cell lines, xenografts and primary cancers. We observed bidirectional relationships between RT and CNAs, with CNAs affecting X-inactivation producing the largest RT shifts. Additionally, we found that clone-specific S-phase enrichment positively correlated with ground-truth proliferation rates in genomically stable but not unstable cells. Together, these results demonstrate robust computational identification of S-phase cells from scWGS data, and highlight the importance of RT and cell cycle properties in studying the genomic evolution of aneuploid tumors.

Epigenetic Regulation is involved in Reversal of Obesity

Epigenetic processes play a crucial role in mediating the impact of environmental energetic challenges, from overconsumption to starvation. Over-nutrition of energy-dense foods and sedentary lifestyles contribute to the development of obesity, characterized by excessive fat storage and impaired metabolic signaling, stemming from disrupted brain signaling. Conversely, dieting and physical activity facilitate body weight rebalancing and trigger adaptive neural responses. These adaptations involve the upregulation of neurogenesis, synaptic plasticity and optimized brain function and energy homeostasis, balanced hormone signaling, normal metabolism, and reduced inflammation. The transition of the brain from a maladaptive to an adaptive state is partially guided by epigenetic mechanisms. While epigenetic mechanisms underlying obesity-related brain changes have been described, their role in mediating the reversal of maladaptation/obesity through lifestyle interventions remains less explored. This review focuses on elucidating epigenetic mechanisms involved in hypothalamic adaptations induced by lifestyle interventions. Given that lifestyle interventions are widely prescribed and accessible approaches for weight loss and maintenance, it is our challenge to uncover epigenetic mechanisms moderating these hypothalamic-functional beneficial changes.

Prediction of the 3D cancer genome from whole-genome sequencing using InfoHiC

The 3D genome prediction in cancer is crucial for uncovering the impact of structural variations (SVs) on tumorigenesis, especially when they are present in noncoding regions. We present InfoHiC, a systemic framework for predicting the 3D cancer genome directly from whole-genome sequencing (WGS). InfoHiC utilizes contig-specific copy number encoding on the SV contig assembly, and performs a contig-to-total Hi-C conversion for the cancer Hi-C prediction from multiple SV contigs. We showed that InfoHiC can predict 3D genome folding from all types of SVs using breast cancer cell line data. We applied it to WGS data of patients with breast cancer and pediatric patients with medulloblastoma, and identified neo topologically associating domains. For breast cancer, we discovered super-enhancer hijacking events associated with oncogenic overexpression and poor survival outcomes. For medulloblastoma, we found SVs in noncoding regions that caused super-enhancer hijacking events of medulloblastoma driver genes (GFI1, GFI1B, and PRDM6). In addition, we provide trained models for cancer Hi-C prediction from WGS at https://github.com/dmcb-gist/InfoHiC , uncovering the impacts of SVs in cancer patients and revealing novel therapeutic targets.

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 2 Mb 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.

Familial severe skeletal Class II malocclusion with gingival hyperplasia caused by a complex structural rearrangement at the KCNJ2-KCNJ16 locus

The aim of this work was to identify the underlying genetic cause in a four-generation family segregating an unusual phenotype comprising a severe form of skeletal Class II malocclusion with gingival hyperplasia. SNP array identified a copy number gain on chromosome 1 (chr1); however, this chromosomal region did not segregate correctly in the extended family. Exome sequencing also failed to identify a candidate causative variant but highlighted co-segregating genetic markers on chr17 and chr19. Short- and long-read genome sequencing allowed us to pinpoint and characterize at nucleotide-level resolution a chromothripsis-like complex rearrangement (CR) inserted into the chr17 co-segregating region at the KCNJ2-SOX9 locus. The CR involved the gain of five different regions from chr1 that are shuffled, chained, and inserted as a single block (∼828 kb) at chr17q24.3. The inserted sequences contain craniofacial enhancers that are predicted to interact with KCNJ2/KCNJ16 through neo-topologically associating domain (TAD) formation to induce ectopic activation. Our findings suggest that the CR inserted at chr17q24.3 is the cause of the severe skeletal Class II malocclusion with gingival hyperplasia in this family and expands the panoply of phenotypes linked to variation at the KCNJ2-SOX9 locus. In addition, we highlight a previously overlooked potential role for misregulation of the KCNJ2/KCNJ16 genes in the pathomechanism of gingival hyperplasia associated with deletions and other rearrangements of the 17q24.2-q24.3 region (MIM 135400).

VUS next in rare diseases? Deciphering genetic determinants of biomolecular condensation

The diagnostic odysseys for rare disease patients are getting shorter as next-generation sequencing becomes more widespread. However, the complex genetic diversity and factors influencing expressivity continue to challenge accurate diagnosis, leaving more than 50% of genetic variants categorized as variants of uncertain significance.Genomic expression intricately hinges on localized interactions among its products. Conventional variant prioritization, biased towards known disease genes and the structure-function paradigm, overlooks the potential impact of variants shaping the composition, location, size, and properties of biomolecular condensates, genuine membraneless organelles swiftly sensing and responding to environmental changes, and modulating expressivity.To address this complexity, we propose to focus on the nexus of genetic variants within biomolecular condensates determinants. Scrutinizing variant effects in these membraneless organelles could refine prioritization, enhance diagnostics, and unveil the molecular underpinnings of rare diseases. Integrating comprehensive genome sequencing, transcriptomics, and computational models can unravel variant pathogenicity and disease mechanisms, enabling precision medicine. This paper presents the rationale driving our proposal and describes a protocol to implement this approach. By fusing state-of-the-art knowledge and methodologies into the clinical practice, we aim to redefine rare diseases diagnosis, leveraging the power of scientific advancement for more informed medical decisions.

Annelid Comparative Genomics and the Evolution of Massive Lineage-Specific Genome Rearrangement in Bilaterians

The organization of genomes into chromosomes is critical for processes such as genetic recombination, environmental adaptation, and speciation. All animals with bilateral symmetry inherited a genome structure from their last common ancestor that has been highly conserved in some taxa but seemingly unconstrained in others. However, the evolutionary forces driving these differences and the processes by which they emerge have remained largely uncharacterized. Here, we analyze genome organization across the phylum Annelida using 23 chromosome-level annelid genomes. We find that while many annelid lineages have maintained the conserved bilaterian genome structure, the Clitellata, a group containing leeches and earthworms, possesses completely scrambled genomes. We develop a rearrangement index to quantify the extent of genome structure evolution and show that, compared to the last common ancestor of bilaterians, leeches and earthworms have among the most highly rearranged genomes of any currently sampled species. We further show that bilaterian genomes can be classified into two distinct categories-high and low rearrangement-largely influenced by the presence or absence, respectively, of chromosome fission events. Our findings demonstrate that animal genome structure can be highly variable within a phylum and reveal that genome rearrangement can occur both in a gradual, stepwise fashion, or rapid, all-encompassing changes over short evolutionary timescales.

Determination of the chromosomal position effects for plug-and-play application in the Myxococcus xanthus chassis cells

The chromosomal position effect can significantly affect the transgene expression, which may provide an efficient strategy for the inauguration of alien genes in new hosts, but has been less explored rationally. The bacterium Myxococcus xanthus harbors a large circular high-GC genome, and the position effect in this chassis may result in a thousand-fold expression variation of alien natural products. In this study, we conducted transposon insertion at TA sites on the M. xanthus genome, and used enrichment and dilution indexes to respectively appraise high and low expression potentials of alien genes at insertion sites. The enrichment sites are characteristically distributed along the genome, and the dilution sites are overlapped well with the horizontal transfer genes. We experimentally demonstrated the enrichment sites as high expression integration sites (HEISs), and the dilution sites unsuitable for gene integration expression. This work highlights that HEISs are the plug-and-play sites for efficient expression of integrated genes.

Rice Promoter Editing: An Efficient Genetic Improvement Strategy

Gene expression levels in rice (Oryza sativa L.) and other plant species are determined by the promoters, which directly control phenotypic characteristics. As essential components of genes, promoters regulate the intensity, location, and timing of gene expression. They contain numerous regulatory elements and serve as binding sites for proteins that modulate transcription, including transcription factors and RNA polymerases. Genome editing can alter promoter sequences, thereby precisely modifying the expression patterns of specific genes, and ultimately affecting the morphology, quality, and resistance of rice. This paper summarizes research on rice promoter editing conducted in recent years, focusing on improvements in yield, heading date, quality, and disease resistance. It is expected to inform the application of promoter editing and encourage further research and development in crop genetic improvement with promote.

Using Transcriptomics to Determine the Mechanism for the Resistance to Fusarium Head Blight of a Wheat-Th. elongatum Translocation Line

Fusarium head blight (FHB), caused by the Fusarium graminearum species complex, is a destructive disease in wheat worldwide. The lack of FHB-resistant germplasm is a barrier in wheat breeding for resistance to FHB. Thinopyrum elongatum is an important relative that has been successfully used for the genetic improvement of wheat. In this study, a translocation line, YNM158, with the YM158 genetic background carrying a fragment of diploid Th. elongatum 7EL chromosome created using 60Co-γ radiation, showed high resistance to FHB under both field and greenhouse conditions. Transcriptome analysis confirmed that the horizontal transfer gene, encoding glutathione S-transferase (GST), is an important contributor to FHB resistance in the pathogen infection stage, whereas the 7EL chromosome fragment carries other genes regulated by F. graminearum during the colonization stage. Introgression of the 7EL fragment affected the expression of wheat genes that were enriched in resistance pathways, including the phosphatidylinositol signaling system, protein processing in the endoplasmic reticulum, plant-pathogen interaction, and the mitogen-activated protein kinase (MAPK) signaling pathway at different stages after F. graminearium infection. This study provides a novel germplasm for wheat resistance to FHB and new insights into the molecular mechanisms of wheat resistance to FHB.

Clonal evolution of the 3D chromatin landscape in patients with relapsed pediatric B-cell acute lymphoblastic leukemia

Relapsed pediatric B-cell acute lymphoblastic leukemia (B-ALL) remains one of the leading causes of cancer mortality in children. We performed Hi-C, ATAC-seq, and RNA-seq on 12 matched diagnosis/relapse pediatric leukemia specimens to uncover dynamic structural variants (SVs) and 3D chromatin rewiring that may contribute to relapse. While translocations are assumed to occur early in leukemogenesis and be maintained throughout progression, we discovered novel, dynamic translocations and confirmed several fusion transcripts, suggesting functional and therapeutic relevance. Genome-wide chromatin remodeling was observed at all organizational levels: A/B compartments, TAD interactivity, and chromatin loops, including some loci shared by 25% of patients. Shared changes were found to drive the expression of genes/pathways previously implicated in resistance as well as novel therapeutic candidates, two of which (ATXN1 and MN1) we functionally validated. Overall, these results demonstrate chromatin reorganization under the selective pressure of therapy and offer the potential for discovery of novel therapeutic interventions.

Unraveling the three-dimensional (3D) genome architecture in Neurodevelopmental Disorders (NDDs)

The human genome, comprising millions of pairs of bases, serves as the blueprint of life, encoding instructions for cellular processes. However, genomes are not merely linear sequences; rather, the complex of DNA and histones, known as chromatin, exhibits complex organization across various levels, which profoundly influence gene expression and cellular function. Central to understanding genome organization is the emerging field of three-dimensional (3D) genome studies. Utilizing advanced techniques such as Hi-C, researchers have unveiled non-random dispositions of genomic elements, highlighting their importance in transcriptional regulation and disease mechanisms. Topologically Associating Domains (TADs), that demarcate regions of chromatin with preferential internal interactions, play crucial roles in gene regulation and are increasingly implicated in various diseases such as cancer and schizophrenia. However, their role in Neurodevelopmental Disorders (NDDs) remains poorly understood. Here, we focus on TADs and 3D conservation across the evolution and between cell types in NDDs. The investigation into genome organization and its impact on disease has led to significant breakthroughs in understanding NDDs etiology such ASD (Autism Spectrum Disorder). By elucidating the wide spectrum of ASD manifestations, researchers aim to uncover the underlying genetic and epigenetic factors contributing to its heterogeneity. Moreover, studies linking TAD disruption to NDDs underscore the importance of spatial genome organization in maintaining proper brain development and function. In summary, this review highlights the intricate interplay between genome organization, transcriptional control, and disease pathology, shedding light on fundamental biological processes and offering insights into the mechanisms underlying NDDs like ASD.

Cooperative insulation of regulatory domains by CTCF-dependent physical insulation and promoter competition

The specificity of gene expression during development requires the insulation of regulatory domains to avoid inappropriate enhancer-gene interactions. In vertebrates, this insulator function is mostly attributed to clusters of CTCF sites located at topologically associating domain (TAD) boundaries. However, TAD boundaries allow some physical crosstalk across regulatory domains, which is at odds with the specific and precise expression of developmental genes. Here we show that developmental genes and nearby clusters of CTCF sites cooperatively foster the robust insulation of regulatory domains. By genetically dissecting a couple of representative loci in mouse embryonic stem cells, we show that CTCF sites prevent undesirable enhancer-gene contacts (i.e. physical insulation), while developmental genes preferentially contribute to regulatory insulation through non-structural mechanisms involving promoter competition rather than enhancer blocking. Overall, our work provides important insights into the insulation of regulatory domains, which in turn might help interpreting the pathological consequences of certain structural variants.

VolcanoSV enables accurate and robust structural variant calling in diploid genomes from single-molecule long read sequencing

Structural variants (SVs) significantly contribute to human genome diversity and play a crucial role in precision medicine. Although advancements in single-molecule long-read sequencing offer a groundbreaking resource for SV detection, identifying SV breakpoints and sequences accurately and robustly remains challenging. We introduce VolcanoSV, an innovative hybrid SV detection pipeline that utilizes both a reference genome and local de novo assembly to generate a phased diploid assembly. VolcanoSV uses phased SNPs and unique k-mer similarity analysis, enabling precise haplotype-resolved SV discovery. VolcanoSV is adept at constructing comprehensive genetic maps encompassing SNPs, small indels, and all types of SVs, making it well-suited for human genomics studies. Our extensive experiments demonstrate that VolcanoSV surpasses state-of-the-art assembly-based tools in the detection of insertion and deletion SVs, exhibiting superior recall, precision, F1 scores, and genotype accuracy across a diverse range of datasets, including low-coverage (10x) datasets. VolcanoSV outperforms assembly-based tools in the identification of complex SVs, including translocations, duplications, and inversions, in both simulated and real cancer data. Moreover, VolcanoSV is robust to various evaluation parameters and accurately identifies breakpoints and SV sequences.

The Ruminant Telomere-to-Telomere (RT2T) Consortium

Telomere-to-telomere (T2T) assemblies reveal new insights into the structure and function of the previously 'invisible' parts of the genome and allow comparative analyses of complete genomes across entire clades. We present here an open collaborative effort, termed the 'Ruminant T2T Consortium' (RT2T), that aims to generate complete diploid assemblies for numerous species of the Artiodactyla suborder Ruminantia to examine chromosomal evolution in the context of natural selection and domestication of species used as livestock.

Genetic variation across and within individuals

Germline variation and somatic mutation are intricately connected and together shape human traits and disease risks. Germline variants are present from conception, but they vary between individuals and accumulate over generations. By contrast, somatic mutations accumulate throughout life in a mosaic manner within an individual due to intrinsic and extrinsic sources of mutations and selection pressures acting on cells. Recent advancements, such as improved detection methods and increased resources for association studies, have drastically expanded our ability to investigate germline and somatic genetic variation and compare underlying mutational processes. A better understanding of the similarities and differences in the types, rates and patterns of germline and somatic variants, as well as their interplay, will help elucidate the mechanisms underlying their distinct yet interlinked roles in human health and biology.

GAGE-seq concurrently profiles multiscale 3D genome organization and gene expression in single cells

The organization of mammalian genomes features a complex, multiscale three-dimensional (3D) architecture, whose functional significance remains elusive because of limited single-cell technologies that can concurrently profile genome organization and transcriptional activities. Here, we introduce genome architecture and gene expression by sequencing (GAGE-seq), a scalable, robust single-cell co-assay measuring 3D genome structure and transcriptome simultaneously within the same cell. Applied to mouse brain cortex and human bone marrow CD34+ cells, GAGE-seq characterized the intricate relationships between 3D genome and gene expression, showing that multiscale 3D genome features inform cell-type-specific gene expression and link regulatory elements to target genes. Integration with spatial transcriptomic data revealed in situ 3D genome variations in mouse cortex. Observations in human hematopoiesis unveiled discordant changes between 3D genome organization and gene expression, underscoring a complex, temporal interplay at the single-cell level. GAGE-seq provides a powerful, cost-effective approach for exploring genome structure and gene expression relationships at the single-cell level across diverse biological contexts.

Structural Variants at the LMNB1 Locus: Deciphering Pathomechanisms in Autosomal Dominant Adult-Onset Demyelinating Leukodystrophy

OBJECTIVES

We aimed to elucidate the pathogenic mechanisms underlying autosomal dominant adult-onset demyelinating leukodystrophy (ADLD), and to understand the genotype/phenotype correlation of structural variants (SVs) in the LMNB1 locus.

BACKGROUND

Since the discovery of 3D genome architectures and topologically associating domains (TADs), new pathomechanisms have been postulated for SVs, regardless of gene dosage changes. ADLD is a rare genetic disease associated with duplications (classical ADLD) or noncoding deletions (atypical ADLD) in the LMNB1 locus.

METHODS

High-throughput chromosome conformation capture, RNA sequencing, histopathological analyses of postmortem brain tissues, and clinical and neuroradiological investigations were performed.

RESULTS

We collected data from >20 families worldwide carrying SVs in the LMNB1 locus and reported strong clinical variability, even among patients carrying duplications of the entire LMNB1 gene, ranging from classical and atypical ADLD to asymptomatic carriers. We showed that patients with classic ADLD always carried intra-TAD duplications, resulting in a simple gene dose gain. Atypical ADLD was caused by LMNB1 forebrain-specific misexpression due to inter-TAD deletions or duplications. The inter-TAD duplication, which extends centromerically and crosses the 2 TAD boundaries, did not cause ADLD. Our results provide evidence that astrocytes are key players in ADLD pathology.

INTERPRETATION

Our study sheds light on the 3D genome and TAD structural changes associated with SVs in the LMNB1 locus, and shows that a duplication encompassing LMNB1 is not sufficient per se to diagnose ADLD, thereby strongly affecting genetic counseling. Our study supports breaking TADs as an emerging pathogenic mechanism that should be considered when studying brain diseases. ANN NEUROL 2024.

Landscape of genomic structural variations in Indian population-based cohorts: Deeper insights into their prevalence and clinical relevance

Structural variations (SV) are large (>50 base pairs) genomic rearrangements comprising deletions, duplications, insertions, inversions, and translocations. Studying SVs is important because they play active and critical roles in regulating gene expression, determining disease predispositions, and identifying population-specific differences among individuals of diverse ancestries. However, SV discoveries in the Indian population using whole-genome sequencing (WGS) have been limited. In this study, using short-read WGS having an average 42X depth of coverage, we identify and characterize 36,210 SVs from 529 individuals enrolled in population-based cohorts in India. These SVs include 24,574 deletions, 2,913 duplications, 8,710 insertions, and 13 inversions; 1.26% (456 out of 36,210) of the identified SVs can potentially impact the coding regions of genes. Furthermore, 56 of these SVs are highly intolerant to loss-of-function changes to the mapped genes, and five SVs impacting ADAMTS17, CCDC40, and RHCE are common in our study individuals. Seven rare SVs significantly impact dosage sensitivity of genes known to be associated with various clinical phenotypes. Most of the SVs in our study are rare and heterozygous. This fine-scale SV discovery in the underrepresented Indian population provides valuable insights that extend beyond Eurocentric human genetic studies.

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.

A practical introduction to holo-omics

Holo-omics refers to the joint study of non-targeted molecular data layers from host-microbiota systems or holobionts, which is increasingly employed to disentangle the complex interactions between the elements that compose them. We navigate through the generation, analysis, and integration of omics data, focusing on the commonalities and main differences to generate and analyze the various types of omics, with a special focus on optimizing data generation and integration. We advocate for careful generation and distillation of data, followed by independent exploration and analyses of the single omic layers to obtain a better understanding of the study system, before the integration of multiple omic layers in a final model is attempted. We highlight critical decision points to achieve this aim and flag the main challenges to address complex biological questions regarding the integrative study of host-microbiota relationships.

Allele-level visualization of transcription and chromatin by high-throughput imaging

The spatial arrangement of the genome within the nucleus is a pivotal aspect of cellular organization and function with implications for gene expression and regulation. While all genome organization features, such as loops, domains, and radial positioning, are nonrandom, they are characterized by a high degree of single-cell variability. Imaging approaches are ideally suited to visualize, measure, and study single-cell heterogeneity in genome organization. Here, we describe two methods for the detection of DNA and RNA of individual gene alleles by fluorescence in situ hybridization (FISH) in a high-throughput format. We have optimized combined DNA/RNA FISH approaches either using simultaneous or sequential detection of DNA and nascent RNA. These optimized DNA and RNA FISH protocols were implemented in a 384-well plate format alongside automated image and data analysis and enable accurate detection of individual gene alleles and their gene expression status across a large cell population. We successfully visualized MYC and EGFR DNA and nascent RNA with allele-level resolution in multiple cell types, and we determined the radial position of active and inactive MYC and EGFR alleles. These optimized DNA/RNA detection approaches are versatile and sensitive tools for mapping of chromatin features and gene activity at the single-allele level and at high throughput.

Facing the challenges to shorten the diagnostic odyssey: first Whole Genome Sequencing experience of a Colombian cohort with suspected rare diseases

Exome and genome sequencing (ES/GS) are routinely used for the diagnosis of genetic diseases in developed countries. However, their implementation is limited in countries from Latin America. We aimed to describe the results of GS in patients with suspected rare genetic diseases in Colombia. We studied 501 patients from 22 healthcare sites from January to December 2022. GS was performed in the index cases using dried blood spots on filtercards. Ancestry analysis was performed under iAdmix. Multiomic testing was performed when needed (biomarker, enzymatic activity, RNA-seq). All tests were performed at an accredited genetic laboratory. Ethnicity prediction data confirmed that 401 patients (80%) were mainly of Amerindian origin. A genetic diagnosis was established for 142 patients with a 28.3% diagnostic yield. The highest diagnostic yield was achieved for pathologies with a metabolic component and syndromic disorders (p < 0.001). Young children had a median of 1 year of diagnostic odyssey, while the median time for adults was significantly longer (15 years). Patients with genetic syndromes have spent more than 75% of their life without a diagnosis, while for patients with neurologic and neuromuscular diseases, the time of the diagnostic odyssey tended to decrease with age. Previous testing, specifically karyotyping or chromosomal microarray were significantly associated with a longer time to reach a definitive diagnosis (p < 0.01). Furthermore, one out of five patients that had an ES before could be diagnosed by GS. The Colombian genome project is the first Latin American study reporting the experience of systematic use of diagnostic GS in rare diseases.

Genetic/epigenetic effects in NF1 microdeletion syndrome: beyond the haploinsufficiency, looking at the contribution of not deleted genes

NF1 microdeletion syndrome, accounting for 5-11% of NF1 patients, is caused by a deletion in the NF1 region and it is generally characterized by a severe phenotype. Although 70% of NF1 microdeletion patients presents the same 1.4 Mb type-I deletion, some patients may show additional clinical features. Therefore, the contribution of several pathogenic mechanisms, besides haploinsufficiency of some genes within the deletion interval, is expected and needs to be defined. We investigated an altered expression of deletion flanking genes by qPCR in patients with type-1 NF1 deletion, compared to healthy donors, possibly contributing to the clinical traits of NF1 microdeletion syndrome. In addition, the 1.4-Mb deletion leads to changes in the 3D chromatin structure in the 17q11.2 region. Specifically, this deletion alters DNA-DNA interactions in the regions flanking the breakpoints, as demonstrated by our 4C-seq analysis. This alteration likely causes position effect on the expression of deletion flanking genes.Interestingly, 4C-seq analysis revealed that in microdeletion patients, an interaction was established between the RHOT1 promoter and the SLC6A4 gene, which showed increased expression. We performed NGS on putative modifier genes, and identified two "likely pathogenic" rare variants in RAS pathway, possibly contributing to incidental phenotypic features.This study provides new insights into understanding the pathogenesis of NF1 microdeletion syndrome and suggests a novel pathomechanism that contributes to the expression phenotype in addition to haploinsufficiency of genes located within the deletion.This is a pivotal approach that can be applied to unravel microdeletion syndromes, improving precision medicine, prognosis and patients' follow-up.

Building Haplotype-Resolved 3D Genome Maps of Chicken Skeletal Muscle

Haplotype-resolved 3D chromatin architecture related to allelic differences in avian skeletal muscle development has not been addressed so far, although chicken husbandry for meat consumption has been prevalent feature of cultures on every continent for more than thousands of years. Here, high-resolution Hi-C diploid maps (1.2-kb maximum resolution) are generated for skeletal muscle tissues in chicken across three developmental stages (embryonic day 15 to day 30 post-hatching). The sequence features governing spatial arrangement of chromosomes and characterize homolog pairing in the nucleus, are identified. Multi-scale characterization of chromatin reorganization between stages from myogenesis in the fetus to myofiber hypertrophy after hatching show concordant changes in transcriptional regulation by relevant signaling pathways. Further interrogation of parent-of-origin-specific chromatin conformation supported that genomic imprinting is absent in birds. This study also reveals promoter-enhancer interaction (PEI) differences between broiler and layer haplotypes in skeletal muscle development-related genes are related to genetic variation between breeds, however, only a minority of breed-specific variations likely contribute to phenotypic divergence in skeletal muscle potentially via allelic PEI rewiring. Beyond defining the haplotype-specific 3D chromatin architecture in chicken, this study provides a rich resource for investigating allelic regulatory divergence among chicken breeds.

How structural variants shape avian phenotypes: Lessons from model systems

Despite receiving significant recent attention, the relevance of structural variation (SV) in driving phenotypic diversity remains understudied, although recent advances in long-read sequencing, bioinformatics and pangenomic approaches have enhanced SV detection. We review the role of SVs in shaping phenotypes in avian model systems, and identify some general patterns in SV type, length and their associated traits. We found that most of the avian SVs so far identified are short indels in chickens, which are frequently associated with changes in body weight and plumage colouration. Overall, we found that relatively short SVs are more frequently detected, likely due to a combination of their prevalence compared to large SVs, and a detection bias, stemming primarily from the widespread use of short-read sequencing and associated analytical methods. SVs most commonly involve non-coding regions, especially introns, and when patterns of inheritance were reported, SVs associated primarily with dominant discrete traits. We summarise several examples of phenotypic convergence across different species, mediated by different SVs in the same or different genes and different types of changes in the same gene that can lead to various phenotypes. Complex rearrangements and supergenes, which can simultaneously affect and link several genes, tend to have pleiotropic phenotypic effects. Additionally, SVs commonly co-occur with single-nucleotide polymorphisms, highlighting the need to consider all types of genetic changes to understand the basis of phenotypic traits. We end by summarising expectations for when long-read technologies become commonly implemented in non-model birds, likely leading to an increase in SV discovery and characterisation. The growing interest in this subject suggests an increase in our understanding of the phenotypic effects of SVs in upcoming years.

Accelerating 3D genomics data analysis with Microcket

The three-dimensional (3D) organization of genome is fundamental to cell biology. To explore 3D genome, emerging high-throughput approaches have produced billions of sequencing reads, which is challenging and time-consuming to analyze. Here we present Microcket, a package for mapping and extracting interacting pairs from 3D genomics data, including Hi-C, Micro-C, and derivant protocols. Microcket utilizes a unique read-stitch strategy that takes advantage of the long read cycles in modern DNA sequencers; benchmark evaluations reveal that Microcket runs much faster than the current tools along with improved mapping efficiency, and thus shows high potential in accelerating and enhancing the biological investigations into 3D genome. Microcket is freely available at https://github.com/hellosunking/Microcket .

msBERT-Promoter: a multi-scale ensemble predictor based on BERT pre-trained model for the two-stage prediction of DNA promoters and their strengths

BACKGROUND

A promoter is a specific sequence in DNA that has transcriptional regulatory functions, playing a role in initiating gene expression. Identifying promoters and their strengths can provide valuable information related to human diseases. In recent years, computational methods have gained prominence as an effective means for identifying promoter, offering a more efficient alternative to labor-intensive biological approaches.

RESULTS

In this study, a two-stage integrated predictor called "msBERT-Promoter" is proposed for identifying promoters and predicting their strengths. The model incorporates multi-scale sequence information through a tokenization strategy and fine-tunes the DNABERT model. Soft voting is then used to fuse the multi-scale information, effectively addressing the issue of insufficient DNA sequence information extraction in traditional models. To the best of our knowledge, this is the first time an integrated approach has been used in the DNABERT model for promoter identification and strength prediction. Our model achieves accuracy rates of 96.2% for promoter identification and 79.8% for promoter strength prediction, significantly outperforming existing methods. Furthermore, through attention mechanism analysis, we demonstrate that our model can effectively combine local and global sequence information, enhancing its interpretability.

CONCLUSIONS

msBERT-Promoter provides an effective tool that successfully captures sequence-related attributes of DNA promoters and can accurately identify promoters and predict their strengths. This work paves a new path for the application of artificial intelligence in traditional biology.

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.

Drought response revealed by chromatin organization variation and transcriptional regulation in cotton

BACKGROUND

Cotton is a major world cash crop and an important source of natural fiber, oil, and protein. Drought stress is becoming a restrictive factor affecting cotton production. To facilitate the development of drought-tolerant cotton varieties, it is necessary to study the molecular mechanism of drought stress response by exploring key drought-resistant genes and related regulatory factors.

RESULTS

In this study, two cotton varieties, ZY007 (drought-sensitive) and ZY168 (drought-tolerant), showing obvious phenotypic differences under drought stress, were selected. A total of 25,898 drought-induced genes were identified, exhibiting significant enrichment in pathways related to plant stress responses. Under drought induction, At subgenome expression bias was observed at the whole-genome level, which may be due to stronger inhibition of Dt subgenome expression. A gene co-expression module that was significantly associated with drought resistance was identified. About 90% of topologically associating domain (TAD) boundaries were stable, and 6613 TAD variation events were identified between the two varieties under drought. We identified 92 genes in ZY007 and 98 in ZY168 related to chromatin 3D structural variation and induced by drought stress. These genes are closely linked to the cotton response to drought stress through canonical hormone-responsive pathways, modulation of kinase and phosphatase activities, facilitation of calcium ion transport, and other related molecular mechanisms.

CONCLUSIONS

These results lay a foundation for elucidating the molecular mechanism of the cotton drought response and provide important regulatory locus and gene resources for the future molecular breeding of drought-resistant cotton varieties.

Pervasive structural heterogeneity rewires glioblastoma chromosomes to sustain patient-specific transcriptional programs

Glioblastoma multiforme (GBM) encompasses brain malignancies marked by phenotypic and transcriptional heterogeneity thought to render these tumors aggressive, resistant to therapy, and inevitably recurrent. However, little is known about how the spatial organization of GBM genomes underlies this heterogeneity and its effects. Here, we compile a cohort of 28 patient-derived glioblastoma stem cell-like lines (GSCs) known to reflect the properties of their tumor-of-origin; six of these were primary-relapse tumor pairs from the same patient. We generate and analyze 5 kbp-resolution chromosome conformation capture (Hi-C) data from all GSCs to systematically map thousands of standalone and complex structural variants (SVs) and the multitude of neoloops arising as a result. By combining Hi-C, histone modification, and gene expression data with chromatin folding simulations, we explain how the pervasive, uneven, and idiosyncratic occurrence of neoloops sustains tumor-specific transcriptional programs via the formation of new enhancer-promoter contacts. We also show how even moderately recurrent neoloops can relate to patient-specific vulnerabilities. Together, our data provide a resource for dissecting GBM biology and heterogeneity, as well as for informing therapeutic approaches.

The probability of chromatin to be at the nuclear lamina has no systematic effect on its transcription level in fruit flies

BACKGROUND

Multiple studies have demonstrated a negative correlation between gene expression and positioning of genes at the nuclear envelope (NE) lined by nuclear lamina, but the exact relationship remains unclear, especially in light of the highly stochastic, transient nature of the gene association with the NE.

RESULTS

In this paper, we ask whether there is a causal, systematic, genome-wide relationship between the expression levels of the groups of genes in topologically associating domains (TADs) of Drosophila nuclei and the probabilities of TADs to be found at the NE. To investigate the nature of this possible relationship, we combine a coarse-grained dynamic model of the entire Drosophila nucleus with genome-wide gene expression data; we analyze the TAD averaged transcription levels of genes against the probabilities of individual TADs to be in contact with the NE in the control and lamins-depleted nuclei. Our findings demonstrate that, within the statistical error margin, the stochastic positioning of Drosophila melanogaster TADs at the NE does not, by itself, systematically affect the mean level of gene expression in these TADs, while the expected negative correlation is confirmed. The correlation is weak and disappears completely for TADs not containing lamina-associated domains (LADs) or TADs containing LADs, considered separately. Verifiable hypotheses regarding the underlying mechanism for the presence of the correlation without causality are discussed. These include the possibility that the epigenetic marks and affinity to the NE of a TAD are determined by various non-mutually exclusive mechanisms and remain relatively stable during interphase.

CONCLUSIONS

At the level of TADs, the probability of chromatin being in contact with the nuclear envelope has no systematic, causal effect on the transcription level in Drosophila. The conclusion is reached by combining model-derived time-evolution of TAD locations within the nucleus with their experimental gene expression levels.

Rare germline structural variants increase risk for pediatric solid tumors

Pediatric solid tumors are rare malignancies that represent a leading cause of death by disease among children in developed countries. The early age-of-onset of these tumors suggests that germline genetic factors are involved, yet conventional germline testing for short coding variants in established predisposition genes only identifies pathogenic events in 10-15% of patients. Here, we examined the role of germline structural variants (SVs)-an underexplored form of germline variation-in pediatric extracranial solid tumors using germline genome sequencing of 1,766 affected children, their 943 unaffected relatives, and 6,665 adult controls. We discovered a sex-biased association between very large (>1 megabase) germline chromosomal abnormalities and a four-fold increased risk of solid tumors in male children. The overall impact of germline SVs was greatest in neuroblastoma, where we revealed burdens of ultra-rare SVs that cause loss-of-function of highly expressed, mutationally intolerant, neurodevelopmental genes, as well as noncoding SVs predicted to disrupt three-dimensional chromatin domains in neural crest-derived tissues. Collectively, our results implicate rare germline SVs as a predisposing factor to pediatric solid tumors that may guide future studies and clinical practice.

Long-read sequencing and structural variant characterization in 1,019 samples from the 1000 Genomes Project

Structural variants (SVs) contribute significantly to human genetic diversity and disease 1-4 . Previously, SVs have remained incompletely resolved by population genomics, with short-read sequencing facing limitations in capturing the whole spectrum of SVs at nucleotide resolution 5-7 . Here we leveraged nanopore sequencing 8 to construct an intermediate coverage resource of 1,019 long-read genomes sampled within 26 human populations from the 1000 Genomes Project. By integrating linear and graph-based approaches for SV analysis via pangenome graph-augmentation, we uncover 167,291 sequence-resolved SVs in these samples, considerably advancing SV characterization compared to population-wide short-read sequencing studies 3,4 . Our analysis details diverse SV classes-deletions, duplications, insertions, and inversions-at population-scale. LINE-1 and SVA retrotransposition activities frequently mediate transductions 9,10 of unique sequences, with both mobile element classes transducing sequences at either the 3'- or 5'-end, depending on the source element locus. Furthermore, analyses of SV breakpoint junctions suggest a continuum of homology-mediated rearrangement processes are integral to SV formation, and highlight evidence for SV recurrence involving repeat sequences. Our open-access dataset underscores the transformative impact of long-read sequencing in advancing the characterisation of polymorphic genomic architectures, and provides a resource for guiding variant prioritisation in future long-read sequencing-based disease studies.

Predicting chromatin conformation contact maps

Over the past 15 years, a variety of next-generation sequencing assays have been developed for measuring the 3D conformation of DNA in the nucleus. Each of these assays gives, for a particular cell or tissue type, a distinct picture of 3D chromatin architecture. Accordingly, making sense of the relationship between genome structure and function requires teasing apart two closely related questions: how does chromatin 3D structure change from one cell type to the next, and how do different measurements of that structure differ from one another, even when the two assays are carried out in the same cell type? In this work, we assemble a collection of chromatin 3D datasets-each represented as a 2D contact map- spanning multiple assay types and cell types. We then build a machine learning model that predicts missing contact maps in this collection. We use the model to systematically explore how genome 3D architecture changes, at the level of compartments, domains, and loops, between cell type and between assay types.

Putative climate adaptation in American pikas (Ochotona princeps) is associated with copy number variation across environmental gradients

Improved understanding of the genetic basis of adaptation to climate change is necessary for maintaining global biodiversity moving forward. Studies to date have largely focused on sequence variation, yet there is growing evidence that suggests that changes in genome structure may be an even more significant source of adaptive potential. The American pika (Ochotona princeps) is an alpine specialist that shows some evidence of adaptation to climate along elevational gradients, but previous work has been limited to single nucleotide polymorphism based analyses within a fraction of the species range. Here, we investigated the role of copy number variation underlying patterns of local adaptation in the American pika using genome-wide data previously collected across the entire species range. We identified 37-193 putative copy number variants (CNVs) associated with environmental variation (temperature, precipitation, solar radiation) within each of the six major American pika lineages, with patterns of divergence largely following elevational and latitudinal gradients. Genes associated (n = 158) with independent annotations across lineages, variables, and/or CNVs had functions related to mitochondrial structure/function, immune response, hypoxia, olfaction, and DNA repair. Some of these genes have been previously linked to putative high elevation and/or climate adaptation in other species, suggesting they may serve as important targets in future studies.

Light regulation of rhodopsin distribution during outer segment renewal in murine rod photoreceptors

Vision under dim light relies on primary cilia elaborated by rod photoreceptors in the retina. This specialized sensory structure, called the rod outer segment (ROS), comprises hundreds of stacked, membranous discs containing the light-sensitive protein rhodopsin, and the incorporation of new discs into the ROS is essential for maintaining the rod's health and function. ROS renewal appears to be primarily regulated by extrinsic factors (light); however, results vary depending on different model organisms. We generated two independent transgenic mouse lines where rhodopsin's fate is tracked by a fluorescently labeled rhodopsin fusion protein (Rho-Timer) and show that rhodopsin incorporation into nascent ROS discs appears to be regulated by both external lighting cues and autonomous retinal clocks. Live-cell imaging of the ROS isolated from mice exposed to six unique lighting conditions demonstrates that ROS formation occurs in a periodic manner in cyclic light, constant darkness, and artificial light/dark cycles. This alternating bright/weak banding of Rho-Timer along the length of the ROS relates to inhomogeneities in rhodopsin density and potential points of structural weakness. In addition, we reveal that prolonged dim ambient light exposure impacts not only the rhodopsin content of new discs but also that of older discs, suggesting a dynamic interchange of material between new and old discs. Furthermore, we show that rhodopsin incorporation into the ROS is greatly altered in two autosomal recessive retinitis pigmentosa mouse models, potentially contributing to the pathogenesis. Our findings provide insights into how extrinsic (light) and intrinsic (retinal clocks and genetic mutation) factors dynamically regulate mammalian ROS renewal.

3C methods in cancer research: recent advances and future prospects

In recent years, Hi-C technology has revolutionized cancer research by elucidating the mystery of three-dimensional chromatin organization and its role in gene regulation. This paper explored the impact of Hi-C advancements on cancer research by delving into high-resolution techniques, such as chromatin loops, structural variants, haplotype phasing, and extrachromosomal DNA (ecDNA). Distant regulatory elements interact with their target genes through chromatin loops. Structural variants contribute to the development and progression of cancer. Haplotype phasing is crucial for understanding allele-specific genomic rearrangements and somatic clonal evolution in cancer. The role of ecDNA in driving oncogene amplification and drug resistance in cancer cells has also been revealed. These innovations offer a deeper understanding of cancer biology and the potential for personalized therapies. Despite these advancements, challenges, such as the accurate mapping of repetitive sequences and precise identification of structural variants, persist. Integrating Hi-C with multiomics data is key to overcoming these challenges and comprehensively understanding complex cancer genomes. Thus, Hi-C is a powerful tool for guiding precision medicine in cancer research and treatment.

Haplotype-resolved 3D chromatin architecture of the hybrid pig

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Functional associations of evolutionarily recent human genes exhibit sensitivity to the 3D genome landscape and disease

Genome organization is intricately tied to regulating genes and associated cell fate decisions. In this study, we examine the positioning and functional significance of human genes, grouped by their evolutionary age, within the 3D organization of the genome. We reveal that genes of different evolutionary origin have distinct positioning relationships with both domains and loop anchors, and remarkably consistent relationships with boundaries across cell types. While the functional associations of each group of genes are primarily cell type-specific, such associations of conserved genes maintain greater stability across 3D genomic features and disease than recently evolved genes. Furthermore, the expression of these genes across various tissues follows an evolutionary progression, such that RNA levels increase from young genes to ancient genes. Thus, the distinct relationships of gene evolutionary age, function, and positioning within 3D genomic features contribute to tissue-specific gene regulation in development and disease.

Unravelling inversions: Technological advances, challenges, and potential impact on crop breeding

Inversions, a type of chromosomal structural variation, significantly influence plant adaptation and gene functions by impacting gene expression and recombination rates. However, compared with other structural variations, their roles in functional biology and crop improvement remain largely unexplored. In this review, we highlight technological and methodological advancements that have allowed a comprehensive understanding of inversion variants through the pangenome framework and machine learning algorithms. Genome editing is an efficient method for inducing or reversing inversion mutations in plants, providing an effective mechanism to modify local recombination rates. Given the potential of inversions in crop breeding, we anticipate increasing attention on inversions from the scientific community in future research and breeding applications.

The impact of HTLV-1 expression on the 3D structure and expression of host chromatin

A typical HTLV-1-infected individual carries >104 different HTLV-1-infected T cell clones, each with a single-copy provirus integrated in a unique genomic site. We previously showed that the HTLV-1 provirus causes aberrant transcription in the flanking host genome and, by binding the chromatin architectural protein CTCF, forms abnormal chromatin loops with the host genome. However, it remained unknown whether these effects were exerted simply by the presence of the provirus or were induced by its transcription. To answer this question, we sorted HTLV-1-infected T-cell clones into cells positive or negative for proviral plus-strand expression, and then quantified host and provirus transcription using RNA-seq, and chromatin looping using quantitative chromosome conformation capture (q4C), in each cell population. We found that proviral plus-strand transcription induces aberrant transcription and splicing in the flanking genome but suppresses aberrant chromatin loop formation with the nearby host chromatin. Reducing provirus-induced host transcription with an inhibitor of transcriptional elongation allows recovery of chromatin loops in the plus-strand-expressing population. We conclude that aberrant host transcription induced by proviral expression causes temporary, reversible disruption of chromatin looping in the vicinity of the provirus.

SpecLoop predicts cell type-specific chromatin loop via transcription factor cooperation

Cell-type-Specific Chromatin Loops (CSCLs) are crucial for gene regulation and cell fate determination. However, the mechanisms governing their establishment remain elusive. Here, we present SpecLoop, a network regularization-based machine learning framework, to investigate the role of transcription factors (TFs) cooperation in CSCL formation. SpecLoop integrates multi-omics data, including gene expression, chromatin accessibility, sequence, protein-protein interaction, and TF binding motif data, to predict CSCLs and identify TF cooperations. Using high resolution Hi-C data as the gold standard, SpecLoop accurately predicts CSCL in GM12878, IMR90, HeLa-S3, K562, HUVEC, HMEC, and NHEK seven cell types, with the AUROC values ranging from 0.8645 to 0.9852 and AUPR values ranging from 0.8654 to 0.9734. Notably SpecLoop demonstrates improved accuracy in predicting long-distance CSCLs and identifies TF complexes with strong predictive ability. Our study systematically explores the TFs and TF pairs associated with CSCL through effective integration of diverse omics data. SpecLoop is freely available at https://github.com/AMSSwanglab/SpecLoop.

PIBF1 regulates multiple gene expression via impeding long-range chromatin interaction to drive the malignant transformation of HPV16 integration epithelial cells

INTRODUCTION

Human papillomavirus (HPV) integration can induce gene expression dysregulation by destroying higher-order chromatin structure in cervical cancer.

OBJECTIVES

We established a 13q22 site-specific HPV16 gene knock-in cell model to interrogate the changes in chromatin structure at the initial stages of host cell malignant transformation.

METHODS

We designed a CRISPR-Cas9 system with sgRNA targeting 13q22 site and constructed the HPV16 gene donor. Cells were cotransfected, screened, and fluorescence sorted. The whole genome sequencing (WGS) was used to confirm the precise HPV16 gene integration site. Western blot and qRT-PCR were used to measure gene expression. In vitro and in vivo analysis were performed to estimate the tumorigenic potential of the HPV16 knock-in cell model. Combined Hi-C, chromatin immunoprecipitation and RNA sequencing analyses revealed correlations between chromatin structure and gene expression. We performed a coimmunoprecipitation assay with anti-PIBF1 antibody to identify endogenous interacting proteins. In vivo analysis was used to determine the role of PIBF1 in the tumor growth of cervical cancer cells.

RESULTS

We successfully established a 13q22 site-specific HPV16 gene knock-in cell model. We found that HPV integration promoted cell proliferation, invasion and stratified growth in vitro, and monoclonal proliferation in vivo. HPV integration divided the affected topologically associated domain (TAD) into two smaller domains, and the progesterone-induced blocking factor 1 (PIBF1) gene near the integration site was upregulated, although PIBF1 was not enriched at the domain boundary by CUT-Tag signal analysis. Moreover, PIBF1 was found to interact with the cohesin complex off chromatin to reduce contact domain formation by disrupting the cohesin ring-shaped structure, causing dysregulation of tumorigenesis-related genes. Xenograft experiments determined the role of PIBF1 in the proliferation in cervical cancer cells.

CONCLUSION

We highlight that PIBF1, a potential chromatin structure regulatory protein, is activated by HPV integration, which provides new insights into HPV integration-driven cervical carcinogenesis.

Comparative study on genomic and epigenomic profiles of retinoblastoma or tuberous sclerosis complex via nanopore sequencing and a joint screening framework

Recurrence and extraocular metastasis in advanced intraocular retinoblastoma (RB) are still major obstacles for successful treatment of Chinese children. Tuberous sclerosis complex (TSC) is a very rare, multisystemic genetic disorder characterized by hamartomatous growth. In this study, we aimed to compare genomic and epigenomic profiles with human RB or TSC using recently developed nanopore sequencing, and to identify disease-associated variations or genes. Peripheral blood samples were collected from either RB or RB/TSC patients plus their normal siblings, followed by nanopore sequencing and identification of disease-specific structural variations (SVs) and differentially methylated regions (DMRs) by a systematic biology strategy named as multiomics-based joint screening framework. In total, 316 RB- and 1295 TSC-unique SVs were identified, as well as 1072 RB- and 1114 TSC-associated DMRs, respectively. We eventually identified 6 key genes for RB for further functional validation. Knockdown of CDK19 with specific siRNAs significantly inhibited Y79 cellular proliferation and increased sensitivity to carboplatin, whereas downregulation of AHNAK2 promoted the cell growth as well as drug resistance. Those two genes might serve as potential diagnostic markers or therapeutic targets of RB. The systematic biology strategy combined with functional validation might be an effective approach for rare pediatric malignances with limited samples and challenging collection process.

Subtelomeric plasticity contributes to gene family expansion in the human parasitic flatworm Schistosoma mansoni

BACKGROUND

The genomic region that lies between the telomere and chromosome body, termed the subtelomere, is heterochromatic, repeat-rich, and frequently undergoes rearrangement. Within this region, large-scale structural changes enable gene diversification, and, as such, large multicopy gene families are often found at the subtelomere. In some parasites, genes associated with proliferation, invasion, and survival are often found in these regions, where they benefit from the subtelomere's highly plastic, rapidly changing nature. The increasing availability of complete (or near complete) parasite genomes provides an opportunity to investigate these typically poorly defined and overlooked genomic regions and potentially reveal relevant gene families necessary for the parasite's lifestyle.

RESULTS

Using the latest chromosome-scale genome assembly and hallmark repeat richness observed at chromosome termini, we have identified and characterised the subtelomeres of Schistosoma mansoni, a metazoan parasitic flatworm that infects over 250 million people worldwide. Approximately 12% of the S. mansoni genome is classified as subtelomeric, and, in line with other organisms, we find these regions to be gene-poor but rich in transposable elements. We find that S. mansoni subtelomeres have undergone extensive interchromosomal recombination and that these sites disproportionately contribute to the 2.3% of the genome derived from segmental duplications. This recombination has led to the expansion of subtelomeric gene clusters containing 103 genes, including the immunomodulatory annexins and other gene families with unknown roles. The largest of these is a 49-copy plexin domain-containing protein cluster, exclusively expressed in the tegument-the tissue located at the host-parasite physical interface-of intramolluscan life stages.

CONCLUSIONS

We propose that subtelomeric regions act as a genomic playground for trial-and-error of gene duplication and subsequent divergence. Owing to the importance of subtelomeric genes in other parasites, gene families implicated in this subtelomeric expansion within S. mansoni warrant further characterisation for a potential role in parasitism.

Alteration of chromosome structure impacts gene expressions implicated in pancreatic ductal adenocarcinoma cells

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy with a five-year survival rate of approximately 10%. Genetic mutations are pivotal drivers in PDAC pathogenesis, but recent investigations also revealed the involvement of non-genetic alterations in the disease development. In this study, we undertook a multi-omics approach, encompassing ATAC-seq, RNA-seq, ChIP-seq, and Hi-C methodologies, to dissect gene expression alterations arising from changes in chromosome accessibility and chromatin three-dimensional interactions in PDAC.

RESULTS

Our findings indicate that chromosomal structural alterations can lead to abnormal expressions on key genes during PDAC development. Notably, overexpression of oncogenes FGFR2, FOXA2, CYP2R1, and CPOX can be attributed to the augmentation of promoter accessibility, coupled with long-range interactions with distal elements. Additionally, our findings indicate that chromosomal structural alterations caused by genomic instability can lead to abnormal expressions in PDACs. As an example, by analyzing chromosomal changes, we identified a putative oncogenic gene, LPAR1, which shows upregulated expression in both PDAC cell lines and clinical samples. The overexpression is correlated with alterations in LPAR1-associated 3D genome structure and chromatin state. We further demonstrated that high LPAR1 activity is required for enhanced PDAC cell migration in vitro.

CONCLUSIONS

Collectively, our findings reveal that the chromosomal conformational alterations, in addition to the well-known genetic mutations, are critical for PDAC tumorigenesis.

Allele-level visualization of transcription and chromatin by high-throughput imaging

The spatial arrangement of the genome within the nucleus is a pivotal aspect of cellular organization and function with implications for gene expression and regulation. While all genome organization features, such as loops, domains, and radial positioning, are non-random, they are characterized by a high degree of single-cell variability. Imaging approaches are ideally suited to visualize, measure, and study single-cell heterogeneity in genome organization. Here, we describe two methods for the detection of DNA and RNA of individual gene alleles by fluorescence in situ hybridization (FISH) in a high-throughput format. We have optimized combined DNA/RNA FISH approaches either using simultaneous or sequential detection. These optimized DNA and RNA FISH protocols, implemented in a 384-well plate format alongside automated image and data analysis, enable accurate detection of chromatin loci and their gene expression status across a large cell population with allele-level resolution. We successfully visualized MYC and EGFR DNA and RNA in multiple cell types, and we determined the radial position of active and inactive MYC and EGFR alleles. These optimized DNA/RNA detection approaches are versatile and sensitive tools for mapping of chromatin features and gene activity at the single-allele level and at high throughput.

Chromatin organization of muscle stem cell

The proper functioning of skeletal muscles is essential throughout life. A crucial crosstalk between the environment and several cellular mechanisms allows striated muscles to perform successfully. Notably, the skeletal muscle tissue reacts to an injury producing a completely functioning tissue. The muscle's robust regenerative capacity relies on the fine coordination between muscle stem cells (MuSCs or "satellite cells") and their specific microenvironment that dictates stem cells' activation, differentiation, and self-renewal. Critical for the muscle stem cell pool is a fine regulation of chromatin organization and gene expression. Acquiring a lineage-specific 3D genome architecture constitutes a crucial modulator of muscle stem cell function during development, in the adult stage, in physiological and pathological conditions. The context-dependent relationship between genome structure, such as accessibility and chromatin compartmentalization, and their functional effects will be analysed considering the improved 3D epigenome knowledge, underlining the intimate liaison between environmental encounters and epigenetics.