Metagenomics characterization of respiratory viral RNA pathogens in children under five years with severe acute respiratory infection in the Free State, South Africa.

Prompt detection of viral respiratory pathogens is crucial in managing respiratory infection including severe acute respiratory infection (SARI). Metagenomics next-generation sequencing (mNGS) and bioinformatics analyses remain reliable strategies for diagnostic and surveillance purposes. This study evaluated the diagnostic utility of mNGS using multiple analysis tools compared with multiplex real-time PCR for the detection of viral respiratory pathogens in children under 5 years with SARI. Nasopharyngeal swabs collected in viral transport media from 84 children admitted with SARI as per the World Health Organization definition between December 2020 and August 2021 in the Free State Province, South Africa, were used in this study. The obtained specimens were subjected to mNGS using the Illumina MiSeq system, and bioinformatics analysis was performed using three web-based analysis tools; Genome Detective, One Codex and Twist Respiratory Viral Research Panel. With average reads of 211323, mNGS detected viral pathogens in 82 (97.6%) of the 84 patients. Viral aetiologies were established in nine previously undetected/missed cases with an additional bacterial aetiology (Neisseria meningitidis) detected in one patient. Furthermore, mNGS enabled the much needed viral genotypic and subtype differentiation and provided significant information on bacterial co-infection despite enrichment for RNA viruses. Sequences of nonhuman viruses, bacteriophages, and endogenous retrovirus K113 (constituting the respiratory virome) were also uncovered. Notably, mNGS had lower detectability rate for severe acute respiratory syndrome coronavirus 2 (missing 18/32 cases). This study suggests that mNGS, combined with multiple/improved bioinformatics tools, is practically feasible for increased viral and bacterial pathogen detection in SARI, especially in cases where no aetiological agent could be identified by available traditional methods.

Metagenomics of gut microbiome for migratory seagulls in Kunming city revealed the potential public risk to human health.

BACKGROUND: Seagull as a migratory wild bird has become most popular species in southwest China since 1980s. Previously, we analyzed the gut microbiota and intestinal pathogenic bacteria configuration for this species by using 16S rRNA sequencing and culture methods. To continue in-depth research on the gut microbiome of migratory seagulls, the metagenomics, DNA virome and RNA virome were both investigated for their gut microbial communities of abundance and diversity in this study. RESULTS: The metagenomics results showed 99.72% of total species was bacteria, followed by viruses, fungi, archaea and eukaryota. In particular, Shigella sonnei, Escherichia albertii, Klebsiella pneumonia, Salmonella enterica and Shigella flexneri were the top distributed taxa at species level. PCoA, NMDS, and statistics indicated some drug resistant genes, such as adeL, evgS, tetA, PmrF, and evgA accumulated as time went by from November to January of the next year, and most of these genes were antibiotic efflux. DNA virome composition demonstrated that Caudovirales was the most abundance virus, followed by Cirlivirales, Geplafuvirales, Petitvirales and Piccovirales. Most of these phages corresponded to Enterobacteriaceae and Campylobacteriaceae bacterial hosts respectively. Caliciviridae, Coronaviridae and Picornaviridae were the top distributed RNA virome at family level of this migratory animal. Phylogenetic analysis indicated the sequences of contigs of Gammacoronavirus and Deltacoronavirus had highly similarity with some coronavirus references. CONCLUSIONS: In general, the characteristics of gut microbiome of migratory seagulls were closely related to human activities, and multiomics still revealed the potential public risk to human health.

Agnostic Sequencing for Detection of Viral Pathogens.

The advent of next-generation sequencing (NGS) technologies has expanded our ability to detect and analyze microbial genomes and has yielded novel molecular approaches for infectious disease diagnostics. While several targeted multiplex PCR and NGS-based assays have been widely used in public health settings in recent years, these targeted approaches are limited in that they still rely on a priori knowledge of a pathogen's genome, and an untargeted or unknown pathogen will not be detected. Recent public health crises have emphasized the need to prepare for a wide and rapid deployment of an agnostic diagnostic assay at the start of an outbreak to ensure an effective response to emerging viral pathogens. Metagenomic techniques can nonspecifically sequence all detectable nucleic acids in a sample and therefore do not rely on prior knowledge of a pathogen's genome. While this technology has been reviewed for bacterial diagnostics and adopted in research settings for the detection and characterization of viruses, viral metagenomics has yet to be widely deployed as a diagnostic tool in clinical laboratories. In this review, we highlight recent improvements to the performance of metagenomic viral sequencing, the current applications of metagenomic sequencing in clinical laboratories, as well as the challenges that impede the widespread adoption of this technology.

