AUTS2 Gene: Keys to Understanding the Pathogenesis of Neurodevelopmental Disorders

Single-cell RNA sequencing analysis of vestibular schwannoma reveals functionally distinct macrophage subsets

BACKGROUND

Vestibular schwannomas (VSs) remain a challenge due to their anatomical location and propensity to growth. Macrophages are present in VS but their roles in VS pathogenesis remains unknown.

OBJECTIVES

The objective was to assess phenotypic and functional profile of macrophages in VS with single-cell RNA sequencing (scRNAseq).

METHODS

scRNAseq was carried out in three VS samples to examine characteristics of macrophages in the tumour. RT-qPCR was carried out on 10 VS samples for CD14, CD68 and CD163 and a panel of macrophage-associated molecules.

RESULTS

scRNAseq revealed macrophages to be a major constituent of VS microenvironment with three distinct subclusters based on gene expression. The subclusters were also defined by expression of CD163, CD68 and IL-1β. AREG and PLAUR were expressed in the CD68+CD163+IL-1β+ subcluster, PLCG2 and NCKAP5 were expressed in CD68+CD163+IL-1β- subcluster and AUTS2 and SPP1 were expressed in the CD68+CD163-IL-1β+ subcluster. RT-qPCR showed expression of several macrophage markers in VS of which CD14, ALOX15, Interleukin-1β, INHBA and Colony Stimulating Factor-1R were found to have a high correlation with tumour volume.

CONCLUSIONS

Macrophages form an important component of VS stroma. scRNAseq reveals three distinct subsets of macrophages in the VS tissue which may have differing roles in the pathogenesis of VS.

H2A monoubiquitination: insights from human genetics and animal models

Metazoan development arises from spatiotemporal control of gene expression, which depends on epigenetic regulators like the polycomb group proteins (PcG) that govern the chromatin landscape. PcG proteins facilitate the addition and removal of histone 2A monoubiquitination at lysine 119 (H2AK119ub1), which regulates gene expression, cell fate decisions, cell cycle progression, and DNA damage repair. Regulation of these processes by PcG proteins is necessary for proper development, as pathogenic variants in these genes are increasingly recognized to underly developmental disorders. Overlapping features of developmental syndromes associated with pathogenic variants in specific PcG genes suggest disruption of central developmental mechanisms; however, unique clinical features observed in each syndrome suggest additional non-redundant functions for each PcG gene. In this review, we describe the clinical manifestations of pathogenic PcG gene variants, review what is known about the molecular functions of these gene products during development, and interpret the clinical data to summarize the current evidence toward an understanding of the genetic and molecular mechanism.

Genome-wide detection of positive and balancing signatures of selection shared by four domesticated rainbow trout populations (Oncorhynchus mykiss)

BACKGROUND

Evolutionary processes leave footprints along the genome over time. Highly homozygous regions may correspond to positive selection of favorable alleles, while maintenance of heterozygous regions may be due to balancing selection phenomena. We analyzed data from 176 fish from four disconnected domestic rainbow trout populations that were genotyped using a high-density Axiom Trout genotyping 665K single nucleotide polymorphism array, including 20 from the US and 156 from three French lines. Using methods based on runs of homozygosity and extended haplotype homozygosity, we detected signatures of selection in these four populations.

RESULTS

Nine genomic regions that included 253 genes were identified as being under positive selection in all four populations Most were located on chromosome 2 but also on chromosomes 12, 15, 16, and 20. In addition, four heterozygous regions that contain 29 genes that are putatively under balancing selection were also shared by the four populations. These were located on chromosomes 10, 13, and 19. Regardless of the homozygous or heterozygous nature of the regions, in each region, we detected several genes that are highly conserved among vertebrates due to their critical roles in cellular and nuclear organization, embryonic development, or immunity. We identified new candidate genes involved in rainbow trout fitness, as well as 17 genes that were previously identified to be under positive selection, 10 of which in other fishes (auts2, atp1b3, zp4, znf135, igf-1α, brd2, col9a2, mrap2, pbx1, and emilin-3).

CONCLUSIONS

Using material from disconnected populations of different origins allowed us to draw a genome-wide map of signatures of positive selection that are shared between these rainbow trout populations, and to identify several regions that are putatively under balancing selection. These results provide a valuable resource for future investigations of the dynamics of genetic diversity and genome evolution during domestication.

