Virus-induced genetics revealed by multidimensional precision medicine transcriptional workflow applicable to COVID-19

ntEmbd: Deep learning embedding for nucleotide sequences

Enabled by the explosion of data and substantial increase in computational power, deep learning has transformed fields such as computer vision and natural language processing (NLP) and it has become a successful method to be applied to many transcriptomic analysis tasks. A core advantage of deep learning is its inherent capability to incorporate feature computation within the machine learning models. This results in a comprehensive and machine-readable representation of sequences, facilitating the downstream classification and clustering tasks. Compared to machine translation problems in NLP, feature embedding is particularly challenging for transcriptomic studies as the sequences are string of thousands of nucleotides in length, which make the long-term dependencies between features from different parts of the sequence even more difficult to capture. This highlights the need for nucleotide sequence embedding methods that are capable of learning input sequence features implicitly. Here we introduce ntEmbd, a deep learning embedding tool that captures dependencies between different features of the sequences and learns a latent representation for given nucleotide sequences. We further provide two sample use cases, describing how learned RNA features can be used in downstream analysis. The first use case demonstrates ntEmbd's utility in classifying coding and noncoding RNA benchmarked against existing tools, and the second one explores the utility of learned representations in identifying adapter sequences in nanopore RNA-seq reads. The tool as well as the trained models are freely available on GitHub at https://github.com/bcgsc/ntEmbd.

The complex, dynamic SpliceOme of the small GTPase transcripts altered by technique, sex, genetics, tissue specificity, and RNA base editing

The small GTPase family is well-studied in cancer and cellular physiology. With 162 annotated human genes, the family has a broad expression throughout cells of the body. Members of the family have multiple exons that require splicing. Yet, the role of splicing within the family has been underexplored. We have studied the splicing dynamics of small GTPases throughout 41,671 samples by integrating Nanopore and Illumina sequencing techniques. Within this work, we have made several discoveries. 1). Using the GTEx long read data of 92 samples, each small GTPase gene averages two transcripts, with 83 genes (51%) expressing two or more isoforms. 2). Cross-tissue analysis of GTEx from 17,382 samples shows 41 genes (25%) expressing two or more protein-coding isoforms. These include protein-changing transcripts in genes such as RHOA, RAB37, RAB40C, RAB4B, RAB5C, RHOC, RAB1A, RAN, RHEB, RAC1, and KRAS. 3). The isolation and library technique of the RNAseq influences the abundance of non-sense-mediated decay and retained intron transcripts of small GTPases, which are observed more often in genes than appreciated. 4). Analysis of 16,243 samples of "Blood PAXgene" identified seven genes (3.7%; RHOA, RAB40C, RAB4B, RAB37, RAB5B, RAB5C, RHOC) with two or more transcripts expressed as the major isoform (75% of the total gene), suggesting a role of genetics in altering splicing. 5). Rare (ARL6, RAB23, ARL13B, HRAS, NRAS) and common variants (GEM, RHOC, MRAS, RAB5B, RERG, ARL16) can influence splicing and have an impact on phenotypes and diseases. 6). Multiple genes (RAB9A, RAP2C, ARL4A, RAB3A, RAB26, RAB3C, RASL10A, RAB40B, and HRAS) have sex differences in transcript expression. 7). Several exons are included or excluded for small GTPase genes (RASEF, KRAS, RAC1, RHEB, ARL4A, RHOA, RAB30, RHOBTB1, ARL16, RAP1A) in one or more forms of cancer. 8). Ten transcripts are altered in hypoxia (SAR1B, IFT27, ARL14, RAB11A, RAB10, RAB38, RAN, RIT1, RAB9A) with RHOA identified to have a transient 3'UTR RNA base editing at a conserved site found in all of its transcripts. Overall, we show a remarkable and dynamic role of splicing within the small GTPase family that requires future explorations.

Autosomal recessive LRP1-related syndrome featuring cardiopulmonary dysfunction, bone dysmorphology, and corneal clouding

