Evidence for a hydrogen sulfide-sensing E3 ligase in yeast In yeast, control of sulfur amino acid metabolism relies upon Met4, a transcription factor which activates the expression of a network of enzymes responsible for the biosynthesis of cysteine and methionine. In times of sulfur abundance, the activity of Met4 is repressed via ubiquitination by the SCFMet30 E3 ubiquitin ligase, but the mechanism by which the F-box protein Met30 senses sulfur status to tune its E3 ligase activity remains unresolved. Herein, we show that Met30 responds to flux through the transsulfuration pathway to regulate the MET gene transcriptional program. In particular, Met30 is responsive to the biological gas hydrogen sulfide, which is sufficient to induce ubiquitination of Met4 in vivo. Additionally, we identify important cysteine residues in Met30s WD-40 repeat region that sense the availability of sulfur in the cell. Our findings reveal how SCFMet30 dynamically senses the flow of sulfur metabolites through the transsulfuration pathway to regulate synthesis of these special amino acids.
Stress-dependent cell stiffening by tardigrade tolerance proteins through reversible formation of cytoskeleton-like filamentous network and gel-transition Tardigrades are able to tolerate almost complete dehydration by entering a reversible ametabolic state called anhydrobiosis and resume their animation upon rehydration. Dehydrated tardigrades are exceptionally stable and withstand various physical extremes. Although trehalose and late embryogenesis abundant (LEA) proteins have been extensively studied as potent protectants against dehydration in other anhydrobiotic organisms, tardigrades produce high amounts of tardigrade-unique protective proteins. Cytoplasmic-abundant heat-soluble (CAHS) proteins are uniquely invented in the lineage of eutardigrades, a major class of the phylum Tardigrada and are essential for their anhydrobiotic survival. However, the precise mechanisms of their action in this protective role are not fully understood. In the present study, we first postulated the presence of tolerance proteins that form protective condensates via phase separation in a stress-dependent manner and searched for tardigrade proteins that reversibly form condensates upon dehydration-like stress. Through comprehensive analysis, we identified 336 such proteins, collectively dubbed "Desolvataion-induced ReversiblY condensing Proteins (DRYPs)". Unexpectedly, we rediscovered CAHS proteins as highly enriched in DRYPs, 3 of which were major components of DRYPs. We revealed that these CAHS proteins reversibly polymerize into many cytoskeleton-like filaments depending on hyperosmotic stress in cultured cells and undergo reversible gel-transition in vitro. CAHS filamentation increases cell stiffness to resist deformation and improves resistance to dehydration-like stress. The conserved putative helical C-terminal region is necessary and sufficient for filament formation by CAHS proteins, and mutations disrupting the secondary structure of this region impaired both the filament formation and the gel transition. On the basis of these results, we propose that CAHS proteins are novel cytoskeletal proteins that form filamentous networks and undergo gel-transition in a stress-dependent manner to provide on-demand physical stabilization of cell integrity against deformative forces during dehydration and also contribute to the exceptional physical stability in a dehydrated state.
Hepatic MIR20B promotes nonalcoholic fatty liver disease by suppressing PPARA BackgroundNon-alcoholic fatty liver disease (NAFLD) is characterized by excessive lipid accumulation and imbalances in lipid metabolism in the liver. Although nuclear receptors (NRs) play a crucial role in hepatic lipid metabolism, the underlying mechanisms of NR regulation in NAFLD remain largely unclear.

MethodsUsing network analysis and RNA-seq to determine the correlation between NRs and microRNA in human NAFLD patients, we revealed that MIR20B specifically targets PPARA. MIR20B mimic and anti-MIR20B were administered to human HepG2 and Huh-7 cells and mouse primary hepatocytes as well as high fat diet (HFD)- or methionine-deficient diet (MCD)-fed mice to verify the specific function of MIR20B in NAFLD. We tested the inhibition of the therapeutic effect of a PPAR agonist, fenofibrate, by MIR20B and the synergic effect of combination of fenofibrate with anti-MIR20B.

