Protein contact map prediction using multiple sequence alignment dropout and consistency learning for sequences with fewer homologs The prediction of protein contact map needs enough normalized number of effective sequence (Nf) in multiple sequence alignment (MSA). When Nf is small, the predicted contact maps are often not satisfactory. To solve this problem, we randomly selected a small part of sequence homologs for proteins with large Nf to generate MSAs with small Nf. From these MSAs, input features were generated and were passed through a consistency learning network, aiming to get the same results when using the features generated from the MSA with large Nf. The results showed that this method effectively improves the prediction accuracy of protein contact maps with small Nf.
A memetic algorithm enables efficient local and global all-atom protein-protein docking with backbone and sidechain flexibility Protein complex formation is encoded by specific interactions at the atomic scale, but the computational cost of modeling proteins at this level often requires the use of simplified energy models and limited conformational flexibility. In particular, the use of all-atom energy functions, backbone and sidechain flexibility results in rugged energy landscapes that are difficult to explore. In this study we develop a protein-protein docking algorithm, EvoDOCK, that combine the strength of a differential evolution algorithm for efficient exploration of the global search space with the benefits of a local optimization method to refine detailed atomic interactions. EvoDOCK enabled accurate and fast local and global protein-protein docking using an all-atom energy function with side-chain flexibility. Comparison with a standard method built on Monte Carlo optimization demonstrated improved accuracy and with increases in computational speed of up to 35 times. The evolutionary algorithm also enabled efficient atomistic docking with backbone flexibility.
Apolipoprotein M attenuates anthracycline cardiotoxicity and lysosomal injury ObjectivesDetermine the role of apolipoprotein M (ApoM) in anthracycline (Dox) cardiotoxicity.

BackgroundApoM binds the cardioprotective sphingolipid sphingosine-1-phosphate (S1P). Circulating ApoM is inversely associated with mortality in human heart failure (HF).

MethodsIn the Penn HF Study (PHFS), we tested the relationship between ApoM and mortality in a subset with anthracycline-induced cardiomyopathy. We measured ApoM in humans and mice treated with Dox and utilized hepatic ApoM transgenic (ApomTG), ApoM knockout (ApomKO), ApoM knock-in mice with impaired S1P binding, and S1P receptor 3 (S1PR3) knockout mice in Dox cardiotoxicity. We assayed autophagy in left ventricular tissue from anthracycline-induced HF patients versus donor controls.

ResultsApoM was inversely associated with mortality in PHFS, and Dox reduced circulating ApoM in mice and breast cancer patients. ApomTG mice were protected from Dox-induced cardiac dysfunction and loss of left ventricular mass. ApomTG attenuated Dox-induced impairment in autophagic flux in vivo and accumulation of insoluble p62, which was also observed in the myocardium of patients with anthracycline-induced HF. In vehicle-treated mice, ApoM negatively regulated transcription factor EB (TFEB), a master regulator of autophagy and lysosomal biogenesis. The effect of ApoM on TFEB required both S1P binding and S1PR3. In the presence of Dox, ApoM preserved TFEB and cardiomyocyte lysosomal abundance assessed as lysosomal associated membrane protein 1 positive structures in vivo, while S1P mimetic pretreatment of cardiomyocytes prevented Dox-induced changes in lysosomal pH.

ConclusionsApoM attenuates Dox cardiotoxicity via the autophagy-lysosome pathway. The association between ApoM and reduced mortality may be explained by its role in sustaining autophagy.

