| Literature DB >> 32778178 |
Margarita Villar1,2, Isabel G Fernández de Mera1, Sara Artigas-Jerónimo1, Marinela Contreras1,3, Christian Gortázar1, José de la Fuente4,5.
Abstract
The coronavirus disease 19 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has affected millions of people worldwide. Recent evidence raised the question about the possibility that cats may be a domestic host for SARS-CoV-2 with unknown implications in disease dissemination. Based on the fact that the domestic cat flea, Ctenocephalides felis, are abundant ectoparasites infesting humans, companion animals and wildlife and that coronavirus-like agents have been identified in the ectoparasite tick vector, Ixodes uriae of seabirds, herein we considered the presence of coronaviruses in general and SARS-CoV-2 in particular in C. felis. We identified coronavirus-derived and cell receptor angiotensin-converting enzyme RNA/proteins in C. felis. Although current evidence suggests that pets are probably dead-end-hosts with small risk of transmission to humans, our results suggested that cat flea may act as biological and/or mechanical vectors of SARS-CoV. Although preliminary, these results indicate a possibility of ectoparasites acting as reservoirs and vectors of SARS-CoV and related beta-coronavirus although with little disease risk due to systemic transmission route, low viremia, virus attenuation or other unknown factors. These results support the need to further study the role of animal SARS-CoV-2 hosts and their ectoparasite vectors in COVID-19 disease spread.Entities:
Keywords: ACE; Arthropod; COVID-19; Cat; Coronavirus; Flea; Proteomics
Mesh:
Substances:
Year: 2020 PMID: 32778178 PMCID: PMC7416815 DOI: 10.1186/s13071-020-04292-y
Source DB: PubMed Journal: Parasit Vectors ISSN: 1756-3305 Impact factor: 3.876
Fig. 1Experimental design and identification of coronavirus-derived RNA and proteins in cat flea. Representative images of a laboratory-reared domestic cat flea and tissues used for extraction of RNA and plasma membrane proteins for RT-PCR and proteomics analysis. Schematic representation of the coronavirus genome organization and virion structure based on SARS-CoV-2 [15]. Coronavirus proteins identified by proteomics analysis included ORF1a, ORF1b, protein M and protein 7b. Genomic RNA (gRNA) serves as mRNA for ORF1a and ORF1b. Other major subgenomic RNAs (sgRNAs) are produced to encode for envelope (e.g. protein M) and accessory (e.g. protein 7b) proteins in addition to the gRNA. Real-time RT-PCR targeting ORF1b identified in coronavirus-derived RNA in laboratory-reared unfed cat flea (sample 2) and feral cat flea (sample 3). MW, molecular weight O’GeneRuler 1 kb Plus DNA Ladder (Thermo Fisher Scientific); sample 1, laboratory-reared unfed cat flea, samples 4 and 5, positive controls (Table 1); samples 6 and 7, nuclease-free water negative controls
RT-PCR targets, oligonucleotide primers and results
| RT-PCR target | Primers: 5′–3′ sequences (amplicon size) | Results (flea sample: Ct) |
|---|---|---|
| Coronavirus generic group ORF1b | 11-FW: TGATGATGSNGTTGTNTGYTAYAA 13-RV: GCATWGTRTGYTGNGARCARAATTC (179 bp) | Fed: na Unfed: 35.24 Feral cat 1: 36.16 Feral cat 2: na |
| SARS-CoV-2, RdRp-IP2 | nCoV_IP2-12669Fw: ATGAGCTTAGTCCTGTTG nCoV_IP2-12759Rv: CTCCCTTTGTTGTGTTGT nCoV_IP2-12696b: AGATGTCTTGTGCTGCCGGTA [5′]Hex [3′]BHQ-1 (108 bp) | Fed: na Unfed: na Feral cat 1: na Feral cat 2: na |
| SARS-CoV-2, RdRp-IP4 | nCoV_IP4-14059Fw: GGTAACTGGTATGATTTCG nCoV_IP4-14146Rv: CTGGTCAAGGTTAATATAGG nCoV_IP4-14084: TCATACAAACCACGCCAGG [5′]Fam [3′]BHQ-1 (107 bp) | Fed: na Unfed: na Feral cat 1: na Feral cat 2: na |
| SARS-CoV-2, E-gene | E_Sarbeco_F1: ACAGGTACGTTAATAGTTAATAGCGT E_Sarbeco_R2: ATATTGCAGCAGTACGCACACA E_Sarbeco_P1: ACACTAGCCATCCTTACTGCGCTTCG [5′]Fam [3′]BHQ-1 (125 bp) | Fed: na Unfed: na Feral cat 1: na Feral cat 2: na |
Notes: RNA samples were extracted using the AllPrep DNA/RNA/Protein Mini Kit (Qiagen, Valencia, CA, USA) from midguts dissected from laboratory-reared unfed and fed cat fleas [3] and from the pools of fleas collected from feral cats. A SYBR green One-Step real-time RT-PCR assay (BioRad, Hercules, CA, USA) targeting the ORF1b was used for generic detection of coronaviruses [16]. Three SARS-CoV-2-specific RT-PCRs targeting the envelop protein E-coding gene and two targets (IP2 and IP4) of RNA-dependent RNA polymerase gene (RdRp) were conducted using the SuperScript III Platinum One-Step qRT-PCR kit (Thermo Fisher Scientific). The protocols used for RT-PCR are included in the WHO guidelines (https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/laboratory-guidance). The positive controls included a positive sample loaned from the University Hospital of Ciudad Real, Spain and an in vitro transcribed RNA derived from the strain BetaCoV_Wuhan_WIV04_2019 (EPI_ISL_402124) loaned by the Pasteur Institute, Paris, France. Real-time RT-PCR was carried out using the CFX96 Touch Real-Time PCR Detection System Thermal Cycler (BioRad, Hercules, CA, USA). Flea samples: Fed and Unfed, laboratory-reared cat flea; Feral cat, pool of feral cat fleas
Abbreviations: na, not amplified; Ct, cycle threshold
Fig. 2Coronavirus-derived sequences with the peptides identified by proteomics analysis and protein coverage. Proteins were analyzed by reverse phase (RP)-liquid chromatography (LC)-mass spectrometry (MS)/MS (RP-LC-MS/MS) using an Easy-nLC II system coupled to an ion trap LCQ Fleet mass spectrometer (Thermo Fisher Scientific) as previously reported [3]. For this study, MS/MS raw files were searched again against a compiled database containing all proteins from Felidae, coronavirus, Drosophila and Ctenocephalidae (78265, 30637, 342897 and 475 Uniprot entries in May 2020, respectively). Amino acid sequences of proteins corresponding to the Uniprot (https://www.uniprot.org) entries shown in parenthesis. The sequences of the peptides identified by proteomics analysis are highlighted in green and the sum of their amino acids reflect protein coverage when compared to the sequence of the identified proteins
Fig. 3Evolutionary conservation of SARS-CoV receptor ACE protein. a Amino acid sequence alignment (Blast E-value = 0.003, max cluster distance = 0.4) was performed with COBALT (https://www.ncbi.nlm.nih.gov/tools/cobalt/cobalt.cgi?LINK_LOC=BlastHomeLink) using protein sequences for ACE in fruit fly (Drosophila melanogaster; Uniprot ID: NP_001285915), water flea (Daphnia pulex; E9GU43), water flea order (D. pulex; A0A162PAD4), body louse (Pediculus humanus corporis; E0VAB8), deer tick (Ixodes scapularis; A0A4D5RPS5), Chinese horseshoe bat (Rhinolophus sinicus; E2DHI7) and human ACE2 (Homo sapiens; Q9BYF1). Conserved regions are highlighted in red. b Prediction of fruit fly ACE and human ACE2 proteins secondary structure using CFSSP: Chou & Fasman Secondary Structure Prediction Server (http://www.biogem.org/tool/chou-fasman/index.php). c Amino acids K31, E35, D38, M82 and K353 identified as involved in the interface between SARS-CoV and human ACE2 [16]. d Slanted cladogram of ACE protein sequences using the Neighbor joining algorithm (max seq difference = 0.85, distance = Grishin protein) at NCBI tree viewer (https://www.ncbi.nlm.nih.gov/tools/treeviewer/)