Microbial Diversity in Children with Gastroenteritis in the Amazon Region of Brazil: Development and Validation of a Molecular Method for Complete Sequencing of Viral Genomes.

INTRODUCTION: Metagenomic sequencing is a powerful tool that is widely used in laboratories worldwide for taxonomic characterization of microorganisms in clinical and environmental samples. In this study, we utilized metagenomics to investigate comprehensively the microbial diversity in fecal samples of children over a four-year period. Our methods were carefully designed to ensure accurate and reliable results. MATERIAL AND METHODS: Validated and analyzed were metagenomic data obtained from sequencing 27 fecal samples from children under 10 years old with gastroenteritis over a four-year period (2012-2016). The fecal specimens were collected from patients who received care at public health facilities in the northern region of Brazil. Sequencing libraries were prepared from cDNA and sequenced on the Illumina HiSeq. Kraken-2 was utilized to classify bacterial taxonomy based on the 16S rRNA gene, using the Silva rRNA database. Additionally, the Diamond program was used for mapping to the non-redundant protein database (NR database). Phylogenomic analyses were conducted using Geneious R10 and MEGA X software, and Bayesian estimation of phylogeny was performed using the MrBayes program. The results indicate significant heterogeneity among norovirus strains, with evidence of recombination and point mutations. This study presents the first complete genome of parechovirus 8 in the region. Additionally, it describes the bacterial populations and bacteriophages present in feces, with a high abundance of Firmicutes and Proteobacteria, including an increased proportion of the Enterobacteriaceae family. The presented data demonstrate the genetic diversity of microbial populations and provide a comprehensive report on viral molecular characterization. These findings are relevant for genomic studies in gastrointestinal infections. The metagenomic approach is a powerful tool for investigating microbial diversity in children with gastroenteritis. However, further studies are imperative to conduct genomic analysis of identified bacterial strains and thoroughly analyze antimicrobial resistance genes.

Diversity of picornaviruses detected in diarrheal samples from children in Belém, Brazilian Amazon (1982-2019).

In this investigation, fecal specimens from children with diarrhea were collected from four community studies conducted between 1982 and 2019 in Belém, Brazilian Amazon. A total of 234 samples were tested by quantitative reverse transcription polymerase chain reaction (RT-qPCR) to detect infections by picornaviruses of the Enterovirus (EV), Parechovirus (HPeV), Cosavirus (HCoSV), Kobuvirus (Aichivirus - AiV) and Salivirus (SalV) genera. The positive samples were subjected to different amplification protocols of the VP1 region of the genome, such as nested PCR or snPCR, and were subsequently genotyped by sequencing VP1 and VP3 of the viral genome. Positivity was observed in 76.5% (179/234) of the samples tested using RT-qPCR for at least one virus, and co-infection was observed in 37.4% (67/179) of the cases. EV was detected in 50.8% (119/234), HPeV in 29.9% (70/234), HCoSV in 27.3% (64/234), and AiV/SalV in 2.1% (5/234) of the specimens tested by RT-qPCR. Using nested PCR and/or snPCR techniques, the positivity rates were 94.11% (112/119) for EV, 72.85% (51/70) for HPeV, and 20.31% (13/64) for HCoSV. It was not possible to amplify the samples that were positive for AiV/SalV. Sequencing revealed 67.2% (80/119) EV, 51.4% (36/70) HPeV, and 20.31% (13/64) HCoSV. Forty-five different types of EV were found among species A, B, and C; HCoSV identified five species, including a possible recombinant strain; all HPeV were identified as belonging to species A, in two samples a possible recombination involving three different strains was verified. This study demonstrated the high circulation and diversity of different types of picornaviruses in fecal samples, including those collected more than 30 years ago. This endorsed the evaluation of important points in the epidemiology of these viruses, such as the presence of co-infection and the possibility of knowing more about these agents, considering that some were recently described; therefore, their detection in older samples can provide more data about their ancestry.

Epidemiologic and Clinical Characteristics of Human Bocavirus Infection in Children with or without Acute Gastroenteritis in Acre, Northern Brazil.

