Delayed but generalist resistance in alternation: a trade-off shaped by the mode of action of antifungals, as revealed by experimental evolution in a plant pathogen. The evolution of resistance to pesticides is a major burden in agriculture. Resistance management involves maximizing selection pressure heterogeneity, particularly by combining active ingredients with different modes of action. We tested the hypothesis that the temporal alternation of active ingredients may delay the build-up of resistance not only by spreading selection pressure over longer periods, but also by decreasing the rate of evolution of resistance to alternated fungicides. Here, we applied an original experimental evolution approach to the economically important crop pathogen Zymoseptoria tritici. We observed the dynamics of Z. tritici resistance in 56 independent lines subjected to 14 continuous or alternation regimes of 3 fungicides contrasting for their mode of action, at their EC95 selection doses. For the first time in a plant pathogen, our results show that alternation is either neutral or slows the evolution of resistance, relative to continuous fungicide use, but results in higher levels of generalism in evolved lines. We demonstrate that the mode of action of resistance of fungicides drivingly underlies this trade-off, more so than the number of fungicides and the frequency of alternation. This trade-off is also dynamic over the course of resistance evolution, as shown by the change in population phenotype structure and the relative impact of selection drivers. These findings open up new possibilities for tailoring resistance management effectively while optimizing smart interplay between alternation components. They also confirm experimental evolution as an untapped but promising approach to dissect adaptation in phytopathogenic fungi.

Author summaryThe efficacy of pesticides has been compromised by the generalization of their use, leading to the rapid and widespread evolution of resistance. This constitutes a major economic and environmental burden in agriculture. The temporal alternation (or cycling) of active ingredients is a management strategy that induces temporal variation of selection of pathogens. Here, we dissected how it can be optimized according to the number and nature of alternated modes of action and to the rhythm of their application and sought to understand how these drivers determine the performance of alternation strategies. We used an approach original in plant pathology, experimental evolution, applied on an economically important fungus, Zymoseptoria tritici. We concluded that alternation can delay the rate of resistance selection with performance depending most likely on the mode of action of alternated antifungals, more so than on other drivers. But we also highlighted that alternation regimes select generalist resistance, that is the ability to resist a large diversity of antifungals. This trade-off was dynamic over time. Our study provides new insights for the informed management of pesticides and the reduction of their side effects.
Discovery of archaeal Fusexins homologous to eukaryotic HAP2/GCS1 gamete fusion proteins Sexual reproduction consists of genome reduction by meiosis and subsequent gamete fusion. Presence of meiotic genes in prokaryotes suggests that DNA repair mechanisms evolved toward meiotic recombination; however, fusogenic proteins resembling those found in eukaryotes were not identified in prokaryotes. Here, we identify archaeal proteins that are homologs of fusexins, a superfamily of fusogens that mediate eukaryotic gamete and somatic cell fusion, as well as virus entry. The crystal structure of a trimeric archaeal Fusexin1 reveals novel features such as a six-helix bundle and an additional globular domain. Ectopically expressed Fusexin1 can fuse mammalian cells, and this process involves the additional domain and a conserved fusion loop. Archaeal fusexin genes exist within integrated mobile elements, potentially linking ancient archaeal gene exchanges and eukaryotic sex.

