Unsupervised clustering analysis of SARS-Cov-2population structure reveals six major subtypes atearly stage across the world Identifying the population structure of the newly emerged coronavirus SARS-CoV-2 has significant potential to inform public health management and diagnosis. As SARS-CoV-2 sequencing data accrued, grouping them into clusters is important for organizing the landscape of the population structure of the virus. Due to the limited prior information on the newly emerged coronavirus, we utilized four different clustering algorithms to group 16,873 SARS-CoV-2 strains, which automatically enables the identification of spatial structure for SARS-CoV-2. A total of six distinct genomic clusters were identified using mutation profiles as input features. Comparison of the clustering results reveals that the four algorithms produced highly consistent results, but the state-of-the-art unsupervised deep learning clustering algorithm performed best and produced the smallest intra-cluster pairwise genetic distances. The varied proportions of the six clusters within different continents revealed specific geographical distributions. In particular, our analysis found that Oceania was the only continent on which the strains were dispersively distributed into six clusters. In summary, this study provides a concrete framework for the use of clustering methods to study the global population structure of SARS-CoV-2. In addition, clustering methods can be used for future studies of variant population structures in specific regions of these fast-growing viruses.
RapidoPGS: A rapid polygenic score calculator for summary GWAS data without a test dataset MotivationPolygenic scores (PGS) aim to genetically predict complex traits at an individual level. PGS are typically trained on genome-wide association summary statistics and require an independent test dataset to tune parameters. More recent methods allow parameters to be tuned on the training data, removing the need for independent test data, but approaches are computationally intensive. Based on fine-mapping principles, we present RapidoPGS, a flexible and fast method to compute PGS requiring summary-level GWAS datasets only, with little computational requirements and no test data required for parameter tuning.

ResultsWe show that RapidoPGS performs slightly less well than two out of three other widely-used PGS methods (LDpred2, PRScs, and SBayesR) for case-control datasets, with median r2 difference: -0.0092, -0.0042, and 0.0064, respectively, but up to 17,000-fold faster with reduced computational requirements. RapidoPGS is implemented in R and can work with user-supplied summary statistics or download them from the GWAS catalog.

Availability and implementationOur method is available with a GPL license as an R package from GitHub.
Visual gene expression reveals a cone to rod developmental progression in deep-sea fishes Vertebrates use cone cells in the retina for colour vision and rod cells to see in dim light. Many deep-sea fishes have adapted to their environment to have only rod cells in the retina, while both rod and cone genes are still preserved in their genomes. As deep-sea fish larvae start their lives in the shallow, and only later submerge to the depth, they have to cope with diverse environmental conditions during ontogeny. Using a comparative transcriptomic approach in 20 deep-sea fish species from eight teleost orders, we report on a developmental cone-to-rod switch. While adults mostly rely on rod opsin (RH1) for vision in dim light, larvae almost exclusively express middle-wavelength-sensitive ("green") cone opsins (RH2) in their retinas. The phototransduction cascade genes follow a similar ontogenetic pattern of cone-followed by rod-specific gene expression in most species, except for the pearleye and sabretooth (Aulopiformes), in which the cone cascade remains dominant throughout development. By inspecting the whole genomes of five deep-sea species (four of them sequenced within this study: Idiacanthus fasciola, Chauliodus sloani; Stomiiformes; Coccorella atlantica, and Scopelarchus michaelsarsi; Aulopiformes), we found that deep-sea fish possess one or two copies of the rod RH1 opsin gene, and up to seven copies of the cone RH2 opsin genes in their genomes, while other cone opsin classes have been mostly lost. Our findings hence provide molecular evidence for a limited opsin gene repertoire and a conserved vertebrate pattern whereby cone photoreceptors develop first and rod photoreceptors are added only at later developmental stages. Loss of histone methyltransferase SETD1B in oogenesis results in the redistribution of genomic histone 3 lysine 4 trimethylation Histone 3 lysine 4 trimethylation (H3K4me3) is an epigenetic mark found at gene promoters and CpG islands. H3K4me3 is essential for mammalian development, yet mechanisms underlying its genomic targeting are poorly understood. H3K4me3 methyltransferases SETD1B and MLL2 are essential for oogenesis. We investigated changes in H3K4me3 in Setd1b conditional knockout (cKO) oocytes using ultra-low input ChIP-seq, with comparisons to DNA methylation and gene expression analyses. H3K4me3 was redistributed in Setd1b cKO oocytes showing losses at active gene promoters associated with downregulated gene expression. Remarkably, many regions also gained H3K4me3, in particular those that were DNA hypomethylated, transcriptionally inactive and CpG-rich, which are hallmarks of MLL2 targets. Consequently, loss of SETD1B disrupts the balance between MLL2 and de novo DNA methyltransferases in determining the epigenetic landscape during oogenesis. Our work reveals two distinct, complementary mechanisms of genomic targeting of H3K4me3 in oogenesis, with SETD1B linked to gene expression and MLL2 to CpG content.

