Genomic language model predicts protein co-regulation and function

Deciphering the relationship between a gene and its genomic context is fundamental to understanding and engineering biological systems. Machine learning has shown promise in learning latent relationships underlying the sequence-structure-function paradigm from massive protein sequence datasets. However, to date, limited attempts have been made in extending this continuum to include higher order genomic context information. Evolutionary processes dictate the specificity of genomic contexts in which a gene is found across phylogenetic distances, and these emergent genomic patterns can be leveraged to uncover functional relationships between gene products. Here, we train a genomic language model (gLM) on millions of metagenomic scaffolds to learn the latent functional and regulatory relationships between genes. gLM learns contextualized protein embeddings that capture the genomic context as well as the protein sequence itself, and encode biologically meaningful and functionally relevant information (e.g. enzymatic function, taxonomy). Our analysis of the attention patterns demonstrates that gLM is learning co-regulated functional modules (i.e. operons). Our findings illustrate that gLM's unsupervised deep learning of the metagenomic corpus is an effective and promising approach to encode functional semantics and regulatory syntax of genes in their genomic contexts and uncover complex relationships between genes in a genomic region.

A genetic screen identifies a role for oprF in Pseudomonas aeruginosa biofilm stimulation by subinhibitory antibiotics

Biofilms are surface-associated communities of bacteria that grow in a self-produced matrix of polysaccharides, proteins, and extracellular DNA (eDNA). Sub-minimal inhibitory concentrations (sub-MIC) of antibiotics induce biofilm formation, potentially as a defensive response to antibiotic stress. However, the mechanisms behind sub-MIC antibiotic-induced biofilm formation are unclear. We show that treatment of Pseudomonas aeruginosa with multiple classes of sub-MIC antibiotics with distinct targets induces biofilm formation. Further, addition of exogenous eDNA or cell lysate failed to increase biofilm formation to the same extent as antibiotics, suggesting that the release of cellular contents by antibiotic-driven bacteriolysis is insufficient. Using a genetic screen for stimulation-deficient mutants, we identified the outer membrane porin OprF and the ECF sigma factor SigX as important. Similarly, loss of OmpA - the Escherichia coli OprF homolog - prevented sub-MIC antibiotic stimulation of E. coli biofilms. Our screen also identified the periplasmic disulfide bond-forming enzyme DsbA and a predicted cyclic-di-GMP phosphodiesterase encoded by PA2200 as essential for biofilm stimulation. The phosphodiesterase activity of PA2200 is likely controlled by a disulfide bond in its regulatory domain, and folding of OprF is influenced by disulfide bond formation, connecting the mutant phenotypes. Addition of reducing agent dithiothreitol prevented sub-MIC antibiotic biofilm stimulation. Finally, activation of a c-di-GMP-responsive promoter follows treatment with sub-MIC antibiotics in the wild-type but not an oprF mutant. Together, these results show that antibiotic-induced biofilm formation is likely driven by a signaling pathway that translates changes in periplasmic redox state into elevated biofilm formation through increases in c-di-GMP.

Identification of microbial pathogens in Neolithic Scandinavian humans

With the Neolithic transition, human lifestyle shifted from hunting and gathering to farming. This change altered subsistence patterns, cultural expression, and population structures as shown by the archaeological/zooarchaeological record, as well as by stable isotope and ancient DNA data. Here, we used metagenomic data to analyse if the transitions also impacted the microbiome composition in 25 Mesolithic and Neolithic hunter-gatherers and 13 Neolithic farmers from several Scandinavian Stone Age cultural contexts. Salmonella enterica, a bacterium that may have been the cause of death for the infected individuals, was found in two Neolithic samples from Battle Axe culture contexts. Several species of the bacterial genus Yersinia were found in Neolithic individuals from Funnel Beaker culture contexts as well as from later Neolithic context. Transmission of e.g. Y. enterocolitica may have been facilitated by the denser populations in agricultural contexts.

Occurrence and mechanisms of tigecycline resistance in carbapenem- and colistin-resistant Klebsiella pneumoniae in Thailand

Tigecycline has been regarded as one of the most important last-resort antibiotics for the treatment of infections caused by extensively drug-resistant (XDR) bacteria, particularly carbapenem- and colistin-resistant Klebsiella pneumoniae (C-C-RKP). However, reports on tigecycline resistance have been growing. Overall, ~ 4000 K. pneumoniae clinical isolates were collected over a five-year period (2017-2021), in which 240 isolates of C-C-RKP were investigated. Most of these isolates (91.7%) were resistant to tigecycline. Notably, a high-risk clone of ST16 was predominantly identified, which was associated with the co-harboring of blaNDM-1 and blaOXA-232 genes. Their major mechanism of tigecycline resistance was the overexpression of efflux pump acrB gene and its regulator RamA, which was caused by mutations in RamR (M184V, Y59C, I141T, A28T, C99/C100 insertion), in RamR binding site (PI) of ramA gene (C139T), in MarR (S82G), and/or in AcrR (L154R, R13Q). Interestingly, four isolates of ST147 carried the mutated tet(A) efflux pump gene. To our knowledge, this is the first report on the prevalence and mechanisms of tigecycline resistance in C-C-RKP isolated from Thailand. The high incidence of tigecycline resistance observed among C-C-RKP in this study reflects an ongoing evolution of XDR bacteria against the last-resort antibiotics, which demands urgent action.

A genome-wide association study reveals the relationship between human genetic variation and the nasal microbiome

The nasal cavity harbors diverse microbiota that contributes to human health and respiratory diseases. However, whether and to what extent the host genome shapes the nasal microbiome remains largely unknown. Here, by dissecting the human genome and nasal metagenome data from 1401 healthy individuals, we demonstrated that the top three host genetic principal components strongly correlated with the nasal microbiota diversity and composition. The genetic association analyses identified 63 genome-wide significant loci affecting the nasal microbial taxa and functions, of which 2 loci reached study-wide significance (p < 1.7 × 10-10): rs73268759 within CAMK2A associated with genus Actinomyces and family Actinomycetaceae; and rs35211877 near POM121L12 with Gemella asaccharolytica. In addition to respiratory-related diseases, the associated loci are mainly implicated in cardiometabolic or neuropsychiatric diseases. Functional analysis showed the associated genes were most significantly expressed in the nasal airway epithelium tissue and enriched in the calcium signaling and hippo signaling pathway. Further observational correlation and Mendelian randomization analyses consistently suggested the causal effects of Serratia grimesii and Yokenella regensburgei on cardiometabolic biomarkers (cystine, glutamic acid, and creatine). This study suggested that the host genome plays an important role in shaping the nasal microbiome.

Mendelian randomization study revealed a gut microbiota-neuromuscular junction axis in myasthenia gravis

A growing number of studies have implicated that gut microbiota abundance is associated with myasthenia gravis (MG). However, the causal relationship underlying the associations is still unclear. Here, we aim to investigate the causal effect of gut microbiota on MG using Mendelian randomization (MR) method. Publicly available Genome-wide association study (GWAS) summary-level data for gut microbiota and for MG were extracted. Inverse variance weighted was used as the main method to analyze causality. The robustness of the results was validated with sensitivity analyses. Our results indicated that genetically predicted increased phylum Lentisphaerae (OR = 1.319, p = 0.026), class Lentisphaerae (OR = 1.306, p = 0.044), order Victivallales (OR = 1.306, p = 0.044), order Mollicutes (OR = 1.424, p = 0.041), and genus Faecalibacterium (OR = 1.763, p = 0.002) were potentially associated with a higher risk of MG; while phylum Actinobacteria (OR = 0.602, p = 0.0124), class Gammaproteobacteria (OR = 0.587, p = 0.036), family Defluviitaleaceae (OR = 0.695, p = 0.047), family Peptococcaceae (OR = 0.698, p = 0.029), and family Family XIII (OR = 0.614, p = 0.017) were related to a lower risk of MG. The present study provides genetic evidence for the causal associations between gut microbiota and MG, thus suggesting novel insights into the gut microbiota-neuromuscular junction axis in the pathogenesis of MG.

