EXPERT: Transfer Learning-enabled context-aware microbial source tracking Microbial source tracking quantifies the potential origin of microbial communities, facilitates better understanding of how the taxonomic structure and community functions were formed and maintained. However, previous methods involve a tradeoff between speed and accuracy, and have encountered difficulty in source tracking under many context-dependent settings. Here, we present EXPERT for context-aware microbial source tracking, in which we adopted a Transfer Learning approach to profoundly elevate and expand the applicability of source tracking, enabling biologically informed novel microbial knowledge discovery. We demonstrate that EXPERT can predict microbial sources with performance superior to other methods in efficiency and accuracy. More importantly, we demonstrate EXPERTs context-aware ability on several applications, including tracking the progression of infant gut microbiome development and monitoring the changes of gut microbiome for colorectal cancer patients. Broadly, transfer learning enables accurate and context-aware microbial source tracking and has the potential for novel microbial knowledge discovery.
The adapted Activity-By-Contact model for enhancer-gene assignment and its application to single-cell data Identifying regulatory regions in the genome is of great interest for understanding the epigenomic landscape in cells. One fundamental challenge in this context is to find the target genes whose expression the regulatory regions affect. A recent successful method is the Activity-By-Contact (ABC) model [10] which scores enhancer-gene interactions based on enhancer activity and the contact frequency of an enhancer to its target gene. It requires two types of assays to measure enhancer activity, which limits the applicability. The ABC score describes regulatory interactions entirely from a genes perspective, and does not account for all the candidate target genes of an enhancer. Moreover, there is no implementation available which would allow for an integration with transcription factor (TF) binding information or an efficient analysis of single-cell data. We show that the ABC-model is comparably accurate with only one assay to measure enhancer activity. Further, we demonstrate that it can yield a higher accuracy by adapting the enhancer activity according to the number of contacts the enhancer has to its candidate target genes. We combined our adapted ABC-model with TF binding information and illustrate an analysis of a single-cell ATAC-seq data set of the human heart. We were able to characterise cell type-specific regulatory interactions as well as to prioritise candidate TFs that drive cell type-specific expression. All executed processing steps are incorporated into our new computational pipeline STARE. The software is available at https://github.com/schulzlab/STARE.
METHYLOME DYNAMICS UPON PROTEASOME INHIBITION BY THE PSEUDOMONAS VIRULENCE FACTOR SYRINGOLIN A DNA methylation is an important epigenetic mark required for proper gene expression and silencing of transposable elements. DNA methylation patterns can be modified by environmental factors such as pathogen infection, where modification of DNA methylation can be associated with plant resistance. To counter the plant defense pathways, pathogens produce effectors molecule, several of which act as proteasome inhibitors. Here we investigated the effect of proteasome inhibition by the bacterial virulence factor Syringolin A on genome-wide DNA methylation. We show that Syringolin A treatment results in an increase of DNA methylation at centromeric and pericentromeric regions of Arabidopsis chromosomes. We identify several CHH DMRs that are enriched in the proximity of transcriptional start sites. Syringolin A treatment does not result in significant changes in small RNA composition. However, significant changes in genome transcriptional activity can be observed, including a strong upregulation of resistance genes that are located on chromosomal arms. We hypothesize that DNA methylation changes could be linked to the upregulation of some atypical members of the de novo DNA methylation pathway: AGO3, AGO9 and DRM1. Our data suggests that modification of genome-wide DNA methylation resulting from an inhibition of the proteasome by bacterial effectors could be part of a epi-genomic arms race against pathogens. The interaction of physical structure and nutrient loading drives ecosystem change in a large tropical lake over 40 years Many lakes across the world are entering novel states and experiencing altered biogeochemical cycling due to local anthropogenic stressors. In the tropics, understanding the drivers of these changes can be difficult due to a lack of documented historic conditions or an absence of continuous monitoring that can distinguish between intra- and inter-annual variation. Over the last forty years (1980-2020), Lake Yojoa (Honduras) has experienced increased watershed development as well as the introduction of a large net-pen Tilapia farm, resulting in a dramatic reduction in seasonal water clarity, increased trophic state and altered nutrient dynamics, shifting Lake Yojoa from an oligotrophic (low productivity) to mesotrophic (moderate productivity) ecosystem. To assess the changes that have occurred in Lake Yojoa as well as putative drivers for those changes, we compared Secchi depth (water clarity), dissolved inorganic nitrogen (DIN), and total phosphorus (TP) concentrations at continuous semi-monthly intervals for the three years between 1979-1983 and again at continuous 16-day intervals for 2018-2020. Between those two periods we observed the loss of a clear water phase that previously occurred in the months when the water column was fully mixed. Seasonal peaks in DIN coincident with mixing suggest that an enhanced accumulation of ammonium in the hypolimnion (the bottom layer of a stratified lake) during stratification, and release to the epilimnion (the top layer of a stratified lake) with mixing maintains high algal abundance and subsequently low secchi depth during what was previously the clear water phase. This interaction of nutrient loading and Lake Yojoas monomictic stratification regime illustrates a key phenomenon in how physical water column structure and nutrients interact in tropical monomictic lakes. This work highlights the need to consider nutrient dynamics of warm anoxic hypolimnions, not just surface water nutrient concentrations, to understand environmental change in these societally important but understudied ecosystems. β-catenin and canonical Wnts in Hydra pattern formation: Insights from mathematical modelling into two distinct systems Formation of the body axes and the subsequent formation of the apical termini are fundamental steps during animal development. Nuclear {beta}-catenin and canonical Wnt3 have been identified as major players in both these processes in Hydra. Based on molecular knowledge of canonical Wnt signalling directly linking nuclear {beta}-catenin and Wnt3 activity, de novo axis formation and head formation have frequently been assumed to be part of the same pattern formation system. In this study, we investigated this relationship by combining new model simulations with available experimental results. We demonstrated that nuclear {beta}-catenin and Wnt3 likely contribute to two distinct de novo pattern formation systems in Hydra, organizing development and differentiation on two different spatial scales. In particular, the results suggested that nuclear {beta}-catenin acts at the scale of the whole body to control axis formation, whereas canonical Wnt3 signalling is involved in a downstream pathway responsible for small-scale patterning of the head. These results suggest that the underlying mechanisms may also be more complex than previously assumed in other animals where axis formation has been ascribed to canonical Wnt signalling. Brain simulation augments machine-learning-based classification of dementia INTRODUCTIONComputational brain network modeling using The Virtual Brain (TVB) simulation platform acts synergistically with machine learning and multi-modal neuroimaging to reveal mechanisms and improve diagnostics in Alzheimers disease.

