Tourmaline: a containerized workflow for rapid and iterable amplicon sequence analysis using QIIME 2 and Snakemake BackgroundAmplicon sequencing (metabarcoding) is a common method to survey diversity of environmental communities whereby a single genetic locus is amplified and sequenced from the DNA of whole or partial organisms, organismal traces (e.g., skin, mucus, feces), or microbes in an environmental sample. Several software packages exist for analyzing amplicon data, among which QIIME 2 has emerged as a popular option because of its broad functionality, plugin architecture, provenance tracking, and interactive visualizations. However, each new analysis requires the user to keep track of input and output file names, parameters, and commands; this lack of automation and standardization is inefficient and creates barriers to meta-analysis and sharing of results.

FindingsWe developed Tourmaline, a Python-based workflow that implements QIIME 2 and is built using the Snakemake workflow management system. Starting from a configuration file that defines parameters and input files--a reference database, a sample metadata file, and a manifest or archive of FASTQ sequences--it uses QIIME 2 to run either the DADA2 or Deblur denoising algorithm, assigns taxonomy to the resulting representative sequences, performs analyses of taxonomic, alpha, and beta diversity, and generates an HTML report summarizing and linking to the output files. Features include support for multiple cores, automatic determination of trimming parameters using quality scores, representative sequence filtering (taxonomy, length, abundance, prevalence, or ID), support for multiple taxonomic classification and sequence alignment methods, outlier detection, and automated initialization of a new analysis using previous settings. The workflow runs natively on Linux and macOS or via a Docker container. We ran Tourmaline on a 16S rRNA amplicon dataset from Lake Erie surface water, showing its utility for parameter optimization and the ability to easily view interactive visualizations through the HTML report, QIIME 2 viewer, and R- and Python-based Jupyter notebooks.

ConclusionsAutomated workflows like Tourmaline enable rapid analysis of environmental and biomedical amplicon data, decreasing the time from data generation to actionable results. Tourmaline is available for download at github.com/aomlomics/tourmaline.
MMMVI: Detecting SARS-CoV-2 Variants of Concern in Metagenomic Wastewater Samples MotivationSARS-CoV-2 is the causative agent of the COVID-19 pandemic. Variants of Concern (VOCs) and Variants of Interest (VOIs) are lineages that represent a greater risk to public health, and can be differentiated from the wildtype virus based on unique profiles of signature mutations. Metagenomic sequence analysis of wastewater represents an emerging form of surveillance that can capture early signal for these variants in a community prior to detection through public health testing or genomic surveillance activities. However, because multiple viral genomes are likely to be present in a metagenomic sample, additional analytical scrutiny of the sequencing reads beyond variant calling is needed to increase confidence in diagnostic determinations.

ResultsWhere multiple signature mutations are present on a given read, they can be used as enhanced biomarkers to confirm the presence of a VOC/VOI in the sample. We have developed MMMVI, a tool to aggregate and report on the likely presence of a VOC/VOI in a set of metagenomic reads based on the detection of reads bearing multiple signature mutations.

AvailabilityMMMVI is implemented in Python, and is available under the MIT licence from https://github.com/dorbarker/voc-identify/

