Regulatory architecture of the RCA gene cluster captures an intragenic TAD boundary, CTCF-Mediated chromatin looping and a long-range intergenic enhancer

The Regulators of Complement Activation (RCA) gene cluster comprises several tandemly arranged genes with shared functions within the immune system. RCA members, such as complement receptor 2 (CR2), are well-established susceptibility genes in complex autoimmune diseases. Altered expression of RCA genes has been demonstrated at both the functional and genetic level, but the mechanisms underlying their regulation are not fully characterised. We aimed to investigate the structural organisation of the RCA gene cluster to identify key regulatory elements that influence the expression of CR2 and other genes in this immunomodulatory region. Using 4C, we captured extensive CTCF-mediated chromatin looping across the RCA gene cluster in B cells and showed these were organised into two topologically associated domains (TADs). Interestingly, an inter-TAD boundary was located within the CR1 gene at a well-characterised segmental duplication. Additionally, we mapped numerous gene-gene and gene-enhancer interactions across the region, revealing extensive co-regulation. Importantly, we identified an intergenic enhancer and functionally demonstrated this element upregulates two RCA members (CR2 and CD55) in B cells. We have uncovered novel, long-range mechanisms whereby autoimmune disease susceptibility may be influenced by genetic variants, thus highlighting the important contribution of chromatin topology to gene regulation and complex genetic disease.This work was supported by the National Institutes of Health [R01 AI24717 to JH], the Australian Government Research Training Program Scholarship at the University of Western Australia [to JC and JSC], the Spanish Government [BFU2016-74961-P to JG-S] and an institutional grant Unidad de Excelencia María de Maeztu [MDM-206-0687 to the Department of Gene Regulation and Morphogenesis, Centro Andaluz de Biología del Desarrol]

Nebulin nemaline myopathy recapitulated in a compound heterozygous mouse model with both a missense and a nonsense mutation in Neb

Nemaline myopathy (NM) caused by mutations in the gene encoding nebulin (NEB) accounts for at least 50% of all NM cases worldwide, representing a significant disease burden. Most NEB-NM patients have autosomal recessive disease due to a compound heterozygous genotype. Of the few murine models developed for NEB-NM, most are Neb knockout models rather than harbouring Neb mutations. Additionally, some models have a very severe phenotype that limits their application for evaluating disease progression and potential therapies. No existing murine models possess compound heterozygous Neb mutations that reflect the genotype and resulting phenotype present in most patients. We aimed to develop a murine model that more closely matched the underlying genetics of NEB-NM, which could assist elucidation of the pathogenetic mechanisms underlying the disease. Here, we have characterised a mouse strain with compound heterozygous Neb mutations; one missense (p.Tyr2303His), affecting a conserved actin-binding site and one nonsense mutation (p.Tyr935*), introducing a premature stop codon early in the protein. Our studies reveal that this compound heterozygous model, Neb(Y2303H, Y935X), has striking skeletal muscle pathology including nemaline bodies. In vitro whole muscle and single myofibre physiology studies also demonstrate functional perturbations. However, no reduction in lifespan was noted. Therefore, Neb(Y2303H,Y935X) mice recapitulate human NEB-NM and are a much needed addition to the NEB-NM mouse model collection. The moderate phenotype also makes this an appropriate model for studying NEB-NM pathogenesis, and could potentially be suitable for testing therapeutic applications.Peer reviewe

Nebulin nemaline myopathy recapitulated in a compound heterozygous mouse model with both a missense and a nonsense mutation in Neb

Nemaline myopathy (NM) caused by mutations in the gene encoding nebulin (NEB) accounts for at least 50% of all NM cases worldwide, representing a significant disease burden. Most NEB-NM patients have autosomal recessive disease due to a compound heterozygous genotype. Of the few murine models developed for NEB-NM, most are Neb knockout models rather than harbouring Neb mutations. Additionally, some models have a very severe phenotype that limits their application for evaluating disease progression and potential therapies. No existing murine models possess compound heterozygous Neb mutations that reflect the genotype and resulting phenotype present in most patients. We aimed to develop a murine model that more closely matched the underlying genetics of NEB-NM, which could assist elucidation of the pathogenetic mechanisms underlying the disease. Here, we have characterised a mouse strain with compound heterozygous Neb mutations; one missense (p.Tyr2303His), affecting a conserved actin-binding site and one nonsense mutation (p.Tyr935*), introducing a premature stop codon early in the protein. Our studies reveal that this compound heterozygous model, Neb(Y2303H, Y935X), has striking skeletal muscle pathology including nemaline bodies. In vitro whole muscle and single myofibre physiology studies also demonstrate functional perturbations. However, no reduction in lifespan was noted. Therefore, Neb(Y2303H,Y935X) mice recapitulate human NEB-NM and are a much needed addition to the NEB-NM mouse model collection. The moderate phenotype also makes this an appropriate model for studying NEB-NM pathogenesis, and could potentially be suitable for testing therapeutic applications.Peer reviewe

A KLHL40 3’ UTR splice-altering variant causes milder NEM8, an under-appreciated disease mechanism

Nemaline myopathy 8 (NEM8) is typically a severe autosomal recessive disorder associated with variants in the kelch-like family member 40 gene (KLHL40). Common features include fetal akinesia, fractures, contractures, dysphagia, respiratory failure, and neonatal death. Here, we describe a man in his 20s with relatively mild NEM8. He presented with hypotonia and bilateral femur fractures at birth, later developing bilateral Achilles’ contractures, scoliosis, and elbow and knee contractures. He had walking difficulties throughout childhood and became wheelchair bound during adolescence after prolonged immobilisation. Muscle MRI during adolescence indicated prominent fat replacement in his pelvic girdle, posterior compartments of thighs, and vastus intermedius. Muscle biopsy revealed nemaline bodies and intranuclear rods. RNA sequencing and western blotting of patient skeletal muscle indicated significant reduction in KLHL40 mRNA and protein respectively. Using gene panel screening, exome sequencing and RNA sequencing, we identified compound heterozygous variants in KLHL40; a truncating 10.9 kb deletion in trans with a likely pathogenic variant (c.*152G>T) in the 3’ untranslated region (UTR). Computational tools SpliceAI and Introme predicted the c.*152G>T variant created a cryptic donor splice site. RNA-seq and in vitro analyses indicated that the c.*152G>T variant induces multiple de novo splicing events that likely provoke nonsense mediated decay of KLHL40 mRNA explaining the loss of mRNA expression and protein abundance in the patient. Analysis of 3’ UTR variants in ClinVar suggests SNPs that introduce aberrant 3’ UTR splicing may be underrecognised in Mendelian disease. We encourage consideration of this mechanism during variant curation.This study was funded by an Australian NHMRC Investigator Grant (APP2007769), Fellowship (APP1117510) and Ideas Grant (APP2002640). The study was supported by an Australian Government Research Training Program (RTP) Scholarship, as well as resources provided by the Pawsey Supercomputing Centre with funding from the Australian Government and the Government of Western Australia.N

A KLHL40 3' UTR splice-altering variant causes milder NEM8, an under-appreciated disease mechanism

Nemaline myopathy 8 (NEM8) is typically a severe autosomal recessive disorder associated with variants in the kelch-like family member 40 gene (KLHL40). Common features include fetal akinesia, fractures, contractures, dysphagia, respiratory failure and neonatal death. Here, we describe a 26-year-old man with relatively mild NEM8. He presented with hypotonia and bilateral femur fractures at birth, later developing bilateral Achilles’ contractures, scoliosis, and elbow and knee contractures. He had walking difficulties throughout childhood and became wheelchair bound from age 13 after prolonged immobilization. Muscle magnetic resonance imaging at age 13 indicated prominent fat replacement in his pelvic girdle, posterior compartments of thighs and vastus intermedius. Muscle biopsy revealed nemaline bodies and intranuclear rods. RNA sequencing and western blotting of patient skeletal muscle indicated significant reduction in KLHL40 mRNA and protein, respectively. Using gene panel screening, exome sequencing and RNA sequencing, we identified compound heterozygous variants in KLHL40; a truncating 10.9 kb deletion in trans with a likely pathogenic variant (c.*152G > T) in the 3′ untranslated region (UTR). Computational tools SpliceAI and Introme predicted the c.*152G > T variant created a cryptic donor splice site. RNA-seq and in vitro analyses indicated that the c.*152G > T variant induces multiple de novo splicing events that likely provoke nonsense mediated decay of KLHL40 mRNA explaining the loss of mRNA expression and protein abundance in the patient. Analysis of 3’ UTR variants in ClinVar suggests variants that introduce aberrant 3’ UTR splicing may be underrecognized in Mendelian disease. We encourage consideration of this mechanism during variant curation.LD is supported by an Australian Government Research Training Program (RTP) Scholarship. GR (Investigator Grant, APP2007769) and NGL (Fellowship APP1117510) are supported by the Australian National Health and Medical Research Council (NHMRC). This work is funded by NHMRC Ideas Grant (APP2002640). This work was supported by resources provided by the Pawsey Supercomputing Centre with funding from the Australian Government and the Government of Western Australia.Peer reviewe