Metagenomic Detection of Divergent Insect- and Bat-Associated Viruses in Plasma from Two African Individuals Enrolled in Blood-Borne Surveillance.

Metagenomic next-generation sequencing (mNGS) has enabled the high-throughput multiplexed identification of sequences from microbes of potential medical relevance. This approach has become indispensable for viral pathogen discovery and broad-based surveillance of emerging or re-emerging pathogens. From 2015 to 2019, plasma was collected from 9586 individuals in Cameroon and the Democratic Republic of the Congo enrolled in a combined hepatitis virus and retrovirus surveillance program. A subset (n = 726) of the patient specimens was analyzed by mNGS to identify viral co-infections. While co-infections from known blood-borne viruses were detected, divergent sequences from nine poorly characterized or previously uncharacterized viruses were also identified in two individuals. These were assigned to the following groups by genomic and phylogenetic analyses: densovirus, nodavirus, jingmenvirus, bastrovirus, dicistrovirus, picornavirus, and cyclovirus. Although of unclear pathogenicity, these viruses were found circulating at high enough concentrations in plasma for genomes to be assembled and were most closely related to those previously associated with bird or bat excrement. Phylogenetic analyses and in silico host predictions suggested that these are invertebrate viruses likely transmitted through feces containing consumed insects or through contaminated shellfish. This study highlights the power of metagenomics and in silico host prediction in characterizing novel viral infections in susceptible individuals, including those who are immunocompromised from hepatitis viruses and retroviruses, or potentially exposed to zoonotic viruses from animal reservoir species.

RAPIDprep: A Simple, Fast Protocol for RNA Metagenomic Sequencing of Clinical Samples.

Emerging infectious disease threats require rapid response tools to inform diagnostics, treatment, and outbreak control. RNA-based metagenomics offers this; however, most approaches are time-consuming and laborious. Here, we present a simple and fast protocol, the RAPIDprep assay, with the aim of providing a cause-agnostic laboratory diagnosis of infection within 24 h of sample collection by sequencing ribosomal RNA-depleted total RNA. The method is based on the synthesis and amplification of double-stranded cDNA followed by short-read sequencing, with minimal handling and clean-up steps to improve processing time. The approach was optimized and applied to a range of clinical respiratory samples to demonstrate diagnostic and quantitative performance. Our results showed robust depletion of both human and microbial rRNA, and library amplification across different sample types, qualities, and extraction kits using a single workflow without input nucleic-acid quantification or quality assessment. Furthermore, we demonstrated the genomic yield of both known and undiagnosed pathogens with complete genomes recovered in most cases to inform molecular epidemiological investigations and vaccine design. The RAPIDprep assay is a simple and effective tool, and representative of an important shift toward the integration of modern genomic techniques with infectious disease investigations.

Bioinformatic Tools for NGS-Based Metagenomics to Improve the Clinical Diagnosis of Emerging, Re-Emerging and New Viruses.

Epidemics and pandemics have occurred since the beginning of time, resulting in millions of deaths. Many such disease outbreaks are caused by viruses. Some viruses, particularly RNA viruses, are characterized by their high genetic variability, and this can affect certain phenotypic features: tropism, antigenicity, and susceptibility to antiviral drugs, vaccines, and the host immune response. The best strategy to face the emergence of new infectious genomes is prompt identification. However, currently available diagnostic tests are often limited for detecting new agents. High-throughput next-generation sequencing technologies based on metagenomics may be the solution to detect new infectious genomes and properly diagnose certain diseases. Metagenomic techniques enable the identification and characterization of disease-causing agents, but they require a large amount of genetic material and involve complex bioinformatic analyses. A wide variety of analytical tools can be used in the quality control and pre-processing of metagenomic data, filtering of untargeted sequences, assembly and quality control of reads, and taxonomic profiling of sequences to identify new viruses and ones that have been sequenced and uploaded to dedicated databases. Although there have been huge advances in the field of metagenomics, there is still a lack of consensus about which of the various approaches should be used for specific data analysis tasks. In this review, we provide some background on the study of viral infections, describe the contribution of metagenomics to this field, and place special emphasis on the bioinformatic tools (with their capabilities and limitations) available for use in metagenomic analyses of viral pathogens.

