Insights into the evolution, virulence and speciation of Babesia MO1 and Babesia divergens through multiomics analyses.

Babesiosis, caused by protozoan parasites of the genus Babesia, is an emerging tick-borne disease of significance for both human and animal health. Babesia parasites infect erythrocytes of vertebrate hosts where they develop and multiply rapidly to cause the pathological symptoms associated with the disease. The identification of new Babesia species underscores the ongoing risk of zoonotic pathogens capable of infecting humans, a concern amplified by anthropogenic activities and environmental changes. One such pathogen, Babesia MO1, previously implicated in severe cases of human babesiosis in the United States, was initially considered a subspecies of B. divergens, the predominant agent of human babesiosis in Europe. Here we report comparative multiomics analyses of B. divergens and B. MO1 that offer insight into their biology and evolution. Our analysis shows that despite their highly similar genomic sequences, substantial genetic and genomic divergence occurred throughout their evolution resulting in major differences in gene functions, expression and regulation, replication rates and susceptibility to antiparasitic drugs. Furthermore, both pathogens have evolved distinct classes of multigene families, crucial for their pathogenicity and adaptation to specific mammalian hosts. Leveraging genomic information for B. MO1, B. divergens, and other members of the Babesiidae family within Apicomplexa provides valuable insights into the evolution, diversity, and virulence of these parasites. This knowledge serves as a critical tool in preemptively addressing the emergence and rapid transmission of more virulent strains.

A comparative genomic study of a hydrocarbon-degrading marine bacterial consortium.

Ocean oil pollution has a large impact on the environment and the health of living organisms. Bioremediation cleaning strategies are promising eco-friendly alternatives for tackling this problem. Previously, we designed and reported a hydrocarbon (HC) degrading microbial consortium of four marine strains belonging to the species Alloalcanivorax xenomutans, Halopseudomonas aestusnigri, Paenarthrobacter sp., and Pseudomonas aeruginosa. However, the knowledge about the metabolic potential of this bacterial consortium for HC bioremediation is not yet well understood. Here, we analyzed the complete genomes of these marine bacterial strains accompanied by a phylogenetic reconstruction along with 138 bacterial strains. Synteny between complete genomes of the same species or genus, revealed high conservation among strains of the same species, covering over 91% of their genomic sequences. Functional predictions highlighted a high abundance of genes related to HC degradation, which may result in functional redundancy within the consortium; however, unique and complete gene clusters linked to aromatic degradation were found in the four genomes, suggesting substrate specialization. Pangenome gain and loss analysis of genes involved in HC degradation provided insights into the evolutionary history of these capabilities, shedding light on the acquisition and loss of relevant genes related to alkane and aromatic degradation. Our work, including comparative genomic analyses, identification of secondary metabolites, and prediction of HC-degrading genes, enhances our understanding of the functional diversity and ecological roles of these marine bacteria in crude oil-contaminated marine environments and contributes to the applied knowledge of bioremediation.

Genomic sequence of three hydrocarbonoclastic pseudomonadaceae strains isolated from marine sediments of the port of Rosarito, Baja California, Mexico.

We report the draft genome sequence of three marine bacteria belonging to Pseudomonas and Stutzerimonas genera, with hydrocarbonoclastic metabolism for oil and monoaromatic hydrocarbon degradation. The genomic information of these organisms contributes to the knowledge of natural and polluted marine environments with ubiquitous presence of hydrocarbons as a selective pressure.

Physiology and comparative genomics of the haloalkalitolerant and hydrocarbonoclastic marine strain Rhodococcus ruber MSA14.

