Dual RNA-Seq analysis of SARS-CoV-2 correlates specific human transcriptional response pathways directly to viral expression.

The SARS-CoV-2 pandemic has challenged humankind's ability to quickly determine the cascade of health effects caused by a novel infection. Even with the unprecedented speed at which vaccines were developed and introduced into society, identifying therapeutic interventions and drug targets for patients infected with the virus remains important as new strains of the virus evolve, or future coronaviruses may emerge that are resistant to current vaccines. The application of transcriptomic RNA sequencing of infected samples may shed new light on the pathways involved in viral mechanisms and host responses. We describe the application of the previously developed "dual RNA-seq" approach to investigate, for the first time, the co-regulation between the human and SARS-CoV-2 transcriptomes. Together with differential expression analysis, we describe the tissue specificity of SARS-CoV-2 expression, an inferred lipopolysaccharide response, and co-regulation of CXCL's, SPRR's, S100's with SARS-CoV-2 expression. Lipopolysaccharide response pathways in particular offer promise for future therapeutic research and the prospect of subgrouping patients based on chemokine expression that may help explain the vastly different reactions patients have to infection. Taken together these findings highlight unappreciated SARS-CoV-2 expression signatures and emphasize new considerations and mechanisms for SARS-CoV-2 therapeutic intervention.

The UCSC Genome Browser database: 2021 update.

For more than two decades, the UCSC Genome Browser database (https://genome.ucsc.edu) has provided high-quality genomics data visualization and genome annotations to the research community. As the field of genomics grows and more data become available, new modes of display are required to accommodate new technologies. New features released this past year include a Hi-C heatmap display, a phased family trio display for VCF files, and various track visualization improvements. Striving to keep data up-to-date, new updates to gene annotations include GENCODE Genes, NCBI RefSeq Genes, and Ensembl Genes. New data tracks added for human and mouse genomes include the ENCODE registry of candidate cis-regulatory elements, promoters from the Eukaryotic Promoter Database, and NCBI RefSeq Select and Matched Annotation from NCBI and EMBL-EBI (MANE). Within weeks of learning about the outbreak of coronavirus, UCSC released a genome browser, with detailed annotation tracks, for the SARS-CoV-2 RNA reference assembly.

Genetic pathways disrupted by ENPP1 deficiency provide insight into mechanisms of osteoporosis, osteomalacia, and paradoxical mineralization.

Ectonucleotide phosphatase/phosphodiesterase 1 (ENPP1) deficiency results in either lethal arterial calcifications ('Generalized Arterial Calcification of Infancy' - GACI), phosphate wasting rickets ('Autosomal Recessive Hypophosphatemic Rickets type 2' - ARHR2), early onset osteoporosis, or progressive spinal rigidity ('Ossification of the Posterior Longitudinal Ligament' - OPLL). As ENPP1 generates a strong endogenous mineralization inhibitor - extracellular pyrophosphate (PPi) - ENPP1 deficiency should not result in reduced bone volume, and therefore the mechanism ENPP1 associated osteoporosis is not apparent given current understanding of the enzyme's function. To investigate genetic pathways driving the skeletal phenotype of ENPP1 deficiency we compared gene expression in Enpp1asj/asj mice and WT sibling pairs by RNAseq and qPCR in whole bones, and in the liver and kidney by qPCR, directly correlating gene expression with measures of bone microarchitectural and biomechanical phenotypes. Unbiased analysis of the differentially expressed genes compared to relevant human disease phenotypes revealed that Enpp1asj/asj mice exhibit strong signatures of osteoporosis, ARHR2 and OPLL. We found that ENPP1 deficient mice exhibited reduced gene transcription of Wnt ligands in whole bone and increased transcription of soluble Wnt inhibitors in the liver and kidney, suggestive of multiorgan inhibition of Wnt activity. Consistent with Wnt suppression in bone, Collagen gene pathways in bone were significantly decreased and Fgf23 was significantly increased, all of which directly correlated with bone microarchitectural defects and fracture risk in Enpp1asj/asj mice. Moreover, the bone findings in 10-week old mice correlated with Enpp1 transcript counts but not plasma [PPi], suggesting that the skeletal phenotype at 10 weeks is driven by catalytically independent ENPP1 function. In contrast, the bone findings in 23-week Enpp1asj/asj mice strongly correlated with plasma PPi, suggesting that chronically low PPi drives the skeletal phenotype in older mice. Finally, correlation between Enpp1 and Fgf23 transcription suggested ENPP1 regulation of Fgf23, which we confirmed by dosing Enpp1asj/asj mice with soluble ENPP1-Fc and observing suppression of intact plasma FGF23 and ALP. In summary, our findings suggest that osteoporosis associated with ENPP1 deficiency involves the suppression of Wnt via catalytically independent Enpp1 pathways, and validates Enpp1asj/asj mice as tools to better understand OPLL and Paradoxical Mineralization Disorders.

Human Heterozygous ENPP1 Deficiency Is Associated With Early Onset Osteoporosis, a Phenotype Recapitulated in a Mouse Model of Enpp1 Deficiency.

Biallelic ENPP1 deficiency in humans induces generalized arterial calcification of infancy (GACI) and/or autosomal recessive hypophosphatemic rickets type 2 (ARHR2). The latter is characterized by markedly increased circulating FGF23 levels and renal phosphate wasting, but aberrant skeletal manifestations associated with heterozygous ENPP1 deficiency are unknown. Here, we report three adult men with early onset osteoporosis who presented with fractures in the thoracic spine and/or left radius, mildly elevated circulating FGF23, and hypophosphatemia. Total hip bone mineral density scans demonstrated osteoporosis (Z-score < -2.5) and HRpQCT demonstrated microarchitectural defects in trabecular and cortical bone. Next-generation sequencing revealed heterozygous loss-of-function mutations in ENPP1 previously observed as biallelic mutations in infants with GACI. In addition, we present bone mass and structure data as well as plasma pyrophosphate (PPi) data of two siblings suffering from ARHR2 in comparison to their heterozygous and wild-type family members indicative of an ENPP1 gene dose effect. The skeletal phenotype in murine Enpp1 deficiency yielded nearly identical findings. Ten-week-old male Enpp1 asj/asj mice exhibited mild elevations in plasma FGF23 and hypophosphatemia, and micro-CT analysis revealed microarchitectural defects in trabecular and cortical bone of similar magnitude to HRpQCT defects observed in humans. Histomorphometry revealed mild osteomalacia and osteopenia at both 10 and 23 weeks. The biomechanical relevance of these findings was demonstrated by increased bone fragility and ductility in Enpp1 asj/asj mice. In summary, ENPP1 exerts a gene dose effect such that humans with heterozygous ENPP1 deficiency exhibit intermediate levels of plasma analytes associated with bone mineralization disturbance resulting in early onset osteoporosis. © 2019 The Authors. Journal of Bone and Mineral Research published by American Society for Bone and Mineral Research.