In search of SARS CoV-2 replication inhibitors: Virtual screening, molecular dynamics simulations and ADMET analysis.

Severe acute respiratory syndrome has relapsed recently as novel coronavirus causing a life threat to the entire world in the absence of an effective therapy. To hamper the replication of the deadly SARS CoV-2 inside the host cells, systematic in silico virtual screening of total 267,324 ligands from Asinex EliteSynergy and BioDesign libraries has been performed using AutoDock Vina against RdRp. The molecular modeling studies revealed the identification of twenty-one macrocyclic hits (2-22) with better binding energy than remdesivir (1), marketed SARS CoV-2 inhibitor. Further, the analysis using rules for drug-likeness and their ADMET profile revealed the candidature of these hits due to superior oral bioavailability and druggability. Further, the MD simulation studies of top two hits (2 and 3) performed using GROMACS 2020.1 for 10 ns revealed their stability into the docked complexes. These results provide an important breakthrough in the design of macrocyclic hits as SARS CoV-2 RNA replicase inhibitor.

Good or not good: Role of miR-18a in cancer biology.

miR-18a is a member of primary transcript called miR-17-92a (C13orf25 or MIR17HG) which also contains five other miRNAs: miR-17, miR-19a, miR-20a, miR-19b and miR-92a. This cluster as a whole shows specific characteristics, where miR-18a seems to be unique. In contrast to the other members, the expression of miR-18a is additionally controlled and probably functions as its own internal controller of the cluster. miR-18a regulates many genes involved in proliferation, cell cycle, apoptosis, response to different kinds of stress, autophagy and differentiation. The disturbances of miR-18a expression are observed in cancer as well as in different diseases or pathological states. The miR-17-92a cluster is commonly described as oncogenic and it is known as 'oncomiR-1', but this statement is a simplification because miR-18a can act both as an oncogene and a suppressor. In this review we summarize the current knowledge about miR-18a focusing on its regulation, role in cancer biology and utility as a potential biomarker.

ADAR1 and MicroRNA; A Hidden Crosstalk in Cancer.

The evolution of cancer cells is believed to be dependent on genetic or epigenetic alterations. However, this concept has recently been challenged by another mode of nucleotide alteration, RNA editing, which is frequently up-regulated in cancer. RNA editing is a biochemical process in which either Adenosine or Cytosine is deaminated by a group of RNA editing enzymes including ADAR (Adenosine deaminase; RNA specific) or APOBEC3B (Apolipoprotein B mRNA Editing Enzyme Catalytic Subunit 3B). The result of RNA editing is usually adenosine to inosine (A-to-I) or cytidine to uridine (C-to-U) transition, which can affect protein coding, RNA stability, splicing and microRNA-target interactions. The functional impact of these alterations is largely unclear and is a subject of extensive research. In the present review, we will specifically focus on the influence of ADARs on carcinogenesis via the regulation of microRNA processing and functioning. This follows a brief review of the current knowledge of properties of ADAR enzyme, RNA editing, and microRNA processing.

Noncoding regions of C. elegans mRNA undergo selective adenosine to inosine deamination and contain a small number of editing sites per transcript.

ADARs (Adenosine deaminases that act on RNA) "edit" RNA by converting adenosines to inosines within double-stranded regions. The primary targets of ADARs are long duplexes present within noncoding regions of mRNAs, such as introns and 3' untranslated regions (UTRs). Because adenosine and inosine have different base-pairing properties, editing within these regions can alter splicing and recognition by small RNAs. However, despite numerous studies identifying multiple editing sites in these genomic regions, little is known about the extent to which editing sites co-occur on individual transcripts or the functional output of these combinatorial editing events. To begin to address these questions, we performed an ultra-deep sequencing analysis of 4 Caenorhabditis elegans 3' UTRs that are known ADAR targets. Synchronous editing events were determined for the long duplexes in vivo. Furthermore, the validity of each editing event was confirmed by sequencing the same regions of mRNA from worms that lack A-to-I editing. This analysis identified a large number of editing sites that can occur within each 3' UTR, but interestingly, each individual transcript contained only a small fraction of these A-to-I editing events. In addition, editing patterns were not random, indicating that an editing event can affect the efficiency of editing at subsequent adenosines. Furthermore, we identified specific sites that can be both positively and negatively correlated with additional sites leading to mutually exclusive editing patterns. These results suggest that editing in noncoding regions is selective and hyper-editing of cellular RNAs is rare.

