Virome comparisons in wild-diseased and healthy captive giant pandas

Diverse papillomaviruses identified from Antarctic fur seals, leopard seals and Weddell seals from the Antarctic

Papillomaviruses (family Papillomaviridae) are non-enveloped, circular, double-stranded DNA viruses known to infect squamous and mucosal epithelial cells. In the family Papillomaviridae there are 53 genera and 133 viral species whose members infect a variety of mammalian, avian, reptilian, and fish species. Within the Antarctic context, papillomaviruses (PVs) have been identified in Adélie penguins (Pygoscelis adeliae, 2 PVs), Weddell seals (Leptonychotes weddellii, 7 PVs), and emerald notothen (Trematomus bernacchii, 1 PV) in McMurdo Sound and Ross Island in eastern Antarctica. Here we identified 13 diverse PVs from buccal swabs of Antarctic fur seals (Arctocephalus gazella, 2 PVs) and leopard seal (Hydrurga leptonyx, 3 PVs) in western Antarctica (Antarctic Peninsula), and vaginal and nasal swabs of Weddell seals (8 PVs) in McMurdo Sound. These PV genomes group into four genera representing 11 new papillomavirus types, of which five are from two Antarctic fur seals and a leopard seal and six from Weddell seals.

Differences between the intestinal microbial communities of healthy dogs from plateau and those of plateau dogs infected with Echinococcus

OBJECTIVE

Cystic echinococcosis (CE) represents a profoundly perilous zoonotic disease. The advent of viral macrogenomics has facilitated the exploration of hitherto uncharted viral territories. In the scope of this investigation, our objective is to scrutinize disparities in the intestinal microbiotic ecosystems of canines dwelling in elevated terrains and those afflicted by Echinococcus infection, employing the tool of viral macrogenomics.

METHODS

In this study, we collected a comprehensive total of 1,970 fecal samples from plateau dogs infected with Echinococcus, as well as healthy control plateau dogs from the Yushu and Guoluo regions in the highland terrain of China. These samples were subjected to viral macrogenomic analysis to investigate the viral community inhabiting the canine gastrointestinal tract.

RESULTS

Our meticulous analysis led to the identification of 136 viral genomic sequences, encompassing eight distinct viral families.

CONCLUSION

The outcomes of this study hold the potential to enhance our comprehension of the intricate interplay between hosts, parasites, and viral communities within the highland canine gut ecosystem. Through the examination of phage presence, it may aid in early detection or assessment of infection severity, providing valuable insights into Echinococcus infection and offering prospects for potential treatment strategies.

Canis Familiaris Papillomavirus Type 26: A Novel Papillomavirus of Dogs and the First Canine Papillomavirus within the Omegapapillomavirus Genus

Domestic dogs are currently recognized as being infected by 25 different canine papillomavirus (CPV) types classified into three genera. A short sequence from a novel CPV type was amplified, along with CPV1, from a papilloma (wart) from the mouth of a dog. The entire 7499 bp genome was amplified, and CPV26 contained putative coding regions that were predicted to produce four early proteins and two late ones. The ORF L1 showed less than 62% similarity for all previously sequenced CPV types but over 69% similarity to multiple Omegapapillomavirus types from a variety of Caniform species including the giant panda, Weddel seal, and polar bear. Phylogenetic analysis confirmed CPV26 clusters within the Omegapapillomavirus genus. Specific primers were used to investigate the presence of CPV26 DNA within a series of 37 canine proliferative lesions. CPV26 DNA was amplified from one lesion, a cutaneous papilloma that also contained CPV6. This is the first time a PV type within the Omegapapillomavirus genus has been detected in a non-domestic species and this provides evidence that the omegapapillomaviruses infected a common ancestor of, and then co-evolved with, the Caniform species. Whether CPV26 causes disease is uncertain, but the absence of an E7 protein may suggest low pathogenicity.

Giant pandas in captivity undergo short-term adaptation in nerve-related pathways

BACKGROUND

Behaviors in captive animals, including changes in appetite, activity level, and social interaction, are often seen as adaptive responses. However, these behaviors may become progressively maladaptive, leading to stress, anxiety, depression, and other negative reactions in animals.

RESULTS

In this study, we investigated the whole-genome sequencing data of 39 giant panda individuals, including 11 in captivity and 28 in the wild. To eliminate the mountain range effect and focus on the factor of captivity only, we first performed a principal component analysis. We then enumerated the 21,474,180 combinations of wild giant pandas (11 chosen from 28) and calculated their distances from the 11 captive individuals. The 11 wild individuals with the closest distances were used for the subsequent analysis. The linkage disequilibrium (LD) patterns demonstrated that the population was almost eliminated. We identified 505 robust selected genomic regions harboring at least one SNP, and the absolute frequency difference was greater than 0.6 between the two populations. GO analysis revealed that genes in these regions were mainly involved in nerve-related pathways. Furthermore, we identified 22 GO terms for which the selection strength significantly differed between the two populations, and there were 10 nerve-related pathways among them. Genes in the differentially abundant regions were involved in nerve-related pathways, indicating that giant pandas in captivity underwent minor genomic selection. Additionally, we investigated the relationship between genetic variation and chromatin conformation structures. We found that nucleotide diversity (θπ) in the captive population was correlated with chromatin conformation structures, which included A/B compartments, topologically associated domains (TADs) and TAD-cliques. For each GO term, we then compared the expression level of genes regulated by the above four factors (AB index, TAD intactness, TAD clique and PEI) with the corresponding genomic background. The retained 10 GO terms were all coordinately regulated by the four factors, and three of them were associated with nerve-related pathways.

CONCLUSIONS

This study revealed that giant pandas in captivity undergo short-term adaptation in nerve-related pathways. Furthermore, it provides new insights into the molecular mechanism of gene expression regulation under short-term adaptation to environmental change.

A survey of fecal virome and bacterial community of the diarrhea-affected cattle in northeast China reveals novel disease-associated ecological risk factors

Limited information on the virome and bacterial community hampers our ability to discern systemic ecological risk factors that cause cattle diarrhea, which has become a pressing issue in the control of disease. A total of 110 viruses, 1,011 bacterial genera, and 322 complete viral genomes were identified from 70 sequencing samples mixed with 1,120 fecal samples from 58 farms in northeast China. For the diarrheic samples, the identified virome and bacterial community varied in terms of composition, abundance, diversity, and geographic distribution in relation to different disease-associated ecological factors; the abundance of identified viruses and bacteria was significantly correlated with the host factors of clinical status, cattle type, and age, and with environmental factors such as aquaculture model and geographical location (P < 0.05); a significant interaction occurred between viruses and viruses, bacteria and bacteria, as well as between bacteria and viruses (P < 0.05). The abundance of SMB53, Butyrivibrio, Facklamia, Trichococcus, and Turicibacter was significantly correlated with the health status of cattle (P < 0.05). The proportion of BRV, BCoV, BKV, BToV, BoNoV, BoNeV, BoAstV, BEV, BoPV, and BVDV in 1,120 fecal samples varied from 1.61% to 12.05%. A series of significant correlations were observed between the prevalence of individual viruses and the disease-associated ecological factors. A genome-based phylogenetic analysis revealed high variability of 10 bovine enteric viruses. The bovine hungarovirus was initially identified in both dairy and beef cattle in China. This study elucidates the fecal virome and bacterial community signatures of cattle affected by diarrhea, and reveals novel disease-associated ecological risk factors, including cattle type, cattle age, aquaculture model, and geographical location.IMPORTANCEThe lack of data on the virome and bacterial community restricts our capability to recognize ecological risk factors for bovine diarrhea disease, thereby hindering our overall comprehension of the disease's cause. In this study, we found that, for the diarrheal samples, the identified virome and bacterial community varied in terms of composition, abundance, diversity, configuration, and geographic distribution in relation to different disease-associated ecological factors. A series of significant correlations were observed between the prevalence of individual viruses and the disease-associated ecological factors. Our study aims to uncover novel ecological risk factors of bovine diarrheal disease by examining the pathogenic microorganism-host-environment disease ecology, thereby providing a new perspective on the control of bovine diarrheal diseases.

Gut microbial community structure and function of Przewalski's horses varied across reintroduced sites in China

