Chromatin Conformation Links Distal Target Genes to CKD Loci

Chromatin accessibility analysis and architectural profiling of human kidneys reveal key cell types and a regulator of diabetic kidney disease

Diabetes is the leading cause of kidney disease that progresses to kidney failure. However, the key molecular and cellular pathways involved in diabetic kidney disease (DKD) pathogenesis are largely unknown. Here, we performed a comparative analysis of adult human kidneys by examining cell type-specific chromatin accessibility by single-nucleus ATAC-seq (snATAC-seq) and analyzing three-dimensional chromatin architecture via high-throughput chromosome conformation capture (Hi-C method) of paired samples. We mapped the cell type-specific and DKD-specific open chromatin landscape and found that genetic variants associated with kidney diseases were significantly enriched in the proximal tubule- (PT) and injured PT-specific open chromatin regions in samples from patients with DKD. BACH1 was identified as a core transcription factor of injured PT cells; its binding target genes were highly associated with fibrosis and inflammation, which were also key features of injured PT cells. Additionally, Hi-C analysis revealed global chromatin architectural changes in DKD, accompanied by changes in local open chromatin patterns. Combining the snATAC-seq and Hi-C data identified direct target genes of BACH1, and indicated that BACH1 binding regions showed increased chromatin contact frequency with promoters of their target genes in DKD. Thus, our multi-omics analysis revealed BACH1 target genes in injured PTs and highlighted the role of BACH1 as a novel regulator of tubular inflammation and fibrosis.

Challenges and considerations for reproducibility of STARR-seq assays

High-throughput methods such as RNA-seq, ChIP-seq, and ATAC-seq have well-established guidelines, commercial kits, and analysis pipelines that enable consistency and wider adoption for understanding genome function and regulation. STARR-seq, a popular assay for directly quantifying the activities of thousands of enhancer sequences simultaneously, has seen limited standardization across studies. The assay is long, with more than 250 steps, and frequent customization of the protocol and variations in bioinformatics methods raise concerns for reproducibility of STARR-seq studies. Here, we assess each step of the protocol and analysis pipelines from published sources and in-house assays, and identify critical steps and quality control (QC) checkpoints necessary for reproducibility of the assay. We also provide guidelines for experimental design, protocol scaling, customization, and analysis pipelines for better adoption of the assay. These resources will allow better optimization of STARR-seq for specific research needs, enable comparisons and integration across studies, and improve the reproducibility of results.

Exploring high-resolution chromatin interaction changes and functional enhancers of myogenic marker genes during myogenic differentiation

Skeletal muscle differentiation (myogenesis) is a complex and highly coordinated biological process regulated by a series of myogenic marker genes. Chromatin interactions between gene's promoters and their enhancers have an important role in transcriptional control. However, the high-resolution chromatin interactions of myogenic genes and their functional enhancers during myogenesis remain largely unclear. Here, we used circularized chromosome conformation capture coupled with next generation sequencing (4C-seq) to investigate eight myogenic marker genes in C2C12 myoblasts (C2C12-MBs) and C2C12 myotubes (C2C12-MTs). We revealed dynamic chromatin interactions of these marker genes during differentiation and identified 163 and 314 significant interaction sites (SISs) in C2C12-MBs and C2C12-MTs, respectively. The interacting genes of SISs in C2C12-MTs were mainly involved in muscle development, and histone modifications of the SISs changed during differentiation. Through functional genomic screening, we also identified 25 and 41 putative active enhancers in C2C12-MBs and C2C12-MTs, respectively. Using luciferase reporter assays for putative enhancers of Myog and Myh3, we identified eight activating enhancers. Furthermore, dCas9-KRAB epigenome editing and RNA-Seq revealed a role for Myog enhancers in the regulation of Myog expression and myogenic differentiation in the native genomic context. Taken together, this study lays the groundwork for understanding 3D chromatin interaction changes of myogenic genes during myogenesis and provides insights that contribute to our understanding of the role of enhancers in regulating myogenesis.

Therapies Targeting Epigenetic Alterations in Acute Kidney Injury-to-Chronic Kidney Disease Transition

