Traces of SARS-CoV-2 RNA in the Blood of COVID-19 Patients
==========================================================

* Ahmed Moustafa
* Ramy K. Aziz

## Abstract

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third virus that caused coronavirus-related outbreaks over the past 20 years. The outbreak was first reported in December 2019 from Wuhan, China, but rapidly progressed into a pandemic of an unprecedented scale since the 1918 flu pandemic. In addition to respiratory complications in COVID-19 patients, clinical characterizations of severe infection cases indicated a number of other comorbidities, including multiple organ failure (liver, kidney, and heart) and septic shock. In an attempt to elucidate COVID-19 pathogenesis in different human organs, we interrogated the presence of the virus in the blood, or any of its components, which might provide a form of trafficking or hiding to the virus. By computationally analyzing high-throughput sequence data from patients with active COVID-19 infection, we found evidence of traces of SARS-CoV-2 RNA in peripheral blood mononuclear cells (PBMC), while the virus RNA was abundant in bronchoalveolar lavage specimens from the same patients. To the best of our knowledge, the presence of SAR-CoV-2 RNA in the PBMC of COVID-19 patients has not been reported before, and this observation could suggest immune presentation, but discounts the possibility of extensive viral infection of lymphocytes or monocytes.

Keywords
*   SARS-CoV-2
*   Coronavirus
*   COVID-19
*   Blood
*   PBMC
*   RNA-Seq

