Highly diverse posaviruses in swine faeces are aquatic in origin

Diversity of Picorna-Like Viruses in the Teltow Canal, Berlin, Germany

The viromes of freshwater bodies are underexplored. The Picornavirales order, with 371 acknowledged species, is one of the most expansive and diverse groups of eukaryotic RNA viruses. In this study, we add 513 picorna-like viruses to the assemblage of more than 2000 unassigned picorna-like viruses. Our set of the aquatic Picornavirales virome of the Teltow Canal in Berlin, Germany, consists of 239 complete and 274 partial genomes. This urban freshwater body is characterized by the predominance of marna-like viruses (30.8%) and dicistro-like viruses (19.1%), whereas picornaviruses, iflaviruses, solinvi-like viruses, polycipi-like viruses, and nora-like viruses are considerably less prevalent. Caliciviruses and secoviruses were absent in our sample. Although presenting characteristic domains of Picornavirales, more than 100 viruses (20.8%) could not be assigned to any of the 9 Picornavirales families. Thirty-three viruses of the Marnaviridae-mostly locarna-like viruses-exhibit a monocistronic genome layout. Besides a wealth of novel virus sequences, viruses with peculiar features are reported. Among these is a clade of untypeable marna-like viruses with dicistronic genomes, but with the capsid protein-encoding open reading frame located at the 5' part of their RNA. A virus with a similar genome layout but clustering with dicistroviruses was also observed. We further detected monocistronic viruses with a polymerase gene related to aparaviruses. The detection of Aichi virus and five novel posa-like viruses indicates a slight burden in municipal wastewater.

Highly diverse ribonucleic acid viruses in the viromes of eukaryotic host species in Yunnan province, China

Background

The diversity in currently documented viruses and their morphological characteristics indicates the need for understanding the evolutionary characteristics of viruses. Notably, further studies are needed to obtain a comprehensive landscape of virome, the virome of host species in Yunnan province, China.

Materials and methods

We implemented the metagenomic next-generation sequencing strategy to investigate the viral diversity, which involved in 465 specimens collected from bats, pangolins, monkeys, and other species. The diverse RNA viruses were analyzed, especially focusing on the genome organization, genetic divergence and phylogenetic relationships.

Results

In this study, we investigated the viral composition of eight libraries from bats, pangolins, monkeys, and other species, and found several diverse RNA viruses, including the Alphacoronavirus from bat specimens. By characterizing the genome organization, genetic divergence, and phylogenetic relationships, we identified five Alphacoronavirus strains, which shared phylogenetic association with Bat-CoV-HKU8-related strains. The pestivirus-like virus related to recently identified Dongyang pangolin virus (DYPV) strains from dead pangolin specimens, suggesting that these viruses are evolving. Some genomes showed higher divergence from known species (e.g., calicivirus CS9-Cali-YN-CHN-2020), and many showed evidence of recombination events with unknown or known strains (e.g., mamastroviruses BF2-astro-YN-CHN-2020 and EV-A122 AKM5-YN-CHN-2020). The newly identified viruses showed extensive changes and could be assigned as new species, or even genus (e.g., calicivirus CS9-Cali-YN-CHN-2020 and iflavirus Ifla-YN-CHN-2020). Moreover, we identified several highly divergent RNA viruses and estimated their evolutionary characteristics among different hosts, providing data for further examination of their evolutionary dynamics.

Conclusion

Overall, our study emphasizes the close association between emerging viruses and infectious diseases, and the need for more comprehensive surveys.

Metagenomic Analysis of RNA Fraction Reveals the Diversity of Swine Oral Virome on South African Backyard Swine Farms in the uMgungundlovu District of KwaZulu-Natal Province

