Viral metagenomics reveals significant viruses in the genital tract of apparently healthy dairy cows

High prevalence of Enterovirus E, Bovine Kobuvirus, and Astrovirus revealed by viral metagenomics in fecal samples from cattle in Central Colombia

Livestock plays a crucial role in ensuring food security and driving the global economy. However, viral infections can have far-reaching consequences beyond economic productivity, affecting the health of cattle, as well as posing risks to human health and other animals. Identifying viruses present in fecal samples, a primary route of pathogen transmission, is essential for developing effective prevention, control, and surveillance strategies. Viral metagenomic approaches offer a broader perspective and hold great potential for detecting previously unknown viruses or uncovering previously undescribed agents. Ubaté Province is Colombia's dairy capital and a key center for livestock production in the country. Therefore, the purpose of this study was to characterize viral communities in fecal samples from cattle in this region. A total of 42 samples were collected from three municipalities in Ubaté Province, located in central Colombia, using a convenient non-probabilistic sampling method. We utilized metagenomic sequencing with Oxford Nanopore Technologies (ONT), combined with diversity and phylogenetic analysis. The findings revealed a consistent and stable viral composition across the municipalities, primarily comprising members of the Picornaviridae family. At the species level, the most frequent viruses were Enterovirus E (EVE) and Bovine Astrovirus (BoAstV). Significantly, this study reported, for the first time in Colombia, the presence of viruses with veterinary importance occurring at notable frequencies: EVE (59%), Bovine Kobuvirus (BKV) (52%), and BoAstV (19%). Additionally, the study confirmed the existence of Circular replicase-encoding single-stranded (CRESS) Virus in animal feces. These sequences were phylogenetically grouped with samples obtained from Asia and Latin America, underscoring the importance of having adequate representation across the continent. The virome of bovine feces in Ubaté Province is characterized by the predominance of potentially pathogenic viruses such as BoAstV and EVE that have been reported with substantial frequency and quantities. Several of these viruses were identified in Colombia for the first time. This study showcases the utility of using metagenomic sequencing techniques in epidemiological surveillance. It also paves the way for further research on the influence of these agents on bovine health and their frecuency across the country.

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.

Comparison of gut viral communities in children under 5 years old and newborns

OBJECTIVES

The gut virome of humans is mainly composed of bacteriophages and their role in shaping the gut microbiome and influencing human health is increasingly recognized. However, little is known about the dynamic changes of the gut virome in children and its role in growth and development. In this study, we collected fecal samples from newborns and children under 5 years old from the same area during the same time period to investigate the gut viral community using viral metagenomic technique.

METHODS

We used viral metagenomics to compare the gut bacteriophage composition between newborns and children under 5 years of age. We collected fecal samples from 45 newborns who were born at the Affiliated Hospital of Jiangsu University and 45 healthy children who were examined at the same hospital. The two groups were classified as the newborn group and the children group.

RESULTS

Our sequencing analysis showed that the number of seqeunce reads of the children group were more than that of the newborn group. The results of alpha diversity and beta diversity both indicated that the diversity of the children group was significantly higher than that of the newborn group and the children group is different from the newborn group. The abundance of gut virome in the children group was also higher than that in the newborn group. The analysis of the genetic characteristics of the viruses showed that the phage genome was scattered and clustered with specificity.

CONCLUSION

Our findings indicate that the gut bacteriophage communities undergo changes over time, presenting diversity and dynamic characteristics. We characterized the composition of gut virome in children and newborns in this region. However, further research is needed to investigate the function of bacteriophages in the ecology of the gastrointestinal tract.

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.

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.

Comparison of gut viral communities in diarrhoea and healthy dairy calves

Calf diarrhoea has been a major cause of economic losses in the global dairy industry. Many factors, including multiple pathogen infections, can directly or indirectly cause calf diarrhoea. This study compared the faecal virome between 15 healthy calves and 15 calves with diarrhoea. Significantly lower diversity of viruses was found in samples from animals with diarrhoea than those in the healthy ones, and this feature may also be related to the age of the calves. Viruses belonging to the families Astroviridae and Caliciviridae that may cause diarrhoea in dairy calves have been characterized, which revealed that reads of caliciviruses and astroviruses in diarrhoea calves were much higher than those in healthy calves. Five complete genomic sequences closely related to Smacoviridae have been identified, which may participate in the regulation of the gut virus community ecology of healthy hosts together with bacteriophages. This research provides a theoretical basis for further understanding of known or potential enteric pathogens related to calf diarrhoea.

Identification and Characterization of Circular Single-Stranded DNA Genomes in Sheep and Goat Milk

