Including random effects in statistical models in ecology: fewer than five levels? As generalized linear mixed-effects models (GLMMs) have become a widespread tool in ecology, the need to guide the use of such tools is increasingly important. One common guideline is that one needs at least five levels of a random effect. Having such few levels makes the estimation of the variance of random effects terms (such as ecological sites, individuals, or populations) difficult, but it need not muddy ones ability to estimate fixed effects terms - which are often of primary interest in ecology. Here, I simulate ecological datasets and fit simple models and show that having too few random effects terms does not influence the parameter estimates or uncertainty around those estimates for fixed effects terms. Thus, it should be acceptable to use fewer levels of random effects if one is not interested in making inference about the random effects terms (i.e. they are  nuisance parameters used to group non-independent data). I also use simulations to assess the potential for pseudoreplication in (generalized) linear models (LMs), when random effects are explicitly ignored and find that LMs do not show increased type-I errors compared to their mixed-effects model counterparts. Instead, LM uncertainty (and p values) appears to be more conservative in an analysis with a real ecological dataset presented here. These results challenge the view that it is never appropriate to model random effects terms with fewer than five levels - specifically when inference is not being made for the random effects, but suggest that in simple cases LMs might be robust to ignored random effects terms. Given the widespread accessibility of GLMMs in ecology and evolution, future simulation studies and further assessments of these statistical methods are necessary to understand the consequences of both violating and blindly following simple guidelines.
Pathways of savannization in a mesic African savanna-forest mosaic O_LIFires in savannas limit tree cover, thereby promoting flammable grass accumulation and fuelling further frequent fires. Meanwhile, forests and thickets form dense canopies that reduce C4-grass fuel loads and creating a humid microclimate, thereby excluding fires under typical climatic conditions.
C_LIO_LIHowever, extreme fires occasionally burn into these closed-canopy systems. Although these rare fires cause substantial tree mortality and can make repeat fires more likely, the long-term consequences of an extreme fire for closed canopy vegetation structure and potential to convert to savanna (hereafter "savannization") remain largely unknown.
C_LIO_LIHere, we analysed whether an extreme fire could, alone, alter species composition, vegetation structure, and fire regimes of closed-canopy ecosystems in an intact savanna-forest-thicket mosaic, or whether successive fires after an initial extreme fire were necessary to trigger a biome transition between from forest to savanna.
C_LIO_LIWe found that forests that only burned once recovered, whereas those that burned again following an initial extreme fire transitioned from closed-canopy forests towards open, grassy savannas.
C_LIO_LIWhile thickets had less tree mortality in fires than forests, repeat fires nonetheless precipitated a transition towards savannas.
C_LIO_LIColonization of the savanna tree community lagged behind the grass community, but also began to transition.
C_LI

