High-resolution genome topology of human retina uncovers super enhancer-promoter interactions at tissue-specific and multifactorial disease loci

Functional genomics in age-related macular degeneration: From genetic associations to understanding disease mechanisms

Genome-wide association studies have been remarkably successful in identifying genetic variants associated with age-related macular degeneration (AMD), demonstrating a strong genetic component largely driven by common variants. However, progress in translating these genetic findings into a deeper understanding of disease mechanisms and new therapies has been slow. Slow progress in this area can be attributed to limited knowledge of the functional impact of AMD-associated non-coding variants on gene function, the molecular mechanisms and cell types underlying disease. This review offers a comprehensive overview of functional genomics approaches to uncover the genetic, epigenetic, cellular and molecular mechanisms underlying AMD and outlines future directions for research.

Single-cell analysis of the epigenome and 3D chromatin architecture in the human retina

Most genetic risk variants linked to ocular diseases are non-protein coding and presumably contribute to disease through dysregulation of gene expression, however, deeper understanding of their mechanisms of action has been impeded by an incomplete annotation of the transcriptional regulatory elements across different retinal cell types. To address this knowledge gap, we carried out single-cell multiomics assays to investigate gene expression, chromatin accessibility, DNA methylome and 3D chromatin architecture in human retina, macula, and retinal pigment epithelium (RPE)/choroid. We identified 420,824 unique candidate regulatory elements and characterized their chromatin states in 23 sub-classes of retinal cells. Comparative analysis of chromatin landscapes between human and mouse retina cells further revealed both evolutionarily conserved and divergent retinal gene-regulatory programs. Leveraging the rapid advancements in deep-learning techniques, we developed sequence-based predictors to interpret non-coding risk variants of retina diseases. Our study establishes retina-wide, single-cell transcriptome, epigenome, and 3D genome atlases, and provides a resource for studying the gene regulatory programs of the human retina and relevant diseases.

Maf-family bZIP transcription factor NRL interacts with RNA-binding proteins and R-loops in retinal photoreceptors

RNA-binding proteins (RBPs) perform diverse functions including the regulation of chromatin dynamics and the coupling of transcription with RNA processing. However, our understanding of their actions in mammalian neurons remains limited. Using affinity purification, yeast-two-hybrid and proximity ligation assays, we identified interactions of multiple RBPs with neural retina leucine (NRL) zipper, a Maf-family transcription factor critical for retinal rod photoreceptor development and function. In addition to splicing, many NRL-interacting RBPs are associated with R-loops, which form during transcription and increase during photoreceptor maturation. Focusing on DHX9 RNA helicase, we demonstrate that its expression is modulated by NRL and that the NRL-DHX9 interaction is positively influenced by R-loops. ssDRIP-Seq analysis reveals both stranded and unstranded R-loops at distinct genomic elements, characterized by active and inactive epigenetic signatures and enriched at neuronal genes. NRL binds to both types of R-loops, suggesting an epigenetically independent function. Our findings suggest additional functions of NRL during transcription and highlight complex interactions among transcription factors, RBPs, and R-loops in regulating photoreceptor gene expression in the mammalian retina.

CTCF-mediated 3D chromatin sets up the gene expression program in the male germline

Spermatogenesis is a unidirectional differentiation process that generates haploid sperm, but how the gene expression program that directs this process is established is largely unknown. Here we determine the high-resolution three-dimensional (3D) chromatin architecture of mouse male germ cells during spermatogenesis and show that CTCF-mediated 3D chromatin dictates the gene expression program required for spermatogenesis. In undifferentiated spermatogonia, CTCF-mediated chromatin interactions between meiosis-specific super-enhancers (SEs) and their target genes precede activation of these SEs on autosomes. These meiotic SEs recruit the master transcription factor A-MYB (MYBL1) in meiotic spermatocytes, which strengthens their 3D contacts and instructs a burst of meiotic gene expression. We also find that at the mitosis-to-meiosis transition, the germline-specific Polycomb protein SCML2 facilitates the resolution of chromatin loops that are specific to mitotic spermatogonia. Moreover, SCML2 and A-MYB help shape the unique 3D chromatin organization of sex chromosomes during meiotic sex chromosome inactivation. We propose that CTCF-mediated 3D chromatin organization regulates epigenetic priming that directs unidirectional differentiation, thereby determining the cellular identity of the male germline.

Multi-omic spatial effects on high-resolution AI-derived retinal thickness

Retinal thickness is a marker of retinal health and more broadly, is seen as a promising biomarker for many systemic diseases. Retinal thickness measurements are procured from optical coherence tomography (OCT) as part of routine clinical eyecare. We processed the UK Biobank OCT images using a convolutional neural network to produce fine-scale retinal thickness measurements across > 29,000 points in the macula, the part of the retina responsible for human central vision. The macula is disproportionately affected by high disease burden retinal disorders such as age-related macular degeneration and diabetic retinopathy, which both involve metabolic dysregulation. Analysis of common genomic variants, metabolomic, blood and immune biomarkers, disease PheCodes and genetic scores across a fine-scale macular thickness grid, reveals multiple novel genetic loci including four on the X chromosome; retinal thinning associated with many systemic disorders including multiple sclerosis; and multiple associations to correlated metabolites that cluster spatially in the retina. We highlight parafoveal thickness to be particularly susceptible to systemic insults. These results demonstrate the gains in discovery power and resolution achievable with AI-leveraged analysis. Results are accessible using a bespoke web interface that gives full control to pursue findings.

