Distinct mutations and lineages of SARS-CoV-2 virus in the early phase of COVID-19 pandemic and subsequent global expansion A novel coronavirus, SARS-CoV-2, has caused over 190 million cases and over 4 million deaths worldwide since it occurred in December 2019 in Wuhan, China. Here we conceptualized the temporospatial evolutionary and expansion dynamics of SARS-CoV-2 by taking a series of cross-sectional view of viral genomes from early outbreak in January 2020 in Wuhan to early phase of global ignition in early April, and finally to the subsequent global expansion by late December 2020. Based on the phylogenetic analysis of the early patients in Wuhan, Wuhan/WH04/2020 is supposed to be a more appropriate reference genome of SARS-CoV-2, instead of the first sequenced genome Wuhan-Hu-1. By scrutinizing the cases from the very early outbreak, we found a viral genotype from the Seafood Market in Wuhan featured with two concurrent mutations (i.e. M type) had become the overwhelmingly dominant genotype (95.3%) of the pandemic one year later. By analyzing 4,013 SARS-CoV-2 genomes from different continents by early April, we were able to interrogate the viral genomic composition dynamics of initial phase of global ignition over a timespan of 14-week. 11 major viral genotypes with unique geographic distributions were also identified. WE1 type, a descendant of M and predominantly witnessed in western Europe, consisted a half of all the cases (50.2%) at the time. The mutations of major genotypes at the same hierarchical level were mutually exclusive, which implying that various genotypes bearing the specific mutations were propagated during human-to-human transmission, not by accumulating hot-spot mutations during the replication of individual viral genomes. As the pandemic was unfolding, we also used the same approach to analyze 261,323 SARS-CoV-2 genomes from the world since the outbreak in Wuhan (i.e. including all the publicly available viral genomes) in order to recapitulate our findings over one-year timespan. By 25 December 2020, 95.3% of global cases were M type and 93.0% of M-type cases were WE1. In fact, at present all the four variants of concern (VOC) are the descendants of WE1 type. This study demonstrates the viral genotypes can be utilized as molecular barcodes in combination with epidemiologic data to monitor the spreading routes of the pandemic and evaluate the effectiveness of control measures. Moreover, the dynamics of viral mutational spectrum in the study may help the early identification of new strains in patients to reduce further spread of infection, guide the development of molecular diagnosis and vaccines against COVID-19, and help assess their accuracy and efficacy in real world at real time.
Migrators Within Migrators: Exploring Transposable Element Landscapes, Activity, and Host Gene Associations in the Monarch Butterfly, Danaus plexippus BackgroundLepidoptera (butterflies and moths) are an important model system in ecology and evolution. A high-quality chromosomal genome assembly is available for the monarch butterfly (Danaus plexippus), but it lacks an in-depth transposable element (TE) annotation, presenting an opportunity to explore monarch TE dynamics and the impact of TEs on shaping the monarch genome.

ResultsWe find 6.21% of the monarch genome is comprised of TEs, a reduction of 6.85% compared to the original TE annotation performed on the draft genome assembly. Monarch TE content is low compared to two closely related species with available genomes, Danaus chrysippus (33.97% TE) and Danaus melanippus (11.87% TE). The biggest TE contributions to genome size in the monarch are LINEs and Penelope-like elements, and three newly identified families, r2-hero_dPle (LINE), penelope-1_dPle (Penelope-like), and hase2-1_dPle (SINE), collectively contribute 34.92% of total TE content. We find evidence of recent TE activity, with two novel Tc1 families rapidly expanding over recent timescales (tc1-1_dPle, tc1- 2_dPle). LINE fragments show signatures of genomic deletions indicating a high rate of TE turnover. We investigate associations between TEs and wing colouration and immune genes and identify a three-fold increase in TE content around immune genes compared to other host genes.

