Insights about genome function from spatial organization of the genome

Artificially inserted strong promoter containing multiple G-quadruplexes induces long-range chromatin modification

Although the role of G-quadruplex (G4) DNA structures has been suggested in chromosomal looping this was not tested directly. Here, to test causal function, an array of G4s, or control sequence that does not form G4s, were inserted within chromatin in cells. In vivo G4 formation of the inserted G4 sequence array, and not the control sequence, was confirmed using G4-selective antibody. Compared to the control insert, we observed a remarkable increase in the number of 3D chromatin looping interactions from the inserted G4 array. This was evident within the immediate topologically associated domain (TAD) and throughout the genome. Locally, recruitment of enhancer histone marks and the transcriptional coactivator p300/Acetylated-p300 increased in the G4-array, but not in the control insertion. Resulting promoter-enhancer interactions and gene activation were clear up to 5 Mb away from the insertion site. Together, these show the causal role of G4s in enhancer function and long-range chromatin interactions. Mechanisms of 3D topology are primarily based on DNA-bound architectural proteins that induce/stabilize long-range interactions. Involvement of the underlying intrinsic DNA sequence/structure in 3D looping shown here therefore throws new light on how long-range chromosomal interactions might be induced or maintained.

A systematic analyses of different bioinformatics pipelines for genomic data and its impact on deep learning models for chromatin loop prediction

Genomic data analysis has witnessed a surge in complexity and volume, primarily driven by the advent of high-throughput technologies. In particular, studying chromatin loops and structures has become pivotal in understanding gene regulation and genome organization. This systematic investigation explores the realm of specialized bioinformatics pipelines designed specifically for the analysis of chromatin loops and structures. Our investigation incorporates two protein (CTCF and Cohesin) factor-specific loop interaction datasets from six distinct pipelines, amassing a comprehensive collection of 36 diverse datasets. Through a meticulous review of existing literature, we offer a holistic perspective on the methodologies, tools and algorithms underpinning the analysis of this multifaceted genomic feature. We illuminate the vast array of approaches deployed, encompassing pivotal aspects such as data preparation pipeline, preprocessing, statistical features and modelling techniques. Beyond this, we rigorously assess the strengths and limitations inherent in these bioinformatics pipelines, shedding light on the interplay between data quality and the performance of deep learning models, ultimately advancing our comprehension of genomic intricacies.

Spontaneous DNA Synapsis by Forming Noncanonical Intermolecular Structures

We report the spontaneous formation of DNA-DNA junctions in solution in the absence of proteins visualised using atomic force microscopy. The synapsis position fits with potential G-quadruplex (G4) sites. In contrast to the Holliday structure, these conjugates have an affinity for G4 antibodies. Molecular modelling was used to elucidate the possible G4/IM-synaptic complex structures. Our results indicate a new role of the intermolecular noncanonical structures in chromatin architecture and genomic rearrangement.

Contribution of 3D genome topological domains to genetic risk of cancers: a genome-wide computational study

Background

Genome-wide association studies have identified statistical associations between various diseases, including cancers, and a large number of single-nucleotide polymorphisms (SNPs). However, they provide no direct explanation of the mechanisms underlying the association. Based on the recent discovery that changes in three-dimensional genome organization may have functional consequences on gene regulation favoring diseases, we investigated systematically the genome-wide distribution of disease-associated SNPs with respect to a specific feature of 3D genome organization: topologically associating domains (TADs) and their borders.

Results

For each of 449 diseases, we tested whether the associated SNPs are present in TAD borders more often than observed by chance, where chance (i.e., the null model in statistical terms) corresponds to the same number of pointwise loci drawn at random either in the entire genome, or in the entire set of disease-associated SNPs listed in the GWAS catalog. Our analysis shows that a fraction of diseases displays such a preferential localization of their risk loci. Moreover, cancers are relatively more frequent among these diseases, and this predominance is generally enhanced when considering only intergenic SNPs. The structure of SNP-based diseasome networks confirms that localization of risk loci in TAD borders differs between cancers and non-cancer diseases. Furthermore, different TAD border enrichments are observed in embryonic stem cells and differentiated cells, consistent with changes in topological domains along embryogenesis and delineating their contribution to disease risk.

