Molecular diversity maintained by long-term balancing selection in mating loci defines multiple mating types in fungi Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of specimens. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci defines mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty specimens of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. Our sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Our analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. The inferred allelic information mirrored the outcome of in vitro crosses, thus allowing us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. The observed allelic classes could potentially generate 14,560 different mating types. We conclude that the genetic diversity of mating in Trichaptum loci is due to long-term balancing selection, with limited recombination and duplication activity. Our large number of sequenced specimens highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Author summaryFungi have complex mating systems, and basidiomycete fungi can encode thousands of mating types. Individuals with the same mating type cannot mate. This sexual system has evolved to facilitate sexual mating, increasing the chances to recombine into advantageous allelic combination and prune deleterious alleles. We explored the genomes of hundred and eighty specimens, combined with experimental mating studies of selected specimens, from a non-model organism (Trichaptum). We characterized the genomic regions controlling sex. The mating ability of the specimens confirmed the role of the mating alleles observed in the genomic data. The detailed analyses of many specimens allowed us to observe gene duplication and rearrangements within the mating loci, increasing the diversity within these loci. We supported previous suggestions of balancing selection in this region, an evolutionary force that maintains genomic diversity. These results supports that our fungal specimens are prone to outcross, which might facilitate the adaptation to new conditions.
An exon-intron split framework to prioritize miRNA-driven post-transcriptional regulatory signals and its application to study energy homeostasis 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.
GPI-anchored FGF directs cytoneme-mediated bidirectional signaling to self-regulate tissue-specific dispersion How signaling proteins generate a multitude of information to organize tissue patterns is critical to understanding morphogenesis. In Drosophila, FGF produced in wing-disc cells regulates the development of the disc-associated air-sac-primordium (ASP). Here, we show that FGF is Glycosylphosphatidylinositol-anchored to the producing cell surface and that this modification both inhibits free FGF secretion and activates target-specific bidirectional FGF-FGFR signaling through cytonemes. FGF-source and recipient ASP cells extend cytonemes that present FGF and FGFR on their surfaces and reciprocally recognize each other over distance by contacting through CAM-like FGF-FGFR binding. Contact-mediated FGF-FGFR binding induces bidirectional signaling, which, in turn, promotes ASP and source cells to polarize cytonemes toward each other and reinforce signaling contacts. Subsequent un-anchoring of FGFR-bound-FGF from the source cell membrane dissociates cytoneme contacts and delivers FGF target-specifically to ASP cytonemes for paracrine functions. Thus, GPI-anchored FGF organizes both source and recipient cells and self-regulates its cytoneme-mediated tissue-specific dispersion and signaling. Swine H1N1 influenza virus variants with enhanced polymerase activity and HA stability promote airborne transmission in ferrets Understanding how animal influenza A viruses (IAVs) acquire airborne transmissibility in humans and ferrets is needed to prepare for and respond to pandemics. Here, we investigated in ferrets the replication and transmission of swine H1N1 isolates P4 and G15, whose majority population had decreased polymerase activity and poor HA stability, respectively. For both isolates, a minor variant was selected and transmitted in ferrets. Polymerase-enhancing variant PA-S321 airborne-transmitted and propagated in one ferret. HA-stabilizing variant HA1-S210 was selected in all G15-inoculated ferrets and was transmitted by contact and airborne routes. With an efficient polymerase and a stable HA, the purified minor variant G15-HA1-S210 had earlier and higher peak titers in inoculated ferrets and was recovered at a higher frequency after airborne transmission than P4 and G15. Overall, HA stabilization played a more prominent role than polymerase enhancement in the replication and transmission of these viruses in ferrets. The results suggest pandemic risk-assessment studies may benefit from deep sequencing to identify minor variants with human-adapted traits.