Machine learning approaches to identify systemic lupus erythematosus in anti-nuclear antibody-positive patients using genomic data and electronic health records

BACKGROUND

Although the 2019 EULAR/ACR classification criteria for systemic lupus erythematosus (SLE) has required at least a positive anti-nuclear antibody (ANA) titer (≥ 1:80), it remains challenging for clinicians to identify patients with SLE. This study aimed to develop a machine learning (ML) approach to assist in the detection of SLE patients using genomic data and electronic health records.

METHODS

Participants with a positive ANA (≥ 1:80) were enrolled from the Taiwan Precision Medicine Initiative cohort. The Taiwan Biobank version 2 array was used to detect single nucleotide polymorphism (SNP) data. Six ML models, Logistic Regression, Random Forest (RF), Support Vector Machine, Light Gradient Boosting Machine, Gradient Tree Boosting, and Extreme Gradient Boosting (XGB), were used to identify SLE patients. The importance of the clinical and genetic features was determined by Shapley Additive Explanation (SHAP) values. A logistic regression model was applied to identify genetic variations associated with SLE in the subset of patients with an ANA equal to or exceeding 1:640.

RESULTS

A total of 946 SLE and 1,892 non-SLE controls were included in this analysis. Among the six ML models, RF and XGB demonstrated superior performance in the differentiation of SLE from non-SLE. The leading features in the SHAP diagram were anti-double strand DNA antibodies, ANA titers, AC4 ANA pattern, polygenic risk scores, complement levels, and SNPs. Additionally, in the subgroup with a high ANA titer (≥ 1:640), six SNPs positively associated with SLE and five SNPs negatively correlated with SLE were discovered.

CONCLUSIONS

ML approaches offer the potential to assist in diagnosing SLE and uncovering novel SNPs in a group of patients with autoimmunity.

Characterization of genomic regions escaping epigenetic reprogramming in sheep

The mammalian genome undergoes two global epigenetic reprogramming events during the establishment of primordial germ cells and in the pre-implantation embryo after fertilization. These events involve the erasure and re-establishment of DNA methylation marks. However, imprinted genes and transposable elements (TEs) maintain their DNA methylation signatures to ensure normal embryonic development and genome stability. Despite extensive research in mice and humans, there is limited knowledge regarding environmentally induced epigenetic marks that escape epigenetic reprogramming in other species. Therefore, the objective of this study was to examine the characteristics and locations of genomic regions that evade epigenetic reprogramming in sheep, as well as to explore the biological functions of the genes within these regions. In a previous study, we identified 107 transgenerationally inherited differentially methylated cytosines (DMCs) in the F1 and F2 generations in response to a paternal methionine-supplemented diet. These DMCs were found in TEs, non-repetitive regions, and imprinted and non-imprinted genes. Our findings suggest that genomic regions, rather than TEs and imprinted genes, have the propensity to escape reprogramming and serve as potential candidates for transgenerational epigenetic inheritance. Notably, 34 transgenerational methylated genes influenced by paternal nutrition escaped reprogramming, impacting growth, development, male fertility, cardiac disorders, and neurodevelopment. Intriguingly, among these genes, 21 have been associated with neural development and brain disorders, such as autism, schizophrenia, bipolar disease, and intellectual disability. This suggests a potential genetic overlap between brain and infertility disorders. Overall, our study supports the concept of transgenerational epigenetic inheritance of environmentally induced marks in mammals.

Proteomics and phosphoproteomics profiling in glutamatergic neurons and microglia in an iPSC model of Jansen de Vries Syndrome

Background

Jansen de Vries Syndrome (JdVS) is a rare neurodevelopmental disorder (NDD) caused by gain-of-function (GOF) truncating mutations in PPM1D exons 5 or 6. PPM1D is a serine/threonine phosphatase that plays an important role in the DNA damage response (DDR) by negatively regulating TP53 (P53). JdVS-associated mutations lead to the formation of a truncated PPM1D protein that retains catalytic activity and has a GOF effect because of reduced degradation. Somatic PPM1D exons 5 and 6 truncating mutations are well-established factors in a number of cancers, due to excessive dephosphorylation and reduced function of P53 and other substrates involved in DDR. Children with JdVS have a variety of neurodevelopmental, psychiatric, and physical problems. In addition, a small fraction has acute neuropsychiatric decompensation apparently triggered by infection or severe non-infectious environmental stress factors.