One-Sentence SummaryCell membrane fusion proteins of viruses and eukaryotes are also present in archaea.
Existence of log-phase Escherichia coli persisters and lasting memory of a starvation pulse The vast majority of a bacterial population is quickly killed when treated with a lethal concentration of antibiotics. The time scale of this killing is often comparable with the bacterial generation time before addition of antibiotics. Yet, a small subpopulation typically survives for an extended period. However, the long-term killing dynamics of bacterial cells has not been fully quantified even in well-controlled laboratory conditions. We constructed a week-long killing assay and followed the survival fraction of Escherichia coli K12 exposed to a high concentration of ciprofloxacin. We found that long-term survivors were formed during exponential growth, with some cells surviving at least 7 days. The long-term dynamics contained at least three timescales, which greatly enhances predictions of the population survival time compared to the biphasic extrapolation from the short term behavior. Furthermore, we observed a surprisingly long memory effect of a brief carbon starvation pulse, which was dependent on the (p)ppGpp synthase relA. Specifically, one hour of carbon starvation prior to antibiotics exposure increased the surviving fraction by nearly 100-fold even after 4 days of ciprofloxacin treatment. Sex determination gene transformer regulates the male-female difference in Drosophila fat storage via the adipokinetic hormone pathway Sex differences in whole-body fat storage exist in many species. For example, Drosophila females store more fat than males. Yet, the mechanisms underlying this sex difference in fat storage remain incompletely understood. Here, we identify a key role for sex determination gene transformer (tra) in regulating the male-female difference in fat storage. Normally, a functional Tra protein is present only in females, where it promotes female sexual development. We show that loss of Tra in females reduced whole-body fat storage, whereas gain of Tra in males augmented fat storage. Tras role in promoting fat storage was largely due to its function in neurons, specifically the Adipokinetic hormone (Akh)-producing cells (APCs). Our analysis of Akh pathway regulation revealed a male bias in APC activity and Akh pathway function, where this sex-biased regulation influenced the sex difference in fat storage by limiting triglyceride accumulation in males. Importantly, Tra loss in females increased Akh pathway activity, and genetically manipulating the Akh pathway rescued Tra-dependent effects on fat storage. This identifies sex-specific regulation of Akh as one mechanism underlying the male-female difference in whole-body triglyceride levels, and provides important insight into the conserved mechanisms underlying sexual dimorphism in whole-body fat storage. SMC protein RecN drives RecA filament translocation and remodelling for in vivo homology search While the molecular repertoire of the homologous recombination pathway is well-studied, the search mechanism that enables recombination between distant homologous regions is poorly understood. Earlier work suggests that the recombinase RecA, an essential component for homology search, forms an elongated filament, nucleating at the break site. How this RecA structure carries out long distance search remains unclear. Here, we follow the dynamics of RecA after induction of a single double-strand break on the Caulobacter chromosome. We find that the RecA-nucleoprotein filament, once formed, rapidly translocates in a directional manner in the cell, undergoing several pole-to-pole traversals, until homology search is complete. Concomitant with translocation, we observe dynamic remodelling of the filament. Importantly in vivo, the RecA filament alone is incapable of such long distance movement; both translocation and dynamic remodelling are contingent on action of SMC-like protein RecN, via its ATPase cycle. We provide a stochastic description of RecN-regulated changes in filament length during translocation via modulation of RecA assembly-disassembly. In summary, we have uncovered the three key elements of homology search driven by RecN: mobility of a finite segment of RecA, filament remodelling and ability to conduct multiple pole-to-pole traversals, which together point to a novel optimal search strategy. Creatine transporter deficient rat model show motor dysfunction, cerebellar alterations and muscle creatine deficiency without muscle atrophy Creatine (Cr) is a nitrogenous organic acid and plays roles as fast phosphate energy buffer to replenish ATP, osmolyte, antioxidant, neuromodulator, and as a compound with anabolic and ergogenic properties in muscle. Cr is taken from the diet or endogenously synthetized by the enzymes AGAT and GAMT, and specifically taken up by the transporter SLC6A8. Loss-of-function mutations in the genes encoding for the enzymes or the transporter cause Cerebral Creatine Deficiency Syndromes (CCDS). CCDS are characterized by brain Cr deficiency, intellectual disability with severe speech delay, behavioral troubles, epilepsy and motor dysfunction. Among CCDS, the X-linked Cr transporter deficiency (CTD) is the most prevalent with no efficient treatment so far. Different animal models of CTD show reduced brain Cr levels, cognitive deficiencies and together they cover other traits similar to those of patients. However, motor function was poorly explored in CTD models and some controversies in the phenotype exist in comparison with CTD patients. Our recently described Slc6a8Y389C knock-in (KI) rat model of CTD showed mild impaired motor function, morphological alterations in cerebellum, reduced muscular mass, Cr deficiency and increased guanidinoacetate content in muscle, although no consistent signs of muscle atrophy. Our results indicate that such motor dysfunction co-occurred with both nervous and muscle dysfunction, suggesting that muscle strength and performance as well as neuronal connectivity might be affected by this Cr deficiency in muscle and brain. ZAF -- The First Open Source Fully Automated Feeder for Aquatic Facilities In the past few decades, aquatic animals have become popular model organisms in biology, spurring a growing need for establishing aquatic facilities. Zebrafish are widely studied and relatively easy to culture using commercial systems. However, a challenging aspect of maintaining aquatic facilities is animal feeding, which is both time- and resource-consuming. We have developed an open-source fully automatic daily feeding system, Zebrafish Automatic Feeder (ZAF). ZAF is reliable, provides a standardized amount of food to every tank, is cost-efficient and easy to build. The advanced version, ZAF+, allows for the precise control of food distribution as a function of fish density per tank, and has a user-friendly interface. Both ZAF and ZAF+ are adaptable to any laboratory environment and facilitate the implementation of aquatic colonies. Here we provide all blueprints and instructions for building the mechanics, electronics, fluidics, as well as to setup the control software and its user-friendly graphical interface. Importantly, the design is modular and can be scaled to meet different user needs. Furthermore, our results show that ZAF and ZAF+ do not adversely affect zebrafish culture, enabling fully automatic feeding for any aquatic facility. An integrative analysis to predict the Active Ingredients and explore Polypharmacological Mechanisms of Orthosiphon stamineus Benth. for human health BackgroundOrthosiphon stamineus Benth. was a dietary supplement or a traditional folk herb with widespread clinical application, but it was lack of comprehensive understanding of its active ingredients and polypharmacological mechanisms. The aim of this work was to systematically study its natural compounds and the molecular mechanisms of the O. stamineus via network pharmacology.

