Cellular Senescence: Defining a Path Forward.

Cellular senescence is a cell state implicated in various physiological processes and a wide spectrum of age-related diseases. Recently, interest in therapeutically targeting senescence to improve healthy aging and age-related disease, otherwise known as senotherapy, has been growing rapidly. Thus, the accurate detection of senescent cells, especially in vivo, is essential. Here, we present a consensus from the International Cell Senescence Association (ICSA), defining and discussing key cellular and molecular features of senescence and offering recommendations on how to use them as biomarkers. We also present a resource tool to facilitate the identification of genes linked with senescence, SeneQuest (available at http://Senequest.net). Lastly, we propose an algorithm to accurately assess and quantify senescence, both in cultured cells and in vivo.

Digital Technology in Diagnostic Breast Pathology and Immunohistochemistry.

Digital technology has been used in the field of diagnostic breast pathology and immunohistochemistry (IHC) for decades. Examples include automated tissue processing and staining, digital data processing, storing and management, voice recognition systems, and digital technology-based production of antibodies and other IHC reagents. However, the recent application of whole slide imaging technology and artificial intelligence (AI)-based tools has attracted a lot of attention. The use of AI tools in breast pathology is discussed briefly as it is covered in other reviews. Here, we present the main application of digital technology in IHC. This includes automation of IHC staining, using image analysis systems and computer vision technology to interpret IHC staining, and the use of AI-based tools to predict marker expression from haematoxylin and eosin-stained digitalized images.

Systemic mechanisms and effects of ionizing radiation: A new 'old' paradigm of how the bystanders and distant can become the players.

Exposure of cells to any form of ionizing radiation (IR) is expected to induce a variety of DNA lesions, including double strand breaks (DSBs), single strand breaks (SSBs) and oxidized bases, as well as loss of bases, i.e., abasic sites. The damaging potential of IR is primarily related to the generation of electrons, which through their interaction with water produce free radicals. In their turn, free radicals attack DNA, proteins and lipids. Damage is induced also through direct deposition of energy. These types of IR interactions with biological materials are collectively called 'targeted effects', since they refer only to the irradiated cells. Earlier and sometimes 'anecdotal' findings were pointing to the possibility of IR actions unrelated to the irradiated cells or area, i.e., a type of systemic response with unknown mechanistic basis. Over the last years, significant experimental evidence has accumulated, showing a variety of radiation effects for 'out-of-field' areas (non-targeted effects-NTE). The NTE involve the release of chemical and biological mediators from the 'in-field' area and thus the communication of the radiation insult via the so called 'danger' signals. The NTE can be separated in two major groups: bystander and distant (systemic). In this review, we have collected a detailed list of proteins implicated in either bystander or systemic effects, including the clinically relevant abscopal phenomenon, using improved text-mining and bioinformatics tools from the literature. We have identified which of these genes belong to the DNA damage response and repair pathway (DDR/R) and made protein-protein interaction (PPi) networks. Our analysis supports that the apoptosis, TLR-like and NOD-like receptor signaling pathways are the main pathways participating in NTE. Based on this analysis, we formulate a biophysical hypothesis for the regulation of NTE, based on DNA damage and apoptosis gradients between the irradiation point and various distances corresponding to bystander (5mm) or distant effects (5cm). Last but not least, in order to provide a more realistic support for our model, we calculate the expected DSB and non-DSB clusters along the central axis of a representative 200.6MeV pencil beam calculated using Monte Carlo DNA damage simulation software (MCDS) based on the actual beam energy-to-depth curves used in therapy.

Proteomics Analysis of Pregnancy in Ewes under Heat Stress Conditions and Melatonin Administration.

Melatonin is an indoleamine with broad spectrum properties that acts as a regulator of antioxidant and immune response in organisms. In our previous studies, melatonin improved redox status and inflammatory response in pregnant ewes under heat stress conditions. In the present study, using proteomics, the proteins regulated by melatonin during different stages of pregnancy and lambing were assessed. Twenty-two ewes equally divided into two groups, the melatonin (M) (n = 11) and control (C) group (n = 11), participated in the study and were exposed to heat stress during the first months of pregnancy. In the M group, melatonin implants were administered throughout pregnancy, every 40 days, until parturition (a total of four implants per ewe). Blood samples were collected at the beginning of the study simultaneously with the administration of the first melatonin implant (blood samples M1, C1), mating (M2, C2), second implant (M3, C3), fourth implant (M4, C4) and parturition (M5, C5), and MALDI-TOF analysis was performed. The results revealed the existence of 42 extra proteins in samples M2, M3 and M4 and 53 in M5 (sample at parturition) that are linked to melatonin. The biological processes of these proteins refer to boosted immune response, the alleviation of oxidative and endoplasmic reticulum stress, energy metabolism, the protection of the maternal organism and embryo development. This proteomics analysis indicates that melatonin regulates protective mechanisms and controls cell proliferation under exogenous or endogenous stressful stimuli during pregnancy and parturition.

Calprotectin: a protein related to cardiovascular risk in adult patients with obstructive sleep apnea.

INTRODUCTION: Increased levels of inflammatory mediators, such as hs-CRP, have been detected in patients with obstructive sleep apnea (OSA) and used as cardiovascular risk and disease outcome predictors. Calprotectin is an inflammatory marker regulating atherogenic processes not investigated in adult OSA patients. The aim of the present study as primary objective was to examine the role of calprotectin as an inflammatory molecule, acting through a distinct pathway to the atherogenic process in adult OSA patients and its associations with hs-CRP and the lipidemic profile of the patients. As a secondary objective was the evaluation of the atherogenic markers post-CPAP treatment. MATERIALS AND METHODS: Seventy-four participants underwent full overnight polysomnography. Blood samples were collected for calprotectin, hs-CRP, total cholesterol, triglycerides, LDL, HDL and glucose levels. Thirty-two OSA patients were reexamined 6 months post-CPAP treatment. RESULTS: Out of 74 participants included in the study, 33 had moderate OSA, 27 had severe OSA and 14 were controls. Calprotectin and hs-CRP were significantly increased in patients with moderate and severe OSA compared to controls (p<0.0001). Calprotectin and hs-CRP levels were positively correlated with apnea-hypopnea index, BMI and total time of sleep with SaO(2)<90% and inversely correlated with SaO(2) minimum and mean values. Calprotectin and hs-CRP levels were significantly improved post-CPAP treatment (p<0.0001). DISCUSSION: Calprotectin may serve as a novel and reliable, biomarker of cardiovascular risk severity in OSA patients. The decrease of calprotectin levels post-CPAP treatment combined with hs-CRP amelioration could provide evidence for reduction of cardiovascular risk post CPAP treatment.

The Role of Early Warning Scoring Systems NEWS and MEWS in the Acute Exacerbation of COPD.

Acute exacerbations of chronic obstructive pulmonary disease (AECOPD) are the most devastating events in the course of the disease. Our aim was to investigate the value of early warning scoring systems: National Early Warning Score (NEWS) and Modified Early Warning Score (MEWS) in AECOPD. This is a prospective observational study of patients with AECOPD who were admitted at hospital. The NEWS and MEWS scores were registered at admission (NEWS-d1, MEWS-d1) and on the second day (NEWS-d2, MEWS-d2). A nasopharyngeal and sputum sample was taken for culture. Follow-up was done at 3 and 6 months after hospitalization. Any possible correlations between NEWS and MEWS and other parameters of COPD were explored. A cohort of 64 patients were included. In-hospital mortality was 4.7% while total mortality at 6 months was 26%. We did not find any significant correlation between in-hospital mortality and any of the scores but we could show a higher mortality and more frequent AECOPD at 6 months of follow-up for those with higher NEWS-d2. NEWS-d2 was associated with higher pCO2 at presentation and a more frequent use of NIV. Higher NEWS-d1 and NEWS-d2 were predictive of a longer hospital stay. The presence of pathogens in the nasopharyngeal sample was related with a higher reduction of both scores on the second day. We therefore support the superiority of NEWS in the evaluation of hospitalized patients with AECOPD. A remaining high NEWS at the second day of hospital stay signals a high risk of hypercapnia and need of NIV but also higher mortality and more frequent exacerbations at 6 months after AECOPD.

Brain proteome response following whole body exposure of mice to mobile phone or wireless DECT base radiation.

The objective of this study was to investigate the effects of two sources of electromagnetic fields (EMFs) on the proteome of cerebellum, hippocampus, and frontal lobe in Balb/c mice following long-term whole body irradiation. Three equally divided groups of animals (6 animals/group) were used; the first group was exposed to a typical mobile phone, at a SAR level range of 0.17-0.37 W/kg for 3 h daily for 8 months, the second group was exposed to a wireless DECT base (Digital Enhanced Cordless Telecommunications/Telephone) at a SAR level range of 0.012-0.028 W/kg for 8 h/day also for 8 months and the third group comprised the sham-exposed animals. Comparative proteomics analysis revealed that long-term irradiation from both EMF sources altered significantly (p < 0.05) the expression of 143 proteins in total (as low as 0.003 fold downregulation up to 114 fold overexpression). Several neural function related proteins (i.e., Glial Fibrillary Acidic Protein (GFAP), Alpha-synuclein, Glia Maturation Factor beta (GMF), and apolipoprotein E (apoE)), heat shock proteins, and cytoskeletal proteins (i.e., Neurofilaments and tropomodulin) are included in this list as well as proteins of the brain metabolism (i.e., Aspartate aminotransferase, Glutamate dehydrogenase) to nearly all brain regions studied. Western blot analysis on selected proteins confirmed the proteomics data. The observed protein expression changes may be related to brain plasticity alterations, indicative of oxidative stress in the nervous system or involved in apoptosis and might potentially explain human health hazards reported so far, such as headaches, sleep disturbance, fatigue, memory deficits, and brain tumor long-term induction under similar exposure conditions.

Analysis of Human and Microbial Salivary Proteomes in Children Offers Insights on the Molecular Pathogenesis of Molar-Incisor Hypomineralization.

Molar incisor hypomineralization is a complex developmental enamel defect that affects the permanent dentition of children with significant functional and aesthetic implications. Saliva is an ideal diagnostic tool and ensures patients' compliance by diminishing the discomfort especially in pediatric population. Lately, salivary proteome analysis has progressively evolved in various biomedical disciplines. As changes in saliva composition are associated with oral diseases, it is reasonable to assume that the saliva proteome of MIH-affected children might be altered compared to healthy children. This study analyzed the human and microbial salivary proteome in children with MIH in order to identify salivary markers indicative of the pathology. The conducted proteomic analysis generated a comprehensive dataset comprising a total of 1515 high confidence identifications and revealed a clear discrimination between the two groups. Statistical comparison identified 142 differentially expressed proteins, while the pathway analysis indicated deregulation of inflammation, immune response mechanisms, and defense response to bacteria in MIH patients. Bacterial proteome analysis showed a lower diversity for the microbial species, which highlights the dysbiotic environment established in the MIH pathology.

Proteomics studies of childhood pilocytic astrocytoma.

Childhood pilocytic astrocytoma is the most frequent brain tumor affecting children. Proteomics analysis is currently considered a powerful tool for global evaluation of protein expression and has been widely applied in the field of cancer research. In the present study, a series of proteomics, genomics, and bioinformatics approaches were employed to identify, classify and characterize the proteome content of low-grade brain tumors as it appears in early childhood. Through bioinformatics database construction, protein profiles generated from pathological tissue samples were compared against profiles of normal brain tissues. Additionally, experiments of comparative genomic hybridization arrays were employed to monitor for genetic aberrations and sustain the interpretation and evaluation of the proteomic data. The current study confirms the dominance of MAPK pathway for the childhood pilocytic astrocytoma occurrence and novel findings regarding the ERK-2 expression are reported.

Proteomic analysis of childhood de novo acute myeloid leukemia and myelodysplastic syndrome/AML: correlation to molecular and cytogenetic analyses.

The aim of this study was to investigate the progression of myelodysplastic syndrome (MDS) to acute myeloid leukemia (AML) and to provide additional data regarding the proteomic analysis of AML. The protein profiles obtained were correlated to cytogenetic and molecular analyses. Bone marrow (BM) and peripheral blood (PB) samples were obtained during MDS diagnosis, at MDS transformation to AML, at de novo AML diagnosis and 3 months following treatment. As controls, non-leukemic pediatric patients were studied. Cytogenetic and molecular analyses were carried out by G banding and polymerase chain reaction followed by sequencing, respectively. Differential proteomic analysis was performed by two-dimensional gel electrophoresis and protein identification by matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry. No significant correlations were noted between protein patterns and cytogenetic or molecular analyses. Certain suppressor genes, metabolic enzymes, immunoglobulins and actin-binding proteins were differentially expressed by BM or PB plasma and cell lysates compared to controls. The obtained data showed that vitamin D and gelsolin played contradicting roles in contributing and restraining leukemogenesis, while MOES, EZRI and AIFM1 could be considered as biomarkers for AML.

Downregulation of Salivary Proteins, Protective against Dental Caries, in Type 1 Diabetes.

Saliva, an essential oral secretion involved in protecting the oral cavity's hard and soft tissues, is readily available and straightforward to collect. Recent studies have analyzed the salivary proteome in children and adolescents with extensive carious lesions to identify diagnostic and prognostic biomarkers. The current study aimed to investigate saliva's diagnostic ability through proteomics to detect the potential differential expression of proteins specific for the occurrence of carious lesions. For this study, we performed bioinformatics and functional analysis of proteomic datasets, previously examined by our group, from samples of adolescents with regulated and unregulated type 1 diabetes, as they compare with healthy controls. Among the differentially expressed proteins relevant to caries pathology, alpha-amylase 2B, beta-defensin 4A, BPI fold containing family B member 2, protein S100-A7, mucin 5B, statherin, salivary proline-rich protein 2, and interleukin 36 gamma were significantly downregulated in poorly-controlled patients compared to healthy subjects. In addition, significant biological pathways (defense response to the bacterium, beta-defensin activity, proline-rich protein activity, oxygen binding, calcium binding, and glycosylation) were deregulated in this comparison, highlighting specific molecular characteristics in the cariogenic process. This analysis contributes to a better understanding of the mechanisms involved in caries vulnerability in adolescents with unregulated diabetes.

Discovery and Pharmacological Evaluation of STEAP4 as a Novel Target for HER2 Overexpressing Breast Cancer.

Breast cancer (BC) is a highly heterogeneous disease encompassing multiple subtypes with different molecular and histopathological features, disease prognosis, and therapeutic responses. Among these, the Triple Negative BC form (TNBC) is an aggressive subtype with poor prognosis and therapeutic outcome. With respect to HER2 overexpressing BC, although advanced targeted therapies have improved the survival of patients, disease relapse and metastasis remains a challenge for therapeutic efficacy. In this study the aim was to identify key membrane-associated proteins which are overexpressed in these aggressive BC subtypes and can serve as potential biomarkers or drug targets. We leveraged on the development of a membrane enrichment protocol in combination with the global profiling GeLC-MS/MS technique, and compared the proteomic profiles of a HER2 overexpressing (HCC-1954) and a TNBC (MDA-MB-231) cell line with that of a benign control breast cell line (MCF-10A). An average of 2300 proteins were identified from each cell line, of which approximately 600 were membrane-associated proteins. Our global proteomic methodology in tandem with invigoration by Western blot and Immunofluorescence analysis, readily detected several previously-established BC receptors like HER2 and EPHA2, but importantly STEAP4 and CD97 emerged as novel potential candidate markers. This is the first time that the mitochondrial iron reductase STEAP4 protein up-regulation is linked to BC (HER2+ subtype), while for CD97, its role in BC has been previously described, but never before by a global proteomic technology in TNBC. STEAP4 was selected for further detailed evaluation by the employment of Immunohistochemical analysis of BC xenografts and clinical tissue microarray studies. Results showed that STEAP4 expression was evident only in malignant breast tissues whereas all the benign breast cases had no detectable levels. A functional role of STEAP4 intervention was established in HER2 overexpressing BC by pharmacological studies, where blockage of the STEAP4 pathway with an iron chelator (Deferiprone) in combination with the HER2 inhibitor Lapatinib led to a significant reduction in cell growth in vitro. Furthermore, siRNA mediated knockdown of STEAP4 also suppressed cell proliferation and enhanced the inhibition of Lapatinib in HER2 overexpressing BC, confirming its potential oncogenic role in BC. In conclusion, STEAP4 may represent a novel BC related biomarker and a potential pharmacological target for the treatment of HER2 overexpressing BC.

Potential biomarkers for Turner in maternal plasma: possibility for noninvasive prenatal diagnosis.

Turner syndrome (TS) is the most common sex chromosome abnormality in females, caused by the complete or partial absence of one X chromosome. To identify biomarkers for TS, we compared the protein composition of maternal plasma samples from pregnant women with normal and TS fetuses, using a proteomic approach consisting of 2D-E separation and MS analysis for the identification of the differentially expressed proteins. Samples were routinely obtained in the second trimester of pregnancy, stored, and used after prenatal determination of the fetal karyotype. Nine proteins (C1S, CO3, CLUS, AFAM, HABP2, IGHA1, HPT, SHBG, and CD5L) were significantly increased in the plasma of women carrying TS fetuses, whereas KNG1, IGJ, and TTHY were decreased. Identified proteins were further evaluated by immunoblot analysis while functional network association was carried out to asses significance. The identification of specific biomarkers may facilitate the development of noninvasive prenatal diagnosis and improve our understanding of the pathology of TS. Nevertheless, testing a larger cohort of pregnant women is necessary to evaluate the relevance of the reported findings.

UniMaP: finding unique mass and peptide signatures in the human proteome.

UNLABELLED: The uniqueness of a measured molecular mass or peptide sequence plays a very important role in the fields of protein identification and peptide/protein-biomarker investigation. We present a publicly available web application that offers information concerning the uniqueness of one or more molecular masses and one or more peptide sequences in the human proteome. When a sequence is found to be unique in humans, the application is able to search across all species querying whether this sequence is unique, not only in humans but also in other species found in the Swiss-Prot Database. The application is also able to search for unique protein fragments derived computationally from enzymatic digestion driven by certain enzymes. Furthermore, the application can list all the unique masses and peptides of a given protein. Through this application, researchers are able to find unique tags, either on a molecular mass level or on a sequence level. These unique tags are remarkably important in research related to protein identification or biomarker discovery and measurements. AVAILABILITY: UniMaP web-application is available at http://bioserver-1.bioacademy.gr/Bioserver/UniMaP/

Modulation of the E2F1-driven cancer cell fate by the DNA damage response machinery and potential novel E2F1 targets in osteosarcomas.

Osteosarcoma is the most common primary bone cancer. Mutations of the RB gene represent the most frequent molecular defect in this malignancy. A major consequence of this alteration is that the activity of the key cell cycle regulator E2F1 is unleashed from the inhibitory effects of pRb. Studies in animal models and in human cancers have shown that deregulated E2F1 overexpression possesses either "oncogenic" or "oncosuppressor" properties, depending on the cellular context. To address this issue in osteosarcomas, we examined the status of E2F1 relative to cell proliferation and apoptosis in a clinical setting of human primary osteosarcomas and in E2F1-inducible osteosarcoma cell line models that are wild-type and deficient for p53. Collectively, our data demonstrated that high E2F1 levels exerted a growth-suppressing effect that relied on the integrity of the DNA damage response network. Surprisingly, induction of p73, an established E2F1 target, was also DNA damage response-dependent. Furthermore, a global proteome analysis associated with bioinformatics revealed novel E2F1-regulated genes and potential E2F1-driven signaling networks that could provide useful targets in challenging this aggressive neoplasm by innovative therapies.

Marked defects in the expression and glycosylation of alpha2-HS glycoprotein/fetuin-A in plasma from neonates with intrauterine growth restriction: proteomics screening and potential clinical implications.

Intrauterine growth restriction (IUGR) has been associated with increased perinatal morbidity and mortality and increased morbidity and metabolic abnormalities later in life. IUGR is characterized as the failure of a fetus to achieve his or her genetic growth potential in utero. Altered protein expression profiles associated with IUGR may be informative on the pathologic mechanisms of this condition and might reveal potential markers for postnatal complications. The aim of this study was to compare protein profiles of umbilical cord plasma from IUGR and appropriate for gestational age full-term neonates. Blood samples from doubly clamped umbilical cord at delivery from 10 IUGR and 10 appropriate for gestational age full-term neonates were analyzed by two-dimensional electrophoresis and MS. Prominent changes of the alpha2-HS glycoprotein/fetuin-A were observed in IUGR cases. Specifically we showed that these changes occur primarily at the level of post-translational modifications of the protein. Using a combination of mass spectrometry and classical biochemical assays, single and heavy chain forms of fetuin-A were found to lack the normally present O-linked sialic acids in IUGR neonates. Fetuin A is a glycoprotein that has been associated with promotion of in vitro cell replication, fetal growth and osteogenesis, and protection from Gram-negative bacterial endotoxins. Prominent defects in glycosylation/sialylation of fetuin-A revealed by our study might be responsible for impaired function of fetuin-A, leading to deficient fetal growth, especially osteogenesis, and/or to the development of complications frequently seen later in the lives of IUGR neonates.

A two-dimensional proteomic profile of Tetrahymena thermophila whole cell lysate.

Tetrahymena thermophila is a unicellular eukaryotic model organism used for a variety of biochemical, molecular and biological studies. According to its macronucleus genome sequence, it is expected to contain more than 27,000 protein-coding genes, although only a small proportion of them have information published specifically about them. Here, we present a reference map for whole cell lysate of T. thermophila obtained using two-dimensional gel electrophoresis (2-DE) combined with mass spectrometry. Although (2-DE) is one of the most efficient techniques for resolving complex protein mixtures and revealing the relative high-abundance proteins, it has not yet been applied generally to ciliates. In order to obtain qualitative protein samples for analysis, an appropriate homogenization method is required. Optimization of the homogenization method led to the analysis of nearly 4500 protein spots, the final identification of 375 different proteins using Mascot software and an additional 258 gene products using a newly developed web service, called Peptide Finder, resulting in a total of 631 different gene products that are considered to constitute the proteomic profile of the whole cell lysate of T. thermophila.

Deep learning: shaping the medicine of tomorrow.

Predicting response to therapy is a major challenge in medicine. Machine learning algorithms are promising tools for assisting this aim. Amongst them, Deep Neural Networks are emerging as the most capable of interrogating across multiple data types. Their further development will lead to sophisticated knowledge extraction, shaping the medicine of tomorrow.

Proteomic advances in salivary diagnostics.

Saliva is identified as functional equivalent to serum, reflecting the physiological state of the body, as well as hormonal, emotional, nutritional and metabolic alterations. The application of mass spectrometry based approaches has allowed a thorough characterization of the saliva proteome and led to the discovery of putative biomarkers. Several salivary biomarkers have been recently explored as potentially useful screening tools in patients diagnosed with metabolic disorders. In this review, we provide an overview of saliva proteomics studies, with a focus on diabetes, and we explore the evidence for the utility of well identified markers for the diagnosis and monitoring of the disease. Emerging approaches in salivary diagnostics that may significantly advance the field of diabetes research are also highlighted.

Expression profiling meta-analysis of ACE2 and TMPRSS2, the putative anti-inflammatory receptor and priming protease of SARS-CoV-2 in human cells, and identification of putative modulators.

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and has resulted in hundreds of thousands of deaths worldwide. While the majority of people with COVID-19 won't require hospitalization, those who do may experience severe life-threatening complications, including acute respiratory distress syndrome. SARS-CoV-2 infects human cells by binding to the cellular surface protein angiotensin-converting enzyme 2 (ACE2); in addition, the cellular transmembrane serine protease 2 (TMPRSS2) is needed for priming of the spike (S) protein of the virus. Virus entry may also depend on the activity of the endosomal/lysosomal cysteine proteases cathepsin B, L (CTSB, CTSL) although their activity is likely dispensable. Given that the uncertainty of how COVID-19 kills, hampers doctors' ability to choose treatments the need for a deep understanding of COVID-19 biology is urgent. Herein, we performed an expression profiling meta-analysis of ACE2, TMPRSS2 and CTSB/L genes (and proteins) in public repository databases and found that all are widely expressed in human tissues; also, the ACE2 and TMPRSS2 genes tend to be co-regulated. The ACE2 and TMPRSS genes expression is (among others) suppressed by TNF, and is induced by pro-inflammatory conditions including obesity, Barrett's esophagus, stomach infection by helicobacter pylori, diabetes, autoimmune diseases and oxidized LDL; by exercise, as well as by growth factors, viruses' infections, cigarette smoke, interferons and androgens. Regarding currently investigated therapies interferon-beta induced ACE2 gene expression in bronchial epithelial cells, while chloroquine tends to upregulate CTSB/L genes. Finally, we analyzed KEGG pathways modulated by ACE2, TMPRSS2 and CTSB/L and probed DrugBank for drugs that target modules of the affected pathways. Our data indicate possible novel high-risk groups for COVID-19; provide a rich resource for future investigations of its pathogenesis and highlight the therapeutic challenges we face.