Pseudomonas aeruginosa accounts for 10–20

The global rise of carbapenem-resistant Acinetobacter baumannii (CRAB) strains poses a critical challenge to healthcare systems due to limited therapeutic options and high mortality rates, especially in intensive care settings. This review explores the epidemiological landscape and molecular mechanisms driving carbapenem resistance, including the production of diverse beta-lactamases (particularly OXA-type enzymes), porin loss, efflux pump overexpression, and mutations in antibiotic targets. Emerging treatment strategies are discussed, such as the use of new beta-lactam–beta-lactamase inhibitor combinations (e.g., sulbactam–durlobactam), siderophore cephalosporins, next-generation polymyxins, as well as novel agents like zosurabalpin and rifabutin (BV100). Alternative approaches—including phage therapy, antimicrobial peptides, CRISPR-based gene editing, and nanoparticle-based delivery systems—are also evaluated for their potential to bypass traditional resistance mechanisms. Furthermore, advances in artificial intelligence and multi-omics integration are highlighted as tools for identifying novel drug targets and predicting resistance profiles. Together, these innovations represent a multifaceted strategy to overcome CRAB infections, yet their successful implementation requires further clinical validation and coordinated surveillance efforts. This analysis highlights the urgent need for continued investment in innovative treatments and effective resistance monitoring to limit the spread of CRAB and protect the effectiveness of last-line antibiotics.

The coronavirus disease 2019 (COVID-19) pandemic has triggered a global health crisis, with over 700 million confirmed cases and at least 7 million deaths reported by early 2024. Children are less vulnerable to severe SARS-CoV-2 infection than adults and typically experience milder respiratory symptoms. However, a rare but significant complication, known as multisystem inflammatory syndrome in children (MIS-C), can develop weeks after infection, characterized by a spectrum of inflammatory symptoms. This study employed whole-exome sequencing and over-representation analysis to identify genetic variants of potential clinical significance related to MIS-C or severe COVID-19 in a group of children with acute respiratory distress syndrome (ARDS), all of whom were unvaccinated for COVID-19. We observed the enrichment of potentially pathogenic genetic variants in genes related to carbohydrate metabolism, particularly glycogen breakdown, in severe COVID-19 pediatric patients, and in genes related to cholesterol and lipoprotein metabolism in MIS-C patients. These findings offer insights into the genetic underpinnings of MIS-C and severe COVID-19, suggesting potential genes and biological pathways for further research.

The coronavirus disease 2019 (COVID-19) pandemic has triggered a global health crisis, with over 700 million confirmed cases and at least 7 million deaths reported by early 2024. Children are less vulnerable to severe SARS-CoV-2 infection than adults and typically experience milder respiratory symptoms. However, a rare but significant complication, known as multisystem inflammatory syndrome in children (MIS-C), can develop weeks after infection, characterized by a spectrum of inflammatory symptoms. This study employed whole-exome sequencing and over-representation analysis to identify distinct pathogenic genetic variants related to MIS-C or severe COVID-19 in a group of children with acute respiratory distress syndrome (ARDS), all of whom were unvaccinated for COVID-19. We observed an enrichment of pathogenic variants in genes related to carbohydrate metabolism, particularly glycogen breakdown, in severe COVID-19 pediatric patients, and in genes related to cholesterol and lipoprotein metabolism in MIS-C patients. These findings offer insights into the genetic underpinnings of MIS-C and severe COVID-19, suggesting potential genes and biological pathways for further research.

Extended-spectrum beta-lactamase producing and ciprofloxacin-non-susceptible Escherichia coli are clinical and environmental issues. We evaluated the susceptibility profile of fosfomycin in non-susceptible E. coli isolated from urine and the environment. We measured the activity of fosfomycin against 319 and 36 E. coli strains from urine and environmental isolates, respectively, collected from rivers. Fosfomycin resistance profiles were investigated using the minimal inhibitory concentration (MIC), according to the Clinical and Laboratory Standards Institute (CLSI) and the European Committee for Antimicrobial Susceptibility Testing (EUCAST) guidelines. Antibiotic susceptibility testing revealed that 5% and 6.6% of urine samples were non-susceptible to fosfomycin according to CLSI and EUCAST guidelines, respectively. The fosfomycin MIC50/90 was 0.5/4 mg/L. Of the 36 E. coli isolates from river water, 11.1% and 13,8% were non-susceptible to fosfomycin according to CLSI and EUCAST, respectively (range ≤0.25 ≥512 mg/L). All the isolates with MIC ≥512 mg/L for fosfomycin showed the fosA3 gene. Fosfomycin resistance was more frequent in the environment than in clinical samples.

