MED27, SLC6A7, and MPPE1 Variants in a Complex Neurodevelopmental Disorder with Severe Dystonia.

BACKGROUND: Despite advances in next generation sequencing technologies, the identification of variants of uncertain significance (VUS) can often hinder definitive diagnosis in patients with complex neurodevelopmental disorders. OBJECTIVE: The objective of this study was to identify and characterize the underlying cause of disease in a family with two children with severe developmental delay associated with generalized dystonia and episodic status dystonicus, chorea, epilepsy, and cataracts. METHODS: Candidate genes identified by autozygosity mapping and whole-exome sequencing were characterized using cellular and vertebrate model systems. RESULTS: Homozygous variants were found in three candidate genes: MED27, SLC6A7, and MPPE1. Although the patients had features of MED27-related disorder, the SLC6A7 and MPPE1 variants were functionally investigated. SLC6A7 variant in vitro overexpression caused decreased proline transport as a result of reduced cell-surface expression, and zebrafish knockdown of slc6a7 exhibited developmental delay and fragile motor neuron morphology that could not be rescued by L-proline transporter-G396S RNA. Lastly, patient fibroblasts displayed reduced cell-surface expression of glycophosphatidylinositol-anchored proteins linked to MPPE1 dysfunction. CONCLUSIONS: We report a family harboring a homozygous MED27 variant with additional loss-of-function SLC6A7 and MPPE1 gene variants, which potentially contribute to a blended phenotype caused by multilocus pathogenic variants. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Hierarchical chromatin organization detected by TADpole.

The rapid development of Chromosome Conformation Capture (3C-based techniques), as well as imaging together with bioinformatics analyses, has been fundamental for unveiling that chromosomes are organized into the so-called topologically associating domains or TADs. While TADs appear as nested patterns in the 3C-based interaction matrices, the vast majority of available TAD callers are based on the hypothesis that TADs are individual and unrelated chromatin structures. Here we introduce TADpole, a computational tool designed to identify and analyze the entire hierarchy of TADs in intra-chromosomal interaction matrices. TADpole combines principal component analysis and constrained hierarchical clustering to provide a set of significant hierarchical chromatin levels in a genomic region of interest. TADpole is robust to data resolution, normalization strategy and sequencing depth. Domain borders defined by TADpole are enriched in main architectural proteins (CTCF and cohesin complex subunits) and in the histone mark H3K4me3, while their domain bodies, depending on their activation-state, are enriched in either H3K36me3 or H3K27me3, highlighting that TADpole is able to distinguish functional TAD units. Additionally, we demonstrate that TADpole's hierarchical annotation, together with the new DiffT score, allows for detecting significant topological differences on Capture Hi-C maps between wild-type and genetically engineered mouse.

Walking along chromosomes with super-resolution imaging, contact maps, and integrative modeling.

Chromosome organization is crucial for genome function. Here, we present a method for visualizing chromosomal DNA at super-resolution and then integrating Hi-C data to produce three-dimensional models of chromosome organization. Using the super-resolution microscopy methods of OligoSTORM and OligoDNA-PAINT, we trace 8 megabases of human chromosome 19, visualizing structures ranging in size from a few kilobases to over a megabase. Focusing on chromosomal regions that contribute to compartments, we discover distinct structures that, in spite of considerable variability, can predict whether such regions correspond to active (A-type) or inactive (B-type) compartments. Imaging through the depths of entire nuclei, we capture pairs of homologous regions in diploid cells, obtaining evidence that maternal and paternal homologous regions can be differentially organized. Finally, using restraint-based modeling to integrate imaging and Hi-C data, we implement a method-integrative modeling of genomic regions (IMGR)-to increase the genomic resolution of our traces to 10 kb.

Conformational changes during pore formation by the perforin-related protein pleurotolysin.

