Metabolic reconstitution by a gnotobiotic microbiota varies over the circadian cycle The capacity of the intestinal microbiota to degrade otherwise indigestible diet components is known to greatly improve the recovery of energy from food. This has led to the hypothesis that increased digestive efficiency may underlie the contribution of the microbiome to obesity. OligoMM12-colonized gnotobiotic mice have a consistently higher fat-mass than germ-free or fully colonized counterparts. We therefore investigated their food intake, digestion efficiency, energy expenditure and respiratory quotient using a novel isolator-housed metabolic cage system which allows long-term measurements without contamination risk. This demonstrated that microbiota- released calories are perfectly balanced by decreased food intake in fully colonized versus gnotobiotic OligoMM12 and germ-free mice fed a standard chow diet, i.e., microbiota-released calories can in fact be well-integrated into appetite control. We also observed no significant difference in energy expenditure per gram lean mass between the different microbiota groups, suggesting that cumulative very small differences in energy balance, or altered energy storage must underlie fat accumulation in OligoMM12 mice. Consistent with altered energy storage, major differences were observed in the type of respiratory substrates used in metabolism over the circadian cycle: in germ-free mice the respiratory exchange ratio was consistently lower than that of fully colonized mice at all times of day, indicative of more reliance on fat and less on glucose metabolism. Intriguingly the RER of OligoMM12-colonized gnotobiotic mice phenocopied fully colonized mice during the dark (active/eating) phase but phenocopied germ-free mice during the light (fasting/resting) phase. Further, OligoMM12-colonized mice showed a germ-free-like drop in liver glycogen storage during the light cycle and both liver and plasma metabolomes of OligoMM12 mice clustered closely with germ-free mice. This implies the existence of microbiota functions that are required to maintain normal host metabolism during the resting/fasting phase of circadian cycle, and which are absent in the OligoMM12 consortium.
Deciphering the potential niche of two novel black yeast fungi from a biological soil crust based on their genomes, phenotypes, and melanin regulation Black yeasts are polyextremotolerant fungi that contain high amounts of melanin in their cell wall and maintain a primarily yeast form. These fungi grow in xeric, nutrient deplete environments which implies that they require highly flexible metabolisms and the ability to form lichen-like mutualisms with nearby algae and bacteria. However, the exact ecological niche and interactions between these fungi and their surrounding community is not well understood. We have isolated two novel black yeast fungi of the genus Exophiala: E. viscosium and E. limosus, which are from dryland biological soil crusts. A combination of whole genome sequencing and various phenotyping experiments have been performed on these isolates to determine their fundamental niches within the biological soil crust consortium. Our results reveal that these Exophiala spp. are capable of utilizing a wide variety of carbon and nitrogen sources potentially from symbiotic microbes, they can withstand many abiotic stresses, and can potentially provide UV resistance to the crust community in the form of secreted melanin. Besides the identification of two novel species within the genus Exophiala, our study also provides new insight into the production and regulation of melanin in extremotolerant fungi.
Protein scaffold-based multimerization of soluble ACE2 efficiently blocks SARS-CoV-2 infection in vitro and in vivo Soluble ACE2 (sACE2) decoy receptors are promising agents to inhibit SARS-CoV-2, as their efficiency is less likely to be affected by common escape mutations in viral proteins. However, their success may be limited by their relatively poor potency. To address this challenge, we developed a large decoy library of sACE2 fusion proteins, generated with several protease inhibitors or multimerization tags. Among these decoys, multimeric sACE2 consisting of SunTag or MoonTag systems, which were originally utilized for signal amplification or gene activation systems, were extremely effective in neutralizing SARS-CoV-2 in pseudoviral systems and in clinical isolates. These novel sACE2 fusion proteins exhibited greater than 100-fold SARS-CoV-2 neutralization efficiency, compared to monomeric sACE2. SunTag or MoonTag in combination with a more potent version of sACE2, which has multiple point mutations for greater binding (v1), achieved near complete neutralization at a sub-nanomolar range, comparable with clinical monoclonal antibodies. Pseudoviruses bearing mutant versions of Spike, alpha, beta, gamma or delta variants, were also neutralized efficiently with SunTag or MoonTag fused sACE2(v1). Finally, therapeutic treatment of sACE2(v1)-MoonTag provided protection against SARS-CoV-2 infection in an in vivo mouse model. Overall, we suggest that the superior activity of the sACE2-SunTag or sACE2-MoonTag fusions is due to the greater occupancy of the multimeric sACE2 receptors on Spike protein as compared to monomeric sACE2. Therefore, these highly potent multimeric sACE2 decoy receptors may offer a promising treatment approach against SARS-CoV-2 infections.