Host-associated microbiota can significantly impact host fitness. Therefore, naturally occurring variations in microbiota may influence the health and persistence of their hosts. This finding is particularly important in reintroduced animals, as they typically experience habitat changes during translocations. However, little is known about how microbiomes are altered in response to conservation translocation. Here, we accessed the gut microbiome of Przewalski's horse (Equus przewalskii) populations in China from three nature reserves (i.e. Xinjiang Kalamaili Nature Reserve, KNR; Dunhuang Xihu National Nature Reserve, DXNNR; and Anxi Extreme-arid Desert Nature Reserve, AENR) using 16s rRNA gene and metagenome sequencing. The results showed that the microbial composition and function differed significantly across locations, while a subset of core taxa was consistently present in most of the samples. The abundance of genes encoding microbe-produced enzymes involved in the metabolism of carbohydrates, especially for glycoside hydrolases, was significantly higher in open-spaced KNR populations than in more confined AENR individuals. This study offers detailed and significant differential characters related to the microbial community and metabolic pathways in various reintroduced sites of Przewalski's horse, which might provide a basis for future microecological and conservation research on endangered reintroduced animals.

Virome of high-altitude canine digestive tract and genetic characterization of novel viruses potentially threatening human health

The majority of currently emerging infectious illnesses are zoonotic infections, which have caused serious public health and economic implications. The development of viral metagenomics has helped us to explore unknown viruses. We collected 1,970 canine feces from Yushu and Guoluo in the plateau region of China for this study to do a metagenomics analysis of the viral community of the canine digestive tract. Our analysis identified 203 novel viruses, classified into 11 known families and 2 unclassified groups. These viruses include the hepatitis E virus, first identified in dogs, and the astrovirus, coronavirus, polyomavirus, and others. The relationship between the newly identified canine viruses and known viruses was investigated through the use of phylogenetic analysis. Furthermore, we demonstrated the cross-species transmission of viruses and predicted new viruses that may cause diseases in both humans and animals, providing technical support for the prevention and control of diseases caused by environmental pollution viruses. IMPORTANCE Most emerging infectious diseases are due to zoonotic disease agents. Because of their effects on the security of human or animal life, agriculture production, and food safety, zoonotic illnesses and livestock diseases are of worldwide significance. Because dogs are closely related to humans and domestic animals, they serve as one of the important links in the transmission of zoonotic and livestock diseases. Canines can contaminate the environment in which humans live such as water and soil through secretions, potentially altering the human gut microbiota or causing diseases. Our study enriched the viral community in the digestive tract microbiome of dogs and found types of viruses that threaten human health, providing technical support for the prevention and control of early warning of diseases caused by environmental contaminant viruses.

Identification of Multiple Novel Viruses in Fecal Samples of Black-Necked Cranes Using Viral Metagenomic Methods

The black-necked crane is the only species of crane that lives in the high-altitude region of the Tibet Plateau. At present, there is little research on viral diseases of the black-necked crane (Grus nigricollis). In this study, a viral metagenomic approach was employed to investigate the fecal virome of black-necked cranes in Saga County, Shigatse City, Tibet, China. The identified virus families carried by black-necked cranes mainly include Genomoviridae, Parvoviridae, and Picornaviridae. The percentages of sequence reads belonging to these three virus families were 1.6%, 3.1%, and 93.7%, respectively. Among them, one genome was characterized as a novel species in the genus Grusopivirus of the family Picornaviridae, four new parvovirus genomes were obtained and classified into four different novel species within the genus Chaphamaparvovirus of the subfamily Hamaparvovirinae, and four novel genomovirus genomes were also acquired and identified as members of three different species, including Gemykroznavirus haeme1, Gemycircularvirus ptero6, and Gemycircularvirus ptero10. All of these viruses were firstly detected in fecal samples of black-necked cranes. This study provides valuable information for understanding the viral community composition in the digestive tract of black-necked cranes in Tibet, which can be used for monitoring, preventing, and treating potential viral diseases in black-necked cranes.

Metagenomic analysis of herbivorous mammalian viral communities in the Northwest Plateau

BACKGROUND

Mammals are potential hosts for many infectious diseases. However, studies on the viral communities of herbivorous mammals in the Northwest Plateau are limited. Here, we studied the viral communities of herbivorous mammals in the Northwest Plateau using virus metagenomic analysis to analyze and compare the viral community composition of seven animal species.

RESULTS

By library construction and next-generation sequencing, contigs and singlets reads with similar viral sequences were classified into 24 viral families. Analyzed from the perspective of sampling areas, the virus community composition was relatively similar in two areas of Wuwei and Jinchang, Gansu Province. Analyzed from the perspective of seven animal species, the viral reads of seven animal species were mostly ssDNA and dominated by CRESS-DNA viruses. Phylogenetic analysis based on viral marker genes indicated that CRESS-DNA viruses and microviruses have high genetic diversity. In addition to DNA viruses, nodaviruses, pepper mild mottle viruses and picornaviruses were RNA viruses that we performed by phylogenetic analysis. The CRESS-DNA viruses and nodaviruses are believed to infect plants and insects, and microviruses can infect bacteria, identifying that they were likely from the diet of herbivorous mammals. Notably, two picornaviruses were identified from red deer and wild horse, showing that the picornavirus found in red deer had the relatively high similarity with human hepatitis A virus, and the picornavirus carried by wild horse could potentially form a new species within the Picornaviridae family.

CONCLUSIONS

This study explored the herbivorous mammalian virus community in the Northwest Plateau and the genetic characteristics of viruses that potentially threaten human health. It reveals the diversity and stability of herbivorous mammalian virus communities in the Northwest Plateau and helps to expand our knowledge of various herbivorous mammalian potentially pathogenic viruses.

Blood virome of patients with traumatic sepsis

Sepsis is one of the possible outcomes of severe trauma, and it poses a dire threat to human life, particularly in immunocompromised people. The most prevalent pathogens are bacteria and fungi, but viruses should not be overlooked. For viral metagenomic analysis, we collected blood samples from eight patients with post-traumatic sepsis before and seven days after treatment. The results demonstrated that Anellovirus predominated the viral community, followed by Siphoviridae and Myoviridae, and that the variations in viral community and viral load before and after treatment were not statistically significant. This study allows us to investigate methods for establishing NGS-based viral diagnostic instruments for detecting viral infections in the blood of sepsis patients so that antiviral therapy can be administered quickly.

Virome Analysis Provides an Insight into the Viral Community of Chinese Mitten Crab Eriocheir sinensis

Recent advances in viromics have led to the discovery of a great diversity of RNA viruses and the identification of a large number of viral pathogens. A systematic exploration of viruses in Chinese mitten crab (Eriocheir sinensis), one of the most important aquatic commercial species, is still lacking. Here, we characterized the RNA viromes of asymptomatic, milky disease (MD)-affected, and hepatopancreatic necrosis syndrome (HPNS)-affected Chinese mitten crabs collected from 3 regions in China. In total, we identified 31 RNA viruses belonging to 11 orders, 22 of which were first reported here. By comparing viral composition between samples, we observed high variation in viral communities across regions, with most of the viral species being region-specific. We proposed to establish several novel viral families or genera based on the phylogenetic relationships and genome structures of viruses discovered in this study, expanding our knowledge of viral diversity in brachyuran crustaceans. IMPORTANCE High-throughput sequencing and meta-transcriptomic analysis provide us with an efficient tool to discover unknown viruses and explore the composition of viral communities in specific species. In this study, we investigated viromes in asymptomatic and diseased Chinese mitten crabs collected from three distant locations. We observed high regional variation in the composition of viral species, highlighting the importance of multi-location sampling. In addition, we classified several novel and ICTV-unclassified viruses based on their genome structures and phylogenetic relationships, providing a new perspective on current viral taxa.

Multiple genotypes infection and molecular characterization of Torque teno neovison virus: A novel Anelloviridae of mink in China

A novel Torque teno neovison virus (TTVs) was identified in specimens collected from dead mink during an outbreak of the Aleutian mink disease virus. Eighteen complete genomic sequences were obtained, ranging from 2109 to 2158 nucleotides in length and consisting of an untranslated region and three open reading frames. The genomic organization of mink TTVs is similar to previously reported anelloviruses. However, the deduced amino acid sequence of its ORF1 protein shows genetic diversity compared to related anelloviruses, suggesting that it represents a putative new species within the Anelloviridae family. This study provides a detailed molecular characterization of the novel mink anelloviruses, including its codon usage pattern, origin, and evolution. Analysis of the viral genomic sequences reveals the existence of multiple genotypes of co-infection. Principal component analysis and phylogenetic trees confirm the coexistence of multiple genotypes. Furthermore, the codon usage analyses indicate that mink TTVs have a genotype-specific codon usage pattern and show a low codon usage bias. Host-specific adaptation analysis suggests that TTVs are less adapted to mink. The possible origin and evolutionary history of mink TTVs were elucidated. Mink TTVs was genetically closely related to giant panda anellovirus, representing a new species. The observed incongruence between the phylogenetic history of TTVs and that of their hosts suggests that the evolution of anellovirus is largely determined by cross-species transmission. The study provides insights into the co-infection and genetic evolution of anellovirus in China.