ResultsWe revealed that MIR20B specifically targets PPARA through miRNA regulatory network analysis of nuclear receptor genes in NAFLD. The expression of MIR20B was upregulated in free fatty acid (FA)- treated hepatocytes and the livers of both obesity-induced mice and NAFLD patients. Overexpression of MIR20B significantly increased hepatic lipid accumulation and triglyceride levels. Furthermore, MIR20B significantly reduced FA oxidation and mitochondrial biogenesis by targeting PPARA. In MIR20B-introduced mice, the effect of fenofibrate to ameliorate hepatic steatosis was significantly suppressed. Finally, inhibition of MIR20B significantly increased FA oxidation and uptake, resulting in improved insulin sensitivity and a decrease in NAFLD progression. Moreover, combination of fenofibrate and anti-MIR20B exhibited the synergic effect on improvement of NAFLD in MCD-fed mice.

ConclusionsTaken together, our results demonstrate that the novel MIR20B targets PPARA, plays a significant role in hepatic lipid metabolism, and present an opportunity for the development of novel therapeutics for NAFLD.

FundingThis research was funded by Korea Mouse Phenotyping Project (2016M3A9D5A01952411), the National Research Foundation of Korea (NRF) grant funded by the Korea government (2020R1F1A1061267, 2018R1A5A1024340, NRF-2021R1I1A2041463, 2020R1I1A1A01074940), and the Future-leading Project Research Fund (1.210034.01) of UNIST. The response network of HSP70 defines vulnerabilities in cancer cells with the inhibited proteasome Human neoplasias are often addicted to the proteasome machinery. However, cancers have evolved efficient response mechanisms to overcome proteasome inhibition with bortezomib and carfilzomib - drugs approved for multiple myeloma treatment. To understand these responses we investigated proteome changes upon the proteasome inhibition with carfilzomib - in multiple myeloma, normal fibroblasts, and cancers of lung, colon, and pancreas. A pathway-oriented siRNA screen based on the proteomics results showed that molecular chaperones, autophagy- and endocytosis-related proteins are cancer-specific vulnerabilities combined with carfilzomib. Targeting of HSPA1A/B (HSP70 family chaperones) most specifically sensitized cancer cells and patient-derived organoids to the proteasome inhibition. A high level of HSPA1A/B mRNA correlated with a low proteasome activity in cancer patient tissues and is a risk factor in cancer patients with a low proteasome expression. Mechanistically, HSPA1A/B governed autophagy, unfolded protein response, endocytic trafficking, and chaperoned the proteasome machinery, suppressing the effect of the proteasome inhibition, but did not control the NRF1/2-driven proteasome subunit transcriptional bounce-back. Consequently, downregulation of NRF1 most specifically decreased the viability of cancer cells with the inhibited proteasome and HSP70. Dihydroceramide desaturase promotes the formation of intraluminal vesicles and inhibits autophagy to increase exosome production Exosomes are important for cell-cell communication. Deficiencies in the human dihydroceramide desaturase gene, DEGS1, increase the dihydroceramide-to-ceramide ratio and causes hypomyelinating leukodystrophy. However, the disease mechanism remains unknown. Here, we developed an in vivo assay with spatially controlled expression of exosome markers in Drosophila eye imaginal discs and showed that the level and activity of the DEGS1 ortholog, ifc, correlated with exosome production. Knocking out ifc decreased the density of the exosome precursor intraluminal vesicles (ILVs) in the multivesicular endosomes and reduced the number of exosomes released. While ifc overexpression and autophagy inhibition both enhanced exosome production, combining the two had no additive effect. Moreover, DEGS1 activity was sufficient to drive ILV formation in vitro. Together, DEGS1/Ifc controls the dihydroceramide-to-ceramide ration and enhances exosome secretion by promoting ILV formation and preventing the autophagic degradation of MVEs.

These findings provide a potential cause for the neuropathy associated with DEGS1-deficient mutations.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=115 SRC="FIGDIR/small/376046v3_ufig1.gif" ALT="Figure 1">
View larger version (43K):
org.highwire.dtl.DTLVardef@c7479corg.highwire.dtl.DTLVardef@12e7f46org.highwire.dtl.DTLVardef@1a5e7f0org.highwire.dtl.DTLVardef@11e7f08_HPS_FORMAT_FIGEXP  M_FIG C_FIG HighlightsO_LIAn in vivo system was developed for observing exosome production in Drosophila.
C_LIO_LIDihydroceramide desaturase (DEGS1/Ifc) promotes exosome production by two means.
C_LIO_LIIfc drives membrane invagination to promote the formation of intraluminal vesicles.
C_LIO_LIIfc inhibits autophagic degradation of MVEs and increases exosome release.
C_LI