Marine hydrocarbonoclastic bacteria can use polycyclic aromatic hydrocarbons as carbon and energy sources, that makes these bacteria highly attractive for bioremediation in oil-polluted waters. However, genomic and metabolic differences between species are still the subject of study to understand the evolution and strategies to degrade PAHs. This study presents Rhodococcus ruber MSA14, an isolated bacterium from marine sediments in Baja California, Mexico, which exhibits adaptability to saline environments, a high level of intrinsic pyrene tolerance (> 5 g L- 1), and efficient degradation of pyrene (0.2 g L- 1) by 30% in 27 days. Additionally, this strain demonstrates versatility by using naphthalene and phenanthrene as individual carbon sources. The genome sequencing of R. ruber MSA14 revealed a genome spanning 5.45 Mbp, a plasmid of 72 kbp, and three putative megaplasmids, lengths between 110 and 470 Kbp. The bioinformatics analysis of the R. ruber MSA14 genome revealed 56 genes that encode enzymes involved in the peripheral and central pathways of aromatic hydrocarbon catabolism, alkane, alkene, and polymer degradation. Within its genome, R. ruber MSA14 possesses genes responsible for salt tolerance and siderophore production. In addition, the genomic analysis of R. ruber MSA14 against 13 reference genomes revealed that all compared strains have at least one gene involved in the alkanes and catechol degradation pathway. Overall, physiological assays and genomic analysis suggest that R. ruber MSA14 is a new haloalkalitolerant and hydrocarbonoclastic strain toward a wide range of hydrocarbons, making it a promising candidate for in-depth characterization studies and bioremediation processes as part of a synthetic microbial consortium, as well as having a better understanding of the catabolic potential and functional diversity among the Rhodococci group.

Exploration of low-frequency allelic variants of SARS-CoV-2 genomes reveals coinfections in Mexico occurred during periods of VOCs turnover.

A total of 14 973 alleles in 29 661 sequenced samples collected between March 2021 and January 2023 by the Mexican Consortium for Genomic Surveillance (CoViGen-Mex) and collaborators were used to construct a thorough map of mutations of the Mexican SARS-CoV-2 genomic landscape containing Intra-Patient Minor Allelic Variants (IPMAVs), which are low-frequency alleles not ordinarily present in a genomic consensus sequence. This additional information proved critical in identifying putative coinfecting variants included alongside the most common variants, B.1.1.222, B.1.1.519, and variants of concern (VOCs) Alpha, Gamma, Delta, and Omicron. A total of 379 coinfection events were recorded in the dataset (a rate of 1.28 %), resulting in the first such catalogue in Mexico. The most common putative coinfections occurred during the spread of Delta or after the introduction of Omicron BA.2 and its descendants. Coinfections occurred constantly during periods of variant turnover when more than one variant shared the same niche and high infection rate was observed, which was dependent on the local variants and time. Coinfections might occur at a higher frequency than customarily reported, but they are often ignored as only the consensus sequence is reported for lineage identification.

Comparative Transcriptome Analysis of Babesia bigemina Attenuated Vaccine and Virulent Strains of Mexican Origin.

Bovine babesiosis, caused by the protozoan Babesia bigemina, is one of the most important hemoparasite diseases of cattle in Mexico and the world. An attenuated B. bigemina strain maintained under in vitro culture conditions has been used as a live attenuated vaccine; however, the biological mechanisms involved in attenuation are unknown. The objective of this study was to identify, through a comparative transcriptomics approach, the components of the B. bigemina virulent parasites that are differentially expressed in vivo, as opposed to those expressed by B. bigemina attenuated vaccine parasites when inoculated into naïve cattle. The biological material under study was obtained by inoculating spleen-intact cattle with infected erythrocytes containing either the attenuated strain or a virulent field strain. After RNA extraction, transcriptomic analysis (RNA-seq) was performed, followed by bioinformatic Differential Expression (DE) analysis and Gene Ontology (GO) term enrichment. The high-throughput sequencing results obtained by analyzing three biological replicates for each parasite strain ranged from 9,504,000 to 9,656,000, and 13,400,000 to 15,750,000 reads for the B. bigemina attenuated and virulent strains, respectively. At least 519 differentially expressed genes were identified in the analyzed strains. In addition, GO analysis revealed both similarities and differences across the three categories: cellular components, biological processes, and molecular functions. The attenuated strain of B. bigemina derived from in vitro culture presents global transcriptomic changes when compared to the virulent strain. Moreover, the obtained data provide insights into the potential molecular mechanisms associated with the attenuation or pathogenicity of each analyzed strain, offering molecular markers that might be associated with virulence or potential vaccine candidates.