Viromics Reveals the High Diversity of Viruses from Fishes of the Tibet Highland

Aquaculture is important for food security and nutrition. The economy has recently been significantly threatened and the risk of zoonoses significantly increased by aquatic diseases, and the ongoing introduction of new aquatic pathogens, particularly viruses, continues to represent a hazard. Yet, our knowledge of the diversity and abundance of fish viruses is still limited. Here, we conducted a metagenomic survey of different species of healthy fishes caught in the Lhasa River, Tibet, China, and sampled intestinal contents, gills, and tissues. To be more precise, by identifying and analyzing viral genomes, we aim to determine the abundance, diversity, and evolutionary relationships of viruses in fish with other potential hosts. Our analysis identified 28 potentially novel viruses, 22 of which may be associated with vertebrates, across seven viral families. During our research, we found several new strains of viruses in fish, including papillomavirus, hepadnavirus, and hepevirus. Additionally, we discovered two viral families, Circoviridae and Parvoviridae, which were prevalent and closely related to viruses that infect mammals. These findings further expand our understanding of highland fish viruses and highlight the emerging view that fish harbor large, unknown viruses. IMPORTANCE The economy and zoonoses have recently been significantly threatened by aquatic diseases. Yet, our knowledge of the diversity and abundance of fish viruses is still limited. We identified the wide genetic diversity of viruses that these fish were harboring. Since there are currently few studies on the virome of fish living in the Tibet highland, our research adds to the body of knowledge. This discovery lays the groundwork for future studies on the virome of fish species and other highland animals, preserving the ecological equilibrium on the plateau.

Evolution of anelloviruses from a circovirus-like ancestor through gradual augmentation of the jelly-roll capsid protein

Anelloviruses are highly prevalent in diverse mammals, including humans, but so far have not been linked to any disease and are considered to be part of the 'healthy virome'. These viruses have small circular single-stranded DNA (ssDNA) genomes and encode several proteins with no detectable sequence similarity to proteins of other known viruses. Thus, anelloviruses are the only family of eukaryotic ssDNA viruses currently not included in the realm Monodnaviria. To gain insights into the provenance of these enigmatic viruses, we sequenced more than 250 complete genomes of anelloviruses from nasal and vaginal swab samples of Weddell seal (Leptonychotes weddellii) from Antarctica and a fecal sample of grizzly bear (Ursus arctos horribilis) from the USA and performed a comprehensive family-wide analysis of the signature anellovirus protein ORF1. Using state-of-the-art remote sequence similarity detection approaches and structural modeling with AlphaFold2, we show that ORF1 orthologs from all Anelloviridae genera adopt a jelly-roll fold typical of viral capsid proteins (CPs), establishing an evolutionary link to other eukaryotic ssDNA viruses, specifically, circoviruses. However, unlike CPs of other ssDNA viruses, ORF1 encoded by anelloviruses from different genera display remarkable variation in size, due to insertions into the jelly-roll domain. In particular, the insertion between β-strands H and I forms a projection domain predicted to face away from the capsid surface and function at the interface of virus-host interactions. Consistent with this prediction and supported by recent experimental evidence, the outermost region of the projection domain is a mutational hotspot, where rapid evolution was likely precipitated by the host immune system. Collectively, our findings further expand the known diversity of anelloviruses and explain how anellovirus ORF1 proteins likely diverged from canonical jelly-roll CPs through gradual augmentation of the projection domain. We suggest assigning Anelloviridae to a new phylum, 'Commensaviricota', and including it into the kingdom Shotokuvirae (realm Monodnaviria), alongside Cressdnaviricota and Cossaviricota.

First Insights into the Occurrence of Circular Single-Stranded DNA Genomes in Asian and African Cattle

Circular replicase-encoding single-stranded (CRESS) DNA viruses and other circular DNA agents are increasingly found in various samples and animals. A specific class of these agents-termed bovine meat and milk factors (BMMF)-has been supposed to act as a factor in indirect carcinogenesis in humans. Initial observations attributed the BMMF to European cattle breeds and foodstuffs produced thereof. In the present study, blood and fecal samples from African and Asian cattle were examined. BMMF molecules and genomoviruses were detected in all bovids under study. The majority (79%) of the 29 circular elements could be assigned to BMMF groups 1 and 2, whereas CRESS viruses of the family Genomoviridae accounted for the smaller part (21%). Two genomoviruses belong to the genus Gemykibivirus and one to the genus Gemykrogvirus. The remaining three might be considered as novel species within the genus Gemycircularvirus. The majority of all isolated molecules originated from fecal samples, whereas only three derived from blood. The results from this study expand our knowledge on the diversity and presence of circular DNA in different ruminants that serve for food production in many countries over the world.

Characterization of a novel papillomavirus identified from a whale (Delphinapterus leucas) pharyngeal metagenomic library

Here, using viral metagenomic method, a novel whale papillomavirus (temporarily named wPV, GenBank accession number OP856597) was discovered in a whale (Delphinapterus leucas) pharyngeal metagenomic library. The complete genome size of wPV is 7179 bp, with GC content of 54.4% and a nucleotide composition of 23.4% A, 22.3% T, 28.4% G, and 25.9% C. The viral genome has a typical papillomavirus organization pattern, and five ORFs were predicted, including two late genes encoding L1 and L2, and three early genes encoding E1, E2, and E6. Pairwise sequence comparison and phylogenetic analysis based on the L1 gene sequence indicated that wPV may be a novel species within genus Dyodeltapapillomavirus. In addition, the E2 region of wPV was predicted to have a potential recombination event. The discovery of this novel papillomavirus increases our understanding of the viral ecology of marine mammals, providing insights into possible future infectious diseases.

Virome of wild rats (Rattus norvegicus) captured far from pig farms in Jiangsu province of China reveals novel porcine circovirus type 2d (PCV2d) sequences

BACKGROUND

Porcine circovirus type 2 (PCV2) has caused great economic losses in the global pig industry. There have been published records of wild rats acting as the reservoirs of PCV2 (only PCV2a and PCV2b), but almost all of which were related to the PCV2-infected swine herds.

RESULTS

In this study, we carried out the detection, amplification, and characterization of novel PCV2 strains in wild rats that were captured far from pig farms. Nested PCR assay demonstrated that the kidney, heart, lung, liver, pancreas, and large and small intestines of rats were screened positive for PCV2. We subsequently sequenced two full genomes of PCV2 in positive sample pools, designated as js2021-Rt001 and js2021-Rt002. Genome sequence analysis indicated that they had the highest similarity to nucleotide sequences of porcine-origin PCV2 isolates in Vietnam. Phylogenetically, js2021-Rt001 and js2021-Rt002 were a part of the PCV2d genotype cluster, which is a predominant genotype circulating worldwide in recent years. The antibody recognition regions, immunodominant decoy epitope, and heparin sulfate binding motif of the two complete genome sequences coincided with those previously reported.

CONCLUSIONS

Our research reported the genomic characterization of two novel PCV2 strains (js2021-Rt001 and js2021-Rt002) and provided the first supported evidence that PCV2d could naturally infect wild rats in China. However, whether the newly identified strains have potential for circulating in nature in vertical and horizontal transmission or inter-species jumping between rats and pigs needs further research.

Diverse DNA virus genomes identified in fecal samples of Mexican free-tailed bats (Tadarida brasiliensis) captured in Chiricahua Mountains of southeast Arizona (USA)

Bats (order Chiroptera) are some of the most abundant mammals on earth and their species ecology strongly influences zoonotic potential. While substantial research has been conducted on bat-associated viruses, particularly on those that can cause disease in humans and/or livestock, globally, limited research has focused on endemic bats in the USA. The southwest region of the US is of particular interest because of its high diversity of bat species. We identified 39 single-stranded DNA virus genomes in the feces of Mexican free-tailed bats (Tadarida brasiliensis) sampled in the Rucker Canyon (Chiricahua Mountains) of southeast Arizona (USA). Twenty-eight of these belong to the virus families Circoviridae (n = 6), Genomoviridae (n = 17), and Microviridae (n = 5). Eleven viruses cluster with other unclassified cressdnaviruses. Most of the viruses identified represent new species. Further research on identification of novel bat-associated cressdnaviruses and microviruses is needed to provide greater insights regarding their co-evolution and ecology relative to bats.