Acute kidney injury (AKI) was previously thought to be a merely transient event; however, recent epidemiological evidence supports the existence of a causal relationship between AKI episodes and subsequent progression to chronic kidney disease (CKD). Although the pathophysiology of this AKI-to-CKD transition is not fully understood, it is mediated by the interplay among multiple components of the kidney including tubular epithelial cells, endothelial cells, pericytes, inflammatory cells, and myofibroblasts. Epigenetic alterations including histone modification, DNA methylation, non-coding RNAs, and chromatin conformational changes, are also expected to be largely involved in the pathophysiology as a “memory” of the initial injury that can persist and predispose to chronic progression of fibrosis. Each epigenetic modification has a great potential as a therapeutic target of AKI-to-CKD transition; timely and target-specific epigenetic interventions to the various temporal stages of AKI-to-CKD transition will be the key to future therapeutic applications in clinical practice. This review elaborates on the latest knowledge of each mechanism and the currently available therapeutic agents that target epigenetic modification in the context of AKI-to-CKD transition. Further studies will elucidate more detailed mechanisms and novel therapeutic targets of AKI-to-CKD transition.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

3D disorganization and rearrangement of genome provide insights into pathogenesis of NAFLD by integrated Hi-C, Nanopore, and RNA sequencing

The three-dimensional (3D) conformation of chromatin is integral to the precise regulation of gene expression. The 3D genome and genomic variations in non-alcoholic fatty liver disease (NAFLD) are largely unknown, despite their key roles in cellular function and physiological processes. High-throughput chromosome conformation capture (Hi-C), Nanopore sequencing, and RNA-sequencing (RNA-seq) assays were performed on the liver of normal and NAFLD mice. A high-resolution 3D chromatin interaction map was generated to examine different 3D genome hierarchies including A/B compartments, topologically associated domains (TADs), and chromatin loops by Hi-C, and whole genome sequencing identifying structural variations (SVs) and copy number variations (CNVs) by Nanopore sequencing. We identified variations in thousands of regions across the genome with respect to 3D chromatin organization and genomic rearrangements, between normal and NAFLD mice, and revealed gene dysregulation frequently accompanied by these variations. Candidate target genes were identified in NAFLD, impacted by genetic rearrangements and spatial organization disruption. Our data provide a high-resolution 3D genome interaction resource for NAFLD investigations, revealed the relationship among genetic rearrangements, spatial organization disruption, and gene regulation, and identified candidate genes associated with these variations implicated in the pathogenesis of NAFLD. The newly findings offer insights into novel mechanisms of NAFLD pathogenesis and can provide a new conceptual framework for NAFLD therapy.

Chromatin Conformation in Development and Disease

Chromatin domains and loops are important elements of chromatin structure and dynamics, but much remains to be learned about their exact biological role and nature. Topological associated domains and functional loops are key to gene expression and hold the answer to many questions regarding developmental decisions and diseases. Here, we discuss new findings, which have linked chromatin conformation with development, differentiation and diseases and hypothesized on various models while integrating all recent findings on how chromatin architecture affects gene expression during development, evolution and disease.

Islet Epigenetic Impacts on β-Cell Identity and Function

The development and maintenance of differentiation is vital to the function of mature cells. Terminal differentiation is achieved by locking in the expression of genes essential for the function of those cells. Gene expression and its memory through generations of cell division is controlled by transcription factors and a host of epigenetic marks. In type 2 diabetes, β cells have altered gene expression compared to controls, accompanied by altered chromatin marks. Mutations, diet, and environment can all disrupt the implementation and preservation of the distinctive β-cell transcriptional signature. Understanding of the full complement of genomic control in β cells is still nascent. This article describes the known effects of histone marks and variants, DNA methylation, how they are regulated in the β cell, and how they affect cell-fate specification, maintenance, and lineage propagation. © 2021 American Physiological Society. Compr Physiol 11:1-18, 2021.

Chromatin architecture reveals cell type-specific target genes for kidney disease risk variants

BACKGROUND

