Optogenetic activation of the neural circuit between dorsal raphe nucleus and the pre-Botzinger complex contributes to inhibition of the seizure-induced respiratory arrest in the DBA/1 mouse SUDEP model Sudden unexpected death in epilepsy (SUDEP) is the leading cause of death among epilepsy patients, occurring even more frequently in cases with anti-epileptic drug resistance. However, the underlying mechanism of SUDEP remains elusive. Our previous study demonstrated that enhancement of serotonin (5-HT) synthesis by intraperitoneal (IP) injection of 5-hydroxytryptophan (5-HTP) significantly reduced the incidence of seizure-induced respiratory arrest (S-IRA) in a DBA/1 mouse SUDEP model. Given that the 5-HT2A receptor (5-HT2AR) plays an important role in mediating the respiration system in the brain, we hypothesized that 5-HT2AR plays a key role in S-IRA and SUDEP. To test this hypothesis, we examined whether the decreased incidence of S-IRA evoked by either acoustic stimulation or pentylenetetrazole (PTZ) injection following 5-HTP administration will be blocked by treatment with ketanserin (KET), a selective antagonist of 5HT2AR, in the DBA/1 mouse SUDEP model. We observed that the reduction in S-IRA by 5-HTP was significantly reversed by IP or intracerebroventricular injection of KET. Considering the localization of 5-HT2AR in the pre-Botzinger complex (PBC), which plays a key role in regulating respiratory rhythm, we next examined whether KET acts on 5-HT2AR in the PBC. To test this hypothesis, we activated the neural circuit between the dorsal raphe nucleus (DR) and PBC using optogenetics technology. We observed that stimulation of TPH2-ChETA-expressing neurons in the DR reduced the incidence of S-IRA evoked by PTZ, and this suppressant effect was significantly reversed by administration of KET in the bilateral PBC with no changes in electroencephalogram activity. The neural circuit between the DR and PBC was confirmed by injection of cholera toxin subunit B555 (CTB-555), a nerve tracer, in the DR or PBC separately. Calcium signaling evoked by PTZ within neurons of the PBC during seizures was significantly reduced by photostimulation of the DR. Taken together, our findings suggest that 5-HT2AR plays a critical role in regulating S-IRA and targeting the serotonergic neural circuit between the DR and PBC is a promising approach to preventing SUDEP.
Pax3 repairs a neural circuit through a program of directed axon outgrowth Repairing damaged or dysfunctional human brain circuits remains an ongoing challenge for biomedical science. While surviving neuronal networks can be reorganised after lesion, for example by neurotrophins, these new connections are disorganised and rarely produce clinical improvement. Here we investigate how to promote axonal growth while retaining correct cellular targeting. We show that, in response to brain-derived neurotrophic factor (BDNF) in target-tissue, potential reinnervating neurons upregulate Pax3. Pax3 in turn increases polysialic acid-neural cell adhesion molecule (PSA-NCAM) on their axon terminals, facilitating their outgrowth and pathfinding, and resulting in correctly-targeted neural circuit repair in the mature nervous system. This is a novel role for Pax3, which we confirmed by showing its expression in afferent neurons is essential for spontaneous and BDNF-induced reinnervation in the developing and mature brains, respectively. Together these results suggest that Pax3 contributes to a repair program, in which axon growth is promoted and direction signaling maintained. These data advance our understanding for accurately rebuilding neural circuits: restricting growth-promotion to potential afferent neurons, as opposed to stimulating the whole circuit, allows axon growth without impairing its guidance.
Experimental evolution in TP53 deficient human gastric organoids recapitulates tumorigenesis The earliest events during human tumor initiation are poorly characterized but may hold clues as to how to detect and prevent malignancy. Here we model occult pre-neoplasia by engineering TP53 deficiency, a common and early event in gastric cancer, into human gastric organoids. By performing experimental evolution in multiple clonally derived cultures over two years we define causal relationships between this common initiating genetic lesion and resulting phenotypes. TP53 loss elicited progressive aneuploidy, including copy number alterations and complex structural variants that are common in gastric cancers and which follow preferred temporal orders. Longitudinal single cell sequencing of TP53 deficient gastric organoids similarly indicates progression towards malignant transcriptional programs. Moreover, lineage tracing with expressed cellular barcodes demonstrates reproducible dynamics whereby initially rare subclones with shared transcriptional programs repeatedly attain clonal dominance. This powerful platform for experimental evolution exposes stringent selection, clonal interference and a striking degree of phenotypic convergence in pre-malignant epithelial organoids, implying that the earliest stages of tumorigenesis may be predictable while illuminating evolutionary constraints and barriers to malignant transformation.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/487529v2_ufig1.gif" ALT="Figure 1">
View larger version (63K):
org.highwire.dtl.DTLVardef@122debborg.highwire.dtl.DTLVardef@4cd2daorg.highwire.dtl.DTLVardef@1a720d6org.highwire.dtl.DTLVardef@15ccaa7_HPS_FORMAT_FIGEXP  M_FIG C_FIG 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. Spatial organization of chromosomes leads to heterogeneous chromatin motion and drives the liquid- or gel-like dynamical behavior of chromatin Chromosome organization and dynamics are involved in regulating many fundamental processes such as gene transcription and DNA repair. Experiments unveiled that chromatin motion is highly heterogeneous inside cell nuclei, ranging from a liquid-like, mobile state to a gel-like, rigid regime. Using polymer modeling, we investigate how these different physical states and dynamical heterogeneities may emerge from the same structural mechanisms. We found that the formation of topologically-associating domains (TADs) is a key driver of chromatin motion heterogeneity. In particular, we demonstrated that the local degree of compaction of the TAD regulates the transition from a weakly compact, fluid state of chromatin to a more compact, gel state exhibiting anomalous diffusion and coherent motion. Our work provides a comprehensive study of chromosome dynamics and a unified view of chromatin motion enabling to interpret the wide variety of dynamical behaviors observed experimentally across different biological conditions, suggesting that the  liquid or  solid behaviour of chromatin are in fact two sides of the same coin. PIN-FORMED1 polarity in the shoot is insensitive to the polarity of neighbouring cells In plants, epidermal cells are planar-polarized along an axis marked by the asymmetric localization patterns of several proteins including PIN-FORMED1 (PIN1), which facilitates the directional efflux of the plant hormone auxin to pattern phyllotaxis (Heisler et al., 2010; Mansfield et al., 2018; Reinhardt et al., 2003). While PIN1 polarity is known to be regulated non-cell autonomously via the MONOPTEROS (MP) transcription factor, how this occurs has not been determined (Bhatia et al., 2016). Here we use mosaic expression of the serine threonine kinase PINOID (PID) to test whether PIN1 polarizes according to the polarity of neighbouring cells. Our findings reveal that PIN1 is insensitive to the polarity of PIN1 in neighbouring cells arguing against auxin flux or extracellular auxin concentrations acting as a polarity cue, in contrast to previous model proposals (Abley et al., 2016; Stoma et al., 2008). Heteologous saRNA-Prime, DNA Dual-Antigen-Boost SARS-CoV-2 Vaccination Elicits Robust Cellular Immunogenicity and Cross-Variant Neutralizing Antibodies We assessed if immune responses are enhanced in CD-1 mice by heterologous vaccination with two different nucleic acid-based COVID-19 vaccines: a next-generation human adenovirus serotype 5 (hAd5)-vectored dual-antigen spike (S) and nucleocapsid (N) vaccine (AdS+N) and a self-amplifying and -adjuvanted S RNA vaccine (SASA S) delivered by a nano-lipid carrier. The AdS+N vaccine encodes S modified with a fusion motif to increase cell-surface expression. The N antigen is modified with an Enhanced T-cell Stimulation Domain (N-ETSD) to direct N to the endosomal/lysosomal compartment and increase MHC class I and II stimulation potential. The S sequence in the SASA S vaccine comprises the D614G mutation, two prolines to stabilize S in the prefusion conformation, and 3 glutamines in the furin cleavage region to increase cross-reactivity across variants. CD-1 mice received vaccination by homologous and heterologous prime > boost combinations. Humoral responses to S were the highest with any regimen including the SASA S vaccine, and IgG bound to wild type and Delta (B.1.617.2) variant S1 at similar levels. An AdS+N boost of an SASA S prime particularly enhanced both CD4+ and CD8+ T-cell responses to both wild type and Delta S peptides relative to all other vaccine regimens. Sera from mice receiving SASA S homologous or heterologous vaccination were found to be highly neutralizing of all pseudovirus strains tested: Wuhan, Beta, Delta, and Omicron strain. The findings here support the clinical testing of heterologous vaccination by an SASA S > AdS+N regimen to provide increased protection against emerging SARS-CoV-2 variants. β-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. Analysis of Drosophila cardiac hypertrophy by micro-computerized tomography for genetic dissection of heart growth mechanisms Heart failure is often preceded by pathological cardiac hypertrophy, a thickening of the heart musculature driven by complex gene regulatory and signaling processes. The Drosophila heart has great potential as a genetic model for deciphering the underlying mechanisms of cardiac hypertrophy. However, current methods for evaluating hypertrophy of the Drosophila heart are laborious and difficult to carry out reproducibly. Here we demonstrate that micro-computerized tomography (microCT) is an accessible, highly reproducible method for non-destructive, quantitative analysis of Drosophila heart morphology and size. To validate our microCT approach for analyzing Drosophila cardiac hypertrophy, we show that expression of constitutively active Ras (Ras85DV12), previously shown to cause hypertrophy of the fly heart, results in significant thickening of both adult and larval heart walls when measured from microCT images. We then show using microCT analysis that genetic upregulation of store-operated Ca2+ entry (SOCE) driven by expression of constitutively active Stim (StimCA) or Orai (OraiCA) proteins also results in significant hypertrophy of the Drosophila heart, through a process that specifically depends on Orai Ca2+ influx channels. Intravital imaging of heart contractility revealed significantly reduced end diastolic dimensions in StimCA and OraiCA expressing hearts, consistent with the hypertrophic phenotype. These results demonstrate that increased SOCE activity is an important driver of hypertrophic cardiomyocyte growth, and demonstrate how microCT analysis combined with tractable genetic tools in Drosophila can be used to delineate molecular signaling processes that underlie cardiac hypertrophy and heart failure.