Characterization of Iflavirus in the Red Flour Beetle, Tribolium castaneum (Coleoptera; Tenebrionidae)

Iflavirus is a group of viruses distributed mainly in arthropod species. We surveyed Tribolium castaneum iflavirus (TcIV) in different laboratory strains and in Sequence Read Archives (SRA) in GenBank. TcIV is highly specific to only T. castaneum and is not found in seven other Tenebrionid species, including the closely related species T. freemani. The same strains from different laboratories and different strains displayed largely different degrees of infections in the examination of 50 different lines by using Taqman-based quantitative PCR. We found that ~63% (27 out of 43 strains) of T. castaneum strains in different laboratories are positive for TcIV PCR with large degrees of variation, in the range of seven orders of magnitude, indicating that the TcIV is highly fluctuating depending on the rearing conditions. The TcIV was prevalent in the nervous system with low levels found in the gonad and gut. The transovarial transmission was supported in the experiment with surface-sterilized eggs. Interestingly, TcIV infection did not show observable pathogenicity. TcIV offers an opportunity to study the interaction between the virus and the immune system of this model beetle species.

Virome diversity of ticks feeding on domestic mammals in China

Ticks are considered the second most common pathogen vectors transmitting a broad range of vital human and veterinary viruses. From 2017 to 2018, 640 ticks were collected in eight different provinces in central and western China. Six species were detected, including H.longicornis, De.everestianus, Rh.microplus, Rh.turanicus, Rh.sanguineous, and Hy.asiaticum. Sixty-four viral metagenomic libraries were constructed on the MiSeq Illumina platform, resulting in 13.44 ​G (5.88 ​× ​10⁷) of 250-bp-end reads, in which 2,437,941 are viral reads. We found 27 nearly complete genome sequences, including 16 genome sequences encoding entire protein-coding regions (lack of 3' or 5' end non-coding regions) and complete viral genomes, distributed in the arboviral family (Chuviridae, Rhabdoviridae, Nairoviridae, Phenuiviridae, Flaviviridae, Iflaviridae) as well as Parvoviridae and Polyomaviridae that cause disease in mammals and even humans. In addition, 13 virus sequences found in Chuviridae, Nairoviridae, Flaviviridae, Iflaviridae, Hepeviridae, Parvoviridae, and Polyomaviridae were identified as belonging to a new virus species in the identified viral genera. Besides, an epidemiological survey shows a high prevalence (9.38% and 15.63%) of two viruses (Ovine Copiparvovirus and Bovine parvovirus 2) in the tick cohort.

Genome of a husavirus from Southern Brazil

New viruses of the Picornavirales order have been discovered with the increase in the number of sequences obtained by high-throughput sequencing, as well as human stool-associated RNA virus (husavirus [HuV]), found in human stool samples. However, there is much to be clarified about HuV. Its cellular host, evolutionary history, and other biological characteristics are still unknown. Therefore, samples collected from human beings and environmental samples in a watershed in Southern Brazil were processed for the metagenomic library. Upon metagenomic analysis, we identified a HuV (husavirus LMM_67754 OP019707) genome with 8,846 bp, which was reported for the first time in Southern Brazil. The new genome presents only 37% of nucleotide identity with Brazilian strains and more than 90% with genomes from China, Vietnam, Venezuela, and the Netherlands. The HuV phylogeny presents significant differences among genomes, probably because multiple introductions of the virus may have occurred. Many questions still need to be answered about HuV. Therefore, more sequences and studies on this virus are necessary to improve the comprehension of the unknown origin of Picornavirales.

Viral metagenomics of the gut virome of diarrheal children with Rotavirus A infection

Diarrhea is a leading cause of morbidity and mortality in children worldwide and represents a major dysbiosis event. Rotavirus has been recognized as a global leading pathogen of diarrhea. This study is aimed at investigating differences in the gut virome between diarrheal children and healthy controls. In 2018, 76 diarrheal fecal samples and 27 healthy fecal samples in Shanghai and 40 diarrheal fecal samples and 19 healthy fecal samples in Taizhou were collected to investigate the composition of the gut virome. Viral metagenomic analyses revealed that the alpha diversity of the diarrheal virome was not significantly different from that of the healthy virome, and the beta diversity had a significant difference between diarrheal and healthy children. The diarrheal virome was mainly dominated by the families Adenoviridae, Astroviridae, Caliciviridae, and Picornaviridae. Meanwhile, the healthy virome also contains phages, including Microviridae and Caudovirales. The high prevalence of diverse enteric viruses in all samples and the little abundance of Microviridae and Caudovirales in diarrheal groups were identified. The study introduced a general overview of the gut virome in diarrheal children, revealed the compositional differences in the gut viral community compared to healthy controls, and provided a reference for efficient treatments and prevention of virus-infectious diarrhea in children.

Virome Analysis of Normal and Growth Retardation Disease-Affected Macrobrachium rosenbergii

The giant freshwater prawn, Macrobrachium rosenbergii, is an important aquaculture species in China. Growth retardation disease (GRD) is a common contagious disease in M. rosenbergii, resulting in slow growth and precocious puberty in prawns, and has caused growing economic losses in the M. rosenbergii industry. To investigate the viral diversity of M. rosenbergii and identify potentially high-risk viruses linked to GRD, virome analysis of the GRD-affected and normal M. rosenbergii was carried out using next-generation sequencing (NGS). A total of 327 contigs (>500 bp) were related to viral sequences belonging to 23 families/orders and a group of unclassified viruses. The majority of the viral contigs in M. rosenbergii belonged to the order Picornavirales, with the Solinviviridae family being the most abundant in both the diseased and normal groups. Furthermore, 16 RNA viral sequences with nearly complete genomes were characterized and phylogenetically analyzed, belonging to the families Solinviviridae, Flaviviridae, Polycipiviridae, Marnaviridae, and Dicistroviridae as well as three new clades of the order Picornavirales. Notably, the cross-species transmission of a picorna-like virus was observed between M. rosenbergii and plants. The "core virome" seemed to be present in the diseased and normal prawns. Still, a clear difference in viral abundance was observed between the two groups. These results showed that the broad diversity of viruses is present in M. rosenbergii and that the association between viruses and disease of M. rosenbergii needs to be further investigated. IMPORTANCE Growth retardation disease (GRD) has seriously affected the development and economic growth of the M. rosenbergii aquaculture industry. Our virome analysis showed that diverse viral sequences were present in M. rosenbergii, significantly expanding our knowledge of viral diversity in M. rosenbergii. Some differences in viral composition were noted between the diseased and normal prawns, indicating that some viruses become more abundant in occurrences or outbreaks of diseases. In the future, more research will be needed to determine which viruses pose a risk for M. rosenbergii. Our study provides important baseline information contributing to disease surveillance and risk assessment in M. rosenbergii aquaculture.

Novel Virus Identification through Metagenomics: A Systematic Review

Metagenomic Next Generation Sequencing (mNGS) allows the evaluation of complex microbial communities, avoiding isolation and cultivation of each microbial species, and does not require prior knowledge of the microbial sequences present in the sample. Applications of mNGS include virome characterization, new virus discovery and full-length viral genome reconstruction, either from virus preparations enriched in culture or directly from clinical and environmental specimens. Here, we systematically reviewed studies that describe novel virus identification through mNGS from samples of different origin (plant, animal and environment). Without imposing time limits to the search, 379 publications were identified that met the search parameters. Sample types, geographical origin, enrichment and nucleic acid extraction methods, sequencing platforms, bioinformatic analytical steps and identified viral families were described. The review highlights mNGS as a feasible method for novel virus discovery from samples of different origins, describes which kind of heterogeneous experimental and analytical protocols are currently used and provides useful information such as the different commercial kits used for the purification of nucleic acids and bioinformatics analytical pipelines.

Viral Metagenomics Reveals a Putative Novel HPV Type in Anogenital Wart Tissues

Viral metagenomics is widely employed to identify novel viruses in biological samples. Recently, although numerous novel human papillomavirus (HPV) types have been identified in clinical samples including anogenital warts (AGWs), many novel HPV sequences remain to be discovered. In this study, a putative novel HPV type designated as HPV-JDFY01 was discovered from library GW05 with 63 sequence reads by the viral metagenomic technique. Its complete genomic sequence was determined by PCR to bridge the gaps between contigs combining Sanger sequencing. The complete genome of HPV-JDFY01 is a 7186 bp encoding 7 open reading frames (ORFs) (E6, E7, E1, E2, E4, L2 and L1) and contains a 487 bp long control region (LCR) between L1 and E6. Sequence and phylogeny analysis indicated that HPV-JDFY01 shared the highest sequence identity of 74.2% with HPV-mSK_244 (MH777383) and well clustered into the genus Gammapapillomavirus. It has the classical genomic organization of Gammapapillomaviruses. Epidemiological investigation showed that one out of the 413 AGW tissue samples was positive for HPV-JDFY01. Further research with large size and different type of samples should be performed to elucidate the epidemiologic status of HPV-JDFY01.