A case for investment in clinical metagenomics in low-income and middle-income countries.

Clinical metagenomics is the diagnostic approach with the broadest capacity to detect both known and novel pathogens. Clinical metagenomics is costly to run and requires infrastructure, but the use of next-generation sequencing for SARS-CoV-2 molecular epidemiology in low-income and middle-income countries (LMICs) offers an opportunity to direct this infrastructure to the establishment of clinical metagenomics programmes. Local implementation of clinical metagenomics is important to create relevant systems and evaluate cost-effective methodologies for its use, as well as to ensure that reference databases and result interpretation tools are appropriate to local epidemiology. Rational implementation, based on the needs of LMICs and the available resources, could ultimately improve individual patient care in instances in which available diagnostics are inadequate and supplement emerging infectious disease surveillance systems to ensure the next pandemic pathogen is quickly identified.

Respiratory metagenomics: route to routine service.

PURPOSE OF REVIEW: The coronavirus disease 2019 pandemic demonstrated broad utility of pathogen sequencing with rapid methodological progress alongside global distribution of sequencing infrastructure. This review considers implications for now moving clinical metagenomics into routine service, with respiratory metagenomics as the exemplar use-case. RECENT FINDINGS: Respiratory metagenomic workflows have completed proof-of-concept, providing organism identification and many genotypic antimicrobial resistance determinants from clinical samples in <6 h. This enables rapid escalation or de-escalation of empiric therapy for patient benefit and reducing selection of antimicrobial resistance, with genomic-typing available in the same time-frame. Attention is now focussed on demonstrating clinical, health-economic, accreditation, and regulatory requirements. More fundamentally, pathogen sequencing challenges the traditional culture-orientated time frame of microbiology laboratories, which through automation and centralisation risks becoming increasingly separated from the clinical setting. It presents an alternative future where infection experts are brought together around a single genetic output in an acute timeframe, aligning the microbiology target operating model with the wider human genomic and digital strategy. SUMMARY: Pathogen sequencing is a transformational proposition for microbiology laboratories and their infectious diseases, infection control, and public health partners. Healthcare systems that link output from routine clinical metagenomic sequencing, with pandemic and antimicrobial resistance surveillance, will create valuable tools for protecting their population against future infectious diseases threats.

Viral metagenomics unveils MW (Malawi) polyomavirus infection in Brazilian pediatric patients with acute respiratory disease.

Viral metagenomics has been extensively applied for the identification of emerging or poorly characterized viruses. In this study, we applied metagenomics for the identification of viral infections among pediatric patients with acute respiratory disease, but who tested negative for SARS-CoV-2. Twelve pools composed of eight nasopharyngeal specimens were submitted to viral metagenomics. Surprisingly, in two of the pools, we identified reads belonging to the poorly characterized Malawi polyomavirus (MWPyV). Then, the samples composing the positive pools were individually tested using quantitative polymerase chain reaction for identification of the MWPyV index cases. MWPyV-positive samples were also submitted to respiratory virus panel testing due to the metagenomic identification of different clinically important viruses. Of note, MWPyV-positive samples tested also positive for respiratory syncytial virus types A and B. In this study, we retrieved two complete MWPyV genome sequences from the index samples that were submitted to phylogenetic inference to investigate their viral origin. Our study represents the first molecular and genomic characterization of MWPyV obtained from pediatric patients in South America. The detection of MWPyV in acutely infected infants suggests that this virus might participate (coparticipate) in cases of respiratory symptoms. Nevertheless, future studies based on testing of a larger number of clinical samples and MWPyV complete genomes appear to be necessary to elucidate if this emerging polyomavirus might be clinically important.

The Case for Studying New Viruses of New Hosts.

Virology has largely focused on viruses that are pathogenic to humans or to the other species that we care most about. There is no doubt that this has been a worthwhile investment. But many transformative advances have been made through the in-depth study of relatively obscure viruses that do not appear on lists of prioritized pathogens. In this review, I highlight the benefits that can accrue from the study of viruses and hosts off the beaten track. I take stock of viral sequence diversity across host taxa as an estimate of the bias that exists in our understanding of host-virus interactions. I describe the gains that have been made through the metagenomic discovery of thousands of new viruses in previously unsampled hosts as well as the limitations of metagenomic surveys. I conclude by suggesting that the study of viruses that naturally infect existing and emerging model organisms represents an opportunity to push virology forward in useful and hard to predict ways.