NEW AND NOTEWORTHYGenetic analysis of cardiac hypertrophy in Drosophila holds immense potential for the discovery of new therapeutic targets to prevent and treat heart failure. However, this potential has been hindered by a lack of rapid and effective methods for analysis of heart size in flies. Here we demonstrate that analysis of the Drosophila heart with micro-computerized tomography yields accurate and highly reproducible heart size measurements that can be used to efficiently analyze heart growth and cardiac hypertrophy in Drosophila. A new method for the sampling and preservation of placental specimens in low-resource settings for the identification of P. falciparum and analysis of nucleic acids Collection, preservation, and shipment of histological specimens in low-resource settings is challenging. We present a novel method that achieved excellent preservation of placental specimens from rural Mali by using formalin fixation, ethanol dehydration, and long-term storage in a solar-powered freezer. Sample preservation success was 92%, permitting evaluation of current and past malaria infection, anemia, placental maturity, and inflammation. Using RNAscope(R) hybridization we were able to visualize cell-specific gene expression patterns in the formalin-fixed paraffin-embedded (FFPE) specimens. Additionally, our method entailed mirrored sampling from the two cut faces of a cotyledon, one for the FFPE workflows and the other for storage in RNAlater and RNA-seq. The mechanism of full activation of tumor suppressor PTEN at the phosphoinositide-enriched membrane Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP3 produced by PI3Ks. Excess PIP3 promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP2 and PIP3. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP3. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP3 allosterically unfolding the N-terminal PIP2 binding domain, positioning it favorably for the polybasic motif interaction with PIP2. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer. Exploring functional protein covariation across single cells using nPOP Many biological processes, such as the cell division cycle, are reflected in protein covariation across single cells. This covariation can be quantified and interpreted by single-cell proteomics with sufficiently high throughput and accuracy. Toward this goal, we developed the nano-ProteOmic sample Preparation (nPOP) method for single-cell proteomics. nPOP uses piezo acoustic dispensing to isolate individual cells in 300 picoliter volumes and performs all subsequent preparation steps in small droplets on a fluorocarbon-coated slide. This design enables simultaneous sample preparation of thousands of single cells, including lysing, digesting, and labeling individual cells in volumes below 20 nl. We used nPOP to prepare 1,888 single cells and 128 negative controls in a single batch. Their analysis enabled quantifying the covariation between thousands of proteins and cell-cycle protein markers. Many protein sets covaried with the cell cycle similarly across all cell types and states, reflecting cell-type independent cell cycle functions. However, the cell cycle covariation of other protein sets differed markedly between cell types, even within subpopulation of melanoma cells expressing markers for drug-resistance priming. The cells expressing these markers accumulated in the G1 phase of the cell cycle and exhibited different covariation of enzymes catabolizing glucose. These results demonstrate that protein covariation across single cells may reveal functionally concerted biological differences between closely related cell states.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=59 SRC="FIGDIR/small/441211v3_ufig1.gif" ALT="Figure 1">
View larger version (53K):
org.highwire.dtl.DTLVardef@1e23b1org.highwire.dtl.DTLVardef@b93faeorg.highwire.dtl.DTLVardef@3c00forg.highwire.dtl.DTLVardef@362b41_HPS_FORMAT_FIGEXP  M_FIG C_FIG Molecular diversity maintained by long-term balancing selection in mating loci defines multiple mating types in fungi Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of specimens. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci defines mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty specimens of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. Our sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Our analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. The inferred allelic information mirrored the outcome of in vitro crosses, thus allowing us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. The observed allelic classes could potentially generate 14,560 different mating types. We conclude that the genetic diversity of mating in Trichaptum loci is due to long-term balancing selection, with limited recombination and duplication activity. Our large number of sequenced specimens highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Author summaryFungi have complex mating systems, and basidiomycete fungi can encode thousands of mating types. Individuals with the same mating type cannot mate. This sexual system has evolved to facilitate sexual mating, increasing the chances to recombine into advantageous allelic combination and prune deleterious alleles. We explored the genomes of hundred and eighty specimens, combined with experimental mating studies of selected specimens, from a non-model organism (Trichaptum). We characterized the genomic regions controlling sex. The mating ability of the specimens confirmed the role of the mating alleles observed in the genomic data. The detailed analyses of many specimens allowed us to observe gene duplication and rearrangements within the mating loci, increasing the diversity within these loci. We supported previous suggestions of balancing selection in this region, an evolutionary force that maintains genomic diversity. These results supports that our fungal specimens are prone to outcross, which might facilitate the adaptation to new conditions. Conserved organ-specific microbial assemblages in different populations of a terrestrial crab Brachyuran crabs originated in the oceans and evolved specific morphological and physiological adaptations to live in freshwater, intertidal and even terrestrial habitats but the role of a selection mechanism involving symbiotic microorganisms long these colonization processes are not known. In this work we investigated the associated microbiota of three populations of a terrestrial brachyuran crab, Chiromantes haematocheir, to find evidence of a conserved crab-specific microbiome unrelated to the population of origin and dissimilar from environmental microbial assemblages. Bacterial 16S rRNA gene and fungal ITS sequences were obtained from selected crab organs and environmental matrices to profile microbial communities. In spite of the presence of truly marine larval stages and no gregarious behaviour, favouring microbiota exchanges, we found common, organ-specific microbiota, associated to the gut and the gills of the crabs (with more than 15% of the genera detected specifically enriched only in one organ). Our results suggest an early establishment of a new common, stable microbiota in the transition from water to land. Comprehensive O-GlcNAc glycoproteomics on NOTCH1 EGF repeats implicated unique Lewis X epitopes in mammals The O-GlcNAc modification of Notch receptors regulates Notch ligand interactions in a manner distinct from other forms of O-glycans on epidermal growth factor-like (EGF) repeats of Notch receptors. Although many proteins, besides Notch receptors, are expected to be O-GlcNAcylated by EGF domain-specific O-GlcNAc transferase (EOGT), only a small number of proteins have been reported to be modified in vivo, and elongated O-GlcNAc glycans have not been extensively explored. To extend our view of the specificity and variety of the glycan modification, we conducted a comprehensive analysis of O-GlcNAc glycans on NOTCH1 in mammals. Mass spectrometric analysis of NOTCH1 fragments expressed in HEK293T cells revealed that several EGF domains with putative O-GlcNAcylation sites were hardly modified with O-GlcNAc. Although amino acid residues before the modification site are preferentially occupied with aromatic residues, Phe and Tyr are preferable to Trp for the apparent modification with O-GlcNAc. Furthermore, a minor form of fucosylated O-GlcNAc glycans was detected in a subset of EGF domains. Fucosylation of O-GlcNAc glycans was enhanced by FUT1, FUT2, or FUT9 expression. The FUT9-dependent Lewis X epitope was confirmed by immunoblotting using an anti-Lewis X antibody. As expected from the similarity in the glycan structures, the Lexis X antigen was detected on O-fucose glycans. Our results refined the putative consensus sequence for the EOGT-dependent extracellular O-GlcNAc modification in mammals and revealed the structural diversity of functional Notch O-glycans. Adult fibroblasts use aggresomes only in distinct cell-states The aggresome is a protein turnover system in which proteins are trafficked along microtubules to the centrosome for degradation. Despite extensive focus on aggresomes in immortalized cell lines, it remains unclear if the aggresome is conserved in all primary cells and all cell-states. Here we examined the aggresome in primary adult mouse dermal fibroblasts shifted into four distinct cell-states. We found that in response to proteasome inhibition, quiescent and immortalized fibroblasts formed aggresomes whereas proliferating and senescent fibroblasts did not. Using this model, we generated a resource to better understand the aggresome, providing a characterization of the proteostasis networks in which the aggresome is used and transcriptomic features associated with the presence or absence of aggresome formation. Using this resource, we validate a previously reported role for p38 MAPK signaling in aggresome formation and identify TAK1 as a novel driver of aggresome formation upstream of p38 MAPKs. Together, our data demonstrate that the aggresome is a non-universal protein degradation system which can be used cell-state specifically and provide a resource for better understanding aggresome formation and function. Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: the Multi-Targeting Drug DREAM Challenge A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ( polypharmacology). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Author SummaryMany modern drugs are developed with the goal of modulating a single cellular pathway or target. However, many drugs are, in fact,  dirty; they target multiple cellular pathways or targets. This phenomenon is known as multi-targeting or polypharmacology. While some strive to develop  cleaner therapeutics that eliminate secondary targets, recent work has shown that multi-targeting therapeutics have key advantages for a variety of diseases. However, while multi-targeting drugs that affect a precisely-defined profile of targets may be more effective, it is difficult to computationally predict which molecules have desirable target profiles. Here, we report the results of a competitive crowdsourcing project (the Multi-Targeting Drug DREAM Challenge), where we challenged participants to predict chemicals that have desired target profiles for cancer and neurodegenerative disease. Reproduction in deep-sea vent shrimps is shaped by diet and phylogeny, with rhythms unlinked to surface production Variations in offspring production according to feeding strategies or food supply have been recognized in many animals from various ecosystems. Despite an unusual trophic structure based on non-photosynthetic primary production, these relationships remain largely under-studied in chemosynthetic ecosystems. Here, we use Rimicaris shrimps from deep-sea hydrothermal vents as a study case to explore relations between reproduction, diets and food supply in these environments. For that, we compared reproductive outputs of three congeneric shrimps differing by their diets. They inhabit vents located under oligotrophic waters of tropical gyres with opposed latitudes, allowing us to also examine the prevalence of phylogenetic vs environmental drivers in their reproductive rhythms. For this we used both our original data and a compilation of published observations on the presence of ovigerous females covering various seasons over the past 35 years. We report distinct egg production trends between Rimicaris species relying solely on chemosymbiosis - R. exoculata and R. kairei - and those relying on mixotrophy - R. chacei - where R. chacei produces broods with higher numbers of smaller eggs. Besides, our data and historical records suggest a reproductive period with substantial proportions of brooding females mostly between January and early April for all examined species whatever the region. Intriguingly, this periodicity does not correspond to seasonal variations in surface production with presence of brooding females during either boreal winter or austral summer. These observations contrast with the long-standing paradigm in deep-sea species for which periodic reproductive patterns have always been attributed to seasonal variations of photosynthetic production sinking from surface. Our results suggest the presence of intrinsic basis for biological rhythms in the deep sea, and bring to light the importance of having year-round observations in order to understand life history of vent animals. A humanized nanobody phage display library yields potent binders of SARS CoV-2 spike Neutralizing antibodies targeting the SARS-CoV-2 spike protein have shown a great preventative/therapeutic potential. Here, we report a rapid and efficient strategy for the development and design of SARS-CoV-2 neutralizing humanized nanobody constructs with sub-nanomolar affinities and nanomolar potencies. CryoEM-based structural analysis of the nanobodies in complex with spike revealed two distinct binding modes. The most potent nanobody, RBD-1-2G(NCATS-BL8125), tolerates the N501Y RBD mutation and remains capable of neutralizing the B.1.1.7 (Alpha) variant. Molecular dynamics simulations provide a structural basis for understanding the neutralization process of nanobodies exclusively focused on the spike-ACE2 interface with and without the N501Y mutation on RBD. A primary human airway air-lung interface (ALI) ex vivo model showed that RBD-1-2G-Fc antibody treatment was effective at reducing viral burden following WA1 and B.1.1.7 SARS-CoV-2 infections. Therefore, this presented strategy will serve as a tool to mitigate the threat of emerging SARS-CoV-2 variants. Relative model selection of evolutionary substitution models can be sensitive to multiple sequence alignment uncertainty BackgroundMultiple sequence alignments (MSAs) represent the fundamental unit of data inputted to most comparative sequence analyses. In phylogenetic analyses in particular, errors in MSA construction have the potential to induce further errors in downstream analyses such as phylogenetic reconstruction itself, ancestral state reconstruction, and divergence time estimation. In addition to providing phylogenetic methods with an MSA to analyze, researchers must also specify a suitable evolutionary model for the given analysis. Most commonly, researchers apply relative model selection to select a model from candidate set and then provide both the MSA and the selected model as input to subsequent analyses. While the influence of MSA errors has been explored for most stages of phylogenetics pipelines, the potential effects of MSA uncertainty on the relative model selection procedure itself have not been explored.

