Patterns of and processes shaping population structure and introgression among recently differentiated Drosophila melanogaster populations Despite a century of genetic analysis, the evolutionary history underlying patterns of exceptional genetic and phenotypic variation in the model organism Drosophila melanogaster remains poorly understood. How genetic and phenotypic variation is partitioned across the range of D. melanogaster, particularly in its putative ancestral range in Subtropical Africa, remains unresolved. Here, we assess patterns of population genetic structure, admixture, mate preference, and genetic incompatibility across a global sample, including 174 new accessions from remote regions within Subtropical Africa. While almost all Out of Africa genomes correspond to a single genetic ancestry, different geographic regions within Africa contain multiple ancestries, with substantial cryptic diversity in Subtropical Africa. Admixture between distinct lineages is prevalent across the range, but admixture rates vary between lineages. Female mate choice within Subtropical Africa is highly polymorphic and behavioral types are not monophyletic. The genetic architecture of mate choice is highly polygenic, including loci associated with neurological development, behavior, olfactory perception, and learning. Finally, we discovered that many segregating putative incompatibilities likely evolved during or after expansion out of Africa. This work contributes to our understanding of the evolutionary history of a key model system, and provides insight into the distribution of polymorphic reproductive barriers.
Constraints on the evolution of toxin-resistant Na,K-ATPases have limited dependence on sequence divergence A growing body of theoretical and experimental evidence suggests that intramolecular epistasis is a major determinant of rates and patterns of protein evolution and imposes a substantial constraint on the evolution of novel protein functions. Here, we examine the role of intramolecular epistasis in the recurrent evolution of resistance to cardiotonic steroids (CTS) across tetrapods, which occurs via specific amino acid substitutions to the -subunit family of Na,K-ATPases (ATP1A). After identifying a series of recurrent substitutions at two key sites of ATP1A that are predicted to confer CTS resistance in diverse tetrapods, we then performed protein engineering experiments to test the functional consequences of introducing these substitutions onto divergent species backgrounds. In line with previous results, we find that substitutions at these sites can have substantial background-dependent effects on CTS resistance. Globally, however, these substitutions also have pleiotropic effects that are consistent with additive rather than background-dependent effects. Moreover, the magnitude of a substitutions effect on activity does not depend on the overall extent of ATP1A sequence divergence between species. Our results suggest that epistatic constraints on the evolution of CTS-resistant forms of Na,K-ATPase likely depend on a small number of sites, with little dependence on overall levels of protein divergence. We propose that dependence on a limited number sites may account for the observation of convergent CTS resistance substitutions observed among taxa with highly divergent Na,K-ATPases.

