title: Insights into G-Quadruplex-Hemin Dynamics Using Atomistic Simulations: Implications for Reactivity and Folding

abstract:
Guanine quadruplex nucleic acids (G4s) are involved in key biological processes such as replication or transcription. Beyond their biological relevance, G4s find applications as biotechnological tools since they readily bind hemin and enhance its peroxidase activity, creating a G4-DNAzyme. The biocatalytic properties of G4-DNAzymes have been thoroughly studied and used for biosensing purposes. Despite hundreds of applications and massive experimental efforts, the atomistic details of the reaction mechanism remain unclear. To help select between the different hypotheses currently under investigation, we use extended explicit-solvent molecular dynamics (MD) simulations to scrutinize the G4/hemin interaction. We find that besides the dominant conformation in which hemin is stacked atop the external G-quartets, hemin can also transiently bind to the loops and be brought to the external G-quartets through diverse delivery mechanisms. The simulations do not support the catalytic mechanism relying on a wobbling guanine. Similarly, catalytic role of the iron-bound water molecule is not in line with our results, however, given the simulation limitations, this observation should be considered with some caution. The simulations rather suggest tentative mechanisms in which the external G-quartet itself could be responsible for the unique H2O2-promoted biocatalytic properties of the G4/hemin complexes. Once stacked atop a terminal G-quartet, hemin rotates about its vertical axis while readily sampling shifted geometries where the iron transiently contacts oxygen atoms of the adjacent G-quartet. This dynamics is not apparent from the ensemble-averaged structure. We also visualize transient interactions between the stacked hemin and the G4 loops. Finally, we investigated interactions between hemin and on-pathway folding intermediates of the parallel-stranded G4 fold. The simulations suggest that hemin drives the folding of parallel-stranded G4s from slip-stranded intermediates, acting as a G4 chaperone. Limitations of the MD technique are briefly discussed.

For Table of Contents Only

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=77 SRC="FIGDIR/small/375691v2_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@384380org.highwire.dtl.DTLVardef@1d7c730org.highwire.dtl.DTLVardef@1562248org.highwire.dtl.DTLVardef@1e4cb9c_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Evaluating and Improving SSU rRNA PCR Primer Coverage for Bacteria, Archaea, and Eukaryotes Using Metagenomes from Global Ocean Surveys

abstract:
Small subunit ribosomal RNA (SSU rRNA) amplicon sequencing can quantitatively and comprehensively profile natural microbiomes, representing a critically important tool for studying diverse global ecosystems. However, results will only be accurate if PCR primers perfectly match the rRNA of all organisms present. To evaluate how well marine microorganisms across all 3 domains are detected by this method, we compared commonly-used primers with > 300 million rRNA gene sequences retrieved from globally-distributed marine metagenomes. The best-performing primers when comparing to 16S rRNA of Bacteria and Archaea were 515Y/926R and 515Y/806RB, which perfectly matched over 96% of all sequences. Considering Cyanobacteria and Chloroplast 16S rRNA, 515Y/926R had the highest coverage (99%), making this set ideal for quantifying marine primary producers. For eukaryotic 18S rRNA sequences, 515Y/926R also performed best (88%), followed by V4R/V4RB (18S rRNA-specific; 82%) - demonstrating that the 515Y/926R combination performs best overall for all 3 domains. Using Atlantic and Pacific Ocean samples, we demonstrate high correspondence between 515Y/926R amplicon abundances (generated for this study) and metagenomic 16S rRNA (median R2=0.98, n=272), indicating amplicons can produce equally accurate community composition data versus shotgun metagenomics. Our analysis also revealed that expected performance of all primer sets could be improved with minor modifications, pointing toward a nearly-completely universal primer set that could accurately quantify biogeochemically-important taxa in ecosystems ranging from the deep-sea to the surface. In addition, our reproducible bioinformatic workflow can guide microbiome researchers studying different ecosystems or human health to similarly improve existing primers and generate more accurate quantitative amplicon data.

Significance StatementPCR amplification and sequencing of marker genes is a low-cost technique for monitoring prokaryotic and eukaryotic microbial communities across space and time, but will only work optimally if environmental organisms match PCR primer sequences exactly. In this study, we evaluated how well primers match globally-distributed short-read oceanic metagenomes. Our results demonstrate primer sets vary widely in performance, and that at least for marine systems, rRNA amplicon data from some primers lack significant biases compared to metagenomes. We also show that it is possible to create a nearly universal primer set for diverse saline environments by defining a specific mixture of a few dozen oligonucleotides and present a software pipeline that can guide rational design of primers for any environment with available metaomic data.

category: ecology


title: iSMNN: Batch Effect Correction for Single-cell RNA-seq data via Iterative Supervised Mutual Nearest Neighbor Refinement

abstract: Batch effect correction is an essential step in the integrative analysis of multiple single cell RNA-seq (scRNA-seq) data. One state-of-the-art strategy for batch effect correction is via unsupervised or supervised detection of mutual nearest neighbors (MNNs). However, both two kinds of methods only detect MNNs across batches on the top of uncorrected data, where the large batch effect may affect the MNN search. To address this issue, we presented iSMNN, a batch effect correction approach via iterative supervised MNN refinement across data after correction. Our benchmarking on both simulation and real datasets showed the advantages of the iterative refinement of MNNs on the performance of correction. Compared to popular alternative methods, our iSMNN is able to better mix the cells of the same cell type across batches. In addition, iSMNN can also facilitate the identification of differentially expression genes (DEGs) relevant to the biological function of certain cell types. These results indicated that iSMNN will be a valuable method for integrating multiple scRNA-seq datasets that can facilitate biological and medical studies at single-cell level.

category: bioinformatics


title: DeepDrug: A general graph-based deep learning framework for drug-drug interactions and drug-target interactions prediction

abstract: Computational approaches for accurate prediction of drug interactions, such as drug-drug interactions (DDIs) and drug-target interactions (DTIs), are highly demanded for biochemical researchers due to the efficiency and cost-effectiveness. Despite the fact that many methods have been proposed and developed to predict DDIs and DTIs respectively, their success is still limited due to a lack of systematic evaluation of the intrinsic properties embedded in the corresponding chemical structure. In this paper, we develop a deep learning framework, named DeepDrug, to overcome the above limitation by using residual graph convolutional networks (RGCNs) and convolutional networks (CNNs) to learn the comprehensive structural and sequential representations of drugs and proteins in order to boost the DDIs and DTIs prediction accuracy. We benchmark our methods in a series of systematic experiments, including binary-class DDIs, multi-class/multi-label DDIs, binary-class DTIs classification and DTIs regression tasks using several datasets. We then demonstrate that DeepDrug outperforms state-of-the-art methods in terms of both accuracy and robustness in predicting DDIs and DTIs with multiple experimental settings. Furthermore, we visualize the structural features learned by DeepDrug RGCN module, which displays compatible and accordant patterns in chemical properties and drug categories, providing additional evidence to support the strong predictive power of DeepDrug. Ultimately, we apply DeepDrug to perform drug repositioning on the whole DrugBank database to discover the potential drug candidates against SARS-CoV-2, where 3 out of 5 top-ranked drugs are reported to be repurposed to potentially treat COVID-19. To sum up, we believe that DeepDrug is an efficient tool in accurate prediction of DDIs and DTIs and provides a promising insight in understanding the underlying mechanism of these biochemical relations. The source code of the DeepDrug can be freely downloaded from https://github.com/wanwenzeng/deepdrug.


title: Exploring protein sequence and functional spaces using adversarial autoencoder

abstract:
MotivationShedding light on the relationships between protein se-quences and functions is a challenging task with many implications in protein evolution, diseases understanding, and protein design. Protein sequence / function space is however hard to comprehend due to its com-plexity. Generative models help to decipher complex systems thanks to their abilities to learn and recreate data specificity. Applied to protein sequences, they can point out relationships between protein positions and functions capture the sequence patterns associated with functions or navigate through uncharted area of molecular evolution.

ResultsIn this study, an unsupervised generative approach based on adversarial auto-encoder (AAE) is proposed to generate and explore new sequences with respect to their functions thanks to the prior distribution allowing a continuous exploration of the latent space. AAEs are tested on three protein families known for their multiple functions. Clustering re-sults on the encoded sequences from the latent space computed by AAEs display high level of homogeneity regarding the protein sequence func-tions. The study also reports and analyzes for the first time two sampling strategies based on latent space interpolation and latent space arithmetic to generate intermediate protein sequences sharing sequential and functional properties of original sequences issued from different families and functions. Generated sequences by interpolation between latent space data points demonstrate the ability of the AAE to generalize and to pro-duce meaningful biological sequences from an evolutionary uncharted area of the biological sequence space. Finally, 3D structure models generated by comparative modelling between different combinations of structures of different sub-families and of generated sequences from latent space or sub-family sequences point out to the ability of the latent space arithmetic to successfully transfer functional properties between sub-families. All in all this study confirms the ability of deep learning frameworks to model biological complexity and bring new tools to explore amino acid sequence and functional spaces.

AvailabilityCode and data used for this study are freely available at https://github.com/T-B-F/aae4seq.

Contacttristan@bitardfeildel.fr

Supplementary informationSupplementary data are available at online.


title: Navigating the amino acid sequence space between functional proteins using a deep learning framework


title: Voltage-gated calcium channels contribute to spontaneous glutamate release directly via nanodomain coupling or indirectly via calmodulin

abstract:
Neurotransmitter release occurs either synchronously with action potentials (evoked release) or spontaneously (spontaneous release). Whether the molecular mechanisms underlying evoked and spontaneous release are identical, especially whether voltage-gated Ca2+ channels (VGCCs) can trigger spontaneous events, is still a matter of debate in glutamatergic synapses. To elucidate this issue, we characterized the VGCC dependence of miniature excitatory postsynaptic currents (mEPSCs) in various synapses with different coupling distances between VGCCs and synaptic vesicles, known as a critical factor in evoked release. We found that most of the extracellular calcium-dependent mEPSCs were attributable to VGCCs in cultured autaptic hippocampal neurons and the mature calyx of Held where VGCCs and vesicles were tightly coupled. Among loosely coupled synapses, mEPSCs were not VGCC-dependent at immature calyx of Held and CA1 pyramidal neuron synapses, whereas VGCCs contribution was significant at CA3 pyramidal neuron synapses. Interestingly, the contribution of VGCCs to spontaneous glutamate release in CA3 pyramidal neurons was abolished by a calmodulin antagonist, calmidazolium. These data suggest that coupling distance between VGCCs and vesicles determines VGCC dependence of spontaneous release at tightly coupled synapses, yet VGCC contribution can be achieved indirectly at loosely coupled synapses.

