Benchmarking of deep learning algorithms for 3D instance segmentation of confocal image datasets Segmenting three dimensional microscopy images is essential for understanding phenomena like morphogenesis, cell division, cellular growth and genetic expression patterns. Recently, deep learning (DL) pipelines have been developed which claim to provide high accuracy segmentation of cellular images and are increasingly considered as the state-of-the-art for image segmentation problems. However, it remains difficult to define their relative performances as the concurrent diversity and lack of uniform evaluation strategies makes it difficult to know how their results compare. In this paper, we first made an inventory of the available DL methods for 3 dimensional (3D) cell segmentation. We next implemented and quantitatively compared a number of representative DL pipelines, alongside a highly efficient non-DL method named MARS. The DL methods were trained on a common dataset of 3D cellular confocal microscopy images. Their segmentation accuracies were also tested in the presence of different image artifacts. A specific method for segmentation quality evaluation was adopted which isolates segmentation errors due to under/over segmentation. This is complemented with a 3D visualization strategy for interactive exploration of segmentation quality. Our analysis shows that the DL pipelines have different levels of accuracy. Two of them, which are end to end 3D and were originally designed for cell boundary detection, show high performance, and offer clear advantages in terms of adaptability to new data.

Author summaryIn recent years a number of deep learning (DL) algorithms based on computational neural networks have been developed which claim to achieve high accuracy and automatic segmentation of 3D microscopy images. Although these algorithms have received considerable attention in the literature, it is difficult to evaluate their relative performances, while it remains unclear whether they really perform better than other, more classical segmentation methods.

