Simultaneous detection of pathogens and antimicrobial resistance genes with the open source, cloud-based, CZ ID pipeline.

Antimicrobial resistant (AMR) pathogens represent urgent threats to human health, and their surveillance is of paramount importance. Metagenomic next generation sequencing (mNGS) has revolutionized such efforts, but remains challenging due to the lack of open-access bioinformatics tools capable of simultaneously analyzing both microbial and AMR gene sequences. To address this need, we developed the CZ ID AMR module, an open-access, cloud-based workflow designed to integrate detection of both microbes and AMR genes in mNGS and whole-genome sequencing (WGS) data. It leverages the Comprehensive Antibiotic Resistance Database and associated Resistance Gene Identifier software, and works synergistically with the CZ ID short-read mNGS module to enable broad detection of both microbes and AMR genes. We highlight diverse applications of the AMR module through analysis of both publicly available and newly generated mNGS and WGS data from four clinical cohort studies and an environmental surveillance project. Through genomic investigations of bacterial sepsis and pneumonia cases, hospital outbreaks, and wastewater surveillance data, we gain a deeper understanding of infectious agents and their resistomes, highlighting the value of integrating microbial identification and AMR profiling for both research and public health. We leverage additional functionalities of the CZ ID mNGS platform to couple resistome profiling with the assessment of phylogenetic relationships between nosocomial pathogens, and further demonstrate the potential to capture the longitudinal dynamics of pathogen and AMR genes in hospital acquired bacterial infections. In sum, the new AMR module advances the capabilities of the open-access CZ ID microbial bioinformatics platform by integrating pathogen detection and AMR profiling from mNGS and WGS data. Its development represents a critical step toward democratizing pathogen genomic analysis and supporting collaborative efforts to combat the growing threat of AMR.

Simultaneous detection of pathogens and antimicrobial resistance genes with the open source, cloud-based, CZ ID pipeline.

Antimicrobial resistant (AMR) pathogens represent urgent threats to human health, and their surveillance is of paramount importance. Metagenomic next generation sequencing (mNGS) has revolutionized such efforts, but remains challenging due to the lack of open-access bioinformatics tools capable of simultaneously analyzing both microbial and AMR gene sequences. To address this need, we developed the Chan Zuckerberg ID (CZ ID) AMR module, an open-access, cloud-based workflow designed to integrate detection of both microbes and AMR genes in mNGS and whole-genome sequencing (WGS) data. It leverages the Comprehensive Antibiotic Resistance Database and associated Resistance Gene Identifier software, and works synergistically with the CZ ID short-read mNGS module to enable broad detection of both microbes and AMR genes. We highlight diverse applications of the AMR module through analysis of both publicly available and newly generated mNGS and WGS data from four clinical cohort studies and an environmental surveillance project. Through genomic investigations of bacterial sepsis and pneumonia cases, hospital outbreaks, and wastewater surveillance data, we gain a deeper understanding of infectious agents and their resistomes, highlighting the value of integrating microbial identification and AMR profiling for both research and public health. We leverage additional functionalities of the CZ ID mNGS platform to couple resistome profiling with the assessment of phylogenetic relationships between nosocomial pathogens, and further demonstrate the potential to capture the longitudinal dynamics of pathogen and AMR genes in hospital acquired bacterial infections. In sum, the new AMR module advances the capabilities of the open-access CZ ID microbial bioinformatics platform by integrating pathogen detection and AMR profiling from mNGS and WGS data. Its development represents a critical step toward democratizing pathogen genomic analysis and supporting collaborative efforts to combat the growing threat of AMR.

Metagenomic Pathogen Sequencing in Resource-Scarce Settings: Lessons Learned and the Road Ahead.

Metagenomic next-generation sequencing (mNGS) is the process of sequencing all genetic material in a biological sample. The technique is growing in popularity with myriad applications including outbreak investigation, biosurveillance, and pathogen detection in clinical samples. However, mNGS programs are costly to build and maintain, and additional obstacles faced by low- and middle-income countries (LMICs) may further widen global inequities in mNGS capacity. Over the past two decades, several important infectious disease outbreaks have highlighted the importance of establishing widespread sequencing capacity to support rapid disease detection and containment at the source. Using lessons learned from the COVID-19 pandemic, LMICs can leverage current momentum to design and build sustainable mNGS programs, which would form part of a global surveillance network crucial to the elimination of infectious diseases.