HIV-1 RRE RNA acts as an RNA silencing suppressor by competing with TRBP-bound siRNAs.

Several proteins and RNAs expressed by mammalian viruses have been reported to interfere with RNA interference (RNAi) activity. We investigated the ability of the HIV-1-encoded RNA elements Trans-Activation Response (TAR) and Rev-Response Element (RRE) to alter RNAi. MicroRNA let7-based assays showed that RRE is a potent suppressor of RNAi activity, while TAR displayed moderate RNAi suppression. We demonstrate that RRE binds to TAR-RNA Binding Protein (TRBP), an essential component of the RNA Induced Silencing Complex (RISC). The binding of TAR and RRE to TRBP displaces small interfering (si)RNAs from binding to TRBP. Several stem-deleted RRE mutants lost their ability to suppress RNAi activity, which correlated with a reduced ability to compete with siRNA-TRBP binding. A lentiviral vector expressing TAR and RRE restricted RNAi, but RNAi was restored when Rev or GagPol were coexpressed. Adenoviruses are restricted by RNAi and encode their own suppressors of RNAi, the Virus-Associated (VA) RNA elements. RRE enhanced the replication of wild-type and VA-deficient adenovirus. Our work describes RRE as a novel suppressor of RNAi that acts by competing with siRNAs rather than by disrupting the RISC. This function is masked in lentiviral vectors co-expressed with viral proteins and thus will not affect their use in gene therapy. The potent RNAi suppressive effects of RRE identified in this study could be used to enhance the expression of RNAi restricted viruses used in oncolysis such as adenoviruses.

Translation attenuation via 3' terminal codon usage in bovine csn1s2 is responsible for the difference in αs2- and ß-casein profile in milk.

The rate of secretion of αs2-casein into bovine milk is approximately 25% of that of ß-casein, yet mammary expression of their respective mRNA transcripts (csn1s2 and csn2) is not different. Our objective was to identify molecular mechanisms that explain the difference in translation efficiency between csn1s2 and csn2. Cell-free translational efficiency of csn2 was 5 times that of csn1s2. Transcripts of csn1s2 distributed into heavier polysomes than csn2 transcripts, indicating an attenuation of elongation and/or termination. Stimulatory and inhibitory effects of the 5' and 3' UTRs on translational efficiency were different with luciferase and casein sequences in the coding regions. Substituting the 5' and 3' UTRs from csn2 into csn1s2 did not improve csn1s2 translation, implicating the coding region itself in the translation difference. Deletion of a 28-codon fragment from the 3' terminus of the csn1s2 coding region, which displays codons with low correlations to cell fitness, increased translation to a par with csn2. We conclude that the usage of the last 28 codons of csn1s2 is the main regulatory element that attenuates its expression and is responsible for the differential translational expression of csn1s2 and csn2.

Identification of miRNAs that modulate glucocerebrosidase activity in Gaucher disease cells.

Gaucher disease is an autosomal recessive disorder caused by deficiency of the enzyme glucocerebrosidase. Although it is a monogenic disease, there is vast phenotypic heterogeneity, even among patients with the same genotype. MicroRNAs (miRNAs) are small non-coding RNAs involved in many biological processes and diseases. To determine whether miRNAs can affect glucocerebrosidase activity, we performed a screen of 875 different miRNA mimics. The screen was performed using Gaucher fibroblasts, and glucocerebrosidase activity was used as the initial outcome parameter. We found several miRNAs that either up- or down-regulated glucocerebrosidase activity. In follow-up assays, we confirmed that one specific miRNA (miR-127-5p) down-regulated both glucocerebrosidase activity and protein levels by down-regulation of LIMP-2, the receptor involved in proper trafficking of glucocerebrosidase from the endoplasmic reticulum to the lysosome. A conditioned media assay demonstrated that cells treated with this miRNA secreted glucocerebrosidase into the extracellular environment, supporting impaired LIMP-2 function. Two other miRNAs, miR-16-5p and miR-195-5p, were found to up-regulate glucocerebrosidase activity by greater than 40% and to enhance expression and protein levels of the enzyme. In conclusion, we show that miRNAs can alter glucocerebrosidase activity in patient cells, indicating that miRNAs can potentially act as modifiers in Gaucher disease.