BACKGROUND:Enzymatic degradation mediated by beta-lactamases constitutes one of the primary mechanisms of resistance to beta-lactam antibiotics in gram-negative bacteria. This enzyme family comprises four molecular classes, categorized into serine beta-lactamases (Classes A, C, and D) and zinc-dependent metallo-beta-lactamases (Class B). Gram-negative bacteria producing beta-lactamase are of significant concern, particularly due to their prevalence in nosocomial infections. A comprehensive understanding of the evolution and dissemination of this enzyme family is essential for effective control of these pathogens. In this study, we conducted the prospecting, phylogenetic analysis, and in silico analysis of beta-lactamases and homologous proteins identified in 1827 bacterial genomes with phenotypic data on beta-lactam resistance. These genomes were distributed among Klebsiella pneumoniae (45%), Acinetobacter baumannii (31%), Pseudomonas aeruginosa (14%), Escherichia coli (6%), and Enterobacter spp. (4%). Using an HMM profile and searching for conserved domains, we mined 2514, 8733, 5424, and 2957 proteins for molecular classes A, B, C, and D, respectively. This set of proteins encompasses canonical subfamilies of beta-lactamases as well as hypothetical proteins and other functional groups. Canonical beta-lactamases were found to be phylogenetically distant from hypothetical proteins, which, in turn, are closer to other representatives of the penicillin-binding-protein (PBP-like) and metallo-beta-lactamase (MBL) families. The catalytic amino acid residues characteristic of beta-lactamases were identified from the sequence alignment and revealed that motifs are less conserved in homologous groups than in beta-lactamases. After comparing the frequency of protein groups in genomes of resistant strains with those of sensitive ones applying Fisher's exact test and relative risk, it was observed that some groups of homologous proteins to classes B and C are more common in the genomes of resistant strains, particularly to carbapenems. We identified the beta-lactamase-like domain widely distributed in gram-negative species of the ESKAPEE group, which highlights its importance in the context of beta-lactam resistance. Some hypothetical homologous proteins have been shown to potentially possess promiscuous activity against beta-lactam antibiotics, however, they do not appear to expressly determine the resistance phenotype. The selective pressure due to the widespread use of antibiotics may favor the optimization of these functions for specialized resistance enzymes.

Toxoplasma gondii is an obligate intracellular parasite of the phylum Apicomplexa and causes toxoplasmosis infections, a disease that affects a quarter of the world’s population and has no effective cure. Epigenetic regulation is one of the mechanisms controlling gene expression and plays an essential role in all organisms. Lysine deacetylases (KDACs) act as epigenetic regulators affecting gene silencing in many eukaryotes. Here, we focus on TgKDAC4, an enzyme unique to apicomplexan parasites, and a class IV KDAC, the least-studied class of deacetylases so far. This enzyme shares only a portion of the specific KDAC domain with other organisms. Phylogenetic analysis from the TgKDAC4 domain shows a putative prokaryotic origin. Surprisingly, TgKDAC4 is located in the apicoplast, making it the only KDAC found in this organelle to date. Transmission electron microscopy assays confirmed the presence of TgKDAC4 in the periphery of the apicoplast. We identified possible targets or/and partners of TgKDAC4 by immunoprecipitation assays followed by mass spectrometry analysis, including TgCPN60 and TgGAPDH2, both located at the apicoplast and containing acetylation sites. Understanding how the protein works could provide new insights into the metabolism of the apicoplast, an essential organelle for parasite survival.