BacPE: a versatile prime-editing platform in bacteria by inhibiting DNA exonucleases

Prime editing allows precise installation of any single base substitution and small insertions and deletions without requiring homologous recombination or double-strand DNA breaks in eukaryotic cells. However, the applications in bacteria are hindered and the underlying mechanisms that impede efficient prime editing remain enigmatic. Here, we report the determination of vital cellular factors that affect prime editing in bacteria. Genetic screening of 129 Escherichia coli transposon mutants identified sbcB, a 3'→5' DNA exonuclease, as a key genetic determinant in impeding prime editing in E. coli, combinational deletions of which with two additional 3'→5' DNA exonucleases, xseA and exoX, drastically enhanced the prime editing efficiency by up to 100-fold. Efficient prime editing in wild-type E. coli can be achieved by simultaneously inhibiting the DNA exonucleases via CRISPRi. Our results pave the way for versatile applications of prime editing for bacterial genome engineering.

Investigating the etiologies of non-malarial febrile illness in Senegal using metagenomic sequencing

The worldwide decline in malaria incidence is revealing the extensive burden of non-malarial febrile illness (NMFI), which remains poorly understood and difficult to diagnose. To characterize NMFI in Senegal, we collected venous blood and clinical metadata in a cross-sectional study of febrile patients and healthy controls in a low malaria burden area. Using 16S and untargeted sequencing, we detected viral, bacterial, or eukaryotic pathogens in 23% (38/163) of NMFI cases. Bacteria were the most common, with relapsing fever Borrelia and spotted fever Rickettsia found in 15.5% and 3.8% of cases, respectively. Four viral pathogens were found in a total of 7 febrile cases (3.5%). Sequencing also detected undiagnosed Plasmodium, including one putative P. ovale infection. We developed a logistic regression model that can distinguish Borrelia from NMFIs with similar presentation based on symptoms and vital signs (F1 score: 0.823). These results highlight the challenge and importance of improved diagnostics, especially for Borrelia, to support diagnosis and surveillance.

Comparison of the DNBSEQ platform and Illumina HiSeq 2000 for bacterial genome assembly

The Illumina HiSeq platform has been a commonly used option for bacterial genome sequencing. Now the BGI DNA nanoball (DNB) nanoarrays platform may provide an alternative platform for sequencing of bacterial genomes. To explore the impact of sequencing platforms on bacterial genome assembly, quality assessment, sequence alignment, functional annotation, mutation detection, and metagenome mapping, we compared genome assemblies based on sequencing of cultured bacterial species using the HiSeq 2000 and BGISEQ-500 platforms. In addition, simulated reads were used to evaluate the impact of insert size on genome assembly. Genome assemblies based on BGISEQ-500 sequencing exhibited higher completeness and fewer N bases in high GC genomes, whereas HiSeq 2000 assemblies exhibited higher N50. The majority of assembly assessment parameters, sequences of 16S rRNA genes and genomes, numbers of single nucleotide variants (SNV), and mapping to metagenome data did not differ significantly between platforms. More insertions were detected in HiSeq 2000 genome assemblies, whereas more deletions were detected in BGISEQ-500 genome assemblies. Insert size had no significant impact on genome assembly. Taken together, our results suggest that DNBSEQ platforms would be a valid substitute for HiSeq 2000 for bacterial genome sequencing.

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

Unraveling the genome of Bacillus velezensis MEP218, a strain producing fengycin homologs with broad antibacterial activity: comprehensive comparative genome analysis

Bacillus sp. MEP218, a soil bacterium with high potential as a source of bioactive molecules, produces mostly C16-C17 fengycin and other cyclic lipopeptides (CLP) when growing under previously optimized culture conditions. This work addressed the elucidation of the genome sequence of MEP218 and its taxonomic classification. The genome comprises 3,944,892 bp, with a total of 3474 coding sequences and a G + C content of 46.59%. Our phylogenetic analysis to determine the taxonomic position demonstrated that the assignment of the MEP218 strain to Bacillus velezensis species provides insights into its evolutionary context and potential functional attributes. The in silico genome analysis revealed eleven gene clusters involved in the synthesis of secondary metabolites, including non-ribosomal CLP (fengycins and surfactin), polyketides, terpenes, and bacteriocins. Furthermore, genes encoding phytase, involved in the release of phytic phosphate for plant and animal nutrition, or other enzymes such as cellulase, xylanase, and alpha 1-4 glucanase were detected. In vitro antagonistic assays against Salmonella typhimurium, Acinetobacter baumanii, Escherichia coli, among others, demonstrated a broad spectrum of C16-C17 fengycin produced by MEP218. MEP218 genome sequence analysis expanded our understanding of the diversity and genetic relationships within the Bacillus genus and updated the Bacillus databases with its unique trait to produce antibacterial fengycins and its potential as a resource of biotechnologically useful enzymes.

Genotyping and phylogenetic analysis of Coxiella burnetii in domestic ruminant and clinical samples in Iran: insights into Q fever epidemiology

Coxiella burnetii, a zoonotic pathogen, is the causative agent of Q fever, an endemic disease in Iran. However, there is currently a lack of available data on the genotypes of C. burnetii in the country. Here, we typed 26 C. burnetii isolates detected in milk, abortion, cotylodon, and cardiac valve samples from various geographical areas and hosts (7 cattle, 8 goats, 10 sheep, and 1 human) using Multilocus Variable Number Tandem Repeat Analysis (MLVA/VNTR) with five loci:ms24, ms27, ms28, ms33, and ms34. As IS1111 was observed to be spontaneously inserted in locus ms23 across all of our examined C. burnetii samples, five loci were employed for MLVA/VNTR genotyping. Among the 26 C. burnetii strains, 22 distinct genotypes (A-V) were identified in the discriminative loci. In silico analysis categorized Iranian C. burnetii strains into five genomic groups along with seven singletons, representing 11 exiting clonal complexes worldwide. Clusters 10 and 11 exclusively consisted of Iranian samples. These findings revealed high genotyping diversity among C. burnetii isolates in Iran. The genotypes circulating in Iran differed significantly from those found in other regions worldwide. To gain a comprehensive understanding of Q fever epidemiology in Iran, it is crucial to conduct large-scale studies that assess the distribution of C. burnetii genotypes across different geographical areas, hosts, and sources.

Investigation of antimicrobial susceptibility and genetic diversity among Staphylococcus pseudintermedius isolated from dogs in Rio de Janeiro