eTOC BlurbThe level and activity of dihydroceramide desaturase (human DEGS1 and Drosophila Ifc) correlate with exosome production. Wu et al. show that DEGS1 drives the formation of intraluminal vesicles in vivo and in vitro. Overexpressing Ifc inhibits autophagy and reduces the degradation of multivesicular endosomes, thus increases exosome release. A flux-based machine learning model to simulate the impact of pathogen metabolic heterogeneity on drug interactions AO_SCPLOWBSTRACTC_SCPLOWDrug combinations are a promising strategy to counter antibiotic resistance. However, current experimental and computational approaches do not account for the entire complexity involved in combination therapy design, such as the effect of pathogen metabolic heterogeneity, changes in the growth environment, drug treatment order and time interval. To address these limitations, we present a comprehensive approach that uses genome-scale metabolic modeling and machine learning to guide combination therapy design. Our mechanistic approach (a) accommodates diverse data types, (b) accounts for time- and order-specific interactions, and (c) accurately predicts drug interactions in various growth conditions and their robustness to pathogen metabolic heterogeneity. Our approach achieved high accuracy (AUC = 0.83 for synergy, AUC = 0.98 for antagonism) in predicting drug interactions for E. coli cultured in 57 metabolic conditions based on experimental validation. The entropy in bacterial metabolic response was predictive of combination therapy outcomes across time scales and growth conditions. Simulation of metabolic heterogeneity using population-FBA identified two sub-populations of E. coli cells defined by regulation of cis- to trans-isomerization of unsaturated fatty acids with distinct effects on drug interactions. Analysis of the vast landscape of condition-specific drug interactions revealed a set of 24 robustly synergistic drug combinations with potential for clinical use.

SO_SCPLOWIGNIFICANCEC_SCPLOWWorldwide, 700,000 people die each year from drug-resistant infections. Drug combinations have great potential to reduce the spread of drug-resistant bacteria. However, their potency is impacted by both the pathogen growth environment and the heterogeneity in pathogen metabolism. The metabolic heterogeneity in a pathogen population allows them to survive antibiotic treatment. Here we present a flexible machine-learning framework that utilizes diverse data types to effectively search through the large design space of both sequential and simultaneous combination therapies across hundreds of simulated growth conditions and pathogen metabolic states. Our approach can serve as a useful guide for the selection of robustly synergistic drug combinations. Strong and localized recurrence controlsdimensionality of neural activity across brainareas The brain contains an astronomical number of neurons, but it is their collective activity that underlies brain function. The number of degrees of freedom that this collective activity explores - its dimensionality - is therefore a fundamental signature of neural dynamics and computation (1-7). However, it is not known what controls this dimensionality in the biological brain - and in particular whether and how local synaptic networks play a role (8-10). Through analysis of high-density Neuropixels recordings (11), we argue that areas across the mouse cortex operate in a sensitive regime that gives these synaptic networks a very strong role in controlling dimensionality. Moreover, we show that this control is expressed through highly tractable features of these synaptic networks. We then analyze these key features via a massive synaptic physiology dataset (12). Quantifying these features in terms of cell-type specific network motifs, we find that the synaptic patterns that impact dimensionality are prevalent in both mouse and human brains. Thus local circuitry scales up systematically to help control the degrees of freedom that brain networks may explore and exploit. X upregulation is not global and extent of upregulation differs between ancestral and acquired X-linked genes Evolution of sex chromosome dosage compensation in mammals remains poorly understood. Ohnos hypothesis state that evolution of dosage compensation occurred through two steps: first, to compensate the dosage imbalance created due to the degradation of Y chromosome in male, upregulation of X-chromosome happened. Subsequently, transmission of X-chromosome upregulation (XCU) in female led to the evolution of X-chromosome inactivation (XCI) to counteract extra dosage of X-linked genes in female cells. Here, we have profiled gene-wise dynamics of XCU in pre-gastrulation mouse embryos at single cell level and found that XCU is dynamically linked with XCI, however, XCU is not global or chromosome-wide like XCI. Therefore, our result raises question whether XCU driven the evolution of XCI. If so, then why XCI is chromosome wide while XCU is not. We propose that XCI might have evolved independent of XCU and therefore refining the current model is necessary. Separately, we show that higher occupancy of different activating factors at upregulated X-linked genes leads to enhanced transcriptional burst frequency and thereby leads to upregulation. On the other hand, our analysis indicates that extent of upregulation, enrichment of different activating marks differs between ancestral and newly acquired X-linked genes. Altogether, our study provides significant insight into the dynamics and mechanistic basis of evolution of sex chromosome dosage compensation.