Dengue virus serotype 2, genotype Cosmopolitan (DENV-2-GII), is one of the most widespread DENV strains globally. In the USA, DENV-2 epidemics have been dominated by DENV-2 genotype Asian-American (DENV-2-GIII), and the first cases of DENV-2-GII were only described in 2019, in Peru, and in 2021 in Brazil. To gain new information about the circulation of DENV-2-GII in Brazil, we sequenced 237 DENV-2 confirmed cases sampled between March 2021 and March 2023 and revealed that DENV-2-GII is already present in all geographic regions of Brazil. The phylogeographic analysis inferred that DENV-2-GII was introduced at least four times in Brazil, between May 2020 and August 2022, generating multiple clades that spread throughout the country with different success. Despite multiple introductions of DENV-2-GII, analysis of the country-wide laboratory surveillance data showed that the Brazilian dengue epidemic in 2022 was dominated by DENV-1 in most states. We hypothesize that massive circulation of DENV-2-GIII in previous years in Brazil might have created a population immune barrier against symptomatic homotypic reinfections by DENV-2-GII, leading to sustained cryptic circulation in asymptomatic cases and localized outbreaks of this new genotype. In summary, our study stresses the importance of arboviral genomic surveillance to close monitoring and better understanding the potential impact of DENV-2-GII in the coming years.

With the emergence of the COVID-19 pandemic, the world has entered into a global health emergency. The increased spread of SARS-CoV-2 has led to reports of severe illness in pediatric patients, simila

Klebsiella pneumoniae is a nosocomial pathogen and an important propagator of multidrug-resistant (MDR) and extensively drug-resistant (XDR) strains. Like other Gram-negative bacteria, they secrete outer membrane vesicles (OMVs) that distribute virulence and resistance factors. Here, we subjected a K. pneumoniae-XDR to subinhibitory concentrations of meropenem, amikacin, polymyxin B, and a combination of these agents to evaluate changes in the protein cargo of OMVs through liquid chromatography-tandem mass spectrometry (LC-MS/MS). Genome sequencing of the clinical isolate K. pneumoniae strain HCD1 (KpHCD1) revealed the presence of 41 resistance genes and 159 virulence factors. We identified 64 proteins in KpHCD1-OMVs modulated with different antibiotic treatments involved in processing genetic information, environmental information, cell envelope formation, energy metabolism, and drug resistance. The OMV proteome expression profile suggests that OMVs may be associated with pathogenicity, survival, stress response, and resistance dissemination.

IntroductionDespite the existing data on the Multisystem Inflammatory Syndrome in Children (MIS-C), the factors that determine these patients evolution remain elusive. Answers may lie, at least in part, in genetics. It is currently under investigation that MIS-C patients may have an underlying innate error of immunity (IEI), whether of monogenic, digenic, or even oligogenic origin.MethodsTo further investigate this hypothesis, 30 patients with MIS-C were submitted to whole exome sequencing. ResultsAnalyses of genes associated with MIS-C, MIS-A, severe covid-19, and Kawasaki disease identified twenty-nine patients with rare potentially damaging variants (50 variants were identified in 38 different genes), including those previously described in IFNA21 and IFIH1 genes, new variants in genes previously described in MIS-C patients (KMT2D, CFB, and PRF1), and variants in genes newly associated to MIS-C such as APOL1, TNFRSF13B, and G6PD. In addition, gene ontology enrichment pointed to the involvement of thirteen major pathways, including complement system, hematopoiesis, immune system development, and type II interferon signaling, that were not yet reported in MIS-C.DiscussionThese data strongly indicate that different gene families may favor MIS- C development. Larger cohort studies with healthy controls and other omics approaches, such as proteomics and RNAseq, will be precious to better understanding the disease dynamics.

The COVID-19 epidemic in Brazil experienced two major country-wide lineage replacements, the first driven by the lineage P.2, formerly classified as variant of interest (VOI) Zeta in late 2020 and the second by the variant of concern (VOC) Gamma in early 2021. To better understand how these SARS-CoV-2 lineage turnovers occurred in Brazil, we analyzed 11,724 high-quality SARS-CoV-2 whole genomes of samples collected in different country regions between September 2020 and April 2021. Our findings indicate that the spatial dispersion of both variants in Brazil was driven by short and long-distance viral transmission. The lineage P.2 harboring Spike mutation E484K probably emerged around late July 2020 in the Rio de Janeiro (RJ) state, which contributed with most (~50%) inter-state viral disseminations, and only became locally established in most Brazilian states by October 2020. The VOC Gamma probably arose in November 2020 in the Amazonas (AM) state, which was responsible for 60-70% of the inter-state viral dissemination, and the earliest timing of community transmission of this VOC in many Brazilian states was already traced to December 2020. We estimate that variant Gamma was 1.56-3.06 more transmissible than variant P.2 co-circulating in RJ and that the median effective reproductive number (Re) of Gamma in RJ and SP states (Re = 1.59-1.91) was lower than in AM (Re = 3.55). In summary, although the epicenter of the lineage P.2 dissemination in Brazil was the heavily interconnected Southeastern region, it displayed a slower rate of spatial spread than the VOC Gamma originated in the more isolated Northern Brazilian region. Our findings also support that the VOC Gamma was more transmissible than lineage P.2, although the viral Re of the VOC varied according to the geographic context.