Human bocavirus (HBoV) is an emerging virus detected around the world that may be associated with cases of acute gastroenteritis (AGE). However, its contribution to AGE has not been elucidated. This study aimed to describe the frequency, clinical features, and HBoV species circulation in children up to 5 years with or without AGE symptoms in Acre, Northern Brazil. A total of 480 stool samples were collected between January and December 2012. Fecal samples were used for extraction, nested PCR amplification, and sequencing for genotyping. Statistical analysis was applied to verify the association between epidemiological and clinical characteristics. Overall, HBoV-positivity was 10% (48/480), with HBoV-positive rates of 8.4% (19/226) and 11.4% (29/254) recorded in diarrheic and non-diarrheic children, respectively. The most affected children were in the age group ranging between 7 and 24 months (50%). HBoV infection was more frequent in children who live in urban areas (85.4%), use water from public networks (56.2%), and live with adequate sewage facilities (50%). Co-detection with other enteric viruses was 16.7% (8/48) and the most prevalent coinfection was RVA+ HBoV (50%, 4/8). HBoV-1 was the most frequent species detected in diarrheic and non-diarrheic children, responsible for 43.8% (21/48) of cases, followed by HBoV-3 (29.2%, 14/48) and HBoV-2 (25%, 12/48). In this study, HBoV infection was not always associated with AGE, as most HBoV cases belonged to the non-diarrheal group. Future studies are warranted in order to determine the role of HBoV in causing acute diarrhea disease.

Detection of G3 human-like rotavirus in institutionalized dogs from Brazil.

Viral gastroenteritis is a common clinical problem in dogs and group A rotavirus (RVA) is one of the agents involved in this etiology. It mainly affects dogs in the first 6 months of life, and these animals are considered an important reservoir and potential transmitters of the virus to other susceptible hosts, such as humans. Among the different types of RVA, G3 is the most detected in dogs, and this genotype is also involved in infections in other animals, including humans. Thus, the present study aims to investigate the presence of RVA in samples of dogs from a public kennel. A total of 64 fecal samples from dogs with diarrhea were analyzed, collected from April 2019 to March 2020, from the kennel of the Zoonosis Control Center, located in Belém, a city in the North of Brazil. The extracted genetic material was subjected to reverse transcription followed by real-time PCR (RT-qPCR); the positives were tested by RT-PCR with a specific primer for the RVA VP7 gene, after nucleotide sequencing and phylogenetic analysis. One sample was subjected to high-performance sequencing. A positivity of 7.8% (5/64) was observed for RVA, all characterized as G3, grouping in the G3-III lineage, with greater similarity to human samples. Different regions of the RVA genome fragments were found. These results emphasize the need for animal health surveillance to better understand the global strain dispersion of RVA and elucidate possible interspecies transmission events, monitoring the genetic diversity of this pathogen.

Genome sequencing of dengue virus serotype 4 in a bat brain sample (Platyrrhinus helleri) from the Brazilian Amazon.

The existence of sylvatic transmission of dengue virus in communities of neotropical bats remains uncertain. In this work we present a near-complete genome of dengue virus serotype 4 obtained from the brain sample of a bat from Platyrrhinus helleri specie collected in the Brazilian Amazon region. The presence of the virus in the brain sample may indicate a possible tropism for the central nervous system in bats, which may justify negative results in previous studies that focused on analysis of other tissues, such as liver and spleen. Besides the duration of dengue virus circulation in the Americas (circa 40 years) may be too short for an implementation of a sylvatic dengue virus cycle. Our findings suggest that continued monitoring is needed to confirm with the neotropical bats could potentially act as a natural reservoir of dengue in the region.

Paramyxoviruses from neotropical bats suggest a novel genus and nephrotropism.

Paramyxoviruses have a broad host range and geographic distribution, including human pathogens transmitted by bats, such as Nipah and Hendra viruses. In this study, we combined high-throughput sequencing and molecular approaches to investigate the presence of paramyxoviruses in neotropical bats (Microchiroptera suborder) in Brazil. We discovered and characterized three novel paramyxoviruses in the kidney tissues of apparently healthy common vampire bats (D. rotundus) and Seba's short-tailed bats (C. perspicillata), which we tentatively named Kanhgág virus (KANV), Boe virus (BOEV), and Guató virus (GUATV). In this study, we classified these viruses as putative species into the Macrojêvirus genus, a newly proposed genus of the Orthoparamyxovirinae subfamily. Using RT-PCR, we detected these viruses in 20.9% (9 out of 43) of bats tested, and viral RNA was detected exclusively in kidney tissues. Attempts to isolate infectious virus were successful for KANV and GUATV. Our results expand the viral diversity, host range, and geographical distribution of the paramyxoviruses.