To clarify these issues, we performed a detailed, quantitative analysis of a number of representative DL pipelines for cell instance segmentation from 3D confocal microscopy image datasets. We developed a protocol for benchmarking the performances of such DL based segmentation pipelines using common training and test datasets, evaluation metrics and visualizations. Using this protocol, we evaluated and compared four different DL pipelines to identify their strengths and limitations. A high performance non-DL method was also included in the evaluation. We show that DL pipelines may show significant differences in their performances depending on their model architecture and pipeline components but overall show excellent adaptability to unseen data. We also show that our benchmarking protocol could be extended to a variety of segmentation pipelines and datasets.
The immunopeptidome from a genomic perspective: Establishing immune-relevant regions for cancer vaccine design A longstanding disconnect between the growing number of MHC Class I immunopeptidomic studies and genomic medicine hinders cancer vaccine design. We develop COD-dipp to genomically map the full spectrum of detected canonical and non-canonical (non-exonic) MHC Class I antigens from 26 cancer studies. We demonstrate that patient mutations in regions overlapping physically identified antigens better predict immunotherapy response when compared to neoantigen predictions. We suggest a vaccine design approach using 140,966 highly immune-visible regions of the genome annotated by their expression and haplotype frequency in the human population. These regions tend to be highly conserved, mutated in cancer and harbor 7.8 times more immunogenicity. Intersecting pan-cancer mutations with these immune surveilled regions revealed a potential to create off-the-shelf multi-epitope vaccines against public neoantigens. Here we release COD-dipp, a cancer vaccine toolkit as a web-application (https://www.proteogenomics.ca/COD-dipp) and open-source high-throughput resource.
An integrative approach leads to the discovery of a novel anti-leukemic peptide from human milk Chemotherapy in childhood leukemia is associated with late morbidity in leukemic survivors, while certain patient subsets are relatively resistant to standard chemotherapy. It is therefore important to identify new agents with sensitivity and selectivity towards leukemic cells, while having less systemic toxicity. Peptide-based therapeutics has gained much attention during the last few years. Here, we used an integrative workflow combining mass spectrometric peptide library construction, in silico anticancer peptide screening, and in vitro leukemic cell studies to discover a novel anti-leukemic peptide having 3+charges and alpha-helical structure, namely HMP-S7, from human breast milk. HMP-S7 showed cytotoxic activity against four distinct leukemic cell lines in a dose-dependent manner but had no effect on solid malignancies or representative normal cells. HMP-S7 induced leukemic cell death by penetrating the plasma membrane to enter the cytoplasm and cause leakage of lactate dehydrogenase, thus acting in a membranolytic manner. Importantly, HMP-S7 exhibited anti-leukemic effect against patient-derived leukemic cells ex vivo. In conclusion, HMP-S7 is a selective anti-leukemic peptide with promise which requires further validation in preclinical and clinical studies.

TeaserIn silico screening of naturally occurring human milk peptides discovers a new anticancer peptide that kills leukemic cells in vitro and ex vivo. Bacterial H-NS contacts DNA at the same irregularly spaced sites in both bridged and hemi-sequestered linear filaments Gene silencing in bacteria is mediated by chromatin proteins, of which Escherichia coli H-NS is a paradigmatic example. H-NS forms nucleoprotein filaments with either one or two DNA duplexes. However, the structures, arrangements of DNA-binding domains (DBDs), and positions of DBD-DNA contacts in linear and bridged filaments are uncertain. To characterize the contacts that silence transcription by RNA polymerase, we combined {middle dot}OH footprinting, molecular dynamics, statistical modeling, and DBD mapping using a chemical nuclease (Fe2+-EDTA) tethered to the DBDs (TEN-map). We find that H-NS DBDs contact DNA at indistinguishable locations in bridged or linear filaments and that the DBDs vary in orientation and position with ~10-bp average spacing. Our results support a hemi-sequestration model of linear-to-bridged H-NS switching in which linear filaments able to inhibit only transcription initiation switch to bridged filaments able to inhibit both initiation and elongation using the same irregularly spaced DNA contact sites.

HighlightsO_LITethered-nuclease mapping (TEN-map) of H-NS DNA-binding domains detects DNA contacts
C_LIO_LIBridged and linear H-NS filaments use the same DNA contact sites
C_LIO_LIH-NS-DNA contacts are unevenly spaced with ~10 bp average separation
C_LIO_LIAT-steps, minor groove width, and electrostatic potential best predict contact sites
C_LI A network approach to genetic circuit designs As genetic circuits become more sophisticated, the size and complexity of data about their designs increases. This data captured goes beyond monolithic genetic sequences and towards circuit modularity and functional details, which are beneficial for analyzing circuit performance and establishing design automation techniques. However, the accessibility, visualisation and usability of design data (and metadata) have received relatively little attention to date. Here, we present a method to turn circuit designs into networks and showcase its potential to enhance the utility of design data. Since networks are dynamic structures, initial graphs can be interactively shaped into sub-networks of relevant information based on requirements such as abstraction, hierarchy and protein interactions. Additionally, several visual changes can be applied, such as colouring or clustering nodes based on types (e.g., genes or promoters), resulting in easier comprehension from a user perspective. This approach allows circuit designs to be coupled to other networks, such as metabolic pathways or implementation protocols captured in graph-like formats. Therefore, we advocate using networks to structure, access and improve synthetic biology information. Electrocommunication signals indicate motivation to compete during dyadic interactions of an electric fish Animals across species compete for limited resources. While in some species competition behavior is solely based on own abilities, others assess their opponents to facilitate these interactions. Using cues and communication signals, contestants gather information about their opponent, adjust their behavior accordingly, and can thereby avoid high costs of escalating fights. We tracked electrocommunication signals, in particular "rises", and agonistic behaviors of the gymnotiform electric fish Apteronotus leptorhynchus in staged competition experiments. A larger body size relative to the opponent was the sole significant predictor for winners. Sex and the frequency of the continuously emitted electric field were only mildly influencing competition outcome. In males, correlations of body-size and winning were stronger than in females and, especially when losing against females, communication and agonistic interactions were enhanced, hinting towards males being more motivated to compete. Fish that lost competitions emitted the majority of rises in dependence on the competitors relative size and sex. The emission of rises was costly since it provoked ritualized biting or chasing behaviors by the other fish. Despite winners being accurately predictable based on rise numbers already after the initial 25 minutes, losers continued to emit rises. The number of rises emitted by losers and the duration of chasing behaviors depended in similar ways on physical attributes of contestants. The detailed evaluation of these correlations hint towards A. leptorhynchus adjusting their competition behavior according to mutual assessment, where rises could signal a losers motivation to continue assessment through ritualized fighting. An adaptive teosinte mexicana introgression modulates phosphatidylcholine levels and is associated with maize flowering time Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp.parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability, and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern US, Canada and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time. Determinants of Murine Skin Microbiota Composition in Homeostasis and Wound Healing Animal microbiota have complex interactions with hosts and environment that determines its composition. Yet the ability of hosts to determine their microbiota composition is less well studied. In this study, to investigate the role host genetics in determining skin microbiota, we used 30 different mouse strains from the recombinant inbred panel, the Collaborative Cross. Murine skin microbiota composition was strongly dependent on murine strain with > 50% of the variation explained by murine strain. In particular, a quantitative trait locus on chromosome 4 associates both with Staphylococcus abundance and principal-component multi-trait analyses. Additionally, excisional wound associated changes in microbiota composition were not uniform across mouse strains and were host-specific, the genetic background accounting for about 40% of the variation in microbiota. Genetic background also had the highest effect on the healing speed of wounds accounting for over 50% of the variation while mouse age and microbiota composition change accounted only for 20% and 5% of the healing speed despite reaching statistical significance. In conclusion, host genetics has a significant impact on the skin microbiota composition during both homeostasis and wound healing. These findings have long reaching implications in our understanding of associations between microbiota dysbiosis and disease. Therapeutic resistance in acute myeloid leukemia cells is mediated by a novel ATM/mTOR pathway regulating oxidative phosphorylation While leukemic cells are susceptible to various therapeutic insults, residence in the bone marrow microenvironment typically confers protection from a wide range of drugs. Thus, understanding the unique molecular changes elicited by the marrow is of critical importance towards improving therapeutic outcomes. In the present study, we demonstrate that aberrant activation of oxidative phosphorylation serves to induce therapeutic resistance in FLT3 mutant human AML cells challenged with FLT3 inhibitor drugs. Importantly, our findings show that AML cells are protected from apoptosis following FLT3 inhibition due to marrow-mediated activation of ATM, which in turn up-regulates oxidative phosphorylation via mTOR signaling. mTOR is required for the bone marrow stroma-dependent maintenance of protein translation, with selective polysome enrichment of oxidative phosphorylation transcripts, despite FLT3 inhibition. To investigate the therapeutic significance of this finding, we tested the mTOR inhibitor everolimus in combination with the FLT3 inhibitor quizartinib in primary human AML xenograft models. While marrow resident AML cells were highly resistant to quizartinib alone, the addition of everolimus induced profound reduction in tumor burden and prevented relapse. Taken together, these data provide a novel mechanistic understanding of marrow-based therapeutic resistance, and a promising strategy for improved treatment of FLT3 mutant AML patients. Dynamics of SARS-CoV-2 mutations reveals regional-specificity and similar trends of N501 and High-Frequency mutation N501Y in different levels of control measures Coronavirus disease 2019 (COVID-19) is a contagious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This disease has spread globally, causing more than 161.5 million cases and 3.3 million deaths to date. Surveillance and monitoring of new mutations in the virus genome are crucial to our understanding of the adaptation of SARS-CoV-2. Moreover, how the temporal dynamics of these mutations is influenced by control measures and non-pharmaceutical interventions (NPIs) is poorly understood. Using 1 058 020 SARS-CoV-2 from sequenced COVID-19 cases from 98 countries (totaling 714 country-month combinations), we perform a normalization by COVID-19 cases to calculate the relative frequency of SARS-CoV-2 mutations and explore their dynamics over time. We found 115 mutations estimated to be present in more than 3 % of global COVID-19 cases and determined three types of mutation dynamics: High-Frequency, Medium-Frequency, and Low-Frequency. Classification of mutations based on temporal dynamics enable us to examine viral adaptation and evaluate the effects of implemented control measures in virus evolution during the pandemic. We showed that Medium-Frequency mutations are characterized by high prevalence in specific regions and/or in constant competition with other mutations in several regions. Finally, taking N501Y mutation as representative of High-Frequency mutations, we showed that level of control measure stringency negatively correlates with the effective reproduction number of SARS-CoV-2 with High-Frequency or not-High-Frequency and both follows similar trends in different levels of stringency. Acetyltransferases GCN5 and PCAF are required for B lymphocyte maturation in mice B lymphocyte development has two DNA recombination processes: V(D)J recombination of the immunoglobulin (Igh) gene variable region and class switching of the Igh constant regions from IgM to IgG, IgA, or IgE. V(D)J recombination is required for successful maturation of B cells from pro-B to pre-B to immature-B and then to mature B cells in the bone marrow. CSR occurs outside of the bone marrow when mature B cells migrate to peripheral lymphoid organs, such as spleen and lymph nodes. Both V(D)J recombination and CSR depend on an open chromatin state that makes DNA accessible to specific enzymes, recombination activating gene (RAG), and activation-induced cytidine deaminase (AID). Acetyltransferases GCN5 and PCAF possess redundant functions acetylating histone H3 lysine 9 (H3K9). Here, we generated a mouse model that lacks both GCN5 and PCAF in B cells. We found that double-deficient mice possess low levels of mature B cells in the bone marrow and peripheral organs, an accumulation of pro-B cells in bone marrow, and reduced CSR levels. We conclude that both GCN5 and PCAF are required for B cell development in vivo. Morphological aspects of opossum (Didelphis aurita) pancreas during intramarsupial development Functionalized nanostructured systems can be used for imaging and drug delivery for anti-tumor therapy, including breast tumors. This is a more efficient approach that offers reduced systemic side effects compared to conventional diagnostic and chemotherapy methods. Multifunctional nanoparticles are potential tools in the diagnosis, location tracing and kill tumor cells through a less invasive manner. Functionalized phosphate-based nanoparticles are capable of encapsulating, or may be associated, with fluorescent probes. In this study, we synthesize a nanoparticle phosphate-based composite (NPC) and functionalize it with poly-ethylene glycol (PEG), hyaluronic acid (HA), the fluorescent probe rhodamin 6G (R6G) and the antimitotic doxorubicin (DOX). We focused on targeting human breast cancer cells reporting the biological effects of functionalized NPC on them. NPC and NPC formulations containing PEG, HA, and R6G did not cause cell viability reduction on MCF-7 and MDA-MB-231 cell lines. The cellular internalization of NPC was quantified by real-time in vitro observation, and confirmed by electron microscopy techniques. Intracellular NPC distribution is detected in the cytoplasm and nucleus of tumor cells by confocal fluorescent images. The percent association of doxorubicin to NPC matrix was approximately 18% and NPC formulations associated with doxorubicin led to a significant reduction in cell viability in MDA-MB-231 and MCF-7 cells. This data suggest the potential use of NPC as a non-cytotoxic platform for association with functional ligands to selective targeting breast cancer cells. NPC use can be also explored in drug delivery to cancer cells. The Skeletal Muscle Circadian Clock Regulates Titin Splicing Through RBM20 Circadian rhythms are maintained by a cell autonomous, transcriptional-translational feedback loop known as the molecular clock. While previous research suggests a role of the molecular clock in regulating skeletal muscle structure and function, no mechanisms have connected the molecular clock to sarcomere filaments. Utilizing inducible, skeletal muscle specific, Bmal1 knockout (iMSBmal1-/-) mice, we showed that knocking out skeletal muscle clock function alters titin isoform expression through both RNAseq and LC-MS. This alteration in titins spring length resulted in sarcomere length heterogeneity. We demonstrate the direct link between altered titin splicing and sarcomere length in vitro using U7 snRNPs that truncate the region of titin altered in iMSBmal1-/- muscle. Finally, we identified a mechanism whereby the skeletal muscle clock regulates titin isoform expression through transcriptional regulation of Rbm20, a potent splicing regulator of titin. Our findings demonstrate the importance of the skeletal muscle circadian clock in maintaining titin isoform through regulation of RBM20 expression. Because circadian rhythm disruption is a feature of many chronic diseases, our results highlight a novel pathway that could be targeted to maintain skeletal muscle structure and function in a range of pathologies. Genetic architecture of natural variation of cardiac performance in flies Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP)1. Genome Wide Associations Studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict Transcription Factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. A connectivity-constrained computational account of topographic organization in primate high-level visual cortex Inferotemporal cortex (IT) in humans and other primates is topo-graphically organized, containing multiple hierarchically-organized areas selective for particular domains, such as faces and scenes. This organization is commonly viewed in terms of evolved domain-specific visual mechanisms. Here, we develop an alternative, domain-general and developmental account of IT cortical organization. The account is instantiated as an Interactive Topographic Network (ITN), a form of computational model in which a hierarchy of model IT areas, subject to connectivity-based constraints, learns high-level visual representations optimized for multiple domains. We find that minimizing a wiring cost on spatially organized feedforward and lateral connections within IT, combined with constraining the feedforward processing to be strictly excitatory, results in a hierarchical, topographic organization. This organization replicates a number of key properties of primate IT cortex, including the presence of domain-selective spatial clusters preferentially involved in the representation of faces, objects, and scenes, columnar responses across separate excitatory and inhibitory units, and generic spatial organization whereby the response correlation of pairs of units falls off with their distance. We thus argue that domain-selectivity is an emergent property of a visual system optimized to maximize behavioral performance while minimizing wiring costs.

