Toward automatic reconstruction of a highly resolved tree of life

Mutagenesis and structural modeling implicate RME-8 IWN domains as conformational control points

After endocytosis, transmembrane cargo is differentially sorted into degradative or recycling pathways. This process is facilitated by recruitment into physically distinct degradative or recycling microdomains on the limiting membrane of individual endosomes. Endosomal sorting complexes required for transport (ESCRT) mark the degradative microdomain, while the recycling domain is marked by the retromer complex and associated proteins RME-8 and SNX-1. The separation of endosomal microdomains is also controlled by RME-8 and SNX-1, at least in part via removal of degradative component HRS/HGRS-1 from the recycling microdomain. This activity is likely due to recruitment and activation of chaperone Hsc70 on the endosome by the RME-8 DNAJ domain. To better understand the mechanism of RME-8 function we performed a new phylogenetic analysis of RME-8 and identified new conserved sequence features. In a complementary approach, we performed structure-function analysis that identified the C-terminus as important for microdomain localization and likely substrate binding, while N-terminal sequences beyond the known single N-terminal PH-like domain are important for endosome recruitment. Random mutagenesis identified IWN4, and by analogy IWN3, to be important for the autoinhibitory DNAJ domain binding, with IWN3 playing a critical role in HRS uncoating activity. Combining AlphaFold structural predictions with in vivo mutation analysis of RME-8, we propose a model whereby SNX-1 and the IWN domains control the conformation of RME-8 and hence the productive exposure of the DNAJ domain. Furthermore, we propose that the activation of RME-8 is cyclical, with SNX-1 acting as an activator and a target of RME-8 uncoating activity.

Impact of environmental factors and bacterial interactions on dust mite allergens in different indoor dust

Indoor dust is the main carrier of indoor pollutants, especially dust mite allergens and bacteria, they can trigger asthma, rhinitis, eczema and other allergic diseases. However, the interactions between dust mite allergens and bacterial communities in different types of indoor dust are not clear. The study focused on particulate and flocculent fibrous dust, explored the concentrations of Der p 1 (Dermatophagoides pteronyssinus) and Der f 1 (D. farinae) in 46 households in Changchun and their environmental influences, characterized the bacterial communities by high-throughput sequencing, and the interactions between Der p 1, Der f 1 and bacterial communities were explored. The results showed that Der p 1 and Der f 1 tended to accumulate more in flocculent fibrous dust, and Der p 1 predominated in the indoor dust samples. The floor height, years of housing occupancy and the living areas all affected the concentrations of dust mite allergens. In bacterial community, Proteobacteria, Firmicutes and Actinobacteria were leading phyla in the two types of dust. Kocuria, Blastococcus and Massilia were dominating genera in particulate dust and Acinetobacter, Lactobacillus, Corynebacterium_1 were dominating genera in flocculent fibrous dust. The overall diversity and species richness of bacteria in particulate dust were significantly higher than those in flocculent dust (p < 0.001). The living area was an important environmental factor affecting the bacterial community in flocculent fibrous dust (p < 0.01). The interaction between the relative abundance of Proteobacteria, Firmicutes and Actinobacteria and dust mite allergen concentrations significantly differed between the two dust types, indicating that bacteria could be used both as food and to establish symbiotic relationships with household dust mites (HDMs) hosts and provide nutrition.

Intraspecies characterization of bacteria via evolutionary modeling of protein domains

The ability to detect and characterize bacteria within a biological sample is crucial for the monitoring of infections and epidemics, as well as for the study of human health and its relationship with commensal microorganisms. To this aim, a commonly used technique is the 16S rRNA gene targeted sequencing. PCR-amplified 16S sequences derived from the sample of interest are usually clustered into the so-called Operational Taxonomic Units (OTUs) based on pairwise similarities. Then, representative OTU sequences are compared with reference (human-made) databases to derive their phylogeny and taxonomic classification. Here, we propose a new reference-free approach to define the phylogenetic distance between bacteria based on protein domains, which are the evolving units of proteins. We extract the protein domain profiles of 3368 bacterial genomes and we use an ecological approach to model their Relative Species Abundance distribution. Based on the model parameters, we then derive a new measurement of phylogenetic distance. Finally, we show that such model-based distance is capable of detecting differences between bacteria in cases in which the 16S rRNA-based method fails, providing a possibly complementary approach , which is particularly promising for the analysis of bacterial populations measured by shotgun sequencing.

Bacterial diversity of a floating vegetation (Phumdi) of Loktak Lake and its extracellular enzymes and bacterial antagonistic property

134 bacterial strains were isolated from phumdis of Loktak Lake. Through 16S rRNA sequencing, Bacillus sp. (23, 17.1%), Staphylococcus sp. (14, 10.4%), Pseudomonas sp. (11, 8.2%) and Acinetobacter sp. (8, 5.9%) were identified as the predominant bacterial taxa of Loktak Lake. B. pumulis (12, 8.9%), S. arlettae (4, 2.9%), P. knackmussii (6, 4.4%) are the leading species of Bacillus, Staphylococcus and Pseudomonas, respectively. Similarly, A. seifertii (2, 1.4%) and A. calcoaceticus (2, 1.4%) are the common species of Acinetobacter. 75 (55.9%) bacterial strains showed the ability to hydrolyze one or more extracellular enzymes tested. Among the extracellular enzymes produced by the bacterial isolates, the presence of elastase activity cannot be underestimated, since the enzyme is involved in the process of bacterial lung infection. Phosphate solubilizing activity could be seen in 11.1% of the bacterial isolates. 27 (20.1%) of the strains shown to have antagonistic activity against one or more tested pathogens. An isolate, MRC 52 showed antagonistic activity against eleven different pathogens including carbapenem resistant E. coli which was further subjected to extraction and identification of the biomolecule exerting the antimicrobial property. Based on GC-MS analysis, the bioactive compound was identified as phenyl ethyl alcohol.

Molecular Evolution and Inheritance Pattern of Sox Gene Family among Bovidae

Sox genes are an evolutionarily conserved family of transcription factors that play important roles in cellular differentiation and numerous complex developmental processes. In vertebrates, Sox proteins are required for cell fate decisions, morphogenesis, and the control of self-renewal in embryonic and adult stem cells. The Sox gene family has been well-studied in multiple species including humans but there has been scanty or no research into Bovidae. In this study, we conducted a detailed evolutionary analysis of this gene family in Bovidae, including their physicochemical properties, biological functions, and patterns of inheritance. We performed a genome-wide cataloguing procedure to explore the Sox gene family using multiple bioinformatics tools. Our analysis revealed a significant inheritance pattern including conserved motifs that are critical to the ability of Sox proteins to interact with the regulatory regions of target genes and orchestrate multiple developmental and physiological processes. Importantly, we report an important conserved motif, EFDQYL/ELDQYL, found in the SoxE and SoxF groups but not in other Sox groups. Further analysis revealed that this motif sequence accounts for the binding and transactivation potential of Sox proteins. The degree of protein–protein interaction showed significant interactions among Sox genes and related genes implicated in embryonic development and the regulation of cell differentiation. We conclude that the Sox gene family uniquely evolved in Bovidae, with a few exhibiting important motifs that drive several developmental and physiological processes.

Small secreted peptides encoded on the wheat (triticum aestivum L.) genome and their potential roles in stress responses

Small secreted peptides (SSPs) are important signals for cell-to-cell communication in plant, involved in a variety of growth and developmental processes, as well as responses to stresses. While a large number of SSPs have been identified and characterized in various plant species, little is known about SSPs in wheat, one of the most important cereal crops. In this study, 4,981 putative SSPs were identified on the wheat genome, among which 1,790 TaSSPs were grouped into 38 known SSP families. The result also suggested that a large number of the putaitive wheat SSPs, Cys-rich peptides in particular, remained to be characterized. Several TaSSP genes were found to encode multiple SSP domains, including CLE, HEVEIN and HAIRPININ domains, and two potentially novel TaSSP family DYY and CRP8CI were identified manually among unpredicted TaSSPs. Analysis on the transcriptomic data showed that a great proportion of TaSSPs were expressed in response to abiotic stresses. Exogenous application of the TaCEPID peptide encoded by TraesCS1D02G130700 enhanced the tolerance of wheat plants to drought and salinity, suggesting porential roles of SSPs in regulating stress responses in wheat.

A revisit to universal single-copy genes in bacterial genomes