Apically-located P4-ATPase1-Lem1 complex internalizes phosphatidylserine and regulates motility-dependent invasion and egress in Toxoplasma gondii.

The membrane asymmetry regulated by P4-ATPases is crucial for the functioning of eukaryotic cells. The underlying spatial translocation or flipping of specific lipids is usually assured by respective P4-ATPases coupled to conforming non-catalytic subunits. Our previous work has identified five P4-ATPases (TgP4-ATPase1-5) and three non-catalytic partner proteins (TgLem1-3) in the intracellular protozoan pathogen, Toxoplasma gondii. However, their flipping activity, physiological relevance and functional coupling remain unknown. Herein, we demonstrate that TgP4-ATPase1 and TgLem1 work together to translocate phosphatidylserine (PtdSer) during the lytic cycle of T. gondii. Both proteins localize in the plasma membrane at the invasive (apical) end of its acutely-infectious tachyzoite stage. The genetic knockout of P4-ATPase1 and conditional depletion of Lem1 in tachyzoites severely disrupt the asexual reproduction and translocation of PtdSer across the plasma membrane. Moreover, the phenotypic analysis of individual mutants revealed a requirement of lipid flipping for the motility, egress and invasion of tachyzoites. Not least, the proximity-dependent biotinylation and reciprocal immunoprecipitation assays demonstrated the physical interaction of P4-ATPase1 and Lem1. Our findings disclose the mechanism and significance of PtdSer flipping during the lytic cycle and identify the P4-ATPase1-Lem1 heterocomplex as a potential drug target in T. gondii.

MiR-590-3p affects the function of adipose-derived stem cells (ADSCs) on the survival of skin flaps by targeting VEGFA.

Introduction: Partial necrosis of skin flaps is still a substantial problem in plastic and reconstructive surgery. In this study, the role of miR-590-3p in adipose-derived stem cells (ADSCs) transplantation in improving the survival of skin flap in a mouse model was delved into. Method: An abdominal perforator flap model was established in mice. The histopathological examination of mice skin tissues after ADSCs transplantation was implemented using Hematoxylin & eosin (H&E) staining. Immunohistochemistry (IHC) or immunofluorescence (IF) staining was utilized to assess the PCNA or CD31 levels. The concentrations of VEGFA in the culture medium were quantified using a VEGFA ELISA kit. Result: The damage of tissue in the skin flap was dramatically relieved by ADSCs transplantation. MiR-590-3p overexpression notably suppressed, while miR-590-3p knockdown facilitated skin flap survival by regulating PCNA, VCAM-1, and VEGFA levels. MiR-590-3p targeted VEGFA to regulate its expression. The knockdown of VEGFA significantly inhibited, while overexpression of VEGFA notably promoted the survival of skin flap. Conclusion: ADSCs transplantation promotes skin flap survival by boosting angiogenesis. The miR-590-3p/VEGFA axis modulates skin flap angiogenesis and survival in ADSCs. These results reveal that interfering with miR-590-3p in ADSCs could potentially be a novel therapeutic target for the improvement of skin flap survival.

Mutational cascade of SARS-CoV-2 leading to evolution and emergence of omicron variant.