Multiomics analysis reveals B. MO1 as a distinct Babesia species and provides insights into its evolution and virulence.

Babesiosis, caused by protozoan parasites of the genus Babesia , is an emerging tick-borne disease of significance for both human and animal health. Babesia parasites infect erythrocytes of vertebrate hosts where they develop and multiply rapidly to cause the pathological symptoms associated with the disease. The identification of various Babesia species underscores the ongoing risk of new zoonotic pathogens capable of infecting humans, a concern amplified by anthropogenic activities and environmental shifts impacting the distribution and transmission dynamics of parasites, their vectors, and reservoir hosts. One such species, Babesia MO1, previously implicated in severe cases of human babesiosis in the midwestern United States, was initially considered closely related to B. divergens , the predominant agent of human babesiosis in Europe. Yet, uncertainties persist regarding whether these pathogens represent distinct variants of the same species or are entirely separate species. We show that although both B. MO1 and B. divergens share similar genome sizes, comprising three nuclear chromosomes, one linear mitochondrial chromosome, and one circular apicoplast chromosome, major differences exist in terms of genomic sequence divergence, gene functions, transcription profiles, replication rates and susceptibility to antiparasitic drugs. Furthermore, both pathogens have evolved distinct classes of multigene families, crucial for their pathogenicity and adaptation to specific mammalian hosts. Leveraging genomic information for B. MO1, B. divergens , and other members of the Babesiidae family within Apicomplexa provides valuable insights into the evolution, diversity, and virulence of these parasites. This knowledge serves as a critical tool in preemptively addressing the emergence and rapid transmission of more virulent strains.

SARS-CoV-2 Omicron variants BA.4 and BA.5 dominated the fifth COVID-19 epidemiological wave in Mexico.

In Mexico, the BA.4 and BA.5 Omicron variants dominated the fifth epidemic wave (summer 2022), superseding BA.2, which had circulated during the inter-wave period. The present study uses genome sequencing and statistical and phylogenetic analyses to examine these variants' abundance, distribution, and genetic diversity in Mexico from April to August 2022. Over 35 % of the sequenced genomes in this period corresponded to the BA.2 variant, 8 % to the BA.4 and 56 % to the BA.5 variant. Multiple subvariants were identified, but the most abundant, BA.2.9, BA.2.12.1, BA.5.1, BA.5.2, BA.5.2.1 and BA.4.1, circulated across the entire country, not forming geographical clusters. Contrastingly, other subvariants exhibited a geographically restricted distribution, most notably in the Southeast region, which showed a distinct subvariant dynamic. This study supports previous results showing that this region may be a significant entry point and contributed to introducing and evolving novel variants in Mexico. Furthermore, a differential distribution was observed for certain subvariants among specific States through time, which may have contributed to the overall increased diversity observed during this wave compared to the previous ones. This study highlights the importance of sustaining genomic surveillance to identify novel variants that may impact public health.

Draft genome sequences of six bacteria isolated from polylactic acid samples from the Mediterranean Sea.

Here, we report the draft genome sequences of six marine strains isolated from plastic samples incubated in the Mediterranean Sea. Genomic analyses place these strains within the Alkalihalobacillus, Bacillus, Halomonas, and Marinobacter genera. Examining the genomes of these non-typical environmental bacteria increases our comprehension of microorganism biology and their potential uses.

Identification of Cry toxin receptor genes homologs in a de novo transcriptome of Premnotrypes vorax (Coleoptera: Curculionidae).

