Zimmermann, Michael T

ABSTRACT The cohesin complex, composed of SMC1, SMC3, RAD21, and STAG1/STAG2, is essential for chromosome cohesion, DNA repair, and transcriptional regulation. Pathogenic variants in cohesin components cause cohesinopathies. The classical characteristics of cohesinopathies include developmental delay (DD), intellectual disability (ID), feeding difficulties, hypotonia, short stature, hearing loss, and dysmorphic features. Here, we present a 5‐year‐old boy with classical cohesinopathy features, including DD/ID and feeding difficulties, along with non‐classical features such as hypoglycemia, recurrent fever, and immunodeficiency. Trio exome sequencing identified a novel de novo missense variant of uncertain significance (NM_005862.3:c.643G>A(p.Val215Ile)) in the STAG1 gene. The variant localizes to the RAD21 interaction interface, and molecular dynamics (MD) simulations revealed conformational changes comparable to other STAG1 variants reported as likely pathogenic in patients, supporting a deleterious effect which may disrupt the STAG1‐RAD21 interaction interface. This case expands the phenotypic and molecular spectrum of STAG1‐related cohesinopathy and advances our understanding of the disease mechanism.

Background/Objectives: Pancreatic ductal adenocarcinoma (PDAC) arises predominantly from activating KRAS mutations, yet individual genetic variants differ markedly in signaling output and clinical impact. G12D, the most prevalent variant, strongly drives oncogenic programs, whereas G12R signals less efficiently through the AKT and ERK pathways and is associated with longer patient survival than G12D-driven PDAC. Methods: To elucidate how these differences influence early cellular transformation, we expressed a panel of KRAS mutants in non-cancerous pancreatic ductal epithelial cells as a model of early PDAC initiation and profiled transcriptional and phospho-proteomic responses. We next examined whether epigenetic differences translate into mutation-specific changes in nuclear organization using quantitative imaging of G12D- and G12R-expressing nuclei at 24 and 48 h. Results: Each variant established a unique regulatory program enriched for chromatin remodelers, histone modifiers, and nuclear structural factors, indicating that variant-specific KRAS signaling rapidly develops divergent epigenetic states. Integrated transcriptomic and phospho-proteomic analyses identified G12D and G12R as the most divergent variants. G12D induced pronounced nuclear remodeling, including increased nuclear size, irregular morphology, and reorganization of the nucleolus and spliceosome, consistent with extensive chromatin and transcriptional reprogramming. In contrast, G12R elicited a weaker response, with minimal or delayed structural changes. Conclusions: Together, these findings demonstrate that KRAS mutational context in pancreatic ductal epithelial cells shapes early transcriptional reprogramming that actively remodels nuclear architecture and nuclear sub-compartments. This work establishes nuclear structural remodeling as a structural state of KRAS-driven epigenetic dysregulation during PDAC initiation.

Background High-impact pathogenic variants in more than a thousand genes are involved in Mendelian forms of neurodevelopmental disorders (NDD). Methods This study describes the molecular and clinical characterisation of 28 probands with NDD harbouring heterozygous AGO1 coding variants, occurring de novo for all those whose transmission could have been verified (26/28). Results A total of 15 unique variants leading to amino acid changes or deletions were identified: 12 missense variants, two in-frame deletions of one codon, and one canonical splice variant leading to a deletion of two amino acid residues. Recurrently identified variants were present in several unrelated individuals: p.(Phe180del), p.(Leu190Pro), p.(Leu190Arg), p.(Gly199Ser), p.(Val254Ile) and p.(Glu376del). AGO1 encodes the Argonaute 1 protein, which functions in gene-silencing pathways mediated by small non-coding RNAs. Three-dimensional protein structure predictions suggest that these variants might alter the flexibility of the AGO1 linker domains, which likely would impair its function in mRNA processing. Affected individuals present with intellectual disability of varying severity, as well as speech and motor delay, autistic behaviour and additional behavioural manifestations. Conclusion Our study establishes that de novo coding variants in AGO1 are involved in a novel monogenic form of NDD, highly similar to the recently reported AGO2 -related NDD.

