Chromosome conformation capture

UBA1-CDK16 : A Sex-Specific Chimeric RNA and Its Role in Immune Sexual Dimorphism

RNA processing mechanisms, such as alternative splicing and RNA editing, have been recognized as critical means to expand the transcriptome. Chimeric RNAs formed by intergenic splicing provide another potential layer of RNA diversification. By analyzing a large set of RNA-Seq data and validating results in over 1,200 blood samples, we identified UBA1-CDK16 , a female-specific chimeric transcript. Intriguingly, both parental genes, are expressed in males and females. Mechanistically, UBA1-CDK16 is produced by cis-splicing between the two adjacent X-linked genes, originating from the inactive X chromosome. A female-specific chromatin loop, formed between the junction sites, facilitates the alternative splicing of its readthrough precursor. This unique chimeric transcript exhibits evolutionary conservation, evolving to be female-specific from non-human primates to humans. Furthermore, our investigation reveals that UBA1-CDK16 is enriched in the myeloid lineage and plays a regulatory role in myeloid differentiation. Notably, female COVID-19 patients who tested negative for this chimeric transcript displayed higher counts of neutrophils, highlighting its potential role in disease pathogenesis. These findings support the notion that chimeric RNAs represent a new repertoire of transcripts that can be regulated independently from the parental genes, and a new class of RNA variance with potential implications in sexual dimorphism and immune responses.

Transgenerational Epigenetic Inheritance of MHC Class I Gene Expression is Regulated by the CCAAT Promoter Element

Transgenerational epigenetic inheritance is defined as the transmission of traits or gene expression patterns across multiple generations that do not derive from DNA alterations. The effect of multiple stress factors or metabolic changes resulting in such inheritance have been documented in plants, worms and flies and mammals. The molecular basis for epigenetic inheritance has been linked to histone and DNA modifications and non-coding RNA. In this study, we show that mutation of a promoter element, the CCAAT box, disrupts stable expression of an MHC Class I transgene, resulting in variegated expression among progeny for at least 4 generations in multiple independently derived transgenic lines. Histone modifications and RNA polII binding correlate with expression, whereas DNA methylation and nucleosome occupancy do not. Mutation of the CCAAT box abrogates NF-Y binding and results in changes to CTCF binding and DNA looping patterns across the gene that correlate with expression status from one generation to the next. These studies identify the CCAAT promoter element as a regulator of stable transgenerational epigenetic inheritance. Considering that the CCAAT box is present in 30% of eukaryotic promoters, this study could provide important insights into how fidelity of gene expression patterns is maintained through multiple generations.

Lymphocyte activation gene-3-associated protein networks are associated with HDL-cholesterol and mortality in the Trans-omics for Precision Medicine program

Deficiency of the immune checkpoint lymphocyte activation gene-3 (LAG3) protein is significantly associated with both elevated HDL-cholesterol (HDL-C) and myocardial infarction risk. We determined the association of genetic variants within ±500 kb of LAG3 with plasma LAG3 and defined LAG3-associated plasma proteins with HDL-C and clinical outcomes. Whole genome sequencing and plasma proteomics were obtained from the Multi-Ethnic Study of Atherosclerosis (MESA) and the Framingham Heart Study (FHS) cohorts as part of the Trans-Omics for Precision Medicine program. In situ Hi-C chromatin capture was performed in EBV-transformed cell lines isolated from four MESA participants. Genetic association analyses were performed in MESA using multivariate regression models, with validation in FHS. A LAG3-associated protein network was tested for association with HDL-C, coronary heart disease, and all-cause mortality. We identify an association between the LAG3 rs3782735 variant and plasma LAG3 protein. Proteomics analysis reveals 183 proteins significantly associated with LAG3 with four proteins associated with HDL-C. Four proteins discovered for association with all-cause mortality in FHS shows nominal associations in MESA. Chromatin capture analysis reveals significant cis interactions between LAG3 and C1S, LRIG3, TNFRSF1A, and trans interactions between LAG3 and B2M. A LAG3-associated protein network has significant associations with HDL-C and mortality.

Rodriguez et al. use whole genome sequencing and plasma proteomics from the Multi-Ethnic Study of Atherosclerosis (MESA) and the Framingham Heart Study (FHS) cohorts of the Trans-Omics for Precision Medicine program and perform in situ Hi-C chromatin capture in cell lines isolated from four MESA participants. They demonstrate that lymphocyte activation gene-3 protein networks are associated with HDL-cholesterol and mortality, which could guide the development of precision medicine.