HighlightsO_LIExtracellular Ca2+ -dependent spontaneous glutamate release is mediated by VGCCs.
C_LIO_LICa2+ cooperativity of spontaneous and evoked release is comparable.
C_LIO_LISynapses at different areas and ages show distinct VGCC dependence on mini events.
C_LIO_LINanodomain coupling between VGCCs and Ca2+ sensors determines VGCC dependence.
C_LI

category: neuroscience


title: Integrative Transcriptomic Analysis Reveals the Immune Mechanism for A CyHV-3-Resistant Common Carp Strain

abstract: Anti-disease breeding is becoming the most promising solution to cyprinid herpesvirus-3 (CyHV-3) infection, the major threat to common carp aquaculture. Virus challenging studies suggested that a breeding strain of common carp is resistant to this disease. This study illustrates the immune mechanisms involved in anti-CyHV-3 ability. An integrative analysis of the protein-coding genes and long non-coding RNAs (lncRNAs) using transcriptomic data was performed. Tissues from the head kidney of common carp were extracted at day 0 (the healthy control) and day 7 after CyHV-3 infection (the survivors), and used to analyze the transcriptome through both Illumina and PacBio sequencing. Following analysis of the GO terms and KEGG pathways involved, the immune-related terms and pathways were merged. In order to dig out details in the immune aspect, The DEGs were filtered using the current common carp immune gene library. Immune gene categories and their corresponding genes in different comparison groups were revealed. Also, the immunological Gene Ontology terms for lncRNA modulation were retained. The weighted gene co-XSexpression network analysis was used to reveal the regulation of immune genes by lncRNA. The results demonstrated that the breeding carp strain develops marked resistance to CyHV-3 through a specific innate immune mechanism. The featured biological processes were autophagy, phagocytosis, cytotoxicity, and virus blockage by lectins and MUC3. Moreover, the immune suppressive signals, such as suppression of IL21R on STAT3, PI3K mediated inhibition on inflammation by dopamine upon infection, as well as the inhibition of NLRC3 on STING during a steady state. Possible susceptible factors for CyHV-3, such as ITGB1, TLR18, and CCL4, were also revealed from the non-breeding strain. The results of this study also suggested that Nramp and PAI regulated by LncRNA could facilitate virus infection and proliferation for infected cells respectively, while T cell leukemia homeobox 3 (TLX3) as well as galectin 3 function by lncRNA may play a role in the resistance mechanism. Therefore, immune factors that are immunogenetically insensitive or susceptible to CyHV-3 infection have been revealed.

category: immunology


title: Dihydroceramide desaturase directs the switch between exosome production and autophagy for neuronal maintenance

abstract:
Exosomes are important for cell-cell communication. Deficiencies in the human dihydroceramide desaturase gene, DEGS1, increase the dihydroceramide-to-ceramide ratio and causes hypomyelinating leukodystrophy. However, the disease mechanism remains unknown. Here, we developed an in vivo assay with spatially controlled expression of exosome markers in Drosophila eye imaginal discs and showed that the level and activity of the DEGS1 ortholog, ifc, correlated with exosome production. Knocking out ifc decreased the density of the exosome precursor intraluminal vesicles (ILVs) in the multivesicular endosomes and reduced the number of exosomes released. While ifc overexpression and autophagy inhibition both enhanced exosome production, combining the two had no additive effect. Moreover, DEGS1 activity was sufficient to drive ILV formation in vitro. Together, DEGS1/Ifc controls the dihydroceramide-to-ceramide ration and enhances exosome secretion by promoting ILV formation and preventing the autophagic degradation of MVEs.

These findings provide a potential cause for the neuropathy associated with DEGS1-deficient mutations.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=115 SRC="FIGDIR/small/376046v3_ufig1.gif" ALT="Figure 1">
View larger version (43K):
org.highwire.dtl.DTLVardef@c7479corg.highwire.dtl.DTLVardef@12e7f46org.highwire.dtl.DTLVardef@1a5e7f0org.highwire.dtl.DTLVardef@11e7f08_HPS_FORMAT_FIGEXP  M_FIG C_FIG HighlightsO_LIAn in vivo system was developed for observing exosome production in Drosophila.
C_LIO_LIDihydroceramide desaturase (DEGS1/Ifc) promotes exosome production by two means.
C_LIO_LIIfc drives membrane invagination to promote the formation of intraluminal vesicles.
C_LIO_LIIfc inhibits autophagic degradation of MVEs and increases exosome release.
C_LI

eTOC BlurbThe level and activity of dihydroceramide desaturase (human DEGS1 and Drosophila Ifc) correlate with exosome production. Wu et al. show that DEGS1 drives the formation of intraluminal vesicles in vivo and in vitro. Overexpressing Ifc inhibits autophagy and reduces the degradation of multivesicular endosomes, thus increases exosome release.

category: cell biology


title: Dihydroceramide desaturase promotes the formation of intraluminal vesicles and inhibits autophagy to increase exosome production


title: Biparatopic sybody constructs neutralize SARS-CoV-2 variants of concern and mitigate emergence of drug resistance