Methods

To understand the molecular basis of JdVS, we developed an induced pluripotent stem cell (iPSC) model system. iPSCs heterozygous for the truncating variant (PPM1D+/tr), were made from a patient, and control lines engineered using CRISPR-Cas9 gene editing. Proteomics and phosphoprotemics analyses were carried out on iPSC-derived glutamatergic neurons and microglia from three control and three PPM1D+/tr iPSC lines. We also analyzed the effect of the TLR4 agonist, lipopolysaccharide, to understand how activation of the innate immune system in microglia could account for acute behavioral decompensation.

Results

One of the major findings was the downregulation of POGZ in unstimulated microglia. Since loss-of-function variants in the POGZ gene are well-known causes of autism spectrum disorder, the decrease in PPM1D+/tr microglia suggests this plays a role in the neurodevelopmental aspects of JdVS. In addition, neurons, baseline, and LPS-stimulated microglia show marked alterations in the expression of several E3 ubiquitin ligases, most notably UBR4, and regulators of innate immunity, chromatin structure, ErbB signaling, and splicing. In addition, pathway analysis points to overlap with neurodegenerative disorders.

Limitations

Owing to the cost and labor-intensive nature of iPSC research, the sample size was small.

Conclusions

Our findings provide insight into the molecular basis of JdVS and can be extrapolated to understand neuropsychiatric decompensation that occurs in subgroups of patients with ASD and other NDDs.

Clinical profiling of MRD48 and functional characterization of two novel pathogenic RAC1 variants

RAC1 is a member of the Rac/Rho GTPase subfamily within the RAS superfamily of small GTP-binding proteins, comprising 3 paralogs playing a critical role in actin cytoskeleton remodeling, cell migration, proliferation and differentiation. De novo missense variants in RAC1 are associated with a rare neurodevelopmental disorder (MRD48) characterized by DD/ID and brain abnormalities coupled with a wide range of additional features. Structural and functional studies have documented either a dominant negative or constitutively active behavior for a subset of mutations. Here, we describe two individuals with previously unreported de novo missense RAC1 variants. We functionally demonstrate their pathogenicity proving a gain-of-function (GoF) effect for both. By reviewing the clinical features of these two individuals and the previously published MRD48 subjects, we further delineate the clinical profile of the disorder, confirming its phenotypic variability. Moreover, we compare the main features of MRD48 with the neurodevelopmental disease caused by GoF variants in the paralog RAC3, highlighting similarities and differences. Finally, we review all previously reported variants in RAC proteins and in the closely related CDC42, providing an updated overview of the spectrum and hotspots of pathogenic variants affecting these functionally related GTPases.

Influence of Auditory Integrative Training on Casein Kinase 2 and Its Impact on Behavioral and Social Interaction in Children with Autism Spectrum Disorder

Considerable disturbances in post-translational protein phosphorylation have recently been discovered in multiple neurological disorders. Casein kinase-2 (CK2) is a tetrameric Ser/Thr protein kinase that phosphorylates a large number of substrates and contributes in several cellular physiological and pathological processes. CK2 is highly expressed in the mammalian brain and catalyzes the phosphorylation of a large number of substrates that are crucial in neuronal or glial homeostasis and inflammatory signaling processes across synapses. In this study, we investigated the impact of auditory integration therapy (AIT) for the treatment of sensory processing abnormalities in autism on plasma CK2 levels. A total of 25 ASD children, aged between 5 and 12 years, were enrolled and participated in the present research study. AIT was performed for two weeks, for a period of 30 min, twice a day, with a 3 h interval between sessions. Before and after AIT, the Childhood Autism Rating Scale (CARS), Social Responsiveness Scale (SRS), and Short Sensory Profile (SSP) scores were calculated, and plasma CK2 levels were assayed using an ELISA test. The CARS and SRS indices of autism severity improved as a result of AIT, which could be related to the decreased level of plasma CK2. However, the mean value of the SSP scores was not significantly increased after AIT. The relationship between CK2 downregulation and glutamate excitotoxicity, neuro-inflammation, and leaky gut, as etiological mechanisms in ASD, was proposed and discussed. Further research, conducted on a larger scale and with a longer study duration, are required to assess whether the cognitive improvement in ASD children after AIT is related to the downregulation of CK2.

CRS: a circadian rhythm score model for predicting prognosis and treatment response in cancer patients

BACKGROUND