Diet-dependent natriuretic peptide receptor C expression in adipose tissue is mediated by PPARγ via long-range distal enhancers

The cardiac natriuretic peptides (NPs) are well established as regulators of blood pressure and fluid volume, but they also stimulate adipocyte lipolysis and control the gene program of nonshivering thermogenesis in brown adipose tissue. The NP "clearance" receptor C (NPRC) functions to clear NPs from the circulation via peptide internalization and degradation and thus is an important regulator of NP signaling and adipocyte metabolism. It is well known that the Nprc gene is highly expressed in adipose tissue and dynamically regulated upon nutrition and environmental changes. However, the molecular basis for how Nprc gene expression is regulated is still poorly understood. Here, we identified the nuclear receptor transcription factor peroxisome proliferator-activated receptor gamma (PPARγ) as a transcriptional regulator of Nprc expression in mouse adipocytes. During 3T3-L1 adipocyte differentiation, levels of Nprc expression increase in parallel with PPARγ induction. Rosiglitazone, a classic PPARγ agonist, increases, whereas siRNA knockdown of PPARγ reduces, Nprc expression in 3T3-L1 adipocytes. By using chromosome conformation capture and luciferase reporter assays, we demonstrate that PPARγ controls Nprc gene expression in adipocytes through its long-range distal enhancers. Furthermore, the induction of Nprc expression in adipose tissue during high-fat diet feeding is found to be associated with increased PPARγ enhancer activity. Our findings define PPARγ as a mediator of adipocyte Nprc gene expression and establish a new connection between PPARγ and the control of adipocyte NP signaling in obesity.

High HDL-Cholesterol Paradox: SCARB1-LAG3-HDL Axis

Purpose of the Review

To evaluate recent studies related to the paradox of high HDL-C with mortality and atherosclerotic cardiovascular disease (ASCVD) risk.

Recent Findings

Two observational studies (Cardiovascular Health in Ambulatory Care Research Team [CANHEART] and Copenhagen City Heart Study and the Copenhagen General Population Study [Copenhagen Heart Studies]) of adults without pre-existing ASCVD have shown a significant U-shaped association of HDL-C with all-cause and cause-specific mortality. Both studies showed that low HDL-C levels consistently increased hazard risk (HR) for all-cause and cause-specific mortality. In the CANHEART study, high HDL-C levels, HDL-C > 90 mg/dL, were associated with increased HR for non-CVD/non-cancer mortality. In the Copenhagen Heart Studies, women with HDL-C ≥ 135 mg/dL showed increased HR for all-cause and CVD mortality, while men with HDL-C > 97 mg/dL showed increased HR for all-cause and CVD mortality. Genetic association studies failed to show that genetic etiologies of high HDL-C significantly reduced risk for myocardial infarction (MI), while hepatocyte nuclear factor-4 (HNF4A) was significantly associated with high HDL-C and increased MI risk. Candidate gene studies have identified scavenger receptor B class I (SCARB1) and lymphocyte activation gene-3 (LAG3) as genes significantly associated with high HDL-C and increased MI risk.

Summary

Low HDL-C remains as a significant factor for increased disease risk while high HDL-C levels are not associated with cardioprotection. Clinical CVD risk calculators need revision.

Conjugated activation of myocardial-specific transcription of Gja5 by a pair of Nkx2-5-Shox2 co-responsive elements

The sinoatrial node (SAN) is the primary pacemaker in the heart. During cardiogenesis, Shox2 and Nkx2-5 are co-expressed in the junction domain of the SAN and regulate pacemaker cell fate through a Shox2-Nkx2-5 antagonism. Cx40 is a marker of working myocardium and an Nkx2-5 transcriptional output antagonized by Shox2, but the underlying regulatory mechanisms remain elusive. Here we characterized a bona fide myocardial-specific Gja5 (coding gene of Cx40) distal enhancer consisting of a pair of Nkx2-5 and Shox2 co-bound elements in the regulatory region of Gja5. Transgenic reporter assays revealed that neither element alone, but the conjugation of both elements together, drives myocardial-specific transcription. Genetic analyses confirmed that the activation of this enhancer depends on Nkx2-5 but is inhibited by Shox2 in vivo, and its presence is essential for Gja5 expression in the myocardium but not the endothelial cells of the heart. Furthermore, chromatin conformation analysis showed an Nkx2-5-dependent loop formation between these two elements and the Gja5 promoter in vivo, indicating that Nkx2-5 bridges the conjugated activation of this enhancer by pairing the two elements to the Gja5 promoter.