Contactdillon.barker@canada.ca
Consequences of alanine-126 mutations in helix-3 on structure and functions of Rad6 E2 ubiquitin-conjugating enzymes Rad6, an E2 ubiquitin-conjugating enzyme conserved from yeast to humans, functions in transcription, genome maintenance and proteostasis. The contributions of many conserved secondary structures of Rad6 and its human homologs UBE2A and UBE2B to their biological functions are not understood. A mutant RAD6 allele with a missense substitution at alanine-126 (A126) of helix-3 that causes defects in telomeric gene silencing, DNA repair and protein degradation was reported over two decades ago. Here, using a combination of genetics, biochemical, biophysical, and computational approaches, we discovered that helix-3 A126 mutations compromise the ability of Rad6 to ubiquitinate target proteins without disrupting interactions with partner E3 ubiquitin-ligases that are required for their various biological functions in vivo. Explaining the defective in vitro or in vivo ubiquitination activities, molecular dynamics simulations and NMR showed that helix-3 A126 mutations cause local disorder of the catalytic pocket of Rad6 in addition to disorganizing the global structure of the protein to decrease its stability in vivo. We also show that helix-3 A126 mutations deform the structures of UBE2A and UBE2B, the human Rad6 homologs, and compromise the in vitro ubiquitination activity and folding of UBE2B. In summary, our studies reveal that the conserved helix-3 is a crucial structural constituent that controls the organization of catalytic pockets, enzymatic activities, and biological functions of the Rad6-family E2 ubiquitin-conjugating enzymes. Network-level modeling of stimulation-induced functional connectivity change: An optogenetic study in non-human primate cortex As aberrant network-level functional connectivity underlies a variety of neural disorders, the ability to induce targeted functional reorganization would be a profound development towards therapies for neural disorders. Brain stimulation has been shown to alter large-scale network-wide functional connectivity, but the mapping from stimulation to the modification is unclear. Here, we leverage advances in neural interfaces, interpretable machine learning, and graph theory to arrive at a model which accurately predicts stimulation-induced network-wide functional reorganization. The model jointly considers the stimulation protocol and the cortical network structure, departing from the standard approach which only considers the stimulation protocol. We validate our approach in the primary sensorimotor cortex of non-human primates using paired optogenetic stimulation through a large-scale optogenetic interface. We observe that the stimulation protocol only predicts a small portion of the induced functional connectivity changes while the network structure predicts much more, indicating that the underlying network is the primary mediator of the response to stimulation. We extract the relationships linking the stimulation and network characteristics to the functional connectivity changes and observe that the mappings diverge over frequency bands and successive stimulations. Finally, we uncover shared processes governing real-time and longer-term effects of stimulation. Our framework represents a paradigm shift for targeted neural stimulation and can be used to interrogate, improve, and develop stimulation-based interventions for neural disorders.