Significance StatementIndividual amino acid residues within a protein work in concert to produce a functionally coherent structure that must be maintained even as orthologous proteins in different species diverge over time. Given this dependence, we expect identical mutations to have more similar effects on protein function in more closely related species. We tested this hypothesis by performing protein-engineering experiments on ATP1A, an enzyme mediating target-site insensitivity to cardiotonic steroids (CTS) in diverse animals. These experiments reveal that the phenotypic effects of substitutions can sometimes be background-dependent, but also that the magnitude of these phenotypic effects does not correlate with overall levels of ATP1A sequence divergence. Our results suggest that epistatic constraints are determined by states at a small number of sites, potentially explaining the frequent convergent CTS resistance substitutions among Na,K-ATPases of highly divergent taxa.
A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Genetic interactions drive heterogeneity in causal variant effect sizes for gene expression and complex traits Despite the growing number of genome-wide association studies (GWAS), it remains unclear to what extent gene-by-gene and gene-by-environment interactions influence complex traits in humans. The magnitude of genetic interactions in complex traits has been difficult to quantify because GWAS are generally underpowered to detect individual interactions of small effect. Here, we develop a method to test for genetic interactions that aggregates information across all trait-associated loci. Specifically, we test whether SNPs in regions of European ancestry shared between European American and admixed African American individuals have the same causal effect sizes. We hypothesize that in African Americans, the presence of genetic interactions will drive the causal effect sizes of SNPs in regions of European ancestry to be more similar to those of SNPs in regions of African ancestry. We apply our method to two traits: gene expression in 296 African Americans and 482 European Americans in the Multi-Ethnic Study of Atherosclerosis (MESA) and low-density lipoprotein cholesterol (LDL-C) in 74K African Americans and 296K European Americans in the Million Veteran Program (MVP). We find significant evidence for genetic interactions in our analysis of gene expression; for LDL-C, we observe a similar point estimate although this is not significant, likely due to lower statistical power. These results suggest that gene-by-gene or gene-by-environment interactions modify the effect sizes of causal variants in human complex traits. ACE2-independent SARS-CoV-2 infection and mouse adaption emerge after passage in cells expressing human and mouse ACE2 Human ACE2 (hACE2) is the key cell attachment and entry receptor for SARS-CoV-2, with the original SARS-CoV-2 isolates unable to use mouse ACE2 (mACE2). Herein we describe a new system for generating mouse-adapted SARS-CoV-2 in vitro by serial passaging virus in co-cultures of cell lines expressing hACE2 and mACE2. Mouse-adapted viruses emerged with up to five amino acid changes in the spike protein, all of which have been seen in human isolates. Mouse-adapted viruses replicated to high titers in C57BL/6J mouse lungs and nasal turbinates, and caused severe lung histopathology. One mouse-adapted virus was also able to replicate efficiently in ACE2-negative cell lines, with ACE2-independent entry by SARS-CoV-2 representing a new biology for SARS-CoV-2 that has potential widespread implications for disease and intervention development. The conserved endocannabinoid anandamide modulates olfactory sensitivity to induce hedonic feeding in C. elegans The ability of cannabis to increase consumption of food has been known for centuries. In addition to producing hyperphagia, cannabinoids can amplify existing preferences for calorically dense, palatable food sources, a phenomenon called hedonic feeding. These effects result from the action of plant-derived cannabinoids on brain receptors where they mimic natural ligands called endocannabinoids. The high degree of conservation of cannabinoid signaling at the molecular level across the animal kingdom suggests hedonic feeding may also be widely conserved. Here we show that exposure of C. elegans to anandamide, an endocannabinoid common to nematodes and mammals, shifts both appetitive and consummatory responses toward nutritionally superior food, an effect analogous to hedonic feeding. We find that anandamides effect on feeding requires the C. elegans cannabinoid receptor NPR-19 but it can also be mediated by the human CB1 cannabinoid receptor, indicating functional conservation between the nematode and mammalian endocannabinoid systems for regulation of food preferences. Furthermore, the effect of anandamide in C. elegans is bidirectional, as it increases appetitive and consummatory responses to superior food but decreases these responses to inferior food. This bidirectionality is mirrored at the cellular level. Anandamides behavioral effects require the AWC chemosensory neurons, and anandamide renders these neurons more sensitive to superior food and less sensitive to inferior food. Our findings reveal a surprising degree of functional conservation in the effects of endocannabinoids on hedonic feeding across species and establish a new system in which to investigate the cellular and molecular basis of endocannabinoid system function in the regulation of food choice. Macroscopic control of synchronous electrical signaling with chemically-excited gene expression Excitable cells can convert electrical signals into chemical outputs to facilitate the active transport of information across larger distances. This electrical-to-chemical conversion requires a tightly regulated expression of ion channels. Alterations of ion channel expression provide landmarks of numerous pathological diseases, such as cardiac arrhythmia, epilepsy, or cancer. Although the activity of ion channels can be locally regulated by external light or chemical stimulus, it remains challenging to coordinate the expression of ion channels on extended spatial-temporal scales in a non-invasive manner. Here, we have engineered yeast S. cerevisiae to read and convert local chemical concentrations into a dynamic electrical field distributed across cell populations. The core mechanism encodes a chemically-excitable dual-feedback gene circuit that precisely tunes the expression domain of potassium channels, globally coordinating cyclic firing of the plasma membrane potential (PMP). We demonstrate that this mechanism leverages an engineered constitutively open bacterial potassium channel KcsA to directly couple chemical stimuli with ion flux through gene expression and it can interface with the host ion channels through the pulsatile production of toxins. Our study provides a robust synthetic transcriptional toolbox underlying the conversion of local chemical environments into spatiotemporally organized electrical impulses for various cellular engineering, synthetic biology, and potential therapeutic applications. Diurnal activity budget and conservation implication of Hippopotamus amphibius at Bui National Park, Ghana Animals apportion time for their daily behavioral activities. We studied the activity budget of Hippopotamus amphibius at the Black Volta River in the Bui National Park, Ghana. We performed instantaneous scan sample using ground count survey between August 2020 to July 2021. We observed that H. amphibius spent most of their day time resting (54.75 %), followed by feeding (22.93 %), walking (19.2 %) and touching (3.12 %). We found significant difference in the time spent between the daily activities (H = 41.67; p < 0.0001). Except activities involving touching, we found no significant seasonal differences in feeding, resting and walking activities by H. amphibius. Understanding the behavioral activities of H. amphibius will enhance management and conservation of the animals. The first multi-tissue diel genome-scale metabolic model of a woody plant highlights the suberin biosynthesis pathways in Quercus suber In the last decade, genome-scale metabolic models have been increasingly used to study plant metabolic behavior at the tissue and multi-tissue level under different environmental conditions. Quercus suber, also known as the cork oak tree, is one of the most important forest communities of the Mediterranean/Iberian region. In this work, we present the genome-scale metabolic model of the Q. suber (iEC7871), the first of a woody plant. The metabolic model comprises 7871 genes, 6231 reactions, and 6481 metabolites across eight compartments. Transcriptomics data was integrated into the model to obtain tissue-specific models for the leaf, inner bark, and phellogen, with specific biomass compositions. The tissue-specific models were merged into a diel multi-tissue metabolic model to predict interactions among the three tissues at the light and dark phases. The metabolic models were also used to analyze the pathways associated with the synthesis of suberin monomers. Nevertheless, the models developed in this work can provide insights into other aspects of the metabolism of Q. suber, such as its secondary metabolism and cork formation. Bioinformatic and mechanistic analysis of the palmerolide PKS-NRPS biosynthetic pathway from the microbiome of an Antarctic ascidian Complex interactions exist between microbiomes and their hosts. Increasingly, defensive metabolites that have been attributed to host biosynthetic capability are now being recognized as products of host-associated microbes. These unique metabolites often have bioactivity targets in human disease and can be purposed as pharmaceuticals. Polyketides are a complex family of natural products that often serve as defensive metabolites for competitive or pro-survival purposes for the producing organism, while demonstrating bioactivity in human diseases as cholesterol lowering agents, anti-infectives, and anti-tumor agents. Marine invertebrates and microbes are a rich source of polyketides. Palmerolide A, a polyketide isolated from the Antarctic ascidian Synoicum adareanum, is a vacuolar-ATPase inhibitor with potent bioactivity against melanoma cell lines. The biosynthetic gene clusters (BGCs) responsible for production of secondary metabolites are encoded in the genomes of the producers as discrete genomic elements. A candidate palmerolide BGC was identified from a S. adareanum microbiome-metagenome based on a high degree of congruence with a chemical structure-based retrobiosynthetic prediction. Protein family homology analysis, conserved domain searches, active site and motif identification were used to identify and propose the function of the [~]75 kb trans-acyltransferase (AT) polyketide synthase-non-ribosomal synthase (PKS-NRPS) domains responsible for the stepwise synthesis of palmerolide A. Though PKS systems often act in a predictable co-linear sequence, this BGC includes multiple trans-acting enzymatic domains, a non-canonical condensation termination domain, a bacterial luciferase-like monooxygenase (LLM), and is found in multiple copies within the metagenome-assembled genome (MAG). Detailed inspection of the five highly similar pal BGC copies suggests the potential for biosynthesis of other members of the palmerolide chemical family. This is the first delineation of a biosynthetic gene cluster from an Antarctic microbial species, recently proposed as Candidatus Synoicohabitans palmerolidicus. These findings have relevance for fundamental knowledge of PKS combinatorial biosynthesis and could enhance drug development efforts of palmerolide A through heterologous gene expression. The structurome of a Clostridium difficile phage and the remarkable accurate prediction of its novel phage receptor-binding protein As natural bacterial predators, bacteriophages have the potential to be developed to tackle antimicrobial resistance, but our exploitation of them is limited by understanding their vast uncharacterised genetic diversity1,2. Fascinatingly, this genetic diversity reflects many ways that phages can make proteins, performing similar functions that together form the familiar phage particle. Critical to infection are phage receptor-binding proteins (RBPs) that bind bacterial  receptors and initiate bacterial entry3. Here we identified and characterised Gp22, a novel RBP for phage CDHS-1 that infects pathogenic C. difficile, but that had no recognisable RBPs. We showed that Gp22 antibodies neutralised CDHS-1 infection and used immunogold-labelling and transmission electron microscopy to identify their location on the capsid. The Gp22 three-dimensional structure was resolved by X-ray crystallography revealing a new RBP class with an N-terminal L-shaped -helical superhelix domain and a C-terminal Mg2+-binding domain. The findings provide novel insights into C. difficile phage biology and phage-host interactions. This will facilitate optimal phage development and future engineering strategies4,5. Furthermore, the AlphaFold2-predicted Gp22 structure, which was strikingly accurate, paves the way for a structurome based transformation and guidance of future phage studies where many proteins lack sequence homology but have recognisable protein structures. Plant-host shift, spatial persistence, and the viability of an invasive insect population Assessing the effects of a plant-host shift is important for monitoring insect populations over long time periods and for interventions in a conservation or pest management framework. In a heterogeneous environment, individuals may disperse between sources and sinks in order to persist. Here we propose a single-species two-patch model that aims to capture the generational movement of an insect that exhibits density-dependent dispersal, to see how shifting between hosts could alter its viability and asymptotic dynamics. We then analyse the stability and persistence properties of the model and further validate it using parameter estimates derived from laboratory experiments. In order to evaluate the potential of this model, we applied it to Drosophila suzukii (Diptera: Drosophilidae), which has become a harmful pest in several countries around the world. Although many studies have investigated the preference and attractiveness of potential hosts on this invasive drosophilid, no studies thus far have investigated whether a shift of fruit host could affect such a species ecological viability or spatiotemporal persistence. The model results show that a shift in host choice can significantly affect the growth potential and fecundity of a species such as D. suzukii, which ultimately could aid such invasive populations in their ability to persist within a changing environment. Neuroinvasion and neurotropism by SARS-CoV-2 variants in the K18-hACE2 mouse Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) not only affects the respiratory tract but also causes neurological symptoms such as loss of smell and taste, headache, fatigue or severe cerebrovascular complications. Using transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) we investigated the spatiotemporal distribution and pathomorphological features in the CNS following intranasal infection with SARS-CoV-2 variants, also after prior influenza A virus infection. Apart from Omicron, we found all variants to frequently spread to and within the CNS. Infection was restricted to neurons and appeared to spread from the olfactory bulb mainly in basally orientated regions in the brain and into the spinal cord, independent of ACE2 expression and without evidence of neuronal cell death, axonal damage or demyelination. However, microglial activation, microgliosis and a mild macrophage and T cell dominated inflammatory response was consistently observed, accompanied by apoptotic death of endothelial, microglial and immune cells, without their apparent infection. Microgliosis and immune cell apoptosis indicate a potential role of microglia for pathogenesis and viral effect in COVID-19 and possible impairment of neurological functions, especially in long COVID. These data may also be informative for the selection of therapeutic candidates, and broadly support investigation of agents with adequate penetration into relevant regions of the CNS. Correlation between speed and turning naturally arises for sparsely sampled cell movements Mechanisms regulating cell movement are not fully understood. One feature of cell movement that determines how far cells displace from an initial position is persistence, the ability to perform movements in a direction similar to the previous movement direction. Persistence is thus determined by turning angles between two sequential displacements and be characterized by an average turning angle or persistence time. Recent studies found that a cells average speed and turning are negatively correlated, suggesting a fundamental cell-intrinsic program whereby cells with a lower turning ability (i.e., larger persistence time) are intrinsically faster (or faster cells turn less). By simulating correlated or persistent random walks (PRWs) using two different frameworks (one based on von Mises-Fisher (vMF) distribution and another based on Ornstein-Uhlenbeck (OU) process) we show that the negative correlation between speed and turning naturally arises when cell trajectories are sub-sampled, i.e., when the frequency of sampling is lower than frequency at which cells make movements. This effect is strongest when the sampling frequency is on the order of magnitude with the typical cell persistence time and when cells vary in persistence time. Both conditions are observed for datasets of T cell movements in vivo that we have analyzed. In simulations the correlation arises due to randomness of cell movements resulting in highly variable persistence times for individual cells that with sub-sampling leads to large variability of average cell speeds. Interestingly, previously suggested methodology of calculating displacement of cell cohorts with different speeds resulted in similar results whether or not there is a cell-intrinsic correlation between cell speed and persistence. For both vMF- and OU-based simulations of PRWs we could find parameter values (distribution of persistence times, speeds, and sampling frequency) that matched experimentally measured correlations between speed and turning for two datasets of T cell movement in vivo suggesting that such simple correlations are not fully informative on the intrinsic link between speed and persistence. Our results thus suggest that sub-sampling may contribute to (and perhaps fully explains) the observed correlation between speed and turning at least for some cell trajectory data and emphasize the role of sampling frequency in inference of critical cellular parameters of cell motility such as speeds.