Graphical AbstractIn oogenesis, SETD1B and CXXC1 target H3K4me3 to actively transcribed gene promoters, while MLL2 targets transcriptionally inactive regions based on underlying CpG composition (upper panel). When SETD1B is ablated, H3K4me3 is lost at a subset of active promoters, resulting in downregulation of transcription (lower panel). Loss of SETD1B alters the activity of MLL2, permitting MLL2 to deposit H3K4me3 at CpG-rich regions, many of which should otherwise be DNA methylated. Thus, it is evident that MLL2 and de novo DNMTs compete for genomic occupancy late in oogenesis, and loss of SETD1B disrupts the balance of these mechanisms.

O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY Climate warming changes synchrony of plants and pollinators Climate warming changes the phenology of many species. When interacting organisms respond differently, climate change may disrupt their interactions and affect the stability of ecosystems. Here, we used GBIF occurrence records to examine phenology trends in plants and their associated insect pollinators in Germany since the 1980s. We found strong phenological advances in plants, but differences in the extent of shifts among pollinator groups. The temporal trends in plant and insect phenologies were generally associated with interannual temperature variation, and thus likely driven by climate change. When examining the synchrony of species-level plant-pollinator interactions, their temporal trends differed among pollinator groups. Overall, plant-pollinator interactions become more synchronized, mainly since the phenology of plants responded more strongly to climate change than that of the pollinators. However, if the observed trends continue, many interactions may become more asynchronous again in the future. Our study suggests that climate change affects the phenologies of both plants and insects, and that it also influences the synchrony of plant-pollinator interactions. The non-specific Lipid Transfer Protein (nsLTP) is involved at early and late stages of symbiosis between Alnus glutinosa and Frankia alni. Alnus glutinosa response to Frankia alni is driven by several sequential physiological modifications that include calcium spiking, root hair deformation, penetration, induction of primordium, formation and growth of nodule. Here, we have conducted a transcriptomic study to analyse plant responses to Frankia alni at early stages of symbiosis establishment.