An Allele-Specific Functional SNP Associated with Two Systemic Autoimmune Diseases Modulates IRF5 Expression by Long-Range Chromatin Loop Formation

Both systemic lupus erythematosus (SLE) and systemic sclerosis (SSc) are autoimmune diseases sharing similar genetic backgrounds. Genome-wide association studies have constantly disclosed numerous genetic variants conferring to both disease risks at 7q32.1, but the functional mechanisms underlying them are still largely unknown. Through a series of bioinformatics and functional analyses, we prioritized a potential independent functional single-nucleotide polymorphism (rs13239597) within TNPO3 promoter region, residing in a putative enhancer element and validated that IRF5 is the distal target gene (∼118 kb) of rs13239597, which is a key regulator involved in pathogenic autoantibody dysregulation, increasing risk of both SLE and SSc. We experimentally validated the long-range chromatin interactions between rs13239597 and IRF5 using chromosome conformation capture assay. We further demonstrated that rs13239597-A acted as an allele-specific enhancer regulating IRF5 expression, independently of TNPO3 by using dual-luciferase reporter assays and CRISPR-Cas9. Particularly, the transcription factor EVI1 could preferentially bind to rs13239597-A allele and increase the enhancer activity to regulate IRF5 expression. Taken together, our results uncovered a mechanistic insight of a noncoding functional variant acting as an allele-specific distal enhancer to directly modulate IRF5 expression, which might obligate in understanding of complex genetic architectures of SLE and SSc pathogenesis.

Dysregulation of a long noncoding RNA reduces leptin leading to a leptin-responsive form of obesity

Quantitative changes in leptin concentration lead to alterations in food intake and body weight, but the regulatory mechanisms that control leptin gene expression are poorly understood. Here we report that fat-specific and quantitative leptin expression is controlled by redundant cis elements and trans factors interacting with the proximal promoter together with a long noncoding RNA (lncOb). Diet-induced obese mice lacking lncOb show increased fat mass with reduced plasma leptin levels and lose weight after leptin treatment, whereas control mice do not. Consistent with this finding, large-scale genetic studies of humans reveal a significant association of single-nucleotide polymorphisms (SNPs) in the region of human lncOb with lower plasma leptin levels and obesity. These results show that reduced leptin gene expression can lead to a hypoleptinemic, leptin-responsive form of obesity and provide a framework for elucidating the pathogenic mechanism in the subset of obese patients with low endogenous leptin levels.

GrapHi-C: graph-based visualization of Hi-C datasets

Objectives

Hi-C is a proximity-based ligation reaction used to detect regions of the genome that are close in 3D space (or “interacting”). Typically, results from Hi-C experiments (contact maps) are visualized as heatmaps or Circos plots. While informative, these visualizations do not directly represent genomic structure and folding, making the interpretation of the underlying 3D genomic organization obscured. Our objective was to generate a graph-based contact map representation that leads to a more intuitive structural visualization.

Results

Normalized contact maps were converted into undirected graphs where each vertex represented a genomic region and each edge represented a detected (intra- and inter-chromosomal) or known (linear) interaction between two regions. Each edge was weighted by the inverse of the linear distance (Hi-C experimental resolution) or the interaction frequency from the contact map. Graphs were generated based on this representation scheme for contact maps from existing fission yeast datasets. Originally, these datasets were used to (1) identify specific principles influencing fission yeast genome organization and (2) uncover changes in fission yeast genome organization during the cell cycle. When compared to the equivalent heatmaps and/or Circos plots, the graph-based visualizations more intuitively depicted the changes in genome organization described in the original studies.

Electronic supplementary material

The online version of this article (10.1186/s13104-018-3507-2) contains supplementary material, which is available to authorized users.

Mapping long-range promoter contacts in human cells with high-resolution capture Hi-C