Significance StatementWe introduce the Interactive Topographic Network, a framework for modeling high-level vision, to demonstrate in computational simulations that the spatial clustering of domains in late stages of the primate visual system may arise from the demands of visual recognition under the constraints of minimal wiring costs and excitatory between-area neuronal communication. The learned organization of the model is highly specialized but not fully modular, capturing many of the properties of organization in primates. Our work is significant for cognitive neuroscience, by providing a domain-general developmental account of topo-graphic functional specialization, and for computational neuroscience, by demonstrating how well-known biological details can be successfully incorporated into neural network models in order to account for critical empirical findings. Fatal COVID-19 ARDS associated with incomplete AEC1 differentiation from the transitional state without senescence or fibrosis ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. Differences in gait stability and acceleration characteristics between healthy young and older females Our aim was to evaluate differences in gait acceleration intensity, variability and stability of feet and trunk between older females and young females using inertial sensors. Twenty older females (OF; mean age 68.4, SD 4.1 years) and eighteen young females (YF; mean age 22.3, SD 1.7 years) were asked to walk straight for 100 meters at their preferred speed, while wearing inertial sensors on heels and lower back. We calculated spatiotemporal measures, foot and trunk acceleration characteristics and their variability, as well as trunk stability using the local divergence exponent (LDE). Two-way analysis of variance (including the factors foot and age), Students t-test, and Mann-Whitney U test were used to compare statistical differences of measures between groups. Cohens d effects were calculated for each variable. Foot maximum vertical acceleration and amplitude, trunk-foot vertical acceleration attenuation, as well as their variability were significantly smaller in OF than in YF. In contrast, trunk mediolateral acceleration amplitude, maximum vertical acceleration, and amplitude, as well as their variability were significantly larger in OF than in YF. Moreover, OF showed lower stability (i.e., higher LDE values) in mediolateral acceleration, mediolateral and vertical angular velocity of the trunk. Even though we measured healthy older females, these participants showed lower vertical foot accelerations with higher vertical trunk acceleration, lower trunk-foot vertical acceleration attenuation, less gait stability, and more variability of the trunk, and hence, were more likely to fall. These findings suggest that instrumented gait measurements may help for early detection of changes or impairments in gait performance, even before this can be observed by clinical eye or gait speed. Evolutionary dynamics of sex-biased genes expressed in cricket brains and gonads Sex-biased gene expression, particularly sex-biased expression in the gonad, has been linked to rates of protein sequence evolution (nonsynonymous to synonymous substitutions, dN/dS) in animals. However, in insects, sex-biased expression studies remain centered on a few holometabolous species, and moreover, other major tissue types such as the brain remain underexplored. Here, we studied sex-biased gene expression and protein evolution in a hemimetabolous insect, the cricket Gryllus bimaculatus. We generated novel male and female RNA-seq data for two sexual tissue types, the gonad and somatic reproductive system, and for two core components of the nervous system, the brain and ventral nerve cord. From a genome-wide analysis, we report several core findings. Firstly, testis-biased genes had accelerated evolution, as compared to ovary-biased and unbiased genes, which was associated with positive selection events. Secondly, while sex-biased brain genes were much less common than for the gonad, they exhibited a striking tendency for rapid protein evolution, an effect that was stronger for the female than male brain. Further, some sex-biased brain genes were linked to sexual functions and mating behaviors, which we suggest may have accelerated their evolution via sexual selection. Thirdly, a tendency for narrow cross-tissue expression breadth, suggesting low pleiotropy, was observed for sex-biased brain genes, suggesting relaxed purifying selection, which we speculate may allow enhanced freedom to evolve adaptive protein functional changes. The findings of rapid evolution of testis-biased genes and male and female-biased brain genes are discussed with respect to pleiotropy, positive selection, and the mating biology of this cricket. Thermodynamic constraints on the assembly and diversity ofmicrobial ecosystems are different near to and far fromequilibrium 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. Probing effects of the SARS-CoV-2 E protein on membrane curvature and intracellular calcium SARS-CoV-2 contains four structural proteins in its genome. These proteins aid in the assembly and budding of new virions at the ER-Golgi intermediate compartment (ERGIC). Current fundamental research efforts largely focus on one of these proteins - the spike (S) protein. Since successful antiviral therapies are likely to target multiple viral components, there is considerable interest in understanding the biophysical role of its other structural proteins, in particular structural membrane proteins. Here, we have focused our efforts on the characterization of the full-length envelope (E) protein from SARS-CoV-2, combining experimental and computational approaches. Recombinant expression of the full-length E protein from SARS-CoV-2 reveals that this membrane protein is capable of independent multimerization, possibly as a tetrameric or smaller species. Fluorescence microscopy shows that the protein localizes intracellularly, and coarse-grained MD simulations indicate it causes bending of the surrounding lipid bilayer, corroborating a potential role for the E protein in viral budding. Although we did not find robust electrophysiological evidence of ion-channel activity, cells transfected with the E protein exhibited reduced intracellular Ca2+, which may further promote viral replication. However, our atomistic MD simulations revealed that previous NMR structures are relatively unstable, and result in models incapable of ion conduction. Our study highlights the importance of using high-resolution structural data obtained from a full-length protein to gain detailed molecular insights, and eventually permitting virtual drug screening. Fine-scale adaptations to environmental variation and growth strategies drive phyllosphere Methylobacterium diversity. Methylobacterium is a prevalent bacterial genus of the phyllosphere. Despite its ubiquity, little is known about the extent to which its diversity reflects neutral processes like migration and drift, versus environmental filtering of life history strategies and adaptations. In two temperate forests, we investigated how phylogenetic diversity within Methylobacterium was structured by biogeography, seasonality, and growth strategies. Using deep, culture-independent barcoded marker gene sequencing coupled with culture-based approaches, we uncovered a considerable diversity of Methylobacterium in the phyllosphere. We cultured different subsets of Methylobacterium lineages depending upon the temperature of isolation and growth (20 {degrees}C or 30 {degrees}C), suggesting long-term adaptation to temperature. To a lesser extent than temperature adaptation, Methylobacterium diversity was also structured across large (>100km; between forests) and small geographical scales (<1.2km within forests), among host tree species, and was dynamic over seasons. By measuring growth of 79 isolates at different temperature treatments, we observed contrasting growth performances, with strong lineage- and season-dependent variations in growth strategies. Finally, we documented a progressive replacement of lineages with a high-yield growth strategy typical of cooperative, structured communities, in favor of those characterized by rapid growth, resulting in convergence and homogenization of community structure at the end of the growing season. Together our results show how Methylobacterium is phylogenetically structured into lineages with distinct growth strategies, which helps explain their differential abundance across regions, host tree species, and time. This works paves the way for further investigation of adaptive strategies and traits within a ubiquitous phyllosphere genus.