ConclusionsWe provide a detailed TE annotation and analysis for the monarch genome, revealing a considerably smaller TE contribution to genome content compared to two closely related Danaus species with available genome assemblies. We identify highly successful novel DNA TE families rapidly expanding over recent timescales, and ongoing signatures of both TE expansion and removal highlight the dynamic nature of repeat content in the monarch genome. Our findings also suggest that insect immune genes are promising candidates for future interrogation of TE-mediated host adaptation.
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. Full-likelihood analysis of genomic data clarifies a complex history of species divergence and introgression in the erato-sara group of Heliconius butterflies AO_SCPLOWBSTRACTC_SCPLOWIntrogressive hybridization plays a key role in adaptive evolution and species diversification in many groups of species. However, frequent hybridization and gene flow between species makes estimation of the species phylogeny and key population parameters challenging. Here, we show that by accounting for phasing and using full-likelihood analysis methods, introgression histories and population parameters can be estimated reliably from whole-genome sequence data. We employ full-likelihood methods under the multispecies coalescent (MSC) model with and without gene flow to analyze the genomic data from six members of the erato-sara clade of Heliconius butterflies and infer the species phylogeny and cross-species introgression events. The methods naturally accommodate random fluctuations in genealogical history across the genome due to deep coalescence. To avoid heterozygote phasing errors in haploid sequences commonly produced by genome assembly methods, we process and compile unphased diploid sequence alignments and use analytical methods to average over uncertainties in heterozygote phase resolution. There is robust evidence for introgression across the genome, both among distantly related species deep in the phylogeny and between sister species in shallow parts of the tree. We obtain chromosome-specific estimates of key population parameters such as introgression directions, times and probabilities, as well as species divergence times and population sizes for modern and ancestral species. We confirm ancestral gene flow between the sara clade and an ancestral population of H. telesiphe, a likely hybrid speciation origin for H. hecalesia, and gene flow between sister species H. erato and H. himera. Inferred introgression among ancestral species also explains the history of two chromosomal inversions deep in the phylogeny of the group. This study illustrates how a full-likelihood approach based on the multispecies coalescent makes it possible to extract rich historical information of species divergence and gene flow from genomic data. Niche differences, not fitness differences, explain coexistence across ecological groups Understanding the drivers of species coexistence is an important objective in ecology. Yet, the multitude of methods to study coexistence hampers cross-community comparisons. Here, we standardized niche and fitness differences (i.e how species limit themselves compared to others and their competitive ability, respectively) across 1018 species pairs to investigate species coexistence across ecological groups and methodological settings (experimental setup, natural co-occurrence, population model used, and growth method). We find that, first, coexistence is driven by large niche differences, not by small fitness differences. Second, species group into clear clusters of coexisting and non-coexisting species along the niche axis. Finally, these clusters are not driven by ecological or methodological settings. This suggests differences between coexisting and non-coexisting communities transcending those measured in our empirical systems. Overall, our results show that species coexistence is mainly influenced by mechanisms acting on niche differences. Differential fate of acellular vascular scaffolds in vivo involves macrophages and T-cell subsets Biological scaffold or implant is a popular choice for the preparation of tissue-engineered organs and has the potential to address donor shortage in the clinics. However, biological scaffolds prepared by physical or chemical agents cause damage to the extracellular matrix by potentially inducing immune responses after implantation. The current study explores an alternative route for the preparation of acellular scaffolds and explores the fate of the prepared scaffolds in a milieu of immune cells following implantation without using immunosuppressant. Using the syngeneic (Lewis male-Lewis female) and allogeneic (Brown Norway male-Lewis female) models and different tissue routes (subcutaneous vs omentum) for transplantation, normal blood vascular scaffolds were implanted which was converted to acellular vascular scaffolds by in vivo natural decellularization at the end of 2 months of observation. We also prepared chemically decellularized acellular scaffolds from normal untreated blood vascular scaffolds using a cocktail of chemicals which was also similarly placed in subcutaneous and omentum sites. Here, we applied in-depth quantitative proteomics along with histology and image analysis to comprehensively describe and compare the proteome of the natural and chemically decellularized scaffold. Our data confirm that site-specific advantages exist in modulating the ECM and regulating the immune responses (macrophage and T cells) following implantation, which possibly led to the production of an acellular scaffold (natural decellularization) under in vivo conditions. The current approach opens up the possibility to create tailor-made acellular scaffolds to build functional blood vessels. In addition, the identification of different tissue sites facilitates differential immune response against the scaffolds. This study provides a rich resource aimed toward an enhanced mechanistic understanding to study immune responses under similar settings in the field of transplantation and regenerative medicine.

Impact statementThe development of a scaffold helps in the preparation of a functional organ in the clinics. In the current study, we prepared an acellular vascular scaffold by utilizing site specific tissue changes and vis-a-vis compared with a conventionally chemically prepared biological scaffold at genomic and protein level, which helped us to identify immunological trigger following implantation. The current study which was carried out without any immunosuppressive agents could help to establish (a) alternative strategies for preparing biological scaffolds as well as (b) implantable sites as potential bioreactors to circumvent any adverse immune reactions for acceptance of the scaffold/implant post implantation. Combinatorial drug-microenvironment interaction mapping reveals cell-extrinsic drug resistance mechanisms and clinically relevant patient subgroups in CLL The tumour microenvironment and genetic alterations collectively influence drug efficacy in cancer, but current evidence is limited to small scale studies and systematic analyses are lacking. We chose Chronic Lymphocytic Leukaemia (CLL), the most common leukaemia in adults, as a model disease to study this complex interplay systematically. We performed a combinatorial assay using 12 drugs individually co-applied with each of 17 microenvironmental stimuli in 192 primary CLL samples, generating a comprehensive map of drug-microenvironment interactions in CLL. This data was combined with whole-exome sequencing, DNA-methylation, RNA-sequencing and copy number variant annotation. Our assay identified four distinct CLL subgroups that differed in their responses to the panel of microenvironmental stimuli. These subgroups were characterized by distinct clinical outcomes independently of known prognostic markers. We investigated the effect of CLL- specific recurrent genetic alterations on microenvironmental responses and identified trisomy 12 as an amplifier of multiple microenvironmental stimuli. We further quantified the impact of microenvironmental stimuli on drug response, confirmed known interactions such as Interleukin (IL) 4 mediated resistance to B cell receptor (BCR) inhibitors, and identified new interactions such as Interferon-{gamma} induced resistance to BCR inhibitors. Finally, we identified interactions which were limited to genetic subgroups. Resistance to chemotherapeutics, such as Fludarabine, induced by Toll-Like Receptor (TLR) agonists could be observed in IGHV unmutated patient samples and IGHV mutated samples with trisomy 12. In-vivo relevance was investigated in CLL-infiltrated lymph nodes, which showed increased IL4 and TLR signalling activity compared to healthy samples (p<0.001). High IL4 activity in lymph nodes correlated with faster disease progression (p=0.038).

