The multi-level phenotypic impact of synonymous substitutions: heterologous gene expression in human cells. Redundancy in the genetic code allows for differences in transcription and/or translation efficiency between sequences carrying synonymous polymorphisms, potentially leading to phenotypic changes. It is commonly admitted that the evolution codon usage bias ("CUB", the over-representation of certain codons in a genome, a gene or in positions along a gene) are driven by a combination of neutral and selective processes, but their relative contribution is a matter of debate, especially in mammals. Particularly, integrative studies quantifying the phenotypic impact of CUB at different molecular and cellular levels are lacking. Here we report a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells.

Six synonymous versions of the shble antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was assessed by means of: i) mRNA-to-DNA and protein-to-mRNA ratios for each shble version; ii) cellular fluorescence, using flow cytometry, as a proxy for single-cell level expression; and iii) real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness.

We show that differences in CUB strongly impact the molecular and the cellular phenotype: i) they result in large differences in mRNA and in protein levels, as well in mRNA-to-protein ratio; ii) they introduce splicing events not predicted by current algorithms; iii) they lead to reproducible phenotypic heterogeneity; iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression.

We interpret that CUB modulate mRNA availability and suitability for translation in human cells, leading to differences in protein levels and eventually eliciting phenotypic differences.

AUTHOR SUMMARYThe genetic code is redundant, with several codons encoding for the same amino acid. These synonymous codons are not used with equal frequencies. Instead, codon usage bias (CUB) varies between species, genes, and even positions along a gene. At each of these levels, CUB are shaped by the overall balance between mutational biases and selection forces. To shed light on the molecular mechanisms underlying molecular and organism phenotypes, integrative studies quantifying the phenotypic impact of CUB at different levels of biological integration are necessary. Here, we monitored the multiscale changes induced by six synonymous versions of an antibiotic resistance gene independently expressed in a human cell line. We show that: 1. both mRNA levels, protein levels, and their ratios are affected by CUB; 2. potential effects on mRNA (splicing and availability for translation) seem to constitute the main level of action; 3. cell fitness is severely impacted by a trade-off between the burden of heterologous expression and the benefits of antibiotic resistance. This integrative study provides new insights on translation regulation and the associated phenotypic impact in human cells, associated to CU.
Activation of the type 3 secretion system of enteropathogenic E. coli leads to remodeling of its membrane composition and function The cell envelope of gram-negative bacteria is a complex structure, essential for bacterial survival and for resistance to many antibiotics. Channels that cross the bacterial envelope and the host cell membrane form secretion systems that are activated upon attachment to host, enabling bacteria to inject effector molecules into the host cell, required for bacterial-host interaction. The type III secretion system (T3SS) is critical for the virulence of several pathogenic bacteria, including enteropathogenic E. coli (EPEC). The EPEC T3SS activation is associated with repression of carbon storage regulator (CsrA), resulting in gene expression remodeling, which is known to affect EPEC central carbon metabolism and contributes to the adaptation to a cell-adherent lifestyle in a poorly understood manner. We reasoned that the changes in bacterial envelope upon attachment to host and the activation of a secretion system may involve a modification of the lipid composition of bacterial envelope. Accordingly, we performed a lipidomics analysis on mutant strains that simulate T3SS activation. We saw a shift in glycerophospholipid metabolism towards the formation of lysophospholipids, attributed to corresponding upregulation of the phospholipase pldA and the acyltransferase ygiH upon T3SS activation in EPEC. We also detected a shift from menaquinones and ubiquinones to undecaprenyl lipids, concomitant to abnormal synthesis of O-antigen. The remodeling of lipid metabolism is mediated by CsrA and associated with increased bacteria cell size and Zeta potential, and a corresponding alteration in EPEC permeability to vancomycin, increasing the sensitivity of T3SS-activated strains and of adherent wild type EPEC to the antibiotic.

ImportanceThe characterization of EPEC membrane lipid metabolism upon attachment to host is an important step towards a better understanding the shift of EPEC, a notable human pathogen, from a planktonic to adherent life style. It may also apply to other pathogenic bacteria that use this secretion system. We predict that upon attachment to host cells the lipid remodeling upon T3SS activation contributes to bacterial fitness and promotes host colonization, and show that it is associated with increased cell permeability and higher sensitivity to vancomycin.