Membrane attack complex/perforin-like (MACPF) proteins comprise the largest superfamily of pore-forming proteins, playing crucial roles in immunity and pathogenesis. Soluble monomers assemble into large transmembrane pores via conformational transitions that remain to be structurally and mechanistically characterised. Here we present an 11 Å resolution cryo-electron microscopy (cryo-EM) structure of the two-part, fungal toxin Pleurotolysin (Ply), together with crystal structures of both components (the lipid binding PlyA protein and the pore-forming MACPF component PlyB). These data reveal a 13-fold pore 80 Å in diameter and 100 Å in height, with each subunit comprised of a PlyB molecule atop a membrane bound dimer of PlyA. The resolution of the EM map, together with biophysical and computational experiments, allowed confident assignment of subdomains in a MACPF pore assembly. The major conformational changes in PlyB are a ∼70° opening of the bent and distorted central ß-sheet of the MACPF domain, accompanied by extrusion and refolding of two α-helical regions into transmembrane ß-hairpins (TMH1 and TMH2). We determined the structures of three different disulphide bond-trapped prepore intermediates. Analysis of these data by molecular modelling and flexible fitting allows us to generate a potential trajectory of ß-sheet unbending. The results suggest that MACPF conformational change is triggered through disruption of the interface between a conserved helix-turn-helix motif and the top of TMH2. Following their release we propose that the transmembrane regions assemble into ß-hairpins via top down zippering of backbone hydrogen bonds to form the membrane-inserted ß-barrel. The intermediate structures of the MACPF domain during refolding into the ß-barrel pore establish a structural paradigm for the transition from soluble monomer to pore, which may be conserved across the whole superfamily. The TMH2 region is critical for the release of both TMH clusters, suggesting why this region is targeted by endogenous inhibitors of MACPF function.

Genome-wide tracing to decipher nuclear organization.

Nuclear organization impacts gene expression activity and cell phenotype. Our current understanding is mainly derived from ensemble-level sequencing studies that reflect the 3D genome structure of millions of cells. These approaches have provided invaluable details on the 3D organizations of the genome and their relation to other nuclear landmarks. However, they mostly lack the ability to provide multimodal information simultaneously at the single-cell level. In recent years, cutting-edge imaging technologies have risen to the challenge of simultaneously describing multiple components of the nuclear space at the single-cell level, paving the way for a deeper understanding of the genome structure-function relationship. This review will focus on the development and utilization of such technologies to gain a multi-component view of the nucleus at single-cell resolution, dissecting the complexity and heterogeneity of nuclear organization.

Chromatin alternates between A and B compartments at kilobase scale for subgenic organization.

Nuclear compartments are prominent features of 3D chromatin organization, but sequencing depth limitations have impeded investigation at ultra fine-scale. CTCF loops are generally studied at a finer scale, but the impact of looping on proximal interactions remains enigmatic. Here, we critically examine nuclear compartments and CTCF loop-proximal interactions using a combination of in situ Hi-C at unparalleled depth, algorithm development, and biophysical modeling. Producing a large Hi-C map with 33 billion contacts in conjunction with an algorithm for performing principal component analysis on sparse, super massive matrices (POSSUMM), we resolve compartments to 500 bp. Our results demonstrate that essentially all active promoters and distal enhancers localize in the A compartment, even when flanking sequences do not. Furthermore, we find that the TSS and TTS of paused genes are often segregated into separate compartments. We then identify diffuse interactions that radiate from CTCF loop anchors, which correlate with strong enhancer-promoter interactions and proximal transcription. We also find that these diffuse interactions depend on CTCF's RNA binding domains. In this work, we demonstrate features of fine-scale chromatin organization consistent with a revised model in which compartments are more precise than commonly thought while CTCF loops are more protracted.

3D reconstruction of genomic regions from sparse interaction data.

Chromosome conformation capture (3C) technologies measure the interaction frequency between pairs of chromatin regions within the nucleus in a cell or a population of cells. Some of these 3C technologies retrieve interactions involving non-contiguous sets of loci, resulting in sparse interaction matrices. One of such 3C technologies is Promoter Capture Hi-C (pcHi-C) that is tailored to probe only interactions involving gene promoters. As such, pcHi-C provides sparse interaction matrices that are suitable to characterize short- and long-range enhancer-promoter interactions. Here, we introduce a new method to reconstruct the chromatin structural (3D) organization from sparse 3C-based datasets such as pcHi-C. Our method allows for data normalization, detection of significant interactions and reconstruction of the full 3D organization of the genomic region despite of the data sparseness. Specifically, it builds, with as low as the 2-3% of the data from the matrix, reliable 3D models of similar accuracy of those based on dense interaction matrices. Furthermore, the method is sensitive enough to detect cell-type-specific 3D organizational features such as the formation of different networks of active gene communities.

Three-dimensional genome organization via triplex-forming RNAs.

An increasing number of long noncoding RNAs (lncRNAs) have been proposed to act as nuclear organization factors during interphase. Direct RNA-DNA interactions can be achieved by the formation of triplex helix structures where a single-stranded RNA molecule hybridizes by complementarity into the major groove of double-stranded DNA. However, whether and how these direct RNA-DNA associations influence genome structure in interphase chromosomes remain poorly understood. Here we theorize that RNA organizes the genome in space via a triplex-forming mechanism. To test this theory, we apply a computational modeling approach of chromosomes that combines restraint-based modeling with polymer physics. Our models suggest that colocalization of triplex hotspots targeted by lncRNAs could contribute to large-scale chromosome compartmentalization cooperating, rather than competing, with architectural transcription factors such as CTCF.