Circadian rhythm regulates complex physiological activities in organisms. A strong link between circadian dysfunction and cancer has been identified. However, the factors of dysregulation and functional significance of circadian rhythm genes in cancer have received little attention.

METHODS

In 18 cancer types from The Cancer Genome Atlas (TCGA), the differential expression and genetic variation of 48 circadian rhythm genes (CRGs) were examined. The circadian rhythm score (CRS) model was created using the ssGSEA method, and patients were divided into high and low groups based on the CRS. The Kaplan-Meier curve was created to assess the patient survival rate. Cibersort and estimate methods were used to identify the infiltration characteristics of immune cells between different CRS subgroups. Gene Expression Omnibus (GEO) dataset is used as verification queue and model stability evaluation queue. The CRS model's ability to predict chemotherapy and immunotherapy was assessed. Wilcoxon rank-sum test was used to compare the differences of CRS among different patients. We use CRS to identify potential "clock-drugs" by the connective map method.

RESULTS

Transcriptomic and genomic analyses of 48 CRGs revealed that most core clock genes are up-regulated, while clock control genes are down-regulated. Furthermore, we show that copy number variation may affect CRGs aberrations. Based on CRS, patients can be classified into two groups with significant differences in survival and immune cell infiltration. Further studies showed that patients with low CRS were more sensitive to chemotherapy and immunotherapy. Additionally, we identified 10 compounds (e.g. flubendazole, MLN-4924, ingenol) that are positively associated with CRS, and have the potential to modulate circadian rhythms.

CONCLUSIONS

CRS can be utilized as a clinical indicator to predict patient prognosis and responsiveness to therapy, and identify potential "clock-drugs".

Networks of placental DNA methylation correlate with maternal serum PCB concentrations and child neurodevelopment

BACKGROUND

Gestational exposure to polychlorinated biphenyls (PCBs) has been associated with elevated risk for neurodevelopmental disorders. Placental epigenetics may serve as a potential mechanism of risk or marker of altered placental function. Prior studies have associated differential placental DNA methylation with maternal PCB exposure or with increased risk of autism spectrum disorder (ASD). However, sequencing-based placental methylomes have not previously been tested for simultaneous associations with maternal PCB levels and child neurodevelopmental outcomes.

OBJECTIVES

We aimed to identify placental DNA methylation patterns associated with maternal PCB levels and child neurodevelopmental outcomes in the high-risk ASD MARBLES cohort.

METHODS

We measured 209 PCB congeners in 104 maternal serum samples collected at delivery. We identified networks of DNA methylation from 147 placenta samples using the Comethyl R package, which performs weighted gene correlation network analysis for whole genome bisulfite sequencing data. We tested placental DNA methylation modules for association with maternal serum PCB levels, child neurodevelopment, and other participant traits.

RESULTS

PCBs 153 + 168, 170, 180 + 193, and 187 were detected in over 50% of maternal serum samples and were highly correlated with one another. Consistent with previous findings, maternal age was the strongest predictor of serum PCB levels, alongside year of sample collection, pre-pregnancy BMI, and polyunsaturated fatty acid levels. Twenty seven modules of placental DNA methylation were identified, including five which significantly correlated with one or more PCBs, and four which correlated with child neurodevelopment. Two modules associated with maternal PCB levels as well as child neurodevelopment, and mapped to CSMD1 and AUTS2, genes previously implicated in ASD and identified as differentially methylated regions in mouse brain and placenta following gestational PCB exposure.

CONCLUSIONS

Placental DNA co-methylation modules were associated with maternal PCBs and child neurodevelopment. Methylation of CSMD1 and AUTS2 could be markers of altered placental function and/or ASD risk following maternal PCB exposure.

DNA methylation in people with anorexia nervosa: Epigenome-wide patterns in actively ill, long-term remitted, and healthy-eater women

OBJECTIVES

Recent studies have reported altered methylation levels at disorder-relevant DNA sites in people who are ill with Anorexia Nervosa (AN) compared to findings in people with no eating disorder (ED) or in whom AN has remitted. The preceding implies state-related influences upon gene expression in people with AN. This study further examined this notion.

METHODS

We measured genome-wide DNA methylation in 145 women with active AN, 49 showing stable one-year remission of AN, and 64 with no ED.

RESULTS