We provide a publicly available resource (www.dietrichlab.de/CLL_Microenvironment/) which uncovers tumour cell extrinsic influences on drug response and disease progression in CLL, and how these interactions are modulated by cell intrinsic molecular features. Prostaglandin EP3 Receptor signaling is required to prevent insulin hypersecretion and metabolic dysfunction in a non-obese mouse model of insulin resistance When homozygous for the LeptinOb mutation (Ob), Black-and-Tan Brachyury (BTBR) mice become morbidly obese and severely insulin resistant, and by 10 weeks of age, frankly diabetic. Previous work has shown Prostaglandin EP3 Receptor (EP3) expression and activity is up-regulated in islets from BTBR-Ob mice as compared to lean controls, actively contributing to their beta-cell dysfunction. In this work, we aimed to test the impact of beta-cell-specific EP3 loss on the BTBR-Ob phenotype by crossing Ptger3 floxed mice with the Rat insulin promoter (RIP)-CreHerr driver strain. Instead, germline recombination of the floxed allele in the founder mouse - an event whose prevalence we identified as directly associated with underlying insulin resistance of the background strain - generated a full-body knockout. Full-body EP3 loss provided no diabetes protection to BTBR-Ob mice, but, unexpectedly, significantly worsened BTBR-lean insulin resistance and glucose tolerance. This in vivo phenotype was not associated with changes in beta-cell fractional area or markers of beta-cell replication ex vivo. Instead, EP3-null BTBR-lean islets had essentially uncontrolled insulin hypersecretion. The selective up-regulation of constitutively-active EP3 splice variants in islets from young, lean BTBR mice as compared to C57BL/6J, where no phenotype of EP3 loss has been observed, provides a potential explanation for the hypersecretion phenotype. In support of this, high islet EP3 expression in Balb/c females vs. Balb/c males was fully consistent with their sexually-dimorphic metabolic phenotype after loss of EP3-coupled Gz protein. Taken together, our findings provide a new dimension to the understanding of EP3 as a critical brake on insulin secretion.

New and NoteworthyIslet Prostaglandin EP3 receptor (EP3) signaling is well-known as up-regulated in the pathophysiological conditions of type 2 diabetes, contributing to beta-cell dysfunction. Unexpected findings in mouse models of non-obese insulin sensitivity and resistance provide a new dimension to our understanding of EP3 as a key modulator of insulin secretion. A previously-unknown relationship between mouse insulin resistance and the penetrance of Rat insulin promoter-driven germline floxed allele recombination is critical to consider when creating beta-cell-specific knockouts.

For Table of Contents Use Only

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=121 SRC="FIGDIR/small/671289v3_ufig1.gif" ALT="Figure 1">
View larger version (37K):
org.highwire.dtl.DTLVardef@15cf1d5org.highwire.dtl.DTLVardef@108f565org.highwire.dtl.DTLVardef@126e479org.highwire.dtl.DTLVardef@61fdf8_HPS_FORMAT_FIGEXP  M_FIG C_FIG Targeted Gene Editing in Plants using CRISPR-Cas9, Single-Stranded DNA Oligonucleotides, and Protoplast Regeneration Versatile genome editing can be facilitated by the insertion of DNA sequences into specific locations. Current protocols involving clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) proteins rely on low efficiency homology-directed repair or non-homologous end joining with modified double-stranded DNA oligonucleotides as donors. Our simple protocol eliminates the need for expensive equipment, chemical and enzymatic donor DNA modification or plasmid construction by using polyethylene glycol-calcium to deliver non-modified single-stranded DNA oligonucleotides and CRISPR-Cas9 ribonucleoprotein into protoplasts. Plants regenerated via edited protoplasts achieved targeted insertion frequencies of up to 50.0% in Nicotiana benthamiana and 13.6% in rapid cycling Brassica oleracea without antibiotic selection. Using a 60-nt donor containing 27 nt in each homologous arm, 6 of 22 regenerated N. benthamiana plants showed targeted insertions, and one contained a precise insertion of a 6-bp HindIII site. The inserted sequences were transmitted to the next generation and invite the possibility of future exploration of versatile genome editing by targeted DNA insertion in plants. β-glucan imprinting remodels macrophage function in response to environmental cues In vitro, exposure of human primary monocytes to the fungal {beta}-glucan enhances their pro-inflammatory responsiveness towards several pathogens. Yet, the role of environmental condition of this process remains unclear. Here we found that {beta}-glucan-induced innate immune memory counteract the anti-inflammatory status of the macrophages. In response to {beta}-glucan imprinting, M-CSF-(M2-like-) macrophages increase their pro-inflammatory responsiveness to secondary stimuli associated with decrease of the M-CSF differentiation hallmarks. In contrast, in GM-CSF-(M1-like-) environment, {beta}-glucan imprinting reduced the pro-inflammatory canonical feature of the macrophages, together with their terminal differentiation marks. Comparing M-CSF and GM-CSF environment, we observed that {beta}-glucan-imprinted macrophages present comparable functions in terms of cytokine responses, phagocytosis, oxidative burst, and angiogenesis. This effect is mediated through Dectin-1 and associated with altered expression of the master regulators of macrophage terminal differentiation, IRF5 and IRF3. {beta}-glucan-induced innate immune memory skews the commitment of the macrophages in complex environment towards similar and less terminally differentiated cells. Together, these observations suggest a potential therapeutic role for {beta}-glucan-induced modulation of innate memory in different pathological contexts. First detection of herpesvirus and frequency of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Substrate specialization in microbes is driven by biochemical constraints of dynamic flux sensing Central carbon metabolism is highly conserved across microbial species, but can catalyze very different pathways depending on the organism and their ecological niche. Here, we study the dynamic re-organization of central metabolism after switches between the two major opposing pathway configurations of central carbon metabolism, glycolysis and gluconeogenesis in Escherichia coli, Pseudomonas aeruginosa and Pseudomonas putida. We combined growth dynamics and dynamic changes of intracellular metabolite levels with a coarse-grained model that integrates fluxes, regulation, protein synthesis and growth and uncovered fundamental limitations of the regulatory network: after nutrient shifts, metabolite concentrations collapse to their equilibrium, rendering the cell unable to sense which direction the flux is supposed to flow through the metabolic network. The cell can partially alleviate this by picking a preferred direction of regulation at the expense of increasing lag times in the opposite direction. Moreover, decreasing both lag times simultaneously comes at the cost of reduced growth rate or higher futile cycling between metabolic enzymes. These three trade-offs can explain why microorganisms specialize for either glycolytic or gluconeogenic substrates and can help elucidate the complex growth patterns exhibited by different microbial species.