abstract: The ongoing COVID-19 pandemic represents an unprecedented global health crisis. Here, we report the identification of a synthetic nanobody (sybody) pair (Sb#15 and Sb#68) that can bind simultaneously to the SARS-CoV-2 spike-RBD and efficiently neutralize pseudotyped and live-viruses by interfering with ACE2 interaction. Two spatially-discrete epitopes identified by cryo-EM translated into the rational design of bispecific and tri-bispecific fusions constructs, exhibiting up to 100- and 1000-fold increase in neutralization potency. Cryo-EM of the sybody-spike complex further revealed a novel up-out RBD conformation. While resistant viruses emerged rapidly in the presence of single binders, no escape variants were observed in presence of the bispecific sybody. The multivalent bispecific constructs further increased the neutralization potency against globally-circulating SARS- CoV-2 variants of concern. Our study illustrates the power of multivalency and biparatopic nanobody fusions for the development of clinically relevant therapeutic strategies that mitigate the emergence of new SARS-CoV-2 escape mutants.

category: biochemistry


title: Soluble SorCS1 binds the insulin receptor to enhance insulin sensitivity.

abstract: Type 2 diabetes has reached endemic proportions and is a substantial burden for the affected patients and the society. Along with lifestyle factors, a number of genetic loci predisposing to type 2 diabetes have been identified, including SORCS1 that encodes the transmembrane receptor SorCS1. The ectodomain of SorCS1 (sol-SorCS1) is shed from plasma membranes but the biological function of this fragment is unknown. Here we show that sol-SorCS1 acts as a high-affinity binding partner for the insulin receptor to stabilize the receptor and increase insulin affinity, protein kinase B activation, and glucose uptake in myocytes. Sol-SorCS1 is liberated from adipocytes, and in diabetic patients the plasma concentration positively correlates with body mass index, but inversely with plasma glucose. In mouse models of insulin resistance, exogenous sol-SorCS1 restored insulin sensitivity. We conclude that sol-SorCS1 increases peripheral insulin sensitivity and propose sol-SorCS1 as a novel insulin sensitizing adipokine and potential antidiabetic agent.

category: physiology


title: EGCG-modified bone graft to modulate the recruitment of M1 macrophage and alleviate the forming of fibrous capsule

abstract: In response to current trends in the modification of guided bone regeneration (GBR) materials, we aimed to build upon our previous studies on epigallocatechin-3-gallate (EGCG) by immersing a commonly used bone graft primarily composed of hydroxyapatite (HA) in EGCG solution, expecting to obtain superior bone-material integration after implantation. Bone grafts are commonly used for bone repair, in which the bone extracellular matrix is stimulated to promote osteogenesis. However, due to its pro-fibrosis effect, this osteoconductive material commonly exhibits implant failure. In addition to providing a basic release profile of EGCG-modified bone graft (E-HA) to clarify the relationship between this material and the environment, we have examined the integration effect via subcutaneous implantation experiments. In this manner, we have assessed the aggregation of pro-inflammatory macrophages, initial angiogenesis, the formation of fibrous capsules, and an enhanced cell viability observed in cultured RAW 264.7 cells. Among these results, we focus on pro-inflammatory macrophages due to their close relationship with fibrosis, which is the most important process in the immune response. Immunofluorescent staining results showed that E-HA substantially compromised the formation of fibrous capsules in hematoxylin-eosin-stained sections, which exhibited less pro-inflammatory macrophage recruitment; meanwhile, the cell viability and primary angiogenesis were improved. This work lays the foundation for future studies on GBR.


title: Streptococcus suis induces expression of cyclooxygenase-2 in porcine lung tissue

abstract: Streptococcus suis is a common pathogen colonising the respiratory tract of pigs. It can cause meningitis, sepsis and pneumonia leading to economic losses in the pig industry worldwide. Cyclooxygenase-2 (COX-2) and its metabolites play an important regulatory role in different biological processes like modulation of inflammation and immune activation. In this report we analysed the induction of COX-2 and the production of its metabolite prostaglandin E2 (PGE2) in a porcine precision-cut lung slice (PCLS) model. Using Western blot analysis, we found a time-dependent induction of COX-2 in the infected tissue resulting in increased PGE2 levels. Immunohistological analysis revealed a strong COX-2 expression in the proximity of the bronchioles between the ciliated epithelial cells and the adjacent alveolar tissue. The morphology, location and vimentin staining suggested that these cells are subepithelial bronchial fibroblasts. Furthermore, we showed that COX-2 expression as well as PGE2 production was detected following infection with two prevalent S. suis serotypes and that the pore-forming toxin suilysin played an important role in this process. Therefore, this study provides new insights in the response of porcine lung cells to S. suis infections and serves as a basis for further studies to define the role of COX-2 and its metabolites in the inflammatory response in porcine lung tissue during infection with S. suis.

category: microbiology


title: SF3B1 promotes tumor malignancy through splicing-independent co-activation of HIF1α

abstract: Mutations in the splicing factor SF3B1 are frequently occurring in various cancers and drive tumor progression through the activation of cryptic splice sites in multiple genes. Recent studies also demonstrate a positive correlation between expression levels of wildtype SF3B1 and tumor malignancy, but underlying mechanisms remain elusive. Here, we report that SF3B1 acts as an activator of HIF signaling through a splicing-independent mechanism. We demonstrate that SF3B1 forms a heterotrimer with HIF1 and HIF1{beta}, facilitating binding of the HIF1 complex to hypoxia response elements (HREs) to activate target gene expression. We further validate the relevance of this mechanism for tumor progression. Monoallelic deletion of Sf3b1 impedes formation and progression of hypoxic pancreatic cancer via impaired HIF signaling, but is well tolerated in normoxic chromophobe renal cell carcinoma. Our work uncovers an essential role of SF3B1 in HIF1 signaling, providing a causal link between high SF3B1 expression and aggressiveness of solid tumors.


title: Biochemical and structural studies of target lectin SapL1 from the emerging opportunistic microfungus Scedosporium apiospermum

abstract: Scedosporium apiospermum is an emerging opportunistic fungal pathogen responsible for life-threatening infections in humans. Host-pathogen interactions often implicate lectins that have become therapeutic targets for the development of carbohydrate mimics for antiadhesive therapy. Here, we present the first report on the identification and characterization of a lectin from S. apiospermum named SapL1. SapL1 was found using bioinformatics as a homolog to the conidial surface lectin FleA from Aspergillus fumigatus known to play a role in the adhesion to host glycoconjugates present in human lung epithelium. In our strategy to obtain recombinant SapL1, we discovered the importance of osmolytes to achieve its expression in soluble form in bacteria. Analysis of glycan arrays indicates specificity for fucosylated oligosaccharides as expected. Submicromolar affinity was measured for fucose using isothermal titration calorimetry. We solved SapL1 crystal structure in complex with -methyl-L-fucoside and analyzed its structural basis for fucose binding. We finally demonstrated that SapL1 binds to bronchial epithelial cells in a fucose-dependent manner. The information gathered here will contribute to the design and development of glycodrugs targeting SapL1.


title: Multi-omics Analysis of Dsup Expressing Human Cells Reveals Open Chromatin Architectural Dynamics Underyling Radioprotection

abstract: Spaceflight has been documented to produce detrimental effects to physiology and genomic stability, partly a result of Galactic Cosmic Radiation (GCR). In recent years, extensive research into extremotolerant organisms has begun to reveal how they survive harsh conditions, such as ionizing radiation. One such organism is the tardigrade (Ramazzottius varieornatus) which can survive up to 5kGy of ionizing radiation and the vacuum of space. In addition to their extensive network of DNA damage response mechanisms, the tardigrade also possesses a unique damage suppressor protein (Dsup) that co-localizes with chromatin in both tardigrade and transduced human cells to protect against DNA damage from reactive oxygen species induced by ionizing radiation. While Dsup has been shown to confer human cells with increased radiotolerance; much of the mechanism of how it does this in the context of human cells remains unknown. Until now there is no knowledge yet of how introduction of Dsup into human cells can perturb molecular networks and if there are any systemic risks associated with foreign gene introduction. Here, we created a stable HEK293 cell line expressing Dsup, validated its radioprotective phenotype, and performed multi-omic analyses across different time points and doses of radiation to delineate molecular mechanism of the radioprotection and assess molecular network pertubations. Dsup expressing human cells showed an enrichment for pathways seen in cells overexpressing HMGN1, a chromosomal architectural protein that has a highly similar nucleosome binding motif. As HMGN1 binding to nucleosomes promotes a less transcriptionally repressed chromatin state, we further explored the hypothesis that Dsup could behave similarly via ATAC-seq analysis and discovered overall selective differential opening and closing of the chromatin landscape. Cut&Run analysis further revealed global increases in histone post translational modifications indicative of open chromatin and global decreases in repressive marks, with Dsup binding preferentially towards promoter regions marked by H3K27ac and H3K4me3. We further validated some of the enriched pathways via in-vitro assays and revealed novel phenotypes that Dsup confers to human cells such as reduction in apoptosis, increased cell proliferation, and increased cell adhesion properties. Our analysis provides evidence that the Dsup protein in the context of HEK293 cells may behave as a chromatin architectural protein and that in addition to its nucleosome shielding effect, may confer radio-resistance via chromatin modulation. These results provide future insight into mitigating some of the major challenges involved with long term spaceflight as well as understanding some of the molecular architectural underpinnings that lead to radioresistant cancer phenotypes back home.

category: synthetic biology


title: Attention and prediction modulations in expected and unexpected visuospatial trajectories.

abstract: Humans are constantly exposed to a rich tapestry of visual information in a potentially changing environment. To cope with the computational burden this engenders, our perceptual system must use prior context to simultaneously prioritise stimuli of importance and suppress irrelevant surroundings. This study investigated the influence of prediction and attention in visual perception by investigating event-related potentials (ERPs) often associated with these processes, N170 and N2pc for prediction and attention, respectively. A contextual trajectory paradigm was used which violated visual predictions and neglected to predetermine areas of spatial interest, to account for the potentially unpredictable nature of a real-life visual scene. Participants (N=36) viewed a visual display of cued and non-cued shapes rotating in a five-step predictable trajectory, with the fifth and final position of either the cued or non-cued shape occurring in a predictable or unpredictable spatial location. To investigate the predictive coding theory of attention we used factors of attention and prediction, whereby attention was manipulated as either cued or non-cued conditions, and prediction manipulated in either predictable or unpredictable conditions. Results showed both enhanced N170 and N2pc amplitudes to unpredictable compared to predictable stimuli. Stimulus cueing status also increased N170 amplitude, but this did not interact with stimulus predictability. The N2pc amplitude was not affected by stimulus cueing status. In accordance with previous research these results suggest the N170 is in part a visual prediction error response with respect to higher-level visual processes, and furthermore the N2pc may index attention reorientation. The results demonstrate prior context influences the sensitivity of the N170 and N2pc electrophysiological responses. These findings add further support to the role of N170 as a prediction error signal and suggest that the N2pc may reflect attentional reorientation in response to unpredicted stimulus locations.


title: The Autism Risk Factor CHD8 Is a Chromatin Activator and Functionally Dependent on the ERK-MAPK Pathway Effector ELK1

abstract: The chromodomain helicase DNA-binding protein CHD8 is among the most frequently found de-novo mutations in autism spectrum disorder (ASD)1-4. Despite its prominent disease involvement, little is known about its molecular function in the human brain. CHD8 is believed to be a chromatin regulator, but mechanisms for its genomic targeting is also unclear. To elucidate the role of CHD8 in human neurons, we generated conditional loss-of-function alleles in pluripotent stem cells. Chromatin accessibility and transcriptional profiling showed that CHD8 is a potent chromatin opener and transcriptional activator of its direct neuronal targets, including a distinct group of ASD genes. We found the chromatin targeting of CHD8 to be highly context-dependent. In human neurons, CHD8 was preferentially bound at promoter sequences which were significantly enriched in ETS motifs. Indeed, the chromatin state of ETS motif-containing promoters was preferentially affected upon loss of CHD8. Among the many ETS transcription factors, we found ELK1 to be the best correlated with CHD8 expression in primary human fetal and adult cortical neurons and most highly expressed in our ES cell-derived neurons. Remarkably, ELK1 was necessary to recruit CHD8 specifically to ETS motifcontaining sites. These findings imply the functional cooperativity between ELK1, a key downstream factor of the MAPK/ERK pathway, and CHD8 on chromatin involvement in human neurons. THEREFORE, the MAPK/ERK/ELK1 axis may also play a role in the pathogenesis caused by CHD8 mutations 5.


title: NyuWa Genome Resource: Deep Whole Genome Sequencing Based Chinese Population Variation Profile and Reference Panel

abstract: The lack of Chinese population specific haplotype reference panel and whole genome sequencing resources has greatly hindered the genetics studies in the worlds largest population. Here we presented the NyuWa genome resource based on deep (26.2X) sequencing of 2,999 Chinese individuals, and constructed NyuWa reference panel of 5,804 haplotypes and 19.3M variants, which is the first publicly available Chinese population specific reference panel with thousands of samples. Compared with other panels, NyuWa reference panel reduces the Han Chinese imputation error rate by the range of 30% to 51%. Population structure and imputation simulation tests supported the applicability of one integrated reference panel for both northern and southern Chinese. In addition, a total of 22,504 loss-of-function variants in coding and noncoding genes were identified, including 11,493 novel variants. These results highlight the value of NyuWa genome resource to facilitate genetics research in Chinese and Asian populations.

category: genomics


title: STOMAGEN-like Regulates Stomatal Generation Causing Increased Water Use Efficiency during C4 Evolution

abstract: C4 photosynthesis optimizes plants carbon and water relations, allowing high photosynthetic rate with low stomatal conductance. Stomata have long been believed as a part of C4 syndrome. However, it remains unclear how stomata traits evolved along the path from C3 to C4. Stomatal patterning was examined in Flaveria genus, a model for studying C4 evolution. Comparative, transgenic and semi-in-vitro experiments were used to study molecular basis that underlies stomatal traits along C4 evolution. Novel results: the evolution from C3 to C4 species through intermediate species is accompanied by a stepwise rather than an abrupt change in the stomatal traits. The initial change occurs near Type II and dramatic change occurs at the C4-like species. On the road to C4, stomata become less in number but bigger in size and changes in stomatal density dominates changes in maximum stomatal conductance (gsmax). The reduction of FSTOMAGEN expression underlies altered gsmax between Flaveria species with different photosynthetic pathways and likely occurs in other C4 lineages. Our study provides insight into the pattern, mechanism and role of stomatal evolution along the road towards C4. This work highlights the stomatal traits in the current C4 evolutionary model and the co-evolution of photosynthetic pathway and stomata.

category: plant biology


title: How do stomata evolve on the road to C4 photosynthesis?


title: The evolution of stomatal traits along the trajectory towards C4 photosynthesis


title: Sucrose promotes D53 accumulation and tillering in rice

abstract:
- Shoot branching, a major component of shoot architecture, is regulated by multiple signals. Previous studies have indicated that sucrose may promote branching through suppressing the inhibitory effect of the hormone strigolactone (SL). However, the molecular mechanisms underlying this effect are unknown.
- Here we used molecular and genetic tools to identify the molecular targets underlying the antagonistic interaction between sucrose and SL.
- We showed that sucrose antagonises the suppressive action of SL on tillering in rice and on the degradation of D53, a major target of SL signalling. Sucrose inhibits the expression of D3, the orthologue of the arabidopsis F-box protein MAX2 required for SL signalling. Over-expression of D3 prevents sucrose from inhibiting D53 degradation and enabled the SL inhibition of tillering under high sucrose. Sucrose also prevents SL-induced degradation of D14, the SL receptor involved in D53 degradation. Interestingly, D14 over-expression enhances D53 protein levels and sucrose-induced tillering.
- Our results show that sucrose inhibits SL perception by targeting key components of SL signalling and, together with previous studies reporting the inhibition of SL synthesis by nitrate and phosphate, demonstrate the central role played by strigolactones in the regulation of plant architecture by nutrients.


title: Theory of temperature-dependent consumer-resource interactions

abstract: Changes in temperature affect consumer-resource interactions, which underpin the functioning of ecosystems. However, existing studies report contrasting predictions regarding the impacts of warming on biological rates and community dynamics. To improve prediction accuracy and comparability, we develop an approach that combines sensitivity analysis and aggregate parameters. The former determines which biological parameters impact the community most strongly. The use of aggregate parameters (i.e., maximal energetic efficiency,{rho} , and interaction strength,{kappa} ), that combine multiple biological parameters, increases explanatory power and reduces the complexity of theoretical analyses. We illustrate the approach using empirically-derived thermal dependence curves of biological rates and applying it to consumer-resource biomass ratio and community stability. Based on our analyses, we generate four predictions: 1) resource growth rate regulates biomass distributions at mild temperatures, 2) interaction strength alone determines the thermal boundaries of the community, 3) warming destabilises dynamics at low and mild temperatures only, 4) interactions strength must decrease faster than maximal energetic efficiency for warming to stabilise dynamics. We argue for the potential benefits of directly working with the aggregate parameters to increase the accuracy of predictions on warming impacts on food webs and promote cross-system comparisons.


