VirusRecom: an information-theory-based method for recombination detection of viral lineages and its application on SARS-CoV-2.

Genomic recombination is an important driving force for viral evolution, and recombination events have been reported for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) during the Coronavirus Disease 2019 pandemic, which significantly alter viral infectivity and transmissibility. However, it is difficult to identify viral recombination, especially for low-divergence viruses such as SARS-CoV-2, since it is hard to distinguish recombination from in situ mutation. Herein, we applied information theory to viral recombination analysis and developed VirusRecom, a program for efficiently screening recombination events on viral genome. In principle, we considered a recombination event as a transmission process of ``information'' and introduced weighted information content (WIC) to quantify the contribution of recombination to a certain region on viral genome; then, we identified the recombination regions by comparing WICs of different regions. In the benchmark using simulated data, VirusRecom showed a good balance between precision and recall compared to two competing tools, RDP5 and 3SEQ. In the detection of SARS-CoV-2 XE, XD and XF recombinants, VirusRecom providing more accurate positions of recombination regions than RDP5 and 3SEQ. In addition, we encapsulated the VirusRecom program into a command-line-interface software for convenient operation by users. In summary, we developed a novel approach based on information theory to identify viral recombination within highly similar sequences, providing a useful tool for monitoring viral evolution and epidemic control.

vsRNAfinder: a novel method for identifying high-confidence viral small RNAs from small RNA-Seq data.

Virus-encoded small RNAs (vsRNA) have been reported to play an important role in viral infection. Unfortunately, there is still a lack of an effective method for vsRNA identification. Herein, we presented vsRNAfinder, a de novo method for identifying high-confidence vsRNAs from small RNA-Seq (sRNA-Seq) data based on peak calling and Poisson distribution and is publicly available at https://github.com/ZenaCai/vsRNAfinder. vsRNAfinder outperformed two widely used methods namely miRDeep2 and ShortStack in identifying viral miRNAs with a significantly improved sensitivity. It can also be used to identify sRNAs in animals and plants with similar performance to miRDeep2 and ShortStack. vsRNAfinder would greatly facilitate effective identification of vsRNAs from sRNA-Seq data.

First detection of Tetraparvovirus ungulate 1 in diseased cattle (Chinese Simmental) from Hunan province, China.

Tetraparvovirus is an emerging parvovirus infecting a variety of mammals and humans, and associated with human diseases including severe acute respiratory infection and acute encephalitis syndrome. In the present study, a Tetraparvovirus ungulate 1 (formerly known as bovine hokovirus) strain HNU-CBY-2023 was identified and characterized from diseased Chinese Simmental from Hunan province, China. The nearly complete genome of HNU-CBY-2023 is 5346 nt in size and showed genomic identities of 85-95.5% to the known Tetraparvovirus ungulate 1 strains from GenBank, indicating a rather genetic variation. Phylogenetic and genetic divergence analyses indicated that Tetraparvovirus ungulate 1 could be divided into two genotypes (I and II), and HNU-CBY-2023 was clustered into genotype II. This study, for the first time, identified Tetraparvovirus ungulate 1 from domestic cattle from mainland China, which will be helpful to understand the prevalence and genetic diversity of Tetraparvovirus ungulate 1.

Identification and characterization of circRNAs encoded by MERS-CoV, SARS-CoV-1 and SARS-CoV-2.