IMPORTANCEDiverse IAVs circulate in animals, yet few acquire the viral traits needed to start a human pandemic. A stabilized HA and mammalian-adapted polymerase have been shown to promote the adaptation of IAVs to humans and ferrets (the gold-standard model for IAV replication, pathogenicity, and transmissibility). Here, we used swine IAV isolates of the gamma lineage a model to investigate the importance of HA stability and polymerase activity in promoting replication and transmission in ferrets. These are emerging viruses that bind to both -2,6- and -2,3-linked receptors. Using isolates containing mixed populations, a stabilized HA was selected within days in inoculated ferrets. An enhanced polymerase was also selected and propagated after airborne transmission to a ferret. Thus, HA stabilization was a stricter requirement, yet both traits promoted transmissibility. Knowing the viral traits needed for pandemic potential, and the relative importance of each, will help identify emerging viruses of greatest concern. Vapor exposure to Δ9-tetrahydrocannabinol (THC) slows locomotion of the Maine Lobster (Homarus americanus) RationaleDespite a long history of use in synaptic physiology, the lobster has been a neglected model for behavioral pharmacology. A restaurateur proposed that exposing lobster to cannabis smoke reduces anxiety and pain during the cooking process. It is unknown if lobster gill respiration in air would result in significant {Delta}9-tetrahydrocannabinol (THC) uptake and whether this would have any detectable behavioral effects.

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

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

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