BACKGROUND: Emergence of new variant of SARS-CoV-2, namely omicron, has posed a global concern because of its high rate of transmissibility and mutations in its genome. Researchers worldwide are trying to understand the evolution and emergence of such variants to understand the mutational cascade events. METHODS: We have considered all omicron genomes (n = 302 genomes) available till 2nd December 2021 in the public repository of GISAID along with representatives of variants of concern (VOC), i.e., alpha, beta, gamma, delta, and omicron; variant of interest (VOI) mu and lambda; and variant under monitoring (VUM). Whole genome-based phylogeny and mutational analysis were performed to understand the evolution of SARS CoV-2 leading to emergence of omicron variant. RESULTS: Whole genome-based phylogeny depicted two phylogroups (PG-I and PG-II) forming variant specific clades except for gamma and VUM GH. Mutational analysis detected 18,261 mutations in the omicron variant, majority of which were non-synonymous mutations in spike (A67, T547K, D614G, H655Y, N679K, P681H, D796Y, N856K, Q954H), followed by RNA dependent RNA polymerase (rdrp) (A1892T, I189V, P314L, K38R, T492I, V57V), ORF6 (M19M) and nucleocapsid protein (RG203KR). CONCLUSION: Delta and omicron have evolutionary diverged into distinct phylogroups and do not share a common ancestry. While, omicron shares common ancestry with VOI lambda and its evolution is mainly derived by the non-synonymous mutations.

Murine model of triosephosphate isomerase deficiency with anemia and severe neuromuscular dysfunction.

Triosephosphate isomerase deficiency (TPI Df) is a rare, aggressive genetic disease that typically affects young children and currently has no established treatment. TPI Df is characterized by hemolytic anemia, progressive neuromuscular degeneration, and a markedly reduced lifespan. The disease has predominately been studied using invertebrate and in vitro models, which lack key aspects of the human disease. While other groups have generated mammalian Tpi1 mutant strains, specifically with the mouse mus musculus, these do not recapitulate key characteristic phenotypes of the human disease. Reported here is the generation of a novel murine model of TPI Df. CRISPR-Cas9 was utilized to engineer the most common human disease-causing mutation, Tpi1 E105D , and Tpi1 null mice were also isolated as a frame-shifting deletion. Tpi1 E105D/null mice experience a markedly shortened lifespan, postural abnormalities consistent with extensive neuromuscular dysfunction, hemolytic anemia, pathological changes in spleen, and decreased body weight. There is a ∼95% reduction in TPI protein levels in Tpi1 E105D/null animals compared to wild-type littermates, consistent with decreased TPI protein stability, a known cause of TPI Df. This work illustrates the capability of Tpi1 E105D/null mice to serve as a mammalian model of human TPI Df. This work will allow for advancement in the study of TPI Df within a model with physiology similar to humans. The development of the model reported here will enable mechanistic studies of disease pathogenesis and, importantly, efficacy testing in a mammalian system for emerging TPI Df treatments.

miRNA expression in COVID-19.

Coronavirus disease 2019 (COVID-19) is regarded as a challenge in health system. Several studies have assessed the immune-related aspect of this disorder to identify the host-related factors that affect the course of COVID-19. microRNAs (miRNAs) as potent regulators of immune responses have gained much attention in this regard. Recent studies have shown aberrant expression of miRNAs in COVID-19 in association with disease course. Differentially expressed miRNAs have been enriched in pathways related with inflammation and antiviral immune response. miRNAs have also been regarded as potential therapeutic targets in COVID-19, particularly for management of pathological consequences of COVID-19. In the current review, we summarize the data about dysregulation of miRNAs in COVID-19.

In vitro and in vivo characterization of erythrosin B and derivatives against Zika virus.

Zika virus (ZIKV) causes significant human diseases without specific therapy. Previously we found erythrosin B, an FDA-approved food additive, inhibited viral NS2B-NS3 interactions, leading to inhibition of ZIKV infection in cell culture. In this study, we performed pharmacokinetic and in vivo studies to demonstrate the efficacy of erythrosin B against ZIKV in 3D mini-brain organoid and mouse models. Our results showed that erythrosin B is very effective in abolishing ZIKV replication in the 3D organoid model. Although pharmacokinetics studies indicated that erythrosin B had a low absorption profile, mice challenged by a lethal dose of ZIKV showed a significantly improved survival rate upon oral administration of erythrosin B, compared to vehicle control. Limited structure-activity relationship studies indicated that most analogs of erythrosin B with modifications on the xanthene ring led to loss or reduction of inhibitory activities towards viral NS2B-NS3 interactions, protease activity and antiviral efficacy. In contrast, introducing chlorine substitutions on the isobenzofuran ring led to slightly increased activities, suggesting that the isobenzofuran ring is well tolerated for modifications. Cytotoxicity studies indicated that all derivatives are nontoxic to human cells. Overall, our studies demonstrated erythrosin B is an effective antiviral against ZIKV both in vitro and in vivo.

