Spatiotemporal transmission of SARS-CoV-2 lineages during 2020-2021 in Pernambuco-Brazil

SARS-CoV-2 virus emerged as a new threat to humans and spread around the world, leaving a large death toll. As of January 2023, Brazil is among the countries with the highest number of registered deaths. Nonpharmacological and pharmacological interventions have been heterogeneously implemented in the country, which, associated with large socioeconomic differences between the country regions, has led to distinct virus spread dynamics. Here, we investigate the spatiotemporal dispersion of SARS-CoV-2 lineages in the Pernambuco state (Northeast Brazil) throughout the distinct epidemiological scenarios that unfolded in the first 2 years of the pandemic. We generated a total of 1,389 new SARS-CoV-2 genomes from June 2020 to August 2021. This sampling captured the arrival, communitary transmission, and the circulation of the B1.1, B.1.1.28, and B.1.1.33 lineages; the emergence of the former variant of interest P.2; and the emergence and fast replacement of all previous variants by the more transmissible variant of concern P.1 (Gamma). Based on the incidence and lineage spread pattern, we observed an East-to-West to inner state pattern of transmission, which is in agreement with the transmission of more populous metropolitan areas to medium- and small-size country-side cities in the state. Such transmission patterns may be partially explained by the main routes of traffic across municipalities in the state. Our results highlight that the fine-grained intrastate analysis of lineages and incidence spread can provide actionable insights for planning future nonpharmacological intervention for air-borne transmissible human pathogens.IMPORTANCEDuring the COVID-19 pandemic, Brazil was one of the most affected countries, mainly due its continental-size, socioeconomic differences among regions, and heterogeneous implementation of intervention methods. In order to investigate SARS-CoV-2 dynamics in the state of Pernambuco, we conducted a spatiotemporal dispersion study, covering the period from June 2020 to August 2021, to comprehend the dynamics of viral transmission during the first 2 years of the pandemic. Throughout this study, we were able to track three significant epidemiological waves of transmission caused by B1.1, B.1.1.28, B.1.1.33, P.2, and P.1 lineages. These analyses provided valuable insights into the evolution of the epidemiological landscape, contributing to a deeper understanding of the dynamics of virus transmission during the early years of the pandemic in the state of Pernambuco.

In the Dawn of an Early Invasion: No Genetic Diversity of Angiostrongylus cantonensis in Ecuador?

The nematode Angiostrongylus cantonensis has been reported worldwide. However, some basic questions remain unanswered about A. cantonensis in Ecuador: (1) Was the invasion of A. cantonensis in Ecuador unique, or did it occur in different waves? (2) Was this invasion as recent as historical records suggest? (3) Did this invasion come from other regions of South America or elsewhere? To address these issues, we assessed the genetic diversity of MT-CO1 gene sequences from isolates obtained in 11 of Ecuador's 24 provinces. Our Bayesian inference phylogenetic tree recovered A. cantonensis as a well-supported monophyletic group. All 11 sequences from Ecuador were identical and identified as AC17a. The haplotype AC17a, found in Ecuador and the USA, formed a cluster with AC17b (USA), AC13 (Thailand), and AC12a-b (Cambodia). Notably, all the samples obtained in Ecuadorian provinces' different geographic and climatic regions had no genetic difference. Despite the lack of genetic information on A. cantonensis in Latin America, except in Brazil, our finding differs from previous studies by its absence of gene diversity in Ecuador. We concluded that the invasion of A. cantonensis in Ecuador may have occurred: (1) as a one-time event, (2) recently, and (3) from Asia via the USA. Further research should include samples from countries neighboring Ecuador to delve deeper into this.

SARS-CoV-2 mutations in Brazil: from genomics to putative clinical conditions

Due to the high rate of transmissibility, Brazil became the new COVID-19 outbreak epicenter and, since then, is being monitored to understand how SARS-CoV-2 mutates and spreads. We combined genomic and structural analysis to evaluate genomes isolated from different regions of Brazil and show that the most prevalent mutations were located in the S, N, ORF3a and ORF6 genes, which are involved in different stages of viral life cycle and its interaction with the host cells. Structural analysis brought to light the positions of these mutations on protein structures, contributing towards studies of selective structure-based drug discovery and vaccine development.