ConclusionsVapor exposure of lobsters, using an e-cigarette based model, produces dose-dependent THC levels in all tissues and reduces locomotor activity. Hot water nociception was temperature dependent, but only minimal anti-nociceptive effect of THC exposure was confirmed. DeCalciOn: A hardware system for real-time decoding of in-vivo calcium imaging data Epifluorescence miniature microscopes ("miniscopes") are widely used for in vivo calcium imaging of neural population activity. Imaging data is usually collected while subjects are engaged in a task and stored for later offline analysis, but emerging techniques for online imaging offer potential for novel real-time experiments in which closed-loop interventions (such as neurostimulation or sensory feedback) are triggered at short latencies in response to neural population activity. Here we introduce DeCalciOn, a plug-and-play hardware device for online population decoding of in vivo calcium signals that can trigger closed-loop feedback at millisecond latencies, and is compatible with miniscopes that use the UCLA Data Acquisition (DAQ) interface. In performance tests, the position of rats (n=13) on a linear track was decoded in real time from hippocampal CA1 population activity by 24 linear classifiers. DeCalciOn required <2.5 ms after each end-of-frame to decode up to 1,024 calcium traces and trigger TTL control outputs. Decoding was most efficient using a  contour-free method of extracting traces from ROIs that were unaligned with neurons in the image, but  contour-based extraction from neuronal ROIs is also supported. DeCalciOn is an easy-to-use system for real-time decoding of calcium fluorescence that enables closed-loop feedback experiments in behaving animals. Library-free BoxCarDIA solves the missing value problem in label-free quantitative proteomics The last decade has seen significant advances in the application of quantitative mass spectrometry-based proteomics technologies to tackle important questions in plant biology. The current standard for quantitative proteomics in plants is the use of data-dependent acquisition (DDA) analysis with or without the use of chemical labels. However, the DDA approach preferentially measures higher abundant proteins, and often requires data imputation due to quantification inconsistency between samples. In this study we systematically benchmarked a recently developed library-free data-independent acquisition (directDIA) method against a state-of-the-art DDA label-free quantitative proteomics workflow for plants. We next developed a novel acquisition approach combining MS1-level BoxCar acquisition with MS2-level directDIA analysis that we call BoxCarDIA. DirectDIA achieves a 33% increase in protein quantification over traditional DDA, and BoxCarDIA a further 8%, without any changes in instrumentation, offline fractionation, or increases in mass-spectrometer acquisition time. BoxCarDIA, especially, offers wholly reproducible quantification of proteins between replicate injections, thereby addressing the long-standing missing-value problem in label-free quantitative proteomics. Further, we find that the gains in dynamic range sampling by directDIA and BoxCarDIA translate to deeper quantification of key, low abundant, functional protein classes (e.g., protein kinases and transcription factors) that are underrepresented in data acquired using DDA. We applied these methods to perform a quantitative proteomic comparison of dark and light grown Arabidopsis cell cultures, providing a critical resource for future plant interactome studies. Our results establish BoxCarDIA as the new method of choice in quantitative proteomics using Orbitrap-type mass-spectrometers, particularly for proteomes with large dynamic range such as that of plants. Plasma proteome analyses in individuals of European and African ancestry identify cis-pQTLs and models for proteome-wide association studies Improved understanding of genetic regulation of proteome can facilitate the identification of causal mechanisms for complex traits. We analyzed data on 4,657 plasma proteins from 7,213 European American (EA) and 1,871 African American (AA) individuals from the ARIC study, and further replicated findings on 467 AA individuals from the AASK study. Here we identified 2,004 proteins in EA and 1,618 in AA, with majority overlapping, which showed associations with common variants in cis-regions. Availability of AA samples led to smaller credible sets and significant number of population-specific cis-pQTLs. Elastic-net produced powerful models for protein prediction in both populations. An application of proteome-wide association studies (PWAS) to serum urate and gout, implicated several proteins, including IL1RN, revealing the promise of the drug anakinra to treat acute gout flares. Our study demonstrates the value of large and diverse ancestry study for genetic mechanisms of molecular phenotypes and their relationship with complex traits. Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector. A flexible electronic strain sensor for the real-time monitoring of tumor regression Assessing the efficacy of cancer therapeutics in mouse models is a critical step in treatment development. However, low resolution measurement tools and small sample sizes make determining drug efficacy in vivo a difficult and time-intensive task. Here we present a commercially scalable wearable electronic sensor that automates the in vivo testing of cancer therapeutics by continuously monitoring micrometer-scale tumor volume progression or regression in real-time. Histology, caliper measurements, and bioluminescence imaging over a one-week treatment period validated the sensors recordings, but the sensors higher time and length scale resolutions provided the only measurement technique capable of continuously monitoring the immediate hour-scale pharmacodynamic response of a given drug. In mice with two subcutaneously implanted tumor models our sensors recorded significant volume reductions in tumors just 5 hours after small molecule or immunotherapy treatment initiation. We anticipate that real-time tumor regression datasets could help expedite and automate the process of screening cancer therapies in vivo. Immune complex solubility and size govern Fc-gamma receptor responses: A scalable cell-based reporter system Fc-gamma receptor (Fc{gamma}R) activation by soluble IgG immune complexes (sICs) represents a major mechanism of inflammation in certain autoimmune diseases such as systemic lupus erythematosus (SLE). A robust and scalable test system allowing for the detection and quantification of sIC bioactivity is missing. Previously described Fc{gamma}R interaction assays are limited to certain Fc{gamma}Rs, lack scalability and flexibility, are not indicative of receptor activation or lack sensitivity towards sIC size. We developed a comprehensive reporter cell panel detecting individual activation of Fc{gamma}Rs from humans and the mouse. The reporter cell lines were integrated into an assay format that provides flexible read-outs enabling the quantification of sIC reactivity via ELISA or a fast detection using flow cytometry. This identified Fc{gamma}RIIA(H) and Fc{gamma}RIIIA as the most sIC-sensitive Fc{gamma}Rs in our test system. Applying the assay we demonstrate that sICs versus immobilized ICs are fundamentally different Fc{gamma}R-ligands with regard to Fc{gamma}R preference and signal strength. Reaching a detection limit in the very low nanomolar range, the assay proved also to be sensitive to sIC stoichiometry and size enabling for the first time a complete reproduction of the Heidelberger-Kendall precipitation curve in terms of immune receptor activation. Analyzing sera from SLE patients and mouse models of lupus and arthritis proved that sIC-dependent Fc{gamma}R activation has predictive capabilities regarding severity of SLE disease. The new methodology provides a sensitive, scalable and comprehensive tool to evaluate the size, amount and bioactivity of sICs in all settings.

One Sentence SummaryIn this study we established a comprehensive Fc{gamma}R reporter cell assay enabling the detection and quantification of soluble immune complexes generated in experimental and clinical settings. Plasmodium DDI1 is a potential dual-therapeutic target and essential chromatin-associated protein DDI1 proteins are involved in a variety of cellular processes, including proteasomal degradation of specific proteins. All DDI1 proteins contain a ubiquitin-like (UBL) domain and a retroviral aspartyl protease (RVP) domain. Some DDI1 proteins also contain a ubiquitin-associated (UBA) domain. The three domains confer distinct activities to DDI1 proteins. The presence of RVP domain makes DDI1 a potential target of HIV protease inhibitors, which also block the development of malaria parasites. Hence, we investigated the DDI1 of malaria parasites to identify its roles during parasite development and potential as a therapeutic target. DDI1 proteins of Plasmodium and other Apicomplexan parasites share the UBL-RVP domain architecture, and some also contain the UBA domain. Plasmodium DDI1 is expressed across all the major life cycle stages and is essential, as conditional depletion of DDI1 protein in the mouse malaria parasite Plasmodium berghei and the human malaria parasite Plasmodium falciparum compromised parasite development. Infection of mice with DDI1 knock-down P. berghei was self-limiting and protected the recovered mice from subsequent infection with homologous as well as heterologous parasites, indicating potential of DDI1 knock-down parasites as a whole organism vaccine. P. falciparum DDI1 (PfDDI1) is associated with chromatin and DNA-protein crosslinks. PfDDI1-depleted parasites accumulated DNA-protein crosslinks and showed enhanced susceptibility to DNA damaging chemicals, indicating a role of PfDDI1 in removal of DNA-protein crosslinks. Knock-down of PfDDI1 increased susceptibility to the retroviral protease inhibitor lopinavir and antimalarial artemisinin, which suggests that simultaneous inhibition of DDI1 could potentiate antimalarial activity of these drugs. Hence, Plasmodium DDI1 could be a dual-therapeutic target for malaria control. 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. 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. Force profile of the two-handed hardstyle kettlebell swing in novice older adults: an exploratory profile BackgroundUnderstanding the force profile of an exercise increases clinical confidence when assessing the benefits and potential risks of a prescribed exercise. This exploratory study presents the force profile of the hardstyle kettlebell swing in novice older adults and compares peak force with kettlebell deadlifts. These data will help inform healthcare providers and coaches who are considering prescribing kettlebell exercises for older adults.

MethodsThirty-five community-dwelling males and females (59-79 years) were recruited, from applicants to participate in the BELL trial. Two-handed hardstyle swings were performed with 8-16 kg. Deadlifts were performed with 8-24 kg and 8-32 kg for females and males, respectively. Ground reaction force was obtained from a floor-mounted force platform. Pairwise comparisons of peak force, forward force, rate of force development, swing cadence, sex, and kettlebell mass, were investigated for the kettlebell swing, with representative force-time curves described. Pairwise comparisons of peak force, sex and kettlebell mass were investigated for the deadlift, with comparisons of peak force, kettlebell mass, and sex, between swings and deadlifts.

ResultsFor kettlebells up to 16 kg, paired samples T-tests show a large exercise effect ({delta} > 1.4) with peak force higher for swings than deadlifts. Data shows: (i) higher peak force during swings than deadlifts ({delta} = 1.77), reaching 4.5 (1.0) N.kg-1, (ii) peak force during an 8 kg swing was greater than a 32kg deadlift, (iii) negligible difference in normalised peak force between males and females performing kettlebell swings, but a moderately large effect size during deadlifts (males > females, {delta} = 0.69), (iv) mean rate of force development of 19.9 (4.7) N.s-1.kg-1 with a very weak, positive correlation with kettlebell mass (y = 14.4 + 0.32x), and trivial effect of sex, (v) mean forward force equal to 5.5% of vertical force during swings, increasing from 3.8 (1.6) % with 8 kg to 7.1 (2.6) % with 16 kg.

ConclusionDuring kettlebell swings, there is negligible difference in normalised net peak force between novice males and females using the same absolute loads. Where ground reaction force is a therapeutic target, kettlebell swings with an 8 kg kettlebell could have similar effects to much heavier deadlifts (>24 kg). Kettlebell swings performed with lighter loads, could provider similar therapeutic value to much heavier deadlifts, and may be a more appealing, affordable, and convenient option for older adults. Decreased Interfacial Dynamics Caused by the N501Y Mutation in the SARS-CoV-2 S1 Spike:ACE2 Complex Corona Virus Disease of 2019 (COVID-19) caused by Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2) has caused a massive health crisis across the globe, with some genetic variants gaining enhanced infectivity and competitive fitness, and thus significantly aggravating the global health concern. In this regard, the recent SARS-CoV-2 alpha variant, B.1.1.7 lineage, reported from the United Kingdom (UK), is of great significance in that it contains several mutations that increase its infection and transmission rates as evident from clinical reports. Specifically, the N501Y mutation in the SARS-CoV-2 spike S1 receptor binding domain (S1-RBD) has been shown to possess an increased affinity for ACE2, although the basis for this is not entirely clear yet. Here, we dissect the mechanism underlying the increased affinity using molecular dynamics (MD) simulations of the available ACE2-S1-RBD complex structure (6M0J) and show a prolonged and stable interfacial interaction of the N501Y mutant S1-RBD with ACE2 compared to the wild type S1-RBD. Additionally, we find that the N501Y mutant S1-RBD displays altered dynamics that likely aids in its enhanced interaction with ACE2. By elucidating a mechanistic basis for the increased affinity of the N501Y mutant S1-RBD for ACE2, we believe that the results presented here will aid in developing therapeutic strategies against SARS-CoV-2 including designing drugs targeting the ACE2-S1-RBD interaction.

SignificanceThe emergence of the new SARS-CoV-2 lineage in the UK in December 2020 has further aggravated the COVID-19 pandemic due to an increased ability of the variant to infect human hosts, likely due to mutations in the viral S1 spike protein including the N501Y S1-RBD mutation that is located at the interface of S1-RBD and ACE2, the host cell receptor for SARS-CoV-2. Given its location at the interface, N501Y S1-RBD mutation can therefore potentially alter the interfacial interaction. Multiple, all-atom, explicit solvent MD simulations of the ACE2-S1-RBD complex carried here indicated a more stable interaction between the N501Y mutant S1-RBD and ACE2 through stabilizing interfacial interactions of residues at one end of the interface that are either sequentially or physically near the mutation site. These mechanistic details will aid in better understanding the mechanism by which the alpha variant has increased infectivity as well as in designing better therapeutics including ACE2-S1 spike protein inhibitors that will, in turn, help thwarting the current and future pandemic.

HighlightsO_LIN501 in the wild type SARS-CoV-2 S1-RBD forms unsustained hydrogen bonds with residues in the ACE2, namely Y41 and K353
C_LIO_LIY501 in the N501Y mutant SARS-CoV-2 S1-RBD is not capable of forming substantial hydrogen bonds with ACE2 within the time span of the current simulation
C_LIO_LIEvidence from analyzing the simulation results suggests that Y501 of S1-RBD could form other types of non-covalent interactions with ACE2, such as van der Waals interactions
C_LIO_LIN501Y S1-RBD mutation stabilizes the position of interfacial residues neighboring to the mutation site, as well as other non-interfacial residues that are distant from the mutation site
C_LIO_LIThese altered dynamics results in more stable interaction of S1-RBD with ACE2 which could be the main reason underlying the reported enhanced affinity of S1-RBD in the SARS-CoV-2 alpha variant (UK B.1.1.7 lineage) to ACE2
C_LI Chromatin accessibility landscape and active transcription factors in primary human invasive lobular and ductal breast carcinomas Background: Invasive lobular breast carcinoma (ILC), the second most prevalent histological subtype of breast cancer, exhibits unique molecular features compared with the more common invasive ductal carcinoma (IDC). While genomic and transcriptomic features of ILC and IDC have been characterized, genome-wide chromatin accessibility pattern differences between ILC and IDC remain largely unexplored. Methods: Here, we characterized tumor-intrinsic chromatin accessibility differences between ILC and IDC using primary tumors from The Cancer Genome Atlas (TCGA) breast cancer assay for transposase-accessible chromatin with sequencing (ATAC-seq) dataset. Results: We identified distinct patterns of genome-wide chromatin accessibility in ILC and IDC. Inferred patient-specific transcription factor (TF) motif activities revealed regulatory differences between and within ILC and IDC tumors. EGR1, RUNX3, TP63, STAT6, SOX family, and TEAD family TFs were higher in ILC, while ATF4, PBX3, SPDEF, PITX family, and FOX family TFs were higher in IDC. Conclusions: This study reveals the distinct epigenomic features of ILC and IDC and the active TFs driving cancer progression that may provide valuable information on patient prognosis. Higher-order effects, continuous species interactions, and trait evolution shape microbial spatial dynamics The assembly and maintenance of microbial diversity in natural communities, despite the abundance of toxin-based antagonistic interactions, presents major challenges for biological understanding. A common framework for investigating such antagonistic interactions involve cyclic dominance games with pairwise interactions. The incorporation of higher-order interactions in such models permits increased levels of microbial diversity, especially in communities where antibiotic producing, sensitive, and resistant strains co-exist. However, most such models involve a small number of discrete species, assume a notion of pure cyclic dominance, and focus on low mutation rate regimes, none of which well represents the highly interlinked, quickly evolving, and continuous nature of microbial phenotypic space. Here, we present an alternative vision of spatial dynamics for microbial communities based on antagonistic interactions--one in which a large number of species interact in continuous phenotypic space, are capable of rapid mutation, and engage in both direct and higher-order interactions mediated by production of and resistance to antibiotics. Focusing on toxin production, vulnerability, and inhibition among species, we observe highly divergent patterns of diversity and spatial community dynamics. We find that species interaction constraints (rather than mobility) best predict spatiotemporal disturbance regimes, whereas community formation time, mobility, and mutation size best explain patterns of diversity. We also report an intriguing relationship among community formation time, spatial disturbance regimes, and diversity dynamics. This relationship, which suggests that both higher-order interactions and rapid evolution are critical for the origin and maintenance of microbial diversity, has broad-ranging links to the maintenance of diversity in other systems.