TeaserBrain stimulation rewires the brain, but the pre-existing network structure of the brain controls the rewiring. 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. Long term analysis of social structure: evidence of age-based consistent associations in male Alpine ibex Despite its recognized importance for understanding the evolution of animal sociality as well as for conservation, long term analysis of social networks of animal populations is still relatively uncommon. We investigated social network dynamics in males of a gregarious mountain ungulate (Alpine ibex, Capra ibex) over ten years focusing on groups, sub-groups and individuals, exploring the dynamics of sociality over different scales. Despite the social structure changing between seasons, the Alpine ibex population was highly cohesive: fission-fusion dynamics lead almost every male in the population to associate with each other male at least once. Nevertheless, we found that male Alpine ibex showed preferential associations that were maintained across seasons and years. Age seemed to be the most important factor driving preferential associations while other characteristics, such as social status, appeared less crucial. We also found that centrality measures were influenced by age and were also related to individual physical condition. The multi-scale and long-term frame of our study helped us show that ecological constrains, such as resource availability, may play a role in shaping associations in a gregarious species, but they cannot solely explain sociality and preferential association that are likely also to be driven by life-history linked physiological and social needs. Our results highlight the importance of long-term studies based on individually recognizable subjects to help us build on our understanding of the evolution of animal sociality. X chromosomes show relaxed selection and complete somatic dosage compensation across Timema stick insect species Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles, and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here we use wholegenome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 11% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution. Shelterin components modulate nucleic acids condensation and phase separation in the context of telomeric DNA Telomeres are nucleoprotein complexes that protect the ends of chromosomes and are essential for chromosome stability in Eukaryotes. In cells, individual telomeres form distinct globules of finite size that appear to be smaller than expected for bare DNA. Moreover, upon changes in their protein composition, telomeres can cluster to form telomere-induced-foci (TIFs) or co-localize with promyelocytic leukemia (PML) nuclear bodies. The physical basis for collapse of individual telomeres and coalescence of multiple ones remains unclear, as does the relationship between these two phenomena. By combining single-molecule measurements, optical microscopy, turbidity assays, and simulations, we show that the telomere scaffolding protein TRF2 can condense individual DNA chains and drives coalescence of multiple DNA molecules, leading to phase separation and the formation of liquid-like droplets. Addition of the TRF2 binding protein hRap1 modulates phase boundaries and tunes the specificity of solution demixing while simultaneously altering the degree of DNA compaction. Our results suggest that the condensation of single telomeres and formation of biomolecular condensates containing multiple telomeres are two different outcomes driven by the same set of molecular interactions. Moreover, binding partners, such as other telomere components, can alter those interactions to promote single-chain DNA compaction over multiple-chain phase separation. RNAi screen of RING/U-box domain ubiquitin ligases identifies critical regulators of tissue regeneration in planarians Regenerative processes depend on the interpretation of signals to coordinate cell behaviors. The role of ubiquitin-mediated signaling is known to be important in many cellular and biological contexts, but its role in regeneration is not well understood. To investigate how ubiquitylation impacts tissue regeneration in vivo, we are studying planarians that are capable of regenerating after nearly any injury using a population of stem cells. Here we used RNAi to screen RING/U-box E3 ubiquitin ligases that are highly expressed in planarian stem cells and stem cell progeny. RNAi screening identified nine genes with functions in regeneration, including the spliceosomal factor prpf19 and histone modifier rnf2; based on their known roles in developmental processes, we further investigated these two genes. We found that prpf19 was required for animal survival but not for stem cell maintenance, suggesting a role in promoting cell differentiation. Because RNF2 is the catalytic subunit of the Polycomb Repressive Complex 1 (PRC1), we also examined other putative members of this complex (CBX and PHC). We observed a striking phenotype of regional tissue misspecification in cbx and phc RNAi planarians. To identify genes regulated by PRC1, we performed RNA-seq after knocking down rnf2 or phc. Although these proteins are predicted to function in the same complex, we found that the set of genes differentially expressed in rnf2 versus phc RNAi were largely non-overlapping. Using in situ hybridization, we showed that rnf2 regulates gene expression levels within a tissue type, whereas phc is necessary for the spatial restriction of gene expression, findings consistent with their respective in vivo phenotypes. This work not only uncovered roles for RING/U-box E3 ligases in stem cell regulation and regeneration, but also identified differential gene targets for two putative PRC1 factors required for maintaining cell-type-specific gene expression in planarians. Establishing combination PAC-1 and TRAIL regimens for treating ovarian cancer based on patient-specific pharmacokinetic profiles using in silico clinical trials Ovarian cancer is commonly diagnosed in its late stages, and new treatment modalities are needed to improve patient outcomes and survival. We have recently established the synergistic effects of combination tumour necrosis factor-related apoptosis-inducing ligand (TRAIL) and procaspase activating compound (PAC-1) therapies in granulosa cell tumours (GCT) of the ovary, a rare form of ovarian cancer, using a mathematical model of the effects of both drugs in a GCT cell line. Here, to understand the mechanisms of combined TRAIL and PAC-1 therapy, study the viability of this treatment strategy, and accelerate preclinical translation, we leveraged our mathematical model in combination with population pharmacokinetics (PopPK) models of both TRAIL and PAC-1 to expand a realistic heterogeneous cohort of virtual patients and optimize treatment schedules. Using this approach, we investigated treatment responses in this virtual cohort and determined optimal therapeutic schedules based on patient-specific pharmacokinetic characteristics. Our results showed that schedules with high initial doses of PAC-1 were required for therapeutic efficacy. Further analysis of individualized regimens revealed two distinct groups of virtual patients within our cohort: one with high PAC-1 elimination, and one with normal PAC-1 elimination. In the high elimination group, high weekly doses of both PAC-1 and TRAIL were necessary for therapeutic efficacy, however virtual patients in this group were predicted to have a worse prognosis when compared to those in the normal elimination group. Thus, PAC-1 pharmacokinetic characteristics, particularly clearance, can be used to identify patients most likely to respond to combined PAC-1 and TRAIL therapy. This work underlines the importance of quantitative approaches in preclinical oncology. 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. 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. Artificial seawater based long-term culture of colonial ascidians Tunicates are highly diverse marine invertebrate filter-feeders that are vertebrates closest relatives. These organisms, despite a drastically different body plan during their adulthood, have a tissue complexity related to that of vertebrates. Ascidians, which compose most of the Tunicata, are benthic sessile hermaphrodites that reproduce sexually through a motile tadpole larval stage. Over half of the known ascidians species are able to reproduce asexually by budding, typically leading to the formation of colonies where animals, called zooids, are interconnected through an external vascular system. In addition, colonial ascidians are established models for important biological processes including allorecognition, immunobiology, aging, angiogenesis and whole-body regeneration. However, the current paucity in breeding infrastructures limits the study of these animals to coastal regions.