Transcriptional control in large genomes often requires looping interactions between distal DNA elements, such as enhancers and target promoters. Current chromosome conformation capture techniques do not offer sufficiently high resolution to interrogate these regulatory interactions on a genomic scale. Here we use Capture Hi-C (CHi-C), an adapted genome conformation assay, to examine the long-range interactions of almost 22,000 promoters in 2 human blood cell types. We identify over 1.6 million shared and cell type-restricted interactions spanning hundreds of kilobases between promoters and distal loci. Transcriptionally active genes contact enhancer-like elements, whereas transcriptionally inactive genes interact with previously uncharacterized elements marked by repressive features that may act as long-range silencers. Finally, we show that interacting loci are enriched for disease-associated SNPs, suggesting how distal mutations may disrupt the regulation of relevant genes. This study provides new insights and accessible tools to dissect the regulatory interactions that underlie normal and aberrant gene regulation.

A rare functional noncoding variant at the GWAS-implicated MIR137/MIR2682 locus might confer risk to schizophrenia and bipolar disorder

Schizophrenia (SZ) genome-wide association studies (GWASs) have identified common risk variants in >100 susceptibility loci; however, the contribution of rare variants at these loci remains largely unexplored. One of the strongly associated loci spans MIR137 (miR137) and MIR2682 (miR2682), two microRNA genes important for neuronal function. We sequenced ∼6.9 kb MIR137/MIR2682 and upstream regulatory sequences in 2,610 SZ cases and 2,611 controls of European ancestry. We identified 133 rare variants with minor allele frequency (MAF) <0.5%. The rare variant burden in promoters and enhancers, but not insulators, was associated with SZ (p = 0.021 for MAF < 0.5%, p = 0.003 for MAF < 0.1%). A rare enhancer SNP, 1:g.98515539A>T, presented exclusively in 11 SZ cases (nominal p = 4.8 × 10(-4)). We further identified its risk allele T in 2 of 2,434 additional SZ cases, 11 of 4,339 bipolar (BP) cases, and 3 of 3,572 SZ/BP study controls and 1,688 population controls; yielding combined p values of 0.0007, 0.0013, and 0.0001 for SZ, BP, and SZ/BP, respectively. The risk allele T of 1:g.98515539A>T reduced enhancer activity of its flanking sequence by >50% in human neuroblastoma cells, predicting lower expression of MIR137/MIR2682. Both empirical and computational analyses showed weaker transcription factor (YY1) binding by the risk allele. Chromatin conformation capture (3C) assay further indicated that 1:g.98515539A>T influenced MIR137/MIR2682, but not the nearby DPYD or LOC729987. Our results suggest that rare noncoding risk variants are associated with SZ and BP at MIR137/MIR2682 locus, with risk alleles decreasing MIR137/MIR2682 expression.

Transcriptional regulatory network shapes the genome structure of Saccharomyces cerevisiae

Among cellular processes gene transcription is central. More and more evidence is mounting that transcription is tightly connected with the spatial organization of the chromosomes. Spatial proximity of genes sharing transcriptional machinery is one of the consequences of this organization. Motivated by information on the physical relationship among genes identified via chromosomal conformation capture methods, we complement the spatial organization with the idea that genes under similar transcription factor control, but possible scattered throughout the genome, might be in physically proximity to facilitate the access of their commonly used transcription factors. Unlike the transcription factory model, "interacting" genes in our "Gene Proximity Model" are not necessarily immediate physical neighbors but are in spatial proximity. Considering the stochastic nature of TF-promoter binding, this local condensation mechanism could serve as a tie to recruit co-regulated genes to guarantee the swiftness of biological reactions. We tested this idea with a simple eukaryotic organism, Saccharomyces cerevisiae. Chromosomal interaction patterns and folding behavior generated by our model re-construct those obtained from experiments. We show that the transcriptional regulatory network has a close linkage with the genome organization in budding yeast, which is fundamental and instrumental to later studies on other more complex eukaryotes.

Peroxisome proliferator-activated receptor α positively regulates complement C3 expression but inhibits tumor necrosis factor α-mediated activation of C3 gene in mammalian hepatic-derived cells