We provide the first study of two siblings with a novel autosomal recessive LRP1-related syndrome identified by rapid genome sequencing and overlapping multiple genetic models. The patients presented with respiratory distress, congenital heart defects, hypotonia, dysmorphology, and unique findings, including corneal clouding and ascites. Both siblings had compound heterozygous damaging variants, c.11420G > C (p.Cys3807Ser) and c.12407T > G (p.Val4136Gly) in LRP1, in which segregation analysis helped dismiss additional variants of interest. LRP1 analysis using multiple human/mouse data sets reveals a correlation to patient phenotypes of Peters plus syndrome with additional severe cardiomyopathy and blood vessel development complications linked to neural crest cells.

Introduction to Personalized Medicine in Pediatrics

Exciting new developments in biomedical and computational sciences provide an extraordinary and unparalleled opportunity to compile, connect, and analyze multiple types of "big data," driving the development of personalized medicine. These insights must begin in early life (ie, pregnancy, neonatal, and infancy) and focus on early prevention, diagnosis, and intervention-areas of medicine where pediatricians are poised to lead the way to a personalized medicine future. The rapid growth of genomics (including pharmacogenomics), transcriptomics, and related "omics" has revolutionized the diagnosis of rare monogenic disorders. It is now clarifying the pathogenesis of complex conditions ranging from autism spectrum disorder to asthma. Collaborations between clinicians and basic scientists integrating multiomics approaches in evaluating children with severe illness are transforming the fields of perinatal, neonatal, and pediatric critical care medicine. Improvements in rapid diagnostic and prognostic information suggest that pediatric personalized medicine is under way and has an exciting future. [Pediatr Ann. 2022;51(10):e381-e386.].

Computational and Experimental Analysis of Genetic Variants

Genomics has grown exponentially over the last decade. Common variants are associated with physiological changes through statistical strategies such as Genome-Wide Association Studies (GWAS) and quantitative trail loci (QTL). Rare variants are associated with diseases through extensive filtering tools, including population genomics and trio-based sequencing (parents and probands). However, the genomic associations require follow-up analyses to narrow causal variants, identify genes that are influenced, and to determine the physiological changes. Large quantities of data exist that can be used to connect variants to gene changes, cell types, protein pathways, clinical phenotypes, and animal models that establish physiological genomics. This data combined with bioinformatics including evolutionary analysis, structural insights, and gene regulation can yield testable hypotheses for mechanisms of genomic variants. Molecular biology, biochemistry, cell culture, CRISPR editing, and animal models can test the hypotheses to give molecular variant mechanisms. Variant characterizations can be a significant component of educating future professionals at the undergraduate, graduate, or medical training programs through teaching the basic concepts and terminology of genetics while learning independent research hypothesis design. This article goes through the computational and experimental analysis strategies of variant characterization and provides examples of these tools applied in publications. © 2022 American Physiological Society. Compr Physiol 12:3303-3336, 2022.

CCR5 and Biological Complexity: The Need for Data Integration and Educational Materials to Address Genetic/Biological Reductionism at the Interface of Ethical, Legal, and Social Implications

In the age of genomics, public understanding of complex scientific knowledge is critical. To combat reductionistic views, it is necessary to generate and organize educational material and data that keep pace with advances in genomics. The view that CCR5 is solely the receptor for HIV gave rise to demand to remove the gene in patients to create host HIV resistance, underestimating the broader roles and complex genetic inheritance of CCR5. A program aimed at providing research projects to undergraduates, known as CODE, has been expanded to build educational material for genes such as CCR5 in a rapid approach, exposing students and trainees to large bioinformatics databases and previous experiments for broader data to challenge commitment to biological reductionism. Our students organize expression databases, query environmental responses, assess genetic factors, generate protein models/dynamics, and profile evolutionary insights into a protein such as CCR5. The knowledgebase generated in the initiative opens the door for public educational information and tools (molecular videos, 3D printed models, and handouts), classroom materials, and strategy for future genetic ideas that can be distributed in formal, semiformal, and informal educational environments. This work highlights that many factors are missing from the reductionist view of CCR5, including the role of missense variants or expression of CCR5 with neurological phenotypes and the role of CCR5 and the delta32 variant in complex critical care patients with sepsis. When connected to genomic stories in the news, these tools offer critically needed Ethical, Legal, and Social Implication (ELSI) education to combat biological reductionism.