Graphical synopsis O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

Standfirst textMicrobes face a series of fundamental trade-offs that limit their ability to optimize simultaneously for both glycolytic and gluconeogenic growth.

Bullet pointsO_LILag times between glycolysis and gluconeogenesis show asymmetry in many microbes: A long lag in one direction, but a short lag in the other.
C_LIO_LILong lag times are caused by an inability to sense fluxes after nutrient shifts.
C_LIO_LIWith existing regulation, lag time asymmetry can only be overcome by reducing either growth rate or increasing futile cycling in metabolism.
C_LI Deficits in seizure threshold and other behaviors in adult mice without gross neuroanatomic injury after late gestation transient prenatal hypoxia Intrauterine hypoxia is a common cause of brain injury in children resulting in a broad spectrum of long-term neurodevelopmental sequela, including life-long disabilities that can occur even in the absence of severe neuroanatomic damage. Postnatal hypoxia-ischemia rodent models are commonly used to understand the effects of ischemia and transient hypoxia on the developing brain. Postnatal models, however, have some limitations. First, they do not test the impact of placental pathologies on outcomes from hypoxia. Second, they primarily recapitulate severe injury because they provoke substantial cell death, which is not seen in children with mild hypoxic injury. Lastly, they do not model preterm hypoxic injury. Prenatal models of hypoxia in mice may allow us to address some of these limitations to expand our understanding of developmental brain injury. The published rodent models of prenatal hypoxia employ multiple days of hypoxic exposure or complicated surgical procedures, making these models challenging to perform consistently in mice. Furthermore, large animal models suggest that transient prenatal hypoxia without ischemia is sufficient to lead to significant functional impairment to the developing brain. However, these large animal studies are resource-intensive and not readily amenable to mechanistic molecular studies. Therefore, here we characterized the effect of late gestation (embryonic day 17.5) transient prenatal hypoxia (5% inspired oxygen) on long-term anatomical and neurodevelopmental outcomes in mice. Late gestation transient prenatal hypoxia increased hypoxia-inducible factor 1 alpha protein levels (a marker of hypoxic exposure) in the fetal brain. Hypoxia exposure predisposed animals to decreased weight at postnatal day 2, which normalized by day 8. However, hypoxia did not affect gestational age at birth, litter size at birth, or pup survival. No differences in fetal brain cell death or long-term gray or white matter changes resulted from hypoxia. Animals exposed to prenatal hypoxia did have several long-term functional consequences, including sex-dichotomous changes. Hypoxia exposure was associated with a decreased seizure threshold and abnormalities in hindlimb strength and repetitive behaviors in males and females. Males exposed to hypoxia had increased anxiety-related deficits, whereas females had deficits in social interaction. Neither sex developed any motor or visual learning deficits. This study demonstrates that late gestation transient prenatal hypoxia in mice is a simple, clinically relevant paradigm for studying putative environmental and genetic modulators of the long-term effects of hypoxia on the developing brain. A General Biphasic-Bodyweight Model for Scaling Basal Metabolic Rate, Glomerular Filtration Rate and Drug Clearance from Birth to Adulthood Understanding the maturation process of human physiology and metabolism has broad medical and pharmaceutical implications. Age and bodyweight are frequently considered as separate variables in modeling the dynamical changes of human organ functions and of drug clearance from birth to adulthood. The objective of this study is to propose a unified, continuous and bodyweight-only equation to quantify the changes of human basal metabolic rate (BMR), glomerular filtration rate (GFR) and drug clearance (CL) from infancy to adulthood. The BMR datasets were retrieved from a comprehensive historical database of male and female subjects (0.02 to 64 years). The CL datasets for 17 drugs and the GFR dataset were generated from age-incorporated maturation-and-growth models with reported parameter values. The model used in the simulation is independent of the proposed model. A statistical approach was used to simulate the model generated CL and GFR data for a hypothetical population with 26 age groups (ranging from 0 to 20 years). Besides, individual CL data for one drug, and sparse PBPK-modeled CL values for two drugs were also included for analysis. A 4-parameter, mixed-allometry equation with two power-law functions of bodyweight was proposed and evaluated as a general model using nonlinear regression and dimensionless analysis. All datasets universally reveal biphasic curves with two distinct linear segments on log-log plots. Compared with simple allometry, the biphasic model fits satisfactorily to all datasets (based on Akaikes Information Criterion and residual plots). The biphasic equation consists of two reciprocal allometric terms that asymptotically determine the overall curvature. The fitting results show a superlinear scaling phase (slope >1; ca. 1.5 - 3.5) below the characteristic bodyweight at the phase transition; and above which, a sublinear scaling phase (slope <1; ca. 0.5 - 0.7) is evident. The phase-transition bodyweight is ranging from 5 to 20 kg (corresponding to 0.5 - 9 years) and the mean value is around 10 kg ([~]2 years) for all data sets. The dimensionless analysis generalizes, and offers quantitative realization of, the maturation and growth process. In conclusion, the proposed mixed-allometry equation is a generic model that quantitatively describes the phase transition occurring in the human maturation process of BMR, GFR and drug CL. Flexibility of Short-chain dehydrogenase is interconnected to its promiscuity for the reduction of multiple ketone intermediates Short-chain dehydrogenases/reductases (SDRs) are a convenient class of enzymes used to synthesize enantiopure alcohols. Several studies describe native or engineered SDRs for converting substrates of interest using cost and time-intensive high-throughput approaches. The classification of SDRs is based on chain length and cofactor binding site. Of these, the shorter  Classical and the longer  Extended enzymes participate in ketoreduction. However, comparative analysis of various modelled SDRs reveals a length independent conserved N-terminal Rossmann fold and a variable C-terminus region for both types. The general hypothesis is that the latter domain influences the enzymes flexibility that may affect the observed promiscuity of the enzyme. We have used a machine learning algorithm on this flexible domain to build a rationale to screen promiscuous candidates. We have built a data set consisting of physicochemical properties derived from the amino-acid composition of enzymes to select closely associated promiscuous mesophilic enzymes. The resulting in vitro studies on pro-pharmaceutical substrates illustrate a direct correlation between the C-terminal lid-loop structure, enzyme melting temperature and the turnover number. We present a walkthrough for exploring promiscuous SDRs for catalyzing enantiopure alcohols of industrial importance. 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. ALLOPURINOL BLOCKS THE AORTIC ANEURYSM IN A MOUSE MODEL OF MARFAN SYNDROME VIA REDUCING AORTIC OXIDATIVE STRESS BackgroundThe pathogenesis and progression of aortic aneurysm in Marfan syndrome (MFS) involves dysregulated TGF-{beta} and nitric oxide signaling, altered hemodynamics, and biomechanical forces. Increasing evidence indicates that redox stress participates in MFS aortopathy, though its contribution is not well established. We reported elevated reactive oxygen species (ROS) formation and NADPH oxidase NOX4 upregulation in MFS patients and mice aortic samples. Here we address the contribution of xanthine oxidoreductase (XOR), which is transcribed as dehydrogenase (XDH) and then post-translationally modified to oxidase (XO), catabolizing mainly purines into uric acid and ROS.