HighlightsO_LICirculating ApoM is inversely associated with survival in human anthracycline-induced cardiomyopathy
C_LIO_LIAnthracycline treatment reduces circulating ApoM in humans and mice
C_LIO_LIIncreasing ApoM attenuates doxorubicin cardiotoxicity, lysosomal injury and preserves myocardial autophagic flux, but does not impact doxorubicin anti-neoplastic efficacy
C_LIO_LIAutophagic impairment is characteristic of human anthracycline cardiomyopathy
C_LI Thermodynamic constraints on the assembly and diversity of microbial ecosystems are different near to and far from equilibrium Non-equilibrium thermodynamics has long been an area of substantial interest to ecologists because most fundamental biological processes, such as protein synthesis and respiration, are inherently energy-consuming. However, most of this interest has focused on developing coarse ecosystem-level maximisation principles, providing little insight into underlying mechanisms that lead to such emergent constraints. Microbial communities are a natural system to decipher this mechanistic basis because their interactions in the form of substrate consumption, metabolite production, and cross-feeding can be described explicitly in thermodynamic terms. Previous work has considered how thermodynamic constraints impact competition between pairs of species, but restrained from analysing how this manifests in complex dynamical systems. To address this gap, we develop a thermodynamic microbial community model with fully reversible reaction kinetics, which allows direct consideration of free-energy dissipation. This also allows species to interact via products rather than just substrates, increasing the dynamical complexity, and allowing a more nuanced classification of interaction types to emerge. Using this model, we find that community diversity increases with substrate lability, because greater free-energy availability allows for faster generation of niches. Thus, more niches are generated in the time frame of community establishment, leading to higher final species diversity. We also find that allowing species to make use of near-to-equilibrium reactions increases diversity in a low free-energy regime. In such a regime, two new thermodynamic interaction types that we identify here reach comparable strengths to the conventional (competition and facilitation) types, emphasising the key role that thermodynamics plays in community dynamics. Our results suggest that accounting for realistic thermodynamic constraints is vital for understanding the dynamics of real-world microbial communities.