Loss of paired immunoglobin-like type 2 receptor B gene associated with age-related macular degeneration impairs photoreceptor function in mouse retina

Genome-wide association studies have uncovered mostly non-coding variants at over 60 genetic loci linked to susceptibility for age-related macular degeneration (AMD). To ascertain the causal gene at the PILRB/PILRA locus, we used a CRISPR strategy to produce germline deletions in the mouse paired immunoglobin-like type 2 receptor (Pilr) genes that encode highly related activating (PILRB) and inhibitory (PILRA) receptors. We show that a combined loss of Pilrb1 and Pilrb2, but not Pilra, leads to an early but relatively stationary defect as the electroretinography (ERG) amplitudes of Pilrb1/2-/- mice exhibit a marked reduction as early as postnatal day 15 and do not show additional significant decrease at 3 and 12-months. No alterations are evident in Müller glia, microglia, bipolar, amacrine and horizontal cells based on immunohistochemistry using cell-type specific markers. PILRB immunostaining is specifically detected at the proximal part of photoreceptor outer segment. Reduced expression of select calcium-regulated phototransduction and synapse-associated proteins, including GCAP1 and 2, PDE6b, AIPL1, PSD95, and CTBP1 indicates dysregulation of calcium homeostasis as a possible mechanism of retinal phenotype in Pilrb1/2-/- mice. Our studies suggest a novel function of PILRB in retinal photoreceptors and an association of PILRB, but not PILRA, with AMD pathogenesis.

EnsembleSE: identification of super-enhancers based on ensemble learning

Super-enhancers (SEs) are typically located in the regulatory regions of genes, driving high-level gene expression. Identifying SEs is crucial for a deeper understanding of gene regulatory networks, disease mechanisms, and the development and physiological processes of organisms, thus exerting a profound impact on research and applications in the life sciences field. Traditional experimental methods for identifying SEs are costly and time-consuming. Existing methods for predicting SEs based solely on sequence data use deep learning for feature representation and have achieved good results. However, they overlook biological features related to physicochemical properties, leading to low interpretability. Additionally, the complex model structure often requires extensive labeled data for training, which limits their further application in biological data. In this paper, we integrate the strengths of different models and proposes an ensemble model based on an integration strategy to enhance the model's generalization ability. It designs a multi-angle feature representation method that combines local structure and global information to extract high-dimensional abstract relationships and key low-dimensional biological features from sequences. This enhances the effectiveness and interpretability of the model's input features, providing technical support for discovering cell-specific and species-specific patterns of SEs. We evaluated the performance on both mouse and human datasets using five metrics, including area under the receiver operating characteristic curve accuracy, and others. Compared to the latest models, EnsembleSE achieved an average improvement of 4.5% in F1 score and an average improvement of 8.05% in recall, demonstrating the robustness and adaptability of the model on a unified test set. Source codes are available at https://github.com/2103374200/EnsembleSE-main.

Maf-family bZIP transcription factor NRL interacts with RNA-binding proteins and R-loops in retinal photoreceptors

RNA-binding proteins (RBPs) perform diverse functions including the regulation of chromatin dynamics and the coupling of transcription with RNA processing. However, our understanding of their actions in mammalian neurons remains limited. Using affinity purification, yeast-two-hybrid and proximity ligation assays, we identified interactions of multiple RBPs with NRL, a Maf-family bZIP transcription factor critical for retinal rod photoreceptor development and function. In addition to splicing, many NRL-interacting RBPs are associated with R-loops, which form during transcription and increase during photoreceptor maturation. Focusing on DHX9 RNA helicase, we demonstrate that its expression is modulated by NRL and that the NRL-DHX9 interaction is positively influenced by R-loops. ssDRIP-Seq analysis reveals both stranded and unstranded R-loops at distinct genomic elements, characterized by active and inactive epigenetic signatures and enriched at neuronal genes. NRL binds to both types of R-loops, suggesting an epigenetically independent function. Our findings suggest additional functions of NRL during transcription and highlight complex interactions among transcription factors, RBPs, and R-loops in regulating photoreceptor gene expression in the mammalian retina.

Structural framework to address variant-gene relationship in primary open-angle glaucoma

