Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Vet. Sci., 16 September 2025

Sec. Veterinary Infectious Diseases

Volume 12 - 2025 | https://doi.org/10.3389/fvets.2025.1677378

First detection of two cycloviruses in cormorant fecal samples in China by high-throughput sequencing technology


Yifei PeiYifei Pei1Shiyu CaiShiyu Cai1Yue XueYue Xue1Yutong Fu,Yutong Fu1,2Jingyue ZhangJingyue Zhang3Quan ShenQuan Shen1Likai JiLikai Ji1Ping WuPing Wu1Hua WangHua Wang1Yan Wang
Yan Wang1*Wen Zhang
Wen Zhang1*Shixing Yang
Shixing Yang1*
  • 1School of Medicine, Jiangsu University, Zhenjiang, Jiangsu, China
  • 2Department of Rehabilitation, Danyang Hospital of Traditional Chinese Medicine, Danyang Jiangsu, China
  • 3Department of Biology, Shenzhen MSU-BIT University, Shenzhen, Guangdong, China

Introduction: The Great Cormorant (Phalacrocorax carbo) is widely distributed across China. As an apex predator in aquatic ecosystems, it plays a tripartite ecological role: acting as a natural host, transmission vector, and indicator species for viruses. Current research confirms that cormorants carry diverse viral pathogens from the families including Flaviviridae, Orthomyxoviridae, Paramyxoviridae, and Polyomaviridae. Significant knowledge gaps persist regarding their virome diversity.

Methods: In this study,46 cormorant fecal swab samples were collected at Xiamen Garden Expo Park, and viralmetagenomics method was conducted to identify two Cycloviruses.

Results: This study identified two novel cycloviruses, Corcyclo-1 (1,856 bp) and Corcyclo-2 (1,831 bp), from cormorant fecal samples using viral metagenomics. Genomic analyses revealed hallmark features of the genus Cyclovirus, including inversely oriented open reading frames (ORFs) encoding the capsid protein (Cap) and replication-associated protein (Rep), as well as a conserved stem-loop sequence TAATACTAT. The Rep gene of Corcyclo-1 contained a 166-bp intron and shared >96.9% amino acid identity with human-, wild boar-, and chicken-derived cyclovirus strains (HaCV-8) from Vietnam and Madagascar, classifying it as a novel strain of HaCV-8. In contrast, Corcyclo-2 harbored a 98-bp intron in its Rep gene and clustered with unclassified cyclovirus strains from bats and mongooses in China and Saint Kitts and Nevis (>97.4% identity), constituting a putative new species. Phylogenetic and pairwise sequence analyses further supported their taxonomic positions. Epidemiological screening demonstrated a high prevalence of Corcyclo-1 (82.6%, 38/46) and Corcyclo-2 (32.6%, 15/46) in cormorant feces. Cross-species surveillance detected Corcyclo-2 in chickens (25.8%, 16/62) and ducks (11.7%, 9/77), whereas Corcyclo-1 was absent in these hosts.

Conclusion: This study represents the first report of cormorant-associated cycloviruses, highlighting their potential for cross-species transmission and providing new insights into the ecological diversity and evolutionary mechanisms of cyclovirus.

Introduction

The Great Cormorant is a widely distributed waterbird across global coastal and inland aquatic ecosystems (1, 2). According to monitoring data from BirdLife International, the global Great Cormorant population is estimated at 1.2 to 2.4 million individuals. However, overfishing in some regions has caused a sharp decline in their food resources. In China, coastal pollution in traditional wintering grounds such as Bohai Bay has resulted in local population declines exceeding 30%. Mass infections by cestode parasites can significantly increase juvenile mortality rates (3). Despite the Great Cormorant's crucial ecological role in marine food chains, systematic research on its virome remains lacking.

The family Circoviridae comprises covalently closed circular single-stranded DNA (ssDNA) viruses, representing the smallest known animal-infecting viral group (4, 5). They are classified within two genera, Circovirus and Cyclovirus (4, 6). Circoviruses exhibit cross-species transmission capabilities among Anseriformes (waterfowl) and various other aquatic birds, causing clinical manifestations including immunosuppression, developmental disorders, and digestive disturbances (79). In poultry, chickens, turkeys, and guinea fowl are highly susceptible to Avian circovirus. Primarily transmitted via the fecal-oral route, this virus causes characteristic pathological lesions such as bursal atrophy and feather dystrophy (10). Wild bird surveillance indicates Circoviridae family members exhibit a 34% positive detection rate in Columbiformes (e.g., the Hill Pigeon, Columba rupestris), with infected individuals experiencing persistent weight loss and diminished flight capacity (11, 12). Although circovirus prevalence in Anseriformes remains unclear, their broad host adaptability suggests a potential zoonotic transmission risk to Bar-headed Goose (Anser indicus) populations (8). Moreover, a previous study reported that a novel human circovirus, HuCV2, was discovered in intravenous drug users in China, showing the highest homology with porcine circovirus type 3 (PCV3) (13). This has heightened concerns about the cross-species transmission of circoviruses.