Differential expression analysis and profiling of hepatic miRNA and isomiRNA in dengue hemorrhagic fever.

Dengue virus causes dengue hemorrhagic fever (DHF) and has been associated to fatal cases worldwide. The liver is one of the most important target tissues in severe cases, due to its intense viral replication and metabolic role. microRNAs role during infection is crucial to understand the regulatory mechanisms of DENV infection and can help in diagnostic and anti-viral therapies development. We sequenced the miRNome of six fatal cases and compared to five controls, to characterize the human microRNAs expression profile in the liver tissue during DHF. Eight microRNAs were differentially expressed, including miR-126-5p, a regulatory molecule of endothelial cells, miR-122-5p, a liver specific homeostasis regulator, and miR-146a-5p, an interferon-regulator. Enrichment analysis with predicted target genes of microRNAs revealed regulatory pathways of apoptosis, involving MAPK, RAS, CDK and FAS. Immune response pathways were related to NF- kB, CC and CX families, IL and TLR. This is the first description of the human microRNA and isomicroRNA profile in liver tissues from DHF cases. The results demonstrated the association of miR-126-5p, miR-122-5p and miR-146a-5p with DHF liver pathogenesis, involving endothelial repair and vascular permeability regulation, control of homeostasis and expression of inflammatory cytokines.

Evaluation of drug susceptibility profile of Mycobacterium tuberculosis Lineage 1 from Brazil based on whole genome sequencing and phenotypic methods.

BACKGROUND: The evaluation of procedures for drug susceptibility prediction of Mycobacterium tuberculosis based on genomic data against the conventional reference method test based on culture is realistic considering the scenario of growing number of tools proposals based on whole-genome sequences (WGS). OBJECTIVES: This study aimed to evaluate drug susceptibility testing (DST) outcome based on WGS tools and the phenotypic methods performed on isolates of M. tuberculosis Lineage 1 from the state of Pará, Brazil, generally associated with low levels of drug resistance. METHODOLOGY: Culture based DST was performed using the Proportion Method in Löwenstein-Jensen medium on 71 isolates that had been submitted to WGS. We analysed the seven main genome sequence-based tools for resistance and lineage prediction applied to M. tuberculosis and for comparison evaluation we have used the Kappa concordance test. FINDINGS: When comparing the WGS-based tools against the DST, we observed the highest level of agreement using TB-profiler. Among the tools, TB-profiler, KvarQ and Mykrobe were those which identified the largest number of TB-MDR cases. Comparing the four most sensitive tools regarding resistance prediction, agreement was observed for 43 genomes. MAIN CONCLUSIONS: Drug resistance profiling using next-generation sequencing offers rapid assessment of resistance-associated mutations, therefore facilitating rapid access to effective treatment.

Near-complete genome of cosavirus A from a child hospitalized with acute gastroenteritis, Brazil.

Acute gastroenteritis (AG) is responsible for 525,000 deaths worldwide in children under-5-years and is caused by the Human Cosavirus (HCoSV; family Picornaviridae, Genus Cosavirus). Although its health importance, a significant percentage of diarrhea cases (≈ 40 %) still of unknown etiology. In Brazil, few studies have reported HCoSV-A sequences analyzing partial 5' UTR. This study characterized the first near-complete genome of a Cosavirus A (strain AM326) from a child hospitalized with AG in Amazonas state, Northern Brazil. High throughput sequencing (HTS) was performed using the HiSeq™ 2500 platform (Illumina) in one fecal specimen collected from the Surveillance of Rotavirus Network of the Evandro Chagas Institute collected in 2017. Sequence reads were assembled by the De Novo approach using three distinct algorithmic (IDBA-UD, Spades, and MegaHit). The final contig was recovered from the HCoSV-AM326 sample revealing 7,735 nt in length (SRA number SRR12535029; GenBank MT023104) and the genetic characterization, as well as phylogenetic analysis demonstrated a new variant strain from Brazil, highlighting the association of HCoSV-A as a possible causative agent of AG. This finding demonstrates the importance of the metagenomic approach to elucidate cases of diarrhea without a defined etiology, as well as providing a better understanding about the virus genetics, evolution and epidemiology.