Primary open-angle glaucoma (POAG) is a complex, multifactorial disease leading to progressive optic neuropathy and irreversible vision loss. Genome-Wide Association Studies (GWAS) have significantly advanced our understanding of the genetic loci associated with POAG. Expanding on these findings, Exome-Wide Association Studies (ExWAS) refine the genetic landscape by identifying rare coding variants with potential functional relevance. Post-GWAS in silico analyses, including fine-mapping, gene-based association testing, and pathway analysis, offer insights into target genes and biological mechanisms underlying POAG. This review aims to provide a comprehensive roadmap for the post-GWAS characterization of POAG genes. We integrate current knowledge from GWAS, ExWAS, and post-GWAS analyses, highlighting key genetic variants and pathways implicated in POAG. Recent advancements in genomics, such as ATAC-seq, CUT&RUN, and Hi-C, are crucial for identifying disease-relevant gene regulatory elements by profiling chromatin accessibility, histone modifications, and three-dimensional chromatin architecture. These approaches help pinpoint regulatory elements that influence gene expression in POAG. Expression Quantitative Trait Loci (eQTL) analysis and Transcriptome-Wide Association Studies (TWAS) elucidate the impact of these elements on gene expression and disease risk, while functional validations like enhancer reporter assays confirm their relevance. The integration of high-resolution genomics with functional assays and the characterization of genes in vivo using animal models provides a robust framework for unraveling the complex genetic architecture of POAG. This roadmap is essential for advancing our understanding and identification of genes and regulatory networks involved in POAG pathogenesis.

Multiscale 3D genome rewiring during PTF1A-mediated somatic cell reprogramming into neural stem cells

The genome is intricately folded into chromatin compartments, topologically associating domains (TADs) and loops unique to each cell type. How this higher-order genome organization regulates cell fate transition remains elusive. Here we show how a single non-neural progenitor transcription factor, PTF1A, reorchestrates the 3D genome during fibroblast transdifferentiation into neural stem cells (NSCs). Multiomics analyses integrating Hi-C data, PTF1A and CTCF DNA-binding profiles, H3K27ac modification, and gene expression, demonstrate that PTF1A binds to subTAD boundaries subsequently associated with elevated CTCF binding and enhanced boundary insulation, and reorganizes chromatin loops, leading to gene expression changes that drive transdifferentiation into NSCs. Moreover, PTF1A activates enhancers and super-enhancers near low-insulation boundaries and modulates H3K27ac deposition, promoting cell fate transitions. Together, our data implicate an involvement of 3D genome in transcriptional and cell fate alterations, and highlight an essential role for PTF1A in gene expression control and multiscale 3D genome remodeling during cell reprogramming.

Overcoming artificial structures in resolution-enhanced Hi-C data by signal decomposition and multi-scale attention

Computational enhancement is an important strategy for inferring high-resolution features from genome-wide chromosome conformation capture (Hi-C) data, which typically have limited resolution. Deep learning has been highly successful in this task but we show that it creates prevalent artificial structures in the enhanced data due to the need to divide the large contact matrix into small patches. In addition, previous deep learning methods largely focus on local patterns, which cannot fully capture the complexity of Hi-C data. Here we propose Smooth, High-resolution, and Accurate Reconstruction of Patterns (SHARP) for enhancing Hi-C data. It uses the novel approach of decomposing the data into three types of signals, due to one-dimensional proximity, contiguous domains, and other fine structures, and applies deep learning only to the third type of signals, such that enhancement of the first two is unaffected by the patches. For the deep learning part, SHARP uses both local and global attention mechanisms to capture multi-scale contextual information. We compare SHARP with state-of-the-art methods extensively, including application to data from new samples and another species, and show that SHARP has superior performance in terms of resolution enhancement accuracy, avoiding creation of artificial structures, identifying significant interactions, and enrichment in chromatin states.

Molecular basis of CRX/DNA recognition and stoichiometry at the Ret4 response element

Two retinal transcription factors, cone-rod homeobox (CRX) and neural retina leucine zipper (NRL), cooperate functionally and physically to control photoreceptor development and homeostasis. Mutations in CRX and NRL cause severe retinal diseases. Despite the roles of NRL and CRX, insight into their functions at the molecular level is lacking. Here, we have solved the crystal structure of the CRX homeodomain in complex with its cognate response element (Ret4) from the rhodopsin proximal promoter region. The structure reveals an unexpected 2:1 stoichiometry of CRX/Ret4 and unique orientation of CRX molecules on DNA, and it explains the mechanisms of pathogenic mutations in CRX. Mutations R41Q and E42K disrupt the CRX protein-protein contacts based on the structure and reduce the CRX/Ret4 binding stoichiometry, suggesting a novel disease mechanism. Furthermore, we show that NRL alters the stoichiometry and increases affinity of CRX binding at the rhodopsin promoter, which may enhance transcription of rod-specific genes and suppress transcription of cone-specific genes.

Genetic Risk of Reticular Pseudodrusen in Age-Related Macular Degeneration: HTRA1 /lncRNA BX842242.1 dominates, with no evidence for Complement Cascade involvement

Age-related macular degeneration (AMD) is a multifactorial retinal disease with a large genetic risk contribution. Reticular pseudodrusen (RPD) is a sub-phenotype of AMD with a high risk of progression to late vision threatening AMD. In a genome-wide association study of 2,165 AMD+/RPD+ and 4,181 AMD+/RPD-compared to 7,660 control participants, both chromosomes 1 ( CFH ) and 10 ( ARMS2/HTRA1 ) major AMD risk loci were reidentified. However association was only detected for the chromosome 10 locus when comparing AMD+/RPD+ to AMD+/RPD-cases. The chromosome 1 locus was notably absent. The chromosome 10 RPD risk region contains a long non-coding RNA (ENSG00000285955/BX842242.1) which colocalizes with genetic markers of retinal thickness. BX842242.1 has a strong retinal eQTL signal, pinpointing the parafoveal photoreceptor outer segment layer. Whole genome sequencing of phenotypically extreme RPD cases identified even stronger enrichment for the chromosome 10 risk genotype.

