The Temporal Pattern of Synaptic Activation Determines the Longevity of Structural Plasticity at Single Dendritic Spines Individual synapses are the points at which information is passed between neurons, yet it is unknown how the diverse patterns of activity that are observed in vivo effect plasticity at the level of single inputs. Here, we aimed to determine what are the structural plasticity consequences of naturalistic patterns of activity at single spines, as these reflect changes in synaptic efficacy. Utilizing two- photon fluorescence imaging and glutamate uncaging, we studied structural plasticity of individual CA1 hippocampal dendritic spines using activation patterns sampled from a Poisson distribution, which resemble endogenous firing patterns from their CA3 inputs. We found that while the majority of inputs initially undergo structural changes, the robustness of this plasticity is determined by the timing structure of the Poisson sampled naturalistic stimulation patterns. Further, we found that structural plasticity elicited by these naturalistic patterns is both NMDAR and protein synthesis dependent, consistent with requirements for other forms of plasticity. Lastly, we found that during the delivery of naturalistic activity patterns, spines underwent rapid and dynamic structural growth that predicted the longevity of plasticity, which was not the case during non-naturalistic stimulation protocols. These data suggest that dendritic spines are able to integrate incoming temporal information and accordingly modulate the longevity of plasticity that is induced.

HighlightsO_LINaturalistic stimulation of single dendritic spines of CA1 hippocampal neurons induces long lasting structural plasticity that depends on the temporal distribution of the synaptic events.
C_LIO_LIStructural plasticity induced by naturalistic stimulation patterns requires NMDA receptor activation and new protein-synthesis.
C_LIO_LIRapid spine structural dynamics during naturalistic activity, but not regular patterns, predict the longevity of subsequent structural plasticity.
C_LI
Understanding transformation tolerant visual object representations in the human brain and convolutional neural networks Forming transformation-tolerant object representations is critical to high-level primate vision. Despite its significance, many details of tolerance in the human brain remain unknown. Likewise, despite the ability of convolutional neural networks (CNNs) to exhibit human-like object categorization performance, whether CNNs form tolerance similar to that of the human brain is unknown. Here we provide the first comprehensive documentation and comparison of three tolerance measures in the human brain and CNNs. We measured fMRI responses from human ventral visual areas to real-world objects across both Euclidean and non-Euclidean feature changes. In single fMRI voxels in higher visual areas, we observed robust object response rank-order preservation across feature changes. This is indicative of functional smoothness in tolerance at the fMRI meso-scale level that has never been reported before. At the voxel population level, we found highly consistent object representational structure across feature changes towards the end of ventral processing. Rank-order preservation, consistency, and a third tolerance measure, cross-decoding success (i.e., a linear classifiers ability to generalize performance across feature changes) showed an overall tight coupling. These tolerance measures were lower for Euclidean than non-Euclidean feature changes in lower visual areas, but increased over the course of ventral processing in most cases. These characteristics of tolerance, however, were absent in eight CNNs pretrained with ImageNet images with varying network architecture, depth, the presence/absence of recurrent processing, or whether a network was pretrained with the original or stylized ImageNet images that encouraged shape processing. Most notably, CNNs do not show increased representational consistency across feature changes at the higher layers. CNNs thus do not appear to develop the same kind of tolerance as the human brain over the course of visual processing.

Significant StatementPerceiving object identity among changes in non-identity features and forming transformation-tolerant object representations is essential to high-level primate vision. Here we provide a comprehensive documentation and comparison of three tolerance measures between the human brain and CNNs pretrained for object classification. While all three measures show increased tolerance in the human brain across four types of feature changes towards the end of ventral visual processing, CNNs fail to develop the same kind of tolerance with visual processing.
Modulations in motor unit discharge are related to changes in fascicle length during isometric contractions The integration of electromyography (EMG) and ultrasound imaging has provided important information about the mechanisms of muscle activation and contraction. Unfortunately, EMG does not allow an accurate assessment of the interplay between the neural drive received by muscles, changes in fascicle length (FL) and the force/torque produced. We aimed to assess the relationship between modulations in tibialis anterior (TA) motor unit (MU) firing rate, FL and dorsiflexion torque (DT) using ultrasound-transparent high-density EMG electrodes. EMG and ultrasound images were recorded simultaneously from TA, using a 32-electrode silicon matrix, while performing isometric dorsiflexion, at diverse ankle joint positions (0{degrees} and 30{degrees} plantar flexion) and torques (20% and 40% of maximum). EMG signals were decomposed into individual MUs and changes in FL were assessed with a fascicle-tracking algorithm. MU firings were converted into a cumulative spike train (CST) that was cross-correlated with DT (CST-DT) and FL (CST-FL). High cross-correlations were found for CST-FL, 0.60 (range: 0.31-0.85) and CST-DT 0.71 (range: 0.31-0.88). Cross-correlation lags revealed that the delay between CST-FL (~75ms) was significantly smaller than CST-DT (~150ms, p<0.001). These delays affected the interpretation of MU recruitment/de-recruitment thresholds, with FL showing similar lengths for both recruitment and de-recruitment. This study is the first to demonstrate that changes in TA FL are closely related to both modulations in MU firing frequency and DT. These relationships allow assessment of the interplay between neural drive, muscle contraction and resultant torque, thereby providing a better understanding of the mechanisms responsible for the generation of muscle force.