De novo transcriptome reveals blood coagulation/antithrombin factors and infection mechanisms in Angiostrongylus cantonensis adult worms

Angiostrongylus cantonensis is the main aetiological agent of eosinophilic meningoencephalitis in humans. Several outbreaks have been documented around the world, cementing its status as an emerging global public health concern. As a result, new strategies for the diagnosis, prophylaxis and treatment of cerebral angiostrongyliasis are urgently needed. In this study, we report on the de novo assembly of the A. cantonensis transcriptome, its full functional annotation and a reconstruction of complete metabolic pathways. All results are available at AngiostrongylusDB (http://angiostrongylus.lad.pucrs.br/admin/welcome). The aim of this study was to identify the active genes and metabolic pathways involved in the mechanisms of infection and survival inside Rattus norvegicus. Among 389 metabolic mapped pathways, the blood coagulation/antithrombin pathways of heparan sulphate/heparin are highlighted. Moreover, we identified genes codified to GP63 (leishmanolysin), CALR (calreticulin), ACE (peptidyl-dipeptidase A), myoglobin and vWD (von Willebrand factor type D domain protein) involved in the infection invasion and survival of the parasite. The large dataset of functional annotations provided and the full-length transcripts identified in this research may facilitate future functional genomics studies and provides a basis for the development of new techniques for the diagnosis, prevention and treatment of cerebral angiostrongyliasis.

In-depth transcriptome reveals the potential biotechnological application of Bothrops jararaca venom gland

Background:

Lack of complete genomic data of Bothrops jararaca impedes molecular biology research focusing on biotechnological applications of venom gland components. Identification of full-length coding regions of genes is crucial for the correct molecular cloning design.

Methods:

RNA was extracted from the venom gland of one adult female specimen of Bothrops jararaca. Deep sequencing of the mRNA library was performed using Illumina NextSeq 500 platform. De novo assembly of B. jararaca transcriptome was done using Trinity. Annotation was performed using Blast2GO. All predicted proteins after clustering step were blasted against non-redundant protein database of NCBI using BLASTP. Metabolic pathways present in the transcriptome were annotated using the KAAS-KEGG Automatic Annotation Server. Toxins were identified in the B. jararaca predicted proteome using BLASTP against all protein sequences obtained from Animal Toxin Annotation Project from Uniprot KB/Swiss-Pro database. Figures and data visualization were performed using ggplot2 package in R language environment.

Results:

We described the in-depth transcriptome analysis of B. jararaca venom gland, in which 76,765 de novo assembled isoforms, 96,044 transcribed genes and 41,196 unique proteins were identified. The most abundant transcript was the zinc metalloproteinase-disintegrin-like jararhagin. Moreover, we identified 78 distinct functional classes of proteins, including toxins, inhibitors and tumor suppressors. Other venom proteins identified were the hemolytic lethal factors stonustoxin and verrucotoxin.

Conclusion:

It is believed that the application of deep sequencing to the analysis of snake venom transcriptomes may represent invaluable insight on their biotechnological potential focusing on candidate molecules.

Identification of cross-reactive markers to strengthen the development of immunodiagnostic methods for angiostrongyliasis and other parasitic infections

Angiostrongylus cantonensis is the main causative agent of eosinophilic meningoencephalitis (EoM) in humans. Molecular diagnostic methods are essential since the identification of larvae in cerebrospinal fluid (CSF) is extremely rare. To date, the detection of a 31 kDa antigen by Western blotting has been the primary immunodiagnostic method for EoM caused by A. cantonensis. However, cross-reactivity with other parasites has been observed. Therefore, we conducted a comparative analysis using sera from individuals with angiostrongyliasis. We also characterized proteins isolated from different cellular sources of A. cantonensis, Toxocara canis, Schistosoma mansoni, and Strongyloides stercoralis with mass spectrometry. A total of 115 cross-reactive proteins were identified. Three of these proteins, heat shock protein, an intermediate filament protein, and galectin 1, represent potential markers for cross-reactivity. In addition, synthetic peptides were generated from previously identified diagnostic targets and tested against sera from individuals infected with several other parasites. As a result, two other markers of cross-reactivity were identified: peptide #4 derived from the 14-3-3 protein and peptide #12 derived from the Lec-5 protein. In contrast, 34 proteins were exclusively present in the Angiostrongylus extracts and represent promising diagnostic molecules for specific identification of A. cantonensis infection. In particular, cytochrome oxidase subunit I is of great interest as a possible immunodiagnostic target for angiostrongyliasis.

Mollusk microbiota shift during Angiostrongylus cantonensis infection in the freshwater snail Biomphalaria glabrata and the terrestrial slug Phillocaulis soleiformis

In the present work, we reported for the first time the microbiome from Phyllocaulis soleiformis and Biomphalaria glabrata assessed using high-throughput DNA sequencing pre- and post-infection with the helminth parasite Angiostrongylus cantonensis. B. glabrata and P. soleiformis were experimentally infected with A. cantonensis. Fecal DNAs from control and infected groups were extracted and subjected to 16S rRNA high-throughput sequencing survey. No significant differences were found in the alpha diversity indexes in Phyllocaulis and Biomphalaria experiments independently. PCoA analysis using the unweighted UniFrac measures showed that both microbiotas behaved differently depending on the host. In Biomphalaria microbiota, control and infected groups were significantly different (p = 0.0219), while Phyllocaulis samples were not (p = 0.5190). The microbiome of P. soleiformis infected with A. cantonensis showed a significant decrease of Sphingobacterium and a substantial increase of Cellvibrio when compared to a control group. The microbiome of B. glabrata infected with A. cantonensis showed a significant decline in the abundance of Flavobacterium, Fluviicola, Nitrospira, Vogesella and an OTU belonging to the family Comamonadaceae, and a significant increase of Uliginosibacterium and an OTU belonging to the family Weeksellaceae when compared to a control group. Overall, the microbiome data reported here provided valuable information with regard to the diversity of bacterial communities that comprise the gut microbiome of gastropods. Furthermore, we report here the effect of the infection of the helminth A. cantonensis in the ratio and distribution of the fecal microbiome of the snails. Further studies are highly valuable in order to better understand those interactions by comparing different microbiome profiles and mollusk models. By now, we anticipate that ecological studies will take significant advantage of these advances, particularly concerning improving our understanding of helminth-microbiome-host interactions.

Microbiota Profile and Impact of Fusobacterium nucleatum in Colorectal Cancer Patients of Barretos Cancer Hospital

Microbial diversity has been pointed out as a major factor in the development and progression of colorectal cancer (CRC). We sought to explore the richness and abundance of the microbial community of a series of colorectal tumor samples treated at Barretos Cancer Hospital, Brazil, through 16S rRNA sequencing. The presence and the impact of Fusobacterium nucleatum (Fn) DNA in CRC prognosis was further evaluated by qPCR in a series of 152 CRC cases. An enrichment for potentially oncogenic bacteria in CRC was observed, with Fusobacterium being the most abundant genus in the tumor tissue. In the validation dataset, Fn was detected in 35/152 (23.0%) of fresh-frozen tumor samples and in 6/57 (10.5%) of paired normal adjacent tissue, with higher levels in the tumor (p = 0.0033). Fn DNA in the tumor tissue was significantly associated with proximal tumors (p = 0.001), higher depth of invasion (p = 0.014), higher clinical stages (p = 0.033), poor differentiation (p = 0.011), MSI-positive status (p < 0.0001), BRAF mutated tumors (p < 0.0001), and the loss of expression of mismatch-repair proteins MLH1 (p < 0.0001), MSH2 (p = 0.003), and PMS2 (p < 0.0001). Moreover, the presence of Fn DNA in CRC tissue was also associated with a worse patient cancer-specific survival (69.9 vs. 82.2% in 5 years; p = 0.028) and overall survival (63.5 vs. 76.5%; p = 0.037). Here we report, for the first time, the association of F. nucleatum presence with important clinical and molecular features in a Brazilian cohort of CRC patients. Tumor detection and classification based on the gut microbiome might provide a promising approach to improve the prediction of patient outcome.

Guarana (Paullinia cupana Mart.) alters gut microbiota and modulates redox status, partially via caffeine in Wistar rats

Microbiota alterations are observed in pathological conditions, and their regulation is a subject of great interest. Gut microbes are affected by diet, and plant polyphenols may have positive effect on gut microbiota; however, plant-derived extracts may have toxic effects. Guarana (Paullinia cupana Mart.) is a nontraditional medicinal plant applied worldwide. Guarana yields an alkaloid and polyphenol-rich seed with antimicrobial, antioxidant, and anti-inflammatory properties, where caffeine is the major compound. We evaluated the effects of guarana seed powder (GSP) and purified caffeine on gut microbial composition and redox and inflammatory parameters in Wistar rats after 21 days of treatment. Fecal microbiota was analyzed utilizing 16S rDNA sequencing. Antioxidant enzymes activities from liver, kidney, and colon, as well as oxidative damage markers, were evaluated. Total nonenzymatic antioxidant potential was also evaluated. Microbiota was altered by both treatments, GSP and caffeine, without loss of diversity. In the liver, the kidney, and the colon, we observed a decrease in the antioxidant enzymes activities in the GSP group with no increase in the expression of oxidative damage markers, although some enzymes were also regulated by caffeine. Taken together, these results suggested that GSP ameliorates redox parameters but negatively affected gut microbiota, partially via caffeine.

Comparison of the nasopharynx microbiome between influenza and non-influenza cases of severe acute respiratory infections: A pilot study

AIMS

Influenza A virus (IAV) can cause severe acute respiratory infection (SARI), and disease outcome may be associated with changes in the microbiome of the nasopharynx. This is a pilot study to characterize the microbiome of the nasopharynx in patients hospitalized with SARI, infected and not infected by IAV.

METHODS AND RESULTS

Using target sequencing of the 16S rRNA gene, we assessed the bacterial community of nasopharyngeal aspirate samples and compared the microbiome of patients infected with IAV with the microbiome of patients who were negative for IAV. We observed differences in the relative abundance of Proteobacteria and Firmicutes between SARI patients, with Streptococcus being enriched and Pseudomonas underrepresented in IAV patients compared with patients who were not infected with IAV.

CONCLUSION

Pseudomonas taxon seems to be in high frequency on the nasopharynx of SARI patients with non-IAV infection and might present a negative association with Streptococcus taxon. Microbial profile appears to be different between SARI patients infected or not infected with IAV.

Mining of potential drug targets through the identification of essential and analogous enzymes in the genomes of pathogens of Glycine max, Zea mays and Solanum lycopersicum

Pesticides are one of the most widely used pest and disease control measures in plant crops and their indiscriminate use poses a direct risk to the health of populations and environment around the world. As a result, there is a great need for the development of new, less toxic molecules to be employed against plant pathogens. In this work, we employed an in silico approach to study the genes coding for enzymes of the genomes of three commercially important plants, soybean (Glycine max), tomato (Solanum lycopersicum) and corn (Zea mays), as well as 15 plant pathogens (4 bacteria and 11 fungi), focusing on revealing a set of essential and non-homologous isofunctional enzymes (NISEs) that could be prioritized as drug targets. By combining sequence and structural data, we obtained an initial set of 568 cases of analogy, of which 97 were validated and further refined, revealing a subset of 29 essential enzymatic activities with a total of 119 different structural forms, most belonging to central metabolic routes, including the carbohydrate metabolism, the metabolism of amino acids, among others. Further, another subset of 26 enzymatic activities possess a tertiary structure specific for the pathogen, not present in plants, men and Apis mellifera, which may be of importance for the development of specific enzymatic inhibitors against plant diseases that are less harmful to humans and the environment.

The Comparative Genomics and Phylogenomics of Leishmania amazonensis Parasite

Leishmaniasis is an infectious disease caused by Leishmania species. Leishmania amazonensis is a New World Leishmania species belonging to the Mexicana complex, which is able to cause all types of leishmaniasis infections. The L. amazonensis reference strain MHOM/BR/1973/M2269 was sequenced identifying 8,802 codifying sequences (CDS), most of them of hypothetical function. Comparative analysis using six Leishmania species showed a core set of 7,016 orthologs. L. amazonensis and Leishmania mexicana share the largest number of distinct orthologs, while Leishmania braziliensis presented the largest number of inparalogs. Additionally, phylogenomic analysis confirmed the taxonomic position for L. amazonensis within the “Mexicana complex”, reinforcing understanding of the split of New and Old World Leishmania. Potential non-homologous isofunctional enzymes (NISE) were identified between L. amazonensis and Homo sapiens that could provide new drug targets for development.