Admixture of evolutionary rates across a hybrid zone Hybridization is a major evolutionary force that can erode genetic differentiation between species, whereas reproductive isolation maintains such differentiation. In studying a hybrid zone between the swallowtail butterflies Papilio syfanius and Papilio maackii, we made the unexpected discovery that genomic substitution rates are unequal between the parental species. This phenomenon creates a novel process in hybridization, where genomic regions most affected by gene flow evolve at similar rates, while genomic regions with greater reproductive isolation evolve at divergent rates. Thus, hybridization mixes evolutionary rates in a way similar to its effect on ancestry. Using coalescent theory, we show that the rate-mixing process provides distinct information about levels of gene flow across different parts of genomes, and that maintenance of divergent substitution rates can be predicted quantitatively from relative sequence divergence (FST) between the hybridizing species at equilibrium. A corollary is that divergent rates will be maintained in regions linked to barrier loci. Overall, we demonstrate that reproductive isolation maintains not only the final outcome of genomic differentiation, but also the rate at which differentiation accumulates. This new information also suggests that the separation of evolutionary rates co-localizes with the separation of gene pools between genomes of incipient species.
Multiple types of genomic variation contribute to adaptive traits in the mustelid subfamily Guloninae Species of the mustelid subfamily Guloninae inhabit diverse habitats on multiple continents, and occupy a variety of ecological niches. They differ in feeding ecologies, reproductive strategies and morphological adaptations. To identify candidate loci associated with adaptations to their respective environments, we generated a de novo assembly of the tayra (Eira barbara), the earliest diverging species in the subfamily, and compared this with the genomes available for the wolverine (Gulo gulo) and the sable (Martes zibellina). Our comparative genomic analyses included searching for signs of positive selection, examining changes in gene family sizes, as well as searching for species-specific structural variants (SVs). Among candidate loci associated with phenotypic traits, we observed many related to diet, body condition and reproduction. For example, for the tayra, which has an atypical gulonine reproductive strategy of aseasonal breeding, we observe species-specific changes in many pregnancy-related genes. For the wolverine, a circumpolar hypercarnivore that must cope with seasonal food scarcity, we observed many changes in genes associated with diet and body condition. All types of genomic variation examined contributed substantially to the identification of candidate loci. This strongly argues for consideration of variation other than single nucleotide polymorphisms in comparative genomics studies aiming to identify loci of adaptive significance.
Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector. Insertions in SARS-CoV-2 genome caused by template switch and duplications give rise to new variants that merit monitoring The appearance of multiple new SARS-CoV-2 variants during the winter of 2020-2021 is a matter of grave concern. Some of these new variants, such as B.1.617.2, B.1.1.7, and B.1.351, manifest higher infectivity and virulence than the earlier SARS-CoV-2 variants, with potential dramatic effects on the course of the COVID-19 pandemic. So far, analysis of new SARS-CoV-2 variants focused primarily on point nucleotide substitutions and short deletions that are readily identifiable by comparison to consensus genome sequences. In contrast, insertions have largely escaped the attention of researchers although the furin site insert in the spike protein is thought to be a determinant of SARS-CoV-2 virulence and other inserts might have contributed to coronavirus pathogenicity as well. Here, we investigate insertions in SARS-CoV-2 genomes and identify 347 unique inserts of different lengths. We present evidence that these inserts reflect actual virus variance rather than sequencing errors. Two principal mechanisms appear to account for the inserts in the SARS-CoV-2 genomes, polymerase slippage and template switch that might be associated with the synthesis of subgenomic RNAs. We show that inserts in the Spike glycoprotein can affect its antigenic properties and thus merit monitoring. At least, three inserts in the N-terminal domain of the Spike (ins245IME, ins246DSWG, and ins248SSLT) that were first detected in 2021 are predicted to lead to escape from neutralizing antibodies, whereas other inserts might result in escape from T-cell immunity. Plant-environment microscopy tracks interactions of Bacillus subtilis with plant roots across the entire rhizosphere Our understanding of plant-microbe interactions in soil is limited by the difficulty of observing processes at the microscopic scale throughout plants large volume of influence. Here, we present the development of 3D live microscopy for resolving plant-microbe interactions across the environment of an entire seedling growing in a transparent soil in tailor-made mesocosms, maintaining physical conditions for the culture of both plants and microorganisms. A tailor made dual-illumination light-sheet system acquired scattering signals from the plant whilst fluorescence signals were captured from transparent soil particles and labelled microorganisms, allowing the generation of quantitative data on samples approximately 3600 mm3 in size with as good as 5 m resolution at a rate of up to one scan every 30 minutes. The system tracked the movement of Bacillus subtilis populations in the rhizosphere of lettuce plants in real time, revealing previously unseen patterns of activity. Motile bacteria favoured small pore spaces over the surface of soil particles, colonising the root in a pulsatile manner. Migrations appeared to be directed towards the root cap, the point "first contact", before subsequent colonisation of mature epidermis cells. Our findings show that microscopes dedicated to live environmental studies present an invaluable tool to understand plant-microbe interactions. 3CAC: improving the classification of phages and plasmids in metagenomic assemblies using assembly graphs MotivationBacteriophages and plasmids usually coexist with their host bacteria in microbial communities and play important roles in microbial evolution. Accurately identifying sequence contigs as phages, plasmids, and bacterial chromosomes in mixed metagenomic assemblies is critical for further unravelling their functions. Many classification tools have been developed for identifying either phages or plasmids in metagenomic assemblies. However, only two classifiers, PPR-Meta and viralVerify, were proposed to simultaneously identify phages and plasmids in mixed metagenomic assemblies. Due to the very high fraction of chromosome contigs in the assemblies, both tools achieve high precision in the classification of chromosomes but perform poorly in classifying phages and plasmids. Short contigs in these assemblies are often wrongly classified or classified as uncertain.