O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY Induced antigen-binding polyreactivity in human serum IgA Previous studies have shown that polyreactive antibodies play an important role in the frontline defense against the dissemination of pathogens in the pre-immune host. Interestingly, antigen-binding polyreactivity can not only be inherent, but also acquired post-translationally. The ability of individual monoclonal IgG and IgE antibodies to acquire polyreactivity following contact with various agents that destabilize protein structure (urea, low pH) or have pro-oxidative potential (heme, ferrous ions) has been studied in detail. However, to the best of our knowledge this property of human IgA has previously been described only cursorily. In the present study pooled human serum IgA and two human monoclonal IgA antibodies were exposed to buffers with acidic pH, to free heme or to ferrous ions, and the antigen-binding behavior of the native and modified IgA to viral and bacterial antigens were compared using immunoblot and ELISA. We observed a dose-dependent increase in reactivity to several bacterial extracts and to pure viral antigens. This newly described property of IgA may have therapeutic potential as has already been shown for pooled IgG with induced polyreactivity. Evaluation of Antibiotic Tolerance in Pseudomonas aeruginosa for Aminoglycosides and its Predicted Gene Regulations Through In-silico Transcriptomic Analysis Pseudomonas aeruginosa causes chronic infections like cystic fibrosis, endocarditis, bacteremia and sepsis, which are life-threatening and difficult to treat. The lack of antibiotic response in P. aeruginosa is due to adaptive resistance mechanism, which prevents the entry of antibiotics into cytosol of the cell to achieve tolerance. Among different groups of antibiotics, aminoglycosides are used as a parental antibiotic for treatment of P. aeruginosa. This study aims to determine the kinetics of antibiotic tolerance and gene expression changes in P. aeruginosa exposed to amikacin, gentamicin, and tobramycin. These antibiotics were exposed to P. aeruginosa at their MICs and the experimental setup was monitored till 72 hours, followed by the measurement of optical density in the interval of every 12 hours. The growth of P. aeruginosa in MICs of antibiotics represents the kinetics of antibiotic tolerance in amikacin, gentamicin, and tobramycin. Transcriptomic profile of antibiotic exposed P. aeruginosa PA14 was taken from Gene Expression Omnibus (GEO), NCBI as microarray datasets. The gene expressions of two datasets were compared by test versus control. Tobramycin exposed P. aeruginosa failed to develop tolerance in MICs 0.5{micro}g/mL, 1{micro}g/mL and 1.5{micro}g/mL. Whereas amikacin and gentamicin treated P. aeruginosa developed tolerance in MICs. This depicts the superior in vitro response of tobramycin over the gentamicin and amikacin. Further, in silico transcriptomic analysis of tobramycin treated P. aeruginosa resulted in low expression of 16s rRNA Methyltransferase E, B & L, alginate biosynthesis genes and several proteins of Type 2 Secretory System (T2SS) and Type 3 Secretory System (T3SS). The Differentially Expressed Genes (DEGs) of alginate biosynthesis, and RNA Methyltransferases suggests increased antibiotic response and low probability of developing resistance. The use of tobramycin as a parental antibiotic with its synergistic combination might combat P. aeruginosa with increased response. Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. Robust and Remote Measurement of Heart Rate Based on A Surveillance Camera Arrhythmia is a marked symptom of many cardiovascular diseases (CVDs). The accurate and timely detection of heart rate can greatly reduce the harm of arrhythmia to people. However, it is still a challenge to robustly and remotely measure heart rate in daily life due to the changing environmental conditions during measurement, such as the varying light intensity, the movement of people, and the uncertain distance between the sensor to people. In this study, we propose a method to accurately measure human heart rate within a distance of 4.5 meters under different light intensities by simply using a surveillance camera. After a 20-second color video of a persons hand is captured by the camera, a method based on the Fast Fourier Transform (FFT) algorithm is designed to extract the blood volume pulse wave to calculate the heart rate. According to the comparison between the real heart rate and results measured by electrocardiography (ECG), the proposed method achieves an accuracy of 95.8% when the measurement is performed within a distance of 4.5 meters and 90% when within 5.0 meters. Our experiments show that when the illuminance varies between 100-1000 lux (lighting level indoor), we still get the correct results. Our experiments also demonstrate that the proposed method accurately obtain heart rate even when the light intensity is below 32 lux (300-500 lux in a workplace environment). The methods strong adaptability to changing environmental conditions makes it applicable to many scenarios, such as homes of the elderly, classrooms, and other public spaces. Low photorespiratory capacity is sufficient to block the induction of photosynthetic capacity in high light The induction of high photosynthetic capacity in high light (HL) is a common response among many herbaceous dicot plants, however, the signals that control this response remain largely unknown. Here, multiple independent lines of evidence are presented in support of the conclusion that low photorespiratory capacity acts a negative signal to limit photosynthetic capacity acclimation in HL in Arabidopsis thaliana. Using a panel of natural accessions, primary nitrogen (N) assimilation and photorespiration rates early after a shift to growth in HL, as well as activities for key enzymes in these pathways, were shown to positively correlate with the magnitude of the subsequent induction of photosynthetic capacity, which occurred several days later. Time-resolved metabolomic data during acclimation to HL were collected using a strongly acclimating ecotype and a weakly acclimating ecotype, revealing in greater detail the differences in N assimilation, photorespiration, and triose-phosphate utilization pathways underlying efficient photosynthetic capacity acclimation. When shifted into HL growth conditions under non-photorespiratory conditions, weakly acclimating ecotypes and even photorespiratory mutants gained the ability to strongly induce high photosynthetic capacity in HL. Thus, a negative, photorespiration-dependent signal early in the HL shift appears to block photosynthetic capacity acclimation in accessions with low photorespiratory capacity, whereas accessions with high photorespiratory capacity are licensed to increase photosynthetic capacity. 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. Stimulus manipulations permit activation of fiber subpopulations in the mouse and rat vagus Afferent and efferent vagal fibers mediate bidirectional communication between the brain and visceral organs. Small, unmyelinated C-afferents constitute the majority of vagal fibers, play critical roles in numerous interoceptive circuits and autonomic reflexes in health and disease and may contribute to the efficacy and safety of vagus nerve stimulation (VNS). Selective engagement of C-afferents with electrical stimuli has not been feasible, due to the default fiber recruitment order: larger fibers first, smaller fibers last. Here, we determine and optimize an electrical stimulus that selectively engages vagal C-afferents. Intermittent KHz-frequency electrical stimulation (KES) activates motor and, preferentially, sensory vagal neurons in the brainstem. During KES, asynchronous activity of C-afferents increases, while that of larger fibers remains largely unchanged. In parallel, KES effectively blocks excitability of larger fibers while moderately suppressing excitability of C-afferents. By compiling selectivity indices in individual animals, we find that optimal KES parameters for C-afferents are >5KHz frequency and 7-10 times engagement threshold (xT) intensity in rats, 15-25xT in mice. These effects can be explained in computational models by how sodium channel responses to KES are shaped by axonal size and myelin. Our results indicate that selective engagement of vagal C-afferents is attainable by intermittent KES. Tetraspanin Cd9b and Cxcl12a/Cxcr4b have a synergistic effect on the control of collective cell migration Collective cell migration is essential for embryonic development and homeostatic processes. During zebrafish development, the posterior lateral line primordium (pLLP) navigates along the embryo flank by collective cell migration. The chemokine receptors, Cxcr4b and Cxcr7b, as well as their cognate ligand, Cxcl12a, are essential for this process. We corroborate that knockdown of the zebrafish cd9 tetraspanin orthologue, cd9b, results in mild pLL abnormalities. Through generation of CRISPR and TALEN mutants, we show that cd9a and cd9b function partially redundantly in pLLP migration, which is delayed in the cd9b single and cd9a; cd9b double mutants. This delay led to a transient reduction in neuromast numbers. Loss of both Cd9a and Cd9b sensitized embryos to reduced Cxcr4b and Cxcl12a levels. Together these results provide evidence that Cd9 modulates collective cell migration of the pLLP during zebrafish development. One interpretation of these observations is that Cd9 contributes to more effective chemokine signalling. Reconstitution of minimal motility system based on Spiroplasma swimming by expressing two bacterial actins in synthetic minimal bacterium Motility is one of the most important features of life, but its evolutionary origin remains unknown. In this study, we focused on Spiroplasma, commensal, or parasitic bacteria. They swim by switching the helicity of a ribbon-like cytoskeleton that comprises six proteins, each of which evolved from a nucleosidase and bacterial actin called MreB. We expressed these proteins in a synthetic, non-motile minimal bacterium, JCVI-syn3.0B, whose reduced genome was computer-designed and chemically synthesized. The synthetic bacterium exhibited swimming motility with features characteristic of Spiroplasma swimming. Moreover, some combinations of the two proteins produced a helical cell shape and swimming, suggesting that the swimming originated from the differentiation and coupling of bacterial actins, and we also obtained a minimal system for motility of the synthetic bacterium.