SynthesisOur results suggest that rare extreme fires, followed by repeated burning can indeed result in savannization in places where savanna and forest represent alternative stable states.
Characterization of methylation profiles in spontaneous preterm birth placental villous tissue Preterm birth is a global public health crisis which results in significant neonatal and maternal mortality. Yet little is known regarding the molecular mechanisms of idiopathic spontaneous preterm birth, and we have few diagnostic markers for adequate assessment of placental development and function. Previous studies of placental pathology and our transcriptomics studies suggest a role for placental maturity in idiopathic spontaneous preterm birth. It is known that placental DNA methylation changes over gestation. We hypothesized that if placental hypermaturity is present in our samples, we would observe a unique idiopathic spontaneous preterm birth DNA methylation profile potentially driving the gene expression differences we previously identified in our placental samples. Our results indicate the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern suggesting hypermaturity. Only seven significant differentially methylated regions fitting the idiopathic spontaneous preterm birth specific (relative to the controls) profile were identified, indicating unusually high similarity in DNA methylation between idiopathic spontaneous preterm birth and term birth samples. We identified an additional 1,718 significantly methylated regions in our gestational age matched controls were the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern, again indicating a striking level of similarity between the idiopathic spontaneous preterm birth and term birth samples. Pathway analysis of these regions revealed differences in genes within the WNT and Cadherin signaling pathways, both of which are essential in placental development and maturation. Taken together, these data demonstrate that the idiopathic spontaneous preterm birth samples are molecularly more mature than expected given their respective gestational age which likely impacts birth timing. First evidence of in vitro cytotoxic effects of marine microlitter on Merluccius merluccius and Mullus barbatus, two Mediterranean commercial fish species Marine litter is composed mainly of plastics and is recognized as a serious threats to marine ecosystems. Ecotoxicological approaches have started elucidating the potential severity of microplastics (MPs) in controlled laboratory studies with pristine materials but no information exist on marine environmental microlitter as a whole. Here, we characterized the litter in the coastal Northern Tyrrhenian sea and in the stomach of two fish species of socio-economic importance, and exposed primary cell cultures of mucosal and lymphoid organs to marine microlitter for evaluating possible cytotoxic effects. An average of 0.30 {+/-} 0.02 microlitter items m-3 was found in water samples. FT-IR analysis revealed that plastic particles, namely HDPE, polyamide and polypropylene were present in 100% and 83.3% of Merluccius merluccius and Mullus barbatus analyzed, which overall ingested 14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual, respectively. Moreover, microlitter was confirmed as a vector of microorganisms. Lastly, the apical end-point of viability was found to be significantly reduced in splenic cells exposed in vitro to two microlitter conditions. Considering the role of the spleen in the mounting of adaptive immune responses, our results warrant more in-depth investigations for clarifying the actual susceptibility of these two species to anthropogenic microlitter.

HighlightsO_LI0.30 {+/-} 0.02 microlitter items m-3 were found at the surface of coastal Northern Tyrrhenian sea
C_LIO_LI14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual were retrieved from the stomach of hakes and mullets
C_LIO_LIThe ingested microlitter contained plastic items
C_LIO_LIMicrolitter was validated as a carrier of bacteria, fungi and flagellates
C_LIO_LISplenic cells exposed to two microlitter conditions for 72 hours suffered cytotoxicity
C_LI CHRM1 is a Druggable Melanoma Target Whose Endogenous Activity is Determined by Inherited Genetic Variation in DOPA Production Melanoma risk is 30 times higher in people with lightly pigmented skin compared to those with darkly pigmented skin. Here we show that this difference results from more than melanin pigment and its ultraviolet radiation (UVR) shielding effect. Using primary human melanocytes representing the full human skin pigment continuum and several preclinical melanoma models, we show that cell-intrinsic differences between dark and light melanocytes regulate melanocyte proliferative capacity, overall cellular differentiation state, and susceptibility to malignant transformation, independently of melanin and UV exposure. We determined that these differences result from dihydroxyphenylalanine (DOPA), a melanin precursor synthesized at higher levels in melanocytes from dark skin. Although DOPA was not previously known to have specific signaling activity, we used both high throughput pharmacologic and genetic in vivo CRISPR screens to determine that DOPA limits melanocyte and melanoma cell proliferation by directly inhibiting the muscarinic acetylcholine receptor M1 (CHRM1), a G Protein-Coupled Receptor (GPCR) not previously known to bind DOPA, nor to affect melanoma pathobiology. Pharmacologic CHRM1 antagonism in melanoma leads to depletion of c-Myc and FOXM1, both of which are proliferation drivers associated with aggressive melanoma. In preclinical mouse melanoma models using both immune deficient and syngeneic immune competent mice, pharmacologic inhibition of CHRM1 or FOXM1 inhibited tumor growth. CHRM1 and FOXM1 may be new therapeutic targets for melanoma. Correlational networking guides the discovery of cryptic natural product biosynthetic enzymes Bacterial natural product biosynthetic genes, canonically clustered, have been increasingly found to rely on cryptic enzymes encoded elsewhere in the genome for completion of biosynthesis. The study and application of lanthipeptides are frequently hindered by cryptic proteases required for final maturation. Here, we establish a global correlation network bridging the gap between lanthipeptide precursors and cryptic proteases. Applying our analysis to 161,954 bacterial genomes, we establish 5,209 correlations between precursors and cryptic proteases, with 91 prioritized. We use network predictions and co-expression analysis to reveal a previously missing protease for the maturation of class I lanthipeptide paenilan. We further discover widely distributed bacterial M16B metallopeptidases of previously unclear biological function as a new family of lanthipeptide proteases. We show the involvement of a pair of bifunctional M16B proteases in the production of novel class III lanthipeptides with high substrate specificity. Together, these results demonstrate the strength of our correlational networking approach to the discovery of cryptic lanthipeptide proteases and potentially other missing enzymes for natural products biosynthesis.