Aging of the eye: Lessons from cataracts and age-related macular degeneration

Aging is the greatest risk factor for chronic human diseases, including many eye diseases. Geroscience aims to understand the effects of the aging process on these diseases, including the genetic, molecular, and cellular mechanisms that underlie the increased risk of disease over the lifetime. Understanding of the aging eye increases general knowledge of the cellular physiology impacted by aging processes at various biological extremes. Two major diseases, age-related cataract and age-related macular degeneration (AMD) are caused by dysfunction of the lens and retina, respectively. Lens transparency and light refraction are mediated by lens fiber cells lacking nuclei and other organelles, which provides a unique opportunity to study a single aging hallmark, i.e., loss of proteostasis, within an environment of limited metabolism. In AMD, local dysfunction of the photoreceptors/retinal pigmented epithelium/Bruch's membrane/choriocapillaris complex in the macula leads to the loss of photoreceptors and eventually loss of central vision, and is driven by nearly all the hallmarks of aging and shares features with Alzheimer's disease, Parkinson's disease, cardiovascular disease, and diabetes. The aging eye can function as a model for studying basic mechanisms of aging and, vice versa, well-defined hallmarks of aging can be used as tools to understand age-related eye disease.

Enhancers on the edge - how the nuclear envelope controls gene regulatory elements

Precise temporal and sequential control of gene expression during development and in response to environmental stimuli requires tight regulation of the physical contact between gene regulatory elements and promoters. Current models describing how the genome folds in 3D space to establish these interactions often ignore the role of the most stable structural nuclear feature - the nuclear envelope. While contributions of 3D folding within/between topologically associated domains (TADs) have been extensively described, mechanical contributions from the nuclear envelope can impact enhancer-promoter interactions both directly and indirectly through influencing intra/inter-TAD interactions. Importantly, these nuclear envelope contributions clearly link this mechanism to development and, when defective, to human disease. Here, we discuss evidence for nuclear envelope regulation of tissue-specific enhancer-promoter pairings, potential mechanisms for this regulation, exciting recent findings that other regulatory elements such as microRNAs and long noncoding RNAs are under nuclear envelope regulation, the possible involvement of condensates, and how disruption of this regulation can lead to disease.

Epigenome-metabolism nexus in the retina: implications for aging and disease

Intimate links between epigenome modifications and metabolites allude to a crucial role of cellular metabolism in transcriptional regulation. Retina, being a highly metabolic tissue, adapts by integrating inputs from genetic, epigenetic, and extracellular signals. Precise global epigenomic signatures guide development and homeostasis of the intricate retinal structure and function. Epigenomic and metabolic realignment are hallmarks of aging and highlight a link of the epigenome-metabolism nexus with aging-associated multifactorial traits affecting the retina, including age-related macular degeneration and glaucoma. Here, we focus on emerging principles of epigenomic and metabolic control of retinal gene regulation, with emphasis on their contribution to human disease. In addition, we discuss potential mitigation strategies involving lifestyle changes that target the epigenome-metabolome relationship for maintaining retinal function.

The chromatin accessibility and transcriptomic landscape of the aging mice cochlea and the identification of potential functional super-enhancers in age-related hearing loss

BACKGROUND

Presbycusis, also referred to as age-related hearing loss (ARHL), is a condition that results from the cumulative effects of aging on an individual's auditory capabilities. Given the limited understanding of epigenetic mechanisms in ARHL, our research focuses on alterations in chromatin-accessible regions.

METHODS

We employed assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) in conjunction with unique identifier (UID) mRNA-seq between young and aging cochleae, and conducted integrated analysis as well as motif/TF-gene prediction. Additionally, the essential role of super-enhancers (SEs) in the development of ARHL was identified by comparative analysis to previous research. Meanwhile, an ARHL mouse model and an aging mimic hair cell (HC) model were established with a comprehensive identification of senescence phenotypes to access the role of SEs in ARHL progression.

RESULTS

The control cochlear tissue exhibited greater chromatin accessibility than cochlear tissue affected by ARHL. Furthermore, the levels of histone 3 lysine 27 acetylation were significantly depressed in both aging cochlea and aging mimic HEI-OC1 cells, highlighting the essential role of SEs in the development of ARHL. The potential senescence-associated super-enhancers (SASEs) of ARHL were identified, most of which exhibited decreased chromatin accessibility. The majority of genes related to the SASEs showed obvious decreases in mRNA expression level in aging HCs and was noticeably altered following treatment with JQ1 (a commonly used SE inhibitor).

CONCLUSION

The chromatin accessibility in control cochlear tissue was higher than that in cochlear tissue affected by ARHL. Potential SEs involved in ARHL were identified, which might provide a basis for future therapeutics targeting SASEs related to ARHL.

Evolutionary conservation of VSX2 super-enhancer modules in retinal development

