Identification and discovery of viruses using next-generation sequencing technology is a fast-developing area with potential wide application in clinical diagnostics, public health monitoring and novel virus discovery. However, tremendous sequence data from NGS study has posed great challenge both in accuracy and velocity for application of NGS study. Here we describe VIP (“Virus Identification Pipeline”), a one-touch computational pipeline for virus identification and discovery from metagenomic NGS data. VIP performs the following steps to achieve its goal: (i) map and filter out background-related reads, (ii) extensive classification of reads on the basis of nucleotide and remote amino acid homology, (iii) multiple k-mer based de novo assembly and phylogenetic analysis to provide evolutionary insight. We validated the feasibility and veracity of this pipeline with sequencing results of various types of clinical samples and public datasets. VIP has also contributed to timely virus diagnosis (~10 min) in acutely ill patients, demonstrating its potential in the performance of unbiased NGS-based clinical studies with demand of short turnaround time. VIP is released under GPLv3 and is available for free download at: https://github.com/keylabivdc/VIP.

On 26 May 2015, an imported Middle East respiratory syndrome coronavirus (MERS-CoV) was identified in Guangdong Province, China, and found to be closely related to the MERS-CoV strain prevalent in South Korea. The full genome of the ChinaGD01 strain was sequenced and analyzed to investigate the epidemiology and evolution of MERS-CoV circulating in South Korea and China.

Background: Mosquitoes host and transmit numerous arthropod-borne viruses (arboviruses) that cause disease in both humans and animals. Effective surveillance of virome profiles in mosquitoes is vital to the prevention and control of mosquito-borne diseases in northwestern China, where epidemics occur frequently. Methods: Mosquitoes were collected in the Shaanxi-Gansu-Ningxia region (Shaanxi Province, Gansu Province, and Ningxia Hui Autonomous Region) of China from June to August 2019. Morphological methods were used for taxonomic identification of mosquito species. High-throughput sequencing and metagenomic analysis were used to characterize mosquito viromes. Results: A total of 22,959 mosquitoes were collected, including Culex pipiens (45.7%), Culex tritaeniorhynchus (40.6%), Anopheles sinensis (8.4%), Aedes (5.2%), and Armigeres subalbatus (0.1%). In total, 3,014,183 (0.95% of clean reads) viral sequences were identified and assigned to 116 viral species (including pathogens such as Japanese encephalitis virus and Getah virus) in 31 viral families, including Flaviviridae, Togaviridae, Phasmaviridae, Phenuiviridae, and some unclassified viruses. Mosquitoes collected in July (86 species in 26 families) showed greater viral diversity than those from June and August. Culex pipiens (69 species in 25 families) and Culex tritaeniorhynchus (73 species in 24 families) carried more viral species than Anopheles sinensis (50 species in 19 families) or Aedes (38 species in 20 families) mosquitoes. Conclusion: Viral diversity and abundance were affected by mosquito species and collection time. The present study elucidates the virome compositions of various mosquito species in northwestern China, improving the understanding of virus transmission dynamics for comparison with those of disease outbreaks.

We isolated a novel virus strain (YN12246) from Culicoides spp. specimens collected at the China-Laos-Myanmar border in southern Yunnan Province. This virus had a cytopathic effect (CPE) on both insect cells (C6/36) and mammalian cells (BHK-21). Electron microscopy revealed the structure of the virions to be spherical with a diameter of 75 nm. Polyacrylamide gel analysis demonstrated that the viral genome consisted of 10 segments of double-stranded RNA (dsRNA), with a distribution pattern of 3-3-3-1. The coding sequences of 9 genome segments of YN12246 (Seg1, Seg3-Seg10) were obtained by high-throughput sequencing and Sanger sequencing. Comparisons of conserved genome segments 1 and 3 (Seg1 and Seg3), encoding the polymerase-VP1 and sub-core T2 protein, respectively, showed that YN12246 groups with the Culicoides-borne orbiviruses. The highest levels of sequence identity were detected between YN12246 and Tibet orbivirus (TIBOV), indicating that they belong to the same virus species (with amino acid identity of 98.8% and 96.4% for the polymerase and T2 protein, respectively). The data presented here confirm that YN12246 is a member of the TIBOV species, which was first isolated from mosquitoes in 2009. This is the first report of the isolation of TIBOV from Culicoides.

High rate of cardiovascular disease (CVD) has been reported among patients with coronavirus disease 2019 (COVID-19). Importantly, CVD, as one of the comorbidities, could also increase the risks of the severity of COVID-19. Here we identified phospholipase A2 group VII (PLA2G7), a well-studied CVD biomarker, as a hub gene in COVID-19 though an integrated hypothesis-free genomic analysis on nasal swabs (n = 486) from patients with COVID-19. PLA2G7 was further found to be predominantly expressed by proinflammatory macrophages in lungs emerging with progression of COVID-19. In the validation stage, RNA level of PLA2G7 was identified in nasal swabs from both COVID-19 and pneumonia patients, other than health individuals. The positive rate of PLA2G7 were correlated with not only viral loads but also severity of pneumonia in non-COVID-19 patients. Serum protein levels of PLA2G7 were found to be elevated and beyond the normal limit in COVID-19 patients, especially among those re-positive patients. We identified and validated PLA2G7, a biomarker for CVD, was abnormally enhanced in COVID-19 at both nucleotide and protein aspects. These findings provided indications into the prevalence of cardiovascular involvements seen in patients with COVID-19. PLA2G7 could be a potential prognostic and therapeutic target in COVID-19.