Conclusions

Our results suggest that, for certain diseases, part of the genetic risk lies in a local genetic variation affecting the genome partitioning in topologically insulated domains. Investigating this possible contribution to genetic risk is particularly relevant in cancers. This study thus opens a way of interpreting genome-wide association studies, by distinguishing two types of disease-associated SNPs: one with an effect on an individual gene, the other acting in interplay with 3D genome organization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-022-00375-2.

Python for gene expression

Genome biology shows substantial progress in its analytical and computational part in the last decades. Differential gene expression is one of many computationally intense areas; it is largely developed under R programming language. Here we explain possible reasons for such dominance of R in gene expression data. Next, we discuss the prospects for Python to become competitive in this area of research in coming years. We indicate that Python can be used already in a field of a single cell differential gene expression. We pinpoint still missing parts in Python and possibilities for improvement.

Python for gene expression

Genome biology shows substantial progress in its analytical and computational part in the last decades. Differential gene expression is one of many computationally intense areas; it is largely developed under R programming language. Here we explain possible reasons for such dominance of R in gene expression data. Next, we discuss the prospects for Python to become competitive in this area of research in coming years. We indicate that Python can be used already in a field of a single cell differential gene expression. We pinpoint still missing parts in Python and possibilities for improvement.

Neutron Scattering Techniques and Complementary Methods for Structural and Functional Studies of Biological Macromolecules and Large Macromolecular Complexes

Abstract

The review describes the application of small-angle scattering (SAS) of neutrons and complementary methods to study the structures of biomacromolecules. Here we cover SAS techniques, such as the contrast variation, the neutron spin-echo, and the solution of direct and inverse problems of three-dimensional reconstruction of the structures of macromolecules from SAS spectra by means of molecular modeling. A special section is devoted to specific objects of research, such as supramolecular complexes, influenza virus nucleoprotein, and chromatin.

Sept8/SEPTIN8 involvement in cellular structure and kidney damage is identified by genetic mapping and a novel human tubule hypoxic model

Chronic kidney disease (CKD), which can ultimately progress to kidney failure, is influenced by genetics and the environment. Genes identified in human genome wide association studies (GWAS) explain only a small proportion of the heritable variation and lack functional validation, indicating the need for additional model systems. Outbred heterogeneous stock (HS) rats have been used for genetic fine-mapping of complex traits, but have not previously been used for CKD traits. We performed GWAS for urinary protein excretion (UPE) and CKD related serum biochemistries in 245 male HS rats. Quantitative trait loci (QTL) were identified using a linear mixed effect model that tested for association with imputed genotypes. Candidate genes were identified using bioinformatics tools and targeted RNAseq followed by testing in a novel in vitro model of human tubule, hypoxia-induced damage. We identified two QTL for UPE and five for serum biochemistries. Protein modeling identified a missense variant within Septin 8 (Sept8) as a candidate for UPE. Sept8/SEPTIN8 expression increased in HS rats with elevated UPE and tubulointerstitial injury and in the in vitro hypoxia model. SEPTIN8 is detected within proximal tubule cells in human kidney samples and localizes with acetyl-alpha tubulin in the culture system. After hypoxia, SEPTIN8 staining becomes diffuse and appears to relocalize with actin. These data suggest a role of SEPTIN8 in cellular organization and structure in response to environmental stress. This study demonstrates that integration of a rat genetic model with an environmentally induced tubule damage system identifies Sept8/SEPTIN8 and informs novel aspects of the complex gene by environmental interactions contributing to CKD risk.

Different class IIa HDACs repressive complexes regulate specific epigenetic responses related to cell survival in leiomyosarcoma cells