The life-threatening coronaviruses MERS-CoV, SARS-CoV-1 and SARS-CoV-2 (SARS-CoV-1/2) have caused and will continue to cause enormous morbidity and mortality to humans. Virus-encoded noncoding RNAs are poorly understood in coronaviruses. Data mining of viral-infection-related RNA-sequencing data has resulted in the identification of 28 754, 720 and 3437 circRNAs encoded by MERS-CoV, SARS-CoV-1 and SARS-CoV-2, respectively. MERS-CoV exhibits much more prominent ability to encode circRNAs in all genomic regions than those of SARS-CoV-1/2. Viral circRNAs typically exhibit low expression levels. Moreover, majority of the viral circRNAs exhibit expressions only in the late stage of viral infection. Analysis of the competitive interactions of viral circRNAs, human miRNAs and mRNAs in MERS-CoV infections reveals that viral circRNAs up-regulated genes related to mRNA splicing and processing in the early stage of viral infection, and regulated genes involved in diverse functions including cancer, metabolism, autophagy, viral infection in the late stage of viral infection. Similar analysis in SARS-CoV-2 infections reveals that its viral circRNAs down-regulated genes associated with metabolic processes of cholesterol, alcohol, fatty acid and up-regulated genes associated with cellular responses to oxidative stress in the late stage of viral infection. A few genes regulated by viral circRNAs from both MERS-CoV and SARS-CoV-2 were enriched in several biological processes such as response to reactive oxygen and centrosome localization. This study provides the first glimpse into viral circRNAs in three deadly coronaviruses and would serve as a valuable resource for further studies of circRNAs in coronaviruses.

An atlas of human viruses provides new insights into diversity and tissue tropism of human viruses.

MOTIVATION: Viruses continue to threaten human health. Yet, the complete viral species carried by humans and their infection characteristics have not been fully revealed. RESULTS: This study curated an atlas of human viruses from public databases and literature, and built the Human Virus Database (HVD). The HVD contains 1131 virus species of 54 viral families which were more than twice the number of the human-infecting virus species reported in previous studies. These viruses were identified in human samples including 68 human tissues, the excreta and body fluid. The viral diversity in humans was age-dependent with a peak in the infant and a valley in the teenager. The tissue tropism of viruses was found to be associated with several factors including the viral group (DNA, RNA or reverse-transcribing viruses), enveloped or not, viral genome length and GC content, viral receptors and the virus-interacting proteins. Finally, the tissue tropism of DNA viruses was predicted using a random-forest algorithm with a middle performance. Overall, the study not only provides a valuable resource for further studies of human viruses but also deepens our understanding toward the diversity and tissue tropism of human viruses. AVAILABILITY AND IMPLEMENTATION: The HVD is available at http://computationalbiology.cn/humanVirusBase/#/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Key mutations on spike protein altering ACE2 receptor utilization and potentially expanding host range of emerging SARS-CoV-2 variants.

Increasing evidence supports inter-species transmission of SARS-CoV-2 variants from humans to domestic or wild animals during the ongoing COVID-19 pandemic, which is posing great challenges to epidemic control. Clarifying the host range of emerging SARS-CoV-2 variants will provide instructive information for the containment of viral spillover. The spike protein (S) of SARS-CoV-2 is the key determinant of receptor utilization, and therefore amino acid mutations on S will probably alter viral host range. Here, to evaluate the impact of S mutations, we tested 27 pseudoviruses of SARS-CoV-2 carrying different spike mutants by infecting Hela cells expressing different angiotensin-converting enzyme 2 (ACE2) orthologs from 20 animals. Of these 27 pseudoviruses, 20 bear single mutation and the other 7 were cloned from emerging SARS-CoV-2 variants, including D614G, Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), Lambda (B.1.429), and Mu (B.1.621). Using pseudoviral reporter assay, we identified that the substitutions of T478I and N501Y enabled the pseudovirus to utilize chicken ACE2, indicating potential infectivity to avian species. Furthermore, the S mutants of real SARS-CoV-2 variants comprising N501Y showed significantly acquired abilities to infect cells expressing mouse ACE2, indicating a critical role of N501Y in expanding SARS-CoV-2 host range. In addition, A262S and T478I significantly enhanced the utilization of various mammal ACE2. In summary, our results indicated that T478I and N501Y substitutions were two S mutations important for receptor adaption of SARS-CoV-2, potentially contributing to the spillover of the virus to many other animal hosts. Therefore, more attention should be paid to SARS-CoV-2 variants with these two mutations.

Key mutations in the spike protein of SARS-CoV-2 affecting neutralization resistance and viral internalization.