MethodsCompounds from O. stamineus were collected by literature retrieval and evaluated by SwissADME with the physicochemical property (ADMET model) and the likelihood for a natural medicine. The connection of active ingredients and target genes was built and confirmed by Cytoscape and AutoDock vina. Then, Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were separately performed to obtain a more in-depth understanding of O. stamineus. Finally, the relationship among active ingredients, targets, and diseases was built to clarify the polypharmacological mechanisms and found relevant active substances for further drug discovery.

ResultsA total of 159 compounds of O. stamineus were collected, and 22 potentially active ingredients were screened out. The Ingredient-Target Interaction Network was built with 12 flavonoids, 3 diterpenes, 3 phenols, and 4 volatile oils, and 65 targets. The Docking analysis indicated that the ingredient-target interaction network was reliable; most ligand-receptor had a strong binding affinity (lowest binding energy: -6.9 kcal/mol). After pathway analysis, 185 significant biological processes and 36 signal pathways were found, and the ingredient-target-disease network of O. stamineus was constructed for polypharmacological mechanisms.

ConclusionOur study clarified the polypharmacological mechanisms via the relationship among active ingredients, targets, and diseases and provided better guidance for subsequent experiments and potential active ingredients for drug discovery or health promotion. An alternatively spliced TREM2 isoform lacking the ligand binding domain is expressed in human brain BackgroundGenetic variants in TREM2 are strongly associated with Alzheimers Disease (AD) risk but alternative splicing in TREM2 transcripts has not been comprehensively described.

ObjectiveRecognizing that alternative splice variants can result in reduced gene expression and/or altered function, we sought to fully characterize splice variation in TREM2.

MethodsHuman blood and anterior cingulate autopsy tissue from 61 donors were used for end-point and quantitative PCR and Western blotting to identify and quantify novel TREM2 isoforms.

ResultsIn addition to previously described transcripts lacking exon 3 or exon 4, or retaining part of intron 3, we identified novel isoforms lacking exon 2, along with isoforms lacking multiple exons. Isoforms lacking exon 2 were predominant at approximately 10% of TREM2 mRNA in the brain. Expression of TREM2 and frequency of exon 2 skipping did not differ between AD samples and non-AD controls (p = 0.1268 and p = 0.4909, respectively). Further, these novel splice isoforms were also observed across multiple tissues with similar frequency (range 5.3 - 13.0%). We found that the exon 2 skipped isoform D2-TREM2 is translated to protein and localizes similarly to full-length TREM2 protein, that both proteins are primarily retained in the Golgi complex, and that D2-TREM2 is expressed in AD and non-AD brain.

