Reference genome-independent taxonomic profiling of microbiomes with mOTUs3 BackgroundTaxonomic profiling is a fundamental task in microbiome research that aims to detect and quantify the relative abundance of microorganisms in biological samples. Available methods using shotgun metagenomic data generally depend on the availability of sequenced and taxonomically annotated reference genomes. However, the majority of microorganisms have not been cultured yet and lack such reference genomes. Thus, a substantial fraction of microbial community members remains unaccounted for during taxonomic profiling of metagenomes, particularly in samples from underexplored environments. To address this issue, we have developed the mOTU profiler, a tool that enables reference genome-independent species-level profiling of metagenomes. As such, it supports the identification and quantification of both "known" and "unknown" species based on a set of select marker genes.

ResultsHere, we present mOTUs3, a command line tool that enables the profiling of metagenomes for >33,000 species-level operational taxonomic units. To achieve this, we leveraged the reconstruction and analysis of >600,000 draft genomes, most of which are metagenome assembled genomes (MAGs), from diverse microbiomes, including soil, freshwater systems, and the gastrointestinal tract of ruminants and other animals, which we found to be greatly underrepresented by reference genomes. Overall, two-thirds of all species-level taxa lacked a reference genome. The cumulative relative abundance of these newly included taxa was low in well-studied microbiomes, such as the human body sites (6-11%). By contrast, they accounted for substantial proportions (ocean, freshwater, soil: 43-63%) or even the vast majority (pig, fish, cattle: 60-80%) of the relative abundance across diverse non-human-associated microbiomes. Using community-developed benchmarks and datasets, we found mOTUs3 to be more accurate than other methods and to be more congruent with 16S rRNA gene-based methods for taxonomic profiling. Furthermore, we demonstrate that mOTUs3 greatly increases the resolution of well-known microbial groups into species-level taxa and helps identify new differentially abundant taxa in comparative metagenomic studies.

ConclusionsWe developed mOTUs3 to enable accurate species-level profiling of metagenomes. Compared to other methods, it provides a more comprehensive view of prokaryotic community diversity, in particular for currently underexplored microbiomes. To facilitate comparative analyses by the research community, it is released with >11,000 precomputed profiles for publicly available metagenomes and is freely available at: https://github.com/motu-tool/mOTUs.
Chromatin accessibility determines cell-type specific intron retention Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including cancer. Despite consistent reports demonstrating intrinsic sequence features that predispose introns to become retained, conflicting findings about cell type-specific IR regulation demand a systematic analysis in a controlled experimental setting. We integrated matched transcriptomics and epigenetics data (including DNA methylation, nucleosome occupancy, histone modifications) from primary human myeloid and lymphoid cells. Using machine learning we trained two complementary models to determine the role of epigenetic factors in the regulation of IR. We show that increased chromatin accessibility contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. With mounting reports linking pathogenic alterations to RNA processing, our findings may have profound implications for the design of therapeutic approaches targeting aberrant splicing.
In Mycobacterium abscessus, the stringent factor Rel regulates metabolism, but is not the only (p)ppGpp synthase The stringent response is a broadly conserved stress response system that exhibits functional variability across bacterial clades. Here, we characterize the role of the stringent factor Rel in the non-tuberculous mycobacterial pathogen, Mycobacterium abscessus (Mab). We found that deletion of rel does not ablate (p)ppGpp synthesis, and that rel does not provide a survival advantage in several stress conditions, or in antibiotic treatment. Transcriptional data show that RelMab is involved in regulating expression of anabolism and growth genes in stationary phase. However, it does not activate transcription of stress response or antibiotic resistance genes, and actually represses transcription of many antibiotic resistance genes. This work shows that there is an unannotated (p)ppGpp synthetase in Mab.