To control the ongoing COVID-19 pandemic, a variety of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccines have been developed. However, the rapid mutations of SARS-CoV-2 spike (S) protein may reduce the protective efficacy of the existing vaccines which is mainly determined by the level of neutralizing antibodies targeting S. In this study, we screened prevalent S mutations and constructed 124 pseudotyped lentiviral particles carrying these mutants. We challenged these pseudoviruses with sera vaccinated by Sinovac CoronaVac and ZF2001 vaccines, two popular vaccines designed for the initial strain of SARS-CoV-2, and then systematically assessed the susceptivity of these SARS-CoV-2 variants to the immune sera of vaccines. As a result, 14 S mutants (H146Y, V320I + S477N, V382L, K444R, L455F + S477N, L452M + F486L, F486L, Y508H, P521R, A626S, S477N + S698L, A701V, S477N + T778I, E1144Q) were found to be significantly resistant to neutralization, indicating reduced protective efficacy of the vaccines against these SARS-CoV-2 variants. In addition, F486L and Y508H significantly enhanced the utilization of human angiotensin-converting enzyme 2, suggesting a potentially elevated infectivity of these two mutants. In conclusion, our results show that some prevalent S mutations of SARS-CoV-2 reduced the protective efficacy of current vaccines and enhance the infectivity of the virus, indicating the necessity of vaccine renewal and providing direction for the development of new vaccines.

Cell membrane proteins with high N-glycosylation, high expression and multiple interaction partners are preferred by mammalian viruses as receptors.

MOTIVATION: Receptor mediated entry is the first step for viral infection. However, the question of how viruses select receptors remains unanswered. RESULTS: Here, by manually curating a high-quality database of 268 pairs of mammalian virus-host receptor interaction, which included 128 unique viral species or sub-species and 119 virus receptors, we found the viral receptors are structurally and functionally diverse, yet they had several common features when compared to other cell membrane proteins: more protein domains, higher level of N-glycosylation, higher ratio of self-interaction and more interaction partners, and higher expression in most tissues of the host. This study could deepen our understanding of virus-receptor interaction. AVAILABILITY AND IMPLEMENTATION: The database of mammalian virus-host receptor interaction is available at http://www.computationalbiology.cn: 5000/viralReceptor. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into the origin of SARS coronavirus.

A large number of SARS-related coronaviruses (SARSr-CoV) have been detected in horseshoe bats since 2005 in different areas of China. However, these bat SARSr-CoVs show sequence differences from SARS coronavirus (SARS-CoV) in different genes (S, ORF8, ORF3, etc) and are considered unlikely to represent the direct progenitor of SARS-CoV. Herein, we report the findings of our 5-year surveillance of SARSr-CoVs in a cave inhabited by multiple species of horseshoe bats in Yunnan Province, China. The full-length genomes of 11 newly discovered SARSr-CoV strains, together with our previous findings, reveals that the SARSr-CoVs circulating in this single location are highly diverse in the S gene, ORF3 and ORF8. Importantly, strains with high genetic similarity to SARS-CoV in the hypervariable N-terminal domain (NTD) and receptor-binding domain (RBD) of the S1 gene, the ORF3 and ORF8 region, respectively, were all discovered in this cave. In addition, we report the first discovery of bat SARSr-CoVs highly similar to human SARS-CoV in ORF3b and in the split ORF8a and 8b. Moreover, SARSr-CoV strains from this cave were more closely related to SARS-CoV in the non-structural protein genes ORF1a and 1b compared with those detected elsewhere. Recombination analysis shows evidence of frequent recombination events within the S gene and around the ORF8 between these SARSr-CoVs. We hypothesize that the direct progenitor of SARS-CoV may have originated after sequential recombination events between the precursors of these SARSr-CoVs. Cell entry studies demonstrated that three newly identified SARSr-CoVs with different S protein sequences are all able to use human ACE2 as the receptor, further exhibiting the close relationship between strains in this cave and SARS-CoV. This work provides new insights into the origin and evolution of SARS-CoV and highlights the necessity of preparedness for future emergence of SARS-like diseases.

Detection and genome characterization of two novel papillomaviruses and a novel polyomavirus in tree shrew (Tupaia belangeri chinensis) in China.

