RASER-FISH: non-denaturing fluorescence in situ hybridization for preservation of three-dimensional interphase chromatin structure.

DNA fluorescence in situ hybridization (FISH) has been a central technique in advancing our understanding of how chromatin is organized within the nucleus. With the increasing resolution offered by super-resolution microscopy, the optimal maintenance of chromatin structure within the nucleus is essential for accuracy in measurements and interpretation of data. However, standard 3D-FISH requires potentially destructive heat denaturation in the presence of chaotropic agents such as formamide to allow access to the DNA strands for labeled FISH probes. To avoid the need to heat-denature, we developed Resolution After Single-strand Exonuclease Resection (RASER)-FISH, which uses exonuclease digestion to generate single-stranded target DNA for efficient probe binding over a 2 d process. Furthermore, RASER-FISH is easily combined with immunostaining of nuclear proteins or the detection of RNAs. Here, we provide detailed procedures for RASER-FISH in mammalian cultured cells to detect single loci, chromatin tracks and topologically associating domains with conventional and super-resolution 3D structured illumination microscopy. Moreover, we provide a validation and characterization of our method, demonstrating excellent preservation of chromatin structure and nuclear integrity, together with improved hybridization efficiency, compared with classic 3D-FISH protocols.

High-resolution targeted 3C interrogation of cis-regulatory element organization at genome-wide scale.

Chromosome conformation capture (3C) provides an adaptable tool for studying diverse biological questions. Current 3C methods generally provide either low-resolution interaction profiles across the entire genome, or high-resolution interaction profiles at limited numbers of loci. Due to technical limitations, generation of reproducible high-resolution interaction profiles has not been achieved at genome-wide scale. Here, to overcome this barrier, we systematically test each step of 3C and report two improvements over current methods. We show that up to 30% of reporter events generated using the popular in situ 3C method arise from ligations between two individual nuclei, but this noise can be almost entirely eliminated by isolating intact nuclei after ligation. Using Nuclear-Titrated Capture-C, we generate reproducible high-resolution genome-wide 3C interaction profiles by targeting 8055 gene promoters in erythroid cells. By pairing high-resolution 3C interaction calls with nascent gene expression we interrogate the role of promoter hubs and super-enhancers in gene regulation.

Genetic and functional insights into CDA-I prevalence and pathogenesis.

BACKGROUND: Congenital dyserythropoietic anaemia type I (CDA-I) is a hereditary anaemia caused by biallelic mutations in the widely expressed genes CDAN1 and C15orf41. Little is understood about either protein and it is unclear in which cellular pathways they participate. METHODS: Genetic analysis of a cohort of patients with CDA-I identifies novel pathogenic variants in both known causative genes. We analyse the mutation distribution and the predicted structural positioning of amino acids affected in Codanin-1, the protein encoded by CDAN1. Using western blotting, immunoprecipitation and immunofluorescence, we determine the effect of particular mutations on both proteins and interrogate protein interaction, stability and subcellular localisation. RESULTS: We identify six novel CDAN1 mutations and one novel mutation in C15orf41 and uncover evidence of further genetic heterogeneity in CDA-I. Additionally, population genetics suggests that CDA-I is more common than currently predicted. Mutations are enriched in six clusters in Codanin-1 and tend to affect buried residues. Many missense and in-frame mutations do not destabilise the entire protein. Rather C15orf41 relies on Codanin-1 for stability and both proteins, which are enriched in the nucleolus, interact to form an obligate complex in cells. CONCLUSION: Stability and interaction data suggest that C15orf41 may be the key determinant of CDA-I and offer insight into the mechanism underlying this disease. Both proteins share a common pathway likely to be present in a wide variety of cell types; however, nucleolar enrichment may provide a clue as to the erythroid specific nature of CDA-I. The surprisingly high predicted incidence of CDA-I suggests that better ascertainment would lead to improved patient care.

