Metabolic network analysis of pre-ASD newborns and 5-year-old children with autism spectrum disorder.

Classical metabolomic and new metabolic network methods were used to study the developmental features of autism spectrum disorder (ASD) in newborns (n = 205) and 5-year-old children (n = 53). Eighty percent of the metabolic impact in ASD was caused by 14 shared biochemical pathways that led to decreased anti-inflammatory and antioxidant defenses, and to increased physiologic stress molecules like lactate, glycerol, cholesterol, and ceramides. CIRCOS plots and a new metabolic network parameter, V ° net, revealed differences in both the kind and degree of network connectivity. Of 50 biochemical pathways and 450 polar and lipid metabolites examined, the developmental regulation of the purine network was most changed. Purine network hub analysis revealed a 17-fold reversal in typically developing children. This purine network reversal did not occur in ASD. These results revealed previously unknown metabolic phenotypes, identified new developmental states of the metabolic correlation network, and underscored the role of mitochondrial functional changes, purine metabolism, and purinergic signaling in autism spectrum disorder.

A validated pangenome-scale metabolic model for the Klebsiella pneumoniae species complex.

The Klebsiella pneumoniae species complex (KpSC) is a major source of nosocomial infections globally with high rates of resistance to antimicrobials. Consequently, there is growing interest in understanding virulence factors and their association with cellular metabolic processes for developing novel anti-KpSC therapeutics. Phenotypic assays have revealed metabolic diversity within the KpSC, but metabolism research has been neglected due to experiments being difficult and cost-intensive. Genome-scale metabolic models (GSMMs) represent a rapid and scalable in silico approach for exploring metabolic diversity, which compile genomic and biochemical data to reconstruct the metabolic network of an organism. Here we use a diverse collection of 507 KpSC isolates, including representatives of globally distributed clinically relevant lineages, to construct the most comprehensive KpSC pan-metabolic model to date, KpSC pan v2. Candidate metabolic reactions were identified using gene orthology to known metabolic genes, prior to manual curation via extensive literature and database searches. The final model comprised a total of 3550 reactions, 2403 genes and can simulate growth on 360 unique substrates. We used KpSC pan v2 as a reference to derive strain-specific GSMMs for all 507 KpSC isolates, and compared these to GSMMs generated using a prior KpSC pan-reference (KpSC pan v1) and two single-strain references. We show that KpSC pan v2 includes a greater proportion of accessory reactions (8.8 %) than KpSC pan v1 (2.5 %). GSMMs derived from KpSC pan v2 also generate more accurate growth predictions, with high median accuracies of 95.4 % (aerobic, n=37 isolates) and 78.8 % (anaerobic, n=36 isolates) for 124 matched carbon substrates. KpSC pan v2 is freely available at https://github.com/kelwyres/KpSC-pan-metabolic-model, representing a valuable resource for the scientific community, both as a source of curated metabolic information and as a reference to derive accurate strain-specific GSMMs. The latter can be used to investigate the relationship between KpSC metabolism and traits of interest, such as reservoirs, epidemiology, drug resistance or virulence, and ultimately to inform novel KpSC control strategies.

Metabolic features of treatment-refractory major depressive disorder with suicidal ideation.

Peripheral blood metabolomics was used to gain chemical insight into the biology of treatment-refractory Major Depressive Disorder with suicidal ideation, and to identify individualized differences for personalized care. The study cohort consisted of 99 patients with treatment-refractory major depressive disorder and suicidal ideation (trMDD-SI n = 52 females and 47 males) and 94 age- and sex-matched healthy controls (n = 48 females and 46 males). The median age was 29 years (IQR 22-42). Targeted, broad-spectrum metabolomics measured 448 metabolites. Fibroblast growth factor 21 (FGF21) and growth differentiation factor 15 (GDF15) were measured as biomarkers of mitochondrial dysfunction. The diagnostic accuracy of plasma metabolomics was over 90% (95%CI: 0.80-1.0) by area under the receiver operator characteristic (AUROC) curve analysis. Over 55% of the metabolic impact in males and 75% in females came from abnormalities in lipids. Modified purines and pyrimidines from tRNA, rRNA, and mRNA turnover were increased in the trMDD-SI group. FGF21 was increased in both males and females. Increased lactate, glutamate, and saccharopine, and decreased cystine provided evidence of reductive stress. Seventy-five percent of the metabolomic abnormalities found were individualized. Personalized deficiencies in CoQ10, flavin adenine dinucleotide (FAD), citrulline, lutein, carnitine, or folate were found. Pathways regulated by mitochondrial function dominated the metabolic signature. Peripheral blood metabolomics identified mitochondrial dysfunction and reductive stress as common denominators in suicidal ideation associated with treatment-refractory major depressive disorder. Individualized metabolic differences were found that may help with personalized management.