Forty-two genes were significantly activated by either with a Frankia culture supernatant or with living cells separated from the roots by a dialysis membrane permitted to identify plant genes which expression changes upon early contact with Frankia. Most of these genes encode biological processes, including oxidative stress and response to stimuli. The most upregulated gene is the non-specific lipid transfer protein (nsLTP) encoding gene with a fold change of 141. Physiological experiments showed that nsLTP increases Frankia nitrogen fixation at sub-lethal concentration. Immunohistochemistry experiments conducted at an early infection stage indicated that nsLTP protein is localized at the deformed root hair region after Frankia inoculation and later in nodules, precisely around bacterial vesicles. Taken together, these results suggest that nsLTP acts at early and late stages of symbiosis, probably by increasing nitrogen uptake by Frankia. 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. Structural basis for Nirmatrelvir in vitro efficacy against SARS-CoV-2 variants The COVID-19 pandemic continues to be a public health threat with emerging variants of SARS-CoV-2. Nirmatrelvir (PF-07321332) is a reversible, covalent inhibitor targeting the main protease (Mpro) of SARS-CoV-2 and the active protease inhibitor in PAXLOVID (nirmatrelvir tablets and ritonavir tablets). We evaluated the in vitro catalytic activity and in vitro potency of nirmatrelvir against the main protease (Mpro) of prevalent variants of concern (VOC) or variants of interest (VOI): Alpha (, B.1.1.7), Beta ({beta}, B.1.351), Delta ({delta}, B1.617.2), Gamma ({gamma}, P.1), Lambda ({lambda}, B.1.1.1.37/C37), Omicron (o, B.1.1.529) as well as the original Washington or wildtype strain. These VOC/VOI carry prevalent mutations at varying frequencies in the Mpro specifically for: , {beta},{gamma} (K90R), {lambda} (G15S) and o (P132H). In vitro biochemical enzymatic assay characterization of the enzyme kinetics of the mutant Mpros demonstrate that they are catalytically comparable to wildtype. Nirmatrelvir has similar potency against each mutant Mpro including P132H that is observed in the Omicron variant with a Ki of 0.635 nM as compared to a Ki of 0.933nM for wildtype. The molecular basis for these observations were provided by solution-phase structural dynamics and structural determination of nirmatrelvir bound to the o, {lambda} and {beta} Mpro at 1.63 - 2.09 [A] resolution. These in vitro data suggest that PAXLOVID has the potential to maintain plasma concentrations of nirmatrelvir many-fold times higher than the amount required to stop the SARS-CoV-2 VOC/VOI, including Omicron, from replicating in cells (1). The Impact of Multivesicular Release on the Transmission of Sensory Information by Ribbon Synapses The statistics of vesicle release determine how information is transferred at the synapse but the classical Poisson model of independent fusion does not hold at the first stages of vision and hearing. There, ribbon synapses encoding analogue signals can coordinate the release of two or more vesicles as a single event. The implications of such multivesicular release (MVR) for spike generation are not known. Here we investigate the hypothesis that MVR might provide advantages over Poisson synapses. We used leaky integrate-and-fire models incorporating the statistics of release measured experimentally from glutamatergic synapses of retinal bipolar cells and compared these with models assuming Poisson inputs constrained to operate at the same average rates. We find that MVR can increase the number of spikes generated per vesicle while reducing interspike intervals and latency to first spike. The combined effect was to increase the efficiency of information transfer (bits per spike) over a range of conditions mimicing retinal ganglion cells of different size. MVR was most advantageous in neurons with short time-constants and reliable synaptic inputs, when a lower degree of convergence was required to trigger spikes. In the special case of a single input driving a neuron with short time-constant, as occurs at the first synapse in hearing, MVR increased information transfer whenever spike generation required more than one vesicle. This study demonstrates how presynaptic integration of vesicles by MVR has the potential to compensate for less effective summation post-synaptically to increase the efficiency with which sensory information is transmitted. Sequence-dependent model of genes with dual σ factor preference Escherichia coli uses {sigma} factors to quickly control large gene cohorts during stress conditions. While most of its genes respond to a single {sigma} factor, approximately 5% of them have dual {sigma} factor preference. The most common are those responsive to both {sigma}70, which controls housekeeping genes, and {sigma}38, which activates genes during stationary growth and stresses. Using RNA-seq and flow-cytometry measurements, we show that  {sigma}70+38 genes are nearly as upregulated in stationary growth as  {sigma}38 genes. Moreover, we find a clear quantitative relationship between their promoter sequence and their response strength to changes in {sigma}38 levels. We then propose and validate a sequence dependent model of {sigma}70+38 genes, with dual sensitivity to {sigma}38 and {sigma}70, that is applicable in the exponential and stationary growth phases, as well in the transient period in between. We further propose a general model, applicable to other stresses and {sigma} factor combinations. Given this, promoters controlling {sigma}70+38 genes (and variants) could become important building blocks of synthetic circuits with predictable, sequence-dependent sensitivity to transitions between the exponential and stationary growth phases. The compound of baicalein, wogonin, and oroxylin-A inhibiting EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. Motor Protein MYO1C is Critical for Photoreceptor Opsin Trafficking and Visual Function Unconventional myosins linked to deafness are also proposed to play a role in retinal cell physiology. However, their direct role in photoreceptor function remains unclear. We demonstrate that systemic loss of the unconventional myosin MYO1C in mice specifically affected opsin trafficking, leading to loss of visual function. Electroretinogram analysis of Myo1c knockout (Myo1c-KO) mice showed a progressive loss of photoreceptor function. Immunohistochemistry and binding assays demonstrated MYO1C localization to photoreceptor inner and outer segments (OS) and identified a direct interaction of rhodopsin with the MYO1C cargo domain. In Myo1c-KO retinas, rhodopsin mislocalized to rod inner segments (IS) and cell bodies, while cone opsins in OS showed punctate staining. In aged mice, the histological and ultrastructural examination of the phenotype of Myo1c-KO retinas showed progressively shorter photoreceptor OS. These results demonstrate that MYO1C is critical for opsin trafficking to the photoreceptor OS and for normal visual function. A multivariate view of cognitive performance reveals positive correlation in the Trinidadian Guppy (Poecilia reticulata) Cognitive variation is common among-individuals and can be consistent across time and context. From an evolutionary perspective, among-individual variation is important as a pre-requisite for natural selection and adaptive evolution. Selection is widely hypothesized to favor high cognitive performance but directional selection should erode variation over time, how then is cognitive variation maintained? As selection does not act on traits in isolation, covariance among specific cognitive traits and/or other aspects of phenotype (e.g. personality) could result in fitness trade-offs that are important in shaping evolutionary dynamics. Here we test this using Trinidadian guppies (Poecilia reticulata), using a multivariate approach by characterising the correlation structure among task-specific cognitive performance measures and a personality trait. We estimate the among-individual correlation matrix (ID) in performance across three cognitive tasks; colour association learning task; motor learning task; reversal learning task, and the personality trait  boldness measured as emergence time from a shelter. We found no support for trade-offs among performance in these tasks. Nor do we find evidence of hypothesised speed-accuracy trade-offs within the association learning task. Rather we find strong positive correlation structure in ID, with 57% of variation explained by the leading eigen vector. While noting that non-cognitive factors and assay composition may affect the structure of ID, we suggest our findings are consistent with the g-model of cognitive performance variation, in which a dominant axis of variation loads positively on all performance measures. Thus, we add to a growing body of support for general variation among individuals in animal cognitive ability. GCN5L1 promotes diastolic dysfunction by inhibiting cardiac pyruvate oxidation Left ventricular diastolic dysfunction is a structural and functional condition that precedes the development of heart failure with preserved ejection fraction (HFpEF). The etiology of diastolic dysfunction includes alterations in fuel substrate metabolism that negatively impact cardiac bioenergetics, and may precipitate the eventual transition to heart failure. To date, the molecular mechanisms that regulate early changes in fuel metabolism leading to diastolic dysfunction remain unclear. In this report, we use a diet-induced obesity model and quantitative acetylproteomics in aged mice to show that inhibitory lysine acetylation of the pyruvate dehydrogenase (PDH) complex promotes energetic deficits and diastolic dysfunction in mouse hearts. Cardiomyocyte-specific deletion of the mitochondrial lysine acetylation regulatory protein GCN5L1 prevented hyperacetylation of the PDH complex subunit PDHA1, allowing aged obese mice to continue using pyruvate as a bioenergetic substrate in the heart. Our findings suggest that changes in mitochondrial protein lysine acetylation represent a key metabolic component of diastolic dysfunction that precedes the development of heart failure. Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable and highly reproducible in terms of microbial composition, function and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation. XRN2-mediated regulation of let-7 microRNAs is of critical pathophysiological significance in Humans MicroRNAs are critical regulators of diverse developmental and physiological processes, and they themselves get regulated both at the level of biogenesis and turnover. We demonstrate that the ribonuclease XRN2 can degrade the mature forms of certain let-7 family members in multiple human cancer cell lines, without affecting their precursors. XRN2 depletion results in a reduction in the expression of a number of oncogenes, and diminishes the proliferative and metastatic potential of cancer cells. The clinical relevance of these observations is also verified in tumour transcriptomics data from public RNA-sequencing datasets, where XRN2 mRNA expression is inversely correlated with the levels of a large number of miRNAs, including let-7 members, and high XRN2 mRNA levels are associated with poor survival in hepatocellular carcinoma, lung adenocarcinoma, and glioblastoma. We demonstrate that the miRNA is released by an as-yet unidentified proteinaceous  miRNA release factor from the grasp of Argonaute before its degradation. Our analyses of the patient-derived transcriptomics data also show that XRN2, via its regulation of let-7, affects multiple pathways in a consistent manner across epithelial and glial cell lineages, and thus, is of critical pathophysiological significance. Bioactive constituents of Verbena officinalis regulate killing efficiency of primary human natural killer cells by accelerating killing processes BackgroundNatural killer (NK) cells play a key role in eliminating tumorigenic and pathogen-infected cells. Verbena officinalis (V. officinalis) has been used as a medical plant in traditional and modern medicine, exhibiting anti-tumor and anti-inflammation activity, but its role in NK cells remains unclear.