Super-enhancers (SEs) are expansive regions of genomic DNA that regulate the expression of genes involved in cell identity and cell fate. We recently identified developmental stage- and cell type-specific modules within the murine Vsx2 SE. Here, we show that the human VSX2 SE modules have similar developmental stage- and cell type-specific activity in reporter gene assays. By inserting the human sequence of one VSX2 SE module into a mouse with microphthalmia, eye size was rescued. To understand the function of these SE modules during human retinal development, we deleted individual modules in human embryonic stem cells and generated retinal organoids. Deleting one module results in small organoids, recapitulating the small-eyed phenotype of mice with microphthalmia, while deletion of the other module led to disruptions in bipolar neuron development. This prototypical SE serves as a model for understanding developmental stage- and cell type-specific effects of neurogenic transcription factors with complex expression patterns. Moreover, by elucidating the gene regulatory mechanisms, we can begin to examine how dysregulation of these mechanisms contributes to phenotypic diversity and disease.

Variants in UBAP1L lead to autosomal recessive rod-cone and cone-rod dystrophy

PURPOSE

Progressive inherited retinal degenerations (IRDs) affecting rods and cones are clinically and genetically heterogeneous and can lead to blindness with limited therapeutic options. The major gene defects have been identified in subjects of European and Asian descent with only few reports of North African descent.

METHODS

Genome, targeted next-generation, and Sanger sequencing was applied to cohort of ∼4000 IRDs cases. Expression analyses were performed including Chip-seq database analyses, on human-derived retinal organoids (ROs), retinal pigment epithelium cells, and zebrafish. Variants' pathogenicity was accessed using 3D-modeling and/or ROs.

RESULTS

Here, we identified a novel gene defect with three distinct pathogenic variants in UBAP1L in 4 independent autosomal recessive IRD cases from Tunisia. UBAP1L is expressed in the retinal pigment epithelium and retina, specifically in rods and cones, in line with the phenotype. It encodes Ubiquitin-associated protein 1-like, containing a solenoid of overlapping ubiquitin-associated domain, predicted to interact with ubiquitin. In silico and in vitro studies, including 3D-modeling and ROs revealed that the solenoid of overlapping ubiquitin-associated domain is truncated and thus ubiquitin binding most likely abolished secondary to all variants identified herein.

CONCLUSION

Biallelic UBAP1L variants are a novel cause of IRDs, most likely enriched in the North African population.

Comparative 3D genome analysis between neural retina and retinal pigment epithelium reveals differential cis-regulatory interactions at retinal disease loci

BACKGROUND

Vision depends on the interplay between photoreceptor cells of the neural retina and the underlying retinal pigment epithelium (RPE). Most genes involved in inherited retinal diseases display specific spatiotemporal expression within these interconnected retinal components through the local recruitment of cis-regulatory elements (CREs) in 3D nuclear space.

RESULTS

To understand the role of differential chromatin architecture in establishing tissue-specific expression at inherited retinal disease loci, we mapped genome-wide chromatin interactions using in situ Hi-C and H3K4me3 HiChIP on neural retina and RPE/choroid from human adult donor eyes. We observed chromatin looping between active promoters and 32,425 and 8060 candidate CREs in the neural retina and RPE/choroid, respectively. A comparative 3D genome analysis between these two retinal tissues revealed that 56% of 290 known inherited retinal disease genes were marked by differential chromatin interactions. One of these was ABCA4, which is implicated in the most common autosomal recessive inherited retinal disease. We zoomed in on retina- and RPE-specific cis-regulatory interactions at the ABCA4 locus using high-resolution UMI-4C. Integration with bulk and single-cell epigenomic datasets and in vivo enhancer assays in zebrafish revealed tissue-specific CREs interacting with ABCA4.

CONCLUSIONS

Through comparative 3D genome mapping, based on genome-wide, promoter-centric, and locus-specific assays of human neural retina and RPE, we have shown that gene regulation at key inherited retinal disease loci is likely mediated by tissue-specific chromatin interactions. These findings do not only provide insight into tissue-specific regulatory landscapes at retinal disease loci, but also delineate the search space for non-coding genomic variation underlying unsolved inherited retinal diseases.

Analysis of long-range chromatin contacts, compartments and looping between mouse embryonic stem cells, lens epithelium and lens fibers

BACKGROUND

Nuclear organization of interphase chromosomes involves individual chromosome territories, "open" and "closed" chromatin compartments, topologically associated domains (TADs) and chromatin loops. The DNA- and RNA-binding transcription factor CTCF together with the cohesin complex serve as major organizers of chromatin architecture. Cellular differentiation is driven by temporally and spatially coordinated gene expression that requires chromatin changes of individual loci of various complexities. Lens differentiation represents an advantageous system to probe transcriptional mechanisms underlying tissue-specific gene expression including high transcriptional outputs of individual crystallin genes until the mature lens fiber cells degrade their nuclei.

RESULTS

