Intense P.1 (Gamma) diversification followed by rapid Delta substitution in Southern Brazil: a SARS-CoV-2 genomic epidemiology study.

The analyses of genetic traits, dispersion patterns and phylogenomics are essential for understanding the evolutionary forces driving SARS-CoV-2 viruses in these three years of COVID-19 pandemics. Brazil is one of the most affected countries in the world and not sufficient genomic studies have been performed. The emergence of P.1 lineage led to one of the most serious public health crises on record. Our study presents the genomic sequencing and characterization of 412 samples from Rio Grande do Sul state, in the Brazilian Southern region, during Gamma and Delta epidemic waves, in 2021. Additionally, molecular evolution tests were performed to identify positively selected sites in Brazil between 2020 and 2022, as well as offer some evolutionary perspective about the maintenance of multiple spike mutations in Omicron lineages. Genomic epidemiology analysis has indicated an intense P.1 (Gamma) diversification followed by rapid Delta substitution in Southern Brazil.

Genomic characterization and molecular evolution of human monkeypox viruses.

Monkeypox virus is a member of the family Poxviridae, as are variola virus and vaccinia virus. It has a linear double-strand DNA genome approximately 197 kb long, containing ~190 non-overlapping ORFs. Comparison of members of the Central and West African clades shows the presence of unique genes that are associated with different disease presentations, depending on the strain. The last smallpox vaccination efforts ended in the mid-1980s, and there is concern about the recent spread of human monkeypox disease around the world. Almost 87,000 human monkeypox cases have been diagnosed in the world, of which more than 10,900 were in Brazil. The aim of this study was to evaluate the epidemiology and molecular evolution of hMpxV. From computational biology analysis of 640 hMpxV genomes from 1962 to 2022, synteny breaks and gene conservation were observed between Central and West clade genomes, and strains belonged with the 2022 outbreak assigned to the West African clade. Evidence was found for diversifying selective pressure at specific sites within protein coding sequences, acting on immunomodulatory processes. The existence of different sites under diversifying and purifying selection in paralog genes indicates adaptive mechanisms underlying the host-pathogen interaction of monkeypox virus in humans.

PCR-RLFP characterization of Leishmania spp. in domestic animals from the south-western border of Brazil.

The aim of this study was to characterize Leishmania spp. from canine and feline samples using Polymerase Chain Reaction (PCR)- Restriction Fragment Length Polymorphism (RFLP). It was conducted in the southern region of Brazil, located at border crossings to Argentina and Uruguay. Samples were collected from 116 dogs (Canis lupus familiaris) and 89 cats (Felis catus). The PCR was performed to screen for an LT1 fragment from kinetoplast DNA (kDNA) target gene, and positive samples were subjected to a second PCR for an internal transcribed spacers (ITS1) region from ribosomal DNA (rDNA) target. RFLP was performed using the Haemophilus aegyptius (HAE III) restriction endonuclease (Fermentas ®). Positive samples by PCR ITS1 were sequenced and deposited in NCBI GenBank, and a phylogenetic analysis was developed. We found that 12.9% (15/116) of the samples from dogs were positive. All the 89 cat samples were negative. Positive samples were tested against Leishmania reference strains presenting different patterns in PCR-RFLP, and these samples showed bands denoting similarity to the standard species of Leishmania infantum, proven through sequencing and phylogenetic analysis. The RFLP technique, alone, was shown to be feasible for practical application and confirmation of the involved Leishmania spp.

Comparative Genomics and Characterization of SARS-CoV-2 P.1 (Gamma) Variant of Concern From Amazonas, Brazil.

BACKGROUND: P.1 lineage (Gamma) was first described in the State of Amazonas, northern Brazil, in the end of 2020, and has emerged as a very important variant of concern (VOC) of SARS-CoV-2 worldwide. P.1 has been linked to increased infectivity, higher mortality, and immune evasion, leading to reinfections and potentially reduced efficacy of vaccines and neutralizing antibodies. METHODS: The samples of 276 patients from the State of Amazonas were sent to a central referral laboratory for sequencing by gold standard techniques, through Illumina MiSeq platform. Both global and regional phylogenetic analyses of the successfully sequenced genomes were conducted through maximum likelihood method. Multiple alignments were obtained including previously obtained unique human SARS-CoV-2 sequences. The evolutionary histories of spike and non-structural proteins from ORF1a of northern genomes were described and their molecular evolution was analyzed for detection of positive (FUBAR, FEL, and MEME) and negative (FEL and SLAC) selective pressures. To further evaluate the possible pathways of evolution leading to the emergence of P.1, we performed specific analysis for copy-choice recombination events. A global phylogenomic analysis with subsampled P.1 and B.1.1.28 genomes was applied to evaluate the relationship among samples. RESULTS: Forty-four samples from the State of Amazonas were successfully sequenced and confirmed as P.1 (Gamma) lineage. In addition to previously described P.1 characteristic mutations, we find evidence of continuous diversification of SARS-CoV-2, as rare and previously unseen P.1 mutations were detected in spike and non-structural protein from ORF1a. No evidence of recombination was found. Several sites were demonstrated to be under positive and negative selection, with various mutations identified mostly in P.1 lineage. According to the Pango assignment, phylogenomic analyses indicate all samples as belonging to the P.1 lineage. CONCLUSION: P.1 has shown continuous evolution after its emergence. The lack of clear evidence for recombination and the positive selection demonstrated for several sites suggest that this lineage emergence resulted mainly from strong evolutionary forces and progressive accumulation of a favorable signature set of mutations.