Significance StatementPersistently diverse microbial communities are one of biologys great puzzles. Using a novel continuous trait space modeling framework that accommodates high mutation rates, elevated species richness, and direct and higher-order antagonistic species interactions, we find that two parameters characterizing mutation size and mobility best explain patterns of microbial diversity. Moreover, community formation time (the duration of the transient phase in community assembly) provides an unexpectedly clear guide to the diversity profiles of the resulting communities. These discoveries showcase how complex, antagonistic interactions mediated by the production of, inhibition of, and vulnerability to toxins (antibiotics) can shape microbial communities, allowing for extraordinarily high levels of diversity and temporal persistence. Polycomb-mediated repression of paternal chromosomes maintains haploid dosage in diploid embryos of Marchantia Complex mechanisms regulate gene dosage throughout eukaryotic life cycles. Mechanisms controlling gene dosage have been extensively studied in animals, however it is unknown how generalizable these mechanisms are to diverse eukaryotes. Here, we use the haploid plant Marchantia polymorpha to assess gene dosage control in its short-lived diploid embryo. We show that throughout embryogenesis, paternal chromosomes are repressed resulting in functional haploidy. The paternal genome is targeted for genomic imprinting by the Polycomb mark H3K27me3 starting at fertilization, rendering the maternal genome in control of embryogenesis. Maintaining haploid gene dosage by this new form of imprinting is essential for embryonic development. Our findings illustrate how haploid-dominant species can regulate gene dosage through paternal chromosome inactivation and initiates the exploration of the link between life cycle history and gene dosage in a broader range of organisms. How grower decision-making is modelled impacts the effect of plant disease control. While the spread of plant disease depends strongly on biological factors driving transmission, it also has a human dimension. Disease control depends on decisions made by individual growers, who are in turn influenced by a broad range of factors. Despite this, human behaviour has rarely been included in plant epidemic models. Considering Cassava Brown Streak Disease, we model how the perceived increase in profit due to disease management influences participation in clean seed systems (CSS). Our models are rooted in game theory, with growers making strategic decisions based on the expected profitability of different control strategies. We find that both the information used by growers to assess profitability and the perception of economic and epidemiological parameters influence long-term participation in the CSS. Over-estimation of infection risk leads to lower participation in the CSS, as growers perceive that paying for the CSS will be futile. Additionally, even though good disease management can be achieved through the implementation of CSS, and a scenario where all controllers use the CSS is achievable when growers base their decision on the average of their entire strategy, CBSD is rarely eliminated from the system. These results are robust to stochastic and spatial effects. Our work highlights the importance of including human behaviour in plant disease models, but also the significance of how that behaviour is included. Conformity leading self-organised rescue from a suboptimal decision-making bias Given the ubiquity of potentially adverse behavioural bias owing to myopic trial-and-error learning, it seems paradoxical that improvements in decision-making performance through conformist social learning, a process widely considered to be bias amplification, still prevail in animal collective behaviour. Here we show, through model analyses and large-scale interactive behavioural experiments with 585 human subjects, that conformist influence can indeed promote favourable risk taking in repeated experience-based decision making, even though many individuals are systematically biased towards adverse risk aversion. Although strong positive feedback conferred by copying the majoritys behaviour could result in unfavourable informational cascades, our differential equation model of collective behavioural dynamics identified a key role for increasing exploration by negative feedback arising when a weak minority influence undermines the inherent behavioural bias. This  collective behavioural rescue, emerging through coordination of positive and negative feedback, highlights a benefit of collective learning in a broader range of environmental conditions than previously assumed and resolves the ostensible paradox of adaptive collective behavioural flexibility under conformist influences. The ecological consequences and evolution of retron-mediated suicide as a way to protect bacteria from being killed by phage In 1984, retrons were described as DNA sequences that code for a reverse transcriptase and a unique single-stranded DNA/RNA hybrid called multicopy single-stranded DNA (msDNA). It would not be until 2020 that a function was shown, when compelling evidence was presented that retrons can protect bacteria from infections by lytic phage via an abortive infection mechanism. When infected with the virulent mutant of the phage lambda, {lambda}VIR, and to a lesser extent other phages, a retron designated Ec48 is activated, the Escherichia coli bearing this element dies and the infecting phage is lost. With the aid of a mathematical model, we explored the a priori conditions under which retrons will protect bacterial populations from invasion by phage and when retron-bearing bacteria will evolve in populations without this element. Using E. coli with and without Ec48, and {lambda}VIR we estimated the parameters of this model and tested the hypotheses generated from our analysis of its properties. Our models and experiments suggest cells expressing retron-mediated abortive infection have no advantage over cells lacking this system as a defense against phage, and we question whether the protection against phage infection is the unique selective force responsible for their evolution in bacteria. 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. Accurate estimation of cell-type resolution transcriptome in bulk tissue through matrix completion Single cell RNA-seq (scRNA-seq) has been widely used to uncover cellular heterogeneity, however, the constraints of cost make it impractical as a routine on large patient cohorts. Here we present ENIGMA, a method that accurately deconvolute bulk tissue RNA-seq into single cell-type resolution given the knowledge gained from scRNA-seq. ENIGMA applies a matrix completion strategy to minimize the distance between mixture transcriptome and weighted combination of cell type-specific expression, allowing quantification of cell type proportions and reconstruction of cell type-specific transcriptome. The superior performance of ENIGMA was validated in simulated and realistic datasets, including disease-related tissues, demonstrating its ability in novel biological findings. 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. NUDT18 catalyzes the hydrolysis of active metabolites of the antivirals Remdesivir, Ribavirin and Molnupiravir Remdesivir and Molnupiravir have gained considerable interest due to their activity against SARS-CoV-2. Cellular hydrolysis of their active triphosphate forms, Remdesivir-TP and Molnupiravir-TP, would decrease drug efficiency. We therefore tested Remdesivir-TP as a substrate against a panel of human hydrolases and found that NUDT18 catalyzes the hydrolysis of Remdesivir-TP. The kcat value of NUDT18 for Remdesivir-TP was determined to 2.6 s-1 and the Km value was 156 M, suggesting that NUDT18 catalyzed hydrolysis of Remdesivir-TP occurs in cells. We demonstrate that the triphosphates of the antivirals Ribavirin and Molnupiravir are hydrolyzed by NUDT18, albeit with a lower efficiency compared to Remdesivir-TP. NUDT18 also hydrolyses the triphosphates of Sofosbuvir and Aciclovir although with significantly lower activity. These results suggest that NUDT18 can act as a cellular sanitizer of modified nucleotides and may influence the antiviral efficacy of Remdesivir, Molnupiravir and Ribavirin. NUDT18 is expressed in respiratory epithelial cells and may limit the antiviral efficacy of Remdesivir and Molnupiravir against SARS-CoV2 replication by decreasing the intracellular concentration of their active metabolites at their intended site of action.