Transcriptional networks supervising class IIa HDAC expression are poorly defined. Here we demonstrate that MEF2D is the key factor controlling HDAC9 transcription. This control, which is part of a negative feed-back loop during muscle differentiation, is hijacked in cancer. In leiomyosarcomas the MEF2D/HDAC9 vicious circuit sustains proliferation and cell survival, through the repression of the death receptor FAS. Comprehensive genome-wide studies demonstrate that HDAC4 and HDAC9 control different genetic programs and show both specific and common genomic binding sites. Although the number of MEF2-target genes commonly regulated is similar, only HDAC4 represses many additional genes that are not MEF2D targets. As expected, HDAC4-/- and HDAC9-/- cells increase H3K27ac levels around the TSS of the respective repressed genes. However, these genes rarely show binding of the HDACs at their promoters. Frequently HDAC4 and HDAC9 bind intergenic regions. We demonstrate that these regions, recognized by MEF2D/HDAC4/HDAC9 repressive complexes, show the features of active enhancers. In these regions HDAC4 and HDAC9 can differentially influence H3K27 acetylation. Our studies describe new layers of class IIa HDACs regulation, including a dominant positional effect, and can contribute to explain the pleiotropic actions of MEF2 TFs.

HOTAIRM1, an enhancer lncRNA, promotes glioma proliferation by regulating long-range chromatin interactions within HOXA cluster genes

The long noncoding RNA HOTAIRM1 reportedly plays important roles in acute myeloid leukemia, gastric cancer and colorectal cancer. Here, we analyzed potential function of HOTAIRM1 in glioma and asked whether it participates in long-range chromatin interactions. We monitored expression of HOTAIRM1 in glioma tissues and correlated levels with patient survival using the TCGA dataset. HOTAIRM1 was highly expressed in glioma tissue, with high levels associated with shortened patient survival time. We then suppressed HOTAIRM1 activity in the human glioblastoma U251 line using CRISPR-cas9 to knock in a truncating polyA fragment. Reporter analysis of these and control cells confirmed that the HOTAIRM1 locus serves as an active enhancer. We then performed Capture-C analysis to identify target genes of that locus and applied RNA antisense purification to assess chromatin interactions between the HOTAIRM1 locus and HOXA cluster genes. HOTAIRM1 knockdown in glioma cells decreased proliferation and reduced expression of HOXA cluster genes. HOTAIRM1 regulates long-range interactions between the HOTAIRM1 locus and HOXA genes. Our work suggests a new mechanism by which HOTAIRM1 regulates glioma progression by regulating high-order chromatin structure and could suggest novel therapeutic targets to treat an intractable cancer.

Tales from topographic oceans: topologically associated domains and cancer

The 3D organization of the genome within the cell nucleus has come into sharp focus over the last decade. This has largely arisen because of the application of genomic approaches that have revealed numerous levels of genomic and chromatin interactions, including topologically associated domains (TADs). The current review examines how these domains were identified, are organized, how their boundaries arise and are regulated, and how genes within TADs are coordinately regulated. There are many examples of the disruption to TAD structure in cancer and the altered regulation, structure and function of TADs are discussed in the context of hormone responsive cancers, including breast, prostate and ovarian cancer. Finally, some aspects of the statistical insight and computational skills required to interrogate TAD organization are considered and future directions discussed.

A tour of 3D genome with a focus on CTCF

The complex three-dimensional (3D) structure of the genome plays critical roles in the maintenance of genome stability, organization, and dynamics and in regulation of gene expression for understanding molecular mechanisms and diseases. Chromatin maintains biological functions and transcriptional activities through long distance interaction and interactions between loops and enhancers-promoters. We firstly overview the architecture and biology of chromatin and loops, topologically associated domains (TADs) and interactions, and compartments and functions. We specifically focus on CCCTC-binding factor (CTCF) in 3D genome organization and function to furthermore understand the significance of CTCF biology, transcriptional regulations, interactions with cohesin, roles in DNA binding, influences of CTCF degradation, and communication with wings-apart like (Wapl) protein. We also summarize the advanced single cell approaches to further monitor dynamics of CTCF functions and structures in the maintenance of 3D genome organization and function at single cell level.

The three-dimensional cancer nucleus