Global pathogenomic analysis identifies known and candidate genetic antimicrobial resistance determinants in twelve species.

Surveillance programs for managing antimicrobial resistance (AMR) have yielded thousands of genomes suited for data-driven mechanism discovery. We present a workflow integrating pangenomics, gene annotation, and machine learning to identify AMR genes at scale. When applied to 12 species, 27,155 genomes, and 69 drugs, we 1) find AMR gene transfer mostly confined within related species, with 925 genes in multiple species but just eight in multiple phylogenetic classes, 2) demonstrate that discovery-oriented support vector machines outperform contemporary methods at recovering known AMR genes, recovering 263 genes compared to 145 by Pyseer, and 3) identify 142 AMR gene candidates. Validation of two candidates in E. coli BW25113 reveals cases of conditional resistance: ΔcycA confers ciprofloxacin resistance in minimal media with D-serine, and frdD V111D confers ampicillin resistance in the presence of ampC by modifying the overlapping promoter. We expect this approach to be adaptable to other species and phenotypes.

Bactabolize is a tool for high-throughput generation of bacterial strain-specific metabolic models.

Metabolic capacity can vary substantially within a bacterial species, leading to ecological niche separation, as well as differences in virulence and antimicrobial susceptibility. Genome-scale metabolic models are useful tools for studying the metabolic potential of individuals, and with the rapid expansion of genomic sequencing there is a wealth of data that can be leveraged for comparative analysis. However, there exist few tools to construct strain-specific metabolic models at scale. Here, we describe Bactabolize, a reference-based tool which rapidly produces strain-specific metabolic models and growth phenotype predictions. We describe a pan reference model for the priority antimicrobial-resistant pathogen, Klebsiella pneumoniae, and a quality control framework for using draft genome assemblies as input for Bactabolize. The Bactabolize-derived model for K. pneumoniae reference strain KPPR1 performed comparatively or better than currently available automated approaches CarveMe and gapseq across 507 substrate and 2317 knockout mutant growth predictions. Novel draft genomes passing our systematically defined quality control criteria resulted in models with a high degree of completeness (≥99% genes and reactions captured compared to models derived from matched complete genomes) and high accuracy (mean 0.97, n=10). We anticipate the tools and framework described herein will facilitate large-scale metabolic modelling analyses that broaden our understanding of diversity within bacterial species and inform novel control strategies for priority pathogens.

Rev1 deficiency induces a metabolic shift in MEFs that can be manipulated by the NAD+ precursor nicotinamide riboside.

Replication stress, caused by Rev1 deficiency, is associated with mitochondrial dysfunction, and metabolic stress. However, the overall metabolic alterations and possible interventions to rescue the deficits due to Rev1 loss remain unclear. Here, we report that loss of Rev1 leads to intense changes in metabolites and that this can be manipulated by NAD + supplementation. Autophagy decreases in Rev1-/- mouse embryonic fibroblasts (MEFs) and can be restored by supplementing the NAD+ precursor nicotinamide riboside (NR). The abnormal mitochondrial morphology in Rev1-/- MEFs can be partially reversed by NR supplementation, which also protects the mitochondrial cristae from rotenone-induced degeneration. In nematodes rev-1 deficiency causes sensitivity to oxidative stress but this cannot be rescued by NR supplementation. In conclusion, Rev1 deficiency leads to metabolic dysregulation of especially lipid and nucleotide metabolism, impaired autophagy, and mitochondrial anomalies, and all of these phenotypes can be improved by NR replenishment in MEFs.

Comprehensive whole genome sequencing with hybrid assembly of multi-drug resistant Candida albicans isolate causing cerebral abscess.

Comprehensive whole genome sequencing (WGS) with hybrid assembly of a multi-drug resistant (MDR) Candida albicans (CA) isolate causing cerebral abscess was performed using Illumina paired end and Oxford Nanopore long read technologies. The innovative technologies utilized here enabled us to resolve fragmented assemblies, and implement comprehensive and detailed genomic analyses involved in antifungal resistance of Candida spp. Functionally important genes (MDR1, CDR2 and SQN2) involved in antifungal resistance were identified and a phylogenetic analysis of the clinical isolate was performed. Additionally, our clinical isolate was found to share 4 single nucleotide polymorphisms with two other sequenced strains of MDR C. auris (381 and 386) including translation elongation factor EF1α and EF3, ATPase activity associated proteins, and the lysine tRNA ligase.

