Joint representations of color and form in mouse visual cortex described by random pooling from rods and cones Spatial transitions in color can aid any visual perception task, and its neural representation - the "integration of color and form" - is thought to begin at primary visual cortex (V1). An integration of color and form is untested in mouse V1, yet studies show that the ventral retina provides the necessary substrate from green-sensitive rods and UV-sensitive cones. Here, we used two-photon imaging in V1 to measure spatial frequency (SF) tuning along four axes of rod and cone contrast space, including luminance and color. We first reveal that V1s sensitivity to color is similar to luminance, yet average SF tuning is significantly shifted lowpass for color. Next, guided by linear models, we used SF tuning along all four color axes to estimate the proportion of neurons that fall into classic models of color opponency - "single-", "double-", and "non-opponent". Few neurons (~6%) fit criteria for double-opponency, which are uniquely tuned for chromatic borders. Most of the population can be described as a unimodal distribution ranging from strongly single-opponent to non-opponent. Consistent with recent studies of the rodent and primate retina, our V1 data is well-described by a simple model in which ON and OFF channels to V1 sample the photoreceptor mosaic randomly.
Hippocampal and medial prefrontal cortices encode structural task representations following progressive and interleaved training schedules Memory generalisations may be underpinned by either encoding- or retrieval-based generalisation mechanisms and different training schedules may bias some learners to favour one of these mechanisms over the other. We used a transitive inference task to investigate whether generalisation is influenced by progressive vs randomly interleaved training, and overnight consolidation. On consecutive days, participants learnt pairwise discriminations from two transitive hierarchies before being tested during fMRI. Inference performance was consistently better following progressive training, and for pairs further apart in the transitive hierarchy. BOLD pattern similarity correlated with hierarchical distances in the left hippocampus (HIP) and medial prefrontal cortex (MPFC) following both training schedules. These results are consistent with the use of structural representations that directly encode hierarchical relationships between task features. However, such effects were only observed in the MPFC for recently learnt relationships and were strongest in participants who were least proficient at the inference task. We conclude that humans preferentially employ encoding-based mechanisms to store map-like relational codes that can be used for memory generalisation. These codes are expressed in the HIP and MPFC following both progressive and interleaved training but are not sufficient for accurate inference.
The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. SPIRO - the automated Petri plate imaging platform designed by biologists, for biologists The imaging of plant seedlings, fungal mycelia and bacterial colonies grown on Petri plates is common in phenotyping assays, and is typically done manually despite the procedures being time-consuming and laborious. The main reason for this is the still limited availability of existing automated phenotyping tools and facilities. Additionally, constructing a custom-made automated solution is a daunting task for most research groups specializing in biology.

Here, we describe SPIRO, the Smart Plate Imaging Robot, an automated platform that acquires time-lapse photos of up to four vertically oriented Petri plates in a single experiment. SPIRO was designed for biologists by biologists; thus, its assembly does not require experience in engineering or programming and its operation is sufficiently intuitive to be carried out without training. SPIRO has a small footprint optimal for fitting into standard incubators for plants and microbes, and is equipped with an LED light source for imaging in the dark, thus allowing acquisition of photos under optimal growth conditions. SPIROs web-based user interface allows setting up experiments and downloading data remotely, without interfering with samples growth. SPIROs excellent image quality is optimal for automated image processing, which we demonstrate with two semi-automated assays for analysis of commonly used phenotypic traits: seed germination and root growth.

Moreover, the robot can be easily customized for a specific use, as all information about SPIRO, including the models for 3D-printed structural components, control software, and scripts for image analysis, are released under permissive open-source licenses. 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. Computation of the mitochondrial age distribution along the axon length We describe a compartmental model of mitochondrial transport in axons, which we apply to compute mitochondrial age at different distances from the soma. The model predicts that the average mitochondrial age at the tip of an axon with a length of 1 cm is approximately 22 hours. The mitochondria are youngest closest to the soma and their age scales approximately linearly with distance from the soma. To the best of the authors knowledge, this is the first attempt to predict the spatial distribution of mitochondrial age within an axon. A sensitivity study of the mean age of mitochondria to various model parameters is also presented. Quantitative analysis of gene expression on the inactive human X chromosome The human X chromosome exists in active (Xa) and inactive (Xi) states. To quantify Xi gene expression and explore its implications, we conducted RNA sequencing in two cell types from individuals with one Xa and zero to three Xis. Linear models fit the data well, indicating equal contributions from each Xi and quantifying significant Xi-driven expression changes for 124 of 406 expressed X-chromosome genes. As revealed by allele-specific expression analyses, these Xi-driven changes reflect both i) Xi transcription and ii) Xi modulation - both positive and negative - of Xa transcript levels. Incorporating population and evolutionary measures of constraint on expression levels, we identified 13 X-chromosome genes most likely to drive the phenotypic consequences of natural variation in Xi copy number. These analytics and the underlying datasets offer a foundation for understanding Turner, Klinefelter, and other sex chromosome aneuploidy syndromes, and sex-chromosome effects on health and disease in euploid males and females. Optimization of avian perching manoeuvres Perching at speed is amongst the most challenging flight behaviours that birds perform1,2, and beyond the capability of current autonomous vehicles. Smaller birds may touchdown by hovering3-8, but larger birds typically swoop upward to perch1,2 - presumably because the adverse scaling of their power margin prohibits slow flapping flight9, and because swooping transfers excess kinetic to potential energy1,2,10. Perching is risky in larger birds6,11, demanding precise control of velocity and pose12-15, but it is unknown how they optimize this challenging manoeuvre. More generally, whereas cruising flight behaviours such as migration and commuting are adapted to minimize cost-of-transport or time-of-flight16, the optimization of unsteady flight manoeuvres remains largely unexplored7,17. Here we show that swooping minimizes neither the time nor energy required to perch safely in Harris hawks Parabuteo unicinctus, but instead minimizes the distance flown under hazardous post-stall conditions. By combining motion capture data from 1,563 flights with flight dynamics modelling, we found that the birds choice of where to transition from powered dive to unpowered climb minimizes the distance from the landing perch over which very high lift coefficients are required. Time and energy are therefore invested to maintain the control authority needed to execute a safe landing, rather than being minimized continuously as in technical applications of autonomous perching under nonlinear feedback control13 and deep reinforcement learning18,19. Naive birds learn this behaviour on-the-fly, so our findings suggest an alternative reward function for reinforcement learning of autonomous perching in air vehicles. Induction of exaggerated cytokine production in human peripheral blood mononuclear cells by a recombinant SARS-CoV-2 spike glycoprotein S1 and its inhibition by dexamethasone An understanding of the pathological inflammatory mechanisms involved in SARS-CoV-2 virus infection is necessary in order to discover new molecular pharmacological targets for SARS-CoV-2 cytokine storm. In this study, the effects of a recombinant SARS-CoV-2 spike glycoprotein S1 was investigated in human peripheral blood mononuclear cells (PBMCs). Stimulation of PBMCs with spike glycoprotein S1 (100 ng/mL) resulted in significant elevation in the production of TNF, IL-6, IL-1{beta} and IL-8. However, pre-treatment with dexamethasone (100 nM) caused significant reduction in the release of these cytokines. Further experiments revealed that S1 stimulation of PBMCs increased phosphorylation of NF-{kappa}B p65 and I{kappa}B, and I{kappa}B degradation. DNA binding of NF-{kappa}B p65 was also significantly increased following stimulation with spike glycoprotein S1. Treatment of PBMCs with dexamethasone (100 nM) or BAY11-7082 (1 {micro}M) resulted in inhibition of spike glycoprotein S1-induced NF-{kappa}B activation. Activation of p38 MAPK by S1 was blocked in the presence of dexamethasone and SKF 86002. CRID3, but not dexamethasone pre-treatment produced significant inhibition of S1-induced activation of NLRP3/caspase-1. Further experiments revealed that S1-induced increase in the production of TNF, IL-6, IL-1{beta} and IL-8 was reduced in the presence of BAY11-7082 and SKF 86002, while CRID3 pre-treatment resulted in the reduction of IL-1{beta} production. These results suggest that SARS-CoV-2 spike glycoprotein S1 stimulated PBMCs to release pro-inflammatory cytokines through mechanisms involving activation of NF-{kappa}B, p38 MAPK and NLRP3 inflammasome. It is proposed that the clinical benefits of dexamethasone in COVID-19 is possibly due to its anti-inflammatory activity in reducing SARS-CoV-2 cytokine storm. Parallel Characterization of cis-Regulatory Elements for Multiple Genes Using CRISPRpath Current pooled CRISPR screens for cis-regulatory elements (CREs) can only accommodate one gene based on its expression level. Here, we describe CRISPRpath, a scalable screening strategy for parallelly characterizing CREs of genes linked to the same biological pathway and converging phenotypes. We demonstrate the ability of CRISPRpath for simultaneously identifying functional enhancers of six genes in the 6-thioguanine-induced DNA mismatch repair pathway using both CRISPR interference (CRISPRi) and CRISPR nuclease (CRISPRn) approaches. 60% of the identified enhancers are known promoters with distinct epigenomic features compared to other active promoters, including increased chromatin accessibility and interactivity. Furthermore, by imposing different levels of selection pressure, CRISPRpath can distinguish enhancers exerting strong impact on gene expression from those exerting weak impact. Our results offer a nuanced view of cis-regulation and demonstrate that CRISPRpath can be leveraged for understanding the complex gene regulatory program beyond transcriptional output at scale. A MALDI-TOF assay identifies nilotinib as an inhibitor of inflammation in acute myeloid leukaemia Inflammatory responses are important in cancer, particularly in the context of monocyte-rich aggressive myeloid neoplasm. We developed a label-free cellular phenotypic drug discovery assay to identify anti-inflammatory drugs in human monocytes derived from acute myeloid leukaemia (AML), by tracking several biological features ionizing from only 2,500 cells using matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry. A proof-of-concept screen showed that the BCR-ABL inhibitor nilotinib, but not the structurally similar imatinib, blocks inflammatory responses. In order to identify the cellular (off-)targets of nilotinib, we performed thermal proteome profiling (TPP). Unlike imatinib, nilotinib and other later generation BCR-ABL inhibitors inhibit the p38-MK2/3 signalling axis which suppressed the expression of inflammatory cytokines, cell adhesion and innate immunity markers in activated human monocytes derived from AML. Thus, our study provides a tool for the discovery of new anti-inflammatory drugs, which could contribute to the treatment of inflammation in myeloid neoplasms and other diseases.

