Functional characteristics of novel pancreatic Pax6 regulatory elements

A bird's eye view on the use of whole exome sequencing in rare congenital ophthalmic diseases

Phenotypic and genotypic heterogeneity in congenital ocular diseases, especially in anterior segment dysgenesis (ASD), have created challenges for proper diagnosis and classification of diseases. Over the last decade, genomic research has indeed boosted our understanding in the molecular basis of ASD and genes associated with both autosomal dominant and recessive patterns of inheritance have been described with a wide range of expressivity. Here we describe the molecular characterization of a cohort of 162 patients displaying isolated or syndromic congenital ocular dysgenesis. Samples were analyzed with diverse techniques, such as direct sequencing, multiplex ligation-dependent probe amplification, and whole exome sequencing (WES), over 20 years. Our data reiterate the notion that PAX6 alterations are primarily associated with ASD, mostly aniridia, since the majority of the cohort (66.7%) has a pathogenic or likely pathogenic variant in the PAX6 locus. Unexpectedly, a high fraction of positive samples (20.3%) displayed deletions involving the 11p13 locus, either partially/totally involving PAX6 coding region or abolishing its critical regulatory region, underlying its significance. Most importantly, the use of WES has allowed us to both assess variants in known ASD genes (i.e., CYP1B1, ITPR1, MAB21L1, PXDN, and PITX2) and to identify rarer phenotypes (i.e., MIDAS, oculogastrointestinal-neurodevelopmental syndrome and Jacobsen syndrome). Our data clearly suggest that WES allows expanding the analytical portfolio of ocular dysgenesis, both isolated and syndromic, and that is pivotal for the differential diagnosis of those conditions in which there may be phenotypic overlaps and in general in ASD.

Analysis of single nuclear chromatin accessibility reveals unique myeloid populations in human pancreatic ductal adenocarcinoma

BACKGROUND

A better understanding of the pancreatic ductal adenocarcinoma (PDAC) immune microenvironment is critical to developing new treatments and improving outcomes. Myeloid cells are of particular importance for PDAC progression; however, the presence of heterogenous subsets with different ontogeny and impact, along with some fluidity between them, (infiltrating monocytes vs. tissue-resident macrophages; M1 vs. M2) makes characterisation of myeloid populations challenging. Recent advances in single cell sequencing technology provide tools for characterisation of immune cell infiltrates, and open chromatin provides source and function data for myeloid cells to assist in more comprehensive characterisation. Thus, we explore single nuclear assay for transposase accessible chromatin (ATAC) sequencing (snATAC-Seq), a method to analyse open gene promoters and transcription factor binding, as an important means for discerning the myeloid composition in human PDAC tumours.

METHODS

Frozen pancreatic tissues (benign or PDAC) were prepared for snATAC-Seq using 10× Chromium technology. Signac was used for preliminary analysis, clustering and differentially accessible chromatin region identification. The genes annotated in promoter regions were used for Gene Ontology (GO) enrichment and cell type annotation. Gene signatures were used for survival analysis with The Cancer Genome Atlas (TCGA)-pancreatic adenocarcinoma (PAAD) dataset.

RESULTS

Myeloid cell transcription factor activities were higher in tumour than benign pancreatic samples, enabling us to further stratify tumour myeloid populations. Subcluster analysis revealed eight distinct myeloid populations. GO enrichment demonstrated unique functions for myeloid populations, including interleukin-1b signalling (recruited monocytes) and intracellular protein transport (dendritic cells). The identified gene signature for dendritic cells influenced survival (hazard ratio = .63, p = .03) in the TCGA-PAAD dataset, which was unique to PDAC.

CONCLUSIONS

These data suggest snATAC-Seq as a method for analysis of frozen human pancreatic tissues to distinguish myeloid populations. An improved understanding of myeloid cell heterogeneity and function is important for developing new treatment targets in PDAC.