Recapitulation of erythropoiesis in congenital dyserythropoietic anaemia type I (CDA-I) identifies defects in differentiation and nucleolar abnormalities.

The investigation of inherited disorders of erythropoiesis has elucidated many of the principles underlying the production of normal red blood cells and how this is perturbed in human disease. Congenital Dyserythropoietic Anaemia type 1 (CDA-I) is a rare form of anaemia caused by mutations in two genes of unknown function: CDAN1 and CDIN1 (previously called C15orf41), whilst in some cases, the underlying genetic abnormality is completely unknown. Consequently, the pathways affected in CDA-I remain to be discovered. To enable detailed analysis of this rare disorder we have validated a culture system which recapitulates all of the cardinal haematological features of CDA-I, including the formation of the pathognomonic 'spongy' heterochromatin seen by electron microscopy. Using a variety of cell and molecular biological approaches we discovered that erythroid cells in this condition show a delay during terminal erythroid differentiation, associated with increased proliferation and widespread changes in chromatin accessibility. We also show that the proteins encoded by CDAN1 and CDIN1 are enriched in nucleoli which are structurally and functionally abnormal in CDA-I. Together these findings provide important pointers to the pathways affected in CDA-I which for the first time can now be pursued in the tractable culture system utilised here.

Dynamics of the 4D genome during in vivo lineage specification and differentiation.

Mammalian gene expression patterns are controlled by regulatory elements, which interact within topologically associating domains (TADs). The relationship between activation of regulatory elements, formation of structural chromatin interactions and gene expression during development is unclear. Here, we present Tiled-C, a low-input chromosome conformation capture (3C) technique. We use this approach to study chromatin architecture at high spatial and temporal resolution through in vivo mouse erythroid differentiation. Integrated analysis of chromatin accessibility and single-cell expression data shows that regulatory elements gradually become accessible within pre-existing TADs during early differentiation. This is followed by structural re-organization within the TAD and formation of specific contacts between enhancers and promoters. Our high-resolution data show that these enhancer-promoter interactions are not established prior to gene expression, but formed gradually during differentiation, concomitant with progressive upregulation of gene activity. Together, these results provide new insight into the close, interdependent relationship between chromatin architecture and gene regulation during development.

A Dynamic Folded Hairpin Conformation Is Associated with α-Globin Activation in Erythroid Cells.

We investigate the three-dimensional (3D) conformations of the α-globin locus at the single-allele level in murine embryonic stem cells (ESCs) and erythroid cells, combining polymer physics models and high-resolution Capture-C data. Model predictions are validated against independent fluorescence in situ hybridization (FISH) data measuring pairwise distances, and Tri-C data identifying three-way contacts. The architecture is rearranged during the transition from ESCs to erythroid cells, associated with the activation of the globin genes. We find that in ESCs, the spatial organization conforms to a highly intermingled 3D structure involving non-specific contacts, whereas in erythroid cells the α-globin genes and their enhancers form a self-contained domain, arranged in a folded hairpin conformation, separated from intermingling flanking regions by a thermodynamic mechanism of micro-phase separation. The flanking regions are rich in convergent CTCF sites, which only marginally participate in the erythroid-specific gene-enhancer contacts, suggesting that beyond the interaction of CTCF sites, multiple molecular mechanisms cooperate to form an interacting domain.

ATR-16 syndrome: mechanisms linking monosomy to phenotype.