Cell type-specific transcriptional programming results from the combinatorial interplay between the repertoire of active regulatory elements. Disease-associated variants disrupt such programming, leading to altered expression of downstream regulated genes and the onset of pathological states. However, due to the non-linear regulatory properties of non-coding elements such as enhancers, which can activate transcription at long distances and in a non-directional way, the identification of causal variants and their target genes remains challenging. Here, we provide a multi-omics analysis to identify regulatory elements associated with functional kidney disease variants, and downstream regulated genes.

RESULTS

In order to understand the genetic risk of kidney diseases, we generated a comprehensive dataset of the chromatin landscape of human kidney tubule cells, including transcription-centered 3D chromatin organization, histone modifications distribution and transcriptome with HiChIP, ChIP-seq and RNA-seq. We identified genome-wide functional elements and thousands of interactions between the distal elements and target genes. The results revealed that risk variants for renal tumor and chronic kidney disease were enriched in kidney tubule cells. We further pinpointed the target genes for the variants and validated two target genes by CRISPR/Cas9 genome editing techniques in zebrafish, demonstrating that SLC34A1 and MTX1 were indispensable genes to maintain kidney function.

CONCLUSIONS

Our results provide a valuable multi-omics resource on the chromatin landscape of human kidney tubule cells and establish a bioinformatic pipeline in dissecting functions of kidney disease-associated variants based on cell type-specific epigenome.

Extracellular Matrix Analysis of Human Renal Arteries in Both Quiescent and Active Vascular State

In vascular tissue engineering strategies, the addition of vascular-specific extracellular matrix (ECM) components may better mimic the in vivo microenvironment and potentially enhance cell-matrix interactions and subsequent tissue growth. For this purpose, the exact composition of the human vascular ECM first needs to be fully characterized. Most research has focused on characterizing ECM components in mature vascular tissue; however, the developing fetal ECM matches the active environment required in vascular tissue engineering more closely. Consequently, we characterized the ECM protein composition of active (fetal) and quiescent (mature) renal arteries using a proteome analysis of decellularized tissue. The obtained human fetal renal artery ECM proteome dataset contains higher levels of 15 ECM proteins versus the mature renal artery ECM proteome, whereas 16 ECM proteins showed higher levels in the mature tissue compared to fetal. Elastic ECM proteins EMILIN1 and FBN1 are significantly enriched in fetal renal arteries and are mainly produced by cells of mesenchymal origin. We functionally tested the role of EMILIN1 and FBN1 by anchoring the ECM secreted by vascular smooth muscle cells (SMCs) to glass coverslips. This ECM layer was depleted from either EMILIN1 or FBN1 by using siRNA targeting of the SMCs. Cultured endothelial cells (ECs) on this modified ECM layer showed alterations on the transcriptome level of multiple pathways, especially the Rho GTPase controlled pathways. However, no significant alterations in adhesion, migration or proliferation were observed when ECs were cultured on EMILIN1- or FNB1-deficient ECM. To conclude, the proteome analysis identified unique ECM proteins involved in the embryonic development of renal arteries. Alterations in transcriptome levels of ECs cultured on EMILIN1- or FBN1-deficient ECM showed that these candidate proteins could affect the endothelial (regenerative) response.

Three-dimensional genome: developmental technologies and applications in precision medicine

In the 20th century, our familiar structure of DNA was the double helix. Due to technical limitations, we do not have a good way to understand the finer structure of the genome, let alone its transcriptional regulation. Until the advent of 3C technologies, we were no longer blind to this one. Three-dimensional (3D) genomics is a new subject, which mainly studies the 3D structure and transcriptional regulation of eukaryotic genomes. Now, this field mainly has Hi-C series and CHIA-PET series technologies. Through 3D genomics, we can understand the basic structure of DNA, understand the growth and development of organisms and the occurrence of diseases, so as to promote human medical and health undertakings. The review introduces the main research techniques of 3D genomics and their characteristics, the latest development of 3D genome structure, the relationship between diseases and 3D genome structure, the applications of 3D genome in precision medicine, and the development of the 4D nucleome project.

Transcriptome analysis reveals microvascular endothelial cell-dependent pericyte differentiation