ImportanceIn this study, we examined the functional roles of the stringent factor Rel in Mycobacterium abscessus (Mab). In most species, stringent factors synthesize the alarmone (p)ppGpp, which globally alters transcription to promote growth arrest and survival under stress and in antibiotic treatment. Our work shows that in Mab, an emerging pathogen which is resistant to many antibiotics, the stringent factor Rel is not solely responsible for synthesizing (p)ppGpp. We find that RelMab downregulates many metabolic genes under stress, but does not upregulate stress response genes and does not promote antibiotic tolerance. This study implies that there is another critical but unannotated (p)ppGpp synthetase in Mab, and suggests that RelMab inhibitors are unlikely to sensitize Mab infections to antibiotic treatment. 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. Multiparametric profiling of cellular lipid metabolism in hypercholesterolemia Systematic insight into cellular dysfunctions can improve understanding of disease etiology, risk assessment and patient stratification. We present a multiparametric high-content imaging platform enabling quantification of low-density lipoprotein (LDL) uptake and lipid storage in cytoplasmic droplets of primary leukocyte subpopulations. We validated this platform with samples from 65 individuals with variable blood LDL-cholesterol (LDL-c) levels, including familial hypercholesterolemia (FH) and non-FH subjects. We integrated lipid storage data into a novel readout, lipid mobilization, measuring the efficiency with which cells deplete lipid reservoirs. Lipid mobilization correlated positively with LDL uptake and negatively with hypercholesterolemia and age, improving differentiation of individuals with normal and elevated LDL-c. Moreover, combination of cell-based readouts with a polygenic risk score for LDL-c explained hypercholesterolemia better than the genetic risk score alone. This platform provides functional insights into cellular lipid trafficking from a few mls of blood and is applicable to dissect metabolic disorders, such as hypercholesterolemia.