Author summaryThere is a growing interest in microbial communities due to their important role in biogeochemical cycling as well as plant and animal health. Although our understanding of thermodynamic constraints on individual cells is rapidly improving, the impact of these constraints on complex microbial communities remains largely unexplored theoretically and empirically. Here, we develop a new microbial community model which allows thermodynamic efficiency and entropy production to be calculated directly. We find that availability of substrates with greater free-energy allows for a faster rate of niche generation, leading to higher final species diversity. We also show that when the free-energy availability is low, species with reactions close to thermodynamic equilibrium are favoured, leading to more diverse and efficient communities. In addition to the conventional interaction types (competition and facilitation), our model reveals the existence of two novel interaction types mediated by products rather than substrates. Though the conventional interactions are generally the strongest, the novel interaction types are significant when free-energy availability is low. Our results suggest that non-equilibrium thermodynamics need to be considered when studying microbial community dynamics. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. A new method for the sampling and preservation of placental specimens in low-resource settings for the identification of P. falciparum and analysis of nucleic acids Collection, preservation, and shipment of histological specimens in low-resource settings is challenging. We present a novel method that achieved excellent preservation of placental specimens from rural Mali by using formalin fixation, ethanol dehydration, and long-term storage in a solar-powered freezer. Sample preservation success was 92%, permitting evaluation of current and past malaria infection, anemia, placental maturity, and inflammation. Using RNAscope(R) hybridization we were able to visualize cell-specific gene expression patterns in the formalin-fixed paraffin-embedded (FFPE) specimens. Additionally, our method entailed mirrored sampling from the two cut faces of a cotyledon, one for the FFPE workflows and the other for storage in RNAlater and RNA-seq. A theoretical model of mitochondrial ATP synthase deficiencies. The role of mitochondrial carriers. The m.8993T>G mutation of the mitochondrial MT-ATP6 gene is associated with NARP syndrome (Neuropathy, Ataxia and Retinitis Pigmentosa). The equivalent point mutation introduced in yeast Saccharomyces cerevisiae mitochondrial DNA considerably reduced the activity of ATP synthase and of cytochrome-c-oxidase preventing yeast growth on oxidative substrates. The over-expression of the mitochondrial oxodicarboxylate carrier (Odc1p) is able to rescue the growth on oxidative substrate by increasing the substrate-level phosphorylation of ADP coupled to conversion of -ketoglutarate (AKG) into succinate with an increase in Complex IV activity. Previous studies showed that equivalent point mutations in ATP synthase behave similarly and can be rescued by Odc1p overexpression and/or uncoupler of OXPHOS from ATP synthesis. In order to better understand the mechanism of ATP synthase mutation bypass, we developed a core model of mitochondrial metabolism based on AKG as respiratory substrate. We describe the different possible metabolite output and the ATP/O ratio values as a function of ATP synthase inhibition. Target Recognition in Tandem WW Domains: Complex Structures for Parallel and Antiparallel Ligand Orientation in h-FBP21 Tandem WW Protein-protein interactions often rely on specialized recognition domains, such as WW domains, which bind to specific proline-rich sequences. The specificity of these protein-protein interactions can be increased by tandem repeats, i.e. two WW domains connected by a linker. With a flexible linker, the WW domains can move freely with respect to each other. Additionally, the tandem WW domains can bind in two different orientations to their target sequences. This makes the elucidation of complex structures of tandem WW domains extremely challenging. Here, we identify and characterize two complex structures of the tandem WW domain of human formin-binding protein 21 and a peptide sequence from its natural binding partner, the core-splicing protein SmB/B. The two structures differ in the ligand orientation, and consequently also in the relative orientation of the two WW domains. We analyze and probe the interactions in the complexes by molecular simulations and NMR experiments. The workflow to identify the complex structures uses molecular simulations, density-based clustering and peptide docking. It is designed to systematically generate possible complex structures for repeats of recognition domains. These structures will help us to understand the synergistic and multivalency effects that generate the astonishing versatility and specificity of protein-protein interactions. Integrative analysis reveals aged clonal B cells, microenvironment and c-Myc activation in the origin of age-related lymphoma While cancer is an age-related disease, most studies focus on genetically engineered younger mouse models. Here, we uncover how cancer develops as a consequence of the naturally aged immune system in mice. B-cell lymphoma frequently occurs in aged mice and is associated with increased cell size, splenomegaly, and a novel clonal B-cell population. Age-emergent B cells clonally expand outside of germinal centers driven by somatic mutations, activated c-Myc and hypermethylated promoters, and both genetically and epigenetically recapitulate human follicular and diffuse-large B-cell lymphomas. Mechanistically, mouse cancerous B cells originate from age-associated B cells, which are atypical memory B cells. Age-associated B cells secrete a spectrum of proinflammatory cytokines and activate paracrinally the expression of c-Myc in surrounding B cells. Although clonal B cells are a product of an aging microenvironment, they evolve being self-sufficient and support malignancy when transferred into young mice. Inhibition of mTOR and c-Myc attenuates premalignant changes in B cells during aging and emerges as a therapeutic strategy to delay the onset of age-related lymphoma. Together, we uncover how aging generates cancerous B cells, characterize a model that captures the origin of spontaneous cancer during aging and identify interventions that may postpone age-associated lymphoma. An interplay between cellular growth and atypical fusion defines morphogenesis of a modular glial niche Neural stem cells (NSCs) are embedded in a multi-layered, intricate cellular microenvironment supporting their activity, the niche. Whilst shape and function are inseparable, the morphogenetic aspects of niche development are poorly understood. Here, we use the formation of the glial network of a NSC niche to investigate acquisition of architectural complexity. Cortex glia (CG) in Drosophila regulate neurogenesis and build a reticular structure around NSCs. We first show that individual CG cells grow tremendously to ensheath several NSC lineages, eventually spanning the entire tissue while partitioning the NSC population. Elaborate proliferative mechanisms convert these cells into syncytia rich in cytoplasmic bridges. Unexpectedly, CG syncytia further undergo homotypic cell-cell fusion, relying on defined molecular players of cell fusion such as cell surface receptors and actin regulators. Exchange of cellular components is however dynamic in space and time, a previously unreported unique mechanism. This atypical cell fusion remodels cellular borders, restructuring the CG syncytia. Ultimately, the coordination of growth and fusion builds the multi-level architecture of the niche, and creates a modular, spatial partition of the NSC population. Our findings provide novel insights into how a niche forms and organises while developing intimate contacts with a stem cell population. Maternal behavioural thermoregulation facilitated evolutionary transitions from egg laying to live birth Live birth is a key innovation that has evolved from egg laying ancestors over 100 times in reptiles. However, egg-laying lizards and snakes often have preferred body temperatures that are lethal to developing embryos and should prevent egg retention: how has viviparity repeatedly evolved in the face of this pervasive mismatch? We resolve this paradox by conducting phylogenetic analyses of adult and embryo thermal preferences across 224 species. Thermal mismatches between mothers and offspring are widespread but resolved by gravid females down-regulating their body temperature towards the thermal optimum of embryos. Importantly, this thermoregulatory behaviour evolved in ancestral egg-laying species before the evolution of live birth. Maternal thermoregulatory behaviour therefore bypasses the constraints imposed by a slowly evolving thermal physiology and is likely to have been a key requirement for repeated transitions to live birth. Cytokinin oxidase/dehydrogenase family genes exhibit functional divergence and overlap in rice growth and development, especially in control of tillering Cytokinins play key roles in plant growth and development; hence, cytokinin biosynthesis and degradation have been extensively studied. Cytokinin oxidase/dehydrogenases (CKXs) are a group of enzymes that regulate oxidative cleavage to maintain cytokinin homeostasis. In rice, 11 OsCKX genes have been identified to date; however, most of their functions remain unknown. Here, we comprehensively analyzed the expression patterns and functions of OsCKX genes. Using CRISPR/Cas9 technology, we constructed mutants of all OsCKX genes to determine the functions of OsCKXs in rice development. The results revealed that the single osckx and higher-order osckx4 osckx9 mutant lines showed functional overlap and subfunctionalization. Notably, the osckx1 osckx2 and osckx4 osckx9 double mutants displayed contrasting phenotypic changes in tiller number and panicle size compared to the wild type. Moreover, we identified several genes with significantly altered expression in osckx4 and osckx9 single and double mutant plants. Many differentially expressed genes were found to be associated with auxin and cytokinin pathways. Additionally, the cytokinins in osckx4 osckx9 mutants were increased compared to the wild type. Overall, our findings provide new insights into the functions of OsCKX genes in rice growth and may be used as a foundation for future studies aimed at improving rice yield. Key words: Cytokinin, expression pattern, OsCKX, panicle, phenotype, rice, tiller