ConclusionSince the TREM2 ligand binding domain is encoded by exon 2, and skipping this exon retains reading frame while conserving localization, we hypothesize that D2-TREM2 acts as an inhibitor of TREM2 and targeting TREM2 splicing may be a novel therapeutic pathway for AD. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Measuring mitochondrial electron transfer complexes in previously frozen cardiac tissue from the offspring of sow: A model to assess exercise-induced mitochondrial bioenergetics changes Measuring the mitochondrial electron transfer complex (ETC) profile from previously frozen heart tissue samples from offspring born to an exercised sow provided descriptive data about exercise-induced mitochondrial biochemical changes in heart tissue from the offspring born to the exercised sow. The hypothesis proposed and tested was that a regular maternal exercise of a sow during pregnancy would increase the mitochondrial efficiency of offspring heart bioenergetics. This hypothesis was tested by isolating mitochondria using a mild-isolation procedure to assess mitochondrial ETC and supercomplex profiles. The procedure described here allowed for the processing of previously frozen archived heart tissues and eliminated the necessity of fresh mitochondria preparation for the assessment of mitochondrial ETC complexes, supercomplexes, and ETC complex activity profiles. This protocol describes the optimal ETC protein complex measurement in multiplexed antibody-based immunoblotting and super complex assessment using blue-native gel electrophoresis.

SUMMARYPreparation of mitochondria-enriched samples from previously frozen archived solid tissues allowed the investigators to perform both functional and analytical assessments of mitochondria in various experimental modalities. This study demonstrates how to prepare mitochondria-enriched preparations from frozen heart tissue and perform analytical assessments of mitochondria. 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. Changes in Behavior and Salivary Serum Amyloid A in cows with Subclinical Mastitis The aim of this study was to identify detailed changes in behavior, and in salivary serum amyloid A (SAA), associated with subclinical mastitis. This included standard sickness behaviors (which form part of an adaptive response to conserve energy, minimize heat loss and facilitate recovery following infection and inflammation) and lesser-studied behaviors, that are biologically predicted to change at low-level infection, and therefore particularly relevant for detecting subclinical mastitis (here labelled  luxury behaviors). SAA is a physiological marker of systemic inflammation, with levels in milk and serum already known to increase during subclinical mastitis. We investigated whether the same was true for SAA in cow saliva. Data were collected for n = 34 commercial barn-housed dairy cows, tested in 17 matched-pairs. Each pair comprised a cow with subclinical mastitis (SCM) and a healthy control (CTRL), identified using somatic cell count (SCC) (SCM: SCC >200 x1000 cells/ml; CTRL: SCC <100 x1000 cells/ml). SCM cows were selected for study ad-hoc, at which point they were paired with a CTRL cow, based upon parity and calving date; consequently, the full data set was accrued over several months. Data were collected for each pair over 3 days: SCC (l4:00-l5:00h) Day 1; behavior (24h from 00:00h) Day 2; salivary serum amyloid-A (SAA) Day 3. We report, for the first time, that an increase in salivary SAA occurs during subclinical mastitis; SAA was higher in SCM cows and demonstrated a positive (weak) correlation with SCC. The behavioral comparisons revealed that SCM cows had reductions in activity (behavioral transitions and distance moved), social exploration, social reactivity (here: likelihood to move away/be displaced following receipt of agonism), performance of social grooming and head butts, and the receipt of agonistic non-contact challenges. In addition, SCM cows received more head swipes, and spent a greater proportion of time lying with their head on their flank than CTRL cows. SCM cows also displayed an altered feeding pattern; they spent a greater proportion of feeding time in direct contact with two conspecifics, and a lower proportion of feeding time at self-locking feed barriers, than CTRL cows. Behavioral measures were found to correlate, albeit loosely, with serum SAA in a direction consistent with predictions for sickness behavior. These included positive correlations with both lying duration and the receipt of all agonistic behavior, and negative correlations with feeding, drinking, the performance of all social and all agonistic behavior, and social reactivity. We conclude that changes in salivary SAA, social behavior, and activity offer potential in the detection of subclinical mastitis and recommend further investigation to substantiate and refine our findings. Hotspots of pest-induced US urban tree death, 2020-2050 O_LIUrban trees are important nature-based solutions for future wellbeing and livability but are at high risk of mortality from insect pests. In the United States (US), 82% of the population live in urban settings and this number is growing, making urban tree mortality a matter of concern for most of its population. Until now, the magnitudes and spatial distributions of risks were unknown.
C_LIO_LIHere, we combine new models of street tree populations in [~]30,000 US communities, species-specific spread predictions for 57 invasive insect species, and estimates of tree death due to insect exposure for 48 host tree genera.
C_LIO_LIWe estimate that 1.4 million street trees will be killed by invasive insects from 2020 through 2050, costing an annualized average of US$ 30M. However, these estimates hide substantial variation: 23% of urban centers will experience 95% of all insect-induced mortality. Further, 90% of all mortality will be due to emerald ash borer (Agrilus planipennis, EAB), which is expected to kill virtually all ash trees (Fraxinus spp.) in >6000 communities.
C_LIO_LIWe define an EAB high-impact zone spanning 902,500km2, largely within the southern and central US, within which we predict the death of 98.8% of all ash trees. "Mortality hotspot cities" include Milwaukee, WI; Chicago, IL; and New York, NY.
C_LIO_LIWe identify Asian wood borers of maple and oak trees as the highest risk future invaders, where a new establishment could cost US$ 4.9B over 30 years.
C_LIO_LIPolicy implications: To plan effective mitigation, managers must know which tree species in which communities will be at the greatest risk, as well as the highest-risk insects. We provide the first country-wide, spatial forecast of urban tree mortality due to invasive insect pests. This framework identifies dominant pest insects and spatial impact hotspots, which can provide the basis for spatial prioritization of spread control efforts such as quarantines and biological control release sites. Our results highlight the need for EAB early-detection efforts as far from current infestations as Seattle, WA. Further, these findings produce a list of biotic and spatiotemporal risk factors for future high-impact US urban forest insect pests.
C_LI Untargeted metabolomic analyses reveal the diversity and plasticity of the specialized metabolome in seeds of different Camelina sativa genotypes Despite the essential role of Specialized Metabolites (SMs) in the interaction of plants with the environment, studying the ability of crop seeds to produce these protective compounds has been neglected. Furthermore, seeds produce a myriad of SMs providing an interesting model to investigate their diversity and plasticity. Camelina sativa gained a lot of interest in the past few years as rustic oil seed crop. A characterization of seed SM landscapes in six camelina genotypes grown in the field and harvested during five growing seasons has been undertaken in this work. This allowed a comprehensive annotation of seed SMs combining analyses that cluster SMs based on their chemical structures and co-accumulation patterns. These data showed broad effects of the environment on the stimulation of the seed-specialized metabolome. Among well annotated compounds, flavonols were identified as the metabolic class characterised by high plasticity, revealing significant variable accumulation according to the year and/or the genotype. Lastly, a deeper characterisation of primary metabolites and lipids in two selected genotypes has been performed. We showed that, in addition to flavonols, alkaloids and glucosinolates displayed a higher phenotypic plasticity with respect to most of the primary metabolites, including some sugars and major storage compounds such as fatty acids, proteins and most lipid classes (e.g. DAG, TAG), but similar plasticity compared to free aminoacids and carboxylic acids. This work highlighted major and unexplored effects of the environment on the seed specialized metabolome demonstrating that seeds exhibit a dynamic and plastic metabolism, with an impact on seed quality.

Significance statementSeeds produce a myriad of Specialized Metabolites (SMs) with an essential role in the interaction of plants with the environment. We characterized SM landscapes, primary metabolites and lipid composition in the seeds of camelina genotypes grown in the open field in five consecutive growing seasons. Our results showed the predominant effect of the environment on the regulation of the seed - specialized metabolome, with a potential impact on seed quality of camelina that may also occur in other oilseed crops. The polyglutamine amyloid nucleus in living cells is a monomer with competing dimensions of order A long-standing goal of the study of amyloids has been to characterize the structural basis of the rate-determining nucleating event. However, the ephemeral nature of that event has made it inaccessible to classical biochemistry, structural biology, and computational approaches. Here, we addressed that limitation by measuring the dependence of amyloid formation on concentration and conformational templates in living cells, whose volumes are sufficiently small to resolve the outcomes of independent nucleation events. We characterized numerous rationally designed sequence variants of polyglutamine (polyQ), a polypeptide that precipitates Huntingtons and other amyloid-associated neurodegenerative diseases when its length exceeds a characteristic threshold. This effort uncovered a pattern of approximately twelve Qs, only for polypeptides exceeding the clinical length threshold, that allow for amyloid nucleation to occur spontaneously within single polypeptides. Nucleation was inhibited by intermolecular phase separation. Using atomistic molecular dynamics simulations, we found that the pattern encodes a minimal steric zipper of interdigitated side chains. Lateral growth of the steric zipper competed with axial growth to produce "pre-amyloid" oligomers. By illuminating the structural mechanism of polyQ amyloid formation in cells, our findings reveal a potential molecular etiology for polyQ diseases, and may provide a roadmap for the design of new therapies. Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable and highly reproducible in terms of microbial composition, function and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation. The first sheep graph-based pan-genome reveals the spectrum of structural variations and their effects on tail phenotypes Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here, we assembled 26 haplotype-resolved genome assemblies from 13 genetically diverse sheep using PacBio HiFi sequencing. We constructed a graph-based ovine pan-genome and discovered 142,422 biallelic insertions and deletions, 7,028 divergent alleles and 13,419 multiallelic variations. We then used a graph-based approach to genotype the biallelic SVs in 684 individuals from 45 domestic breeds and two wild species. Integration with RNA-seq data allows to identify candidate expression-associated SVs. We demonstrate a direct link of SVs and phenotypes by localizing the putative causative insertion in HOXB13 gene responsible for the long-tail trait and identifying multiple large SVs associated with the fat-tail. Beyond generating a benchmark resource for ovine structural variants, our study highlights that animal genetic research will greatly benefit from using a pan-genome graph rather than a single reference genome. Circulating multimeric immune complexes drive immunopathology in COVID-19 A dysregulated immune response with high levels of SARS-CoV-2 specific IgG antibodies characterizes patients with severe or critical COVID-19. Although a robust IgG response is traditionally considered to be protective, excessive triggering of activating Fc-gamma-receptors (Fc{gamma}Rs) could be detrimental and cause immunopathology. Here, we document that patients who develop soluble circulating IgG immune complexes (sICs) during infection are subject to enhanced immunopathology driven by Fc{gamma}R activation. Utilizing cell-based reporter systems we provide evidence that sICs are predominantly formed prior to a specific humoral response against SARS-CoV-2. sIC formation, together with increased afucosylation of SARS-CoV-2 specific IgG eventually leads to an enhanced CD16 (Fc{gamma}RIII) activation of immune cells reaching activation levels comparable active systemic lupus erythematosus (SLE) disease. Our data suggest a vicious cycle of escalating immunopathology driven by an early formation of sICs in predisposed patients. These findings reconcile the seemingly paradoxical findings of high antiviral IgG responses and systemic immune dysregulation in severe COVID-19.

