Comparison of Capture Hi-C Analytical Pipelines

Using paired-end read orientations to assess technical biases in capture Hi-C

Hi-C and capture Hi-C (CHi-C) both leverage paired-end sequencing of chimeric fragments to gauge the strength of interactions based on the total number of paired-end reads mapped to a common pair of restriction fragments. Mapped paired-end reads can have four relative orientations, depending on the genomic positions and strands of the two reads. We assigned one paired-end read orientation to each of the four possible re-ligations that can occur between two given restriction fragments. In a large hematopoietic cell dataset, we determined the read pair counts of interactions separately for each orientation. Interactions with imbalances in the counts occur much more often than expected by chance for both Hi-C and CHi-C. Based on such imbalances, we identified target restriction fragments enriched at only one instead of both ends. By matching them to the baits used for the experiments, we confirmed our assignment of paired-end read orientations and gained insights that can inform bait design. An analysis of unbaited fragments shows that, beyond bait effects, other known types of technical biases are reflected in count imbalances. Taking advantage of distance-dependent contact frequencies, we assessed the impact of such biases. Our results have the potential to improve the design and interpretation of CHi-C experiments.

Melanoma-specific mutation hotspots in distal, non-coding, promoter-interacting regions implicate novel candidate driver genes

BACKGROUND

To develop targeted treatments, it is crucial to identify the full spectrum of genetic drivers in melanoma, including those in non-coding regions. However, recent efforts to explore non-coding regions have primarily focused on gene-adjacent elements such as promoters and non-coding RNAs, leaving intergenic distal regulatory elements largely unexplored.

METHODS

We used Hi-C chromatin contact data from melanoma cells to map distal, non-coding, promoter-interacting regulatory elements genome-wide in melanoma. Using this "promoter-interaction network", alongside whole-genome sequence and gene expression data from the Pan Cancer Analysis of Whole Genomes, we developed multivariate linear regression models to identify distal somatic mutation hotspots that affect promoter activity.

RESULTS

We identified eight recurrently mutated hotspots that are novel, melanoma-specific, located in promoter-interacting distal regulatory elements, alter transcription factor binding motifs, and affect the expression of genes (e.g., HSPB7, CLDN1, ADCY9 and FDXR) previously implicated as tumour suppressors/oncogenes in various cancers.

CONCLUSIONS

Our study suggests additional non-coding drivers beyond the well-characterised TERT promoter in melanoma, offering new insights into the disruption of complex regulatory networks by non-coding mutations that may contribute to melanoma development. Furthermore, our study provides a framework for integrating multiple levels of biological data to uncover cancer-specific non-coding drivers.

Functional Genomics in Psoriasis

Psoriasis is an autoimmune cutaneous condition that significantly impacts quality of life and represents a burden on society due to its prevalence. Genome-wide association studies (GWASs) have pinpointed several psoriasis-related risk loci, underlining the disease's complexity. Functional genomics is paramount to unveiling the role of such loci in psoriasis and disentangling its complex nature. In this review, we aim to elucidate the main findings in this field and integrate our discussion with gold-standard techniques in molecular biology-i.e., Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-and high-throughput technologies. These tools are vital to understanding how disease risk loci affect gene expression in psoriasis, which is crucial in identifying new targets for personalized treatments in advanced precision medicine.

Revisiting Assessment of Computational Methods for Hi-C Data Analysis

The performances of algorithms for Hi-C data preprocessing, the identification of topologically associating domains, and the detection of chromatin interactions and promoter-enhancer interactions have been mostly evaluated using semi-quantitative or synthetic data approaches, without utilizing the most recent methods, since 2017. In this study, we comprehensively evaluated 24 popular state-of-the-art methods for the complete end-to-end pipeline of Hi-C data analysis, using manually curated or experimentally validated benchmark datasets, including a CRISPR dataset for promoter-enhancer interaction validation. Our results indicate that, although no single method exhibited superior performance in all situations, HiC-Pro, DomainCaller, and Fit-Hi-C2 showed relatively balanced performances of most evaluation metrics for preprocessing, topologically associating domain identification, and chromatin interaction/promoter-enhancer interaction detection, respectively. The comprehensive comparison presented in this manuscript provides a reference for researchers to choose Hi-C analysis tools that best suit their needs.

3D genome organization links non-coding disease-associated variants to genes

Genome sequencing has revealed over 300 million genetic variations in human populations. Over 90% of variants are single nucleotide polymorphisms (SNPs), the remainder include short deletions or insertions, and small numbers of structural variants. Hundreds of thousands of these variants have been associated with specific phenotypic traits and diseases through genome wide association studies which link significant differences in variant frequencies with specific phenotypes among large groups of individuals. Only 5% of disease-associated SNPs are located in gene coding sequences, with the potential to disrupt gene expression or alter of the function of encoded proteins. The remaining 95% of disease-associated SNPs are located in non-coding DNA sequences which make up 98% of the genome. The role of non-coding, disease-associated SNPs, many of which are located at considerable distances from any gene, was at first a mystery until the discovery that gene promoters regularly interact with distal regulatory elements to control gene expression. Disease-associated SNPs are enriched at the millions of gene regulatory elements that are dispersed throughout the non-coding sequences of the genome, suggesting they function as gene regulation variants. Assigning specific regulatory elements to the genes they control is not straightforward since they can be millions of base pairs apart. In this review we describe how understanding 3D genome organization can identify specific interactions between gene promoters and distal regulatory elements and how 3D genomics can link disease-associated SNPs to their target genes. Understanding which gene or genes contribute to a specific disease is the first step in designing rational therapeutic interventions.