Viruses in unexplained encephalitis cases in American black bears (Ursus americanus).

Viral infections were investigated in American black bears (Ursus americanus) from Nevada and northern California with and without idiopathic encephalitis. Metagenomics analyses of tissue pools revealed novel viruses in the genera Circoviridae, Parvoviridae, Anelloviridae, Polyomaviridae, and Papillomaviridae. The circovirus and parvovirus were of particular interest due to their potential importance as pathogens. We characterized the genomes of these viruses and subsequently screened bears by PCR to determine their prevalence. The circovirus (Ursus americanus circovirus, UaCV) was detected at a high prevalence (10/16, 67%), and the chaphamaparvovirus (Ursus americanus parvovirus, UaPV) was found in a single bear. We showed that UaCV is present in liver, spleen/lymph node, and brain tissue of selected cases by in situ hybridization (ISH) and PCR. Infections were detected in cases of idiopathic encephalitis and in cases without inflammatory brain lesions. Infection status was not clearly correlated with disease, and the significance of these infections remains unclear. Given the known pathogenicity of a closely related mammalian circovirus, and the complex manifestations of circovirus-associated diseases, we suggest that UaCV warrants further study as a possible cause or contributor to disease in American black bears.

High prevalence, genetic diversity and a potentially novel genotype of Sapelovirus A (Picornaviridae) in enteric and respiratory samples in Hungarian swine farms.

All of the known porcine sapeloviruses (PSVs) currently belong to a single genotype in the genus Sapelovirus (family Picornaviridae). Here, the complete genome of a second, possibly recombinant, genotype of PSV strain SZ1M-F/PSV/HUN2013 (MN807752) from a faecal sample of a paraplegic pig in Hungary was characterized using viral metagenomics and RT-PCR. This sapelovirus strain showed only 64 % nucleotide identity in the VP1 region to its closest PSV-1 relative. Complete VP1 sequence-based epidemiological investigations of PSVs circulating in Hungary showed the presence of diverse strains found in high prevalence in enteric and respiratory samples collected from both asymptomatic and paraplegic pigs from 12 swine farms. Virus isolation attempts using PK-15 cell cultures were successful in 3/8 cases for the classic but not the novel PSV genotype. Sequence comparisons of faeces and isolate strains derived VP1 showed that cultured PSV strains not always represent the dominant PSVs found in vivo.

Plasma virome of 781 Brazilians with unexplained symptoms of arbovirus infection include a novel parvovirus and densovirus.

Plasma from patients with dengue-like symptoms was collected in 2013 to 2016 from the Brazilian states of Tocantins and Amapa. 781 samples testing negative for IgM against Dengue, Zika, and Chikungunya viruses and for flaviviruses, alphaviruses and enteroviruses RNA using RT-PCRs were analyzed using viral metagenomics. Viral particles-associated nucleic acids were enriched, randomly amplified, and deep sequenced in 102 mini-pools generating over 2 billion reads. Sequence data was analyzed for the presence of known and novel eukaryotic viral reads. Anelloviruses were detected in 80%, human pegivirus 1 in 19%, and parvovirus B19 in 17% of plasma pools. HIV and enteroviruses were detected in two pools each. Previously uncharacterized viral genomes were also identified, and their presence in single plasma samples confirmed by PCR. Chapparvovirus and ambidensovirus genomes, both in the Parvoviridae family, were partially characterized showing 33% and 34% identity in their NS1 sequences to their closest relative. Molecular surveillance using pre-existing plasma from febrile patients provides a readily scalable approach for the detection of novel, potentially emerging, viruses.

High Diversity and Novel Enteric Viruses in Fecal Viromes of Healthy Wild and Captive Thai Cynomolgus Macaques (Macaca fascicularis).

