Learning interpretable cellular responses tocomplex perturbations in high-throughput screens Recent advances in multiplexed single-cell transcriptomics experiments are facilitating the high-throughput study of drug and genetic perturbations. However, an exhaustive exploration of the combinatorial perturbation space is experimentally unfeasible, so computational methods are needed to predict, interpret, and prioritize perturbations. Here, we present the compositional perturbation autoencoder (CPA), which combines the interpretability of linear models with the flexibility of deep-learning approaches for single-cell response modeling. CPA encodes and learns transcriptional drug responses across different cell type, dose, and drug combinations. The model produces easy-to-interpret embeddings for drugs and cell types, which enables drug similarity analysis and predictions for unseen dosage and drug combinations. We show that CPA accurately models single-cell perturbations across compounds, doses, species, and time. We further demonstrate that CPA predicts combinatorial genetic interactions of several types, implying that it captures features that distinguish different interaction programs. Finally, we demonstrate that CPA can generate in-silico 5,329 missing genetic combination perturbations (97.6% of all possibilities) with diverse genetic interactions. We envision our model will facilitate efficient experimental design and hypothesis generation by enabling in-silico response prediction at the single-cell level, and thus accelerate therapeutic applications using single-cell technologies.
The Clade Displacement Index: how to detect horizontal gene transfers in unrooted gene trees While most genes of any organism are inherited vertically - i.e. from its parent organisms - sometimes they can be exchanged between unrelated species in a process known as the horizontal gene transfer (HGT). Studies of HGT contribute to our knowledge about the mechanisms of evolution, including the emergence of new pathogens, and a great deal of effort has been put into different methods of finding transferred genes. The golden standard of HGT detection is the analysis of the incongruence between the gene and the species trees. Those methods typically require rooted trees, in which the direction of evolution is known. Gene trees are typically unrooted, and rooting them is yet another step in HGT analysis, prone to errors which may lead to wrong conclusions. A natural question arises: can HGTs be detected in gene trees without rooting them at all?

It turns out that, for a particular, yet broad, class of transfers, the answer to this question is: yes. It also turns out that the same methodology can be applied to complement the bootstrap support in assessing the stability of gene tree topology. In this article, we present the Clade Displacement Index, a measure of shift of a given clades location between two trees. We derive algorithms to compute it and give several examples of its applications to HGT detection and gene tree stability analysis. We finish by pointing out directions for further studies and an example that shows that not all HGTs are detectable without knowing the location of the root of the gene tree.

A Jupyter Notebook with the implementation and applications of CDI described in this paper is available at https://github.com/mciach/CDI
Mitofusin 1 and 2 regulation of mitochondrial DNA content is a critical determinant of glucose homeostasis The dynamin-like GTPases Mitofusin 1 and 2 (Mfn1 and Mfn2) are essential for mitochondrial function, which has been principally attributed to their regulation of fission/fusion dynamics. Here, we report that Mfn1 and 2 are critical for glucose-stimulated insulin secretion (GSIS) primarily through control of mitochondrial DNA (mtDNA) content. Whereas Mfn1 and Mfn2 individually were dispensable for glucose homeostasis, combined Mfn1/2 deletion in {beta}-cells reduced mtDNA content, impaired mitochondrial morphology and networking, and decreased respiratory function, ultimately resulting in severe glucose intolerance. Importantly, gene dosage studies unexpectedly revealed that Mfn1/2 control of glucose homeostasis was dependent on maintenance of mtDNA content, rather than mitochondrial structure. Mfn1/2 maintain mtDNA content by regulating the expression of the crucial mitochondrial transcription factor Tfam, as Tfam overexpression ameliorated the reduction in mtDNA content and GSIS in Mfn1/2-deficient {beta}-cells. Thus, the primary physiologic role of Mfn1 and 2 in {beta}-cells is coupled to the preservation of mtDNA content rather than mitochondrial architecture, and Mfn1 and 2 may be promising targets to overcome mitochondrial dysfunction and restore glucose control in diabetes. Task-independent metrics of computational hardness predict human cognitive performance The survival of human organisms depends on our ability to solve complex tasks in the face of limited cognitive resources. However, little is known about the factors that drive the complexity of those tasks. Here, building on insights from computational complexity theory, we quantify the computational hardness of cognitive tasks using a set of task-independent metrics related to the computational resource requirements of individual instances of a task. We then examine the relation between those metrics and human behavior and find that they predict both time spent on a task as well as accuracy in three canonical cognitive tasks. Our findings demonstrate that performance in cognitive tasks can be predicted based on generic metrics of their inherent computational hardness.