HighlightThe osckx4 osckx9 double mutant had a significantly greater number of tillers, whereas the osckx1 osckx2 double mutant showed the opposite phenotypic change, compared to the wild type Royal knifefish generate powerful suction feeding through large neurocranial elevation and high epaxial muscle power Suction feeding in ray-finned fishes involves powerful buccal cavity expansion to accelerate water and food into the mouth. Previous XROMM studies in largemouth bass (Micropterus salmoides), bluegill sunfish (Lepomis macrochirus), and channel catfish (Ictalurus punctatus) have shown that more than 90% of suction power in high performance strikes comes from the axial musculature. Thus, the shape of the axial muscles and skeleton may impact suction feeding mechanics. Royal knifefish (Chitala blanci) have an unusual postcranial morphology, with a ventrally flexed vertebral column and relatively large mass of epaxial muscle. Based on their body shape, we hypothesized that royal knifefish would generate high power strikes by utilizing large neurocranial elevation, vertebral column extension, and epaxial shortening. As predicted, C. blanci generated high suction expansion power compared to the other three species studied to date (up to 160 W), which was achieved by increasing both the rate of volume change and the intraoral subambient pressure. The large epaxial muscle (25% of body mass) shortened at high velocities to produce large neurocranial elevation and vertebral extension (up to 41 deg, combined), as well as high muscle mass-specific power (up to 800 W kg-1). For the highest power strikes, axial muscles generated 95% of the power, and 64% of the axial muscle mass consisted of the epaxial muscles. The epaxial-dominated suction expansion of royal knifefish supports our hypothesis that postcranial morphology may be a strong predictor of suction feeding biomechanics.