The Density of Regulatory Information Is a Major Determinant of Evolutionary Constraint on Noncoding DNA in Drosophila

Evolutionary analyses have estimated that ∼60% of nucleotides in intergenic regions of the Drosophila melanogaster genome are functionally relevant, suggesting that regulatory information may be encoded more densely in intergenic regions than has been revealed by most functional dissections of regulatory DNA. Here, we approached this issue through a functional dissection of the regulatory region of the gene shavenbaby (svb). Most of the ∼90 kb of this large regulatory region is highly conserved in the genus Drosophila, though characterized enhancers occupy a small fraction of this region. By analyzing the regulation of svb in different contexts of Drosophila development, we found that the regulatory information that drives svb expression in the abdominal pupal epidermis is organized in a different way than the elements that drive svb expression in the embryonic epidermis. While in the embryonic epidermis svb is activated by compact enhancers separated by large inactive DNA regions, svb expression in the pupal epidermis is driven by regulatory information distributed over broader regions of svb cis-regulatory DNA. In the same vein, we observed that other developmental genes also display a dense distribution of putative regulatory elements in their regulatory regions. Furthermore, we found that a large percentage of conserved noncoding DNA of the Drosophila genome is contained within regions of open chromatin. These results suggest that part of the evolutionary constraint on noncoding DNA of Drosophila is explained by the density of regulatory information, which may be greater than previously appreciated.

AT7867 promotes pancreatic progenitor differentiation of human iPSCs

Generation of pure pancreatic progenitor (PP) cells is critical for clinical translation of stem cell-derived islets. Herein, we performed PP differentiation with and without AKT/P70 inhibitor AT7867 and characterized the resulting cells at protein and transcript level in vitro and in vivo upon transplantation into diabetic mice. AT7867 treatment increased the percentage of PDX1+NKX6.1+ (-AT7867: 50.9% [IQR 48.9%-53.8%]; +AT7867: 90.8% [IQR 88.9%-93.7%]; p = 0.0021) and PDX1+GP2+ PP cells (-AT7867: 39.22% [IQR 36.7%-44.1%]; +AT7867: 90.0% [IQR 88.2%-93.6%]; p = 0.0021). Transcriptionally, AT7867 treatment significantly upregulated PDX1 (p = 0.0001), NKX6.1 (p = 0.0005), and GP2 (p = 0.002) expression compared with controls, while off-target markers PODXL (p < 0.0001) and TBX2 (p < 0.0001) were significantly downregulated. Transplantation of AT7867-treated PPs resulted in faster hyperglycemia reversal in diabetic mice compared with controls (time and group: p < 0.0001). Overall, our data show that AT7867 enhances PP cell differentiation leading to accelerated diabetes reversal.

Future directions in managing aniridia-associated keratopathy

Congenital aniridia is a panocular disorder that is typically characterized by iris hypoplasia and aniridia-associated keratopathy (AAK). AAK results in the progressive loss of corneal transparency and thereby loss of vision. Currently, there is no approved therapy to delay or prevent its progression, and clinical management is challenging because of phenotypic variability and high risk of complications after interventions; however, new insights into the molecular pathogenesis of AAK may help improve its management. Here, we review the current understanding about the pathogenesis and management of AAK. We highlight the biological mechanisms involved in AAK development with the aim to develop future treatment options, including surgical, pharmacological, cell therapies, and gene therapies.

Transcription modulates chromatin dynamics and locus configuration sampling

In living cells, the 3D structure of gene loci is dynamic, but this is not revealed by 3C and FISH experiments in fixed samples, leaving a notable gap in our understanding. To overcome these limitations, we applied the highly predictive heteromorphic polymer (HiP-HoP) model to determine chromatin fiber mobility at the Pax6 locus in three mouse cell lines with different transcription states. While transcriptional activity minimally affects movement of 40-kbp regions, we observed that motion of smaller 1-kbp regions depends strongly on local disruption to chromatin fiber structure marked by H3K27 acetylation. This also substantially influenced locus configuration dynamics by modulating protein-mediated promoter-enhancer loops. Importantly, these simulations indicate that chromatin dynamics are sufficiently fast to sample all possible locus conformations within minutes, generating wide dynamic variability within single cells. This combination of simulation and experimental validation provides insight into how transcriptional activity influences chromatin structure and gene dynamics.