The dynamic nature of the coronavirus receptor, angiotensin-converting enzyme 2 (ACE2) in differentiating airway epithelia.

Once inhaled, SARS-CoV-2 particles enter respiratory ciliated cells by interacting with angiotensin converting enzyme 2 (ACE2). Understanding the nature of ACE2 within airway tissue has become a recent focus particularly in light of the COVID-19 pandemic. Airway mucociliary tissue was generated in-vitro using primary human nasal epithelial cells and the air-liquid interface (ALI) model of differentiation. Using ALI tissue, three distinct transcript variants of ACE2 were identified. One transcript encodes the documented full-length ACE2 protein. The other two transcripts are unique truncated isoforms, that until recently had only been predicted to exist via sequence analysis software. Quantitative PCR revealed that all three transcript variants are expressed throughout differentiation of airway mucociliary epithelia. Immunofluorescence analysis of individual ACE2 protein isoforms exogenously expressed in cell-lines revealed similar abilities to localize in the plasma membrane and interact with the SARS CoV 2 spike receptor binding domain. Immunohistochemistry on differentiated ALI tissue using antibodies to either the N-term or C-term of ACE2 revealed both overlapping and distinct signals in cells, most notably only the ACE2 C-term antibody displayed plasma-membrane localization. We also demonstrate that ACE2 protein shedding is different in ALI Tissue compared to ACE2-transfected cell lines, and that ACE2 is released from both the apical and basal surfaces of ALI tissue. Together, our data highlights various facets of ACE2 transcripts and protein in airway mucociliary tissue that may represent variables which impact an individual's susceptibility to SARS-CoV-2 infection, or the severity of Covid-19.

RNAi-mediated siRNA sequences to combat the COVID-19 pandemic with the inhibition of SARS-CoV2.

The outbreak of the COVID-19 pandemic has cost five million lives to date, and was caused by a positive-sense RNA virus named SARS-CoV2. The lack of drugs specific to SARS-CoV2, leads us to search for an effective and specific therapeutic approach. Small interfering RNA (siRNA) is able to activate the RNA interference (RNAi) pathway to silence the specific targeted gene and inhibit the viral replication, and it has not yet attracted enough attention as a SARS-CoV2 antiviral agent. It could be a potential weapon to combat this pandemic until the completion of full scale, effective mass vaccination. For this study, specific siRNAs were designed using a web-based bioinformatics tool (siDirect2.0) against 14 target sequences. These might have a high probability of silencing the essential proteins of SARS-CoV2. such as: 3CLpro/Mpro/nsp5, nsp7, Rd-Rp/nsp12, ZD, NTPase/HEL or nsp13, PLpro/nsp3, envelope protein (E), spike glycoprotein (S), nucleocapsid phosphoprotein (N), membrane glycoprotein (M), ORF8, ORF3a, nsp2, and its respective 5' and 3'-UTR. Among these potential drug targets, the majority of them contain highly conserved sequences; the rest are chosen on the basis of their role in viral replication and survival. The traditional vaccine development technology using SARS-CoV2 protein takes 6-8 months; meanwhile the virus undergoes several mutations in the candidate protein chosen for vaccine development. By the time the protein-based vaccine reaches the market, the virus would have undergone several mutations, such that the antibodies against the viral sequence may not be effective in restricting the newly mutated viruses. However, siRNA technology can make sequences based on real time viral mutation status. This has the potential for suppressing SARS-CoV2 viral replication, through RNAi technology.

Emergence, evolution, and vaccine production approaches of SARS-CoV-2 virus: Benefits of getting vaccinated and common questions.