Evolutionary and Molecular Analysis of Complete Genome Sequences of Norovirus From Brazil: Emerging Recombinant Strain GII.P16/GII.4.

Noroviruses (NoVs) are enteric viruses that cause acute gastroenteritis, and the pandemic GII.4 genotype is spreading and evolving rapidly. The recombinant GII.P16/GII.4_Sydney strain emerged in 2016, replacing GII.P31/GII.4_Sydney (GII.P31 formerly known as GII.Pe) in some countries. We analyzed the complete genome of 20 NoV strains (17 GII.P31/GII.4_ Sydney and 3 GII.P16/GII.4_Sydney) from Belém and Manaus, Brazil, collected from 2012 to 2016. Phylogenetic trees were constructed by maximum likelihood method from 191 full NoV-VP1 sequences, demonstrated segregation of the Sydney lineage in two larger clades, suggesting that GII.4 strains associated with GII.P16 already have modifications compared with GII.P31/GII.4. Additionally, the Bayesian Markov Chain Monte Carlo method was used to reconstruct a time-scaled phylogenetic tree formed by GII.P16 ORF1 sequences (n = 117) and three complete GII.P16 sequences from Belém. The phylogenetic tree indicated the presence of six clades classified into different capsid genotypes and locations. Evolutionary rates of the ORF1 gene of GII.P16 strains was estimated at 2.01 × 10-3 substitutions/site/year, and the most recent common ancestors were estimated in 2011 (2011-2012, 95% HPD). Comparing the amino acid (AA) sequence coding for ORF1 with the prototype strain GII.P16/GII.4, 36 AA changes were observed, mainly in the non-structural proteins p48, p22, and RdRp. GII.P16/GII.4 strains of this study presented changes in amino acids 310, 333, 373, and 393 of the antigenic sites in the P2 subdomain, and ML tree indicating the division within the Sydney lineage according to the GII.P16 and GII.P31 polymerases. Notably, as noroviruses have high recombination rates and the GII.4 genotype was prevalent for a long time in several locations, additional and continuous evolutionary analyses of this new genotype should be needed in the future.

High Genomic Variability in Equine Infectious Anemia Virus Obtained from Naturally Infected Horses in Pantanal, Brazil: An Endemic Region Case.

Equine infectious anemia virus (EIAV) is a persistent lentivirus that causes equine infectious anemia (EIA). In Brazil, EIAV is endemic in the Pantanal region, and euthanasia is not mandatory in this area. All of the complete genomic sequences from field viruses are from North America, Asia, and Europe, and only proviral genomic sequences are available. Sequences from Brazilian EIAV are currently available only for gag and LTR regions. Thus, the present study aimed for the first time to sequence the entire EIAV genomic RNA in naturally infected horses from an endemic area in Brazil. RNA in plasma from naturally infected horses was used for next-generation sequencing (NGS), and gaps were filled using Sanger sequencing methodology. Complete viral genomes of EIAV from two horses were obtained and annotated (Access Number: MN560970 and MN560971). Putative genes were analyzed and compared with previously described genes, showing conservation in gag and pol genes and high variations in LTR and env sequences. Amino acid changes were identified in the p26 protein, one of the most common targets used for diagnosis, and p26 molecular modelling showed surface amino acid alterations in some epitopes. Brazilian genome sequences presented 88.6% nucleotide identity with one another and 75.8 to 77.3% with main field strains, such as EIAV Liaoning, Wyoming, Ireland, and Italy isolates. Furthermore, phylogenetic analysis suggested that this Brazilian strain comprises a separate monophyletic group. These results may help to better characterize EIAV and to overcome the challenges of diagnosing and controlling EIA in endemic regions.

Evaluation of drug susceptibility profile of Mycobacterium tuberculosis Lineage 1 from Brazil based on whole genome sequencing and phenotypic methods