METHODSWe enhance large-scale whole-brain simulation in TVB with a cause-and-effect model linking local Amyloid {beta} PET with altered excitability. We use PET and MRI data from 33 participants of Alzheimers Disease Neuroimaging Initiative (ADNI3) combined with frequency compositions of TVB-simulated local field potentials (LFP) for machine-learning classification.

RESULTSThe combination of empirical neuroimaging features and simulated LFPs significantly outperformed the classification accuracy of empirical data alone by about 10% (weighted F1-score empirical 64.34% vs. combined 74.28%). Informative features showed high biological plausibility regarding the Alzheimers-typical spatial distribution.

DISCUSSIONThe cause-and-effect implementation of local hyperexcitation caused by Amyloid {beta} can improve the machine-learning-driven classification of Alzheimers and demonstrates TVBs ability to decode information in empirical data employing connectivity-based brain simulation.

RESEARCH IN CONTEXTO_LISYSTEMATIC REVIEW. Machine-learning has been proven to augment diagnostics of dementia in several ways. Imaging-based approaches enable early diagnostic predictions. However, individual projections of long-term outcome as well as differential diagnosis remain difficult, as the mechanisms behind the used classifying features often remain unclear. Mechanistic whole-brain models in synergy with powerful machine learning aim to close this gap.
C_LIO_LIINTERPRETATION. Our work demonstrates that multi-scale brain simulations considering Amyloid {beta} distributions and cause-and-effect regulatory cascades reveal hidden electrophysiological processes that are not readily accessible through measurements in humans. We demonstrate that these simulation-inferred features hold the potential to improve diagnostic classification of Alzheimers disease.
C_LIO_LIFUTURE DIRECTIONS. The simulation-based classification model needs to be tested for clinical usability in a larger cohort with an independent test set, either with another imaging database or a prospective study to assess its capability for long-term disease trajectories.
C_LI Cytosolic aggregation of mitochondrial proteins disrupts cellular homeostasis by stimulating the aggregation of other proteins Mitochondria are organelles with their own genomes, but they rely on the import of nuclear-encoded proteins that are translated by cytosolic ribosomes. Therefore, it is important to understand whether failures in the mitochondrial uptake of these nuclear-encoded proteins can cause proteotoxic stress and identify response mechanisms that may counteract it. Here, we report that upon impairments in mitochondrial protein import, high-risk precursor and immature forms of mitochondrial proteins form aberrant deposits in the cytosol. These deposits then cause further cytosolic accumulation and consequently aggregation of other mitochondrial proteins and disease-related proteins, including -synuclein and amyloid {beta}. This aggregation triggers a cytosolic protein homeostasis imbalance that is accompanied by specific molecular chaperone responses at both the transcriptomic and protein levels. Altogether, our results provide evidence that mitochondrial dysfunction, specifically protein import defects, contributes to impairments in protein homeostasis, thus revealing a possible molecular mechanism by which mitochondria are involved in neurodegenerative diseases. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. Disruption of the microphysiological niche alters matrix deposition and causes loss of mitochondrial Ca2+ control in skeletal muscle fibers Cells rapidly lose their physiological phenotype upon disruption of their extracellular matrix (ECM)-intracellular cytoskeleton interactions. Here, we investigated acute effects of ECM disruption on cellular and mitochondrial morphology, transcriptomic signatures, and Ca2+ handling in adult mouse skeletal muscle fibers. Adult skeletal muscle fibers were isolated from mouse toe muscle either by collagenase-induced dissociation of the ECM or by mechanical dissection that leaves the contiguous ECM intact. Experiments were generally performed four hours after cell isolation. At this time, there were striking differences in the gene expression patterns between fibers isolated with the two methods; 24h after cell isolation, enzymatically dissociated fibers had transcriptomic signatures resembling dystrophic phenotypes. Mitochondrial appearance was grossly similar in the two groups, but 3D electron microscopy revealed shorter and less branched mitochondria in enzymatically dissociated than in mechanically dissected fibers. Similar increases in free cytosolic [Ca2+] during repeated tetanic stimulation were accompanied by marked mitochondrial Ca2+ uptake only in enzymatically dissociated muscle fibers. The aberrant mitochondrial Ca2+ uptake was partially prevented by the mitochondrial Ca2+ uniporter inhibitor Ru360 and by cyclosporine A and NV556, which inhibit the mitochondrial protein Ppif (also called cyclophilin D). Importantly, inhibition of Ppif with NV556 significantly improved survival of mice with mitochondrial myopathy in which muscle mitochondria take up excessive amounts of Ca2+ even with an intact ECM. In conclusion, skeletal muscle fibers isolated by collagenase-induced dissociation of the ECM display aberrant mitochondrial Ca2+ uptake, which involves a Ppif-dependent mitochondrial Ca2+ influx resembling that observed in mitochondrial myopathies. Deep learning can predict multi-omic biomarkers from routine pathology images: A systematic large-scale study We assessed the pan-cancer predictability of multi-omic biomarkers from haematoxylin and eosin (H&E)-stained whole slide images (WSI) using deep learning (DL) throughout a systematic study. A total of 13,443 DL models predicting 4,481 multi-omic biomarkers across 32 cancer types were trained and validated. The investigated biomarkers included a broad range of genetic, transcriptomic, proteomic, and metabolic alterations. Furthermore, established markers relevant for prognosis, molecular subtypes and clinical outcomes were included. Overall, we established the general feasibility of predicting multi-omic markers directly from routine histology images with DL across solid cancer types, where 50% of the models could perform at an area under the curve (AUC) of more than 0.633 (with 25% of the models having an AUC larger than 0.711). A wide range of biomarkers were detectable from routine histology images across all investigated cancer types, with a mean AUC of at least 0.62 in almost all malignancies. Strikingly, we observed that biomarker predictability was mostly consistent and not dependent on sample size and class ratio, suggesting a degree of true predictability inherent in histomorphology. Together, the results of our study show the potential of DL to predict a multitude of biomarkers across the omics spectrum using only routine slides. This paves the way for accelerating diagnosis and developing more precise treatments for cancer patients. An hidden Markov model to estimate homozygous-by-descent probabilities associated with nested layers of ancestors Inbreeding results from the mating of related individuals and has negative consequences because it brings together deleterious variants in one individual. Genomic estimates of the inbreeding coefficients are preferred to pedigree-based estimators as they measure the realized inbreeding levels and they are more robust to pedigree errors. Several methods identifying homozygous-by-descent (HBD) segments with hidden Markov models (HMM) have been recently developed and are particularly valuable when the information is degraded or heterogeneous (e.g., low-fold sequencing, low marker density, heterogeneous genotype quality or variable marker spacing). We previously developed a multiple HBD class HMM where HBD segments are classified in different groups based on their length (e.g., recent versus old HBD segments) but we recently observed that for high inbreeding levels with many HBD segments, the estimated contributions might be biased towards more recent classes (i.e., associated with large HBD segments) although the overall estimated level of inbreeding remained unbiased. We herein propose a new model in which the HBD classification is modeled in successive nested levels with decreasing expected HBD segment lengths, the underlying exponential rates being directly related to the number of generations to the common ancestor. The non-HBD classes are now modeled as a mixture of HBD segments from later generations and shorter non-HBD segments (i.e., both with higher rates). The new model has improved statistical properties and performs better on simulated data compared to our previous version. We also show that the parameters of the model are easier to interpret and that the model is more robust to the choice of the number of classes. Overall, the new model results in an improved partitioning of inbreeding in different HBD classes and should be preferred. In-depth characterization of denitrifier communities across different soil ecosystems in the tundra BackgroundIn contrast to earlier assumptions, there is now mounting evidence for the role of tundra soils as important sources of the greenhouse gas nitrous oxide (N2O). However, the microorganisms involved in the cycling of N2O in this system remain largely uncharacterized. Since tundra soils are variable sources and sinks of N2O, we aimed at investigating differences in community structure across different soil ecosystems in the tundra.