Methods and ResultsWe report that in MFS patients, XOR protein levels increased in aortic samples. In MFS mice (Fbn1C1041G/+), XOR mRNA transcripts augmented in the aorta of young but not older MFS mice. Moreover, the enzymatic activity of XO form increased with respect to XDH. The administration of the XOR inhibitor allopurinol (ALO) blocked the progression of aortic root aneurysm in MFS mice. ALO was also protective when administrated before the appearance of aneurysm. MFS-induced elastic fiber fragmentation, fibrotic remodeling, nuclear translocation of pNRF2, and increased 3-nitrotyrosine levels in the aortic tunica media as well as endothelial dysfunction were all prevented by allopurinol. The large overproduction of H2O2 in MFS aorta did not occur in allopurinol-treated mice whose uric acid plasma levels remained unaltered.

ConclusionsThis study strengthens the concept that redox stress is an important determinant of aortic aneurysm formation and progression in MFS and supports a clinical trial for ALO in the pharmacological treatment of MFS aortopathy.

GRAPHICAL ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=146 SRC="FIGDIR/small/464182v3_ufig1.gif" ALT="Figure 1">
View larger version (31K):
org.highwire.dtl.DTLVardef@eab02aorg.highwire.dtl.DTLVardef@1d3b39eorg.highwire.dtl.DTLVardef@1b54d0forg.highwire.dtl.DTLVardef@1d61ad0_HPS_FORMAT_FIGEXP  M_FIG C_FIG Cryptic pathogen-sugar interactions revealed by universal saturation transfer analysis Many host pathogen interactions such as human viruses (including non-SARS-coronaviruses) rely on attachment to host cell-surface glycans. There are conflicting reports about whether the Spike protein of SARS-CoV-2 binds to sialic acid commonly found on host cell-surface N-linked glycans. In the absence of a biochemical assay, the ability to analyze the binding of glycans to heavily- modified proteins and resolve this issue is limited. Classical Saturation Transfer Difference (STD) NMR can be confounded by overlapping sugar resonances that compound with known experimental constraints. Here we present  universal saturation transfer analysis (uSTA), an NMR method that builds on existing approaches to provide a general and automated workflow for studying protein-ligand interactions. uSTA reveals that B-origin-lineage-SARS-CoV-2 spike trimer binds sialoside sugars in an  end on manner and modelling guided by uSTA localises binding to the spike N-terminal domain (NTD). The sialylated-polylactosamine motif is found on tetraantennary human N-linked-glycoproteins in deeper lung and may have played a role in zoonosis. Provocatively, sialic acid binding is abolished by mutations in some subsequent SARS- CoV-2 variants-of-concern. A very high resolution cryo-EM structure confirms the NTD location and  end on mode; it rationalises the effect of NTD mutations and the structure-activity relationship of sialic acid analogues. uSTA is demonstrated to be a robust, rapid and quantitative tool for analysis of binding, even in the most demanding systems.