Predominance of the SARS-CoV-2 Lineage P.1 and Its Sublineage P.1.2 in Patients from the Metropolitan Region of Porto Alegre, Southern Brazil in March 2021.

Almost a year after the COVID-19 pandemic had begun, new lineages (B.1.1.7, B.1.351, P.1, and B.1.617.2) associated with enhanced transmissibility, immunity evasion, and mortality were identified in the United Kingdom, South Africa, and Brazil. The previous most prevalent lineages in the state of Rio Grande do Sul (RS, Southern Brazil), B.1.1.28 and B.1.1.33, were rapidly replaced by P.1 and P.2, two B.1.1.28-derived lineages harboring the E484K mutation. To perform a genomic characterization from the metropolitan region of Porto Alegre, we sequenced viral samples to: (i) identify the prevalence of SARS-CoV-2 lineages in the region, the state, and bordering countries/regions; (ii) characterize the mutation spectra; (iii) hypothesize viral dispersal routes by using phylogenetic and phylogeographic approaches. We found that 96.4% of the samples belonged to the P.1 lineage and approximately 20% of them were assigned as the novel P.1.2, a P.1-derived sublineage harboring signature substitutions recently described in other Brazilian states and foreign countries. Moreover, sequences from this study were allocated in distinct branches of the P.1 phylogeny, suggesting multiple introductions in RS and placing this state as a potential diffusion core of P.1-derived clades and the emergence of P.1.2. It is uncertain whether the emergence of P.1.2 and other P.1 clades is related to clinical or epidemiological consequences. However, the clear signs of molecular diversity from the recently introduced P.1 warrant further genomic surveillance.

Mutation hotspots and spatiotemporal distribution of SARS-CoV-2 lineages in Brazil, February 2020-2021.

The COVID-19 pandemic has already reached more than 110 million people and is associated with 2.5 million deaths worldwide. Brazil is the third worst-hit country, with approximately 10.2 million cases and 250 thousand deaths. International efforts have been established to share information about Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemiology and evolution to support the development of effective strategies for public health and disease management. We aimed to analyze the high-quality genome sequences from Brazil from February 2020-2021 to identify mutation hotspots, geographical and temporal distribution of SARS-CoV-2 lineages by using phylogenetics and phylodynamics analyses. We describe heterogeneous sequencing efforts, the progression of the different lineages along time, evaluating mutational spectra and frequency oscillations derived from the prevalence of specific lineages across different Brazilian regions. We found at least seven major (1-7) and two minor clades related to the six most prevalent lineages in the country and described its spatial distribution and dynamics. The emergence and recent frequency shift of lineages (P.1 and P.2) carrying mutations of concern in the spike protein (e. g., E484K, N501Y) draws attention due to their association with immune evasion and enhanced receptor binding affinity. Improvements in genomic surveillance are of paramount importance and should be extended in Brazil to better inform policy makers about better decisions to fight the COVID-19 pandemic.

E484K as an innovative phylogenetic event for viral evolution: Genomic analysis of the E484K spike mutation in SARS-CoV-2 lineages from Brazil.

The COVID-19 pandemic caused by SARS-CoV-2 has affected millions of people since its beginning in 2019. The propagation of new lineages and the discovery of key mechanisms adopted by the virus to overlap the immune system are central topics for the entire public health policies, research and disease management. Since the second semester of 2020, the mutation E484K has been progressively found in the Brazilian territory, composing different lineages over time. It brought multiple concerns related to the risk of reinfection and the effectiveness of new preventive and treatment strategies due to the possibility of escaping from neutralizing antibodies. To better characterize the current scenario we performed genomic and phylogenetic analyses of the E484K mutated genomes sequenced from Brazilian samples in 2020. From October 2020, more than 40% of the sequenced genomes present the E484K mutation, which was identified in three different lineages (P.1, P.2 and B.1.1.33 - posteriorly renamed as N.9) in four Brazilian regions. We also evaluated the presence of E484K associated mutations and identified selective pressures acting on the spike protein, leading us to some insights about adaptive and purifying selection driving the virus evolution.

Genomic epidemiology of SARS-CoV-2 in Esteio, Rio Grande do Sul, Brazil.

