Next Article in Journal
Impact of Predator Exclusion and Habitat on Seroprevalence of New World Orthohantavirus Harbored by Two Sympatric Rodents within the Interior Atlantic Forest
Next Article in Special Issue
Genomic Characterization of a Novel Alphacoronavirus Isolated from Bats, Korea, 2020
Previous Article in Journal
Alveolar-like Macrophages Attenuate Respiratory Syncytial Virus Infection
Previous Article in Special Issue
Evolutionary Signatures Governing the Codon Usage Bias in Coronaviruses and Their Implications for Viruses Infecting Various Bat Species
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 13
Issue 10
10.3390/v13101962
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic Characterization of Diverse Bat Coronavirus HKU10 in Hipposideros Bats

by
Ning Wang
1,2,3,*,†,
Chu-Ming Luo
1,2,†,
Xing-Lou Yang
1,
Hai-Zhou Liu
1,
Li-Biao Zhang
4,
Wei Zhang
1,
Bei Li
1,
Yan Zhu
1,
Cheng Peng
1,
Zheng-Li Shi
1 and
Ben Hu
1,*
1
CAS Key Laboratory of Special Pathogens and Biosafety, Wuhan Institute of Virology, Chinese Academy of Sciences, Wuhan 430071, China
2
University of Chinese Academy of Sciences, Beijing 100864, China
3
Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
4
Guangdong Key Laboratory of Animal Conservation and Resource Utilization, Guangdong Public Laboratory of Wild Animal Conservation and Utilization, Guangdong Institute of Applied Biological Resources, Guangzhou 510260, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2021, 13(10), 1962; https://doi.org/10.3390/v13101962
Submission received: 30 July 2021 / Revised: 19 September 2021 / Accepted: 23 September 2021 / Published: 29 September 2021
(This article belongs to the Special Issue Bats and Coronaviruses)
Browse Figures
Versions Notes

Abstract

:
Bats have been identified as natural reservoirs of a variety of coronaviruses. They harbor at least 19 of the 33 defined species of alpha- and betacoronaviruses. Previously, the bat coronavirus HKU10 was found in two bat species of different suborders, Rousettus leschenaultia and Hipposideros pomona, in south China. However, its geographic distribution and evolution history are not fully investigated. Here, we screened this viral species by a nested reverse transcriptase PCR in our archived samples collected over 10 years from 25 provinces of China and one province of Laos. From 8004 bat fecal samples, 26 were found to be positive for bat coronavirus HKU10 (BtCoV HKU10). New habitats of BtCoV HKU10 were found in the Yunnan, Guangxi, and Hainan Provinces of China, and Louang Namtha Province in Laos. In addition to H. pomona, BtCoV HKU10 variants were found circulating in Aselliscus stoliczkanus and Hipposideros larvatus. We sequenced full-length genomes of 17 newly discovered BtCoV HKU10 strains and compared them with previously published sequences. Our results revealed a much higher genetic diversity of BtCoV HKU10, particularly in spike genes and accessory genes. Besides the two previously reported lineages, we found six novel lineages in their new habitats, three of which were located in Yunnan province. The genotypes of these viruses are closely related to sampling locations based on polyproteins, and correlated to bat species based on spike genes. Combining phylogenetic analysis, selective pressure, and molecular-clock calculation, we demonstrated that Yunnan bats harbor a gene pool of BtCoV HKU10, with H. pomona as a natural reservoir. The cell tropism test using spike-pseudotyped lentivirus system showed that BtCoV HKU10 could enter cells from human and bat, suggesting a potential interspecies spillover. Continuous studies on these bat coronaviruses will expand our understanding of the evolution and genetic diversity of coronaviruses, and provide a prewarning of potential zoonotic diseases from bats.

1. Introduction

Coronaviruses (CoVs) are enveloped, non-segmented positive-strand RNA viruses, with a genome ranging from 26–32 kb [1]. According to the nomenclature and taxonomy recently released by the International Committee for Taxonomy of Viruses (ICTV), CoVs are classified into four genera, Alpha-, Beta-, Gamma- and Deltacoronavirus, within the Orthocoronavirinae subfamily, Coronaviridae family, Cornidovirineae suborder, and Nidovirales order [2]. All coronaviruses share a similar genome organization and discontinuous transcription mechanism [3]. Two overlapping open reading frames, (ORF)1a/b, occupying 2/3 genome at the 5′ end, encode polyproteins which are cleaved into 15 or 16 nonstructural proteins (nsp1 or nsp2 through nsp16) in host cells. The structural proteins spike (S), envelope (E), membrane (M) and nucleocapsid (N) are expressed from a nested set of 3′ co-terminal subgenomic messenger RNAs (sg mRNAs).
CoVs infect various animals, including but not limited to pig, cat, bird, bat and human. While mammal-infecting CoVs can be found in all four genera, bat CoVs are considered as the genetic source of Alphacoronavirus and Betacoronavirus, and avian CoVs are considered as the genetic source of Gammacoronavirus and Deltacoronavirus. [4]. The seven human coronaviruses discovered to date all belong to Alphacoronavirus and Betacoronavirus, and five of them are phylogenetically related to, or even originate from CoVs circulating in bats [5,6,7,8,9,10,11].
CoVs are known to possess high mutation rates and recombination rates, which may allow them to adapt to new hosts and ecological niches [12]. The spike protein recognizes and binds to cellular receptor and then mediate CoV entry. It determines the host specificity of CoVs together with the cellular receptor and is essential for crossing the species barrier. Interspecies transmission of CoVs causes emerging zoonotic diseases that pose a threat to public health. Notably, three of them, including severe acute respiratory syndrome coronavirus (SARS-CoV) [13,14,15], Middle East respiratory syndrome coronavirus (MERS-CoV) [16] and SARS-CoV-2 [11], have caused severe respiratory diseases with high mortality in humans. Previous studies revealed that SARS-CoV originates from bats [5,6,17], while MERS-CoV is phylogenetically associated with the varieties of coronaviruses found in vesper bats [7,10,18], and SARS-CoV-2 likely originates in wildlife, including bats and pangolins. In addition, the swine acute diarrhea syndrome coronavirus, which caused heavy losses on pig farms in recent years, was the result of the spillover of a bat alphacoronavirus [11,19,20,21]. These data indicate that bat coronaviruses probably pose high spillover potentials.
Previous studies identified a bat coronavirus species, bat coronavirus HKU10 (BtCoV HKU10), from Guangdong Province and Hong Kong of China, and southeast Asian countries, such as Laos and Thailand [22,23,24,25]. These viruses were from two bat species belonging to different suborders, Hipposideros pomona and Rousettus leschenaultia. Analysis based on positive selection and molecular-clock calculation suggested that BtCoV HKU10 was probably a recent interspecies transmission from R. leschenaultia to H. pomona [25]. In China, Rousettus bats are distributed only in southern regions, while Hipposideros bats are widely found throughout China. In order to better understand the geographic distribution, evolution, and interspecies transmission potential of these viruses, we performed surveillance of this virus in wider geographic areas in our archived samples collected from 25 provinces of China and one province of Laos. Diverse BtCoV HKU10 strains were discovered in south and southwest China and Louang Namtha province of Laos, but not in the northern China. We sequenced the full-length genomes of 17 newly detected HKU10 strains. Phylogenetic analysis, positive selection and molecular-clock calculation demonstrated that interspecies transmission of BtCoV HKU10 may have occurred from H. pomona to R. leschenaultia, which is opposite to previous report. Our study supports the theory that H. pomona is the natural reservoir of BtCoV HKU10, and Yunnan is the geographical source of these viruses.

2. Materials and Methods

2.1. Ethics Statement

All sampling procedures were performed by veterinarians with approval from Animal Ethics Committee of the Wuhan Institute of Virology (WIVH05210201, approved on 9 July 2012). The study was conducted in accordance with the Guide for the Care and Use of Wild Mammals in Research of the People’s Republic of China.

2.2. Sampling

Bat samplings were conducted from September 2006 to June 2016, as described previously [17,26]. Bat fecal swab and pellet samples were collected at different seasons in 25 provinces in China and one province in northern Laos adjacent to Yunnan province of China.