Cystic fibrosis (CF) is a genetic disease caused by mutations in the cystic fibrosis transmembrane conductance regulator gene that leads to respiratory complications and mortality. Studies have shown shifts in the respiratory microbiota during disease progression in individuals with CF. In addition, CF patients experience short cycles of acute intermittent aggravations of symptoms called pulmonary exacerbations, which may be characterized by a decrease in lung function and weight loss. The resident microbiota become imbalanced, promoting biofilm formation, and reducing the effectiveness of therapy. The aim of this study was to monitor patients aged 8-23 years with CF to evaluate their lower respiratory microbiota using 16S rRNA sequencing. The most predominant pathogens observed in microbiota, Staphylococcus (Staph) and Pseudomonas (Pseud) were correlated with clinical variables, and the in vitro capacity of biofilm formation for these pathogens was tested. A group of 34 patients was followed up for 84 days, and 306 sputum samples were collected and sequenced. Clustering of microbiota by predominant pathogen showed that children with more Staph had reduced forced expiratory volume in one second (FEV1) and forced vital capacity (FVC) compared to children with Pseud. Furthermore, the patients' clinical condition was consistent with the results of pulmonary function. More patients with pulmonary exacerbation were observed in the Staph group than in the Pseud group, as confirmed by lower body mass index and pulmonary function. Additionally, prediction of bacterial functional profiles identified genes encoding key enzymes involved in virulence pathways in the Pseud group. Importantly, this study is the first Brazilian study to assess the lower respiratory microbiota in a significant group of young CF patients. In this sense, the data collected for this study on the microbiota of children in Brazil with CF provide a valuable contribution to the knowledge in the field.

Antimicrobial resistance is an old and silent pandemic. Resistant organisms emerge in parallel with new antibiotics, leading to a major global public health crisis over time. Antibiotic resistance may be due to different mechanisms and against different classes of drugs. These mechanisms are usually found in the same organism, giving rise to multidrug-resistant (MDR) and extensively drug-resistant (XDR) bacteria. One resistance mechanism that is closely associated with the emergence of MDR and XDR bacteria is the efflux of drugs since the same pump can transport different classes of drugs. In Gram-negative bacteria, efflux pumps are present in two configurations: a transmembrane protein anchored in the inner membrane and a complex formed by three proteins. The tripartite complex has a transmembrane protein present in the inner membrane, a periplasmic protein, and a porin associated with the outer membrane. In Pseudomonas aeruginosa, one of the main pathogens associated with respiratory tract infections, four main sets of efflux pumps have been associated with antibiotic resistance: MexAB-OprM, MexXY, MexCD-OprJ, and MexEF-OprN. In this review, the function, structure, and regulation of these efflux pumps in P. aeruginosa and their actions as resistance mechanisms are discussed. Finally, a brief discussion on the potential of efflux pumps in P. aeruginosa as a target for new drugs is presented.

Amid the SARS-CoV-2 continuously changing epidemic profile, this study details the space-time dynamics of the emergence of the Delta lineage across Brazilian territories, pointing out its multiple introductions in the country and its most prevalent sublineages. Some of these sublineages have their emergence, alongside their genomic composition and geographic distribution, detailed here for the first time.