ResultsHere we present 3CAC, a new three-class classifier that improves the precision of phage and plasmid classification. 3CAC starts with an initial three-class classification generated by existing classifiers and improves the classification of short contigs and contigs with low confidence classification by using proximity in the assembly graph. Evaluation on simulated metagenomes and on real human gut microbiome samples showed that 3CAC outperformed PPR-Meta and viralVerify in both precision and recall, and increased F1-score by 10-60 percentage points.

AvailabilityThe 3CAC software is available on https://github.com/Shamir-Lab/3CAC.

Contactrshamir@tau.ac.il

Supplementary informationSupplementary data are available at Bioinformatics online. Conformity leading self-organised rescue from a suboptimal decision-making bias Given the ubiquity of potentially adverse behavioural bias owing to myopic trial-and-error learning, it seems paradoxical that improvements in decision-making performance through conformist social learning, a process widely considered to be bias amplification, still prevail in animal collective behaviour. Here we show, through model analyses and large-scale interactive behavioural experiments with 585 human subjects, that conformist influence can indeed promote favourable risk taking in repeated experience-based decision making, even though many individuals are systematically biased towards adverse risk aversion. Although strong positive feedback conferred by copying the majoritys behaviour could result in unfavourable informational cascades, our differential equation model of collective behavioural dynamics identified a key role for increasing exploration by negative feedback arising when a weak minority influence undermines the inherent behavioural bias. This  collective behavioural rescue, emerging through coordination of positive and negative feedback, highlights a benefit of collective learning in a broader range of environmental conditions than previously assumed and resolves the ostensible paradox of adaptive collective behavioural flexibility under conformist influences. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. Intron-mediated induction of phenotypic heterogeneity Introns are universally present in the nuclear genomes of eukaryotes1. The budding yeast, an otherwise intron-poor species, preserves two sets of ribosomal protein (RP) genes differing primarily in their introns2-4. Despite recent findings on the role of RP introns under stress and starvation5-7, understanding the contribution of introns to ribosome regulation remains challenging. Here, combining isogrowth profiling8 with single-cell protein measurements9, we found that introns can mediate inducible phenotypic heterogeneity conferring a clear fitness advantage. Osmotic stress leads to bimodal expression of the small ribosomal subunit protein Rps22B, mediated by an intron in the 5 untranslated region of its transcript. The two resulting yeast subpopulations differ in their ability to cope with starvation. Low Rps22B protein levels resulted in prolonged survival under sustained starvation, while high Rps22B levels enabled cells to grow faster after transient starvation. Further, yeast growing at high sugar concentrations - similar to those in ripe grapes - exhibit bimodal Rps22B expression when approaching stationary phase. Differential intron-mediated regulation of RP genes thus provides a way to diversify the population when starvation looms in natural environments. Our findings reveal a new role for introns in inducing phenotypic heterogeneity in changing environments and suggest that duplicated RP genes in yeast contribute to resolving the evolutionary conflict between precise expression control and environmental responsiveness10. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Cargo-mediated mechanisms reduce inter-motor mechanical interference, promote load-sharing and enhance processivity in teams of molecular motors In cells, multiple molecular motors work together as teams to carry cargoes such as vesicles and organelles over long distances to their destinations by stepping along a network of cytoskeletal filaments. How motors that typically mechanically interfere with each other, work as teams is unclear. Here we explored the possibility that purely physical mechanisms may potentially enhance teamwork, both at the single motor and cargo level. To explore these mechanisms, we developed a 3D dynamical simulation of cargo transport along microtubules by teams of canonically non-cooperative kinesin-1 motors. We accounted for cargo membrane fluidity by explicitly simulating the Brownian dynamics of motors on the cargo surface and considered both the load and ATP dependence of single motor functioning. We showed explicitly, for the first time, that surface fluidity leads to the reduction of negative mechanical interference between kinesins, enhancing load sharing thereby decreasing single motor off-rates and increasing processivity. Remarkably, we also showed that, independent of fluidity, increasing numbers of bound motors pull the cargo closer to the microtubule, increasing the on-rates of individual unbound motors, resulting in a cooperative increase in bound motor numbers that depends on 3D cargo geometry. At the cargo level, surface fluidity makes more motors available for binding, though this effect is significant only at low ATP or high motor density. Interestingly, we find that the fluidity induced reduction in mechanical interference dominates over the increased availability of motors for typical physiological and in vitro conditions, this allowing us to reconcile different experimental results in different regimes. Finally, we show that these effects can altogether result in enhanced mechanical efficiency and tunable cargo run-lengths for teams of molecular motors over physiological ranges of fluidity with implications for new experimental validation efforts, transport in vivo, artificial cargo design and multi-motor-driven mechanical processes in general.