Key PointsLabel-free cell-based assay identifies new anti-inflammatory drugs using MALDI-TOF MS. Nilotinib reduces inflammation by inhibition of MAPK14-MK2/3 signalling axis in AML. HYBRID REACTION-DIFFUSION AND CLOCK-AND-WAVEFRONT MODEL FOR THE ARREST OF OSCILLATIONS IN THE SOMITOGENESIS SEGMENTATION CLOCK AO_SCPLOWBSTRACTC_SCPLOWThe clock and wavefront paradigm is arguably the most widely accepted model for explaining the embryonic process of somitogenesis. According to this model, somitogenesis is based upon the interaction between a genetic oscillator, known as segmentation clock, and a differentiation wavefront, which provides the positional information indicating where each pair of somites is formed. Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation. Cotterell et al. called it a progressive oscillatory reaction-diffusion (PORD) model. In the Meinhardt-PORD model, somitogenesis is driven by short-range interactions and the posterior movement of the front is a local, emergent phenomenon, which is not controlled by global positional information. With this model, it is possible to explain some experimental observations that are incompatible with the clock and wavefront model. However the Meinhardt-PORD model has some important disadvantages of its own. Namely, it is quite sensitive to fluctuations and depends on very specific initial conditions (which are not biologically realistic). In this work, we propose an equivalent Meinhardt-PORD model, and then amend it to couple it with a wavefront consisting of a receding morphogen gradient. By doing so, we get a hybrid model between the Meinhardt-PORD and the clock-and-wavefront ones, which overcomes most of the deficiencies of the two originating models.