Numerous RNA viruses have been reported in backyard swine populations in various countries. In the absence of active disease surveillance, a persistent knowledge gap exists on the diversity of RNA viruses in South African backyard swine populations. This is the first study investigating the diversity of oral RNA virome of the backyard swine in South Africa. We used three samples of backyard swine oral secretion (saliva) collected from three distantly located backyard swine farms (BSFs) in the uMgungundlovu District, KwaZulu-Natal, South Africa. Total viral RNA was extracted and used for the library preparation for deep sequencing using the Illumina HiSeq X instrument. The FASTQ files containing paired-end reads were analyzed using Genome Detective v 1.135. The assembled nucleotide sequences were analyzed using the PhyML phylogenetic tree. The genome sequence analysis identified a high diversity of swine enteric viruses in the saliva samples obtained from BSF2 and BSF3, while only a few viruses were identified in the saliva obtained from BSF1. The swine enteric viruses belonged to various animal virus families; however, two fungal viruses, four plant viruses, and five unclassified RNA viruses were also identified. Specifically, viruses of the family Astroviridae, according to the number of reads, were the most prevalent. Of note, the genome sequences of Rotavirus A (RVA) and Rotavirus C (RVC) at BSF2 and RVC and Hepatitis E virus (HEV) at BSF3 were also obtained. The occurrence of various swine enteric viruses in swine saliva suggests a high risk of diarrhoeic diseases in the backyard swine. Of note, zoonotic viruses in swine saliva, such as RVA, RVC, and HEV, indicate a risk of zoonotic spillover to the exposed human populations. We recommend the implementation of biosecurity to ensure sustainable backyard swine farming while safeguarding public health.

Application of shotgun metagenomics sequencing and targeted sequence capture to detect circulating porcine viruses in the Dutch-German border region

Porcine viruses have been emerging in recent decades, threatening animal and human health, as well as economic stability for pig farmers worldwide. Next‐generation sequencing (NGS) can detect and characterize known and unknown viruses but has limited sensitivity when an unbiased approach, such as shotgun metagenomics sequencing, is used. To increase the sensitivity of NGS for the detection of viruses, we applied and evaluated a broad viral targeted sequence capture (TSC) panel and compared it to an unbiased shotgun metagenomic approach. A cohort of 36 pooled porcine nasal swab and blood serum samples collected from both sides of the Dutch–German border region were evaluated. Overall, we detected 46 different viral species using TSC, compared to 40 viral species with a shotgun metagenomics approach. Furthermore, we performed phylogenetic analysis on recovered influenza A virus (FLUAV) genomes from Germany and revealed a close similarity to a zoonotic influenza strain previously detected in the Netherlands. Although TSC introduced coverage bias within the detected viruses, it improved sensitivity, genome sequence depth and contig length. In‐depth characterization of the swine virome, coupled with developing new enrichment techniques, can play a crucial role in the surveillance of circulating porcine viruses and emerging zoonotic pathogens.

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

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

Multiple clades of Husavirus in South America revealed by next generation sequencing

Husavirus (HuV) is an unclassified virus of the order Picornavirales that has already been identified worldwide in various locations. The genetic, epidemiological, and pathogenic characteristics are, however, little understood. In children with acute gastroenteritis, this study used next-generation sequencing to recognize unknown sources of viruses. In particular, 251 fecal samples obtained from individuals were sequenced in southern, northeastern, and northern Brazil. all samples were also analyzed using culture methods and parasitological tests to classify other enteric pathogens such as bacteria, parasites, and viruses. 1.9% of the samples tested positive for HuV, for a total of 5 positive children, with a mean age of 2 year, with three males and two females. Detailed molecular characterization of full genomes showed that Brazilian HuVs’ nucleotide divergence is less than 11%. The genetic gap between Brazilian sequences and the closest HuV reported previously, on the other hand, is 18%. The study showed that Brazilian sequences are closely related to the HuV defined in Viet Nam in 2013, further characterization based on phylogenetics. At least two divergent clades of HuV in South America were also seen in the phylogenetic study.

The Husavirus Posa-Like Viruses in China, and a New Group of Picornavirales