Graphical abstract O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY Ineffectual AEC1 differentiation from KRT8hi transitional state without fibrosis is associated with fatal COVID-19 ARDS 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. Differed cell division angle, position of cell proliferative area, and localization of ANGUSTIFOLIA3 in lateral organs Leaf meristem is a cell proliferative zone present in the lateral organ primordia. In this study, we investigated how the proliferative zone affects the final morphology of the lateral organs. We examined how cell proliferative zones differ in the primordia of planar floral organs and polar auxin transport inhibitor (PATI)-treated leaves from normal foliage leaf primordia of Arabidopsis thaliana with a focus on the spatial accumulation pattern of ANGUSTIFOLIA3 (AN3), a key element for leaf meristem positioning. We found that organ shape changes by PATI treatment were correlated to the angle of the cell division plane relative to the leaf primordia axis in the leaf meristem (cell division angle), but not with leaf-meristem positioning, size of the leaf meristem, or the localization pattern of AN3 protein. In contrast, different shapes between sepals and petals compared with foliage leaves were associated with both altered meristem position associated with altered AN3 expression patterns and different distributions of cell division angles. These results suggest that lateral organ shapes are regulated via two aspects: position of meristem and cell division angles

Summary statementDifferent lateral organs with different morphology possess different properties of meristems; cell division angles, position of cell proliferative area and AN3 localization patterns. Detecting DNA of novel fungal pathogens using ResNets and a curated fungi-hosts data collection BackgroundEmerging pathogens are a growing threat, but large data collections and approaches for predicting the risk associated with novel agents are limited to bacteria and viruses. Pathogenic fungi, which also pose a constant threat to public health, remain understudied. Relevant data remains comparatively scarce and scattered among many different sources, hindering the development of sequencing-based detection workflows for novel fungal pathogens. No prediction method working for agents across all three groups is available, even though the cause of an infection is often difficult to identify from symptoms alone.

ResultsWe present a curated collection of fungal host range data, comprising records on human, animal and plant pathogens, as well as other plant-associated fungi, linked to publicly available genomes. We show that it can be used to predict the pathogenic potential of novel fungal species directly from DNA sequences with either sequence homology or deep learning. We develop learned, numerical representations of the collected genomes and visualize the landscape of fungal pathogenicity. Finally, we train multi-class models predicting if next-generation sequencing reads originate from novel fungal, bacterial or viral threats.

ConclusionsThe neural networks trained using our data collection enable accurate detection of novel fungal pathogens. A curated set of over 1,400 genomes with host and pathogenicity metadata supports training of machine learning models and sequence comparison, not limited to the pathogen detection task.