Genomic Analysis of the Suspicious SARS-CoV-2 Sequences in the Public Sequencing Database.

SARS-CoV-2 has infected more than 600 million people. However, the origin of the virus is still unclear; knowing where the virus came from could help us prevent future zoonotic epidemics. Sequencing data, particularly metagenomic data, can profile the genomes of all species in the sample, including those not recognized at the time, thus allowing for the identification of the progenitor of SARS-CoV-2 in samples collected before the pandemic. We analyzed the data from 5,196 SARS-CoV-2-positive sequencing runs in the NCBI's SRA database with collection dates prior to 2020 or unknown. We found that the mutation patterns obtained from these suspicious SARS-CoV-2 reads did not match the genome characteristics of an unknown progenitor of the virus, suggesting that they may derive from circulating SARS-CoV-2 variants or other coronaviruses. Despite a negative result for tracking the progenitor of SARS-CoV-2, the methods developed in the study could assist in pinpointing the origin of various pathogens in the future. IMPORTANCE Sequences that are homologous to the SARS-CoV-2 genome were found in numerous sequencing runs that were not associated with the SARS-CoV-2 studies in the public database. It is unclear whether they are derived from the possible progenitor of SARS-CoV-2 or contamination of more recent SARS-CoV-2 variants circulated in the population due to the lack of information on the collection, library preparation, and sequencing processes. We have developed a computational framework to infer the evolutionary relationship between sequences based on the comparison of mutations, which enabled us to rule out the possibility that these suspicious sequences originate from unknown progenitors of SARS-CoV-2.

Application of CRISPR-Based Human and Bacterial Ribosomal RNA Depletion for SARS-CoV-2 Shotgun Metagenomic Sequencing.

OBJECTIVES: The aim of this study is to evaluate the effectiveness of a CRISPR-based human and bacterial ribosomal RNA (rRNA) depletion kit (JUMPCODE Genomics) on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) shotgun metagenomic sequencing in weakly positive respiratory samples. METHODS: Shotgun metagenomics was performed on 40 respiratory specimens collected from solid organ transplant patients and deceased intensive care unit patients at UCLA Medical Center in late 2020 to early 2021. Human and bacterial rRNA depletion was performed on remnant library pools prior to sequencing by Illumina MiSeq. Data quality was analyzed using Geneious Prime, whereas the identification of SARS-CoV-2 variants and lineages was determined by Pangolin. RESULTS: The average genome coverage of the rRNA-depleted respiratory specimens increased from 72.55% to 93.71% in overall samples and from 29.3% to 83.3% in 15 samples that failed to achieve sufficient genome coverage using the standard method. Moreover, rRNA depletion enhanced genome coverage to over 85% in 11 (73.3%) of 15 low viral load samples with cycle threshold values up to 35, resulting in the identification of genotypes. CONCLUSION: The CRISPR-based human and bacterial rRNA depletion enhanced the sensitivity of SARS-CoV-2 shotgun metagenomic sequencing, especially in low viral load samples.

Mobilome-driven segregation of the resistome in biological wastewater treatment.

Biological wastewater treatment plants (BWWTP) are considered to be hotspots for the evolution and subsequent spread of antimicrobial resistance (AMR). Mobile genetic elements (MGEs) promote the mobilization and dissemination of antimicrobial resistance genes (ARGs) and are thereby critical mediators of AMR within the BWWTP microbial community. At present, it is unclear whether specific AMR categories are differentially disseminated via bacteriophages (phages) or plasmids. To understand the segregation of AMR in relation to MGEs, we analyzed meta-omic (metagenomic, metatranscriptomic and metaproteomic) data systematically collected over 1.5 years from a BWWTP. Our results showed a core group of 15 AMR categories which were found across all timepoints. Some of these AMR categories were disseminated exclusively (bacitracin) or primarily (aminoglycoside, MLS and sulfonamide) via plasmids or phages (fosfomycin and peptide), whereas others were disseminated equally by both. Combined and timepoint-specific analyses of gene, transcript and protein abundances further demonstrated that aminoglycoside, bacitracin and sulfonamide resistance genes were expressed more by plasmids, in contrast to fosfomycin and peptide AMR expression by phages, thereby validating our genomic findings. In the analyzed communities, the dominant taxon Candidatus Microthrix parvicella was a major contributor to several AMR categories whereby its plasmids primarily mediated aminoglycoside resistance. Importantly, we also found AMR associated with ESKAPEE pathogens within the BWWTP, and here MGEs also contributed differentially to the dissemination of the corresponding ARGs. Collectively our findings pave the way toward understanding the segmentation of AMR within MGEs, thereby shedding new light on resistome populations and their mediators, essential elements that are of immediate relevance to human health.