SUMMARY STATEMENTRoyal knifefish rely on their distinct postcranial morphology--with a curved vertebral column and large dorsal body muscles--to produce large neurocranial elevation and powerful suction feeding. Immunogenicity of SARS-CoV-2 trimetric spike protein associated to Poly(I:C) plus Alum The SARS-CoV-2 pandemic has had a social and economic impact worldwide, and vaccination is an efficient strategy for diminishing those damages. New adjuvant formulations are required for the high vaccine demands, especially adjuvant formulations that induce a Th1 phenotype. Herein we assess a vaccination strategy using a combination of Alum and polyinosinic:polycytidylic acid (Poly(I:C)) adjuvants plus the SARS-CoV-2 spike protein in a prefusion trimeric conformation by an intradermal (ID) route. We found high levels of IgG anti-spike antibodies in the serum by enzyme linked immunosorbent assay (ELISA) and high neutralizing titers against SARS-CoV-2 in vitro by neutralization assay, after one or two boosts. By evaluating the production of IgG subtypes, as expected, we found that formulations containing Poly(I:C) induced IgG2a whereas Alum did not. The combination of these two adjuvants induced high levels of both IgG1 and IgG2a. In addition, cellular immune responses of CD4+ and CD8+ T cells producing interferon-gamma were equivalent, demonstrating that the Alum + Poly(I:C) combination supported a Th1 profile. Based on the high neutralizing titers, we evaluated B cells in the germinal centers, which are specific for receptor-binding domain (RBD) and spike, and observed that more positive B cells were induced upon the Alum + Poly(I:C) combination. Moreover, these B cells produced antibodies against both RBD and non-RBD sites. We also studied the impact of this vaccination preparation (spike protein with Alum + Poly(I:C)) in the lungs of mice challenged with inactivated SARS-CoV-2 virus. We found a production of IgG, but not IgA, and a reduction in neutrophil recruitment in the bronchoalveolar lavage fluid (BALF) of mice, suggesting that our immunization scheme reduced lung inflammation. Altogether, our data suggest that Alum and Poly(I:C) together is a possible adjuvant combination for vaccines against SARS-CoV-2 by the intradermal route. NECAB2 participates in an endosomal pathway of mitochondrial stress response at striatal synapses Synaptic signaling depends on ATP generated by mitochondria. Due to extensive connectivity, the striatum is especially vulnerable to mitochondrial dysfunction and thus requires efficient mitochondrial quality control and repair. We found that global knockout of the neuronal calcium-binding protein 2 (NECAB2) in the mouse causes loss of striatal synapses and behavioral phenotypes related to striatal dysfunction such as reduced motivation and altered sensory gating. Striatal mitochondria from Necab2 knockout mice are more abundant and smaller. They are characterized by increased respiration and superoxide production resulting in oxidative stress. This accumulation of dysfunctional mitochondria is caused by a defective assembly of mitochondria with early endosomes in a pathway that involves the small GTPase Rab5 and its guanine nucleotide exchange factor Alsin/ALS2. NECAB2 therefore participates in an endosomal pathway of mitochondrial stress response and repair important for striatal function. Regulation of the PKD2 channel function by TACAN Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in membrane receptor PKD1 or cation channel PKD2. TACAN (also named TMEM120A), recently reported as an ion channel in neuron cells for mechano and pain sensing, is also distributed in diverse non-neuronal tissues such as kidney, heart and intestine, suggesting its involvement in other functions. In this study, we found that TACAN is in complex with PKD2 in native renal cell lines. Using the two-electrode voltage clamp in Xenopus oocytes we found that TACAN inhibited the channel activity of PKD2 gain-of-function mutant F604P. The first and last transmembrane domains of TACAN were found to interact with the PKD2 C-and N-terminal portions, respectively. We showed that the TACAN N-terminus acted as a blocking peptide and that TACAN inhibits the PKD2 function through the PKD2/TACAN binding. By patch clamping in mammalian cells, we found that TACAN inhibits both the single channel conductance and open probability of PKD2 and mutant F604P. PKD2 co-expressed with TACAN, but not PKD2 alone, exhibited pressure sensitivity. Furthermore, we also found that TACAN aggravates PKD2-dependent tail curvature and pronephric cysts in larval zebrafish, in support of the in vitro inhibitory effects of TACAN. In summary, this study revealed that TACAN acts as a PKD2 inhibitor and mediates mechano sensitivity of the PKD2/TACAN channel complex. The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. The extracellular loops of OmpA control the slow rate of in vitro folding Outer membrane proteins are all beta barrels and these barrels have a variety of well-documented loop conformations. Here we test the effect of three different loop types on outer membrane protein A (OmpA) folding. We designed twelve 5-residue loops and experimentally tested the effect of replacing the long loops of outer membrane protein OmpA with the designed loops. Our studies succeeded in creating the smallest known outer membrane barrel. We find that significant changes in OmpA loops do not have a strong overall effect on OmpA folding. However, when decomposing folding into a fast rate and a slow rate we find that changes in loops strongly affect the slow rate of OmpA folding. Extracellular loop types with higher levels of hydrogen bonds had more instances of increasing the slow folding rate and extracellular loop types with low levels of hydrogen bonds had more instances of decreasing the slow folding rate. Having the slow rate affected by loop composition is consistent with the slow rate being associated with the insertion step of outer membrane protein folding. Optimization of avian perching manoeuvres Perching at speed is amongst the most challenging flight behaviours that birds perform1,2, and beyond the capability of current autonomous vehicles. Smaller birds may touchdown by hovering3-8, but larger birds typically swoop upward to perch1,2 - presumably because the adverse scaling of their power margin prohibits slow flapping flight9, and because swooping transfers excess kinetic to potential energy1,2,10. Perching is risky in larger birds6,11, demanding precise control of velocity and pose12-15, but it is unknown how they optimize this challenging manoeuvre. More generally, whereas cruising flight behaviours such as migration and commuting are adapted to minimize cost-of-transport or time-of-flight16, the optimization of unsteady flight manoeuvres remains largely unexplored7,17. Here we show that swooping minimizes neither the time nor energy required to perch safely in Harris hawks Parabuteo unicinctus, but instead minimizes the distance flown under hazardous post-stall conditions. By combining motion capture data from 1,563 flights with flight dynamics modelling, we found that the birds choice of where to transition from powered dive to unpowered climb minimizes the distance from the landing perch over which very high lift coefficients are required. Time and energy are therefore invested to maintain the control authority needed to execute a safe landing, rather than being minimized continuously as in technical applications of autonomous perching under nonlinear feedback control13 and deep reinforcement learning18,19. Naive birds learn this behaviour on-the-fly, so our findings suggest an alternative reward function for reinforcement learning of autonomous perching in air vehicles. Cluster-specific gene markers enhance Shigella and Enteroinvasive Escherichia coli in silico serotyping Shigella and enteroinvasive Escherichia coli (EIEC) cause human bacillary dysentery with similar invasion mechanisms and share similar physiological, biochemical and genetic characteristics. The ability to differentiate Shigella and EIEC from each other is important for clinical diagnostic and epidemiologic investigations. The existing genetic signatures may not discriminate between Shigella and EIEC. However, phylogenetically, Shigella and EIEC strains are composed of multiple clusters and are different forms of E. coli. In this study, we identified 10 Shigella clusters, 7 EIEC clusters and 53 sporadic types of EIEC by examining over 17,000 publicly available Shigella/EIEC genomes. We compared Shigella and EIEC accessory genomes to identify the cluster-specific gene markers or marker sets for the 17 clusters and 53 sporadic types. The gene markers showed 99.63% accuracy and more than 97.02% specificity. In addition, we developed a freely available in silico serotyping pipeline named Shigella EIEC Cluster Enhanced Serotype Finder (ShigEiFinder) by incorporating the cluster-specific gene markers and established Shigella/EIEC serotype specific O antigen genes and modification genes into typing. ShigEiFinder can process either paired end Illumina sequencing reads or assembled genomes and almost perfectly differentiated Shigella from EIEC with 99.70% and 99.81% cluster assignment accuracy for the assembled genomes and mapped reads respectively. ShigEiFinder was able to serotype over 59 Shigella serotypes and 22 EIEC serotypes and provided a high specificity with 99.40% for assembled genomes and 99.38% for mapped reads for serotyping. The cluster markers and our new serotyping tool, ShigEiFinder (https://github.com/LanLab/ShigEiFinder), will be useful for epidemiologic and diagnostic investigations.

Impact statementThe differentiation of Shigella strains from enteroinvasive E. coli (EIEC) is important for clinical diagnosis and public health epidemiologic investigations. The similarities between Shigella and EIEC strains make this differentiation very difficult as both share common ancestries within E. coli. However, Shigella and EIEC are phylogenetically separated into multiple clusters, making high resolution separation using cluster specific genomic markers possible. In this study, we identified 17 Shigella or EIEC clusters including five that were newly identified through examination of over 17,000 publicly available Shigella and EIEC genomes. We further identified an individual or a set of cluster-specific gene markers for each cluster using comparative genomic analysis. These markers can then be used to classify isolates into clusters and were used to develop an in silico pipeline, ShigEiFinder (https://github.com/LanLab/ShigEiFinder) for accurate differentiation, cluster typing and serotyping of Shigella and EIEC from Illumina sequencing reads or assembled genomes. This study will have broad application from understanding the evolution of Shigella/EIEC to diagnosis and epidemiology.

Data summarySequencing data have been deposited at the National Center for Biotechnology Information under BioProject number PRJNA692536.

RepositoriesRaw sequence data are available from NCBI under the BioProject number PRJNA692536. Inhibitory proteins block substrate access by occupying the active site cleft of Bacillus subtilis intramembrane metalloprotease SpoIVFB Intramembrane proteases function in numerous signaling pathways that impact health, but knowledge about regulation of intramembrane proteolysis is limited. We examined inhibition of intramembrane metalloprotease SpoIVFB by proteins BofA and SpoIVFA. We found that BofA residues in and near a predicted transmembrane segment are required for SpoIVFB inhibition, and cross-linking experiments indicated that this transmembrane segment occupies the SpoIVFB active site cleft. SpoIVFA is also required for SpoIVFB inhibition. The inhibitory proteins block access of the substrate N-terminal Proregion to the membrane-embedded SpoIVFB active site, based on additional cross-linking experiments; however, the inhibitory proteins did not prevent interaction between the substrate C-terminal region and the SpoIVFB soluble domain. A structural model was built of SpoIVFB in complex with BofA and parts of SpoIVFA and substrate, using partial homology and constraints from cross-linking and co-evolutionary analyses. The model predicts that conserved BofA residues interact to stabilize a transmembrane segment and a membrane-embedded C-terminal region. SpoIVFA is predicted to bridge the BofA C-terminal region and SpoIVFB, forming a membrane-embedded inhibition complex. Our results reveal a novel mechanism of intramembrane protease inhibition with clear implications for relief from inhibition in vivo and design of inhibitors as potential therapeutics. "Adopt a Microorganism" methodology enhances the discourse richness of high school students in both hybrid learning and remote education The Internet has changed the way teachers and students access information and build knowledge. The recent COVID-19 pandemic has created challenges for both teachers and students and a demand for new methodologies of remote learning. In the life sciences, mixing online content with practical activities represents an even greater challenge. In microbiology, the implementation of an active teaching methodology, the #Adopt project, based on the social network Facebook(R), represents an excellent option for connecting remote education with classroom activities. In 2020, the version applied in high school, "Adopt a Microorganism", was adapted to meet the demands of emergency remote education owing to the suppression of face-to-face activities caused by the pandemic. In the present study, we assessed how the change in methodology impacted the discourse richness of students from high school integrated with technical education in the Business Administration program of the Federal Institute of Sao Paulo, Sorocaba Campus. Three questionnaires related to the groups of microorganisms (Archaea, Bacteria, Virus, Fungi, and Protozoan) were applied. The students responses in the 2019 and 2020 classes were compared concerning content richness and multiplicity of concepts through the application of the Shannon diversity index, an approach that is generally used to assess biodiversity in different environments. The observed results suggest that remote learning provided students with a conceptual basis and richness of content equivalent to that achieved by students subjected to the hybrid teaching model. In conclusion, this study suggests that the #Adopt project methodology increases students discourse richness in microbiology even without face-to-face traditional classes. A novel sodium signaling complex regulates uterine activity Depolarization of the myometrial smooth muscle cell (MSMC) resting membrane potential is necessary for the transition of the uterus from a quiescent state to a contractile state. The molecular mechanisms involved in this transition are not completely understood. Here, we report a novel coupled system between the Na+-activated K+ channel (SLO2.1) and the non-selective Na+ leak channel (NALCN) which determines the MSMC membrane potential. We show that SLO2.1 currents are activated by an inward Na+ leak current carried by the NALCN channel leading to MSMC hyperpolarization. These results show an unanticipated role for the Na+ leak currents in activating a negative feedback system countering the excitable effects of Na+ currents. This is a novel role for the NALCN channel in which Na+ acts as an intracellular signaling molecule. In fact, we report here that the net effect of Na+ entry through NALCN channels is a hyperpolarization of the MSMCs plasma membrane because of the activation of SLO2.1 K channel. Importantly, we also report that a decrease in NALCN/SLO2.1 activity triggers both Ca2+ entries through VDCCs, promoting myometrial contraction. Consistently, with a functional coupling, our data show that NALCN and SLO2.1 are in proximity to one another in human MSMCs. We propose that the spatial arrangement of SLO2.1 and NALCN permits these channels to functionally interact in order to regulate human MSMC membrane potential and cell excitability to modulate uterine contractile activity. Spatio-temporal, optogenetic control of gene expression in organoids Organoids derived from stem cells become increasingly important to study human development and to model disease. However, methods are needed to control and study spatio-temporal patterns of gene expression in organoids. To this aim, we combined optogenetics and gene perturbation technologies to activate or knock-down RNA of target genes, at single-cell resolution and in programmable spatio-temporal patterns. To illustrate the usefulness of our approach, we locally activated Sonic Hedgehog (SHH) signaling in an organoid model for human neurodevelopment. High-resolution spatial transcriptomic and single-cell analyses showed that this local induction was sufficient to generate stereotypically patterned organoids in three dimensions and revealed new insights into SHHs contribution to gene regulation in neurodevelopment.

With this study, we propose optogenetic perturbations in combination with spatial transcriptomics as a powerful technology to reprogram and study cell fates and tissue patterning in organoids. Frameshift and wild-type proteins are highly similar because the genetic code and genomes were optimized for frameshift tolerance Frameshift protein sequences encoded by alternative reading frames of coding genes have been considered meaningless, and frameshift mutations have been considered of little importance for the molecular evolution of coding genes and proteins. However, functional frameshifts have been found widely existing. It was puzzling how a frameshift protein kept its structure and functionality while its amino-acid sequence was changed substantially. Here we show that frame similarities between frameshifts and wild types are higher than random similarities and are defined at the genetic code, gene, and genome levels. In the standard genetic code, frameshift codon substitutions are more conservative than random substitutions. The frameshift tolerability of the standard genetic code ranks in the top 2.0-3.5% of alternative genetic codes, showing that the genetic code is nearly optimal for frameshift tolerance. Furthermore, frameshift-resistant codons (codon pairs) appear more frequently than expected in many genes and certain genomes, showing that the frameshift optimality is reflected not only in the genetic code but more importantly, in its allowance of further optimizing the frameshift tolerance of a particular gene or genome, which shed light on the role of frameshift mutations in molecular and genomic evolution. Evaluation of methods for estimating coalescence times using ancestral recombination graphs The ancestral recombination graph (ARG) is a structure that describes the joint genealogies of sampled DNA sequences along the genome. Recent computational methods have made impressive progress towards scalably estimating whole-genome genealogies. In addition to inferring the ARG, some of these methods can also provide ARGs sampled from a defined posterior distribution. Obtaining good samples of ARGs is crucial for quantifying statistical uncertainty and for estimating population genetic parameters such as effective population size, mutation rate, and allele age. Here, we use standard neutral coalescent simulations to benchmark the estimates of pairwise coalescence times from three popular ARG inference programs: ARGweaver, Relate, and tsinfer+tsdate. We compare 1) the true coalescence times to the inferred times at each locus; 2) the distribution of coalescence times across all loci to the expected exponential distribution; 3) whether the sampled coalescence times have the properties expected of a valid posterior distribution. We find that inferred coalescence times at each locus are most accurate in ARGweaver, and often more accurate in Relate than in tsinfer+tsdate. However, all three methods tend to overestimate small coalescence times and underestimate large ones. Lastly, the posterior distribution of ARGweaver is closer to the expected posterior distribution than Relates, but this higher accuracy comes at a substantial trade-off in scalability. The best choice of method will depend on the number and length of input sequences and on the goal of downstream analyses, and we provide guidelines for the best practices.