Cynomolgus macaques are common across South East Asian countries including Thailand. The National Primate Research Center of Thailand, Chulalongkorn University (NPRCT-CU) captures wild-borne cynomolgus macaque for research use. Limited information is available on the enteric viruses and possible zoonotic infections into or from cynomolgus macaques. We characterized and compare the fecal virome of two populations; healthy wild-originated captive cynomolgus macaques (n = 43) reared in NPRCT-CU and healthy wild cynomolgus macaques (n = 35). Over 90% of recognized viral sequence reads amplified from feces were from bacterial viruses. Viruses from seven families of mammalian viruses were also detected (Parvoviridae, Anelloviridae, Picornaviridae, Adenoviridae, Papillomaviridae, Herpesviridae, and Caliciviridae). The genomes of a member of a new picornavirus genus we named Mafapivirus, a primate chapparvovirus, and a circular Rep-encoding single-strand (CRESS) DNA virus were also characterized. Higher abundance of CRESS DNA viruses of unknown tropism and invertebrate-tropic ambidensovirus were detected in wild versus captive macaques likely reflecting dietary differences. Short term rearing in captivity did not have a pronounced effect on the diversity of mammalian viruses of wild cynomolgus macaques. This study is the first report of the fecal virome of cynomolgus macaques, non-human primates frequently used in biomedical research and vaccination studies.

Diverse picornaviruses are prevalent among free-living and laboratory rats (Rattus norvegicus) in Hungary and can cause disseminated infections.

In this study, the full length genomes of three phylogenetically distant picornaviruses (family Picornaviridae) belonging to the genus Rosavirus (rat08/rRoB/HUN, MN116648), Kobuvirus (rat08/rAiA/HUN, MN116647), and Cardiovirus (rat08/rCaB/HUN, MN116646) were obtained from a single faecal sample of a free-living Norway rat (Rattus norvegicus) in Hungary using viral metagenomics and RT-PCR/Sanger sequencing. The acquired complete genomes were in silico analyzed in detail revealing the presence of a second minor open reading frame encoding an alternative Leader peptide (L*) in rat08/rCaB/HUN and a ca. 222 nt-long sequence repeat with compact secondary RNA structure in the 3' UTR of rat08/rRoB/HUN. The studied rat picornaviruses were frequently detectable by RT-PCR with relatively high viral loads ranged between 8.99E+02 and 1.29E+06 copies/ml in rat faecal samples collected from five geographically distant locations throughout Hungary. The VP1 sequence-based phylogenetic analyses show the presence of multiple, mostly location-specific lineages for all three picornaviruses. Rat rosavirus and rat cardiovirus were identified in spleen while rat cardiovirus was also detected in liver, muscle and kidney samples with variable copy numbers (6.42E+01-1.90E+05 copies/µg total RNA) suggesting extra-intestinal dissemination. Both viruses were also prevalent (70.0% and 18.2%) among two populations of laboratory rats ("Wistar-type" and "hooded-type") held in different, isolated laboratory animal units.

Lyon-IARC Polyomavirus DNA in Feces of Diarrheic Cats.

A viral metagenomic analysis of feces from an unexplained outbreak of feline diarrhea revealed the presence of Lyon-IARC polyomavirus (LIPyV) DNA. LIPyV, whose genome was originally sequenced from swabs of human skin, was fecally shed by three out of five diarrheic cats.

Chapparvovirus DNA Found in 4% of Dogs with Diarrhea.

Feces from dogs in an unexplained outbreak of diarrhea were analyzed by viral metagenomics revealing the genome of a novel parvovirus. The parvovirus was named cachavirus and was classified within the proposed Chapparvovirus genus. Using PCR, cachavirus DNA was detected in two of nine tested dogs from that outbreak. In order to begin to elucidate the clinical impact of this virus, 2,053 canine fecal samples were screened using real-time PCR. Stool samples from 203 healthy dogs were positive for cachavirus DNA at a rate of 1.47%, while 802 diarrhea samples collected in 2017 and 964 samples collected in 2018 were positive at rates of 4.0% and 4.66% frequencies, respectively (healthy versus 2017-2018 combined diarrhea p-value of 0.05). None of 83 bloody diarrhea samples tested positive. Viral loads were generally low with average real-time PCR Ct values of 36 in all three positive groups. The species tropism and pathogenicity of cachavirus, the first chapparvovirus reported in feces of a placental carnivore, remains to be fully determined.

Genomes of ubiquitous marine and hypersaline Hydrogenovibrio, Thiomicrorhabdus and Thiomicrospira spp. encode a diversity of mechanisms to sustain chemolithoautotrophy in heterogeneous environments.