Antimicrobial resistance is a public health burden with worldwide impacts and was recently identified as one of the major causes of death in 2019. Fosfomycin is an antibiotic commonly used to treat urinary tract infections, and resistance to it in Enterobacteriaceae is mainly due to the metalloenzyme FosA3 encoded by the fosA3 gene. In this work, we adapted a CRISPR-Cas9 system named pRE-FOSA3 to restore the sensitivity of a fosA3+ Escherichia coli strain. The fosA3+ E. coli strain was generated by transforming synthetic fosA3 into a nonpathogenic E. coli TOP10. To mediate the fosA3 disruption, two guide RNAs (gRNAs) were selected that used conserved regions within the fosA3 sequence of more than 700 fosA3+ E. coli isolates, and the resensitization plasmid pRE-FOSA3 was assembled by cloning the gRNA into pCas9. gRNA_195 exhibited 100% efficiency in resensitizing the bacteria to fosfomycin. Additionally, the edited strain lost the ampicillin resistance encoded in the same plasmid containing the synthetic fosA3 gene, despite not being the CRISPR-Cas9 target, indicating plasmid clearance. The in vitro analysis presented here points to a path that can be explored to assist the development of effective alternative methods of treatment against fosA3+ bacteria.

The genus Herbaspirillum is widely studied for its ability to associate with grasses and to perform biological nitrogen fixation. However, the bacteria of the Herbaspirillum genus have frequently been isolated from clinical samples. Understanding the genomic characteristics that allow these bacteria to switch environments and become able to colonize human hosts is essential for monitoring emerging pathogens and predicting outbreaks. In this work, we describe the sequencing, assembly, and annotation of the genome of H. frisingense AU14559 isolated from the sputum of patients with cystic fibrosis, and its comparison with the genomes of the uropathogenic strain VT-16–41 and the environmental strains GSF30, BH-1, IAC152, and SG826. The genes responsible for biological nitrogen fixation were absent from all strains except for GSF30. On the other hand, genes encoding virulence and host interaction factors were mostly shared with environmental strains. We also identified a large set of intrinsic antibiotic resistance genes that were shared across all strains. Unlike other strains, in addition to unique genomic islands, AU14559 has a mutation that renders the biosynthesis of rhamnose and its incorporation into the exopolysaccharide unfeasible. These data suggest that H. frisingense has characteristics that provide it with the metabolic diversity needed to infect and colonize human hosts.

Abstract Background We report a genomic surveillance of SARS-CoV-2 lineages circulating in Paraná, southern Brazil, from March 2020 to April 2021. Our analysis, based on 333 genomes, revealed that the first variants detected in the state of Paraná in March 2020 were the B.1.1.33 and B.1.1.28 variants. The variants B.1.1.28 and B.1.1.33 were predominant throughout 2020 until the introduction of the variant P.2 in August 2020 and a variant of concern (VOC), Gamma (P.1), in January 2021. The VOC Gamma, a ramification of the B.1.1.28 lineage first detected in Manaus (northern Brazil), has grown rapidly since December 2020 and was thought to be responsible for the deadly second wave of COVID-19 throughout Brazil. Methods The 333 genomic sequences of SARS-CoV-2 from March 2020 to April 2021 were generated as part of the genomic surveillance carried out by Fiocruz in Brazil Genomahcov Fiocruz. SARS-CoV-2 sequencing was performed using representative samples from all geographic areas of Paraná. Phylogenetic analyses were performed using the 333 genomes also included other SARS-CoV-2 genomes from the state of Paraná and other states in Brazil that were deposited in the GISAID. In addition, the time-scaled phylogenetic tree was constructed with up to 3 random sequences of the Gamma variant from each state in Brazil in each month of 2021. In this analysis we also added the sequences identified as the B.1.1.28 lineage of the Amazonas state and and the Gamma-like-II (P.1-like-II) lineage identified in different regions of Brazil. Results Phylogenetic analyses of the SARS-CoV-2 genomes that were previously classified as the VOC Gamma lineage by WHO/PANGO showed that some genomes from February to April 2021 branched in a monophyletic clade and that these samples grouped together with genomes recently described with the lineage Gamma-like-II. Additionally, a new mutation (E661D) in the spike (S) protein has been identified in nearly 10% of the genomes classified as the VOC Gamma from Paraná in March and April 2021.Finally, we analyzed the correlation between the lineage and the Gamma variant frequency, age group (patients younger or older than 60 years old) and the clinical data of 86 cases from the state of Paraná. Conclusions Our results provided a reliable picture of the evolution of the SARS-CoV-2 pandemic in the state of Paraná characterized by the dominance of the Gamma strain, as well as a high frequencies of the Gamma-like-II lineage and the S:E661D mutation. Epidemiological and genomic surveillance efforts should be continued to unveil the biological relevance of the novel mutations detected in the VOC Gamma in Paraná.

