PMC:3713832 / 1489-8665 JSONTXT

{"target":"https://pubannotation.org/docs/sourcedb/PMC/sourceid/3713832","sourcedb":"PMC","sourceid":"3713832","source_url":"https://www.ncbi.nlm.nih.gov/pmc/3713832","text":"The Study\nWe obtained pharyngeal and anal swab specimens of 414 insectivorous bats. Samples of each species were pooled and then processed with a viral particle–protected nucleic acid purification method (6). The extracted RNA and DNA were amplified by sequence-independent PCR. The amplified viral nucleic acid libraries of the bat species were then sequenced with the Illumina/Solexa GAII sequencer (Illumina, San Diego, CA, USA). Those reads generated by the Illumina/Solexa GAII with length of 80 bases were directly aligned to the protein sequences in the National Center for Biotechnology Information nonredundant protein database by the blastx program in the BLAST software package, version 2.2.22 (www.ncbi.nlm.nih.gov/blast) with parameters “-e 1e-5 -F T -b 10 -v 10.” No assembly was performed before alignment. Sequence similarity–based taxonomic assignments were conducted as described (7). We found 1,075 reads of betacoronavirus in Rhinolophus pusillus bats in Shaanxi and 92 reads of betacoronavirus in Chaerephon plicata bats in Yunnan.\nWe estimated the approximate locations of those reads on the CoV genome and their relative distances on the basis of alignment results exported with MEGAN 4–MetaGenome Analyzer (http://ab.inf.uni-tuebingen.de/software/megan/). The located reads were then used for reads-based nested PCR to identify genomic sequences. We established the complete genome sequences of 2 betacoronaviruses (Bat Rp-coronavirus/Shaanxi2011 and Bat Cp-coronavirus/Yunnan2011), which are 29,484 nt and 29,452 nt, respectively. The G+C content of Bat Rp-coronavirus/Shaanxi2011 and Bat Cp-coronavirus/Yunnan2011 is 41.6% and 40.9%, respectively.\nWe conducted complete genome comparison and phylogenetic analysis on the basis of polymerase and spike protein. Pairwise genome sequence alignment was conducted by using EMBOSS Needle software (www.ebi.ac.uk/Tools/psa/emboss_needle/) with default parameters. The overall nucleotide sequences between Bat Rp-coronavirus/Shaanxi2011 and Bat Cp-coronavirus/Yunnan2011 indicated 88.7% nt identity. They shared 87.4%–89.5% nt identity with SARS-CoV, 88%–89.9% nt identity with the bat SARS-like CoV (bat SARS-CoV Rm1), and 87.6%–89.6% nt identity with the civet SARS-like CoV (civet SARS-CoV SZ16). On the other hand, comparison between the betacoronavirus genomes and human betacoronavirus (HCoV-OC43) showed only 49.9%–50.4% nt overall identity, whereas the betacoronavirus genomes and HCoV-EMC showed 52.1% nt overall identity.\nThe RNA-dependent RNA polymerase (RdRp, the 12th nonstructural protein codified to open reading frame 1a,b) is a highly conserved gene of CoVs, which is frequently used for phylogenetic comparison (8,9). MEGA5.0 (www.megasoftware.net) was used to construct the phylogenetic trees on the basis of the nucleotide sequences and deduced amino acid sequences. First, we used the MUSCLE package and default parameters (www.megasoftware.net/) to construct the alignment. The best substitution model was then evaluated with the Model Selection package implemented in MEGA5. Finally, we used the maximum-likelihood method with an appropriate model to process the phylogenetic analysis with 1,000 bootstrap replicates. We constructed a phylogenetic tree based on the nucleotide sequences of the RdRp gene to show the evolutionary relationship between these 2 betacoronaviruses and other CoVs (Figure 1). Reference CoV genome sequences were downloaded from GenBank and aligned with the fragments of the newly discovered CoVs. The RdRp genes of Bat Rp-coronavirus/Shaanxi2011 and Bat Cp-coronavirus/Yunnan2011 were highly similar, sharing 93.1% nt identity. The phylogenetic analysis demonstrated that betacoronaviruses and the bat SARS-like CoVs in our study are clustered (93.1%–93.4% nt identity) and are close in distance to SARS-CoVs (92.9%–94.8% nt identity) and civet SARS-like CoVs (93.1%–94.8% nt identity) but that bat CoV (BtCoV-HKU9) and HCoV-OC43 are placed among the relatively distant groups (65.8%–65.9% and 62.9%–63.5% nt identities with the betacoronaviruses, respectively). Therefore, collectively we called these betacoronaviruses and bat SARS-like CoVs the bat SARS-like cluster of CoVs. Bat Rp-coronavirus/Shaanxi2011 and Bat Cp-coronavirus/Yunnan2011 showed little genetic similarity (\u003c66.2%–67.3% nt identity) to HCoV-EMC.\nFigure 1 Phylogenetic tree of novel betacoronaviruses based on the nucleotide sequence of the RdRp gene. The following coronaviruses (CoVs) and GenBank accession numbers were used: bat severe acute respiratory syndrome CoV Rm1 (bat SARS-CoV Rm1; DQ412043), bat SARS-CoV Rp3 (DQ071615), bat SARS-CoV Rf1 (DQ412042), bat SARS-CoV HKU3 (DQ022305),SARS-CoV isolate Tor2/FP1–10895 (SARS-CoV Tor2; JX163925), SARS-CoV BJ182–12 (SARS-CoV BJ182; EU371564), SARS-CoV (NC004718), civet SARS-CoV SZ3 (AY304486), civet SARS-CoV SZ16 (AY304488), bat CoV HKU9 (BtCoV-HKU9; EF065513), bat CoV HKU4 (BtCoV-HKU4; EF065505), bat CoV HKU5 (BtCoV-HKU5; EF065509), human betacoronvirus 2c EMC/2012 (HCoV-EMC; JX869059), human CoV OC43 (HCoV-OC43; NC005147), HCoV-HKU1 (NC006577), bat coronavirus HKU2 (BtCoV-HKU2; NC009988), bat coronavirus 1A (BtCoV-1A; NC010437), HCoV-229E (NC002645), HCoV-NL63 (NC005831), bat CoV HKU8 (BtCoV-HKU8; NC010438), scotophilus bat CoV 512 (BtCoV-512; NC009657), avian infectious bronchitis virus (IBV; NC001451), beluga whale CoV SW1 (BWCoV; NC010646). Scale bar indicates genetic distance estimated by using TN93+G+I model implemented in MEGA5 (www.megasoftware.net). The spike proteins of CoVs are responsible for receptor binding and host species adaptation, and their genes therefore constitute one of the most variable regions within CoV genomes (10,11). The phylogenetic tree based on the amino acid sequences of spike protein (Figure 2) suggests that the selected betacoronaviruses were mainly divided into 5 clusters: SARS cluster; bat SARS-like cluster; civet SARS-like cluster; human betacoronavirus cluster; and EMC cluster. Bat Rp-coronavirus/Shaanxi2011 and Bat Cp-coronavirus/Yunnan2011 shared 89.4% aa identity in spike proteins, which consisted of 1,240 aa and 1,241 aa, respectively. The spike proteins of the CoVs in our analysis have 89.8%–92.7% aa identity with those of bat SARS-like CoVs, with substantial similarity in the receptor-binding domain. The close relationship also was observed with the SARS-CoVs (79.2%–79.4% aa identity) and civet SARS-like CoVs (78.9%–79.1% aa identity). In contrast, the human betacoronaviruses and EMC cluster formed separate clusters distinct from SARS-related CoVs that showed only 27.8%–29.4% aa and 28.8%–30.5% aa identities with the betacoronaviruses, respectively, in our analysis. The genome sequences reported here have been deposited into GenBank (accession nos. JX993987–JX993988).\nFigure 2 Phylogenetic tree of novel betacoronaviruses based on the deduced amino acid sequence of spike protein. SARS, severe acute respiratory syndrome; CoV, coronavirus; HCoV, human CoV; BtCoV, bat CoV; BWCoV, beluga whale CoV; IBV, avian infectious bronchitis. Scale bar indicates genetic distance estimated by using WAG+G+I+F model implemented in MEGA5 (www.megasoftware.net).","divisions":[{"label":"title","span":{"begin":0,"end":9}},{"label":"p","span":{"begin":10,"end":1052}},{"label":"p","span":{"begin":1053,"end":1673}},{"label":"p","span":{"begin":1674,"end":2499}},{"label":"p","span":{"begin":2500,"end":4334}},{"label":"figure","span":{"begin":4335,"end":5515}},{"label":"label","span":{"begin":4335,"end":4343}},{"label":"caption","span":{"begin":4345,"end":5515}},{"label":"p","span":{"begin":4345,"end":5515}},{"label":"p","span":{"begin":5516,"end":6794}},{"label":"label","span":{"begin":6795,"end":6803}}],"tracks":[{"project":"2_test","denotations":[{"id":"23739658-22855479-27697034","span":{"begin":205,"end":206},"obj":"22855479"},{"id":"23739658-21813714-27697035","span":{"begin":899,"end":900},"obj":"21813714"},{"id":"23739658-15613317-27697036","span":{"begin":2698,"end":2699},"obj":"15613317"},{"id":"23739658-22933277-27697037","span":{"begin":2700,"end":2701},"obj":"22933277"},{"id":"23739658-16195424-27697038","span":{"begin":5699,"end":5701},"obj":"16195424"},{"id":"23739658-17617433-27697039","span":{"begin":5702,"end":5704},"obj":"17617433"}],"attributes":[{"subj":"23739658-22855479-27697034","pred":"source","obj":"2_test"},{"subj":"23739658-21813714-27697035","pred":"source","obj":"2_test"},{"subj":"23739658-15613317-27697036","pred":"source","obj":"2_test"},{"subj":"23739658-22933277-27697037","pred":"source","obj":"2_test"},{"subj":"23739658-16195424-27697038","pred":"source","obj":"2_test"},{"subj":"23739658-17617433-27697039","pred":"source","obj":"2_test"}]},{"project":"MyTest","denotations":[{"id":"23739658-22855479-27697034","span":{"begin":205,"end":206},"obj":"22855479"},{"id":"23739658-21813714-27697035","span":{"begin":899,"end":900},"obj":"21813714"},{"id":"23739658-15613317-27697036","span":{"begin":2698,"end":2699},"obj":"15613317"},{"id":"23739658-22933277-27697037","span":{"begin":2700,"end":2701},"obj":"22933277"},{"id":"23739658-16195424-27697038","span":{"begin":5699,"end":5701},"obj":"16195424"},{"id":"23739658-17617433-27697039","span":{"begin":5702,"end":5704},"obj":"17617433"}],"namespaces":[{"prefix":"_base","uri":"https://www.uniprot.org/uniprot/testbase"},{"prefix":"UniProtKB","uri":"https://www.uniprot.org/uniprot/"},{"prefix":"uniprot","uri":"https://www.uniprot.org/uniprotkb/"}],"attributes":[{"subj":"23739658-22855479-27697034","pred":"source","obj":"MyTest"},{"subj":"23739658-21813714-27697035","pred":"source","obj":"MyTest"},{"subj":"23739658-15613317-27697036","pred":"source","obj":"MyTest"},{"subj":"23739658-22933277-27697037","pred":"source","obj":"MyTest"},{"subj":"23739658-16195424-27697038","pred":"source","obj":"MyTest"},{"subj":"23739658-17617433-27697039","pred":"source","obj":"MyTest"}]}],"config":{"attribute types":[{"pred":"source","value type":"selection","values":[{"id":"2_test","color":"#ec93c1","default":true},{"id":"MyTest","color":"#93dbec"}]}]}}