Research into the three-dimensional (3D) organization of the cancer cell genome started over 100 years ago. We follow an exciting avenue of research in this field, from Hansemann's early observations of aberrant mitoses and nuclei in cancer cells in the late 19th century to Boveri's theory of the cancer cell in the early 20th century, to current views of nuclear organization and its changes in cancer. Molecular and imaging methods go hand in hand with providing us with a better understanding of the spatial nature of the cancer cell genome. This has led to the concept that the structural order of the nucleus can be used as cancer cell biomarker.

The telomere bouquet is a hub where meiotic double-strand breaks, synapsis, and stable homolog juxtaposition are coordinated in the zebrafish, Danio rerio

Meiosis is a cellular program that generates haploid gametes for sexual reproduction. While chromosome events that contribute to reducing ploidy (homologous chromosome pairing, synapsis, and recombination) are well conserved, their execution varies across species and even between sexes of the same species. The telomere bouquet is a conserved feature of meiosis that was first described nearly a century ago, yet its role is still debated. Here we took advantage of the prominent telomere bouquet in zebrafish, Danio rerio, and super-resolution microscopy to show that axis morphogenesis, synapsis, and the formation of double-strand breaks (DSBs) all take place within the immediate vicinity of telomeres. We established a coherent timeline of events and tested the dependence of each event on the formation of Spo11-induced DSBs. First, we found that the axis protein Sycp3 loads adjacent to telomeres and extends inward, suggesting a specific feature common to all telomeres seeds the development of the axis. Second, we found that newly formed axes near telomeres engage in presynaptic co-alignment by a mechanism that depends on DSBs, even when stable juxtaposition of homologous chromosomes at interstitial regions is not yet evident. Third, we were surprised to discover that ~30% of telomeres in early prophase I engage in associations between two or more chromosome ends and these interactions decrease in later stages. Finally, while pairing and synapsis were disrupted in both spo11 males and females, their reproductive phenotypes were starkly different; spo11 mutant males failed to produce sperm while females produced offspring with severe developmental defects. Our results support zebrafish as an important vertebrate model for meiosis with implications for differences in fertility and genetically derived birth defects in males and females.

Inherent to reproduction is the transmission of genetic information from one generation to the next. In sexually reproducing organisms, each parent contributes an equal amount of genetic information, packaged in chromosomes, to the offspring. Diploid organisms, like humans, have two copies of every chromosome, while their haploid gametes (e.g. eggs and sperm) have only one. This reduction in ploidy depends on the segregation of chromosomes during meiosis, resulting in gametes with one copy of each chromosome. Missegregation of the chromosomes in the parents leads to abnormal chromosome numbers in the offspring, which is usually lethal or has detrimental developmental effects. While it has been known for over a century that homologous chromosomes pair and recombine to facilitate proper segregation, how homologs find their partners has remained elusive. A structure that has been central to the discussion of homolog pairing is the bouquet, or the dynamic clustering of telomeres during early stages of meiosis. Here we use zebrafish to show that the telomere bouquet is the site where key events leading to homologous chromosome pairing are coordinated. Furthermore, we show that deletion of spo11, a gene required for proper recombination in most studied organisms, resulted in very different effects in males and females where males were sterile while females produced deformed progeny.

Using Chromatin Architecture to Understand the Genetics and Transcriptomics of Juvenile Idiopathic Arthritis

The presence of abnormal gene expression signatures is a well-described feature of the oligoarticular and polyarticular forms of juvenile idiopathic arthritis. In this review, we discuss how new insights into genetic risk for JIA and the three dimensional architecture of the genome may be used to develop a better understanding of the mechanisms driving these gene expression patterns.

Chromosome positioning and male infertility: it comes with the territory

The production of functional spermatozoa through spermatogenesis requires a spatially and temporally highly regulated gene expression pattern, which in case of alterations, leads to male infertility. Changes of gene expression by chromosome anomalies, gene variants, and epigenetic alterations have been described as the main genetic causes of male infertility. Recent molecular and cytogenetic approaches have revealed that higher order chromosome positioning is essential for basic genome functions, including gene expression. This review addresses this issue by exposing well-founded evidences which support that alterations on the chromosome topology in spermatogenetic cells leads to defective sperm function and could be considered as an additional genetic cause of male infertility.