MotivationWe have limited information on how cellular lipid uptake and processing differ between individuals and influence the development of metabolic diseases, such as hypercholesterolemia. Available assays are labor intensive, require skilled personnel and are difficult to scale to higher throughput, making it challenging to obtain systematic functional cell-based data from individuals. To overcome this problem, we established a scalable automated analysis pipeline enabling reliable quantification of multiple cellular readouts, including lipid uptake, storage and mobilization, from different white blood cell populations. This approach provides new personalized insights into the cellular basis of hypercholesterolemia and obesity.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/440471v4_ufig1.gif" ALT="Figure 1">
View larger version (84K):
org.highwire.dtl.DTLVardef@9222c1org.highwire.dtl.DTLVardef@27e289org.highwire.dtl.DTLVardef@89ba82org.highwire.dtl.DTLVardef@33c9cd_HPS_FORMAT_FIGEXP  M_FIG C_FIG A Human-Immune-System (HIS) humanized mouse model (DRAGA: HLA-A2. HLA-DR4. Rag1 KO.IL-2Rγc KO. NOD) for COVID-19 We report the first Human Immune System (HIS)-humanized mouse model ("DRAGA": HLA-A2.HLA-DR4.Rag1KO.IL-2R{gamma}cKO.NOD) for COVID-19 research. This mouse is reconstituted with human cord blood-derived, HLA-matched hematopoietic stem cells. It engrafts human epi/endothelial cells expressing the human ACE2 receptor for SARS-CoV-2 and TMPRSS2 serine protease co-localized on lung epithelia. HIS-DRAGA mice sustained SARS-CoV-2 infection, showing deteriorated clinical condition, replicating virus in the lungs, and human-like lung immunopathology including T-cell infiltrates, microthrombi and pulmonary sequelae. Among T-cell infiltrates, lung-resident (CD103+) CD8+ T cells were sequestered in epithelial (CD326+) lung niches and secreted granzyme B and perforin, indicating cytotoxic potential. Infected mice also developed antibodies against the SARS-CoV-2 viral proteins. Hence, HIS-DRAGA mice showed unique advantages as a surrogate in vivo human model for studying SARS-CoV-2 immunopathology and for testing the safety and efficacy of candidate vaccines and therapeutics. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. A generic hidden Markov model for multi-parent populations A common step in the analysis of multi-parent populations is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of multi-parent populations. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly-available data on Diversity Outbred and Collaborative Cross mice Expression and subcellular localization of USH1C/harmonin in the human retina provides insights into pathophysiology and therapy of Usher syndrome Usher syndrome (USH) is the most common form of hereditary deafness-blindness in humans. USH is a complex genetic disorder, assigned to three clinical subtypes differing in onset, course and severity, with USH1 being the most severe. Rodent USH1 models do not reflect the ocular phenotype observed in human patients to date; hence, little is known about the pathophysiology of USH1 in the human eye. One of the USH1 genes, USH1C, exhibits extensive alternative splicing and encodes numerous harmonin protein isoforms that function as scaffolds for organizing the USH interactome. RNA-seq analysis of human retina uncovered harmonin a1 as the most abundant transcript of USH1C. Bulk RNA-seq analysis and immunoblotting showed abundant expression of harmonin in Muller glia cells (MGCs) and retinal neurons. Furthermore, harmonin was localized in the terminal endfeet and apical microvilli of MGCs, presynaptic region (pedicle) of cones, and outer segments of rods as well as at adhesive junctions of MGCs and photoreceptors in the outer limiting membrane (OLM). Our data provide evidence for the interaction of harmonin with OLM molecules in MGCs and photoreceptors and rhodopsin in photoreceptors. Subcellular expression and colocalization of harmonin correlates with the clinical phenotype observed in USH1C patients. In addition, primary cilia defects in USH1C patient-derived fibroblasts could be reverted by delivery of harmonin a1 transcript isoform. Our data provide novel insights into photoreceptor cell biology, USH1C pathophysiology and for developing gene therapy treatment. 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. Cell-type specific circadian bioluminescence rhythms in Dbp reporter mice Circadian rhythms are endogenously generated physiological and molecular rhythms with a cycle length of about 24 h. Bioluminescent reporters have been exceptionally useful for studying circadian rhythms in numerous species. Here, we report development of a reporter mouse generated by modification of a widely expressed and highly rhythmic gene encoding D-site albumin promoter binding protein (Dbp). In this line of mice, firefly luciferase is expressed from the Dbp locus in a Cre-recombinase-dependent manner, allowing assessment of bioluminescence rhythms in specific cellular populations. A mouse line in which luciferase expression was Cre-independent was also generated. The Dbp reporter alleles do not alter Dbp gene expression rhythms in liver or circadian locomotor activity rhythms. In vivo and ex vivo studies show the utility of the reporter alleles for monitoring rhythmicity. Our studies reveal cell-type specific characteristics of rhythms among neuronal populations within the suprachiasmatic nuclei ex vivo. In vivo studies show Dbp-driven bioluminescence rhythms in the liver of Albumin-Cre;DbpKI/+ "liver reporter" mice. After a shift of the lighting schedule, locomotor activity achieved the proper phase relationship with the new lighting cycle more rapidly than hepatic bioluminescence did. As previously shown, restricting food access to the daytime altered the phase of hepatic rhythmicity. Our model allowed assessment of the rate of recovery from misalignment once animals were provided with food ad libitum. These studies confirm the previously demonstrated circadian misalignment following environmental perturbations and reveal the utility of this model for minimally invasive, longitudinal monitoring of rhythmicity from specific mouse tissues. An automated workflow for multiplexed single-cell proteomics sample preparation at unprecedented sensitivity. The analysis of single-cell proteomes has recently become a viable complement to transcriptomics and genomics studies. Proteins are the primary driver of cellular functionality and mRNA levels are often an unreliable proxy of such. Therefore, the global analysis of the proteome is essential to study cellular identities. Multiplexed and label-free mass spectrometry-based approaches with single-cell resolution have lately attributed surprising heterogeneity to presumed homogenous cell populations. Even though specialized experimental designs and instrumentation have demonstrated remarkable advances, the efficient sample preparation of single cells still lag. Here, we introduce the proteoCHIP, a universal option for single-cell proteomics sample preparation at surprising sensitivity and throughput. The automated processing using a commercial system combining single-cell isolation and picoliter dispensing, the cellenONE(R), reduces final sample volumes to low nanoliters submerged in a hexadecane layer simultaneously eliminating error-prone manual sample handling and overcoming evaporation. The specialized proteoCHIP design allows direct injection of single cells via a standard autosampler resulting in around 1,500 protein groups per analytical run at remarkable reporter ion signal to noise while reducing or eliminating the carrier proteome. We identified close to 2,600 proteins across 170 multiplexed single cells from two highly similar human cell types. This dedicated loss-less workflow allows distinguishing in vitro co-differentiated cell types of self-organizing cardiac organoids based on indicative markers across 150 single cells. In-depth characterization revealed enhanced cellular motility of cardiac endothelial cells and sarcomere organization in cardiomyocytes. Our versatile and automated sample preparation has not only proven to be easily adaptable but is also sufficiently sensitive to drive biological applications of single-cell proteomics. 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. Diurnal activity budget and conservation implication of Hippopotamus amphibius at Bui National Park, Ghana Animals apportion time for their daily behavioral activities. We studied the activity budget of Hippopotamus amphibius at the Black Volta River in the Bui National Park, Ghana. We performed instantaneous scan sample using ground count survey between August 2020 to July 2021. We observed that H. amphibius spent most of their day time resting (54.75 %), followed by feeding (22.93 %), walking (19.2 %) and touching (3.12 %). We found significant difference in the time spent between the daily activities (H = 41.67; p < 0.0001). Except activities involving touching, we found no significant seasonal differences in feeding, resting and walking activities by H. amphibius. Understanding the behavioral activities of H. amphibius will enhance management and conservation of the animals. Life history adaptations to fluctuating environments: Combined effects of demographic buffering and lability of demographic parameters Demographic buffering and lability have been identified as adaptive strategies to optimise fitness in a fluctuating environment. These are not mutually exclusive, however we lack efficient methods to measure their relative importance for a given life history. Here, we decompose the stochastic growth rate (fitness) into components arising from nonlinear responses and variance-covariance of demographic parameters to an environmental driver, which allows studying joint effects of buffering and lability. We apply this decomposition for 154 animal matrix population models under different scenarios, to explore how these main fitness components vary across life histories. Faster-living species appear more responsive to environmental fluctuations, either positively or negatively. They have the highest potential for strong adaptive demographic lability, while demographic buffering is a main strategy in slow-living species. Our decomposition provides a comprehensive framework to study how organisms adapt to variability through buffering and lability, and to predict species responses to climate change. 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 Fatal neuroinvasion and SARS-CoV-2 tropism in K18-hACE2 mice is partially independent on hACE2 expression Animal models recapitulating distinctive features of severe COVID-19 are critical to enhance our understanding of SARS-CoV-2 pathogenesis. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. The precise mechanisms of lethality in this mouse model remain unclear. Here, we evaluated the spatiotemporal dynamics of SARS-CoV-2 infection for up to 14 days post-infection. Despite infection and moderate pneumonia, rapid clinical decline or death of mice was invariably associated with viral neuroinvasion and direct neuronal injury (including brain and spinal neurons). Neuroinvasion was observed as early as 4 dpi, with virus initially restricted to the olfactory bulb supporting axonal transport via the olfactory neuroepithelium as the earliest portal of entry. No evidence of viremia was detected suggesting neuroinvasion occurs independently of entry across the blood brain barrier. SARS-CoV-2 tropism was not restricted to ACE2-expressing cells (e.g., AT1 pneumocytes), and some ACE2-positive lineages were not associated with the presence of viral antigen (e.g., bronchiolar epithelium and brain capillaries). Detectable ACE2 expression was not observed in neurons, supporting overexpression of ACE2 in the nasal passages and neuroepithelium as more likely determinants of neuroinvasion in the K18-hACE2 model. Although our work incites caution in the utility of the K18-hACE2 model to study global aspects of SARS-CoV-2 pathogenesis, it underscores this model as a unique platform for exploring the mechanisms of SARS-CoV-2 neuropathogenesis that may have clinical relevance acknowledging the growing body of evidence that suggests COVID-19 may result in long-standing neurologic consequences.