The emergence of coronavirus disease 2019 (COVID-19) pandemic in Wuhan city, China at the end of 2019 made it urgent to identify the origin of the causal pathogen and its molecular evolution, to appropriately design an effective vaccine. This study analyzes the evolutionary background of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2 or SARS-2) in accordance with its close relative SARS-CoV (SARS-1), which was emerged in 2002. A comparative genomic and proteomic study was conducted on SARS-2, SARS-1, and Middle East respiratory syndrome coronavirus (MERS), which was emerged in 2012. In silico analysis inferred the genetic variability among the tested viruses. The SARS-1 genome harbored 11 genes encoding 12 proteins, while SARS-2 genome contained only 10 genes encoding for 10 proteins. MERS genome contained 11 genes encoding 11 proteins. The analysis also revealed a slight variation in the whole genome size of SARS-2 comparing to its siblings resulting from sequential insertions and deletions (indels) throughout the viral genome particularly ORF1AB, spike, ORF10 and ORF8. The effective indels were observed in the gene encoding the spike protein that is responsible for viral attachment to the angiotensin-converting enzyme 2 (ACE2) cell receptor and initiating infection. These indels are responsible for the newly emerging COVID-19 variants αCoV, ßCoV, γCoV and δCoV. Nowadays, few effective COVID-19 vaccines developed based on spike (S) glycoprotein were approved and become available worldwide. Currently available vaccines can relatively prevent the spread of COVID-19 and suppress the disease. The traditional (killed or attenuated virus vaccine and antibody-based vaccine) and innovated vaccine production technologies (RNA- and DNA-based vaccines and viral vectors) are summarized in this review. We finally highlight the most common questions related to COVID-19 disease and the benefits of getting vaccinated.

Novel role of macrophage TXNIP-mediated CYLD-NRF2-OASL1 axis in stress-induced liver inflammation and cell death.

Background & Aims: The stimulator of interferon genes (STING)/TANK-binding kinase 1 (TBK1) pathway is vital in mediating innate immune and inflammatory responses during oxidative/endoplasmic reticulum (ER) stress. However, it remains unknown whether macrophage thioredoxin-interacting protein (TXNIP) may regulate TBK1 function and cell death pathways during oxidative/ER stress. Methods: A mouse model of hepatic ischaemia/reperfusion injury (IRI), the primary hepatocytes, and bone marrow-derived macrophages were used in the myeloid-specific TXNIP knockout (TXNIPM-KO) and TXNIP-proficient (TXNIPFL/FL) mice. Results: The TXNIPM-KO mice were resistant to ischaemia/reperfusion (IR) stress-induced liver damage with reduced serum alanine aminotransferase (ALT)/aspartate aminotransferase (AST) levels, macrophage/neutrophil infiltration, and pro-inflammatory mediators compared with the TXNIPFL/FL controls. IR stress increased TXNIP, p-STING, and p-TBK1 expression in ischaemic livers. However, TXNIPM-KO inhibited STING, TBK1, interferon regulatory factor 3 (IRF3), and NF-κB activation with interferon-ß (IFN-ß) expression. Interestingly, TXNIPM-KO augmented nuclear factor (erythroid-derived 2)-like 2 (NRF2) activity, increased antioxidant gene expression, and reduced macrophage reactive oxygen species (ROS) production and hepatic apoptosis/necroptosis in IR-stressed livers. Mechanistically, macrophage TXNIP deficiency promoted cylindromatosis (CYLD), which colocalised and interacted with NADPH oxidase 4 (NOX4) to enhance NRF2 activity by deubiquitinating NOX4. Disruption of macrophage NRF2 or its target gene 2',5' oligoadenylate synthetase-like 1 (OASL1) enhanced Ras GTPase-activating protein-binding protein 1 (G3BP1) and TBK1-mediated inflammatory response. Notably, macrophage OASL1 deficiency induced hepatocyte apoptotic peptidase activating factor 1 (APAF1), cytochrome c, and caspase-9 activation, leading to increased caspase-3-initiated apoptosis and receptor-interacting serine/threonine-protein kinase 3 (RIPK3)-mediated necroptosis. Conclusions: Macrophage TXNIP deficiency enhances CYLD activity and activates the NRF2-OASL1 signalling, controlling IR stress-induced liver injury. The target gene OASL1 regulated by NRF2 is crucial for modulating STING-mediated TBK1 activation and Apaf1/cytochrome c/caspase-9-triggered apoptotic/necroptotic cell death pathway. Our findings underscore a novel role of macrophage TXNIP-mediated CYLD-NRF2-OASL1 axis in stress-induced liver inflammation and cell death, implying the potential therapeutic targets in liver inflammatory diseases. Lay summary: Liver inflammation and injury induced by ischaemia and reperfusion (the absence of blood flow to the liver tissue followed by the resupply of blood) is a significant cause of hepatic dysfunction and failure following liver transplantation, resection, and haemorrhagic shock. Herein, we uncover an underlying mechanism that contributes to liver inflammation and cell death in this setting and could be a therapeutic target in stress-induced liver inflammatory injury.