Comparisons revealed 205 differentially methylated sites between active and no ED groups, and 162 differentially methylated sites between active and remitted groups (Q < 0.01). Probes tended to map onto genes relevant to psychiatric, metabolic and immune functions. Notably, several of the genes identified here as being differentially methylated in people with AN (e.g. SYNJ2, PRKAG2, STAT3, CSGALNACT1, NEGR1, NR1H3) have figured in previous studies on AN. Effects also associated illness chronicity and lower BMI with more pronounced DNA methylation alterations, and remission of AN with normalisation of DNA methylation.

CONCLUSIONS

Findings corroborate earlier results suggesting reversible DNA methylation alterations in AN, and point to particular genes at which epigenetic mechanisms may act to shape AN phenomenology.

Atrophin-1 Function and Dysfunction in Dentatorubral-Pallidoluysian Atrophy

Dentatorubral-pallidoluysian atrophy (DRPLA) is a rare, incurable genetic disease that belongs to the group of polyglutamine (polyQ) diseases. DRPLA is the most common in the Japanese population; however, its global prevalence is also increasing due to better clinical recognition. It is characterized by cerebellar ataxia, myoclonus, epilepsy, dementia, and chorea. DRPLA is caused by dynamic mutation of CAG repeat expansion in ATN1 gene encoding the atrophin-1 protein. In the cascade of molecular disturbances, the pathological form of atrophin-1 is the initial factor, which has not been precisely characterized so far. Reports indicate that DRPLA is associated with disrupted protein-protein interactions (in which an expanded polyQ tract plays a crucial role), as well as gene expression deregulation. There is a great need to design efficient therapy that would address the underlying neurodegenerative process and thus prevent or alleviate DRPLA symptoms. An in-depth understanding of the normal atrophin-1 function and mutant atrophin-1 dysfunction is crucial for this purpose. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

The altered expression of cytoskeletal and synaptic remodeling proteins during epilepsy

The cytoskeleton plays an important role in epilepsy; however, the mechanism is unknown. Therefore, this study aimed to reveal the mechanism of cytoskeletal proteins in epilepsy by investigating the expression of cytoskeletal proteins and synaptophysin (SYP) in mice at 0, 3, 6, and 24 h, 3 days, and 7 days in a kainic acid (KA)-induced epileptic model. Our results demonstrated that the expression of F-actin decreased significantly between 3 and 6 h, 6 and 24 h, and 24 h and 3 days (P < 0.05). Meanwhile, the expression of the neurofilament light chain, neurofilament medium chain, and neurofilament heavy chain subunits was significantly decreased (P < 0.001) at 3 h after the KA injection compared to the KA 0 h group, followed by an elevation at 6 h and a further decrease at 24 h compared to at 6 h. SYP expression was significantly decreased between 0 and 3 h as well as between 3 and 6 h (P < 0.05). At 24 h, the level was increased compared to at 6 h and continued to increase at 3 days after the KA injection. Thus, we propose that cytoskeletal proteins may be involved in the pathogenesis of epilepsy.

Circulating cell-free DNA methylation mirrors alterations in cerebral patterns in epilepsy

BACKGROUND

DNA methylation profiling of circulating cell-free DNA (cfDNA) has rapidly become a promising strategy for biomarker identification and development. The cell-type-specific nature of DNA methylation patterns and the direct relationship between cfDNA and apoptosis can potentially be used non-invasively to predict local alterations. In addition, direct detection of altered DNA methylation patterns performs well as a biomarker. In a previous study, we demonstrated marked DNA methylation alterations in brain tissue from patients with mesial temporal lobe epilepsy with hippocampal sclerosis (MTLE-HS).

RESULTS

We performed DNA methylation profiling in cfDNA isolated from the serum of MTLE patients and healthy controls using BeadChip arrays followed by systematic bioinformatic analysis including deconvolution analysis and integration with DNase accessibility data sets. Differential cfDNA methylation analysis showed an overrepresentation of gene ontology terms and transcription factors related to central nervous system function and regulation. Deconvolution analysis of the DNA methylation data sets ruled out the possibility that the observed differences were due to changes in the proportional contribution of cortical neurons in cfDNA. Moreover, we found no overrepresentation of neuron- or glia-specific patterns in the described cfDNA methylation patterns. However, the MTLE-HS cfDNA methylation patterns featured a significant overrepresentation of the epileptic DNA methylation alterations previously observed in the hippocampus.

CONCLUSIONS