Expanding known viral diversity in plants: virome of 161 species alongside an ancient canal

BACKGROUND

Since viral metagenomic approach was applied to discover plant viruses for the first time in 2006, many plant viruses had been identified from cultivated and non-cultivated plants. These previous researches exposed that the viral communities (virome) of plants have still largely uncharacterized. Here, we investigated the virome in 161 species belonging to 38 plant orders found in a riverside ecosystem.

RESULTS

We identified 245 distinct plant-associated virus genomes (88 DNA and 157 RNA viruses) belonging to 27 known viral families, orders, or unclassified virus groups. Some viral genomes were sufficiently divergent to comprise new species, genera, families, or even orders. Some groups of viruses were detected that currently are only known to infect organisms other than plants. It indicates a wider host range for members of these clades than previously recognized theoretically. We cannot rule out that some viruses could be from plant contaminating organisms, although some methods were taken to get rid of them as much as possible. The same viral species could be found in different plants and co-infections were common.

CONCLUSIONS

Our data describe a complex viral community within a single plant ecosystem and expand our understanding of plant-associated viral diversity and their possible host ranges.

Complete Mitogenomes of Ticks Ixodes acutitarsus and Ixodes ovatus Parasitizing Giant Panda: Deep Insights into the Comparative Mitogenomic and Phylogenetic Relationship of Ixodidae Species

Ticks rank second in the world as vectors of disease. Tick infestation is one of the factors threatening the health and survival of giant pandas. Here, we describe the mitogenomes of Ixodes acutitarsus and Ixodes ovatus parasitizing giant pandas, and perform comparative and phylogenetic genomic analyses on the newly sequenced and other available mitogenomes of hard ticks. All six newly determined mitogenomes contain a typical gene component and share an ancient Arthropoda gene arrangement pattern. Our study suggests that I. ovatus is a species complex with high genetic divergence, indicating that different clades of I. ovatus represent distinct species. Comparative mitogenomic analyses show that the average A + T content of Ixodidae mitogenomes is 78.08%, their GC-skews are strongly negative, while AT-skews fluctuate around 0. A large number of microsatellites are detected in Ixodidae mitogenomes, and the main microsatellite motifs are mononucleotide A and trinucleotide AAT. We summarize five gene arrangement types, and identify the trnY-COX1-trnS1-COX2-trnK-ATP8-ATP6-COX3-trnG fragment is the most conserved region, whereas the region near the control region is the rearrangement hotspot in Ixodidae mitogenomes. The phylogenetic trees based on 15 genes provide a very convincing relationship (Ixodes + (Robertsicus + ((Bothriocroton + Haemaphysalis) + (Amblyomma + (Dermacentor + (Rhipicentor + (Hyalomma + Rhipicephalus))))))) with very strong supports. Remarkably, Archaeocroton sphenodonti is embedded in the Haemaphysalis clade with strong supports, resulting in paraphyly of the Haemaphysalis genus, so in-depth morphological and molecular studies are essential to determine the taxonomic status of A. sphenodonti and its closely related species. Our results provide new insights into the molecular phylogeny and evolution of hard ticks, as well as basic data for population genetics assessment and efficient surveillance and control for the giant panda-infesting ticks.

Viral metagenomics of pharyngeal secretions from children with acute respiratory diseases with unknown etiology revealed diverse viruses

Acute respiratory tract infections are a major public health problem and the leading cause of morbidity in children younger than 5 years old. This study investigated the potential reasons of unexplained acute respiratory infections in children in Xuzhou and its environs during 2018-2019.We collected pharyngeal swab samples from 411 children under the age of five who presented with symptoms of unexplained acute respiratory infection and were negative for bacteria, mycoplasma, and influenza viruses. Using viral metagenomic techniques, viral nucleic acids were extracted, enriched, and sequenced from the samples. Results indicated that Picornaviridae, Parvoviridae, Paramyxoviridae, Coronaviridae, and Anelloviridae were the five virus families with the highest relative content of sequence reads. And we detected 35 HBoV-positive and 12 HEV-positive samples out of 411 samples by the polymerase chain reaction (PCR). Partial or nearly complete genome sequences of viruses belonging to the families Picornaviridae, Parvoviridae, and Anelloviridae were characterized, and phylogenetic trees were constructed based on the nucleic acid or amino acid sequences of the predicted viral open reading frames (ORFs), as well as genotyping of the viruses. In addition, we observed recombination events in the Saffold virus and Coxsackievirus A9 by analyzing the genetic characteristics of the viruses revealed in this study. This study provides vital information for the prevention and treatment of acute respiratory infections in children younger than five years old.

Viromic analysis of feces from laboratory rabbits reveals a new Circovirus

BACKGROUND

Members of the genus Circovirus with the family Circoviridae are responsible for fatal diseases that can affect mammals and birds. Beak and feather disease virus (BFDV) is responsible for fatal diseases that could affect birds, causing the psittacine beak and feather disease. The current study discovered a new Circovirus from feces of laboratory rabbits and name it RabCV, which shows close relationship to BFDVs.

RESULTS

We investigated the feces virome of 10 laboratory rabbits using the viral metagenomic method. In these samples, we detected a new rabbit-associated Circovirus (RabCV) and performed phylogenetic analysis based on replication-associated (Rep) protein. The result showed that the RabCV was closely clustered with BFDVs, sharing the identity of 56.7%-57.2% with them based on the whole genome sequence. PCR screening in a cohort of 38 laboratory rabbits showed that 3 out of the 38 rabbits were positive for this new rabbit-associated Circovirus.

CONCLUSION

A new Circovirus was discovered from feces of rabbits, which showed low prevalence in the healthy laboratory rabbits. BFDV is responsible for fatal diseases that could affect birds, which suggested that the potential threat of the new rabbit-associated Circovirus to the health of laboratory rabbits needs further study.

Virome of Giant Panda-Infesting Ticks Reveals Novel Bunyaviruses and Other Viruses That Are Genetically Close to Those from Giant Pandas

Tick infestations have been reported as one of the factors threatening the health of giant pandas, but studies of viral pathogens carried by ticks feeding on the blood of giant pandas are limited. To assess whether blood-sucking ticks of giant pandas can carry viral pathogens and if so, whether the viruses in ticks are associated with those previously detected in giant panda hosts, we determined the viromes of ticks detached from giant pandas in a field stocking area in Sichuan Province, southwest China. Using viral metagenomics we identified 32 viral species in ticks, half of which (including anellovirus [n = 9], circovirus [n = 3], and gemycircularvirus [n = 4]) showed homology to viruses carried by giant pandas and their associated host species (such as red pandas and mosquitoes) in the same living domain. Remarkably, several viruses in this study phylogenetically assigned as bunyavirus, hepe-like virus, and circovirus were detected with relatively high abundance, but whether these newly identified tick-associated viruses can replicate in ticks and then transmit to host animals during a blood meal will require further investigation. These findings further expand our understanding of the role of giant panda-infesting ticks in the local ecosystem, especially related to viral acquisition and transmission, and lay a foundation to assess the risk for giant panda exposure to tick-borne viruses. IMPORTANCE Ticks rank only second to mosquitoes as blood-feeding arthropods, capable of spreading pathogens (including viruses, bacteria, and parasites) to hosts during a blood meal. To better understand the relationship between viruses carried by ticks and viruses that have been reported in giant pandas, it is necessary to analyze the viromes of giant panda-parasitic blood-sucking ticks. This study collected 421 ticks on the body surface of giant pandas in Sichuan Province, China. We characterized the extensive genetic diversity of viruses harbored by these ticks and reported frequent communication of viruses between giant pandas and their ticks. While most of the virome discovered here are nonpathogenic viruses from giant pandas and potentially tick-specific viruses, we revealed some possible tick-borne viruses, represented by novel bunyaviruses. This research contributes to the literature because currently there are few studies on the virome of giant panda-infesting ticks.

Dynamic alterations of the mice gut virome after Coxsackievirus B3 infection