Staphylococcus pseudintermedius is an opportunistic pathogen causing a variety of infections that are difficult to treat, especially because of the development of antimicrobial resistance. It has a clonal distribution around the world. To have a better understanding of the MRSP population, we search the presence of MRSP in colonized or infected dogs. Samples from 99 dogs with infections and 35 from asymptomatic dogs were collected. Isolates were identified by mass spectrometry and Multiplex-PCR. The mecA gene was confirmed by conventional PCR. MRSP strains were analyzed by whole-genome sequencing. 75 S. pseudintermedius were identified, most from infection cases. The species were isolated from 70 out of the 135 dogs. Penicillin and Trimethoprim/Sulfamethoxazole presented higher resistance rates. Forty-seven strains were classified as multi-drug resistant (MDR), and were more isolated from dogs with infection (P < 0.05). Eighteen samples were classified as MRSP, representing 24.0% of the population. Six of 16 MRSP sequenced samples belonged to the world spread clone ST71; others belonged to unknown clones. Most samples carried the SCCmec type IIIA. Twenty-one different genetic resistance determinants were found among MRPS strains. MRSP is circulating among infected and colonized dogs in Rio de Janeiro, Brazil.

Synthetic genetic oscillators demonstrate the functional importance of phenotypic variation in pneumococcal-host interactions

Phenotypic variation is the phenomenon in which clonal cells display different traits even under identical environmental conditions. This plasticity is thought to be important for processes including bacterial virulence, but direct evidence for its relevance is often lacking. For instance, variation in capsule production in the human pathogen Streptococcus pneumoniae has been linked to different clinical outcomes, but the exact relationship between variation and pathogenesis is not well understood due to complex natural regulation. In this study, we use synthetic oscillatory gene regulatory networks (GRNs) based on CRISPR interference (CRISPRi) together with live cell imaging and cell tracking within microfluidics devices to mimic and test the biological function of bacterial phenotypic variation. We provide a universally applicable approach for engineering intricate GRNs using only two components: dCas9 and extended sgRNAs (ext-sgRNAs). Our findings demonstrate that variation in capsule production is beneficial for pneumococcal fitness in traits associated with pathogenesis providing conclusive evidence for this longstanding question.

Systematic investigation of recipient cell genetic requirements reveals important surface receptors for conjugative transfer of IncI2 plasmids

Bacterial conjugation is a major horizontal gene transfer mechanism. While the functions encoded by many conjugative plasmids have been intensively studied, the contribution of recipient chromosome-encoded genes remains largely unknown. Here, we analyzed the genetic requirement of recipient cells for conjugation of IncI2 plasmid TP114, which was recently shown to transfer at high rates in the gut microbiota. We performed transfer assays with ~4,000 single-gene deletion mutants of Escherichia coli. When conjugation occurs on a solid medium, we observed that recipient genes impairing transfer rates were not associated with a specific cellular function. Conversely, transfer assays performed in broth were largely dependent on the lipopolysaccharide biosynthesis pathway. We further identified specific structures in lipopolysaccharides used as recipient cell surface receptors by PilV adhesins associated with the type IVb accessory pilus of TP114. Our strategy is applicable to study other mobile genetic elements and understand important host cell factors for their dissemination.

Gut microbial structural variation associates with immune checkpoint inhibitor response

The gut microbiota may have an effect on the therapeutic resistance and toxicity of immune checkpoint inhibitors (ICIs). However, the associations between the highly variable genomes of gut bacteria and the effectiveness of ICIs remain unclear, despite the fact that merely a few gene mutations between similar bacterial strains may cause significant phenotypic variations. Here, using datasets from the gut microbiome of 996 patients from seven clinical trials, we systematically identify microbial genomic structural variants (SVs) using SGV-Finder. The associations between SVs and response, progression-free survival, overall survival, and immune-related adverse events are systematically explored by metagenome-wide association analysis and replicated in different cohorts. Associated SVs are located in multiple species, including Akkermansia muciniphila, Dorea formicigenerans, and Bacteroides caccae. We find genes that encode enzymes that participate in glucose metabolism be harbored in these associated regions. This work uncovers a nascent layer of gut microbiome heterogeneity that is correlated with hosts' prognosis following ICI treatment and represents an advance in our knowledge of the intricate relationships between microbiota and tumor immunotherapy.

Fungal secondary metabolism is governed by an RNA-binding protein CsdA/RsdA complex

Production of secondary metabolites is controlled by a complicated regulatory network in eukaryotic cells. Several layers of regulators are involved in this process, ranging from pathway-specific regulation, to epigenetic control, to global regulation. Here, we discover that interaction of an RNA-binding protein CsdA with a regulator RsdA coordinates fungal secondary metabolism. Employing a genetic deletion approach and transcriptome analysis as well as metabolomics analysis, we reveal that CsdA and RsdA synergistically regulate fungal secondary metabolism comprehensively. Mechanistically, comprehensive genetic and biochemical studies prove that RsdA and CsdA co-localize in the nucleus and physically interact to achieve their functions. In particular, we demonstrate that CsdA mediates rsdA expression by binding specific motif "GUCGGUAU" of its pre-mRNA at a post-transcriptional level. We thus uncover a mechanism in which RNA-binding protein physically interacts with, and controls the expression level of, the RsdA to coordinate fungal secondary metabolism.

Development of a portable on-site applicable metagenomic data generation workflow for enhanced pathogen and antimicrobial resistance surveillance

Rapid, accurate and comprehensive diagnostics are essential for outbreak prevention and pathogen surveillance. Real-time, on-site metagenomics on miniaturized devices, such as Oxford Nanopore Technologies MinION sequencing, could provide a promising approach. However, current sample preparation protocols often require substantial equipment and dedicated laboratories, limiting their use. In this study, we developed a rapid on-site applicable DNA extraction and library preparation approach for nanopore sequencing, using portable devices. The optimized method consists of a portable mechanical lysis approach followed by magnetic bead-based DNA purification and automated sequencing library preparation, and resulted in a throughput comparable to a current optimal, laboratory-based protocol using enzymatic digestion to lyse cells. By using spike-in reference communities, we compared the on-site method with other workflows, and demonstrated reliable taxonomic profiling, despite method-specific biases. We also demonstrated the added value of long-read sequencing by recovering reads containing full-length antimicrobial resistance genes, and attributing them to a host species based on the additional genomic information they contain. Our method may provide a rapid, widely-applicable approach for microbial detection and surveillance in a variety of on-site settings.

Genomic analysis and assessment of pathogenic (toxicogenic) potential of Staphylococcus haemolyticus and Bacillus paranthracis consortia isolated from bovine mastitis in Russia

Three stable microbial consortia, each composed of Bacillus paranthracis and Staphylococcus haemolyticus strains, were isolated from milk of cows diagnosed with mastitis in three geographically remote regions of Russia. The composition of these consortia remained stable following multiple passages on culture media. Apparently, this stability is due to the structure of the microbial biofilms formed by the communities. The virulence of the consortia depended on the B. paranthracis strains. It seems plausible that the ability of the consortia to cause mastitis in cattle was affected by mutations of the cytK gene of B. paranthracis.

Native microbiome dominates over host factors in shaping the probiotic genetic evolution in the gut

Probiotics often acquire potentially adaptive mutations in vivo, gaining new functional traits through gut selection. While both the host and microbiome can contribute to probiotics' genetic evolution, separating the microbiome and the host's contribution to such selective pressures remains challenging. Here, we introduced germ-free (GF) and specific pathogen-free (SPF) mouse models to track how probiotic strains, i.e., Lactiplantibacillus plantarum HNU082 (Lp082) and Bifidobacterium animalis subsp. lactis V9 (BV9), genetically evolved under selection pressures derived from host factors alone and both host and microbial ecological factors. Notably, compared to the genome of a probiotic strain before consumption, the host only elicited <15 probiotic mutations in probiotic genomes that emerged in the luminal environment of GF mice, while a total of 840 mutations in Lp082 mutants and 21,579 mutations in BV9 were found in SPF mice, <0.25% of those derived from both factors that were never captured by other experimental evolution studies, indicating that keen microbial competitions exhibited the predominant evolutionary force in shaping probiotic genetic composition (>99.75%). For a given probiotic, functional genes occurring in potentially adaptive mutations induced by hosts (GF mice) were all shared with those found in mutants of SPF mice. Collectively, the native microbiome consistently drove a more rapid and divergent genetic evolution of probiotic strains in seven days of colonization than host factors did. Our study further laid a theoretical foundation for genetically engineering probiotics for better gut adaptation through in vitro artificial gut ecosystems without the selection pressures derived from host factors.