AvailabilityThe data, models and code are hosted at https://zenodo.org/record/5846345, https://zenodo.org/record/5711877, and https://gitlab.com/dacs-hpi/deepac. 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. New Segmentation and Feature Extraction Algorithm for Classification of White Blood Cells in Peripheral Smear Images This article addresses a new method for the classification of white blood cells (WBCs) using image processing techniques and machine learning methods. The proposed method consists of three steps: detecting the nucleus and cytoplasm, extracting features, and classification. At first, a new algorithm is designed to segment the nucleus. For the cytoplasm to be detected, only a part of it located inside the convex hull of the nucleus is involved in the process. This attitude helps us overcome the difficulties of segmenting the cytoplasm. In the second phase, three shapes and four novel color features are devised and extracted. Finally, by using an SVM model, the WBCs are classified. The segmentation algorithm can detect the nucleus with a dice similarity coefficient of 0.9675. The proposed method can categorize WBCs in Raabin-WBC, LISC, and BCCD datasets with accuracies of 94.65 %, 92.21 %, and 94.20 %, respectively. It is worth mentioning that the hyperparameters of the classifier are fixed only with the Raabin-WBC dataset, and these parameters are not readjusted for LISC and BCCD datasets. The obtained results demonstrate that the proposed method is robust, fast, and accurate. A computational modeling approach for predicting multicell patterns based on signaling-induced differential adhesion Differential adhesion within cell populations enables the emergence of unique patterns in heterogeneous multicellular systems. This process has been previously explored using synthetically engineered heterogenous multicell spheroid systems, in which cell subpopulations engage in bidirectional intercellular signaling to regulate the expression of different cadherins. While engineered cell systems provide excellent experimental tools to observe pattern formation in cell populations, computational models may be leveraged to more systematically explore the key parameters that drive the emergence of different patterns. We developed and validated two- and three-dimensional agent-based models (ABMs) of spheroid patterning for cells engineered with a bidirectional signaling circuit that regulates N- and P-cadherin expression. The model was used to predict how varying initial cell seedings, cadherin induction probabilities, or homotypic adhesion strengths between cells impact spheroid patterning, and unsupervised machine learning techniques were used to map system parameters to unique spheroid patterns. The model was then deployed to design new synthetic cell signaling circuits based on a desired final multicell pattern. Force redistribution during clathrin-mediated endocytosis revealed by new phase-separating force sensors Forces are central to countless cellular processes, yet in vivo force measurement at the molecular scale remains difficult if not impossible. During clathrin-mediated endocytosis, forces produced by the actin cytoskeleton are transmitted to the plasma membrane by a multi-protein coat for membrane deformation. However, the magnitudes of these forces remain unknown. Here, we present new in vivo force sensors that phase separate under force. We measured the forces on the fission yeast HIP1R homologue End4p, a protein that links the membrane to the actin cytoskeleton. End4p is under [~]15 pN force near the actin cytoskeleton, [~] 11 pN near the clathrin lattice, and [~]9 pN near the plasma membrane. Our results demonstrate that forces are collected and redistributed across the endocytic machinery.

One-Sentence SummaryPhase-separating coiled-coil force sensors reveal an unexpected mode of force transduction during endocytosis in vivo. Embryonic exposure to cannabidiol disrupts active neural circuits, an effect increased by Δ-9-tetrahydrocannabinol and involving CB1R and CB2R in zebrafish In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. Ensovibep, a novel trispecific DARPin candidate that protects against SARS-CoV-2 variants SARS-CoV-2 has infected millions of people globally and continues to undergo evolution. Emerging variants can be partially resistant to vaccine induced immunity and therapeutic antibodies, emphasizing the urgent need for accessible, broad-spectrum therapeutics. Here, we report a comprehensive study of ensovibep, the first trispecific clinical DARPin candidate, that can simultaneously engage all three units of the spike protein trimer to potently inhibit ACE2 interaction, as revealed by structural analyses. The cooperative binding of the individual modules enables ensovibep to retain inhibitory potency against all frequent SARS-CoV-2 variants, including Omicron BA.1 and BA.2, as of February 2022. Moreover, viral passaging experiments show that ensovibep, when used as a single agent, can prevent development of escape mutations comparably to a cocktail of monoclonal antibodies (mAb). Finally, we demonstrate that the very high in vitro antiviral potency also translates into significant therapeutic protection and reduction of pathogenesis in Roborovski dwarf hamsters infected with either the SARS-CoV-2 wild-type or the Alpha variant. In this model, ensovibep prevents fatality and provides substantial protection equivalent to the standard of care mAb cocktail. These results support further clinical evaluation and indicate that ensovibep could be a valuable alternative to mAb cocktails and other treatments for COVID-19. 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. Citation needed?Wikipedia bibliometrics during the first waveof the COVID-19 pandemic. BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. Force correction analysis method for derivation of multidimensional free energy landscapes from adaptively-biased replica simulations A methodology is proposed for the calculation of multidimensional free-energy landscapes of molecular systems, based on analysis of multiple Molecular Dynamics trajectories wherein adaptive biases have been applied to enhance the sampling of different collective variables. In this approach, which we refer to as Force Correction Analysis Method (FCAM), local averages of the total and biasing forces are evaluated post-hoc, and the latter are subtracted from the former to obtain unbiased estimates of the mean force across collective-variable space. Multidimensional free-energy surfaces and minimum free-energy pathways are then derived from integration of the mean force landscape through kinetic Monte Carlo algorithm. To evaluate the proposed method, a series of numerical tests and comparisons with existing approaches were carried out for small molecules, peptides, and proteins, based on all-atom trajectories generated with standard, concurrent and replica-exchange Metadynamics in collective-variable spaces ranging from one- to six-dimensional. The tests confirm the correctness of the FCAM formulation and demonstrate that calculated mean forces and free energies converge rapidly and accurately, outperforming other methods used to unbias this kind of simulation data.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=102 SRC="FIGDIR/small/431654v5_ufig1.gif" ALT="Figure 1">
View larger version (45K):
org.highwire.dtl.DTLVardef@1412289org.highwire.dtl.DTLVardef@1425d53org.highwire.dtl.DTLVardef@1fcf5d0org.highwire.dtl.DTLVardef@6e6941_HPS_FORMAT_FIGEXP  M_FIG TOC/Abstract Graphic

C_FIG SufB intein splicing in Mycobacterium tuberculosis is influenced by two remote conserved N-extein Histidines Inteins are auto-processing domains that implement a multi-step biochemical reaction termed protein splicing, marked by cleavage and formation of peptide bonds. They excise from a precursor protein, generating a functional protein via covalent bonding of flanking exteins. We report the kinetic study of splicing and cleavage reaction in [Fe-S] cluster assembly protein SufB from Mycobacterium tuberculosis. Although it follows a canonical intein splicing pathway, distinct features are added by extein residues present in the active site. Sequence analysis identified two conserved histidines in the N-extein region; His-5 and His-38. Kinetic analyses of His-5Ala and His-38Ala SufB mutants exhibited significant reductions in splicing and cleavage rates relative to the SufB wild-type precursor protein. Structural analysis and molecular dynamics simulations suggested that Mtu SufB displays a unique mechanism where two remote histidines work concurrently to facilitate N-terminal cleavage reaction. His-38 is stabilized by the solvent-exposed His-5, and can impact N-S acyl shift by direct interaction with the catalytic Cys1. Development of inteins as biotechnological tools or as pathogen specific novel antimicrobial targets requires a more complete understanding of such unexpected roles of conserved extein residues in protein splicing. Antibiotic inhibition of the Plasmodium apicoplast decreases haemoglobin degradation and antagonises dihydroartemisinin action The World Health Organisation (WHO) recommends artemisinin (ART) combinations for treatment of uncomplicated Plasmodium falciparum malaria. Understanding the interaction between co-administered drugs within combination therapies is clinically important to prevent unintended consequences. The WHO guidelines recommend second line treatments that combine artesunate with tetracycline, doxycycline, or clindamycin--antibiotics that target the Plasmodium relict plastid, the apicoplast. In addition, antibiotics can be used simultaneously against other infectious diseases, leading to their inadvertent combination with ARTs. One consequence of apicoplast inhibition is a perturbation to haemoglobin uptake and trafficking--a pathway required for activation of ART derivatives. Here, we show that apicoplast-targeting antibiotics reduce the abundance of the catalyst of ART activation (free haem) in P. falciparum, likely through diminished haemoglobin digestion. We demonstrate antagonism between ART and these antibiotics, suggesting that apicoplast inhibitors reduce ART activation. These data have potential clinical implications due to the reliance on--and widespread use of--both ARTs and these antibiotics in malaria endemic regions. Oct1 recruits the histone lysine demethylase Utx to promote lineage specification The pathways used by cells to transition between undifferentiated, pluripotent state and tissue-specific states are incompletely understood. Here we show that the widely-expressed transcription factor Oct1/Pou2f1 activates silent, developmental lineage-appropriate genes to "canalize" developmental progression. Using Oct1 inducible knockout embryonic stem cells, we show that that Oct1 deficiency impairs mesodermal and terminal muscle differentiation in a manner that can be rescued by Oct1 retroviral expression. Using bulk RNA-seq, we show that mesoderm-specific genes are not correctly induced early in the differentiation timecourse. Single-cell gene expression profiling reveals that Oct1-deficient cells lose coherence in temporal induction of lineage programs, and show inappropriate developmental lineage branching resulting in poorly differentiated cells state with epithelial characteristics and hallmarks of oxidative stress. In embryonic stem cells, Oct1 co-binds with Oct4 to genes critical for mesoderm induction. The Utx/Kdm6a histone lysine demethylase also binds to many of these genes, and using a prototypic Pax3 gene we show that Oct1 recruits Utx to remove inhibitory H3K27me3 marks and activate expression. The specificity of the ubiquitous Oct1 protein for mesodermal genes can be explained by cooperative interactions with lineage-driving Smad transcription factors, as we show that Smad and Oct binding sites frequently coexist mesoderm-specific genes, that Oct1 and Smad3 interact, and that the sites and factors act cooperatively at the Myog enhancer. Overall, these results identify Oct1 as a key mediator of the induction of mesoderm lineage-specific genes. RapidoPGS: A rapid polygenic score calculator for summary GWAS data without a test dataset MotivationPolygenic scores (PGS) aim to genetically predict complex traits at an individual level. PGS are typically trained on genome-wide association summary statistics and require an independent test dataset to tune parameters. More recent methods allow parameters to be tuned on the training data, removing the need for independent test data, but approaches are computationally intensive. Based on fine-mapping principles, we present RapidoPGS, a flexible and fast method to compute PGS requiring summary-level GWAS datasets only, with little computational requirements and no test data required for parameter tuning.

ResultsWe show that RapidoPGS performs slightly less well than two out of three other widely-used PGS methods (LDpred2, PRScs, and SBayesR) for case-control datasets, with median r2 difference: -0.0092, -0.0042, and 0.0064, respectively, but up to 17,000-fold faster with reduced computational requirements. RapidoPGS is implemented in R and can work with user-supplied summary statistics or download them from the GWAS catalog.

Availability and implementationOur method is available with a GPL license as an R package from GitHub. Bulk-assembly of monodisperse coacervates and giant unilamellar vesicles with programmable hierarchical complexity Compartmentalization is an essential step for the emergence of protocellular life on the early earth and the endeavor to bottom-up construct synthetic cells in the lab. Among manifold natural strategies for compartmentalization, assembly of lipid vesicles or coacervate droplets with resemblance to cell envelope or cytoplasm have emerged as two dominant paradigms. However, the spontaneous assembly of lipid vesicles or coacervates generally results in polydisperse microcompartments and lacks the ability to efficiently integrate functional building blocks towards defined complexity, which give rise to the difficulty to achieve the order and complexity towards the scope of that of life via spontaneous assembly. Herein, we show that the interplay of coacervates and colloidal particles provides a bulk-assembly approach to form monodisperse coacervates and giant unilamellar vesicles (GUVs) with programmable hierarchical complexity. The hierarchical complexity of the synthetic cell microcompartments can be delicately engineered from the control of the spatial or temporal organization of components in entity unit to manifold forms of interconnected synthetic cell consortia, achieving a compartment hierarchy analogous to the organelle-cell-tissue structure. The increase in complexity gives rise to emergent properties, for instance, collective morphology evolution of coacervate assemblages, gated permeability of GUVs without sophisticated protein machinery, and remarkable structural and functional stability of GUVs at extreme conditions. This work paves an unprecedented step to form monodisperse, hierarchical coacervates and GUVs via bulk-assembly, which provides a prebiotically plausible approach for the proper integration of inanimate matter towards the order and complexity of life and monodisperse, hierarchical microdroplets for follow-up scientific and technological applications. Fast, high-throughput production of improved rabies viral vectors for specific, efficient and versatile transsynaptic retrograde labeling To understand the function of neuronal circuits it is crucial to disentangle the connectivity patterns within the network. However, most tools currently used to explore connectivity have low-throughput, low-selectivity, or limited accessibility. Here, we report the development of an improved packaging system for the production of the highly neurotrophic RVdGenvA-CVS-N2c rabies viral vectors, yielding titers orders of magnitude higher with no background contamination, at a fraction of the production time, while preserving the efficiency of transsynaptic labeling. Along with the production pipeline, we developed suites of "starter" AAV and bicistronic RVdG-CVS-N2c vectors, enabling retrograde labeling from a wide range of neuronal populations, tailored for diverse experimental requirements. We demonstrate the power and flexibility of the new system by uncovering hidden local and distal inhibitory connections in the hippocampal formation and by imaging the functional properties of a cortical microcircuit across weeks. Our novel production pipeline provides a convenient approach to generate new rabies vectors, while our toolkit flexibly and efficiently expands the current capacity to label, manipulate and image the neuronal activity of interconnected neuronal circuits in vitro and in vivo. Detecting phylodiversity-dependent diviersification with a general phylogenetic inference framework AO_SCPLOWBSTRACTC_SCPLOWDiversity-dependent diversification models have been extensively used to study the effect of ecological limits and feedback of community structure on species diversification processes, such as speciation and extinction. Current diversity-dependent diversification models characterise ecological limits by carrying capacities for species richness. Such ecological limits have been justified by niche filling arguments: as species diversity increases, the number of available niches for diversification decreases.