ResultsWe assessed the consistency of relative model selection when presented with multiple perturbed versions of a given MSA. We find that while relative model selection is mostly robust to MSA uncertainty, in a substantial proportion of circumstances, relative model selection identifies distinct best-fitting models from different MSAs created from the same set of sequences. We find that this issue is more pervasive for nucleotide data compared to amino-acid data. However, we also find that it is challenging to predict whether relative model selection will be robust or sensitive to uncertainty in a given MSA.

ConclusionsWe find that that MSA uncertainty can affect virtually all steps of phylogenetic analysis pipelines to a greater extent than has previously been recognized, including relative model selection. Characterising ontogeny of quantity discrimination in zebrafish A sense of non-symbolic numerical magnitudes is widespread in the animal kingdom and has been documented in adult zebrafish. Here we investigated the ontogeny of this ability using a group size preference task in juvenile zebrafish. Fish showed group size preference from 21 days post fertilization (dpf) and reliably chose the larger group when presented with discriminations of between 1 vs. 3, 2 vs. 5 and 2 vs. 3 conspecifics but not 2 vs. 4 conspecifics. When the ratio between the number of conspecifics in each group was maintained at 1:2, fish could discriminate between 1 vs. 2 individuals and 3 vs. 6, but again, not when given a choice between 2 vs. 4 individuals. These findings are in agreement with studies in other species suggesting the systems involved in quantity representation do not operate separately from other cognitive mechanisms. Rather they suggest quantity processing in fish may be the result of an interplay between attentional, cognitive and memory-related mechanisms as in humans and other animals. Our results emphasise the potential of the use of zebrafish to explore the genetic and neural processes underlying the ontogeny and function of number cognition. In vitro maturation of Toxoplasma gondii bradyzoites in human myotubes and their metabolomic characterization The apicomplexan parasite Toxoplasma gondii forms bradyzoite-containing tissue cysts that cause chronic and drug-tolerant infections. However, current in vitro models do not allow long-term culture of these cysts to maturity. Here, we developed a human myotube-based in vitro culture model of functionally mature tissue cysts that are orally infectious to mice and tolerate exposure to a range of antibiotics and temperature stresses. Metabolomic characterization of purified cysts reveals global changes that comprise increased levels of amino acids and decreased abundance of nucleobase- and tricarboxylic acid cycle-associated metabolites. In contrast to fast replicating tachyzoite forms of T. gondii these tissue cysts tolerate exposure to the aconitase inhibitor sodium fluoroacetate. Direct access to persistent stages of T. gondii under defined cell culture conditions will be essential for the dissection of functionally important host-parasite interactions and drug evasion mechanisms. It will also facilitate the identification of new strategies for therapeutic intervention. Bias correction for inverse variance weighting Mendelian randomization BackgroundInverse-variance weighted two-sample Mendelian randomization (IVW-MR) is the most widely used approach that utilizes genome-wide association studies (GWAS) summary statistics to infer the existence and the strength of the causal effect between an exposure and an outcome. Estimates from this approach can be subject to different biases due to the use of weak instruments and winners curse, which can change as a function of the overlap between the exposure and outcome samples.