2.3. RNA Extraction, PCR Screening and Sequencing

Viral RNA was extracted from 200 μL of fecal swab or pellet samples with High Pure Viral RNA Kit (Roche Diagnostics GmbH, Mannheim, Germany) as per the manufacturer’s instructions. RNA was eluted in 50 μL of Elution buffer, aliquoted, and stored at −80 °C. A one-step heminested RT-PCR (Invitrogen, San Diego, CA, USA) (Table S1), targeting a 440 nucleotide (nt) fragment of RNA-dependent RNA polymerase (RdRp), was employed to detect the presence of coronavirus sequences, as described previously [27]. PCR products were gel purified and sequenced with an ABI Prism 3730 DNA analyzer (Applied Biosystems, Foster City, CA, USA). Alternatively, the PCR products were cloned into pGEM-T Easy Vector (Promega, Madison, WI, USA) for sequencing. The positive samples in this study were classified using the abbreviated name of sampling location followed by sample ID (e.g., YN3723).
To confirm the bat species of an individual sample, we performed PCR to amplify the cytochrome b (Cytob) or NADH dehydrogenase subunit 1 (ND1) gene using DNA extracted from the feces or swabs [28,29]. The gene sequences were assembled excluding the primer sequences; BLASTN was used to identify host species based on the most closely related sequences with highest query coverage and a minimum identity of 95%.

2.4. Sequencing of Full-Length Genomes

Full genome sequences were determined as previously reported [26]. Briefly, sequences were determined by one-step PCR (Invitrogen, San Diego, CA, USA) amplification with degenerate primers designed on the basis of multiple alignment of available alphaCoV sequences deposited in GenBank, and then amplified with SuperScript IV Reverse Transcriptase (Invitrogen, San Diego, CA, USA) and Expand Long Template PCR System (Roche Diagnostics GmbH, Mannheim, Germany) with specific primers (Figure S1 and Table S1). For the reverse transcription, the reaction mix contained 1 μL of 2 μM gene-specific reverse primer, 1 μL of 10 mM dNTP mix, 5 μL of RNA and 6 μL of nuclease-free water, The mixture annealed at 65 °C for 5 min and then incubated on ice for 2 min. Then, we added 4 μL of 5× SSIV buffer, 1 μL of 100 mM DTT, 1 μL of RNase inhibitor and 1 μL of SuperScript IV reverse transcriptase. Mixed and incubated the mixture at 50 °C for 10 min and inactivated at 80 °C for 10 min. For long template PCT, the 50 μL reaction mix contained 5 μL PCR buffer 2, 500 μM of dNTP, 300 nM of each primer, 0.75 μL of Expand long template enzyme mix, and 5 μL of reverse transcription product. The amplification was performed as follows: 94 °C 5 min followed by 10 cycles consisting of 94 °C 10 s, 50 °C 30 s, 68 °C 10 min, then 25 cycles of 94 °C 10 s, 50 °C 30 s, 68 °C 10 min (+20 s for each cycle), and final elongation at 68 °C for 7 min. PCR products were gel purified and sequenced. Sequences of the 5′ and 3′ genomic ends were obtained by 5′ and 3′ RACE (SMARTer RACE 5′/3′ Kit; Clontech, Mountain View, CA, USA), respectively.

2.5. Genome Analysis