Secondary AbstractMeasurement of cell movements often results in a negative correlation between average speed and average turning angle suggesting an existence of a universal, cell-intrinsic movement program. We show that such a negative correlation may arise if cells in the population differ in their ability for persistent movement when the movement data are sub-sampled. We show that sub-sampling of cell trajectories generated using two different frameworks of persistent random walk can match the experimentally observed correlation between speed and turning for T cell movements in vivo. The pH Dependence of Niclosamide Solubility, Dissolution, and Morphology Motivates Potentially Universal Mucin-Penetrating Nasal and Throat Sprays for COVID19, its Contagious Variants, and Other Respiratory Viral Infections MotivationWith the coronavirus pandemic still raging, prophylactic nasal and early treatment throat sprays could help prevent infection and reduce viral load. Niclosamide has the potential to treat a broad range of viral infections if local bioavailability is optimized as mucin-penetrating solutions as opposed to more traditional microparticle-based sprays that cannot penetrate the mucin.

ExperimentalpH-dependence of supernatant concentrations and dissolution rates of niclosamide were measured in buffered solutions by Nanodrop-UV/Vis-spectroscopy for niclosamide from different suppliers, as precipitated material, and as cosolvates. Data was compared to predictions from Henderson Hasselbalch and precipitation pH models. Optimal microscopy was used to observe the morphologies of precipitated and converted niclosamide.