MethodsWe developed a method (MRlap) that simultaneously considers weak instrument bias and winners curse, while accounting for potential sample overlap. Assuming spike-and-slab genomic architecture and leveraging LD-score regression and other techniques, we could analytically derive, reliably estimate, and hence correct for the bias of IVW-MR using association summary statistics only.

ResultsWe tested our approach using simulated data for a wide range of realistic settings. In all the explored scenarios, our correction reduced the bias, in some situations by as much as 30 folds. When applied to real data on obesity-related exposures, we observed significant differences between IVW-based and corrected effects, both for non-overlapping and fully overlapping samples. Additionally, our results are consistent with the fact that the strength of the biases will decrease as the sample size increases. Finally, we also showed that the overall bias is also dependent on the genetic architecture of the exposure, and traits with low heritability and/or high polygenicity are more strongly affected.

ConclusionsOur method not only reduces bias in causal effect estimation but also enables the use of much larger GWAS sample sizes, by allowing for potentially overlapping samples. 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. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Defining An Expanded RAS Conformational Landscape Based on Over 700 Experimentally Determined Structures of KRAS, NRAS, and HRAS For many human cancers and tumor-associated diseases, mutations in the RAS isoforms (KRAS, NRAS, and HRAS) are the most common oncogenic alterations, making these proteins high-priority therapeutic targets. Effectively targeting the RAS isoforms requires an exact understanding of their active, inactive, and druggable conformations. However, there is no structure-guided catalogue of RAS conformations to guide therapeutic targeting or examining the structural impact of RAS mutations. We present an expanded classification of RAS conformations based on analyzing their catalytic switch 1 (SW1) and switch 2 (SW2) loops. From all 721 available human KRAS, NRAS, and HRAS structures in the Protein Data Bank (PDB) (206 RAS-protein complexes, 190 inhibitor-bound, and 325 unbound, including 204 WT and 517 mutated structures), we created a broad conformational classification based on the spatial positions of residue Y32 in SW1 and residue Y71 in SW2. Subsequently, we defined additional conformational subsets (some previously undescribed) by clustering all well modeled SW1 and SW2 loops using a density-based machine learning algorithm with a backbone dihedral-based distance metric. In all, we identified three SW1 conformations and nine SW2 conformations, each which are associated with different nucleotide states (GTP-bound, nucleotide-free, and GDP-bound) and specific bound proteins or inhibitor sites. The GTP-bound SW1 conformation can be further subdivided based on the hydrogen (H)-bond type made between residue Y32 and the GTP {gamma}-phosphate: water-mediated, direct, or no H-bond. Further analyzing these structures clarified the catalytic impact of the G12D and G12V RAS mutations, and the inhibitor chemistries that bind to each druggable RAS conformation. To facilitate future RAS structural analyses, we have created a web database, called Rascore, presenting an updated and searchable dataset of human KRAS, NRAS, and HRAS structures in the PDB, and which includes a page for analyzing user uploaded RAS structures by our algorithm (http://dunbrack.fccc.edu/rascore/).

SignificanceAnalyzing >700 experimentally determined RAS structures helped define an expanded landscape of active, inactive and druggable RAS conformations, the structural impact of common RAS mutations, and previously uncharacterized RAS-inhibitor binding modes.