The genus Cyclovirus, classified within the family Circoviridae, comprises small, non-enveloped viruses with covalently closed circular single-stranded DNA (ssDNA) genomes (1.7–2.1 kb) (4). First identified in 2010, cycloviruses share structural similarities with circoviruses, including bidirectional open reading frames (ORFs) encoding the capsid (Cap) and replication-associated (Rep) proteins, as well as a conserved stem-loop motif (TAATACTAT) near the origin of replication (14). However, cycloviruses exhibit distinct genomic features such as frequent introns in the Rep gene and lower overall nucleotide identity with circoviruses (<70%), justifying their classification as a separate genus (15).

Cycloviruses have been detected in diverse hosts, including humans, livestock, wildlife (e.g., bats, rodents, and birds), and environmental samples, suggesting broad host tropism and potential zoonotic transmission (16, 17). Notably, human-associated cycloviruses (e.g., HaCV-8) have been linked to neurological and respiratory infections, though their pathogenicity remains poorly understood (14). In animals, cyclovirus infections are often subclinical, but their role in immunosuppression or coinfections warrants further investigation.

Despite their ecological prevalence, cyclovirus research lags behind circoviruses due to challenges in culturing and the lack of standardized detection protocols. Wild birds, particularly aquatic species, are underrepresented in cyclovirus surveillance, leaving gaps in understanding their evolutionary dynamics and cross-species transmission risks. This study addresses this gap by characterizing novel cycloviruses in Great Cormorants, a keystone species in aquatic ecosystems.

Viral metagenomics technology has been successfully applied for pathogen screening in various waterfowl species, offering a viable approach to study viral diversity and potential health risks in this species (18, 19). This study employed viral metagenomics to analyze 46 fecal samples from wild Great Cormorants collected in Xiamen, Fujian Province during 2021. This study reports the first full-genome characterization of a great cormorant-associated circovirus, revealing molecular evidence that wetland waterbirds harbor viral genes homologous to human/mammalian viruses. The complete viral genome sequences obtained provide essential data for establishing molecular detection standards targeting the Cap gene for avian circoviruses in southeastern coastal China. This has significance for monitoring viral diseases in Great Cormorants and assessing zoonotic transmission risk.

Materials and methods

Sample collection and preparation

In 2021, 46 cormorant fecal swab samples were collected at Xiamen Garden Expo Park, with species identification confirmed by ornithologists during collection (Figure 1). The swabs were soaked in 1 mL of Dulbecco's phosphate-buffered saline (DPBS) and vortexed for 10 min. Following centrifugation, the supernatant was moved to a 1.5 mL tube and kept at −80 °C for later analysis. The research was authorized by Jiangsu University's Ethics Committee (Approval No. 2018UJS18023), and sampling followed China's Wildlife Protection Law. The following experimental procedures were conducted following the methods described by Yijie Sun et al. (8).

Figure 1
Map showing Fujian Province with an enlarged section of Xiamen City highlighted. A red dotted line connects the province to Xiamen City, where the Xiamen Garden Expo Park is marked.

Figure 1. The geographical position of the sample site (red dot) in Xiamen, Fujian, China.

Viral nucleic acid extraction and library construction

A pool of 9–10 swab extracts was used to construct each library, generating five libraries (swab39, swab40, swab41, swab42, and swab43). Each library used 500 μL of swab suspension (50–55.6 μL supernatant per sample), which was filtered using a 0.45-μm filter (Merck Millipore, USA) to exclude bacteria and eukaryotic cells. The filtrate underwent nuclease treatment (Turbo DNase from Thermo Fisher Scientific, USA, BaselineZero DNase from the US epicenter, Benzoase nuclease from Novagen Corporation, USA, and RNase A from Thermo Fisher Scientific) at 37 °C for 90 min to digest exposed nucleic acids. The QIAamp MinElute Virus Spin Kit (Qiagen, Germany) was employed for viral nucleic acid extraction, followed by concentration measurement with a Qubit 4 fluorometer (Invitrogen, USA). Purified nucleic acids were stored at −80 °C.

The combined viral nucleic acids (DNA/RNA) were reverse transcribed using SuperScript III Reverse Transcriptase (Invitrogen, USA) with 100 pmol random hexamers. The RT protocol included: 25 °C for 10 min, 50 °C for 60 min, 85 °C for 5 min, and 95 °C for 2 min. Products were rapidly cooled and kept on ice for 2 min. cDNA synthesis was completed with Klenow fragment (NEB, USA) at 37 °C for 60 min and 75 °C for 20 min.