title: Sustained ErbB activation causes demyelination and hypomyelination by driving necroptosis of mature oligodendrocytes and apoptosis of oligodendrocyte precursor cells

abstract:
Oligodendrocytes are vulnerable to genetic and environmental insults and its injury leads to demyelinating diseases. The roles of ErbB receptors in maintaining the CNS myelin integrity are largely unknown. Here we overactivate ErbB receptors that mediate signaling of either neuregulin or EGF family growth factors and found their synergistic activation caused deleterious outcomes in white matter. Sustained ErbB activation induced by the tetracycline-dependent mouse tool Plp-tTA resulted in demyelination, axonal degeneration, oligodendrocyte precursor cell (OPC) proliferation, astrogliosis, and microgliosis in white matter. Moreover, there was hypermyelination prior to these inflammatory pathological events. In contrast, sustained ErbB activation induced by another tetracycline-dependent mouse tool Sox10+/rtTA caused hypomyelination in the corpus callosum and optic nerve, which appeared to be a developmental deficit and did not associate with OPC regeneration, astrogliosis, or microgliosis. By tracing the differentiation states of cells expressing tTA/rtTA-dependent transgene or pulse-labeled reporter proteins in vitro and in vivo, we found that Plp-tTA targeted mainly mature oligodendrocytes (MOs), whereas Sox10+/rtTA targeted OPCs and newly-formed oligodendrocytes. The distinct phenotypes of mice with ErbB overactivation induced by Plp-tTA and Sox10+/rtTA consolidated their non-overlapping targeting preferences in the oligodendrocyte lineage, and enabled us to demonstrate that ErbB overactivation in MOs induced necroptosis that caused inflammatory demyelination, whereas in OPCs induced apoptosis that caused non-inflammatory hypomyelination. Early interference with aberrant ErbB activation ceased oligodendrocyte deaths and restored myelin development in both mice. This study suggests that aberrant ErbB activation is an upstream pathogenetic mechanism of demyelinating diseases, providing a potential therapeutic target.

Significance statementPrimary oligodendropathy is one of the etiological mechanisms for multiple sclerosis, and oligodendrocyte necroptosis is a pathological hallmark in the disease. Moreover, the demyelinating disease is now a broad concept that embraces schizophrenia, in which white matter lesions are an emerging feature. ErbB overactivation has been implicated in schizophrenia by genetic analysis and postmortem studies. This study suggests the etiological implications of ErbB overactivation in myelin pathogenesis and elucidates the pathogenetic mechanisms.


title: Analysis of Hippocampal Subfields in Sickle Cell Disease Using Ultrahigh Field MRI

abstract: Sickle cell disease (SCD) is an inherited hemoglobinopathy that causes organ dysfunction, including cerebral vasculopathy and neurological complications. Hippocampal segmentation with newer and advanced 7 Tesla (7T) MRI protocols has revealed atrophy in specific subregions in other neurodegenerative and neuroinflammatory diseases, however, there is limited evidence of hippocampal involvement in SCD. Thus, we explored whether SCD may be also associated with abnormalities in hippocampal subregions. We conducted 7T MRI imaging in individuals with SCD, including the HbSS, HbSC and HbS/beta thalassemia genotypes (n=37), and healthy race and age-matched controls (n=40), using a customized head coil. Both T1 and T2 weighted images were used for automatic segmentation of the hippocampus subfields. Individuals with SCD had significantly smaller volume of the Dentate Gyrus and Cornu Ammonis (CA) 2 and 3 as compared to the control group. Other hippocampal subregions also showed a trend towards smaller volumes in the SCD group. These findings support previous reports of reduced volume in the temporal lobe in SCD patients. Further studies are necessary to investigate the mechanisms that lead to structural changes in the hippocampus subfields and their relationship with cognitive performance in SCD patients.


title: Hippocampal-cortical encoding activity predicts the precision of episodic memory

abstract: Our recollections of past experiences can vary both in the number of specific event details accessible from memory and the precision with which such details are reconstructed. Prior neuroimaging evidence suggests the success and precision of episodic recollection to rely on distinct neural substrates during memory retrieval. In contrast, the specific encoding mechanisms supporting later memory precision, and whether they differ from those underlying successful memory formation in general, are currently unknown. Here, we combined continuous measures of memory retrieval with model-based analyses of behavioural and neuroimaging data to tease apart the encoding correlates of successful memory formation and mnemonic precision. In the MRI scanner, participants encoded object-scene displays, and later reconstructed features of studied objects using a continuous scale. We observed overlapping encoding activity in inferior prefrontal and posterior perceptual regions to predict both which object features were later remembered versus forgotten, and the precision with which they were reconstructed from memory. In contrast, hippocampal encoding activity significantly predicted the precision, but not overall success, of subsequent memory retrieval. The current results identify a hippocampal-cortical encoding basis for episodic memory precision, and suggest a contribution of shared cortical encoding mechanisms to the formation of both accessible and precise memory representations.


title: Sustained pupil responses are modulated by predictability of auditory sequences

abstract:
The brain is highly sensitive to auditory regularities and exploits the predictable order of sounds in many situations, from parsing complex auditory scenes, to the acquisition of language. To understand the impact of stimulus predictability on perception, it is important to determine how the detection of predictable structure influences processing and attention. Here we use pupillometry to gain insight into the effect of sensory regularity on arousal. Pupillometry is a commonly used measure of salience and processing effort, with more perceptually salient or perceptually demanding stimuli consistently associated with larger pupil diameters.

In two experiments we tracked human listeners pupil dynamics while they listened to sequences of 50ms tone pips of different frequencies. The order of the tone pips was either random, contained deterministic (fully predictable) regularities (experiment 1, n = 18, 11 female) or had a probabilistic regularity structure (experiment 2, n = 20, 17 female). The sequences were rapid, preventing conscious tracking of sequence structure thus allowing us to focus on the automatic extraction of different types of regularities. We hypothesized that if regularity facilitates processing by reducing processing demands, a smaller pupil diameter would be seen in response to regular relative to random patterns. Conversely, if regularity is associated with heightened arousal and attention (i.e. engages processing resources) the opposite pattern would be expected. In both experiments we observed a smaller sustained (tonic) pupil diameter for regular compared with random sequences, consistent with the former hypothesis and confirming that predictability facilitates sequence processing.

Significance statementThe brain is highly sensitive to auditory regularities. To appreciate the impact that detecting predictability has on perception, we need to better understand how a predictable structure influences processing and attention. We recorded listeners pupil responses to sequences of tones that followed either a predictable or unpredictable pattern, as the pupil can be used to implicitly tap into these different cognitive processes. We found that the pupil showed a smaller sustained diameter to predictable sequences, indicating that predictability eased processing rather than boosted attention. The findings suggest that the pupil response can be used to study the automatic extraction of regularities, and that the effects are most consistent with predictability helping the listener to efficiently process upcoming sounds.


title: SANS (USH1G) regulates pre-mRNA splicing by mediating the intra-nuclear transfer of tri-snRNP complexes