NEW AND NOTEWORTHYBy employing ultrasound-transparent high-density surface EMG electrodes, we show that modulations in tibialis anterior motor unit discharge rate were closely related to both changes in its fascicle length and resultant torque. These relationships allowed quantifying delays between neural drive and muscle shortening as well as muscle shortening and torque during submaximal isometric contractions, providing an accurate estimation of the time required to generate muscle force and subsequent production of torque via the tendon. Human genetic analyses of organelles highlight the nucleus, but not the mitochondrion, in age-related trait heritability Most age-related human diseases are accompanied by a decline in cellular organelle integrity, including impaired lysosomal proteostasis and defective mitochondrial oxidative phosphorylation. An open question, however, is the degree to which inherited variation in or near genes encoding each organelle contributes to age-related disease pathogenesis. Here, we evaluate if genetic loci encoding organelle proteomes confer greater-than-expected age-related disease risk. As mitochondrial dysfunction is a "hallmark" of aging, we begin by assessing nuclear and mitochondrial DNA loci near genes encoding the mitochondrial proteome and surprisingly observe a lack of enrichment across 24 age-related traits. Within nine other organelles, we find no enrichment with one exception: the nucleus, where enrichment emanates from nuclear transcription factors. In agreement, we find that genes encoding several organelles tend to be "haplosufficient," while we observe strong purifying selection against heterozygous protein-truncating variants impacting the nucleus. Our work identifies common variation near transcription factors as having outsize influence on age-related trait risk, motivating future efforts to determine if and how this inherited variation then contributes to observed age-related organelle deterioration. Variation in white spruce needle respiration at the species range limits: a potential impediment to northern expansion White spruce (Picea glauca) spans a massive range from arctic treeline to temperate forests, yet the variability in respiratory physiology and related implications for tree carbon balance at the extremes of this distribution remain as enigmas. Working at both the most northern and southern extents of the white spruce distribution range more than 5000 km apart, we measured the short- term temperature response of dark respiration (R/T) at upper and lower canopy positions. R/T curves were fit to both polynomial and thermodynamic models so that model parameters could be compared among locations, canopy positions, and with previously published data. Respiration measured at 25{degrees}C (R25) was 68% lower at the southern location than at the northern location (0.73{+/-}0.15 vs. 2.27{+/-}0.02 mol m-2 s-1), resulting in a significantly lower (p< 0.01) intercept in R/T response in temperate trees. Only at the southern location did upper canopy leaves have a steeper temperature response than lower canopy leaves, likely reflecting steeper canopy gradients in light. No differences were observed in the maximum temperature of respiration. At the northern range limit respiration is nearly twice that of the average R25 reported in a global leaf respiration database. This large carbon cost likely challenges tree survival and contributes to restricting the location of the northern treeline. We predict that without significant thermal acclimation, foliage respiration will increase with projected end-of-the-century warming and will likely constrain the future range limits of this important boreal species.