BACKGROUND: Deletions removing 100s-1000s kb of DNA, and variable numbers of poorly characterised genes, are often found in patients with a wide range of developmental abnormalities. In such cases, understanding the contribution of the deletion to an individual's clinical phenotype is challenging. METHODS: Here, as an example of this common phenomenon, we analysed 41 patients with simple deletions of ~177 to ~2000 kb affecting one allele of the well-characterised, gene dense, distal region of chromosome 16 (16p13.3), referred to as ATR-16 syndrome. We characterised deletion extents and screened for genetic background effects, telomere position effect and compensatory upregulation of hemizygous genes. RESULTS: We find the risk of developmental and neurological abnormalities arises from much smaller distal chromosome 16 deletions (~400 kb) than previously reported. Beyond this, the severity of ATR-16 syndrome increases with deletion size, but there is no evidence that critical regions determine the developmental abnormalities associated with this disorder. Surprisingly, we find no evidence of telomere position effect or compensatory upregulation of hemizygous genes; however, genetic background effects substantially modify phenotypic abnormalities. CONCLUSIONS: Using ATR-16 as a general model of disorders caused by CNVs, we show the degree to which individuals with contiguous gene syndromes are affected is not simply related to the number of genes deleted but depends on their genetic background. We also show there is no critical region defining the degree of phenotypic abnormalities in ATR-16 syndrome and this has important implications for genetic counselling.

Cohesin Disrupts Polycomb-Dependent Chromosome Interactions in Embryonic Stem Cells.

How chromosome organization is related to genome function remains poorly understood. Cohesin, loop extrusion, and CCCTC-binding factor (CTCF) have been proposed to create topologically associating domains (TADs) to regulate gene expression. Here, we examine chromosome conformation in embryonic stem cells lacking cohesin and find, as in other cell types, that cohesin is required to create TADs and regulate A/B compartmentalization. However, in the absence of cohesin, we identify a series of long-range chromosomal interactions that persist. These correspond to regions of the genome occupied by the polycomb repressive system and are dependent on PRC1. Importantly, we discover that cohesin counteracts these polycomb-dependent interactions, but not interactions between super-enhancers. This disruptive activity is independent of CTCF and insulation and appears to modulate gene repression by the polycomb system. Therefore, we discover that cohesin disrupts polycomb-dependent chromosome interactions to modulate gene expression in embryonic stem cells.

Single-allele chromatin interactions identify regulatory hubs in dynamic compartmentalized domains.

The promoters of mammalian genes are commonly regulated by multiple distal enhancers, which physically interact within discrete chromatin domains. How such domains form and how the regulatory elements within them interact in single cells is not understood. To address this we developed Tri-C, a new chromosome conformation capture (3C) approach, to characterize concurrent chromatin interactions at individual alleles. Analysis by Tri-C identifies heterogeneous patterns of single-allele interactions between CTCF boundary elements, indicating that the formation of chromatin domains likely results from a dynamic process. Within these domains, we observe specific higher-order structures that involve simultaneous interactions between multiple enhancers and promoters. Such regulatory hubs provide a structural basis for understanding how multiple cis-regulatory elements act together to establish robust regulation of gene expression.

A tissue-specific self-interacting chromatin domain forms independently of enhancer-promoter interactions.

Self-interacting chromatin domains encompass genes and their cis-regulatory elements; however, the three-dimensional form a domain takes, whether this relies on enhancer-promoter interactions, and the processes necessary to mediate the formation and maintenance of such domains, remain unclear. To examine these questions, here we use a combination of high-resolution chromosome conformation capture, a non-denaturing form of fluorescence in situ hybridisation and super-resolution imaging to study a 70 kb domain encompassing the mouse α-globin regulatory locus. We show that this region forms an erythroid-specific, decompacted, self-interacting domain, delimited by frequently apposed CTCF/cohesin binding sites early in terminal erythroid differentiation, and does not require transcriptional elongation for maintenance of the domain structure. Formation of this domain does not rely on interactions between the α-globin genes and their major enhancers, suggesting a transcription-independent mechanism for establishment of the domain. However, absence of the major enhancers does alter internal domain interactions. Formation of a loop domain therefore appears to be a mechanistic process that occurs irrespective of the specific interactions within.

Robust CRISPR/Cas9 Genome Editing of the HUDEP-2 Erythroid Precursor Line Using Plasmids and Single-Stranded Oligonucleotide Donors.