Our results support the use of cfDNA methylation profiling as a rational approach to seeking non-invasive and reproducible epilepsy biomarkers.

Pathophysiological Mechanism of Neurodevelopmental Disorders-Overview

Technological advancements in next-generation DNA sequencing have enabled elucidation of many genetic causes of neurodevelopmental disorders (NDDs) over the last two decades [...].

Disruption of the gene regulatory programme in neurodevelopmental disorders

Cortical development consists of a series of synchronised events, including fate transition of cortical progenitors, neuronal migration, specification and connectivity. It is becoming clear that gene expression programs governing these events rely on the interplay between signalling molecules, transcription factors and epigenetic mechanisms. When genetic or environmental factors disrupt expression of genes involved in important brain development processes, neurodevelopmental disorders can occur. This review aims to highlight how recent advances in technologies have helped uncover and imitate the gene regulatory mechanisms commonly disrupted in neurodevelopmental disorders.

auts2 Features and Expression Are Highly Conserved during Evolution Despite Different Evolutionary Fates Following Whole Genome Duplication

The AUTS2 gene plays major roles during brain development and is associated with various neuropathologies including autism. Data in non-mammalian species are scarce, and the aim of our study was to provide a comprehensive analysis of auts2 evolution in teleost fish, which are widely used for in vivo functional analysis and biomedical purposes. Comparative genomics in 78 species showed that auts2a and auts2b originate from the teleost-specific whole genome duplication (TGD). auts2a, which is highly similar to human AUTS2, was almost systematically retained following TGD. In contrast, auts2b, which encodes for a shorter protein similar to a short human AUTS2 isoform, was lost more frequently and independently during evolution. RNA-seq analysis in 10 species revealed a highly conserved profile with predominant expression of both genes in the embryo, brain, and gonads. Based on protein length, conserved domains, and expression profiles, we speculate that the long human isoform functions were retained by auts2a, while the short isoform functions were retained by auts2a and/or auts2b, depending on the lineage/species. auts2a showed a burst in expression during medaka brain formation, where it was expressed in areas of the brain associated with neurodevelopmental disorders. Together, our data suggest a strong conservation of auts2 functions in vertebrates despite different evolutionary scenarios in teleosts.

[Whole-transcriptome sequencing analysis of placental differential miRNA expression profile in Down syndrome]

OBJECTIVE

To identify new biomarkers and molecular pathogenesis of Down syndrome (DS) by analyzing differentially expressed miRNAs in the placentas and their biological pathways.

METHODS

Whole transcriptome sequencing was used to identify the differentially expressed miRNAs in DS (n=3) and normal placental samples (n=3) diagnosed by prenatal diagnosis. The target genes were predicted using miRWalk, Targetscan and miRDB, and GO and KEGG pathway analyses were performed for gene enrichment studies.

RESULTS

We identified a total of 82 differentially expressed miRNAs in the placental tissues of DS, including 29 up-regulated miRNAs (fold change ≥2, P < 0.05) and 15 down-regulated miRNAs (fold change ≥2, P < 0.05), among which 10 miRNAs with relatively high expression abundance were selected for further analysis, including 4 up-regulated and 6 down-regulated miRNAs. These selected miRNAs shared the common target genes BTBD3 and AUTS2, both of which were associated with neurodevelopment. GO analysis showed that the target genes of the selected miRNAs were mainly enriched in protein binding, hydrolytic enzymes, metal ion binding protein combining, transferase activity, nucleotide, cytoplasmic constituents, nucleus composition, transcriptional regulation, RNA metabolism regulation, DNA-dependent RNA polymerase Ⅱ promoter transcriptional regulation, eye development, and sensory organ development. KEGG enrichment analysis showed that the target genes of these differentially expressed miRNAs were involved in the signaling pathways including tumor-related signaling pathway, PI3K-Akt signaling pathway, Ras signaling pathway, Rap1 signaling pathway, cytoskeletal regulatory signaling pathway, purine metabolization-related signaling pathway and P53 signaling pathway.

CONCLUSION

The differentially expressed miRNAs may play important roles in placental damage and pregnancy pathology in DS and provide clues for the prevention and treatment of mental retardation-related diseases.

Thrombosis-related circulating miR-16-5p is associated with disease severity in patients hospitalised for COVID-19