abstract: Splicing is catalyzed by the spliceosome, a compositionally dynamic complex assembled stepwise on pre-mRNA. We reveal links between splicing machinery components with the intrinsically disordered ciliopathy protein SANS. Pathogenic mutations in SANS/USH1G lead to Usher syndrome - the most common cause of deaf-blindness. SANS functions have been associated with cytoplasmic processes so far. Here, we reveal molecular links between SANS and pre-mRNA splicing catalyzed by the spliceosome, a compositionally dynamic complex assembled stepwise on pre-mRNA in the nucleus. Here, we show SANS associates with Cajal bodies and nuclear speckles, where SANS interacts with components of spliceosomal sub-complexes such as SF3B1 and the large splicing cofactor SON but also with PRPFs and snRNAs related to the tri-snRNP complex. SANS is required for the transfer of tri-snRNPs from Cajal bodies to nuclear speckles. SANS depletion alters the kinetics of spliceosome assembly, leading to accumulation of Complex A. SANS deficiency and USH1G pathogenic mutations affects splicing of genes related to cell proliferation and USH. Thus, we provide the first evidence that splicing deregulation may participate in the pathophysiology of Usher syndrome.


title: MicroRNA-24-3p promotes skeletal muscle differentiation and regeneration by regulating HMGA1

abstract: Numerous studies have established the critical roles of microRNAs in regulating post-transcriptional gene expression in diverse biological processes. Here, we report on the role and mechanism of miR-24-3p in skeletal muscle differentiation and regeneration. miR-24-3p promotes myoblast differentiation and skeletal muscle regeneration by directly targeting high mobility group AT-hook 1 (HMGA1) and regulating it and its direct downstream target, the inhibitor of differentiation 3 (ID3). miR-24-3p knockdown in neonatal mice increases PAX7-positive proliferating muscle stem cells (MuSCs) by derepressing Hmga1 and Id3. Similarly, inhibition of miR-24-3p in the tibialis anterior muscle prevents Hmga1 and Id3 downregulation and impairs regeneration. These findings provide evidence that the miR-24-3p/HMGA1/ID3 axis is required for MuSC differentiation and skeletal muscle regeneration in vivo.

category: developmental biology


title: Conflict over fertilization underlies the transient evolution of reinforcement

abstract: When two species meet in secondary contact, the production of low fitness hybrids may be prevented by the adaptive evolution of increased prezygotic isolation, a process known as reinforcement. Theoretical challenges to the evolution of reinforcement are generally cast as a coordination problem--e.g., recombination may break the association between trait and preference loci. However, the evolution of reinforcement also poses a potential conflict between mates. For example, the opportunity costs to hybridization may differ between the sexes or species. This is particularly likely for reinforcement based on postmating prezygotic (PMPZ) incompatibilities, as the ability to fertilize both conspecific and heterospecific eggs is beneficial to male gametes, but heterospecific mating may incur a cost for female gametes. We develop a population genetic model of interspecific conflict over reinforcement inspired by "gametophytic factors", which act as PMPZ barriers among Zea mays subspecies. We demonstrate that this conflict results in the transient evolution of reinforcement--after females adaptively evolve to reject gametes lacking a signal common in conspecific gametes, this gamete signal adaptively introgresses into the other population. Ultimately the male gamete signal fixes in both species, and prezygotic isolation returns to pre-reinforcement levels. We interpret geographic patterns of isolation among Z. mays subspecies in light of these findings, and suggest when and how this conflict can be resolved. Our results suggest that sexual conflict over fertilization may pose an understudied obstacle to the evolution of reinforcement via PMPZ incompatibilities.

category: evolutionary biology


title: The SpliZ generalizes "Percent Spliced In" to reveal regulated splicing at single-cell resolution

abstract: To date, detecting robust single-cell-regulated splicing is viewed as out of reach from droplet based technologies such as 10x Chromium. This prevents the discovery of single-cell-regulated splicing in rare cell types or those that are difficult or impossible to sequence deeply. Here, we introduce a novel, robust, and computationally efficient set of statistics, the Splicing Z Score (SpliZ) and SpliZVD, to detect regulated splicing in single cell RNA-seq including 10x Chromium. The SpliZ(VD) provides annotation-free detection of differentially regulated, complex alternative splicing events. The SpliZ generalizes and increases statistical power compared to the Percent Spliced In (PSI) and mathematically reduces to PSI for simple exon-skipping. We applied the SpliZ to primary human lung cells to discover hundreds of genes with new regulated cell-type-specific splicing. The SpliZ has wide application to enable biological discovery of genes predicted to have functionally significant splicing programs including those regulated in development.


title: Building blocks of protein structures - Physics meets Biology

abstract: The native state structures of globular proteins are stable and well-packed indicating that self-interactions are favored over protein-solvent interactions under folding conditions. We use this as a guiding principle to derive the geometry of the building blocks of protein structures - -helices and strands assembled into {beta}-sheets - with no adjustable parameters, no amino acid sequence information, and no chemistry. There is an almost perfect fit between the dictates of mathematics and physics and the rules of quantum chemistry. Our theory establishes an energy landscape that channels protein evolution by providing sequence-independent platforms for elaborating sequence-dependent functional diversity. Our work highlights the vital role of discreteness in life and has implications for the creation of artificial life and on the nature of life elsewhere in the cosmos.

category: biophysics


title: Functional network plasticity of the human social brain

abstract: AO_SCPLOWBSTRACTC_SCPLOWDespite the importance of our ability to understand the thoughts and feelings of others, the social brain remains incompletely understood. Here, we studied the plasticity of social brain function in healthy adults following the targeted training of attention-mindfulness, socio-affective, and socio-cognitive skills for 9 months. All participants were followed with repeated multimodal neuroimaging and behavioral testing. Longitudinal analyses of functional networks indicated marked and specific reorganization following mental training. Socio-cognitive training resulting in an increased integration of multiple demand and default mode regions whereas attention-mindfulness resulted in their segregation. Socio-affective training resulted in an increased functional integration of ventral attention network with these regions. Changes in functional network organization were robust after varying analysis parameters, and predictive of change in behavioral markers of compassion and perspective-taking. Our results advance the understanding of the social brain, describing its intrinsic functional organization and reorganization following mental training.


title: Changing the social brain: plasticity along macro-scale axes of functional connectivity following social mental training.


title: Characterization and structural basis of a lethal mouse-adapted SARS-CoV-2

abstract:
The ongoing SARS-CoV-2 pandemic has brought an urgent need for animal models to study the pathogenicity of the virus. Herein, we generated and characterized a novel mouse-adapted SARS-CoV-2 strain, named MASCp36, that causes severe acute respiratory symptoms and mortality in standard laboratory mice. Particularly, this model exhibits age and gender related skewed distribution of mortality akin to severe COVID-19, and the 50% lethal dose (LD50) of MASCp36 was 58 PFU in 9-month-old, male BALB/c mice. Deep sequencing identified three amino acid substitutions, N501Y, Q493H, and K417N, subsequently emerged at the receptor binding domain (RBD) of MASCp36, during in vivo passaging. All three mutations in RBD significantly enhanced the binding affinity to its endogenous receptor, mouse ACE2 (mACE2). Cryo-electron microscopy (cryo-EM) analysis of human ACE2 (hACE2) or mACE2 in complex with the RBD of MASCp36 at 3.1 to 3.7 angstrom resolution elucidates molecular basis for the receptor-binding switch driven by specific amino acid substitutions. Interestingly, N501Y and Q493H enhanced the binding affinity to human ACE2 (hACE2); while triple mutations N501Y/Q493H/K417N decreased affinity to hACE2, thus led to the reduced infectivity of MASCp36 to human cells. Our study not only provides a robust platform for studying the pathogenesis of severe COVID-19 and rapid evaluation of coutermeasures against SARS-CoV-2, but also unveils the molecular mechanism for the rapid adaption and evolution of SARS-CoV-2 in human and animals.

One sentence summaryA mouse adapted SARS-CoV-2 strain that harbored specific amino acid substitutions in the RBD of S protein showed 100% mortality in aged, male BALB/c mice.


title: Nsp1 of SARS-CoV-2 Stimulates Host Translation Termination

abstract: The Nsp1 protein of SARS-CoV-2 regulates the translation of host and viral mRNAs in cells. Nsp1 inhibits host translation initiation by occluding the entry channel of the 40S ribosome subunit. The structural study of SARS-CoV-2 Nsp1-ribosomal complexes reported post-termination 80S complex containing Nsp1 and the eRF1 and ABCE1 proteins. Considering the presence of Nsp1 in the post-termination 80S ribosomal complex simultaneously with eRF1, we hypothesized that Nsp1 may be involved in translation termination. Using a cell-free translation system and reconstituted in vitro translation system, we show that Nsp1 stimulates translation termination in the stop codon recognition stage at all three stop codons. This stimulation targets the release factor 1 (eRF1) and does not affect the release factor 3 (eRF3). The activity of Nsp1 in translation termination is provided by its N-terminal domain and the minimal required part of eRF1 is NM domain. We assume that biological meaning of Nsp1 activity in translation termination is binding with the 80S ribosomes translating host mRNAs and removal them from the pool of the active ribosomes.

category: molecular biology


title: Hand2 delineates mesothelium progenitors and is reactivated in mesothelioma

abstract: The mesothelium forms epithelial membranes that line the bodies cavities and surround the internal organs. Mesothelia widely contribute to organ homeostasis and regeneration, and their dysregulation can result in congenital anomalies of the viscera, ventral wall defects, and mesothelioma tumors. Nonetheless, the embryonic ontogeny and developmental regulation of mesothelium formation has remained uncharted. Here, we combine genetic lineage tracing, in toto live imaging, and single-cell transcriptomics in zebrafish to track mesothelial progenitor origins from the lateral plate mesoderm (LPM). Our single-cell analysis uncovers a post-gastrulation gene expression signature centered on hand2 that delineates distinct progenitor populations within the forming LPM. Combining gene expression analysis and imaging of transgenic reporter zebrafish embryos, we chart the origin of mesothelial progenitors to the lateral-most, hand2-expressing LPM and confirm evolutionary conservation in mouse. Our time-lapse imaging of transgenic hand2 reporter embryos captures zebrafish mesothelium formation, documenting the coordinated cell movements that form pericardium and visceral and parietal peritoneum. We establish that the primordial germ cells migrate associated with the forming mesothelium as ventral migration boundary. Functionally, hand2 mutants fail to close the ventral mesothelium due to perturbed migration of mesothelium progenitors. Analyzing mouse and human mesothelioma tumors hypothesized to emerge from transformed mesothelium, we find de novo expression of LPM-associated transcription factors, and in particular of Hand2, indicating the re-initiation of a developmental transcriptional program in mesothelioma. Taken together, our work outlines a genetic and developmental signature of mesothelial origins centered around Hand2, contributing to our understanding of mesothelial pathologies and mesothelioma.