BACKGROUND The evaluation of procedures for drug susceptibility prediction of Mycobacterium tuberculosis based on genomic data against the conventional reference method test based on culture is realistic considering the scenario of growing number of tools proposals based on whole-genome sequences (WGS). OBJECTIVES This study aimed to evaluate drug susceptibility testing (DST) outcome based on WGS tools and the phenotypic methods performed on isolates of M. tuberculosis Lineage 1 from the state of Pará, Brazil, generally associated with low levels of drug resistance. METHODOLOGY Culture based DST was performed using the Proportion Method in Löwenstein-Jensen medium on 71 isolates that had been submitted to WGS. We analysed the seven main genome sequence-based tools for resistance and lineage prediction applied to M. tuberculosis and for comparison evaluation we have used the Kappa concordance test. FINDINGS When comparing the WGS-based tools against the DST, we observed the highest level of agreement using TB-profiler. Among the tools, TB-profiler, KvarQ and Mykrobe were those which identified the largest number of TB-MDR cases. Comparing the four most sensitive tools regarding resistance prediction, agreement was observed for 43 genomes. MAIN CONCLUSIONS Drug resistance profiling using next-generation sequencing offers rapid assessment of resistance-associated mutations, therefore facilitating rapid access to effective treatment.

Genomic characterization of orthobunyavirus of veterinary importance in America.

During 2013, in Argentina, three new isolates of serogroup Bunyamwera virus (genus Orthobunyavirus, family Peribunyaviridae) were recovered from two horses with encephalitis, and from an aborted equine fetus. In the present study, we report the complete genome sequence, genetic characterization, and phylogenetic analysis of three new strains isolated in Argentina to clarifying their relationship within the Bunyamwera serogroup virus and to investigate the evolutionary history of viruses with segmented genomes.

Comparative analysis of the transcriptome of the Amazonian fish species Colossoma macropomum (tambaqui) and hybrid tambacu by next generation sequencing.

BACKGROUND: The C. macropomum is a characiform fish from the Amazon basin that has been hybridized with other pacu species to produce commercial hybrids, such as the tambacu. However, little is known of the functional genomics of the parental species or these hybrid forms. The transcriptome of C. macropomum and tambacu were sequenced using 454 Roche platform (pyrosequencing) techniques to characterize the domains of Gene Ontology (GO) and to evaluate the levels of gene expression in the two organisms. RESULTS: The 8,188,945 reads were assembled into 400,845 contigs. A total of 58,322 contigs were annotated with a predominance of biological processes for both organisms, as determined by Gene Ontology (GO). Similar numbers of metabolic pathways were identified in both the C. macropomum and the tambacu, with the metabolism category presenting the largest number of transcripts. The BUSCO analysis indicated that our assembly was more than 40% complete. We identified 21,986 genes for the two fishes. The P and Log2FC values indicated significant differences in the levels of gene expression, with a total of 600 up-regulated genes. CONCLUSION: In spite of the lack of a reference genome, the functional annotation was successful, and confirmed a considerable difference in the specificity and levels of gene expression between the two organisms. This report provides a comprehensive baseline for the genetic management of these commercially important fishes, in particular for the identification of specific genes that may represent markers involved in the immunity, growth, and fertility of these organisms, with potential practical applications in aquaculture management.

Metagenomic Analysis of Samples from Three Bat Species Collected in the Amazon Rain Forest.

We report here the sequencing of five microbiome samples collected from different bat species in the Amazon rain forest. All contigs matching virus sequences were assigned to members of the Retroviridae family, while the bacterial contigs matched several bacterial species mostly belonging to the Proteobacteria phylum.

Complete Endogenous Retrovirus Genome Sequence from a Brazilian Vampire Bat (Desmodus rotundus).

The strain Desmodus rotundus endogenous retrovirus (DrERV) QR09 was obtained from a bat tissue sample collected from Desmodus rotundus in the Brazilian rain forest. The complete genome was sequenced using the next-generation sequencing strategy. The full-length genome of DrERV QR09 is 8,256 nucleotides in length and showed high similarity with other DrERVs.

First isolation of West Nile virus in Brazil.