Novel posa-like viral genomes were first identified in swine fecal samples using metagenomics and were designated as unclassified viruses in the order Picornavirales. In the present study, nine husavirus strains were identified in China. Their genomes share 94.1–99.9% similarity, and alignment of these nine husavirus strains identified 697 nucleotide polymorphism sites across their full-length genomes. These nine strains were directly clustered with the Husavirus 1 lineage, and their genomic arrangement showed similar characteristics. These posa-like viruses have undergone a complex evolutionary process, and have a wide geographic distribution, complex host spectrum, deep phylogenetic divergence, and diverse genomic organizations. The clade of posa-like viruses forms a single group, which is evolutionarily distinct from other known families and could represent a distinct family within the Picornavirales. The genomic arrangement of Picornavirales and the new posa-like viruses are different, whereas the posa-like viruses have genomic modules similar to the families Dicistroviridae and Marnaviridae. The present study provides valuable genetic evidence of husaviruses in China, and clarifies the phylogenetic dynamics and the evolutionary characteristics of Picornavirales.

Increased viral read counts and metagenomic full genome characterization of porcine astrovirus 4 and Posavirus 1 in sows in a swine farm with unexplained neonatal piglet diarrhea

Neonatal diarrhea in piglets may cause major losses in affected pig herds. The present study used random high-throughput RNA sequencing (metagenomic next generation sequencing, mNGS) to investigate the virome of sows from a farm with persistent neonatal piglet diarrhea in comparison to two control farms without diarrhea problems. A variety of known swine gastrointestinal viruses was detected in the control farms as well as in the problem farm (Mamastrovirus, Enterovirus, Picobirnavirus, Posavirus 1, Kobuvirus, Proprismacovirus). A substantial increase in normalized viral read counts was observed in the affected farm compared to the control farms. The increase was attributable to a single viral species in each of the sampled sows (porcine astrovirus 4 and Posavirus 1). The complete genomes of a porcine astrovirus 4 and two co-infecting Posavirus 1 were de novo assembled and characterized. The 6734 nt single-stranded RNA genome of porcine astrovirus 4 (PoAstV-4) strain Belgium/2019 contains three overlapping open reading frames (nonstructural protein 1ab, nonstructural protein 1a, capsid protein). Posavirus 1 strains Belgium/01/2019 and Belgium/02/2019 have a 9814 nt single-stranded positive-sense RNA genome encoding a single open reading frame (polyprotein precursor) containing the five expected Picornavirales-conserved protein domains. The study highlights the potential of mNGS workflows to study unexplained neonatal diarrhea in piglets and contributes to the scarce availability of both PoAstV-4 and Posavirus-1 whole genome sequences from Western Europe.

Molecular-based investigation and genetic characterization of porcine stool-associated RNA virus (posavirus) lineages 1 to 3 in pigs in South Korea from 2017 to 2019

Recent results on the detection and genetic characterization of stool-associated RNA viruses from different species have increased the knowledge about the extreme genetic diversity of picornaviruses. This study aimed to investigate the presence of unclassified porcine stool-associated RNA viruses (posaviruses) in South Korea and to elucidate the molecular evolution of the viruses. By RT-PCR, posaviruses 1 and 3 were exclusively found in fecal samples and consistently detected in three consecutive years in six of eight provinces, with 148/697 (21.2%) and 33/84 (39.3%) positive samples and farms, respectively. Every age group but the older age groups (finisher, sow) had significantly higher positive rates of posavirus 1 than posavirus 3. An analysis of the RNA-dependent RNA polymerase sequences by likelihood mapping and maximum-likelihood-based phylogenetic analysis revealed that stool-associated RNA viruses formed four supergroups that were well separated from all recognized families of the order Picornavirales. Five genomes of Korean posaviruses generated in this study were phylogenetically grouped with posavirus 1 and posavirus 3 and were predicted to have the typical genome organization of picornaviruses.

Phylogenetic analysis of novel posaviruses detected in feces of Japanese pigs with posaviruses and posa-like viruses of vertebrates and invertebrates