Chromatin organization between mouse embryonic stem (ES) cells, newborn (P0.5) lens epithelium and fiber cells were analyzed using Hi-C. Localization of CTCF in both lens chromatins was determined by ChIP-seq and compared with ES cells. Quantitative analyses show major differences between number and size of TADs and chromatin loop size between these three cell types. In depth analyses show similarities between lens samples exemplified by overlaps between compartments A and B. Lens epithelium-specific CTCF peaks are found in mostly methylated genomic regions while lens fiber-specific and shared peaks occur mostly within unmethylated DNA regions. Major differences in TADs and loops are illustrated at the ~ 500 kb Pax6 locus, encoding the critical lens regulatory transcription factor and within a larger ~ 15 Mb WAGR locus, containing Pax6 and other loci linked to human congenital diseases. Lens and ES cell Hi-C data (TADs and loops) together with ATAC-seq, CTCF, H3K27ac, H3K27me3 and ENCODE cis-regulatory sites are shown in detail for the Pax6, Sox1 and Hif1a loci, multiple crystallin genes and other important loci required for lens morphogenesis. The majority of crystallin loci are marked by unexpectedly high CTCF-binding across their transcribed regions.

CONCLUSIONS

Our study has generated the first data on 3-dimensional (3D) nuclear organization in lens epithelium and lens fibers and directly compared these data with ES cells. These findings generate novel insights into lens-specific transcriptional gene control, open new research avenues to study transcriptional condensates in lens fiber cells, and enable studies of non-coding genetic variants linked to cataract and other lens and ocular abnormalities.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation provides a crucial epigenetic mark linking genetic variations to environmental influence. We have analyzed array-based DNA methylation profiles of 160 human retinas with co-measured RNA-seq and >8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 methylation quantitative trait loci and 12,505 expression quantitative trait loci) and 13,747 DNA methylation loci affecting gene expression, with over one-third specific to the retina. Methylation and expression quantitative trait loci show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration. Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of macular degeneration pathology by genotype-environment interaction in retina.

Multi-omics analysis in human retina uncovers ultraconserved cis-regulatory elements at rare eye disease loci

Cross-species genome comparisons have revealed a substantial number of ultraconserved non-coding elements (UCNEs). Several of these elements have proved to be essential tissue- and cell type-specific cis-regulators of developmental gene expression. Here, we characterize a set of UCNEs as candidate CREs (cCREs) during retinal development and evaluate the contribution of their genomic variation to rare eye diseases, for which pathogenic non-coding variants are emerging. Integration of bulk and single-cell retinal multi-omics data reveals 594 genes under potential cis-regulatory control of UCNEs, of which 45 are implicated in rare eye disease. Mining of candidate cis-regulatory UCNEs in WGS data derived from the rare eye disease cohort of Genomics England reveals 178 ultrarare variants within 84 UCNEs associated with 29 disease genes. Overall, we provide a comprehensive annotation of ultraconserved non-coding regions acting as cCREs during retinal development which can be targets of non-coding variation underlying rare eye diseases.

Whole genome sequencing of 4,787 individuals identifies gene-based rare variants in age-related macular degeneration

Genome-wide association studies have contributed extensively to the discovery of disease-associated common variants. However, the genetic contribution to complex traits is still largely difficult to interpret. We report a genome-wide association study of 2394 cases and 2393 controls for age-related macular degeneration (AMD) via whole-genome sequencing, with 46.9 million genetic variants. Our study reveals significant single-variant association signals at four loci and independent gene-based signals in CFH, C2, C3, and NRTN. Using data from the Exome Aggregation Consortium (ExAC) for a gene-based test, we demonstrate an enrichment of predicted rare loss-of-function variants in CFH, CFI, and an as-yet unreported gene in AMD, ORMDL2. Our method of using a large variant list without individual-level genotypes as an external reference provides a flexible and convenient approach to leverage the publicly available variant datasets to augment the search for rare variant associations, which can explain additional disease risk in AMD.

Integrating genetic regulation and single-cell expression with GWAS prioritizes causal genes and cell types for glaucoma

Primary open-angle glaucoma (POAG), characterized by retinal ganglion cell death, is a leading cause of irreversible blindness worldwide. However, its molecular and cellular causes are not well understood. Elevated intraocular pressure (IOP) is a major risk factor, but many patients have normal IOP. Colocalization and Mendelian randomization analysis of >240 POAG and IOP genome-wide association study (GWAS) loci and overlapping expression and splicing quantitative trait loci (e/sQTLs) in 49 GTEx tissues and retina prioritizes causal genes for 60% of loci. These genes are enriched in pathways implicated in extracellular matrix organization, cell adhesion, and vascular development. Analysis of single-nucleus RNA-seq of glaucoma-relevant eye tissues reveals that the POAG and IOP colocalizing genes and genome-wide associations are enriched in specific cell types in the aqueous outflow pathways, retina, optic nerve head, peripapillary sclera, and choroid. This study nominates IOP-dependent and independent regulatory mechanisms, genes, and cell types that may contribute to POAG pathogenesis.

Stage-specific dynamic reorganization of genome topology shapes transcriptional neighborhoods in developing human retinal organoids