IMPORTANCECOVID-19 is predominantly a respiratory disease caused by SARS-CoV-2 that has infected more than 191 million people with over 4 million fatalities (2021-07-20). The development of animal models recapitulating distinctive features of severe COVID-19 is critical to enhancing our understanding of SARS-CoV-2 pathogenesis and in the evaluation of vaccine and therapeutic efficacy. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. Here, we show lethality of this model is invariably associated with viral neuroinvasion linked with viral replication and assembly. Importantly, pneumonia albeit invariably present was generally moderate with the absence of culturable infectious virus at peak neuroinvasion. The dynamics of viral neuroinvasion and pneumonia were only partially dependent on hACE2. Overall, this study provides an in-depth sequential characterization of the K18-hACE2 model following SARS-CoV-2 infection, highlighting its significance to further study the mechanisms of SARS-CoV-2 neuropathogenesis. 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. 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. Analysis of diverse eukaryotes suggests the existence of an ancestral 1 mitochondrial apparatus derived from the bacterial type II secretion system The type 2 secretion system (T2SS) is present in some Gram-negative eubacteria and used to secrete proteins across the outer membrane. Here we report that certain representative heteroloboseans, jakobids, malawimonads and hemimastigotes unexpectedly possess homologues of core T2SS components. We show that at least some of them are present in mitochondria, and their behaviour in biochemical assays is consistent with the presence of a mitochondrial T2SS-derived system (miT2SS). We additionally identified 23 protein families co-occurring with miT2SS in eukaryotes. Seven of these proteins could be directly linked to the core miT2SS by functional data and/or sequence features, whereas others may represent different parts of a broader functional pathway, possibly also involving the peroxisome. Its distribution in eukaryotes and phylogenetic evidence together indicate that the miT2SS-centred pathway is an ancestral eukaryotic trait. Our findings thus have direct implications for the functional properties of the early mitochondrion. 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. Disrupting short-term memory in premotor cortex affects serial dependence in visuomotor integration Human behavior is biased by past experience. For example, when intercepting a moving target, the speed of previous targets will bias responses in future trials. Neural mechanisms underlying this so-called serial dependence are still under debate. Here, we tested the hypothesis that the previous trial leaves a neural trace in brain regions associated with encoding task-relevant information in visual and/or motor regions. We reasoned that injecting noise by means of transcranial magnetic stimulation (TMS) over premotor and visual areas would degrade such memory traces and hence reduce serial dependence. To test this hypothesis, we applied bursts of TMS pulses to right visual motion processing region hV5/MT+ and to left dorsal premotor cortex during inter-trial intervals of a coincident timing task performed by twenty healthy human participants (15 female). Without TMS, participants presented a bias towards the speed of the previous trial when intercepting moving targets. TMS over dorsal premotor cortex decreased serial dependence in comparison to the control Vertex stimulation, whereas TMS applied over hV5/MT+ did not. In addition, TMS seems to have specifically affected the memory trace that leads to serial dependence, as we found no evidence that participants behavior worsened after applying TMS. These results provide causal evidence that an implicit short-term memory mechanism in premotor cortex keeps information from one trial to the next, and that this information is blended with current trial information so that it biases behavior in a visuomotor integration task with moving objects.