Interspecific competition can drive plasmid loss from a focal species in a microbial community

Plasmids are key disseminators of antimicrobial resistance genes and virulence factors, and it is therefore critical to predict and reduce plasmid spread within microbial communities. The cost of plasmid carriage is a key metric that can be used to predict plasmids' ecological fate, and it is unclear whether plasmid costs are affected by growth partners in a microbial community. We carried out competition experiments and tracked plasmid maintenance using a model system consisting of a synthetic and stable five-species community and a broad host-range plasmid, engineered to carry different payloads. We report that both the cost of plasmid carriage and its long-term maintenance in a focal strain depended on the presence of competitors, and that these interactions were species specific. Addition of growth partners increased the cost of a high-payload plasmid to a focal strain, and accordingly, plasmid loss from the focal species occurred over a shorter time frame. We propose that the destabilising effect of interspecific competition on plasmid maintenance may be leveraged in clinical and natural environments to cure plasmids from focal strains.

Genomic insights into Bacillus subtilis MBB3B9 mediated aluminium stress mitigation for enhanced rice growth

Aluminium (Al) toxicity in acid soil ecosystems is a major impediment to crop production as it drastically affects plant root growth, thereby acquisition of nutrients from the soil. Plant growth-promoting bacteria offers an interesting avenue for promoting plant growth under an Al-phytotoxic environment. Here, we report the plant growth-promoting activities of an acid-tolerant isolate of Bacillus subtilis that could ameliorate acid-induced Al-stress in rice (Oryza sativa L.). The whole genome sequence data identified the major genes and genetic pathways in B. subtilis MBB3B9, which contribute to the plant growth promotion in acidic pH. Genetic pathways for organic acid production, denitrification, urea metabolism, indole-3-acetic acid (IAA) production, and cytokinin biosynthesis were identified as major genetic machinery for plant growth promotion and mitigation of Al-stress in plants. The in-vitro analyses revealed the production of siderophores and organic acid production as primary mechanisms for mitigation of Al-toxicity. Other plant growth-promoting properties such as phosphate solubilization, zinc solubilization, and IAA production were also detected in significant levels. Pot experiments involving rice under acidic pH and elevated concentrations of aluminium chloride (AlCl3) suggested that soil treatment with bacterial isolate MBB3B9 could enhance plant growth and productivity compared to untreated plants. A significant increase in plant growth and productivity was recorded in terms of plant height, chlorophyll content, tiller number, panicle number, grain yield, root growth, and root biomass production.

Multi-breed host rumen epithelium transcriptome and microbiome associations and their relationship with beef cattle feed efficiency

Understanding host-microbial interactions in the rumen and its influence on desirable production traits may lead to potential microbiota manipulation or genetic selection for improved cattle feed efficiency. This study investigated the host transcriptome and its correlation with the rumen archaea and bacteria differential abundance of two pure beef cattle breeds (Angus and Charolais) and one composite beef hybrid (Kinsella) divergent for residual feed intake (RFI; low-RFI vs. high-RFI). Using RNA-Sequencing of rumen tissue and 16S rRNA gene amplicon sequencing, differentially expressed genes (FDR ≤ 0.05, |log2(Fold-change) >|2) and differentially abundant (p-value < 0.05) archaea and bacteria amplicon sequence variants (ASV) were determined. Significant correlations between gene expression and ASVs (p-value < 0.05) were determine using Spearman correlation. Interesting associations with muscle contraction and the modulation of the immune system were observed for the genes correlated with bacterial ASVs. Potential functional candidate genes for feed efficiency status were identified for Angus (CCL17, CCR3, and CXCL10), Charolais (KCNK9, GGT1 and IL6), and Kinsella breed (ESR2). The results obtained here provide more insights regarding the applicability of target host and rumen microbial traits for the selection and breeding of more feed efficient beef cattle.

Probiotics alleviate constipation and inflammation in late gestating and lactating sows

Constipation and systemic inflammation are common in late pregnant and lactating sows, which cause health problems like uteritis, mastitis, dystocia, or even stillbirth, further influencing piglets' survival and growth. Probiotic supplementation can improve such issues, but the beneficial mechanism of relieving constipation and enhancing gut motility remains underexplored. This study aimed to investigate the effects and mechanism of probiotic supplementation in drinking water to late pregnant sows on constipation, inflammation, and piglets' growth performance. Seventy-four sows were randomly allocated to probiotic (n = 36) and control (n = 38) groups. Probiotic treatment significantly relieved sow constipation, enhanced serum IL-4 and IL-10 levels while reducing serum IL-1β, IL-12p40, and TNF-α levels, and increased piglet daily gain and weaning weight. Furthermore, probiotic administration reshaped the sow gut bacteriome and phageome structure/diversity, accompanied by increases in some potentially beneficial bacteria. At 113 days of gestation, the probiotic group was enriched in several gut microbial bioactive metabolites, multiple carbohydrate-active enzymes that degrade pectin and starch, fecal butyrate and acetate, and some serum metabolites involved in vitamin and amino acid metabolism. Our integrated correlation network analysis revealed that the alleviation of constipation and inflammation was associated with changes in the sow gut bacteriome, phageome, bioactive metabolic potential, and metabolism.

Negative chemotaxis of Ligilactobacillus agilis BKN88 against gut-derived substances

Ligilactobacillus agilis is a motile lactic acid bacterium found in the gastrointestinal tracts of animals. The findings of our previous study suggest that the motility of L. agilis BKN88 enables gut colonization in murine models. However, the chemotactic abilities of motile lactobacilli remain unknown. This study aimed to identify the gut-derived chemoeffectors and their corresponding chemoreceptors in L. agilis BKN88. Chemotaxis assays with chemotactic and non-chemotactic (ΔcheA) L. agilis strains revealed that low pH, organic acids, and bile salts served as repellents. L. agilis BKN88 was more sensitive to bile and acid than the gut-derived non-motile lactobacilli, implying that L. agilis might utilize motility and chemotaxis instead of exhibiting stress tolerance/resistance. L. agilis BKN88 contains five putative chemoreceptor genes (mcp1-mcp5). Chemotaxis assays using a series of chemoreceptor mutants revealed that each of the five chemoreceptors could sense multiple chemoeffectors and that these chemoreceptors were functionally redundant. Mcp2 and Mcp3 sensed all tested chemoeffectors. This study provides further insights into the interactions between chemoreceptors and ligands of motile lactobacilli and the unique ecological and evolutionary features of motile lactobacilli, which may be distinct from those of non-motile lactobacilli.

The minimal SUF system is not required for Fe-S cluster biogenesis in the methanogenic archaeon Methanosarcina acetivorans