BACKGROUND: Papillomaviruses (PVs) and polyomaviruses (PyVs) infect diverse vertebrates including human and cause a broad spectrum of outcomes from asymptomatic infection to severe disease. There has been no PV and only one PyV detected in tree shrews, though the genomic properties of tree shrews are highly similar to those of the primates. METHODS: Swab and organ samples of tree shrews collected in the Yunnan Province of China, were tested by viral metagenomic analysis and random PCR to detect the presence of PVs and PyVs. By PCR amplification using specific primers, cloning, sequencing and assembling, genomes of two PVs and one PyV were identified in the samples. RESULTS: Two novel PVs and a novel PyV, named tree shrew papillomavirus 1 and 2 (TbelPV1 and TbelPV2) and polyomavirus 1 (TbelPyV1) were characterized in the Chinese tree shrew (Tupaia belangeri chinensis). The genomes of TbelPV1, TbelPV2, and TbelPyV1 are 7410 bp, 7526 bp, and 4982 bp in size, respectively. The TbelPV1 genome contains 7 putative open-reading frames (ORFs) coding for viral proteins E1, E2, E4, E6, E7, L1, and L2; the TbelPV2 genome contains 6 ORFs coding for viral proteins E1, E2, E6, E7, L1, and L2; and the TbelPyV1 genome codes for the typical small and large T antigens of PyV, as well as the VP1, VP2, and VP3 capsid proteins. Genomic comparison and phylogenetic analysis indicated that TbelPV1 and TbelPV2 represented 2 novel PV genera of Papillomaviridae, and TbelPyV1 represented a new species of genus Alphapolyomavirus. Our epidemiologic study indicated that TbelPV1 and TbelPV2 were both detected in oral swabs, while TbelPyV1 was detected in oral swabs and spleens. CONCLUSION: Two novel PVs (TbelPV1 and TbelPV2) and a novel PyV (TbelPyV) were discovered in tree shrews and their genomes were characterized. TbelPV1, TbelPV2, and TbelPyV1 have the highest similarity to Human papillomavirus type 63, Ursus maritimus papillomavirus 1, and Human polyomavirus 9, respectively. TbelPV1 and TbelPV2 only showed oral tropism, while TbelPyV1 showed oral and spleen tropism.

Genomic characterization and pathogenicity of porcine deltacoronavirus strain CHN-HG-2017 from China.

Porcine deltacoronavirus (PDCoV) was first detected in Hong Kong and has recently spread to many countries around the world. PDCoV causes acute diarrhea and vomiting in pigs, resulting in significant economic losses in the global pork industry. In this study, a Chinese PDCoV strain, designated CHN-HG-2017, was isolated from feces of a suckling piglet with severe watery diarrhea on a farm located in central China. Subsequently, the virus was identified by an indirect immunofluorescence assay and electron microscopy. A nucleotide sequence alignment showed that the whole genome of CHN-HG-2017 is 97.6%-99.1% identical to other PDCoV strains. Analysis of potential recombination sites showed that CHN-HG-2017 is a possible recombinant originating from the strains CH/SXD1/2015 and Vietnam/HaNoi6/2015. Furthermore, the pathogenicity of this recombinant PDCoV strain was investigated in 5-day-old piglets by oral inoculation. The challenged piglets developed typical symptoms, such as vomiting, anorexia, diarrhea and lethargy, from 1 to 7 days post-inoculation (DPI). Viral shedding was detected in rectal swabs until 14 DPI in the challenged piglets. Interestingly, high titers of virus-neutralizing antibodies in sera were detected at 21 DPI. Tissues of small intestines from CHN-HG-2017-infected piglets at 4 DPI displayed significant macroscopic and microscopic lesions with clear viral antigen expression. Our analysis of the full genome sequence of a recombinant PDCoV and its virulence in suckling piglets might provide new insights into the pathogenesis of PDCoV and facilitate further investigation of this newly emerged pathogen.

Genomic characterization of a novel astrovirus identified in Amur tigers from a zoo in China.