Posaviruses and posa-like viruses are unclassified viruses with sequence similarity to viruses of the order Picornavirales. They have been reported in various vertebrates and invertebrates. We identified 11 posavirus-like sequences in porcine feces and performed phylogenic analysis. Previously reported Japanese posaviruses and those identified in this study clustered with posavirus 1, 4, and 7 and husavirus 1, while five viruses branched into three independent lineages, tentatively named posavirus 10, 11, and 12. Interestingly, posaviruses, except for posavirus 8 and 9, husaviruses, and rasaviruses, formed a cluster consisting of viruses only from pigs, humans, and rats, while posavirus 8 and 9, fisavirus, and basaviruses clustered with posa-like viruses from invertebrates.

First demonstration of the circulation of a pneumovirus in French pigs by detection of anti-swine orthopneumovirus nucleoprotein antibodies

The presence of pneumoviruses in pigs is poorly documented. In this study, we used the published sequence of the nucleoprotein (N) of the recently identified Swine Orthopneumovirus (SOV) to express and purify SOV N as a recombinant protein in Escherichia coli. This protein was purified as nanorings and used to set up an enzyme-linked immunosorbent assay, which was used to analyse the presence of anti-pneumovirus N antibodies in swine sera. Sera collected from different pig farms in the West of France and from specific pathogen free piglets before colostrum uptake showed indirectly that a pneumovirus is circulating in pig populations with some variations between animals. Piglets before colostrum uptake were sero-negative for anti-pneumovirus antibodies while most of the other pigs showed positivity. Interestingly, in two farms presenting respiratory clinical signs and negative or under control for some common respiratory pathogens, pigs were detected positive for anti-pneumovirus antibodies. Globally, anti-pneumovirus N antibody concentrations were variable between and within farms. Further studies will aim to isolate the circulating virus and determine its potential pathogenicity. SOV could potentially become a new member of the porcine respiratory complex, important on its own or in association with other viral and bacterial micro-organisms.

Origins and Evolution of the Global RNA Virome

Viruses with RNA genomes dominate the eukaryotic virome, reaching enormous diversity in animals and plants. The recent advances of metaviromics prompted us to perform a detailed phylogenomic reconstruction of the evolution of the dramatically expanded global RNA virome. The only universal gene among RNA viruses is the gene encoding the RNA-dependent RNA polymerase (RdRp). We developed an iterative computational procedure that alternates the RdRp phylogenetic tree construction with refinement of the underlying multiple-sequence alignments. The resulting tree encompasses 4,617 RNA virus RdRps and consists of 5 major branches; 2 of the branches include positive-sense RNA viruses, 1 is a mix of positive-sense (+) RNA and double-stranded RNA (dsRNA) viruses, and 2 consist of dsRNA and negative-sense (-) RNA viruses, respectively. This tree topology implies that dsRNA viruses evolved from +RNA viruses on at least two independent occasions, whereas -RNA viruses evolved from dsRNA viruses. Reconstruction of RNA virus evolution using the RdRp tree as the scaffold suggests that the last common ancestors of the major branches of +RNA viruses encoded only the RdRp and a single jelly-roll capsid protein. Subsequent evolution involved independent capture of additional genes, in particular, those encoding distinct RNA helicases, enabling replication of larger RNA genomes and facilitating virus genome expression and virus-host interactions. Phylogenomic analysis reveals extensive gene module exchange among diverse viruses and horizontal virus transfer between distantly related hosts. Although the network of evolutionary relationships within the RNA virome is bound to further expand, the present results call for a thorough reevaluation of the RNA virus taxonomy.IMPORTANCE The majority of the diverse viruses infecting eukaryotes have RNA genomes, including numerous human, animal, and plant pathogens. Recent advances of metagenomics have led to the discovery of many new groups of RNA viruses in a wide range of hosts. These findings enable a far more complete reconstruction of the evolution of RNA viruses than was attainable previously. This reconstruction reveals the relationships between different Baltimore classes of viruses and indicates extensive transfer of viruses between distantly related hosts, such as plants and animals. These results call for a major revision of the existing taxonomy of RNA viruses.

Complex virome in feces from Amerindian children in isolated Amazonian villages