Complement C3 is a pivotal component of three cascades of complement activation. The liver is the main source of C3 in circulation and expression and secretion of C3 by hepatocytes is increased during acute inflammation. However, the mechanism of the regulation of the C3 gene in hepatocytes is not well elucidated. We showed that the C3 gene is the direct target for peroxisome proliferator-activated receptor α (PPARα) in human hepatoma HepG2 cells and mouse liver. Using PPARα siRNA and synthetic PPARα agonist WY-14643 and antagonist MK886 we showed that activation of PPARα results in up-regulation of C3 gene expression and protein secretion by HepG2 cells. The PPAR response element (PPRE), which is able to bind PPARα in vitro and in vivo, was found in the human C3 promoter. PPRE is conserved between human and mouse, and WY-14643 stimulates mouse C3 expression in the liver. TNFα increases C3 gene via NF-κB and, to a lesser extent, MEK1/2 signaling pathways, whereas TNFα-mediated stimulation of C3 protein secretion depends on activation of MEK1/2, p38, and JNK in HepG2 cells. Activation of PPARα abolishes TNFα-mediated up-regulation of C3 gene expression and protein secretion due to interference with NF-κB via PPRE-dependent mechanism in HepG2 cells. TNFα decreases PPARα protein content via NF-κB and MEK1/2 signaling pathways and inhibits PPARα binding with the human C3 promoter in HepG2 cells. These results suggest novel mechanism controlling C3 expression in hepatocytes during acute phase inflammation and demonstrate a crosstalk between PPARα and TNFα in the regulation of complement system.

Pairing of homologous regions in the mouse genome is associated with transcription but not imprinting status

Although somatic homologous pairing is common in Drosophila it is not generally observed in mammalian cells. However, a number of regions have recently been shown to come into close proximity with their homologous allele, and it has been proposed that pairing might be involved in the establishment or maintenance of monoallelic expression. Here, we investigate the pairing properties of various imprinted and non-imprinted regions in mouse tissues and ES cells. We find by allele-specific 4C-Seq and DNA FISH that the Kcnq1 imprinted region displays frequent pairing but that this is not dependent on monoallelic expression. We demonstrate that pairing involves larger chromosomal regions and that the two chromosome territories come close together. Frequent pairing is not associated with imprinted status or DNA repair, but is influenced by chromosomal location and transcription. We propose that homologous pairing is not exclusive to specialised regions or specific functional events, and speculate that it provides the cell with the opportunity of trans-allelic effects on gene regulation.

Long-range DNA looping and gene expression analyses identify DEXI as an autoimmune disease candidate gene

The chromosome 16p13 region has been associated with several autoimmune diseases, including type 1 diabetes (T1D) and multiple sclerosis (MS). CLEC16A has been reported as the most likely candidate gene in the region, since it contains the most disease-associated single-nucleotide polymorphisms (SNPs), as well as an imunoreceptor tyrosine-based activation motif. However, here we report that intron 19 of CLEC16A, containing the most autoimmune disease-associated SNPs, appears to behave as a regulatory sequence, affecting the expression of a neighbouring gene, DEXI. The CLEC16A alleles that are protective from T1D and MS are associated with increased expression of DEXI, and no other genes in the region, in two independent monocyte gene expression data sets. Critically, using chromosome conformation capture (3C), we identified physical proximity between the DEXI promoter region and intron 19 of CLEC16A, separated by a loop of >150 kb. In reciprocal experiments, a 20 kb fragment of intron 19 of CLEC16A, containing SNPs associated with T1D and MS, as well as with DEXI expression, interacted with the promotor region of DEXI but not with candidate DNA fragments containing other potential causal genes in the region, including CLEC16A. Intron 19 of CLEC16A is highly enriched for transcription-factor-binding events and markers associated with enhancer activity. Taken together, these data indicate that although the causal variants in the 16p13 region lie within CLEC16A, DEXI is an unappreciated autoimmune disease candidate gene, and illustrate the power of the 3C approach in progressing from genome-wide association studies results to candidate causal genes.

Preferential associations between co-regulated genes reveal a transcriptional interactome in erythroid cells

The discovery of interchromosomal interactions in higher eukaryotes points to a functional interplay between genome architecture and gene expression, challenging the view of transcription as a one-dimensional process. However, the extent of interchromosomal interactions and the underlying mechanisms are unknown. Here we present the first genome-wide analysis of transcriptional interactions using the mouse globin genes in erythroid tissues. Our results show that the active globin genes associate with hundreds of other transcribed genes, revealing extensive and preferential intra- and interchromosomal transcription interactomes. We show that the transcription factor Klf1 mediates preferential co-associations of Klf1-regulated genes at a limited number of specialized transcription factories. Our results establish a new gene expression paradigm, implying that active co-regulated genes and their regulatory factors cooperate to create specialized nuclear hot spots optimized for efficient and coordinated transcriptional control.