Astroviruses (AstVs) have a very wide range of hosts and are associated with enteric and extra-enteric disease in mammals and birds. Cross-species transmission of AstVs has been observed frequently. In the present study, the genome of a novel astrovirus from Amur tigers (Panthera tigris) from a zoo in China was characterized and was found to have the typical genomic features of other mammal AstVs. It showed the highest nucleotide sequence similarity (46.1-87.3% identity) to AstVs from cats, indicating a close phylogenetic relationship and possible cross-species transmission between them. To our knowledge, this is the first identification and characterization of AstV from tigers, and this virus is the third astrovirus identified in hosts of the family Felidae. The results of this study will be helpful for understanding the origin, genetic diversity, and cross-species transmission of AstV.

Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor.

The 2002-3 pandemic caused by severe acute respiratory syndrome coronavirus (SARS-CoV) was one of the most significant public health events in recent history. An ongoing outbreak of Middle East respiratory syndrome coronavirus suggests that this group of viruses remains a key threat and that their distribution is wider than previously recognized. Although bats have been suggested to be the natural reservoirs of both viruses, attempts to isolate the progenitor virus of SARS-CoV from bats have been unsuccessful. Diverse SARS-like coronaviruses (SL-CoVs) have now been reported from bats in China, Europe and Africa, but none is considered a direct progenitor of SARS-CoV because of their phylogenetic disparity from this virus and the inability of their spike proteins to use the SARS-CoV cellular receptor molecule, the human angiotensin converting enzyme II (ACE2). Here we report whole-genome sequences of two novel bat coronaviruses from Chinese horseshoe bats (family: Rhinolophidae) in Yunnan, China: RsSHC014 and Rs3367. These viruses are far more closely related to SARS-CoV than any previously identified bat coronaviruses, particularly in the receptor binding domain of the spike protein. Most importantly, we report the first recorded isolation of a live SL-CoV (bat SL-CoV-WIV1) from bat faecal samples in Vero E6 cells, which has typical coronavirus morphology, 99.9% sequence identity to Rs3367 and uses ACE2 from humans, civets and Chinese horseshoe bats for cell entry. Preliminary in vitro testing indicates that WIV1 also has a broad species tropism. Our results provide the strongest evidence to date that Chinese horseshoe bats are natural reservoirs of SARS-CoV, and that intermediate hosts may not be necessary for direct human infection by some bat SL-CoVs. They also highlight the importance of pathogen-discovery programs targeting high-risk wildlife groups in emerging disease hotspots as a strategy for pandemic preparedness.

Characterization of a Minimal Type of Promoter Containing the -10 Element and a Guanine at the -14 or -13 Position in Mycobacteria.

Three key promoter elements, i.e., -10, -35, and T-15G-14N, are recognized by the σ subunit of RNA polymerase. Among them, promoters with the -10 element and either -35 or T-15G-14N are known to initiate transcription efficiently, but recent systematic analyses have identified a large group of promoters in Mycobacterium tuberculosis that contain only a -10 consensus. How these promoters initiate transcription remains poorly understood. Here, we show that promoters containing the -10 element and an upstream G located at the -14 or -13 position can successfully initiate transcription in mycobacteria. Importantly, this new type of promoter is active in the absence of other promoter consensuses, suggesting that it is a minimal promoter type. Mutation of the upstream G in promoters decreased the efficiencies of their binding with RNA polymerase and their abilities to initiate transcription in both in vitro and in vivo analyses. A glutamic acid in σ region 3.0 is essential for recognizing G-14 and G-13 and is conserved in both principal and principal-like σ factors in mycobacteria, indicating that recognition of this minimal type of promoter might be a common mechanism for transcription initiation. Consistently, more than 70% of the identified promoters in M. tuberculosis contained G-14 or G-13 upstream of the conserved -10 element, and thousands of promoters in representative mycobacterial species have been predicted using the -10 consensus and G-14 or G-13 Altogether, our study presents a universal mechanism for transcription initiation from a minimal promoter in mycobacteria, which might also be applicable to other bacteria.IMPORTANCE In contrast to the detailed information for recognizing classic promoters in the model organism Escherichia coli, very little is known about how transcription is initiated in the human pathogen Mycobacterium tuberculosis In this study, we characterized a new type of promoter in mycobacteria that requires only a -10 consensus and an upstream G-14 or G-13 Residues important for recognizing the -10 element and the upstream G are conserved in σA and σB from mycobacterial species. According to such features, thousands of promoters in mycobacteria can be predicted using the -10 consensus and G-14 or G-13, which suggests that transcription from this new type of promoter might be widespread. Our findings provide insightful information for characterizing promoters in mycobacteria.