ResultsSupernatant-concentrations of niclosamide increased with increasing pH, from 1.77uM at pH 3.66 to 30uM at pH 8, and more rapidly from 90uM at pH8.5 to 300uM at pH 9.1, reaching 641uM at pH 9.5. Logarithmic rates for dissolution increased by [~]3x for pHs 8.62 to 9.44. However, when precipitated from supersaturated solution, niclosamide equilibrated to much lower final supernatant concentrations, reflective of more stable polymorphs at each pH that were also apparent for niclosamide from other suppliers and cosolvates.

ConclusionsGiven niclosamides activity against COVID19, its more contagious variants, and other respiratory viral infections, these niclosamide solutions, that put the virus in lockdown, could represent universal prophylactic nasal and early treatment throat sprays. As solutions they would be the simplest and potentially most effective formulations from both an efficacy standpoint as well as manufacturing and distribution, with no cold chain. They now just need testing. 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. Glutathione regulates transcriptional activation of iron transporters via S-nitrosylation of bHLH factors to modulate subcellular iron homeostasis While glutathione (GSH) is known to regulate iron (Fe) deficiency response in plants, its involvement in modulating subcellular Fe homeostasis remained elusive. In this study, we report that the GSH depleted mutants, cad2-1 and pad2-1 displayed increased sensitivity to Fe deficiency with significant down-regulation of the vacuolar Fe exporters, AtNRAMP3 and AtNRAMP4 and the chloroplast Fe importer, AtPIC1. Moreover, the pad2-1 mutant accumulated higher Fe content in vacuoles and lower in chloroplasts compared with Col-0 under Fe limited condition. Exogenous GSH treatment could enhance the chloroplast Fe content in Col-0 plants but failed to do so in the nramp3nramp4 double mutant demonstrating the role of GSH in modulating subcellular Fe homeostasis. Pharmacological experiments, mutant analysis and promoter assay revealed that this regulation involved transcriptional activation of the transporter genes by a GSH-GSNO module. The Fe responsive bHLH transcription factors (TFs), AtbHLH29, AtbHLH38 and AtbHLH101 were identified to interact with the promoters of these genes which were in turn activated via S-nitrosylation. Together, the present study delineates the role of GSH-GSNO module in regulating subcellular Fe homeostasis by transcriptional activation of the Fe transporters AtNRAMP3, AtNRAMP4 and AtPIC1 via S-nitrosylation of the bHLH TFs during Fe deficiency.