Somitogenesis, the process by which somites are formed, is an essential developmental stage in many vertebrates. This process occurs with a strikingly regular periodicity, that is preserved among embryos of a single species. The clock and wavefront paradigm is arguably the most widely accepted model for explaining somitogenesis. However, it is incapable of explaining some experimental facts, like the appearance of somites in the absence of an external wavefront (i.e. a receding morphogen gradient). Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation, and called it a progressive oscillatory reaction-diffusion (PORD) model. The Meinhardt-PORD model tackles some of the deficiencies of the clock and wave-front models, but it has some serious issues of its own. In the present work, we introduce an equivalent Meinhardt-PORD model, and then amend it to couple it with a receding morphogen gradient. By doing so, we get a hybrid model that incorporates characteristics of the Meinhardt-PORD and clock-and-wavefront models. We show that this hybrid model under-goes a bifurcation, from a stable to an unstable limit cycle, as the value of the parameter accounting for a background regulatory input (associated to the receding morphogen gradient) decreases. This bifurcation allows the model to explain why somites can form in the absence of an external wavefront, reassesses the role of the receding morphogen gradient as a conductor for somitogenesis, and makes the model behavior robust to random fluctuations, as well as independent from specific initial conditions (the latter, are two of the weak points of the Meinhardt-PORD model). We argue that this findings provide convincing evidence that reaction-diffusion and positional information (receding morphogen gradient) mechanisms could work together in somitogenesis. Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, LK-DFBA, that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this frameworks success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that the optimal constraint approach was generally robust to local perturbations of the system, indicating that just a single wild-type dataset could allow identification of the ideal constraint for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems. Elucidation of SARS-CoV-2 budding mechanisms through molecular dynamics simulations of M and E protein complexes SARS-CoV-2 and other coronaviruses pose major threats to global health, yet computational efforts to understand them have largely overlooked the process of budding, a key part of the coronavirus life cycle. When expressed together, coronavirus M and E proteins are sufficient to facilitate budding into the ER-Golgi intermediate compartment (ERGIC). To help elucidate budding, we ran atomistic molecular dynamics (MD) simulations using the Feig laboratorys refined structural models of the SARS-CoV-2 M protein dimer and E protein pentamer. Our MD simulations consisted of M protein dimers and E protein pentamers in patches of membrane. By examining where these proteins induced membrane curvature in silico, we obtained insights around how the budding process may occur. Multiple M protein dimers acted together to induce global membrane curvature through protein-lipid interactions while E protein pentamers kept the membrane planar. These results could eventually help guide development of antiviral therapeutics which inhibit coronavirus budding. 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. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

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

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

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. Muscle activation strategies of the vastus lateralis according to sex AimDespite men exhibiting greater muscle strength and fatigibility than women, it remains unclear if there are sex-based differences in muscle recruitment strategies e.g. motor unit (MU) recruitment and modulation of firing rate (FR) at normalised forces and during progressive increases in force.

MethodsTwenty-nine healthy male and thirty-one healthy female participants (18-35 years) were studied. Intramuscular electromyography was used to record individual motor unit potentials (MUPs) and near fibre MUPs from the vastus lateralis (VL) during 10% and 25% maximum isometric voluntary contractions (MVC), and spike-triggered averaging was used to obtain motor unit number estimates (MUNE) of the VL. Multilevel mixed-effects linear regression models were used to investigate the effects of sex at each contraction level.

ResultsMen exhibited greater muscle strength (p<0.001) and size (p<0.001) than women, with no difference in force steadiness at 10% or 25% MVC. Women had smaller MUs and higher FR at 10% MVC (both p<0.02), similar to that at 25% MVC in MU size (p=0.062) and FR (p=0.031). However, both sexes showed similar increases in MU size and FR when moving from low-to mid-level contractions. There were no sex differences in any near fibre MUP parameters or in MUNE.

ConclusionIn the vastus lateralis, women produce muscle force via different neuromuscular recruitment strategies to men which is characterised by smaller MUs discharging at higher rates. However, similar strategies are employed to increase force production from low to moderate contractions. These findings of similar proportional increases between sexes support the use of mixed sex cohorts in studies of this nature.