Iron-sulfur (Fe-S) proteins are essential for the ability of methanogens to carry out methanogenesis and biological nitrogen fixation (diazotrophy). Nonetheless, the factors involved in Fe-S cluster biogenesis in methanogens remain largely unknown. The minimal SUF Fe-S cluster biogenesis system (i.e., SufBC) is postulated to serve as the primary system in methanogens. Here, the role of SufBC in Methanosarcina acetivorans, which contains two sufCB gene clusters, was investigated. The CRISPRi-dCas9 and CRISPR-Cas9 systems were utilized to repress or delete sufC1B1 and sufC2B2, respectively. Neither the dual repression of sufC1B1 and sufC2B2 nor the deletion of both sufC1B1 and sufC2B2 affected the growth of M. acetivorans under any conditions tested, including diazotrophy. Interestingly, deletion of only sufC1B1 led to a delayed-growth phenotype under all growth conditions, suggesting that the deletion of sufC2B2 acts as a suppressor mutation in the absence of sufC1B1. In addition, the deletion of sufC1B1 and/or sufC2B2 did not affect the total Fe-S cluster content in M. acetivorans cells. Overall, these results reveal that the minimal SUF system is not required for Fe-S cluster biogenesis in M. acetivorans and challenge the universal role of SufBC in Fe-S cluster biogenesis in methanogens.

Conjugative transfer of streptococcal prophages harboring antibiotic resistance and virulence genes

Prophages play important roles in the transduction of various functional traits, including virulence factors, but remain debatable in harboring and transmitting antimicrobial resistance genes (ARGs). Herein we characterize a prevalent family of prophages in Streptococcus, designated SMphages, which harbor twenty-five ARGs that collectively confer resistance to ten antimicrobial classes, including vanG-type vancomycin resistance locus and oxazolidinone resistance gene optrA. SMphages integrate into four chromosome attachment sites by utilizing three types of integration modules and undergo excision in response to phage induction. Moreover, we characterize four subtypes of Alp-related surface proteins within SMphages, the lethal effects of which are extensively validated in cell and animal models. SMphages transfer via high-frequency conjugation that is facilitated by integrative and conjugative elements from either donors or recipients. Our findings explain the widespread of SMphages and the rapid dissemination of ARGs observed in members of the Streptococcus genus.

Identifying and tracking mobile elements in evolving compost communities yields insights into the nanobiome

Microbial evolution is driven by rapid changes in gene content mediated by horizontal gene transfer (HGT). While mobile genetic elements (MGEs) are important drivers of gene flux, the nanobiome-the zoo of Darwinian replicators that depend on microbial hosts-remains poorly characterised. New approaches are necessary to increase our understanding beyond MGEs shaping individual populations, towards their impacts on complex microbial communities. A bioinformatic pipeline (xenoseq) was developed to cross-compare metagenomic samples from microbial consortia evolving in parallel, aimed at identifying MGE dissemination, which was applied to compost communities which underwent periodic mixing of MGEs. We show that xenoseq can distinguish movement of MGEs from demographic changes in community composition that otherwise confounds identification, and furthermore demonstrate the discovery of various unexpected entities. Of particular interest was a nanobacterium of the candidate phylum radiation (CPR) which is closely related to a species identified in groundwater ecosystems (Candidatus Saccharibacterium), and appears to have a parasitic lifestyle. We also highlight another prolific mobile element, a 313 kb plasmid hosted by a Cellvibrio lineage. The host was predicted to be capable of nitrogen fixation, and acquisition of the plasmid coincides with increased ammonia production. Taken together, our data show that new experimental strategies combined with bioinformatic analyses of metagenomic data stand to provide insight into the nanobiome as a driver of microbial community evolution.

Microdiversity of the vaginal microbiome is associated with preterm birth

Preterm birth (PTB) is the leading cause of neonatal morbidity and mortality. The vaginal microbiome has been associated with PTB, yet the mechanisms underlying this association are not fully understood. Understanding microbial genetic adaptations to selective pressures, especially those related to the host, may yield insights into these associations. Here, we analyze metagenomic data from 705 vaginal samples collected during pregnancy from 40 women who delivered preterm spontaneously and 135 term controls from the Multi-Omic Microbiome Study-Pregnancy Initiative. We find that the vaginal microbiome of pregnancies that ended preterm exhibited unique genetic profiles. It was more genetically diverse at the species level, a result which we validate in an additional cohort, and harbored a higher richness and diversity of antimicrobial resistance genes, likely promoted by transduction. Interestingly, we find that Gardnerella species drove this higher genetic diversity, particularly during the first half of the pregnancy. We further present evidence that Gardnerella spp. underwent more frequent recombination and stronger purifying selection in genes involved in lipid metabolism. Overall, our population genetics analyses reveal associations between the vaginal microbiome and PTB and suggest that evolutionary processes acting on vaginal microbes may play a role in adverse pregnancy outcomes such as PTB.

Social demographics determinants for resistome and microbiome variation of a multiethnic community in Southern Malaysia

The prevalence of antibiotic-resistant bacteria in Southeast Asia is a significant concern, yet there is limited research on the gut resistome and its correlation with lifestyle and environmental factors in the region. This study aimed to profile the gut resistome of 200 individuals in Malaysia using shotgun metagenomic sequencing and investigate its association with questionnaire data comprising demographic and lifestyle variables. A total of 1038 antibiotic resistance genes from 26 classes were detected with a mean carriage rate of 1.74 ± 1.18 gene copies per cell per person. Correlation analysis identified 14 environmental factors, including hygiene habits, health parameters, and intestinal colonization, that were significantly associated with the resistome (adjusted multivariate PERMANOVA, p < 0.05). Notably, individuals with positive yeast cultures exhibited a reduced copy number of 15 antibiotic resistance genes. Network analysis highlighted Escherichia coli as a major resistome network hub, with a positive correlation to 36 antibiotic-resistance genes. Our findings suggest that E. coli may play a pivotal role in shaping the resistome dynamics in Segamat, Malaysia, and its abundance is strongly associated with the community's health and lifestyle habits. Furthermore, the presence of yeast appears to be associated with the suppression of antibiotic-resistance genes.

Prevalence and epidemiological investigation of mgrB-dependent colistin resistance in extensively drug resistant Klebsiella pneumoniae in Iran

Carbapenemases-producing K. pneumoniae are challenging antimicrobial therapy of hospitalised patients, which is further complicated by colistin resistance. The aim of this study was to investigate the molecular epidemiological insights into carbapenemases-producing and colistin-resistant clinical K. pneumoniaeA total of 162 colistin resistant clinical strains of K. pneumoniae were collected during 2017-2019. Antimicrobial susceptibility and the colistin minimum inhibitory concentration were determined. Using PCR assay, the prevalence of resistance-associated genes including blaKPC, blaIMP, blaVIM, blaOXA-48, blaNDM-1 and mcr-1 to -9 was examined. Additionally, a PCR assay was used to examine the mgrB gene in colistin-resistant bacteria. 94.4% of the tested strains were resistant to imipenem and 96.3% were resistant to meropenem. Colistin resistance (MIC > 4 µg/L) was observed in 161 isolates (99.4%) by Colistin Broth Disk Elution method. The KPC enzyme was the most common carbapenemase and was identified in 95 strains (58.6%), followed by the IMP, VIM and OXA-48 detected in 47 (29%), 23 (14.2%) and 12 (7.4%) isolates, respectively. However, no NDM-1 gene was detected. Additionally, none of the studied isolates harbored mcr variants, while mgrB gene was observed in 152 (92.6%) isolates. Colistin resistance of K. pneumoniae isolates may be associated with mgrB gene mutation. To stop the spread of resistant K. pneumoniae, surveillance must be improved, infection prevention protocols must be followed, and antibiotic stewardship must be practised.

Small RNA sequencing of field Culex mosquitoes identifies patterns of viral infection and the mosquito immune response

