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.
Cluster-specific gene markers enhance Shigella and Enteroinvasive Escherichia coli in silico serotyping Shigella and enteroinvasive Escherichia coli (EIEC) cause human bacillary dysentery with similar invasion mechanisms and share similar physiological, biochemical and genetic characteristics. The ability to differentiate Shigella and EIEC from each other is important for clinical diagnostic and epidemiologic investigations. The existing genetic signatures may not discriminate between Shigella and EIEC. However, phylogenetically, Shigella and EIEC strains are composed of multiple clusters and are different forms of E. coli. In this study, we identified 10 Shigella clusters, 7 EIEC clusters and 53 sporadic types of EIEC by examining over 17,000 publicly available Shigella/EIEC genomes. We compared Shigella and EIEC accessory genomes to identify the cluster-specific gene markers or marker sets for the 17 clusters and 53 sporadic types. The gene markers showed 99.63% accuracy and more than 97.02% specificity. In addition, we developed a freely available in silico serotyping pipeline named Shigella EIEC Cluster Enhanced Serotype Finder (ShigEiFinder) by incorporating the cluster-specific gene markers and established Shigella/EIEC serotype specific O antigen genes and modification genes into typing. ShigEiFinder can process either paired end Illumina sequencing reads or assembled genomes and almost perfectly differentiated Shigella from EIEC with 99.70% and 99.81% cluster assignment accuracy for the assembled genomes and mapped reads respectively. ShigEiFinder was able to serotype over 59 Shigella serotypes and 22 EIEC serotypes and provided a high specificity with 99.40% for assembled genomes and 99.38% for mapped reads for serotyping. The cluster markers and our new serotyping tool, ShigEiFinder (https://github.com/LanLab/ShigEiFinder), will be useful for epidemiologic and diagnostic investigations.

Impact statementThe differentiation of Shigella strains from enteroinvasive E. coli (EIEC) is important for clinical diagnosis and public health epidemiologic investigations. The similarities between Shigella and EIEC strains make this differentiation very difficult as both share common ancestries within E. coli. However, Shigella and EIEC are phylogenetically separated into multiple clusters, making high resolution separation using cluster specific genomic markers possible. In this study, we identified 17 Shigella or EIEC clusters including five that were newly identified through examination of over 17,000 publicly available Shigella and EIEC genomes. We further identified an individual or a set of cluster-specific gene markers for each cluster using comparative genomic analysis. These markers can then be used to classify isolates into clusters and were used to develop an in silico pipeline, ShigEiFinder (https://github.com/LanLab/ShigEiFinder) for accurate differentiation, cluster typing and serotyping of Shigella and EIEC from Illumina sequencing reads or assembled genomes. This study will have broad application from understanding the evolution of Shigella/EIEC to diagnosis and epidemiology.

Data summarySequencing data have been deposited at the National Center for Biotechnology Information under BioProject number PRJNA692536.

RepositoriesRaw sequence data are available from NCBI under the BioProject number PRJNA692536.
High-throughput submicron-resolution microscopyof entire C. elegans populations under strongimmobilization by cooling cultivation plates Despite its profound impact on biology, most high-resolution in vivo microscopy approaches remain low throughput because current immobilization techniques require significant manual effort. We greatly accelerate imaging of the nematode Caenorhabditis elegans by implementing a simple cooling approach to easily immobilize entire populations directly on their cultivation plates. We optimize and characterize cooling immobilization. Counterintuitively, relatively warmer temperatures immobilize animals significantly more effectively than colder temperatures utilized in prior studies. This enhanced immobilization enables clear submicron-resolution fluorescence imaging, which is challenging to achieve with most current immobilization techniques. We demonstrate 64x magnification 3D imaging and timelapse recording of neurons in adults and embryos without motion blur. Compared to standard azide immobilization, cooling immobilization reduces the animal preparation and recovery time by >98%, significantly increasing experimental speed. By obviating individual animal manipulation, our approach could also empower automated imaging of large C. elegans populations within standard experimental setups and workflows. A terrain treadmill to study animal locomotion through large obstacles A major challenge to understanding locomotion in complex 3-D terrain with large obstacles is to create tools for controlled, systematic lab experiments. Existing terrain arenas only allow observations at small spatiotemporal scales ([~]10 body length, [~]10 stride cycles). Here, we create a terrain treadmill to enable high-resolution observations of animal locomotion through large obstacles over large spatiotemporal scales. An animal moves through modular obstacles on an inner sphere, while a rigidly-attached, concentric, transparent outer sphere rotated with the opposite velocity via closed-loop feedback to keep the animal on top. During sustained locomotion, a discoid cockroach moved through pillar obstacles for 25 minutes ({approx} 2500 strides) over 67 m ({approx} 1500 body lengths), and was contained within a radius of 4 cm (0.9 body length) for 83% of the duration, even at speeds of up to 10 body length/s. The treadmill enabled observation of diverse locomotor behaviors and quantification of animal-obstacle interaction.

Summary StatementA device keeps a small animal on top of a sphere while it traverses obstacles, creating a "terrain treadmill" to study locomotion over a long time and distance at high-resolution. The Heterochromatin protein 1 is a master regulator in RNA splicing precision deficient in ulcerative colitis Defects in RNA splicing have been linked to numerous human disorders, but remain poorly explored in inflammatory bowel disease (IBD). Here, we report that, in the gut epithelium of patients with ulcerative colitis (UC), the expression of the chromatin and alternative splicing regulator HP1{gamma} is strongly reduced. Accordingly, inactivation of the HP1{gamma} gene in the mouse gut triggered several IBD-like traits, including inflammation and dysbiosis. In parallel, we discovered that its loss of function broadly increased splicing noise, reducing requirement for canonical splicing consensus sequences, and favoring the usage of cryptic splice sites at numerous genes with key functions in gut biology. This notably resulted in the production of progerin, a noncanonical toxic splice variant of prelamin A mRNA, responsible for the Hutchinson Gilford Progeria Syndrome (HGPS) of premature aging. Likewise, production of progerin transcript was found to be a signature of colonic cells from UC patients. Thus, our study identifies HP1{gamma} as a regulator of RNA metabolism in vivo, providing a unique mechanism linking anti-inflammation and accuracy of RNA splicing in the gut epithelium. HP1 defect may confer a general disturbance in RNA splicing precision to scrutinize in IBD and more generally in accelerating aging diseases. Intestine-enriched apolipoprotein b orthologs are required for stem cell differentiation and regeneration in planarians Lipid metabolism plays an instructive role in regulating stem cell state and differentiation. However, the roles of lipid mobilization and utilization in stem cell-driven regeneration are unclear. Planarian flatworms readily restore missing tissue due to injury-induced activation of pluripotent somatic stem cells called neoblasts. Here, we identify two intestine-enriched orthologs of apolipoprotein b, apob-1 and apob-2, which mediate transport of neutral lipid stores from the intestine to target tissues including neoblasts, and are required for tissue homeostasis and regeneration. Inhibition of apob function by RNAi causes head regression and lysis in uninjured animals, and delays body axis re-establishment and regeneration of multiple organs in amputated fragments. Furthermore, apob RNAi causes expansion of the population of differentiating neoblast progeny and dysregulates expression of genes enriched in differentiating and mature cells in eight major cell type lineages. We conclude that intestine-derived lipids serve as a source of metabolites required for neoblast differentiation. Pseudomonas aeruginosa C-terminal processing protease CtpA assembles into a hexameric structure that requires activation by a spiral-shaped lipoprotein binding partner Pseudomonas aeruginosa CtpA is a carboxyl terminal-processing protease that partners with the outer membrane lipoprotein LbcA to degrade at least five cell wall-associated proteins, four of which are cell wall hydrolases. This activity plays an important role in supporting P. aeruginosa virulence in a mouse model of acute pneumonia. However, almost nothing is known about the molecular mechanisms underlying CtpA and LbcA function. Here, we used structural analysis to show that CtpA alone assembles into an inactive hexamer comprising a trimer of dimers, which limits its substrate access and prevents nonspecific degradation. The adaptor protein LbcA is a right-handed open spiral with 11 tetratricopeptide repeats, which might wrap around a substrate to deliver it to CtpA for degradation. By structure-guided mutagenesis and functional assays, we also showed that the interfaces of the CtpA trimer-of-dimers, and an N-terminal helix of LbcA, are important for LbcA-mediated substrate degradation by CtpA both in vitro and in vivo. This work improves our understanding of the molecular mechanism of a CTP within the C-terminal processing peptidase-3 group.

IMPORTANCECarboxyl-terminal processing proteases (CTPs) are found in all three domains of life. In bacteria, some CTPs have been associated with virulence, raising the possibility that they could be theraputic targets. However, relatively little is known about their molecular mechanisms of action. In Pseudomonas aeruginosa, CtpA supports virulence by working in complex with the outer membrane lipoprotein LbcA to degrade cell wall hydrolases. Here, we report structure-function analyses of CtpA and LbcA, which reveals that CtpA assembles into an inactive hexamer comprising a trimer of dimers. LbcA is monomeric, with an N-terminal region important for binding to and activating CtpA, followed by a spiral structure composed of 11 tetratricopetide repeats, which could wrap around a substrate for delivery to CtpA. This work provides the first structure of a CTP-3 group member, revealing a unique mutimeric arrangement and insight into how this important proteolytic system functions. Modeling microRNA-driven post-transcriptional regulation by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Cell-specific Bioorthogonal Tagging of Glycoproteins Altered glycosylation is an undisputed corollary of cancer development. Understanding these alterations is paramount but hampered by limitations underlying cellular model systems. For instance, the intricate interactions between tumour and host cannot be adequately recapitulated in monoculture of tumour-derived cell lines. More complex co-culture models usually rely on sorting procedures for proteome analyses and rarely capture the details of protein glycosylation. Here, we report a strategy termed Bio-Orthogonal Cell line-specific Tagging of Glycoproteins (BOCTAG). Cells are equipped by transfection with an artificial biosynthetic pathway that transforms bioorthogonally tagged sugars into the corresponding nucleotide-sugars. Only transfected cells incorporate bioorthogonal tags into glycoproteins in the presence of non-transfected cells. We employ BOCTAG as an imaging technique and to annotate cell-specific glycosylation sites in mass spectrometry-glycoproteomics. We demonstrate application in co-culture and mouse models, allowing for profiling of the glycoproteome as an important modulator of cellular function. 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. Propionate ameliorates diabetes-induced neurological dysfunction through regulating the PI3K/Akt/eNOS signaling pathway A large body of research has established diabetes-related cognitive deterioration, sometimes known as "diabetic encephalopathy." Current evidence supports that oxidative stress, neuronal apoptosis, and cerebral microcirculation weakness are associated with cognition deficits induced by diabetes. The present study explores the effect of propionate on neurological deficits, cerebral blood flow, and oxidative stress in diabetic mice. Propionate can markedly improve neurological function, which is correlated with its capabilities of stimulating nitric oxide (NO) production, increasing cerebral microcirculation, suppressing oxidative stress and reducing neuron loss in the hippocampus. In addition, the results of Western Blotting indicated that the brain-protective function of propionate in streptozocin (STZ)-induced insulin-dependent diabetes mellitus (IDDM) mice is related to phosphoinositide 3-kinase (PI3K)/serine-threonine protein kinase (Akt)/endothelial nitric oxide synthase(eNOS) signaling pathway. In a diabetic mouse model, propionate reduces cerebral microcirculation, hippocampus apoptosis, and neurological impairment. Thus, propionate, now employed as a food preservative, may also help slow diabetes-induced cognitive loss.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=110 SRC="FIGDIR/small/464493v2_ufig1.gif" ALT="Figure 1">
View larger version (27K):
org.highwire.dtl.DTLVardef@81fa09org.highwire.dtl.DTLVardef@1469451org.highwire.dtl.DTLVardef@b37fa1org.highwire.dtl.DTLVardef@8037e7_HPS_FORMAT_FIGEXP  M_FIG C_FIG HighlightsO_LIPropionate ameliorates diabetic-induced depression-like behavior, spatial learning and memory deficits.
C_LIO_LIPropionate increases cerebral microcirculation by enhancing NO production in diabetic mice.
C_LIO_LIPropionate suppresses oxidative stress and reduces neuron loss in the hippocampus.
C_LIO_LIThe brain-protective function of propionate is related to PI3K/Akt/eNOS signaling pathway.
C_LI Define protein variant functions with high-complexity mutagenesis libraries and enhanced mutation detection software ASMv1.0 Open reading frame (ORF) variant libraries have advanced our ability to query the functions of a large number of variants of a protein simultaneously in a single experiment. Variant libraries targeting full-length ORFs typically consists of all possible single-amino-acid substitutions and a stop codon at each amino-acid position. Because a variant differs from the template ORF by a single codon variation, variant quantification presents the most profound challenge to this technology. Efforts such as dividing a library into sub-libraries for direct sequencing, or tag-directed subassembly are practical only for small ORFs. Our approach, however, features an enhanced variant-calling tool that processes massive shotgun sequencing data, allowing one to screen a single variant library for ORFs sized up to [~]3600 bases. The variant-calling tool detects variants reliably, and also presents variant annotations in datafiles enabling analyses that have reshaped our strategies governing library design, screen deconvolution, sequencing and sequencing data analysis. 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. Sarcoplasmic reticulum calcium handling in young DBA/2J vs C57BL/10 mdx models of Duchenne muscular dystrophy The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. Identification of Antimalarial Compounds that Inhibit Apicomplexan AP2 Proteins in the Human Malaria Parasite Plasmodium falciparum Plasmodium parasites are reliant on the Apicomplexan AP2 (ApiAP2) transcription factor family to regulate gene expression programs. AP2 DNA binding domains have no homologs in the human or mosquito host genomes, making them potential antimalarial drug targets. Using an in-silico screen to dock thousands of small molecules into the crystal structure of the AP2-EXP (Pf3D7_1466400) AP2 domain (PDB:3IGM), we identified compounds that interact with this domain. Four compounds were found to compete for DNA binding with AP2-EXP and at least one additional ApiAP2 protein. Our top ApiAP2 competitor compound perturbs the transcriptome of P. falciparum trophozoites and results in a decrease in abundance of log2 fold change > 2 for 50% (46/93) of AP2-EXP target genes. Additionally, two ApiAP2 competitor compounds have multi-stage anti-Plasmodium activity against blood and mosquito stage parasites. In summary, we describe a novel set of antimalarial compounds that are targeted against the ApiAP2 family of proteins. These compounds may be used for future chemical genetic interrogation of ApiAP2 proteins or serve as starting points for a new class of antimalarial therapeutics.

Author SummaryPlasmodium parasites are the causative agent of malaria, which resulted in over 600,000 deaths in 2021. Due to resistance arising for every antimalarial therapeutic deployed to date, new drug targets and druggable pathways must be explored. To address this concern, we used a molecular docking screen to predict competitors of DNA binding by the parasite specific family of Apicomplexan AP2 (ApiAP2) transcription factor proteins for testing in vitro and in vivo. We find that ApiAP2 competing compounds have antimalarial activity consistent with the disruption of gene regulation. This work will further our understanding of both the biological role and targetability of parasite transcriptional regulation. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Estimating the fundamental niche: accounting for the uneven availability of existing environmental conditions in the calibration area Studies that question important conceptual and methodological aspects of the field of ecological niche modeling (and species distribution modeling) have cast doubts on whether it is possible to estimate the fundamental niche of a species using presence-only data. The main limitation in niche estimation is that presence data come from the realized niche, which is only a subset of the fundamental niche. Most of the existing methods lack the ability to overcome this limitation and therefore fit niches that resemble the realized niche. To obtain a more accurate estimate of the fundamental niche, we propose using the geographic region that is accessible to a species (based on its dispersal ability) to determine a sampling distribution, in environmental space, from which we can quantify the likelihood of observing a particular environmental combination in a sample of presence points. We incorporate this sampling distribution into a multivariate normal model (i.e., a Mahalanobis distance model) by creating a weight function that modifies the probabilities of observing different environmental combinations in a sample of presences. This modification accounts for the uneven availability of environmental conditions. We show that the parameters of this weighted-normal model can be approximated with a maximum likelihood estimation approach; and then used to draw ellipsoids (confidence regions) that represent the fundamental niche of the species. We illustrate the application of our model with two worked examples. First, we use presence data for an invasive species and an accessible area that includes only its native range to evaluate whether the fitted model predicts confirmed establishments of the species outside its native range. Second, we use presence data for closely related species with known accessible areas to demonstrate how the different dispersal abilities of the species constrain a classic Mahalanobis distance model. Overall, we show that accounting for the distribution of environmental conditions that are accessible to a species indeed affects the estimation of the ellipsoids used to model its fundamental niche. A Machine Learning Approach to Identify Small Molecule Inhibitors of Secondary Nucleation inalpha-Synuclein Aggregation Drug development is an increasingly active area of application of machine learning methods, due to the need to overcome the high attrition rates of conventional drug discovery pipelines. This issue is especially pressing for neurodegenerative diseases where very few disease-modifying drugs have been approved. To address this problem, we describe a machine learning approach to identify specific inhibitors of the proliferation of -synuclein aggregates through secondary nucleation, a process that has been implicated in Parkinsons disease and related synucleinopathies. We use a combination of docking simulations followed by machine learning to first identify initial hit compounds and then explore the chemical space around these compounds. Our results demonstrate that this approach leads to the identification of novel chemical matter with an improved hit rate and potency over conventional similarity search approaches. Double trouble: having two hosts reduces infection prevalence in vectored trypanosomatids Trypanosomatids are a diverse family of protozoan parasites, some of which cause devastating human and livestock diseases. There are two distinct infection life-cycles in trypanosomatids; some species complete their entire life-cycle in a single host (monoxenous) while others infect two hosts (dixenous). Dixenous trypanosomatids are mostly vectored by insects, and the human trypanosomatid diseases are caused mainly by vectored parasites. While infection prevalence has been described for subsets of hosts and trypanosomatids, little is known about whether monoxenous and dixenous trypanosomatids differ in infection prevalence. Here, we use meta-analyses to synthesise all published evidence of trypanosomatid infection prevalence for the last two decades, encompassing 931 unique host-trypansomatid systems. In examining 581 studies that describe infection prevalence, we find, strikingly, that monoxenous species are two-fold more prevalent than dixenous species across all hosts. We also find that dixenous trypanosomatids have significantly lower infection prevalence in insects than their non-insect hosts. Within monoxenous trypanosomatids, genera infecting bees are most prevalent and infection prevalence does not vary between wild and managed bees. To our knowledge, these results reveal for the first time, a fundamental difference in infection prevalence according to host specificity where vectored species suffer from lower infection prevalence as a result of a  jack of all trades, master of none style trade-off between the vector and subsequent hosts. Intramolecular carbon isotope signals reflect metabolite allocation in plants Stable isotopes at natural abundance are key tools to study physiological processes occurring outside the temporal scope of manipulation and monitoring experiments. Whole-molecule carbon isotope ratios (13C/12C) enable assessments of plant carbon uptake yet conceal information about carbon allocation. Here, we identify an intramolecular 13C/12C signal at treering glucose C-5 and C-6 and develop experimentally testable theories on its origin. More specifically, we assess the potential of processes within C3 metabolism for signal introduction based (inter alia) on constraints on signal propagation posed by metabolic networks. We propose that the intramolecular signal reports carbon allocation into major metabolic pathways in actively photosynthesising leaf cells including the anaplerotic, shikimate, and non-mevalonate pathway. We support our theoretical framework by linking it to previously reported whole-molecule 13C/12C increases in cellulose of ozone-treated Betula pendula and a highly significant relationship between the intramolecular signal and tropospheric ozone concentration. Our theory postulates a pronounced preference of leaf-cytosolic triose-phosphate isomerase to catalyse the forward reaction in vivo (dihydroxyacetone phosphate to glyceraldehyde 3-phosphate). In conclusion, intramolecular 13C/12C analysis resolves information about carbon uptake and allocation enabling more comprehensive assessments of carbon metabolism than whole-molecule 13C/12C analysis.

HighlightIntramolecular 13C/12C analysis resolves information about carbon uptake and allocation (and associated environmental controls) enabling more comprehensive assessments of carbon metabolism, plant-environment interactions, and environmental variability than whole-molecule 13C/12C analysis. LinearTurboFold: Linear-Time Global Prediction of Conserved Structures for RNA Homologs with Applications to SARS-CoV-2 The constant emergence of COVID-19 variants reduces the effectiveness of existing vaccines and test kits. Therefore, it is critical to identify conserved structures in SARS-CoV-2 genomes as potential targets for variant-proof diagnostics and therapeutics. However, the algorithms to predict these conserved structures, which simultaneously fold and align multiple RNA homologs, scale at best cubically with sequence length, and are thus infeasible for coronaviruses, which possess the longest genomes ([~]30,000 nt) among RNA viruses. As a result, existing efforts on modeling SARS-CoV-2 structures resort to single sequence folding as well as local folding methods with short window sizes, which inevitably neglect long-range interactions that are crucial in RNA functions. Here we present LinearTurboFold, an efficient algorithm for folding RNA homologs that scales linearly with sequence length, enabling unprecedented global structural analysis on SARS-CoV-2. Surprisingly, on a group of SARS-CoV-2 and SARS-related genomes, LinearTurbo-Folds purely in silico prediction not only is close to experimentally-guided models for local structures, but also goes far beyond them by capturing the end-to-end pairs between 5 and 3 UTRs ([~]29,800 nt apart) that match perfectly with a purely experimental work. Furthermore, LinearTurboFold identifies novel conserved structures and conserved accessible regions as potential targets for designing efficient and mutation-insensitive small-molecule drugs, antisense oligonucleotides, siRNAs, CRISPR-Cas13 guide RNAs and RT-PCR primers. LinearTurboFold is a general technique that can also be applied to other RNA viruses and full-length genome studies, and will be a useful tool in fighting the current and future pandemics.

Significance StatementConserved RNA structures are critical for designing diagnostic and therapeutic tools for many diseases including COVID-19. However, existing algorithms are much too slow to model the global structures of full-length RNA viral genomes. We present LinearTurboFold, a linear-time algorithm that is orders of magnitude faster, making it the first method to simultaneously fold and align whole genomes of SARS-CoV-2 variants, the longest known RNA virus ([~]30 kilobases). Our work enables unprecedented global structural analysis and captures long-range interactions that are out of reach for existing algorithms but crucial for RNA functions. LinearTurboFold is a general technique for full-length genome studies and can help fight the current and future pandemics. Cortical granularity shapes the organization of afferent paths to the amygdala and its striatal targets in nonhuman primate The prefrontal cortex (PFC) and insula, amygdala, and striatum form interconnected networks that drive motivated behaviors. We previously found a connectional trend in which granularity of the ventromedial and orbital PFC/insula predicted connections to the amygdala, and also the breadth of amygdalo-striatal efferents, including projections beyond the classic ventral striatum. To further interrogate connectional relationships among the cortex, amygdala, and striatum, and to further define the limbic (amygdala-recipient) striatum, we conducted tract tracing studies in two cohorts of Macaques (Male n = 14, Female n = 1). We focused on the cortico-amygdalo-striatal (indirect) and cortico- limbic striatal (direct) paths originating in the entire PFC and insula. Larger data sets and a quantitative approach revealed cortical rules in which cortical granularity predicts the complexity and location of projections to both the basal nucleus of the amygdala and striatum. Remarkably, projections from cortical-like basal nucleus to the striatum followed similar patterns. In both direct and indirect paths to the limbic striatum, agranular cortices formed a foundational, broad projection, and were joined by inputs from progressively more differentiated cortices. In amygdalo-striatal paths, the ventral basal nucleus was the  foundational input, with progressively more dorsal basal nucleus regions gradually adding inputs as the limbic striatum extended caudally. Together, the  indirect and  direct paths followed consistent principles in which cortical granularity dictated the strength and complexity of projections at their targets. Cluster analyses independently confirmed these connectional trends, and also highlighted connectional features that predicted termination in specific subregions of the basal nucleus and limbic striatum.

Significance StatementThe limbic system broadly refers to brain circuits that coordinate emotional responses. Here, we investigate circuits of the amygdala, which are involved in coding the emotional value of external cues, and their influence on the striatum. Regions of prefrontal cortex and insula form gradients of overlapping inputs to the amygdalas basal nucleus, which feed forward to the striatum. Direct cortical inputs to these amygdala-recipient striatal areas are surprisingly organized according to similar principles, but subtly shift from the classic ventral striatum to the caudal ventral striatum. Together, these distinct subsystems--cortico-amygdala-striatal circuits and direct cortico-striatal circuits-- provide substantial opportunity for different levels of internal, sensory, and external experiences to be integrated within the striatum, a major motor-behavioral interface. Characterization of methylation profiles in spontaneous preterm birth placental villous tissue Preterm birth is a global public health crisis which results in significant neonatal and maternal mortality. Yet little is known regarding the molecular mechanisms of idiopathic spontaneous preterm birth, and we have few diagnostic markers for adequate assessment of placental development and function. Previous studies of placental pathology and our transcriptomics studies suggest a role for placental maturity in idiopathic spontaneous preterm birth. It is known that placental DNA methylation changes over gestation. We hypothesized that if placental hypermaturity is present in our samples, we would observe a unique idiopathic spontaneous preterm birth DNA methylation profile potentially driving the gene expression differences we previously identified in our placental samples. Our results indicate the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern suggesting hypermaturity. Only seven significant differentially methylated regions fitting the idiopathic spontaneous preterm birth specific (relative to the controls) profile were identified, indicating unusually high similarity in DNA methylation between idiopathic spontaneous preterm birth and term birth samples. We identified an additional 1,718 significantly methylated regions in our gestational age matched controls were the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern, again indicating a striking level of similarity between the idiopathic spontaneous preterm birth and term birth samples. Pathway analysis of these regions revealed differences in genes within the WNT and Cadherin signaling pathways, both of which are essential in placental development and maturation. Taken together, these data demonstrate that the idiopathic spontaneous preterm birth samples are molecularly more mature than expected given their respective gestational age which likely impacts birth timing. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Changes to the mtDNA copy number during yeast culture growth We show that the mtDNA copy number in growing cultures of the yeast Saccharomyces cerevisiae increases by a factor of up to 4, being lowest ([~]10 per haploid genome) and stable during rapid fermentative growth, and highest at the end of the respiratory phase. When yeast are grown on glucose, the onset of the mtDNA copy number increase coincides with the early stages of the diauxic shift, and the increase continues through respiration. A lesser yet still substantial copy number increase occurs when yeast are grown on a nonfermentable carbon source, i.e. when there is no diauxic shift. The mtDNA copy number increase during and for some time after the diauxic shift is not driven by an increase in cell size. The copy number increase occurs in both haploid and diploid strains, but is markedly attenuated in a diploid wild isolate that is a ready sporulator. Strain-to-strain differences in mtDNA copy number are least apparent in fermentation and most apparent in late respiration or stationary phase. While changes in mitochondrial morphology and function were previously known to accompany changes in physiological state, it had not been previously shown that the mtDNA copy number changes substantially over time in a clonal growing culture. The mtDNA copy number in yeast is therefore a highly dynamic phenotype. LRR protein RNH1 dampens the inflammasome activation and is associated with adverse clinical outcomes in COVID-19 patients Inflammasomes are cytosolic innate immune sensors of pathogen infection and cellular damage that induce caspase-1 mediated inflammation upon activation. Although inflammation is protective, uncontrolled excessive inflammation can cause inflammatory diseases and can be detrimental, such as in COVID-19. However, the underlying mechanisms that control inflammasome activation are incompletely understood. Here we report that the leucine rich repeat (LRR) protein Ribonuclease inhibitor (RNH1), which shares homology with LRRs of NLRP proteins, attenuates inflammasome activation. Deletion of RNH1 in macrophages increases IL-1{beta} production and caspase-1 activation for inflammasome stimuli. Mechanistically, RNH1 decreases pro-IL-1{beta} expression and induces proteasome-mediated caspase-1 degradation. Corroborating this, mouse models of monosodium urate (MSU)-induced peritonitis and LPS-induced endotoxemia, which are dependent on caspase-1, respectively show increased neutrophil infiltration and lethality in Rnh1-/- mice compared to WT mice. Furthermore, RNH1 protein levels are negatively correlated with inflammation and disease severity in hospitalized COVID-19 patients. We propose that RNH1 is a new inflammasome regulator with relevance to COVID-19 severity.