TeaserThe ability of humans to solve cognitive tasks is affected by generic mathematical properties of problems related to their computational complexity. 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. P5A-ATPase controls the ER translocation of Wnt in neuronal migration Wnt family are conserved secretory proteins required for developmental patterning and tissue homeostasis. The mechanisms underlying intracellular maturation and intercellular signal transduction of Wnt proteins have been extensively studied. However, how Wnt is targeted to the endoplasmic reticulum (ER) for processing and secretion remains elusive. Here we report that CATP-8/P5A-ATPase directs neuronal migration non-cell autonomously in C. elegans by regulating EGL-20/Wnt biogenesis. CATP-8 functions as a translocase to translocate EGL-20/Wnt nascent polypeptide into the ER by interacting with the hydrophobic core region of EGL-20 signal sequence. Such regulation of Wnt biogenesis by P5A-ATPase is conserved in human cells. These findings reveal physiological roles of P5A-ATPase in neural development and identify Wnt proteins as direct substrates of P5A-ATPase for ER translocation. 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. Muscle activation strategies of the vastus lateralis according to sex AimDespite men exhibiting greater muscle strength and fatigibility than women, it remains unclear if there are sex-based differences in muscle recruitment strategies e.g. motor unit (MU) recruitment and modulation of firing rate (FR) at normalised forces and during progressive increases in force.

MethodsTwenty-nine healthy male and thirty-one healthy female participants (18-35 years) were studied. Intramuscular electromyography was used to record individual motor unit potentials (MUPs) and near fibre MUPs from the vastus lateralis (VL) during 10% and 25% maximum isometric voluntary contractions (MVC), and spike-triggered averaging was used to obtain motor unit number estimates (MUNE) of the VL. Multilevel mixed-effects linear regression models were used to investigate the effects of sex at each contraction level.

ResultsMen exhibited greater muscle strength (p<0.001) and size (p<0.001) than women, with no difference in force steadiness at 10% or 25% MVC. Women had smaller MUs and higher FR at 10% MVC (both p<0.02), similar to that at 25% MVC in MU size (p=0.062) and FR (p=0.031). However, both sexes showed similar increases in MU size and FR when moving from low-to mid-level contractions. There were no sex differences in any near fibre MUP parameters or in MUNE.

ConclusionIn the vastus lateralis, women produce muscle force via different neuromuscular recruitment strategies to men which is characterised by smaller MUs discharging at higher rates. However, similar strategies are employed to increase force production from low to moderate contractions. These findings of similar proportional increases between sexes support the use of mixed sex cohorts in studies of this nature.

Key pointsO_LIIncreases in muscle force production are mediated by motor unit (MU) recruitment, and MU firing rate (FR).
C_LIO_LIWomen are underrepresented in studies of human neuromuscular research and markedly differ to men in a number of aspects of neuromuscular function, yet little is known of the recruitment strategies of each.
C_LIO_LIHere we demonstrate men and women have similar vastus lateralis MU number estimates, yet women recruit smaller MUs with higher FR than men at normalised contraction levels. However, increases in force are achieved via similar trajectories of MU recruitment and MU FR in men and women.
C_LIO_LIAlthough men and women exhibit divergent neuromuscular recruitment strategies to achieve normalised forces, increases in force are achived similarly and support the inclusion of mixed sex cohorts in studies of this nature.
C_LI 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. An integrative analysis to predict the Active Ingredients and explore Polypharmacological Mechanisms of Orthosiphon stamineus Benth. for human health BackgroundOrthosiphon stamineus Benth. was a dietary supplement or a traditional folk herb with widespread clinical application, but it was lack of comprehensive understanding of its active ingredients and polypharmacological mechanisms. The aim of this work was to systematically study its natural compounds and the molecular mechanisms of the O. stamineus via network pharmacology.

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

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