BACKGROUND: Brazil is the third country most affected by Coronavirus disease-2019 (COVID-19), but viral evolution in municipality resolution is still poorly understood in Brazil and it is crucial to understand the epidemiology of viral spread. We aimed to track molecular evolution and spread of Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Esteio (Southern Brazil) using phylogenetics and phylodynamics inferences from 21 new genomes in global and regional context. Importantly, the case fatality rate (CFR) in Esteio (3.26%) is slightly higher compared to the Rio Grande do Sul (RS) state (2.56%) and the entire Brazil (2.74%). RESULTS: We provided a comprehensive view of mutations from a representative sampling from May to October 2020, highlighting two frequent mutations in spike glycoprotein (D614G and V1176F), an emergent mutation (E484K) in spike Receptor Binding Domain (RBD) characteristic of the B.1.351 and P.1 lineages, and the adjacent replacement of 2 amino acids in Nucleocapsid phosphoprotein (R203K and G204R). E484K was found in two genomes from mid-October, which is the earliest description of this mutation in Southern Brazil. Lineages containing this substitution must be subject of intense surveillance due to its association with immune evasion. We also found two epidemiologically-related clusters, including one from patients of the same neighborhood. Phylogenetics and phylodynamics analysis demonstrates multiple introductions of the Brazilian most prevalent lineages (B.1.1.33 and B.1.1.248) and the establishment of Brazilian lineages ignited from the Southeast to other Brazilian regions. CONCLUSIONS: Our data show the value of correlating clinical, epidemiological and genomic information for the understanding of viral evolution and its spatial distribution over time. This is of paramount importance to better inform policy making strategies to fight COVID-19.

Zrg1, a cryptococcal protein associated with regulation of growth in nutrient deprivation conditions.

Cryptococcus gattii is one of the causes of cryptococcosis, a life-threatening disease generally characterized by pneumonia and/or meningitis. Zinc is an essential element for life, being required for the activity of many proteins with catalytic and structural roles. Here, we characterize ZRG1 (zinc-related gene 1), which codes a product involved in zinc metabolism. Transcriptional profiling revealed that zinc availability regulated the expression of ZRG1, and its null mutants demonstrated impaired growth in zinc- and nitrogen-limiting conditions. Moreover, zrg1 strains displayed alterations in the expression of the zinc homeostasis-related genes ZAP1 and ZIP1. Notably, cryptococcal cells lacking Zrg1 displayed upregulation of autophagy-like phenotypes. Despite no differences were detected in the classical virulence-associated traits; cryptococcal cells lacking ZRG1 displayed decreased capacity for survival inside macrophages and attenuated virulence in an invertebrate model. Together, these results indicate that ZRG1 plays an important role in proper zinc metabolism, and is necessary for cryptococcal fitness and virulence.

Transcriptional Analysis Allows Genome Reannotation and Reveals that Cryptococcus gattii VGII Undergoes Nutrient Restriction during Infection.

Cryptococcus gattii is a human and animal pathogen that infects healthy hosts and caused the Pacific Northwest outbreak of cryptococcosis. The inhalation of infectious propagules can lead to internalization of cryptococcal cells by alveolar macrophages, a niche in which C. gattii cells can survive and proliferate. Although the nutrient composition of macrophages is relatively unknown, the high induction of amino acid transporter genes inside the phagosome indicates a preference for amino acid uptake instead of synthesis. However, the presence of countable errors in the R265 genome annotation indicates significant inhibition of transcriptomic analysis in this hypervirulent strain. Thus, we analyzed RNA-Seq data from in vivo and in vitro cultures of C. gattii R265 to perform the reannotation of the genome. In addition, based on in vivo transcriptomic data, we identified highly expressed genes and pathways of amino acid metabolism that would enable C. gattii to survive and proliferate in vivo. Importantly, we identified high expression in three APC amino acid transporters as well as the GABA permease. The use of amino acids as carbon and nitrogen sources, releasing ammonium and generating carbohydrate metabolism intermediaries, also explains the high expression of components of several degradative pathways, since glucose starvation is an important host defense mechanism.

sRNAs as possible regulators of retrotransposon activity in Cryptococcus gattii VGII.

BACKGROUND: The absence of Argonaute genes in the fungal pathogen Cryptococcus gattii R265 and other VGII strains indicates that yeasts of this genotype cannot have a functional RNAi pathway, an evolutionarily conserved gene silencing mechanism performed by small RNAs. The success of the R265 strain as a pathogen that caused the Pacific Northwest and Vancouver Island outbreaks may imply that RNAi machinery loss could be beneficial under certain circumstances during evolution. As a result, a hypermutant phenotype would be created with high rates of genome retrotransposition, for instance. This study therefore aimed to evaluate in silicio the effect of retrotransposons and their control mechanisms by small RNAs on genomic stability and synteny loss of C. gattii R265 through retrotransposons sequence comparison and orthology analysis with other 16 C. gattii genomic sequences available. RESULTS: Retrotransposon mining identified a higher sequence count to VGI genotype compared to VGII, VGIII, and VGIV. However, despite the lower retrotransposon number, VGII exhibited increased synteny loss and genome rearrangement events. RNA-Seq analysis indicated highly expressed retrotransposons as well as sRNA production. CONCLUSIONS: Genome rearrangement and synteny loss may suggest a greater retrotransposon mobilization caused by RNAi pathway absence, but the effective presence of sRNAs that matches retrotransposon sequences means that an alternative retrotransposon silencing mechanism could be active in genomic integrity maintenance of C. gattii VGII strains.