ResultsWe analysed 1.4 Tb of metagenomic data from soils in northern Finland covering a range of ecosystems from dry upland soils to water-logged fens and obtained 796 manually binned and curated metagenome-assembled genomes (MAGs). We then searched for MAGs harbouring genes involved in denitrification, an important process driving N2O emissions. Communities of potential denitrifiers were dominated by microorganisms with truncated denitrification pathways (i.e., lacking one or more denitrification genes) and differed across soil ecosystems. Upland soils showed a strong N2O sink potential and were dominated by members of the Alphaproteobacteria such as Bradyrhizobium and Reyranella. Fens, which had in general net-zero N2O fluxes, had a high abundance of poorly characterized taxa affiliated with the Chloroflexota lineage Ellin6529 and the Acidobacteriota subdivision Gp23.

ConclusionsBy coupling an in-depth characterization of microbial communities with in situ measurements of N2O fluxes, our results suggest that the observed spatial patterns of N2O fluxes in the tundra are related to differences in the composition of denitrifier communities. Failure to account for behavioral variability significantly compromises accuracy in indirect population monitoring Indirect wildlife population surveying largely depends upon counts of artefacts of behavior (e.g., nests or dung). Likelihood to encounter these artefacts is derived from both artefact production and decay, and variability in artefact production behavior is considered to contribute minimally to inaccuracy in wildlife estimation. Here we demonstrate how ignoring behavioral variability contributes to significant population misestimation, using an example of an endangered ape, the bonobo (Pan paniscus). Until now, a single estimate of (i.e., nest construction) has been used to extrapolate bonobo densities, assumed to be broadly representative of bonobo sign production behavior. We estimated nest construction rate across seasons and social groups at the Kokolopori Bonobo Reserve, DRC and find that nest construction rates in bonobos to be highly variable across populations as well as seasonal. Failure to account for this variability led to potentially severe degradation in accuracy of bonobo population estimates of abundance, accounting for a likely overestimation of bonobo numbers by 34%, and at worst as high as 80%. With this example, we demonstrate that failure to account for inter- and intra-population behavioral variation compromises the ability to estimate both relative and absolute wildlife abundances. We argue that variation in sign production is but one of several potential ways that behavioral variability can affect conservation monitoring, should be measured across contexts whenever possible, and must be considered in population estimation confidence intervals. With increasing attention to behavioral variability as a potential tool for conservation, conservationists must also account for the impact that behavioral variability can play upon wildlife population estimation. Our results underline the importance of observational research to wildlife monitoring schemes as a critical but overlooked component of conservation management. We extensively discuss the various avenues through which behavioral variability is likely to impact wildlife monitoring accuracy and precision and propose potential avenues for accounting for behavioral variability in wildlife monitoring. An alignment free approach confirms semantic properties of species proteomes Alignment-based methods dominate molecular biology. However, by primarily allowing one-to-one comparisons, these methods are focused on a gene-centered viewpoint and lack the broad context essential to analyze how complex biological systems function and evolve. In actuality, a gene is part of genome where more than one sequence contributes to the functional network and evolutionary trajectory of the cell. The need for conservation of established interactions, is arguably more important to the evolutionary success of species than conservation of individual function. To test whether such contextual information exists, a distributional semantics method - Latent Semantic Analysis (LSA), was applied to thousands of species proteomes. Using natural language processing, Latent Taxonomic Signatures (LTSs) were identified that outperformed existing alignment-based BLAST methods when random protein sequences were being mapped to annotated taxonomy according to GenBank. LTSs are a novel proteome distributed feature, suggesting the existence of evolutionary constraints imposed on individual proteins by their proteome context. Even orphan proteins are exhibiting LTSs, which makes their uniqueness linked to a specific taxonomic level questionable. Unlike more simple bias, LTSs represent a self-similarity pattern, where random sets of species proteins show the same statistical properties of a complete proteome at many scales. Natural language processing and machine learning provide insights not easily discernable using alignment based methods suggestive there is more to species related differences than just translational optimization. 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. Extracellular Vesicle (EV) Mitochondria and EV/HSP27 Mixtures Protect the Ischemic Brain Endothelium Ischemic stroke causes brain endothelial cell (BEC) death and damages tight junction integrity of the blood-brain barrier (BBB). We harnessed the innate mitochondrial load of endothelial cell-derived extracellular vesicles (EVs) and utilized mixtures of EV/exogenous heat shock protein 27 (HSP27) as a one-two punch strategy to increase BEC survival (via EV mitochondria) and preserve their tight junction integrity (via HSP27 effects). We demonstrated that the medium-to-large (m/lEV) but not small EVs (sEV) transferred their mitochondrial load, which subsequently colocalized with the mitochondrial network of the recipient primary human BECs. BECs treated with m/lEVs increased relative ATP levels and displayed superior mitochondrial function. m/lEVs isolated from rotenone (mitochondrial complex I inhibitor) exposed BECs (RTN-m/lEVs) did not increase BECs ATP levels compared to naive m/lEVs. In contrast, RTN-sEVs functionality was minimally affected compared to naive sEVs. Compared to EV/HSP27 mixtures, PEG-DET/EV/HSP27 ternary mixtures increased HSP27 protein loading. (PEG-DET/HSP27)/EV ternary and EV/HSP27 binary mixtures decreased the paracellular permeability of small and large molecular mass fluorescent tracers in oxygen glucose-deprived primary human BECs. This one-two punch approach to increase BEC metabolic function and tight junction integrity is a promising strategy for BBB protection and prevention of long-term neurological dysfunction post-ischemic stroke.