ConclusionOur study clarified the polypharmacological mechanisms via the relationship among active ingredients, targets, and diseases and provided better guidance for subsequent experiments and potential active ingredients for drug discovery or health promotion. TCF11 Has a Potent Tumor-Repressing Effect than Its Prototypic Nrf1α by Definition of both Similar yet Different Regulatory Profiles, with a Striking Disparity from Nrf2 Nrf1 and Nrf2, as two principal CNC-bZIP transcription factors, regulate similar but different targets involved in a variety of biological functions for maintaining cell homeostasis and organ integrity. Of note, the unique topobiological behavior of Nrf1 makes its functions more complicated than Nrf2, because it is allowed for alternatively transcribing and selectively splicing to yield multiple isoforms (e.g., TCF11, Nrf1). In order to gain a better understanding of their similarities and differences in distinct regulatory profiles, all four distinct cell models for stably expressing TCF11, TCF11{Delta}N, Nrf1 or Nrf2 have been herein established by an Flp-In T-REx-293 system and then identified by transcriptomic sequencing. Further analysis revealed that Nrf1 and TCF11 have similar yet different regulatory profiles, although both contribute basically to positive regulation of their co-targets, which are disparate from those regulated by Nrf2. Such disparity in those gene regulation by Nrf1 and Nrf2 was further corroborated by scrutinizing comprehensive functional annotation of their specific and/or common target genes. Conversely, the mutant TCF11{Delta}N, resulting from a deletion of the N-terminal amino acids 2-156 from TCF11, resembles Nrf2 with the largely consistent structure and function. Interestingly, our further experimental evidence demonstrates that TCF11 acts as a potent tumor-repressor relative to Nrf1, albeit both isoforms possess a congruous capability to prevent malignant growth of tumor and upregulate those genes critical for improving the survival of patients with hepatocellular carcinoma. Genome-wide analysis of the non-coding RNA synthetic genetic network reveals extensive plasticity and unique phenotypes in yeast Recent studies identified non-coding RNAs (ncRNAs) with unknown function that are responsible for major fitness changes in yeast. To understand ncRNA interplay and aid their functional assignment, the synthetic genetic array (SGA) methodology was employed to create >15,000 double mutants and to score their epistasis in different environments. Unlike the protein network, ncRNAs mostly displayed positive epistasis in rich medium. Interestingly, the negative interactions significantly increased under stressors, showing environmental-dependent functions for ncRNAs. No correlation was found between the network of ncRNAs and that of their neighbouring genes, suggesting functional independence. The U3 paralogs, SNR17A and SNR17B, share the majority of genetic interactions in rich medium as expected. For example, SUT480 interacted with both paralogs and its function was linked to 18S rRNA processing. However, under stressors, a large number of unique epistatic interactions were observed, supporting the notion that SNR17A and SNR17B have diverged and sub-functionalised after genome duplication. Major orchestration of shikimate, early phenylpropanoid and stilbenoid pathways by Subgroup 2 R2R3-MYBs in grapevine The stilbenoid pathway is responsible for the production of resveratrol and its derivatives in grapevine. A few transcription factors (TFs) have been previously identified as regulators of this pathway but the extent of this control is yet to be fully understood. Here we demonstrate how DNA affinity purification sequencing (DAP-Seq) allows for genome-wide TF binding site interrogation in a non-model species. We obtained 5,190 and 4,443 binding events assigned to 4,041 and 3,626 genes for MYB14 and MYB15, respectively (around 40% of peaks being located within -10kb of transcription start sites). DAP-Seq of MYB14 and MYB15 was combined with aggregate gene centred co-expression networks built from more than 1,400 transcriptomic datasets from leaves, fruits and flowers to narrow down bound genes to a set of high confidence targets. The analysis of MYB14, MYB15 and MYB13, a third uncharacterised member of Subgroup 2 (S2), showed that in addition to the few previously known stilbene synthase (STS ) targets, these three regulators bind to 30 out of 47 STS family genes. Moreover all three MYBs bind to several PAL, C4H and 4CL genes, in addition to shikimate pathway genes, the WRKY03 stilbenoid co-regulator and novel resveratrol-modifying gene candidates amongst which ROMT2 -3 were validated enzymatically. A high proportion of DAP-Seq bound genes was induced in the activated transcriptomes of transient MYB15 -overexpressing stilbenoid-producing grapevine leaves, validating our methodological approach for identifying gene regulatory networks of specialised metabolism. Overall, MYB genes from Subgroup 2 appear to play a key role in binding and directly regulating several primary and secondary metabolic steps leading to an increased flux towards stilbenoid production. Delayed induction of type I and III interferons mediates nasal epithelial cell permissiveness to SARS-CoV-2 The nasal epithelium is a plausible entry point for SARS-CoV-2, a site of pathogenesis and transmission, and may initiate the host response to SARS-CoV-2. Antiviral interferon (IFN) responses are critical to outcome of SARS-CoV-2. Yet little is known about the interaction between SARS-CoV-2 and innate immunity in this tissue. Here we applied single-cell RNA sequencing and proteomics to a primary cell model of human nasal epithelium differentiated at air-liquid interface. SARS-CoV-2 demonstrated widespread tropism for nasal epithelial cell types. The host response was dominated by type I and III IFNs and interferon-stimulated gene products. This response was notably delayed in onset relative to viral gene expression and compared to other respiratory viruses. Nevertheless, once established, the paracrine IFN response began to impact on SARS-CoV-2 replication. When provided prior to infection, recombinant IFN{beta} or IFN{lambda}1 induced an efficient antiviral state that potently restricted SARS-CoV-2 viral replication, preserving epithelial barrier integrity. These data suggest that the IFN-I/III response to SARS-CoV-2 initiates in the nasal airway and suggest nasal delivery of recombinant IFNs to be a potential chemoprophylactic strategy. Reconstitution of minimal motility system based on Spiroplasma swimming by expressing two bacterial actins in synthetic minimal bacterium Motility is one of the most important features of life, but its evolutionary origin remains unknown. In this study, we focused on Spiroplasma, commensal, or parasitic bacteria. They swim by switching the helicity of a ribbon-like cytoskeleton that comprises six proteins, each of which evolved from a nucleosidase and bacterial actin called MreB. We expressed these proteins in a synthetic, non-motile minimal bacterium, JCVI-syn3.0B, whose reduced genome was computer-designed and chemically synthesized. The synthetic bacterium exhibited swimming motility with features characteristic of Spiroplasma swimming. Moreover, some combinations of the two proteins produced a helical cell shape and swimming, suggesting that the swimming originated from the differentiation and coupling of bacterial actins, and we also obtained a minimal system for motility of the synthetic bacterium.