title: Identification of enriched hyperthermophilic microbial communities from a deep-sea hydrothermal vent chimney under electrolithoautotrophic culture conditions.

abstract: Deep-sea hydrothermal vents are extreme and complex ecosystems based on a trophic chain. We are still unsure of the identities of the first colonizers of these environments and their metabolism, but they are thought to be (hyper)thermophilic autotrophs. Here we investigate whether the electric potential observed across hydrothermal chimneys could serve as an energy source for these first colonizers. Experiments were performed in a two-chamber microbial electrochemical system inoculated with deep-sea hydrothermal chimney samples, with a cathode as sole electron donor, CO2 as sole carbon source, and nitrate, sulfate, or oxygen as electron acceptors. After a few days of culture, all three experiments showed growth of electrotrophic biofilms consuming the electrons (directly or indirectly) and producing organic compounds including acetate, glycerol, and pyruvate. Within the biofilms, the only known autotroph species retrieved were members of Archaeoglobales. Various heterotrophic phyla also grew through trophic interactions, with Thermococcales growing in all three experiments as well as other bacterial groups specific to each electron acceptor. This electrotrophic metabolism as energy source driving initial microbial colonization of conductive hydrothermal chimneys is discussed.


title: Cortical network underlying speech production during delayed auditory feedback

abstract: Hearing ones own voice is critical for fluent speech production as it allows for the detection and correction of vocalization errors in real-time. This behavior known as the auditory feedback control of speech is impaired in various neurological disorders ranging from stuttering to aphasia, however the underlying neural mechanisms are still poorly understood. Computational models of speech motor control suggest that, during speech production, the brain uses an efference copy of the motor command to generate an internal estimate of the speech output. When actual feedback differs from this internal estimate, an error signal is generated to correct the estimate and subsequent motor commands to produce intended speech. We were able to localize these neural markers using electrocorticographic recordings from neurosurgical subjects during a delayed auditory feedback (DAF) paradigm. In this task, subjects hear their voice with a time delay as they produced words and sentences (similar to an echo on a conference call), which is well known to disrupt fluency by causing slow and stutter-like speech in humans. We observed a significant response enhancement in auditory cortex that scaled with the duration of feedback delay indicating an auditory speech error signal. Immediately following auditory cortex, dorsal precentral gyrus (dPreCG), a region that has not been implicated in auditory feedback processing before, exhibited a markedly similar response enhancement suggesting a tight coupling between the two regions. Critically, response enhancement in dPreCG occurred only when subjects profoundly slowed down their speech during articulation of long utterances due to a continuous mismatch between produced speech and reafferent feedback. These results suggest that dPreCG plays an essential role in updating the internal speech estimates to maintain fluency as well as coordinating the efference copy and auditory error signals during speech production.


title: Structure and variability of optogenetic responses identify the operating regime of cortex

abstract: Predicting the response of the cortical microcircuit to perturbations is a prerequisite to determine the mechanisms that mediate its response to stimulus; yet, an encompassing perspective that describes the full ensemble of the networks response in models that accurately recapitulate recorded data is still lacking. Here we develop a class of mathematically tractable models that exactly describe the modulation of the distribution of cell-type-specific calcium-imaging activity with the contrast of a visual stimulus. The inferred parameters recover signatures of the connectivity structure found in mouse visual cortex. Analysis of this structure subsequently reveals parameter-independent relations between the responses of different cell types to perturbations and each interneurons role in circuit-stabilization. Leveraging recent theoretical approaches, we derive explicit expressions for the distribution of responses to partial perturbations which reveal a novel, counter intuitive effect in the sign of response functions. Finally applying the theory to inferring feedback to V1 during locomotion, we find that it is predominantly mediated by both SOM and VIP modulation.


title: Chromosome-length genome assembly and structural variations of the primal Basenji dog (Canis lupus familiaris) genome

abstract:
BackgroundBasenjis are considered an ancient dog breed of central African origins that still live and hunt with tribesmen in the African Congo. Nicknamed the barkless dog, Basenjis possess unique phylogeny, geographical origins and traits, making their genome structure of great interest. The increasing number of available canid reference genomes allows us to examine the impact the choice of reference genome makes with regard to reference genome quality and breed relatedness.

ResultsHere, we report two high quality de novo Basenji genome assemblies: a female, China (CanFam_Bas), and a male, Wags. We conduct pairwise comparisons and report structural variations between assembled genomes of three dog breeds: Basenji (CanFam_Bas), Boxer (CanFam3.1) and German Shepherd Dog (GSD) (CanFam_GSD). CanFam_Bas is superior to CanFam3.1 in terms of genome contiguity and comparable overall to the high quality CanFam_GSD assembly. By aligning short read data from 58 representative dog breeds to three reference genomes, we demonstrate how the choice of reference genome significantly impacts both read mapping and variant detection.

ConclusionsThe growing number of high-quality canid reference genomes means the choice of reference genome is an increasingly critical decision in subsequent canid variant analyses. The basal position of the Basenji makes it suitable for variant analysis for targeted applications of specific dog breeds. However, we believe more comprehensive analyses across the entire family of canids is more suited to a pangenome approach. Collectively this work highlights the importance the choice of reference genome makes in all variation studies.


title: FcγR responses to soluble immune complexes of varying size: A scalable cell-based reporter system

abstract:
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.


title: FcγR responses to soluble immune complexes are governed by solubility and size


title: Immune complex solubility and size govern individual Fc-gamma receptor responses: Novel biomarkers for autoimmune diseases


title: Immune complex solubility and size govern Fc-gamma receptor responses


title: Immune complex solubility and size govern Fc-gamma receptor responses: A scalable cell-based reporter system


title: Detection and functional resolution of soluble multimeric immune complexes by a comprehensive FcγR reporter cell panel


title: snpXplorer: a web application to explore SNP-associations and annotate SNP-sets

abstract:
Genetic association studies are frequently used to study the genetic basis of numerous human phenotypes. However, the rapid interrogation of how well a certain genomic region associates across traits as well as the interpretation of genetic associations is often complex and requires the integration of multiple sources of annotation, which involves advanced bioinformatic skills. We developed snpXplorer, an easy-to-use web-server application for exploring Single Nucleotide Polymorphisms (SNP) association statistics and to functionally annotate sets of SNPs. snpXplorer can superimpose association statistics from multiple studies, and displays regional information including SNP associations, structural variations, recombination rates, eQTL, linkage disequilibrium patterns, genes and gene-expressions per tissue. By overlaying multiple GWAS studies, snpXplorer can be used to compare levels of association across different traits, which may help the interpretation of variant consequences. Given a list of SNPs, snpXplorer can also be used to perform variant-to-gene mapping and gene-set enrichment analysis to identify molecular pathways that are overrepresented in the list of input SNPs. snpXplorer is freely available at https://snpxplorer.net. Source code, documentation, example files and tutorial videos are available within the Help section of snpXplorer and at https://github.com/TesiNicco/snpXplorer.

Key points: O_LIsnpXplorer shows GWAS summary statistics, regional information and helps deciphering GWAS outcomes
C_LIO_LIsnpXplorer interactively compares association levels of a genomic region across phenotypes
C_LIO_LIsnpXplorer performs variant-to-gene mapping and gene-set enrichment analysis
C_LI

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=151 SRC="FIGDIR/small/377879v2_ufig1.gif" ALT="Figure 1">
View larger version (47K):
org.highwire.dtl.DTLVardef@13a9ba5org.highwire.dtl.DTLVardef@c07b9aorg.highwire.dtl.DTLVardef@f2ce4forg.highwire.dtl.DTLVardef@c6d228_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Functional correlates of immediate early gene expression in mouse visual cortex

abstract: During visual development, response properties of layer 2/3 neurons in visual cortex are shaped by experience. Both visual and visuomotor experience are necessary to coordinate the integration of bottom-up visual input and top-down motor-related input. Whether visual and visuomotor experience engage different plasticity mechanisms, possibly associated with the two separate input pathways, is still unclear. To begin addressing this, we measured the expression level of three different immediate early genes (IEG) (c-fos, egr1 or Arc) and neuronal activity in layer 2/3 neurons of visual cortex before and after a mouses first visual exposure in life, and subsequent visuomotor learning. We found that expression levels of all three IEGs correlated positively with neuronal activity, but that first visual and first visuomotor exposure resulted in differential changes in IEG expression patterns. In addition, IEG expression levels differed depending on whether neurons exhibited primarily visually driven or motor-related activity. Neurons with strong motor-related activity preferentially expressed EGR1, while neurons that developed strong visually driven activity preferentially expressed Arc. Our findings are consistent with the interpretation that bottom-up visual input and top-down motor-related input are associated with different IEG expression patterns and hence possibly also with different plasticity pathways.


title: Dynamic Domain Specificity In Human Ventral Temporal Cortex

abstract: An enduring neuroscientific debate concerns the extent to which neural representation is restricted to neural populations specialized for particular domains of perceptual input, or distributed outside of highly selective populations as well. A critical level for this debate is the neural representation of the identity of individual images, such as individual-level face or written word recognition. Here, intracranial recordings throughout ventral temporal cortex across 17 human subjects were used to assess the spatiotemporal dynamics of individual word and face processing within and outside regions strongly selective for these categories of visual information. Individual faces and words were first discriminable primarily only in strongly selective areas, beginning at about 150 milliseconds after word or face onset, and then discriminable both within and outside strongly selective areas approximately 170 milliseconds later. Regions of ventral temporal cortex that were and were not strongly selective both contributed non-redundant information to the discrimination of individual images. These results can reconcile previous results endorsing disparate poles of the domain specificity debate by highlighting the temporally segregated contributions of different functionally defined cortical areas to individual level representations. This work supports a dynamic model of neural representation characterized by successive domain-specific and distributed processing stages.


title: THUNDER: A reference-free deconvolution method to infer cell type proportions from bulk Hi-C data