Mosquito-borne disease remains a significant burden on global health. In the United States, the major threat posed by mosquitoes is transmission of arboviruses, including West Nile virus by mosquitoes of the Culex genus. Virus metagenomic analysis of mosquito small RNA using deep sequencing and advanced bioinformatic tools enables the rapid detection of viruses and other infecting organisms, both pathogenic and non-pathogenic to humans, without any precedent knowledge. In this study, we sequenced small RNA samples from over 60 pools of Culex mosquitoes from two major areas of Southern California from 2017 to 2019 to elucidate the virome and immune responses of Culex. Our results demonstrated that small RNAs not only allowed the detection of viruses but also revealed distinct patterns of viral infection based on location, Culex species, and time. We also identified miRNAs that are most likely involved in Culex immune responses to viruses and Wolbachia bacteria, and show the utility of using small RNA to detect antiviral immune pathways including piRNAs against some pathogens. Collectively, these findings show that deep sequencing of small RNA can be used for virus discovery and surveillance. One could also conceive that such work could be accomplished in various locations across the world and over time to better understand patterns of mosquito infection and immune response to many vector-borne diseases in field samples.

Genetic diversity of Helicobacter pylori type IV secretion system cagI and cagN genes and their association with clinical diseases

A number of cagPAI genes in the Helicobacter pylori genome are considered the most evolved genes under a diversifying selection and evolutionary pressure. Among them, cagI and cagN are described as a part of the two different-operon of cagPAI that are involved in the T4SS machinery, but the definite association of these factors with clinical manifestations is still unclear. A total of 70 H. pylori isolates were obtained from different gastroduodenal patients. All isolates were examined for the presence of primary H. pylori virulence genes by PCR analysis. Direct DNA sequence analysis was performed for the cagI and cagN genes. The results were compared with the reference strain. The cagI, cagN, cagA, cagL, vacA s1m1, vacA s1m2, vacA s2m2, babA2, sabA, and dupA genotypes were detected in 80, 91.4, 84, 91.4, 32.8, 42.8, 24.4, 97.1, 84.3, and 84.3% of the total isolates, respectively. The most variable codon usage in cagI was observed at residues 20-25, 55-60, 94, 181-199, 213-221, 241-268, and 319-320, while the most variable codon usage in CagN hypervariable motif (CagNHM) was observed at residues 53 to 63. Sequencing data analysis of cagN revealed a hypothetical hexapeptide motif (EAKDEN/K) in residues of 278-283 among six H. pylori isolates, which needs further studies to evaluate its putative function. The present study demonstrated a high prevalence of cagI and cagN genes among Iranian H. pylori isolates with gastroduodenal diseases. Furthermore, no significant correlation between cagI and cagN variants and clinical diseases was observed in the present study. However, all patients had a high prevalence of cagPAI genes including cagI, cagN, cagA, and cagL, which indicates more potential role of these genes in disease outcome.

Host-adaptive traits in the plant-colonizing Pseudomonas donghuensis P482 revealed by transcriptomic responses to exudates of tomato and maize

Pseudomonads are metabolically flexible and can thrive on different plant hosts. However, the metabolic adaptations required for host promiscuity are unknown. Here, we addressed this knowledge gap by employing RNAseq and comparing transcriptomic responses of Pseudomonas donghuensis P482 to root exudates of two plant hosts: tomato and maize. Our main goal was to identify the differences and the common points between these two responses. Pathways upregulated only by tomato exudates included nitric oxide detoxification, repair of iron-sulfur clusters, respiration through the cyanide-insensitive cytochrome bd, and catabolism of amino and/or fatty acids. The first two indicate the presence of NO donors in the exudates of the test plants. Maize specifically induced the activity of MexE RND-type efflux pump and copper tolerance. Genes associated with motility were induced by maize but repressed by tomato. The shared response to exudates seemed to be affected both by compounds originating from the plants and those from their growth environment: arsenic resistance and bacterioferritin synthesis were upregulated, while sulfur assimilation, sensing of ferric citrate and/or other iron carriers, heme acquisition, and transport of polar amino acids were downregulated. Our results provide directions to explore mechanisms of host adaptation in plant-associated microorganisms.

Metabolic cooperation between conspecific genotypic groups contributes to bacterial fitness

Microbial interactions are important for the survival of species and the stability of the microbial ecosystem. Although bacteria have diverse conspecific genotypes in the natural microbial ecosystem, little is known about whether wild-type strains within species would interact with each other and how the intraspecific interaction influences the growth of the species. In this work, using Lactobacillus acetotolerans, a dominant species with diverse conspecific genotypes in natural food fermentation ecosystems as a case, we studied the interactions between different genotypic groups of this species. In interspecific and intraspecific pairwise cocultures, the growth of L. acetotolerans decreased, but the increase of the phylogenetic similarity would reduce this negative effect, indicating a potential intraspecific interaction of this species. Meanwhile, the strain classification method affected the analysis of intraspecific interactions, which can be efficiently demonstrated using 99.5% average nucleotide identity (ANI) as the strain-level classification method. Using this ANI classification method, we revealed the population fitness significantly increased in cocultures of different genotypic groups. Facilitation involving 11 amino acids was identified between different ANI genotypic groups, which was beneficial for increasing population fitness. This work revealed that wild-type conspecific strains could interact with each other via cooperative metabolic changes and benefit each other to increase fitness. It shed new light on the survival and stability of species in natural microbial ecosystems.

Convergent Within-Host Adaptation of Pseudomonas aeruginosa through the Transcriptional Regulatory Network

Bacteria adapt to their host by mutating specific genes and by reprogramming their gene expression. Different strains of a bacterial species often mutate the same genes during infection, demonstrating convergent genetic adaptation. However, there is limited evidence for convergent adaptation at the transcriptional level. To this end, we utilize genomic data of 114 Pseudomonas aeruginosa strains, derived from patients with chronic pulmonary infection, and the P. aeruginosa transcriptional regulatory network. Relying on loss-of-function mutations in genes encoding transcriptional regulators and predicting their effects through the network, we demonstrate predicted expression changes of the same genes in different strains through different paths in the network, implying convergent transcriptional adaptation. Furthermore, through the transcription lens we associate yet-unknown processes, such as ethanol oxidation and glycine betaine catabolism, with P. aeruginosa host adaptation. We also find that known adaptive phenotypes, including antibiotic resistance, which were identified before as achieved by specific mutations, are achieved also through transcriptional changes. Our study has revealed novel interplay between the genetic and transcriptional levels in host adaptation, demonstrating the versatility of the adaptive arsenal of bacterial pathogens and their ability to adapt to the host conditions in a myriad of ways. IMPORTANCE Pseudomonas aeruginosa causes significant morbidity and mortality. The pathogen's remarkable ability to establish chronic infections greatly depends on its adaptation to the host environment. Here, we use the transcriptional regulatory network to predict expression changes during adaptation. We expand the processes and functions known to be involved in host adaptation. We show that the pathogen modulates the activity of genes during adaptation, including genes implicated in antibiotic resistance, both directly via genomic mutations and indirectly via mutations in transcriptional regulators. Furthermore, we detect a subgroup of genes whose predicted changes in expression are associated with mucoid strains, a major adaptive phenotype in chronic infections. We propose that these genes constitute the transcriptional arm of the mucoid adaptive strategy. Identification of different adaptive strategies utilized by pathogens during chronic infection has major promise in the treatment of persistent infections and opens the door to personalized tailored antibiotic treatment in the future.