Preliminary sequences were assembled using DNAStar lasergene V7 (DNAStar, Madison, WI, USA). Putative open reading frames (ORFs) and deduced amino acid (aa) sequences were predicted using the NCBI’s ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder/ accessed on 29 September 2017) with a minimal ORF length of 150 nt, followed by manual inspection. The sequences of the 5′ untranslated region (5′-UTR) and 3′ untranslated region (3′-UTR) were defined, and the leader sequence, the leader and body TRSs were illustrated as previous described [30]. Phylogenetic trees based on nt or aa sequences were constructed using the Maximum Likelihood algorithm with bootstrap values determined by 1000 replicates in the MEGA 6 software package [31]. Full-length genome sequences detected in this study were aligned with those of previously reported BtCoV HKU10 using MUSCLE [31]. The aligned sequences were scanned for recombination events by using Recombination Detection Program (RDP) [32]. The potential recombination events suggested by strong p values (<10–20) were further confirmed using similarity plot and bootscan analyses implemented in Simplot 3.5.1 [33]. The number of synonymous substitutions per synonymous site, Ks, and the number of nonsynonymous substitutions per nonsynonymous site, Ka, for each coding region were calculated using the Ka/Ks calculation tool of Norwegian Bioinformatics Platform (http://services.cbu.uib.no/tools/kaks accessed on 22 October 2017) with default parameters [34].

2.6. Virus Isolation

The virus isolation was performed as previously described [5]. Briefly, Vero E6 monolayer was maintained in DMEM medium supplemented with 10% fetal bovine serum (FBS). Fecal supernatant was acquired via gradient centrifuge and then added to Vero E6 cells, and 1:10 diluted in DMEM. Cells were incubated at 37 °C for 1 h and the inoculum was then replaced by fresh DMEM containing 2% FBS and the antibiotic–antimycotic (Gibco, Grand Island, NY, USA). Cells were incubated for another 5 days and checked daily for cytopathic effect. Three passages were carried. Both culture supernatant and cell pallet were examined for CoV by RT-PCR [35].

2.7. Estimation of Divergence Time

Divergence time was calculated using RdRp gene sequence and the Bayesian Markov chain Monte Carlo (MCMC) approach as implemented in BEAST (version 2.4.4) as described previously [25]. One parametric model (constant size) and one nonparametric model (Bayesian Skyline with five groups) for tree priors were used for the inference. Analyses were performed with the SRD06 substitution model using both strict and relaxed (uncorrelated lognormal (Ucld) and uncorrelated exponential (Uced)) molecular clock. The MCMC run was 1 × 109 steps long, with sampling every 1000 steps. Convergence was assessed on the basis of the effective sampling size after a 10% run-in using Tracer software version 1.6. The mean time of the most recent common ancestor (tMRCA) and the highest posterior density regions at 95% (HPD) were calculated. Bayes Skyline under a relaxed-clock model with Uced was adopted for making inferences, as Bayes factor analysis showed it fit the data better than other models tested. The trees were summarized in a target tree by the Tree Annotator program included in the BEAST package by choosing the tree with the maximum sum of posterior probabilities (maximum clade credibility) after a 10% burn-in.

2.8. Cell Lines

Bat primary or immortalized cells (Rhinolophus sinicus kidney immortalized cells, RsKT; Rhinolophus sinicus lung primary cells, RsLu4323; Rhinolophus sinicus brain immortalized cells, RsBrT; Rhinolophus affinis kidney primary cells, RaK4324; Rousettus leschenaultia Kidney immortalized cells, RlKT; Hipposideros pratti lung immortalized cells, HpLuT) generated in our laboratory were all cultured in DMEM/F12 with 15% FBS [26]. Pteropus alecto kidney cells (Paki) was maintained in DMEM/F12 supplemented with 10% FBS [36]. Human cells (Hep-2, A549, Calu-3, H292, Caco-2, Huh-7, Hela, HEK-293T, RD, CCF-STTG1, THP-1), bat lung epithelial cells (Tb1-Lu), and other mammalian cells (LLC-ML2, Vero, BHK21, MDCK, FK, SIEC, NIH-3T3, V79) originated from American Type Culture Collection (ATCC, www.atcc.org accessed on 18 September 2021) were maintained according to the recommendations.

2.9. BtCoV HKU10 Spike-Mediated Pseudovirus Cell Tropism Screening

Retroviruses pseudotyped with the spike from Ro-BtCoV HKU10, YN7560 and MERS-CoV were used to infect human, bat or other mammalian cells in 96-well plates. The pseudovirus particles were confirmed with Western blotting and negative-staining electron microscopy. Pseudovirus containing supernatants were analyzed and transferred to PVDF membrane. HA tag antibodies (Proteintech, Wuhan, China) and P24 antibodies (Keyuananbo, Wuhan, China) were used to detect recombinant S proteins with a C-terminal HA tag and HIV P24 protein, respectively. Pseudovirus containing supernatants were 0.45 μm filtered and centrifuge at 128,000 g for 2 h, and then resuspend the pellet. We loaded 5 μL of purified pseudoviruses into the grids and incubated for 3 min, and 5 μL of 1% phosphotungstic acid was applied for negative staining. Pseudovirus particles were examined in a H-7000FA transmission electron microscope. The production process, measurements of infection and luciferase activity were conducted as described previously [7,37].

2.10. Nucleotide Sequence Accession Numbers

The sequences of viral genomes determined in this study have been submitted to GenBank sequence database under accession numbers MN477899 to MN477915.

3. Results

3.1. Prevalence of BtCoV HKU10

A total of 8004 fecal specimens were collected from 25 provinces in China and Louang Namtha province in Laos (Figure 1A). These bats belonged to 69 bat species of 6 bat families according to morphological or molecular identification. By RT-PCR and sequencing, we found 26 were positive for BtCoV HKU10 (Figure 1B and Supplementary Table S2). All positive samples were from three bat species: Aselliscus stoliczkanus (1/165), Hipposideros larvatus (8/196), and H. pomona (17/186) collected in southern and southwestern provinces of China (Yunnan, Guangxi, Guangdong and Hainan) and Louang Namtha province of Laos. We did not find any BtCoV-HKU10 or related viruses from 144 R. leschenaultia samples. Neither cytopathic effect nor viral replication was detected, indicating the failure of virus isolation.

3.2. Genomic Characterization of Different BtCoV HKU10 Strains

We sequenced the complete RdRp genes of all 26 strains. These sequences shared high identity to known BtCoV HKU10 (nt: 81.8–99.2%, deduced amino acid: 91.4–99.7%). Phylogenetic analysis based on the RdRp sequences indicated that the newly detected viruses, together with the known BtCoV HKU10 strains, formed eight lineages closely associated with the sampling locations (Figure 2). Lineages 1 and 2 were previously reported from Hong Kong and Guangdong, respectively, while the other six were newly discovered. Three of them were from Yunnan province and one each was from Hainan, Guangxi and Laos, respectively (Figure 2).
We then sequenced the complete genomes of 17 BtCoV HKU10 strains, considering the variations in RdRp sequences, bat species, and sampling location (Figure 3). These genomes ranged between 28,021 to 28,540 nt, with 38.19 to 39.81% G + C contents. All of them possessed the same putative transcription regulatory sequence (TRS) core motif as known BtCoV HKU10 except the TRS preceding ORF3 in YN2714/2727 (Table S3). The genome organization were similar to known BtCoV HKU10 (Figure 3). The major differences were located in the regions encoding accessory protein ORF7 downstream the N gene. In total, four types of ORF7: (1) GD141391, GD141402, GX160935, and HN140937 had similar accessory genes to known BtCoV-HKU10, which included ORF7a, ORF7b and ORF7c; (2) Strain YN3723 and YN3740 had a truncated ORF7c as a result of stop-codon substitution; (3) Strain LA7496 lost ORF7b and ORF7c; (4) the rest lost ORF7c.
The size of the 5′ untranslated regions of BtCoV HKU10 strains in H. pomona, H. larvatus, and A. stoliczkanus were 302, 299, and 302 nt, respectively. The replicase ORF1ab occupies 20.3 to 20.4 kb and encoded 16 putative proteins. The genomes of these viruses shared >73.3% nt sequence identities to known BtCoV HKU10. A separate deduced amino acid sequence comparison of seven conserved ORF1ab domains was performed as suggested by the International Committee on Taxonomy of Viruses (ICTV) for formal CoV species delineation. These domains possessed >90% aa identities to the prototype of BtCoV HKU10, suggesting that these viruses are BtCoV HKU10 variants.
We compared the genome nt identities, structural protein and non-structural protein aa identities among these viruses (Supplementary Table S4). Based on the full-length genomic sequences, viruses from A. stoliczkanus and H. pomona shared higher similarity (>82.6% and >85.0% nt identities, respectively) to known Hi-BatCoV HKU10 than those from H. lavartus (<72.0% nt identities). However, viruses from the same species were more closely related based on spike protein sequences. The spike proteins of BatCoV HKU10 from H. pomona and A. stoliczkanus were more similar (>79.8% aa identities) than the ones from H. larvatus and R. leschenaultia (<70.7% aa identities). For other structural and non-structural proteins, E, M, N and NS3 were conserved, while the NS7a, NS7b and NS7c possessed low aa identities among different strains (as low as 33.3%, 21.1% and 41.3%, respectively, for NS7a, b, and c).
We further analyzed the deduced amino acid sequence of polyprotein 1ab. The putative mature nonstructural proteins (NSPs) in the ORF1ab encoding the replicase complex was predicted as previously described. Five putative cleavage sites recognized by 3C-like proteinase (3CLpro) and papainlike (PLpro) in YN2714 strain showed difference to other BtCoV-HKU10 strains. Different putative cleavage motifs were also discovered between NSP1/2 and NSP5/6 in various strains of BtCoV HKU10s (Table 1).

3.3. Phylogenetic Analyses of Polyprotein 1ab and Structural Proteins

We constructed the phylogenetic trees using the deduced aa sequences of pp1ab or nucleotide sequences of S, E, M, N genes and genomes, respectively (Figure 4, Figure S3). All the trees except the S showed similar topology and they revealed that the BtCoV HKU10s were likely phylogeographically associated. The virus strains from south China (Guangdong, Guangxi, Hainan, and Hong Kong) were always clustered together despite the fact that they were from two different bat species. However, two strains detected in Yunnan province (YN2714 and YN2727) were distantly related to other strains in all five trees, suggesting a long independent evolutionary history in H. larvatus. The five HKU10 strains from H. pomona in Yunnan province were separated into two lineages with different topology, indicating virus genetic diversity in Yunnan. LA7496 from Laos formed a relatively distinct lineage from other HKU10s detected in Laos from its different hosts.
Remarkably, a different phylogenetic pattern was observed in the S protein tree, in which the HKU10s formed four distinct clusters which were correlated to their host species rather than sampling locations. The S proteins from H. Pomona were much closer to that in A. stoliczkanus than those in H. larvatus, though H. pomona and H. larvatus belonged to the same bat genus. Upon the phylogenetic tree, S protein of HKU10s in H. pomona could be further divided into two branches, and the virus strains in Yunnan were distributed in both branches. The viruses in each of the two branches shared high sequence similarities.

3.4. Recombination Analysis

The full-length genome sequences of BtCoV HKU10s were screened for potential recombination events as previously described. Briefly, the sequences were scanned sequentially by a bootscan algorithm and similarity plot analysis. Multiple potential recombination events were observed at ORF1a and S gene (Figure 5). YN3723 was likely to be a recombinant strain from three HKU10s (YN3740, YN4996, and YN7345) discovered in the same cave in Yunnan province, with strong p-value (<10−32). Breakpoints were identified in the genome of YN3723 at nt 20,275 and 23,265, the region between which encoded the C-terminal of pp1ab, S1 subunit and N-terminal of S2 subunit of S protein. In this region, the strain YN3723 was highly similar (95% nt identity) to the strain YN3740. Meanwhile, in the upstream ORF1ab sequence from nt 20,275, YN3723 displayed the highest genetic similarity (99% nt identity) to that of the strain YN4996.
Other BtCoV HKU10 were also screened for evidence of potential recombination events. However, no significant recombination breakpoint among these viruses and other HKU10 strains.

3.5. Estmation of Synonymous and Nonsynonymous Substitution Rates

High Ka/Ks ratios and substantial changes in the spike proteins of coronaviruses may reflect rapid viral evolution soon after introduction into a new animal host, as demonstrated in previous studies [38]. We examined the Ka/Ks ratios among BtCoV HKU10 in different hosts. A set of low ratios were observed in polymerase and structural genes of YN2714 and YN2727 (Table 2), providing another piece of evidence for the long-independent evolution history of BtCoV-HKU10 in H. lavatus bats. BtCoV HKU10 in H. pomona, R. leschenaultii and A. stoliczkanus possessed highly similar genome sequences (>80% sequence identity), indicating that interspecies transmission among these three bat species occurred recently, with subsequent viral adaption in the new host species [25]. To test the hypothesis, we then determined the Ka/Ks ratio for the various coding regions in HKU10s of different host species (Table 2). Higher Ka/Ks ratios were observed in Ro-BtCoV HKU10 S (0.382) than in other HKU10s (0.229–0.273). Higher Ka/Ks ratios were also observed in ORF7a, 7b, and 7c genes. The Ka/Ks ratio of 7b genes from H. pomona, H. larvatus and all HKU10s was 1.037, 1.161, and 0.951, respectively, which suggested that the 7b gene was likely under positive selection in these two clades.

3.6. Estimation of Divergence Date

We used relaxed clock model with UCed on RdRp sequences to calculate the divergence dates of BtCoV HKU10. The most recent common ancestor (MRCA) of BtCoV HKU10 was estimated at 1783 (HPDs, 1410.72-1981.43), approximately 238 years ago. For HKU10s in different host species, the MRCA of in H. pomona, R. leschenaultia, H. larvatus and A. stoliczkanus were estimated at 1851 (HPDs, 1681.82-1982.84), 1996 (HPDs, 1979.87-2004.09), 1995 (HPDs, 1965.80-2009.92), and 1995 (HPDs, 1958.63-2012.47), respectively (Figure 6). The estimated mean substitution rate of the RdRp data set was 3.705 × 10−4 substitution per site per year, which is comparable to previous estimations for other CoVs [25,39,40].

3.7. BtCoV HKU10 Spike-Mediated Pseudovirus Entry

To understand the infectivity of BtCoV HU10, we constructed pseudotyped lentiviruses carrying spike proteins from Ro-BtCoV HKU10 (NC018871) and LA7560, respectively. A total of 27 cell lines from human, bat, and other mammals were challenged. Huh-7 (human lung) and RD (human muscle) cells are susceptible to these viruses with and without treatment of trypsin. Pseudovirus could enter Tb1-Lu (bat lung) cell only when cells were treated with trypsin. HKU10 seemed unable to infect respiratory cells. As previously reported, the pseudovirus assay for MERS-CoV spike could enter A549, Calu-3, Caco-2, Huh-7, HEK-293T, LLC-MK2, and Vero cells [41,42] (Figure 7 and Table S5).

4. Discussion

In this study, we screened for BtCoV HKU10 from our archived 8004 bat fecal samples collected in 25 provinces of China. We found twenty-six (4.7%) BtCoV HKU10 strains in three bat species (547 bats) from more habitats including Guangdong, Guangxi, Hainan and Yunnan Province of China, and the Louang Namtha Province of Laos. Combining previous reports in Guangdong Province and Hong Kong in China, Thailand and Laos, we further demonstrate that BtCoV HKU10 are widely spread in south and southwest China and southeast Asia. Two bat species, A. stoliczkanus and H. larvatus, were firstly identified as the host of BtCoV HKU10. Our results expand the host ranges and geographic distributions of BtCoV HKU10.
The genomes of BtCoV HKU10, both previously and newly discovered are highly similar except in the S genes and accessory genes ORF7. Analyses based on the most conserved regions (ORF1ab and RdRp genes) indicate that the newly discovered viral strains together with previously described strains can be classified into eight lineages which are closely associated with their sampling locations. The previously reported BtCoV HKU10 strains from Guangdong and Hong Kong are closely related, although the two strains were found in two members of different bat suborders: H. pomona and R. leschenaultia. In our study, viruses from Yunnan belong to three lineages and the one found in H. larvatus is distantly related to BtCoV HKU10 found in H. pomona, suggesting an another interspecies transmission event of BtCoV HKU10. Previous studies suggested that the BtCoV HKU10 transmitted from R. leschenaulitia to H. pomona. However, in our extensive surveillances, we only found BtCoV HKU10 among Hipposideros bats, especially in H. pomona. We suggest that H. pomona, not Rousettus bats, are the primary reservoirs of BtCoV HKU10 based on the following reasons: first, the richest divergence, wider geographic distribution and longer independent evolution history of BatCoV HKU10 were observed in Hipposideros hosts; second, the Ka/Ks ratio of the previously reported S gene in R. leschenaultia (0.382) is higher than that from other hosts (0.229–0.273), which indicate a higher overall selection pressure and changing in Ro-BatCoV HKU10 spike; third, with previously described method, the molecular-clock analysis of all available RdRp genes dated the MRCA of all BatCoV HKU10 strains at 1783 (HPDs, 1410.72–1981.43). The MRCA of H. pomona and R. leschenaultia were estimated at around 1851 (HPDs, 1681.82–1982.84) and 1996 (HPDs, 1979.87–2004.09), respectively, indicating a recent interspecies transmission from H. Pomona to R. leschenaultia.
The major differences of BtCoV HKU10 from different hosts and locations were found in their S proteins. As we failed to cultivate BtCoV HKU10 in vitro, we used spike-pseudotyped lentivirus system for cell-entry and host adaptation studies. We screened 27 cells form human, bat and other mammals and found that the tested BtCoV HKU10 could infect multiples cell lines including human and bat species, providing evidence of potential interspecies transmission at cell-entry stage. Coronaviruses demonstrate a complex pattern for receptor recognition as both protein and sugar at cell surface could be used as cellular receptor mediating virus entry. The sugar-receptor-binding region is located at the N-terminal domain (NTD) of S1 region [43,44,45], while the protein-receptor-binding domain is mainly located at the C-terminal domain (CTD) of S1 domain [43,46,47,48]. Protease cleavage triggers for membrane fusion by coronavirus spike protein, which is another determining factors for coronavirus entry. The host proteases, including proprotein convertases, extracellular proteases, cell surface proteases, and lysosomal proteases cleave coronaviruses at different stages during virus maturation and entry [49,50]. In this study, we detected the spike protein of the pseudoviruses produced by 293T cells with Western blotting and found no evidence for human proprotein convertases cleavage (Figure S2). The mature pseudo-typed virions with intact spike provide an evidence that human proprotein convertases may not participate in the BtCoV HKU10 spike maturation just as SARS-CoV [51].
Bats harbor a vast number of coronaviruses. The ability to fly, their worldwide distribution and intimate contact with human and other domestic animals provide bats with opportunities for virus spillover from their natural reservoirs to human society, as demonstrated by the SARS and MERS outbreaks. With increasing genetic information about bat coronaviruses, bats are considered as the natural reservoirs of most alpha- and beta-coronaviruses. However, little information was known regarding the risk of their interspecies transmission. Further investigations will be focused on virus isolation, receptor usage and the pathogenesis of these CoVs, which will provide further information for predicting the potential spillover of these viruses.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v13101962/s1, Figure S1: Strategy of BatCoV HKU10 genomes sequencing. Figure S2: Characterization of BtCoV HKU10 spike mediated pseudovirus. Figure S3: Phylogenetic tree of BatCoV HKU10. Table S1: Primer pairs for HKU10 detection and genome sequencing. Table S2: Detection of HKU10 in bats by RT-PCR. Table S3: Genome size, coding potential, and putative transcription regulatory sequences of 11 BtCoV HKU10 strains. Table S4: Comparison of whole genomes (nt) and coding regions (aa) of representation strains of BtCoV HKU10. Table S5: Susceptibility test of Ro-BtCoV HKU10 and BtCoV HKU10 LA7560.

Author Contributions

Conceptualization, C.-M.L., N.W. and B.H.; methodology, N.W. and C.-M.L.; formal analysis, N.W., C.-M.L., B.H., X.-L.Y. and H.-Z.L.; investigation, N.W. and C.-M.L.; resources, W.Z., B.L., L.-B.Z. and C.P.; writing—original draft preparation, N.W. and C.-M.L.; writing—review and editing, Z.-L.S. and B.H.; visualization, N.W.; supervision, Z.-L.S. and B.H.; project administration, Y.Z.; funding acquisition, Z.-L.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was jointly funded by the strategic priority research program of the Chinese academy of sciences (XDB29010101), the National Natural Science Foundation of China (U2002218, 81290341), Scientific and technological basis special project (2013FY113500) from the Ministry of Science and Technology of the People’s Republic of China; United States Agency for International Development (USAID) Emerging Pandemic Threats PREDICT project grant (Cooperative Agreement no. AID-OAA-A-14-00102).

Institutional Review Board Statement

The study was conducted according with the Guide for the Care and Use of Wild Mammals in Re-search of the People’s Republic of China, and approved by Animal Ethics Committee of the Wuhan Institute of Virology (WIVH05210201).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available under accession numbers MN477899 to MN477915.

Acknowledgments

We thank Yunzhi Zhang from Dali University and Yunnan Institute of Endemic Diseases Control and Prevention for his contribution to part of the fieldwork in this study.

Conflicts of Interest

All the authors declare no conflict of interest.

References

  1. Knipe, D.M.; Howley, P.M.; Cohen, J.I.; Griffin, D.E.; Lamb, R.A.; Martin, M.A.; Racaniello, V.R.; Roizman, B. Fields Virology, 6th ed.; Wolters Kluwer/Lippincott Williams and Wilkins: Philadelphia, PA, USA, 2013; Volume 1. [Google Scholar]
  2. International Committee on Taxonomy of Viruses. ICTV Master Species List 2018a v1. 2018. Available online: http://www.ictvonline.org/ (accessed on 7 April 2019).
  3. Miller, W.A.; Koev, G. Synthesis of subgenomic RNAs by positive-strand RNA viruses. Virology 2000, 273, 1–8. [Google Scholar] [CrossRef] [Green Version]
  4. Woo, P.C.Y.; Lau, S.K.P.; Lam, C.S.F.; Lau, C.C.Y.; Tsang, A.K.L.; Lau, J.H.N.; Bai, R.; Teng, J.L.L.; Tsang, C.C.C.; Wang, M.; et al. Discovery of Seven Novel Mammalian and Avian Coronaviruses in the Genus Deltacoronavirus Supports Bat Coronaviruses as the Gene Source of Alphacoronavirus and Betacoronavirus and Avian Coronaviruses as the Gene Source of Gammacoronavirus and Deltacoronavirus. J. Virol. 2012, 86, 3995–4008. [Google Scholar] [PubMed] [Green Version]
  5. Ge, X.Y.; Li, J.L.; Yang, X.L.; Chmura, A.A.; Zhu, G.J.; Epstein, J.H.; Mazet, J.K.; Hu, B.; Zhang, W.; Peng, C.; et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature 2013, 503, 535. [Google Scholar] [CrossRef] [PubMed]
  6. Hu, B.; Zeng, L.P.; Yang, X.L.; Ge, X.Y.; Zhang, W.; Li, B.; Xie, J.Z.; Shen, X.R.; Zhang, Y.Z.; Wang, N.; et al. Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into the origin of SARS coronavirus. PLoS Pathog. 2017, 13, e1006698. [Google Scholar] [CrossRef]
  7. Luo, C.-M.; Wang, N.; Yang, X.-L.; Liu, H.-Z.; Zhang, W.; Li, B.; Hu, B.; Peng, C.; Geng, Q.-B.; Zhu, G.-J.; et al. Discovery of novel bat coronaviruses in south China that use the same receptor as MERS coronavirus. J. Virol. 2018, 92, e00116-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Tao, Y.; Shi, M.; Chommanard, C.; Queen, K.; Zhang, J.; Markotter, W.; Kuzmin, I.V.; Holmes, E.C.; Tong, S. Surveillance of Bat Coronaviruses in Kenya Identifies Relatives of Human Coronaviruses NL63 and 229E and Their Recombination History. J. Virol. 2017, 91, e01953-16. [Google Scholar] [CrossRef] [Green Version]
  9. Corman, V.M.; Baldwin, H.J.; Tateno, A.F.; Zerbinati, R.M.; Annan, A.; Owusu, M.; Nkrumah, E.E.; Maganga, G.D.; Oppong, S.; Adu-Sarkodie, Y.; et al. Evidence for an Ancestral Association of Human Coronavirus 229E with Bats. J. Virol. 2015, 89, 11858–11870. [Google Scholar] [CrossRef] [Green Version]
  10. Lau, S.K.; Li, K.S.; Tsang, A.K.; Lam, C.S.; Ahmed, S.; Chen, H.; Chan, K.H.; Woo, P.C.; Yuen, K.Y. Genetic characterization of Betacoronavirus lineage C viruses in bats reveals marked sequence divergence in the spike protein of pipistrellus bat coronavirus HKU5 in Japanese pipistrelle: Implications for the origin of the novel Middle East respiratory syndrome coronavirus. J. Virol. 2013, 87, 8638–8650. [Google Scholar]
  11. Zhou, P.; Yang, X.L.; Wang, X.G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.R.; Zhu, Y.; Li, B.; Huang, C.L.; et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Lai, M.M.C. Recombination in large RNA viruses: Coronaviruses. Semin. Virol. 1996, 7, 381–388. [Google Scholar] [CrossRef]
  13. Guan, Y.; Zheng, B.J.; He, Y.Q.; Liu, X.L.; Zhuang, Z.X.; Cheung, C.L.; Luo, S.W.; Li, P.H.; Zhang, L.J.; Guan, Y.J.; et al. Isolation and characterization of viruses related to the SARS coronavirus from animals in Southern China. Science 2003, 302, 276–278. [Google Scholar] [CrossRef] [Green Version]
  14. Ksiazek, T.G.; Erdman, D.; Goldsmith, C.S.; Zaki, S.R.; Peret, T.; Emery, S.; Tong, S.X.; Urbani, C.; Comer, J.A.; Lim, W.; et al. A novel coronavirus associated with severe acute respiratory syndrome. N. Engl. J. Med. 2003, 348, 1953–1966. [Google Scholar] [CrossRef] [PubMed]
  15. Drosten, C.; Gunther, S.; Preiser, W.; van der Werf, S.; Brodt, H.R.; Becker, S.; Rabenau, H.; Panning, M.; Kolesnikova, L.; Fouchier, R.A.M.; et al. Identification of a novel coronavirus in patients with severe acute respiratory syndrome. N. Engl. J. Med. 2003, 348, 1967–1976. [Google Scholar] [CrossRef] [PubMed]
  16. Zaki, A.M.; van Boheemen, S.; Bestebroer, T.M.; Osterhaus, A.; Fouchier, R.A.M. Isolation of a Novel Coronavirus from a Man with Pneumonia in Saudi Arabia. N. Engl. J. Med. 2012, 367, 1814–1820. [Google Scholar] [CrossRef] [PubMed]
  17. Li, W.D.; Shi, Z.L.; Yu, M.; Ren, W.Z.; Smith, C.; Epstein, J.H.; Wang, H.Z.; Crameri, G.; Hu, Z.H.; Zhang, H.J.; et al. Bats are natural reservoirs of SARS-like coronaviruses. Science 2005, 310, 676–679. [Google Scholar] [CrossRef] [PubMed]
  18. Anthony, S.J.; Gilardi, K.; Menachery, V.D.; Goldstein, T.; Ssebide, B.; Mbabazi, R.; Navarrete-Macias, I.; Liang, E.; Wells, H.; Hicks, A.; et al. Further Evidence for Bats as the Evolutionary Source of Middle East Respiratory Syndrome Coronavirus. mBio 2017, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Zhou, P.; Fan, H.; Lan, T.; Yang, X.-L.; Shi, W.; Zhang, W.; Zhu, Y.; Zhang, Y.; Xie, Q.-M.; Mani, S.; et al. Fatal swine acute diarrhea syndrome caused by an HKU2-related coronavirus of bat origin. Nature 2018, 556, 255–258. [Google Scholar] [CrossRef] [PubMed]
  20. Pan, Y.F.; Tian, X.Y.; Qin, P.; Wang, B.; Zhao, P.W.; Yang, Y.L.; Wang, L.X.; Wang, D.D.; Song, Y.H.; Zhang, X.B.; et al. Discovery of a novel swine enteric alphacoronavirus (SeACoV) in southern China. Vet. Microbiol. 2017, 211, 15–21. [Google Scholar] [CrossRef]
  21. Gong, L.; Li, J.; Zhou, Q.F.; Xu, Z.C.; Chen, L.; Zhang, Y.; Xue, C.Y.; Wen, Z.F.; Cao, Y.C. A New Bat-HKU2-like Coronavirus in Swine, China, 2017. Emerg. Infect. Dis 2017, 23, 1607–1609. [Google Scholar] [CrossRef] [Green Version]
  22. Ge, X.Y.; Wang, N.; Zhang, W.; Hu, B.; Li, B.; Zhang, Y.Z.; Zhou, J.H.; Luo, C.M.; Yang, X.L.; Wu, L.J.; et al. Coexistence of multiple coronaviruses in several bat colonies in an abandoned mineshaft. Virol. Sin. 2016, 31, 31–40. [Google Scholar] [CrossRef] [Green Version]
  23. Wacharapluesadee, S.; Duengkae, P.; Rodpan, A.; Kaewpom, T.; Maneeorn, P.; Kanchanasaka, B.; Yingsakmongkon, S.; Sittidetboripat, N.; Chareesaen, C.; Khlangsap, N.; et al. Diversity of coronavirus in bats from Eastern Thailand. Virol. J. 2015, 12, 57. [Google Scholar] [CrossRef] [Green Version]
  24. Lacroix, A.; Duong, V.; Hul, V.; San, S.; Davun, H.; Omaliss, K.; Chea, S.; Hassanin, A.; Theppangna, W.; Silithammavong, S.; et al. Genetic diversity of coronaviruses in bats in Lao PDR and Cambodia. Infect. Genet. Evol. 2017, 48, 10–18. [Google Scholar] [CrossRef]
  25. Lau, S.K.P.; Li, K.S.M.; Tsang, A.K.L.; Shek, C.T.; Wang, M.; Choi, G.K.Y.; Guo, R.T.; Wong, B.H.L.; Poon, R.W.S.; Lam, C.S.F.; et al. Recent Transmission of a Novel Alphacoronavirus, Bat Coronavirus HKU10, from Leschenault’s Rousettes to Pomona Leaf-Nosed Bats: First Evidence of Interspecies Transmission of Coronavirus between Bats of Different Suborders. J. Virol. 2012, 86, 11906–11918. [Google Scholar] [CrossRef] [Green Version]
  26. Wang, N.; Luo, C.M.; Liu, H.Z.; Yang, X.L.; Hu, B.; Zhang, W.; Li, B.; Zhu, Y.; Zhu, G.J.; Shen, X.R.; et al. Characterization of a New Member of Alphacoronavirus with Unique Genomic Features in Rhinolophus Bats. Viruses 2019, 11, 379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Watanabe, S.; Masangkay, J.S.; Nagata, N.; Morikawa, S.; Mizutani, T.; Fukushi, S.; Alviola, P.; Omatsu, T.; Ueda, N.; Iha, K.; et al. Bat coronaviruses and experimental infection of bats, the Philippines. Emerg Infect Dis 2010, 16, 1217–1223. [Google Scholar] [CrossRef] [PubMed]
  28. Irwin, D.M.; Kocher, T.D.; Wilson, A.C. Evolution of the Cytochrome-B Gene of Mammals. J. Mol. Evol. 1991, 32, 128–144. [Google Scholar] [CrossRef] [PubMed]
  29. Mayer, F.; von Helversen, O. Cryptic diversity in European bats. Proc. R. Soc. B Biol. Sci. 2001, 268, 1825–1832. [Google Scholar] [CrossRef] [Green Version]
  30. Yang, D.; Leibowitz, J.L. The structure and functions of coronavirus genomic 3′ and 5′ ends. Virus Res. 2015, 206, 120–133. [Google Scholar] [CrossRef]
  31. Hall, B.G. Building Phylogenetic Trees from Molecular Data with MEGA. Mol. Biol. Evol. 2013, 30, 1229–1235. [Google Scholar] [CrossRef] [Green Version]
  32. Martin, D.P.; Lemey, P.; Lott, M.; Moulton, V.; Posada, D.; Lefeuvre, P. RDP3: A flexible and fast computer program for analyzing recombination. Bioinformatics 2010, 26, 2462–2463. [Google Scholar] [CrossRef]
  33. Lole, K.S.; Bollinger, R.C.; Paranjape, R.S.; Gadkari, D.; Kulkarni, S.S.; Novak, N.G.; Ingersoll, R.; Sheppard, H.W.; Ray, S.C. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J. Virol. 1999, 73, 152–160. [Google Scholar] [CrossRef] [Green Version]
  34. Siltberg, J.; Liberles, D.A. A simple covarion-based approach to analyse nucleotide substitution rates. J. Evol. Biol. 2002, 15, 588–594. [Google Scholar] [CrossRef] [Green Version]
  35. Luna, L.K.D.; Heiser, V.; Regamey, N.; Panning, M.; Drexler, J.F.; Mulangu, S.; Poon, L.; Baumgarte, S.; Haijema, B.J.; Kaiser, L.; et al. Generic detection of coronaviruses and differentiation at the prototype strain level by reverse transcription-PCR and nonfluorescent low-density microarray. J. Clin. Microbiol. 2007, 45, 1049–1052. [Google Scholar] [CrossRef] [Green Version]
  36. Crameri, G.; Todd, S.; Grimley, S.; McEachern, J.A.; Marsh, G.A.; Smith, C.; Tachedjian, M.; De Jong, C.; Virtue, E.R.; Yu, M.; et al. Establishment, immortalisation and characterisation of pteropid bat cell lines. PLoS ONE 2009, 4, e8266. [Google Scholar] [CrossRef]
  37. Yang, Y.; Du, L.Y.; Liu, C.; Wang, L.L.; Ma, C.Q.; Tang, J.; Baric, R.S.; Jiang, S.; Li, F. Receptor usage and cell entry of bat coronavirus HKU4 provide insight into bat-to-human transmission of MERS coronavirus. Proc. Natl. Acad. Sci. USA 2014, 111, 12516–12521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Song, H.D.; Tu, C.C.; Zhang, G.W.; Wang, S.Y.; Zheng, K.; Lei, L.C.; Chen, Q.X.; Gao, Y.W.; Zhou, H.Q.; Xiang, H.; et al. Cross-host evolution of severe acute respiratory syndrome coronavirus in palm civet and human. Proc. Natl. Acad. Sci. USA 2005, 102, 2430–2435. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Pyrc, K.; Dijkman, R.; Deng, L.; Jebbink, M.F.; Ross, H.A.; Berkhout, B.; Van der Hoek, L. Mosaic structure of human coronavirus NL63, one thousand years of evolution. J. Mol. Biol. 2006, 364, 964–973. [Google Scholar] [CrossRef] [PubMed]
  40. Vijgen, L.; Keyaerts, E.; Moes, E.; Thoelen, I.; Wollants, E.; Lemey, P.; Vandamme, A.M.; Van Ranst, M. Complete genomic sequence of human coronavirus OC43: Molecular clock analysis suggests a relatively recent zoonotic coronavirus transmission event. J. Virol. 2005, 79, 1595–1604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Raj, V.S.; Mou, H.H.; Smits, S.L.; Dekkers, D.H.W.; Muller, M.A.; Dijkman, R.; Muth, D.; Demmers, J.A.A.; Zaki, A.; Fouchier, R.A.M.; et al. Dipeptidyl peptidase 4 is a functional receptor for the emerging human coronavirus-EMC. Nature 2013, 495, 251–254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Chan, J.F.; Chan, K.H.; Choi, G.K.; To, K.K.; Tse, H.; Cai, J.P.; Yeung, M.L.; Cheng, V.C.; Chen, H.; Che, X.Y.; et al. Differential cell line susceptibility to the emerging novel human betacoronavirus 2c EMC/2012: Implications for disease pathogenesis and clinical manifestation. J. Infect. Dis. 2013, 207, 1743–1752. [Google Scholar] [CrossRef] [Green Version]
  43. Liu, C.; Tang, J.; Ma, Y.M.; Liang, X.Y.; Yang, Y.; Peng, G.Q.; Qi, Q.Q.; Jiang, S.B.; Li, J.R.; Du, L.Y.; et al. Receptor Usage and Cell Entry of Porcine Epidemic Diarrhea Coronavirus. J. Virol. 2015, 89, 6121–6125. [Google Scholar] [CrossRef] [Green Version]
  44. Li, W.T.; Hulswit, R.J.G.; Widjaja, I.; Raj, V.S.; McBride, R.; Peng, W.J.; Widagdo, W.; Tortorici, M.A.; van Dieren, B.; Lang, Y.; et al. Identification of sialic acid-binding function for the Middle East respiratory syndrome coronavirus spike glycoprotein. Proc. Natl. Acad. Sci. USA 2017, 114, E8508–E8517. [Google Scholar] [CrossRef] [Green Version]
  45. Hulswit, R.J.G.; Lang, Y.F.; Bakkers, M.J.G.; Li, W.T.; Li, Z.S.; Schouten, A.; Ophorst, B.; van Kuppeveld, F.J.M.; Boons, G.J.; Bosch, B.J.; et al. Human coronaviruses OC43 and HKU1 bind to 9-O-acetylated sialic acids via a conserved receptor-binding site in spike protein domain A. Proc. Natl. Acad. Sci. USA 2019, 116, 2681–2690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Lu, G.W.; Hu, Y.W.; Wang, Q.H.; Qi, J.X.; Gao, F.; Li, Y.; Zhang, Y.F.; Zhang, W.; Yuan, Y.; Bao, J.K.; et al. Molecular basis of binding between novel human coronavirus MERS-CoV and its receptor CD26. Nature 2013, 500, 227. [Google Scholar] [CrossRef] [Green Version]
  47. Wong, S.K.; Li, W.H.; Moore, M.J.; Choe, H.; Farzan, M. A 193-amino acid fragment of the SARS coronavirus S protein efficiently binds angiotensin-converting enzyme 2. J. Biol. Chem. 2004, 279, 3197–3201. [Google Scholar] [CrossRef] [Green Version]
  48. Hofmann, H.; Simmons, G.; Rennekamp, A.J.; Chaipan, C.; Gramberg, T.; Heck, E.; Geier, M.; Wegele, A.; Marzi, A.; Bates, P.; et al. Highly conserved regions within the spike proteins of human coronaviruses 229E and NL63 determine recognition of their respective cellular receptors. J. Virol. 2006, 80, 8639–8652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Zheng, Y.; Shang, J.; Yang, Y.; Liu, C.; Wan, Y.S.; Geng, Q.B.; Wang, M.; Baric, R.; Li, F. Lysosomal Proteases Are a Determinant of Coronavirus Tropism. J. Virol. 2018, 92, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Millet, J.K.; Whittaker, G.R. Host cell proteases: Critical determinants of coronavirus tropism and pathogenesis. Virus Res. 2015, 202, 120–134. [Google Scholar] [CrossRef]
  51. Li, F. Structure, Function, and Evolution of Coronavirus Spike Proteins. Annu. Rev. Virol. 2016, 3, 237–261. [Google Scholar] [CrossRef] [Green Version]
Viruses 13 01962 g001 550
Figure 1. Locations of sampling (A) and BtCoV HKU10 positives (B). Sampling locations are in gray and bat species are listed in color. BtCoV HKU10 positives are marked in square (A) and in red dot (B).
Figure 1. Locations of sampling (A) and BtCoV HKU10 positives (B). Sampling locations are in gray and bat species are listed in color. BtCoV HKU10 positives are marked in square (A) and in red dot (B).
Viruses 13 01962 g001
Viruses 13 01962 g002 550
Figure 2. Phylogenetic analysis of BtCoV HKU10 detected in bats. Full-length RdRp genes of viruses detected in this study were aligned with those of published BtCoV HKU10 and representative CoV strains. The tree was constructed by the maximum-likelihood method with bootstrap values determined with 1000 replicates. The scale bar indicates the estimated number of substitutions per 20 nucleotides. Filled diamonds indicate the CoVs detected in this study. FIPV, feline infectious peritonitis virus; PEDV, porcine epidemic diarrhea virus.
Figure 2. Phylogenetic analysis of BtCoV HKU10 detected in bats. Full-length RdRp genes of viruses detected in this study were aligned with those of published BtCoV HKU10 and representative CoV strains. The tree was constructed by the maximum-likelihood method with bootstrap values determined with 1000 replicates. The scale bar indicates the estimated number of substitutions per 20 nucleotides. Filled diamonds indicate the CoVs detected in this study. FIPV, feline infectious peritonitis virus; PEDV, porcine epidemic diarrhea virus.
Viruses 13 01962 g002
Viruses 13 01962 g003 550
Figure 3. Schematic diagram of genomic organization of BtCoV HKU10 strains. The genomic regions or ORFs are described and compared in Supplementary Tables S3 and S4. Solid bars indicate conserved genes and grey letters indicate group-specific genes. Upper letters indicate structural proteins and lower letters indicate nonstructural proteins (p1a and p1b) and accessory proteins.
Figure 3. Schematic diagram of genomic organization of BtCoV HKU10 strains. The genomic regions or ORFs are described and compared in Supplementary Tables S3 and S4. Solid bars indicate conserved genes and grey letters indicate group-specific genes. Upper letters indicate structural proteins and lower letters indicate nonstructural proteins (p1a and p1b) and accessory proteins.
Viruses 13 01962 g003
Viruses 13 01962 g004a 550Viruses 13 01962 g004b 550Viruses 13 01962 g004c 550
Figure 4. Phylogenetic trees of pp1ab (aa sequences), (S, N, M, E) (nt sequences). Trees were constructed by the maximum likelihood method using the Poisson model with bootstrap value determined with 1000 replicates. The scale bar indicates the estimated number of substitutions per 10 amino acids or nucleotides. Filled diamonds indicate viruses detected in this study.
Figure 4. Phylogenetic trees of pp1ab (aa sequences), (S, N, M, E) (nt sequences). Trees were constructed by the maximum likelihood method using the Poisson model with bootstrap value determined with 1000 replicates. The scale bar indicates the estimated number of substitutions per 10 amino acids or nucleotides. Filled diamonds indicate viruses detected in this study.
Viruses 13 01962 g004aViruses 13 01962 g004bViruses 13 01962 g004c
Viruses 13 01962 g005 550
Figure 5. Evidence of recombination in BtCoV HKU10s. Similarity (A) (window of 40 nt, step size of 40 nt) and recombination (B) (window of 1500 nt, step size of 150 nt) plots were generated using Simplot (V3.5.1) with default settings. Full-length genome sequence of YN3723 was used as query sequence and YN3740, YN4996 and YN7345 as reference sequences. All analyses were performed with Kimura model, a window size of 1500 base pairs, and a step size of 150 base pairs. The map of query genome sequences (C) are used to position breakpoints.
Figure 5. Evidence of recombination in BtCoV HKU10s. Similarity (A) (window of 40 nt, step size of 40 nt) and recombination (B) (window of 1500 nt, step size of 150 nt) plots were generated using Simplot (V3.5.1) with default settings. Full-length genome sequence of YN3723 was used as query sequence and YN3740, YN4996 and YN7345 as reference sequences. All analyses were performed with Kimura model, a window size of 1500 base pairs, and a step size of 150 base pairs. The map of query genome sequences (C) are used to position breakpoints.
Viruses 13 01962 g005
Viruses 13 01962 g006 550
Figure 6. Estimation of the tMRCA of BtCoV HKU10. The time-scaled phylogeny was summarized from all MCMC phylogenies of the RdRp gene data set analyzed under the relaxed-clock model with an exponential distribution (Uced) in BEAST version 2.4.4. Filled diamonds indicate viruses detected in this study.
Figure 6. Estimation of the tMRCA of BtCoV HKU10. The time-scaled phylogeny was summarized from all MCMC phylogenies of the RdRp gene data set analyzed under the relaxed-clock model with an exponential distribution (Uced) in BEAST version 2.4.4. Filled diamonds indicate viruses detected in this study.
Viruses 13 01962 g006
Viruses 13 01962 g007 550
Figure 7. HKU10 spike-mediated pseudovirus entry. Ro-BtCoV HKU10, LA7560, and MERS-CoV spike-mediated pseudovirus were produced to screen their cell tropism. For Trypsin treatment, pseudoviruses were incubated with 50 μg/mL TPCK-treated trypsin at room temperature for 10 min, and then add 50 μg/mL soybean trypsin inhibitor to stop the reactions. Trypsin treated or untreated pseudoviruses were added to cells and incubated. Then, 48 h post infection, cells were lysed and luciferase were detected.
Figure 7. HKU10 spike-mediated pseudovirus entry. Ro-BtCoV HKU10, LA7560, and MERS-CoV spike-mediated pseudovirus were produced to screen their cell tropism. For Trypsin treatment, pseudoviruses were incubated with 50 μg/mL TPCK-treated trypsin at room temperature for 10 min, and then add 50 μg/mL soybean trypsin inhibitor to stop the reactions. Trypsin treated or untreated pseudoviruses were added to cells and incubated. Then, 48 h post infection, cells were lysed and luciferase were detected.
Viruses 13 01962 g007
Table
Table 1. Characteristics of predicted nonstructural proteins of ORF1ab in different strains of BtCoV HKU10 in this study.
Table 1. Characteristics of predicted nonstructural proteins of ORF1ab in different strains of BtCoV HKU10 in this study.
NSPPutative Functional Domain(s)YN2714LA7496YN887GD141391
Amino Acids Position in ORF1abPredicted Size (aa of Protein)C-End Predicted Cleavage SiteAmino Acids Position in ORF1abPredicted Size (aa of Protein)C-End Predicted Cleavage SiteAmino Acids Position in ORF1abPredicted Size (aa of Protein)C-End Predicted Cleavage SiteAmino Acids Position in ORF1abPredicted Size (aa of Protein)C-End Predicted Cleavage Site
NSP1UnknownM1-A195195NA|GPM1-A195195VA|KPM1-A195195VA|KPM1-A195195VA|KV
NSP2UnknownG196-G888693TG|GGK196-G888693RG|SGK196-G888693RG|SGK196-G888693RG|SG
NSP3ADRP, PL2 proG889-G24631575CG|SGS889-G25291641CG|SGS889-G25231635CG|SGS889-G25181630CG|SG
NSP4Hydrophobid domainS2464-Q2941478LQ|AGS2530-Q3007478LQ|SGS2524-Q3001478LQ|SGS2519-Q2996478LQ|SG
NSP53CL proA2942-Q3243302LQ|STS3008-Q3309302LQ|STS3002-Q3303302LQ|SNS2997-Q3298302LQ|SN
NSP6Hydrophobid domainS3244-Q3519276VQ|SKS3310-Q3585276VQ|SKS3304-Q3579276VQ|SKS3299-Q3574276VQ|SK
NSP7ReplicaseS3520-Q360283LQ|SVS3586-Q366883LQ|SVS3580-Q366283LQ|SVS3575-Q365783LQ|SV
NSP8ReplicaseS3603-Q3797195LQ|NNS3669-Q3863195LQ|NNS3663-Q3857195LQ|NNS3658-Q3852195LQ|NN
NSP9ReplicaseN3798-Q3905108LQ|AGN3864-Q3971108LQ|AGN3858-Q3965108LQ|AGN3853-Q3960108LQ|AG
NSP10RNA synthesis proteinA3906-Q4041136VQ|AFA3972-Q4108137MQ|AFA3966-Q4102137MQ|AFA3961-Q4097137MQ|AF
NSP11Unknown (short peptide at the end of ORF1a)A4042-N405817 A4109-N412517 A4103-N411917 A4098-N411417
NSP12RdRpA4042-Q4968927LQ|AAA4109-Q5035927LQ|SAA4103-Q5029927LQ|SAA4098-Q5024927LQ|SA
NSP13Hel, NTPaseA4969-Q5565597LQ|AGS5036-Q5632597LQ|AGS5030-Q5626597LQ|AGS5025-Q5621597LQ|AG
NSP14ExoN, NMTA5566-Q6083518LQ|SLA5633-Q6150518LQ|SLA5627-Q6144518LQ|SLA5622-Q6139518LQ|SL
NSP15NeudoUS6084-Q6422339LQ|SAS6151-Q6487337LQ|SAS6145-Q6483339LQ|SAS6140-Q6476337LQ|SA
NSP162′-O-MTS6423-K6724302 S6488-C6799312 S6484-K6785302 S6477-H6783307
Table
Table 2. Estimation of nonsynonymous and synonymous substitution rates in the genomes of BtCoV HKU10 in different host species.
Table 2. Estimation of nonsynonymous and synonymous substitution rates in the genomes of BtCoV HKU10 in different host species.
GenesKa/Ks Ratio
BtCoV HKU10 BtCoV HKU10BtCoV HKU10 BtCoV HKU10 BtCoV HKU10
(H. pomona, 20 Strains)(A. stoliczkanus, 1 Strain)(H. lavatus, 2 Strains)(R. leschenaultia, 2 Strains)(All, 25 Strains)
NSP10.3540.2140.2370.1500.297
NSP20.3090.1310.3910.1710.403
NSP30.3590.1900.3030.1390.370
NSP40.1890.1040.2250.0590.274
NSP50.1760.0330.2050.0030.221
NSP60.3170.1670.0970.0530.477
NSP70.265Ka = 0, Ks = 0.035Ka = 0, Ks = 0Ka = 0, Ks = 00.166
NSP80.1740.5000.0100.0120.287
NSP90.198Ka = 0, Ks = 00.0730.0310.598
NSP100.200Ka = 0, Ks = 0.0320.1200.0560.138
NSP11Ka = 0, Ks = 0.066Ka = 0, Ks = 0.131Ka = 0, Ks = 0.066Ka = 0, Ks = 01.176
NSP120.1230.0530.0920.0120.067
NSP130.1120.1610.0870.0080.0738
NSP140.0810.0520.1420.0340.109
NSP150.1270.1140.0860.0260.167
NSP160.1030.1170.1260.1140.105
S0.2400.2430.2720.3820.259
ORF3a0.0690.169Ka = 0, Ks = 00.0510.274
E00.2040.227Ka = 0, Ks = 00.218
M0.1920.0240.1470.2750.066
N0.1510.1880.2180.1460.207
ORF7a0.3390.4030.4420.3360.662
ORF7b1.0368NA1.1610.1360.951
ORF7c0.630NANA0.6090.4433
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, N.; Luo, C.-M.; Yang, X.-L.; Liu, H.-Z.; Zhang, L.-B.; Zhang, W.; Li, B.; Zhu, Y.; Peng, C.; Shi, Z.-L.; et al. Genomic Characterization of Diverse Bat Coronavirus HKU10 in Hipposideros Bats. Viruses 2021, 13, 1962. https://doi.org/10.3390/v13101962

AMA Style

Wang N, Luo C-M, Yang X-L, Liu H-Z, Zhang L-B, Zhang W, Li B, Zhu Y, Peng C, Shi Z-L, et al. Genomic Characterization of Diverse Bat Coronavirus HKU10 in Hipposideros Bats. Viruses. 2021; 13(10):1962. https://doi.org/10.3390/v13101962

Chicago/Turabian Style

Wang, Ning, Chu-Ming Luo, Xing-Lou Yang, Hai-Zhou Liu, Li-Biao Zhang, Wei Zhang, Bei Li, Yan Zhu, Cheng Peng, Zheng-Li Shi, and et al. 2021. "Genomic Characterization of Diverse Bat Coronavirus HKU10 in Hipposideros Bats" Viruses 13, no. 10: 1962. https://doi.org/10.3390/v13101962

APA Style

Wang, N., Luo, C. -M., Yang, X. -L., Liu, H. -Z., Zhang, L. -B., Zhang, W., Li, B., Zhu, Y., Peng, C., Shi, Z. -L., & Hu, B. (2021). Genomic Characterization of Diverse Bat Coronavirus HKU10 in Hipposideros Bats. Viruses, 13(10), 1962. https://doi.org/10.3390/v13101962

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop