Individual beliefs about temporal continuity explain variation of perceptual biases Perception of magnitudes such as duration or distance is often found to be systematically biased. The biases, which result from incorporating prior knowledge in the perceptual process, can vary considerably between individuals. The variations are commonly attributed to differences in sensory precision and reliance on priors. However, another factor not considered so far is the implicit belief about how successive sensory stimuli are generated: independently from each other or with certain temporal continuity. The main types of explanatory models proposed so far - static or iterative - mirror this distinction but cannot adequately explain individual biases. Here we propose a new unifying model that explains individual variation as combination of sensory precision and beliefs about temporal continuity and predicts the experimentally found changes in biases when altering temporal continuity. Thus, according to the model, individual differences in perception depend on beliefs about how stimuli are generated in the world.
Trial and error: a hierarchical modeling approach to test-retest assessment The concept of test-retest reliability indexes the consistency of a measurement across time. High reliability is critical for any scientific study, but specifically for the study of individual differences. Evidence of poor reliability of commonly used behavioral and functional neuroimaging tasks is mounting. Reports on low reliability of task-based fMRI have called into question the adequacy of using even the most common, well-characterized cognitive tasks with robust population-level effects, to measure individual differences. Here, we lay out a hierarchical framework that estimates reliability as a correlation divorced from trial-level variability, and show that reliability estimates tend to be higher compared to the conventional framework that adopts condition-level modeling and ignores across-trial variability. We examine how estimates from the two frameworks diverge and assess how different factors (e.g., trial and subject sample sizes, relative magnitude of cross-trial variability) impact reliability estimates. We also show that, under specific circumstances, the two statistical frameworks converge. Results from the two approaches are approximately equivalent if (a) the trial sample size is sufficiently large, or (b) cross-trial variability is in the same order of magnitude as, or less than, cross-subject variability. As empirical data indicate that cross-trial variability is large in most tasks, this work highlights that a large number of trials (e.g., greater than 100) may be required to achieve precise reliability estimates. We reference the tools TRR and 3dLMEr for the community to apply trial-level models to behavior and neuroimaging data and discuss how to make these new measurements most useful for current studies.
Somatic mutations of activated signaling genes, transcription factors, or tumor suppressors are a precondition for leukemic transformation from myelodysplastic syndromes: a sequencing analysis of 64 paired samples The transformation biology of secondary AML from MDS is still not fully understood. Here, we performed a large cohort of paired sequences including target, whole-exome and single cell sequencing to search AML transformation-related mutations (TRM). The results showed that fifty-five out of the 64 (85.9%) patients presented presumptive TRM involving activated signaling, transcription factors, or tumor suppressors. Most of TRM (63.6%, 35 cases) emerged at the leukemia transformation point. All five of the remaining nine patients analyzed by paired whole exome sequencing showed TRM which are not included in the reference targets. Single-cell sequencing indicated that the activated cell signaling route was related to TRM which take place prior to phenotypic development. Of note, defined TRM was limited to a small set of genes (less than ten, in the order: NRAS/KRAS, CEBPA, TP53, FLT3, RUNX1, CBL, PTPN11 and WT1, accounted for 91.0% of the mutations). In conclusion, somatic mutations involving in activated signaling, transcription factors, or tumor suppressors appeared to be a precondition for AML transformation from myelodysplastic syndromes. The TRM may be considered as new therapy targets. An epigenomic shift in amygdala marks the transition to maternal behaviors in alloparenting virgin female mice In many species, adults care for young offspring that are not their own, a phenomenon called alloparenting. However, most nonparental adults must be sensitized by repeated or extended exposures to newborns before they will robustly display parental-like behaviors. To capture neurogenomic events underlying the transition to active parental caring behaviors, we analyzed brain gene expression and chromatin profiles of virgin female mice co-housed with mothers during pregnancy and after birth. After an initial display of antagonistic behaviors and a surge of defense-related gene expression, we observed a dramatic shift in the chromatin landscape specifically in amygdala of the pup-exposed virgin females, accompanied by a dampening of anxiety-related gene expression. This epigenetic shift coincided with hypothalamic expression of the oxytocin gene and the emergence of behaviors and gene expression patterns classically associated with maternal care. The results outline a neurogenomic program associated with dramatic behavioral changes and suggest molecular networks relevant to human postpartum mental health. Using machine learning to design adeno-associated virus capsids with high likelihood of viral assembly We study the application of machine learning in designing adeno-associated virus (AAV) capsid sequences with high likelihood of viral assembly, i.e. capsid viability. Specifically, we design and implement Origami, a model-based optimization algorithm, to identify highly viable capsid sequences within the vast space of 2033 possibilities. Our evaluation shows that Origami performs well in terms of optimality and diversity of model-designed sequences. Moreover, these sequences are ranked according to their viability score. This helps designing experiments given budget constraint. 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. A Neurogenic Signature Involving Monoamine Oxidase-A controls Human Thermogenic Adipose Tissue Development Mechanisms that control "beige/brite" thermogenic adipose tissue development may be harnessed to improve human metabolic health. To define these mechanisms, we developed a species-hybrid model in which human mesenchymal progenitor cells were used to develop white or thermogenic/beige adipose tissue in mice. The hybrid adipose tissue developed distinctive features of human adipose tissue, such as larger adipocyte size, despite its neurovascular architecture being entirely of murine origin. Thermogenic adipose tissue recruited a denser, qualitatively distinct vascular network, differing in genes mapping to circadian rhythm pathways, and denser sympathetic innervation. The enhanced thermogenic neurovascular network was associated with human adipocyte expression of THBS4, TNC, NTRK3 and SPARCL1, which enhance neurogenesis, and decreased expression of MAOA and ACHE, which control neurotransmitter tone. Systemic inhibition of MAOA, which is present in human but absent in mouse adipocytes, induced browning of human but not mouse adipose tissue, revealing the physiological relevance of this pathway. Our results reveal species-specific cell type dependencies controlling the development of thermogenic adipose tissue and point to human adipocyte MAOA as a potential target for metabolic disease therapy. Priority conservation areas and a global population estimate for the Critically Endangered Philippine Eagle derived from modelled range metrics using remote sensing habitat characteristics Many range-restricted taxa are currently experiencing population declines yet lack fundamental information regarding distribution and population size. Establishing baseline estimates for both these key biological parameters is however critical for directing conservation planning for at-risk range-restricted species. The International Union for the Conservation of Nature (IUCN) Red List uses three range metrics that define species distributions and inform extinction risk assessments: extent of occurrence (EOO), area of occupancy (AOO) and area of habitat (AOH). However, calculating all three metrics using standard IUCN approaches relies on a geographically representative sample of locations, which for rare species is often spatially biased. Here, we apply model-based interpolation using Species Distribution Models (SDMs), correlating occurrences with remote-sensing covariates, to calculate IUCN range metrics, protected area coverage and a global population estimate for the Critically Endangered Philippine Eagle (Pithecophaga jefferyi). Our final range wide continuous SDM had high predictive accuracy (Continuous Boyce Index = 0.927) and when converted to a binary model estimated an AOH = 23,185 km2, a maximum EOO = 605,759 km2, a minimum EOO = 272,272 km2, with an AOO = 53,867 km2. Based on inferred habitat from the AOH metric, we estimate a global population of 318 breeding pairs (range: 258-362 pairs), or 636 mature individuals, across the Philippine Eagle global range. Protected areas covered 34 % of AOH, 15 % less than the target representation, with the continuous model identifying key habitat as priority conservation areas. We demonstrate that even when occurrences are geographically biased, robust habitat models can be built that enable quantification of baseline IUCN range metrics, protected area coverage, and a population size estimate. In the absence of adequate location data for many rare and threatened taxa, our method is a promising spatial modelling tool with widespread applications, in particular for island endemics facing high extinction risk. SERCA function is impaired in skeletal and cardiac muscles from young DBA/2J mdx mice The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. Single-dose AAV-based vaccine induces a high level of neutralizing antibodies against SARS-CoV-2 in rhesus macaques Coronavirus disease 2019 (COVID-19), which is triggered by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, continues to threaten global public health. Developing a vaccine that only requires single immunization but provides long-term protection for the prevention and control of COVID-19 is important. Here, we developed an adeno-associated virus (AAV)-based vaccine expressing a stable receptor-binding domain (SRBD) protein. The vaccine requires only a single shot but provides effective neutralizing antibodies (NAbs) over 598 days in rhesus macaques (Macaca mulatta). Importantly, our results showed that the NAbs were kept in high level and long lasting against authentic wild-type SARS-CoV-2, Beta, Delta and Omicron variants using plaque reduction neutralization test. Of note, although we detected pre-existing AAV2/9 antibodies before immunization, the vaccine still induced high and effective NAbs against COVID-19 in rhesus macaques. AAV-SRBD immune serum also efficiently inhibited the binding of ACE2 with RBD in the SARS-CoV-2 B.1.1.7 (Alpha), B.1.351 (Beta), P.1/P.2 (Gamma), B.1.617.2 (Delta), B.1.617.1/3(Kappa), and C.37 (Lambda) variants. Thus, these data suggest that the vaccine has great potential to prevent the spread of SARS-CoV-2. A novel 5 prime DNA binding site in RFC facilitates PCNA loading for gap DNA repair RFC uses ATP to assemble PCNA onto primed sites for replicative DNA polymerases {delta} and {gamma}. The RFC pentamer forms a central chamber that binds 3' ss/ds DNA junctions to load PCNA onto DNA during replication. We show here five structures revealing a novel 2nd DNA binding site in RFC that binds a 5' duplex. This 5' DNA site is located in the N-terminal BRCT domain and AAA+ module of the large Rfc1 subunit. It appears 5' DNA binds after 3' DNA binding to RFC. Our structures reveal ideal binding to a 7-nt gap, which includes 2 bp unwound by the clamp loader. Biochemical studies show enhanced binding to 5 and 10 nt gaps, consistent with the structural results. We propose that the 5' site facilitates RFCs PCNA loading activity at a DNA damage-induced gap to recruit gap-filling polymerases. These findings are consistent with genetic studies showing that base excision repair of gaps greater than 1 base requires PCNA and involves the 5' DNA binding domain of Rfc1. We further observe that a 5' end facilitates PCNA loading at an RPA coated 50-nt gap, suggesting a potential role of the RFC 5' DNA site in longer gaps including during lagging strand DNA synthesis. A revised view on the evolution of glutamine synthetase isoenzymes in plants Glutamine synthetase (GS) is a key enzyme responsible for the incorporation of inorganic nitrogen in the form of ammonium into the amino acid glutamine. The genes encoding GS are among the oldest existing genes in living organisms. In plants, two groups of functional GS enzymes are found: eubacterial GSIIb (GLN2) and eukaryotic GSIIe (GLN1/GS). Phylogenetic analyses have shown that the GLN2 group originated from bacteria following horizontal gene transfer. Only GLN1/GS genes are found in vascular plants, which suggests that they are involved in the final adaptation of plants to terrestrial life. The present phylogenetic study reclassifies the different GS of seed plants into three clusters: GS1a, GS1b and GS2. The presence of genes encoding GS2 has been expanded to Cycadopsida gymnosperms, which suggests the origin of this gene in a common ancestor of Cycadopsida, Ginkgoopsida and angiosperms. GS1a genes have been identified in all gymnosperms, basal angiosperms and some Magnoliidae species. Previous studies in conifers and the gene expression profiles obtained in ginkgo and magnolia in the present work could explain the absence of GS1a in more recent angiosperm species (e.g., monocots and eudicots) due to the redundant roles of GS1a and GS2 in photosynthetic cells. Altogether, the results provide a better understanding of the evolution of plant GS isoenzymes and their physiological roles, which is valuable for improving crop nitrogen use efficiency and productivity. Intragenomic rearrangements of SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses Variation of the betacoronavirus SARS-CoV-2 has been the bane of COVID-19 control. Documented variation includes point mutations, deletions, insertions, and recombination among closely or distantly related coronaviruses. Here, we describe yet another aspect of genome variation by beta- and alphacoronaviruses. Specifically, we report numerous genomic insertions of 5-untranslated region sequences into coding regions of SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses. To our knowledge this is the first systematic description of such insertions. In many cases, these insertions change viral protein sequences and further foster genomic flexibility and viral adaptability through insertion of transcription regulatory sequences in novel positions within the genome. Among human Embecorivus betacoronaviruses, for instance, from 65% to all of the surveyed sequences in publicly available databases contain 5-UTR-derived inserted sequences. In limited instances, there is mounting evidence that these insertions alter the fundamental biological properties of mutant viruses. Intragenomic rearrangements add to our appreciation of how variants of SARS-CoV-2 and other beta- and alphacoronaviruses may arise.