To promote a wider scientific spreading and popularity of colonial ascidians, we have developed a flexible recirculating husbandry setup for their long-term in-lab culture. Our system is inspired both by the flow-through aquariums used by coastal ascidian labs, as well as by the recirculating in-lab systems used for zebrafish research. Our hybrid system thus combines colony breeding, water filtering and food culturing in a semi-automated system where specimens develop on hanging microscopy glass slides. Temperature, light/dark cycles, flow speed and feeding rates can be controlled independently in four different breeding environments to provide room for species-specific optimization as well as for running experiments. This setup is complemented with a quarantine for the acclimatization of wild isolates.

Herein we present our success in breeding Botrylloides diegensis, a species of colonial ascidians, for more than 3 years in recirculating artificial seawater over 600 km away from their natural habitat. We show that colonies adapt well to in-lab culturing provided that a suitable marine microbiome is present, and that a specific strain can be isolated, propagated and efficiently used for research over prolonged periods of time. The flexible and modular structure of our system can be scaled and adapted to the needs of specific species, such as Botryllus schlosseri, as well as of particular laboratory spaces. Overall, we show that Botrylloides diegensis can be proficiently bred in-land and suggest that our results can be extended to other species of colonial ascidians to promote research on these fascinating animals.

HighlightsO_LIFirst artificial seawater based recirculating aquaculture for colonial ascidians
C_LIO_LIOver 3 years of continuous breeding
C_LIO_LISemi-automated setup with minimized maintenance
C_LIO_LIGood biomass production for strain propagation
C_LIO_LI4 different culture conditions for optimized breeding for species of interest
C_LI Vapor exposure to Δ9-tetrahydrocannabinol (THC) slows locomotion of the Maine Lobster (Homarus americanus) RationaleDespite a long history of use in synaptic physiology, the lobster has been a neglected model for behavioral pharmacology. A restaurateur proposed that exposing lobster to cannabis smoke reduces anxiety and pain during the cooking process. It is unknown if lobster gill respiration in air would result in significant {Delta}9-tetrahydrocannabinol (THC) uptake and whether this would have any detectable behavioral effects.

ObjectiveThe primary goal was to determine tissue THC levels in the lobster after exposure to THC vapor. Secondary goals were to determine if THC vapor altered locomotor behavior or nociception.

MethodsTissue samples were collected (including muscle, brain and hemolymph) from Homarus americanus (N=3 per group) following 30 or 60 minutes of exposure to vapor generated by an e-cigarette device using THC (100 mg/mL in a propylene glycol vehicle). Separate experiments assessed locomotor behavior and hot water nociceptive responses following THC vapor exposure.

ResultsTHC vapor produced duration-related THC levels in all tissues examined. Locomotor activity was decreased (distance, speed, time-mobile) by 30 min inhalation of THC. Lobsters exhibit a temperature-dependent withdrawal response to immersion of tail, antennae or claws in warm water; this is novel evidence of thermal nociception for this species. THC exposure for 60 minutes had only marginal effect on nociception under the conditions assessed.