The gut microbiome plays an essential role in the human health and dysbiosis has been implicated in numerous diseases. Coxsackievirus B3 infects millions of humans yearly and yet limited research has explored dynamic alterations of the gut virome after infection. Here, we established the mouse model of Coxsackievirus B3 infection and collected fecal samples at several time points to investigate alterations of the gut virome using viral metagenomic analysis. We found that the mice virome was dominated by Caudovirales and Microviridae, and phylogenetic analyses showed that both Caudovirales and Microviridae had high diversity. The gut virome had significant variations with the increase of Caudovirales and the decrease of Microviridae after infection. We proposed that Caudovirales and Microviridae may be biomarkers for the Coxsackievirus infection process. This study provides a reference for the dynamic changes of the gut virome after human Enterovirus infection, which may help guide the rational drug use in clinical treatment and provide new ideas for preventing Enterovirus infection.

Genome Analysis of Goose-Origin Astroviruses Causing Fatal Gout in Shanghai, China Reveals One of Them Belonging to a Novel Type Is a Recombinant Strain

Since 2014, a goose-origin astroviruses disease, which is characterized by urate precipitation in viscera, has rapidly spread to major commercial goose provinces leading to huge economic losses in the poultry industry of China. In March 2020, a goose farm locates in Shanghai, China, where there was no goose astroviruses (GAstVs) infection reported before, experienced an outbreak of gout disease in geese. The etiological investigation was carried out by virus metagenomics and bacterial culture and two GAstVs strains, designated as CHSH01 and CHSH02, were determined. Their complete genomes were measured to 7,154 and 7,330 nt in length, excludingthe poly(A) tail, respectively, and had different genomic features and classifications. CHSH01 shared a very low sequence identity with other strains in terms of not only the complete genome but also different ORFs. Phylogenetic analysis showed CHSH02 belonged to GAstV-2, which was the predominant species in the geese with gout in China according to the previous study. Meanwhile, CHSH01 strain displayed low identity with other AstVs, and phylogenetic and recombination analysis suggested that CHSH01 belonging to a novel type was a recombinant strain, one parent strain of which was an AstV determined from a bar-headed goose (a kind of migrant bird). Moreover, the primary epidemiological investigation showed that the two strains were prevalent in the same goose farm and co-infection occurred. These findings arise the potential cross-species transmission of CHSH01 between domestic and wild fowl.

Detection of Equus caballus papillomavirus-2 in equine penile/preputial papillomas and squamous cell carcinomas in southern Brazil

For approximately one decade, a novel papillomavirus termed Equus caballus papillomavirus-2 (EcPV-2) has been associated with equine penile/preputial papillomas and squamous cell carcinomas (SCCs). It is currently believed that the virus has a carcinogenic activity, being able to induce such neoplastic lesions. After being first described, EcPV-2 has been detected in many countries from North America, Europe, and Asia; however, to date, it has not been reported in Brazil. The aim of this research was to investigate the presence of EcPV-2 in penile/preputial papillomas and SCCs of Brazilian horses. Forty samples diagnosed as equine penile and/or preputial papillomas, carcinomas in situ (CIS), or SCCs in two veterinary anatomic pathology services from southern Brazil were investigated. Histologic evaluation and immunohistochemistry (IHC) using a BPV-1 antibody were performed. Posteriorly, the samples were submitted to polymerase chain reaction using two broad-spectrum (MY09/11 and FAP) and one EcPV-2-specific primer sets. Positive samples were sequenced. PV antigen expression was detected in one papilloma, one CIS, and two SCCs by IHC. Five SCCs, one papilloma, and one CIS were PV-positive on PCR. Sequencing of the seven PCR products revealed homology with EcPV-2. This study confirms the occurrence of EcPV-2 infection in Brazilian horses. Moreover, the results presented here provide useful information concerning the phylogeny from the viruses detected in our samples. We hope to encourage further studies on this novel agent, contributing to its characterization, and, possibly, to the eventual development of preventive measurements, including a possible vaccine.

Viral metagenomics unveiled extensive communications of viruses within giant pandas and their associated organisms in the same ecosystem

Cross-species transmission events were commonplace, with numerous cases of host-switching during the viral evolutionary history, but relatively little evidence for onward transmission in different species living in the same ecosystem. For understanding the communications of viruses in giant pandas (Ailuropoda melanoleuca) and their associated organisms, based on a large size of samples (N = 2305) collected between 2015 and 2020 from giant panda (N = 776) and other four giant panda-associated organisms in the same ecosystem, red pandas (N = 700), stray cats (N = 32), wild rats (N = 42), and mosquitoes (N = 755), viromics was used for the virus identification and subsequent virus traceability. The results showed that a feline panleukopenia virus (FPV) was found in giant pandas with clinical signs of vomiting and mild diarrhea. Meanwhile, the same FPV strain was also prevalent in the healthy red panda (Ailurus fulgens) population. From the viromes of the five different organisms, 250 virus genomes were determined. Our data revealed that besides FPV, other putative pathogenic viruses, such as red panda amdoparvoviruses (RPAVs) and Getah viruses (GETVs) were responsible for previous disease or death of some red pandas. We also demonstrated that a number of viruses were involved in potential interspecies jumping events between giant pandas and their associated species. Collectively, our results shed light on the genetic diversity and relationship of diverse viral pathogens in 'Giant pandas-Associated animals-Arthropods' and report some cases of possible viral host-switching among these host species living in the same ecosystem.

Virome in the cloaca of wild and breeding birds revealed a diversity of significant viruses

BACKGROUND

Wild birds may harbor and transmit viruses that are potentially pathogenic to humans, domestic animals, and other wildlife.

RESULTS

Using the viral metagenomic approach, we investigated the virome of cloacal swab specimens collected from 3182 birds (the majority of them wild species) consisting of > 87 different species in 10 different orders within the Aves classes. The virus diversity in wild birds was higher than that in breeding birds. We acquired 707 viral genomes from 18 defined families and 4 unclassified virus groups, with 265 virus genomes sharing < 60% protein sequence identities with their best matches in GenBank comprising new virus families, genera, or species. RNA viruses containing the conserved RdRp domain with no phylogenetic affinity to currently defined virus families existed in different bird species. Genomes of the astrovirus, picornavirus, coronavirus, calicivirus, parvovirus, circovirus, retrovirus, and adenovirus families which include known avian pathogens were fully characterized. Putative cross-species transmissions were observed with viruses in wild birds showing > 95% amino acid sequence identity to previously reported viruses in domestic poultry. Genomic recombination was observed for some genomes showing discordant phylogenies based on structural and non-structural regions. Mapping the next-generation sequencing (NGS) data respectively against the 707 genomes revealed that these viruses showed distribution pattern differences among birds with different habitats (breeding or wild), orders, and sampling sites but no significant differences between birds with different behavioral features (migratory and resident).

CONCLUSIONS

The existence of a highly diverse virome highlights the challenges in elucidating the evolution, etiology, and ecology of viruses in wild birds. Video Abstract.

The virome of the white-winged vampire bat Diaemus youngi is rich in circular DNA viruses

In the Neotropical region, the white-winged vampire bat (Diaemus youngi) is the rarest of the three species of vampire bats. This bat species feeds preferentially on bird blood, and there is limited information on the viruses infecting D. youngi. Hence, this study aimed to expand the knowledge about the viral diversity associated with D. youngi by sampling and pooling the lungs, liver, kidneys, heart, and intestines of all animals using high-throughput sequencing (HTS) on the Illumina MiSeq platform. A total of three complete and 10 nearly complete circular virus genomes were closely related to gemykrogvirus (Genomoviridae family), smacovirus (Smacoviridae family), and torque teno viruses (TTVs) (Anelloviridae family). In addition, three sequences of bat paramyxovirus were detected and found to be closely related to viruses reported in Pomona roundleaf bats and rodents. The present study provides a snapshot of the viral diversity associated with white-winged vampire bats and provides a baseline for comparison to viruses detected in future outbreaks.

Gut Virome of the World's Highest-Elevation Lizard Species (Phrynocephalus erythrurus and Phrynocephalus theobaldi) Reveals Versatile Commensal Viruses