Membrane-localized expression, production and assembly of Vibrio parahaemolyticus T3SS2 provides evidence for transertion

It has been proposed that bacterial membrane proteins may be synthesized and inserted into the membrane by a process known as transertion, which involves membrane association of their encoding genes, followed by coupled transcription, translation and membrane insertion. Here, we provide evidence supporting that the pathogen Vibrio parahaemolyticus uses transertion to assemble its type III secretion system (T3SS2), to inject virulence factors into host cells. We propose a two-step transertion process where the membrane-bound co-component receptor (VtrA/VtrC) is first activated by bile acids, leading to membrane association and expression of its target gene, vtrB, located in the T3SS2 pathogenicity island. VtrB, the transmembrane transcriptional activator of T3SS2, then induces the localized expression and membrane assembly of the T3SS2 structural components and its effectors. We hypothesize that the proposed transertion process may be used by other enteric bacteria for efficient assembly of membrane-bound molecular complexes in response to extracellular signals.

Effect of in vitro cultivation on human gut microbiota composition using 16S rDNA amplicon sequencing and metabolomics approach

Gut microbiota (GM) plays many key functions and helps maintain the host's health. Consequently, the development of GM cultivation under in vitro stimulating physiological conditions has gained extreme interest in different fields. In this study, we evaluated the impact of four culture media: Gut Microbiota Medium (GMM), Schaedler Broth (SM), Fermentation Medium (FM), and Carbohydrate Free Basal Medium (CFBM) on preserving the biodiversity and metabolic activity of human GM in batch in vitro cultures using PMA treatment coupled with 16S rDNA sequencing (PMA-seq) and LC-HR-MS/MS untargeted metabolomics supplemented with GC-MS SCFA profiling. Before the experiments, we determined the possibility of using the pooled faecal samples (MIX) from healthy donors (n = 15) as inoculum to reduce the number of variables and ensure the reproducibility of in vitro cultivation tests. Results showed the suitability of pooling faecal samples for in vitro cultivation study. Non-cultured MIX inoculum was characterized by higher α-diversity (Shannon effective count, and Effective microbial richness) compared to inocula from individual donors. After 24 h of cultivation, a significant effect of culture media composition on GM taxonomic and metabolomic profiles was observed. The SM and GMM had the highest α-diversity (Shannon effective count). The highest number of core ASVs (125) shared with non-cultured MIX inoculum and total SCFAs production was observed in the SM. These results might contribute to the development of standardized protocols for human GM in vitro cultivation by preventing methodological bias in the data.

Portable and cost-effective genetic detection and characterization of Plasmodium falciparum hrp2 using the MinION sequencer

The prevalence of Plasmodium falciparum hrp2 (pfhrp2)-deleted parasites threatens the efficacy of the most used and sensitive malaria rapid diagnostic tests and highlights the need for continued surveillance for this gene deletion. While PCR methods are adequate for determining pfhrp2 presence or absence, they offer a limited view of its genetic diversity. Here, we present a portable sequencing method using the MinION. Pfhrp2 amplicons were generated from individual samples, barcoded, and pooled for sequencing. To overcome potential crosstalk between barcodes, we implemented a coverage-based threshold for pfhrp2 deletion confirmation. Amino acid repeat types were then counted and visualized with custom Python scripts following de novo assembly. We evaluated this assay using well-characterized reference strains and 152 field isolates with and without pfhrp2 deletions, of which 38 were also sequenced on the PacBio platform to provide a standard for comparison. Of 152 field samples, 93 surpassed the positivity threshold, and of those samples, 62/93 had a dominant pfhrp2 repeat type. PacBio-sequenced samples with a dominant repeat-type profile from the MinION sequencing data matched the PacBio profile. This field-deployable assay can be used alone for surveilling pfhrp2 diversity or as a sequencing-based addition to the World Health Organization's existing deletion surveillance protocol.

AcrAB efflux pump impacts on the survival of adherent-invasive Escherichia coli strain LF82 inside macrophages

The tripartite complex AcrAB-TolC is the major RND pump in Escherichia coli and other Enterobacteriaceae. It consists of the AcrB transporter, which is embedded in the inner membrane, the AcrA adapter located in the periplasm, and the channel protein TolC responsible for the transport of substrates towards the extracellular environment. Besides conferring resistance to many classes of antibiotics, AcrAB plays a role in the pathogenesis and virulence of several bacterial pathogens. Here we report that the AcrAB pump heavily affects the infection process of the LF82 strain, the prototype of Adherent-Invasive Escherichia coli (AIEC) which are highly abundant in the ileal mucosa of Chron disease patients. We found that the deletion of genes encoding AcrA and/or AcrB leads to decreased survival of LF82 within macrophages. Ectopic AcrAB expression in a acrAB defective mutant restores the wild type condition. Furthermore, we demonstrate that inhibition of AcrB and replacement of the transporter with an unfunctional AcrB also interfere with bacterial viability inside macrophages. Overall, these data suggest a pivotal role of the AcrAB efflux pump in bacteria-host cell interactions also in AIEC.

Regulation of pseurotin A biosynthesis by GliZ and zinc in Aspergillus fumigatus

Recently, we reported that zinc regulates gliotoxin biosynthesis via ZafA, which is a zinc-responsive transcriptional activator. From an HPLC analysis of culture media of Aspergillus fumigatus, we found a trend of decreasing gliotoxin production but increasing pseurotin A and fumagillin production in proportion to the zinc concentration. The expression of the genes involved in pseurotin A biosynthesis was upregulated under high zinc concentrations. Furthermore, upregulated expression of pseurotin A biosynthetic genes and higher production of pseurotin A were observed in the zafA deletion strain. Interestingly, the deletion of gliZ, a transcriptional activator of gliotoxin biosynthesis genes, resulted in upregulated expression of pseurotin A biosynthetic genes and increased production of pseurotin A. We detected upregulation of fumR expression in the gliZ and zafA deletion mutants. The overexpression of gliZ observed in the zafA deletion mutant resulted in the failure of the mutant to increase pseurotin A production, which is a phenotype of the zafA deletion mutant. These results suggest that ZafA sequentially regulates pseurotin A biosynthesis through GliZ. Finally, we found through a murine virulence test that the gliZ and fumR double-deletion mutants showed a delayed death rate compared with the single-deletion mutants of either gliZ or fumR. Taken together, these results suggested that the biosynthesis of gliotoxin and pseurotin A are regulated in opposite ways by zinc utilization and that each secondary metabolite is synthesized when the synthesis of another secondary metabolite fails to protect it against the defense system of the host.

Illuminating the oral microbiome and its host interactions: tools and approaches for molecular ecological studies

A more comprehensive understanding of oral diseases like caries and periodontitis is dependent on an intimate understanding of the microbial ecological processes that are responsible for disease development. With this review, we provide a comprehensive overview of relevant molecular ecology techniques that have played critical roles in the current understanding of human oral biofilm development, interspecies interactions, and microbiome biogeography. The primary focus is on relevant technologies and examples available in the oral microbiology literature. However, most, if not all, of the described technologies should be readily adaptable for studies of microbiomes from other mucosal sites in the body. Therefore, this review is intended to serve as a reference guide used by microbiome researchers as they inevitably transition into molecular mechanistic studies of the many significant phenotypes observed clinically.

Large-scale genomic analysis of Mycobacterium tuberculosis reveals extent of target and compensatory mutations linked to multi-drug resistant tuberculosis