abstract: Hi-C data provide population averaged estimates of three-dimensional chromatin contacts across cell types and states in bulk samples. Effective analysis of Hi-C data entails controlling for the potential confounding factor of differential cell type proportions across heterogeneous bulk samples. We propose a novel unsupervised deconvolution method for inferring cell type composition from bulk Hi-C data, the Two-step Hi-c UNsupervised DEconvolution appRoach (THUNDER). We conducted extensive simulations to test THUNDER based on combining two published single-cell Hi-C (scHi-C) datasets. THUNDER more accurately estimates the underlying cell type proportions compared to supervised and unsupervised methods (e.g., MuSiC, TOAST, and NMF). We further demonstrate the practical utility of THUNDER to estimate cell type proportions and identify cell-type-specific interactions in Hi-C data from adult human cortex tissue samples. THUNDER will be a useful tool in adjusting for varying cell type composition in population samples, facilitating valid and more powerful downstream analysis such as differential chromatin organization studies. Additionally, THUNDER estimated contact profiles provide a useful exploratory framework to investigate cell-type-specificity of the chromatin interactome while experimental data is still rare.


title: Tracing the evolution of aneuploid cancers from multiregional sequencing

abstract: Clonal deconvolution of mutational landscapes is crucial to understand the evolutionary dynamics of cancer. Two limiting factors for clonal deconvolution that have remained unresolved are variation in tumor cell fraction (TCF) and chromosomal copy number across different samples of the same tumor. We developed a semi-supervised algorithm that tracks variant calls through multi-sample spatiotemporal tumor data. While normalizing allele frequencies based on TCF it also adjusts for copy number changes at clonal deconvolution. Absent a priori copy number data, it renders in silico copy number estimations from bulk sequences. Using published and simulated tumor sequences, we reliably segregated clonal/subclonal variants even at low sequencing depth ([~]50x). Given at least one pure tumor sample (>70% TCF), we could normalize and deconvolve paired samples down to a TCF of 40%. This renders a reliable clonal reconstruction well adapted to multi-regionally sampled solid tumor, which are often aneuploid and contaminated by non-cancer cells.


title: Neural evidence for narrative-based processing of evidence and bias in juror decision making

abstract:
Efforts to explain complex human decisions have focused on competing theories emphasizing utility and narrative mechanisms. These are difficult to distinguish using behavior alone. Both narrative and utility theories have been proposed to explain juror decisions, which are among the most consequential complex decisions made in a modern society. Here, we use patterns of brain activation derived from large neuroimaging databases to look for signatures of the cognitive processes associated with different models of juror decision making. We asked jury-eligible subjects to rate the strength of a series of criminal cases while recording the resulting patterns of brain activation. We identified cognitive processes associated with the cumulative evidence presented. Evidence accumulation correlated with multiple narrative processes, including reading and recall. Of the cognitive processes traditionally viewed as components of utility, activation patterns associated with uncertainty, but not value, were more active with stronger evidence. Independent of utility and narrative, activations linked to reasoning and relational logic also correlated with increasing evidence. Formal comparison of narrative and utility models supported a central role for narrative in integrating evidence in complex decisions.

Significance StatementThe last decade has seen a sharply increased interest in narrative as a central cognitive process in human decision making and as an important factor in the evolution of human societies. However, the roles of narrative versus utility models of decision making remain hotly debated. While available models frequently produce similar behavioral predictions, they rely on different cognitive processes and so their roles can be separated using the right neural tests. Here, we use brain imaging during mock juror decisions to show that cognitive processes associated with narrative, and to a lesser extent utility, were engaged while subjects evaluated evidence. These results support a role for narrative in complex decisions.


title: Neural support for contributions of utility and narrative processing of evidence in juror decision making


title: A model of egocentric to allocentric understanding in mammalian brains

abstract: In the mammalian brain, allocentric representations support efficient self-location and flexible navigation. A number of distinct populations of these spatial responses have been identified but no unified function has been shown to account for their emergence. Here we developed a network, trained with a simple predictive objective, that was capable of mapping egocentric information into an allocentric spatial reference frame. The prediction of visual inputs was sufficient to drive the appearance of spatial representations resembling those observed in rodents: head direction, boundary vector, and place cells, along with the recently discovered egocentric boundary cells, suggesting predictive coding as a principle for their emergence in animals. Strikingly, the network learned a solution for head direction tracking and stabilisation convergent with known biological connectivity. Moreover, like mammalian representations, responses were robust to environmental manipulations, including exposure to novel settings. In contrast to existing reinforcement learning approaches, agents equipped with this network were able to flexibly reuse learnt behaviours --adapting rapidly to unfamiliar environments. Thus, our results indicate that these representations, derived from a simple egocentric predictive framework, form an efficient basis-set for cognitive mapping.


title: Synaptic integration by co-localized G protein-coupled receptors in presynaptic terminals

abstract:
In presynaptic terminals, membrane-delimited Gi/o-mediated presynaptic inhibition is ubiquitous and acts through G{beta}{gamma} to inhibit Ca2+ entry, or directly at SNARE complexes to inhibit Ca2+-dependent synaptotagmin-SNARE complex interactions. At CA1-subicular presynaptic terminals 5-HT1B and GABAB receptors colocalize. GABAB receptors inhibit Ca2+ entry, whereas 5-HT1B receptors target SNARE complexes. We demonstrate in male and female rats that GABAB receptors receptors alter Pr, whereas 5-HT1B receptors reduce evoked cleft glutamate concentrations allowing differential inhibition of AMPA and NMDA receptor EPSCs. This reduction in cleft glutamate concentration was confirmed by imaging glutamate release using a genetic sensor (iGluSnFR).Simulations of glutamate release and postsynaptic glutamate receptor currents were made. We tested effects of changes in vesicle numbers undergoing fusion at single synapses, relative placement of fusing vesicles and postsynaptic receptors, and the rate of release of glutamate from a fusion pore. Experimental effects of Pr changes, consistent with GABAB receptor effects, were straightforwardly represented by changes in numbers of synapses. The effects of 5-HT1B receptor-mediated inhibition are well-fit by simulated modulation of the release rate of glutamate into the cleft. Colocalization of different actions of GPCRs provide synaptic integration within presynaptic terminals. Train-dependent presynaptic Ca2+ accumulation forces frequency-dependent recovery of neurotransmission during 5-HT1B receptor activation. This is consistent with competition between Ca2+-synaptotagmin and G{beta}{gamma} at SNARE complexes. Thus, stimulus trains in 5-HT1B receptor agonist unveil dynamic synaptic modulation and a sophisticated hippocampal output filter that itself is modulated by colocalized GABAB receptors which alter presynaptic Ca2+. In combination these pathways allow complex presynaptic integration.

Significance StatementTwo G protein coupled receptors colocalize at presynaptic sites, to mediate presynaptic modulation by G{beta}{gamma}, but one - a GABAB receptor inhibits Ca2+ entry while another - a 5-HT1B receptor competes with Ca2+-synaptotagmin binding to the synaptic vesicle machinery. We have investigated downstream effects of signaling and integrative properties of these receptors. Their effects are profoundly different. GABAB receptors alter Pr leaving synaptic properties unchanged, while 5-HT1B receptors fundamentally change properties of synaptic transmission, modifying AMPA receptor but sparing NMDA receptor responses. Co-activation of these receptors allows synaptic integration because of convergence of GABAB receptor alteration on Ca2+ and the effect of this altered Ca2+ signal on 5-HT1B receptor signaling. This presynaptic convergence provides a novel form of synaptic integration.


title: Synaptic integration of subquantal neurotransmission by co-localized G protein coupled receptors in presynaptic terminals


title: A "How-to" Guide for Interpreting Parameters in Habitat-Selection Analyses

abstract:
O_LIHabitat-selection analyses allow researchers to link animals to their environment via habitat-selection or step-selection functions, and are commonly used to address questions related to wildlife management and conservation efforts. Habitat-selection analyses that incorporate movement characteristics, referred to as integrated step-selection analyses, are particularly appealing because they allow modeling of both movement and habitat-selection processes.
C_LIO_LIDespite their popularity, many users struggle with interpreting parameters in habitat-selection and step-selection functions. Integrated step-selection analyses also require several additional steps to translate model parameters into a full-fledged movement model, and the mathematics supporting this approach can be challenging for many to understand.
C_LIO_LIUsing simple examples, we demonstrate how weighted distribution theory and the inhomogeneous Poisson point-process can facilitate parameter interpretation in habitat-selection analyses. Further, we provide a "how to" guide illustrating the steps required to implement integrated step-selection analyses using the amt package.
C_LIO_LIBy providing clear examples with open-source code, we hope to make habitat-selection analyses more understandable and accessible to end users.
C_LI


title: Bedshift: perturbation of genomic interval sets

abstract:
Functional genomics experiments, like ChIP-Seq or ATAC-Seq, produce results that are summarized as a region set. Many tools have been developed to analyze region sets, including computing similarity metrics to compare them. However, there is no way to objectively evaluate the effectiveness of region set similarity metrics. In this paper we present Bedshift, a command-line tool and Python API to generate new BED files by making random perturbations to an original BED file. Perturbed files have known similarity to the original file and are therefore useful to benchmark similarity metrics. To demonstrate, we used Bedshift to create an evaluation dataset of hundreds of perturbed files generated by shifting, adding, and dropping regions from a reference BED file. Then, we compared four similarity metrics: Jaccard score, coverage score, Euclidean distance, and cosine similarity. Our results highlight differences in behavior among these metrics, such as that Jaccard score is most sensitive to added or dropped regions, while coverage score is most sensitive to shifted regions. Together, we show that Bedshift is a useful tool for creating randomized region sets for a variety of uses.

AvailabilityBSD2-licensed source code and documentation can be found at https://bedshift.databio.org.


title: TIGA: Target illumination GWAS analytics