## Introduction

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third virus that caused coronavirus-related outbreaks over the past 20 years. The first outbreak occurred in Asia in 2002-2003 causing Severe Acute Respiratory Syndrome (SARS), hence, the name, SARS-CoV, which back then was not related to any of the known viruses (Marra et al. 2003; Rota et al. 2003). Between 2002 and 2003, 8,098 people became sick with SARS, and of those 774 died (i.e., a mortality rate of 9.5%). Since 2004, there have been no more reports of SARS cases (via [NHS](https://www.nhs.uk/conditions/sars/), [WHO](https://www.cdc.gov/sars/about/faq.html), and [CDC](https://www.cdc.gov/sars/about/faq.html)).

The second coronavirus-related outbreak started in the Arabian Peninsula in 2012 (Zaki et al. 2012) causing a more fetal disease, Middle East respiratory syndrome (MERS), with a significantly higher mortality rate of 40% of the cases infected by MERS-CoV virus (Zumla, Hui, and Perlman 2015).

More recently, in December 2019, the third coronavirus-related outbreak was first reported in Wuhan, China, by the emergence of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), initially dubbed “the 2019 novel coronavirus” (2019-nCoV). The spread of the virus led to a pandemic of an unprecedented scale since the 1918 flu pandemic.

As of May 10, almost 4 million confirmed COVID-19 cases globally and 27,000 deaths have been reported by WHO ([WHO Dashboard](https://covid19.who.int/), continuously updated). In addition to the respiratory complications in COVID-19 patients, clinical characterizations of severe infection cases indicated further comorbidities, including multiple organ failure (liver, kidney, and heart) and septic shock (Poston, Patel, and Davis 2020; Cascella et al. 2020; Hui Li et al. 2020).

The genome sequence of SARS-COV-2 has been determined and made public (Lu et al. 2020), and since then thousands (more than 17,000) of genomes have been sequenced from all around the world (Shu and McCauley 2017). The availability of those genomic sequences allows the rapid screening of viral RNA in human tissues as well as environmental samples (e.g., sewage (Bibby and Peccia 2013)) using multi-omic wet lab technologies, as well as in silico screening tools, for publicly available metatranscriptomic samples.

In an attempt to elucidate COVID-19 pathogenesis in different human organs, we conducted this study to interrogate the presence of the virus in the blood, or any of its components, which might provide a form of trafficking or hiding to the virus, notably that some precarious studies reported the ability of the virus to infect lymphocytes (X. Wang et al. 2020) whereas others went as far as to suggest the virus exerts its pathogenesis through “attacking hemoglobin,” although this hypothesis has been heavily criticized (Read 2020). Moreover, the virus was reported to be found in the plasma of COVID-19 patients (Huang et al. 2020). Finally, peripheral blood mononuclear cells (PBMC) were shown to harbor other infectious viruses, such as HIV, HCV, and HBV (W.-K. Wang et al. 2002; Z. Li, Hou, and Cao 2015).

For the aforementioned reasons, we computationally analyzed high-throughput sequence data from patients with active COVID-19 infection, and found evidence of traces of SARS-CoV-2 RNA in their PMBCs, while their bronchoalveolar lavage samples had large amounts of the viral RNA.

## Methods

### Dataset

Publicly accessible raw RNA-Seq FASTQ sequences published by Xiong et al (Xiong et al. 2020) were obtained from the Genome Sequence Archive (Y. Wang et al. 2017) (GSA accession [CRA002390](https://bigd.big.ac.cn/gsa/browse/CRA002390)).

### Comparison of RNA Abundance and Gene Expression Profiles

Filtered RNA-Seq reads were searched (with blastx (Altschul et al. 1990)) against the RefSeq protein database (O’Leary et al. 2016) (release 99) using DIAMOND (Buchfink, Xie, and Huson 2015) with an e-value cutoff < 1e-10. The counts of the matching RefSeq accessions were normalized to the total number of reads per library. Principal Component Analysis (PCA) was performed on the normalized abundance of the RefSeq proteins.

### Detection of Viral RNA

Raw sequences were processed for quality control using fastp (Chen et al. 2018). Filtered FASTQ sequences were aligned to the SAR-CoV-2 reference genome (GenBank accession [NC_045512](https://www.ncbi.nlm.nih.gov/nuccore/1798174254)) using BWA (Heng Li and Durbin 2009). Generated BAM files were filtered for mapped sequences with a quality score > 40 and alignment score > 90 using Sambamba (Tarasov et al. 2015). Identified SAR-CoV-2 matching sequences were manually inspected and searched against NCBI “nt” using BLAST (blastn) for verification (Altschul et al. 1990).

## Results

We analyzed the RNA-Seq dataset published by Xiong et al (Xiong et al. 2020). **Table 1** describes the processed samples and the sequence reads statistics of each sample:

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/05/18/2020.05.10.20097055/T1)

Table 1. Summary of the samples included in the study
Nine RNA-Seq samples (Xiong et al. 2020) were analyzed, including 3 PBMC from healthy donors, 3 PBMC from patients, and 2 BALF from patients with an additional replicate each. The numbers are paired-end reads.

Based on the global gene expression in PBMC, PCA showed a strong, expected separation between BALF and PBMC samples, and further slight separation in the PBMC from healthy controls and COVID-19 patients (**Figure 1**). However, one healthy control was placed within the COVID-19 cluster.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/18/2020.05.10.20097055/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/05/18/2020.05.10.20097055/F1)

Figure 1. Principal Component Analysis based on RNA abundance data
Samples are color-coded, red for healthy donors, and green for COVID-19 patients. PC1 explains 63% of the variance and PC2 explains 30% of the variance.

We found viral sequences in all of the BALF samples with a median of abundance 2.15% of the total reads. We also identified 2 paired-ended viral reads in a PBMC sample of a COVID-19 patient, matching the SAR-CoV-2 polyprotein (pp1ab) (accession [NP_828849](https://www.ncbi.nlm.nih.gov/protein/30124074)) and SAR-CoV-2 surface glycoprotein (accession [YP_009724390](https://www.ncbi.nlm.nih.gov/protein/1796318598)) (**Figures 2, 3**).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/18/2020.05.10.20097055/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/05/18/2020.05.10.20097055/F2)

Figure 2. Mapping of SARS-CoV-2 RNA against the SARS-CoV-2 genome

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/18/2020.05.10.20097055/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/05/18/2020.05.10.20097055/F3)

Figure 3. Alignment of the detected viral RNA sequences to SARS-CoV-2 genome

Viral RNA sequences from the COVID-19 patients were mapped against the SARS-CoV-2 reference genome (GenBank accession [NC_045512](https://www.ncbi.nlm.nih.gov/nuccore/1798174254)). The x-axis is the nucleotide position on the virus genome. The y-axis is the coverage (in reads; not normalized) of the genomic position by RNA-Seq viral sequences. Panel "A” represents the BALF samples. Panel "B” represents the PBMC samples.

The two paired-end viral reads from a COVID-19 PBMC sample aligned to the SARS-CoV-2 reference genome (GenBank accession [NC_045512](https://www.ncbi.nlm.nih.gov/nuccore/1798174254)). Panel A is the alignment first paired-end read. Panel B is the alignment second paired-end read. Perfect identities between the reference and the two overlapping mates are indicated by “*”.

## Discussion

Coronavirus-related infections have been reported to be associated with hematological changes including lymphopenia, thrombocytopenia, and leukopenia through infecting blood cells, bone marrow stromal cells, or inducing autoantibodies (Yang et al. 2003).

In an early study characterizing the clinical features of COVID-19 patients, Huang and coworkers indicated that using RT-PCR allowed them to detect coronavirus in plasma isolated samples from the patients (Huang et al. 2020). In their report, they preferred to use the term "RNAaemia,” rather than "viraemia,” which they defined as the presence of a virus in the blood because they did not perform tests to confirm the presence of an infectious SAR-CoV-2 virus in the blood of the patients.

On the other hand, there have been no reports of detecting viral RNA in blood cells, notably PBMC. On the contrary, in their preliminary analysis of RNA isolated from PBMC, [Corley et al](https://doi.org/10.1101/2020.04.13.039263) confirmed that they did not detect viral sequences (currently, a preprint [10.1101/2020.04.13.039263v1](https://www.biorxiv.org/content/10.1101/2020.04.13.039263v1)).

Using high-throughput sequencing has been repeatedly demonstrated to be an effective approach for the identification and quantification of viruses in the blood (Moustafa et al. 2017) following similar methods as those used in viral metagenomics and uncultivated viral genomics Breitbart 2003 (Breitbart et al. 2003), Aziz et al (Aziz et al. 2015) and Roux et al (Breitbart et al. 2003).

Therefore, we planned to exhaustively look for SAR-CoV-2 RNA sequences in publicly available RNA-Seq PBMC datasets. As of the writing of this report, there has been only one publicly available RNA-Seq PBMC dataset published by (Xiong et al. 2020), in which the group profiled global gene expression in BALF and PBMC specimens of COVID-19. Predictably, we detected viral RNA in all BALF samples (2 patients, 2 replicates each) with an average abundance of 2.15% of the total RNA in those samples, which included human RNA. However, we identified two paired-end reads in RNA isolated from PBMC from only one (accession [CRR119891](https://bigd.big.ac.cn/gsa/browse/CRA002390/CRR119891)) out of three patients, which aligned to the SAR-CoV-2 genome. Expectedly, no SARS-CoV-2 RNA was detected in the PBMC from healthy controls (3 donors). These RNA traces are certainly quite rare; however, they confidently and specifically belong to SARS-COV-2. One viral RNA read translates into polyprotein (pp1ab, accession [NP_828849](https://www.ncbi.nlm.nih.gov/protein/30124074)), which is the largest protein of coronaviruses and involved in the replication and transcription of the viral genome. The other viral RNA translates into surface (spike) glycoprotein (accession [YP_009724390](https://www.ncbi.nlm.nih.gov/protein/1796318598)), which mediates the entry of the viral into the human cells expressing human angiotensin-converting enzyme 2 (hACE2) (Ou et al. 2020; Walls et al. 2020).

Although we are not rejecting the possibility of cross-contamination or barcode bleeding (Mitra et al. 2015; Kircher, Sawyer, and Meyer 2012) for detecting the viral RNA in one PBMC RNA-Seq sample, such possibility is unlikely, given that control samples had zero hits to SARS-CoV-2 RNA. We are also considering the possibility of SAR-CoV-2 being sampled by antigen-presenting cells (most likely dendritic cells) or presented to T lymphocytes, which are in the PBMC population.

One more possibility, which we believe needs many more samples to consider, is that SARS-CoV-2 may be specifically or coincidentally internalized by one of the mononuclear cell types, which may suggest a mechanism for the chronicity of the SAR-CoV-2 infection. This hypothesis requires further testing. However, in the light of our data, it is hard to support the early reports that SARS-COV-2 targets T-lymphocytes *in vivo*, as suggested earlier, in a correspondence, based on cell culture experiment with pseudotyped viruses (X. Wang et al. 2020).

With more data becoming publicly available, it will be possible to revisit this hypothesis and others to improve our understanding of the progression and replication of SAR-CoV-2 in infected individuals.

## Data Availability

The data was published by Xiong et al. 2020 and is publicly available through the Genome Sequence Archive with accession CRA002390.

[https://bigd.big.ac.cn/gsa/browse/CRA002390](https://bigd.big.ac.cn/gsa/browse/CRA002390) 

*   Received May 10, 2020.
*   Revision received May 10, 2020.
*   Accepted May 18, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  Altschul, S. F.,  W. Gish,  W. Miller,  E. W. Myers, and  D. J. Lipman. 1990. “Basic Local Alignment Search Tool.” Journal of Molecular Biology 215 (3): 403–10.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/jmbi.1990.9999&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2231712&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F18%2F2020.05.10.20097055.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990ED16700008&link_type=ISI) 

2.  Aziz, Ramy K.,  Bhakti Dwivedi,  Sajia Akhter,  Mya Breitbart, and  Robert A. Edwards. 2015. “Multidimensional Metrics for Estimating Phage Abundance, Distribution, Gene Density, and Sequence Coverage in Metagenomes.” Frontiers in Microbiology 6 (May): 381.
    
    

3.  Bibby, Kyle, and  Jordan Peccia. 2013. “Identification of Viral Pathogen Diversity in Sewage Sludge by Metagenome Analysis.” Environmental Science & Technology 47 (4): 1945–51.
    
    

4.  Breitbart, Mya,  Ian Hewson,  Ben Felts,  Joseph M. Mahaffy,  James Nulton,  Peter Salamon, and  Forest Rohwer. 2003. “Metagenomic Analyses of an Uncultured Viral Community from Human Feces.” Journal of Bacteriology 185 (20): 6220–23.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MjoiamIiO3M6NToicmVzaWQiO3M6MTE6IjE4NS8yMC82MjIwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMTgvMjAyMC4wNS4xMC4yMDA5NzA1NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

5.  Buchfink, Benjamin,  Chao Xie, and  Daniel H. Huson. 2015. “Fast and Sensitive Protein Alignment Using DIAMOND.” Nature Methods 12 (1): 59–60.
    
    

6.  Cascella, Marco,  Michael Rajnik,  Arturo Cuomo,  Scott C. Dulebohn, and  Raffaela Di Napoli. 2020. “Features, Evaluation and Treatment Coronavirus (COVID-19).” In StatPearls. Treasure Island (FL): StatPearls Publishing.
    
    

7.  Chen, Shifu,  Yanqing Zhou,  Yaru Chen, and  Jia Gu. 2018. “Fastp: An Ultra-Fast All-in-One FASTQ Preprocessor.” Bioinformatics 34 (17): i884–90.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bty560&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30423086&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F18%2F2020.05.10.20097055.atom) 

8.  Huang, Chaolin,  Yeming Wang, Xingwang Li,  Lili Ren,  Jianping Zhao,  Yi Hu,  Li Zhang, et al. 2020. “Clinical Features of Patients Infected with 2019 Novel Coronavirus in Wuhan, China.” The Lancet 395 (10223): 497–506.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F18%2F2020.05.10.20097055.atom) 

9.  Kircher, Martin,  Susanna Sawyer, and  Matthias Meyer. 2012. “Double Indexing Overcomes Inaccuracies in Multiplex Sequencing on the Illumina Platform.” Nucleic Acids Research 40 (1): e3.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkr771&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22021376&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F18%2F2020.05.10.20097055.atom) 

10. Li, Heng, and  Richard Durbin. 2009. “Fast and Accurate Short Read Alignment with Burrows-Wheeler Transform.” Bioinformatics 25 (14): 1754–60.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btp324&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19451168&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F18%2F2020.05.10.20097055.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000267665900006&link_type=ISI) 

11. Li, Hui,  Liang Liu,  Dingyu Zhang,  Jiuyang Xu,  Huaping Dai,  Nan Tang,  Xiao Su, and  Bin Cao. 2020. “SARS-CoV-2 and Viral Sepsis: Observations and Hypotheses.” The Lancet, April. [https://doi.org/10.1016/S0140-6736(20)30920-X](https://doi.org/10.1016/S0140-6736(20)30920-X).
    
    

12. Li, Zixiong,  Xiaomei Hou, and  Guangwen Cao. 2015. “Is Mother-to-Infant Transmission the Most Important Factor for Persistent HBV Infection?” Emerging Microbes & Infections 4 (5): e30.
    
    

13. Lu, Roujian,  Xiang Zhao,  Juan Li,  Peihua Niu,  Bo Yang,  Honglong Wu,  Wenling Wang, et al. 2020. “Genomic Characterisation and Epidemiology of 2019 Novel Coronavirus: Implications for Virus Origins and Receptor Binding.” The Lancet 395 (10224): 565–74.
    
    

14. Marra, Marco A.,  Steven J. M. Jones,  Caroline R. Astell,  Robert A. Holt,  Angela Brooks-Wilson,  Yaron S. N. Butterfield,  Jaswinder Khattra, et al. 2003. “The Genome Sequence of the SARS-Associated Coronavirus.” Science 300 (5624): 1399–1404.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMDAvNTYyNC8xMzk5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMTgvMjAyMC4wNS4xMC4yMDA5NzA1NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

15. Mitra, Abhishek,  Magdalena Skrzypczak,  Krzysztof Ginalski, and  Maga Rowicka. 2015. “Strategies for Achieving High Sequencing Accuracy for Low Diversity Samples and Avoiding Sample Bleeding Using Illumina Platform.” PloS One 10 (4): e0120520.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0120520&link_type=DOI) 

16. Moustafa, Ahmed,  Chao Xie,  Ewen Kirkness,  William Biggs,  Emily Wong,  Yaron Turpaz,  Kenneth Bloom, et al. 2017. “The Blood DNA Virome in 8,000 Humans.” PLoS Pathogens 13 (3): e1006292.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1006292&link_type=DOI) 

17. O’Leary, Nuala A.,  Mathew W. Wright,  J. Rodney Brister,  Stacy Ciufo,  Diana Haddad,  Rich McVeigh,  Bhanu Rajput, et al. 2016. “Reference Sequence (RefSeq) Database at NCBI: Current Status, Taxonomic Expansion, and Functional Annotation.” Nucleic Acids Research 44 (D1): D733–45.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkv1189&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26553804&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F18%2F2020.05.10.20097055.atom) 

18. Ou, Xiuyuan,  Yan Liu,  Xiaobo Lei,  Pei Li,  Dan Mi,  Lili Ren,  Li Guo, et al. 2020. “Characterization of Spike Glycoprotein of SARS-CoV-2 on Virus Entry and Its Immune Cross-Reactivity with SARS-CoV.” Nature Communications 11 (1): 1620.
    
    

19. Poston, Jason T.,  Bhakti K. Patel, and  Andrew M. Davis. 2020. “Management of Critically Ill Adults With COVID-19.” JAMA: The Journal of the American Medical Association, March. [https://doi.org/10.1001/jama.2020.4914](https://doi.org/10.1001/jama.2020.4914).
    
    

20. Read, Randy. 2020. “Flawed Methods in ‘COVID-19: Attacks the 1-Beta Chain of Hemoglobin and Captures the Porphyrin to Inhibit Human Heme Metabolism.’”
    
    

21. Rota, Paul A.,  M. Steven Oberste,  Stephan S. Monroe,  W. Allan Nix,  Ray Campagnoli,  Joseph P. Icenogle,  Silvia Peñaranda, et al. 2003. “Characterization of a Novel Coronavirus Associated with Severe Acute Respiratory Syndrome.” Science 300 (5624): 1394–99.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMDAvNTYyNC8xMzk0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMTgvMjAyMC4wNS4xMC4yMDA5NzA1NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

22. Shu, Yuelong, and  John McCauley. 2017. “GISAID: Global Initiative on Sharing All Influenza Data - from Vision to Reality.” Euro Surveillance: Bulletin Europeen Sur Les Maladies Transmissibles = European Communicable Disease Bulletin 22 (13). [https://doi.org/10.2807/1560-7917.ES.2017.22.13.30494](https://doi.org/10.2807/1560-7917.ES.2017.22.13.30494).
    
    

23. Tarasov, Artem,  Albert J. Vilella,  Edwin Cuppen,  Isaac J. Nijman, and  Pjotr Prins. 2015. “Sambamba: Fast Processing of NGS Alignment Formats.” Bioinformatics 31 (12): 2032–34.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btv098&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25697820&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F18%2F2020.05.10.20097055.atom) 

24. Walls, Alexandra C.,  Young-Jun Park,  M. Alejandra Tortorici,  Abigail Wall,  Andrew T. McGuire, and  David Veesler. 2020. “Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein.” Cell 181 (2): 281–92.e6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=doi.org/10.1016/j.cell.2020.02.058&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F18%2F2020.05.10.20097055.atom) 

25. Wang, Wei-Kung,  Tzu-Ling Sung,  Yu-Chen Tsai,  Chuan-Liang Kao,  Shu-Mei Chang, and  Chwan-Chuen King. 2002. “Detection of Dengue Virus Replication in Peripheral Blood Mononuclear Cells from Dengue Virus Type 2-Infected Patients by a Reverse Transcription-Real-Time PCR Assay.” Journal of Clinical Microbiology 40 (12): 4472–78.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjEwOiI0MC8xMi80NDcyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMTgvMjAyMC4wNS4xMC4yMDA5NzA1NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

26. Wang, Xinling,  Wei Xu,  Gaowei Hu,  Shuai Xia,  Zhiping Sun,  Zezhong Liu,  Youhua Xie,  Rong Zhang,  Shibo Jiang, and  Lu Lu. 2020. “SARS-CoV-2 Infects T Lymphocytes through Its Spike Protein-Mediated Membrane Fusion.” Cellular & Molecular Immunology, April. [https://doi.org/10.1038/s41423-020-0424-9](https://doi.org/10.1038/s41423-020-0424-9).
    
    

27. Wang, Yanqing,  Fuhai Song,  Junwei Zhu,  Sisi Zhang,  Yadong Yang,  Tingting Chen,  Bixia Tang, et al. 2017. “GSA: Genome Sequence Archive.” Genomics, Proteomics & Bioinformatics 15 (1): 14–18.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.gpb.2017.01.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28387199&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F18%2F2020.05.10.20097055.atom) 

28. Xiong, Yong,  Yuan Liu,  Liu Cao,  Dehe Wang,  Ming Guo,  Ao Jiang,  Dong Guo, et al. 2020. “Transcriptomic Characteristics of Bronchoalveolar Lavage Fluid and Peripheral Blood Mononuclear Cells in COVID-19 Patients.” Emerging Microbes & Infections 9 (1): 761–70.
    
    

29. Yang, Mo,  Kam-Lun E. Hon,  Karen Li,  Tai-Fai Fok, and  Chi-Kong Li. 2003. “The Effect of SARS Coronavirus on Blood System: Its Clinical Findings and the Pathophysiologic Hypothesis.” Zhongguo Shi Yan Xue Ye Xue Za Zhi / Zhongguo Bing Li Sheng Li Xue Hui = Journal of Experimental Hematology / Chinese Association of Pathophysiology 11 (3): 217–21.
    
    

30. Zaki, Ali M.,  Sander van Boheemen,  Theo M. Bestebroer,  Albert D. M. E. Osterhaus, and  Ron A. M. Fouchier. 2012. “Isolation of a Novel Coronavirus from a Man with Pneumonia in Saudi Arabia.” The New England Journal of Medicine 367 (19): 1814–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1211721&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23075143&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F18%2F2020.05.10.20097055.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000310773200008&link_type=ISI) 

31. Zumla, Alimuddin,  David S. Hui, and  Stanley Perlman. 2015. “Middle East Respiratory Syndrome.” The Lancet 386 (9997): 995–1007.