Cell population-based framework of genetic epidemiology in the single-cell omics era

Genetic epidemiology is a rapidly advancing field due to the recent availability of large amounts of omics data. In recent years, it has become possible to obtain omics information at the single-cell level, so genetic epidemiological models need to be updated to integrate with single-cell expression data. In this perspective paper, we propose a cell population-based framework for genetic epidemiology in the single-cell era. In this framework, genetic diversity influences phenotypic diversity through the diversity of cell population profiles, which are defined as high-dimensional probability distributions of the state spaces of biomolecules of each omics layer. We discuss how biomolecular experimental measurement data can capture the different properties of this distribution. In particular, single-cell data constitute a sample from this population distribution where only some coordinate values are observable. From a data analysis standpoint, we introduce methodology for feature extraction from cell population profiles. Finally, we discuss how this framework can be applied not only to genetic epidemiology but also to systems biology.

High-Density Blood Transcriptomics Reveals Precision Immune Signatures of SARS-CoV-2 Infection in Hospitalized Individuals

The immune response to COVID-19 infection is variable. How COVID-19 influences clinical outcomes in hospitalized patients needs to be understood through readily obtainable biological materials, such as blood. We hypothesized that a high-density analysis of host (and pathogen) blood RNA in hospitalized patients with SARS-CoV-2 would provide mechanistic insights into the heterogeneity of response amongst COVID-19 patients when combined with advanced multidimensional bioinformatics for RNA. We enrolled 36 hospitalized COVID-19 patients (11 died) and 15 controls, collecting 74 blood PAXgene RNA tubes at multiple timepoints, one early and in 23 patients after treatment with various therapies. Total RNAseq was performed at high-density, with >160 million paired-end, 150 base pair reads per sample, representing the most sequenced bases per sample for any publicly deposited blood PAXgene tube study. There are 770 genes significantly altered in the blood of COVID-19 patients associated with antiviral defense, mitotic cell cycle, type I interferon signaling, and severe viral infections. Immune genes activated include those associated with neutrophil mechanisms, secretory granules, and neutrophil extracellular traps (NETs), along with decreased gene expression in lymphocytes and clonal expansion of the acquired immune response. Therapies such as convalescent serum and dexamethasone reduced many of the blood expression signatures of COVID-19. Severely ill or deceased patients are marked by various secondary infections, unique gene patterns, dysregulated innate response, and peripheral organ damage not otherwise found in the cohort. High-density transcriptomic data offers shared gene expression signatures, providing unique insights into the immune system and individualized signatures of patients that could be used to understand the patient's clinical condition. Whole blood transcriptomics provides patient-level insights for immune activation, immune repertoire, and secondary infections that can further guide precision treatment.

Hexosylceramides and Glycerophosphatidylcholine GPC(36:1) Increase in Multi-Organ Dysfunction Syndrome Patients with Pediatric Intensive Care Unit Admission over 8-Day Hospitalization