abstract: Genome wide association studies (GWAS) can reveal important genotype-phenotype associations, however, data quality and interpretability issues must be addressed. For drug discovery scientists seeking to prioritize targets based on the available evidence, these issues go beyond the single study. Here, we describe rational ranking, filtering and interpretation of inferred gene-trait associations and data aggregation across studies by leveraging existing curation and harmonization efforts. Each gene-trait association is evaluated for confidence, with scores derived solely from aggregated statistics, linking a protein-coding gene and phenotype. We propose a method for assessing confidence in gene-trait associations from evidence aggregated across studies, including a bibliometric assessment of scientific consensus based on the iCite Relative Citation Ratio, and meanRank scores, to aggregate multivariate evidence. This method, intended for drug target hypothesis generation, scoring and ranking, has been implemented as an analytical pipeline, available as open source, with public datasets of results, and a web application designed for usability by drug discovery scientists, at https://unmtid-shinyapps.net/tiga/.


title: A Markov Random Field Model for Network-based Differential Expression Analysis of Single-cell RNA-seq Data

abstract:
BackgroundRecent development of single cell sequencing technologies has made it possible to identify genes with different expression (DE) levels at the cell type level between different groups of samples. In this article, we propose to borrow information through known biological networks to increase statistical power to identify differentially expressed genes (DEGs).

ResultsWe develop MRFscRNAseq, which is based on a Markov Random Field (MRF) model to appropriately accommodate gene network information as well as dependencies among cell types to identify cell-type specific DEGs. We implement an Expectation-Maximization (EM) algorithm with mean field-like approximation to estimate model parameters and a Gibbs sampler to infer DE status. Simulation study shows that our method has better power to detect cell-type specific DEGs than conventional methods while appropriately controlling type I error rate. The usefulness of our method is demonstrated through its application to study the pathogenesis and biological processes of idiopathic pulmonary fibrosis (IPF) using a single-cell RNA-sequencing (scRNA-seq) data set, which contains 18,150 protein-coding genes across 38 cell types on lung tissues from 32 IPF patients and 28 normal controls.

ConclusionsThe proposed MRF model is implemented in the R package MRFscRNAseq available on GitHub. By utilizing gene-gene and cell-cell networks, our method increases statistical power to detect differentially expressed genes from scRNA-seq data.


title: Costs of choosiness can promote reproductive isolation in parapatry

abstract: Species often replace each other spatially along contact zones, giving scope for parapatric speciation. In models of parapatric speciation driven by assortative mating, costs of female choosiness have so far be ignored. Yet, if females encounter only a limited number of males, those that are too choosy may remain unmated, and this should create direct sexual selection against choosiness. In our individual-based model of parapatric ecological speciation, disruptive viability selection leads to divergence of an ecological trait. Additionally, female choosiness (following a  matching mating rule based on the same ecological trait) can evolve at the risk of remaining unmated, and can limit gene flow between diverging populations. In line with previous litterature, out of the contact zone, the evolution of cost-free choosiness stops at intermediate values due to indirect selection against strong choosiness. Here we show that a weak cost of choosiness, by modifying genotypic frequencies on which viability selection acts, reduces this selection pressure, thus permitting the evolution of stronger choosiness than in the absence of costs. In strong contrast to sympatric models, costs of choosiness can therefore promote reproductive isolation in parapatry.


title: ACE2-Targeting Monoclonal Antibody as Potent and Broad-Spectrum Coronavirus Blocker

abstract: The evolution of coronaviruses, such as SARS-CoV-2, makes broad-spectrum coronavirus preventional or therapeutical strategies highly sought after. Here we report a human angiotensin-converting enzyme 2 (ACE2)-targeting monoclonal antibody, 3E8, blocked the S1-subunits and pseudo-typed virus constructs from multiple coronaviruses including SARS-CoV-2, SARS-CoV-2 mutant variants (SARS-CoV-2-D614G, B.1.1.7, B.1.351, B.1.617.1 and P.1), SARS-CoV and HCoV-NL63, without markedly affecting the physiological activities of ACE2 or causing severe toxicity in ACE2 "knock-in" mice. 3E8 also blocked live SARS-CoV-2 infection in vitro and in a prophylactic mouse model of COVID-19. Cryo-EM and "alanine walk" studies revealed the key binding residues on ACE2 interacting with the CDR3 domain of 3E8 heavy chain. Although full evaluation of safety in non-human primates is necessary before clinical development of 3E8, we provided a potentially potent and "broad-spectrum" management strategy against all coronaviruses that utilize ACE2 as entry receptors and disclosed an anti-coronavirus epitope on human ACE2.


title: Neuraminidase inhibitors rewire neutrophil function in vivo in murine sepsis and ex vivo in COVID-19

abstract:
Neutrophil overstimulation plays a crucial role in tissue damage during severe infections. Neuraminidase (NEU)-mediated cleavage of surface sialic acid has been demonstrated to regulate leukocyte responses. Here, we report that antiviral NEU inhibitors constrain host NEU activity, surface sialic acid release, ROS production, and NETs released by microbial-activated human neutrophils. In vivo, treatment with Oseltamivir results in infection control and host survival in peritonitis and pneumonia models of sepsis. Single-cell RNA sequencing re-analysis of publicly data sets of respiratory tract samples from critical COVID-19 patients revealed an overexpression of NEU1 in infiltrated neutrophils. Moreover, Oseltamivir or Zanamivir treatment of whole blood cells from severe COVID-19 patients reduces host NEU-mediated shedding of cell surface sialic acid and neutrophil overactivation. These findings suggest that neuraminidase inhibitors can serve as host-directed interventions to dampen neutrophil dysfunction in severe infections.

At a GlanceIn a severe systemic inflammatory response, such as sepsis and COVID-19, neutrophils play a central role in organ damage. Thus, finding new ways to inhibit the exacerbated response of these cells is greatly needed. Here, we demonstrate that in vitro treatment of whole blood with the viral neuraminidase inhibitors Oseltamivir or Zanamivir, inhibits the activity of human neuraminidases as well as the exacerbated response of neutrophils. In experimental models of severe sepsis, oseltamivir decreased neutrophil activation and increased the survival rate of mice. Moreover, Oseltamivir or Zanamivir ex vivo treatment of whole blood cells from severe COVID-19 patients rewire neutrophil function.


title: Identifying and prioritizing potential human-infecting viruses from their genome sequences

abstract: Determining which animal viruses may be capable of infecting humans is currently intractable at the time of their discovery, precluding prioritization of high-risk viruses for early investigation and outbreak preparedness. Given the increasing use of genomics in virus discovery and the otherwise sparse knowledge of the biology of newly-discovered viruses, we developed machine learning models that identify candidate zoonoses solely using signatures of host range encoded in viral genomes. Within a dataset of 861 viral species with known zoonotic status, our approach outperformed models based on the phylogenetic relatedness of viruses to known human-infecting viruses (AUC = 0.773), distinguishing high-risk viruses within families that contain a minority of human-infecting species and identifying putatively undetected or so far unrealized zoonoses. Analyses of the underpinnings of model predictions suggested the existence of generalisable features of viral genomes that are independent of virus taxonomic relationships and that may preadapt viruses to infect humans. Our model reduced a second set of 645 animal-associated viruses that were excluded from training to 272 high and 41 very high-risk candidate zoonoses and showed significantly elevated predicted zoonotic risk in viruses from non-human primates, but not other mammalian or avian host groups. A second application showed that our models could have identified SARS-CoV-2 as a relatively high-risk coronavirus strain and that this prediction required no prior knowledge of zoonotic SARS-related coronaviruses. Genome-based zoonotic risk assessment provides a rapid, low-cost approach to enable evidence-driven virus surveillance and increases the feasibility of downstream biological and ecological characterisation of viruses.


title: Shapes and Genescapes: Mapping Multivariate Phenotype-Biological Process Associations for Craniofacial Shape

abstract: Realistic mappings of genes to morphology are inherently multivariate on both sides of the equation. The importance of coordinated gene effects on morphological phenotypes is clear from the intertwining of gene actions in signaling pathways, gene regulatory networks, and developmental processes underlying the development of shape and size. Yet, current approaches tend to focus on identifying and localizing the effects of individual genes and rarely leverage the information content of high dimensional phenotypes. Here, we explicitly model the joint effects of biologically coherent collections of genes on a multivariate trait--craniofacial shape -- in a sample of n = 1,145 mice from the Diversity Outbred (DO) experimental line. We use biological process gene ontology (GO) annotations to select skeletal and facial development gene sets and solve for the axis of shape variation that maximally covaries with gene set marker variation. We use our process-centered, multivariate genotype-phenotype (MGP) approach to determine the overall contributions to craniofacial variation of genes involved in relevant processes and how variation in different processes corresponds to multivariate axes of shape variation. Further, we compare the directions of effect in phenotype space of mutations to the primary axis of shape variation associated with broader pathways within which they are thought to function. Finally, we leverage the relationship between mutational and pathway-level effects to predict phenotypic effects beyond craniofacial shape in specific mutants. We also introduce an online application which provides users the means to customize their own process-centered craniofacial shape analyses in the DO. The process-centered approach is generally applicable to any continuously varying phenotype and thus has wide-reaching implications for complex-trait genetics.


title: Local topology, bifurcations and mutation hot-spots in proteins with SARS-CoV-2 spike protein as an example

abstract: Novel topological methods are introduced to protein research. The aim is to identify hot-spot sites where a bifurcation can change the local topology of the protein backbone. Since the shape of a protein is intimately related to its biological function, a mutation that takes place at such a bifurcation hot-spot has an enhanced capacity to change the proteins biological function. The methodology applies to any protein but it is developed with the SARS-CoV-2 spike protein as a timely example. First, topological criteria are introduced to identify and classify potential mutation hot-spot sites along the protein backbone. Then, the expected outcome of a substitution mutation is estimated for a general class of hot-spots, by a comparative analysis of the backbone segments that surround the hot-spot sites. This analysis employs the statistics of commensurable amino acid fragments in the Protein Data Bank, in combination with general stereochemical considerations. It is observed that the notorious D614G substitution of the spike protein is a good example of such a mutation hot-spot. Several topologically similar examples are then analyzed in detail, some of them are even better candidates for a mutation hot-spot than D614G. The local topology of the recently observed N501Y mutation is also inspected, and it is found that this site is prone to a different kind of local topology changing bifurcation.