SignificanceUnderstanding mechanisms of variation in coronaviruses is vital to control of their associated diseases. Beyond point mutations, insertions, deletions and recombination, we here describe for the first time intragenomic rearrangements and their relevance to changes in transmissibility, immune escape and/or virulence documented during the SARS-CoV-2 pandemic. 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. Towards a mathematical understanding of colonization resistance in multispecies microbial communities Microbial community composition and dynamics are key to health and disease. Explaining the forces generating and shaping diversity in the microbial consortia making up our bodys defenses is a major aim of current research in microbiology. For this, tractable models are needed, that bridge the gap between observations of patterns and underlying mechanisms. While most microbial dynamics models are based on the Lotka-Volterra framework, we still do not have an analytic quantity for colonization resistance, by which a microbial systems fitness as a whole can be understood. Here, inspired by an epidemiological perspective, we propose a rather general modeling framework whereby colonization resistance can be clearly mathematically defined and studied. In our model, N similar species interact with each other through a co-colonization interaction network encompassing pairwise competition and cooperation, abstractly mirroring how organisms effectively modify their micro-scale environment in relation to others. This formulation relies on a generic notion of shared resources between members of a consortium, yielding explicit frequency-dependent dynamics among N species, in the form of a replicator equation, and offering a precise definition of colonization resistance. We demonstrate that colonization resistance arises and evolves naturally in a multispecies system as a collective quadratic term in a replicator equation, describing dynamic mean invasion fitness. Each pairwise invasion growth rate between two ecological partners, [Formula], is derived explicitly from species asymmetries and mean traits. This makes the systemic colonization resistance [Formula] also an emergent function of global mean-field parameters and trait variation architecture, weighted by the evolving relative abundances among species. In particular, if the underlying invasion fitness matrix {Lambda} displays species-specific  invasiveness or  invasibility structure, colonization resistance will be insensitive to mean micro-scale cooperation or competition. However, in general, colonization resistance depends on and may undergo critical transitions with changes in mean  environment, e.g. cooperation and growth level in a community. We illustrate several key links between our proposed measure of colonization resistance and invader success, including sensitivity to timing, and to the intrinsic pairwise invasion architecture of the resident community. Our simulations reveal that symmetric and invader-driven mutual invasion among resident species tend to maximize systemic colonization resistance to outsiders, when compared to resident-driven, anti-symmetric, almost anti-symmetric and random {Lambda} structures. We contend this modeling approach is a powerful new avenue to study, test and validate interaction networks and invasion topologies in diverse microbial consortia, and quantify analytically their role in colonization resistance, system function, and invasibility. Effects of substrate on shrimp growth, water quality and bacterial community in the biofloc system nursing Penaeus vannamei postlarvae This study aimed to investigate the effects of substrate on water quality, shrimp growth and bacterial community in the biofloc system with a salinity of 5{per thousand}. Two treatments, biofloc system with (sB) or without (nB) addition of elastic solid packing filler (nylon) as substrate, were set up. Penaeus vannamei postlarvae (PL, [~] stage 15) were stocked at a density of 4000 PL m-3 in each treatment with triplicates for a 28-days culture experiment, taking glucose as carbon source (C:N 15:1). Results showed that the survival rate (96.3{+/-}3.6%), FCR (0.76{+/-}0.06) and productivity (1.54{+/-}0.12 kg m-3) in sB treatment were significantly better than those in nB treatment (81.0{+/-}7.1%, 0.98{+/-}0.08 and 1.14{+/-}0.09 kg m-3, P<0.05). All water parameters were in the recommended ranges. Substrate showed significant effect on TAN, TSS, turbidity, biofloc volume, pH and carbonate alkalinity (P < 0.05). Actinobacteria (4.0-22.7%), Bacteroidetes (10.4-33.5%), Firmicutes (0.2-11.2%), Planctomycetes (4.0-14.9%) and Proteobacteria (29.4-59.0%) were the most dominant phyla for both treatments. However, the bacterial community in sB treatment showed to be significantly different from that in nB treatment (Jaccard distance 0.94{+/-}0.01, P=0.001). Substrate showed significant effects on Shannon, Heip, Pielou and Simpson index, as well as relative abundances of Actinobacteria, Bacteroidetes and Planctomycetes (P < 0.05). The results suggested that addition of substrate affected the shrimp growth, water quality and bacterial community in the biofloc system nursing P. vannamei PL with a 5{per thousand} salinity. 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. Coupled yawing and orbital rotation in kinesin motility revealed by polarization measurement of gold nanorods Kinesin motor domains generate impulses of force and movement that have both translational and rotational components, raising the question of how the rotational component contributes to motor functions. We used a new assay in which kinesin-coated gold nanorods (kinesin-GNRs) move on suspended microtubules, for two plus-end-directed and weakly-processive kinesins: single-headed KIF1A, dimeric ZEN-4. Polarization of the light scattered by two types of kinesin-GNRs periodically oscillated as they orbited the microtubule along a left-handed helical trajectory. Our analyses revealed that each kinesin-GNR unidirectionally rotates about its yaw axis as it translocates, and that the period of this yaw-axis rotation corresponds to two periods of its left-handed helical orbit around the microtubule axis. Theoretical analyses suggest that the yaw-axis rotation enhances biased lateral displacement of the kinesin team. Our study reveals biaxial rotation as a new mode of motility in kinesin teams that may help the team to sidestep roadblocks. A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Applied tutorial for the design and fabrication of biomicrofluidic devices by resin 3D printing Stereolithographic (SL) 3D printing, especially digital light processing (DLP) printing, is a promising rapid fabrication method for bio-microfluidic applications such as clinical tests, lab-on-a-chip devices, and sensor integrated devices. The benefits of 3D printing lead many to believe this fabrication method will accelerate the use of microfluidics, but there are a number of potential obstacles to overcome for bioanalytical labs to fully utilize this technology. For commercially available printing materials, this includes challenges in producing prints with the print resolution and mechanical stability required for a particular design, along with cytotoxic components within many SL resins and low optical compatibility for imaging experiments. Potential solutions to these problems are scattered throughout the literature and rarely available in head-to-head comparisons. Therefore, we present here a concise guide to the principles of resin 3D printing most relevant for fabrication of bioanalytical microfluidic devices. Intended to quickly orient labs that are new to 3D printing, the tutorial includes the results of selected systematic tests to inform resin selection, strategies for design optimization, and improvement of biocompatibility of resin 3D printed bio-microfluidic devices. The nucleus bypasses obstacles by deforming like a drop with surface tension mediated by lamin A/C Migrating cells must deform their stiff cell nucleus to move through pores and fibers in tissue. Lamin A/C is known to hinder cell migration by limiting nuclear deformation and passage through confining channels, but its role in nuclear deformation and passage through fibrous environments is less clear. We studied cell and nuclear migration through discrete, closely spaced, slender obstacles which mimic the mechanical properties of collagen fibers. Nuclei bypassed slender obstacles while preserving their overall morphology by deforming around them with deep local invaginations of little resisting force. The obstacles did not impede the nuclear trajectory or cause a rupture of the nuclear envelope. Nuclei likewise deformed around single collagen fibers in cells migrating in 3D collagen gels. In contrast to its limiting role in nuclear passage through confining channels, lamin A/C facilitated nuclear deformation and passage through fibrous environments; nuclei in lamin-null (Lmna-/-) cells lost their overall morphology and became entangled on the obstacles. Analogous to surface tension-mediated deformation of a liquid drop, lamin A/C imparts a surface tension on the nucleus that allows nuclear invaginations with little mechanical resistance, preventing nuclear entanglement and allowing nuclear passage through fibrous environments. Native American Ancestry and Pigmentation Allele Contributions to Skin Color in a Caribbean Population Interest in the genetic basis of variation in skin pigmentation in Native American populations led us to seek indigenous populations of the Western Hemisphere with African and minimal European admixture to study the effect of Native American ancestry on skin color. Admixture analysis from DNA collected from 458 individuals in the Kalinago territory of the Commonwealth of Dominica showed shared ancestry with East Asians at K=3 and 55% Native American, 32% African, and 11% European ancestry at K=6, the highest Native American ancestry of Caribbean populations. Skin pigmentation was 20 to 80 melanin units, averaging 46. Three albino individuals were homozygous for multi-nucleotide polymorphism OCA2NW273KV of African origin, whose population allele frequency was 0.03 and single allele effect size was -8 melanin units. Hypopigmenting allele frequencies for SLC24A5A111T and SLC45A2L374F were 0.14 and 0.05, whose single allele effect sizes were -6 and -3, respectively. Skin color plots of individuals lacking known hypopigmenting alleles suggests that Native American Ancestry reduced pigmentation by more than 20 melanin units (low and high estimates 21.8 and 28.5). Shared ancestry with East Asians at K=3 suggests potential sharing of one or more pigmentation alleles. Muscle Activation Profile During Perturbed Walking is Modulated According to Body State During an unexpected loss of balance, avoiding a fall requires people to readjust their footing rapidly and effectively. We investigated the activation patterns of the ankle and knee muscles, and muscle fiber type recruitment resulting from unannounced, medio-lateral (i.e., right/left) horizontal-surface walking perturbations in twenty healthy adults (27.00{+/-}2.79 years, 10 females). Surface electromyography (sEMG) total spectral power for specific frequency bands (40-60Hz, 60-150Hz, 150-250Hz, 250-400Hz and 400-1000Hz), from tibialis anterior (TA) and vastus lateralis (VL) muscles were analyzed. Compared to non-perturbed walking, we found a significant increase in the total spectral power of lower-extremity muscles during the first 3 seconds after perturbation. When two feet were on the ground in time of perturbation we found a different muscle fiber type recruitment pattern between VL and TA muscles. This was not significant for perturbations implemented when one foot was on the ground. Our findings suggest that muscle operating frequency is modulated in real time to fit body state and somatosensory input in the context of functional goal requirements such as a rapid change of footing in response to unexpected loss of balance in single and double-support phases of gait.