The study of cellular processes and gene regulation in terminal erythroid development has been greatly facilitated by the generation of an immortalised erythroid cell line derived from Human Umbilical Derived Erythroid Precursors, termed HUDEP-2 cells. The ability to efficiently genome edit HUDEP-2 cells and make clonal lines hugely expands their utility as the insertion of clinically relevant mutations allows study of potentially every genetic disease affecting red blood cell development. Additionally, insertion of sequences encoding short protein tags such as Strep, FLAG and Myc permits study of protein behaviour in the normal and disease state. This approach is useful to augment the analysis of patient cells as large cell numbers are obtainable with the additional benefit that the need for specific antibodies may be circumvented. This approach is likely to lead to insights into disease mechanisms and provide reagents to allow drug discovery. HUDEP-2 cells provide a favourable alternative to the existing immortalised erythroleukemia lines as their karyotype is much less abnormal. These cells also provide sufficient material for a broad range of analyses as it is possible to generate in vitro-differentiated erythroblasts in numbers 4-7 fold higher than starting cell numbers within 9-12 days of culture. Here we describe an efficient, robust and reproducible plasmid-based methodology to introduce short (<20 bp) DNA sequences into the genome of HUDEP-2 cells using the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 Cas9 system combined with single-stranded oligodeoxynucleotide (ssODN) donors. This protocol produces genetically modified lines in ~30 days and could also be used to generate knock-out and knock-in mutations.

Predicting the three-dimensional folding of cis-regulatory regions in mammalian genomes using bioinformatic data and polymer models.

The three-dimensional (3D) organization of chromosomes can be probed using methods like Capture-C. However, it is unclear how such population-level data relate to the organization within a single cell, and the mechanisms leading to the observed interactions are still largely obscure. We present a polymer modeling scheme based on the assumption that chromosome architecture is maintained by protein bridges, which form chromatin loops. To test the model, we perform FISH experiments and compare with Capture-C data. Starting merely from the locations of protein binding sites, our model accurately predicts the experimentally observed chromatin interactions, revealing a population of 3D conformations.

Factors influencing success of clinical genome sequencing across a broad spectrum of disorders.

To assess factors influencing the success of whole-genome sequencing for mainstream clinical diagnosis, we sequenced 217 individuals from 156 independent cases or families across a broad spectrum of disorders in whom previous screening had identified no pathogenic variants. We quantified the number of candidate variants identified using different strategies for variant calling, filtering, annotation and prioritization. We found that jointly calling variants across samples, filtering against both local and external databases, deploying multiple annotation tools and using familial transmission above biological plausibility contributed to accuracy. Overall, we identified disease-causing variants in 21% of cases, with the proportion increasing to 34% (23/68) for mendelian disorders and 57% (8/14) in family trios. We also discovered 32 potentially clinically actionable variants in 18 genes unrelated to the referral disorder, although only 4 were ultimately considered reportable. Our results demonstrate the value of genome sequencing for routine clinical diagnosis but also highlight many outstanding challenges.

Expanded GAA repeats impair FXN gene expression and reposition the FXN locus to the nuclear lamina in single cells.

Abnormally expanded DNA repeats are associated with several neurodegenerative diseases. In Friedreich's ataxia (FRDA), expanded GAA repeats in intron 1 of the frataxin gene (FXN) reduce FXN mRNA levels in averaged cell samples through a poorly understood mechanism. By visualizing FXN expression and nuclear localization in single cells, we show that GAA-expanded repeats decrease the number of FXN mRNA molecules, slow transcription, and increase FXN localization at the nuclear lamina (NL). Restoring histone acetylation reverses NL positioning. Expanded GAA-FXN loci in FRDA patient cells show increased NL localization with increased silencing of alleles and reduced transcription from alleles positioned peripherally. We also demonstrate inefficiencies in transcription initiation and elongation from the expanded GAA-FXN locus at single-cell resolution. We suggest that repressive epigenetic modifications at the expanded GAA-FXN locus may lead to NL relocation, where further repression may occur.