Multi-kingdom gut microbiota analyses define COVID-19 severity and post-acute COVID-19 syndrome.

Our knowledge of the role of the gut microbiome in acute coronavirus disease 2019 (COVID-19) and post-acute COVID-19 is rapidly increasing, whereas little is known regarding the contribution of multi-kingdom microbiota and host-microbial interactions to COVID-19 severity and consequences. Herein, we perform an integrated analysis using 296 fecal metagenomes, 79 fecal metabolomics, viral load in 1378 respiratory tract samples, and clinical features of 133 COVID-19 patients prospectively followed for up to 6 months. Metagenomic-based clustering identifies two robust ecological clusters (hereafter referred to as Clusters 1 and 2), of which Cluster 1 is significantly associated with severe COVID-19 and the development of post-acute COVID-19 syndrome. Significant differences between clusters could be explained by both multi-kingdom ecological drivers (bacteria, fungi, and viruses) and host factors with a good predictive value and an area under the curve (AUC) of 0.98. A model combining host and microbial factors could predict the duration of respiratory viral shedding with 82.1% accuracy (error ± 3 days). These results highlight the potential utility of host phenotype and multi-kingdom microbiota profiling as a prognostic tool for patients with COVID-19.

The gut virome in health and disease: new insights and associations.

PURPOSE OF REVIEW: Recent years have seen great strides made in the field of viral metagenomics. Many studies have reported alterations in the virome in different disease states. The vast majority of the human intestinal virome consists of bacteriophages, viruses that infect bacteria. The dynamic relationship between gut bacterial populations and bacteriophages is influenced by environmental factors that also impact host health and disease. In this review, we focus on studies highlighting the dynamics of the gut virome and fluctuations associated with disease states. RECENT FINDINGS: Novel correlations have been identified between the human gut virome and diseases such as obesity, necrotizing enterocolitis and severe acute respiratory syndrome coronavirus 2 infection. Further associations between the virome and cognition, diet and geography highlight the complexity of factors that can influence the dynamic relationship between gut bacteria, bacteriophages and health. SUMMARY: Here, we highlight some novel associations between the virome and health that will be the foundation for future studies in this field. The future development of microbiome-based interventions, identification of biomarkers, and novel therapeutics will require a thorough understanding of the gut virome and its dynamics.

Detection and discovery of plant viruses in soybean by metagenomic sequencing.

BACKGROUND: Viruses negatively impact soybean production by causing diseases that affect yield and seed quality. Newly emerging or re-emerging viruses can also threaten soybean production because current control measures may not be effective against them. Furthermore, detection and characterization of new plant viruses requires major efforts when no sequence or antibody-based resources are available. METHODS: In this study, soybean fields were scouted for virus-like disease symptoms during the 2016-2019 growing seasons. Total RNA was extracted from symptomatic soybean parts, cDNA libraries were prepared, and RNA sequencing was performed using high-throughput sequencing (HTS). A custom bioinformatic workflow was used to identify and assemble known and unknown virus genomes. RESULTS: Several viruses were identified in single or mixed infections. Full- or nearly full-length genomes were generated for tobacco streak virus (TSV), alfalfa mosaic virus (AMV), tobacco ringspot virus (TRSV), soybean dwarf virus (SbDV), bean pod mottle virus (BPMV), soybean vein necrosis virus (SVNV), clover yellow vein virus (ClYVV), and a novel virus named soybean ilarvirus 1 (SIlV1). Two distinct ClYVV isolates were recovered, and their biological properties were investigated in Nicotiana benthamiana, broad bean, and soybean. In addition to infections by individual viruses, we also found that mixed viral infections in various combinations were quite common. CONCLUSIONS: Taken together, the results of this study showed that HTS-based technology is a valuable diagnostic tool for the identification of several viruses in field-grown soybean and can provide rapid information about expected viruses as well as viruses that were previously not detected in soybean.