Historical biomonitoring of pollution trends in the North Pacific using archived samples from the Continuous Plankton Recorder Survey.

First started in 1931, the Continuous Plankton Recorder (CPR) Survey is the longest-running and most geographically extensive marine plankton sampling program in the world. This pilot study investigates the feasibility of biomonitoring the spatiotemporal trends of marine pollution using archived CPR samples from the North Pacific. We selected specimens collected from three different locations (British Columbia Shelf, Northern Gulf of Alaska, and Aleutian Shelf) in the North Pacific between 2002 and 2020. Comprehensive profiling of the plankton chemical exposome was conducted using liquid and gas chromatography coupled with tandem mass spectrometry (LC-MS/MS and GC-MS/MS). Our results show that phthalates, plasticizers, persistent organic pollutants (POPs), pesticides, pharmaceuticals, and personal care products were present in the plankton exposome, and that many of these pollutants have decreased in amount over the last two decades, which was most pronounced for tri-n-butyl phosphate. In addition, the plankton exposome differed significantly by regional human activities, with the most polluted samples coming from the nearshore area. Exposome-wide association analysis revealed that bioaccumulation of environmental pollutants was highly correlated with the biomass of different plankton taxa. Overall, this study demonstrates that exposomic analysis of archived samples from the CPR Survey is effective for long-term biomonitoring of the spatial and temporal trends of environmental pollutants in the marine environment.

Metabolomic and exposomic biomarkers of risk of future neurodevelopmental delay in human milk.

BACKGROUND: The chemical composition of human milk has long-lasting effects on brain development. We examined the prognostic value of the human milk metabolome and exposome in children with the risk of neurodevelopmental delay (NDD). METHODS: This retrospective cohort study included 82 mother-infant pairs (40 male and 42 female infants). A total of 59 milk samples were from mothers with typically developing children and 23 samples were from mothers of children at risk. Milk samples were collected before 9 months of age (4.6 ± 2.5 months, mean ± SD). Neurocognitive development was assessed by maternal report at 14.2 ± 3.1 months using the Ages and Stages Questionnaires-2. RESULTS: Metabolome and exposome profiling identified 453 metabolites and 61 environmental chemicals in milk. Machine learning tools identified changes in deoxysphingolipids, phospholipids, glycosphingolipids, plasmalogens, and acylcarnitines in the milk of mothers with children at risk for future delay. A predictive classifier had a diagnostic accuracy of 0.81 (95% CI: 0.66-0.96) for females and 0.79 (95% CI: 0.62-0.94) for males. CONCLUSIONS: Once validated in larger studies, the chemical analysis of human milk might be added as an option in well-baby checks to help identify children at risk of NDD before the first symptoms appear. IMPACT: Maternal milk for infants sampled before 9 months of age contained sex-specific differences in deoxysphingolipids, sphingomyelins, plasmalogens, phospholipids, and acylcarnitines that predicted the risk of neurodevelopmental delay at 14.2 months of age. Once validated, this early biosignature in human milk might be incorporated into well-baby checks and help to identify infants at risk so early interventions might be instituted before the first symptoms appear.

High-quality genome-scale metabolic network reconstruction of probiotic bacterium Escherichia coli Nissle 1917.

BACKGROUND: Escherichia coli Nissle 1917 (EcN) is a probiotic bacterium used to treat various gastrointestinal diseases. EcN is increasingly being used as a chassis for the engineering of advanced microbiome therapeutics. To aid in future engineering efforts, our aim was to construct an updated metabolic model of EcN with extended secondary metabolite representation. RESULTS: An updated high-quality genome-scale metabolic model of EcN, iHM1533, was developed based on comparison with 55 E. coli/Shigella reference GEMs and manual curation, including expanded secondary metabolite pathways (enterobactin, salmochelins, aerobactin, yersiniabactin, and colibactin). The model was validated and improved using phenotype microarray data, resulting in an 82.3% accuracy in predicting growth phenotypes on various nutrition sources. Flux variability analysis with previously published 13C fluxomics data validated prediction of the internal central carbon fluxes. A standardised test suite called Memote assessed the quality of iHM1533 to have an overall score of 89%. The model was applied by using constraint-based flux analysis to predict targets for optimisation of secondary metabolite production. Modelling predicted design targets from across amino acid metabolism, carbon metabolism, and other subsystems that are common or unique for influencing the production of various secondary metabolites. CONCLUSION: iHM1533 represents a well-annotated metabolic model of EcN with extended secondary metabolite representation. Phenotype characterisation and the iHM1533 model provide a better understanding of the metabolic capabilities of EcN and will help future metabolic engineering efforts.