The number of viruses circulating in small isolated human populations may be reduced by viral extinctions and rare introductions. Here we used viral metagenomics to characterize the eukaryotic virome in feces from healthy children from a large urban center and from three Amerindian villages with minimal outside contact. Numerous human enteric viruses, mainly from the Picornaviridae and Caliciviridae families, were sequenced from each of the sites. Multiple children from the same villages shed closely related viruses reflecting frequent transmission clusters. Feces of isolated villagers also contained multiple viral genomes of unknown cellular origin from the Picornavirales order and CRESS-DNA group and higher levels of nematode and protozoan DNA. Despite cultural and geographic isolation, the diversity of enteric human viruses was therefore not reduced in these Amazonian villages. Frequent viral introductions and/or increased susceptibility to enteric infections may account for the complex fecal virome of Amerindian children in isolated villages.

Enteric virome of Ethiopian children participating in a clean water intervention trial

BACKGROUND

The enteric viruses shed by different populations can be influenced by multiple factors including access to clean drinking water. We describe here the eukaryotic viral genomes in the feces of Ethiopian children participating in a clean water intervention trial.

METHODOLOGY/PRINCIPAL FINDINGS

Fecal samples from 269 children with a mean age of 2.7 years were collected from 14 villages in the Amhara region of Ethiopia, half of which received a new hand-dug water well. Feces from these villages were then analyzed in 29 sample pools using viral metagenomics. A total of 127 different viruses belonging to 3 RNA and 3 DNA viral families were detected. Picornaviridae family sequence reads were the most commonly found, originating from 14 enterovirus and 6 parechovirus genotypes plus multiple members of four other picornavirus genera (cosaviruses, saliviruses, kobuviruses, and hepatoviruses). Picornaviruses with nearly identical capsid VP1 were detected in different pools reflecting recent spread of these viral strains. Next in read frequencies and positive pools were sequences from the Caliciviridae family including noroviruses GI and GII and sapoviruses. DNA viruses from multiple genera of the Parvoviridae family were detected (bocaviruses 1-4, bufavirus 3, and dependoparvoviruses), together with four species of adenoviruses and common anelloviruses shedding. RNA in the order Picornavirales and CRESS-DNA viral genomes, possibly originating from intestinal parasites or dietary sources, were also characterized. No significant difference was observed between the number of mammalian viruses shed from children from villages with and without a new water well.

CONCLUSIONS

We describe an approach to estimate the efficacy of potentially virus transmission-reducing interventions and the first complete (DNA and RNA viruses) description of the enteric viromes of East African children. A wide diversity of human enteric viruses was found in both intervention and control groups. Mammalian enteric virome diversity was not reduced in children from villages with a new water well. This population-based sampling also provides a baseline of the enteric viruses present in Northern Ethiopia against which to compare future viromes.

Detection and complete genome characteristics of Posavirus 1 from pigs in China

Porcine stool-associated RNA virus 1 (Posavirus 1) is a novel member of picornaviruses and first identified from fecal samples of 30-day-old pigs with diarrhea in USA in 2011. To evaluate the existence of Posavirus 1 in swine herds, 118 clinical samples from diarrheal pigs and 31 fecal swabs from healthy pigs were collected and detected by reverse transcription-polymerase chain reaction (RT-PCR) using Posavirus 1-specific primers. Only five fecal samples from diarrheal pigs on two swine farms were positive for Posavirus 1. The complete genome sequences [excluding poly (A) tail] of two representative isolates SDQD-25 and HBTS-11 are determined and consist of 9840 and 9819 nucleotides in length, and encode one putative polyprotein of 3070 and 2952 amino acids, respectively. They share 90.3% homology with each other and 81.3-95.4% homologies with American Posavirus 1 isolates or strains at the nucleotide sequence level. The phylogenetic analysis based on the entire genomes of reference picornavirus strains or isolates showed SDQD-25, HBTS-11 cluster together with American Posavirus 1 isolates or strains, yet are clearly distant from the other picornaviruses. The complete genome sequences of Chinese Posavirus 1 isolates will enrich the information of Posavirus 1 sequence database and further expedite posavirus research on the genetic diversity, epidemiology, and evolution in China.

Characterization of Posa and Posa-like virus genomes in fecal samples from humans, pigs, rats, and bats collected from a single location in Vietnam