Deep Metagenomic Sequencing for Endophthalmitis Pathogen Detection Using a Nanopore Platform.

PURPOSE: To evaluate the utility of nanopore sequencing for identifying potential causative pathogens in endophthalmitis, comparing culture results against full-length 16S rRNA nanopore sequencing (16S Nanopore), whole genome nanopore sequencing (Nanopore WGS), and Illumina (Illumina WGS). DESIGN: Cross-sectional diagnostic comparison. METHODS: Patients with clinically suspected endophthalmitis underwent intraocular vitreous biopsy as per standard care. Clinical samples were cultured by conventional methods, together with full-length 16S rRNA and WGS using nanopore and Illumina sequencing platforms. RESULTS: Of 23 patients (median age 68.5 years [range 47-88]; 14 males [61%]), 18 cases were culture-positive. Nanopore sequencing identified the same cultured organism in all of the culture-positive cases and identified potential pathogens in two culture-negative cases (40%). Nanopore WGS was able to additionally detect the presence of bacteriophages in three samples. The agreements at genus level between culture and 16S Nanopore, Nanopore WGS, and Illumina WGS were 75%, 100%, and 78%, respectively. CONCLUSIONS: Whole genome sequencing has higher sensitivity and provides a viable alternative to culture and 16S sequencing for detecting potential pathogens in endophthalmitis. Moreover, WGS has the ability to detect other potential pathogens in culture-negative cases. Whilst Nanopore and Illumina WGS provide comparable data, nanopore sequencing provides potential for cost-effective point-of-care diagnostics.

Combination of Whole Genome Sequencing and Metagenomics for Microbiological Diagnostics.

Whole genome sequencing (WGS) provides the highest resolution for genome-based species identification and can provide insight into the antimicrobial resistance and virulence potential of a single microbiological isolate during the diagnostic process. In contrast, metagenomic sequencing allows the analysis of DNA segments from multiple microorganisms within a community, either using an amplicon- or shotgun-based approach. However, WGS and shotgun metagenomic data are rarely combined, although such an approach may generate additive or synergistic information, critical for, e.g., patient management, infection control, and pathogen surveillance. To produce a combined workflow with actionable outputs, we need to understand the pre-to-post analytical process of both technologies. This will require specific databases storing interlinked sequencing and metadata, and also involves customized bioinformatic analytical pipelines. This review article will provide an overview of the critical steps and potential clinical application of combining WGS and metagenomics together for microbiological diagnosis.

Pathogenic and metagenomic evaluations reveal the correlations of porcine epidemic diarrhea virus, porcine kobuvirus and porcine astroviruses with neonatal piglet diarrhea.

Porcine epidemic diarrhea virus (PEDV) frequently causes diarrhea outbreaks. However, whether newly discovered enteric viruses such as porcine kobuvirus (PKV) and porcine astroviruses (PAstVs) are also correlated with diarrhea is still unclear. Diarrhea outbreaks were reported in a PEDV-vaccinated pig farm in Xinjiang Uygur Autonomous Region of China from 2019 to 2020. PEDV was a common pathogen detected in fecal samples by routine RT-PCR assays. The PEDV positive fecal sample was used for pathogenic analysis due to the failure isolation of PEDV. The challenged neonatal piglets appeared watery diarrhea within one day post infection (dpi) and all died within 6 dpi. Histopathological and immunohistochemical examinations supported that PEDV is a major pathogen causing intestinal lesions. To further explore enteric viruses associated with neonatal piglet diarrhea, metagenomics sequencing was performed for the diarrheic piglets. Remarkably, PKV was the most abundant virus (58.33%) followed by PEDV (34.45%) and PAstVs (7.22%), which were also confirmed by real-time RT-PCR assays. Significant in vivo replications of PEDV and PKV could only be observed in challenged piglets whilst PAstVs maintained similar virus loads in both challenged and mock infected piglets. Overall, this study provides first pathogenic and metagenomic evidence that significant proliferations of PEDV and PKV are closely associated with severe diarrhea in neonatal piglets, while PAstVs likely play limited roles in neonatal piglet diarrhea.