Chemolithoautotrophic bacteria from the genera Hydrogenovibrio, Thiomicrorhabdus and Thiomicrospira are common, sometimes dominant, isolates from sulfidic habitats including hydrothermal vents, soda and salt lakes and marine sediments. Their genome sequences confirm their membership in a deeply branching clade of the Gammaproteobacteria. Several adaptations to heterogeneous habitats are apparent. Their genomes include large numbers of genes for sensing and responding to their environment (EAL- and GGDEF-domain proteins and methyl-accepting chemotaxis proteins) despite their small sizes (2.1-3.1 Mbp). An array of sulfur-oxidizing complexes are encoded, likely to facilitate these organisms' use of multiple forms of reduced sulfur as electron donors. Hydrogenase genes are present in some taxa, including group 1d and 2b hydrogenases in Hydrogenovibrio marinus and H. thermophilus MA2-6, acquired via horizontal gene transfer. In addition to high-affinity cbb3 cytochrome c oxidase, some also encode cytochrome bd-type quinol oxidase or ba3 -type cytochrome c oxidase, which could facilitate growth under different oxygen tensions, or maintain redox balance. Carboxysome operons are present in most, with genes downstream encoding transporters from four evolutionarily distinct families, which may act with the carboxysomes to form CO2 concentrating mechanisms. These adaptations to habitat variability likely contribute to the cosmopolitan distribution of these organisms.

Genomic evolution, recombination, and inter-strain diversity of chelonid alphaherpesvirus 5 from Florida and Hawaii green sea turtles with fibropapillomatosis.

Chelonid alphaherpesvirus 5 (ChHV5) is a herpesvirus associated with fibropapillomatosis (FP) in sea turtles worldwide. Single-locus typing has previously shown differentiation between Atlantic and Pacific strains of this virus, with low variation within each geographic clade. However, a lack of multi-locus genomic sequence data hinders understanding of the rate and mechanisms of ChHV5 evolutionary divergence, as well as how these genomic changes may contribute to differences in disease manifestation. To assess genomic variation in ChHV5 among five Hawaii and three Florida green sea turtles, we used high-throughput short-read sequencing of long-range PCR products amplified from tumor tissue using primers designed from the single available ChHV5 reference genome from a Hawaii green sea turtle. This strategy recovered sequence data from both geographic regions for approximately 75% of the predicted ChHV5 coding sequences. The average nucleotide divergence between geographic populations was 1.5%; most of the substitutions were fixed differences between regions. Protein divergence was generally low (average 0.08%), and ranged between 0 and 5.3%. Several atypical genes originally identified and annotated in the reference genome were confirmed in ChHV5 genomes from both geographic locations. Unambiguous recombination events between geographic regions were identified, and clustering of private alleles suggests the prevalence of recombination in the evolutionary history of ChHV5. This study significantly increased the amount of sequence data available from ChHV5 strains, enabling informed selection of loci for future population genetic and natural history studies, and suggesting the (possibly latent) co-infection of individuals by well-differentiated geographic variants.

Diverse and highly recombinant anelloviruses associated with Weddell seals in Antarctica.

The viruses circulating among Antarctic wildlife remain largely unknown. In an effort to identify viruses associated with Weddell seals (Leptonychotes weddellii) inhabiting the Ross Sea, vaginal and nasal swabs, and faecal samples were collected between November 2014 and February 2015. In addition, a Weddell seal kidney and South Polar skua (Stercorarius maccormicki) faeces were opportunistically sampled. Using high throughput sequencing, we identified and recovered 152 anellovirus genomes that share 63-70% genome-wide identities with other pinniped anelloviruses. Genome-wide pairwise comparisons coupled with phylogenetic analysis revealed two novel anellovirus species, tentatively named torque teno Leptonychotes weddellii virus (TTLwV) -1 and -2. TTLwV-1 (n = 133, genomes encompassing 40 genotypes) is highly recombinant, whereas TTLwV-2 (n = 19, genomes encompassing three genotypes) is relatively less recombinant. This study documents ubiquitous TTLwVs among Weddell seals in Antarctica with frequent co-infection by multiple genotypes, however, the role these anelloviruses play in seal health remains unknown.

Discovery of a novel circular DNA virus in the Forbes sea star, Asterias forbesi.