Universal single-copy genes (USCGs) are widely used for species classification and taxonomic profiling. Despite many studies on USCGs, our understanding of USCGs in bacterial genomes might be out of date, especially how different the USCGs are in different studies, how well a set of USCGs can distinguish two bacterial species, whether USCGs can separate different strains of a bacterial species, to name a few. To fill the void, we studied USCGs in the most updated complete bacterial genomes. We showed that different USCG sets are quite different while coming from highly similar functional categories. We also found that although USCGs occur once in almost all bacterial genomes, each USCG does occur multiple times in certain genomes. We demonstrated that USCGs are reliable markers to distinguish different species while they cannot distinguish different strains of most bacterial species. Our study sheds new light on the usage and limitations of USCGs, which will facilitate their applications in evolutionary, phylogenomic, and metagenomic studies.

Genomic features of Mycoplasma bovis subtypes currently circulating in France

Background

Mycoplasma (M.) bovis is a major etiological agent of bovine respiratory disease, which is the most economically costly disease of cattle worldwide. Cattle disease surveillance on M. bovis is increasingly using gene-based techniques, such as multilocus sequence typing (MLST), or genome-based techniques such as core genome MLST that both require only partial genomic data. However, accurate up-to-date surveillance also demands complete, circular genomes that can be used as reference to track the evolution of the different lineages. Yet, in France, two of the main subtypes currently circulating still have no representing genome in public databases. Here, to address this gap, we provide and compare three new complete M. bovis genomes obtained from recent clinical isolates that represent major subtypes circulating in France and Europe.

Results

Genomes were obtained using a hybrid assembly strategy (Illumina and Nanopore) with fine-tuning of settings and inputs used in the Unicycler assembly pipeline, such as size selection of reads and quality trimming of the FASTQ files. The main characteristics and synteny of the genomes were compared. The three genomes mainly differed by their content in terms of mobile genetic elements, i.e. integrative conjugative elements (ICE) and insertion sequences (IS), a feature that impacts their structure. For instance, strain L15527, representing subtype3 (st3), harbours an exceptionally high number of ICEs, which results in a bigger-sized genome than all those previously described and could be associated with the propensity of st3 to gain and fix mutations through chromosomal transfer mechanisms. In contrast, strain F9160, of st1, is very close to the PG45 type strain isolated in 1961 in the USA, and harbours a huge number of IS. These features may be associated with an evolution towards a host-restricted state or in a “closed” host or environment reservoir until a recent re-emergence.

Conclusions

Whole-genome comparison of the three French M. bovis subtypes provides valuable resources for future studies combining epidemiology, phylogenetic data, and phylodynamic methods.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08818-9.

When LUCA met gnomAD: genetic constraints on universal genes in humans

LUCA, the last universal common ancestor, is the hypothetical most recent common ancestor of the three domains of life which share the universal genes (UG). It seems interesting to evaluate whether the UG phylogeny has had an impact on current Human gene constraints. A list of human homologs of UG was retrieved from the eggNOG database. We analyzed this LUCA gene (LG) group, and a random sample of 500 genes from the gnomAD database (RG group). Gene constraint metrics were retrieved from gnomAD and associations with Mendelian diseases and modes of inheritance were retrieved from OMIM. The LG group consisted of 277 genes and the RG group, 492 (8 genes were in LG group). 38.6% of the genes in the LG group and 25.2% of the genes in the RG group were associated with a Mendelian disease (p < 0.0001). The mode of inheritance was more often autosomal recessive (69.0 vs. 50.5%), and less often autosomal dominant (19.0 vs. 31.3%), or mixed (6.0 vs. 12.1%) for those associated with the LG group (p = 0.048). The LG group was significantly more constrained for missense variants (MOEUF, 0.919 vs. 0.997, p < 0.0001) and was borderline significantly more constrained for loss-of-function variants (LOEUF, 0.872 vs. 0.947, p = 0.051). These results suggest that the UG in humans differs from the rest of the genome in terms of constraints and associated Mendelian diseases. It suggests that phylogenic data can explain some of the characteristics of human genes and could help in interpreting variants.

Population Structure of German Cockroaches (Blattodea: Ectobiidae) in an Urban Environment Based on Single Nucleotide Polymorphisms

German cockroaches (Blattella germanica L.) harbor and disperse medically important pathogens and are a source of allergens that impact human health and wellbeing. Management of this pest requires an understanding of their distribution and dispersal. In this study, we collected German cockroaches from three apartment buildings in New Jersey, USA. We identified single-nucleotide polymorphisms (SNPs) from DNA extractions using next generation sequencing. We analyzed the SNPs and characterized cockroach population genetic structure using Fst, principal component, phylogenetic, and STRUCTURE analyses. We found significant differences in German cockroach population structure among the buildings. Within buildings, we found variable population structure that may be evidence for multiple colonization events. This study shows that SNPs derived from next generation sequencing provide a powerful tool for analyzing the genetic population structure of these medically important pests.

Expression of fatty acid synthase genes and their role in development and arboviral infection of Aedes aegypti

BACKGROUND

Fatty acids are the building blocks of complex lipids essential for living organisms. In mosquitoes, fatty acids are involved in cell membrane production, energy conservation and expenditure, innate immunity, development and reproduction. Fatty acids are synthesized by a multifunctional enzyme complex called fatty acid synthase (FAS). Several paralogues of FAS were found in the Aedes aegypti mosquito. However, the molecular characteristics and expression of some of these paralogues have not been investigated.

METHODS

Genome assemblies of Ae. aegypti were analyzed, and orthologues of human FAS was identified. Phylogenetic analysis and in silico molecular characterization were performed to identify the functional domains of the Ae. aegypti FAS (AaFAS). Quantitative analysis and loss-of-function experiments were performed to determine the significance of different AaFAS transcripts in various stages of development, expression following different diets and the impact of AaFAS on dengue virus, serotype 2 (DENV2) infection and transmission.

RESULTS

We identified seven putative FAS genes in the Ae. aegypti genome assembly, based on nucleotide similarity to the FAS proteins (tBLASTn) of humans, other mosquitoes and invertebrates. Bioinformatics and molecular analyses suggested that only five of the AaFAS genes produce mRNA and therefore represent complete gene models. Expression levels of AaFAS varied among developmental stages and between male and female Ae. aegypti. Quantitative analyses revealed that expression of AaFAS1, the putative orthologue of the human FAS, was highest in adult females. Transient knockdown (KD) of AaFAS1 did not induce a complete compensation by other AaFAS genes but limited DENV2 infection of Aag2 cells in culture and the midgut of the mosquito.

CONCLUSION

AaFAS1 is the predominant AaFAS in adult mosquitoes. It has the highest amino acid similarity to human FAS and contains all enzymatic domains typical of human FAS. AaFAS1 also facilitated DENV2 replication in both cell culture and in mosquito midguts. Our data suggest that AaFAS1 may play a role in transmission of dengue viruses and could represent a target for intervention strategies.

Analysis of Microbial Community Structure and Diversity in Burial Soil of Yangguanzhai Cemetery

As one of the common physical remains in archaeological discoveries, human bones are important bases for studying the history of human development, which is of great significance for exploring the evolution law of ancient human, reconstructing ancient human society, and tracking the development of human civilization. However, in the process of human bone burial, in addition to being affected by physical and chemical factors, it will also be affected by microorganisms in the buried soil, resulting in a variety of diseases. According to the determination and analysis of the microbial community structure and diversity in the burial soil of Yangguanzhai Site in Gaoling District in Xi’an city, Shaanxi Province, this paper attempts to explore the influence of microorganisms in the burial environment on human bones, in order to provide scientific proof for the microbial prevention and control of bone relics in the archaeological excavation site. For the first time, Illumina NovaSeq high-throughput sequencing technology was used to analyze the microbial community structure in the burial soil. At the phylum level, there were 8 dominant bacteria species in the soil samples of tombs, which were Firmicutes, Actinobacteriota, Actinobacteria, Proteobacteria, Acidobacteriota, Methylomirabilota, Chloroflexi, Bacteroidota. At the genus level, there were 12 dominant species in the soil samples of tombs, including MIZ17, MND1, Gaiella, oc32, Kroppenstedtia, Halomonas, Bacteroides, Dongia, Faecalibacterium, Nocardioides, Pseudomonas, Pseudonocardia. The overall microorganisms in the soil of Yangguanzhai Cemetery were relatively well-distributed, and the microbial community structure near human bones is the most abundant and diverse. Therefore, it is necessary to take some measures to control microorganisms and protect human bones.

Phylogenomics of the Phylum Proteobacteria: Resolving the Complex Relationships