One-Sentence SummaryThe minimal system comprised two bacterial actins that provided cell helicity and swimming to the synthetic minimal bacterium. Nascent alt-protein chemoproteomics reveals a repressor of ribosome biogenesis Many unannotated microproteins and alternative proteins (alt-proteins) have recently been found to be co-encoded with canonical proteins, but few of their functions are known. Motivated by the hypothesis that alt-proteins undergoing active or stress-induced synthesis could play important cellular roles, here, we developed a chemoproteomic pipeline to identify nascent alt-proteins in human cells. We identified 22 actively translated unannotated alt-proteins, one of which is upregulated after DNA damage stress. We further defined MINAS-60 (MIcroprotein that Negatively regulates ASsembly of the pre-60S ribosomal subunit), a nucleolar localized alt-protein co-encoded with human RBM10.Depletion of MINAS-60 increases the amount of the mature 60S ribosomal subunit, consequently upregulating global protein synthesis and cell proliferation by repressing late-stage pre-60S assembly and export of the 60S ribosome subunit to the cytoplasm. Together, these results implicate MINAS-60 as a repressor of ribosome biogenesis, and demonstrate that chemoproteomics can enable generation of functional hypotheses for uncharacterized alt-proteins. A phenomics approach for in vitro antiviral drug discovery BackgroundThe current COVID-19 pandemic has highlighted the need for new and fast methods to identify novel or repurposed therapeutic drugs. Here we present a method for untargeted phenotypic drug screening of virus-infected cells, combining Cell Painting with antibody-based detection of viral infection in a single assay. We designed an image analysis pipeline for segmentation and classification of virus-infected and non-infected cells, followed by extraction of morphological properties.

ResultsWe show that the methodology can successfully capture virus-induced phenotypic signatures of MRC-5 human lung fibroblasts infected with Human coronavirus 229E (CoV-229E). Moreover, we demonstrate that our method can be used in phenotypic drug screening using a panel of nine host- and virus-targeting antivirals. Treatment with effective antiviral compounds reversed the morphological profile of the host cells towards a non-infected state.