Extended AbstractThe surface proteins found on both pathogens and host cells mediate entry (and exit) and influence disease progression and transmission. Both types can bear host-generated post- translational modifications such as glycosylation that are essential for function but can confound biophysical methods used for dissecting key interactions. Several human viruses (including non- SARS-coronaviruses) attach to host cell-surface N-linked glycans that include forms of sialic acid (sialosides). There remains, however, conflicting evidence as to if or how SARS-associated coronaviruses might use such a mechanism. Here, we demonstrate quantitative extension of  saturation transfer protein NMR methods to a complete mathematical model of the magnetization transfer caused by interactions between protein and ligand. The method couples objective resonance-identification via a deconvolution algorithm with Bloch-McConnell analysis to enable a structural, kinetic and thermodynamic analysis of ligand binding beyond previously-perceived limits of exchange rates, concentration or system. Using an automated and openly available workflow this  universal saturation transfer analysis (uSTA) can be readily-applied in a range of even heavily-modified systems in a general manner to now obtain quantitative binding interaction parameters (KD, kEx). uSTA proved critical in mapping direct interactions between natural sialoside sugar ligands and relevant virus-surface attachment glycoproteins - SARS-CoV-2-spike and influenza-H1N1-haemagglutinin variants - by quantitating ligand signal in spectral regions otherwise occluded by resonances from mobile protein glycans (that also include sialosides). In B- origin-lineage-SARS-CoV-2 spike trimer  end on-binding to sialoside sugars was revealed contrasting with  extended surface-binding for heparin sugar ligands; uSTA-derived constraints used in structural modelling suggested sialoside-glycan binding sites in a beta-sheet-rich region of spike N-terminal domain (NTD). Consistent with this, uSTA-glycan binding was minimally- perturbed by antibodies that neutralize the ACE2-binding domain (RBD) but strongly disrupted in spike from the B1.1.7/alpha and B1.351/beta variants-of-concern, which possess hotspot mutations in the NTD. Sialoside binding in B-origin-lineage-NTD was unequivocally pinpointed by cryo-EM to a site that is created from residues that are notably deleted in variants (e.g. H69,V70,Y145 in alpha). An analysis of beneficial genetic variances in cohorts of patients from early 2020 suggests a model in which this site in the NTD of B-origin-lineage-SARS-CoV-2 (but not in alpha/beta-variants) may have exploited a specific sialylated-polylactosamine motif found on tetraantennary human N-linked-glycoproteins in deeper lung. Together these confirm a novel binding mode mediated by the unusual NTD of SARS-CoV-2 and suggest how it may drive virulence and/or zoonosis via modulation of glycan attachment. Since cell-surface glycans are widely relevant to biology and pathology, uSTA can now provide ready, quantitative, widespread analysis of complex, host-derived and post-translationally modified proteins with putative ligands relevant to disease even in previously confounding complex systems. Kinetoplastid-specific X2-family kinesins interact with a kinesin-like pleckstrin homology domain protein that localizes to the trypanosomal microtubule quartet Kinesins are motor proteins found in all eukaryotic lineages that move along microtubules to mediate cellular processes such as mitosis and intracellular transport. In trypanosomatids, the kinesin superfamily has undergone a prominent expansion, resulting in one of the most diverse kinesin repertoires that includes the two kinetoplastid-restricted families X1 and X2. Here, we characterize in Trypanosoma brucei TbKifX2A, an orphaned X2 kinesin. TbKifX2A tightly interacts with TbPH1, a kinesin-like protein with a likely inactive motor domain, a rarely reported occurrence. Both TbKifX2A and TbPH1 localize to the microtubule quartet (MtQ), a characteristic but poorly understood cytoskeletal structure that wraps around the flagellar pocket as it extends to the cell body anterior. The proximal proteome of TbPH1 revealed two other interacting proteins, the flagellar pocket protein FP45 and intriguingly another X2 kinesin, TbKifX2C. Simultaneous ablation of TbKifX2A/TbPH1 results in the depletion of FP45 and TbKifX2C and also an expansion of the flagellar pocket, among other morphological defects. TbKifX2A is the first motor protein to be localized to the MtQ. The observation that TbKifX2C also associates with the MtQ suggests that the X2 kinesin family may have co-evolved with the MtQ, both kinetoplastid-specific traits. Reduced glutathione decreases cell adhesion and increases cell volume. Glutathione is the most abundant thiol in animal cells. Reduced glutathione (GSH) is a major intracellular antioxidant neutralizing free radicals and detoxifying electrophiles. It plays important roles in many cellular processes, including cell differentiation, proliferation, and apoptosis. In the present study we demonstrate that extracellular concentration of reduced glutathione markedly increases cell volume within few hours, in a dose-response manner. Pre-incubation of cells with BSO, the inhibitor of {gamma}-glutamylcysteine synthetase, responsible for the first step in intracellular glutathione synthesis did not change the effect of reduced glutathione on cell volume suggesting a mechanism limited to the interaction of extracellular reduced glutathione on cell membrane. Similarly, inhibition of {gamma}-glutamylcyclotransferase involved in intracellular glutamate production had no effect on the action of reduced glutathione. Oxidized glutathione exerted no effect on cell volume. Results show that reduced GSH decreases cell adhesion resulting in an increased cell volume. Since many cell types are able to export GSH, the present results suggest that this could be a fundamental self-regulation of cell volume, giving the cells a self-control on their adhesion proteins. Recovery of high-quality assembled genomes via metagenome binning guided with single-cell amplified genomes BackgroundObtaining high-quality (HQ) reference genomes from microbial communities is crucial for understanding the phylogeny and function of uncultured microbes in complex microbial ecosystems. Despite the improved bioinformatic approaches to generate curated metagenome-assembled genomes (MAGs), existing metagenomic binners often fail to obtain reliable MAGs, and thus, they are nowhere comparable to genomes sequenced from isolates in terms of strain level resolution. Here, we present a single-cell genome-guided metagenome binning (MetaSAG) to reconstruct the strain-resolved genomes from microbial communities at once.