Summary StatementWhite spruce (Picea glauca) needle respiration at the northern limit of the species range is three times higher than at the southern range limit (when measured at 25 {degrees}C). This high carbon cost likely challenges tree survival and contributes to the location of the northern treeline. HuR-dependent SOD2 protein synthesis is an early adaptation to anchorage independence During metastasis, cancer cells must adapt to survive a loss of anchorage and evade anoikis. An important pro-survival adaptation is the ability of metastatic tumor cells to increase their antioxidant capacity and restore cellular redox balance. Although much is known about the transcriptional regulation of antioxidant enzymes in response to stress, how cells rapidly adapt to alter antioxidant enzyme levels is less well understood. Using ovarian cancer cells as a model, we demonstrate that an increase in protein expression of the mitochondrial superoxide dismutase SOD2 is a very early event initiated in response to cellular detachment. SOD2 protein synthesis is rapidly induced within 0.5-2 hours of matrix detachment, and polyribosome profiling demonstrates an increase in the number of ribosomes bound to SOD2 mRNA, indicating an increase in SOD2 translation in response to anchorage-independence. Mechanistically, we find that anchorage-independence specifically induces cytosolic accumulation of the RNA binding protein HuR/ELAVL1 and leads to increased HuR binding to SOD2 mRNA. Using HuR siRNA-mediated knock-down, we show that the presence of HuR is necessary for the increase in SOD2 mRNA association with the heavy polyribosome fraction and SOD2 protein synthesis observed in anchorage-independence. Cellular detachment activates the stress-response protein kinase p38 MAPK, which is necessary for HuR-SOD2 mRNA binding and optimal increases in SOD2 protein expression. These findings illustrate a novel post-transcriptional regulatory mechanisms of SOD2, enabling cells to rapidly increase their mitochondrial antioxidant capacity as an acute response to anchorage-independence. Deep Heterogeneous Dilation of LSTM for Transient-phase Gesture Prediction through High-density Electromyography: Towards Application in Neurorobotics Deep networks have been recently proposed to estimate motor intention using conventional bipolar surface electromyography (sEMG) signals for myoelectric control of neurorobots. In this regard, Deepnets are generally challenged by long training times (affecting practicality and calibration), complex model architectures (affecting the predictability of the outcomes), and a large number of trainable parameters (increasing the need for big data). Capitalizing on our recent work on homogeneous temporal dilation in a Recurrent Neural Network (RNN) model, this paper proposes, for the first time, heterogeneous temporal dilation in an LSTM model and applies that to high-density surface electromyography (HD-sEMG), allowing for the decoding of dynamic temporal dependencies with tunable temporal foci. In this paper, a 128-channel HD-sEMG signal space is considered due to the potential for enhancing the spatiotemporal resolution of human-robot interfaces. Accordingly, this paper addresses a challenging motor intention decoding problem of neurorobots, namely, transient intention identification. Our approach uses only the dynamic and transient phase of gesture movements when the signals are not stabilized or plateaued, which can significantly enhance the temporal resolution of human-robot interfaces. This would eventually enhance seamless real-time implementations. Additionally, this paper introduces the concept of "dilation foci" to modulate the modeling of temporal variation in transient phases. In this work a high number (e.g., 65) of gestures is included, which adds to the complexity and significance of the understudied problem. Our results show state-of-the-art performance for gesture prediction in terms of accuracy, training time, and model convergence. Social integration predicts survival in female white-faced capuchin monkeys Across multiple species of social mammals, a growing number of studies have found that that individual sociality is associated with survival. In long-lived species, like primates, lifespan is one of the main components of fitness. We used 18 years of data from the Lomas Barbudal Monkey project to quantify social integration in 11 capuchin (Cebus capucinus) groups and tested whether female survivorship was associated with females tendencies to interact with three types of partners: (1) all group members, (2) adult females, (3) adult males. We found strong evidence that females who engaged more with other females in affiliative interactions and foraged in close proximity experienced increased survivorship. We found some weak evidence that females might also benefit from engaging in more support in agonistic contexts with other females. These benefits were evident in models that account for the females rank and group size. Female interactions with all group members also increased survival, but the estimates of the effects were more uncertain. In interactions with adult males, only females who provided more grooming to males survived longer. The results presented here suggest that social integration may result in survival-related benefits. Females might enjoy these benefits through exchanging grooming for other currencies, such as coalitionary support or tolerance. Host species identity shapes the diversity and structure of insect microbiota As for most of the life that inhabits our planet, microorganisms play an essential role in the fitness of insects, including nutrition, reproduction, defence, and many other functions. More recently, we assisted in an exponential growth of studies describing the taxonomical composition of bacterial communities across insects phylogeny. However, there is still an outstanding question that needs to be answered: which factors contribute most to shape insects microbiomes? This study tries to find an answer to this question by taking advantage of publicly available sequencing data and reanalysing over 4,000 samples of insect-associated bacterial communities under a common framework. Results suggest that insect taxonomy has a wider impact on the structure and diversity of their associated microbial communities than the other factors considered (diet, sex, life stage, sample origin, and treatment). However, when specifically testing for signatures of co-diversification of insect species and their microbiota, analyses found weak support for this, suggesting that while insect species strongly drive the structure and diversity of insect microbiota, the diversification of those microbial communities did not follow their hosts phylogeny. Furthermore, a parallel survey of the literature highlights several methodological limitations that need to be considered in future research endeavours. Sperm Histone H3 Lysine 4 tri-methylation serves as a metabolic sensor of paternal obesity and is associated with the inheritance of metabolic dysfunction ObjectiveParental environmental exposures can strongly influence descendant risks for adult disease. How paternal obesity changes the sperm chromatin leading to the acquisition of metabolic disease in offspring remains controversial and ill-defined. The objective of this study was to assess: (1) whether obesity induced by a high-fat diet alters sperm histone methylation; (2) whether paternal obesity can induce metabolic disturbances across generations; (3) whether there could be cumulative damage to the sperm epigenome leading to enhanced metabolic dysfunction in descendants; and (4) whether obesity-sensitive regions associate with embryonic epigenetic and transcriptomic profiles. Using a genetic mouse model of epigenetic inheritance, we investigated the role of histone H3 lysine 4 methylation (H3K4me3) in the paternal transmission of metabolic dysfunction. This transgenic mouse overexpresses the histone demethylase enzyme KDM1A in the developing germline and has an altered sperm epigenome at the level of histone H3K4 methylation. We hypothesized that challenging transgenic sires with a high-fat diet would further erode the sperm epigenome and lead to enhanced metabolic disturbances in the next generations.

MethodsTo assess whether paternal obesity can have inter- or transgenerational impacts, and if so, to identify potential mechanisms of this non-genetic inheritance, we used wildtype C57BL/6NCrl and transgenic males with a pre-existing altered sperm epigenome. To induce obesity, sires were fed either a control or high-fat diet (10% or 60% kcal fat, respectively) for 10-12 weeks, then bred to wildtype C57BL/6NCrl female fed a regular diet. F1 and F2 descendants were characterized for metabolic phenotypes by examining the effects of paternal obesity by sex, on body weight, fat mass distribution, the liver transcriptome, intraperitoneal glucose and insulin tolerance tests. To determine whether obesity altered the F0 sperm chromatin, native chromatin immunoprecipitation-sequencing targeting H3K4me3 was performed. To gain insight into mechanisms of paternal transmission, we compared our sperm H3K4me3 profiles with embryonic and placental chromatin states, histone modification and gene expression profiles.

ResultsObesity-induced alterations in H3K4me3 occurred at genes implicated in metabolic, inflammatory, and developmental processes. These processes were associated with offspring metabolic dysfunction and corresponded to genes enriched for H3K4me3 in embryos, and overlapped embryonic and placenta gene expression profiles. Transgenerational susceptibility to metabolic disease was only observed when obese F0 had a pre-existing modified sperm epigenome. This coincided with increased H3K4me3 alterations in sperm and more severe phenotypes affecting their offspring.

ConclusionsOur data suggest sperm H3K4me3 might serve as a metabolic sensor that connects paternal diet with offspring phenotypes via the placenta. This non-DNA based knowledge of inheritance has the potential to improve our understanding of how environment shapes heritability and may lead to novel routes for the prevention of disease. This study highlights the need to further study the connection between the sperm epigenome, placental development and childrens health. The isoforms of pyruvate kinase act as nutrient sensors for the β-cell KATP channel Insulin release from pancreatic {beta}-cells is triggered by the metabolism-dependent closure of ATP-sensitive K+ (KATP) channels. Here, we used isoform-specific deletion of pyruvate kinase (PK) to identify the unique roles of PKm1 and PKm2 in the {beta}-cell nutrient response. In the absence of PK activity, cytosolic ATP/ADP was insufficient to close KATP channels. Both PK isoforms were found on the plasma membrane, where they endow KATP channels with nutrient sensitivity and specificity. Although PKm1 accounts for >90% of cellular activity, PKm2 controls glucosedependent KATP closure. Amino acids, which generate mitochondrial phosphoenolpyruvate (PEP) without allosterically activating PKm2, depend on PKm1 to raise ATP/ADP and close KATP channels. Shifting the {beta}-cell from PKm2 to PKm1 correspondingly increased amino acid sensitivity and boosted secretory output. In addition to these membrane depolarizing "on-switches," {beta}-cell deletion of PCK2 exposed a mitochondrial PEP-dependent "off-switch" that facilitates Ca2+ extrusion. Restricting PCK2 prolonged Ca2+ influx and increased insulin secretion. In summary, it is not the KATP channel alone, but the PK-KATP microcompartment that confers nutrient sensitivity to the {beta}-cell plasma membrane. The differential response of the PK isoforms to the glycolytic and mitochondrial sources of PEP influences {beta}-cell nutrient preference, and controls the oscillatory cycle regulating insulin secretion. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Agonist efficiency from concentration-response curves: structural implications and applications Agonists are evaluated by a concentration-response curve (CRC), with a midpoint (EC50) that indicates potency, a high-concentration asymptote that indicates efficacy and a low-concentration asymptote that indicates constitutive activity. A third agonist attribute, efficiency ({eta}), is the fraction of binding energy that is applied to the conformational change that activates the receptor. We show that {eta} can be calculated from EC50 and the asymptotes of a CRC derived from either single-channel or whole-cell responses. For 20 agonists of skeletal muscle nicotinic receptors, the distribution of {eta} values is bimodal with population means at 51% (including ACh, nornicotine and DMPP) and 40% (including epibatidine, varenicline and cytisine). The value of {eta} is related inversely to the size of the agonists head-group, with high-versus low-efficiency ligands having an average volume of 70 [A]3 versus 102 [A]3. Most binding site mutations have only a small effect on ACh efficiency except for Y190A (35%), W149A (60%) and those at G153 (42%). If {eta} is known, the midpoint and high-concentration asymptote can be calculated from each other. Hence, an entire CRC can be estimated from the response to a single agonist concentration, and efficacy can be estimated from EC50 of a CRC that has been normalized to 1. Given {eta}, the level of constitutive activity can be estimated from a single CRC.