Characterization of stem-cell-derived islets during differentiation and after implantation

Recapitulation of embryonic pancreatic development has enabled development of methods for in vitro islet cell differentiation using human pluripotent stem cells (hPSCs), which have the potential to cure diabetes. Advanced methods for optimal generation of stem-cell-derived islets (SC-islets) has enabled successful diabetes reversal in rodents and shown promising early clinical trial outcomes. The main impediment for use of SC-islets is concern about safety because of off-target growth resulting from contaminated residual cells. In this review, we summarize the different endocrine and non-endocrine cell populations that have been described to emerge throughout β cell differentiation and after transplantation. We discuss the most recent approaches to enrich endocrine populations and remove off-target cells. Finally, we discuss the critical quality control and release criteria testing that we anticipate will be required prior to transplantation to ensure product safety.

Single-cell RNA sequencing of the Strongylocentrotus purpuratus larva reveals the blueprint of major cell types and nervous system of a non-chordate deuterostome

Identifying the molecular fingerprint of organismal cell types is key for understanding their function and evolution. Here, we use single-cell RNA sequencing (scRNA-seq) to survey the cell types of the sea urchin early pluteus larva, representing an important developmental transition from non-feeding to feeding larva. We identify 21 distinct cell clusters, representing cells of the digestive, skeletal, immune, and nervous systems. Further subclustering of these reveal a highly detailed portrait of cell diversity across the larva, including the identification of neuronal cell types. We then validate important gene regulatory networks driving sea urchin development and reveal new domains of activity within the larval body. Focusing on neurons that co-express Pdx-1 and Brn1/2/4, we identify an unprecedented number of genes shared by this population of neurons in sea urchin and vertebrate endocrine pancreatic cells. Using differential expression results from Pdx-1 knockdown experiments, we show that Pdx1 is necessary for the acquisition of the neuronal identity of these cells. We hypothesize that a network similar to the one orchestrated by Pdx1 in the sea urchin neurons was active in an ancestral cell type and then inherited by neuronal and pancreatic developmental lineages in sea urchins and vertebrates.

The Spectrum of PAX6 Mutations and Genotype-Phenotype Correlations in the Eye

The transcription factor PAX6 is essential in ocular development in vertebrates, being considered the master regulator of the eye. During eye development, it is essential for the correct patterning and formation of the multi-layered optic cup and it is involved in the developing lens and corneal epithelium. In adulthood, it is mostly expressed in cornea, iris, and lens. PAX6 is a dosage-sensitive gene and it is highly regulated by several elements located upstream, downstream, and within the gene. There are more than 500 different mutations described to affect PAX6 and its regulatory regions, the majority of which lead to PAX6 haploinsufficiency, causing several ocular and systemic abnormalities. Aniridia is an autosomal dominant disorder that is marked by the complete or partial absence of the iris, foveal hypoplasia, and nystagmus, and is caused by heterozygous PAX6 mutations. Other ocular abnormalities have also been associated with PAX6 changes, and genotype-phenotype correlations are emerging. This review will cover recent advancements in PAX6 regulation, particularly the role of several enhancers that are known to regulate PAX6 during eye development and disease. We will also present an updated overview of the mutation spectrum, where an increasing number of mutations in the non-coding regions have been reported. Novel genotype-phenotype correlations will also be discussed.

Tales from topographic oceans: topologically associated domains and cancer

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

BAT Hi-C maps global chromatin interactions in an efficient and economical way