PurposeThe impact of bioactive constituents of V. officinalis on immune responses remains largely elusive. In this work we investigated the potential targets of V. officinalis and focused on killing efficiency and related functions of NK cells regulated by bioactive constituents of V. officinalis.

Study design/methodsWe used primary human NK cells from peripheral blood mononuclear cells. Potential regulatory roles of selected compounds in NK killing function-related genes were first analyzed by network pharmacology approaches. Potential engagement between compounds and NK receptors was examined by molecular docking assay. Killing efficiency of NK cells was determined with real-time killing assay and live-cell imaging. Proliferation was examined by CFSE staining. Expression of cytotoxic proteins was analyzed using flow cytometry. Lytic granule release was quantified by CD107a degranulation assay. Contact time required for killing and determination of serial killers were analyzed using live cell imaging results.

ResultsResults from network pharmacology approaches and molecular docking analyses predicted that the five bioactive compounds (acteoside, apigenin, kaempferol, verbenalin and hastatoside) from V. officinalis can regulated various key genes essential for NK effector functions and can bind to NK inhibitory receptors as well as activating receptors, indicating their potential role in regulating NK killing capacity. Analyses of NK killing efficiency revealed that out of these five compounds verbenalin exhibited an encouraging potency to enhance NK killing efficiency. Further investigation showed that verbenalin did not affect proliferation, expression of cytotoxic proteins, or lytic granule degranulation, but rather reduced contact time required for killing, enhanced total killing events per NK cell, and elevated the fraction of killing competent NK cells.