De novo and rare inherited mutations implicate the transcriptional coregulator TCF20/SPBP in autism spectrum disorder.

BACKGROUND: Autism spectrum disorders (ASDs) are common and have a strong genetic basis, yet the cause of ∼70-80% ASDs remains unknown. By clinical cytogenetic testing, we identified a family in which two brothers had ASD, mild intellectual disability and a chromosome 22 pericentric inversion, not detected in either parent, indicating de novo mutation with parental germinal mosaicism. We hypothesised that the rearrangement was causative of their ASD and localised the chromosome 22 breakpoints. METHODS: The rearrangement was characterised using fluorescence in situ hybridisation, Southern blotting, inverse PCR and dideoxy-sequencing. Open reading frames and intron/exon boundaries of the two physically disrupted genes identified, TCF20 and TNRC6B, were sequenced in 342 families (260 multiplex and 82 simplex) ascertained by the International Molecular Genetic Study of Autism Consortium (IMGSAC). RESULTS: IMGSAC family screening identified a de novo missense mutation of TCF20 in a single case and significant association of a different missense mutation of TCF20 with ASD in three further families. Through exome sequencing in another project, we independently identified a de novo frameshifting mutation of TCF20 in a woman with ASD and moderate intellectual disability. We did not identify a significant association of TNRC6B mutations with ASD. CONCLUSIONS: TCF20 encodes a transcriptional coregulator (also termed SPBP) that is structurally and functionally related to RAI1, the critical dosage-sensitive protein implicated in the behavioural phenotypes of the Smith-Magenis and Potocki-Lupski 17p11.2 deletion/duplication syndromes, in which ASD is frequently diagnosed. This study provides the first evidence that mutations in TCF20 are also associated with ASD.

Homozygous mutations in a predicted endonuclease are a novel cause of congenital dyserythropoietic anemia type I.

The congenital dyserythropoietic anemias are a heterogeneous group of rare disorders primarily affecting erythropoiesis with characteristic morphological abnormalities and a block in erythroid maturation. Mutations in the CDAN1 gene, which encodes Codanin-1, underlie the majority of congenital dyserythropoietic anemia type I cases. However, no likely pathogenic CDAN1 mutation has been detected in approximately 20% of cases, suggesting the presence of at least one other locus. We used whole genome sequencing and segregation analysis to identify a homozygous T to A transversion (c.533T>A), predicted to lead to a p.L178Q missense substitution in C15ORF41, a gene of unknown function, in a consanguineous pedigree of Middle-Eastern origin. Sequencing C15ORF41 in other CDAN1 mutation-negative congenital dyserythropoietic anemia type I pedigrees identified a homozygous transition (c.281A>G), predicted to lead to a p.Y94C substitution, in two further pedigrees of SouthEast Asian origin. The haplotype surrounding the c.281A>G change suggests a founder effect for this mutation in Pakistan. Detailed sequence similarity searches indicate that C15ORF41 encodes a novel restriction endonuclease that is a member of the Holliday junction resolvase family of proteins.

A biophysical model for transcription factories.

SUMMARY: Transcription factories are nuclear domains where gene transcription takes place although the molecular basis for their formation and maintenance are unknown. In this study, we explored how the properties of chromatin as a polymer may contribute to the structure of transcription factories. We found that transcriptional active chromatin contains modifications like histone H4 acetylated at Lysine 16 (H4K16ac). Single fibre analysis showed that this modification spans the entire body of the gene. Furthermore, H4K16ac genes cluster in regions up to 500 Kb alternating active and inactive chromatin. The introduction of H4K16ac in chromatin induces stiffness in the chromatin fibre. The result of this change in flexibility is that chromatin could behave like a multi-block copolymer with repetitions of stiff-flexible (active-inactive chromatin) components. Copolymers with such structure self-organize through spontaneous phase separation into microdomains. Consistent with such model H4K16ac chromatin form foci that associates with nascent transcripts. We propose that transcription factories are the result of the spontaneous concentration of H4K16ac chromatin that are in proximity, mainly in cis.