Analysis, Modeling, and Visualization of Chromosome Conformation Capture Experiments.

Chromatin Conformation Capture techniques have unveiled several layers of chromosome organization such as the segregation in compartments, the folding in topologically associating domains (TADs), and site-specific looping interactions. The discovery of this genome hierarchical organization emerged from the computational analysis of chromatin capture data. With the increasing availability of such data, automatic pipelines for the robust comparison, grouping, and classification of multiple experiments are needed. Here we present a pipeline based on the TADbit framework that emphasizes reproducibility, automation, quality check, and statistical robustness. This comprehensive modular pipeline covers all the steps from the sequencing products to the visualization of reconstructed 3D models of the chromatin.

Dynamics of genome architecture and chromatin function during human B cell differentiation and neoplastic transformation.

To investigate the three-dimensional (3D) genome architecture across normal B cell differentiation and in neoplastic cells from different subtypes of chronic lymphocytic leukemia and mantle cell lymphoma patients, here we integrate in situ Hi-C and nine additional omics layers. Beyond conventional active (A) and inactive (B) compartments, we uncover a highly-dynamic intermediate compartment enriched in poised and polycomb-repressed chromatin. During B cell development, 28% of the compartments change, mostly involving a widespread chromatin activation from naive to germinal center B cells and a reversal to the naive state upon further maturation into memory B cells. B cell neoplasms are characterized by both entity and subtype-specific alterations in 3D genome organization, including large chromatin blocks spanning key disease-specific genes. This study indicates that 3D genome interactions are extensively modulated during normal B cell differentiation and that the genome of B cell neoplasias acquires a tumor-specific 3D genome architecture.

Transcriptional activation during cell reprogramming correlates with the formation of 3D open chromatin hubs.

Chromosome structure is a crucial regulatory factor for a wide range of nuclear processes. Chromosome conformation capture (3C)-based experiments combined with computational modelling are pivotal for unveiling 3D chromosome structure. Here, we introduce TADdyn, a tool that integrates time-course 3C data, restraint-based modelling, and molecular dynamics to simulate the structural rearrangements of genomic loci in a completely data-driven way. We apply TADdyn on in situ Hi-C time-course experiments studying the reprogramming of murine B cells to pluripotent cells, and characterize the structural rearrangements that take place upon changes in the transcriptional state of 21 genomic loci of diverse expression dynamics. By measuring various structural and dynamical properties, we find that during gene activation, the transcription starting site contacts with open and active regions in 3D chromatin domains. We propose that these 3D hubs of open and active chromatin may constitute a general feature to trigger and maintain gene transcription.

RNA proximity sequencing reveals the spatial organization of the transcriptome in the nucleus.

The global, three-dimensional organization of RNA molecules in the nucleus is difficult to determine using existing methods. Here we introduce Proximity RNA-seq, which identifies colocalization preferences for pairs or groups of nascent and fully transcribed RNAs in the nucleus. Proximity RNA-seq is based on massive-throughput RNA barcoding of subnuclear particles in water-in-oil emulsion droplets, followed by cDNA sequencing. Our results show RNAs of varying tissue-specificity of expression, speed of RNA polymerase elongation and extent of alternative splicing positioned at varying distances from nucleoli. The simultaneous detection of multiple RNAs in proximity to each other distinguishes RNA-dense from sparse compartments. Application of Proximity RNA-seq will facilitate study of the spatial organization of transcripts in the nucleus, including non-coding RNAs, and its functional relevance.

Human pancreatic islet three-dimensional chromatin architecture provides insights into the genetics of type 2 diabetes.

Genetic studies promise to provide insight into the molecular mechanisms underlying type 2 diabetes (T2D). Variants associated with T2D are often located in tissue-specific enhancer clusters or super-enhancers. So far, such domains have been defined through clustering of enhancers in linear genome maps rather than in three-dimensional (3D) space. Furthermore, their target genes are often unknown. We have created promoter capture Hi-C maps in human pancreatic islets. This linked diabetes-associated enhancers to their target genes, often located hundreds of kilobases away. It also revealed >1,300 groups of islet enhancers, super-enhancers and active promoters that form 3D hubs, some of which show coordinated glucose-dependent activity. We demonstrate that genetic variation in hubs impacts insulin secretion heritability, and show that hub annotations can be used for polygenic scores that predict T2D risk driven by islet regulatory variants. Human islet 3D chromatin architecture, therefore, provides a framework for interpretation of T2D genome-wide association study (GWAS) signals.

TEMPy: a Python library for assessment of three-dimensional electron microscopy density fits.