In recent years, a variety of circular replicase-encoding single-stranded (CRESS) DNA viruses and unclassified virus-like DNA elements have been discovered in a broad range of animal species and environmental samples. Key questions to be answered concern their presence in the human diet and their potential impact on disease emergence. Especially DNA elements termed bovine meat and milk factors (BMMF) are suspected to act as co-factors in the development of colon and breast cancer. To expand our knowledge on the occurrence of these potential pathogens in human nutrition, a total of 73 sheep and 40 goat milk samples were assayed by combining rolling circle amplification (RCA), PCR and Sanger sequencing. The present study further includes retail milk from the aforementioned species. We recovered 15 single stranded (ss) circular genomes. Of those, nine belong to the family Genomoviridae and six are members of the unclassified group of BMMF. Thus, dairy sheep and goats add to dispersal of CRESS viruses and circular ssDNA elements, which enter the food chain via milk. The presence of these entities is therefore more widespread in Bovidae than initially assumed and seems to be part of the common human nutrition.

The genetic diversity of "papillomavirome" in bovine teat papilloma lesions

Background

Papillomaviruses are small nonenveloped, circular double-stranded DNA viruses that belong to the Papillomaviridae family. To date, 29 Bos taurus papillomavirus (BPV) types have been described. Studies involving mixed BPV infections have rarely been reported in contrast to human papillomavirus (HPV), which is commonly described in numerous studies showing coinfections. Moreover, previous studies had shown that HPV coinfections increase the risk of carcinogenesis. In the present study, we used rolling-circle amplification followed by a high-throughput sequencing (RCA-HTS) approach in 23 teat papillomas from southern Brazil.

Results

Eleven well-characterized BPV types and 14 putative new BPV types were genetically characterized into the Xi, Epsilon and Dyoxipapillomavirus genera according to phylogenetic analysis of the L1 gene, which expands the previous 29 BPV types to 43. Moreover, BPV coinfections were detected in the majority (56.3%) of the papilloma lesions analyzed, suggesting a genetic diverse “papillomavirome” in bovine teat warts.

Conclusions

The data generated in this study support the possibility that a wide range of BPV is probably underdetected by conventional molecular detection tools, and that BPV coinfections are underestimated and probably genetic diverse. Additionally, 14 new BPV types were characterized, increasing the knowledge regarding BPV genetic diversity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00114-3.

Complete genome and phylogenetic analysis of bovine papillomavirus type 15 in Southern Xinjiang dairy cow

Background

Bovine papilloma is a neoplastic disease caused by bovine papillomaviruses (BPVs), which were recently divided into 5 genera and at least 24 genotypes.

Objectives

The complete genome sequence of BPV type 15 (BPV Aks-02), a novel putative BPV type from skin samples from infected cows in Southern Xinjiang China, was determined by collecting warty lesions, followed by DNA extraction and amplicon sequencing.

Methods

DNA was analyzed initially by polymerase chain reaction (PCR) using the degenerate primers FAP59 and FAP64. The complete genome sequences of the BPV Aks-02 were amplified by PCR using the amplification primers and sequencing primers. Sequence analysis and phylogenetic analysis were performed using bio-informatic software.

Results

The nucleotide sequence of the L1 open reading frame (ORF) of BPV Aks-02 was 75% identity to the L1 ORF of BPV-9 reference strain from GenBank. The complete genome consisted of 7,189 base pairs (G + C content of 42.50%) that encoded 5 early (E8, E7, E1, E2, and E4) and 2 late (L1 and L2) genes. The E7 protein contained a consensus CX₂CX₂₉CX₂C zinc-binding domain and a LxCxE motif. Among the different members of this group, the percentages of the complete genome and ORFs (including 5 early and 2 late ORFs) sequence identity of BPV Aks-02 were closer to the genus Xipapillomavirus 1 of the Xipapillomavirus genus. Phylogenetic analysis and sequence similarities based on the L1 ORF of BPV Aks-02 revealed the same cluster.

Conclusions

The results suggest that BPV type (BPV Aks-02) clustered with members of the Xipapillomavirus genus as BPV 15 and were closely related to Xipapillomavirus 1.

Interrogating the bovine reproductive tract metagenomes using culture-independent approaches: a systematic review

Undesirable microbial infiltration into the female bovine reproductive tracts, for example during calving or mating, is likely to disturb the commensal microflora. Persistent establishment and overgrowth of certain pathogens induce reproductive diseases, render the female bovine reproductive tract unfavourable for pregnancy or can result in transmission to the foetus, leading to death and abortion or birth abnormalities. This review of culture-independent metagenomics studies revealed that normal microflora in the female bovine reproductive tract is reasonably consistently dominated by bacteria from the phyla Bacteroidetes, Firmicutes, Proteobacteria, following by Actinobacteria, Fusobacteria and Tenericutes. Reproductive disease development in the female bovine reproductive tract was demonstrated across multiple studies to be associated with high relative abundances of bacteria from the phyla Bacteroidetes and Fusobacteria. Reduced bacterial diversity in the reproductive tract microbiome in some studies of cows diagnosed with reproductive diseases also indicated an association between dysbiosis and bovine reproductive health. Nonetheless, the bovine genital tract microbiome remains underexplored, and this is especially true for the male genital tract. Future research should focus on the functional aspects of the bovine reproductive tract microbiomes, for example their contributions to cattle fertility and susceptibility towards reproductive diseases.