Intragenic enhancers act as alternative promoters.

A substantial amount of organismal complexity is thought to be encoded by enhancers which specify the location, timing, and levels of gene expression. In mammals there are more enhancers than promoters which are distributed both between and within genes. Here we show that activated, intragenic enhancers frequently act as alternative tissue-specific promoters producing a class of abundant, spliced, multiexonic poly(A)(+) RNAs (meRNAs) which reflect the host gene's structure. meRNAs make a substantial and unanticipated contribution to the complexity of the transcriptome, appearing as alternative isoforms of the host gene. The low protein-coding potential of meRNAs suggests that many meRNAs may be byproducts of enhancer activation or underlie as-yet-unidentified RNA-encoded functions. Distinguishing between meRNAs and mRNAs will transform our interpretation of dynamic changes in transcription both at the level of individual genes and of the genome as a whole.

Codanin-1 mutations in congenital dyserythropoietic anemia type 1 affect HP1{alpha} localization in erythroblasts.

Congenital dyserythropoietic anemia type 1 (CDA-1), a rare inborn anemia characterized by abnormal chromatin ultrastructure in erythroblasts, is caused by abnormalities in codanin-1, a highly conserved protein of unknown function. We have produced 3 monoclonal antibodies to codanin-1 that demonstrate its distribution in both nucleus and cytoplasm by immunofluorescence and allow quantitative measurements of patient and normal material by Western blot. A detailed analysis of chromatin structure in CDA-1 erythroblasts shows no abnormalities in overall histone composition, and the genome-wide epigenetic landscape of several histone modifications is maintained. However, immunofluorescence analysis of intermediate erythroblasts from patients with CDA-1 reveals abnormal accumulation of HP1α in the Golgi apparatus. A link between mutant codanin-1 and the aberrant localization of HP1α is supported by the finding that codanin-1 can be coimmunoprecipitated by anti-HP1α antibodies. Furthermore, we show colocalization of codanin-1 with Sec23B, the protein defective in CDA-2 suggesting that the CDAs might be linked at the molecular level. Mice containing a gene-trapped Cdan1 locus demonstrate its widespread expression during development. Cdan1(gt/gt) homozygotes die in utero before the onset of primitive erythropoiesis, suggesting that Cdan1 has other critical roles during embryogenesis.

Global gene expression analysis of human erythroid progenitors.

Understanding the pattern of gene expression during erythropoiesis is crucial for a synthesis of erythroid developmental biology. Here, we isolated 4 distinct populations at successive erythropoietin-dependent stages of erythropoiesis, including the terminal, pyknotic stage. The transcriptome was determined using Affymetrix arrays. First, we demonstrated the importance of using defined cell populations to identify lineage and temporally specific patterns of gene expression. Cells sorted by surface expression profile not only express significantly fewer genes than unsorted cells but also demonstrate significantly greater differences in the expression levels of particular genes between stages than unsorted cells. Second, using standard software, we identified more than 1000 transcripts not previously observed to be differentially expressed during erythroid maturation, 13 of which are highly significantly terminally regulated, including RFXAP and SMARCA4. Third, using matched filtering, we identified 12 transcripts not previously reported to be continuously up-regulated in maturing human primary erythroblasts. Finally, using transcription factor binding site analysis, we identified potential transcription factors that may regulate gene expression during terminal erythropoiesis. Our stringent lists of differentially regulated and continuously expressed transcripts containing many genes with undiscovered functions in erythroblasts are a resource for future functional studies of erythropoiesis. Our Human Erythroid Maturation database is available at https://cellline.molbiol.ox.ac.uk/eryth/index.html. [corrected].