The white potato worm Premnotrypes vorax (Hustache) (Coleoptera: Curculionidae) is one of the most destructive insect pests of potato crops in South America. Like many coleopteran insects, P. vorax shows low susceptibility to Cry insecticidal proteins produced by the bacterium Bacillus thuringiensis (Bt). However, the presence of Cry toxin receptors in the midgut of this this insect has never been studied. The main Cry-binding proteins described in other insect species are cadherin (CAD), aminopeptidase N (APN), alkaline phosphatase (ALP) and ATP-binding cassette (ABC) transporters. In this study, we analyzed and validated a de novo assembled transcriptome of Illumina sequencing data to identify and to characterize homologs of Cry toxin receptors. We identified the protein sequences in P. vorax that show high identity with their orthologous sequences of the Cry toxin binding proteins in other coleopteran larvae such as APN, ALP, CAD and ABC transporter. This study provides preliminary identification of putative receptor genes of Cry proteins that would be useful for future studies involving biocontrol of this important potato crop pest.

Multiple approaches of cellular metabolism define the bacterial ancestry of mitochondria.

We breathe at the molecular level when mitochondria in our cells consume oxygen to extract energy from nutrients. Mitochondria are characteristic cellular organelles that derive from aerobic bacteria and carry out oxidative phosphorylation and other key metabolic pathways in eukaryotic cells. The precise bacterial origin of mitochondria and, consequently, the ancestry of the aerobic metabolism of our cells remain controversial despite the vast genomic information that is now available. Here, we use multiple approaches to define the most likely living relatives of the ancestral bacteria from which mitochondria originated. These bacteria live in marine environments and exhibit the highest frequency of aerobic traits and genes for the metabolism of fundamental lipids that are present in the membranes of eukaryotes, sphingolipids, and cardiolipin.

Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico.

Over 200 different SARS-CoV-2 lineages have been observed in Mexico by November 2021. To investigate lineage replacement dynamics, we applied a phylodynamic approach and explored the evolutionary trajectories of five dominant lineages that circulated during the first year of local transmission. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in Mexico. Lineages B.1.1.222 and B.1.1.519 exhibited similar dynamics, constituting clades that likely originated in Mexico and persisted for >12 months. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. For the largest B.1.617.2 clades, we further explored viral lineage movements across Mexico. Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

Babesia duncani multi-omics identifies virulence factors and drug targets.

Babesiosis is a malaria-like disease in humans and animals that is caused by Babesia species, which are tick-transmitted apicomplexan pathogens. Babesia duncani causes severe to lethal infection in humans, but despite the risk that this parasite poses as an emerging pathogen, little is known about its biology, metabolic requirements or pathogenesis. Unlike other apicomplexan parasites that infect red blood cells, B. duncani can be continuously cultured in vitro in human erythrocytes and can infect mice resulting in fulminant babesiosis and death. We report comprehensive, detailed molecular, genomic, transcriptomic and epigenetic analyses to gain insights into the biology of B. duncani. We completed the assembly, 3D structure and annotation of its nuclear genome, and analysed its transcriptomic and epigenetics profiles during its asexual life cycle stages in human erythrocytes. We used RNA-seq data to produce an atlas of parasite metabolism during its intraerythrocytic life cycle. Characterization of the B. duncani genome, epigenome and transcriptome identified classes of candidate virulence factors, antigens for diagnosis of active infection and several attractive drug targets. Furthermore, metabolic reconstitutions from genome annotation and in vitro efficacy studies identified antifolates, pyrimethamine and WR-99210 as potent inhibitors of B. duncani to establish a pipeline of small molecules that could be developed as effective therapies for the treatment of human babesiosis.

SARS-CoV-2 BW lineage, a fast-growing Omicron variant from southeast Mexico bearing relevant escape mutations.