Summary statementGlutathione regulates subcellular iron homeostasis under iron deficiency via GSNO dependent transcriptional activation of AtNRAMP3, AtNRAMP4 and AtPIC1 genes by S-nitrosylation of the iron responsive bHLH transcription factors, AtbHLH29, AtbHLH38 and AtbHLH101. Royal knifefish generate powerful suction feeding through large neurocranial elevation and high epaxial muscle power Suction feeding in ray-finned fishes involves powerful buccal cavity expansion to accelerate water and food into the mouth. Previous XROMM studies in largemouth bass (Micropterus salmoides), bluegill sunfish (Lepomis macrochirus), and channel catfish (Ictalurus punctatus) have shown that more than 90% of suction power in high performance strikes comes from the axial musculature. Thus, the shape of the axial muscles and skeleton may impact suction feeding mechanics. Royal knifefish (Chitala blanci) have an unusual postcranial morphology, with a ventrally flexed vertebral column and relatively large mass of epaxial muscle. Based on their body shape, we hypothesized that royal knifefish would generate high power strikes by utilizing large neurocranial elevation, vertebral column extension, and epaxial shortening. As predicted, C. blanci generated high suction expansion power compared to the other three species studied to date (up to 160 W), which was achieved by increasing both the rate of volume change and the intraoral subambient pressure. The large epaxial muscle (25% of body mass) shortened at high velocities to produce large neurocranial elevation and vertebral extension (up to 41 deg, combined), as well as high muscle mass-specific power (up to 800 W kg-1). For the highest power strikes, axial muscles generated 95% of the power, and 64% of the axial muscle mass consisted of the epaxial muscles. The epaxial-dominated suction expansion of royal knifefish supports our hypothesis that postcranial morphology may be a strong predictor of suction feeding biomechanics.