However, as species diversify they may diverge from one another phenotypically, which may open new niches for new species. Alternatively, this phenotypic divergence may not affect the species diversification process or even inhibit further diversification. Hence, it seems natural to explore the consequences of phylogenetic diversity-dependent (or phylodiversity-dependent) diversification. Current likelihood methods for estimating diversity-dependent diversification parameters cannot be used for this, as phylodiversity is continuously changing as time progresses and species form and become extinct.

Here, we present a new method based on Monte Carlo Expectation-Maximization (MCEM), designed to perform statistical inference on a general class of species diversification models and implemented in the R package emphasis. We use the method to fit phylodiversity-dependent diversification models to 14 phylogenies, and compare the results to the fit of a richness-dependent diversification model. We find that in a number of phylogenies, phylogenetic divergence indeed spurs speciation even though species richness reduces it. Not only do we thus shine a new light on diversity-dependent diversification, we also argue that our inference framework can handle a large class of diversification models for which currently no inference method exists. Rattus Economicus: A neural model for rational regulation of water-seeking effort in rodents In the laboratory, animals motivation to work tends to be positively correlated with reward magnitude. But in nature, rewards earned by work are essential to survival (e.g., working to find water), and the payoff of that work can vary on long timescales (e.g., seasonally). Under these constraints, the strategy of working less when rewards are small could be fatal. We found that instead, rats in a closed economy did more work for water rewards when the rewards were stably smaller, a phenomenon also observed in human labor supply curves. Like human consumers, rats also showed elasticity of demand, consuming far more water per day when its price in effort was lower. The neural mechanisms underlying such "rational" market behaviors remain largely unexplored. We propose a dynamic utility maximization model that can account for the dependence of rat labor supply (trials/day) on the wage rate (ml/trial), and also predict the temporal dynamics of when rats work. Based on data from mice, we hypothesize that SFOGLUT neurons in lamina terminalis continuously compute the instantaneous marginal utility of voluntary work for water reward, and causally determine the amount and timing of work.