PURPOSE: The swift expansion of the BW.1 SARS-CoV-2 variant coincided with a rapid increase of COVID-19 cases occurring in Southeast Mexico in October, 2022, which marked the start of Mexico's sixth epidemiological wave. In Yucatan, up to 92% (58 of 73) of weekly sequenced genomes between epidemiological week 42 and 47 were identified as either BW.1 or its descendant, BW.1.1 in the region, during the last trimester of 2022. In the current study, a comprehensive genomic comparison was carried out to characterize the evolutionary history of the BW lineage, identifying its origins and its most important mutations. METHODS: An alignment of all the genomes of the BW lineage and its parental BA.5.6.2 variant was carried out to identify their mutations. A phylogenetic and ancestral sequence reconstruction analysis with geographical inference, as well as a longitudinal analysis of point mutations, were performed to trace back their origin and contrast them with key RBD mutations in variant BQ.1, one of the fastest-growing lineages to date. RESULTS: Our ancestral reconstruction analysis portrayed Mexico as the most probable origin of the BW.1 and BW.1.1 variants. Two synonymous substitutions, T7666C and C14599T, support their Mexican origin, whereas other two mutations are specific to BW.1: S:N460K and ORF1a:V627I. Two additional substitutions and a deletion are found in its descending subvariant, BW.1.1. Mutations found in the receptor binding domain, S:K444T, S:L452R, S:N460K, and S:F486V in BW.1 have been reported to be relevant for immune escape and are also key mutations in the BQ.1 lineage. CONCLUSIONS: BW.1 appears to have arisen in the Yucatan Peninsula in Southeast Mexico sometime around July 2022 during the fifth COVID-19 wave. Its rapid growth may be in part explained by the relevant escape mutations also found in BQ.1.

Omicron-BA.1 Dispersion Rates in Mexico Varied According to the Regional Epidemic Patterns and the Diversity of Local Delta Subvariants.

PURPOSE: The Omicron subvariant BA.1 of SARS-CoV-2 was first detected in November 2021 and quickly spread worldwide, displacing the Delta variant. In this work, a characterization of the spread of this variant in Mexico is presented. METHODS: The time to fixation of BA.1, the diversity of Delta sublineages, the population density, and the level of virus circulation during the inter-wave interval were determined to analyze differences in BA.1 spread. RESULTS: BA.1 began spreading during the first week of December 2021 and became dominant in the next three weeks, causing the fourth COVID-19 epidemiological surge in Mexico. Unlike previous variants, BA.1 did not exhibit a geographically distinct circulation pattern. However, a regional difference in the speed of the replacement of the Delta variant was observed. CONCLUSIONS: Viral diversity and the relative abundance of the virus in a particular area around the time of the introduction of a new lineage seem to have influenced the spread dynamics, in addition to population density. Nonetheless, if there is a significant difference in the fitness of the variants, or if the time allowed for the competition is sufficiently long, it seems the fitter virus will eventually become dominant, as observed in the eventual dominance of the BA.1.x variant in Mexico.

Argonaute5 and its associated small RNAs modulate the transcriptional response during the rhizobia-Phaseolus vulgaris symbiosis.

Both plant- and rhizobia-derived small RNAs play an essential role in regulating the root nodule symbiosis in legumes. Small RNAs, in association with Argonaute proteins, tune the expression of genes participating in nodule development and rhizobial infection. However, the role of Argonaute proteins in this symbiosis has been overlooked. In this study, we provide transcriptional evidence showing that Argonaute5 (AGO5) is a determinant genetic component in the root nodule symbiosis in Phaseolus vulgaris. A spatio-temporal transcriptional analysis revealed that the promoter of PvAGO5 is active in lateral root primordia, root hairs from rhizobia-inoculated roots, nodule primordia, and mature nodules. Transcriptional analysis by RNA sequencing revealed that gene silencing of PvAGO5 affected the expression of genes involved in the biosynthesis of the cell wall and phytohormones participating in the rhizobial infection process and nodule development. PvAGO5 immunoprecipitation coupled to small RNA sequencing revealed the small RNAs bound to PvAGO5 during the root nodule symbiosis. Identification of small RNAs associated to PvAGO5 revealed miRNAs previously known to participate in this symbiotic process, further supporting a role for AGO5 in this process. Overall, the data presented shed light on the roles that PvAGO5 plays during the root nodule symbiosis in P. vulgaris.