Key pointsO_LIIncreases in muscle force production are mediated by motor unit (MU) recruitment, and MU firing rate (FR).
C_LIO_LIWomen are underrepresented in studies of human neuromuscular research and markedly differ to men in a number of aspects of neuromuscular function, yet little is known of the recruitment strategies of each.
C_LIO_LIHere we demonstrate men and women have similar vastus lateralis MU number estimates, yet women recruit smaller MUs with higher FR than men at normalised contraction levels. However, increases in force are achieved via similar trajectories of MU recruitment and MU FR in men and women.
C_LIO_LIAlthough men and women exhibit divergent neuromuscular recruitment strategies to achieve normalised forces, increases in force are achived similarly and support the inclusion of mixed sex cohorts in studies of this nature.
C_LI Population bottlenecks constrain host microbiome diversity and genetic variation impeding fitness It is becoming increasingly clear that microbial symbionts influence key aspects of their hosts fitness, and vice versa. This may fundamentally change our thinking about how microbes and hosts interact in influencing fitness and adaptation to changing environments. Here we explore how reductions in population size commonly experienced by threatened species influence microbiome diversity. Consequences of such reductions are normally interpreted in terms of a loss of genetic variation, increased inbreeding and associated inbreeding depression. However, fitness effects of population bottlenecks might also be mediated through microbiome diversity, such as through loss of functionally important microbes. Here we utilise 50 Drosophila melanogaster lines with different histories of population bottlenecks to explore these questions. The lines were phenotyped for egg-to-adult viability and their genomes sequenced to estimate genetic variation. The bacterial 16S rRNA gene was amplified in these lines to investigate microbial diversity. We found that 1) host population bottlenecks constrained microbiome richness and diversity, 2) core microbiomes of hosts with low genetic variation were constituted from subsets of microbiomes found in flies with higher genetic variation, 3) both microbiome diversity and host genetic variation contributed to host population fitness, 4) connectivity and robustness of bacterial networks was low in the inbred lines regardless of host genetic variation, 5) reduced microbial diversity was associated with weaker evolutionary responses of hosts in stressful environments, and 6) these effects were unrelated to Wolbachia density. These findings suggest that population bottlenecks reduce hologenomic variation (combined host and microbial genetic variation). Thus, while the current biodiversity crisis focuses on population sizes and genetic variation of eukaryotes, an additional focal point should be the microbial diversity carried by the eukaryotes, which in turn may influence host fitness and adaptability with consequences for the persistence of populations.

Author summaryIt is becoming increasingly clear that organisms and the microbes that live on or in them - their microbiome - affect each other in profound ways that we are just beginning to understand. For instance, a diverse microbiome can help maintain metabolic functions or fight pathogens causing diseases. A disrupted microbiome may be especially critical for animals and plants that occur in low numbers because of threats from e.g. human exploitation or climate change, as they may already suffer from genetic challenges such as inbreeding and reduced evolutionary potential. The importance of such a reduction in population size, called a bottleneck, on the microbial diversity and the potential interactive effects on host health remains unexplored. Here we experimentally test these associations by investigating the microbiomes of 50 inbred or non-inbred populations of vinegar flies. We found that restricting the population size constrain the hosts genetic variation and simultaneously decreases the diversity of the microbiome that they harbor, and that both effects were detrimental to host fitness. The microbial communities in inbred host populations were less robust than in their non-inbred counterparts, suggesting that we should increasingly consider the microbiome diversity, which may ultimately influence the health and persistence of threatened species. Is Vaccination a Viable Method to Control Johne's Disease Caused by Mycobacterium avium subsp. paratuberculosis? Data from 12 Million Ovine Vaccinations and 7.6 Million Carcass Examinations in New South Wales, Australia from 1999-2009 BackgroundMycobacterium avium subsp. paratuberculosis (MAP) causes Johnes disease (or paratuberculosis), a chronic wasting disease of ruminants and other animals resulting from granulomatous enteritis. There are increasing concerns that MAP is zoonotic. The prevalence of Johnes disease is increasing worldwide. In an attempt to control an epidemic of ovine Johnes disease (OJD) in New South Wales (NSW), a government/industry sponsored voluntary vaccination/on-farm management program commenced in 2000. We report herein an observational study of changes in disease prevalence as vaccination progressed, based on abattoir surveillance data for OJD from 1999 to 2009. We also discuss the epidemiological, policy, regulatory, research, economic and sociological elements that contributed to the development of a mature control program, whose aim was to halt the epidemic spread of OJD in a naive sheep population.

MethodsNSW was divided into areas of "High" (HPA), "Medium" (MPA) and "Low" (LPA) OJD prevalence. A killed whole cell vaccine (Gudair(R)) was administered to sheep from 2000 to 2009. Trained examiners evaluated the viscera of adult sheep carcasses at slaughter for gross evidence of OJD. MAP infection was confirmed by histopathology.