Statement of significanceReceptors are molecular machines that convert chemical energy from agonist binding into mechanical energy of a global conformational change that generates a cell response. Agonists are distinguished by their potency (proportional to affinity) and efficacy, but also by the efficiency at which their binding energy is applied to receptor activation. Here we show that agonist efficiency can be estimated from a single concentrationresponse curve, and estimate efficiencies of 20 nicotinic receptor agonists. These have a bimodal distribution with larger agonists belonging to the lower-efficiency population. We further show that mutations of some binding site residues alter efficiency, and that knowledge of agonist efficiency simplifies and extends dose-response curve analysis. Age-associated B cells are long-lasting effectors that restrain reactivation of latent γHV68 Age-associated B cells (ABCs; CD19+CD11c+T-bet+) are a unique population that are increased in an array of viral infections, though their role during latent infection is largely unexplored. Here, we use murine gammaherpesvirus 68 ({gamma}HV68) to demonstrate that ABCs remain elevated long-term during latent infection and express IFN{gamma} and TNF. Using a strain of {gamma}HV68 that is cleared following acute infection, we show that ABCs persist in the absence of latent virus, though their expression of IFN{gamma} and TNF is decreased. With a fluorescent virus we demonstrate that ABCs are infected with {gamma}HV68 at similar rates to other previously activated B cells. We find that mice without ABCs display defects in anti-viral IgG2a/c antibodies and are more susceptible to {gamma}HV68 reactivation when challenged with heterologous infection. Together, these results indicate that ABCs are a persistent effector subset during latent viral infection that restrains {gamma}HV68 reactivation. Infection and transmission of ancestral SARS-CoV-2 and its alpha variant in pregnant white-tailed deer SARS-CoV-2, a novel Betacoronavirus, was first reported circulating in human populations in December 2019 and has since become a global pandemic. Recent history involving SARS-like coronavirus outbreaks (SARS-CoV and MERS-CoV) have demonstrated the significant role of intermediate and reservoir hosts in viral maintenance and transmission cycles. Evidence of SARS-CoV-2 natural infection and experimental infections of a wide variety of animal species has been demonstrated, and in silico and in vitro studies have indicated that deer are susceptible to SARS-CoV-2 infection. White-tailed deer (Odocoileus virginianus) are amongst the most abundant, densely populated, and geographically widespread wild ruminant species in the United States. Human interaction with white-tailed deer has resulted in the occurrence of disease in human populations in the past. Recently, white-tailed deer fawns were shown to be susceptible to SARS-CoV-2. In the present study, we investigated the susceptibility and transmission of SARS-CoV-2 in adult white-tailed deer. In addition, we examined the competition of two SARS-CoV-2 isolates, representatives of the ancestral lineage A (SARS-CoV-2/human/USA/WA1/2020) and the alpha variant of concern (VOC) B.1.1.7 (SARS-CoV-2/human/USA/CA_CDC_5574/2020), through co-infection of white-tailed deer. Next-generation sequencing was used to determine the presence and transmission of each strain in the co-infected and contact sentinel animals. Our results demonstrate that adult white-tailed deer are highly susceptible to SARS-CoV-2 infection and can transmit the virus through direct contact as well as vertically from doe to fetus. Additionally, we determined that the alpha VOC B.1.1.7 isolate of SARS-CoV-2 outcompetes the ancestral lineage A isolate in white-tailed deer, as demonstrated by the genome of the virus shed from nasal and oral cavities from principal infected and contact animals, and from virus present in tissues of principal infected deer, fetuses and contact animals. Neuroinvasion and neurotropism by SARS-CoV-2 variants in the K18-hACE2 mouse Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) not only affects the respiratory tract but also causes neurological symptoms such as loss of smell and taste, headache, fatigue or severe cerebrovascular complications. Using transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) we investigated the spatiotemporal distribution and pathomorphological features in the CNS following intranasal infection with SARS-CoV-2 variants, also after prior influenza A virus infection. Apart from Omicron, we found all variants to frequently spread to and within the CNS. Infection was restricted to neurons and appeared to spread from the olfactory bulb mainly in basally orientated regions in the brain and into the spinal cord, independent of ACE2 expression and without evidence of neuronal cell death, axonal damage or demyelination. However, microglial activation, microgliosis and a mild macrophage and T cell dominated inflammatory response was consistently observed, accompanied by apoptotic death of endothelial, microglial and immune cells, without their apparent infection. Microgliosis and immune cell apoptosis indicate a potential role of microglia for pathogenesis and viral effect in COVID-19 and possible impairment of neurological functions, especially in long COVID. These data may also be informative for the selection of therapeutic candidates, and broadly support investigation of agents with adequate penetration into relevant regions of the CNS. Unbiased integration of single cell transcriptome replicates Single cell transcriptomic approaches are becoming mainstream, with replicate experiments commonly performed with the same single cell technology. Methods that enable integration of these datasets by removing batch effects while preserving biological information are required for unbiased data interpretation. Here we introduce Canek for this purpose. Canek leverages information from mutual nearest neighbor to combine local linear corrections with cell-specific non-linear corrections within a fuzzy logic framework. Using a combination of real and synthetic datasets, we show that Canek corrects batch effects while introducing the least amount of bias compared with competing methods. Canek is computationally efficient and can easily integrate thousands of single-cell transcriptomes from replicated experiments. Trade-offs in biosensor optimization for dynamic pathway engineering Recent progress in synthetic biology allows the construction of dynamic control circuits for metabolic engineering. This technology promises to overcome many challenges encountered in traditional pathway engineering, thanks to their ability to self-regulate gene expression in response to bioreactor perturbations. The central components in these control circuits are metabolite biosensors that read out pathway signals and actuate enzyme expression. However, the construction of metabolite biosensors is a major bottleneck for strain design, and a key challenge is to understand the relation between biosensor dose-response curves and pathway performance. Here we employ multiobjective optimization to quantify performance trade-offs that arise in the design and calibration of metabolite biosensors. Our approach reveals strategies for tuning dose-response curves along an optimal trade-off between production flux and the cost of an increased expression burden on the host. We explore properties of control architectures built in the literature, and identify their advantages and caveats in terms of performance and robustness to growth conditions and leaky promoters. We demonstrate the optimality of a control circuit for glucaric acid production in Escherichia coli, which has been shown to increase titer by 2.5-fold as compared to static designs. Our results lay the groundwork for the automated design of control circuits for pathway engineering, with applications in the food, energy and pharmaceutical sectors. Correlational networking guides the discovery of cryptic natural product biosynthetic enzymes Bacterial natural product biosynthetic genes, canonically clustered, have been increasingly found to rely on cryptic enzymes encoded elsewhere in the genome for completion of biosynthesis. The study and application of lanthipeptides are frequently hindered by cryptic proteases required for final maturation. Here, we establish a global correlation network bridging the gap between lanthipeptide precursors and cryptic proteases. Applying our analysis to 161,954 bacterial genomes, we establish 5,209 correlations between precursors and cryptic proteases, with 91 prioritized. We use network predictions and co-expression analysis to reveal a previously missing protease for the maturation of class I lanthipeptide paenilan. We further discover widely distributed bacterial M16B metallopeptidases of previously unclear biological function as a new family of lanthipeptide proteases. We show the involvement of a pair of bifunctional M16B proteases in the production of novel class III lanthipeptides with high substrate specificity. Together, these results demonstrate the strength of our correlational networking approach to the discovery of cryptic lanthipeptide proteases and potentially other missing enzymes for natural products biosynthesis.