Abstract Motivation Genomic data are prevalent, leading to frequent encounters with uninterpreted variants or mutations with unknown mechanisms of effect. Researchers must manually aggregate data from multiple sources and across related proteins, mentally translating effects between the genome and proteome, to attempt to understand mechanisms. Materials and methods P2T2 presents diverse data and annotation types in a unified protein-centric view, facilitating the interpretation of coding variants and hypothesis generation. Information from primary sequence, domain, motif, and structural levels are presented and also organized into the first Paralog Annotation Analysis across the human proteome. Results Our tool assists research efforts to interpret genomic variation by aggregating diverse, relevant, and proteome-wide information into a unified interactive web-based interface. Additionally, we provide a REST API enabling automated data queries, or repurposing data for other studies. Conclusion The unified protein-centric interface presented in P2T2 will help researchers interpret novel variants identified through next-generation sequencing. Code and server link available at github.com/GenomicInterpretation/p2t2.

ABSTRACT Inactivating mutations of the ENPP1 gene are associated with generalized arterial calcification of infancy (GACI) and less often autosomal-recessive hypophosphatemic rickets type 2 (ARHR2). We aimed to investigate the spectrum of phenotypes in a family with monoallelic and biallelic mutations of ENPP1 after identification through whole exome sequencing of a 54-year-old female with biallelic mutation of ENPP1, c.323G > T; p.Cys108Phe and c.1441C > T; p.Arg481Trp. Including the proband, 2 subjects had biallelic mutations, 5 had monoallelic mutations, and 2 had no mutation of ENPP1. The maternal mutation, a known pathogenic variant associated with GACI, was found in 3 subjects with monoallelic mutations, while the paternal mutation, which was not previously reported, was present in 2 subjects with monoallelic mutations. Both subjects with biallelic mutations had bowing of bilateral femurs, periarticular mineral deposition, normocalcemic primary hyperparathyroidism with multigland parathyroidectomy, increased carotid intima-media thickness, and enthesopathy was also noted in one subject. Intact FGF23 was elevated in both subjects with biallelic mutations, while C-terminal FGF23 was only elevated in one and PPi was reduced in one. Subjects with monoallelic mutations did not have periarticular calcifications or bone deformities. To conclude, patients with biallelic GACI causing mutations can survive well into adulthood, and despite the same biallelic ENPP1 pathogenic variants, clinical and biochemical manifestations can significantly differ, and include enthesopathy and primary hyperparathyroidism, which have not been previously described. Although carriers of monoallelic ENPP1 variants appear unaffected by classic disease manifestations, there may be subtle biochemical and clinical findings that warrant further investigation. © 2019 American Society for Bone and Mineral Research.

Background Identifying the biologic functions of groups of genes identified in high-throughput studies currently requires considerable time and/or bioinformatics experience. This is due in part to each resource housed within separate databases, requiring users to know about them, and integrate across them. Time consuming and often repeated for each study, integrating across resources and merging with data under study is an increasingly common bioinformatics task. Methods We developed an open-source R software package for assisting researchers in annotating their genesets with functions, pathways, and their interconnectivity across a diversity of network resources. Results We present rapid integration of term annotation and network resources (RITAN) for the rapid and comprehensive annotation of a list of genes using functional term and pathway resources and their relationships among each other using multiple network biology resources. Currently, and to comply with data redistribution policies, RITAN allows rapid access to 16 term annotations spanning gene ontology, biologic pathways, and immunologic modules, and nine network biology resources, with support for user-supplied resources; we provide recommendations for additional resources and scripts to facilitate their addition to RITAN. Having the resources together in the same system allows users to derive novel combinations. RITAN has a growing set of tools to explore the relationships within resources themselves. These tools allow users to merge resources together such that the merged annotations have a minimal overlap with one another. Because we index both function annotation and network interactions, the combination allows users to expand small groups of genes using links from biologic networks—either by adding all neighboring genes or by identifying genes that efficiently connect among input genes—followed by term enrichment to identify functions. That is, users can start from a core set of genes, identify interacting genes from biologic networks, and then identify the functions to which the expanded list of genes contribute. Conclusion We believe RITAN fills the important niche of bridging the results of high-throughput experiments with the ever-growing corpus of functional annotations and network biology resources. Availability Rapid integration of term annotation and network resources is available as an R package at github.com/MTZimmer/RITAN and BioConductor.org .