Novel bat adenoviruses with low G+C content shed new light on the evolution of adenoviruses.

Bats have been reported to carry diverse adenoviruses. However, most bat adenoviruses have been identified on the basis of partial genome sequences, and knowledge on the evolution of bat adenoviruses remains limited. In this study, we isolated and characterized four novel adenoviruses from two distinct bat species, and their full-length genomes were sequenced. Sequence analysis revealed that these isolates represented three distinct species of the genus Mastadenovirus. However, all isolates had an exceptionally low G+C content and relatively short genomes compared with other known mastadenoviruses. We further analysed the relationships among the G+C content, 5'-C-phosphate-G-3' (CpG) representation and genome size in the family Adenoviridae. Our results revealed that the CpG representation in adenoviral genomes depends primarily on the level of methylation, and the genome size displayed significant positive correlations with both G+C content and CpG representation. Since ancestral adenoviruses are believed to have contained short genomes, those probably had a low G+C content, similar to the genomes of these bat strains. Our results suggest that bats are important natural reservoirs for adenoviruses and play important roles in the evolution of adenoviruses.

Bat Severe Acute Respiratory Syndrome-Like Coronavirus WIV1 Encodes an Extra Accessory Protein, ORFX, Involved in Modulation of the Host Immune Response.

UNLABELLED: Bats harbor severe acute respiratory syndrome (SARS)-like coronaviruses (SL-CoVs) from which the causative agent of the 2002-2003 SARS pandemic is thought to have originated. However, despite the fact that a large number of genetically diverse SL-CoV sequences have been detected in bats, only two strains (named WIV1 and WIV16) have been successfully cultured in vitro These two strains differ from SARS-CoV only in containing an extra open reading frame (ORF) (named ORFX), between ORF6 and ORF7, which has no homology to any known protein sequences. In this study, we constructed a full-length cDNA clone of SL-CoV WIV1 (rWIV1), an ORFX deletion mutant (rWIV1-ΔX), and a green fluorescent protein (GFP)-expressing mutant (rWIV1-GFP-ΔX). Northern blotting and fluorescence microscopy indicate that ORFX was expressed during WIV1 infection. A virus infection assay showed that rWIV1-ΔX replicated as efficiently as rWIV1 in Vero E6, Calu-3, and HeLa-hACE2 cells. Further study showed that ORFX could inhibit interferon production and activate NF-κB. Our results demonstrate for the first time that the unique ORFX in the WIV1 strain is a functional gene involving modulation of the host immune response but is not essential for in vitro viral replication. IMPORTANCE: Bats harbor genetically diverse SARS-like coronaviruses (SL-CoVs), and some of them have the potential for interspecies transmission. A unique open reading frame (ORFX) was identified in the genomes of two recently isolated bat SL-CoV strains (WIV1 and -16). It will therefore be critical to clarify whether and how this protein contributes to virulence during viral infection. Here we revealed that the unique ORFX is a functional gene that is involved in the modulation of the host immune response but is not essential for in vitro viral replication. Our results provide important information for further exploration of the ORFX function in the future. Moreover, the reverse genetics system we constructed will be helpful for study of the pathogenesis of this group of viruses and to develop therapeutics for future control of emerging SARS-like infections.

Detection of alpha- and betacoronaviruses in rodents from Yunnan, China.