Circular Rep-Encoding Single-Stranded DNA Sequences in Milk from Water Buffaloes (Bubalus arnee f. bubalis)

Isolation and characterization of circular replicase-encoding single-stranded (ss) DNA from animal, plant and environmental samples are rapidly evolving in virology. We detected 21 circular DNA elements, including one genomoviral sequence, in individual milk samples from domesticated Asian water buffaloes (Bubalus arnee f. bubalis). Most of the obtained genomes are related to Sphinx 1.76 and Sphinx 2.36 sequences and share a high degree of similarity to recently published circular DNAs—named BMMF (bovine meat and milk factors)—that have been isolated from commercial milk, as well as from bovine serum. Characteristic features such as rep genes, tandem repeats and inverted repeats were detected. These BMMF have recently been found to be present in taurine-type dairy cattle breeds descending from the aurochs (Bos primigenius). Importantly, the occurrence of BMMF has been linked to the higher incidence of colorectal and breast cancer in North America and Western Europe compared with Asia. This is the first report of circular ssDNA detected in milk from the domesticated form of the wild Asian water buffalo (B. arnee) belonging to the subfamily Bovinae. This novelty should be taken into account in view of the above-mentioned cancer hypothesis.

Identification and Characterization of a Novel Recombinant Porcine Astrovirus from Pigs in Anhui, China

Porcine astroviruses (PAstVs) have wide distribution in swine herds worldwide. At present, five porcine astrovirus genotypes have been identified. In this study, using viral metagenomics, a novel PAstV strain (designated as Ahast) was identified in fecal samples from pigs in Anhui of China, and the complete genomic sequence of Ahast was obtained by assembling and PCR amplification. Genomic structural analysis indicated that Ahast had a typical ribosomal frameshifting signal, and some conserve amino acid motifs were also found in virally encoded proteins. Phylogenetic analysis and sequence comparison indicated that this virus belonged to porcine astrovirus genotype 4 (PAstV4), which formed a clade clustered with other PAstV4. Multiple recombinant events were confirmed by recombination analysis and indicated that Ahast was a potential recombinant. Epidemiological investigation indicated that PAstV4 has a 10.7% prevalence in this pig farm. The new recombinant identified in this study will be beneficial to comprehend the origin, genetic diversity, and evolution of porcine astroviruses in Anhui of China.

Anogenital-Associated Papillomaviruses in Animals: Focusing on Bos taurus Papillomaviruses

In contrast to the diverse studies on human papillomaviruses (HPVs), information on animal PVs associated with anogenital lesions is limited. In the animal kingdom, papillomas occur more commonly in cattle than in any other animals, and diverse types of Bos taurus papillomaviruses (BPVs) exist, including the very recently discovered BPV type 29 (BPV29). From this perspective, we will review previous studies describing PV types associated with anogenitals in animals, with a focus on BPVs. To date, two classical BPV types, classified into Deltapapillomavirus (BPV1 and BPV2) and Dyokappapapillomavirus (BPV22), and two novel Xipapillomaviruses (BPV28 and BPV29) have been identified from anogenital lesions and tissues of the domestic cow. Due to the limited reports describing anogenital-associated PVs in animals, the relationships between their phylogenetic and pathogenetic properties are still undiscovered. Animal studies are valuable not only for the veterinary field but also for human medicine, as animal diseases have been shown to mimic human diseases. Studies of anogenital-associated PVs in animals have a positive impact on various research fields.

Fecal, oral, blood and skin virome of laboratory rabbits

Here, we investigated the fecal, oral, blood, and skin virome of 10 laboratory rabbits using a viral metagenomic method. In the oral samples, we detected a novel polyomavirus (RabPyV), and phylogenetic analysis based on the large T antigen, VP1 and VP2 regions indicated that the novel strain might have undergone a recombination event. Recombination analysis based on related genomes confirmed that RabPyV is a multiple recombinant between rodent-like and avian-like polyomaviruses. In fecal samples, three partial or complete genome sequences of viruses belonging to the families Picobirnaviridae, Parvoviridae, Microviridae and Coronaviridae were characterized, and phylogenetic trees were constructed based on the predicted amino acid sequences of viral proteins. This study increases the amount of genetic information on viruses present in laboratory rabbits.