Proteobacteria is one of the largest and phenotypically most diverse divisions within the domain bacteria. Due to the economic importance, this phylum demands an urgent need for a clear and scientifically sound classification system to streamline their characterization. The goal of our study was to carefully reevaluate the current system of classification and suggest changes wherein necessary. Phylogenetic trees of 84 Proteobacteria were constructed using single gene-based phylogeny involving 16S rRNA genes and protein sequences of 85 conserved genes, whole genome-based phylogenetic tree using CVtree3.0, amino acid Identity matrix tree, and concatenated tree with aforementioned conserved genes. The results of our study confirm the polyphyletic relationship between Desulfurella acetivorans, a Deltaproteobacteria with Epsilonproteobacteria. The group Syntrophobacterales was found to be polyphyletic with respect to Desulfarculus baarsii and the group Thiotrichales was found to be splitting in different phylogenetic trees. Placement of phylogenetic groups belonging to Rhodocyclales, Oceonospirilalles, and Chromatiales is controversial and requires further study and revisions. Based on our analysis, we strongly support reclassification of Magnetococcales as a separate class Etaproteobacteria. From our results, we conclude that concatenated trees of conserved proteins are a more accurate method for phylogenetic analysis, as compared to other methods used.

The human physiology of well-being: A systematic review on the association between neurotransmitters, hormones, inflammatory markers, the microbiome and well-being

To understand the pathways through which well-being contributes to health, we performed a systematic review according to the Preferred Reporting Items for Systematic Review and Meta-Analysis (PRISMA) guidelines on the association between well-being and physiological markers in four categories, neurotransmitters, hormones, inflammatory markers, and microbiome. We identified 91 studies. Neurotransmitter studies (k_{number of studies}=9) reported only a possible positive association between serotonin and well-being. For the hormone studies (k = 48), a lower momentary cortisol level was related to higher well-being (meta-analytic r = -0.06), and a steeper diurnal slope of cortisol levels. Inflammatory marker studies (k = 36) reported negative or non-significant relations with well-being, with meta-analytic estimates of respectively r = -0.07 and r = -0.05 for C-reactive protein and interleukin-6. Microbiome studies (k = 4) reported inconsistent associations between different bacteria abundance and well-being. The results indicate possible but small roles of serotonin, cortisol, and inflammatory markers in explaining differences in well-being. The inconsistent and limited results for other markers and microbiome require further research. Future directions for a complete picture of the physiological factors underlying well-being are proposed.

A consensus view of the proteome of the last universal common ancestor

The availability of genomic and proteomic data from across the tree of life has made it possible to infer features of the genome and proteome of the last universal common ancestor (LUCA). A number of studies have done so, all using a unique set of methods and bioinformatics databases. Here, we compare predictions across eight such studies and measure both their agreement with one another and with the consensus predictions among them. We find that some LUCA genome studies show a strong agreement with the consensus predictions of the others, but that no individual study shares a high or even moderate degree of similarity with any other individual study. From these observations, we conclude that the consensus among studies provides a more accurate depiction of the core proteome of the LUCA and its functional repertoire. The set of consensus LUCA protein family predictions between all of these studies portrays a LUCA genome that, at minimum, encoded functions related to protein synthesis, amino acid metabolism, nucleotide metabolism, and the use of common, nucleotide‐derived organic cofactors.

A Standardized Framework for Better Understanding of Phenotypic Differences within Bacterial Phyla Based on Protein Domain

We propose a standardized framework to classify target species based on their protein domains, which can be utilized in different contexts, like eukaryotes and prokaryotes. In this study, by applying the framework to the bacterial kingdom as an implementation example and comparing the results with the current taxonomy standards at the phylum level, we came to the conclusion that the sequence of domains rather than the content of domains in a protein and the presence of one domain rather than the number of occurrences of one domain play more important roles in deciding bacterial phenotypes as well as matching the current taxonomy. In addition, the comparison also helps us to better focus on the species that conflict with the current phylum category, as well as to further investigate their phenotypic or genotypic differences. IMPORTANCE A 3-step framework was designed which can be applied to clustering species based on their protein domains, and different candidate models are proposed in each step for better adaptation of various scenarios. We show its implementation for the bacterial kingdom as an example, which helps us to find the most appropriate model combination that will best reflect the relationship between domains and phenotypes in this context. In addition, identifying species that are distant in the results but should be closely related phylogenetically can help us to focus on the mismatch for better understanding of their key phenotypic or genotypic differences.

Protein Network Analysis of Whole Exome Sequencing of Severe Preeclampsia

Preeclampsia is a hypertensive disorder of pregnancy, which complicates up to 15% of US deliveries. It is an idiopathic disorder associated with several different phenotypes. We sought to determine if the genetic architecture of preeclampsia can be described by clusters of patients with variants in genes in shared protein interaction networks. We performed a case-control study using whole exome sequencing on early onset preeclamptic mothers with severe clinical features and control mothers with uncomplicated pregnancies between 2016 and 2020. A total of 143 patients were enrolled, 61 women with early onset preeclampsia with severe features based on ACOG criteria, and 82 control women at term, matched for race and ethnicity. A network analysis and visualization tool, Proteinarium, was used to confirm there are clusters of patients with shared gene networks associated with severe preeclampsia. The majority of the sequenced patients appear in two significant clusters. We identified one case dominant and one control dominant cluster. Thirteen genes were unique to the case dominated cluster. Among these genes, LAMB2, PTK2, RAC1, QSOX1, FN1, and VCAM1 have known associations with the pathogenic mechanisms of preeclampsia. Using bioinformatic analysis, we were able to identify subsets of patients with shared protein interaction networks, thus confirming our hypothesis about the genetic architecture of preeclampsia.

Diurnal and eating-associated microbial patterns revealed via high-frequency saliva sampling

The oral microbiome is linked to oral and systemic health, but its fluctuation under frequent daily activities remains elusive. Here, we sampled saliva at 10- to 60-min intervals to track the high-resolution microbiome dynamics during the course of human activities. This dense time series data showed that eating activity markedly perturbed the salivary microbiota, with tongue-specific Campylobacter concisus and Oribacterium sinus and dental plaque-specific Lautropia mirabilis, Rothia aeria, and Neisseria oralis increased after every meal in a temporal order. The observation was reproducible in multiple subjects and across an 11-mo period. The microbiome composition showed significant diurnal oscillation patterns at different taxonomy levels with Prevotella/Alloprevotella increased at night and Bergeyella HMT 206/Haemophilus slowly increased during the daytime. We also identified microbial co-occurring patterns in saliva that are associated with the intricate biogeography of the oral microbiome. Microbial source tracking analysis showed that the contributions of distinct oral niches to the salivary microbiome were dynamically affected by daily activities, reflecting the role of saliva in exchanging microbes with other oral sites. Collectively, our study provides insights into the temporal microbiome variation in saliva and highlights the need to consider daily activities and diurnal factors in design of oral microbiome studies.

Conformational spread drives the evolution of the calcium-calmodulin protein kinase II

The calcium calmodulin (Ca²⁺/CaM) dependent protein kinase II (CaMKII) decodes Ca²⁺ frequency oscillations. The CaMKIIα isoform is predominantly expressed in the brain and has a central role in learning. I matched residue and organismal evolution with collective motions deduced from the atomic structure of the human CaMKIIα holoenzyme to learn how its ring architecture abets function. Protein dynamic simulations showed its peripheral kinase domains (KDs) are conformationally coupled via lateral spread along the central hub. The underlying β-sheet motions in the hub or association domain (AD) were deconvolved into dynamic couplings based on mutual information. They mapped onto a coevolved residue network to partition the AD into two distinct sectors. A second, energetically stressed sector was added to ancient bacterial enzyme dimers for assembly of the ringed hub. The continued evolution of the holoenzyme after AD–KD fusion targeted the sector’s ring contacts coupled to the KD. Among isoforms, the α isoform emerged last and, it alone, mutated rapidly after the poikilotherm–homeotherm jump to match the evolution of memory. The correlation between dynamics and evolution of the CaMKII AD argues single residue substitutions fine-tune hub conformational spread. The fine-tuning could increase CaMKIIα Ca²⁺ frequency response range for complex learning functions.

Mutational Pattern Induced by 5-Fluorouracil and Oxaliplatin in the Gut Microbiome

