Since the identification of SARS-CoV-2 in December 2019 a large number of SARS-CoV-2 genomes has been sequenced with unprecedented speed around the world and deposited in several databases. This marks a unique opportunity to study how a virus spread and evolve in a worldwide context. However, currently there is not a useful haplotype classification system to help tracking the virus evolution. Here we identified eleven mutations with 10 % or more frequency in a data set of 7848 genomes. Using these mutations, we identified 6 SARS-CoV-2 haplotypes or OTUs (Operational Taxonomic Unit) that correlate well with a phylogenomic tree. After that, we analyzed the geographical and temporal distribution of these OTUs, as well as their correlation with patient status. Our geographical analysis showed different OTUs prevalence between continents and the temporal distribution analysis revealed an evolution-like pattern in SARS-CoV-2. Finally, we observed a homogenous distribution of OTUs in mild and severe patients and a great prevalence of OTU 2 in asymptomatic patients. However, genomes in the asymptomatic category, comes from isolates on three consecutive days in February (15 to 17), weakening this observation and highlighting the need to increase genomic analyzes in asymptomatic and severe patients. Our classification system is phylogenetically consistent and allows us to easily track geographic and temporal distribution of important mutations around the world. In the next months, it could be updated using similar steps that we used here.

Since the identification of SARS-CoV-2, a large number of genomes have been sequenced with unprecedented speed around the world. This marks a unique opportunity to analyze virus spreading and evolution in a worldwide context. Currently, there is not a useful haplotype description to help to track important and globally scattered mutations. Also, differences in the number of sequenced genomes between countries and/or months make it difficult to identify the emergence of haplotypes in regions where few genomes are sequenced but a large number of cases are reported. We propose an approach based on the normalization by COVID-19 cases of relative frequencies of mutations using all the available data to identify major haplotypes. Furthermore, we can use a similar normalization approach to tracking the temporal and geographic distribution of haplotypes in the world. Using 171,461 genomes, we identify five major haplotypes or operational taxonomic units (OTUs) based on nine high-frequency mutations. OTU_3 characterized by mutations R203K and G204R is currently the most frequent haplotype circulating in four of the six continents analyzed (South America, North America, Europe, Asia, Africa, and Oceania). On the other hand, during almost all months analyzed, OTU_5 characterized by the mutation T85I in nsp2 is the most frequent in North America. Recently (since September), OTU_2 has been established as the most frequent in Europe. OTU_1, the ancestor haplotype, is near to extinction showed by its low number of isolations since May. Also, we analyzed whether age, gender, or patient status is more related to a specific OTU. We did not find OTU’s preference for any age group, gender, or patient status. Finally, we discuss structural and functional hypotheses in the most frequently identified mutations, none of those mutations show a clear effect on the transmissibility or pathogenicity.

In this study, the feasibility of using a biohydrometallurgical technique for selective metals recovery from electronic waste (e‐waste) by bacterial bioleaching was investigated. Acidithiobacillus was identified in coal mining acid mine drainage (AMD). The microorganism was studied using specific sequencing of a 16s rDNA fragment. The potential for the dissolution of copper from waste printed wire boards (PWBs) using the isolated Acidithiobacillus ferrooxidans (A. ferroxidans) was evaluated. The bioleaching experiments were performed in an orbital shaker at 30 °C and 170 rpm, with 10 % (v/v) inoculum and a pulp density of 30 g/L. The copper concentration was determined by energy dispersive x‐ray fluorescence (XRF). The result shows that copper recovery from PWBs using our A. ferrooxidans strain was 95 % after 8 days, which showed the feasibility of this process.

The second messenger cyclic diguanylate monophosphate (c-di-GMP) is a central regulator of bacterial lifestyle, controlling several behaviors, including the switch between sessile and motile states. The c-di-GMP levels are controlled by the interplay between diguanylate cyclases (DGCs) and phosphodiesterases, which synthesize and hydrolyze this second messenger, respectively. These enzymes often contain additional domains that regulate activity via binding of small molecules, covalent modification, or protein-protein interactions. A major challenge remains to understand how DGC activity is regulated by these additional domains or interaction partners in specific signaling pathways. Here, we identified a pair of co-transcribed genes ( and ) in the phytopathogenic, Gram-negative bacterium subsp. (Xac), whose mutations resulted in opposing motility phenotypes. We show that the periplasmic cache domain of XAC2382, a membrane-associated DGC, interacts with XAC2383, a periplasmic binding protein, and we provide evidence that this interaction regulates XAC2382 DGC activity. Moreover, we solved the crystal structure of XAC2383 with different ligands, indicating a preference for negatively charged phosphate-containing compounds. We propose that XAC2383 acts as a periplasmic sensor that, upon binding its ligand, inhibits the DGC activity of XAC2382. Of note, we also found that this previously uncharacterized signal transduction system is present in several other bacterial phyla, including Gram-positive bacteria. Phylogenetic analysis of homologs of the XAC2382-XAC2383 pair supports several independent origins that created new combinations of XAC2382 homologs with a conserved periplasmic cache domain with different cytoplasmic output module architectures.