Microvascular homeostasis is strictly regulated, requiring close interaction between endothelial cells and pericytes. Here, we aimed to improve our understanding of how microvascular crosstalk affects pericytes. Human-derived pericytes, cultured in absence, or presence of human endothelial cells, were studied by RNA sequencing. Compared with mono-cultured pericytes, a total of 6704 genes were differentially expressed in co-cultured pericytes. Direct endothelial contact induced transcriptome profiles associated with pericyte maturation, suppression of extracellular matrix production, proliferation, and morphological adaptation. In vitro studies confirmed enhanced pericyte proliferation mediated by endothelial-derived PDGFB and pericyte-derived HB-EGF and FGF2. Endothelial-induced PLXNA2 and ACTR3 upregulation also triggered pericyte morphological adaptation. Pathway analysis predicted a key role for TGFβ signaling in endothelial-induced pericyte differentiation, whereas the effect of signaling via gap- and adherens junctions was limited. We demonstrate that endothelial cells have a major impact on the transcriptional profile of pericytes, regulating endothelial-induced maturation, proliferation, and suppression of ECM production.

QTL mapping of rat blood pressure loci on RNO1 within a homologous region linked to human hypertension on HSA15

Fine-mapping of regions linked to the inheritance of hypertension is accomplished by genetic dissection of blood pressure quantitative trait loci (BP QTLs) in rats. The goal of the current study was to further fine-map two genomic regions on rat chromosome 1 with opposing blood pressure effects (BP QTL1b1 and BP QTL1b1a), the homologous region of which on human chromosome 15 harbors BP QTLs. Two new substrains were constructed and studied from the previously reported BP QTL1b1, one having significantly lower systolic BP by 17 mmHg than that of the salt-sensitive (S) rat (P = 0.007). The new limits of BP QTL1b1 were between 134.09 Mb and 135.40 Mb with a 43% improvement from the previous 2.31 Mb to the current 1.31 Mb interval containing 4 protein-coding genes (Rgma, Chd2, Fam174b, and St8sia2), 2 predicted miRNAs, and 4 lncRNAs. One new substrain was constructed and studied from the previously reported BPQTL1b1a having a significantly higher systolic BP by 22 mmHg (P = 0.006) than that of the S rat. The new limits of BPQTL1b1a were between 133.53 Mb and 134.52 Mb with a 32% improvement from the previous1.45 Mb to the current 990.21 Kb interval containing 1 protein-coding gene, Mctp2, and a lncRNA. The congenic segments of these two BP QTLs overlapped between 134.09 Mb and 134.52 Mb. No exonic variants were detected in any of the genes. These findings reiterate complexity of genetic regulation of BP within QTL regions, where elements beyond protein-coding sequences could be factors in controlling BP.

Indoxyl Sulfate Stimulates Angiogenesis by Regulating Reactive Oxygen Species Production via CYP1B1

Indoxyl sulfate (IS) is an accumulative protein-bound uremic toxin found in patients with kidney disease. It is reported that IS impairs the vascular endothelium, but a comprehensive overview of all mechanisms active in IS-injury currently remains lacking. Here we performed RNA sequencing in human umbilical vein endothelial cells (HUVECs) after IS or control medium treatment and identified 1293 genes that were affected in a IS-induced response. Gene enrichment analysis highlighted pathways involved in altered vascular formation and cell metabolism. We confirmed these transcriptome profiles at the functional level by demonstrating decreased viability and increased cell senescence in response to IS treatment. In line with the additional pathways highlighted by the transcriptome analysis, we further could demonstrate that IS exposure of HUVECs promoted tubule formation as shown by the increase in total tubule length in a 3D HUVECs/pericytes co-culture assay. Notably, the pro-angiogenic response of IS and increased ROS production were abolished when CYP1B1, one of the main target genes that was highly upregulated by IS, was silenced. This observation indicates IS-induced ROS in endothelial cells is CYP1B1-dependent. Taken together, our findings demonstrate that IS promotes angiogenesis and CYP1B1 is an important factor in IS-activated angiogenic response.

Understanding the 3D genome: Emerging impacts on human disease