One-Sentence SummaryThe minimal system comprised two bacterial actins that provided cell helicity and swimming to the synthetic minimal bacterium. High-quality Arabidopsis thaliana Genome Assembly with Nanopore and HiFi Long Reads Arabidopsis thaliana is an important and long-established model species for plant molecular biology, genetics, epigenetics, and genomics. However, the latest version of reference genome still contains significant number of missing segments. Here, we report a high-quality and almost complete Col-0 genome assembly with two gaps (Col-XJTU) using combination of Oxford Nanopore Technology ultra-long reads, PacBio high-fidelity long reads, and Hi-C data. The total genome assembly size is 133,725,193 bp, introducing 14.6 Mb of novel sequences compared to the TAIR10.1 reference genome. All five chromosomes of Col-XJTU assembly are highly accurate with consensus quality (QV) scores > 60 (ranging from 62 to 68), which are higher than those of TAIR10.1 reference (QV scores ranging from 45 to 52). We have completely resolved chromosome (Chr) 3 and Chr5 in a telomere-to-telomere manner. Chr4 has been completely resolved except the nucleolar organizing regions, which comprise long repetitive DNA fragments. The Chr1 centromere (CEN1), reportedly around 9 Mb in length, is particularly challenging to assemble due to the presence of tens of thousands of CEN180 satellite repeats. Using the cutting-edge sequencing data and novel computational approaches, we assembled about 4 Mb of sequence for CEN1 and a 3.5-Mb-long CEN2. We investigated the structure and epigenetics of centromeres. We detected four clusters of CEN180 monomers, and found that the centromere-specific histone H3-like protein (CENH3) exhibits a strong preference for CEN180 cluster 3. Moreover, we observed hypomethylation patterns in CENH3-enriched regions. We believe that this high-quality genome assembly, Col-XJTU, would serve as a valuable reference to better understand the global pattern of centromeric polymorphisms, as well as genetic and epigenetic features in plants. Rossmann-toolbox: a deep learning-based protocol for the prediction and design of cofactor specificity in Rossmann-fold proteins The Rossmann fold enzymes are involved in essential biochemical pathways such as nucleotide and amino acid metabolism. Their functioning relies on interaction with cofactors, small nucleoside-based compounds specifically recognized by a conserved {beta}{beta} motif shared by all Rossmann fold proteins. While Rossmann methyltransferases recognize only a single cofactor type, the S-Adenosylmethionine (SAM), the oxidoreductases, depending on the family, bind nicotinamide (NAD, NADP) or flavin-based (FAD) cofactors. In this study, we show that despite its short length, the {beta}{beta} motif unambiguously defines the specificity towards the cofactor. Following this observation, we trained two complementary deep learning models for the prediction of the cofactor specificity based on the sequence and structural features of the {beta}{beta} motif. A benchmark on two independent test sets, one containing {beta}{beta} motifs bearing no resemblance to those of the training set, and the other comprising 38 experimentally confirmed cases of rational design of the cofactor specificity, revealed the nearly perfect performance of the two methods. The Rossmann-toolbox protocols can be accessed via the webserver at https://lbs.cent.uw.edu.pl/rossmann-toolbox and are available as a Python package at https://github.com/labstructbioinf/rossmann-toolbox.