Glycero- and sphingo-lipids are important in plasma membrane structure, caloric storage and signaling. An un-targeted lipidomics approach for a cohort of critically ill pediatric intensive care unit (PICU) patients undergoing multi-organ dysfunction syndrome (MODS) was compared to sedation controls. After IRB approval, patients meeting the criteria for MODS were screened, consented (n = 24), and blood samples were collected from the PICU at HDVCH, Michigan; eight patients needed veno-arterial extracorporeal membrane oxygenation (VA ECMO). Sedation controls were presenting for routine sedation (n = 4). Plasma lipid profiles were determined by nano-electrospray (nESI) direct infusion high resolution/accurate mass spectrometry (MS) and tandem mass spectrometry (MS/MS). Biostatistics analysis was performed using R v 3.6.0. Sixty-one patient samples over three time points revealed a ceramide metabolite, hexosylceramide (Hex-Cer) was high across all time points (mean 1.63-3.19%; vs. controls 0.22%). Fourteen species statistically differentiated from sedation controls (p-value ≤ 0.05); sphingomyelin (SM) [SM(d18:1/23:0), SM(d18:1/22:0), SM(d18:1/23:1), SM(d18:1/21:0), SM(d18:1/24:0)]; and glycerophosphotidylcholine (GPC) [GPC(36:01), GPC(18:00), GPC(O:34:02), GPC(18:02), GPC(38:05), GPC(O:34:03), GPC(16:00), GPC(40:05), GPC(O:36:03)]. Hex-Cer has been shown to be involved in viral infection and may be at play during acute illness. GPC(36:01) was elevated in all MODS patients at all time points and is associated with inflammation and brain injury.

Balancing precision versus cohort transcriptomic analysis of acute and recovery phase of viral bronchiolitis

Viral infections affecting the lower respiratory tract place enormous burdens on hospitals. As neither vaccines nor antiviral agents exist for many viruses, understanding risk factors and outcomes in each patient using minimally invasive analysis, such as blood, can lead to improved health care delivery. A cohort of PAXgene RNA sequencing of infants admitted with moderate or severe acute bronchiolitis and respiratory syncytial virus were compared with case-control statistical analysis and cohort-based outlier mapping for precision transcriptomics. Patients with severe bronchiolitis had signatures connected to the immune system, interferon signaling, and cytokine signaling, with marked sex differences in XIST, RPS4Y1, KDM5D, and LINC00278 for severity. Several patients had unique secondary infections, cytokine activation, immune responses, biological pathways, and immune cell activation, highlighting the need for defining patient-level transcriptomic signatures. Balancing relative contributions of cohort-based biomarker discoveries with patient's biological responses is needed to understand the totality of mechanisms of adverse outcomes in viral bronchiolitis.

Pediatric Multi-Organ Dysfunction Syndrome: Analysis by an Untargeted "Shotgun" Lipidomic Approach Reveals Low-Abundance Plasma Phospholipids and Dynamic Recovery over 8-Day Period, a Single-Center Observational Study

Lipids are molecules involved in metabolism and inflammation. This study investigates the plasma lipidome for markers of severity and nutritional status in critically ill children. Children with multi-organ dysfunction syndrome (MODS) (n = 24) are analyzed at three time-points and cross-referenced to sedation controls (n = 4) for a total of N = 28. Eight of the patients with MODS, needed veno-arterial extracorporeal membrane oxygenation (VA ECMO) support to survive. Blood plasma lipid profiles are quantified by nano-electrospray (nESI), direct infusion high resolution/accurate mass spectrometry (MS), and tandem mass spectrometry (MS/MS), and compared to nutritional profiles and pediatric logistic organ dysfunction (PELOD) scores. Our results show that PELOD scores were not significantly different between MODS and ECMO cases across time-points (p = 0.66). Lipid profiling provides stratification between sedation controls and all MODS patients for total lysophosphatidylserine (lysoPS) (p-value = 0.004), total phosphatidylserine (PS) (p-value = 0.015), and total ether-linked phosphatidylethanolamine (ether-PE) (p-value = 0.03) after adjusting for sex and age. Nutrition intake over time did not correlate with changes in lipid profiles, as measured by caloric and protein intake. Lipid measurement in the intensive care environment shows dynamic changes over an 8-day pediatric intensive care unit (PICU) course, suggesting novel metabolic indicators for defining critically ill children.