Porcine stool-associated RNA virus (posavirus), and Human stool-associated RNA virus (husavirus) are viruses in the order Picornavirales recently described in porcine and human fecal samples. The tentative group (Posa and Posa-like viruses: PPLVs) also includes fish stool-associated RNA virus (fisavirus) as well as members detected in insects (Drosophila subobscura and Anopheles sinensis) and parasites (Ascaris suum). As part of an agnostic deep sequencing survey of animal and human viruses in Vietnam, we detected three husaviruses in human fecal samples, two of which share 97–98% amino acid identity to Dutch husavirus strains and one highly divergent husavirus with only 25% amino acid identity to known husaviruses. In addition, the current study found forty-seven complete posavirus genomes from pigs, ten novel rat stool-associated RNA virus genomes (tentatively named rasavirus), and sixteen novel bat stool-associated RNA virus genomes (tentatively named basavirus). The five expected Picornavirales protein domains (helicase, 3C-protease, RNA-dependent RNA polymerase, and two Picornavirus capsid domain) were found to be encoded by all PPLV genomes. In addition, a nucleotide composition analysis revealed that the PPLVs shared compositional properties with arthropod viruses and predicted non-mammalian hosts for all PPLV lineages. The study adds seventy-six genomes to the twenty-nine PPLV genomes currently available and greatly extends our sequence knowledge of this group of viruses within the Picornavirales order.

Virome comparisons in wild-diseased and healthy captive giant pandas

BACKGROUND

The giant panda (Ailuropoda melanoleuca) is a vulnerable mammal herbivore living wild in central China. Viral infections have become a potential threat to the health of these endangered animals, but limited information related to these infections is available.

METHODS

Using a viral metagenomic approach, we surveyed viruses in the feces, nasopharyngeal secretions, blood, and different tissues from a wild giant panda that died from an unknown disease, a healthy wild giant panda, and 46 healthy captive animals.

RESULTS

The previously uncharacterized complete or near complete genomes of four viruses from three genera in Papillomaviridae family, six viruses in a proposed new Picornaviridae genus (Aimelvirus), two unclassified viruses related to posaviruses in Picornavirales order, 19 anelloviruses in four different clades of Anelloviridae family, four putative circoviruses, and 15 viruses belonging to the recently described Genomoviridae family were sequenced. Reflecting the diet of giant pandas, numerous insect virus sequences related to the families Iflaviridae, Dicistroviridae, Iridoviridae, Baculoviridae, Polydnaviridae, and subfamily Densovirinae and plant viruses sequences related to the families Tombusviridae, Partitiviridae, Secoviridae, Geminiviridae, Luteoviridae, Virgaviridae, and Rhabdoviridae; genus Umbravirus, Alphaflexiviridae, and Phycodnaviridae were also detected in fecal samples. A small number of insect virus sequences were also detected in the nasopharyngeal secretions of healthy giant pandas and lung tissues from the dead wild giant panda. Although the viral families present in the sick giant panda were also detected in the healthy ones, a higher proportion of papillomaviruses, picornaviruses, and anelloviruses reads were detected in the diseased panda.

CONCLUSION

This viral survey increases our understanding of eukaryotic viruses in giant pandas and provides a baseline for comparison to viruses detected in future infectious disease outbreaks. The similar viral families detected in sick and healthy giant pandas indicate that these viruses result in commensal infections in most immuno-competent animals.

Identification of further diversity among posaviruses

Recently, there have been reports of new members of posavirus-like viruses in the order Picornavirales. In this study, using a metagenomics approach, 11 posavirus-like sequences (>7,000 nucleotides) were detected in 155 porcine fecal samples. Phylogenetic analysis revealed that the newly identified virus sequences, together with other posavirus-like viruses, form distinct clusters within the order Picornavirales, composed of eight genogroups and unassigned sequences based on amino acid sequences of the helicase and RNA-dependent RNA polymerase regions, with <40 % and <50 % sequence identity, respectively. We propose further classifications of highly diverse posavirus populations based on newly identified sequences from Japanese pig feces.