ConclusionsVapor exposure of lobsters, using an e-cigarette based model, produces dose-dependent THC levels in all tissues and reduces locomotor activity. Hot water nociception was temperature dependent, but only minimal anti-nociceptive effect of THC exposure was confirmed. The cardiac TRPA1 channel drives calcium-mediated mechano-arrhythmogenesis SUMMARY PARAGRAPHMaintenance of cardiac function involves a regulatory loop in which electrical excitation causes the heart to contract through excitation-contraction coupling (ECC),1 and the mechanical state of the heart directly affects its electrical activity through mechano-electric coupling (MEC).2 However, in pathological states such as acute ischaemia that alter early or late electro-mechanical coordination (i.e., disturbances in ECC or repolarisation-relaxation coupling, RRC), MEC may contribute to the initiation and / or sustenance of arrhythmias (mechano-arrhythmogenesis).3 The molecular identity of specific factor(s) underlying mechano-arrhythmogenesis in acute ischaemia, however, remain undefined.4 By rapid stretch of rabbit single left ventricular cardiomyocytes, we show that upon ATP-sensitive potassium channel-induced alterations of RRC, overall vulnerability to mechano-arrhythmogenesis is increased, with mechano-sensitive5-11 transient receptor potential kinase ankyrin 1 (TRPA1) channels12 acting as the molecular driver through a Ca2+-mediated mechanism. Specifically, TRPA1 activation drives stretch-induced excitation and creates a substrate for self-sustained arrhythmias, which are maintained by increased cytosolic free Ca2+ concentration ([Ca2+]i) and spontaneous [Ca2+]i oscillations. This TRPA1-dependent mechano-arrhythmogenesis involves microtubules, and can be prevented by block of TRPA1 or buffering of [Ca2+]i. Thus, in cardiac pathologies with disturbed RRC dynamics and / or augmented TRPA1 activity, TRPA1 may represent an anti-arrhythmic target with untapped therapeutic potential.13-17 High-throughput variant detection using a color-mixing strategy Many diseases are related to multiple genetic alterations along a single gene. Probing for highly multiple (>10) variants in a single qPCR tube is not possible due to a limited number of fluorescence channels and one variant per channel, so many more tubes are needed. Here, we experimentally validate our novel color-mixing strategy that uses fluorescence combinations as digital color codes to probe multiple variants simultaneously. The color-mixing strategy relies on a simple intra-tube assay that can probe for 15 variants as part of an inter-tube assay that can probe for an exponentially increased number of variants. The color-mixing strategy is achieved using multiplex double-stranded toehold probes modified with fluorophores and quenchers; the probes are designed to be quenched or luminous after binding to wildtype or variant templates. We used the color-mixing strategy to probe for 21 pathogenic mutations in thalassemia and to distinguish between heterozygous and homozygous variants in 6 tubes, with a specificity of 99% and a sensitivity of 94%. To support tuberculosis diagnosis, we used the same strategy to simultaneously probe in Mycobacterium tuberculosis for rifampicin-resistance mutations occurring within one 81-bp region and one 48-bp region in rpoB gene, plus five isoniazid-resistance mutations in inhA and katG genes. 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. Identification of COVID-19 prognostic markers and therapeutic targets through meta-analysis and validation of Omics data from nasopharyngeal samples While our battle with the COVID-19 pandemic continues, a multitude of Omics data has been generated from patient samples in various studies, which remains to be translated. We conducted a meta-analysis of published transcriptome and proteome profiles of nasal swab and bronchioalveolar lavage fluid (BALF) samples of COVID-19 patients, to shortlist high confidence upregulated host factors. Subsequently, mRNA overexpression of selected genes was validated in nasal swab/BALF samples from a cohort of COVID-19 positive/negative, symptomatic/asymptomatic individuals. Analysis of these data revealed S100 family genes (S100A6, S100A8, S100A9, and S100P) as prognostic markers of COVID-19 disease. Furthermore, Thioredoxin gene (TXN) was identified as a significant upregulated host factor in our overlap analysis. An FDA-approved drug Auranofin, which inhibits Thioredoxin reduction, was found to mitigate SARS-CoV-2 replication in vitro and in vivo in the hamster challenge model. Overall, this study translates COVID-19 host response Big Data into potential clinical interventions. Plant inositol-phosphate-glycans and a fucosylated xyloglucan oligosaccharide are accumulated upon Arabidopsis thaliana/ Botrytis cinerea infection In mammals, insulin is involved in controlling blood glucose levels and its role in modulating immunity is being more and more documented. This hormone promotes the release of inositolphosphate glycans (IPG) which act as mediators. In plants, one IG has already been identified in plant culture cells (Smith and Fry, 1999; Smith et al., 1999) but, to our knowledge, no IPG have been yet identified. Here, we discovered 7 IPG that are accumulated upon Arabidopsis thaliana-Botrytis cinerea interaction, concomitantly with oligogalacturonides and a fucosylated xyloglucan oligosaccharide. Further structural characterization showed that they come from the hydrolysis of polar heads of Serie A to H glycosylinositol phosphorylceramides presumably via a phospholipase C activity. Taken together with the emerging role of insulin as immune regulator, these results question the role of IPG as damage associated molecular pattern both in animal and plant kingdoms. Protein cost minimization promotes the emergence of coenzyme redundancy Coenzymes distribute a variety of chemical moieties throughout cellular metabolism, participating in group (e.g., phosphate, acyl) and electron transfer. For a variety of reactions requiring acceptors or donors of specific resources, there often exist degenerate sets of molecules (e.g., NAD(H) and NADP(H)) that carry out similar functions. Although the physiological roles of various coenzyme systems are well established, it is unclear what selective pressures may have driven the emergence of coenzyme redundancy. Here we use genome-wide metabolic modeling approaches to decompose the selective pressures driving enzymatic specificity for either NAD(H) or NADP(H) in the metabolic network of Escherichia coli. We found that few enzymes are thermodynamically constrained to using a single coenzyme, and in principle, a metabolic network relying on only NAD(H) is feasible. However, structural and sequence analyses revealed widespread conservation of residues that retain selectivity for either NAD(H) or NADP(H), suggesting that additional forces may shape specificity. Using a model accounting for the cost of oxidoreductase enzyme expression, we found that coenzyme redundancy universally reduces the minimal amount of protein required to catalyze coenzyme-coupled reactions, inducing individual reactions to strongly prefer one coenzyme over another when reactions are near thermodynamic equilibrium. We propose that protein minimization generically promotes coenzyme redundancy, and that coenzymes typically thought to exist in a single pool (e.g., CoA) may exist in more than one form (e.g., dephospho-CoA).