Key pointsO_LIThe Rossmann fold encompasses a multitude of diverse enzymes involved in most of the essential cellular pathways
C_LIO_LIProteins belonging to the Rossmann fold co-evolved with their nucleoside-based cofactors and require them for the functioning
C_LIO_LIManipulating the cofactor specificity is an important step in the process of enzyme engineering
C_LIO_LIWe developed an end-to-end pipeline for the prediction and design of the cofactor specificity of the Rossmann fold proteins
C_LIO_LIOwing to the utilization of deep learning approaches the pipeline achieved nearly perfect accuracy
C_LI No impact of biocontrol agent's predation cues on development time or size of surviving Aedes albopictus under optimal nutritional availability Cyclopoid copepods have been applied successfully to limit populations of highly invasive Aedes albopictus mosquitoes that can transmit diseases of public health importance. However, there is concern that changes in certain mosquito traits, induced by exposure to copepod predation, might increase the risk of disease transmission. In this study, third instar Ae. albopictus larvae, hereafter referred to as "focal individuals", were exposed to Megacyclops viridis predator cues associated with both consumption of newly-hatched mosquito larvae and attacks on focal individuals. The number of newly-hatched larvae surrounding each focal larva was held constant to control for density effects on size, and the focal individuals day of pupation and wing length were recorded for each replicate.