Clinical implicationsThe identification of sICs as drivers of an escalating immunopathology in predisposed patients opens new avenues regarding intervention strategies to alleviate critical COVID-19 progression.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=78 SRC="FIGDIR/small/449893v4_ufig1.gif" ALT="Figure 1">
View larger version (17K):
org.highwire.dtl.DTLVardef@b4616corg.highwire.dtl.DTLVardef@682d1aorg.highwire.dtl.DTLVardef@16946cborg.highwire.dtl.DTLVardef@a6ef7d_HPS_FORMAT_FIGEXP  M_FIG C_FIG A vicious cycle of immunopathology in COVID-19 patients is driven by soluble multimeric immune complexes (sICs). SARS-CoV-2 infection triggers sIC formation in prone individuals. Activation of Fc{gamma}RIII/CD16 expressing immune cells by sICs precedes a humoral response to SARS-CoV2 infection. sICs and infection add to IgG afucosylation, further enhancing Fc{gamma}RIII/CD16 activation by opsonized targets. High inflammation induces further sIC mediated immune cell activation ultimately leading to an escalating immunopathology. Impaired activation of Transposable Elements in SARS-CoV-2 infection Transposable element (TE) transcription is induced in response to viral infections. TE induction triggers a robust and durable interferon (IFN) response, providing a host defense mechanism. Still, the connection between SARS-CoV-2 IFN response and TEs remains largely unknown. Here, we analyzed TE expression changes in response to SARS-CoV-2 infection in different human cellular models. We find that compared to other viruses, which cause global upregulation of TEs, SARS-CoV-2 infection results in a significantly milder TE response in both primary lung epithelial cells and in iPSC-derived lung alveolar type 2 cells. TE activation precedes, and correlates with, the induction of IFN-related genes, suggesting that the limited activation of TEs following SARS-CoV-2 infection may be the reason for the weak IFN response. Diminished TE activation was not observed in lung cancer cell lines with very high viral load. Moreover, we identify two variables which explain most of the observed diverseness in immune responses: basal expression levels of TEs in the pre-infected cells, and the viral load. Finally, analyzing the SARS-CoV-2 interactome, as well as the epigenetic landscape around the TEs that are activated following infection, we identify SARS-CoV-2 interacting proteins, which may regulate chromatin structure and TE transcription in response to a high viral load. This work provides a functional explanation for SARS-CoV-2s success in its fight against the host immune system, and suggests that TEs could be used as sensors and serve as potential drug targets for COVID-19.