Author summaryIn cells, multiple molecular motors work together as teams to carry cargoes such as vesicles and organelles over long distances to their destinations by stepping along a network of protein filaments. How do intrinsically non-cooperative motors function as teams within cells? In this paper, we show, using computer simulations, that the fluid surfaces of cellular cargo reduce the mechanical interference between motors allowing better load sharing and decreasing their detachment rates, thereby increasing the distance over which they can carry cargo. We additionally found that motors pull cargo closer to the filament increasing the attachment rates of other unbound motors. These effects act synergistically with an increased availability of motors due to fluidity to further increase travel distances. Our simulation model accounted for the dynamics of the motors on the cargo surface and their load and ATP dependent kinetics, allowing us to connect single motor properties to overall transport. Our work on understanding how teamwork arises in mechanically coupled motors sheds new light on cellular processes, reconciles existing observations, encourages new experimental validation efforts and can also suggest new ways of improving mechanical efficiency and transport of artificial cargo powered by motor teams. Nutrient metabolism regulates insulin granule formation in the pancreatic islet β-cell via ER redox homeostasis Defects in the pancreatic {beta}-cells secretion system are well-described in Type 2 diabetes (T2D) and include impaired proinsulin processing and a deficit in mature insulin-containing secretory granules; however, the cellular mechanisms underlying these defects and the contribution of hyperglycemia to this process remain poorly understood. Here, we used an in situ fluorescent pulse-chase strategy and proximity labeling-based quantitative proteomics analysis to study proinsulin trafficking and demonstrate a direct link to glucose metabolism via the production of redox intermediates that facilitate proinsulin export from the ER. We show that ER export of proinsulin is delayed in T2D models resulting in decreased insulin granule formation and further demonstrate this process can be regulated by NADPH and reducing equivalent availability. Together, these data highlight a critical role for nutrient metabolism and mitochondrial dysfunction in the maladaptive remodeling of the {beta}-cells secretory pathway in the decline of {beta}-cell function in T2D. Modulation of RNA polymerase processivity affects double-strand break repair in the presence of a DNA end-binding protein Homologous recombination is the predominant double-strand DNA break repair pathway in Escherichia coli. A recent report of non-homologous end joining involving the Ku-like Gam protein of bacteriophage Mu (MuGam) has sparked interest in understanding the functions of DNA end-binding proteins in bacteria. MuGam binds to DNA ends, but how it interferes with DNA repair or transcription in live bacteria remains elusive. In E. coli, RNA polymerase secondary channel interactors, such as the Gre factors, play a role in coordinating transcription with DNA replication and break repair. Here we show that MuGam inhibits break repair by slowing down DNA resection and impeding recombination in living cells. Loss of GreA restores break repair in the presence of MuGam by allowing for increased DNA resection due to a potential release of MuGam from the DNA. Using MuGam as a DNA break sensor, we found that breaks are generated when translation is inhibited, more so in the presence of GreA, supporting the model where transcription-translation uncoupling increases transcription/replication collisions. Significantly, this work reveals that modulation of RNA polymerase can impact DNA break repair in presence of a Ku-like protein. Loss of Nuclear DNA Ligase III Reverts PARP Inhibitor Resistance in BRCA1/53BP1 Double-deficient Cells by Exposing ssDNA Gaps Inhibitors of poly(ADP-ribose) (PAR) polymerase (PARPi) have entered the clinic for the treatment of homologous recombination (HR)-deficient cancers. Despite the success of this approach, preclinical and clinical research with PARPi has revealed multiple resistance mechanisms, highlighting the need for identification of novel functional biomarkers and combination treatment strategies. Functional genetic screens performed in cells and organoids that acquired resistance to PARPi by loss of 53BP1 identified loss of LIG3 as an enhancer of PARPi toxicity in BRCA1-deficient cells. Enhancement of PARPi toxicity by LIG3 depletion is dependent on BRCA1 deficiency but independent of the loss of 53BP1 pathway. Mechanistically, we show that LIG3 loss promotes formation of MRE11-mediated post-replicative ssDNA gaps in BRCA1-deficient and BRCA1/53BP1 double-deficient cells exposed to PARPi, leading to an accumulation of chromosomal abnormalities. LIG3 depletion also enhances efficacy of PARPi against BRCA1-deficient mammary tumors in mice, suggesting LIG3 as a potential therapeutic target. A chemical screen based on an interruption of zebrafish gastrulation identifies the HTR2C inhibitor Pizotifen as a suppressor of EMT-mediated metastasis Metastasis is responsible for approximately 90% of cancer-associated mortality but few models exist that allow for rapid and effective screening of anti-metastasis drugs. Current mouse models of metastasis are too expensive and time consuming to use for rapid and high-throughput screening. Therefore, we created a unique screening concept utilizing conserved mechanisms between zebrafish gastrulation and cancer metastasis for identification of potential anti-metastatic drugs. We hypothesized that small chemicals that interrupt zebrafish gastrulation might also suppress metastatic progression of cancer cells and developed a phenotype-based chemical screen to test the hypothesis. The screen used epiboly, the first morphogenetic movement in gastrulation, as a marker and enabled 100 chemicals to be tested in five hours. The screen tested 1280 FDA-approved drugs and identified Pizotifen, an antagonist for serotonin receptor 2C (HTR2C) as an epiboly-interrupting drug. Pharmacologic and genetic inhibition of HTR2C suppressed metastatic progression in a mouse model. Blocking HTR2C with Pizotifen restored epithelial properties to metastatic cells through inhibition of Wnt-signaling. In contrast, HTR2C induced epithelial to mesenchymal transition (EMT) through activation of Wnt-signaling and promoted metastatic dissemination of human cancer cells in a zebrafish xenotransplantation model. Taken together, our concept offers a novel platform for discovery of anti-metastasis drugs. Honeycomb crystallography: comb formation under geometric frustrations 1As honeybees build their nests in pre-existing tree cavities, they must deal with dealing with the presence of geometric constraints, resulting in non-regular hexagons and topological defects in the comb. In this work, we study how bees adapt to their environment in order to regulate the comb structure. Specifically, we identify the irregularities in honeycomb structure in the presence of various geometric frustrations. We 3D-print experimental frames with a variety of constraints imposed on the imprinted foundations. The combs constructed by the bees show clear evidence of reoccurring patterns built by bees in response to specific geometric frustrations on these starter frames. Furthermore, using an experimental-modeling framework, we demonstrate that these patterns can be successfully modeled and replicated through a simulated annealing process, in which the minimized potential is a variation of the Lennard-Jones potential that only considers first-neighbor interactions according to a Delaunay triangulation. Our simulation results not only confirm the connection between honeycomb structures and other crystal systems such as graphene, but also show that irregularities in the honeycomb structure can be explained as the result of local interactions between honeybees and their immediate surroundings, leading to emergent global order. Additionally, our computational model can be used to describe specific strategies that bees use to effectively solve each geometric mismatch problem while minimizing cost of comb building. Embryonic exposure to cannabidiol disrupts active neural circuits, an effect increased by Δ-9-tetrahydrocannabinol and involving CB1R and CB2R in zebrafish In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. Small proline-rich proteins (SPRRs) are epidermally produced antimicrobial proteins that defend the cutaneous barrier by direct bacterial membrane disruption Human skin functions as a physical barrier, preventing the entry of foreign pathogens while also accommodating a myriad of commensal microorganisms. A key contributor to the skin landscape is the sebaceous gland. Mice devoid of sebocytes are prone to skin infection, yet our understanding of how sebocytes function in host defense is incomplete. Here we show that the small proline-rich proteins, SPRR1 and SPRR2 are bactericidal in skin. SPRR1B and SPPR2A were induced in human sebocytes by exposure to the bacterial cell wall component lipopolysaccharide (LPS). Further, LPS injected into mouse skin triggered the expression of the mouse SPRR orthologous genes, Sprr1a and Sprr2a, through stimulation of MYD88. Both mouse and human SPRR proteins displayed potent bactericidal activity against MRSA (methicillin-resistant Staphylococcus aureus), Pseudomonas aeruginosa and skin commensals. Thus, Sprr1a-/-;Sprr2a-/- mice are more susceptible to MRSA and Pseudomonas aeruginosa skin infection. Lastly, mechanistic studies demonstrate that SPRR proteins exert their bactericidal activity through binding and disruption of the bacterial membrane. Taken together, these findings provide insight into the regulation and antimicrobial function of SPRR proteins in skin and how the skin defends the host against systemic infection. Numerical Simulation of Concussive-generated Cortical Spreading Depolarization to Optimize DC-EEG Electrode Spacing for Non-invasive Visual Detection BackgroundCortical Spreading Depolarization (SD) is a propagating depolarization wave of neurons and glial cells in the cerebral gray matter. SD occurs in all forms of severe acute brain injury as documented using invasive detection methods. Based on many experimental studies of mechanical brain deformation and concussion, the occurrence of SDs in human concussion has often been hypothesized. However, this hypothesis cannot be confirmed in humans as SDs can only be detected with invasive detection methods that would require either a craniotomy or a burr hole to be performed on athletes. Typical electroencephalography (EEG) electrodes, placed on the scalp, can detect the possible presence of SD but have not been able to accurately and reliably identify SDs.