ConclusionsThe method can be used in drug discovery for morphological profiling of novel antiviral compounds on both infected and non-infected cells. Multiple spillovers and onward transmission of SARS-Cov-2 in free-living and captive White-tailed deer Many animal species are susceptible to SARS-CoV-2 and could potentially act as reservoirs, yet transmission of the virus in non-human free-living animals has not been documented. White-tailed deer (Odocoileus virginianus), the predominant cervid in North America, are susceptible to SARS-CoV-2 infection, and experimentally infected fawns can transmit the virus. To test the hypothesis that SARS-CoV-2 may be circulating in deer, we tested 283 retropharyngeal lymph node (RPLN) samples collected from 151 free-living and 132 captive deer in Iowa from April 2020 through December of 2020 for the presence of SARS-CoV-2 RNA. Ninety-four of the 283 deer (33.2%; 95% CI: 28, 38.9) samples were positive for SARS-CoV-2 RNA as assessed by RT-PCR. Notably, between November 23, 2020 and January 10, 2021, 80 of 97 (82.5%; 95% CI 73.7, 88.8) RPLN samples had detectable SARS-CoV-2 RNA by RT-PCR. Whole genome sequencing of the 94 positive RPLN samples identified 12 SARS-CoV-2 lineages, with B.1.2 (n = 51; 54.5%), and B.1.311 (n = 19; 20%) accounting for ~75% of all samples. The geographic distribution and nesting of clusters of deer and human lineages strongly suggest multiple zooanthroponotic spillover events and deer-to-deer transmission. The discovery of sylvatic and enzootic SARS-CoV-2 transmission in deer has important implications for the ecology and long-term persistence, as well as the potential for spillover to other animals and spillback into humans. These findings highlight an urgent need for a robust and proactive "One Health" approach to obtaining a better understanding of the ecology and evolution of SARS-CoV-2.

One-Sentence SummarySARS-CoV-2 was detected in one-third of sampled white-tailed deer in Iowa between September 2020 and January of 2021 that likely resulted from multiple human-to-deer spillover and deer-to-deer transmission events. Orientation dependent proton transverse relaxation in the human brain white matter: The magic angle effect on a cylindrical helix PurposeTo overcome some limitations of prior orientation-dependent proton transverse relaxation formalisms in white matter (WM) with a novel framework based on the generalized magic angle effect function.

MethodsA cylindrical helix model was developed embracing anisotropic rotational and translational diffusion of restricted molecules in human brain WM, with the former characterized by an axially symmetric system. Transverse relaxation rates R2 and [Formula] were divided into isotropic [Formula] and anisotropic parts, [Formula], with  denoting an open angle and{varepsilon} 0 an orientation ({Phi}) offset from DTI-derived primary diffusivity direction. The proposed framework (Fit A) was compared with prior models without{varepsilon} 0 on previously published water and methylene proton transverse relaxation rates from developing, healthy, and pathological WM at 3T. Goodness of fit was represented by root-mean-square error (RMSE). F-test and linear correlation were used with statistical significance set to P [&le;] 0.05.

ResultsFit A significantly (P<0.01) outperformed prior models as demonstrated by reduced RMSEs, e.g., 0.349 vs. 0.724 in myelin water. Fitted{varepsilon} 0 was in good agreement with calculated{varepsilon} 0 from directional diffusivities. Compared with those from healthy adult, the fitted [Formula], and  from neonates were substantially reduced but{varepsilon} 0 increased, consistent with incomplete myelination. Significant positive [Formula] and negative ( and [Formula]) correlations were found with aging (demyelination) in elderly.

ConclusionThe developed framework can better characterize orientation dependences from a wide range of proton transverse relaxation measurements in human brain WM, shedding new light on myelin microstructural alterations at the molecular level. The Evolution of Protein Folds by Creative Destruction Mechanisms by which new protein folds emerge and diverge pose central questions in biological sciences. Incremental mutation and step-wise adaptation explain relationships between topologically similar protein folds. However, the universe of folds is diverse and riotous, suggesting roles of more potent and creative forces. Sequence and structure similarity are observed between topologically distinct folds, indicating that proteins with distinct folds may share common ancestry.