Chromosome Conformation Capture (3C)-based technologies, such as Hi-C, have represented a significant breakthrough in investigating the structure and function of higher-order genome architecture. However, the mapping of global chromatin interactions remains challenging across many biological conditions due to high background noise and financial constraints, especially for small laboratories. Here, we describe the Bridge linker-Alul-Tn5 Hi-C (BAT Hi-C) method, which is a simple and efficient method for delineating chromatin conformational features of mouse embryonic stem (mES) cells and uncover DNA loops. This protocol combines Alul fragmentation and biotinylated linker-mediated proximity ligation to obtain kilobase (kb) resolution with a marked increase in the amount of unique read pairs. The protocol also includes chromatin isolation to reduce background noise and Tn5 tagmentation to cut down on preparation time. Importantly, with only one-third sequencing depth, our method revealed the same spectrum of chromatin contacts as in situ Hi-C. BAT Hi-C is an economical (i.e., approximately $40 for library preparation) and straightforward (total hands-on time of 3 days) tool that is ideal for the in-depth analysis of long-range chromatin looping events in a genome-wide fashion.

Heterogeneous Loop Model to Infer 3D Chromosome Structures from Hi-C

Adapting a well-established formalism in polymer physics, we develop a minimalist approach to infer three-dimensional folding of chromatin from Hi-C data. The three-dimensional chromosome structures generated from our heterogeneous loop model (HLM) are used to visualize chromosome organizations that can substantiate the measurements from fluorescence in situ hybridization, chromatin interaction analysis by paired-end tag sequencing, and RNA-seq signals. We demonstrate the utility of the HLM with several case studies. Specifically, the HLM-generated chromosome structures, which reproduce the spatial distribution of topologically associated domains from fluorescence in situ hybridization measurement, show the phase segregation between two types of topologically associated domains explicitly. We discuss the origin of cell-type-dependent gene-expression level by modeling the chromatin globules of α-globin and SOX2 gene loci for two different cell lines. We also use the HLM to discuss how the chromatin folding and gene-expression level of Pax6 loci, associated with mouse neural development, are modulated by interactions with two enhancers. Finally, HLM-generated structures of chromosome 19 of mouse embryonic stem cells, based on single-cell Hi-C data collected over each cell-cycle phase, visualize changes in chromosome conformation along the cell-cycle. Given a contact frequency map between chromatic loci supplied from Hi-C, HLM is a computationally efficient and versatile modeling tool to generate chromosome structures that can complement interpreting other experimental data.

Polymer Simulations of Heteromorphic Chromatin Predict the 3D Folding of Complex Genomic Loci

Chromatin folded into 3D macromolecular structures is often analyzed by chromosome conformation capture (3C) and fluorescence in situ hybridization (FISH) techniques, but these frequently provide contradictory results. Chromatin can be modeled as a simple polymer composed of a connected chain of units. By embedding data for epigenetic marks (H3K27ac), chromatin accessibility (assay for transposase-accessible chromatin using sequencing [ATAC-seq]), and structural anchors (CCCTC-binding factor [CTCF]), we developed a highly predictive heteromorphic polymer (HiP-HoP) model, where the chromatin fiber varied along its length; combined with diffusing protein bridges and loop extrusion, this model predicted the 3D organization of genomic loci at a population and single-cell level. The model was validated at several gene loci, including the complex Pax6 gene, and was able to determine locus conformations across cell types with varying levels of transcriptional activity and explain different mechanisms of enhancer use. Minimal a priori knowledge of epigenetic marks is sufficient to recapitulate complex genomic loci in 3D and enable predictions of chromatin folding paths.

•

HiP-HoP: highly predictive heteromorphic polymer model to analyze chromatin structure

•

Polymer simulations use widely available epigenetic and protein binding data as input

•

Validate HiP-HoP model at complex loci using 3D FISH and Capture-C

•

Simulations uncover striking conformational variability in chromatin fiber folding