Significance StatementHuman perception and action are biased by the recent past. The origin of such serial bias is still not fully understood, but a few components seem to be fundamental for its emergence: the brain needs to keep previous trial information in short-term memory and blend it with incoming information. Here, we present evidence that a premotor area has a potential role in storing previous trial information in short-term memory in a visuomotor task, and that this information is responsible for biasing ongoing behavior. These results corroborate the perspective that areas associated with processing information of a stimulus or task also participate in maintaining that information in short-term memory even when this information is no longer relevant for current behavior. Quantifying ERK-activity in response to inhibition of the BRAFV600E-MEK-ERK cascade using mathematical modelling BackgroundSimultaneous inhibition of multiple components of the BRAF-MEK-ERK cascade (vertical inhibition) has become a standard of care for treating BRAF-mutant melanoma. However, the molecular mechanism of how vertical inhibition synergistically suppresses intracellular ERK-activity, and consequently cell proliferation, are yet to be fully elucidated.

MethodsWe develop a mechanistic mathematical model that describes how the mutant BRAF-inhibitor, dabrafenib, and the MEK-inhibitor, trametinib, affect BRAFV600E-MEK-ERK signalling. The model is based upon a system of chemical reactions that describes cascade signalling dynamics. Using mass action kinetics, the chemical reactions are re-expressed as ordinary differential equations which are parameterised by in vitro data and solved numerically to obtain the temporal evolution of cascade component concentrations.

ResultsThe model provides a quantitative method to compute how dabrafenib and trametinib can be used in combination to synergistically inhibit ERK-activity in BRAFV600E-mutant melanoma cells. The model elucidates molecular mechanisms of vertical inhibition of the BRAFV600E-MEK-ERK cascade and delineates how elevated BRAF concentrations generate drug resistance to dabrafenib and trametinib. The computational simulations further suggest that elevated ATP levels could be a factor in drug resistance to dabrafenib.