The Nextera XT DNA Sample Preparation Kit (Illumina, USA) was used for library preparation following standard protocols, including adapter ligation and 15 cycles of limited amplification. Personalbio conducted sequencing on the Illumina NovaSeq platform, generating 250-bp paired-end reads with dual-indexing.

Bioinformatic processing pipeline

The NovaSeq platform produced 250-bp paired-end reads, which were demultiplexed using Illumina's software and analyzed through a custom bioinformatics pipeline running on a 32-node Linux cluster. Adapters were removed using VecScreen with default settings. Sequence ends with Phred scores below 10 were removed. Bowtie2 v2.2.4 aligned reads against the BLAST NT database to eliminate bacterial sequences (20). Remaining reads were reconstructed with SOAPdenovo2 (k-mer = 63, default settings) (21). Contigs and singlets were matched to a custom viral proteome database through BLASTx analysis (E-value cutoff <10−5). This viral database was curated from NCBI's viral reference proteome FASTA files (https://ftp.ncbi.nih.gov/refseq/release/viral/) using taxonomy annotation from the viral kingdom. Potential viral sequences underwent additional screening against an NVNR protein database to exclude false identifications.

Sequence alignment and ORF prediction

Pairwise alignment of viral amino acid sequences was performed using SDTv1.2 with MUSCLE parameters and Geneious 11.1.2 with MUSCLE alignment. ORFs were identified using Geneious 11.1.2 and NCBI ORF Finder.

Phylogenetic analysis

Phylogenetic analysis examined evolutionary relationships using amino acid sequences of circoviruses identified in this study, along with top BLASTx matches and representative viral strains (22). Sequence alignments were performed with MUSCLE in MEGA X. MrBayes v3.2.7 generated trees with mixed amino acid models (“prset aamodelpr = mixed”). The analysis ran for 10 million generations (sampling every 50), discarding initial 25% as burn-in. Convergence was achieved at split frequency SD < 0.01, with bootstrap values assigned to nodes.

Primer development and cyclovirus detection by PCR

Two sets of primers were designed based on the Cap gene of Corcyclo-1 and Corcyclo-2. Using these primers, we performed PCR on DNA extracted from fecal samples of cormorants, chickens, and ducks. PCR reactions (25 μL total volume) contained 3 μL DNA template, 10 μL HotStarTaq Master Mix, 2 μL each primer (20 pmol), and ddH2O. Controls included positive and negative samples (extraction/water controls). Cycling parameters: 95 °C for 2 min; 30 cycles of 94 °C (20 s), 55 °C (30 s), 72 °C (30 s); final extension at 72 °C for 5 min. Amplicons were bidirectionally sequenced ≥2 times via Sanger sequencing (Sangon Biotech, China).

Nucleotide sequence accession number

The complete viral genome sequences identified in this study were deposited in GenBank under the accession numbers VP731248 and VP731249.

Results

Viral metagenomic overview

Illumina sequencing generated 13,877,028 raw reads across five libraries (Supplementary Table S1). After bioinformatics analysis, a total of 568,343 sequences reads matching viral proteins, accounting for 4.1% of the total number of raw reads. The majority of the eukaryotic viruses identifies in this study including the families of Circoviridae (48.5% of the total analyzed virus reads), Picornaviridae (31.1%), Astroviridae (12.8%), Parvoviridae (5.2%), Caliciviridae (1.7%), Genomoviridae (0.4%), and Smacoviridae (0.2%) (Figure 2, Supplementary Table S2). Considering the potential pathogenic role, and the fact that no circovirus has been detected in cormorants to date, we chose the circoviruses for future comprehensive study.

Figure 2
Pie chart showing the distribution of virus families. Circoviridae holds the largest portion at 48.5%, followed by Picornaviridae at 31.1%, Astroviridae at 12.8%, Parvoviridae at 5.2%, Caliciviridae at 1.7%, Genomoviridae at 0.4%, and Smacoviridae at 0.2%.

Figure 2. Taxonomic analysis of the fecal virome detected in cormorants at the family level.

The genomic analysis of two novel cycloviruses

In this study, De novo assemble yielded two complete cyclovirus genomes: Corcyclo-1 (1,856 bp) and Corcyclo-2 (1,831 bp). Both genomes exhibited hallmark features of the genus Cyclovirus including two inversely oriented open reading frames encoding the Cap protein and Rep protein (23). The GC contents of Corcyclo-1 and Corcyclo-2 are 47.2 % and 47.5 %, respectively. For Corcyclo-1, ORF1 is 669 bp in length encoding a 222-aa Cap protein, while ORF2 is 879 bp in length encoding a 292-aa Rep protein. For Corcyclo-2, Cap protein encoded by ORF1 is 232 aa in length, while Rep protein is 277 aa (Figures 3A, C). The Rep genes of both cycloviruses were predicted to be interrupted by a 166-bp or 98-bp intron with a typical splice donor (GT) and splice acceptor site (AG) respectively. The conserved motifs of cyclovirus (WWDGY, DDFYGW, DRYP, CSK) were presented in the Rep proteins of both Corcyclo-1 and Corcyclo-2. The typical nonamer sequence of the stem loop structure (TAATACTAT) was conserved in the genomes of Corcyclo-1 and Corcyclo-2 (Figures 3B, D). Despite structural similarities, Corcyclo-2 lacked ORF overlap due to a premature termination codon in its Rep gene.

Figure 3
Circular genomic maps and structures of two plasmids. Panel A shows Corcyclo-1 with a size of 1,856 base pairs, featuring capsid and replication regions. Panel B is a hairpin loop diagram for Corcyclo-1's DNA. Panel C displays Corcyclo-2, measuring 1,831 base pairs, also highlighting capsid and replication regions. Panel D is a hairpin loop diagram for Corcyclo-2's DNA. Both diagrams illustrate genetic features with base pair positions and labels.

Figure 3. The genomic organization and stem-loop structures of cycloviruses identified in cormorants. (A) The genomic organization of Corcyclo-1 identified in cormorants; (B) the stem-loop structure of Corcyclo-1 identified in cormorants; (C) the genomic organization of Corcyclo-2 identified in cormorants; (D) the stem-loop structure of Corcyclo-2 identified in cormorants.

Phylogenetic analysis

Two phylogenetic trees were constructed using Cap and Rep protein sequences from Corcyclo-1 and Corcyclo-2, and reference strains from the Cyclovirus species (Figures 4A, B). For Corcyclo-1, Cap phylogeny clustered it with other cyclovirus strains isolated from cerebrospinal fluid and feces samples from human (GenBank no. KF031465 and KM392288), wild boar (GenBank no. KF031470), and chicken (GenBank no. KF031471). Both of them belong to Human associated cyclovirus 8 (HaCV-8). Rep phylogeny further supported this classification, suggesting cross-species transmission potential. For Corcyclo-2, both Cap and Rep phylogenies grouped it with unclassified cyclovirus strains from feces samples of bat and mongoose.

Figure 4
Two phylogenetic trees labeled A and B, showing the evolutionary relationships of various strains or species. Each tree includes branching diagrams with colored gradient bars representing bootstrap percentages. Names of strains are listed alongside each branch. Tree A and B display different groupings, with some branch names highlighted in red. The image likely represents a comparison of genetic data.

Figure 4. Phylogenetic analysis of cycloviruses identified in cormorants. (A) phylogenetic tree constructed based on Cap proteins; (B) phylogenetic tree constructed based on Rep proteins. Viruses identified in this study are marked with red.

Pairwise alignment of amino acid sequences and nucleotide sequences

In order to determine the classification of these cycloviruses detected in this study, the viral encoded proteins and whole genome sequences were compared with the corresponding sequences of known cycloviruses. For Corcyclo-1, Cap protein alignment revealed >96.9% amino acid identity with cyclovirus strains from humans, wild boars, and chickens in Vietnam and Madagascar, and <92% identity with other representative strains (Figure 5A). Rep alignment showed >97.2% identity with strains from humans, wild boars, and chickens in Vietnam and Madagascar, but <92% identity with other strains (Figure 5B). Whole-genome alignment exhibited >96.9% nucleotide identity with strains from humans, wild boars, and chickens in Vietnam and Madagascar, and <86.3% identity with others strains (Figure 5C). For Corcyclo-2, Cap alignment indicated >97.4% amino acid identity with unclassified cyclovirus strains from bat and mongoose in China and Saint Kitts and Nevis, but <53.5% identity with other representatives (Figures 5A, B). Rep alignment showed >99.3% identity with unclassified cyclovirus strains from bat and mongoose in China and Saint Kitts and Nevis, while <91.4% identity with other strains. Whole-genome alignment exhibited >96.9% nucleotide identity with Cyclovirus strains from bat and mongoose in China and Saint Kitts and Nevis, and <88.9% with others strains (Figure 5C). Viral taxonomy standards require ≤ 80% whole-genome nucleotide similarity and consistent topological positioning of the replication origin (ori) relative to Rep/Cap coding regions for species classification within a genus (24). Based on this classification standard, Corcyclo-1 should be classified as Human associated cyclovirus 8 of the Hepatovirus genus, while Corcyclo-2 is classified as an unclassified cyclovirus.

Figure 5
Triangular heatmaps labeled A, B, and C show pairwise data correlations. Each heatmap features a color gradient from blue to red, representing low to high correlation. Annotations include text surrounding each heatmap.

Figure 5. Pairwise comparison of nucleotide and amino acid sequences of cycloviruses identified in cormorants with representative strains of different species and genera in the Circoviridae family. (A) pairwise comparison based on Cap proteins; (B) pairwise comparison based on Rep proteins; (C) pairwise comparison based on nucleotide sequences of full genomes.

Cyclovirus prevalence estimation

To investigate the prevalence of Corcyclo-1 and Corcyclo-2, viral DNA genomes were extracted from all 46 fecal samples. Two sets of primers, designed based on the Cap gene of Corcyclo-1 and Corcyclo-2, were respectively used to perform PCR screening (Supplementary Table S3). The results indicated that 38 samples were positive for Corcyclo-1 (38/46, 82.6% positive rate), while only 15 samples were positive for Corcyclo-2 (15/46, 32.6% positive rate) (Supplementary Figures S1, S2). Since both cycloviruses have previously been detected in other animal species, we investigated their cross-species prevalence by screening 62 chicken and 77 duck fecal samples from surrounding areas using the above primers. Results showed that Corcyclo-2 was detected in both sample types, with positive rates of 25.8% (16/62) in chickens and 11.7% (9/77) in ducks (Supplementary Figures S3, S4), in contrast, Corcyclo-1 was not detected in fecal samples from chickens and ducks.

Discussion

Birds harbor diverse viromes, and the viral diversity in healthy individuals suggests their role as natural viral reservoirs (18). The Great Cormorant is protected by China's National Forestry and Grassland Administration and widely distributed. Previous studies identified multiple pathogenic viruses infecting this species: Nazerian et al. (25) first isolated a novel herpesvirus from diseased chick blood; Bányai et al. (25) discovered a novel gammapolyomavirus in deceased individuals at Budapest Zoo, Hungary; Varsani et al. (26) identified multiple viruses (including Cressdnavirus, Microvirus, Gyrovirus, and Caudovirus) in fecal samples of Double-crested Cormorants in the United States (27). Currently, systematic surveillance data on viral carriage or infection in Chinese cormorant populations are absent.

The family Circoviridae currently comprises approximately 155 species classified into two genera (Circovirus and Cyclovirus) (6). Breakthroughs in PCR and metagenomic sequencing have significantly expanded our understanding of their diversity (28, 29). Current data reveal that members of the genus Circovirus are widely distributed in human clinical samples (positivity rate 2.9%−7.2%), livestock hosts (15%−34% positivity in domestic pigs), and aquatic environments (89% positivity in wastewater treatment plants), demonstrating remarkable environmental adaptability (3033).

This study identified and characterized two novel cycloviruses in cormorant fecal samples. The discovery of Corcyclo-1 and Corcyclo-2 significantly expands knowledge regarding the genetic diversity and host range of Circoviridae. Full genome sequencing and phylogenetic analysis indicate that Corcyclo-1 belongs to HaCV-8. In this study, HaCV-8 was found to have a high prevalence in cormorants, but it was not detected in the fecal samples of chickens and ducks from the same area. Further epidemiological investigations are needed to determine its pathogenicity in cormorants and assess its potential for transmission to humans. In contrast, evolutionary analysis and sequence comparison showed that Corcyclo-2 discovered in this study is closely related to the unclassified cyclovirus in the feces of wild bats in China and the cyclovirus in the feces of mongooses in Saint Kitts and Nevis, and it is classified as an unknown novel cyclovirus. Epidemiological surveys showed that Corcyclo-2 was also detected in the fecal samples of chickens and ducks from the same area, indicating that it can infect multiple species. However, whether it is pathogenic to the infected animals needs further detailed research. Since Corcyclo-2 was detected in cormorants, chickens, and ducks in this study, and genetically similar strains have previously been identified in bats, we hypothesize that the virus is capable of cross-species transmission. Transmission is likely mediated by the fecal–oral route, and contact transmission cannot be excluded. Therefore, continuous surveillance of Corcyclo-2 infection in wild birds and domestic poultry is essential for the early detection and control of its potential pathogenic threat to both animal and public health.

To date, novel cyclovirus-like viruses have been isolated from diverse species, indicating broad host infectivity. Although no significant avian clinical signs were observed during sampling, undetected clinical disease due to insufficient sampling cannot be excluded. Consequently, further experimental and epidemiological studies are needed to elucidate their pathogenicity and host range. While this study advances understanding of cormorant-associated viruses, limitations remain. First, reliance solely on fecal samples precludes comprehensive viral distribution assessment, as these viruses may inhabit other tissues or biological fluids. Second, despite discovering Corcyclo-1 and Corcyclo-2, functional validation of their pathogenicity and transmission modes under natural conditions is lacking. Finally, sampling spatiotemporal constraints may limit representation of population-wide viral carriage, especially given potential geographical and seasonal influences on viral diversity and abundance.

Conclusion

This study reports the first identification of two novel viruses, Corcyclo-1 and Corcyclo-2, substantially enhancing our knowledge of viral diversity carried by wild cormorants. Both Corcyclo-1 and Corcyclo-2 are members of the family Circoviridae. Corcyclo-1 exhibits significant homology and structural similarity to Human associated cyclovirus 8 (HaCV-8), a virus whose nucleic acid has been detected in human samples. However, the clinical significance and pathogenicity of HaCV-8 in humans remain entirely unknown. The high prevalence of this virus in cormorants warrants further investigation to determine its true host range and to assess the potential risk of zoonotic transmission. Corcyclo-2 is closely related to several unclassified cycloviruses, forming an independent branch and is designated as an unclassified cyclovirus. Epidemiological investigations showed that Corcyclo-2 could also infect chickens and ducks in the same region, while Corcyclo-1 was not detected in fecal samples from chickens and ducks. These findings contribute to understanding the genetic diversity and host range of these viruses. Future research should better expand the cormorant sample size and characterize the epidemiology, pathogenicity, and potential for cross-species transmission of these newly discovered viruses.

Data availability statement

The complete viral genome sequences identified in this study were deposited in GenBank under the accession numbers PV731248 and PV731249.

Ethics statement

The animal study was approved by Jiangsu University Ethics Committee. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

YP: Data curation, Investigation, Writing – original draft, Writing – review & editing. SC: Investigation, Writing – review & editing. YX: Investigation, Writing – review & editing. YF: Writing – review & editing. JZ: Investigation, Writing – review & editing. QS: Data curation, Writing – review & editing. LJ: Data curation, Writing – review & editing. PW: Data curation, Writing – review & editing. HW: Data curation, Writing – review & editing. WZ: Conceptualization, Funding acquisition, Supervision, Writing – review & editing. YW: Conceptualization, Supervision, Writing – review & editing. SY: Conceptualization, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the National Key Research and Development Programs of China 299 (2023YFD1801300) and the National Natural Science Foundation of China (82341106).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2025.1677378/full#supplementary-material

Supplementary Table S1 | The summary of library information of cormorant.

Supplementary Table S2 | The sequence numbers of eukaryotic viruses in family level in different libraries.

Supplementary Table S3 | The summary of primer information of Corcyclo-1 and Corcyclo-2.

Supplementary Figure S1 | PCR screening results for Corcyclo-1 in cormorant fecal samples.

Supplementary Figure S2 | PCR screening results for Corcyclo-2 in cormorant fecal samples.

Supplementary Figure S3 | PCR screening results for Corcyclo-2 in chicken fecal samples.

Supplementary Figure S4 | PCR screening results for Corcyclo-2 in duck fecal samples.

References

1. Rusconi A, Prati P, Bragoni R, Castelli M, Postiglione U, Rigamonti S, et al. Occurrence of Eustrongylides excisus (Nematoda:Dioctophymatidae) in European perch (Perca fluviatilis) and great cormorant (Phalacrocorax carbo) in lake annone, Northern Italy. J Parasitol. (2022) 108:209–16. doi: 10.1645/20-175

PubMed Abstract | Crossref Full Text | Google Scholar

2. van Leeuwen CHA, Lovas-Kiss Á, Ovegård M, Green AJ. Great cormorants reveal overlooked secondary dispersal of plants and invertebrates by piscivorous waterbirds. Biol Lett. (2017) 13:10. doi: 10.1098/rsbl.2017.0406

PubMed Abstract | Crossref Full Text | Google Scholar

3. Davidovich N, Yasur-Landau D, Behar A, Pretto T, Scholz T. Invasive parasites and global change: Evidence for the recent and rapid spillover of a potential pathogen of tilapias with a complex, three-host life cycle. Heliyon. (2023) 9:e18831. doi: 10.1016/j.heliyon.2023.e18831

PubMed Abstract | Crossref Full Text | Google Scholar

4. Breitbart M, Delwart E, Rosario K, Segalés J, Varsani A, Ictv Report C, et al. Virus taxonomy profile: circoviridae. J Gen Virol. (2017) 98:1997–8. doi: 10.1099/jgv.0.000871

PubMed Abstract | Crossref Full Text | Google Scholar

5. Rosario K, Breitbart M, Harrach B, Segalés J, Delwart E, Biagini P, et al. Revisiting the taxonomy of the family Circoviridae: establishment of the genus Cyclovirus and removal of the genus Gyrovirus. Arch Virol. (2017) 162:1447–63. doi: 10.1007/s00705-017-3247-y

PubMed Abstract | Crossref Full Text | Google Scholar

6. Varsani A, Harrach B, Roumagnac P, Benko M, Breitbart M, Delwart E, et al. 2024 taxonomy update for the family Circoviridae. Arch Virol. (2024) 169:176. doi: 10.1007/s00705-024-06107-2

PubMed Abstract | Crossref Full Text | Google Scholar

7. Todd D. Circoviruses: immunosuppressive threats to avian species: a review. Avian Pathol. (2000) 29:373–94. doi: 10.1080/030794500750047126

PubMed Abstract | Crossref Full Text | Google Scholar

8. Sun Y, Wang Y, Ji L, Pei Y, Sun X, Ji L, et al. Identification of multiple novel viruses in bar-headed goose feces from Tibet of China. Front Vet Sci. (2024) 11:1485631. doi: 10.3389/fvets.2024.1485631

PubMed Abstract | Crossref Full Text | Google Scholar

9. Hong YT, Kang M, Jang HK. Pathogenesis of duck circovirus genotype 1 in experimentally infected Pekin ducks. Poult Sci. (2018) 97:3050–7. doi: 10.3382/ps/pey177

PubMed Abstract | Crossref Full Text | Google Scholar

10. Li G, Zhou D, Zhao M, Liu Q, Hao X, Yan T, et al. Kinetic analysis of pathogenicity and tissue tropism of gyrovirus 3 in experimentally infected chickens. Vet Res. (2021) 52:120. doi: 10.1186/s13567-021-00990-2

PubMed Abstract | Crossref Full Text | Google Scholar

11. Nath BK, Das T, Peters A, Gupta SD, Sarker S, Forwood JK, et al. Australasian pigeon circoviruses demonstrate natural spillover infection. Viruses. (2023) 15:2025. doi: 10.3390/v15102025

PubMed Abstract | Crossref Full Text | Google Scholar

12. Zhang Z, Dai W, Wang S, Dai D. Epidemiology and genetic characteristics of pigeon circovirus (PiCV) in eastern China. Arch Virol. (2015) 160:199–206. doi: 10.1007/s00705-014-2255-4

PubMed Abstract | Crossref Full Text | Google Scholar

13. Li Y, Zhang P, Ye M, Tian RR, Li N, Cao L, et al. Novel circovirus in blood from intravenous drug users, Yunnan, China. Emerg Infect Dis. (2023) 29:1015–9. doi: 10.3201/eid2905.221617

PubMed Abstract | Crossref Full Text | Google Scholar

14. Phan TG, Luchsinger V, Avendaño LF, Deng X, Delwart E. Cyclovirus in nasopharyngeal aspirates of Chilean children with respiratory infections. J Gen Virol. (2014) 95:922–7. doi: 10.1099/vir.0.061143-0

PubMed Abstract | Crossref Full Text | Google Scholar

15. Li L, Kapoor A, Slikas B, Bamidele OS, Wang C, Shaukat S, et al. Multiple diverse circoviruses infect farm animals and are commonly found in human and chimpanzee feces. J Virol. (2010) 84:1674–82. doi: 10.1128/JVI.02109-09

PubMed Abstract | Crossref Full Text | Google Scholar

16. Opriessnig T, Xiao CT, Mueller NJ, Denner J. Emergence of novel circoviruses in humans and pigs and their possible importance for xenotransplantation and blood transfusions. Xenotransplantation. (2024) 31:e12842. doi: 10.1111/xen.12842

PubMed Abstract | Crossref Full Text | Google Scholar

17. Gainor K, Becker A, Malik YS, Ghosh S. Detection and complete genome analysis of circoviruses and cycloviruses in the small Indian mongoose (Urva auropunctata): identification of novel species. Viruses. (2021) 13:1700. doi: 10.3390/v13091700

PubMed Abstract | Crossref Full Text | Google Scholar

18. Shan T, Yang S, Wang H, Wang H, Zhang J, Gong G, et al. Virome in the cloaca of wild and breeding birds revealed a diversity of significant viruses. Microbiome. (2022) 10:60. doi: 10.1186/s40168-022-01246-7

PubMed Abstract | Crossref Full Text | Google Scholar

19. Lu X, Ji L, Wang H, Zhang Q, Wang X, Liu Y, et al. Highly diverse RNA viruses and phage sequences concealed within birds. Microbiol Spectr. (2024) 12:e0080224. doi: 10.1128/spectrum.00802-24

PubMed Abstract | Crossref Full Text | Google Scholar

20. Schäffer AA, Nawrocki EP, Choi Y, Kitts PA, Karsch-Mizrachi I, McVeigh R. VecScreen_plus_taxonomy: imposing a tax(onomy) increase on vector contamination screening. Bioinformatics. (2018) 34:755–9. doi: 10.1093/bioinformatics/btx669

PubMed Abstract | Crossref Full Text | Google Scholar

21. Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience. (2012) 1:18. doi: 10.1186/2047-217X-1-18

PubMed Abstract | Crossref Full Text | Google Scholar

22. Kapli P, Yang Z, Telford MJ. Phylogenetic tree building in the genomic age. Nat Rev Genet. (2020) 21:428–44. doi: 10.1038/s41576-020-0233-0

PubMed Abstract | Crossref Full Text | Google Scholar

23. Nath BK, Das S, Roby JA, Sarker S, Luque D, Raidal SR, et al. Structural perspectives of beak and feather disease virus and porcine circovirus proteins. Viral Immunol. (2021) 34:49–59. doi: 10.1089/vim.2020.0097

PubMed Abstract | Crossref Full Text | Google Scholar

24. Aiyar A. The use of CLUSTAL W and CLUSTAL X for multiple sequence alignment. Methods Mol Biol. (2000) 132:221–41. doi: 10.1385/1-59259-192-2:221

PubMed Abstract | Crossref Full Text | Google Scholar

25. French EL, Purchase HG, Nazerian K, A. new herpesvirus isolated from a nestling cormorant (Phalacrocorax melanoleucos). Avian Pathol. (1973) 2:3–15. doi: 10.1080/03079457308418059

PubMed Abstract | Crossref Full Text | Google Scholar

26. Fehér E, Kaszab E, Bali K, Hoitsy M, Sós E, Bányai K, et al. Novel gammapolyomavirus in a great cormorant (Phalacrocorax carbo). Arch Virol. (2022) 167:1721–4. doi: 10.1007/s00705-022-05478-8

PubMed Abstract | Crossref Full Text | Google Scholar

27. Khalifeh A, Sahnan S, Kale S, Olivo D, Lund MC, Kraberger S, et al. Genomes of cressdnaviruses, microviruses, a gyrovirus, and a caudovirus identified in the fecal samples of a mallard and double-crested cormorant. Microbiol Resour Announc. (2024) 13:e0033224. doi: 10.1128/mra.00332-24

PubMed Abstract | Crossref Full Text | Google Scholar

28. Liu Y, Zhang X, Han X, Liu J, Yao L. Development of a droplet digital PCR method for detection of porcine circovirus 4. BMC Vet Res. (2023) 19:129. doi: 10.1186/s12917-023-03690-5

PubMed Abstract | Crossref Full Text | Google Scholar

29. Alex CE, Fahsbender E, Altan E, Bildfell R, Wolff P, Jin L, et al. Viruses in unexplained encephalitis cases in American black bears (Ursus americanus). PLoS ONE. (2020) 15:e0244056. doi: 10.1371/journal.pone.0244056

PubMed Abstract | Crossref Full Text | Google Scholar

30. Yang Y, Xu T, Wen J, Yang L, Lai S, Sun X, et al. Prevalence and phylogenetic analysis of porcine circovirus type 2 (PCV2) and type 3 (PCV3) in the Southwest of China during 2020-2022. Front Vet Sci. (2022) 9:1042792. doi: 10.3389/fvets.2022.1042792

PubMed Abstract | Crossref Full Text | Google Scholar

31. Wang Y, Guo C, He J. Aquatic circoviruses: emerging pathogens in global aquaculture - from discovery to disease management. J Virol. (2025) 99:e0173724. doi: 10.1128/jvi.01737-24

PubMed Abstract | Crossref Full Text | Google Scholar

32. Wu S, Yip CC, Situ J, Li Z, Ho SS, Cai J, et al. Human circovirus in patients with hepatitis, Hong Kong. Emerg Infect Dis. (2024) 30:2521–31. doi: 10.3201/eid3012.241114

PubMed Abstract | Crossref Full Text | Google Scholar

33. Liu Y, Gong QL, Nie LB, Wang Q, Ge GY, Li DL, et al. Prevalence of porcine circovirus 2 throughout China in 2015-2019: a systematic review and meta-analysis. Microb Pathog. (2020) 149:104490. doi: 10.1016/j.micpath.2020.104490

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: viral metagenomics, circovirus, sequence alignment, phylogenetic analysis, genomic structure

Citation: Pei Y, Cai S, Xue Y, Fu Y, Zhang J, Shen Q, Ji L, Wu P, Wang H, Wang Y, Zhang W and Yang S (2025) First detection of two cycloviruses in cormorant fecal samples in China by high-throughput sequencing technology. Front. Vet. Sci. 12:1677378. doi: 10.3389/fvets.2025.1677378

Received: 31 July 2025; Accepted: 25 August 2025;
Published: 16 September 2025.

Edited by:

Mengmeng Zhao, Foshan University, China

Reviewed by:

Jun Luo, South China Agricultural University, China
Gen Li, South China Agricultural University, China
Han Gao, Foshan University, China

Copyright © 2025 Pei, Cai, Xue, Fu, Zhang, Shen, Ji, Wu, Wang, Wang, Zhang and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Shixing Yang, am9obnNvbnlhbmcxOTc5QDE2My5jb20=; Wen Zhang, ejAyMTZ3ZW5AeWFob28uY29t; Yan Wang, d2FuZ3lhbl9qdHVAMTI2LmNvbQ==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.