One Sentence SummaryMultimerization of sACE2 markedly enhanced the neutralization of SARS-CoV-2 by blocking multiple viral spike proteins simultaneously. 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. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. Common coupling map advances GPCR-G protein selectivity Two-thirds of human hormones and one-third of clinical drugs act on membrane receptors that couple to G proteins to achieve appropriate functional responses. While G protein transducers from literature are annotated in the Guide to Pharmacology database, two recent large-scale datasets now expand the receptor-G protein  couplome. However, these three datasets differ in scope and reported G protein couplings giving different coverage and conclusions on GPCR-G protein signaling. Here, we report a common coupling map uncovering novel couplings supported by both large-scale studies, the selectivity/promiscuity of GPCRs and G proteins, and how the co-coupling and co-expression of G proteins compare to the families from phylogenetic relationships. The coupling map and insights on GPCR-G protein selectivity will catalyze advances in receptor research and cellular signaling towards the exploitation of G protein signaling bias in design of safer drugs. GWAS-associated variants, non-genetic factors, and transient transcriptome in Multiple Sclerosis etiopathogenesis: a colocalization analysis A clinically actionable understanding of multiple sclerosis (MS) etiology goes through GWAS interpretation, prompting research on new gene regulatory models. We previously suggested a stochastic etiologic model where small-scale random perturbations could reach a threshold for MS development. The recently described mapping of the transient transcriptome (TT), including intergenic and intronic RNAs, seems appropriate to verify this model through a rigorous colocalization analysis. We show that genomic regions coding for the TT were significantly enriched for MS-associated GWAS variants and DNA binding sites for molecular transducers mediating putative, non-genetic, etiopathogenetic factors for MS (e.g., vitamin D deficiency, Epstein Barr virus latent infection, B cell dysfunction). These results suggest a model whereby TT-coding regions are hotspots of convergence between genetic ad non-genetic factors of risk/protection for MS, and plausibly for other complex disorders. Our colocalization analysis also provides a freely available data resource (www.mscoloc.com) for future research on MS transcriptional regulation. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. Two Types of Cinnamoyl-CoA Reductase Function Divergently in Tissue Lignification, Phenylpropanoids Flux Control, and Inter-pathway Cross-talk with Glucosinolates as Revealed in Brassica napus Cinnamoyl-CoA reductase (CCR) is the entry point of lignin pathway and a crucial locus in dissection and manipulation of associated traits, but its functional dissection in Brassicaceae plants is largely lagged behind though Arabidopsis thaliana CCR1 has been characterized to certain extent. Here, 16 CCR genes are identified from Brassica napus and its parental species B. rapa and B. oleracea. Brassicaceae CCR genes are divided into CCR1 subfamily and CCR2 subfamily with divergent organ-specificity, yellow-seed trait participation and stresses responsiveness. CCR1 is preferential in G- and H-lignins biosynthesis and vascular development, while CCR2 has a deviation to S-lignin biosynthesis and interfascicular fiber development. CCR1 has stronger effects on lignification-related development, lodging resistance, phenylpropanoid flux control and seed coat pigmentation, whereas CCR2 controls sinapates levels. CCR1 upregulation could delay bolting and flowering time, while CCR2 upregulation weakens vascular system in leaf due to suppressed G lignin accumulation. Besides, CCR1 and CCR2 are deeply but almost oppositely linked with glucosinolates metabolism through inter-pathway crosstalk. Strangely, upregulation of both CCR1 and CCR2 could not enhance resistance to UV-B and S. sclerotiorum though CCR2 is sharply induced by them. These results provide systemic dissection on Brassica CCRs and CCR1-CCR2 divergence in Brassicaceae.