The gut virome is a reservoir of diverse symbiotic and pathogenic viruses coevolving with their hosts, and yet limited research has explored the gut viromes of highland-dwelling rare species. Using viral metagenomic analysis, the viral communities of the Phrynocephalus lizards living in the Qinghai-Tibet Plateau were investigated. Phage-encoded functional genes and antibiotic resistance genes (ARGs) were analyzed. The viral communities of different lizard species were all predominated by bacteriophages, especially the Caudovirales order. The virome of Phrynocephalus erythrurus living around the Namtso Lake possessed a unique structure, with the greatest abundance of the Parvoviridae family and the highest number of exclusive viral species. Several vertebrate-infecting viruses were discovered, including caliciviruses, astroviruses, and parvoviruses. Phylogenetic analyses demonstrated that the virus hallmark genes of bacteriophages possessed high genetic diversity. After functional annotation, the majority of phage-associated functional genes were classified in the energy metabolism category. In addition, plenty of ARGs belonging to the multidrug category were discovered, and five ARGs were exclusive to the virome from Phrynocephalus theobaldi. This study provided the first insight into the structure and function of the virome in highland lizards, contributing to the protection of threatened lizard species. Also, our research is of exemplary significance for the gut virome research of lizard species and other cold-blooded and highland-dwelling animals, prompting a better understanding of the interspecific differences and transmission of commensal viruses. IMPORTANCE The Phrynocephalus lizards inhabiting the Qinghai-Tibet Plateau (QTP) are considered to be the highest-altitude lizard species in the world, and they have been added to the IUCN list of threatened species. Living in the QTP with hypoxic, arid, natural conditions, the lizards presented a unique pattern of gut virome, which could provide both positive and negative effects, such as the enrichment of functional genes and the dissemination of antibiotic resistance genes (ARGs). This work provides the foundation for further research on the gut virome in these endangered lizard species and other cold-blooded and highland-dwelling animals, contributing to the maintenance of ecological balance on the plateau.

Within-host quantitation of anellovirus genome complexity from clinical samples

Anellovirus (AV) is a ubiquitous and diverse virus in the human population. An individual can be infected with multiple AV genera and species that form a heterogeneous repertoire, called the anellome. Due to its exceptional genetic diversity, efficient evaluation of anellome complexity remains a methodological challenge. In the current study, AV genome was first enriched from patient serum samples through two-phase rolling circle amplification. Following Illumina sequencing, anellome was analyzed with an advanced bioinformatics pipeline, including read extraction at three similarity levels, de novo assembly, species assignment, and determination of relative abundance among AV variants. The method was validated in the mock sample and then applied to 21 hepatitis C virus (HCV) patients with and without hepatocellular carcinoma (HCC). Overall, there was a large variance regarding AV richness, ranging from 2 to 51 AV species. In contrast to HCV patients without HCC, HCC incidence was associated with reduced richness (12.6 ± 14.4 vs. 35.4 ± 13.6, p = 0.001) and Shannon entropy (0.4 ± 0.34 vs. 0.61 ± 0.12, p = 0.095) at the AV species level. Interestingly, AV genus beta and gamma expanded in the anellome in 7 of 10 HCC patients. These observations shed light on the potential association between anellome and HCC incidence in patients with chronic HCV infection. The method presented here represents a valuable tool to investigate the role of anellome in human health and disease.

Viral metagenomics reveals diverse viruses in the fecal samples of children with diarrhea

Diarrhea is the third leading cause of death in developing countries in children under the age of five. About half a million children die of diarrhea every year, most of which in developing countries. Viruses are the main pathogen of diarrhea. In China, the fecal virome of children with diarrhea has been rarely studied. Using an unbiased viral metagenomics approach, we analyzed the fecal virome in children with diarrhea. Many DNA or RNA viruses associated with diarrhea identified in those fecal samples were mainly from six families of Adenoviridae, Astroviridae, Caliciviridae, Parvoviridae, Picornaviridae, and Reoviridae. Among them, the family of Caliciviridae accounts for the largest proportion of 78.42%, following with Adenoviridae (8.94%) and Picornaviridae (8.36%). In addition to those diarrhea-related viruses that have already been confirmed to cause human diarrhea, the viruses not associated with diarrhea were also identified including anellovirus and picobirnavirus. This study increased our understanding of diarrheic children fecal virome and provided valuable information for the prevention and treatment of viral diarrhea in this area.

•

Many DNA or RNA viruses associated with diarrhea were identified in this study.

•

Viruses belonging to the family of Caliciviridae were the most main pathogen that induced children diarrhea.

•

In addition to those diarrhea-related viruses, the viruses not associated with diarrhea were also identified.

The virome of bronchoalveolar lavage fluid from patients with fever of unknown origin

Background: Viral metagenomics, a high-throughput sequencing combined with virus sequence-independent amplification by random PCR, allows for unbiased detection of virtually any viruses present in samples. Materials & methods: In order to investigate the virome of bronchoalveolar lavage fluid from patients with fever of unknown origin, 58 samples collected from diseased patients were characterized and compared. Results: Some representatives of Anelloviridae were identified, we found the torque teno virus (TTV) accounts for the majority of virus communities and were more prevalent in the specimens of febrile patients. Phylogenetic analysis suggested that these anellovirus isolates were close to the previous TTV available in GenBank®. Conclusion: All these data indicate that the human anellovirus species TTV may associated with fever of unknown origin.

Metagenomic Analysis of Viral Community in the Yangtze River Expands Known Eukaryotic and Prokaryotic Virus Diversity in Freshwater

Viruses in aquatic ecosystems are characterized by extraordinary abundance and diversity. Thus far, there have been limited studies focused on viral communities in river water systems. Here, we investigated the virome of the Yangtze River Delta using viral metagenomic analysis. The compositions of viral communities from six sampling sites were analyzed and compared. By using library construction and next generation sequencing, contigs and singlet reads similar to viral sequences were classified into 17 viral families, including nine dsDNA viral families, four ssDNA viral families and four RNA viral families. Statistical analysis using Friedman test suggested that there was no significant difference among the six sampling sites (P ​> ​0.05). The viromes in this study were all dominated by the order Caudovirales, and a group of Freshwater phage uvFW species were particularly prevalent among all the samples. The virome from Nanjing presented a unique pattern of viral community composition with a relatively high abundance of family Parvoviridae. Phylogenetic analyses based on virus hallmark genes showed that the Caudovirales order and CRESS-DNA viruses presented high genetic diversity, while viruses in the Microviridae and Parvoviridae families and the Riboviria realm were relatively conservative. Our study provides the first insight into viral community composition in large river ecosystem, revealing the diversity and stability of river water virome, contributing to the proper utilization of freshwater resource.

•

First insight into viral community composition in large river ecosystem.

•

Virus hallmark genes present both diverse and conservative characteristics.

•

The composition of viral communities is similar on the whole.

•

Slight regional variation of virome is existed in individual areas.

A variety of highly divergent eukaryotic ssDNA viruses in sera of pigs

Over the last decade, viral metagenomics has been established as a non-targeted approach for identifying viruses in stock animals, including pigs. This has led to the identification of a vast diversity of small circular ssDNA viruses. The present study focuses on the investigation of eukaryotic circular Rep-encoding single-stranded (CRESS) DNA viral genomes present in serum of commercially reared pigs from southern Brazil. Several CRESS DNA viral genomes were detected, including representatives of the families Smacoviridae (n=5), Genomoviridae (n=3), Redondoviridae (n=1), Nenyaviridae (n=1) and other yet unclassified genomes (n=9), plus a circular DNA molecule, which probably belongs to the phylum Cressdnaviricota. A novel genus within the family Smacoviridae, tentatively named 'Suismacovirus', comprising 21 potential new species, is proposed. Although the reported genomes were recovered from pigs with clinical signs of respiratory disease, further studies should examine their potential role as pathogens. Nonetheless, these findings highlight the diversity of circular ssDNA viruses in serum of domestic pigs, expand the knowledge on CRESS DNA viruses' genetic diversity and distribution and contribute to the global picture of the virome of commercially reared pigs.

Viral metagenomics reveals two novel anelloviruses in feces of experimental rats

BACKGROUND

Rodents are widely distributed and are the natural reservoirs of a diverse group of zoonotic viruses. Thus, analyzing the viral diversity harbored by rodents could assist efforts to predict and reduce the risk of future emergence of zoonotic viral diseases. Rodents are commonly used in animal testing, particularly mice and rats. Experimental rats are important animal models, and a history of pathogenic infections in these animals will directly affect the animal trial results. The pathogenicity of Anellovirus (AV) remains poorly understood due to the lack of a suitable model cell line or animal to support the viral cycle. This study aimed to discover possible anelloviruses from the virome in feces of experimental rats by viral metagenomic technique.

METHODS

Fecal samples were collected from 10 commercial SD rats and pooled into a sample pool and then subjected to libraries construction which was then sequenced on Illumina MiSeq platform. The sequenced reads were analyzed using viral metagenomic analysis pipeline and two novel anelloviruses (AVs) were identified from fecal sample of experimental rats. The prevalence of these two viruses was investigated by conventional PCR.

RESULTS