Chemotherapeutic agents, such as 5-fluorouracil (5-FU) and oxaliplatin (Oxi), can not only kill the cancer cell but also influence the proliferation of gut microbiota; however, the interaction between these drugs and gut microbiota remains poorly understood. In this study, we developed a powerful framework for taxonomy composition and genomic variation analysis to investigate the mutagenesis effect and proliferation influence of chemotherapeutic agents, such as 5-FU and Oxi, on gut microbiota and the interaction between these drugs and gut microbiota during chemotherapy. Using the gut microbiome data, we detected 1.45 million variations among the chemotherapy groups and found the drugs significantly affected mutation signatures of gut microbiota. Oxi notably increased transversion rate, whereas 5-FU reduced the rate. Traits related to cell division and nutrient mobilization showed evidence of strong selection pressure from chemotherapeutic agents. In addition, drug-associated bacteriome shift patterns and functional alterations were found: the metabolism changes in the 5-FU group implied that gut microbiota could provide additional nicotinamide adenine dinucleotide (NAD⁺) to inhibit cancer cell autophagy; in the Oxi group, the ribosome and lysine biosynthesis genes were obviously enriched. Our study provides a blueprint for characterizing the role of microbes and drug–microbe interaction in the gut microbiota response to chemotherapy.

Ancient Origins of Cytoskeletal Crosstalk: Spectraplakin-like Proteins Precede the Emergence of Cortical Microtubule Stabilization Complexes as Crosslinkers

Adhesion between cells and the extracellular matrix (ECM) is one of the prerequisites for multicellularity, motility, and tissue specialization. Focal adhesions (FAs) are defined as protein complexes that mediate signals from the ECM to major components of the cytoskeleton (microtubules, actin, and intermediate filaments), and their mutual communication determines a variety of cellular processes. In this study, human cytoskeletal crosstalk proteins were identified by comparing datasets with experimentally determined cytoskeletal proteins. The spectraplakin dystonin was the only protein found in all datasets. Other proteins (FAK, RAC1, septin 9, MISP, and ezrin) were detected at the intersections of FAs, microtubules, and actin cytoskeleton. Homology searches for human crosstalk proteins as queries were performed against a predefined dataset of proteomes. This analysis highlighted the importance of FA communication with the actin and microtubule cytoskeleton, as these crosstalk proteins exhibit the highest degree of evolutionary conservation. Finally, phylogenetic analyses elucidated the early evolutionary history of spectraplakins and cortical microtubule stabilization complexes (CMSCs) as model representatives of the human cytoskeletal crosstalk. While spectraplakins probably arose at the onset of opisthokont evolution, the crosstalk between FAs and microtubules is associated with the emergence of metazoans. The multiprotein complexes contributing to cytoskeletal crosstalk in animals gradually gained in complexity from the onset of metazoan evolution.

ThermoBase: A database of the phylogeny and physiology of thermophilic and hyperthermophilic organisms