Resistance to isoniazid (INH) and rifampicin (RIF) first-line drugs in Mycobacterium tuberculosis (Mtb), together called multi-drug resistance, threatens tuberculosis control. Resistance mutations in katG (for INH) and rpoB (RIF) genes often come with fitness costs. To overcome these costs, Mtb compensatory mutations have arisen in rpoC/rpoA (RIF) and ahpC (INH) loci. By leveraging the presence of known compensatory mutations, we aimed to detect novel resistance mutations occurring in INH and RIF target genes. Across ~ 32 k Mtb isolates with whole genome sequencing (WGS) data, there were 6262 (35.7%) with INH and 5435 (30.7%) with RIF phenotypic resistance. Known mutations in katG and rpoB explained ~ 99% of resistance. However, 188 (0.6%) isolates had ahpC compensatory mutations with no known resistance mutations in katG, leading to the identification of 31 putative resistance mutations in katG, each observed in at least 3 isolates. These putative katG mutations can co-occur with other INH variants (e.g., katG-Ser315Thr, fabG1 mutations). For RIF, there were no isolates with rpoC/rpoA compensatory mutations and unknown resistance mutations. Overall, using WGS data we identified putative resistance markers for INH that could be used for genotypic drug-resistance profiling. Establishing the complete repertoire of Mtb resistance mutations will assist the clinical management of tuberculosis.

Breast cancer patients from the Midwest region of the United States have reduced levels of short-chain fatty acid-producing gut bacteria

As geographical location can impact the gut microbiome, it is important to study region-specific microbiome signatures of various diseases. Therefore, we profiled the gut microbiome of breast cancer (BC) patients of the Midwestern region of the United States. The bacterial component of the gut microbiome was profiled utilizing 16S ribosomal RNA sequencing. Additionally, a gene pathway analysis was performed to assess the functional capabilities of the bacterial microbiome. Alpha diversity was not significantly different between BC and healthy controls (HC), however beta diversity revealed distinct clustering between the two groups at the species and genera level. Wilcoxon Rank Sum test revealed modulation of several gut bacteria in BC specifically reduced abundance of those linked with beneficial effects such as Faecalibacterium prausnitzii. Machine learning analysis confirmed the significance of several of the modulated bacteria found by the univariate analysis. The functional analysis showed a decreased abundance of SCFA (propionate) production in BC compared to HC. In conclusion, we observed gut dysbiosis in BC with the depletion of SCFA-producing gut bacteria suggesting their role in the pathobiology of breast cancer. Mechanistic understanding of gut bacterial dysbiosis in breast cancer could lead to refined prevention and treatment.

Whole genome sequence characterization of Aspergillus terreus ATCC 20541 and genome comparison of the fungi A. terreus

Aspergillus terreus is well-known for lovastatin and itaconic acid production with biomedical and commercial importance. The mechanisms of metabolite formation have been extensively studied to improve their yield through genetic engineering. However, the combined repertoire of carbohydrate-active enzymes (CAZymes), cytochrome P450s (CYP) enzymes, and secondary metabolites (SMs) in the different A. terreus strains has not been well studied yet, especially with respect to the presence of biosynthetic gene clusters (BGCs). Here we present a 30 Mb whole genome sequence of A. terreus ATCC 20541 in which we predicted 10,410 protein-coding genes. We compared the CAZymes, CYPs enzyme, and SMs across eleven A. terreus strains, and the results indicate that all strains have rich pectin degradation enzyme and CYP52 families. The lovastatin BGC of lovI was linked with lovF in A. terreus ATCC 20541, and the phenomenon was not found in the other strains. A. terreus ATCC 20541 lacked a non-ribosomal peptide synthetase (AnaPS) participating in acetylaszonalenin production, which was a conserved protein in the ten other strains. Our results present a comprehensive analysis of CAZymes, CYPs enzyme, and SM diversities in A. terreus strains and will facilitate further research in the function of BGCs associated with valuable SMs.

Spatio-temporal patterns of Synechococcus oligotypes in Moroccan lagoonal environments

Synechococcus are unicellular cyanobacteria susceptible to environmental fluctuations and can be used as bioindicators of eutrophication in marine ecosystems. We examined their distribution in two Moroccan lagoons, Marchica on the Mediterranean coast and Oualidia on the Atlantic, in the summers of 2014 and 2015 using 16S rRNA amplicon oligotyping. Synechococcus representatives recruited a higher number of reads from the 16S rRNA in Marchica in comparison to Oualidia. We identified 31 Synechococcus oligotypes that clustered into 10 clades with different distribution patterns. The Synechococcus community was mainly represented by oligotype 1 (clade III) in Marchica. Cooccurring clades IV and I had an important relative abundance in Marchica in the summer of 2014, which is unusual, as these clades are widespread in cold waters. Moreover, Clades VII and subcluster "5.3" formed a sizeable percentage of the Synechococcus community in Marchica. Notably, we found low Synechococcus sequence counts in the Atlantic Lagoon. These results showed that the relative abundance of Synechococcus reads is not constant over space and time and that rare members of the Synechococcus community did not follow a consistent pattern. Further studies are required to decipher Synechococcus dynamics and the impact of environmental parameters on their spatial and temporal distributions.

Genomic epidemiology on the move

This Genome Watch article explores how portable sequencing technology is re-shaping the field of genomic epidemiology.

Undecaprenyl phosphate translocases confer conditional microbial fitness

The microbial cell wall is essential for maintenance of cell shape and resistance to external stressors¹. The primary structural component of the cell wall is peptidoglycan, a glycopolymer with peptide crosslinks located outside of the cell membrane¹. Peptidoglycan biosynthesis and structure are responsive to shifting environmental conditions such as pH and salinity^2-6, but the mechanisms underlying such adaptations are incompletely understood. Precursors of peptidoglycan and other cell surface glycopolymers are synthesized in the cytoplasm and then delivered across the cell membrane bound to the recyclable lipid carrier undecaprenyl phosphate⁷ (C55-P, also known as UndP). Here we identify the DUF368-containing and DedA transmembrane protein families as candidate C55-P translocases, filling a critical gap in knowledge of the proteins required for the biogenesis of microbial cell surface polymers. Gram-negative and Gram-positive bacteria lacking their cognate DUF368-containing protein exhibited alkaline-dependent cell wall and viability defects, along with increased cell surface C55-P levels. pH-dependent synthetic genetic interactions between DUF368-containing proteins and DedA family members suggest that C55-P transporter usage is dynamic and modulated by environmental inputs. C55-P transporter activity was required by the cholera pathogen for growth and cell shape maintenance in the intestine. We propose that conditional transporter reliance provides resilience in lipid carrier recycling, bolstering microbial fitness both inside and outside the host.

SARS-CoV-2 diagnostic testing rates determine the sensitivity of genomic surveillance programs

The first step in SARS-CoV-2 genomic surveillance is testing to identify people who are infected. However, global testing rates are falling as we emerge from the acute health emergency and remain low in many low- and middle-income countries (mean = 27 tests per 100,000 people per day). We simulated COVID-19 epidemics in a prototypical low- and middle-income country to investigate how testing rates, sampling strategies and sequencing proportions jointly impact surveillance outcomes, and showed that low testing rates and spatiotemporal biases delay time to detection of new variants by weeks to months and can lead to unreliable estimates of variant prevalence, even when the proportion of samples sequenced is increased. Accordingly, investments in wider access to diagnostics to support testing rates of approximately 100 tests per 100,000 people per day could enable more timely detection of new variants and reliable estimates of variant prevalence. The performance of global SARS-CoV-2 genomic surveillance programs is fundamentally limited by access to diagnostic testing.