Highlights[tpltrtarr] Medium-to-large extracellular vesicles (m/lEVs), not small EVs contain mitochondria
[tpltrtarr]m/lEVs increased ATP and mitochondrial function in brain endothelial cells (BECs)
[tpltrtarr]Donor m/lEV mitochondria colocalized with host mitochondria in recipient BECs
[tpltrtarr]m/lEVs from rotenone-exposed BECs did not increase recipient BEC ATP levels
[tpltrtarr]PEG-DET/EV/HSP27 ternary mixtures increase exogenous HSP27 protein loading


Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=133 SRC="FIGDIR/small/466491v4_ufig1.gif" ALT="Figure 1">
View larger version (43K):
org.highwire.dtl.DTLVardef@155f618org.highwire.dtl.DTLVardef@8fa778org.highwire.dtl.DTLVardef@13bd44aorg.highwire.dtl.DTLVardef@ee771e_HPS_FORMAT_FIGEXP  M_FIG C_FIG A distinct brain expression profile for genes within large introgression deserts and under positive selection in Homo sapiens Analyses of ancient DNA from extinct hominins have provided unique insights into the complex evolutionary history of Homo sapiens, intricately related to that of the Neanderthals and the Denisovans as revealed by several instances of admixture events. These analyses have also allowed the identification of introgression deserts: genomic regions in our species that are depleted of  archaic haplotypes. The presence of genes like FOXP2 in these deserts has been taken to be suggestive of brain-related functional differences between Homo species. Here, we seek a deeper characterization of these regions, taking into account signals of positive selection in our lineage. Analyzing publicly available transcriptomic data from the human brain at different developmental stages, we found that structures outside the cerebral neocortex, in particular the cerebellum, the striatum and the mediodorsal nucleus of the thalamus show the most divergent transcriptomic profiles when considering genes within large introgression deserts and under positive selection. PySTACHIO: Python Single-molecule TrAcking stoiCHiometry Intensity and simulatiOn, a flexible, extensible, beginner-friendly and optimized program for analysis of single-molecule microscopy As camera pixel arrays have grown larger and faster, and optical microscopy techniques ever more refined, there has been an explosion in the quantity of data acquired during routine light microcopy. At the single-molecule level, analysis involves multiple steps and can rapidly become computationally expensive, in some cases intractable on office workstations. Complex bespoke software can present high activation barriers to entry for new users. Here, we redevelop our quantitative single-molecule analysis routines into an optimized and extensible Python program, with GUI and command-line implementations to facilitate use on local machines and remote clusters, by beginners and advanced users alike. We demonstrate that its performance is on par with previous MATLAB implementations but runs an order of magnitude faster. We tested it against challenge data and demonstrate its performance is comparable to state-of-the-art analysis platforms. We show the code can extract fluorescence intensity values for single reporter dye molecules and, using these, estimate molecular stoichiometries and cellular copy numbers of fluorescently-labeled biomolecules. It can evaluate 2D diffusion coefficients for the characteristically short single-particle tracking data. To facilitate benchmarking we include data simulation routines to compare different analysis programs. Finally, we show that it works with 2-color data and enables colocalization analysis based on overlap integration, to infer interactions between differently labelled biomolecules. By making this freely available we aim to make complex light microscopy single-molecule analysis more democratized.