Integrative RNA profiling of TBEV-infected neurons and astrocytes reveals potential pathogenic effectors.

Tick-borne encephalitis virus (TBEV), the most medically relevant tick-transmitted flavivirus in Eurasia, targets the host central nervous system and frequently causes severe encephalitis. The severity of TBEV-induced neuropathogenesis is highly cell-type specific and the exact mechanism responsible for such differences has not been fully described yet. Thus, we performed a comprehensive analysis of alterations in host poly-(A)/miRNA/lncRNA expression upon TBEV infection in vitro in human primary neurons (high cytopathic effect) and astrocytes (low cytopathic effect). Infection with severe but not mild TBEV strain resulted in a high neuronal death rate. In comparison, infection with either of TBEV strains in human astrocytes did not. Differential expression and splicing analyses with an in silico prediction of miRNA/mRNA/lncRNA/vd-sRNA networks found significant changes in inflammatory and immune response pathways, nervous system development and regulation of mitosis in TBEV Hypr-infected neurons. Candidate mechanisms responsible for the aforementioned phenomena include specific regulation of host mRNA levels via differentially expressed miRNAs/lncRNAs or vd-sRNAs mimicking endogenous miRNAs and virus-driven modulation of host pre-mRNA splicing. We suggest that these factors are responsible for the observed differences in the virulence manifestation of both TBEV strains in different cell lines. This work brings the first complex overview of alterations in the transcriptome of human astrocytes and neurons during the infection by two TBEV strains of different virulence. The resulting data could serve as a starting point for further studies dealing with the mechanism of TBEV-host interactions and the related processes of TBEV pathogenesis.

Characterization of Single Nucleotide Variants of OPN3 Gene in Melanocytic Nevi and Melanoma.

In this study, we examined single nucleotide variants (SNVs) of the OPN3 gene in malignant melanoma and melanocytic nevi. A total of 20 variants of SNVs were detected. Of these variants, five nonsynonymous mutations of OPN3 were identified, including c.T152C, c.T401C, c.G547A, c.G768A, and c.G992A. Three prediction tools, MutationTaster2, Polymorphism Phenotyping version 2, and PROVEAN (Protein Variation Effect Analyzer), which predict possible impact of an amino acid substitution, suggested that the mutations could be deleterious. Nine SNVs occurred in 3' untranslated regions, whereas two were observed in 5' untranslated regions. In all cases, four intronic variants were identified. In addition, we identified nine 3' untranslated region SNVs in OPN3; one of them (OPN3[NM_014322:c.∗83C>T]) is predicted to disrupt a conserved microRNA (has-miR-376c-3p) target site, located in position 86-93 of OPN3 3' untranslated region. Our findings suggest that there is a strong possibility that OPN3 SNVs play a role in the pathogenesis of melanocytic tumors via prediction of functional phenotype.

MLb-LDLr: A Machine Learning Model for Predicting the Pathogenicity of LDLr Missense Variants.

Untreated familial hypercholesterolemia (FH) leads to atherosclerosis and early cardiovascular disease. Mutations in the low-density lipoprotein receptor (LDLr) gene constitute the major cause of FH, and the high number of mutations already described in the LDLr makes necessary cascade screening or in vitro functional characterization to provide a definitive diagnosis. Implementation of high-predicting capacity software constitutes a valuable approach for assessing pathogenicity of LDLr variants to help in the early diagnosis and management of FH disease. This work provides a reliable machine learning model to accurately predict the pathogenicity of LDLr missense variants with specificity of 92.5% and sensitivity of 91.6%.