Coronavirus disease 2019 (COVID-19) is a contagious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This disease has spread globally, causing more than 161.5 million cases and 3.3 million deaths to date. Surveillance and monitoring of new mutations in the virus’ genome are crucial to our understanding of the adaptation of SARS-CoV-2. Moreover, how the temporal dynamics of these mutations is influenced by control measures and non-pharmaceutical interventions (NPIs) is poorly understood. Using 1,058,020 SARS-CoV-2 from sequenced COVID-19 cases from 98 countries (totaling 714 country-month combinations), we perform a normalization by COVID-19 cases to calculate the relative frequency of SARS-CoV-2 mutations and explore their dynamics over time. We found 115 mutations estimated to be present in more than 3% of global COVID-19 cases and determined three types of mutation dynamics: high-frequency, medium-frequency, and low-frequency. Classification of mutations based on temporal dynamics enable us to examine viral adaptation and evaluate the effects of implemented control measures in virus evolution during the pandemic. We showed that medium-frequency mutations are characterized by high prevalence in specific regions and/or in constant competition with other mutations in several regions. Finally, taking N501Y mutation as representative of high-frequency mutations, we showed that level of control measure stringency negatively correlates with the effective reproduction number of SARS-CoV-2 with high-frequency or not-high-frequency and both follows similar trends in different levels of stringency.

The recent outbreak of yellow fever (YF) in São Paulo during 2016–2019 has been one of the most severe in the last decades, spreading to areas with low vaccine coverage. The aim of this study was to assess the genetic diversity of the yellow fever virus (YFV) from São Paulo 2016–2019 outbreak, integrating the available genomic data with new genomes from patients from the Hospital das Clínicas da Faculdade de Medicina da Universidade de São Paulo (HCFMUSP). Using phylodynamics, we proposed the existence of new IE subclades, described their sequence signatures, and determined their locations and time of origin. Plasma or urine samples from acute severe YF cases (n = 56) with polymerase chain reaction (PCR) positive to YFV were submitted to viral genome amplification using 12 sets of primers. Thirty-nine amplified genomes were subsequently sequenced using next-generation sequencing (NGS). These 39 sequences, together with all the complete genomes publicly available, were aligned and used to determine nucleotide/amino acids substitutions and perform phylogenetic and phylodynamic analysis. All YFV genomes generated in this study belonged to the genotype South American I subgroup E. Twenty-one non-synonymous substitutions were identified among the new generated genomes. We analyzed two major clades of the genotypes IE, IE1, and IE2 and proposed the existence of subclades based on their sequence signatures. Also, we described the location and time of origin of these subclades. Overall, our findings provide an overview of YFV genomic characterization and phylodynamics of the 2016–2019 outbreak contributing to future virological and epidemiological studies.

Hepatitis B virus (HBV) spreads efficiently among all human populations worldwide. HBV is classified into ten genotypes (A to J) with their geographic distribution and clinical features. In Mexico, HBV genotype H is the leading cause of hepatitis B and has been detected in indigenous populations, suggesting that HBV genotype H may be native to Mexico. However, little is known about the evolutionary history of HBV genotype H. Thus, we aimed to determine the age of HBV genotype H in Mexico using molecular dating techniques. Ninety-two HBV sequences of the reverse transcriptase (RT) domain of the polymerase gene (~1,251 bp) were analyzed; 48 were genotype H, 43 were genotype F, and the oldest HBV sequence from America was included as the root. All sequences were aligned, and the most recent common ancestor (TMRCA) time was calculated using the Bayesian Skyline Evolutionary Analysis. Our results estimate a TMRCA for the genotype H in Mexico of 2070.9 (667.5–4489.2) years before the present (YBP). We identified four major diversification events in genotype H, named H1, H2, H3, and H4. The TMRCA of H1 was 1213.0 (253.3–2638.3) YBP, followed by H2 1175.5 (557.5–2424.2) YBP, H3 949.6 (279.3–2105.0) YBP, and H4 1230.5 (336.3, 2756.7) YBP. We estimated that genotype H diverged from its sister genotype F around 8140.8 (1867.5–18012.8) YBP. In conclusion, this study found that genotype H in Mexico has an estimated age of 2070.9 (667.5–4489.2) YBP and has experienced at least four major diversification events since then.

At over 0.6% of the population, Peru has one of the highest SARS-CoV-2 mortality rate in the world. Much effort to sequence genomes has been done in this country since mid-2020. However, an adequate analysis of the dynamics of the variants of concern and interest (VOCIs) is missing. We investigated the dynamics of the COVID-19 pandemic in Peru with a focus on the second wave, which had the greatest case fatality rate. The second wave in Peru was dominated by Lambda and Gamma. Analysis of the origin of Lambda shows that it most likely emerged in Peru before the second wave (June–November, 2020). After its emergence it reached Argentina and Chile from Peru where it was locally transmitted. During the second wave in Peru, we identify the coexistence of two Lambda and three Gamma sublineages. Lambda sublineages emerged in the center of Peru whereas the Gamma sublineages more likely originated in the north-east and mid-east. Importantly, it is observed that the center of Peru played a prominent role in transmitting SARS-CoV-2 to other regions within Peru.