HighlightBrassicaceae contains two types of Cinnamoyl-CoA reductase. As revealed in Brassica napus, they are divergently involved in lignin monomer biosynthesis, tissue lignification, phenylpropanoid flux control, and inter-pathway crosstalk with glucosinolates. A mitochondria-specific mutational signature of aging: increased rate of A>G substitutions on a heavy chain. The mutational spectrum of the mitochondrial DNA (mtDNA) does not resemble any of the known mutational signatures of the nuclear genome and variation in mtDNA mutational spectra between different organisms is still incomprehensible. Since mitochondria is tightly involved in aerobic energy production, it is expected that mtDNA mutational spectra is affected by the oxidative damage. Assuming that oxidative damage increases with age, we analyze mtDNA mutagenesis of different species. Analysing (i) dozens thousands of somatic mtDNA mutations in samples of different age (ii) 70053 polymorphic synonymous mtDNA substitutions, reconstructed in 424 mammalian species with different generation length and (iii) synonymous nucleotide content of 650 complete mitochondrial genomes of mammalian species we observed that the frequency of AH>GH substitutions (H - heavy chain notation) is twice higher in species with high versus low generation length making their mtDNA more AH poor and GH rich. Considering that AH>GH substitutions are also sensitive to the time spent single stranded (TSSS) during asynchroniuos mtDNA replication we demonstrated that AH>GH substitution rate is a function of both species-specific generation length and position specific TSSS. We propose that AH>GH is a mitochondria-specific signature of oxidative damage associated with both aging and TSSS. Selection towards lighter seeds in Danish grasslands over an eight-year period Time-series analysis of cover data from 236 Danish grassland sites demonstrated community selection on seed mass. Across four grassland habitat types and during an eight-year period there was significant community selection towards lighter seeds. In the analysis, the continuous seed size variable was used to group plant species into functional types. This method was chosen in order to account for the sampling process of the abundance data, and is in contrast to most other analyses of trait selection, where the community weighted mean of the traits is used as the dependent variable. Generally, little is known on the underlying causes of possible plant trait community selection in semi-natural grasslands that are not subjected to successional processes. Inter-tissue convergence of gene expression during ageing suggests age-related loss of tissue and cellular identity Developmental trajectories of gene expression may reverse in their direction during ageing, a phenomenon previously linked to cellular identity loss. Our analysis of cerebral cortex, lung, liver and muscle transcriptomes of 16 mice, covering development and ageing intervals, revealed widespread but tissue-specific ageing-associated expression reversals. Cumulatively, these reversals create a unique phenomenon: mammalian tissue transcriptomes diverge from each other during postnatal development, but during ageing, they tend to converge towards similar expression levels, a process we term Divergence followed by Convergence, or DiCo. We found that DiCo was most prevalent among tissue-specific genes and associated with loss of tissue identity, which is confirmed using data from independent mouse and human datasets. Further, using publicly available single-cell transcriptome data, we showed that DiCo could be driven both by alterations in tissue cell type composition and also by cell-autonomous expression changes within particular cell types. 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. New multiscale characterisation methodology for effective determination of isolationstructure-function relationship of extracellular vesicles. Extracellular vesicles (EVs) have been lauded as next generation medicines, but very few EV-based therapeutics have progressed to clinical use. Limited clinical translation is largely due to technical barriers that hamper our ability to mass-produce EVs, i.e. to isolate, purify and characterise them effectively. Technical limitations in comprehensive characterisation of EVs leads to unpredicted biological effects of EVs. Here, using a range of optical and non-optical techniques, we showed that the differences in molecular composition of EVs isolated using two isolation methods correlated with the differences in their biological function. Our results demonstrated that the isolation method determines the composition of isolated EVs at single and sub-population levels. Besides the composition, we measured for the first time the dry mass and predicted sedimentation of EVs. These parameters were shown to correlate well with the biological and functional effects of EVs on single cell and cell cultures. We anticipate that our new multiscale characterisation approach, which goes beyond traditional experimental methodology, will support fundamental understanding of EVs as well as elucidate the functional effects of EVs in in vitro and in vivo studies. Our findings and methodology will be pivotal for developing optimal isolation methods and establishing EVs as mainstream therapeutics and diagnostics. This innovative approach is applicable to a wide range of sectors including biopharma and biotechnology as well as to regulatory agencies. RIG-I immunotherapy overcomes radioresistance in p53-positive malignant melanoma Radiation therapy induces cytotoxic DNA damage, which results in cell-cycle arrest and activation of cell-intrinsic death pathways, but its application has been limited by the radioresistance of tumors, such as in malignant melanoma. RIG-I is a cytosolic immune receptor expressed in all somatic cells, including tumor cells, with a key role in sensing viral RNA. RIG-I specific oligonucleotide ligands elicit a robust cell-intrinsic antiviral response and immunogenic cell death in tumor cells and are being tested in clinical trials. Nonetheless, their potential to overcome radioresistance has not yet been explored. Here, we demonstrate that activation of RIG-I enhances the extent and immunogenicity of irradiation-induced tumor cell death in human and murine melanoma cell lines in vitro and improved survival in the murine B16 melanoma model. Pathway analysis of transcriptomic data revealed a central role for p53 downstream of the combination treatment, which was corroborated using p53-/- B16 cells. In vivo, the effect of irradiation on immune-cell activation and inhibition of tumor growth was absent in mice carrying p53-/- B16 tumors, while the response to RIG-I stimulation in those mice was maintained. Our results identify p53 as pivotal for the synergistic antitumoral effect of RIG-I and irradiation, resulting in potent induction of immunogenic tumor-cell death. Thus, the administration of RIG-I ligands in combination with radiotherapy is a promising therapeutic approach to treating radioresistant tumors with a functional p53 pathway, such as malignant melanoma. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. NUDT18 catalyzes the hydrolysis of active metabolites of the antivirals Remdesivir, Ribavirin and Molnupiravir Remdesivir and Molnupiravir have gained considerable interest due to their activity against SARS-CoV-2. Cellular hydrolysis of their active triphosphate forms, Remdesivir-TP and Molnupiravir-TP, would decrease drug efficiency. We therefore tested Remdesivir-TP as a substrate against a panel of human hydrolases and found that NUDT18 catalyzes the hydrolysis of Remdesivir-TP. The kcat value of NUDT18 for Remdesivir-TP was determined to 2.6 s-1 and the Km value was 156 M, suggesting that NUDT18 catalyzed hydrolysis of Remdesivir-TP occurs in cells. We demonstrate that the triphosphates of the antivirals Ribavirin and Molnupiravir are hydrolyzed by NUDT18, albeit with a lower efficiency compared to Remdesivir-TP. NUDT18 also hydrolyses the triphosphates of Sofosbuvir and Aciclovir although with significantly lower activity. These results suggest that NUDT18 can act as a cellular sanitizer of modified nucleotides and may influence the antiviral efficacy of Remdesivir, Molnupiravir and Ribavirin. NUDT18 is expressed in respiratory epithelial cells and may limit the antiviral efficacy of Remdesivir and Molnupiravir against SARS-CoV2 replication by decreasing the intracellular concentration of their active metabolites at their intended site of action. Phenotypic heterogeneity driven by plasticity of the intermediate EMT state governs disease progression and metastasis in breast cancer The Epithelial-to-Mesenchymal Transition (EMT) is a developmental cellular program frequently coopted by cancer cells1 and is a key contributor to both heterogeneity in solid tumors2-4 and later stage chemo-resistance and metastasis5,6. Rather than being a switch from an epithelial to a mesenchymal state, increasing evidence points to the existence of intermediate EMT states, wherein cells co-express both epithelial and mesenchymal traits7-13. Multiple stable intermediate EMT states possessing unique characteristics exist across the EMT spectrum7,8,14,15, contributing to the complex heterogeneity of tumors to promote metastasis16. While much work has been carried out identifying and characterizing EMT-inducing transcription factors17-20, the transcriptional and epigenetic networks responsible for the stability and maintenance of the midpoints along the EMT spectrum are poorly defined. In addition, there are currently no approaches to identifying and quantifying intermediate EMT subpopulations within patient tumors to evaluate their prognostic significance. Using clonally isolated derivatives of the SUM149PT breast cancer cell line, we systematically interrogate how each EMT state independently contributes to heterogeneity and influences metastatic progression, uncovering the role of RUNX2 in stabilizing certain intermediate states. Using SUM149PT-derived tumors as a training set, we develop an entropy-based model to quantify phenotypic heterogeneity and EMT status. Remarkably, the majority of cell states captured in the SUM149PT model are represented in a selection of patient tumors, laying the foundation for quantification of epithelial-mesenchymal heterogeneity and understanding the role of the intermediate EMT state in tumor progression. ConVarT: a search engine for orthologous variants and functional inference of human genetic variants The availability of genetic variants, together with phenotypic annotations from model organisms, facilitates comparing these variants with equivalent variants in humans. However, existing databases and search tools do not make it easy to scan for equivalent variants, namely "matching variants" (MatchVars) between humans and other organisms. Therefore, we developed an integrated search engine called ConVarT (http://www.convart.org/) for matching variants between humans, mice, and C. elegans. ConVarT incorporates annotations (including phenotypic and pathogenic) into variants, and these previously unexploited phenotypic MatchVars from mice and C. elegans can give clues about the functional consequence of human genetic variants. Our analysis shows that many phenotypic variants in different genes from mice and C. elegans, so far, have no counterparts in humans, and thus, can be useful resources when evaluating a relationship between a new human mutation and a disease. Diffusion and distal linkages govern interchromosomal dynamics during meiotic prophase The pairing of homologous chromosomes (homologs) in meiosis is essential for distributing the correct numbers of chromosomes into haploid gametes. In budding yeast, pairing depends on the formation of 150-200 Spo11-mediated double-strand breaks (DSBs) that are distributed among 16 homolog pairs, but it is not known if all, or only a subset of these DSBs, contribute to the close juxtaposition of homologs. Having established a system to measure the position of fluorescently tagged chromosomal loci in 3D space over time, we analyzed locus trajectories to determine how frequently, and how long, loci spend colocalized or apart. Continuous imaging revealed highly heterogeneous cell-to-cell behavior of foci, with the majority of cells exhibiting a "mixed" phenotype where foci move into and out of proximity, even at late stages of prophase, suggesting that the axial structures of the synaptonemal complex may be more dynamic than anticipated. The observed plateaus of the mean-squared change in distance (MSCD) between foci informed the development of a biophysical model of two diffusing polymers that captures the loss of centromere linkages as cells enter meiosis, nuclear confinement, and the formation of Spo11-dependent linkages. The predicted number of linkages per chromosome in our theoretical model closely approximates the small number (~2-4) of estimated synapsis-initiation sites, suggesting that excess DSBs have negligible effects on the overall juxtaposition of homologs. These insights into the dynamic in-terchromosomal behavior displayed during homolog pairing demonstrate the power of combining time-resolved in vivo analysis with modeling at the granular level.

Significance StatementEssential for sexual reproduction, meiosis is a specialized cell division required for the production of haploid gametes. Critical to this process is the pairing, recombination, and segregation of homologous chromosomes (homologs). While pairing and recombination are linked, it is not known how many linkages are sufficient to hold homologs in proximity. Here, we reveal that random diffusion and the placement of a small number of linkages are sufficient to establish the apparent "pairing" of homologs. We also show that colocalization between any two loci is more dynamic than anticipated. Our study is the first to provide observations of live interchromosomal dynamics during meiosis and illustrates the power of combining single-cell measurements with theoretical polymer modeling. A chemical screen based on an interruption of zebrafish gastrulation identifies the HTR2C inhibitor Pizotifen as a suppressor of EMT-mediated metastasis Metastasis is responsible for approximately 90% of cancer-associated mortality but few models exist that allow for rapid and effective screening of anti-metastasis drugs. Current mouse models of metastasis are too expensive and time consuming to use for rapid and high-throughput screening. Therefore, we created a unique screening concept utilizing conserved mechanisms between zebrafish gastrulation and cancer metastasis for identification of potential anti-metastatic drugs. We hypothesized that small chemicals that interrupt zebrafish gastrulation might also suppress metastatic progression of cancer cells and developed a phenotype-based chemical screen to test the hypothesis. The screen used epiboly, the first morphogenetic movement in gastrulation, as a marker and enabled 100 chemicals to be tested in five hours. The screen tested 1280 FDA-approved drugs and identified Pizotifen, an antagonist for serotonin receptor 2C (HTR2C) as an epiboly-interrupting drug. Pharmacologic and genetic inhibition of HTR2C suppressed metastatic progression in a mouse model. Blocking HTR2C with Pizotifen restored epithelial properties to metastatic cells through inhibition of Wnt-signaling. In contrast, HTR2C induced epithelial to mesenchymal transition (EMT) through activation of Wnt-signaling and promoted metastatic dissemination of human cancer cells in a zebrafish xenotransplantation model. Taken together, our concept offers a novel platform for discovery of anti-metastasis drugs. Define protein variant functions with high-complexity mutagenesis libraries and enhanced mutation detection software ASMv1.0 Open reading frame (ORF) variant libraries have advanced our ability to query the functions of a large number of variants of a protein simultaneously in a single experiment. Variant libraries targeting full-length ORFs typically consists of all possible single-amino-acid substitutions and a stop codon at each amino-acid position. Because a variant differs from the template ORF by a single codon variation, variant quantification presents the most profound challenge to this technology. Efforts such as dividing a library into sub-libraries for direct sequencing, or tag-directed subassembly are practical only for small ORFs. Our approach, however, features an enhanced variant-calling tool that processes massive shotgun sequencing data, allowing one to screen a single variant library for ORFs sized up to [~]3600 bases. The variant-calling tool detects variants reliably, and also presents variant annotations in datafiles enabling analyses that have reshaped our strategies governing library design, screen deconvolution, sequencing and sequencing data analysis. 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. Impact of circadian time of dosing on cardiomyocyte-autonomous effects of glucocorticoids Mitochondrial capacity is critical to adapt the high energy demand of the heart to circadian oscillations and diseased states. Glucocorticoids regulate the circadian cycle of energy metabolism, but little is known about how circadian timing of exogenous glucocorticoid dosing directly regulates heart metabolism through cardiomyocyte-autonomous mechanisms. While chronic oncedaily intake of glucocorticoids promotes metabolic stress and heart failure, we recently discovered that intermittent once-weekly dosing of exogenous glucocorticoids promoted muscle metabolism in normal and obese skeletal muscle. However, the effects of glucocorticoid intermittence on heart metabolism and heart failure remain unknown. Here we investigated the extent to which circadian time of dosing regulates the effects of the glucocorticoid prednisone in heart metabolism and function in conditions of single pulse or chronic intermittent dosing. In WT mice, we found that prednisone improved cardiac content of NAD+ and ATP with light-phase dosing (ZT0), while the effects were blocked by dark-phase dosing (ZT12). The drug effects on mitochondrial function were cardiomyocyte-autonomous, as shown by inducible cardiomyocyte-restricted glucocorticoid receptor (GR) ablation, and depended on an intact cardiomyocyte clock, as shown by inducible cardiomyocyte-restricted ablation of Brain and Muscle ARNT-like 1 (BMAL1). Conjugating time-of-dosing with chronic intermittence, we found that once-weekly prednisone improved metabolism and function in heart after myocardial injury dependent on circadian time of intake, i.e. with lightphase but not dark-phase dosing. Our study identifies cardiac-autonomous mechanisms through which circadian-specific intermittent dosing reconverts glucocorticoid drugs to metabolic boosters for the heart. Stimulus manipulations permit activation of fiber subpopulations in the mouse and rat vagus Afferent and efferent vagal fibers mediate bidirectional communication between the brain and visceral organs. Small, unmyelinated C-afferents constitute the majority of vagal fibers, play critical roles in numerous interoceptive circuits and autonomic reflexes in health and disease and may contribute to the efficacy and safety of vagus nerve stimulation (VNS). Selective engagement of C-afferents with electrical stimuli has not been feasible, due to the default fiber recruitment order: larger fibers first, smaller fibers last. Here, we determine and optimize an electrical stimulus that selectively engages vagal C-afferents. Intermittent KHz-frequency electrical stimulation (KES) activates motor and, preferentially, sensory vagal neurons in the brainstem. During KES, asynchronous activity of C-afferents increases, while that of larger fibers remains largely unchanged. In parallel, KES effectively blocks excitability of larger fibers while moderately suppressing excitability of C-afferents. By compiling selectivity indices in individual animals, we find that optimal KES parameters for C-afferents are >5KHz frequency and 7-10 times engagement threshold (xT) intensity in rats, 15-25xT in mice. These effects can be explained in computational models by how sodium channel responses to KES are shaped by axonal size and myelin. Our results indicate that selective engagement of vagal C-afferents is attainable by intermittent KES. Parallel introgression, not recurrent emergence, explains apparent elevational ecotypes of polyploid Himalayan snowtrout The recurrence of similar evolutionary patterns within different habitats often reflects parallel selective pressures acting upon either standing or independently occurring genetic variation to produce a convergence of phenotypes. This interpretation (i.e. parallel divergences within adjacent streams) has been hypothesized for drainage-specific morphological  ecotypes observed in polyploid snowtrout (Cyprinidae: Schizothorax). However, parallel patterns of differential introgression during secondary contact is a viable alternative hypothesis. Here, we utilized ddRADseq (N=35,319 de-novo and N=10,884 transcriptome-aligned SNPs), as derived from Nepali/Bhutanese samples (N=48 each), to test these competing hypotheses. We first employed genome-wide allelic depths to derive appropriate ploidy models, then a Bayesian approach to yield genotypes statistically consistent under the inferred expectations. Elevational  ecotypes were consistent in geometric-morphometric space, but with phylogenetic relationships at the drainage-level, sustaining an hypothesis of independent emergence. However, partitioned analyses of phylogeny and admixture identified subsets of loci under selection that retained genealogical concordance with morphology, suggesting instead that apparent patterns of morphological/phylogenetic discordance are driven by widespread genomic homogenization. Here, admixture occurring in secondary contact effectively  masks previous isolation. Our results underscore two salient factors:1) Morphological adaptations are retained despite hybridization; and 2) The degree of admixture varies across tributaries, presumably concomitant with underlying environmental or anthropogenic factors.