Significance statementMetabolism relies on a small class of molecules (coenzymes) that serve as universal donors and acceptors of key chemical groups and electrons. Although metabolic networks crucially depend on structurally redundant coenzymes (e.g., NAD(H) and NADP(H)) associated with different enzymes, the criteria that led to the emergence of this redundancy remain poorly understood. Our combination of modeling, and structural and sequence analysis indicates that coenzyme redundancy is not essential for metabolism, but rather an evolved strategy promoting efficient usage of enzymes when biochemical reactions are near equilibrium. Our work suggests that early metabolism may have operated with fewer coenzymes, and that adaptation for metabolic efficiency may have driven the rise of coenzyme diversity in living systems. A human antibody reveals a conserved site on beta-coronavirus spike proteins and confers protection against SARS-CoV-2 infection Broadly neutralizing antibodies (bnAbs) to coronaviruses (CoVs) are valuable in their own right as prophylactic and therapeutic reagents to treat diverse CoVs and, importantly, as templates for rational pan-CoV vaccine design. We recently described a bnAb, CC40.8, from a coronavirus disease 2019 (COVID-19)-convalescent donor that exhibits broad reactivity with human beta-coronaviruses ({beta}-CoVs). Here, we showed that CC40.8 targets the conserved S2 stem-helix region of the coronavirus spike fusion machinery. We determined a crystal structure of CC40.8 Fab with a SARS-CoV-2 S2 stem-peptide at 1.6 [A] resolution and found that the peptide adopted a mainly helical structure. Conserved residues in {beta}-CoVs interacted with CC40.8 antibody, thereby providing a molecular basis for its broad reactivity. CC40.8 exhibited in vivo protective efficacy against SARS-CoV-2 challenge in two animal models. In both models, CC40.8-treated animals exhibited less weight loss and reduced lung viral titers compared to controls. Furthermore, we noted CC40.8-like bnAbs are relatively rare in human COVID-19 infection and therefore their elicitation may require rational structure-based vaccine design strategies. Overall, our study describes a target on {beta}-CoV spike proteins for protective antibodies that may facilitate the development of pan-{beta}-CoV vaccines.

SUMMARYA human mAb isolated from a COVID-19 donor defines a protective cross-neutralizing epitope for pan-{beta}-CoV vaccine design strategies Cell type-specific inference of differential expression in spatial transcriptomics Spatial transcriptomics enables spatially resolved gene expression measurements at near single-cell resolution. There is a pressing need for computational tools to enable the detection of genes that are differentially expressed (DE) within specific cell types across tissue context. We show that current approaches cannot learn cell type-specific DE due to changes in cell type composition across space and the fact that measurement units often detect transcripts from more than one cell type. Here, we introduce a statistical method, Cell type-Specific Inference of Differential Expression (C-SIDE), that identifies cell type-specific patterns of differential gene expression while accounting for localization of other cell types. We model spatial transcriptomics gene expression as an additive mixture across cell types of general log-linear cell type-specific expression functions. This approach provides a unified framework for defining and identifying gene expression changes in a wide-range of relevant contexts: changes due to pathology, anatomical regions, physical proximity to specific cell types, and cellular microenvironment. Furthermore, our approach enables statistical inference across multiple samples and replicates when such data is available. We demonstrate, through simulations and validation experiments on Slide-seq and MER-FISH datasets, that our approach accurately identifies cell type-specific differential gene expression and provides valid uncertainty quantification. Lastly, we apply our method to characterize spatially-localized tissue changes in the context of disease. In an Alzheimers mouse model Slide-seq dataset, we identify plaque-dependent patterns of cellular immune activity. We also find a putative interaction between tumor cells and myeloid immune cells in a Slide-seq tumor dataset. We make our C-SIDE method publicly available as part of the open source R package https://github.com/dmcable/spacexr. Computation of the mitochondrial age distribution along the axon length We describe a compartmental model of mitochondrial transport in axons, which we apply to compute mitochondrial age at different distances from the soma. The model predicts that the average mitochondrial age at the tip of an axon with a length of 1 cm is approximately 22 hours. The mitochondria are youngest closest to the soma and their age scales approximately linearly with distance from the soma. To the best of the authors knowledge, this is the first attempt to predict the spatial distribution of mitochondrial age within an axon. A sensitivity study of the mean age of mitochondria to various model parameters is also presented. SARS-CoV-2 infects lung epithelial cells and induces senescence and an inflammatory response in patients with severe COVID-19 RationaleSARS-CoV-2 infection of the respiratory system can progress to a life threatening multi-systemic disease, mediated via an excess of cytokines ("cytokine storm"), but the molecular mechanisms are poorly understood.