SUMMARY STATEMENTRoyal knifefish rely on their distinct postcranial morphology--with a curved vertebral column and large dorsal body muscles--to produce large neurocranial elevation and powerful suction feeding. Mitochondrial genomes of two parasitic Cuscuta species lack clear evidence of horizontal gene transfer and retain unusually fragmented ccmFC genes BackgroundThe intimate association between parasitic plants and their hosts favours the exchange of genetic material, potentially leading to horizontal gene transfer (HGT) between plants. With the recent publication of several parasitic plant nuclear genomes, there has been considerable focus on such non-sexual exchange of genes. To enhance the picture on HGT events in a widely distributed parasitic genus, Cuscuta (dodders), we assembled and analyzed the organellar genomes of two recently sequenced species, C. australis and C. campestris, making this the first account of complete mitochondrial genomes (mitogenomes) for this genus.

ResultsThe mitogenomes are 265,696 and 275,898 bp in length and contain a typical set of mitochondrial genes, with ten missing or pseudogenized genes often lost from angiosperm mitogenomes. Each mitogenome also possesses a structurally unusual ccmFC gene, which exhibits splitting of one exon and a shift to trans-splicing of its intron. Based on phylogenetic analysis of mitochondrial genes from across angiosperms and similarity-based searches, there is little to no indication of HGT into the Cuscuta mitogenomes. A few candidate regions for plastome-to-mitogenome transfer were identified, with one suggestive of possible HGT.

ConclusionsThe lack of HGT is surprising given examples from the nuclear genomes, and may be due in part to the relatively small size of the Cuscuta mitogenomes, limiting the capacity to integrate foreign sequences. GUIdEStaR (G-quadruplex, uORF, IRES, Epigenetics, Small RNA, Repeats), the integrated metadatabase in conjunction with neural network methods GUIdEStaR integrates existing databases of various types of G-quadruplex, upstream Open Reading Frame (uORF), Internal Ribosome Entry Site (IRES), methylation to RNA and histone protein, small RNA, and repeats. GUIdEStaR consists of approx. 40,000 genes and 320,000 transcripts. An mRNA transcript is divided into 5 regions (5UTR, 3UTR, exon, intron, and biological region) where each region contains presence-absence data of 169 different types of elements. Recently, artificial intelligence (AI) based analysis of sequencing data has been gaining popularity in the area of bioinformatics. GUIdEStaR generates datasets that can be used as inputs to AI methods. At the GUIdEStaR homepage, users submit gene symbols by clicking a "Send" button, and shortly result files in CSV format are available for download at the result website. Users have an option to send the result files to their email addresses. Additionally, the entire database and the example Java codes are also freely available for download. Here, we demonstrate the database usage with three neural network classification studies-1) small RNA study for classifying transcription factor (TF) genes into either one of TF mediated by small RNA originated from SARS-CoV-2 or by human microRNA (miRNA), 2) cell membrane receptor study for classifying receptor genes as either with virus interaction or without one, and 3) nonsense mediated mRNA decay (NMD) study for classifying cell membrane and nuclear receptors as either NMD target or non-target. GUIdEStaR is available for access to the easy-to-use web-based database at www.guidestar.kr and for download at https://sourceforge.net/projects/guidestar. A platform of patient-derived microtumors identifies therapeutic vulnerabilities in ovarian cancer. BackgroundIn light of the frequent development of therapeutic resistance in cancer treatment, there is a strong need for personalized model systems representing patient tumor heterogeneity, while enabling parallel drug testing and identification of appropriate treatment responses in individual patients. Using ovarian cancer as a prime example of a heterogeneous tumor disease with complex microenvironment and high recurrence rates, we developed a 3D preclinical ovarian cancer model comprised of patient-derived microtumors (PDM) and autologous tumor-infiltrating lymphocytes (TILs) for identification of treatment vulnerabilities and validation of anti-cancer drug efficacy using immunohistochemistry, immune cell phenotyping, functional assays and protein profiling analyses.