Recent burst of new technologies that allow for quantitatively delineating chromatin structure has greatly expanded our understanding of how the genome is organized in the three-dimensional (3D) space of the nucleus. It is now clear that the hierarchical organization of the eukaryotic genome critically impacts nuclear activities such as transcription, replication, as well as cellular and developmental events such as cell cycle, cell fate decision and embryonic development. In this review, we discuss new insights into how the structural features of the 3D genome hierarchy are established and maintained, how this hierarchy undergoes dynamic rearrangement during normal development and how its perturbation will lead to human disease, highlighting the accumulating evidence that links the diverse 3D genome architecture components to a multitude of human diseases and the emerging mechanisms by which 3D genome derangement causes disease phenotypes.

Non-coding DNA in IBD: from sequence variation in DNA regulatory elements to novel therapeutic potential

Genome-wide association studies have identified over 200 loci associated with IBD. We and others have recently shown that, in addition to variants in protein-coding genes, the majority of the associated loci are related to DNA regulatory elements (DREs). These findings add a dimension to the already complex genetic background of IBD. In this review we summarise the existing evidence on the role of DREs in IBD. We discuss how epigenetic research can be used in candidate gene approaches that take non-coding variants into account and can help to pinpoint the essential pathways and cell types in the pathogenesis of IBD. Despite the increased level of genetic complexity, these findings can contribute to novel therapeutic options that target transcription factor binding and enhancer activity. Finally, we summarise the future directions and challenges of this emerging field.

Integrated Functional Genomic Analysis Enables Annotation of Kidney Genome-Wide Association Study Loci

BACKGROUND

Linking genetic risk loci identified by genome-wide association studies (GWAS) to their causal genes remains a major challenge. Disease-associated genetic variants are concentrated in regions containing regulatory DNA elements, such as promoters and enhancers. Although researchers have previously published DNA maps of these regulatory regions for kidney tubule cells and glomerular endothelial cells, maps for podocytes and mesangial cells have not been available.

METHODS

We generated regulatory DNA maps (DNase-seq) and paired gene expression profiles (RNA-seq) from primary outgrowth cultures of human glomeruli that were composed mainly of podocytes and mesangial cells. We generated similar datasets from renal cortex cultures, to compare with those of the glomerular cultures. Because regulatory DNA elements can act on target genes across large genomic distances, we also generated a chromatin conformation map from freshly isolated human glomeruli.

RESULTS

We identified thousands of unique regulatory DNA elements, many located close to transcription factor genes, which the glomerular and cortex samples expressed at different levels. We found that genetic variants associated with kidney diseases (GWAS) and kidney expression quantitative trait loci were enriched in regulatory DNA regions. By combining GWAS, epigenomic, and chromatin conformation data, we functionally annotated 46 kidney disease genes.

CONCLUSIONS

We demonstrate a powerful approach to functionally connect kidney disease-/trait-associated loci to their target genes by leveraging unique regulatory DNA maps and integrated epigenomic and genetic analysis. This process can be applied to other kidney cell types and will enhance our understanding of genome regulation and its effects on gene expression in kidney disease.

Genome-Wide Maps of Transcription Regulatory Elements and Transcription Enhancers in Development and Disease

Gene expression is regulated by numerous elements including enhancers, insulators, transcription factors, and architectural proteins. Regions of DNA distal to the transcriptional start site, called enhancers, play a central role in the temporal and tissue-specific regulation of gene expression through RNA polymerase II. The identification of enhancers and other cis regulatory elements has largely been possible due to advances in next generation sequencing technologies. Enhancers regulate gene expression through chromatin loops mediated by architectural proteins such as YY1, CTCF, the cohesin complex, and LDB1. Additionally, enhancers can be transcribed to produce noncoding RNAs termed enhancer RNAs that likely participate in transcriptional regulation. The central role of enhancers in regulating gene expression implicates them in both normal physiology but also many disease states. The importance of enhancers is evident by the suggested role of SNPs, duplications, and other alterations of enhancer function in many diseases, ranging from cancer to atherosclerosis to chronic kidney disease. Although much progress has been made in recent years, the field of enhancer biology and our knowledge of the cis regulome remains a work in progress. This review will highlight recent seminal studies which demonstrate the role of enhancers in normal physiology and disease pathogenesis. © 2019 American Physiological Society. Compr Physiol 9:439-455, 2019.