Three-dimensional electron microscopy is currently one of the most promising techniques used to study macromolecular assemblies. Rigid and flexible fitting of atomic models into density maps is often essential to gain further insights into the assemblies they represent. Currently, tools that facilitate the assessment of fitted atomic models and maps are needed. TEMPy (template and electron microscopy comparison using Python) is a toolkit designed for this purpose. The library includes a set of methods to assess density fits in intermediate-to-low resolution maps, both globally and locally. It also provides procedures for single-fit assessment, ensemble generation of fits, clustering, and multiple and consensus scoring, as well as plots and output files for visualization purposes to help the user in analysing rigid and flexible fits. The modular nature of TEMPy helps the integration of scoring and assessment of fits into large pipelines, making it a tool suitable for both novice and expert structural biologists.

Mutations in SLC12A5 in epilepsy of infancy with migrating focal seizures.

The potassium-chloride co-transporter KCC2, encoded by SLC12A5, plays a fundamental role in fast synaptic inhibition by maintaining a hyperpolarizing gradient for chloride ions. KCC2 dysfunction has been implicated in human epilepsy, but to date, no monogenic KCC2-related epilepsy disorders have been described. Here we show recessive loss-of-function SLC12A5 mutations in patients with a severe infantile-onset pharmacoresistant epilepsy syndrome, epilepsy of infancy with migrating focal seizures (EIMFS). Decreased KCC2 surface expression, reduced protein glycosylation and impaired chloride extrusion contribute to loss of KCC2 activity, thereby impairing normal synaptic inhibition and promoting neuronal excitability in this early-onset epileptic encephalopathy.

Conserved mechanisms of microtubule-stimulated ADP release, ATP binding, and force generation in transport kinesins.

Kinesins are a superfamily of microtubule-based ATP-powered motors, important for multiple, essential cellular functions. How microtubule binding stimulates their ATPase and controls force generation is not understood. To address this fundamental question, we visualized microtubule-bound kinesin-1 and kinesin-3 motor domains at multiple steps in their ATPase cycles--including their nucleotide-free states--at ∼ 7 Šresolution using cryo-electron microscopy. In both motors, microtubule binding promotes ordered conformations of conserved loops that stimulate ADP release, enhance microtubule affinity and prime the catalytic site for ATP binding. ATP binding causes only small shifts of these nucleotide-coordinating loops but induces large conformational changes elsewhere that allow force generation and neck linker docking towards the microtubule plus end. Family-specific differences across the kinesin-microtubule interface account for the distinctive properties of each motor. Our data thus provide evidence for a conserved ATP-driven mechanism for kinesins and reveal the critical mechanistic contribution of the microtubule interface.

Allosteric signalling in the outer membrane translocation domain of PapC usher.

PapC ushers are outer-membrane proteins enabling assembly and secretion of P pili in uropathogenic E. coli. Their translocation domain is a large ß-barrel occluded by a plug domain, which is displaced to allow the translocation of pilus subunits across the membrane. Previous studies suggested that this gating mechanism is controlled by a ß-hairpin and an α-helix. To investigate the role of these elements in allosteric signal communication, we developed a method combining evolutionary and molecular dynamics studies of the native translocation domain and mutants lacking the ß-hairpin and/or the α-helix. Analysis of a hybrid residue interaction network suggests distinct regions (residue 'communities') within the translocation domain (especially around ß12-ß14) linking these elements, thereby modulating PapC gating. Antibiotic sensitivity and electrophysiology experiments on a set of alanine-substitution mutants confirmed functional roles for four of these communities. This study illuminates the gating mechanism of PapC ushers and its importance in maintaining outer-membrane permeability.

Stepwise visualization of membrane pore formation by suilysin, a bacterial cholesterol-dependent cytolysin.

Membrane attack complex/perforin/cholesterol-dependent cytolysin (MACPF/CDC) proteins constitute a major superfamily of pore-forming proteins that act as bacterial virulence factors and effectors in immune defence. Upon binding to the membrane, they convert from the soluble monomeric form to oligomeric, membrane-inserted pores. Using real-time atomic force microscopy (AFM), electron microscopy (EM), and atomic structure fitting, we have mapped the structure and assembly pathways of a bacterial CDC in unprecedented detail and accuracy, focussing on suilysin from Streptococcus suis. We show that suilysin assembly is a noncooperative process that is terminated before the protein inserts into the membrane. The resulting ring-shaped pores and kinetically trapped arc-shaped assemblies are all seen to perforate the membrane, as also visible by the ejection of its lipids. Membrane insertion requires a concerted conformational change of the monomeric subunits, with a marked expansion in pore diameter due to large changes in subunit structure and packing.