Thermophiles and hyperthermophiles are those organisms which grow at high temperature (> 40°C). The unusual properties of these organisms have received interest in multiple fields of biological research, and have found applications in biotechnology, especially in industrial processes. However, there are few listings of thermophilic and hyperthermophilic organisms and their relevant environmental and physiological data. Such repositories can be used to standardize definitions of thermophile and hyperthermophile limits and tolerances and would mitigate the need for extracting organism data from diverse literature sources across multiple, sometimes loosely related, research fields. Therefore, we have developed ThermoBase, a web-based and freely available database which currently houses comprehensive descriptions for 1238 thermophilic or hyperthermophilic organisms. ThermoBase reports taxonomic, metabolic, environmental, experimental, and physiological information in addition to literature resources. This includes parameters such as coupling ions for chemiosmosis, optimal pH and range, optimal temperature and range, optimal pressure, and optimal salinity. The database interface allows for search features and sorting of parameters. As such, it is the goal of ThermoBase to facilitate and expedite hypothesis generation, literature research, and understanding relating to thermophiles and hyperthermophiles within the scientific community in an accessible and centralized repository. ThermoBase is freely available online at the Astrobiology Habitable Environments Database (AHED; https://ahed.nasa.gov), at the Database Center for Life Science (TogoDB; http://togodb.org/db/thermobase), and in the S1 File.

Reclassification of 11 Members of the Family Rhodobacteraceae at Genus and Species Levels and Proposal of Pseudogemmobacter hezensis sp. nov

A novel Gram-stain-negative, aerobic, motile bacterial strain, D13-10-4-6^T, was isolated from the bark sample of Populus × euramericana. The strain could grow at 15-35°C, at pH 6-10 and in 0-4% (w/v) NaCl, and the strain tested positive for oxidase and catalase activities. The main polar lipids were phosphatidylmonomethylethanolamine, diphosphatidylglycerol, phosphatidylglycerol, and phosphatidylethanolamine. The main respiratory quinone was Q-10, and the predominant fatty acid was C_18:1 ω7c. The phylogenetic analyses showed that the strain belonged to the genus Pseudogemmobacter of the family Rhodobacteraceae. The family Rhodobacteraceae is an ecologically diverse group that includes bacteria from aquatic to terrestrial ecosystems. As a consequence, the classification of the family Rhodobacteraceae is difficult, not least when the early taxonomy work relied heavily on 16S rRNA gene analysis. Recently, the taxonomic status of many members of the family has been revised based on the genome analysis; however, there are still some classification conflicts due to the lack of genome sequences and parallel publication time. In this study, phylogenetic trees based on 16S rRNA gene, gyrB gene, and 120 concatenated proteins, the average amino acid identity (AAI) and percentage of conserved proteins (POCP) have been used for the analysis of strain D13-10-4-6^T and other members of 15 genera within the family to further clarify their taxonomic relationships. For the data of phylogeny, AAI, and POCP, the taxonomic proposals are (1) reclassification of Rhodobacter tardus as the type species of a novel genus, Stagnihabitans gen. nov., as Stagnihabitans tardus comb. nov.; (2) reclassification of Tabrizicola alkalilacus, Tabrizicola sediminis, Tabrizicola algicola into a novel genus, Pseudotabrizicola gen. nov., as Pseudotabrizicola alkalilacus comb. nov., Pseudotabrizicola sediminis comb. nov., Pseudotabrizicola algicola comb. nov.; (3) reclassification of Rhodobacter sediminicola into the genus Cereibacter as Cereibacter sediminicola comb. nov.; (4) reclassification of Rhodobacter flagellatus, Rhodobacter thermarum, and Xinfangfangia soli into the genus Tabrizicola as Tabrizicola flagellatus comb. nov., Tabrizicola thermarum comb. Nov., and Tabrizicola soli comb. nov.; (5) reclassification of Xinfangfangia humi into the genus Pseudogemmobacter as Pseudogemmobacter humicola comb. nov.; (6) classification of strain D13-10-4-6^T as a novel species of the genus Pseudogemmobacter, for which the name P. hezensis sp. nov. is proposed, the type strain is D13-10-4-6^T (= CFCC 12033^T = KCTC 82215^T).

Integrating remote sensing with ecology and evolution to advance biodiversity conservation

Remote sensing has transformed the monitoring of life on Earth by revealing spatial and temporal dimensions of biological diversity through structural, compositional and functional measurements of ecosystems. Yet, many aspects of Earth's biodiversity are not directly quantified by reflected or emitted photons. Inclusive integration of remote sensing with field-based ecology and evolution is needed to fully understand and preserve Earth's biodiversity. In this Perspective, we argue that multiple data types are necessary for almost all draft targets set by the Convention on Biological Diversity. We examine five key topics in biodiversity science that can be advanced by integrating remote sensing with in situ data collection from field sampling, experiments and laboratory studies to benefit conservation. Lowering the barriers for bringing these approaches together will require global-scale collaboration.

Omicron-included mutation-induced changes in epitopes of SARS-CoV-2 spike protein and effectiveness assessments of current antibodies

The COVID-19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is spreading globally and continues to rage, posing a serious threat to human health and life quality. Antibody therapy and vaccines both have shown great efficacy in the prevention and treatment of COVID-19, whose development progress and adaptation range have attracted wide attention. However, with the emergence of variant strains of SARS-CoV-2, the neutralization activity of therapeutic or vaccine-induced antibodies may be reduced, requiring long-term virus monitoring and drug upgrade in response to its evolution. In this paper, conformational changes including continuous epitopes (CPs), discontinuous epitopes (DPs) and recognition interfaces of the three representative SARS-CoV-2 spike protein (SP) mutants (i.e., the Delta (B.1.617.2), Mu (B.1.621) and Omicron (B.1.1.529) strains), were analyzed to evaluate the effectiveness of current mainstream antibodies. The results showed that the conformation of SP wild type (WT) and mutants both remained stable, while the local antigenic epitopes underwent significant changes. Sufficient flexibility of SP CPs is critical for effective antibody recognition. The DPs of Delta, Mu and Omicron variants have showed stronger binding to human angiotensin converting enzyme-2 (hACE2) than WT; the possible drug resistance mechanisms of antibodies against three different epitopes (i.e., NTD_DP, RBD1_DP and RBD2_DP) were also proposed, respectively; the RBD2 of Delta, NTD of Mu, NTD and RBD2 of Omicron are deserve more attention in the subsequent design of next-generation vaccines. The simulation results not only revealed structural characteristics of SP antigenic epitopes, but also provided guidance for antibody modification, vaccine design and effectiveness evaluation.

Landscape of mobile genetic elements and their antibiotic resistance cargo in prokaryotic genomes

Prokaryotic Mobile Genetic Elements (MGEs) such as transposons, integrons, phages and plasmids, play important roles in prokaryotic evolution and in the dispersal of cargo functions like antibiotic resistance. However, each of these MGE types is usually annotated and analysed individually, hampering a global understanding of phylogenetic and environmental patterns of MGE dispersal. We thus developed a computational framework that captures diverse MGE types, their cargos and MGE-mediated horizontal transfer events, using recombinases as ubiquitous MGE marker genes and pangenome information for MGE boundary estimation. Applied to ∼84k genomes with habitat annotation, we mapped 2.8 million MGE-specific recombinases to six operational MGE types, which together contain on average 13% of all the genes in a genome. Transposable elements (TEs) dominated across all taxa (∼1.7 million occurrences), outnumbering phages and phage-like elements (<0.4 million). We recorded numerous MGE-mediated horizontal transfer events across diverse phyla and habitats involving all MGE types, disentangled and quantified the extent of hitchhiking of TEs (17%) and integrons (63%) with other MGE categories, and established TEs as dominant carriers of antibiotic resistance genes. We integrated all these findings into a resource (proMGE.embl.de), which should facilitate future studies on the large mobile part of genomes and its horizontal dispersal.

The Giardia ventrolateral flange is a lamellar membrane protrusion that supports attachment

Attachment to the intestinal epithelium is critical to the lifestyle of the ubiquitous parasite Giardia lamblia. The ventrolateral flange is a sheet-like membrane protrusion at the interface between parasites and attached surfaces. This structure has been implicated in attachment, but its role has been poorly defined. Here, we identified a novel actin associated protein with putative WH2-like actin binding domains we named Flangin. Flangin complexes with Giardia actin (GlActin) and is enriched in the ventrolateral flange making it a valuable marker for studying the flanges' role in Giardia biology. Live imaging revealed that the flange grows to around 1 μm in width after cytokinesis, then remains uniform in size during interphase, grows in mitosis, and is resorbed during cytokinesis. A flangin truncation mutant stabilizes the flange and blocks cytokinesis, indicating that flange disassembly is necessary for rapid myosin-independent cytokinesis in Giardia. Rho family GTPases are important regulators of membrane protrusions and GlRac, the sole Rho family GTPase in Giardia, was localized to the flange. Knockdown of Flangin, GlActin, and GlRac result in flange formation defects. This indicates a conserved role for GlRac and GlActin in forming membrane protrusions, despite the absence of canonical actin binding proteins that link Rho GTPase signaling to lamellipodia formation. Flangin-depleted parasites had reduced surface contact and when challenged with fluid shear force in flow chambers they had a reduced ability to remain attached, confirming a role for the flange in attachment. This secondary attachment mechanism complements the microtubule based adhesive ventral disc, a feature that may be particularly important during mitosis when the parental ventral disc disassembles in preparation for cytokinesis. This work supports the emerging view that Giardia's unconventional actin cytoskeleton has an important role in supporting parasite attachment.

Deep distributed computing to reconstruct extremely large lineage trees

Phylogeny estimation (the reconstruction of evolutionary trees) has recently been applied to CRISPR-based cell lineage tracing, allowing the developmental history of an individual tissue or organism to be inferred from a large number of mutated sequences in somatic cells. However, current computational methods are not able to construct phylogenetic trees from extremely large numbers of input sequences. Here, we present a deep distributed computing framework to comprehensively trace accurate large lineages (FRACTAL) that substantially enhances the scalability of current lineage estimation software tools. FRACTAL first reconstructs only an upstream lineage of the input sequences and recursively iterates the same produce for its downstream lineages using independent computing nodes. We demonstrate the utility of FRACTAL by reconstructing lineages from >235 million simulated sequences and from >16 million cells from a simulated experiment with a CRISPR system that accumulates mutations during cell proliferation. We also successfully applied FRACTAL to evolutionary tree reconstructions and to an experiment using error-prone PCR (EP-PCR) for large-scale sequence diversification.

Strain-level characterization of broad host range mobile genetic elements transferring antibiotic resistance from the human microbiome

Mobile genetic elements (MGEs) carrying antibiotic resistance genes (ARGs) disseminate ARGs when they mobilise into new bacterial hosts. The nature of such horizontal gene transfer (HGT) events between human gut commensals and pathogens remain poorly characterised. Here, we compare 1354 cultured commensal strains (540 species) to 45,403 pathogen strains (12 species) and find 64,188 MGE-mediated ARG transfer events between the two groups using established methods. Among the 5931 MGEs, we find 15 broad host range elements predicted to have crossed different bacterial phyla while also occurring in animal and environmental microbiomes. We experimentally demonstrate that predicted broad host range MGEs can mobilise from commensals Dorea longicatena and Hungatella hathewayi to pathogen Klebsiella oxytoca, crossing phyla simultaneously. Our work establishes the MGE-mediated ARG dissemination network between human gut commensals and pathogens and highlights broad host range MGEs as targets for future ARG dissemination management.

Evolutionary genomic relationships and coupling in MK-STYX and STYX pseudophosphatases

The dual specificity phosphatase (DUSP) family has catalytically inactive members, called pseudophosphatases. They have mutations in their catalytic motifs that render them enzymatically inactive. This study analyzes the significance of two pseudophosphatases, MK-STYX [MAPK (mitogen-activated protein kinase phosphoserine/threonine/tyrosine-binding protein]) and STYX (serine/threonine/tyrosine-interacting protein), throughout their evolution and provides measurements and comparison of their evolutionary conservation. Phylogenetic trees were constructed to show any deviation from various species evolutionary paths. Data was collected on a large set of proteins that have either one of the two domains of MK-STYX, the DUSP domain or the cdc-25 homology (CH2) /rhodanese-like domain. The distance between species pairs for MK-STYX or STYX and Ka/Ks ratio were calculated. In addition, both pseudophosphatases were ranked among a large set of related proteins, including the active homologs of MK-STYX, MKP (MAPK phosphatase)-1 and MKP-3. MK-STYX had one of the highest species-species protein distances and was under weaker purifying selection pressure than most proteins with its domains. In contrast, the protein distances of STYX were lower than 82% of the DUSP-containing proteins and was under one of the strongest purifying selection pressures. However, there was similar selection pressure on the N-terminal sequences of MK-STYX, STYX, MKP-1, and MKP-3. We next perform statistical coupling analysis, a process that reveals interconnected regions within the proteins. We find that while MKP-1,-3, and STYX all have 2 functional units (sectors), MK-STYX only has one, and that MK-STYX is similar to MKP-3 in the evolutionary coupling of the active site and KIM domain. Within those two domains, the mean coupling is also most similar for MK-STYX and MKP-3. This study reveals striking distinctions between the evolutionary patterns of MK-STYX and STYX, suggesting a very specific role for each pseudophosphatase, further highlighting the relevance of these atypical members of DUSP as signaling regulators. Therefore, our study provides computational evidence and evolutionary reasons to further explore the properties of pseudophosphatases, in particular MK-STYX and STYX.

Isolation and Characterization of Lytic Bacteriophages Targeting Diverse Enterobacter spp. Clinical Isolates

Background: Enterobacter spp. are opportunistic pathogens that cause nosocomial infections. Bacteriophages could be used to treat antibiotic-resistant Enterobacter infections. Materials and Methods: We used 10 genetically diverse clinical Enterobacter spp. isolates to identify lytic bacteriophages in hospital and municipal wastewater. Comparative genomics was performed on host bacterial isolates and isolated phages. Activity of each phage against all 10 host isolates was determined. We also tested phage activity against paired isolates from two patients who developed ceftazidime-avibactam resistance. Results: Bacteria belonged to three Enterobacter species and Klebsiella aerogenes. We isolated 12 bacteriophages, most of which belonged to the Myoviridae and Autographiviridae families. Most phages were able to lyse multiple bacterial isolates, and many lysed isolates of different species. Ceftazidime-avibactam-resistant isolates were still phage susceptible, and one isolate showed increased susceptibility compared with the parent isolate. Conclusion: The phages we isolated expand the diversity of Enterobacter-targeting phages, and could be useful for treating antibiotic-resistant Enterobacter infections.

Expression and distribution of synaptotagmin family members in the zebrafish retina

Synaptotagmins belong to a large family of proteins. Although various synaptotagmins have been implicated as Ca²⁺ sensors for vesicle replenishment and release at conventional synapses, their roles at retinal ribbon synapses remain incompletely understood. Zebrafish is a widely used experimental model for retinal research. We therefore investigated the homology between human, rat, mouse, and zebrafish synaptotagmins 1-10 using a bioinformatics approach. We also characterized the expression and distribution of various synaptotagmin (syt) genes in the zebrafish retina using RT-PCR, qPCR, and in situhybridization, focusing on the family members whose products likely underlie Ca²⁺ -dependent exocytosis in the central nervous system (synaptotagmins 1, 2, 5, and 7). Most zebrafish synaptotagmins are well conserved and can be grouped in the same classes as mammalian synaptotagmins, based on crucial amino acid residues needed for coordinating Ca²⁺ binding and determining phospholipid binding affinity. The only exception is synaptotagmin 1b, which lacks 34 amino acid residues in the C2B domain and is therefore unlikely to bind Ca²⁺ there. Additionally, the products of zebrafish syt5a and syt5b genes share identity with mammalian class 1 and 5 synaptotagmins. Zebrafish syt1, syt2, syt5, and syt7 paralogues are found in the zebrafish brain, eye, and retina, excepting syt1b, which is only present in the brain. The complementary expression pattern of the remaining paralogues in the retina suggests that syt1a and syt5a may underlie synchronous release and syt7a and syt7b may mediate asynchronous release or other Ca²⁺ -dependent processes in different retinal neurons.

An estimate of the deepest branches of the tree of life from ancient vertically-evolving genes

Core gene phylogenies provide a window into early evolution, but different gene sets and analytical methods have yielded substantially different views of the tree of life. Trees inferred from a small set of universal core genes have typically supported a long branch separating the archaeal and bacterial domains. By contrast, recent analyses of a broader set of non-ribosomal genes have suggested that Archaea may be less divergent from Bacteria, and that estimates of inter-domain distance are inflated due to accelerated evolution of ribosomal proteins along the inter-domain branch. Resolving this debate is key to determining the diversity of the archaeal and bacterial domains, the shape of the tree of life, and our understanding of the early course of cellular evolution. Here, we investigate the evolutionary history of the marker genes key to the debate. We show that estimates of a reduced Archaea-Bacteria (AB) branch length result from inter-domain gene transfers and hidden paralogy in the expanded marker gene set. By contrast, analysis of a broad range of manually curated marker gene datasets from an evenly sampled set of 700 Archaea and Bacteria reveals that current methods likely underestimate the AB branch length due to substitutional saturation and poor model fit; that the best-performing phylogenetic markers tend to support longer inter-domain branch lengths; and that the AB branch lengths of ribosomal and non-ribosomal marker genes are statistically indistinguishable. Furthermore, our phylogeny inferred from the 27 highest-ranked marker genes recovers a clade of DPANN at the base of the Archaea and places the bacterial Candidate Phyla Radiation (CPR) within Bacteria as the sister group to the Chloroflexota.

Curtobacterium spp. and Curtobacterium flaccumfaciens: Phylogeny, Genomics-Based Taxonomy, Pathogenicity, and Diagnostics

The genus of Curtobacterium, belonging to the Microbacteriaceae family of the Actinomycetales order, includes economically significant pathogenic bacteria of soybeans and other agricultural crops. Thorough phylogenetic and full-genome analysis using the latest genomic data has demonstrated a complex and contradictory taxonomic picture within the group of organisms classified as the Curtobacterium species. Based on these data, it is possible to delineate about 50 new species and to reclassify a substantial part of the Curtobacterium strains. It is suggested that 53 strains, including most of the Curtobacterium flaccumfaciens pathovars, can compose a monophyletic group classified as C. flaccumfaciens. A genomic analysis using the most recent inventory of bacterial chromosomal and plasmid genomes deposited to GenBank confirmed the possible role of Microbacteriaceae plasmids in pathogenicity and demonstrated the existence of a group of related plasmids carrying virulence factors and possessing a gene distantly related to DNA polymerase found in bacteriophages and archaeal and eukaryotic viruses. A PCR diagnostic assay specific to the genus Curtobacterium was developed and tested. The presented results assist in the understanding of the evolutionary relations within the genus and can lay the foundation for further taxonomic updates.

Computational analyses of bacterial strains from shotgun reads

Shotgun sequencing is routinely employed to study bacteria in microbial communities. With the vast amount of shotgun sequencing reads generated in a metagenomic project, it is crucial to determine the microbial composition at the strain level. This study investigated 20 computational tools that attempt to infer bacterial strain genomes from shotgun reads. For the first time, we discussed the methodology behind these tools. We also systematically evaluated six novel-strain-targeting tools on the same datasets and found that BHap, mixtureS and StrainFinder performed better than other tools. Because the performance of the best tools is still suboptimal, we discussed future directions that may address the limitations.

A robust approach to estimate relative phytoplankton cell abundances from metagenomes

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.

Transmission of and Infection With COVID-19 Among Vaccinated and Unvaccinated Attendees of an Indoor Wedding Reception in Minnesota

IMPORTANCE

Characterizing rates of SARS-CoV-2 infection among vaccinated and unvaccinated persons with the same exposure is critical to understanding the association of vaccination with the risk of infection with the Delta variant. Additionally, evidence of Delta variant transmission by children to vaccinated adults has important public health implications.

OBJECTIVE

To characterize transmission and infection of SARS-CoV-2 among vaccinated and unvaccinated attendees of an indoor wedding reception.

DESIGN, SETTING, AND PARTICIPANTS

This cohort study included attendees at an indoor wedding reception in Minnesota in July 2021. Data were collected from REDCap surveys and routine surveillance interviews. The full list of attendees and a partial list of emails were obtained. Fifty-seven attendees completed the emailed survey. Eighteen additional attendees were identified from the state health department COVID-19 surveillance database.

EXPOSURES

Attendance at an indoor event.

MAIN OUTCOMES AND MEASURES

Risk of SARS-CoV-2 infection among vaccinated and unvaccinated attendees, identification of an index case, whole genome sequencing (WGS) to identify the COVID-19 variant, understanding of transmission patterns, and assessment of secondary transmission. The primary case definition was an individual with a positive SARS-CoV-2 test who attended the wedding in the 14 days prior to their illness.

RESULTS

Data were gathered for 75 attendees (mean [SE] age, 37.5 [13.7] years; 57 [76%] female individuals), of whom 56 (75%) were fully vaccinated, 4 (5%) were partially vaccinated, and 15 (20%) were unvaccinated. Of 62 attendees who were tested, 29 (47%) tested positive, including 16 of 46 fully vaccinated attendees (35%), 2 of 4 partially vaccinated attendees (50%), and 11 of 12 unvaccinated attendees (92%). Being unvaccinated was associated with a higher risk of infection compared with being vaccinated (risk ratio, 2.64; 95% CI, 1.71-4.06; P = .001). One unvaccinated adult required hospitalization. An unvaccinated child who was symptomatic on the event date was identified as the index case. Eleven specimens were available for WGS. All sequenced specimens were closely related and were identified as the Delta variant. WGS supported secondary transmission from a vaccinated individual with SARS-CoV-2.

CONCLUSIONS AND RELEVANCE

This cohort study identified a COVID-19 Delta variant outbreak at an indoor event despite a high proportion of vaccinated attendees. It found that vaccination was associated with a reduced risk of infection.

Population genomics of Bacillus anthracis from an anthrax hyperendemic area reveals transmission processes across spatial scales and unexpected within-host diversity

Genomic sequencing has revolutionized our understanding of bacterial disease epidemiology, but remains underutilized for zoonotic pathogens in remote endemic settings. Anthrax, caused by the spore-forming bacterium Bacillus anthracis, remains a threat to human and animal health and rural livelihoods in low- and middle-income countries. While the global genomic diversity of B. anthracis has been well-characterized, there is limited information on how its populations are genetically structured at the scale at which transmission occurs, critical for understanding the pathogen's evolution and transmission dynamics. Using a uniquely rich dataset, we quantified genome-wide SNPs among 73 B. anthracis isolates derived from 33 livestock carcasses sampled over 1 year throughout the Ngorongoro Conservation Area, Tanzania, a region hyperendemic for anthrax. Genome-wide SNPs distinguished 22 unique B. anthracis genotypes (i.e. SNP profiles) within the study area. However, phylogeographical structure was lacking, as identical SNP profiles were found throughout the study area, likely the result of the long and variable periods of spore dormancy and long-distance livestock movements. Significantly, divergent genotypes were obtained from spatio-temporally linked cases and even individual carcasses. The high number of SNPs distinguishing isolates from the same host is unlikely to have arisen during infection, as supported by our simulation models. This points to an unexpectedly wide transmission bottleneck for B. anthracis, with an inoculum comprising multiple variants being the norm. Our work highlights that inferring transmission patterns of B. anthracis from genomic data will require analytical approaches that account for extended and variable environmental persistence, as well as co-infection.

Evidence for the Rapid and Divergent Evolution of Mycoplasmas: Structural and Phylogenetic Analysis of Enolases

Mycoplasmas are a group of prokaryotes without cell walls that have evolved through several rounds of degenerative evolution. With a low cell DNA G + C content and definitively long genetic lineages, mycoplasmas are thought to be in a state of rapid evolution. However, little associated evidence has been provided. Enolase is a key enzyme in glycolysis that is widely found in all species from the three domains, and it is evolutionarily conserved. In our previous studies, enolase acted as a virulence factor and participated in cell-surface adhesion in Mycoplasma hyopneumoniae. Furthermore, unique loop regions were first found in the crystal structure of Mhp Eno. Here, enolase structures from Mycoplasma pneumoniae and Mycoplasma bovis were determined. An extra helix 7 is specific and conservatively found in almost all mycoplasma enolases, as confirmed by crystal structures and sequence alignment. Particular motifs for helix 7, which is composed of F-K/G-K-L/F-K-X-A-I, have been proposed and could be regarded as molecular markers. To our surprise, the genetic distances between any two mycoplasma enolases were obviously longer than those between the two corresponding species themselves, indicating divergent evolution of mycoplasma enolases, whereas no horizontal gene transfer was detected in mycoplasma enolase genens. Furthermore, different evolutionary patterns were adopted by different loop regions of mycoplasma enolase. Enolases from different Mycoplasma species also showed different affinities for PLG and fibronectin. Our results indicate the rapid and divergent evolution of mycoplasma enolase and mycoplasmas. This study will also aid understanding the independent evolution of Mycoplasma species after separation from their common ancestor.

Phylogenetic systematics of Yphthimoides Forster, 1964 and related taxa, with notes on the biogeographical history of Yphthimoides species

Species losses are increasing and may have an impact on our understanding of patterns of evolutionary pathways and phylogenetic relationships among the groups being lost. The knowledge of such patterns can contribute to preventing future losses by identifying which lineages have higher or lower diversification rates, thus informing conservation strategies. Recent years have seen a significant growth in studies of butterfly systematics, allowing a better understanding of evolutionary relationships among most groups and revealing significant taxonomic chaos in several groups. One of the latter groups is the nymphalid subtribe Euptychiina (Satyrinae), which has been shown to include a number of non-monophyletic genera based on recent molecular phylogenetic analyses. Among others, these genera include Yphthimoides, which is widespread throughout the Neotropical region but particularly diverse in the southeastern Neotropics, and a pair of related genera, Pharneuptychia Forster, 1964 and Moneuptychia Forster, 1964. Using molecular data, this study scope and aims was to provide a phylogenetic hypothesis that corroborates Yphthimoides as presently conceived being non-monophyletic, a result reinforced by a comparative study of the male genitalic morphology. Our results also show that Pharneuptychia and Moneuptychia, plus a species misplaced elsewhere in the Euptychiina, Euptychoides castrensis (Schaus, 1902), form a well supported clade, and that the latter 'species' is a complex of cryptic species. We therefore propose a number of taxonomic rearrangements in the present work to resolve these issues: Yphthimoides eriphule (A. Butler, 1867) will be moved to a new genus; Y. affinis (A. Butler, 1867), Y. maepius (Godart, [1824]), Y. mimula (Hayward, 1954), Y. neomaenas (Hayward, 1967) and Y. mythra (Weymer, 1911) are being transferred to Malaveria Viloria & Benmesbah, 2021; Pharneuptychia innocentia (Godart, [1824]) will be moved to another genus to be described; and Euptychoides castrensis, Pharneuptychia romanina (Bryk, 1953) and Yphthimoides viviana (Romieux, 1927) are being moved to Moneuptychia. The dating of divergences points to a split between the ancestral lineage of Yphthimoides and its sister group, Carminda Ebert and Dias, inDias 1998, during the last half of the Miocene, around 11.86 Mya, and to the diversification of the Pharneuptychia during the same time 11.35 (± 3.52) Mya. Biogeographic analysis showed that the most recent common ancestor of Yphthimoides started to diversify either in the the Brazilian Cerrado savannas or in a combined area of Cerrado and South Atlantic Forest, with a possible change in the ancestral habitat of Carminda. Furthermore, ancestral character mapping favors a savanna origin hypothesis over a forest origin hypothesis.

Longitudinal, virological, and serological assessment of hospitalized COVID-19 patients

Here we described the virological and serological assessment of 23 COVID-19 patients hospitalized and followed up in Milan, Italy, during the first wave of COVID-19 pandemic. Nasopharyngeal (NPS), anal swabs, and blood samples were collected from 23 COVID-19 patients, at hospital admission, and periodically up to discharge, for a median time of 20 days (3-83 days). RNA was isolated and tested for SARS-CoV-2 by qRT-PCR; anti-SARS-CoV-2 IgM and IgG antibody titers were evaluated in serum samples by ELISA. SARS-CoV-2 genome was detected in the NPS swabs of the 23 patients, at the admission, and 8/19 (42.1%) were still positive at the discharge. Anal swabs were positive to SARS-CoV-2 RNA detection in 20/23 (86.9%) patients; 6/19 (31.6%) were still positive at discharge. The mean time of RNA negative conversion was 17 days (4-36 days) and 33 days (4-77 days), for NPS and anal swabs, respectively. SARS-CoV-2-RNA was detected in the blood of 6/23 (26.1%) patients. Thirteen/23 (56.5%) and 17/23 (73.9%) patients were seropositive for IgM and IgG, respectively, at the admission, and the median IgM and IgG levels significantly (p < 0.05) increased after 13 days. Although the limited cohort size, our report provides evidence that SARS-CoV-2 is shed through multiple routes, with important implications in healthcare settings.

Deciphering Microbial Gene Family Evolution Using Duplication-Transfer-Loss Reconciliation and RANGER-DTL

Phylogenetic reconciliation has emerged as a principled, highly effective technique for investigating the origin, spread, and evolutionary history of microbial gene families. Proper application of phylogenetic reconciliation requires a clear understanding of potential pitfalls and sources of error, and knowledge of the most effective reconciliation-based tools and protocols to use to maximize accuracy. In this book chapter, we provide a brief overview of Duplication-Transfer-Loss (DTL) reconciliation, the standard reconciliation model used to study microbial gene families and provide a step-by-step computational protocol to maximize the accuracy of DTL reconciliation and minimize false-positive evolutionary inferences.

Phylogenetic Signal, Congruence, and Uncertainty across Bacteria and Archaea

Reconstruction of the Tree of Life is a central goal in biology. Although numerous novel phyla of bacteria and archaea have recently been discovered, inconsistent phylogenetic relationships are routinely reported, and many inter-phylum and inter-domain evolutionary relationships remain unclear. Here, we benchmark different marker genes often used in constructing multidomain phylogenetic trees of bacteria and archaea and present a set of marker genes that perform best for multidomain trees constructed from concatenated alignments. We use recently-developed Tree Certainty metrics to assess the confidence of our results and to obviate the complications of traditional bootstrap-based metrics. Given the vastly disparate number of genomes available for different phyla of bacteria and archaea, we also assessed the impact of taxon sampling on multidomain tree construction. Our results demonstrate that biases between the representation of different taxonomic groups can dramatically impact the topology of resulting trees. Inspection of our highest-quality tree supports the division of most bacteria into Terrabacteria and Gracilicutes, with Thermatogota and Synergistota branching earlier from these superphyla. This tree also supports the inclusion of the Patescibacteria within the Terrabacteria as a sister group to the Chloroflexota instead of as a basal-branching lineage. For the Archaea, our tree supports three monophyletic lineages (DPANN, Euryarchaeota, and TACK/Asgard), although we note the basal placement of the DPANN may still represent an artifact caused by biased sequence composition. Our findings provide a robust and standardized framework for multidomain phylogenetic reconstruction that can be used to evaluate inter-phylum relationships and assess uncertainty in conflicting topologies of the Tree of Life.

Solid-State NMR Investigations of Extracellular Matrixes and Cell Walls of Algae, Bacteria, Fungi, and Plants

Extracellular matrixes (ECMs), such as the cell walls and biofilms, are important for supporting cell integrity and function and regulating intercellular communication. These biomaterials are also of significant interest to the production of biofuels and the development of antimicrobial treatment. Solid-state nuclear magnetic resonance (ssNMR) and magic-angle spinning-dynamic nuclear polarization (MAS-DNP) are uniquely powerful for understanding the conformational structure, dynamical characteristics, and supramolecular assemblies of carbohydrates and other biomolecules in ECMs. This review highlights the recent high-resolution investigations of intact ECMs and native cells in many organisms spanning across plants, bacteria, fungi, and algae. We spotlight the structural principles identified in ECMs, discuss the current technical limitation and underexplored biochemical topics, and point out the promising opportunities enabled by the recent advances of the rapidly evolving ssNMR technology.

Real-time kinetic studies of Mycobacterium tuberculosis LexA-DNA interaction

Transcriptional repressor, LexA, regulates the ‘SOS’ response, an indispensable bacterial DNA damage repair machinery. Compared with its Escherichia coli ortholog, LexA from Mycobacterium tuberculosis (Mtb) possesses a unique N-terminal extension of additional 24 amino acids in its DNA-binding domain (DBD) and 18 amino acids insertion at its hinge region that connects the DBD to the C-terminal dimerization/autoproteolysis domain. Despite the importance of LexA in ‘SOS’ regulation, Mtb LexA remains poorly characterized and the functional importance of its additional amino acids remained elusive. In addition, the lack of data on kinetic parameters of Mtb LexA–DNA interaction prompted us to perform kinetic analyses of Mtb LexA and its deletion variants using Bio-layer Interferometry (BLI). Mtb LexA is seen to bind to different ‘SOS’ boxes, DNA sequences present in the operator regions of damage-inducible genes, with comparable nanomolar affinity. Deletion of 18 amino acids from the linker region is found to affect DNA binding unlike the deletion of the N-terminal stretch of extra 24 amino acids. The conserved RKG motif has been found to be critical for DNA binding. Overall, the present study provides insights into the kinetics of the interaction between Mtb LexA and its target ‘SOS’ boxes. The kinetic parameters obtained for DNA binding of Mtb LexA would be instrumental to clearly understand the mechanism of ‘SOS’ regulation and activation in Mtb.

The Wzi outer membrane protein mediates assembly of a tight capsular polysaccharide layer on the Acinetobacter baumannii cell surface

Identification of novel therapeutic targets is required for developing alternate strategies to treat infections caused by the extensively drug-resistant bacterial pathogen, Acinetobacter baumannii. As capsular polysaccharide (CPS) is a prime virulence determinant required for evasion of host immune defenses, understanding the pathways for synthesis and assembly of this discrete cell-surface barrier is important. In this study, we assess cell-bound and cell-free CPS material from A. baumannii AB5075 wildtype and transposon library mutants and demonstrate that the Wzi outer membrane protein is required for the proper assembly of the CPS layer on the cell surface. Loss of Wzi resulted in an estimated 4.4-fold reduction in cell-associated CPS with a reciprocal increase in CPS material shed in the extracellular surrounds. Transmission electron microscopy revealed a disrupted CPS layer with sparse patches of CPS on the external face of the outer membrane when Wzi function was lost. However, this genotype did not have a significant effect on biofilm formation. Genetic analysis demonstrated that the wzi gene is ubiquitous in the species, though the nucleotide sequences were surprisingly diverse. Though divergence was not concomitant with variation at the CPS biosynthesis K locus, an association between wzi type and the first sugar of the CPS representing the base of the structure most likely to interact with Wzi was observed.

N. S, Pachamuthu K, Chenna S, Nair A, Shivaprasad PV. Essential role of γ-clade RNA-dependent RNA polymerases in rice development and yield-related traits is linked to their atypical polymerase activities regulating specific genomic regions

RNA-dependent RNA polymerases (RDR) generate double-stranded (ds)RNA triggers for RNA silencing across eukaryotes. Among the three clades, α-clade and β-clade members are key components of RNA silencing and mediators of stress responses across eukaryotes. However, γ-clade members are unusual in that they are represented in phylogenetically distant plants and fungi, and their functions are unknown. Using genetic, bioinformatic and biochemical methods, we show that γ-clade RDRs from Oryza sativa L. are involved in plant development as well as regulation of expression of coding and noncoding RNAs. Overexpression of γ-clade RDRs in transgenic rice and tobacco plants resulted in robust growth phenotype, whereas their silencing in rice displayed strong inhibition of growth. Small (s)RNA and RNA-seq analysis of OsRDR3 mis-expression lines suggested that it is specifically involved in the regulation of repeat-rich regions in the genome. Biochemical analysis confirmed that OsRDR3 has robust polymerase activities on both single stranded (ss)RNA and ssDNA templates similar to the activities reported for α-clade RDRs such as AtRDR6. Our results provide the first evidence of the importance of γ-clade RDRs in plant development, their atypical biochemical activities and their contribution to the regulation of gene expression.

A portable library of phosphate-depletion based synthetic promoters for customable and automata control of gene expression in bacteria

Industrial biotechnology gene expression systems relay on constitutive promoters compromising cellular growth from the start of the bioprocess, or on inducible devices, which require manual addition of cognate inducers. To overcome this shortcoming, we engineered an automata regulatory system based on cell-stress mechanisms. Specifically, we engineered a synthetic and highly portable phosphate-depletion library of promoters inspired by bacterial PHO starvation system (Pliar promoters). Furthermore, we fully characterized 10 synthetic promoters within the background of two well-known bacterial workhorses such as E. coli W and P. putida KT2440. The promoters displayed an interesting host-dependent performance and a wide strength spectrum ranging from 0.4- to 1.3-fold when compared to the wild-type phosphatase alkaline promoter (PphoA). By comparing with available gene expression systems, we proved the suitability of this new library for the automata and effective decoupling of growth from production in P. putida. Growth phase-dependent expression of these promoters could therefore be activated by fine tuning the initial concentration of phosphate in the medium. Finally, the Pliar library was implemented in the SEVA platform in a ready-to-use mode allowing its broad use by the scientific community.

Prevalence and Specificity of Chemoreceptor Profiles in Plant-Associated Bacteria

Chemosensory pathways are among the most abundant prokaryotic signal transduction systems, allowing bacteria to sense and respond to environmental stimuli. Signaling is typically initiated by the binding of specific molecules to the ligand binding domain (LBD) of chemoreceptor proteins (CRs). Although CRs play a central role in plant-microbiome interactions such as colonization and infection, little is known about their phylogenetic and ecological specificity. Here, we analyzed 82,277 CR sequences from 11,806 representative microbial species covering the whole prokaryotic phylogeny, and we classified them according to their LBD type using a de novo homology clustering method. Through phylogenomic analysis, we identified hundreds of LBDs that are found predominantly in plant-associated bacteria, including several LBDs specific to phytopathogens and plant symbionts. Functional annotation of our catalogue showed that many of the LBD clusters identified might constitute unknown types of LBDs. Moreover, we found that the taxonomic distribution of most LBD types that are specific to plant-associated bacteria is only partially explained by phylogeny, suggesting that lifestyle and niche adaptation are important factors in their selection. Finally, our results show that the profile of LBD types in a given genome is related to the lifestyle specialization, with plant symbionts and phytopathogens showing the highest number of niche-specific LBDs. The LBD catalogue and information on how to profile novel genomes are available at https://github.com/compgenomicslab/CRs.

IMPORTANCE Considering the enormous variety of LBDs at sensor proteins, an important question resides in establishing the forces that have driven their evolution and selection. We present here the first clear demonstration that environmental factors play an important role in the selection and evolution of LBDs. We were able to demonstrate the existence of LBD families that are highly enriched in plant-associated bacteria but show a wide phylogenetic spread. These findings offer a number of research opportunities in the field of single transduction, such as the exploration of similar relationships in chemoreceptors of bacteria with a different lifestyle, like those inhabiting or infecting the human intestine. Similarly, our results raise the question whether similar LBD types might be shared by members of different sensor protein families. Lastly, we provide a comprehensive catalogue of CRs classified by their LBD region that includes a large number of putative new LBD types.