Experimental Evolution Reveals Unifying Systems-Level Adaptations but Diversity in Driving Genotypes.

Genotype-fitness maps of evolution have been well characterized for biological components, such as RNA and proteins, but remain less clear for systems-level properties, such as those of metabolic and transcriptional regulatory networks. Here, we take multi-omics measurements of 6 different E. coli strains throughout adaptive laboratory evolution (ALE) to maximal growth fitness. The results show the following: (i) convergence in most overall phenotypic measures across all strains, with the notable exception of divergence in NADPH production mechanisms; (ii) conserved transcriptomic adaptations, describing increased expression of growth promoting genes but decreased expression of stress response and structural components; (iii) four groups of regulatory trade-offs underlying the adjustment of transcriptome composition; and (iv) correlates that link causal mutations to systems-level adaptations, including mutation-pathway flux correlates and mutation-transcriptome composition correlates. We thus show that fitness landscapes for ALE can be described with two layers of causation: one based on system-level properties (continuous variables) and the other based on mutations (discrete variables). IMPORTANCE Understanding the mechanisms of microbial adaptation will help combat the evolution of drug-resistant microbes and enable predictive genome design. Although experimental evolution allows us to identify the causal mutations underlying microbial adaptation, it remains unclear how causal mutations enable increased fitness and is often explained in terms of individual components (i.e., enzyme rate) as opposed to biological systems (i.e., pathways). Here, we find that causal mutations in E. coli are linked to systems-level changes in NADPH balance and expression of stress response genes. These systems-level adaptation patterns are conserved across diverse E. coli strains and thus identify cofactor balance and proteome reallocation as dominant constraints governing microbial adaptation.

Genome-scale metabolic network reconstructions of diverse Escherichia strains reveal strain-specific adaptations.

Bottom-up approaches to systems biology rely on constructing a mechanistic basis for the biochemical and genetic processes that underlie cellular functions. Genome-scale network reconstructions of metabolism are built from all known metabolic reactions and metabolic genes in a target organism. A network reconstruction can be converted into a mathematical format and thus lend itself to mathematical analysis. Genome-scale models (GEMs) of metabolism enable a systems approach to characterize the pan and core metabolic capabilities of the Escherichia genus. In this work, GEMs were constructed for 222 representative strains of Escherichia across HC1100 levels spanning the known Escherichia phylogeny. The models were used to study Escherichia metabolic diversity and speciation on a large scale. The results show that unique strain-specific metabolic capabilities correspond to different species and nutrient niches. This work is a first step towards a curated reconstruction of pan-Escherichia metabolism. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Estimated Roles of the Carrier and the Bacterial Strain When Methicillin-Resistant Staphylococcus aureus Decolonization Fails: a Case-Control Study.