Virus Metagenomics in Farm Animals: A Systematic Review

A majority of emerging infectious diseases are of zoonotic origin. Metagenomic Next-Generation Sequencing (mNGS) has been employed to identify uncommon and novel infectious etiologies and characterize virus diversity in human, animal, and environmental samples. Here, we systematically reviewed studies that performed viral mNGS in common livestock (cattle, small ruminants, poultry, and pigs). We identified 2481 records and 120 records were ultimately included after a first and second screening. Pigs were the most frequently studied livestock and the virus diversity found in samples from poultry was the highest. Known animal viruses, zoonotic viruses, and novel viruses were reported in available literature, demonstrating the capacity of mNGS to identify both known and novel viruses. However, the coverage of metagenomic studies was patchy, with few data on the virome of small ruminants and respiratory virome of studied livestock. Essential metadata such as age of livestock and farm types were rarely mentioned in available literature, and only 10.8% of the datasets were publicly available. Developing a deeper understanding of livestock virome is crucial for detection of potential zoonotic and animal pathogens and One Health preparedness. Metagenomic studies can provide this background but only when combined with essential metadata and following the “FAIR” (Findable, Accessible, Interoperable, and Reusable) data principles.

FUNDC1-mediated mitophagy in bovine papillomavirus-infected urothelial cells

E5 protein, the major oncoprotein of the bovine Deltapapillomavirus genus, has been detected in 17 of the 19 urothelial cancers by molecular and morphological procedures. In 10 urothelial cancers, the oxygen sensitive subunit HIF-1α, which is upregulated by hypoxia, was overexpressed. Mitophagy, the selective autophagic removal of dysfunctional mitochondria, was upregulated in hypoxic neoplastic cells infected by BPVs which was mediated by FUNDC1, a mitochondrial outer-membrane protein. The FUNDC1 receptor was amplified by PCR, and amplicon sequencing showed a 100% homology with bovine FUNDC1 sequences deposited in GenBank (accession number: NM_001104982). Both transcripts and protein levels of FUNDC1 were significantly decreased in hypoxic neoplastic cells relative to healthy, non-neoplastic cells. FUNDC1 interacted with the LC3 protein, a marker of autophagosome (mitophagosome) membrane, the Hsc70/Hsp70 chaperone, and Bag3 co-chaperone. Bag3 may play a role in mitophagosome formation together with the Synpo2 protein, and may be involved in the degradation of Hsc70/Hsp70-bound CHIP-ubiquitinated cargoes, in association with its chaperone. Ultrastructural findings revealed the presence of mitochondria exhibiting severe fragmentation and loss of cristae, as well as numerous mitochondria-containing autophagosomes. Total and phosphorylated GTPase dynamin-related protein 1 (DRP1), which plays a crucial role in mitochondrial fission, a pre-requisite for mitophagy, was overexpressed at the mitochondrial level. Total and phosphorylated mitochondrial fission factor (Mff), mitochondrial fission protein 1 (Fis1), mitochondrial dynamics 51 (MiD51), and MiD49, which are DRP1 receptors responsible and/or co-responsible for its mitochondrial recruitment were overexpressed.

Bovine papillomavirus E5 oncoprotein expression and its association with an interactor network in aggresome-autophagy pathway

E5 protein, the major oncoprotein of bovine Deltapapillomavirus (BPV), was found to be expressed in 18 of 21 examined urothelial cancers of cattle. E5 oncoprotein was found to interact with p62 which was degraded through the autophagosome-lysosome pathway as well as LC3-II and appeared to be involved in the phosphorylation of the α-subunit of eukaryotic initiation factor 2 (eIF2α). Autophagy was morphologically documented by transmission electron microscope (TEM) through the detection of double-membrane autophagosomes and autolysosomes. Overexpression of Bag3 known to mediate selective autophagy was also demonstrated. Furthermore, Bag3 and BPV E5 oncoprotein were seen to co-localize with dynein and 14-3-3γ, which suggested that Bag3 could be involved in inducing the retrograde transport of BPV E5 along microtubules to aggresomes, perinuclear sites with high autophagic flux. Electron dense perinuclear structures consistent with aggresomes were also documented by TEM in urothelial cancer cells. Finally, Bag3 was found to also interact with synaptopodin 2 (Synpo2), which would seem to contribute to cargo degradation as it has been shown to facilitate autophagosome formation. This study provides mechanistic insights into the potential role(s) of autophagy in BPV disease, which can help to develop future treatment and control measures for BPV infection. Activation of autophagy correlates positively with BPV infection and may play a role in biological behavior of bladder cancer as urothelial carcinomas of cattle are known to be characterized by a relatively low rate of metastasis.