ConclusionsIn this work, we have identified verbenalin as the most potent bioactive constituent from V. officinalis to enhance NK cell killing efficiency via accelerating killing processes and converting non-killers to killing competent NK cells without affecting the lytic granule pathway. Dimerization of kringle 1 domain from hepatocyte growth factor/scatter factor provides a potent minimal MET receptor agonist. Hepatocyte Growth Factor/Scatter Factor (HGF/SF) and its cognate receptor MET play several essential roles in embryogenesis and regeneration in post-natal life of epithelial organs such as liver, kidney, lung, and pancreas, prompting a strong interest in harnessing HGF/SF-MET signalling for regeneration of epithelial organs after acute or chronic damage. The limited stability and tissue diffusion of native HGF/SF, however, which reflect the tightly controlled, local mechanism of action of the morphogen, have led to a major search for HGF/SF mimics for therapy. In this work, we describe the rational design, production and characterisation of K1K1, a novel minimal MET agonist consisting of two copies of the kringle 1 domain of HGF/SF placed in tandem. K1K1 is highly stable and displays biological activities equivalent or superior to native HGF/SF in a variety of in vitro assay systems and in a mouse model of liver disease. These data suggest that this engineered ligand may find wide applications in acute and chronic diseases of the liver and other epithelial organs dependent on MET activation. The extracellular loops of OmpA control the slow rate of in vitro folding Outer membrane proteins are all beta barrels and these barrels have a variety of well-documented loop conformations. Here we test the effect of three different loop types on outer membrane protein A (OmpA) folding. We designed twelve 5-residue loops and experimentally tested the effect of replacing the long loops of outer membrane protein OmpA with the designed loops. Our studies succeeded in creating the smallest known outer membrane barrel. We find that significant changes in OmpA loops do not have a strong overall effect on OmpA folding. However, when decomposing folding into a fast rate and a slow rate we find that changes in loops strongly affect the slow rate of OmpA folding. Extracellular loop types with higher levels of hydrogen bonds had more instances of increasing the slow folding rate and extracellular loop types with low levels of hydrogen bonds had more instances of decreasing the slow folding rate. Having the slow rate affected by loop composition is consistent with the slow rate being associated with the insertion step of outer membrane protein folding. MT-MAG: Accurate and interpretable machine learning for complete or partial taxonomic assignments of metagenome-assembled genomes We propose MT-MAG, a novel machine learning-based software tool for the complete or partial hierarchically-structured taxonomic classification of metagenome-assembled genomes (MAGs). MT-MAG is capable of classifying large and diverse metagenomic datasets: a total of 245.68 Gbp in the training sets, and 9.6 Gbp in the test sets analyzed in this study. MT-MAG is, to the best of our knowledge, the first machine learning method for taxonomic assignment of metagenomic data that offers a "partial classification" option. MT-MAG outputs complete or partial classification paths, and interpretable numerical classification confidences of its classifications, at all taxonomic ranks. To assess the performance of MT-MAG, we define a "weighted classification accuracy," with a weighting scheme reflecting the fact that partial classifications at different ranks are not equally informative. For the two benchmarking datasets analyzed (genomes from human gut microbiome species, and bacterial and archaeal genomes assembled from cow rumen metagenomic sequences), MT-MAG achieves an average of 80.13% in weighted classification accuracy. At the Species level, MT-MAG outperforms DeepMicrobes, the only other comparable software tool, by an average of 35.75% in weighted classification accuracy. In addition, MT-MAG is able to completely classify an average of 67.7% of the sequences at the Species level, compared with DeepMicrobes which only classifies 47.45%. Moreover, MT-MAG provides additional information for sequences that are not completely classified at the Species level, resulting in the partial or complete classification of 95.15%, of the genomes in the datasets analyzed. In particular, MT-MAG classifies, on average, 88.84% of the test sequences to the Phylum level, 88.39% to the Class level, 86.81% to the Order level, 81.17% to the Family level, and 71.13% to the Genus level. Interferon alpha-based combinations suppress SARS-CoV-2 infection in vitro and in vivo There is an urgent need for new antivirals with powerful therapeutic potential and tolerable side effects. In the present study, we found that recombinant human interferon-alpha (IFNa) triggers intrinsic and extrinsic cellular antiviral responses, as well as reduces replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in vitro. Although IFNa alone was insufficient to completely abolish SARS-CoV-2 replication, combinations of IFNa with remdesivir or other antiviral agents (EIDD-2801, camostat, cycloheximide, or convalescent serum) showed strong synergy and effectively inhibited SARS-CoV-2 infection in human lung epithelial Calu-3 cells. Furthermore, we showed that the IFNa-remdesivir combination suppressed virus replication in human lung organoids, and that its single prophylactic dose attenuated SARS-CoV-2 infection in lungs of Syrian hamsters. Transcriptome and metabolomic analyses showed that the combination of IFNa-remdesivir suppressed virus-mediated changes in infected cells, although it affected the homeostasis of uninfected cells. We also demonstrated synergistic antiviral activity of IFNa2a-based combinations against other virus infections in vitro. Altogether, our results indicate that IFNa2a-based combination therapies can achieve higher efficacy while requiring lower dosage compared to monotherapies, making them attractive targets for further pre-clinical and clinical development. Regulation of DNA damage response by trimeric G-protein Signaling Upon sensing DNA double-strand breaks (DSBs), eukaryotic cells either die or repair DSBs via one of two competing pathways, i.e., non-homologous end-joining (NHEJ) or homologous recombination (HR). We show that cell fate after DNA damage hinges on the guanine nucleotide-exchange modulator of heterotrimeric G-protein, Gi*{beta}{gamma}, GIV/Girdin. GIV suppresses HR by binding and sequestering BRCA1, a key coordinator of multiple steps within the HR pathway, away from DSBs; it does so using a C-terminal motif that binds BRCA1s BRCT-modules via both phospho-dependent and -independent mechanisms. GIV promotes NHEJ, and binds and activates Gi and enhances the  free G{beta}{gamma}[-&gt;]PI-3-kinase[-&gt;]Akt pathway, thus revealing the enigmatic origin of prosurvival Akt signals during dsDNA repair. Absence of GIV, or the loss of either of its two functions impaired DNA repair, and induced cell death when challenged with numerous cytotoxic agents. That GIV selectively binds few other BRCT-containing proteins suggests convergent signaling such that heterotrimeric G-proteins may finetune sensing, repair, and outcome after DNA damage.