Metagenomic analysis reveals differences in the co-occurrence and abundance of viral species in SARS-CoV-2 patients with different severity of disease.

BACKGROUND: SARS-CoV-2 infections have a wide spectrum of clinical manifestations whose causes are not completely understood. Some human conditions predispose to severe outcome, like old age or the presence of comorbidities, but many other facets, including coinfections with other viruses, remain poorly characterized. METHODS: In this study, the eukaryotic fraction of the respiratory virome of 120 COVID-19 patients was characterized through whole metagenomic sequencing. RESULTS: Genetic material from respiratory viruses was detected in 25% of all samples, whereas human viruses other than SARS-CoV-2 were found in 80% of them. Samples from hospitalized and deceased patients presented a higher prevalence of different viruses when compared to ambulatory individuals. Small circular DNA viruses from the Anneloviridae (Torque teno midi virus 8, TTV-like mini virus 19 and 26) and Cycloviridae families (Human associated cyclovirus 10), Human betaherpesvirus 6, were found to be significantly more abundant in samples from deceased and hospitalized patients compared to samples from ambulatory individuals. Similarly, Rotavirus A, Measles morbillivirus and Alphapapilomavirus 10 were significantly more prevalent in deceased patients compared to hospitalized and ambulatory individuals. CONCLUSIONS: Results show the suitability of using metagenomics to characterize a broader peripheric virological landscape of the eukaryotic virome in SARS-CoV-2 infected patients with distinct disease outcomes. Identified prevalent viruses in hospitalized and deceased patients may prove important for the targeted exploration of coinfections that may impact prognosis.

Superspreading event of Covid-19 in adolescents: is there a difference between the vaccinated and the unvaccinated?

Not available.

The New SARS-CoV-2 Variants and Their Epidemiological Impact in Mexico.

The COVID-19 disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus started its deadly journey into a global pandemic in Wuhan, China, in December 2019, where it was first isolated. Subsequently, multiple variants of the virus have been identified worldwide. In this review, we discuss the overall landscape of the pandemic in Mexico, including its most prevalent variants, their surveillance at a genomic level, and how they impacted the epidemiology of the disease. We also evaluate the heterologous vaccination in Mexico and how it may have influenced group immunity and helped mitigate the pandemic. Finally, we present an integrated view that could help us to understand the pandemic and serve as an example of the situation in Latin America.

Dominance of Three Sublineages of the SARS-CoV-2 Delta Variant in Mexico.

In this study, we analyzed the sequences of SARS-CoV-2 isolates of the Delta variant in Mexico, which has completely replaced other previously circulating variants in the country due to its transmission advantage. Among all the Delta sublineages that were detected, 81.5 % were classified as AY.20, AY.26, and AY.100. According to publicly available data, these only reached a world prevalence of less than 1%, suggesting a possible Mexican origin. The signature mutations of these sublineages are described herein, and phylogenetic analyses and haplotype networks are used to track their spread across the country. Other frequently detected sublineages include AY.3, AY.62, AY.103, and AY.113. Over time, the main sublineages showed different geographical distributions, with AY.20 predominant in Central Mexico, AY.26 in the North, and AY.100 in the Northwest and South/Southeast. This work describes the circulation, from May to November 2021, of the primary sublineages of the Delta variant associated with the third wave of the COVID-19 pandemic in Mexico and highlights the importance of SARS-CoV-2 genomic surveillance for the timely identification of emerging variants that may impact public health.