Key points O_LIUnlike other viruses, SARS-CoV-2 invokes a weak and inefficient transposable element (TE) response
C_LIO_LITE induction precedes and predicts IFN response
C_LIO_LIBasal TE expression and viral load explain immune responses
C_LIO_LIDistinct chromatin and enhancer binding factors occupancy on TEs induced by SARS-CoV-2
C_LI Spatio-temporal, optogenetic control of gene expression in organoids Organoids derived from stem cells become increasingly important to study human development and to model disease. However, methods are needed to control and study spatio-temporal patterns of gene expression in organoids. To this aim, we combined optogenetics and gene perturbation technologies to activate or knock-down RNA of target genes, at single-cell resolution and in programmable spatio-temporal patterns. To illustrate the usefulness of our approach, we locally activated Sonic Hedgehog (SHH) signaling in an organoid model for human neurodevelopment. High-resolution spatial transcriptomic and single-cell analyses showed that this local induction was sufficient to generate stereotypically patterned organoids in three dimensions and revealed new insights into SHHs contribution to gene regulation in neurodevelopment.

With this study, we propose optogenetic perturbations in combination with spatial transcriptomics as a powerful technology to reprogram and study cell fates and tissue patterning in organoids. Multivariate Genomewide Association Analysis by Iterative Hard Thresholding 1In genome-wide association studies (GWAS), analyzing multiple correlated traits is potentially superior to conducting multiple univariate analyses. Standard methods for multivariate GWAS operate marker-by-marker and are computationally intensive. We present a penalized regression algorithm for multivariate GWAS based on iterative hard thresholding (IHT) and implement it in a convenient Julia package MendelIHT.jl (https://github.com/OpenMendel/MendelIHT.jl). In simulation studies with up to 100 traits, IHT exhibits similar true positive rates, smaller false positive rates, and faster execution times than GEMMAs linear mixed models and mv-PLINKs canonical correlation analysis. On UK Biobank data, our IHT software completed a 3-trait joint analysis in 20 hours and an 18-trait joint analysis in 53 hours, requiring up to 80GB of computer memory. In short, our software enables geneticists to fit a single regression model that simultaneously considers the effect of all SNPs and dozens of traits. Cellular-level phenotyping of tumor-immune microenvironment (TiME) in patients in vivo reveals distinct inflammation and endothelial anergy signatures Immunotherapies have shown unprecedented clinical benefits in several malignancies1-3. However, clinical responses remain variable and unpredictable, indicating the need to develop predictive platforms that can improve patient stratification4. Phenotyping of tumors into hot, altered, or cold5 based on assessment of only T-cell infiltration in static tumor biopsies provides suboptimal prediction of immunotherapy response6,7. In vivo dynamic mechanisms within the tumor microenvironment such as tumor angiogenesis and leukocyte trafficking5,8,9 also play a central role in modulating anti-tumor immunity and therefore immunotherapy response. Here, we report novel tumor immune microenvironment (TiME) phenotyping in vivo in patients with non-invasive spatially-resolved cellular-level imaging based on endogenous contrast. Investigating skin cancers as a model, with reflectance confocal microscopy (RCM) imaging10, we determined four major phenotypes with variable prevalence of vasculature (Vasc) and inflammation (Inf) features: VaschiInfhi, VaschiInflo, VascloInfhi and Vascmed/hiInflo. The VaschiInfhi phenotype correlates with high immune activation, exhaustion, and vascular signatures while VaschiInflo with endothelial anergy and immune exclusion. Automated quantification of TiME features demonstrates moderate-high accuracy and correlation with corresponding gene expression. Prospectively analyzed response to topical immunotherapy show highest response in VascloInfhi, and reveals the added value of vascular features in predicting treatment response. Our novel in vivo cellular-level imaging and phenotyping approach can potentially advance our fundamental understanding of TiME, develop robust predictors for immunotherapy outcomes and identify novel targetable pathways in future. Game of theories: an agent-based model of how variation in contest cost-benefits determines the evolution of contest resolution rules Game-theory based models are used to understand rules that animals use to settle contests over indivisible resources. However, the empirical literature of contests indicates controversial support to models, with some species supporting different models and other species showing no support to any model. Since strategies used to resolve contests may have different associated costs, it is possible that different conditions have determined the evolution of distinct assessment strategies used by animals. We used an agent-based model to explore the importance of the following conditions: resource availability, probability of reproduction with resource, and damage costs on evolution of assessment strategies. We used self- and mutual-assessment models as a heurist framework to build agents with different assessment strategies. In our model, agents competed for resources in scenarios with different combinations of resource availability, probability of reproduction with resource, and damage costs. We found that agents following the self-assessment with damage strategy were prevalent in scenarios with no probability of reproduction without the resource, independently of other variables. We also found that agents following the non-aggressive strategy occurred in all scenarios. However, agents using the non-aggressive strategy were prevalent only in scenarios with probability of reproduction with the resource. Finally, we observed that agents using mutual-assessment occurred only in a scenario with high risk of damage, low availability of resources, and with probability of reproduction without the resource. These results indicate that agents following the self-assessment with damage and non-aggressive strategies may be able to stay at most scenarios. Inhibitors of ERp44, PDIA1, and AGR2 induce disulfide-mediated oligomerization of Death Receptors 4 and 5 and cancer cell death Breast cancer mortality remains unacceptably high, indicating a need for safer and more effective therapeutic agents. Disulfide bond Disrupting Agents (DDAs) were previously identified as a novel class of anticancer compounds that selectively kill cancers that overexpress the Epidermal Growth Factor Receptor (EGFR) or its family member HER2. DDAs kill EGFR+ and HER2+ cancer cells via the parallel downregulation of EGFR, HER2, and HER3 and activation/oligomerization of Death Receptors 4 and 5 (DR4/5). However, the mechanisms by which DDAs mediate these effects are unknown. Affinity purification analyses employing biotinylated-DDAs reveal that the Protein Disulfide Isomerase (PDI) family members AGR2, PDIA1, and ERp44 are DDA target proteins. Further analyses demonstrate that shRNA-mediated knockdown of AGR2 and ERp44, or expression of ERp44 mutants, enhance basal and DDA-induced DR5 oligomerization. DDA treatment of breast cancer cells disrupts PDIA1 and ERp44 mixed disulfide bonds with their client proteins. Together, the results herein reveal DDAs as the first small molecule, active site inhibitors of AGR2 and ERp44, and demonstrate roles for AGR2 and ERp44 in regulating the activity, stability, and localization of DR4 and DR5, and activation of Caspase 8. Mouse frontal cortex mediates additive multisensory decisions To interpret the world and make accurate perceptual decisions, the brain must combine information across sensory modalities. For instance, it must combine vision and hearing to localize objects based on their image and sound. Probability theory suggests that evidence from multiple independent cues should be combined additively, and several studies support this hypothesis1-14. However, other studies suggest that mammals, particularly mice, can exhibit sensory bias during audiovisual integration15-19; furthermore, the cortical substrates of multisensory integration remain uncertain20. Here we show that to localize a stimulus, mice combine auditory and visual spatial cues additively, a computation supported by additive multisensory integration in frontal cortex. We developed an audiovisual localization task where mice turn a wheel to indicate the joint position of an image and a sound. Scanning optogenetic inactivation21-24 during stimulus presentation showed that visual cortex contributes unisensory information whereas frontal cortex contributes multisensory information to the mouses decision. Neuropixels recordings of >10,000 neurons indicated that neural activity in frontal area MOs (secondary motor cortex) reflects an additive combination of visual and auditory signals. An accumulator model2,25,26 applied to the sensory representations of MOs neurons reproduced both the observed choices and reaction times. Sensory signals in MOs were specific to trained mice, and indeed the same accumulator applied to naive mice failed to perform the task. These results suggest that frontal cortex adapts through learning to integrate information across sensory cortices, providing a signal that is then transformed into a binary decision by a downstream accumulator.