ConclusionsThe model can be used to systematically motivate which dabrafenib-trametinib dose-combinations, for treating BRAFV600E-mutated melanoma, warrant experimental investigation. A humanized nanobody phage display library yields potent binders of SARS CoV-2 spike Neutralizing antibodies targeting the SARS-CoV-2 spike protein have shown a great preventative/therapeutic potential. Here, we report a rapid and efficient strategy for the development and design of SARS-CoV-2 neutralizing humanized nanobody constructs with sub-nanomolar affinities and nanomolar potencies. CryoEM-based structural analysis of the nanobodies in complex with spike revealed two distinct binding modes. The most potent nanobody, RBD-1-2G(NCATS-BL8125), tolerates the N501Y RBD mutation and remains capable of neutralizing the B.1.1.7 (Alpha) variant. Molecular dynamics simulations provide a structural basis for understanding the neutralization process of nanobodies exclusively focused on the spike-ACE2 interface with and without the N501Y mutation on RBD. A primary human airway air-lung interface (ALI) ex vivo model showed that RBD-1-2G-Fc antibody treatment was effective at reducing viral burden following WA1 and B.1.1.7 SARS-CoV-2 infections. Therefore, this presented strategy will serve as a tool to mitigate the threat of emerging SARS-CoV-2 variants. Text has been revised, image definition has been revised. Maize is a staple food all over the world. Models for climate change suggest that, in the future, cloud formation might be reduced in the tropics increasing the exposure to Ultraviolet-B (UV-B) radiation, a DNA-damaging agent. UV-B (290 to 320 nm) has been shown to affect yield in maize. In this project we have determined the differences in DNA repair efficiencies between U.S. inbred lines B73 and Mo17, and Central American purple landraces from Guatemala and Costa Rica. Our results from single cell electrophoresis experiments (Comet Assay) suggest that the landrace Pujagua Santa Cruz (P1, Costa Rican) was resistant to damage caused by the radiomimetic agent zeocin (24h/100 micrograms per mL), while landrace Pujagua La Cruz (P2, Costa Rican) was able to repair DNA damage after one hour. On the other hand, line Mo17 (Missouri, USA) was unable to repair the damage, while B73 (Iowa, USA) and the landraces Jocopilas (Guatemalan), Orotina Congo (Costa Rican) and Talamanca (Costa Rican) were partially able to repair DNA damage. High Resolution Melting (HRM) curve analysis of putative homologous DNA repair gene ZeaATM1 showed that both P1 and P2 had differences in the melting temperatures for this gene compared to B73 and Mo17, while P1 showed additional differences in ZeaSOG1, ZeaRAD51 and ZeaBRCA1, suggesting that in this landrace the presence of polymorphisms may be common among key genes for this pathway. Taken together our results suggest that key adaptive differences in DNA repair efficiencies exist between inbred lines and landraces of maize and that some Central America landraces could be used as a valuable pool of alleles for plant breeding aiming to increase tolerance to radiation. A coarse-grained NADH redox model enables inference of subcellular metabolic fluxes from fluorescence lifetime imaging Mitochondrial metabolism is of central importance to diverse aspects of cell and developmental biology. Defects in mitochondria are associated with many diseases, including cancer, neuropathology, and infertility. Our understanding of mitochondrial metabolism in situ and dysfunction in diseases are limited by the lack of techniques to measure mitochondrial metabolic fluxes with sufficient spatiotemporal resolution. Herein, we developed a new method to infer mitochondrial metabolic fluxes in living cells with subcellular resolution from fluorescence lifetime imaging of NADH. This result is based on the use of a generic coarse-grained NADH redox model. We tested the model in mouse oocytes and human tissue culture cells subject to a wide variety of perturbations by comparing predicted fluxes through the electron transport chain (ETC) to direct measurements of oxygen consumption rate. Interpreting the FLIM measurements of NADH using this model, we discovered a homeostasis of ETC flux in mouse oocytes: perturbations of nutrient supply and energy demand do not change ETC flux despite significantly impacting NADH metabolic state. Furthermore, we observed a subcellular spatial gradient of ETC flux in mouse oocytes and found that this gradient is primarily a result of heterogeneous mitochondrial proton leak. We concluded from these observations that ETC flux in mouse oocytes is not controlled by energy demand or supply, but by the intrinsic rates of mitochondrial respiration.