We report a genomic surveillance of SARS-CoV-2 lineages circulating in Paraná, Southern Brazil, from March 2020 to April 2021. Our analysis, based on 333 genomes, revealed that the first variants detected in the state of Paraná in March 2020 were the B.1.1.33 and B.1.1.28 variants. The variants B.1.1.28 and B.1.1.33 were predominant throughout 2020 until the introduction of the variant P.2 in August 2020 and a variant of concern (VOC), P.1, in January 2021. Phylogenetic analyses of the SARS-CoV-2 genomes that were previously classified as the VOC P.1 lineage by PANGO showed that some genomes from February to April 2021 branched in a monophyletic clade and that these samples grouped together with genomes recently described with the lineage P.1-like-II. An extended phylogenetic analysis, including SARS-CoV-2 genomes from all over Brazil, showed that the P.1-like-II lineage appears at a high frequency in the southern region of the country. The P.1-like-II lineage genomes share some, but not all, defining mutations of the VOC P.1. For instance, it has the previously described ORF1a:D2980H and N:P383 L unique mutations and the newly detected ORF1a:P1213 L and ORF1b:K2340N mutations. Additionally, a new mutation (E661D) in the spike (S) protein has been identified in nearly 10% of the genomes classified as the VOC P.1 from Paraná in March and April 2021. We also report the identification of the S:W152C mutation in one genome from Paraná, classified as the N.10 variant. Finally, we analyzed the correlation between the lineage and the P.1 variant frequency, age group (patients younger or older than 60 years old) and the clinical data of 86 cases from the state of Paraná. This analysis does not support an association between the P.1 variant prevalence and COVID-19 severity or age strata. Our results provided a reliable picture of the evolution of the SARS-CoV-2 pandemic in the state of Paraná characterized by the dominance of the P.1 strain, as well as a high frequencies of the P.1-like-II lineage and the S:E661D mutations. Epidemiological and genomic surveillance efforts should be continued to unveil the biological relevance of the novel mutations detected in the VOC P.1 in Paraná.

Nonulosonic acids (NulOs) are a group of nine-carbon monosaccharides with different functions in nature. N-acetylneuraminic acid (Neu5Ac) is the most common NulO. It covers the membrane surface of all human cells and is a central molecule in the process of self-recognition via SIGLECS receptors. Some pathogenic bacteria escape the immune system by copying the sialylation of the host cell membrane. Neu5Ac production in these bacteria is catalysed by the enzyme NeuB. Some bacteria can also produce other NulOs named pseudaminic and legionaminic acids, through the NeuB homologues PseI and LegI, respectively. In Opisthokonta eukaryotes, the biosynthesis of Neu5Ac is catalysed by the enzyme NanS. In this study, we used publicly available data of sequences of NulOs synthases to investigate its distribution within the three domains of life and its relationship with pathogenic bacteria. We mined the KEGG database and found 425 NeuB sequences. Most NeuB sequences (58.74 %) from the KEGG orthology database were classified as from environmental bacteria; however, sequences from pathogenic bacteria showed higher conservation and prevalence of a specific domain named SAF. Using the HMM profile we identified 13 941 NulO synthase sequences in UniProt. Phylogenetic analysis of these sequences showed that the synthases were divided into three main groups that can be related to the lifestyle of these bacteria: (I) predominantly environmental, (II) intermediate and (III) predominantly pathogenic. NeuB was widely distributed in the groups. However, LegI and PseI were more concentrated in groups II and III, respectively. We also found that PseI appeared later in the evolutionary process, derived from NeuB. We use this same methodology to retrieve sialic acid synthase sequences from Archaea and Eukarya. A large-scale phylogenetic analysis showed that while the Archaea sequences are spread across the tree, the eukaryotic NanS sequences were grouped in a specific branch in group II. None of the bacterial NanS sequences grouped with the eukaryotic branch. The analysis of conserved residues showed that the synthases of Archaea and Eukarya present a mutation in one of the three catalytic residues, an E134D change, related to a Neisseria meningitidis reference sequence. We also found that the conservation profile is higher between NeuB of pathogenic bacteria and NanS of eukaryotes than between NeuB of environmental bacteria and NanS of eukaryotes. Our large-scale analysis brings new perspectives on the evolution of NulOs synthases, suggesting their presence in the last common universal ancestor.