The complete genomic sequence of these two AVs were determined and fully characterized, with strain name ratane153-zj1 and ratane153-zj2. The circular genomes of ratane153-zj1 and ratane153-zj2 are 2785 nt and 1930 nt in length, respectively, and both include three ORFs. Ratane153-zj1 closely clustered with members within the genus Wawtorquevirus and formed a separate branch based on the phylogenetic tree constructed over the amino acid sequence of ORF1 of the two AVs identified in this study and other related AVs. While the complete amino acid sequences of ORF1 of ratane153-zj2 (nt 335 to 1390) had the highest sequence identity with an unclassified AV (GenBank No. ATY37438) from Chinchilla lanigera, and they clustered with one AV (GenBank No. QYD02305) belonging to the genus Etatorquevirus from Lynx rufus. Conventional PCR with two sets of specific primers designed based on the two genomes, respectively, showed that they were detectable at a low frequency in cohorts of experimental rats.

CONCLUSION

Our study expanded the genome diversity of AVs and provided genetic background information of viruses existed in experimental rats.

Virome in Fecal Samples From Wild Giant Pandas (Ailuropoda Melanoleuca)

The giant panda (Ailuropoda melanoleuca) is one of the most endangered mammals in the world; anthropogenic habitat loss and poaching still threaten the survival of wild pandas. Viral infection has become one of the potential threats to the health of these animals, but the available information related to these infections is still limited. In order to detect possible vertebrate viruses, the virome in the fecal samples of seven wild giant pandas from Qinling Mountains was investigated by using the method of viral metagenomics. From the fecal virome of wild giant pandas, we determined six nearly complete genomes belonging to the order Picornavirales, two of which may be qualified as a novel virus family or genus. In addition, four complete genomes belonging to the Genomoviridae family were also fully characterized. This virological investigation has increased our understanding of the gut viral community in giant pandas. Whether these viruses detected in fecal samples can really infect giant panda needs further research.

Identification of a novel circovirus in blood sample of giant pandas (Ailuropoda melanoleuca)

The members of the family Circoviridae are considered to be one of the smallest autonomously replicating viruses that are classified into two genera, Circovirus and Cyclovirus. Circoviruses have been found in a variety of vertebrates, but whether they infect endangered protected animals has not been studied in much detail. Here, viral metagenomics and PCR methods were used to detect and verify viral nucleic acid in the blood sample from giant pandas. According to these methods, the complete genome sequence of a novel circovirus, the giant panda associated circovirus (GPCV) from the blood sample of three giant pandas was identified. The GPCV genome is 2090 bp in size and reveals two putative ambisense open-reading frames, encoding the major structural capsid protein and the replication associated protein, respectively, the latter having two predicted introns. Pairwise sequence comparison and phylogenetic analyses indicated GPCV was a putative new species within genus Circovirus based on the species demarcation criteria of the International Committee on the Taxonomy of Viruses. It is the first time that circovirus has been identified from blood sample of giant pandas. These efforts will contribute to future analyses to illuminate the evolutionary relationships between classified and newly identified members of the family Circoviridae.

Taxonomic update for mammalian anelloviruses (family Anelloviridae)

Anelloviruses are small negative-sense single-stranded DNA viruses with genomes ranging in size from 1.6 to 3.9 kb. The family Anelloviridae comprised 14 genera before the present changes. However, in the last five years, a large number of diverse anelloviruses have been identified in various organisms. Here, we undertake a global analysis of mammalian anelloviruses whose full genome sequences have been determined and have an intact open reading frame 1 (ORF1). We established new criteria for the classification of anelloviruses, and, based on our analyses, we establish new genera and species to accommodate the unclassified anelloviruses. We also note that based on the updated species demarcation criteria, some previously assigned species (n = 10) merge with other species. Given the rate at which virus sequence data are accumulating, and with the identification of diverse anelloviruses, we acknowledge that the taxonomy will have to be dynamic and continuously evolve to accommodate new members.

Virome analysis of ticks in a forest region of Liaoning, China: characterization of a novel hepe-like virus sequence

Background

Ticks (class Arachnida, subclass Acari) are vectors of transmitting a broad range of pathogenic microorganisms, protozoa, and viruses affecting humans and animals. Liaoning Province is rich in forests where different animals and, abundant Haemaphysalis longicornis ticks exist.

Methods

Using viral metagenomics, we analyzed the virome in 300 Haemaphysalis longicornis ticks collected from June to August 2015 in the forested region of Liaoning Province, China.

Results

From the 300 ticks, 1,218,388 high-quality reads were generated, of which 5643 (0.463%) reads showed significant sequence identity to known viruses. Sequence and phylogenetic analysis revealed that viral sequences showing a close relationship with Dabieshan tick virus, Aleutian mink disease virus, adeno-associated virus, Gokushovirus, avian gyrovirus 2 were present in the virome of these ticks. However, the significance of these viruses to human and animal health requires further investigation. Notably, an hepe-like virus, named tick-borne hepe-like virus sequence, was obtained and was highly prevalent in these ticks with a rate of 50%. Nevertheless, one constraint of our study was the limited geographical distribution of the sampled ticks.

Conclusion

Our study offers an overview of the virome in ticks from a forest region of Liaoning Province and provides further awareness of the viral diversity of ticks.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01632-x.

Family Genomoviridae: 2021 taxonomy update

The family Genomoviridae (phylum Cressdnaviricota, class Repensiviricetes, order Geplafuvirales) includes viruses with circular single-stranded DNA genomes encoding two proteins, the capsid protein and the rolling-circle replication initiation protein. The genomes of the vast majority of members in this family have been sequenced directly from diverse environmental or animal- and plant-associated samples, but two genomoviruses have been identified infecting fungi. Since the last taxonomic update of the Genomoviridae, a number of new members of this family have been sequenced. Here, we report on the most recent taxonomic update, including the creation of one new genus, Gemytripvirus, and classification of ~420 new genomoviruses into 164 new species. We also announce the adoption of the "Genus + freeform epithet" binomial system for the naming of all 236 officially recognized species in the family Genomoviridae. The updated taxonomy presented in this article has been accepted by the International Committee on Taxonomy of Viruses (ICTV).

Complex evolutionary history of felid anelloviruses

Anellovirus infections are highly prevalent in mammals, however, prior to this study only a handful of anellovirus genomes had been identified in members of the Felidae family. Here we characterise anelloviruses in pumas (Puma concolor), bobcats (Lynx rufus), Canada lynx (Lynx canadensis), caracals (Caracal caracal) and domestic cats (Felis catus). The complete anellovirus genomes (n = 220) recovered from 149 individuals were diverse. ORF1 protein sequence similarity network analysis coupled with phylogenetic analysis, revealed two distinct clusters that are populated by felid-derived anellovirus sequences, a pattern mirroring that observed for the porcine anelloviruses. Of the two-felid dominant anellovirus groups, one includes sequences from bobcats, pumas, domestic cats and an ocelot, and the other includes sequences from caracals, Canada lynx, domestic cats and pumas. Coinfections of diverse anelloviruses appear to be common among the felids. Evidence of recombination, both within and between felid-specific anellovirus groups, supports a long coevolution history between host and virus.

Monsavirus in monkey rectal swab and throat swab specimens in China: Proposal for Posaliviridae as a new family in Picornavirales

Posa-like viruses have been detected in the fecal samples of several host species and are considered unclassified members of Picornavirales. Here, we identified genomic fragments of novel posa-like viruses (monsaviruses) in monkey specimens through next generation sequencing and obtained 11 full-length genomes. This monsavirus shared 88.5-89.2% nucleotide similarity with the Tottori-HG1 strain (GenBank accession LC123275). In total, 713 nucleotide polymorphism sites were identified, indicating their persistent evolution during circulation. The genomic organization and phylogenetic relationship of monsavirus were determined. Subsequent phylogenetic analysis of the conserved replication block of Hel-Pro-RdRp and core RNA-dependent RNA polymerase domain-based analysis of posa-like viruses showed significant separation compared with other known families. Further, posa-like virus genomes possessed the classical replication block of picornavirus in the 5' part of genome and picorna-like capsid domains at the structural coding region of 3' part of genome. Based on these results, we proposed the new family Posaliviridae, within Picornavirales. Four genera, which showed 68.6-75.5% amino acid distances but similar genomic organization including the conserved replication block of Hel-Pro-RdRp, the same order of the genomic coding region, and picorna-like capsid domains, were identified. The flexible genomic organization strategy and a large evolutionary scale of Posaliviridae was explicit. This study provides novel information on monsaviruses and important taxonomic data for the family Posaliviridae.