MethodsTo explore the possibility of a non-invasive method to resolve this hurdle, we developed a finite element numerical model that simulates scalp voltage changes that are induced by a brain-surface SD. We then compared our simulation results with retrospectively evaluated data in aneurysmal subarachnoid hemorrhage (aSAH) patients from Drenckhahn et al. (Brain 135:853, 2012).

ResultsThe ratio of peak scalp to simulated peak cortical voltage, Vscalp/Vcortex, was 0.0735, whereas the ratio from the retrospectively evaluated data was 0.0316 (0.0221, 0.0527) [median (1st quartile, 3rd quartile), n = 161, p < 0.001, one sample Wilcoxon signed rank test]. These differing values provide validation because their differences can be attributed to differences in shape between concussive- and aSAH-SDs, as well as the inherent limitations in human study voltage measurements. This simulated scalp surface potential was used to design a virtual scalp detection array. Error analysis and visual reconstruction showed that 1 cm is the optimal electrode spacing to visually identify the propagating scalp voltage from a cortical SD. Electrode spacings of 2 cm and above produce distorted images and high errors in the reconstructed image.

ConclusionOur analysis suggests that concussive (and other) SDs can be detected from the scalp, which could confirm SD occurrence in human concussion, provide concussion diagnosis based on an underlying physiological mechanism, and lead to non-invasive SD detection in the setting of severe acute brain injury. Regulation of the PKD2 channel function by TACAN Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in membrane receptor PKD1 or cation channel PKD2. TACAN (also named TMEM120A), recently reported as an ion channel in neuron cells for mechano and pain sensing, is also distributed in diverse non-neuronal tissues such as kidney, heart and intestine, suggesting its involvement in other functions. In this study, we found that TACAN is in complex with PKD2 in native renal cell lines. Using the two-electrode voltage clamp in Xenopus oocytes we found that TACAN inhibited the channel activity of PKD2 gain-of-function mutant F604P. The first and last transmembrane domains of TACAN were found to interact with the PKD2 C-and N-terminal portions, respectively. We showed that the TACAN N-terminus acted as a blocking peptide and that TACAN inhibits the PKD2 function through the PKD2/TACAN binding. By patch clamping in mammalian cells, we found that TACAN inhibits both the single channel conductance and open probability of PKD2 and mutant F604P. PKD2 co-expressed with TACAN, but not PKD2 alone, exhibited pressure sensitivity. Furthermore, we also found that TACAN aggravates PKD2-dependent tail curvature and pronephric cysts in larval zebrafish, in support of the in vitro inhibitory effects of TACAN. In summary, this study revealed that TACAN acts as a PKD2 inhibitor and mediates mechano sensitivity of the PKD2/TACAN channel complex. Microorganisms Characterization in Anaerobic Solventogenic Processing of Volatile Fatty Acids 1This work evaluates four different physicochemical pretreatments (acidic, thermal, acidic-thermal and thermal-acidic) on an anaerobic inoculum used for alcohol production from acetate and butyrate. All experiments were conducted in single batches using acetate and butyrate as substrates at 30 {degrees}C and with a pressurized headspace of pure H2 at 2.15 atm (218.2 MPa). Thermal and acidic-thermal pretreatments lead to higher production of both ethanol and butanol. Mathematical modelling shows that the highest attainable concentrations of ethanol and butanol produced were 122 mg L-1 and 97 mg L-1 for the thermal pretreatment (after 17.5 days) and 87 mg L-1 and 143 mg L-1 for the acidic-thermal pretreatment (after 18.9 days). Acetate was produced in all assays. Thermodynamic data indicated that a high H2 partial pressure favoured solventogenic metabolic pathways. Finally, sequencing data showed that both thermal and acidic-thermal pretreatments selected mainly the bacterial genera Pseudomonas, Brevundimonas and Clostridium. The acidic-thermal pretreatment selected a bacterial community more adapted to the conversion of acetate and butyrate into ethanol and butanol, respectively. Thermal-acidic pretreatment was unstable, showing significant variability between replicates. Acidic pretreatment showed the lowest alcohol production. 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. Genetic Loci and Metabolic States Associated With Murine Epigenetic Aging Changes in DNA methylation (DNAm) are linked to aging. Here, we profile highly conserved CpGs in 339 predominantly female mice belonging to the BXD family for which we have deep longevity and genomic data. We use a  pan-mammalian microarray that provides a common platform for assaying the methylome across mammalian clades. We computed epigenetic clocks and tested associations with DNAm entropy, diet, weight, metabolic traits, and genetic variation. We describe the multifactorial variance of methylation at these CpGs, and show that high fat diet augments the age-associated changes. Entropy increases with age. The progression to disorder, particularly at CpGs that gain methylation over time, was predictive of genotype-dependent life expectancy. The longer-lived BXD strains had comparatively lower entropy at a given age. We identified two genetic loci that modulate rates of epigenetic age acceleration (EAA): one on chromosome (Chr) 11 that encompasses the Erbb2/Her2 oncogenic region, and a second on Chr19 that contains a cytochrome P450 cluster. Both loci harbor genes associated with EAA in humans including STXBP4, NKX2-3, and CUTC. Transcriptome and proteome analyses revealed associations with oxidation-reduction, metabolic, and immune response pathways. Our results highlight concordant loci for EAA in humans and mice, and demonstrate a tight coupling between the metabolic state and epigenetic aging. CryoEM of endogenous mammalian V-ATPase interacting with the TLDc protein mEAK-7 V-ATPases are rotary proton pumps that serve as signaling hubs with numerous proposed binding partners in cells. We used cryoEM combined with exhaustive focused classification to detect endogenous proteins that associate with V-ATPase from porcine kidney. A super-stoichiometric copy of subunit C was found in ~3% of complexes, while an additional ~1.6% of complexes bound mEAK-7, a protein with proposed roles in dauer formation in nematodes and mTOR signaling in mammals. High-resolution cryoEM of porcine kidney V-ATPase with recombinant mEAK-7 shows that mEAK-7s TLDc domain interacts with V-ATPases stator while its C-terminal  helix binds V-ATPases rotor. This crosslink would be expected to inhibit rotary catalysis. However, unlike inhibition of yeast V-ATPase by the TLDc protein Oxr1p, exogenous mEAK-7 does not inhibit V-ATPase and mEAK-7 overexpression in cells does not alter lysosomal or phagosomal pH. Instead, cryoEM suggests that interaction of mEAK-7 with V-ATPase is disrupted by ATP-induced rotation of the rotor. Comparison of Oxr1p and mEAK-7 binding explains this difference. Together, these results show that differences in V-ATPase binding by TLDc domain-containing proteins can lead to effects ranging from strong inhibition to formation of labile interactions that are sensitive to the enzymes activity. 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.