MethodsPDM and TILs were isolated from fresh primary ovarian cancer tissue specimen using mechanical disruption and limited enzymatic digestion and were subsequently cultured in suspension in defined media in the absence of serum. The heterogeneous cellular composition of isolated PDM as well as autologous TILs was analyzed by FFPE immunohistochemistry and multi-color flow cytometry, respectively. For in-depth protein profiling of PDM we established Reverse Phase Protein Array (RPPA) analyses of >110 total and phospho-proteins. Treatment efficacy in response to chemotherapeutics as well as immunotherapeutic compounds was assessed in PDM and PDM-TIL co-cultures using a functional viability assay in microplate format.

ResultsThe enzymatic digestion of primary ovarian cancer tissue and suspension culture in defined serum-free media allowed fast and efficient recovery of patient-derived microtumors (PDM). Immunohistochemical analyses demonstrated histopathological comparability of ovarian cancer PDM with corresponding patient tumor tissue. Reverse Phase Protein Array (RPPA)-based analyses of >110 total and phospho-proteins enabled the identification of patient-specific sensitivities to standard, platinum-based therapy and thereby the prediction of potential treatment-responders. Finally, combining PDM and autologous TILs for individual efficacy testing of immune checkpoint inhibitor treatment demonstrated the potential for patient-specific enhancement of cytotoxic TIL activity by this therapeutic approach.

Conclusion3D patient-derived ovarian cancer microtumors represent a preclinical, ex vivo tumor model that reflects intertumoral heterogeneity and represent the cellular complexity of individual patient tumors. Combining protein pathway analysis and anti-cancer drug efficacy testing of PDM enables drug mode-of-action analyses and therapeutic sensitivity prediction within a clinically relevant time frame after surgery. Follow-up studies in larger cohorts are currently under way to further evaluate the applicability of this platform to support clinical decision-making and personalizing cancer treatment. Co-option of immune effectors by the hormonal signalling system triggering metamorphosis in Drosophila melanogaster Insect metamorphosis is regulated by the production, secretion and degradation of two peripheral hormones: 20-hydroxyecdysone (ecdysone) and juvenile hormone (JH). In addition to their roles in developmental regulation, increasing evidence suggests that these hormones are involved in innate immunity processes, such as phagocytosis and the induction of antimicrobial peptide (AMP) production. AMP regulation includes systemic responses as well as local responses at surface epithelia that contact with the external environment. At pupariation, Drosophila melanogaster increases dramatically the expression of three AMP genes, drosomycin (drs), drosomycin-like 2 (drsl2) and drosomycin-like 5 (drsl5). We show that the systemic action of drs at pupariation is dependent on ecdysone signalling in the fat body and operates via the ecdysone downstream target, Broad-Z4. In parallel, ecdysone also regulates local responses, specifically through the activation of drsl2 expression in the gut. Finally, we confirm the relevance of this ecdysone dependent AMP expression for the control of bacterial persistence by showing that flies lacking drs expression in the fat body have higher bacterial persistence over metamorphosis. Together, our data emphasize the importance of the association between ecdysone signalling and immunity using in vivo studies and establish a new role for ecdysone at pupariation, which impacts developmental success by regulating the immune system in a stage-dependent manner. We speculate that this co-option of immune effectors by the hormonal system may constitute a pre-emptive mechanism to control bacterial numbers in the pupa at the core of metamorphosis evolution. P2X2 receptor subunit interfaces are missense variant hotspots where mutations tend to increase apparent ATP affinity Background and PurposeP2X receptors (P2XRs) are trimeric ligand-gated ion channels (LGICs) that open a cation-selective pore in response to ATP binding to their large extracellular domain (ECD). The seven known P2XR subtypes can assemble as homo- or heterotrimeric complexes and contribute to numerous physiological functions, including nociception, inflammation and hearing. The overall structure of P2XRs is well established, but little is known about the spectrum and prevalence of human genetic variations and the functional implications in specific domains.

Experimental ApproachHere we examine the impact of P2X2 receptor (P2X2R) inter-subunit interface missense variants identified in the human population or through structural predictions. We test both single and double mutants through electrophysiological and biochemical approaches.

Key resultsWe demonstrate that predicted ECD inter-subunit interfaces display a higher-than-expected density of missense variations and that the majority of mutations that disrupt putative inter-subunit interactions result in channels with higher apparent ATP affinity. Lastly, we show that double mutants at the subunit interface show significant energetic coupling, especially if located in close proximity.