ImportanceMethylobacterium is a bacterial group tied to plants. Despite its ubiquity and importance to their hosts, little is known about the processes driving Methylobacterium community dynamics. By combining traditional culture-dependent and independent (metabarcoding) approaches, we monitored Methylobacterium diversity in two temperate forests over a growing season. On the surface of tree leaves, we discovered remarkably diverse and dynamic Methylobacterium communities over short temporal (from June to October) and spatial scales (within 1.2 km). Because we cultured different subsets of Methylobacterium diversity depending on the temperature of incubation, we suspected that these dynamics partly reflected climatic adaptation. By culturing strains in lab conditions mimicking seasonal variations, we found that diversity and environmental variations were indeed good predictors of Methylobacterium growth performances. Our findings suggest that Methylobacterium community dynamics at the surface of tree leaves results from the succession of strains with contrasted growth strategies in response to environmental variations. 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. Single-Cell Multiomics Reveals Distinct Cell States at the Top of the Human Hematopoietic Hierarchy The advent of single cell (Sc) genomics has challenged the dogma of haematopoiesis as a tree-like structure of stepwise lineage commitment through distinct and increasingly restricted progenitor populations. Instead, analysis of ScRNA-seq has proposed that the earliest events in human hematopoietic stem cell (HSC) differentiation are characterized by only subtle molecular changes, with hematopoietic stem and progenitor cells (HSPCs) existing as a continuum of low-primed cell-states that gradually transition into a specific lineage (CLOUD-HSPCs). Here, we combine ScRNA-seq, ScATAC-seq and cell surface proteomics to dissect the heterogeneity of CLOUD-HSPCs at different stages of human life. Within CLOUD-HSPCs, pseudotime ordering of both mRNA and chromatin data revealed a bifurcation of megakaryocyte/erythroid and lympho/myeloid trajectories immediately downstream a subpopulation with an HSC-specific enhancer signature. Importantly, both HSCs and lineage-restricted progenitor populations could be prospectively isolated based on correlation of their molecular signatures with CD35 and CD11A expression, respectively. Moreover, we describe the changes that occur in this heterogeneity as hematopoiesis develops from neonatal to aged bone marrow, including an increase of HSCs and depletion of lympho-myeloid biased MPPs. Thus, this study dissects the heterogeneity of human CLOUD-HSPCs revealing distinct HSPC-states of relevance in homeostatic settings such as ageing. Novel MOG analogues to explore the MCT2 pharmacophore, α-ketoglutarate biology and cellular effects of N-oxalylglycine -ketoglutarate (KG) is a central metabolic node with a broad influence on cellular physiology. The KG analogue N-oxalylglycine (NOG) and its membrane-permeable pro-drug derivative dimethyl-oxalylglycine (DMOG) have been extensively used as tools to study prolyl hydroxylases (PHDs) and other KG-dependent processes. In cell culture media, DMOG is rapidly converted to MOG, which enters cells through monocarboxylate transporter MCT2, leading to intracellular NOG concentrations that are sufficiently high to inhibit glutaminolysis enzymes and cause cytotoxicity. Therefore, the degree of (D)MOG instability together with MCT2 expression levels determine the intracellular targets NOG engages with and, ultimately, its effects on cell viability. Here we designed and characterised a series of MOG analogues with the aims of improving stability and exploring the functional requirements for interaction with MCT2, a relatively understudied member of the SLC16 family. We report MOG analogues that maintain ability to enter cells via MCT2, and identify compounds that do not inhibit glutaminolysis or cause cytotoxicity but can still inhibit PHDs. We use these analogues to show that glutaminolysis-induced activation of mTORC1 can be uncoupled from PHD activity. Therefore, these new compounds can help deconvolute cellular effects that result from the polypharmacological action of NOG. First detection of herpesvirus and prevalence of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity.