key: cord-0708593-cxpb3a38 authors: Pierzynowska, Karolina; Gaffke, Lidia; Węgrzyn, Grzegorz title: Transcriptomic analyses suggest that mucopolysaccharidosis patients may be less susceptible to COVID‐19 date: 2020-08-15 journal: FEBS Lett DOI: 10.1002/1873-3468.13908 sha: e7bb1dc45546e463818657031eaf3786b4fef57e doc_id: 708593 cord_uid: cxpb3a38 We used transcriptomic (RNA‐seq) analyses to determine if patients suffering from all types and subtypes of mucopolysaccharidosis (MPS), a severe inherited metabolic disease, may be more susceptible to COVID‐19. The expression levels of genes encoding proteins potentially involved in SARS‐CoV‐2 development were estimated in MPS cell lines. Four genes (GTF2F2, RAB18, TMEM97, PDE4DIP) coding for proteins potentially facilitating virus development were down‐regulated, while two genes (FBN1, MFGE8), the products of which potentially interfere with virus propagation, were up‐regulated in most MPS types. Although narrowing of respiratory tract and occurrence of thick mucus, characteristic of MPS, are risk factors for COVID‐19, transcriptomic analyses suggest that MPS cells might be less, rather than more, susceptible to SARS‐CoV‐2 infection. Recently started and still ongoing outbreak of coronavirus disease 2019, or COVID-19 (as officially named by World Health Organization), is caused by SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) [1] . Over 16 million patients have been diagnosed by the date of submission of this report, including over 650,000 fatal cases, and this pandemic spread of the disease caused serious disturbance not only in medical care but also in socio-economic aspects of human life [2] [3] [4] . Apart from general understanding of the mechanisms of the disease, which should lead to development of effective therapies and prevention strategies, including production of vaccines, identification of groups of people either especially susceptible or resistant to SARS-CoV-2 infection appears to be crucial for appropriate management of the epidemy. On the other hand, despite thousands of articles published on general or specific aspects of COVID-19, only a small fraction of publications concerns problems of patients suffering from other diseases as risk groups. Recently published reports include papers analyzing problems linked to SARS-CoV-2mediated infections of patients with as different disorders as dental problems [5] , cardiovascular diseases [6, 7] , dermatologic diseases [8] , or gastrointestinal dysfunctions [9] , to list just a few. However, patients suffering from genetic diseases may be a special risk group, not only as potentially affected persons but also from the epidemiological point of view, due to severity of the primary conditions and requirement for constant care by others. On the other hand, primary genetic defects often cause multiple changes in expression of a battery of genes which may potentially affect, either positively or negatively, viral development. Mucopolysaccharidoses (MPS) are a group of 11 genetic diseases caused by mutations in one of genes coding for lysosomal enzymes involved in degradation of glycosaminoglycans (GAGs), complex carbohydrates [10] . GAG storage leads to development of severe symptoms which progressively appear in the course of the disease development, and the clinical patient's state worsens gradually, causing a death usually by the end of the second decade [11] . Most organs are affected in MPS patients, including heart, gastrointestinal tract, respiratory system, skin, bones, This article is protected by copyright. All rights reserved muscles, central nervous system (in some types) and others [11] . Thus, one might suspect that such patients might be a high risk group for COVID-19, however, no reports were published to date on COVID-19 in MPS. On the other hand, recent studies indicated that in cells of MPS patients, there are hundreds of genes which expression is significantly changed relative to healthy persons [12] [13] [14] . This might potentially modify susceptibility of these patients to infection by SARS-CoV-2. Therefore, in this report, we have analyzed transcriptomes of cells derived from patients suffering from all MPS types relative to control cell (derived from a healthy person) in the light of potential susceptibility to COVID-19. Cultures of fibroblast of all MPS types, as well as control fibroblasts (HDFa line), were This article is protected by copyright. All rights reserved RNA isolation was conducted in 4 biological repeats for each cell line (i.e. 4 independent cultures from different passages). The procedure was performed exactly as described previously [12] . RNA-seq analyses were performed as described previously [12] . Raw data are deposited at NCBI Sequence Read Archive (SRA), under accession no. PRJNA562649. Fibroblasts (6 x 10 5 cells) were passaged on 10 cm-diameter plates, and allowed to attach overnight, and the cultivation was continued for 48 h. A solution composed of 1% Triton X-100, 0.5 mM EDTA, 150 mM NaCl, 50 mM Tris, pH 7.5, and a mixture of protease and phosphatase inhibitors (Roche Applied Science, Penzberg, Germany) was used for cell lysis. One-way ANOVA on log 2 (1 + x) values revealing normal continuous distribution was employed to determine statistical significance. The Benjamini-Hochberg method was used to estimate the false discovery rate (FDR). Post hoc Student's t-test with Bonferroni correction was used for comparisons between two groups. The R software v3.4.3 was used for all statistical analyses. This article is protected by copyright. All rights reserved We have analyzed transcriptomes of fibroblasts derived from patients suffering from all known MPS types (I, II, IIIA, IIIB, IIIC, IIID, IVA, IVB, VI, VII, and IX) , obtained in the RNAseq analysis, relative to control fibroblasts (HDFa cell line) (SRA accession no. PRJNA562649). The quality of the RNA-seq data was confirmed previously and assessed to be reliable on the basis of RT-qPCR analyses of selected genes [12] [13] [14] . Because of the severity of the disease, fibroblasts are most often used type of cells in genetic and genomic studies on MPS; in fact availability of other biological material is technically and ethically poor, due to conditions of patients. We have considered human genes coding for proteins which have been reported recently as those interacting with polypeptides encoded by SARS-CoV-2 [15, 16] . The newest database lists 332 physical interactions between proteins encoded by human and the virus [16] . We found that among genes coding for human proteins involved in interactions with SARS-CoV-2's proteins, there were between 13 and 41 transcripts (depending on MPS type) which levels were significantly up-or down-regulated in MPS fibroblasts relative to control cells ( Figure 1 ). The full list of such transcripts, which expression was significantly changed in fibroblasts of particular MPS types, is included in Supplementary material (Table S1 ). It is worth noting that vast majority of changes in expression of almost every single gene is the same for the direction of the change, irrespective of the MPS type, i.e. every gene is either down-or up-regulated in all MPS types. These results may suggest that the response of MPS patients to SARS-CoV-2 infection may considerably differ from that of general population due to differences in intracellular virus development. Interestingly, in most MPS types, number of down-regulated genes was higher than up-regulated genes ( Figure 1 , Table S1 ). Unfortunately, expressions of two human genes crucial for the virus adsorption and intracellular penetration, ACE2 and TMPRSS2, coding for the SARS-CoV-2 receptor and specific protease, respectively [1] [2] [3] [4] , were at very low levels in fibroblasts (FPKM values <1), precluding a reliable analysis. Nevertheless, it might be worth mentioning that while low mRNA levels of ACE2 could be detected in HDFa cells, the determined corresponding values for vast majority of MPS cell lines were zero, except MPS IVA and VI, where negligible amounts of the transcripts could be found. If down-regulation of ACE2 occurs in MPS cells indeed, the patients might be less susceptible for COVID-2, however, it should be underlined that such an assumption is highly This article is protected by copyright. All rights reserved speculative. Despite these difficulties, we aimed to detect ACE2 and TMPRSS2 proteins by high sensitive Western-blotting (using the WES system). Although we failed to visualize ACE2 protein (data not shown), we were able to detect TMPRSS2 ( Figure 2E ). Quantitative analysis indicated that the level of this protein was decreased in MPS I, II and IIIB fibroblasts while increased in other MPS types relative to control cells ( Figure 2F ). This might suggest differential efficiency of adsorption of SARS-CoV-2 on the cells, depending on the MPS type. When considering genes which were up-or down-regulated in most MPS types, we could find 6 genes which expression was changed in at least 7 types relative to control cells. Expression analyses of these genes are presented as heat-maps (Figure 2A) , and fold-changes of levels of particular transcripts are demonstrated in Table 1 . Four of these genes (GTF2F2, RAB18, TMEM97, and PDE4DIP) were down-regulated, and two (FBN1 and MFGE8) were up-regulated in all MPS types. We asked if transcriptomic data could be confirmed at the level of abundance of corresponding gene products. Therefore, using Western-blotting, we tested levels of selected proteins, encoded by genes which revealed significantly changed expression at the mRNA level in at least 7 MPS types. The results presented in Figure 2B , 2C and 2D indicated that levels of the PDE4DIP protein were decreased in fibroblast of all MPS types but MPS IIIC, and levels of FBN1 were increased in all MPS types but MPS IIIC and IIID, relative to control cells. Therefore, we concluded that the presented transcriptomic data are reliable and observed mRNA levels generally correspond to levels of gene products. Again, the Western-blotting results confirmed that expressions of tested genes are similarly regulated (up-or down-, depending on the gene) in most MPS types. Such uniform distribution of up-and down-regulated genes among MPS types indicates that the observed changes are not coincidental, but rather they represent a general tendency in MPS cells (this conclusion is corroborated by results for other genes, presented in Table S1 ). The GTF2F2 gene codes for the general transcription factor IIF subunit 2, known as TFIIF [17] , thus, it is likely involved in the viral gene expression positive control. RAB18 is a Rasrelated small GTPase involved in regulation of membrane trafficking and vesicular transport [18] . TMEM97 is a transmembrane protein which regulates cholesterol levels [19] . The PDE4DIP gene codes for the phosphodiesterase 4D-interacting protein which is involved in Golgi apparatus functions [20] . One may predict that all these proteins are required for effective SARS-CoV-2 development, thus, down-regulation of their genes might cause lower susceptibility to COVID-2. This article is protected by copyright. All rights reserved The FBN1 gene codes for fibrillin-1, a protein which is proteolytically cleaved to generate an extracellular matrix glycoprotein and asprosin, involved in glucose homeostasis [21] . MFGE8 also codes for a polyprotein which after cleavage in converted to several polypeptides, among which lactahedrin is the major one; lactahedrin is a membrane glycoprotein promoting phagocytosis of apoptotic cells [22] . Therefore, one might speculate that up-regulation of these two genes in MPS cells could promote anti-viral defense cellular mechanisms, again lowering susceptibility of the patients to COVID-2. We have also analyzed transcripts which levels were particularly significantly changed (log 2 fold change (FC)>2) in certain MPS types. The results are shown in Table 2 . In individual MPS types, following genes were down-regulated: CAV1, UBE2I, HOXC6, TMEM97, NUP88, and ARL6IP6. On the other hand, MFGE8 and ATF5 were up-regulated. CAV1 codes for caveolin-1, a protein involved in promoting cell cycle progression [23] . UBE2I is a ubiquitin conjugating enzyme E2I, involved in protein degradation mechanism, while it can also stimulate transcription [24] . HOXC6 encodes a transcription factor involved in morphogenesis [25] . The NUP88 gene product is one of nucleoporins, required for communication of the nucleus with cytoplasm, particularly for mRNA transport [26] . ARL6IP6 codes for ADP ribosylation factor like GTPase 6 interacting protein 6, a polypeptide involved in the protein modification processes [27] . One may predict that, together with TMEM97 (described above), In summary, results of transcriptomic analyses, presented in this report, might suggest that cells of patients suffering from MPS can reveal lower susceptibility to infection by SARS-CoV-2 due to specific changes in expression of genes coding for proteins involved in interactions with viral proteins, despite increased levels of TMPRSS2 in some MPS types. On the other hand, it is necessary to take into consideration the fact that MPS patients are generally more sensitive to various infectious agents due to secondary anatomical and physiological changes, appearing as a This article is protected by copyright. All rights reserved result of the primary GAG storage [10, 11] . These changes include narrowing of respiratory tract and occurrence of thick mucus which can be significant risk factors when considering COVID-19. Thus, one should still recommend particular caution when protecting MPS patients against SARS-CoV-2 infection. Moreover, in MPS types in which enzyme replacement therapy is available, potential interruption of the treatment (due to strategies of some hospitals, reluctant to the whole day procedure of intravenous administration of the drug to severely affected patients in the clinic) may have fatal consequences, as published reports clearly indicated that cessation of this therapy results in dramatic worsening of conditions of patients [29] [30] [31] which, in turn, would make them especially endangered by COVID-19. In fact, frequent and severe respiratory tract infections were reported as being among characteristic symptoms of enzyme replacement therapy interruption in MPS patients [29] [30] [31] . The question what is the specific cause of changes in expression of tested genes in MPS cells remains unanswered. As indicated in recently published reports [12] [13] [14] , expression of hundreds of genes is changed in MPS cells relative to controls. Specific mechanisms of these regulations remain to be elucidated, however, it was proposed that they reflect secondary and tertiary changes in cellular processes, after primary GAG storage [32] . On the other hand, it has been indicated that severity of MPS symptoms in patients does not simply correlate with GAG levels in organisms of such patients [33] . Therefore, specific regulations of cellular processes, including control of gene expression, in MPS cells, is a complex process, involving a complicated Accepted Article Table 1 . Human genes (which products are potentially involved in interaction with SARS-CoV-2 proteins) with expression significantly changed (at FDR<0.1; p<0.1) in most (7 or more) MPS types relative to control cells (HDFa). Down-regulated genes are marked in blue, and up-regulated genes are marked in red. Non-colored values indicated results in which no statistically significant differences between MPS and HDFa were determined. Transcript I II III A III B III C III D IVA IVB VI VII Transcript I II III A III B III C III D IVA IVB VI VII IX CAV1 ↓ - - - - - - ↓ - - - Accepted Article UBE2I - ↓ - - - - - - ↓ - - HOXC6 - - - - ↓ - - - - - - HOXC6 - - - - ↓ - - ↓ - - - TMEM97 - - - - - - - - - - ↓ NUP88 - - - - - ↓ - - - - - ARL6IP6 - ↓ - ↓ ↓ - - - - - - MFGE8 ↑ - ↑ ↑ - - - ↑ - - - MFGE8 ↑ - ↑ ↑ - - - ↑ - - - MFGE8 ↑ - ↑ ↑ - - - ↑ - - - MFGE8 ↑ - ↑ ↑ - - - - - - - ATF5 - - - - - - - - - - ↑ ATF5 - - - - - - - - - - ↑ FIGURE LEGENDS This article is protected by copyright. All rights reserved Western-blotting of two proteins, PDE4DIP and FBN1 (with GAPDH as a loading control), and quantitative analysis (from three independent experiments) is presented in panels C and D, respectively (with asterisks indicating statistically significant differences, at p<0.05, relative to control). Panel E shows representative results of Western-blotting of the TMPRSS2 protein (with GAPDH as a loading control), and quantitative analysis (from three independent experiments) is presented in panel F (with asterisks indicating statistically significant differences, at p<0.05, relative to control). Supplementary material - Table S1 . Human genes (which products are potentially involved in interaction with SARS-CoV-2 proteins) with changed expression (at FDR<0.1; p<0.1) in specific MPS type relative to control cells (HDFa). Up-arrowed and down-arrowed symbols indicate upand down-regulation, respectively, while minus symbols indicate no significant changes. For some genes, more than one transcript were identified. COVID-19 outbreak: an overview Epidemiological and clinical aspects of COVID-19; a narrative review Coronavirus disease 2019: a bibliometric analysis and review Coronavirus: update related to the current outbreak of COVID-19 Coronavirus disease 19 (COVID-19): implications for clinical dental care Cardiac and arrhythmic complications in Covid-19 patients COVID-19 and the heart The effect of quarantine and isolation for COVID-19 in general population and dermatologic treatments COVID-19 and the gastrointestinal tract: more than meets the eye Pathogenesis of mucopolysaccharidoses, an update A basic understanding of mucopolysaccharidosis: Incidence, clinical features, diagnosis, and management Underestimated aspect of mucopolysaccharidosis pathogenesis: global changes in cellular processes revealed by transcriptomic studies Genetic base of behavioral disorders in mucopolysaccharidoses: transcriptomic studies Expression of genes involved in apoptosis is dysregulated in mucopolysaccharidoses as revealed by pilot transcriptomic analyses Master regulator analysis of the SARS-CoV-2/human interactome A SARS-CoV-2 protein interaction map reveals targets for drug repurposing Accepted Article This article is protected by copyright. All rights reserved Genetic characterization of general transcription factors TFIIF and TFIIB of Homo sapiens sapiens Rab18: new insights into the function of an essential protein Small-molecule modulators of sigma1 and sigma2/TMEM97 in the context of cancer: foundational concepts and emerging themes Molecule of the month Fibrillin-1 and fibrillin-1-derived asprosin in adipose tissue function and metabolic disorders Lactadherin: an unappreciated haemostasis regulator and potential therapeutic agent Caveolin-1, cancer and therapy resistance UBE2I (UBC9), a SUMO-conjugating enzyme, localizes to nuclear speckles and stimulates transcription in mouse oocytes The HOXC6 homeodomain-containing proteins The human homologue of yeast CRM1 is in a dynamic subcomplex with CAN/Nup214 and a novel nuclear pore component Nup88 Identification of new transmembrane proteins concentrated at the nuclear envelope using organellar proteomics of mesenchymal cells The transcription factor ATF5: role in cellular differentiation, stress responses, and cancer Rapid deterioration of a patient with mucopolysaccharidosis type I during interruption of enzyme replacement therapy Effect of rapid cessation of enzyme replacement therapy: a report of 5 cases and a review of the literature Effect of rapid cessation of enzyme replacement therapy: a report of 5 more cases Changes in cellular processes occurring in mucopolysaccharidoses as underestimated pathomechanisms of these diseases Correlation between severity of mucopolysaccharidoses and combination of the residual enzyme activity and efficiency of glycosaminoglycan synthesis