A single-stranded DNA (ssDNA) virus, Asterias forbesi-associated circular virus (AfaCV), was discovered in a Forbes sea star displaying symptoms of sea star wasting disease (SSWD). The AfaCV genome organization is typical of circular Rep-encoding ssDNA (CRESS-DNA) viruses and is similar to that of members of the family Circoviridae. PCR-based surveys indicate that AfaCV is not clearly associated with SSWD, whereas the sea star-associated densovirus (SSaDV), recently implicated in SSWD in the Pacific, was prevalent in symptomatic specimens. AfaCV represents the first CRESS-DNA virus detected in echinoderms, adding to the growing diversity of these viruses recently recovered from invertebrates.

Development of a Serological Assay for the Sea Lion (Zalophus californianus) Anellovirus, ZcAV.

New diseases in marine animals are emerging at an increasing rate, yet methodological limitations hinder characterization of viral infections. Viral metagenomics is an effective method for identifying novel viruses in diseased animals; however, determining virus pathogenesis remains a challenge. A novel anellovirus (Zalophus californianus anellovirus, ZcAV) was recently reported in the lungs of captive California sea lions involved in a mortality event. ZcAV was not detected by PCR in the blood of these animals, creating the inability to assess the prevalence of ZcAV in live sea lions. This study developed an enzyme-linked immunosorbent assay (ELISA) to detect antibodies to ZcAV in sea lion serum. To assess ZcAV prevalence, paired serum and lung samples (n = 96) from wild sea lions that stranded along the California coast were tested through ELISA and PCR, respectively. Over 50% of the samples tested positive for ZcAV by ELISA (34%), PCR (29%), or both (11%) assays. ZcAV is prevalent in stranded wild sea lion populations and results suggest that PCR assays alone may grossly underestimate ZcAV exposure. This ELISA provides a tool for testing live sea lions for ZcAV exposure and is valuable for subsequent studies evaluating the potential pathogenicity of this anellovirus.

Densovirus associated with sea-star wasting disease and mass mortality.

Populations of at least 20 asteroid species on the Northeast Pacific Coast have recently experienced an extensive outbreak of sea-star (asteroid) wasting disease (SSWD). The disease leads to behavioral changes, lesions, loss of turgor, limb autotomy, and death characterized by rapid degradation ("melting"). Here, we present evidence from experimental challenge studies and field observations that link the mass mortalities to a densovirus (Parvoviridae). Virus-sized material (i.e., <0.2 µm) from symptomatic tissues that was inoculated into asymptomatic asteroids consistently resulted in SSWD signs whereas animals receiving heat-killed (i.e., control) virus-sized inoculum remained asymptomatic. Viral metagenomic investigations revealed the sea star-associated densovirus (SSaDV) as the most likely candidate virus associated with tissues from symptomatic asteroids. Quantification of SSaDV during transmission trials indicated that progression of SSWD paralleled increased SSaDV load. In field surveys, SSaDV loads were more abundant in symptomatic than in asymptomatic asteroids. SSaDV could be detected in plankton, sediments and in nonasteroid echinoderms, providing a possible mechanism for viral spread. SSaDV was detected in museum specimens of asteroids from 1942, suggesting that it has been present on the North American Pacific Coast for at least 72 y. SSaDV is therefore the most promising candidate disease agent responsible for asteroid mass mortality.

Divergent astrovirus associated with neurologic disease in cattle.

Using viral metagenomics of brain tissue from a young adult crossbreed steer with acute onset of neurologic disease, we sequenced the complete genome of a novel astrovirus (BoAstV-NeuroS1) that was phylogenetically related to an ovine astrovirus. In a retrospective analysis of 32 cases of bovine encephalitides of unknown etiology, 3 other infected animals were detected by using PCR and in situ hybridization for viral RNA. Viral RNA was restricted to the nervous system and detected in the cytoplasm of affected neurons within the spinal cord, brainstem, and cerebellum. Microscopically, the lesions were of widespread neuronal necrosis, microgliosis, and perivascular cuffing preferentially distributed in gray matter and most severe in the cerebellum and brainstem, with increasing intensity caudally down the spinal cord. These results suggest that infection with BoAstV-NeuroS1 is a potential cause of neurologic disease in cattle.