BACKGROUND: Rodents represent the most diverse mammals on the planet and are important reservoirs of human pathogens. Coronaviruses infect various animals, but to date, relatively few coronaviruses have been identified in rodents worldwide. The evolution and ecology of coronaviruses in rodent have not been fully investigated. RESULTS: In this study, we collected 177 intestinal samples from thress species of rodents in Jianchuan County, Yunnan Province, China. Alphacoronavirus and betacoronavirus were detected in 23 rodent samples from three species, namely Apodemus chevrieri (21/98), Eothenomys fidelis (1/62), and Apodemus ilex (1/17). We further characterized the full-length genome of an alphacoronavirus from the A. chevrieri rat and named it as AcCoV-JC34. The AcCoV-JC34 genome was 27,649 nucleotides long and showed a structure similar to the HKU2 bat coronavirus. Comparing the normal transcription regulatory sequence (TRS), 3 variant TRS sequences upstream the spike (S), ORF3, and ORF8 genes were found in the genome of AcCoV-JC34. In the conserved replicase domains, AcCoV-JC34 was most closely related to Rattus norvegicus coronavirus LNRV but diverged from other alphacoronaviruses, indicating that AcCoV-JC34 and LNRV may represent a novel alphacoronavirus species. However, the S and nucleocapsid proteins showed low similarity to those of LRNV, with 66.5 and 77.4% identities, respectively. Phylogenetic analysis revealed that the S genes of AcCoV-JC34, LRNV, and HKU2 formed a distinct lineage with all known coronaviruses. CONCLUSIONS: Both alphacoronaviruses and betacoronaviruses were detected in Apodemus chevrieri in the Yunnan Province of China, indicating that Apodemus chevrieri is an important host for coronavirus. Several new features were identified in the genome of an Apodemus chevrieri coronavirus. The phylogenetic distance to other coronaviruses suggests a variable origin and evolutionary route of the S genes of AcCoV-JC34, LRNV, and HKU2. These results indicate that the diversity of rodent coronaviruses is much higher than previously expected. Further surveillance and functional studies of these coronaviruses will help to better understand the importance of rodent as host for coronaviruses.

Detection and genome characterization of four novel bat hepadnaviruses and a hepevirus in China.

BACKGROUND: In recent years, novel hepadnaviruses, hepeviruses, hepatoviruses, and hepaciviruses have been discovered in various species of bat around the world, indicating that bats may act as natural reservoirs for these hepatitis viruses. In order to further assess the distribution of hepatitis viruses in bat populations in China, we tested the presence of these hepatitis viruses in our archived bat liver samples that originated from several bat species and various geographical regions in China. METHODS: A total of 78 bat liver samples (involving two families, five genera, and 17 species of bat) were examined using nested or heminested reverse transcription PCR (RT-PCR) with degenerate primers. Full-length genomic sequences of two virus strains were sequenced followed by phylogenetic analyses. RESULTS: Four samples were positive for hepadnavirus, only one was positive for hepevirus, and none of the samples were positive for hepatovirus or hepacivirus. The hepadnaviruses were discovered in the horseshoe bats, Rhinolophus sinicus and Rhinolophus affinis, and the hepevirus was found in the whiskered bat Myotis davidii. The full-length genomic sequences were determined for one of the two hepadnaviruses identified in R. sinicus (designated BtHBVRs3364) and the hepevirus (designated BtHEVMd2350). A sequence identity analysis indicated that BtHBVRs3364 had the highest degree of identity with a previously reported hepadnavirus from the roundleaf bat, Hipposideros pomona, from China, and BtHEVMd2350 had the highest degree of identity with a hepevirus found in the serotine bat, Eptesicus serotinus, from Germany, but it exhibited high levels of divergence at both the nucleotide and the amino acid levels. CONCLUSIONS: This is the first study to report that the Chinese horseshoe bat and the Chinese whiskered bat have been found to carry novel hepadnaviruses and a novel hepevirus, respectively. The discovery of BtHBVRs3364 further supports the significance of host switches evolution while opposing the co-evolutionary theory associated with hepadnaviruses. According to the latest criterion of the International Committee on Taxonomy of Viruses (ICTV), we hypothesize that BtHEVMd2350 represents an independent genotype within the species Orthohepevirus D of the family Hepeviridae.

Novel bat adenoviruses with an extremely large E3 gene.