We have generated a high-resolution Hi-C map of developing human retinal organoids to elucidate spatiotemporal dynamics of genomic architecture and its relationship with gene expression patterns. We demonstrate progressive stage-specific alterations in DNA topology and correlate these changes with transcription of cell-type-restricted gene markers during retinal differentiation. Temporal Hi-C reveals a shift toward A compartment for protein-coding genes and B compartment for non-coding RNAs, displaying high and low expression, respectively. Notably, retina-enriched genes are clustered near lost boundaries of topologically associated domains (TADs), and higher-order assemblages (i.e., TAD cliques) localize in active chromatin regions with binding sites for eye-field transcription factors. These genes gain chromatin contacts at their transcription start site as organoid differentiation proceeds. Our study provides a global view of chromatin architecture dynamics associated with diversification of cell types during retinal development and serves as a foundational resource for in-depth functional investigations of retinal developmental traits.

CTCF-mediated 3D chromatin predetermines the gene expression program in the male germline

Spermatogenesis is a unidirectional differentiation process that generates haploid sperm, but how the gene expression program that directs this process is established is largely unknown. Here we determine the high-resolution 3D chromatin architecture of male germ cells during spermatogenesis and show that CTCF-mediated 3D chromatin predetermines the gene expression program required for spermatogenesis. In undifferentiated spermatogonia, CTCF-mediated chromatin contacts on autosomes pre-establish meiosis-specific super-enhancers (SE). These meiotic SE recruit the master transcription factor A-MYB in meiotic spermatocytes, which strengthens their 3D contacts and instructs a burst of meiotic gene expression. We also find that at the mitosis-to-meiosis transition, the germline-specific Polycomb protein SCML2 resolves chromatin loops that are specific to mitotic spermatogonia. Moreover, SCML2 and A-MYB establish the unique 3D chromatin organization of sex chromosomes during meiotic sex chromosome inactivation. We propose that CTCF-mediated 3D chromatin organization enforces epigenetic priming that directs unidirectional differentiation, thereby determining the cellular identity of the male germline.

Exploring the contribution of ARMS2 and HTRA1 genetic risk factors in age-related macular degeneration

Age-related macular degeneration (AMD) is the leading cause of severe irreversible central vision loss in individuals over 65 years old. Genome-wide association studies (GWASs) have shown that the region at chromosome 10q26, where the age-related maculopathy susceptibility (ARMS2/LOC387715) and HtrA serine peptidase 1 (HTRA1) genes are located, represents one of the strongest associated loci for AMD. However, the underlying biological mechanism of this genetic association has remained elusive. In this article, we extensively review the literature by us and others regarding the ARMS2/HTRA1 risk alleles and their functional significance. We also review the literature regarding the presumed function of the ARMS2 protein and the molecular processes of the HTRA1 protein in AMD pathogenesis in vitro and in vivo, including those of transgenic mice overexpressing HtrA1/HTRA1 which developed Bruch's membrane (BM) damage, choroidal neovascularization (CNV), and polypoidal choroidal vasculopathy (PCV), similar to human AMD patients. The elucidation of the molecular mechanisms of the ARMS2 and HTRA1 susceptibility loci has begun to untangle the complex biological pathways underlying AMD pathophysiology, pointing to new testable paradigms for treatment.

Identification of Evolutionarily Conserved VSX2 Enhancers in Retinal Development

Super-enhancers (SEs) are expansive regions of genomic DNA that regulate the expression of genes involved in cell identity and cell fate. Recently, we found that distinct modules within a murine SE regulate gene expression of master regulatory transcription factor Vsx2 in a developmental stage- and cell-type specific manner. Vsx2 is expressed in retinal progenitor cells as well as differentiated bipolar neurons and Müller glia. Mutations in VSX2 in humans and mice lead to microphthalmia due to a defect in retinal progenitor cell proliferation. Deletion of a single module within the Vsx2 SE leads to microphthalmia. Deletion of a separate module within the SE leads to a complete loss of bipolar neurons, yet the remainder of the retina develops normally. Furthermore, the Vsx2 SE is evolutionarily conserved in vertebrates, suggesting that these modules are important for retinal development across species. In the present study, we examine the ability of these modules to drive retinal development between species. By inserting the human build of one Vsx2 SE module into a mouse with microphthalmia, eye size was rescued. To understand the implications of these SE modules in a model of human development, we generated human retinal organoids. Deleting one module results in small organoids, recapitulating the small-eyed phenotype of mice with microphthalmia, while deletion of the other module leads to a complete loss of ON cone bipolar neurons. This prototypical SE serves as a model for uncoupling developmental stage- and cell-type specific effects of neurogenic transcription factors with complex expression patterns. Moreover, by elucidating the gene regulatory mechanisms, we can begin to examine how dysregulation of these mechanisms contributes to phenotypic diversity and disease.

10q26 - The enigma in age-related macular degeneration