Principal findingsFrom 2000-2009, 12 million vaccine doses were administered in NSW (91%; 10.9 million in the HPA). Many of the vaccinated flocks were suffering > 5% annual mortality in adult sheep, with some individual flocks with 10-15% losses attributable to OJD. A total of 7.6 million carcasses were examined (38%; 2.9 million from the HPA). Overall, 16% of slaughter consignments (sheep consigned to the abattoir from a single vendor) were positive for OJD, of which 94% were from the HPA. In the HPA, the percentage of animals with lesions attributable to OJD at slaughter fell progressively from 2.4% (10,406/432,860) at commencement of vaccination in 2000 to 0.8% (1,573/189,564) by 2009. Herd immunity from vaccination in the HPA was estimated at 70% by 2009. This coincided with a progressive decrease in reports of clinical disease and mortalities in vaccinated flocks.

SignificanceWe show a decrease in the prevalence of lesions attributable to OJD in NSW concomitant with initiation of voluntary vaccination, on-farm management plans, abattoir monitoring and feedback of animal prevalence data to sheep producers. We conclude that a target of [&le;] 1% regional prevalence of OJD affected sheep at slaughter is achievable using these interventions. Defining An Expanded RAS Conformational Landscape Based on Over 700 Experimentally Determined Structures of KRAS, NRAS, and HRAS For many human cancers and tumor-associated diseases, mutations in the RAS isoforms (KRAS, NRAS, and HRAS) are the most common oncogenic alterations, making these proteins high-priority therapeutic targets. Effectively targeting the RAS isoforms requires an exact understanding of their active, inactive, and druggable conformations. However, there is no structure-guided catalogue of RAS conformations to guide therapeutic targeting or examining the structural impact of RAS mutations. We present an expanded classification of RAS conformations based on analyzing their catalytic switch 1 (SW1) and switch 2 (SW2) loops. From all 721 available human KRAS, NRAS, and HRAS structures in the Protein Data Bank (PDB) (206 RAS-protein complexes, 190 inhibitor-bound, and 325 unbound, including 204 WT and 517 mutated structures), we created a broad conformational classification based on the spatial positions of residue Y32 in SW1 and residue Y71 in SW2. Subsequently, we defined additional conformational subsets (some previously undescribed) by clustering all well modeled SW1 and SW2 loops using a density-based machine learning algorithm with a backbone dihedral-based distance metric. In all, we identified three SW1 conformations and nine SW2 conformations, each which are associated with different nucleotide states (GTP-bound, nucleotide-free, and GDP-bound) and specific bound proteins or inhibitor sites. The GTP-bound SW1 conformation can be further subdivided based on the hydrogen (H)-bond type made between residue Y32 and the GTP {gamma}-phosphate: water-mediated, direct, or no H-bond. Further analyzing these structures clarified the catalytic impact of the G12D and G12V RAS mutations, and the inhibitor chemistries that bind to each druggable RAS conformation. To facilitate future RAS structural analyses, we have created a web database, called Rascore, presenting an updated and searchable dataset of human KRAS, NRAS, and HRAS structures in the PDB, and which includes a page for analyzing user uploaded RAS structures by our algorithm (http://dunbrack.fccc.edu/rascore/).

SignificanceAnalyzing >700 experimentally determined RAS structures helped define an expanded landscape of active, inactive and druggable RAS conformations, the structural impact of common RAS mutations, and previously uncharacterized RAS-inhibitor binding modes. SkewDB: A comprehensive database of GC and 10 other skews for over 28,000 chromosomes and plasmids GC skew denotes the relative excess of G nucleotides over C nucleotides on the leading versus the lagging replication strand of eubacteria. While the effect is small, typically around 2.5%, it is robust and pervasive. GC skew and the analogous TA skew are a localized deviation from Chargaffs second parity rule, which states that G and C, and T and A occur with (mostly) equal frequency even within a strand.

Most bacteria also show the analogous TA skew. Different phyla show different kinds of skew and differing relations between TA and GC skew.

This article introduces an open access database (https://skewdb.org) of GC and 10 other skews for over 28,000 chromosomes and plasmids. Further details like codon bias, strand bias, strand lengths and taxonomic data are also included. The SkewDB database can be used to generate or verify hypotheses. Since the origins of both the second parity rule, as well as GC skew itself, are not yet satisfactorily explained, such a database may enhance our understanding of microbial DNA. A nickase Cas9 gene-drive system promotes super-Mendelian inheritance in Drosophila CRISPR-based gene-drive systems have been proposed for managing insect populations, including disease-transmitting mosquitoes, due to their ability to bias their inheritance towards super-Mendelian rates (>50%). Current technologies employ a Cas9 that introduces DNA double-strand breaks into the opposing wildtype allele to replace it with a copy of the gene drive allele via DNA homology-directed repair. Yet, the use of different Cas9s versions is unexplored, and alternative approaches could increase the available toolkit for gene-drive designs. Here, we report a novel gene-drive approach that relies on Cas9 nickases that generate staggered paired nicks in DNA to propagate the engineered gene-drive cassette. We show that generating 5 overhangs in the system yields efficient allelic conversion. The nickase gene-drive arrangement produces large, stereotyped deletions that are advantageous for targeting essential genes. Indeed, our nickase approach should expand the repertoire for gene-drive designs aimed at applications in mosquitoes and beyond. Relationships between aboveground plant traits and carbon cycling in tundra plant communities O_LIThe functional composition and diversity of plant communities are globally applicable predictors of ecosystem functioning. Yet, it is unclear how traits influence carbon cycling. This is an important question in the tundra where vegetation shifts are occurring across the entire biome, and where soil organic carbon stocks are large and vulnerable to environmental change.
C_LIO_LITo study how traits affect carbon cycling in the tundra, we built a model that explained carbon cycling (above-ground and soil organic carbon stocks, and photosynthetic and respiratory fluxes) with abiotic conditions (air temperature and soil moisture), plant community functional composition (average plant height, leaf dry matter content (LDMC) and specific leaf area (SLA)), and functional diversity (weighted standard deviations of the traits). Data was collected from an observational study setting from northern Finland.
C_LIO_LIThe explanatory power of the models was relatively high, but a large part of variation in soil organic carbon stocks remained unexplained. Plant height was the strongest predictor of all carbon cycling variables except soil carbon stocks. Communities of larger plants were associated with larger CO2 fluxes and above-ground carbon stocks. Communities with fast leaf economics (i.e. high SLA and low LDMC) had higher photosynthesis, ecosystem respiration, and soil organic carbon stocks.
C_LIO_LIWithin-community variability in plant height, SLA, and LDMC affected ecosystem functions differently. SLA and LDMC diversity increased CO2 fluxes and soil organic carbon stocks, while height diversity increased the above-ground carbon stock. The contributions of functional diversity metrics to ecosystem functioning were about as important as those of average SLA and LDMC traits.
C_LIO_LISynthesis: Plant height, SLA, and LDMC have clear effects on tundra carbon cycling. The importance of functional diversity highlights a potentially important mechanism controlling the vast tundra carbon pools that should be better recognized. More research on root traits and decomposer communities is needed to understand the below-ground mechanisms regulating carbon cycling in the tundra.
C_LI 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 Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. CTP and parS coordinate ParB partition complex dynamics and ParA-ATPase activation for ParABS-mediated DNA partitioning ParABS partition systems, comprising the centromere-like DNA sequence parS, the parS-binding ParB-CTPase and the nucleoid-binding ParA-ATPase, ensure faithful segregation of bacterial chromosomes and low-copy-number plasmids. F-plasmid partition complexes containing ParBF and parSF move by generating and following a local concentration gradient of nucleoid-bound ParAF. However, the process through which ParBF activates ParAF-ATPase has not been defined. We studied CTP- and parSF-modulated ParAF--ParBF complex assembly, in which DNA-bound ParAF-ATP dimers are activated for ATP hydrolysis by interacting with two ParBF N-terminal domains. CTP or parSF enhances the ATPase rate without significantly accelerating ParAF--ParBF complex assembly. Together, parSF and CTP accelerate ParAF--ParBF assembly without further significant increase in ATPase rate. Magnetic-tweezers experiments showed that CTP promotes multiple ParBF loading onto parSF-containing DNA, generating condensed partition complex-like assemblies. We propose that ParBF in the partition complex adopts a conformation that enhances ParBF--ParBF and ParAF--ParBF interactions promoting efficient partitioning.