Conclusions and ImplicationsWe provide the first structural mapping of the mutational burden across the human population in a LGIC and show that the density of missense mutations is constrained between protein domains, indicating evolutionary selection at the domain level. Our data may indicate that, unlike other LGICs, P2X2Rs have evolved an intrinsically high threshold for activation, possibly to allow for additional modulation or as a cellular protection mechanism against overstimulation.

Bullet point summary What is already known: O_LIP2X2 receptors are ATP-activated ion channels implicated in hearing and nociceptice pathways
C_LI
 What this study adds: O_LIA structural mapping of missense variants observed in the human population
C_LIO_LIWe identify the intersubunit-interface as a variant hotspot and decipher functional impact of mutations
C_LI
 Clinical significance: O_LIThe development of both inhibitors and activators of P2X2 receptor function may be required
C_LI Itaconate and its derivatives repress C2C12 myogenesis A Krebs cycle intermediate metabolite, itaconate, has gained attention as a potential antimicrobial and autoimmune disease treatment due to its anti-inflammatory effects. While itaconate and its derivatives pose an attractive therapeutic option for the treatment of inflammatory diseases, the effects outside the immune system still remain limited, particularly in the muscle. Therefore, we endeavored to determine if itaconate signaling impacts muscle differentiation. Utilizing the well-established C2C12 model of in vitro myogenesis, we evaluated the effects of itaconate and its derivatives on transcriptional and protein markers of muscle differentiation as well as mitochondrial function. We found itaconate and the derivatives dimethyl itaconate and 4-octyl itaconate disrupt differentiation media-induced myogenesis. A primary biological effect of itaconate is a succinate dehydrogenase (SDH) inhibitor. We find the SDH inhibitors dimethyl malonate and harzianopyridone phenocopie the anti-myogenic effects of itaconate. Furthermore, we find treatment with exogenous succinate results in blunted myogenesis. Together our data indicate itaconate and its derivatives interfere with in vitro myogenesis, potentially through inhibition of SDH and subsequent succinate accumulation. More importantly, our findings suggest the therapeutic potential of itaconate and its derivatives could be limited due to deleterious effects on myogenesis. Discovery and functional assessment of a novel adipocyte population driven by intracellular Wnt/β-catenin signaling in mammals Wnt/{beta}-catenin signaling has been well established as a potent inhibitor of adipogenesis. Here, we identified a population of adipocytes that exhibit persistent activity of Wnt/{beta}-catenin signaling, as revealed by the Tcf/Lef-Gfp reporter allele, in embryonic and adult mouse fat depots, named as Wnt+ adipocytes. We showed that the {beta}-catenin-mediated signaling activation in these cells is Wnt ligand- and receptor-independent but relies on AKT/mTOR pathway and is essential for cell survival. Such adipocytes are distinct from classical ones in transcriptomic and genomic signatures and can be induced from various sources of mesenchymal stromal cells including human cells. Genetic lineage-tracing and targeted cell ablation studies revealed that these adipocytes convert into beige adipocytes directly and are also required for beige fat recruitment under thermal challenge, demonstrating both cell autonomous and non-cell autonomous roles in adaptive thermogenesis. Furthermore, mice bearing targeted ablation of these adipocytes exhibited glucose intolerance, while mice receiving exogenously supplied such cells manifested enhanced glucose utilization. Our studies uncover a unique adipocyte population in regulating beiging in adipose tissues and systemic glucose homeostasis. Biophysical characterization of the oligomeric states of recombinant Immunoglobulins type-M and their C1q binding kinetics by Biolayer Interferometry The Immunoglobulins type-M (IgMs) are one of the first antibody classes mobilized during immune responses against pathogens and tumor cells. Binding to specific target antigens enables the interaction with the C1 complex which strongly activates the classical complement pathway. This biological function is the basis for the huge therapeutic potential of IgMs but due to their high oligomeric complexity, in vitro production, biochemical and biophysical characterizations are challenging. In the present study, we present recombinant production of two IgM models (IgM617 and IgM012) in pentameric and hexameric states and the evaluation of their polymer distribution using different biophysical methods (AUC, SEC-MALLS, Mass Photometry and Transmission Electron Microscopy). Each IgM oligomer has individual specific expression pattern and yield with different protein quality likely due to intrinsic IgM properties and patterning. Nevertheless, the purified recombinant IgMs retain their ability to activate complement in a C1q dependent manner. And more importantly, a new method to evaluate their functional quality attribute by characterizing the kinetics of C1q binding to recombinant IgM has been developed using BioLayer Interferometry (BLI). We show that recombinant IgMs possess similar C1q binding properties as IgMs purified from human plasma.