New & NoteworthyTo study muscle spectral profiles in response to loss of balance, we investigated the dynamics of muscle spectral power changes, across different frequency bands after unannounced mechanical perturbations during walking. We showed increased activation of high-frequency motor units of the lower-limb muscles, subside 3 seconds after perturbation. Differences in power increase of specific frequency bands suggest that muscle activation is modulated in real time to fit body state in the context of functional goal requirements. Vesalius: high-resolution in silico anatomization of Spatial Transcriptomic data using Image Analysis Characterization of tissue architecture promises to deliver insights into development, cell communication and disease. In silico spatial domain retrieval methods have been developed for spatial transcriptomics (ST) data assuming transcriptional similarity of neighboring barcodes. However, domain retrieval approaches with this assumption cannot work in complex tissues composed of multiple cell types. This task becomes especially challenging in cellular resolution ST methods. We developed Vesalius to decipher tissue anatomy from ST data by applying image processing technology. Vesalius uniquely detected territories composed of multiple cell types and successfully recovered tissue structures in high-resolution ST data including in mouse brain, embryo, liver, and colon. Utilizing this tissue architecture, Vesalius identified tissue morphology specific gene expression and regional specific gene expression changes for astrocytes, interneuron, oligodendrocytes, and entorhinal cells in the mouse brain. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. A genomic and morphometric analysis of alpine bumblebees: ongoing reductions in tongue length but no clear genetic component Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyze population structure in these populations. We generated a genome assembly of B. balteatus. We then analyzed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species ("incognitus"). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait likely has a complex genetic and environmental basis.