In silico Analyses of Immune System Protein Interactome Network, Single-Cell RNA Sequencing of Human Tissues, and Artificial Neural Networks Reveal Potential Therapeutic Targets for Drug Repurposing Against COVID-19

Background: There is pressing urgency to identify therapeutic targets and drugs that allow treating COVID-19 patients effectively. Methods: We performed in silico analyses of immune system protein interactome network, single-cell RNA sequencing of human tissues, and artificial neural networks to reveal potential therapeutic targets for drug repurposing against COVID-19. Results: We screened 1,584 high-confidence immune system proteins in ACE2 and TMPRSS2 co-expressing cells, finding 25 potential therapeutic targets significantly overexpressed in nasal goblet secretory cells, lung type II pneumocytes, and ileal absorptive enterocytes of patients with several immunopathologies. Then, we performed fully connected deep neural networks to find the best multitask classification model to predict the activity of 10,672 drugs, obtaining several approved drugs, compounds under investigation, and experimental compounds with the highest area under the receiver operating characteristics. Conclusion: After being effectively analyzed in clinical trials, these drugs can be considered for treatment of severe COVID-19 patients. Scripts can be downloaded at https://github.com/muntisa/immuno-drug-repurposing-COVID-19.

The Feasibility of Studying Metabolites in PICU Multi-Organ Dysfunction Syndrome Patients over an 8-Day Course Using an Untargeted Approach

Metabolites are generated from critical biological functions and metabolism. This pediatric study reviewed plasma metabolites in patients suffering from multi-organ dysfunction syndrome (MODS) in the pediatric intensive care unit (PICU) using an untargeted metabolomics approach. Patients meeting the criteria for MODS were screened for eligibility and consented (n = 24), and blood samples were collected at baseline, 72 h, and 8 days; control patients (n = 4) presented for routine sedation in an outpatient setting. A subset of MODS patients (n = 8) required additional support with veno-atrial extracorporeal membrane oxygenation (VA-ECMO) therapy. Metabolites from thawed blood plasma were determined from ion pairing reversed-phase liquid chromatography–mass spectrometry (LC-MS) analysis. Chromatographic peak alignment, identification, relative quantitation, and statistical and bioinformatics evaluation were performed using MAVEN and MetaboAnalyst 4.0. Metabolite analysis revealed 115 peaks per sample. From the partial least squares-discriminant analysis (PLS-DA) with variance of importance (VIP) scores above ≥2.0, 7 dynamic metabolites emerged over the three time points: tauro-chenodeoxycholic acid (TCDCA), hexose, p-hydroxybenzoate, hydroxyphenylacetic acid (HPLA), 2_3-dihydroxybenzoic acid, 2-keto-isovalerate, and deoxyribose phosphate. After Bonferroni adjustment for repeated measures, hexose and p-hydroxybenzoate were significant at one time point or more. Kendall’s tau-b test was used for internal validation of creatinine. Metabolites may be benign or significant in describing a patient’s pathophysiology and require operator interpretation.

Precision Medicine and Public Health: New Challenges for Effective and Sustainable Health

The development of high-throughput omics technologies represents an unmissable opportunity for evidence-based prevention of adverse effects on human health. However, the applicability and access to multi-omics tests are limited. In Italy, this is due to the rapid increase of knowledge and the high levels of skill and economic investment initially necessary. The fields of human genetics and public health have highlighted the relevance of an implementation strategy at a national level in Italy, including integration in sanitary regulations and governance instruments. In this review, the emerging field of public health genomics is discussed, including the polygenic scores approach, epigenetic modulation, nutrigenomics, and microbiomes implications. Moreover, the Italian state of implementation is presented. The omics sciences have important implications for the prevention of both communicable and noncommunicable diseases, especially because they can be used to assess the health status during the whole course of life. An effective population health gain is possible if omics tools are implemented for each person after a preliminary assessment of effectiveness in the medium to long term.

SARS-CoV-2 infection: molecular mechanisms of severe outcomes to suggest therapeutics