Graphical abstract O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY Herbivore life histories are altered by drought stress in their hosts plants Drought associated with climate change can stress plants, altering their interactions with phytophagous arthropods. Dryness not only impacts cultivated plants but also their parasites, which in some cases are favoured by drought. Herbivorous arthropods feeding on drought-stressed plants typically produce bigger offspring and develop faster. However, how much responses to abiotic factors differ among populations of a species remains poorly documented. Here, we evaluate variability among populations of a major agricultural pest, the two spotted spider mite, Tetranychus urticae, in response to drought-stress. We compare key life-history parameters of twelve populations that originate from climates ranging from wet and cool Atlantic locations to medium to dry hot Mediterranean locations. We evaluated how plant drought stress affects four life history traits: development time, fecundity, sex-ratio and emigration rate in an experiment comparing well-watered and drought-stressed bean plants. Mites feeding on drought-stressed plants developed faster and attempted to leave leaves less often, and young females were more fecund. The mites from wet temperate climates exhibited greater plasticity between the two water regimes than mites originating from dryer and hot climates. Thus, climate in the area of origin influences mite response to drought. 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. Plasma proteome analyses in individuals of European and African ancestry identify cis-pQTLs and models for proteome-wide association studies Improved understanding of genetic regulation of proteome can facilitate the identification of causal mechanisms for complex traits. We analyzed data on 4,657 plasma proteins from 7,213 European American (EA) and 1,871 African American (AA) individuals from the ARIC study, and further replicated findings on 467 AA individuals from the AASK study. Here we identified 2,004 proteins in EA and 1,618 in AA, with majority overlapping, which showed associations with common variants in cis-regions. Availability of AA samples led to smaller credible sets and significant number of population-specific cis-pQTLs. Elastic-net produced powerful models for protein prediction in both populations. An application of proteome-wide association studies (PWAS) to serum urate and gout, implicated several proteins, including IL1RN, revealing the promise of the drug anakinra to treat acute gout flares. Our study demonstrates the value of large and diverse ancestry study for genetic mechanisms of molecular phenotypes and their relationship with complex traits. SufB intein splicing in Mycobacterium tuberculosis is influenced by two remote conserved N-extein Histidines Inteins are auto-processing domains that implement a multi-step biochemical reaction termed protein splicing, marked by cleavage and formation of peptide bonds. They excise from a precursor protein, generating a functional protein via covalent bonding of flanking exteins. We report the kinetic study of splicing and cleavage reaction in [Fe-S] cluster assembly protein SufB from Mycobacterium tuberculosis. Although it follows a canonical intein splicing pathway, distinct features are added by extein residues present in the active site. Sequence analysis identified two conserved histidines in the N-extein region; His-5 and His-38. Kinetic analyses of His-5Ala and His-38Ala SufB mutants exhibited significant reductions in splicing and cleavage rates relative to the SufB wild-type precursor protein. Structural analysis and molecular dynamics simulations suggested that Mtu SufB displays a unique mechanism where two remote histidines work concurrently to facilitate N-terminal cleavage reaction. His-38 is stabilized by the solvent-exposed His-5, and can impact N-S acyl shift by direct interaction with the catalytic Cys1. Development of inteins as biotechnological tools or as pathogen specific novel antimicrobial targets requires a more complete understanding of such unexpected roles of conserved extein residues in protein splicing. A computational modeling approach for predicting multicell patterns based on signaling-induced differential adhesion Differential adhesion within cell populations enables the emergence of unique patterns in heterogeneous multicellular systems. This process has been previously explored using synthetically engineered heterogenous multicell spheroid systems, in which cell subpopulations engage in bidirectional intercellular signaling to regulate the expression of different cadherins. While engineered cell systems provide excellent experimental tools to observe pattern formation in cell populations, computational models may be leveraged to more systematically explore the key parameters that drive the emergence of different patterns. We developed and validated two- and three-dimensional agent-based models (ABMs) of spheroid patterning for cells engineered with a bidirectional signaling circuit that regulates N- and P-cadherin expression. The model was used to predict how varying initial cell seedings, cadherin induction probabilities, or homotypic adhesion strengths between cells impact spheroid patterning, and unsupervised machine learning techniques were used to map system parameters to unique spheroid patterns. The model was then deployed to design new synthetic cell signaling circuits based on a desired final multicell pattern. Structure-Functional-Selectivity Relationship Studies on A-86929 Analogs and Small biaryl Fragments toward Discovery of biased D1 agonists Dopamine regulates normal functions such as movement, reinforcement learning, and cognition, and its dysfunction has been implicated in multiple psychiatric and neurological disorders. Dopamine acts through the D1- (D1R and D5R) and D2-class (D2R, D3R and D4R) of seven transmembrane receptors, and activates both G-protein- and {beta}-arrestin-dependent signaling pathways, to mediate its physiological effects. Current dopamine receptor-based therapies are used to ameliorate motor deficits in Parkinsons disease, or as antipsychotic medications for schizophrenia. These drugs show efficacy for ameliorating only some symptoms caused by dopamine dysfunction and are plagued by debilitating side-effects. Studies in primates and rodents have shown that shifting the balance of dopamine receptor signaling towards the arrestin pathway can be beneficial for inducing normal movement, while reducing motor side-effects such as dyskinesias, and can be efficacious at enhancing cognitive function compared to balanced agonists. Several structure-activity-relationship (SAR) studies have embarked on discovering {beta}-arrestin-biased dopamine agonists, focused on D2 partial agonists, non-catechol D1 agonists, and mixed D1/D2R dopamine receptor agonists. Here, we describe an SAR study to identify novel D1R {beta}-arrestin biased ligands using A-86929, a high-affinity D1R catechol agonist, as a core scaffold. Previously described and novel analogs of A-86929 were synthesized and screened in vitro for structure-functional-selectivity relationships (SFSR) studies to identify chemical motifs responsible for {beta}-arrestin biased activity at both D1 and D2Rs. Most of the A-86929 analogs screened were G protein biased but none of them were exclusively arrestin-biased. Additionally, various catechol aryl fragments were designed and synthesized. Other compounds surveyed included hydroxyl and chloro analogs of dopamine to test for hydrogen bonding and ionic interactions. Some of these small molecular probes displayed weak bias towards the {beta}-arrestin pathway. Continued in-depth SFSR studies informed by structure determination, molecular modeling, and mutagenesis studies will facilitate discovery of potent and efficacious arrestin-biased dopamine receptor ligands. Analyses of symbiotic bacterial communities in the plant pest Bemisia tabaci reveal high prevalence of Candidatus Hemipteriphilus asiaticus on the African continent AO_SCPLOWBSTRACTC_SCPLOWMicrobial symbionts are widespread in insects and some of them have been associated to adaptive changes. Primary symbionts (P-symbionts) have a nutritional role that allows their hosts to feed on unbalanced diets (plant sap, wood, blood). Most of them have undergone genome reduction, but their genomes still retain genes involved in pathways that are necessary to synthesize the nutrients that their hosts need. However, in some P-symbionts, essential pathways are incomplete and secondary symbionts (S-symbionts) are required to complete parts of their degenerated functions. The P-symbiont of the phloem sap-feeder Bemisia tabaci, Candidatus Portiera aleyrodidarium, lacks genes involved in the synthesis of vitamins, cofactors, and also of some essential amino-acids. Seven S-symbionts have been detected in the B. tabaci species complex. Phenotypic and genomic analyses have revealed various effects, from reproductive manipulation to fitness benefits, notably some of them have complementary metabolic capabilities to Candidatus Portiera aleyrodidarium, suggesting that their presence may be obligatory. In order to get the full picture of the symbiotic community of this pest, we investigated, through metabarcoding approaches, the symbiont content of individuals from Burkina Faso, a West African country where B. tabaci induces severe crop damage. While no new putative B. tabaci S-symbiont was identified, Candidatus Hemipteriphilus asiaticus, a symbiont only described in B. tabaci populations from Asia, was detected for the first time on this continent. Phylogenetic analyses however reveal that it is a different strain than the reference found in Asia. Specific diagnostic PCRs showed a high prevalence of these S-symbionts and especially of Candidatus Hemipteriphilus asiaticus in different genetic groups. These results suggest that Candidatus Hemipteriphilus asiaticus may affect the biology of B. tabaci and provide fitness advantage in some B. tabaci populations.