Methicillin-resistant Staphylococcus aureus (MRSA) is a common bacterial pathogen that frequently colonizes healthy individuals, with potential to cause invasive infection. In Denmark, to keep the prevalence low, MRSA carriers are recommended to undergo decolonization treatments, but achieving decolonization is challenging. Knowledge about the factors contributing to decolonization is scarce. We aimed to identify bacterial genome and clinical factors influencing MRSA decolonization. We identified all new MRSA patients above 2 years of age within the Hvidovre catchment area, Copenhagen, Denmark, in 2017 and 2018. Carriers were defined as chronic carriers (cases) if they were MRSA positive after two or more treatments and as nonchronic carriers (controls) if they were MRSA free after the first or second treatment. Using whole-genome sequencing (WGS), we constructed a pangenome of bacterial strains. With the incorporation of bacterial genome and clinical patient data, machine learning and multivariate analyses were performed to determine the factors associated with decolonization. A total of 477 MRSA carriers were included. An age of ≥13 years was significantly associated with nonchronic carriage. We identified 278 bacterial genetic features that were statistically significantly associated with chronic carriage (P < 0.05 by Fisher's exact test). Chronic MRSA carriage was predicted with 68% accuracy using a combination of bacterial genome data and patient clinical data. Decolonization success is multifactorial. Apart from the 68% predicted accuracy found in this study, we estimate that the remaining 32% is a result of host factors and microbiome composition. IMPORTANCE Carriage of methicillin-resistant Staphylococcus aureus (MRSA) and other multiresistant bacteria is a prerequisite for infection and transmission. Successful decolonization treatment removes these risks. We aimed to identify bacterial genome and host clinical factors that influence MRSA decolonization to estimate the roles of the carrier and the bacterial strain, respectively, when decolonization fails. The long-term goal, beyond this study, is to optimize decolonization success, minimize MRSA transmission, and, ultimately, improve the quality of life of MRSA carriers.

Pangenome analysis of Enterobacteria reveals richness of secondary metabolite gene clusters and their associated gene sets.

In silico genome mining provides easy access to secondary metabolite biosynthetic gene clusters (BGCs) encoding the biosynthesis of many bioactive compounds, which are the basis for many important drugs used in human medicine. However, the association between BGCs and other functions encoded in the genomes of producers have remained elusive. Here, we present a systems biology workflow that integrates genome mining with a detailed pangenome analysis for detecting genes associated with a particular BGC. We analyzed 3,889 enterobacterial genomes and found 13,266 BGCs, represented by 252 distinct BGC families and 347 additional singletons. A pangenome analysis revealed 88 genes putatively associated with a specific BGC coding for the colon cancer-related colibactin that code for diverse metabolic and regulatory functions. The presented workflow opens up the possibility to discover novel secondary metabolites, better understand their physiological roles, and provides a guide to identify and analyze BGC associated gene sets.

Cerebrospinal fluid and plasma metabolomics of acute endurance exercise.

Metabolomics has emerged as a powerful new tool in precision medicine. No studies have yet been published on the metabolomic changes in cerebrospinal fluid (CSF) produced by acute endurance exercise. CSF and plasma were collected from 19 young active adults (13 males and 6 females) before and 60 min after a 90-min monitored outdoor run. The median age, BMI, and VO2 max of subjects was 25 years (IQR 22-31), 23.2 kg/m2 (IQR 21.7-24.5), and 47 ml/kg/min (IQR 38-51), respectively. Targeted, broad-spectrum metabolomics was performed by liquid chromatography, tandem mass spectrometry (LC-MS/MS). In the CSF, purines and pyrimidines accounted for 32% of the metabolic impact after acute endurance exercise. Branch chain amino acids, amino acid neurotransmitters, fatty acid oxidation, phospholipids, and Krebs cycle metabolites traceable to mitochondrial function accounted for another 52% of the changes. A narrow but important channel of metabolic communication was identified between the brain and body by correlation network analysis. By comparing these results to previous work in experimental animal models, we found that over 80% of the changes in the CSF correlated with a cascade of mitochondrial and metabolic changes produced by ATP signaling. ATP is released as a co-neurotransmitter and neuromodulator at every synapse studied to date. By regulating brain mitochondrial function, ATP release was identified as an early step in the kinetic cascade of layered benefits produced by endurance exercise.

Transmission of Klebsiella strains and plasmids within and between grey-headed flying fox colonies.

The grey-headed flying fox (Pteropus poliocephalus) is an endemic Australian fruit bat, known to carry zoonotic pathogens. We recently showed they harbour bacterial pathogen Klebsiella pneumoniae and closely related species in the K. pneumoniae species complex (KpSC); however, the dynamics of KpSC transmission and gene flow within flying fox colonies are poorly understood. High-resolution genome comparisons of 39 KpSC isolates from grey-headed flying foxes identified five putative strain transmission clusters (four intra- and one inter-colony). The instance of inter-colony strain transmission of K. africana was found between two flying fox populations within flying distance, indicating either direct or indirect transmission through a common food/water source. All 11 plasmids identified within the KpSC isolates showed 73% coverage (mean) and ≥95% identity to human-associated KpSC plasmids, indicating gene flow between human clinical and grey-headed flying fox isolates. Along with strain transmission, inter-species horizontal plasmid transmission between K. pneumoniae and Klebsiella africana was also identified within a flying fox colony. Finally, genome-scale metabolic models were generated to predict and compare substrate usage to previously published KpSC models, from human and environmental sources. These models indicated no distinction on the basis of metabolic capabilities. Instead, metabolic capabilities were consistent with population structure and ST/lineage.

Systems biology approach to functionally assess the Clostridioides difficile pangenome reveals genetic diversity with discriminatory power.

Combatting Clostridioides difficile infections, a dominant cause of hospital-associated infections with incidence and resulting deaths increasing worldwide, is complicated by the frequent emergence of new virulent strains. Here, we employ whole-genome sequencing, high-throughput phenotypic screenings, and genome-scale models of metabolism to evaluate the genetic diversity of 451 strains of C. difficile. Constructing the C. difficile pangenome based on this set revealed 9,924 distinct gene clusters, of which 2,899 (29%) are defined as core, 2,968 (30%) are defined as unique, and the remaining 4,057 (41%) are defined as accessory. We develop a strain typing method, sequence typing by accessory genome (STAG), that identifies 176 genetically distinct groups of strains and allows for explicit interrogation of accessory gene content. Thirty-five strains representative of the overall set were experimentally profiled on 95 different nutrient sources, revealing 26 distinct growth profiles and unique nutrient preferences; 451 strain-specific genome scale models of metabolism were constructed, allowing us to computationally probe phenotypic diversity in 28,864 unique conditions. The models create a mechanistic link between the observed phenotypes and strain-specific genetic differences and exhibit an ability to correctly predict growth in 76% of measured cases. The typing and model predictions are used to identify and contextualize discriminating genetic features and phenotypes that may contribute to the emergence of new problematic strains.

Mathematical models to study the biology of pathogens and the infectious diseases they cause.

Mathematical models have many applications in infectious diseases: epidemiologists use them to forecast outbreaks and design containment strategies; systems biologists use them to study complex processes sustaining pathogens, from the metabolic networks empowering microbial cells to ecological networks in the microbiome that protects its host. Here, we (1) review important models relevant to infectious diseases, (2) draw parallels among models ranging widely in scale. We end by discussing a minimal set of information for a model to promote its use by others and to enable predictions that help us better fight pathogens and the diseases they cause.

A curated collection of Klebsiella metabolic models reveals variable substrate usage and gene essentiality.

The Klebsiella pneumoniae species complex (KpSC) is a set of seven Klebsiella taxa that are found in a variety of niches and are an important cause of opportunistic health care-associated infections in humans. Because of increasing rates of multi-drug resistance within the KpSC, there is a growing interest in better understanding the biology and metabolism of these organisms to inform novel control strategies. We collated 37 sequenced KpSC isolates isolated from a variety of niches, representing all seven taxa. We generated strain-specific genome-scale metabolic models (GEMs) for all 37 isolates and simulated growth phenotypes on 511 distinct carbon, nitrogen, sulfur, and phosphorus substrates. Models were curated and their accuracy was assessed using matched phenotypic growth data for 94 substrates (median accuracy of 96%). We explored species-specific growth capabilities and examined the impact of all possible single gene deletions using growth simulations in 145 core carbon substrates. These analyses revealed multiple strain-specific differences, within and between species, and highlight the importance of selecting a diverse range of strains when exploring KpSC metabolism. This diverse set of highly accurate GEMs could be used to inform novel drug design, enhance genomic analyses, and identify novel virulence and resistance determinants. We envisage that these 37 curated strain-specific GEMs, covering all seven taxa of the KpSC, provide a valuable resource to the Klebsiella research community.

Comparative pangenomics: analysis of 12 microbial pathogen pangenomes reveals conserved global structures of genetic and functional diversity.

BACKGROUND: With the exponential growth of publicly available genome sequences, pangenome analyses have provided increasingly complete pictures of genetic diversity for many microbial species. However, relatively few studies have scaled beyond single pangenomes to compare global genetic diversity both within and across different species. We present here several methods for "comparative pangenomics" that can be used to contextualize multi-pangenome scale genetic diversity with gene function for multiple species at multiple resolutions: pangenome shape, genes, sequence variants, and positions within variants. RESULTS: Applied to 12,676 genomes across 12 microbial pathogenic species, we observed several shared resolution-specific patterns of genetic diversity: First, pangenome openness is associated with species' phylogenetic placement. Second, relationships between gene function and frequency are conserved across species, with core genomes enriched for metabolic and ribosomal genes and accessory genomes for trafficking, secretion, and defense-associated genes. Third, genes in core genomes with the highest sequence diversity are functionally diverse. Finally, certain protein domains are consistently mutation enriched across multiple species, especially among aminoacyl-tRNA synthetases where the extent of a domain's mutation enrichment is strongly function-dependent. CONCLUSIONS: These results illustrate the value of each resolution at uncovering distinct aspects in the relationship between genetic and functional diversity across multiple species. With the continued growth of the number of sequenced genomes, these methods will reveal additional universal patterns of genetic diversity at the pangenome scale.