Exposing late instar Ae. albopictus to predation decreased their chances of surviving to adulthood, and three focal larvae that died in the predator treatment showed signs of melanization, indicative of wounding. Among surviving focal Ae. albopictus, no significant difference in either pupation day or wing length was observed due to copepod predation. The absence of a significant sublethal impact from M. viridis copepod predation on surviving later-stage larvae in this analysis supports the use of M. viridis as a biocontrol agent against Ae. albopictus.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=106 SRC="FIGDIR/small/460697v2_ufig1.gif" ALT="Figure 1">
View larger version (9K):
org.highwire.dtl.DTLVardef@1abaabborg.highwire.dtl.DTLVardef@1561651org.highwire.dtl.DTLVardef@d325a7org.highwire.dtl.DTLVardef@4d1830_HPS_FORMAT_FIGEXP  M_FIG C_FIG Social behavioral profiling by unsupervised deep learning reveals a stimulative effect of dopamine D3 agonists on zebrafish sociality It has been a major challenge to systematically evaluate and compare how pharmacological perturbations influence social behavioral outcomes. Although some pharmacological agents are known to alter social behavior, precise description and quantification of such effects have proven difficult. The complexity of brain functions regulating sociality makes it challenging to predict drug effects on social behavior without testing in live animals, and most existing behavioral assays are low-throughput and provide only unidimensional readouts of social function. To achieve richer characterization of drug effects on sociality, we developed a scalable social behavioral assay for zebrafish named ZeChat based on unsupervised deep learning. High-dimensional and dynamic social behavioral phenotypes are automatically classified using this method. By screening a neuroactive compound library, we found that different classes of chemicals evoke distinct patterns of social behavioral fingerprints. By examining these patterns, we discovered that dopamine D3 agonists possess a social stimulative effect on zebrafish. The D3 agonists pramipexole, piribedil, and 7-hydroxy-DPAT-HBr rescued social deficits in a valproic acid-induced zebrafish autism model. The ZeChat platform provides a promising approach for dissecting the pharmacology of social behavior and discovering novel social-modulatory compounds. Comparative toxicogenomics of glyphosate and Roundup herbicides by mammalian stem cell-based genotoxicity assays and molecular profiling in Sprague-Dawley rats Whether glyphosate-based herbicides (GBHs) are more potent than glyphosate alone at activating cellular mechanisms, which drive carcinogenesis remains controversial. As GBHs are more cytotoxic that glyphosate, we reasoned they may also be more capable of activating carcinogenic pathways. We tested this hypothesis by comparing the effects of glyphosate with Roundup GBHs both in vitro and in vivo. First, glyphosate was compared with representative GBHs namely MON 52276 (EU), MON 76473 (UK) and MON 76207 (USA) using the mammalian stem cell-based ToxTracker system. Here, MON 52276 and MON 76473, but not glyphosate and MON 76207, activated oxidative stress and unfolded protein responses. Second, molecular profiling of liver was performed in female Sprague-Dawley rats exposed to glyphosate or MON 52276 (both at 0.5, 50, 175 mg/kg bw/day glyphosate) for 90 days. MON 52276 but not glyphosate increased hepatic steatosis and necrosis. MON 52276 and glyphosate altered the expression of genes in liver reflecting TP53 activation by DNA damage and circadian rhythm regulation. Genes most affected in liver were similarly altered in kidneys. Small RNA profiling in liver showed decreased amounts of miR-22 and miR-17 from MON 52276 ingestion. Glyphosate decreased mir-30 while miR-10 levels were increased. DNA methylation profiling of liver revealed 5,727 and 4,496 differentially methylated CpG sites between the control and glyphosate and MON 52276 exposed animals respectively. Apurinic/apyrimidinic DNA damage formation in liver was increased with glyphosate exposure. Altogether, our results show that Roundup formulations cause more biological changes linked with carcinogenesis than glyphosate. 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. 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. Negative DNA supercoiling makes protein-mediated looping deterministic and ergodic within the bacterial doubling time Protein-mediated DNA looping is fundamental to gene regulation and such loops occur stochastically in purified systems. Additional proteins increase the probability of looping, but these probabilities maintain a broad distribution. For example, the probability of lac repressor-mediated looping in individual molecules ranged 0-100%, and individual molecules exhibited representative behavior only in observations lasting an hour or more. Titrating with HU protein progressively compacted the DNA without narrowing the 0-100% distribution. Increased negative supercoiling produced an ensemble of molecules in which all individual molecules more closely resembled the average. Furthermore, in only twelve minutes of observation, well within the doubling time of the bacterium, most molecules exhibited the looping probability of the ensemble. DNA supercoiling, an inherent feature of all genomes, appears to impose time-constrained, emergent behavior on otherwise random molecular activity. On the incongruence of genotype-phenotype and fitness landscapes The mapping from genotype to phenotype to fitness typically involves multiple nonlinearities that can transform the effects of mutations. For example, mutations may contribute additively to a phenotype, but their effects on fitness may combine non-additively because selection favors a low or intermediate value of that phenotype. This can cause incongruence between the topographical properties of a fitness landscape and its underlying genotype-phenotype landscape. Yet, genotype-phenotype landscapes are often used as a proxy for fitness landscapes to study the dynamics and predictability of evolution. Here, we use theoretical models and empirical data on transcription factor-DNA interactions to systematically study the incongruence of genotype-phenotype and fitness landscapes when selection favors a low or intermediate phenotypic value. Using the theoretical models, we prove a number of fundamental results. For example, selection for low or intermediate phenotypic values does not change simple sign epistasis into reciprocal sign epistasis, yet it changes reciprocal sign epistasis into simple sign epistasis and no sign epistasis with equal probability. More broadly, we show that such selection tends to create fitness landscapes that are more rugged than the underlying genotype-phenotype landscape, but this increased ruggedness typically does not frustrate adaptive evolution because the local adaptive peaks in the fitness landscape tend to be nearly as tall as the global peak. Many of these results carry forward to the empirical genotype-phenotype landscapes, which may help to explain why low- and intermediate-affinity transcription factor-DNA interactions are so prevalent in eukaryotic gene regulation.

AUTHOR SUMMARYHow do mutations change phenotypic traits and organismal fitness? This question is often addressed in the context of a classic metaphor of evolutionary theory -- the fitness landscape. A fitness landscape is akin to a physical landscape, in which genotypes define spatial coordinates, and fitness defines the elevation of each coordinate. Evolution then acts like a hill-climbing process, in which populations ascend fitness peaks as a consequence of mutation and selection. It is becoming increasingly common to construct such landscapes using experimental data from high-throughput sequencing technologies and phenotypic assays, in systems such as macromolecules and gene regulatory circuits. Although these landscapes are typically defined by molecular phenotypes, and are therefore more appropriately referred to as genotype-phenotype landscapes, they are often used to study evolutionary dynamics. This requires the assumption that the molecular phenotype is a reasonable proxy for fitness, which need not be the case. For example, selection may favor a low or intermediate phenotypic value, causing incongruence between a fitness landscape and its underlying genotype-phenotype landscape. Here, we study such incongruence using a diversity of theoretical models and experimental data from gene regulatory systems. We regularly find incongruence, in that fitness landscapes tend to comprise more peaks than their underlying genotype-phenotype landscapes. However, using evolutionary simulations, we show that this increased ruggedness need not impede adaptation.