Despite comprehensive research efforts over the last decades, the pathomechanisms of age-related macular degeneration (AMD) remain far from being understood. Large-scale genome wide association studies (GWAS) were able to provide a defined set of genetic aberrations which contribute to disease risk, with the strongest contributors mapping to distinct regions on chromosome 1 and 10. While the chromosome 1 locus comprises factors of the complement system with well-known functions, the role of the 10q26-locus in AMD-pathophysiology remains enigmatic. 10q26 harbors a cluster of three functional genes, namely PLEKHA1, ARMS2 and HTRA1, with most of the AMD-associated genetic variants mapping to the latter two genes. High linkage disequilibrium between ARMS2 and HTRA1 has kept association studies from reliably defining the risk-causing gene for long and only very recently the genetic risk region has been narrowed to ARMS2, suggesting that this is the true AMD gene at this locus. However, genetic associations alone do not suffice to prove causality and one or more of the 14 SNPs on this haplotype may be involved in long-range control of gene expression, leaving HTRA1 and PLEKHA1 still suspects in the pathogenic pathway. Both, ARMS2 and HTRA1 have been linked to extracellular matrix homeostasis, yet their exact molecular function as well as their role in AMD pathogenesis remains to be uncovered. The transcriptional regulation of the 10q26 locus adds an additional level of complexity, given, that gene-regulatory as well as epigenetic alterations may influence expression levels from 10q26 in diseased individuals. Here, we provide a comprehensive overview on the 10q26 locus and its three gene products on various levels of biological complexity and discuss current and future research strategies to shed light on one of the remaining enigmatic spots in the AMD landscape.

The CUT&RUN suspect list of problematic regions of the genome

BACKGROUND

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is an increasingly popular technique to map genome-wide binding profiles of histone modifications, transcription factors, and co-factors. The ENCODE project and others have compiled blacklists for ChIP-seq which have been widely adopted: these lists contain regions of high and unstructured signal, regardless of cell type or protein target, indicating that these are false positives. While CUT&RUN obtains similar results to ChIP-seq, its biochemistry and subsequent data analyses are different. We found that this results in a CUT&RUN-specific set of undesired high-signal regions.

RESULTS

We compile suspect lists based on CUT&RUN data for the human and mouse genomes, identifying regions consistently called as peaks in negative controls. Using published CUT&RUN data from our and other labs, we show that the CUT&RUN suspect regions can persist even when peak calling is performed with SEACR or MACS2 against a negative control and after ENCODE blacklist removal. Moreover, we experimentally validate the CUT&RUN suspect lists by performing reiterative negative control experiments in which no specific protein is targeted, showing that they capture more than 80% of the peaks identified.

CONCLUSIONS

We propose that removing these problematic regions can substantially improve peak calling in CUT&RUN experiments, resulting in more reliable datasets.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation (DNAm) provides a crucial epigenetic mark linking genetic variations to environmental influence. We analyzed array-based DNAm profiles of 160 human retinas with co-measured RNA-seq and > 8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 mQTLs and 12,505 eQTLs) and 13,747 eQTMs (DNAm loci affecting gene expression), with over one-third specific to the retina. mQTLs and eQTMs show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration (AMD). Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of AMD pathology by genotype-environment interaction in retina.

Super-enhancer landscape rewiring in cancer: The epigenetic control at distal sites

Super-enhancers evolve as elements at the top of the hierarchical control of gene expression. They are important end-gatherers of signaling pathways that control stemness, differentiation or adaptive responses. Many epigenetic regulations focus on these regions, and not surprisingly, during the process of tumorigenesis, various alterations can account for their dysfunction. Super-enhancers are emerging as key drivers of the aberrant gene expression landscape that sustain the aggressiveness of cancer cells. In this review, we will describe and discuss about the structure of super-enhancers, their epigenetic regulation, and the major changes affecting their functionality in cancer.

Superenhancers as master gene regulators and novel therapeutic targets in brain tumors

Transcriptional deregulation, a cancer cell hallmark, is driven by epigenetic abnormalities in the majority of brain tumors, including adult glioblastoma and pediatric brain tumors. Epigenetic abnormalities can activate epigenetic regulatory elements to regulate the expression of oncogenes. Superenhancers (SEs), identified as novel epigenetic regulatory elements, are clusters of enhancers with cell-type specificity that can drive the aberrant transcription of oncogenes and promote tumor initiation and progression. As gene regulators, SEs are involved in tumorigenesis in a variety of tumors, including brain tumors. SEs are susceptible to inhibition by their key components, such as bromodomain protein 4 and cyclin-dependent kinase 7, providing new opportunities for antitumor therapy. In this review, we summarized the characteristics and identification, unique organizational structures, and activation mechanisms of SEs in tumors, as well as the clinical applications related to SEs in tumor therapy and prognostication. Based on a review of the literature, we discussed the relationship between SEs and different brain tumors and potential therapeutic targets, focusing on glioblastoma.

Mini-review: Gene regulatory network benefits from three-dimensional chromatin conformation and structural biology

Gene regulatory networks are now at the forefront of precision biology, which can help researchers better understand how genes and regulatory elements interact to control cellular gene expression, offering a more promising molecular mechanism in biological research. Interactions between the genes and regulatory elements involve different promoters, enhancers, transcription factors, silencers, insulators, and long-range regulatory elements, which occur at a ∼10 µm nucleus in a spatiotemporal manner. In this way, three-dimensional chromatin conformation and structural biology are critical for interpreting the biological effects and the gene regulatory networks. In the review, we have briefly summarized the latest processes in three-dimensional chromatin conformation, microscopic imaging, and bioinformatics, and we have presented the outlook and future directions for these three aspects.