Bats carry diverse RNA viruses, some of which are responsible for human diseases. Compared to bat-borne RNA viruses, relatively little information is known regarding bat-borne DNA viruses. In this study, we isolated and characterized three novel bat adenoviruses (BtAdV WIV9-11) from Rhinolophus sinicus. Their genomes, which are highly similar to each other but distinct from those of previously sequenced adenoviruses (AdVs), are 37 545, 37 566 and 38 073 bp in size, respectively. An unusually large E3 gene was identified in their genomes. Phylogenetic and taxonomic analyses suggested that these isolates represent a distinct species of the genus Mastadenovirus. Cell susceptibility assays revealed a broad cell tropism for these isolates, indicating that they have a potentially wide host range. Our results expand the understanding of genetic diversity of bat AdVs.

Coronavirus nsp10/nsp16 Methyltransferase Can Be Targeted by nsp10-Derived Peptide In Vitro and In Vivo To Reduce Replication and Pathogenesis.

UNLABELLED: The 5' cap structures of eukaryotic mRNAs are important for RNA stability and protein translation. Many viruses that replicate in the cytoplasm of eukaryotes have evolved 2'-O-methyltransferases (2'-O-MTase) to autonomously modify their mRNAs and carry a cap-1 structure (m7GpppNm) at the 5' end, thereby facilitating viral replication and escaping innate immune recognition in host cells. Previous studies showed that the 2'-O-MTase activity of severe acute respiratory syndrome coronavirus (SARS-CoV) nonstructural protein 16 (nsp16) needs to be activated by nsp10, whereas nsp16 of feline coronavirus (FCoV) alone possesses 2'-O-MTase activity (E. Decroly et al., J Virol 82:8071-8084, 2008, http://dx.doi.org/10.1128/JVI.00407-08; M. Bouvet et al., PLoS Pathog 6:e1000863, 2010, http://dx.doi.org/10.1371/journal.ppat.1000863; E. Decroly et al., PLoS Pathog 7:e1002059, 2011, http://dx.doi.org/10.1371/journal.ppat.1002059; Y. Chen et al., PLoS Pathog 7:e1002294, 2011, http://dx.doi.org/10.1371/journal.ppat.1002294) . In this study, we demonstrate that stimulation of nsp16 2'-O-MTase activity by nsp10 is a universal and conserved mechanism in coronaviruses, including FCoV, and that nsp10 is functionally interchangeable in the stimulation of nsp16 of different coronaviruses. Based on our current and previous studies, we designed a peptide (TP29) from the sequence of the interaction interface of mouse hepatitis virus (MHV) nsp10 and demonstrated that the peptide inhibits the 2'-O-MTase activity of different coronaviruses in biochemical assays and the viral replication in MHV infection and SARS-CoV replicon models. Interestingly, the peptide TP29 exerted robust inhibitory effects in vivo in MHV-infected mice by impairing MHV virulence and pathogenesis through suppressing virus replication and enhancing type I interferon production at an early stage of infection. Therefore, as a proof of principle, the current results indicate that coronavirus 2'-O-MTase activity can be targeted in vitro and in vivo. IMPORTANCE: Coronaviruses are important pathogens of animals and human with high zoonotic potential. SARS-CoV encodes the 2'-O-MTase that is composed of the catalytic subunit nsp16 and the stimulatory subunit nsp10 and plays an important role in virus genome replication and evasion from innate immunity. Our current results demonstrate that stimulation of nsp16 2'-O-MTase activity by nsp10 is a common mechanism for coronaviruses, and nsp10 is functionally interchangeable in the stimulation of nsp16 among different coronaviruses, which underlies the rationale for developing inhibitory peptides. We demonstrate that a peptide derived from the nsp16-interacting domain of MHV nsp10 could inhibit 2'-O-MTase activity of different coronaviruses in vitro and viral replication of MHV and SARS-CoV replicon in cell culture, and it could strongly inhibit virus replication and pathogenesis in MHV-infected mice. This work makes it possible to develop broad-spectrum peptide inhibitors by targeting the nsp16/nsp10 2'-O-MTase of coronaviruses.