ResultsMetaSAG employs single-cell amplified genomes (SAGs) generated with microfluidic technology as binning guides to recover improved draft genomes with the metagenomic data. To assess the performance of reconstructing genomes from various microbial communities, we compared MetaSAG with four conventional metagenomic binners using a cell mock community, human gut microbiota, and skin microbiota samples. MetaSAG showed precise contig binning and higher recovery rates (>97%) of rRNA and plasmids compared to conventional binners in genome reconstruction from the cell mock community. In human microbiota samples, MetaSAG recovered the largest number of genomes with a total of 103 gut microbial genomes (21 HQ and 65 showed >90% completeness) and 45 skin microbial genomes (10 HQ and 40 showed >90% completeness), respectively. Conventional binners recovered one Staphylococcus hominis genome, whereas MetaSAG recovered two S. hominis genomes from the identical skin microbiota sample. Single-cell sequencing indicated that these S. hominis genomes clearly derived from two distinct strains harboring specifically different plasmids. We found that all conventional S. hominis MAGs had substantial lack or excess of the genome sequences and contamination of other Staphylococcus bacteria (S. epidermidis).

ConclusionsMetaSAG enabled us to obtain the strain-resolved genomes in the mock community and human microbiota samples by assigning metagenomic sequences correctly and covering both highly conserved genes such as rRNA genes and unique extrachromosomal elements, including plasmids. MetaSAG will provide HQ genomes that are difficult to obtain with metagenomic analyses alone and will facilitate the understanding of microbial ecosystems by elucidating detailed metabolic pathways and horizontal gene transfer networks through microbial genomes. MetaSAG is available at https://github.com/kojiari/metasag. Decoding nonspecific interactions between human nuclear transport proteins: A computational study The nuclear protein transport between the nucleus and cytosol can be considered a core process of cell regulation. Specially designed proteins in nature such as importins, exportins, and some other transporters facilitate this transport in the cell and control the cellular processes. Transient and weak protein-protein interactions are basis of these various biomolecular processes. Prior to cargo transports, the transport proteins recognize the Nuclear localization signals (NLSs) and Nuclear export signals (NESs) of cargo proteins and, bind to the RanGTP. Also, these proteins bind with other similar protein subunits along with RanGTP to transport cargos. Cell is enormously crowded place where DNA, RNA, proteins, lipids and small molecules cooperatively facilitate numerous cellular processes. In such environment, existence of nonspecific interactions between proteins is quite obvious. Considering this hypothesis, in this study, protein-protein docking approach was applied to determine the binding affinities of 12 human nuclear transport proteins. Results showed that KPNA1, TNPO1 and TNPO3 have greater affinity to bind with other transport proteins. Also, among 78 complexes (12 homodimers and 66 heterodimers), KPNA1-KPNB1, KPNA1-TNPO1 and KPNA1-TNPO3 complexes have the highest stability.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=156 SRC="FIGDIR/small/436462v2_ufig1.gif" ALT="Figure 1">
View larger version (40K):
org.highwire.dtl.DTLVardef@d81d48org.highwire.dtl.DTLVardef@6a929dorg.highwire.dtl.DTLVardef@bfa16eorg.highwire.dtl.DTLVardef@ff610b_HPS_FORMAT_FIGEXP  M_FIG C_FIG Initially, 12 human nuclear transport proteins PDB structures were retrieved from the 1. Protein data bank (PDB). These proteins had some missing terminals and residues thus, we used 2. SWISS-MODEL and 3. MODELLER v.10.1 to model those regions in these proteins. Next, we used widely popular web server, 4. ClusPro v.2.0 for protein-protein docking analysis among 12 proteins. Then, we employed 5. PRODIGY web server to calculate the binding affinities of 78 complexes (12 homodimers & 66 heterodimers). Finally, we utilised three web tools, 6. Arpeggio, 7. PIMA and 8. PDBePISA to analyse top-three complexes (KPNA1-KPNB1, KPNA1-TNPO1 & TNPO3) for in-depth interactions and energetics. The obstacles and potential clues of prime editing applications in tomato, a dicot plant Precision genome editing is highly desired for crop improvement. The recently emerged CRISPR/Cas technology offers great potential applications in precision plant genome engineering. A prime editing (PE) approach combining a reverse transcriptase (RT) with a Cas9 nickase and a "priming" extended guide RNA has shown a high frequency for precise genome modification in mammalian cells and several plant species. Nevertheless, the applications of the PE approach in dicot plants are still limited and inefficient. We designed and tested prime editors for precision editing of a synthetic sequence in a transient assay and for desirable alleles of 10 loci in tomato by stable transformation. Our data obtained by targeted deep sequencing also revealed only low PE efficiencies in both the tobacco and tomato systems. Further assessment of the activities of the PE components uncovered potential reasons for the inefficiency of the PE complexes. Modifying the pegRNA sequences by shortening or introducing mismatches to the primer binding sequences (PBS) in order to reduce their melting temperatures (Tm) did not enhance the PE efficiency at the SlBMP21, SlALC and SlALS1 loci. Our data show challenges of PE approach in tomato, indicating a further improvement of the PE system for the successful applications such as use of improved expression systems. Our work provides an important clue for the successful application of the PE approach in crop improvement. Dynamic gain decomposition reveals functional effects of dendrites, ion channels and input statistics in population coding The information processing capabilities of large neuronal circuits form the basis of higher-level cognition. It is unknown how the constituent cells properties interact with the external stimuli to shape the networks collective activity patterns and the associated computations. The dynamic gain function, a spectrally resolved approximation of a populations coding capability, is sensitive to the properties of stimuli and cells and could provide a mechanistic link from cell-level to network-level performance. However, the dynamic gains shape is co-determined in a complex manner by axonal and somatodendritic parameters and the populations operating regime, which so far precluded an understanding of any individual parameters impact. Here, we decomposed the dynamic gain function into three components corresponding to separate signal transformations. This allowed attribution of network-level encoding features to specific cell-level parameters. Applying the method to data from real neurons and biophysically plausible models, we found: 1. The encoding bandwidth of real neurons, approximately 400 Hz, is constrained by the voltage dependence of axonal currents during an early stage of action potential initiation. 2. The state-of-the-art models only achieve encoding bandwidths around 100 Hz and are mostly limited by subthreshold processes instead. 3. Large dendrites and low-threshold potassium currents modulate the bandwidth by shaping the subthreshold stimulus-to-voltage transformation. Dynamic gain decomposition can physiologically interpret changes in the dynamic gain curve, which is particularly desirable in studies of cellular effects in spectrinopathies and neurodegeneration. By pinpointing shortcomings of current models, it also guides inference of neuron models best suited for large-scale network simulations.

ImpactThe dynamic gain function quantifies how neurons can engage in collective, network-level activity, shaping brain rhythms and corresponding information encoding capabilities. Its shape results from a complex interaction between properties of different molecules (ion channels) and cell compartments (morphology, resistance), and is so far only understood for the simplest neuron models. Here we provide an interpretable analysis, decomposing dynamic gain based on the stimulus transformation steps in individual neurons. It is applicable to models of any complexity and real neurons, and allows attribution of dynamic gain changes to specific sub- and suprathreshold processes. Using this decomposition method, we reveal the relevance of subthreshold potassium channels for ultrafast information encoding and expose the shortcomings of even the state-of-the-art neuron models. Naming Human Diseases: Ethical Principles of Curating Exclusive Substitute for Inopportune Nosology BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. Teosinte mexicana introgression modulates phosphatidylcholine levels and maize flowering time Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp.parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability, and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern US, Canada and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time.