To the best of our knowledge, this is the first demonstration of a bacterial lipid remodeling due to activation of a secretion system.
An enhanced inverted encoding model for neural reconstructions Here we present a more interpretable and flexible approach for reconstructing the contents of perception, attention, and memory from neuroimaging data. Our enhanced inverted encoding model (eIEM) incorporates methodological improvements including proper accounting of population-level tuning functions and a trial-by-trial prediction error-based metric where reconstruction quality is measured in meaningful units. Improved flexibility is further gained via eIEMs novel goodness-of-fit feature: for trial-by-trial reconstructions, goodness-of-fits are obtained independently (non-circularly) to prediction error and can be applied to any IEM procedure or decoding metric, resulting in improved reconstruction quality and brain-behavior correlations. We validate eIEM from methodological principles, simulated neuroimaging datasets, and three pre-existing fMRI datasets spanning perception, attention, and working memory. Notably, eIEM is easy to apply and broadly accessible - our Python package (https://pypi.org/project/inverted-encoding) implements eIEM in one line of code - and is easily modifiable to compare performance metrics and/or scale up to more complex models. Towards a cumulative science of vocal markers of autism: a cross-linguistic meta-analysis-based investigation of acoustic markers in American and Danish autistic children Acoustic atypicalities in speech production are argued to be potential markers of clinical features in Autism Spectrum Disorder (ASD). A recent meta-analysis highlighted shortcomings in the field, in particular small sample sizes and study heterogeneity (Fusaroli et al., 2017). We showcase a cumulative (i.e., explicitly building on previous studies both conceptually and statistically) yet self-correcting (i.e., critically assessing the impact of cumulative statistical techniques) approach to prosody in ASD to overcome these issues.

We relied on the recommendations contained in the meta-analysis to build and analyze a cross-linguistic corpus of multiple speech productions in 77 autistic and 72 neurotypical children and adolescents (>1000 recordings in Danish and US English). We used meta-analytically informed and skeptical priors, with informed priors leading to more generalizable inference. We replicated findings of a minimal cross-linguistically reliable distinctive acoustic profile for ASD (higher pitch and longer pauses) with moderate effect sizes. We identified novel reliable differences between the two groups for normalized amplitude quotient, maxima dispersion quotient, and creakiness. However, the differences were small, and there is likely no one acoustic profile characterizing all autistic individuals. We identified reliable relations of acoustic features with individual differences (age, gender), and clinical features (speech rate and ADOS sub-scores).

Besides cumulatively building our understanding of acoustic atypicalities in ASD, the study shows how to use systematic reviews and meta-analyses to guide the design and analysis of follow-up studies. We indicate future directions: larger and more diverse cross-linguistic datasets, focus on heterogeneity, self-critical cumulative approaches and open science.

Lay SummaryAutistic individuals are reported to speak in distinctive ways. Distinctive vocal production can affect social interactions and social development and could represent a noninvasive way to support the assessment of ASD. We systematically checked whether acoustic atypicalities highlighted in previous articles could be actually found across multiple recordings and two languages. We find a minimal acoustic profile of ASD: higher pitch, longer pauses, increased hoarseness and creakiness of the voice. However, there is much individual variability (by age, sex, language, and clinical characteristics). This suggests that the search for one common "autistic voice" might be naive and more fine-grained approaches are needed. GCN5L1 promotes diastolic dysfunction by inhibiting cardiac pyruvate oxidation Left ventricular diastolic dysfunction is a structural and functional condition that precedes the development of heart failure with preserved ejection fraction (HFpEF). The etiology of diastolic dysfunction includes alterations in fuel substrate metabolism that negatively impact cardiac bioenergetics, and may precipitate the eventual transition to heart failure. To date, the molecular mechanisms that regulate early changes in fuel metabolism leading to diastolic dysfunction remain unclear. In this report, we use a diet-induced obesity model and quantitative acetylproteomics in aged mice to show that inhibitory lysine acetylation of the pyruvate dehydrogenase (PDH) complex promotes energetic deficits and diastolic dysfunction in mouse hearts. Cardiomyocyte-specific deletion of the mitochondrial lysine acetylation regulatory protein GCN5L1 prevented hyperacetylation of the PDH complex subunit PDHA1, allowing aged obese mice to continue using pyruvate as a bioenergetic substrate in the heart. Our findings suggest that changes in mitochondrial protein lysine acetylation represent a key metabolic component of diastolic dysfunction that precedes the development of heart failure. Inability to switch from ARID1A-BAF to ARID1B-BAF impairs exit from pluripotency and commitment towards neural crest formation in ARID1B-related neurodevelopmental disorders The BAF complex modulates chromatin accessibility. Specific BAF configurations have functional consequences, and subunit switches are essential for cell differentiation. ARID1B and its paralog ARID1A encode for mutually exclusive BAF subunits. De novo ARID1B haploinsufficient mutations cause a neurodevelopmental disorder spectrum, including Coffin-Siris syndrome, which is characterized by neurological and craniofacial features. Here, we reprogrammed ARID1B+/- Coffin-Siris patient-derived skin fibroblasts into iPSCs and modeled cranial neural crest cell (CNCC) formation. We discovered that ARID1B is active only during the first stage of this process, coinciding with neuroectoderm specification, where it is part of a lineage-specific BAF configuration (ARID1B-BAF), which includes SMARCA4 and nine additional subunits. ARID1B-BAF acts as a gatekeeper, ensuring exit from pluripotency and lineage commitment, by attenuating NANOG, SOX2 and thousands of enhancers directly regulated by these two pluripotency factors at the iPSC stage.

In iPSCs, these enhancers are maintained active by an ARID1A-containing BAF. At the onset of differentiation, cells transition from ARID1A-BAF to ARID1B-BAF, eliciting attenuation of the NANOG/SOX2 networks, and triggering pluripotency exit. Coffin-Siris patient cells fail to perform the ARID1A/ARID1B switch and maintain ARID1A-BAF at pluripotency enhancers throughout all stages of CNCC formation. This leads to a persistent and aberrant SOX2 and NANOG activity, which impairs CNCC formation. In fact, despite showing the typical neural crest signature (TFAP2A+, SOX9+), ARID1B-haploinsufficient CNCCs are also NANOG/OCT4-positive, in stark contrast with the ARID1B-wt CNCCs, which are NANOG/OCT4-negative.

These findings suggest a connection between ARID1B mutations, neuroectoderm formation, and a pathogenic mechanism for Coffin-Siris syndrome. Inter-chromosomal linkage disequilibrium and linked fitness cost loci associated with selection for herbicide resistance The adaptation of weedy plants to herbicide is both a significant problem in agriculture and a model for the study of rapid adaptation under regimes of strong selection. Despite recent advances in our understanding of simple genetic changes that lead to resistance, a significant gap remains in our knowledge of resistance controlled by many loci and the evolutionary factors that influence the maintenance of resistance over time. Here, we perform a multi-level analysis involving whole genome sequencing and assembly, resequencing and gene expression analysis to both uncover putative loci involved in nontarget herbicide resistance and to examine evolutionary forces underlying the maintenance of resistance in natural populations. We found loci involved in herbicide detoxification, stress sensing, and alterations in the shikimate acid pathway to be under selection, and confirmed that detoxification is responsible for glyphosate resistance using a functional assay. Furthermore, we found interchromosomal linkage disequilibrium (ILD), most likely associated with epistatic selection, to influence NTSR loci found on separate chromosomes thus potentially mediating resistance through generations. Additionally, by combining the selection screen, differential expression and LD analysis, we identified fitness cost loci that are strongly linked to resistance alleles, indicating the role of genetic hitchhiking in maintaining the cost. Overall, our work strongly suggests that NTSR glyphosate resistance in I. purpurea is conferred by multiple genes which are maintained through generations via ILD, and that the fitness cost associated with resistance in this species is a by-product of genetic-hitchhiking. Evidence for a hydrogen sulfide-sensing E3 ligase in yeast In yeast, control of sulfur amino acid metabolism relies upon Met4, a transcription factor which activates the expression of a network of enzymes responsible for the biosynthesis of cysteine and methionine. In times of sulfur abundance, the activity of Met4 is repressed via ubiquitination by the SCFMet30 E3 ubiquitin ligase, but the mechanism by which the F-box protein Met30 senses sulfur status to tune its E3 ligase activity remains unresolved. Herein, we show that Met30 responds to flux through the transsulfuration pathway to regulate the MET gene transcriptional program. In particular, Met30 is responsive to the biological gas hydrogen sulfide, which is sufficient to induce ubiquitination of Met4 in vivo. Additionally, we identify important cysteine residues in Met30s WD-40 repeat region that sense the availability of sulfur in the cell. Our findings reveal how SCFMet30 dynamically senses the flow of sulfur metabolites through the transsulfuration pathway to regulate synthesis of these special amino acids. Cytosolic aggregation of mitochondrial proteins disrupts cellular homeostasis by stimulating the aggregation of other proteins Mitochondria are organelles with their own genomes, but they rely on the import of nuclear-encoded proteins that are translated by cytosolic ribosomes. Therefore, it is important to understand whether failures in the mitochondrial uptake of these nuclear-encoded proteins can cause proteotoxic stress and identify response mechanisms that may counteract it. Here, we report that upon impairments in mitochondrial protein import, high-risk precursor and immature forms of mitochondrial proteins form aberrant deposits in the cytosol. These deposits then cause further cytosolic accumulation and consequently aggregation of other mitochondrial proteins and disease-related proteins, including -synuclein and amyloid {beta}. This aggregation triggers a cytosolic protein homeostasis imbalance that is accompanied by specific molecular chaperone responses at both the transcriptomic and protein levels. Altogether, our results provide evidence that mitochondrial dysfunction, specifically protein import defects, contributes to impairments in protein homeostasis, thus revealing a possible molecular mechanism by which mitochondria are involved in neurodegenerative diseases. Conformist social learning leads to self-organised prevention against adverse bias in risky decision making Given the ubiquity of potentially adverse behavioural bias owing to myopic trial-and-error learning, it seems paradoxical that improvements in decision-making performance through conformist social learning, a process widely considered to be bias amplification, still prevail in animal collective behaviour. Here we show, through model analyses and large-scale interactive behavioural experiments with 585 human subjects, that conformist influence can indeed promote favourable risk taking in repeated experience-based decision making, even though many individuals are systematically biased towards adverse risk aversion. Although strong positive feedback conferred by copying the majoritys behaviour could result in unfavourable informational cascades, our differential equation model of collective behavioural dynamics identified a key role for increasing exploration by negative feedback arising when a weak minority influence undermines the inherent behavioural bias. This  collective behavioural rescue, emerging through coordination of positive and negative feedback, highlights a benefit of collective learning in a broader range of environmental conditions than previously assumed and resolves the ostensible paradox of adaptive collective behavioural flexibility under conformist influences. Clonal evolution in serially passaged Cryptococcus neoformans x deneoformans hybrids reveals a heterogenous landscape of genomic change Cryptococcus neoformans x deneoformans hybrids (also known as serotype AD hybrids) are basidiomycete yeasts that are common in a clinical setting. Like many hybrids, the AD hybrids are largely locked at the F1 stage and are mostly unable to undergo normal meiotic reproduction. However, these F1 hybrids, which display a high ([~]10%) sequence divergence are known to genetically diversify through mitotic recombination and aneuploidy, and this diversification may be adaptive. In this study, we evolved a single AD hybrid genotype in six diverse environments by serial passaging and then used genome resequencing of evolved clones to determine evolutionary mechanisms of adaptation. The evolved clones generally increased fitness after passaging, accompanied by an average of 3.3 point mutations, 2.9 loss of heterozygosity (LOH) events, and 0.7 trisomic chromosomes per clone. LOH occurred through nondisjunction of chromosomes, crossing over consistent with break-induced replication, and gene conversion, in that order of prevalence. The breakpoints of these recombination events were significantly associated with regions of the genome with lower sequence divergence between the parents and clustered in subtelomeric regions, notably in regions that had undergone introgression between the two parental species. Parallel evolution was observed, particularly through repeated homozygosity via nondisjunction, yet there was little evidence of environment-specific parallel change for either LOH, aneuploidy, or mutations. These data show that AD hybrids have both a remarkable genomic plasticity and yet are challenged in the ability to recombine through sequence divergence and chromosomal rearrangements, a scenario likely limiting the precision of adaptive evolution to novel environments. Quantifying isolation-by-resistance and connectivity in dendritic ecological networks Statistical methods that quantify individual movements have limited applicability within dendritic ecological networks (DENs), to include riverscapes. We promote an approach herein by deriving an analytical framework (=DO_SCPLOWISTC_SCPLOWNO_SCPLOWETC_SCPLOW and RO_SCPLOWESISTC_SCPLOWNO_SCPLOWETC_SCPLOW) that more appropriately quantifies differentiation by meshing individual DEN segments with fitted distances and annotated environmental features. We first explored the spatial arrangement of diversity in a river network by fitting pairwise distances against graph edges. We then derived effective resistance models (e.g., as a composite of an array of possible environmental covariates) by using a genetic algorithm for heuristic model selection, invoking circuit theory (CO_SCPLOWIRCUITC_SCPLOWSO_SCPLOWCAPEC_SCPLOW) to form a graph-based complement to RO_SCPLOWESISTANCEC_SCPLOWGA (a landscape ecology package). Although numerous pairwise distance metrics can be employed, we utilized genetic distances encompassing 13,218 ddRAD loci from N=762 Speckled Dace (Cyprinidae: Rhinichthys osculus) sampled at 78 Colorado River localities ({micro}=9.77/site). When distances for each locus were fitted into the network, we found intraspecific divergences due to stream reaches were greater than expected when distances alone. In fitting model-averaged  effective resistance models to the same network, we showed the manner by which a host of anthropogenic and environmental variables contribute to patterns of connectivity underlying riverscape genetic differentiation. Our framework represents a contemporary approach to deriving metapopulation/ metacommunity structure within DENs, in that it allows: (a) The extension of projections into unsampled temporal/spatial components; (b) Comparisons to be made among species and/or drainages, both of which benefit multi-species management; and (c) the quantification of locus-specific patterns, from which adaptive hypotheses could then be derived. O_QD"The intimate relation of hydrological conditions to fish life is proverbial; to be hopelessly or dangerously maladjusted to ones environment is to be a  fish out of water" (Hubbs 1941).

C_QD Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Protection induced by a human monoclonal antibody recognizing two different epitopes in a conserved region of streptococcal M proteins Group A streptococci have evolved multiple strategies to evade human antibodies, making it challenging to create effective vaccines or antibody treatments. Here, we have generated antibodies derived from the memory B cells of an individual who had successfully cleared a group A streptococcal infection. The antibodies bind with high affinity in the central region of the surface-bound M protein. Such antibodies are typically non-opsonic. However, one antibody could effectively promote vital immune functions, including phagocytosis and in vivo protection. Remarkably, this antibody primarily interacts through a bivalent dual-Fab cis mode, where the Fabs bind to two distinct epitopes in the M protein. The dual-Fab cis binding phenomenon is conserved across different groups of M types. In contrast, other antibodies binding with normal single-Fab mode to the same region can not bypass the M proteins virulent effects. A broadly binding, protective monoclonal antibody could be a candidate for anti-streptococcal therapy. Our findings highlight the concept of dual-Fab cis binding as a means to access conserved, and normally non-opsonic regions, for protective antibody targeting. IRGM/Irgm1 Aggravates Progression of Atherosclerosis by Inducing Macrophage Apoptosis through the MAPK Signaling Pathway AimsAtherosclerosis underlies most cardiovascular diseases, among which acute coronary syndrome (ACS) caused by plaque rupture (PR) often leads to death. Immune-related GTPases (IRGM/Irgm1) have been extensively studied in inflammatory diseases, but their role in atherosclerosis is unclear. Determining how IRGM/Irgm1 promotes atherosclerotic plaque vulnerability will provide information for new biomarkers and/or therapeutic targets.

Methods and resultsWe identified ruptured and unruptured plaques by optical coherence tomography, and found that serum IRGM was highly expressed in patients with ST-segment elevation myocardial infarction in PR. We used ApoE-/-Irgm1+/+, ApoE-/-Irgm1+/- mice and chimeric mice to establish a model of advanced atherosclerosis. The results of pathological experiments showed that Irgm1 caused plaque necrosis. The ratio of neutral lipids and cholesterol crystals increases, while the content of collagen fibers decreases, aggravating the destabilization of atherosclerotic plaques. In vitro, we used multiple approaches to confirm that Irgm1 promotes macrophage apoptosis by promoting the production of reactive oxygen species and activating the MAPK signaling pathway.

ConclusionsIRGM may be a potential risk factor for PR. Mechanistic studies have shown that IRGM/Irgm1 contributes to the formation and rupture of fragile plaques. This is partly mediated by the induction of macrophage apoptosis via the MAPK signaling pathway. IRGM may offer new strategies for early treatment of ACS.

Translation viewOur findings indicate that IRGM/Irgm1 contributes to formation and rupture of vulnerable plaques. It suggests that IRGM may provide a new target for the early treatment of ACS. Multi-omics investigation of tacrolimus nephrotoxicity BackgroundTacrolimus, an immunosuppressive drug prescribed to a majority of transplanted patients is nephrotoxic, through still unclear mechanisms. This study aims to evaluate the impact of tacrolimus on a lineage of proximal tubular cells using a multi-omics approach.

MethodsLLC-PK1 cells were exposed to 5 M of tacrolimus for 24h. Intracellular proteins and metabolites, and extracellular metabolites were extracted and analysed by LC-MS/MS. The transcriptional expression of the dysregulated proteins PCK-1, FBP1 and FBP2 was measured using RT-qPCR.

ResultsIn our cell model, tacrolimus impacted different metabolic pathways including those of arginine (e.g., citrulline, ornithine) (p < 0.0001), amino acids (e.g., valine, isoleucine, aspartic acid) (p < 0.0001) and pyrimidine (p<0.01). In addition, it induced oxidative stress (p < 0.01) as shown by a decrease in total cell glutathione quantity. It impacted cell energy through an increase in Krebs cycle intermediates (e.g., citrate, aconitate, fumarate) (p < 0.01) and down-regulation of PCK-1 (p < 0.05) and FPB1 (p < 0.01), which are key enzymes in gluconeogenesis. Apart from glucose synthesis, gluconeogenesis is an important process in kidney mediated acid-base balance control.

ConclusionThe variations found using this multi-omics approach clearly point towards a dysregulation of energy production in epithelial cells of the renal tubule, and potentially of their functions, that may be implicated in tacrolimus nephrotoxicity in the clinics. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. The emergence of highly fit SARS-CoV-2 variants accelerated by epistasis and caused by recombination The SARS-CoV-2 pandemic recently entered an alarming new phase with the emergence of the variants of concern (VOC) and understanding their biology is paramount to predicting future ones. Current efforts mainly focus on mutations in the spike glycoprotein (S), but changes in other regions of the viral proteome are likely key. We analyzed more than 900,000 SARS-CoV-2 genomes with a computational systems biology approach including a haplotype network and protein structural analyses to reveal lineage-defining mutations and their critical functional attributes. Our results indicate that increased transmission is promoted by epistasis, i.e., combinations of mutations in S and other viral proteins. Mutations in the non-S proteins involve immune-antagonism and replication performance, suggesting convergent evolution. Furthermore, adaptive mutations appear in geographically disparate locations, suggesting that either independent, repeat mutation events or recombination among different strains are generating VOC. We demonstrate that recombination is a stronger hypothesis, and may be accelerating the emergence of VOC by bringing together cooperative mutations. This emphasizes the importance of a global response to stop the COVID-19 pandemic. OrganoID: a versatile deep learning platform for tracking and analysis of single-organoid dynamics Organoids have immense potential as ex vivo disease models for drug discovery and personalized drug screening. Dynamic changes in individual organoid morphology, number, and size can indicate important drug responses, however these metrics are difficult and labor-intensive to obtain for high-throughput image datasets. Here, we present OrganoID, a robust image analysis platform that automatically recognizes, labels, and tracks single organoids, pixel-by-pixel, in brightfield and phase-contrast microscopy experiments. The platform was trained on images of pancreatic cancer organoids and validated on separate images of pancreatic, lung, colon, and adenoid cystic carcinoma organoids, which showed excellent agreement with manual measurements of organoid count (96%) and size (95%) without any parameter adjustments. Single-organoid tracking accuracy remained above 89% over a four-day time-lapse microscopy study. Automated single-organoid morphology analysis of a chemotherapy dose-response experiment identified decreased organoid circularity as an important morphological feature reflecting drug response. OrganoID enables straightforward, detailed, and accurate image analysis to accelerate the use of organoids in high-throughput, data-intensive biomedical applications. ATP dependent polymerization dynamics of bacterial actin proteins involved in Spiroplasma swimming MreB is a bacterial protein belonging to the actin superfamily. It polymerises into an antiparallel double-stranded filament that generally functions in cell shape determination by maintaining cell wall synthesis. Spiroplasma eriocheiris, a helical wall-less bacterium, has five classes of MreB homologs (SpeMreB1-5) that are likely to be involved in swimming motility. Here, we investigated the structure, ATPase activity, and polymerisation dynamics of SpeMreB3 and SpeMreB5. SpeMreB3 polymerised into an antiparallel double-stranded filament, and SpeMreB5 formed sheets, including the antiparallel filament, upon the binding of a nucleotide. SpeMreB3 showed slow Pi release owing to the lack of an amino acid motif conserved in the catalytic centre of MreB family proteins. Our crystal structures of SpeMreB3 and analyses of the mutant variants showed that the amino acid motif most likely plays a role in eliminating the proton of the nucleophilic water for ATP hydrolysis. Our sedimentation assay suggests that SpeMreB3 has a lower polymerisation activity than SpeMreB5, while their polymerisation dynamics are qualitatively similar to those of other actin superfamily proteins, in which ATP hydrolysis stabilises the filament, and Pi release leads to depolymerisation. 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. Visualization of loop extrusion by DNA nanoscale tracing in single human cells The spatial organization of the genome is essential for its functions, including gene expression, DNA replication and repair, as well as chromosome segregation1. Biomolecular condensation and loop extrusion have been proposed as the principal driving forces that underlie the formation of non-random chromatin structures such as topologically associating domains2, 3. However, whether the actual 3D-folding of DNA in single cells is consistent with these mechanisms has been difficult to address in situ. Here, we developed LoopTrace, a fluorescence imaging workflow for high-resolution reconstruction of 3D genome architecture that conserves chromatin structure at the nanoscale and can resolve the 3D-fold of chromosomal DNA with better than 5-kb precision in single human cells. Our results show that the chromatin fibre behaves as a random coil up to the megabase scale and is further structured by contacts between sites that anchor loops. Our single cell folds reveal that such looping interactions are sparse and lead to a large heterogeneity of folds with one or two dynamically positioned loop bases as the main reproducible feature of megabase scale chromosomal regions. Clustering folds by their 3D conformations revealed a series of structures consistent with progressive loop extrusion between major anchor sites. Consistently, the looping interactions and their non-random positioning depend on the presence of the loop extrusion enzyme cohesin and its anchor protein CTCF, respectively. Our approach is scalable and will be instrumental to image the functional 3D architecture of the genome directly at the nanoscale. Single amino-acid mutation in a Drosophila melanogaster ribosomal protein: an insight in uL11 transcriptional activity The ribosomal protein uL11 is located at the basis of the ribosome P-stalk and plays a paramount role in translational efficiency. In addition, no mutant for uL11 is available suggesting that this gene is haplo-insufficient as many other Ribosomal Protein Genes (RPGs). We have previously shown that overexpression of Drosophila melanogaster uL11 enhances the transcription of many RPGs and Ribosomal Biogenesis genes (RiBis) suggesting that uL11 might globally regulate the level of translation through its transcriptional activity. Moreover, uL11 trimethylated on lysine 3 (uL11K3me3) interacts with the chromodomain of the Enhancer of Polycomb and Trithorax Corto, and both proteins co- localize with RNA Polymerase II at many sites on polytene chromosomes. These data have led to the hypothesis that the N-terminal end of uL11, and more particularly the trimethylation of lysine 3, supports the extra-ribosomal activity of uL11 in transcription. To address this question, we mutated the lysine 3 codon using a CRISPR/Cas9 strategy and obtained several lysine 3 mutants. We describe here the first mutants of D. melanogaster uL11. Unexpectedly, the uL11K3A mutant, in which the lysine 3 codon is replaced by an alanine, displays a genuine Minute phenotype known to be characteristic of RPG deletions (longer development, low fertility, high lethality, thin and short bristles) whereas the uL11K3Y mutant, in which the lysine 3 codon is replaced by a tyrosine, is unaffected. In agreement, the rate of translation decreases in uL11K3A but not in uL11K3Y. Co-immunoprecipitation experiments show that the interaction between uL11 and the Corto chromodomain is impaired by both mutations. RNA- seq analysis from wing imaginal discs shows enrichment in the GO categories "glutathione metabolism" for up-regulated genes in both uL11K3A and uL11K3Y mutants and "regulation of transcription" for down-regulated genes in uL11K3A only. Analysis of the cis-regulatory sequences of these genes suggests that uL11 might regulate transcription of target genes in concert with the couple of transcription factors Mad/Med that mediate response to the Bone Morphogenetic Protein (BMP) signaling pathway. From fitting the average to fitting the individual: A cautionary tale for mathematicalmodelers An outstanding challenge in the clinical care of cancer is moving from a one-size-fits-all approach that relies on population-level statistics towards personalized therapeutic design. Mathematical modeling is a powerful tool in treatment personalization, as it allows for the incorporation of patient-specific data so that treatment can be tailor-designed to the individual. Herein, we work with a mathematical model of murine cancer immunotherapy that has been previously-validated against the average of an experimental dataset. We ask the question: what happens if we try to use this same model to perform personalized fits, and therefore make individualized treatment recommendations? Typically, this would be done by choosing a single fitting methodology, and a single cost function, identifying the individualized best-fit parameters, and extrapolating from there to make personalized treatment recommendations. Our analyses show the potentially problematic nature of this approach, as predicted personalized treatment response proved to be sensitive to the fitting methodology utilized. We also demonstrate how a small amount of the right additional experimental measurements could go a long way to improve consistency in personalized fits. Finally, we show how quantifying the robustness of the average response could also help improve confidence in personalized treatment recommendations.

Author summaryAs we enter the era of healthcare where personalized medicine becomes a more common approach to treating cancer patients, harnessing the power of mathematical models will only become more essential. Using a preclinical dataset on cancer immunotherapy, we explore the challenges and limitations that arise when trying to move from fitting and making predictions for the population-level average, to fitting and making predictions for an individual. We find that the standard of approach of picking a single fitting methodology and a single cost function may end up having limited predictive value when applied to individual data. We also show how having a small amount of the right additional experimental data, and establishing the robustness of average treatment response, can help improve confidence in personalized model predictions. First evidence of in vitro cytotoxic effects of marine microlitter on Merluccius merluccius and Mullus barbatus, two Mediterranean commercial fish species Marine litter is composed mainly of plastics and is recognized as a serious threats to marine ecosystems. Ecotoxicological approaches have started elucidating the potential severity of microplastics (MPs) in controlled laboratory studies with pristine materials but no information exist on marine environmental microlitter as a whole. Here, we characterized the litter in the coastal Northern Tyrrhenian sea and in the stomach of two fish species of socio-economic importance, and exposed primary cell cultures of mucosal and lymphoid organs to marine microlitter for evaluating possible cytotoxic effects. An average of 0.30 {+/-} 0.02 microlitter items m-3 was found in water samples. FT-IR analysis revealed that plastic particles, namely HDPE, polyamide and polypropylene were present in 100% and 83.3% of Merluccius merluccius and Mullus barbatus analyzed, which overall ingested 14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual, respectively. Moreover, microlitter was confirmed as a vector of microorganisms. Lastly, the apical end-point of viability was found to be significantly reduced in splenic cells exposed in vitro to two microlitter conditions. Considering the role of the spleen in the mounting of adaptive immune responses, our results warrant more in-depth investigations for clarifying the actual susceptibility of these two species to anthropogenic microlitter.

HighlightsO_LI0.30 {+/-} 0.02 microlitter items m-3 were found at the surface of coastal Northern Tyrrhenian sea
C_LIO_LI14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual were retrieved from the stomach of hakes and mullets
C_LIO_LIThe ingested microlitter contained plastic items
C_LIO_LIMicrolitter was validated as a carrier of bacteria, fungi and flagellates
C_LIO_LISplenic cells exposed to two microlitter conditions for 72 hours suffered cytotoxicity
C_LI Structural and molecular determinants for the interaction of ExbB from Serratia marcescens and HasB, a TonB paralog ExbB and ExbD are cytoplasmic membrane proteins that associate with TonB to convey the energy of the proton-motive force to outer membrane receptors in Gram-negative bacteria for iron uptake. The opportunistic pathogen Serratia marcescens (Sm) possesses both TonB and a heme-specific TonB paralog, HasB. ExbBSm has a long periplasmic extension absent in other bacteria such as E. coli (Ec). Long ExbBs are found in several genera of Alphaproteobacteria, most often in correlation with a hasB gene. We investigated specificity determinants of ExbBSm and HasB. We determined the cryo-EM structures of ExbBSm and of the ExbB-ExbDSm complex from S. marcescens. ExbBSm alone is a stable pentamer, and its complex includes two ExbD monomers. We showed that ExbBSm extension interacts with HasB and is involved in heme acquisition and we identified key residues in the membrane domain of ExbBSm and ExbBEc, essential for function and likely involved in the interaction with TonB/HasB. Our results shed light on the new class of inner membrane energy machinery formed by ExbB,ExbD and HasB.