SARS-CoV-2 tropism for the ACE2 receptor, along with the multifaceted inflammatory reaction, is likely to drive the generalized hypercoagulable and thrombotic state seen in patients with COVID-19. Using the original bioinformatic workflow and network medicine approaches we reanalysed four coronavirus-related expression datasets and performed co-expression analysis focused on thrombosis and ACE2 related genes. We identified microRNAs (miRNAs) which play role in ACE2-related thrombosis in coronavirus infection and further, we validated the expressions of precisely selected miRNAs-related to thrombosis (miR-16-5p, miR-27a-3p, let-7b-5p and miR-155-5p) in 79 hospitalized COVID-19 patients and 32 healthy volunteers by qRT-PCR. Consequently, we aimed to unravel whether bioinformatic prioritization could guide selection of miRNAs with a potential of diagnostic and prognostic biomarkers associated with disease severity in patients hospitalized for COVID-19. In bioinformatic analysis, we identified EGFR, HSP90AA1, APP, TP53, PTEN, UBC, FN1, ELAVL1 and CALM1 as regulatory genes which could play a pivotal role in COVID-19 related thrombosis. We also found miR-16-5p, miR-27a-3p, let-7b-5p and miR-155-5p as regulators in the coagulation and thrombosis process. In silico predictions were further confirmed in patients hospitalized for COVID-19. The expression levels of miR-16-5p and let-7b in COVID-19 patients were lower at baseline, 7-days and 21-day after admission compared to the healthy controls (p < 0.0001 for all time points for both miRNAs). The expression levels of miR-27a-3p and miR-155-5p in COVID-19 patients were higher at day 21 compared to the healthy controls (p = 0.007 and p < 0.001, respectively). A low baseline miR-16-5p expression presents predictive utility in assessment of the hospital length of stay or death in follow-up as a composite endpoint (AUC:0.810, 95% CI, 0.71-0.91, p < 0.0001) and low baseline expression of miR-16-5p and diabetes mellitus are independent predictors of increased length of stay or death according to a multivariate analysis (OR: 9.417; 95% CI, 2.647-33.506; p = 0.0005 and OR: 6.257; 95% CI, 1.049-37.316; p = 0.044, respectively). This study enabled us to better characterize changes in gene expression and signalling pathways related to hypercoagulable and thrombotic conditions in COVID-19. In this study we identified and validated miRNAs which could serve as novel, thrombosis-related predictive biomarkers of the COVID-19 complications, and can be used for early stratification of patients and prediction of severity of infection development in an individual.Abbreviations: ACE2, angiotensin-converting enzyme 2AF, atrial fibrillationAPP, Amyloid Beta Precursor ProteinaPTT, activated partial thromboplastin timeAUC, Area under the curveAβ, amyloid betaBMI, body mass indexCAD, coronary artery diseaseCALM1, Calmodulin 1 geneCaM, calmodulinCCND1, Cyclin D1CI, confidence intervalCOPD, chronic obstructive pulmonary diseaseCOVID-19, Coronavirus disease 2019CRP, C-reactive proteinCV, CardiovascularCVDs, cardiovascular diseasesDE, differentially expressedDM, diabetes mellitusEGFR, Epithelial growth factor receptorELAVL1, ELAV Like RNA Binding Protein 1FLNA, Filamin AFN1, Fibronectin 1GEO, Gene Expression OmnibushiPSC-CMs, Human induced pluripotent stem cell-derived cardiomyocytesHSP90AA1, Heat Shock Protein 90 Alpha Family Class A Member 1Hsp90α, heat shock protein 90αICU, intensive care unitIL, interleukinIQR, interquartile rangelncRNAs, long non-coding RNAsMI, myocardial infarctionMiRNA, MiR, microRNAmRNA, messenger RNAncRNA, non-coding RNANERI, network-medicine based integrative approachNF-kB, nuclear factor kappa-light-chain-enhancer of activated B cellsNPV, negative predictive valueNXF, nuclear export factorPBMCs, Peripheral blood mononuclear cellsPCT, procalcitoninPPI, Protein-protein interactionsPPV, positive predictive valuePTEN, phosphatase and tensin homologqPCR, quantitative polymerase chain reactionROC, receiver operating characteristicSARS-CoV-2, severe acute respiratory syndrome coronavirus 2SD, standard deviationTLR4, Toll-like receptor 4TM, thrombomodulinTP53, Tumour protein P53UBC, Ubiquitin CWBC, white blood cells.