GRAPHIC ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=190 SRC="FIGDIR/small/452842v3_ufig1.gif" ALT="Figure 1">
View larger version (66K):
org.highwire.dtl.DTLVardef@3d5202org.highwire.dtl.DTLVardef@1d32d89org.highwire.dtl.DTLVardef@152a85dorg.highwire.dtl.DTLVardef@1ab9953_HPS_FORMAT_FIGEXP  M_FIG C_FIG HIGHLIGHTSO_LINon-receptor G protein modulator, GIV/Girdin binds BRCA1
C_LIO_LIBinding occurs in both canonical and non-canonical modes
C_LIO_LIGIV sequesters BRCA1 away from dsDNA breaks, suppresses HR
C_LIO_LIActivation of Gi by GIV enhances Akt signals, favors NHEJ
C_LI

IN BRIEFIn this work, the authors show that heterotrimeric G protein signaling that is triggered by non-receptor GEF, GIV/Girdin, in response to double-stranded DNA breaks is critical for decisive signaling events which favor non-homologous end-joining (NHEJ) and inhibit homologous recombination (HR). Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. Illuminating the plant transposon insertion landscape in real time using Cas9-targeted Nanopore sequencing and a novel pipeline Transposable elements (TEs) are a major source of genomic novelty, contributing to plant adaptation, speciation and new cultivar production. The often large, complex genomes of plants make identifying TE insertions from short reads challenging, while whole-genome sequencing remains expensive. In addition, rare somatic insertions reflecting mobilome dynamic are difficult to track using short reads. To overcome these challenges, we applied the Cas9-targeted Nanopore sequencing (CANS) approach and a novel pipeline, NanoCasTE, to trace both genetically inherited and somatic TE insertions in plants. We performed CANS of EVD retrotransposon in wild-type Arabidopsis thaliana and obtained up to 40x coverage of the targets after only a few hours of sequencing on a MinION sequencer. To estimate the ability to detect new TE and T-DNA insertions, we exploited the A. thaliana ddm1 mutant, which has elevated TE activity. Using CANS, we detect the hemizygous T-DNA insertion site and genetically inherited EVD insertions. CANS also detected >800 somatic TE insertions allowing monitoring of TE insertions in real time. By this, we built a high-density physical map of EVD insertion preferences in the A. thaliana ddm1 and found that EVD insertions are biased toward pericentromeric regions in ddm1 that differ from wild-type A. thaliana ecotypes suggesting a pivotal role of DNA methylation in the shaping of EVD insertion preferences. Using CANS, we were able to detect transposition events of ONSEN transposon family that occurred in response to heat stress and found hotspots of ONSEN insertions. Our results indicate that CANS and NanoCasTE are effective tools for the detection of TEIs in real time, opening an avenue for obtaining a deeper understanding of mobilome organization in response to stress and in different genetic backgrounds. The ability of CANS to detect genetically inherited TEIs and T-DNA integration sites make it useful for screening insertion mutants and transgenic plants. IRGM/Irgm1 Aggravates Progression of Atherosclerosis by Inducing Macrophage Apoptosis through the MAPK Signaling Pathway AimsAtherosclerosis underlies most cardiovascular diseases, among which acute coronary syndrome (ACS) caused by plaque rupture (PR) often leads to death. Immune-related GTPases (IRGM/Irgm1) have been extensively studied in inflammatory diseases, but their role in atherosclerosis is unclear. Determining how IRGM/Irgm1 promotes atherosclerotic plaque vulnerability will provide information for new biomarkers and/or therapeutic targets.

Methods and resultsWe identified ruptured and unruptured plaques by optical coherence tomography, and found that serum IRGM was highly expressed in patients with ST-segment elevation myocardial infarction in PR. We used ApoE-/-Irgm1+/+, ApoE-/-Irgm1+/- mice and chimeric mice to establish a model of advanced atherosclerosis. The results of pathological experiments showed that Irgm1 caused plaque necrosis. The ratio of neutral lipids and cholesterol crystals increases, while the content of collagen fibers decreases, aggravating the destabilization of atherosclerotic plaques. In vitro, we used multiple approaches to confirm that Irgm1 promotes macrophage apoptosis by promoting the production of reactive oxygen species and activating the MAPK signaling pathway.

ConclusionsIRGM may be a potential risk factor for PR. Mechanistic studies have shown that IRGM/Irgm1 contributes to the formation and rupture of fragile plaques. This is partly mediated by the induction of macrophage apoptosis via the MAPK signaling pathway. IRGM may offer new strategies for early treatment of ACS.

Translation viewOur findings indicate that IRGM/Irgm1 contributes to formation and rupture of vulnerable plaques. It suggests that IRGM may provide a new target for the early treatment of ACS. 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.