BACKGROUND: Serological evidence of West Nile virus (WNV) infection has been reported in different regions of Brazil from equine and human hosts but the virus had never been isolated in the country. OBJECTIVES: We sought to identify the viral etiology of equine encephalitis in Espírito Santo state. METHODS: We performed viral culture in C6/36 cells, molecular detection of WNV genome, histopathology and immunohistochemistry from horse cerebral tissue. We also carried out sequencing, phylogenetic analysis and molecular clock. FINDINGS: Histopathologic analysis from horse cerebral tissue showed injury related to encephalitis and WNV infection was confirmed by immunohistochemistry. The virus was detected by reverse transcription quantitative polymerase chain reaction (RT-qPCR) from brain tissue and subsequently isolated in C6/36 cells. WNV full-length genome was sequenced showing the isolated strain belongs to lineage 1a. The molecular clock indicated that Brazilian WNV strain share the same common ancestor that were circulating in US during 2002-2005. MAIN CONCLUSIONS: Here we report the first isolation of WNV in Brazil from a horse with neurologic disease, which was clustered into lineage 1a with others US WNV strains isolated in beginning of 2000's decade.

Description and phylogeny of the mitochondrial genome of Sabethes chloropterus, Sabethes glaucodaemon and Sabethes belisarioi (Diptera: Culicidae).

Mosquitoes as Sabethes chloropterus, Sabethes glaucodaemon, Sabethes belisarioi are species of medical and epidemiological importance for arboviruses transmission such as yellow fever and St. Louis encephalitis. Despite this, no information about these three species mitochondrial DNA has been found in literature. Our study presents a mitochondrial genome description, including identity, SNPs, mutation rate, and phylogeny analysis using COX1, COX2, NADH4, NADH5, CYOB genes. The Sa. chloropterus, Sa. glaucodaemon and Sa. belisaroi mitochondrial genome sizes 15.609 bp, 15.620 bp, 15.907 bp, respectively, with 37 functional genes, presenting about 4.982 single nucleotide polymorphisms and 13.291 identical sites between them, besides all genes with dN/dS < 1 ratio, and also a greater approximation between Sa. glaucodaemon and Sa. chloropterus than with Sa. belisarioi. Due to the importance of mitochondrial DNA for population structure studies, evolution, and others, we expect that this data can contribute to other studies related to these mosquitoes and their viruses.

Characterization of Triniti virus supports its reclassification in the family Peribunyaviridae.

Triniti virus (TNTV) has been isolated in Trinidad and Tobago and in Brazil. To date little is known about this virus, which is classified as an ungrouped virus within the family Togaviridae. Here, three isolates of TNTV were characterized both genetically and antigenically. The genome was shown to contain three RNA segments: small (S), medium (M) and large (L). Genome organization, protein sizes and protein motifs were similar to those of viruses in the genus Orthobunyavirus, family Peribunyaviridae. Antigenic reactivity revealed the three TNTV isolates to be closely related, but no serologic cross-reaction with other orthobunyaviruses. Morphological observation by transmission electron microscopy indicated that virus size and symmetry were compatible with those of viruses in the family Peribunyaviridae. Our serological, morphological and molecular results support the taxonomic reclassification of TNTV as a member of the genus Orthobunyavirus, family Peribunyaviridae.

Oropouche orthobunyavirus: Genetic characterization of full-length genomes and development of molecular methods to discriminate natural reassortments.

Oropouche orthobunyavirus (OROV) has significant impact in public health in Amazon region. This arbovirus is one of the most common causes of febrile illness in Brazil, and is responsible for several epidemics since 1960's. In this study, we sequenced and characterized the complete coding sequences (S-, M- and L-RNA) of 35 OROV isolates from Brazil. Here, we classified 20 strains in genotype I from Pará and Maranhão states, nine as genotype II from Pará and Rondônia states confirmed, four classified into genotype III from Acre, Maranhão, Minas Gerais and Rondônia states and two genotype IV from Amazonas State. Also, we did not observe reassortment events involving the OROV isolates. In addition, we developed novel RT-PCR tools to identify reassortment events among OROV strains. These data will be useful to better understand the molecular epidemiology and diagnostic of OROV infections.