We found evidence of common ancestry between three distinct {beta}-barrel folds: OB, SH3 and cradle loop barrel (CLB). The data suggest a mechanism of fold evolution that interconverts SH3, OB and CLB. This mechanism, which we call creative destruction, can be generalized to explain other examples of fold evolution including circular permutation. In creative destruction, an open reading frame duplicates or otherwise merges with another to produce a fused polypeptide. A merger forces two ancestral domains into a new sequence and spatial context. The fused polypeptide can explore folding landscapes that are inaccessible to either of the independent ancestral domains. However, the folding landscapes of the fused polypeptide are not fully independent of those of the ancestral domains. Creative destruction is thus partially conservative in that a daughter fold would inherit some motifs from the ancestral folds. After a merger and refolding, adaptive processes such as mutation and loss of extraneous segments optimize the new daughter fold.

SignificanceMechanisms of emergence and early diversification of structured proteins present deep and difficult problems in evolutionary biology. Here we excavate the deepest evolutionary history, found within the translation machinery, which is an ancient molecular fossil and the birthplace of all proteins. We demonstrate common origins of some of the simplest, oldest and most common protein folds. Furthermore, the data suggest a mechanism, that we call creative destruction, that explains at molecular level how simple folds spawn more complex folds. In this mechanism, new folds emerge from old folds via gene duplication, expression, exploration of new folding landscapes and adaptation. Creative destruction explains the facile emergence of complex from simple architectures in a funneled exploration. TemplateFlow: FAIR-sharing of multi-scale, multi-species brain models Reference anatomies of the brain and corresponding atlases play a central role in experimental neuroimaging workflows and are the foundation for reporting standardized results. The choice of such references --i.e., templates-- and atlases is one relevant source of methodological variability across studies, which has recently been brought to attention as an important challenge to reproducibility in neuroscience. TemplateFlow is a publicly available framework for human and nonhuman brain models. The framework combines an open database with software for access, management, and vetting, allowing scientists to distribute their resources under FAIR --findable, accessible, interoperable, reusable-- principles. TemplateFlow supports a multifaceted insight into brains across species, and enables multiverse analyses testing whether results generalize across standard references, scales, and in the long term, species, thereby contributing to increasing the reliability of neuroimaging results. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. Development of a yeast-based assay for bioavailable phosphorus Preventing eutrophication of inland freshwater ecosystems requires quantifying the phosphorus (P) content of the streams and rivers that feed them. Typical methods for measuring P assess soluble reactive P (SRP) or total P (TP) and require expensive analytical techniques that produce hazardous waste. Here we present a novel method for measuring the more relevant bioavailable P (BAP); this assay utilizes the growth of familiar bakers yeast, avoids production of hazardous waste, and reduces cost relative to measurements of SRP and TP. The yeast BAP (yBAP) assay takes advantage of the observation that yeast density at saturating growth increases linearly with provided P. We show that this relationship can be used to measure P in freshwater in concentration ranges relevant to eutrophication. In addition, we measured yBAP in water containing known amounts of fertilizer and in samples from agricultural waterways. We observed that the majority of yBAP values were between those obtained from standard SRP and TP measurements, demonstrating that the assay is compatible with real-world settings. The cost-effective and nonhazardous nature of the yeast-based assay suggests that it could have utility in a range of settings, offering added insight to identify water systems at risk of eutrophication from excess phosphorus. Identification of a conserved neutralizing epitope present on spike proteins from all highly pathogenic coronaviruses Three pathogenic human coronaviruses have emerged within the last 20 years, with SARS-CoV-2 causing a global pandemic. Although therapeutic antibodies targeting the SARS-CoV-2 spike currently focus on the poorly conserved receptor-binding domain, targeting essential neutralizing epitopes on the more conserved S2 domain may provide broader protection. We report an antibody binding an epitope conserved in the pre-fusion core of MERS-CoV, SARS-CoV and SARS-CoV-2 spike S2 domains. Antibody 3A3 binds a conformational epitope with ~2.5 nM affinity and neutralizes spike from SARS-CoV, SARS-CoV-2 and variants of concern in in vitro pseudovirus assays. Hydrogen-deuterium exchange mass spectrometry identified residues 980-1006 in the flexible hinge region at the S2 apex as the 3A3 epitope, suggesting 3A3 prevents the S2 conformational rearrangements required for conversion to the spike post-fusion state and virus-host cell fusion. This work defines a conserved vulnerable site on the SARS-CoV-2 S2 domain and guides the design of pan-protective spike immunogens. Modeling miRNA-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.