Bi-allelic loss-of-function OBSCN variants predispose individuals to severe recurrent rhabdomyolysis

Rhabdomyolysis is the acute breakdown of skeletal myofibres in response to an initiating factor, most commonly toxins and over exertion. A variety of genetic disorders predispose to rhabdomyolysis through different pathogenic mechanisms, particularly in patients with recurrent episodes. However, most cases remain without a genetic diagnosis. Here we present six patients who presented with severe and recurrent rhabdomyolysis, usually with onset in the teenage years; other features included a history of myalgia and muscle cramps. We identified 10 bi-allelic loss-of-function variants in the gene encoding obscurin (OBSCN) predisposing individuals to recurrent rhabdomyolysis. We show reduced expression of OBSCN and loss of obscurin protein in patient muscle. Obscurin is proposed to be involved in sarcoplasmic reticulum function and Ca2+ handling. Patient cultured myoblasts appear more susceptible to starvation as evidenced by a greater decreased in sarcoplasmic reticulum Ca2+ content compared to control myoblasts. This likely reflects a lower efficiency when pumping Ca2+ back into the sarcoplasmic reticulum and/or a decrease in Ca2+ sarcoplasmic reticulum storage ability when metabolism is diminished. OSBCN variants have previously been associated with cardiomyopathies. None of the patients presented with a cardiomyopathy and cardiac examinations were normal in all cases in which cardiac function was assessed. There was also no history of cardiomyopathy in first degree relatives, in particular in any of the carrier parents. This cohort is relatively young, thus follow-up studies and the identification of additional cases with bi-allelic null OBSCN variants will further delineate OBSCN-related disease and the clinical course of disease. Cabrera-Serrano et al. show that biallelic loss-of-function variants in the gene encoding obscurin (OBSCN) predispose individuals to recurrent and severe episodes of rhabdomyolysis, typically with onset in the teenage years.Peer reviewe

Histo-Blood Group Gene Polymorphisms as Potential Genetic Modifiers of Infection and Cystic Fibrosis Lung Disease Severity

The pulmonary phenotype in cystic fibrosis (CF) is variable; thus, environmental and genetic factors likely contribute to clinical heterogeneity. We hypothesized that genetically determined ABO histo-blood group antigen (ABH) differences in glycosylation may lead to differences in microbial binding by airway mucus, and thus predispose to early lung infection and more severe lung disease in a subset of patients with CF. infection in the severe or mild groups. Multivariate analyses of other clinical phenotypes, including gender, asthma, and meconium ileus demonstrated no differences between groups based on ABH type. infection, nor was there any association with other clinical phenotypes in a group of 808 patients homozygous for the ΔF508 mutation

Mucin Variable Number Tandem Repeat Polymorphisms and Severity of Cystic Fibrosis Lung Disease: Significant Association with MUC5AC

Variability in cystic fibrosis (CF) lung disease is partially due to non-CFTR genetic modifiers. Mucin genes are very polymorphic, and mucins play a key role in the pathogenesis of CF lung disease; therefore, mucin genes are strong candidates as genetic modifiers. DNA from CF patients recruited for extremes of lung phenotype was analyzed by Southern blot or PCR to define variable number tandem repeat (VNTR) length polymorphisms for MUC1, MUC2, MUC5AC, and MUC7. VNTR length polymorphisms were tested for association with lung disease severity and for linkage disequilibrium (LD) with flanking single nucleotide polymorphisms (SNPs). No strong associations were found for MUC1, MUC2, or MUC7. A significant association was found between the overall distribution of MUC5AC VNTR length and CF lung disease severity (p = 0.025; n = 468 patients); plus, there was robust association of the specific 6.4 kb HinfI VNTR fragment with severity of lung disease (p = 6.2 x 10(-4) after Bonferroni correction). There was strong LD between MUC5AC VNTR length modes and flanking SNPs. The severity-associated 6.4 kb VNTR allele of MUC5AC was confirmed to be genetically distinct from the 6.3 kb allele, as it showed significantly stronger association with nearby SNPs. These data provide detailed respiratory mucin gene VNTR allele distributions in CF patients. Our data also show a novel link between the MUC5AC 6.4 kb VNTR allele and severity of CF lung disease. The LD pattern with surrounding SNPs suggests that the 6.4 kb allele contains, or is linked to, important functional genetic variation

A novel lung disease phenotype adjusted for mortality attrition for cystic fibrosis Genetic modifier studies

Genetic studies of lung disease in Cystic Fibrosis are hampered by the lack of a severity measure that accounts for chronic disease progression and mortality attrition. Further, combining analyses across studies requires common phenotypes that are robust to study design and patient ascertainment