HighlightsO_LIWe present PySTACHIO, a refined version of our spot tracking algorithm
C_LIO_LIWe demonstrate highly improved performance over previous MATLAB versions
C_LIO_LIPySTACHIO can accurately estimate stoichiometries and 2D diffusion coefficients
C_LIO_LIPerformance is comparable to state-of-the-art packages on challenge data
C_LIO_LIPySTACHIO has both GUI and command line interfaces and can be hosted as a web app
C_LI

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=77 SRC="FIGDIR/small/435952v3_ufig1.gif" ALT="Figure 1">
View larger version (27K):
org.highwire.dtl.DTLVardef@17cc45dorg.highwire.dtl.DTLVardef@3fb66eorg.highwire.dtl.DTLVardef@e20360org.highwire.dtl.DTLVardef@4ba47d_HPS_FORMAT_FIGEXP  M_FIG C_FIG The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Measuring transcription factor binding and gene expression using barcoded self-reporting transposon calling cards and transcriptomes Calling cards technology using self-reporting transposons enables the identification of DNA-protein interactions through RNA sequencing. Although immensely powerful, current implementations of calling cards in bulk experiments on populations of cells are technically cumbersome and require many replicates to identify independent insertions into the same genomic locus. Here, we have drastically reduced the cost and labor requirements of calling card experiments in bulk populations of cells by introducing a DNA barcode into the calling card itself. An additional barcode incorporated during reverse transcription enables simultaneous transcriptome measurement in a facile and affordable protocol. We demonstrate that barcoded self-reporting transposons recover in vitro binding sites for four basic helix-loop-helix transcription factors with important roles in cell fate specification: ASCL1, MYOD1, NEUROD2, and NGN1. Further, simultaneous calling cards and transcriptional profiling during transcription factor overexpression identified both binding sites and gene expression changes for two of these factors. Lastly, we demonstrated barcoded calling cards can record binding in vivo in the mouse brain. In sum, RNA-based identification of transcription factor binding sites and gene expression through barcoded self-reporting transposon calling cards and transcriptomes is an efficient and powerful method to infer gene regulatory networks in a population of cells. Nuclear Progestin Receptor Mediated Linkage of Blood Coagulation and Ovulation Ovulation is a dramatic remodeling process that includes rupture of blood capillaries and clotting, but coagulation is not thought to directly regulate this process. Herein, we report remarkable increases of coagulation factors V (f5, [~]3145-fold) and tissue factor (f3a, [~]120-fold) in zebrafish ovarian follicle cells during ovulation. This increase was mediated through the nuclear progestin receptor (Pgr), which is essential for ovulation in zebrafish, and was totally abolished in ovarian follicular cells from pgr-/- mutants. In addition, promoter activities of f5 and f3a were significantly enhanced by progestin (DHP) via Pgr. Similar regulation of human F5 promoter activity was induced via human PGRB, suggesting a conserved mechanism. Site-directed mutagenesis of the zebrafish f5 promoter further demonstrated a direct regulation of coagulation factors via progestin response elements. Moreover, a stark increase of erythrocytes occurred in capillaries meshed in wildtype preovulatory follicles but was absent in pgr-/- mutants. Interestingly, anticoagulants significantly inhibited ovulation both in vitro and in vivo, respectively. Furthermore, reduced fecundity was observed in f5+/- female zebrafish. Taken together, our study provides plausible evidence for steroid regulation of coagulation factors, and a new hypothesis for blood clotting triggered ovulation in vertebrates. 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. Targeted isolation of panels of diverse human protective broadly neutralizing antibodies against SARS-like viruses The emergence of current SARS-CoV-2 variants of concern (VOCs) and potential future spillovers of SARS-like coronaviruses into humans pose a major threat to human health and the global economy 1-7. Development of broadly effective coronavirus vaccines that can mitigate these threats is needed 8, 9. Notably, several recent studies have revealed that vaccination of recovered COVID-19 donors results in enhanced nAb responses compared to SARS-CoV-2 infection or vaccination alone 10-13. Here, we utilized a targeted donor selection strategy to isolate a large panel of broadly neutralizing antibodies (bnAbs) to sarbecoviruses from two such donors. Many of the bnAbs are remarkably effective in neutralization against sarbecoviruses that use ACE2 for viral entry and a substantial fraction also show notable binding to non-ACE2-using sarbecoviruses. The bnAbs are equally effective against most SARS-CoV-2 VOCs and many neutralize the Omicron variant. Neutralization breadth is achieved by bnAb binding to epitopes on a relatively conserved face of the receptor binding domain (RBD) as opposed to strain-specific nAbs to the receptor binding site that are commonly elicited in SARS-CoV-2 infection and vaccination 14-18. Consistent with targeting of conserved sites, select RBD bnAbs exhibited in vivo protective efficacy against diverse SARS-like coronaviruses in a prophylaxis challenge model. The generation of a large panel of potent bnAbs provides new opportunities and choices for next-generation antibody prophylactic and therapeutic applications and, importantly, provides a molecular basis for effective design of pan-sarbecovirus vaccines. Induction of exaggerated cytokine production in human peripheral blood mononuclear cells by a recombinant SARS-CoV-2 spike glycoprotein S1 and its inhibition by dexamethasone An understanding of the pathological inflammatory mechanisms involved in SARS-CoV-2 virus infection is necessary in order to discover new molecular pharmacological targets for SARS-CoV-2 cytokine storm. In this study, the effects of a recombinant SARS-CoV-2 spike glycoprotein S1 was investigated in human peripheral blood mononuclear cells (PBMCs). Stimulation of PBMCs with spike glycoprotein S1 (100 ng/mL) resulted in significant elevation in the production of TNF, IL-6, IL-1{beta} and IL-8. However, pre-treatment with dexamethasone (100 nM) caused significant reduction in the release of these cytokines. Further experiments revealed that S1 stimulation of PBMCs increased phosphorylation of NF-{kappa}B p65 and I{kappa}B, and I{kappa}B degradation. DNA binding of NF-{kappa}B p65 was also significantly increased following stimulation with spike glycoprotein S1. Treatment of PBMCs with dexamethasone (100 nM) or BAY11-7082 (1 {micro}M) resulted in inhibition of spike glycoprotein S1-induced NF-{kappa}B activation. Activation of p38 MAPK by S1 was blocked in the presence of dexamethasone and SKF 86002. CRID3, but not dexamethasone pre-treatment produced significant inhibition of S1-induced activation of NLRP3/caspase-1. Further experiments revealed that S1-induced increase in the production of TNF, IL-6, IL-1{beta} and IL-8 was reduced in the presence of BAY11-7082 and SKF 86002, while CRID3 pre-treatment resulted in the reduction of IL-1{beta} production. These results suggest that SARS-CoV-2 spike glycoprotein S1 stimulated PBMCs to release pro-inflammatory cytokines through mechanisms involving activation of NF-{kappa}B, p38 MAPK and NLRP3 inflammasome. It is proposed that the clinical benefits of dexamethasone in COVID-19 is possibly due to its anti-inflammatory activity in reducing SARS-CoV-2 cytokine storm. Lipid membranes modulate the activity of RNA through sequence-dependent interactions RNA is a ubiquitous biomolecule that can serve as both catalyst and information carrier. Understanding how RNA bioactivity is controlled is crucial for elucidating its physiological roles and potential applications in synthetic biology. Here we show that lipid membranes can act as RNA organization platforms, introducing a novel mechanism for ribo-regulation. The activity of R3C ribozyme can be modified by the presence of lipid membranes, with direct RNA-lipid interactions dependent on RNA nucleotide content, base pairing and length. In particular, the presence of guanine in short RNAs is crucial for RNA-lipid interactions, and G-quadruplex formation further promotes lipid binding. Lastly, by artificially modifying the R3C substrate sequence to enhance membrane binding we generated a lipid-sensitive ribozyme reaction with riboswitch-like behavior. These findings introduce RNA-lipid interactions as a tool for developing synthetic riboswitches and novel RNA-based lipid biosensors, and bear significant implications for RNA World scenarios for the origin of life.