ObjectivesTo investigate whether SARS-CoV-2 may induce cellular senescence in lung epithelial cells, leading to secretion of inflammatory cytokines, known as the senescence-associated secretory phenotype (SASP).

MethodsAutopsy lung tissue samples from eleven COVID-19 patients and sixty age-matched non-infected controls were analysed by immunohistochemistry for SARS-CoV-2 and markers of cellular senescence (SenTraGor, p16INK4A) and key SASP cytokines (interleukin-1{beta}, interleukin-6). We also investigated whether SARS-CoV-2 infection of an epithelial cell line induces senescence and cytokine secretion.

Measurements and Main ResultsSARS-CoV-2 was detected by immunocytochemistry and electron microscopy predominantly in alveolar type-2 (AT2) cells, which also expressed the angiotensin-converting-enzyme 2 (ACE2), a critical entry receptor for this virus. In COVID-19 samples, AT2 cells displayed increased markers of senescence [p16INK4A, SenTraGor staining positivity in 12{+/-}1.2% of cells compared to 1.7{+/-}0.13% in non-infected controls (p<0.001)], with markedly increased expression of interleukin-1{beta} and interleukin-6 (p<0.001). Infection of epithelial cells (Vero E6) with SARS-CoV-2 in-vitro induced senescence and DNA damage (increased SenTraGor and {gamma}-H2AX), and reduced proliferation (Ki67) compared to uninfected control cells (p<0.01).

ConclusionsWe demonstrate that in severe COVID-19 patients, AT2 cells are infected with SARS-CoV-2 and show senescence and expression of proinflammatory cytokines. We also show that SARS-CoV-2 infection of epithelial cells may induce senescence and inflammation, indicating that cellular senescence may be an important molecular mechanism of severe COVID-19. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Joint interactions with humans may pose a higher risk of zoonotic outbreaks than interactions with conspecifics among wildlife populations at human-wildlife interfaces O_LIPandemics caused by wildlife-origin pathogens, like COVID-19, highlight the importance of understanding the ecology of zoonotic transmission and outbreaks among wildlife populations at human-wildlife interfaces. To-date, the relative effects of human-wildlife and wildlife-wildlife interactions on the likelihood of such outbreaks remain unclear.
C_LIO_LIIn this study, we used social network analysis and epidemiological Susceptible Infected Recovered (SIR) models, to track zoonotic outbreaks through wild animals joint propensities to engage in social-ecological co-interactions with humans, and their social grooming interactions with conspecifics.
C_LIO_LIWe collected behavioral and demographic data on 10 groups of macaques (Macaca spp.) living in (peri)urban environments across Asia. Outbreak sizes predicted by the SIR models were related to structural features of the social networks, and particular properties of individual animals connectivity within those networks.
C_LIO_LIOutbreak sizes were larger when the first-infected animal was highly central, in both types of networks. Across host-species, particularly for rhesus and bonnet macaques, the effects of network centrality on outbreak sizes were stronger through macaques human co-interaction networks compared to grooming networks.
C_LIO_LIOur findings, independent of pathogen-transmissibility, suggest that wildlife populations in the Anthropocene are vulnerable to zoonosis more so due to their propensities to aggregate around anthropogenic factors than their gregariousness with conspecifics. Thus, the costs of zoonotic outbreaks may outweigh the potential/perceived benefits of jointly interacting with humans to procure anthropogenic food. From One Health perspectives, animals that consistently interact with both humans and conspecifics across time and space are useful targets for disease spillover assessments and control.
C_LI