Introduction:The year 2020 was defined by the 29,903 base pairs of RNA that codes for the SARS-CoV-2 genome. SARS-CoV-2 infects humans to cause COVID-19, spreading from patient-to-patient yet impacts patients very divergently.Areas covered: Within this review, we address the known molecular mechanisms and supporting data for COVID-19 clinical course and pathology, clinical risk factors and molecular signatures, therapeutics of severe COVID-19, and reinfection/vaccination. Literature and published datasets were reviewed using PubMed, Google Scholar, and NCBI SRA tools. The combination of exaggerated cytokine signaling, pneumonia, NETosis, pyroptosis, thrombocytopathy, endotheliopathy, multiple organ dysfunction syndrome (MODS), and acute respiratory distress syndrome (ARDS) create a positive feedback loop of severe damage in patients with COVID-19 that impacts the entire body and may persist for months following infection. Understanding the molecular pathways of severe COVID-19 opens the door for novel therapeutic design. We summarize the current insights into pathology, risk factors, secondary infections, genetics, omics, and drugs being tested to treat severe COVID-19.Expert opinion: A growing level of support suggests the need for stronger integration of biomarkers and precision medicine to guide treatment strategies of severe COVID-19, where each patient has unique outcomes and thus require guided treatment.

Progress and challenge for computational quantification of tissue immune cells

Tissue immune cells have long been recognized as important regulators for the maintenance of balance in the body system. Quantification of the abundance of different immune cells will provide enhanced understanding of the correlation between immune cells and normal or abnormal situations. Currently, computational methods to predict tissue immune cell compositions from bulk transcriptomes have been largely developed. Therefore, summarizing the advantages and disadvantages is appropriate. In addition, an examination of the challenges and possible solutions for these computational models will assist the development of this field. The common hypothesis of these models is that the expression of signature genes for immune cell types might represent the proportion of immune cells that contribute to the tissue transcriptome. In general, we grouped all reported tools into three groups, including reference-free, reference-based scoring and reference-based deconvolution methods. In this review, a summary of all the currently reported computational immune cell quantification tools and their applications, limitations, and perspectives are presented. Furthermore, some critical problems are found that have limited the performance and application of these models, including inadequate immune cell type, the collinearity problem, the impact of the tissue environment on the immune cell expression level, and the deficiency of standard datasets for model validation. To address these issues, tissue specific training datasets that include all known immune cells, a hierarchical computational framework, and benchmark datasets including both tissue expression profiles and the abundances of all the immune cells are proposed to further promote the development of this field.

SARS-CoV-2-Encoded Proteome and Human Genetics: From Interaction-Based to Ribosomal Biology Impact on Disease and Risk Processes

SARS-CoV-2 (COVID-19) has infected millions of people worldwide, with lethality in hundreds of thousands. The rapid publication of information, both regarding the clinical course and the viral biology, has yielded incredible knowledge of the virus. In this review, we address the insights gained for the SARS-CoV-2 proteome, which we have integrated into the Viral Integrated Structural Evolution Dynamic Database, a publicly available resource. Integrating evolutionary, structural, and interaction data with human proteins, we present how the SARS-CoV-2 proteome interacts with human disorders and risk factors ranging from cytokine storm, hyperferritinemic septic, coagulopathic, cardiac, immune, and rare disease-based genetics. The most noteworthy human genetic potential of SARS-CoV-2 is that of the nucleocapsid protein, where it is known to contribute to the inhibition of the biological process known as nonsense-mediated decay. This inhibition has the potential to not only regulate about 10% of all biological transcripts through altered ribosomal biology but also associate with viral-induced genetics, where suppressed human variants are activated to drive dominant, negative outcomes within cells. As we understand more of the dynamic and complex biological pathways that the proteome of SARS-CoV-2 utilizes for entry into cells, for replication, and for release from human cells, we can understand more risk factors for severe/lethal outcomes in patients and novel pharmaceutical interventions that may mitigate future pandemics.