AbstractDiffuse idiopathic skeletal hyperostosis (DISH) is a disorder principally characterized by calcification and ossification of spinal ligaments and entheses. Fibrodysplasia ossificans progressiva (FOP) is a rare autosomal dominant disabling disorder characterized by progressive ossification of skeletal muscle, fascia, tendons, and ligaments. These conditions manifest phenotypic overlap in the ossification of tendons and ligaments. We describe herein a patient with DISH, exhibiting heterotopic ossification of the posterior longitudinal ligament where clinical whole exome sequencing identified a variant within ACVR1, a gene implicated in FOP. This variant, p.K400E, is a novel variant, not identified previously, and occurs in a highly conserved region across orthologs. We used sequence‐based predicative algorithms, molecular modeling, and molecular dynamics simulations, to test the potential for p.K400E to alter the structure and dynamics of ACVR1. We applied the same modeling and simulation methods to established FOP variants, to identify the detailed effects that they have on the ACVR1 protein, as well as to act as positive controls against which the effects of p.K400E could be evaluated. Our in silico molecular analyses support p.K400E as altering the behavior of ACVR1. In addition, functional testing to measure the effect of this variant on BMP‐pSMAD 1/5/8 target genes was carried out which revealed this variant to cause increased ID1 and Msx2 expression compared with the wild‐type receptor. This analysis supports the potential for the variant of uncertain significance to contribute to the patient's phenotype.

AbstractHeterochromatin protein 1α (HP1α) is a protein that mediates cancer‐associated processes in the cell nucleus. Proteomic experiments, reported here, demonstrate that HP1α complexes with importin α (IMPα), a protein necessary for its nuclear transport. This data is congruent with Simple Linear Motif (SLiM) analyses that identify an IMPα‐binding motif within the linker that joins the two globular domains of this protein. Using molecular modeling and dynamics simulations, we develop a model of the IMPα‐HP1α complex and investigate the impact of phosphorylation and genomic variants on their interaction. We demonstrate that phosphorylation of the HP1α linker likely regulates its association with IMPα, which has implications for HP1α access to the nucleus, where it functions. Cancer‐associated genomic variants do not abolish the interaction of HP1α but instead lead to rearrangements where the variant proteins maintain interaction with IMPα, but with less specificity. Combined, this new mechanistic insight bears biochemical, cell biological, and biomedical relevance.

AbstractAnaplastic large cell lymphomas (ALCLs) represent a relatively common group of T-cell non-Hodgkin lymphomas (T-NHLs) that are unified by similar pathologic features but demonstrate marked genetic heterogeneity. ALCLs are broadly classified as being anaplastic lymphoma kinase (ALK)+ or ALK−, based on the presence or absence of ALK rearrangements. Exome sequencing of 62 T-NHLs identified a previously unreported recurrent mutation in the musculin gene, MSCE116K, exclusively in ALK− ALCLs. Additional sequencing for a total of 238 T-NHLs confirmed the specificity of MSCE116K for ALK− ALCL and further demonstrated that 14 of 15 mutated cases (93%) had coexisting DUSP22 rearrangements. Musculin is a basic helix-loop-helix (bHLH) transcription factor that heterodimerizes with other bHLH proteins to regulate lymphocyte development. The E116K mutation localized to the DNA binding domain of musculin and permitted formation of musculin–bHLH heterodimers but prevented their binding to authentic target sequence. Functional analysis showed MSCE116K acted in a dominant-negative fashion, reversing wild-type musculin-induced repression of MYC and cell cycle inhibition. Chromatin immunoprecipitation–sequencing and transcriptome analysis identified the cell cycle regulatory gene E2F2 as a direct transcriptional target of musculin. MSCE116K reversed E2F2-induced cell cycle arrest and promoted expression of the CD30–IRF4–MYC axis, whereas its expression was reciprocally induced by binding of IRF4 to the MSC promoter. Finally, ALCL cells expressing MSCE116K were preferentially targeted by the BET inhibitor JQ1. These findings identify a novel recurrent MSC mutation as a key driver of the CD30–IRF4–MYC axis and cell cycle progression in a unique subset of ALCLs.