title: A Chitin-Binding Protein Ultra-highly expressed in the Outer Fold of Mantle is Related to Shell Colour in Pacific Oyster Crassostrea gigas

abstract: Molluscs constitute the second largest phylum of animals in the world, and shell colour and stripes are one of their most important phenotypic characteristics. Studies on the mechanism of shell pigmentation help understand the evolutionary and ecological significance of shell colour and serve as the basis for shell colour breeding. In this study, a matched-pair design was used in comparing the black- and white-striped mantles of the same oyster. The result showed that the stripes of shell surface are corresponding to the stripes of the mantle edge. Transcriptomic analysis identified an uncharacterized protein gene (we named it as Crassostrea gigas chitin-binding protein, CgCBP) highly expressed in the black mantles. Three folds (inner fold, middle fold and outer fold) were found on the mantle edge of oyster, but only the outer fold has the same colour as the shell. Transcriptomic comparison indicated that the three folds of the mantle edge are functionally differentiated, and the CgCBP gene is ultra-highly expressed only in the outer fold. We obtained new black and white shell periostraca from the shell notching experiment and found their structural differences by scanning electron microscope. Chitin was successfully extracted and identified from the shell periostraca. Proteomic analysis revealed differences in protein composition between the two shell periostraca. Particularly, the black shell periostraca have more proteins related to melanin biosynthesis and chitin binding compared with the white ones, and melanin particles were observed in the black mantle edge using transmission electron microscope. Magnetic bead binding and Western blot experiments showed that the CgCBP protein can specifically bind to chitin and in vivo RNAi screening indicated that CgCBP knockdown can change the structure of the shell periostracum and reduce its pigmentation. All these results suggest that the outer fold of mantle may have more important roles in shell pigmentation, shell periostracum structure is functionally correlated with shell pigmentation, and the CgCBP gene ultra-highly expressed in the outer fold may influence shell pigmentation by affecting its periostracum structure.


title: Structure-guided investigations into HMG-CoA reductase as a herbicide target

abstract: HMGR (3-hydroxy-3-methylglutaryl-coenzyme A reductase) is the rate limiting enzyme of the eukaryotic mevalonate pathway successfully targeted by statins to treat hypercholesterolemia in humans. The first statins to market have also been shown to be herbicidal yet their development as herbicides has been hindered by the risk of off-target effects. Here we present the crystal structure of a HMGR from Arabidopsis thaliana (AtHMG1) that exhibits a wider active site than all know structures from non-plant species. This enabled the rational design of statin derivatives with plant-specificity. By sequence analysis of natural statin fungal gene clusters, we could also engineer statin-resistant transgenic plants. These results suggest HMGR is a viable herbicide target, modifiable to provide a tolerance trait.

category: plant biology


title: An rRNA fragment in extracellular vesicles secreted by human airway epithelial cells increases the fluoroquinolone sensitivity of P. aeruginosa

abstract:
Lung infection by antibiotic resistant strains of Pseudomonas aeruginosa is a well-known concern for immunocompromised hosts including people with lung diseases such as cystic fibrosis. We have previously demonstrated that extracellular vesicles (EVs) secreted by primary human airway epithelial cells (AEC) deliver miRNA let-7b-5p to P. aeruginosa where it suppresses biofilm formation and increases sensitivity to beta-lactam antibiotics. In this study we used RNA-seq to characterize the small RNA (sRNA) content of EVs secreted by AEC and demonstrate transfer of multiple distinct RNA fragments from EVs to P. aeruginosa. Bioinformatic predictions reveal that several sRNAs may target all three subunits of the fluoroquinolone efflux pump MexHI-OpmD, an effect predicted to increase antibiotic sensitivity to fluoroquinolone antibiotics. Exposure of P. aeruginosa to EVs resulted in a significant reduction in the protein levels of MexH (-48%), MexI (-50%) and OpmD (-35%). Moreover, EVs reduced planktonic growth of P. aeruginosa in the presence of the fluoroquinolone antibiotic ciprofloxacin by 20%. A mexGHI-opmD deletion mutant of P. aeruginosa phenocopied this increased sensitivity to ciprofloxacin. Finally, we found that a fragment of an 18S rRNA external transcribed spacer that was transferred to P. aeruginosa by EVs was sufficient to reduce planktonic growth of P. aeruginosa in the presence of ciprofloxacin, to reduce the minimum inhibitory concentration (MIC) of P. aeruginosa for ciprofloxacin by over 50%, and to significantly reduce protein levels of MexH and OpmD. In conclusion, an rRNA fragment secreted by AEC in EVs increases the ciprofloxacin sensitivity of P. aeruginosa by targeting and down-regulating the fluoroquinolone efflux pump MexHI-OpmD. A combination of rRNA fragments and ciprofloxacin packaged in nanoparticles or EVs may benefit patients with antibiotic-resistant P. aeruginosa infections.

Author SummaryAccording to the World Health Organization and the U.S. Centers for Disease Control and Prevention the development of antibiotic resistant strains of bacteria, including Pseudomonas aeruginosa, are a significant global threat to human health. Thus, development of new approaches to eliminate antibiotic resistant infections is required. In this study, we report that lung epithelial cells secrete extracellular vesicles (EVs) that fuse with and deliver small rRNAs to P. aeruginosa, and that the rRNAs increase the sensitivity of P. aeruginosa to the antibiotic ciprofloxacin by reducing protein levels of the drug efflux pump MexHI-OpmD. We identified one rRNA fragment that by itself significantly reduced the protein levels of MexH and OpmD and increased the ability of ciprofloxacin to kill P. aeruginosa. We propose that developing synthetic vesicles containing a combination of the rRNA that inhibits antibiotic efflux pumps and ciprofloxacin would benefit patients with antibiotic resistant P. aeruginosa infections.

category: microbiology


title: The baseline immunological and hygienic status of pigs impact disease severity of African swine fever

abstract: African Swine Fever virus (ASFV) is a large double-enveloped DNA virus of the Asfarviridae family that causes a lethal hemorrhagic disease in domestic pigs and wild boars. Since 2007, a highly virulent genotype II strain has emerged and spread in Europe and South-East Asia, where millions of animals succumbed to the disease. Field- and laboratory-attenuated strains of ASFV cause highly variable clinical disease severity and survival, and mechanisms remain unclear. We hypothesized that the immunological and hygienic status of pigs is a determinant of ASF disease course. Here we compared the immunological profile at baseline and in response to ASFV infection in specific pathogen-free (SPF) and farm-raised Large White domestic pigs. At steady state, SPF pigs showed lower white blood cell counts and a lower basal inflammatory and antiviral transcriptomic profile compared to farm pigs associated with profound differences in gut microbiome composition. After inoculation with a highly virulent ASFV genotype II strain (Armenia 2008), severe clinical signs, viremia and pro-inflammatory cytokines appeared sooner in SPF pigs, indicating a reduced capacity to control early virus replication. In contrast, during infection with an attenuated field isolate (Estonia 2014), SPF pigs presented a milder and shorter clinical disease with full recovery, whereas farm pigs presented severe protracted disease with 50% lethality. Interestingly, farm pigs showed higher production of inflammatory cytokines, whereas SPF pigs produced more anti-inflammatory IL-1ra early after infection, and presented a stronger expansion of leukocytes in the recovery phase. Altogether, our data indicate that the hygiene-dependent innate immune status has a double-edge sword impact on immune responses in ASF pathogenesis. While the higher baseline innate immune activity helps the host in reducing initial virus replication, it promotes immunopathological cytokine responses, and delays lymphocyte proliferation after infection with an attenuated strain. Such effects should be considered for live vaccine development and vigilance.


title: Incorporating selection pressures into risk predictions of chemical resistance evolution in pest organisms

abstract: The evolution of chemical resistance in pests is a global threat to food security and human health, yet resistance issues are mostly dealt with after they occur. Predictive models of resistance risk provide an opportunity for field practitioners to implement proactive strategies to manage pest species more effectively. These models require identification of variables that are associated with resistance risk. Phagy and voltinism, for example, have been previously shown as useful predictors of resistance evolution. Yet despite the importance of chemical selection pressure in resistance evolution, this predictor has been difficult to incorporate into predictive models due to incomplete and (or) poor datasets on chemical usage patterns. Here we test whether chemical registrations of pesticide products can act as a proxy for selection pressures in predictive models. We obtained pesticide product registrations from publicly available chemical databases for 427 agricultural arthropod pest species in the USA and 209 Australian species. We constructed Bayesian logistic regressions of resistance status as a function of the number of pesticide product registrations for 42 chemical Mode of Action (MoA) groups, as well as phagy and voltinism. These models were well supported by our data and cross-validation analyses and could discriminate between observations of resistance and susceptibility in both USA and Australian species sets. Moreover, we found strong support for a positive correlation between pesticide product registrations and resistance status. We propose a "registration-selection-resistance" hypothesis to explain this positive association, involving a positive feedback loop between marketing and social factors of pesticide use with the evolutionary processes of selection and adaptation. Pesticide product registrations provide a predictive variable that is specific to both pest species and chemical MoA groups. Our work aids the development of predictive models for proactive resistance management, combining ecological and operational factors, with broad applicability across a suite of agricultural pests, including weeds and fungi.

category: evolutionary biology


title: Whole-Genome Profiling of Ethyl Methanesulfonate Mutagenesis in Tomato

abstract: The induced mutations accelerate crop improvement by providing novel disease resistance and yield alleles. However, the alleles with no perceptible phenotype but having an altered function remain hidden in mutagenized plants. The whole-genome sequencing (WGS) of mutagenized individuals uncovers the complete spectrum of mutations in the genome. We sequenced 132 doubly ethyl methanesulfonate (EMS)-mutagenized lines of tomato and detected ca. 41 million SNPs and 5.5 million short-INDELs. We found a very high average density of mutations 1/3.05 Kb compared to other species. About 97% of the genome had mutations, including the genes, promoters, UTRs, and introns. More than 1/3rd of genes in the mutagenized population had one or more deleterious mutations predicted by SIFT. Nearly 1/4th of deleterious genes mapped on tomato metabolic pathways modulating multiple pathway steps. Contrary to the reported GC>AT transition bias for EMS, we found EMS also produced nearly equal AT>GC transitions. Comparing mutation frequency among synonymous codons revealed that the most preferred codon is least mutagenic towards EMS. The reduction in carotenoids in {zeta}-carotene isomerase mutant fruits and chloroplasts relocation loss in phototropin1 mutant validated the mutation discovery pipeline. Our study makes a large repertoire of mutations accessible to genetic studies and the breeding of tomato.


title: Functional Recovery from Human Induced Pluripotent Stem Cell-Derived Dopamine Neuron Grafts is Dependent on Neurite Outgrowth

abstract: Transplantation of human pluripotent stem cell-derived dopaminergic (DA) neurons is an emerging therapeutic strategy for Parkinsons disease (PD). In this study, we demonstrate, for the first time, that that functional recovery after DA neuron transplant is critically dependent on graft-host integration, and not dependent on graft size or DA neuron content. Specifically, in anticipation of an autologous DA neuron transplant strategy, we studied two human induced pluripotent stem cell lines derived from people with idiopathic PD. We confirmed the cells ability to differentiate into mature DA neurons in vitro by assessing electrophysiology and depolarization-induced neurotransmitter release. To evaluate efficacy, we transplanted DA neuron precursors into a hemiparkinsonian rat model. We differentiated the cell lines for either 18 (d18) or 25 days (d25) to investigate the effect of cell maturity on efficacy. We found comparable graft survival in d18 and d25 groups from both cell lines, but behavioral analysis revealed that only d18 preparations, from both cell lines, resulted in recovery of motor impairments. Immunohistochemical analysis did not reveal any DA neuron-associated markers that correlated with efficacy. However, we found striking and consistent differences in graft neurite outgrowth between the two culture timepoints. The functional grafts from d18 cells had more outgrowth than non-functional d25 grafts. A time course of gene expression during differentiation of the cell lines revealed differences in genes associated with DA neuron development and neurite plasticity. These results are the first to demonstrate that graft-host integration, and not DA neuron content, is key to graft-induced functional recovery. The gene expression profiling may offer insight into the optimal developmental stage for graft efficacy.

category: neuroscience


title: Peer presence elicits task-independent changes within and beyond the mentalizing network across children and adults

abstract:
There is ample behavioral evidence that others mere presence can affect any behavior in human and non-human animals, generally facilitating the expression of mastered responses while impairing the acquisition of novel ones. Much less is known about i) how the brain orchestrates the modulation of such a wide array of behaviors by others presence and ii) when these neural underpinnings mature during development. To address these issues, fMRI data were collected in children and adults alternately observed and unobserved by a familiar peer. Subjects performed two tasks. One, numerosity comparison, depends on number-processing brain areas, the other, phonological comparison, on language-processing areas. Consistently with previous behavioral findings, peer observation facilitated both tasks, and childrens improvement was comparable to adults. Regarding brain activation, we found virtually no evidence of observation-driven changes within the number- or language-related areas specific to each task. Rather, we observed the same task-independent changes for both numerosity and phonological comparisons. This unique neural signature encompassed a large brain network of domain-general areas involved in social cognition, especially mentalizing, attention, and reward. It was also largely shared by children and adults. The one exception was childrens right temporoparietal junction, which failed to show the observation- driven lesser deactivation seen in adults. These findings indicate that social facilitation of some human education-related skills is i) primarily orchestrated by domain-general brain networks, rather than by task-selective substrates, and ii) relatively mature early in the course of education, thus having a protracted impact on academic achievements that may have heretofore been underestimated.

HighlightO_LIBasic math and reading skills were measured in children and adults.
C_LIO_LIParticipants were alternately observed and unobserved by a familiar peer.
C_LIO_LIBehavior showed that children were as facilitated by peer observation as adults.
C_LIO_LIfMRI data showed task-independent, observation-driven changes in mentalizing, attention, and reward brain regions.
C_LIO_LIAll adults neural changes were also found in children, except one located in right temporo- parietal junction.
C_LI


title: Mechanism of dopamine traveling waves in the striatum: theory and experiment

abstract: Striatal dopamine (DA) encodes reward, with recent work showing that DA release occurs in spatiotemporal waves. However, the mechanism of DA waves is unknown. Here we report that the striatal cholinergic neuropil also exhibits wave-like activity, and that the spatial scale of striatal DA release is extended by nicotinic receptors. Based on these findings we hypothesized that the local reciprocal interaction between cholinergic interneurons (CIN) and DA axons suffices to drive endogenous traveling waves. We show that the morphological and physiological properties of the CIN-DA interaction can be modeled as a reaction-diffusion system that gives rise to traveling waves. Analytically-tractable versions of the model show that the structure and the nature of propagation of CIN and DA traveling waves depend on their coupling, and that traveling waves can give rise to empirically observed correlations between these signals. Our model provides a biophysical mechanism for wave formation and predicts that the observed DA and CIN waves are strongly coupled phenomena.


title: Analog of Kynurenic Acid Decreases Tau Pathology by Modulating Astrogliosis in Rat Model for Tauopathy.

abstract:
Background and purposeKynurenines have immunomodulatory and neuroactive properties and can influence the central nervous system. Previous studies showed the involvement of the kynurenines in the pathogenesis and progression of neurodegenerative disease. In neurodegenerative disorders, including tauopathies, the tryptophan metabolism is shifted toward neurotoxic agents and the reduction of neuroprotectant products. Astrocyte-derived kynurenic acid serves as a neuroprotectant. However, systemic administration of kynurenic acid is not effective because of low permeability across the blood-brain barrier (BBB).

Experimental ApproachWe used a kynurenic acid analog with similar biological activity but higher brain permeability to overcome BBB limitations. In the present study, we used amide derivate of kynurenic acid N-(2-N, N-dimethylaminoethyl)-4-oxo-1H-quinoline-2-carboxamid (KYNA-1). We administered KYNA-1 for three months to tau transgenic rats SHR-24 and analyzed the effect on tau pathology and activation of glial cells. Primary glial cell cultures were applied to identify the mechanism of the KYNA-1 effect.

Key resultsKYNA-1 was not toxic to rats after chronic three-month administration. When chronically administered, KYNA-1 reduced hyperphosphorylation of insoluble tau in the brain of transgenic rats. Noteworthily, the plasma total tau was also reduced. We determined that the effect of KYNA-1 on tau pathology was induced through the modulation of glial activation. KYNA-1 inhibited LPS induced activation of astrocytes and induced transformation of microglia to M2 phenotype.

Conclusion and ImplicationsWe identified that the administration of KYNA-1 reduced tau hyperphosphorylation and neuroinflammation. KYNA-1 may serve as a promising treatment for tauopathies.

What is already known?O_LIStudies showed tryptophan-kynurenine pathway changes in neurodegenerative disorders including tauopathies
C_LIO_LIKynurenines exert immunomodulatory and neuroactive properties and have influence on the central nervous system
C_LI

What does this study add?O_LIChronic administration of synthetic analog of kynurenic acid (KYNA-1) reduces tau phosphorylation and astrogliosis in a transgenic rat model for tauopathies
C_LIO_LIThe analog reversed LPS-induced inflammatory changes in glial cell cultures
C_LI

What is the clinical significance?O_LIAdministration of KYNA-1 analog shifted the tryptophan metabolism in the neuroprotectant direction
C_LIO_LINeuroprotective analogs KYNA-1 can serve as a new and effective potential therapeutic approach for tauopathies
C_LI


title: Cross-modal alterations of Alpha Activity do not reflect inhibition of early sensory processing: A frequency tagging study

abstract: Selective attention refers to the prioritization of relevant and suppression of irrelevant sensory input. Oscillatory alpha-band activity ([~]10 Hz) has been suggested to facilitate this process by functional inhibition of early sensory regions. The objective of the current EEG study was to compare this theory to a  gating function at later stages of the processing stream. Visual cues indicated whether an upcoming target required visual orientation or more demanding auditory pitch discrimination. To gauge inhibition, we probed the visual and auditory systems with frequency-tagging during the cue-to-target interval, flickering the fixation cross at 36 Hz and simultaneously playing amplitude-modulated white noise at 40 Hz. In line with previous work, we observed a post-cue increase in posterior alpha activity when anticipating auditory targets. Surprisingly, the amplitude of visual and auditory responses increased in anticipation of auditory targets, correlated with alpha activity. This indicates that responsiveness of the visual stream was not inhibited when attention was directed to auditory processing. Furthermore, anticipation of demanding auditory targets heightened the vigilance of the sensory system, independent of sensory modality. Our findings suggest that alpha modulation does not reflect  gain control in early sensory areas, but rather exhibits its inhibitory effects at later stages.


title: Clonal spreading of tumor-infiltrating T cells underlies the robust antitumor immune responses

abstract: The repertoire of tumor-infiltrating T cells is an emerging perspective for characterizing effective antitumor T-cell responses. Oligoclonal expansion of tumor T-cell repertoire has been evaluated; however, their association with antitumor effects is unclear. We demonstrated that the polyclonal fraction of the tumor-reactive T-cell repertoire consisting of relatively minor clones increased in tumor-bearing mice treated with anti-PD-L1 or anti-CD4 monoclonal antibody, which was correlated with antitumor effect. Meanwhile, the size of the oligoclonal fraction consisting of major clones remained unchanged. Moreover, the polyclonal fraction was enriched in progenitor exhausted T cells, which are essential for a durable antitumor response, and was more dependent on CCR7+ migratory dendritic cells, which are responsible for priming tumor-reactive T cells in the tumor-draining lymph nodes. These results suggest that the expansion of diverse tumor-reactive clones ("clonal spreading") is an important mechanism by which anti-PD-L1 and anti-CD4 treatments induce robust and durable antitumor T-cell responses.

category: immunology


title: Hidden Talents: Silent Gene Clusters Encoding Magnetic Organelle Biosynthesis in a Non-Magnetotactic Phototrophic Bacterium

abstract: Horizontal gene transfer is a powerful source of innovations in prokaryotes that can affect almost any cellular system, including microbial organelles. However, typically rapid loss of non-functional gene acquisitions obscures the mechanisms determining the fate of horizontally transferable genes. Here, we report the first discovery of a horizontally inherited gene cluster encoding biosynthesis of magnetosomes, the organelles used by magnetotactic bacteria for navigation, in a non-magnetotactic phototrophic bacterium Rhodovastum atsumiense. We show that these clusters were inactivated through transcriptional silencing and antisense RNA regulation, but retain functionality, as several genes were able to complement the orthologous deletions in a remotely related magnetotactic bacterium. The laboratory transfer of foreign magnetosome genes to R. atsumiense was found to endow the strain with magnetosome biosynthesis, but strong negative selection led to rapid loss of this trait upon subcultivation. Our results provide insight into the horizontal dissemination of gene clusters encoding complex prokaryotic organelles and illuminate the potential mechanisms of their genomic preservation in a dormant state.


title: Effects of abiotic environment on invertebrate herbivory depend on plant community context in a montane grassland

abstract: Invertebrate herbivores are important and diverse, and their abundance and impacts are expected to undergo unprecedented shifts under climate change. Yet, past studies of invertebrate herbivory have documented a wide variety of responses to changing temperature, making it challenging to predict the direction and magnitude of these shifts. One explanation for these idiosyncratic responses is that changing environmental conditions may drive concurrent changes in plant communities and herbivore traits. Thus, the impacts of changing temperature on herbivory might depend on how temperature combines and interacts with characteristics of plant communities and the herbivores that occupy them. Here, we test this hypothesis by surveying invertebrate herbivory in 220, 0.5 meter-diameter herbaceous plant communities along a 1101-meter elevational gradient. Our results suggest that increasing temperature can drive community-level herbivory via at least three overlapping mechanisms: increasing temperature directly reduced herbivory, indirectly affected herbivory by reducing phylogenetic diversity of the plant community, and indirectly affected herbivory by altering the effects of functional and phylogenetic diversity on herbivory. Consequently, increasing functional diversity of plant communities had a negative effect on herbivory, but only in colder environments while a positive effect of increasing phylogenetic diversity was observed in warmer environments. Moreover, accounting for differences among herbivore feeding guilds considerably improved model fit, because different herbivore feeding guilds varied in their response to temperature and plant community composition. Together, these results highlight the importance of considering both plant and herbivore community context in order to predict how climate change will alter invertebrate herbivory.

category: ecology


title: Spatial and social structure of rewilded laboratory mice

abstract:
As an essential biomedical model organism, house mice have been studied intensely under laboratory conditions, yet they evolved to survive and reproduce in complex and dynamic environments. There has been recent interest in the study of  rewilded mice reared in complex outdoor environments, particularly for understanding the brain and behavior. Yet little work has examined lab mouse behavior under free-living conditions. Here, we characterize the emergent spatial and social structure of replicated populations of C57BL/6J (C57) mice over 10 days in large outdoor field enclosures and compare them to populations of recently wild-derived outbred house mice under the same conditions. We observed shared aspects of space use and social structure across all trials but found that C57 societies differed from those emerging from outbred mice across multiple dimensions. Males of both genotypes rapidly established and then defended territories. Female C57 mice spent more time with other individuals and explored more space relative to all other groups. These behavioral differences resulted in C57 mice rapidly forming less stable, but more densely connected, social networks than outbred wild-derived mice. These data suggest that laboratory domestication has had larger effects on female mouse social organization than their male counterparts. Importantly, this work demonstrates that C57 mice recapitulate many, but not all, aspects of social structures generated by wild mice in outdoor conditions. Rewilding allows for tractable, replicable, and ecologically realistic approaches to studying mouse behavior and can facilitate the study of the biological basis of higher order social organization.

HIGHLIGHTSO_LIWe describe emergent spatial and social structures of rewilded C57BL/6J (C57) lab mice across replicated trials in outdoor field enclosures and compare them to wild-derived outbred mice
C_LIO_LIBoth C57 and outbred males rapidly establish and maintain territories
C_LIO_LIC57 females explore the field enclosures substantially more than any other group
C_LIO_LIWith the exception of C57 females, most mice spent the majority of their recorded time alone
C_LIO_LIThe resulting societies formed by C57 mice are less modular, more densely connected, and less stable than those formed by wild-derived outbred mice
C_LI

category: animal behavior and cognition


title: Predicting A/B compartments from histone modifications using deep learning

abstract: Genomes in 3D are folded into organizational units that can influence critical biological functions. In particular, the organization of chromatin into A and B compartments segregates its active regions from the inactive regions. Compartments, evident in Hi-C contact matrices, have been used to describe cell-type specific changes in A/B organization. However, obtaining Hi-C data for all cell and tissue types of interest is prohibitively expensive, which has limited the widespread consideration of compartment status. We present a prediction tool called Compartment prediction using Recurrent Neural Network (CoRNN) that models the relationship between the compartmental organization of the genome and histone modification enrichment. Our model predicts A/B compartments, in a cross-cell type setting, with an average area under the ROC score of 90.9%. Our cell type specific compartment predictions show high overlap with known functional elements. We investigate our predictions by systematically removing combinations of histone marks and find that H3K27ac and H3K36me3 are the most predictive marks. We then perform a detailed investigation of loci where compartment status cannot be accurately predicted from these marks.These regions represent chromatin with ambiguous compartmental status, likely due to variations in status within the population of cells. As such these ambiguous loci also show highly variable compartmental status between biological replicates in the same GM12878 cell type. Our software and trained model are publicly available at https://github.com/rsinghlab/CoRNN.

category: bioinformatics


title: fastman: A fast algorithm for visualizing GWAS results using Manhattan and Q-Q plots

abstract: Visualization of GWAS summary statistics, specifically P-values, as Manhattan plots is widespread in GWAS publications, and many popular software tools are available, such as the R package qqman. But there is substantial need for further development, such as the handling of non-human data. We provide a new R package, fastman, with major additional capabilities. It handles genomes of non-model organisms, even those at a draft stage, i.e. contigs that havent been compiled to chromosomes. Non-numeric chromosome IDs are supported. It supports plotting of other genetic scores, such as FST, D statistics, selection statistics such as PBS, or other kinds of GWAS statistics such as beta. Importantly, negative or two-tailed values are supported in this package. We implement a heuristic algorithm that drastically reduces plotting time for huge datasets without any loss of visual precision, while allowing for many different data types and missing data. We provide substantial additional flexibility in highlighting and annotation. In summary, we have developed a package fastman in R for fast and efficient visualization of GWAS results and other genomewide scores using Manhattan and Q-Q plots. The package can create plots directly from association outputs by PLINK. Alternatively, it can produce plots from any R data frame with custom columns and is equipped to handle big datasets with fast plot generation. It is available for public use on https://github.com/kaustubhad/fastman.


title: An objective criterion to evaluate sequence-similarity networks helps in dividing the protein family sequence space

abstract:
The deluge of genomic data raises various challenges for computational protein annotation. The definition of superfamilies, based on conserved folds, or of families, showing more recent homology signatures, allow a first categorization of the sequence space. However, for precise functional annotation or the identification of the unexplored parts within a family, a division into subfamilies is essential. As curators of an expert database, the Carbohydrate Active Enzymes database (CAZy), we began, more than 15 years ago, to manually define subfamilies based on phylogeny reconstruction. However, facing the increasing amount of sequence and functional data, we required more scalable and reproducible methods. The recently popularized sequence similarity networks (SSNs), allows to cope with very large families and computation of many subfamily schemes. Still, the choice of the optimal SSN subfamily scheme only relies on expert knowledge so far, without any data-driven guidance from within the network. In this study, we therefore decided to investigate several network properties to determine a criterion which can be used by curators to evaluate the quality of subfamily assignments. The performance of the closeness centrality criterion, a network property to indicate the connectedness within the network, shows high similarity to the decisions of expert curators from eight distinct protein families. Closeness centrality also suggests that in some cases multiple levels of subfamilies could be possible, depending on the granularity of the research question, while it indicates when no subfamily emerged in some family evolution. We finally used closeness centrality to create subfamilies in four families of the CAZy database, providing a finer functional annotation and highlighting subfamilies without biochemically characterized members for potential future discoveries.

Author SummaryProteins perform a lot of functions within living cells. To determining their broad function, we group similar amino-acid sequences into families as their shared ancestry argue for shared functionality. Thats what we do in the CAZy database, which covers >300 Cazbohydrate-Active enZyme families nowadays. However, we need to divide families into subfamilies to provide finer readibility into (meta)genomes and guide biochemists towards unexplored regions of the sequence space. We recently used the popularized Sequence Similarity Networks (SSN) to delineate subfamilies in the large GH16 family, but had to entirely rely on expert knowledge to evaluate and take the final decision until now, which is not scalable, not eough automated and less reproducible. To accelerate the construction of protein subfamilies from sequence similarity networks, we present here an investigation of different network properties, to use as indicators for optimal subfamily divisions. The closeness centrality criterion performed best on artificial data, and recapitulates the decisions of expert curators. We used this criterion to divide four more CAZy families into subfamilies, showed that for others no subfamilies exist.

We are therefore able to create new protein subfamilies faster and with more reliability.


title: Finding heterogeneously methylated genomicregions using ONT reads

abstract:
SummaryNanopore reads encode information on the methylation status of cytosines in CpG dinucleotides. The length of the reads makes it comparatively easy to look at patterns consisting of multiple loci; here we exploit this property to look for regions where one can define subpopulations of cells based on methylation patterns. As a benchmark we run our clustering algorithm on known imprinted genes and show that the clustering based on methylation is consistent with the phasing of the genome; we then scan chromosome 15 looking for windows corresponding to heterogeneous methylation. We can also compute the covariance of methylation across these regions while keeping into account the mixture of different types of reads.

Availabilityhttps://github.com/EmanueleRaineri/releases

Contactemanuele.raineri@cnag.crg.eu, simon.heath@cnag.crg.eu

Supplementary informationTables, figures, and some further explanations of the algorithms are available as online supplementary information.


title: cvlr: Finding heterogeneously methylated genomicregions using ONT reads


title: ISM-FLUX: single-step MINFLUX with an array detector

abstract: Single-molecule localization based on the concept of MINFLUX allows for molecular resolution imaging and tracking. However, MINFLUX has a limited field-of-view and therefore requires a precise and photon-costly pre-localization step. We propose ISM-FLUX, a localization technique that combines structured illumination with structured detection. We show through simulations that replacing the point-detector with a small (e.g. 5 x 5 pixels) single-photon camera and sequentially exciting the sample with four spatially separated doughnut-shaped beams, a localization precision between 1 and 15 nm can be obtained over a field-of-view (FOV) of more than 800 nm with 100 photons. The large FOV, induced by the camera, makes the pre-localization step obsolete. We calculate the effect of different parameters, such as the position of the doughnut beams and the number of pixels, on the localization precision. We predict that the combination of a good localization precision and the experimental simplicity of ISM-FLUX will help the wide adoption of MINFLUX based microscopy.

category: biophysics


title: Low Level Whole-Brain Radiation Enhances Theranostic Potential Of Single Domain Antibody Fragments For HER2-Positive Brain Metastases

abstract:
BackgroundSingle-domain antibody fragments (aka VHH, [~]13 kDa) are promising delivery systems for brain tumor theranostics; however, achieving efficient delivery of VHH to intracranial lesions remains challenging due to the tumor-brain barrier. Here, we evaluate low-dose whole-brain irradiation as a strategy to increase the delivery of an anti-HER2 VHH to breast cancer-derived intracranial tumors in mice.

MethodsMice with intracranial HER2-positive BT474BrM3 tumors received 10-Gy fractionated cranial irradiation and evaluated using non-invasive imaging methods. The anti-HER2 VHH 5F7 was labeled with 18F, administered intravenously to irradiated mice and controls, and PET/CT imaging was conducted at various intervals after irradiation. Tumor uptake of 18F-labeled 5F7 in irradiated and control mice was compared by PET/CT image analysis and correlated with tumor volumes. In addition, longitudinal dynamic contrast-enhanced MRI (DCE-MRI) was conducted to visualize and quantify the potential effects of radiation on tumor perfusion and permeability.

ResultsIncreased 18F-labeled 5F7 intracranial tumor uptake was observed with PET in mice that received cranial irradiation, with maximum tumor accumulation seen approximately 12 days post initial radiation treatment. No radiation-induced changes in HER2 expression were detected by Western blot, flow cytometry, or on tissue sections. DCE-MRI imaging demonstrated transiently increased tumor perfusion and permeability after irradiation, consistent with the higher tumor uptake of 18F-labeled anti-HER2 5F7 in irradiated mice.

ConclusionLow-level brain irradiation induces dynamic changes in tumor vasculature that increase the intracranial tumor delivery of an anti-HER2 VHH, which could facilitate the use of labeled sdAb to detect, monitor, and treat HER2-expressing brain metastases.

Key pointsO_LILow-level radiation enhances uptake of HER2-specific VHH in intracranial tumors.
C_LIO_LIXRT + radiolabeled VHH shows promise as a treatment strategy for breast cancer brain metastases.
C_LI

Importance of the StudyImproving the detection and treatment of brain metastases (BM) that overexpress human epidermal growth factor receptor type 2 (HER2) is an urgent medical need. Drug delivery to BM is confounded by their tumor vasculature, which is more restrictive than in GBM. Single domain antibody fragments, about one-tenth the size of antibodies, could be promising theranostic vectors for BM provided sufficient BM uptake could be achieved. In this study, we utilized longitudinal PET imaging to demonstrate that low-dose whole-brain irradiation (WBRT) significantly increased 18F-labeled HER2-specific 5F7 VHH uptake in intracranial HER2-positive tumors in mice. Combining low dose WBRT with 5F7 VHH labeled with -or {beta}-particle emitting radionuclides could provide an effective and specific targeted radiotherapeutic strategy for treating patients with HER2-expressing BM.

category: cancer biology


title: helixCAM: A Platform for Programmable Cellular Assembly in Bacteria and Human Cells

abstract: Interactions between cells are indispensable for creating structure and signaling. The ability to direct precise cell-cell interactions would be powerful for engineering tissues, understanding signaling pathways, and directing immune cell targeting. In humans, intercellular interactions are mediated by cell adhesion molecules (CAMs). However, CAMs are natively expressed by many cells and have cross-reactivity, making them unsuitable for programming specific interactions. Here, we showcase "helixCAM," a platform for engineering novel CAMs by presenting coiled-coil peptides on the cell surface. helixCAMs were able to create specific cell-cell interactions and direct patterned aggregate formation in bacteria and human cells. We built rationally designed helixCAM libraries and discovered novel high-performance helixCAM pairs. High-affinity helixCAMs were then used for multicellular engineering applications, such as spherical layering, adherent cell targeting, and surface patterning. The helixCAMs are an expandable platform for directing complex multicellular assemblies, and we foresee its utility in tissue engineering, immunology, and developmental biology.

category: synthetic biology


title: More than just hitchhikers: a survey of bacterial communities associated with diatoms originating from marine reptiles

abstract: Diatoms and bacteria are known for being the first colonizers of submerged surfaces including the skin of marine reptiles. Sea turtle carapace and skin harbour diverse prokaryotic and eukaryotic microbial taxa, including several epizoic diatoms. However, the importance of diatom-bacteria associations is hardly investigated in biofilms associated with animal hosts. This study provides a detailed inventory of diatoms, bacteria, and diatom-associated bacteria originating from several loggerhead sea turtles using a combination of metabarcoding and culturing approaches. Carapace and skin samples rbcL and 16S rRNA amplicon sequencing showed a high diversity of diatoms and bacteria, respectively. Cultures of putative epizoic and non-epizoic diatom strains contained from 18 to 101 bacterial amplicon sequence variants (ASVs) and their bacterial assemblages strongly reflected those of their source host. Diatom strains allowed for enrichment and isolation of rare-in-source bacterial families such as Marinobacteraceae, Alteromonadaceae, and Alcanivoracaceae. When accounting for phylogenetic relationships between bacterial ASVs, we observed related diatom genera might retain related microbial taxa in culture, regardless of the source environment. These data provide deeper insights into several levels of sea turtle epizoic diatom and bacterial communities, and reveal the potential of epizoic biofilms as a source of novel microbes and possibly important diatom-bacteria associations.


title: Wolbachia inhibits ovarian formation and increases blood avidity in female Aedes aegypti-a novel time-specific reproductive alteration

abstract: Wolbachia, a gram-negative endosymbiotic bacterium widespread in arthropods, is well-known for changing the reproduction of its host in ways that increase its rate of spread, but there are also costs to hosts that can reduce this. Here we investigated a novel reproductive alteration of Wolbachia wAlbB on its host Aedes aegypti, based on studies of mosquito life history traits, ovarian dissections and gene expression assays. We found that an extended period of the larval stage as well as the egg stage as previously shown can increase the proportion of Wolbachia-infected females that become infertile; the impact of Wolbachia on infertility therefore accumulated before pupation. We found that ovarian formation was blocked in infertile females, even though infected females had relatively lower Wolbachia densities than fertile infected females. Infertile females also showed a higher frequency of blood feeding following a prior blood meal, indicating that they do not enter a complete gonotrophic cycle. Treatments leading to infertility decreased expression of genes related to reproduction. Our results demonstrate effects related to the causes and consequences of infertile wAlbB-infected Ae. aegypti females with implications for Wolbachia releases and have evolutionary implications for Wolbachia infections in novel hosts.


title: Gene regulation by synergism between mRNA decay and gene transcription

abstract: It has become increasingly clear in the last few years that gene expression in eukaryotes is not a linear process starting from transcription in the nucleus to the cytoplasm, but a circular one where the mRNA level is controlled by crosstalk between nuclear transcription and cytoplasmic decay pathways. One of the possible purposes of this crosstalk is to keep the mRNA level approximately constant. This is called mRNA buffering and happens when transcription and mRNA degradation act at compensatory rates. However, if transcription and mRNA degradation act synergistically, enhanced gene expression regulation would occur. In this work we mathematically modeled the effects of RNA binding proteins (RBP) when they have positive or negative effects on mRNA synthesis and decay rates. We found that they can buffer or enhance gene expression responses depending on their respective effects on transcription and mRNA stability. Then we analyzed new and previously published genomic datasets obtained for several yeast mutants related to either transcription or mRNA decay, but they are not known to possess activity in the other process. We show that some of them, which were presumed only transcription factors (Sfp1) or only decay factors (Puf3, Upf2/3), may represent examples of RBPs that make specific synergistic crosstalk to enhance the control of the mRNA levels of their target genes by combining antagonistic effects on transcriptional and mRNA stability.

category: molecular biology


title: Deciphering the roadmap of in vivo reprogramming towards pluripotency

abstract: Differentiated cells can be converted to pluripotent stem cells (iPSCs) upon ectopic expression of transcription factors OCT4, SOX2, KLF4 and MYC (OSKM) in a process known as reprogramming. Great efforts have been made to dissect intermediate states of in vitro reprogramming and how they are affected by culture conditions, while the roadmap of in vivo reprogramming remains unexplored. Here, we use single cell RNA sequencing to capture cells undergoing reprogramming in the adult pancreas. We identify markers along the trajectory from acinar identity to pluripotency, which allow in situ visualization of the intermediate states of reprogramming. Importantly, different tissues expressing OSKM, such as pancreas, stomach and colon, share markers of intermediate reprogramming, suggesting a conserved in vivo reprogramming path. Our in vivo roadmap defines landmarks along in vivo reprogramming that could be useful for applications in tissue regeneration and cellular rejuvenation based on intermediate reprogramming states.

category: developmental biology


title: Identification and implication of tissue-enriched ligands in epithelial-endothelial crosstalk during pancreas development

abstract: Development of the pancreas is driven by an intrinsic program coordinated with signals from other cell types in the epithelial environment. These intercellular communications have been so far challenging to study because of the low concentration, localized production and diversity of the signals released. Here, we combined scRNAseq data with a computational interactomic approach to identify signals involved in the reciprocal interactions between the various cell types of the developing pancreas. This in silico approach yielded 40,607 potential ligand-target interactions between the different main pancreatic cell types. Among this vast network of interactions, we focused on three ligands potentially involved in communications between epithelial and endothelial cells. Bmp7 and Wnt7b, expressed by pancreatic epithelial cells and predicted to target endothelial cells, and Sema6d, involved in the reverse interaction. In situ hybridization confirmed the localized expression of Bmp7 in the pancreatic epithelial tip cells and of Wnt7b in the trunk cells. On the contrary, Sema6d was enriched in endothelial cells. Functional experiments on ex vivo cultured pancreatic explants indicated that tip cell-produced Bmp7 restrained development of endothelial cells. This work identified ligands with a restricted tissular and cellular distribution and highlighted the role of Bmp7 in the intercellular communications shaping vessel development during pancreas organogenesis.


title: Change in prey genotype frequency rescues predator from extinction

abstract: Indirect evolutionary rescue (IER) is a mechanism where a non-evolving population is saved from extinction in an otherwise lethal environment by evolution in an interacting population. This process has been described in a predator-prey model, where extinction of the predator is prevented by a shift in the frequency of defended towards undefended prey when reduced predator densities lower selection for defended prey. We test here how increased mortality and the initial frequencies of the prey types affect IER. Combining the analysis of model simulations and experiments with rotifers feeding on an algal population we show IER in the presence of increased predator mortality. We found that IER was dependent on the ability of the prey population to evolve as well as on the frequency of the defended prey. High initial frequencies of defended prey resulted in predator extinction despite the possibility for prey evolution, as the increase in undefended prey was delayed too much to allow predator rescue. This frequency dependency for IER was more pronounced for higher predator mortalities. Our findings can help informing the development of conservation and management strategies that consider evolutionary responses in communities to environmental changes.


title: MALE-MALE SOCIAL BONDING, COALITIONARY SUPPORT AND REPRODUCTIVE SUCCESS IN WILD GUINEA BABOONS

abstract: Male-male bonds may confer substantial fitness benefits. The adaptive value of these relationships is often attributed to coalitionary support, which aids in rank ascension and female defence, ultimately resulting in greater reproductive success. We investigated the link between male-male sociality and both coalitionary support and reproductive success in wild Guinea baboons. This species lives in a tolerant multi-level society with reproductive units comprising a male and 1-6 females at the core. Males are philopatric, form differentiated, stable, and equitable affiliative relationships ( strong bonds) with other males, and lack a clear rank hierarchy. Here, we analysed behavioural and paternity data for 30 males and 50 infants collected over four years in the Niokolo-Koba National Park, Senegal. Strongly bonded males supported each other more frequently during conflicts, but strong bonds did not promote reproductive success. Instead, males that spent less time socializing with other males were associated with a higher number of females and sired more offspring. Notably, reproductively active males still maintained bonds with other males, but adjusted their social investment in relation to life-history stage. Long-term data will be needed to test if the adaptive value of male bonding lies in longer male tenure and/or in promoting group cohesion.


title: Sex-dependent transcriptional control of cardiac electrophysiology by histone acetylation modifiers based on the GTEx database

abstract: Development of safer drugs based on epigenetic modifiers, e.g. histone deacetylase inhibitors (HDACi), requires better understanding of their effects on cardiac electrophysiology. Using RNAseq data from the genotype-tissue-expression database (GTEx), we created models that link the abundance of chromatin modifiers, such as histone acetylation enzymes (HDAC/SIRT/HATs), and the gene expression of ion channels (IC) via select cardiac transcription factors (TFs) in male and female adult human hearts (left ventricle, LV). Gene expression data (transcripts per million, TPM) from GTEx donors (21 to 70 y.o.) were filtered, normalized and transformed to Euclidian space to allow quantitative comparisons in 84 female and 158 male LVs. Sex-specific partial least-square (PLS) regression models, linking gene expression data for HDAC/SIRT/HATs to TFs and to ICs gene expression, revealed tight co-regulation of cardiac ion channels by HDAC/SIRT/HATs, with stronger clustering in the male LV. Co-regulation of genes encoding excitatory and inhibitory processes in cardiac tissue by the histone modifiers may help their predominantly net-neutral effects on cardiac electrophysiology. ATP1A1, encoding for the Na/K pump, represented an outlier - with orthogonal regulation by the histone modifiers to most of the ICs. The HDAC/SIRT/HAT effects were mediated by strong (+) TF regulators of ICs, e.g. MEF2A and TBX5, in both sexes. Furthermore, for male hearts, PLS models revealed a stronger (+)/(-) mediatory role on ICs for NKX25 and TGF1B/KLF4, respectively, while RUNX1 exhibited larger (-) TF effects on ICs in females. Male-trained PLS models of HDAC/SIRT/HAT effects on ICs underestimated the effects on some ICs in females. Insights from the GTEx dataset about the co-expression and transcriptional co-regulation of histone-modifying enzymes, transcription factors and key cardiac ion channels in a sex-specific manner can help inform safer drug design.

category: bioengineering


title: CPGL: Prediction of compound-protein interaction by integrating graph attention network with long short-term memory neural network

abstract: Recent advancements of artificial intelligence based on deep learning algorithms have made it possible to computationally predict compound-protein interaction (CPI) without conducting laboratory experiments. In this manuscript, we integrated a graph attention network (GAT) for compounds and a long short-term memory neural network (LSTM) for proteins, used end-to-end representation learning for both compounds and proteins, and proposed a deep learning algorithm, CPGL (CPI with GAT and LSTM) to optimize the feature extraction from compounds and proteins and to improve the model robustness and generalizability. CPGL demonstrated an excellent predictive performance and outperforms recently reported deep learning models. Based on 3 public CPI datasets, C.elegans, Human and BindingDB, CPGL represented 1 - 5% improvement compared to existing deep-learning models. Our method also achieves excellent results on datasets with imbalanced positive and negative proportions constructed based on the above two datasets. More importantly, using 2 label reversal datasets, GPCR and Kinase, CPGL showed superior performance compared to other existing deep learning models. The AUC were substantially improved by 15% to 50% on the Kinase dataset, indicative of the robustness and generalizability of CPGL.


title: scAEGAN: Unification of Single-Cell Genomics Data by Adversarial Learning of Latent Space Correspondences

abstract: Recent progress in Single-Cell Genomics have produced different library protocols and techniques for profiling of one or more data modalities in individual cells. Machine learning methods have separately addressed specific integration challenges (libraries, samples, paired-unpaired data modalities). We formulate an unifying data-driven methodology addressing all these challenges. To this end, we design a hybrid architecture using an autoencoder (AE) network together with adversarial learning by a cycleGAN (cGAN) network, jointly referred to as scAEGAN. The AE learns a low-dimensional embedding of each condition, whereas the cGAN learns a non-linear mapping between the AE representations. The core insight is that the AE respects each samples uniqueness, whereas the cGAN exploits the distributional data similarity in the latent space. We evaluate scAEGAN using simulated data and real datasets of a single-modality (scRNA-seq), different library preparations (Fluidigm C1, CelSeq, CelSeq2, SmartSeq), and several data modalities such as paired scRNA-seq and scATAC-seq. We find that scAEGAN outperforms Seurat3 in library integration, is more robust against data sparsity, and beats Seurat 4 in integrating paired data from the same cell. Furthermore, in predicting one data modality from another, scAEGAN outperforms Babel. We conclude scAEGAN surpasses current state-of-the-art methods across several seemingly different integration challenges.


title: SiCoDEA: a simple, fast and complete app for analyzing the effect of individual drugs and their combinations

abstract: The administration of combinations of drugs is a method widely used in the treatment of different pathologies as it can lead to an increase in the therapeutic effect and a reduction in the dose compared to the administration of the single drugs. For these reasons, it is of interest to study combinations of drugs and in particular to determine whether a specific combination has a synergistic, antagonistic or additive effect, i.e greater, less than or equal to the effect expected by the sum of the individual drugs. For this purpose, various mathematical models have been developed, which use different methods to evaluate the synergy of a combination of drugs. Most of these methods are based on the Loewe Additivity Principle (Loewe et al., 1953), which has its key step in choosing the model used for predicting the effect of individual drugs. Creating a model for this purpose and calculating its parameters, however, requires a certain level of mathematical and programming knowledge or the use of commercial software. For this purpose, therefore, we have developed an open access and easy to use app that allows to explore different models and to choose the most fitting for the specific experimental data: SiCoDEA (Single and Combined Drug Effect Analysis, https://sicodea.shinyapps.io/shiny/). The data used to test SiCoDEA comes from cell line samples treated with different drug combinations and analyzed through metabolic or viability assays. SiCoDEA is developed through a Shiny interactive and easy-to-use interface (R based). There are five models taken into consideration for the analysis of single drugs and the calculation of combination index. The first is one of the most used, that of the median-effect; while the others are different forms of the log-logistic equation, with two, three and four parameters. The purpose of SiCoDEA is, on the one hand, to provide an easy-to-use tool for analyzing drug combination data and, on the other hand, also to have a view of the various steps and to offer different results based on the model chosen. An important prerequisite in analyzing drug combinations is in fact the dose-response curve calculated for individual drugs. For this purpose, SiCoDEA allows you to view the plots of the individual drugs both to evaluate the distribution of the calculated points and therefore identify any outliers, and to view the curve of the different models taken into consideration and evaluate which one best fits the data. A table showing all the R2 values for the five different models is created with the curve plot. In addition to the type of model, it is also possible to choose between two different normalization methods, one based on the maximum or minimum value and the other on the value calculated at drug concentrations equal to zero. For the chosen options, a plot is then created that shows the trend of the combination index for the different drug combinations and consequently whether it is synergy, antagonism or additivity. Finally, it is possible to export the results in single png files or in a summary report in pdf. SiCoDEA is an open-source app among the most complete and offers more functions even than the famous CompuSyn (Chou et al., 2010), as it allows you to analyze drug curves with different models, rather than just one, and it also allows the analysis of single drug curves.


title: Estimation of microtubule-generated forces using a DNA origami nanospring

abstract: Microtubules are dynamic cytoskeletal filaments that can generate forces when polymerizing and depolymerizing. Proteins that follow growing or shortening microtubule ends and couple forces to cargo movement are important for a wide range of cellular processes. Quantifying these forces and the composition of protein complexes at dynamic microtubule ends is challenging and requires sophisticated instrumentation. Here we present an experimental approach to estimate microtubule-generated forces through the extension of a fluorescent spring-shaped DNA origami molecule. Optical readout of the spring extension enables recording of force production simultaneously with single-molecule fluorescence of proteins getting recruited to the site of force generation. DNA nanosprings enable multiplexing of force measurements and only require a fluorescence microscope and basic laboratory equipment. We validate the performance of DNA nanosprings against results obtained using optical trapping. Finally, we demonstrate the use of the nanospring to study proteins that couple microtubule growth and shortening to force generation.


title: Serotonin, Etonogestrel and breathing activity in murine Congenital Central Hypoventilation Syndrome

abstract: Congenital Central Hypoventilation Syndrome, a rare disease caused by PHOX2B mutation, is associated with absent or blunted CO2/H+ chemosensitivity due to the dysfunction of PHOX2B neurons of the retrotrapezoid nucleus. No pharmacological treatment is available. Clinical observations have reported non-systematic CO2/H+ chemosensitivity recovery under desogestrel. Here, we used a preclinical model of Congenital Central Hypoventilation Syndrome, the retrotrapezoid nucleus conditional Phox2b mutant mouse, to investigate whether etonogestrel, the active metabolite of desogestrel, led to a restoration of chemosensitivity by acting on serotonin neurons known to be sensitive to etonogestrel, or retrotrapezoid nucleus PHOX2B residual cells that persist despite the mutation. The effect of etonogestrel, alone or combined with serotonin drugs, on the respiratory rhythm of medullary-spinal cord preparations from Phox2b mutants and wildtype mice was analyzed under metabolic acidosis. c-FOS, serotonin and PHOX2B were immunodetected. Serotonin metabolic pathways were characterized by ultra-high-performance liquid chromatography. We observed etonogestrel restored chemosensitivity in Phox2b mutants in a non-systematic way. Histological differences between Phox2b mutants with restored chemosensitivity and others indicated greater activation of serotonin neurons of the raphe obscurus nucleus but no effect on retrotrapezoid nucleus PHOX2B residual cells. Pharmacology of serotonin systems modulated the respiratory effect of etonogestrel differently according to serotonin metabolic pathways. Etonogestrel induced a restoration of chemosensitivity in Phox2b mutants by acting on serotonin neurons. Our work thus highlights that the state of serotonin systems was critically important for the occurrence of an etonogestrel-restoration, an element to consider in potential therapeutic intervention in Congenital Central Hypoventilation Syndrome patients.

category: physiology


title: Autophagy modulates apical growth and development in the moss Physcomitrium patens

abstract:
Different to root hairs and pollen tubes, Physcomitrium patens apical growing protonemal cells have the singularity that they continue to undergo cell divisions as the plant develops, allowing to study autophagy in the context of a multicellular apical growing tissue coupled to development. Herein, we showed that the core autophagy machinery is present in the moss P. patens, and deeply characterized the growth and development of wild-type, atg5 and atg7 loss-of-function mutants under optimal and nutrient-deprived conditions. Our results showed that the growth of the different morphological and functional protonemata apical growing cells, chloronema and caulonema, is differentially modulated by this process. These differences depend on the protonema cell type and position along the protonemal filament, and growth condition. As a global plant response, the absence of autophagy triggers the spread of the colony through protonemata growth at the expense of a reduction in buds and gametophore development, and thus the adult gametophytic and reproductive phases. Altogether this study provides valuable information indicating that autophagy has roles during apical growth with differential responses within the cell types of the same tissue and contributes to life cycle progression and thus the development of the 2D and 3D tissues of P. patens.

HIGHLIGHTAutophagy is differentially induced in protonemal cells, and contributes to apical growth, life cycle progression, and thus the development of the 2D and 3D tissues of P. patens.


title: Extensive de novo activity stabilizes epigenetic inheritance of CG methylation in Arabidopsis transposons

abstract: Cytosine methylation within CG dinucleotides (mCG) can be epigenetically inherited over many generations. Such inheritance is thought to be mediated by a semiconservative mechanism that produces binary present/absent methylation patterns. However, we show here that in Arabidopsis thaliana h1ddm1 mutants, intermediate heterochromatic mCG is stably inherited across many generations and is quantitatively associated with transposon expression. We develop a mathematical model that estimates the rates of semiconservative maintenance failure and de novo methylation at each transposon, demonstrating that mCG can be stably inherited at any level via a dynamic balance of these activities. We find that DRM2 - the core methyltransferase of the RNA-directed DNA methylation pathway - catalyzes most of the heterochromatic de novo mCG, with de novo rates orders of magnitude higher than previously thought, whereas chromomethylases make smaller contributions. Our results demonstrate that stable epigenetic inheritance of mCG in plant heterochromatin is enabled by extensive de novo methylation.

category: genetics


title: Modular safe-harbor Transgene Insertion (MosTI) for targeted single-copy and extrachromosomal array integration in C. elegans.

abstract: Efficient and reproducible transgenesis facilitates and accelerates research using genetic model organisms. Here we describe a modular safe harbor transgene insertion (MosTI) for use in C. elegans which improves targeted insertion of single-copy transgenes by homology directed repair and targeted integration of extrachromosomal arrays by non-homologous end-joining. MosTI allows easy conversion between selection markers at insertion site and a collection of universal targeting vectors with commonly used promoters and fluorophores. Insertions are targeted at three permissive safe-harbor intergenic locations and transgenes are reproducibly expressed in somatic and germ cells. Chromosomal integration is mediated by CRISPR/Cas9, and positive selection is based on a set of split markers (unc-119, hygroR, and gfp) where only animals with chromosomal insertions are rescued, resistant to antibiotics, or fluorescent, respectively. Single-copy insertion is efficient using either constitutive or heat-shock inducible Cas9 expression (25 - 75%) and insertions can be generated from a multiplexed injection mix. Extrachromosomal array integration is also efficient (7 - 44%) at MosTI landing sites or at the endogenous unc-119 locus. We use short-read sequencing to estimate the plasmid copy numbers for eight integrated arrays (6 to 37 copies) and long-read Nanopore sequencing to determine the structure and size (5.4 Mb) of one array. Using universal targeting vectors, standardized insertion strains, and optimized protocols, it is possible to construct complex transgenic strains which should facilitate the study of increasingly complex biological problems in C. elegans.


title: 3D printed SAXS chamber for controlled in-situ dialysis and optical characterization

abstract: 3D printing changes the scope of how samples can be mounted for small angle X-ray scattering (SAXS). In this paper we present a 3D printed X-ray chamber, which allows for in-situ exchange of buffer and in-situ optical transmission spectroscopy. The chamber is made of cyclic olefin copolymers (COC), including COC X-ray windows providing ultra low SAXS background. The design integrates a membrane insert for in-situ dialysis of the 100 {micro}l sample volume against a reservoir, which enables measurements of the same sample under multiple conditions using an in-house X-ray setup equipped with a 17.4 keV molybdenum source. We demonstrate the designs capabilities by measuring reversible structural changes in lipid and polymer systems as a function of salt concentration and pH. In the same chambers optical light transmission spectroscopy was carried out measuring optical turbidity of the mesophases and local pH values using pH-responsive dyes. Microfluidic exchange and optical spectroscopy combined with in-situ X-ray scattering enables vast applications for the study of responsive materials.


title: Vγ9Vδ2 T cells are potent inhibitors of SARS-CoV-2 replication and exert effector phenotypes in COVID-19 patients

abstract:
V{gamma}9V{delta}2 T cells play a key role in the innate immune response to viral infections, including SARS-CoV-1 and 2, and are activated through butyrophilin (BTN)-3A. Here, the objectives were to: 1) characterize the effects of SARS-CoV-2 infection on the number, phenotype, and activation of V{gamma}9V{delta}2 T cells in infected patients, and 2) assess the effects of in vitro SARS-CoV-2 infection on the expression of BTN3A and its impact on the activation and response of V{gamma}9V{delta}2 T cells to an anti-BTN3A antibody. Blood V{gamma}9V{delta}2 T cells decreased in clinically mild SARS-CoV-2 infections compared to healthy volunteers (HV). This decrease was maintained up to 28 days and in the recovery period. Terminally differentiated V{gamma}9V{delta}2 T cells tend to be enriched on the day of diagnosis, 28 days after and during the recovery period compared to HV. Furthermore, these cells showed cytotoxic and inflammatory activities as shown by TNF, IFN{gamma} and CD107a/b increase following anti-BTN3A activation. Moreover, BTN3A upregulation and V{gamma}9V{delta}2 T cell infiltration were observed in a lung biopsy from a fatal SARS-CoV-2 infection, as compared to HV. In vitro, SARS-CoV-2 infection significantly increased BTN3A expression in macrophages and lung cell lines. The activation via BTN3A enhanced the anti-SARS-CoV-2 V{gamma}9V{delta}2 T cells cytotoxicity and IFN-{gamma} and TNF in SARS-CoV-2 infected patient. Increasing concentrations of anti-BTN3A were accompanied by an inhibition of viral replication. Altogether, these data suggest that V{gamma}9V{delta}2 T cells are important in the immune response against SARS-CoV-2 infection and that activation by an anti-BTN3A antibody may enhance their response.

KEY POINTSO_LISARS-CoV-2 mediates upregulation of the key receptor of V{gamma}9V{delta}2 T cells BTN3A on lung tissues and cell lines as well as monocytes
C_LIO_LIDuring SARS-CoV-2 infection, V{gamma}9V{delta}2 are differentiated and efficiently degranulate and secrete cytokines upon activation with BTN3A mAb
C_LI


title: Elucidating host cell response pathways and repurposing therapeutics for SARS-CoV-2 and other coronaviruses using gene expression profiles of chemical and genetic perturbations

abstract: COVID-19 is an ongoing pandemic that has been causing devastation across the globe for over 2 years. Although there are multiple vaccines that can prevent severe symptoms, effective COVID-19 therapeutics are still of importance. Using our proprietary in silico SMarTR engine, we screened more than 22,000 unique compounds represented by over half a million gene expression profiles to uncover compounds that can be repurposed for SARS-CoV-2 and other coronaviruses in a timely and cost-efficient manner. We then tested 13 compounds in vitro and found three with potency against SARS-CoV-2 with reasonable cytotoxicity. Bortezomib and homoharringtonine are some of the most promising hits with IC50 of 1.39 M and 0.16 M, respectively for SARS-CoV-2. Tanespimycin and homoharringtonine were effective against the common cold coronaviruses. In-depth analysis highlighted proteasome, ribosome, and heat shock pathways as key targets in modulating host responses during viral infection. Further studies of these pathways and compounds have provided novel and impactful insights into SARS-CoV-2 biology and host responses that could be further leveraged for COVID-19 therapeutics development.

category: genomics


title: Redox regulation of the SARS-CoV-2 main protease provides new opportunities for drug design

abstract: AO_SCPLOWBSTRACTC_SCPLOWBesides vaccines, the development of antiviral drugs targeting SARS-CoV-2 is critical for stopping the current COVID-19 pandemic and preventing future outbreaks. The SARS-CoV-2 main protease (Mpro), a cysteine protease with essential functions in viral replication, has been validated as an effective drug target. Here, we show that Mpro is subject to redox regulation and reversibly switches between the enzymatically active dimer and the functionally dormant monomer through redox modifications of cysteine residues. These include sulfenylation, disulfide formation between the catalytic cysteine and a proximal cysteine, and generation of an allosteric lysine-cysteine SONOS bridge that is required for structural stability under oxidative stress conditions, such as those exerted by the innate immune system. We identify homo- and heterobifunctional reagents that mimic the redox switching and possess antiviral activity. The discovered redox switches are conserved in main proteases from other coronaviruses, e.g. MERS and SARS-CoV, indicating their potential as common druggable sites.

category: biochemistry


title: Persistent cross-species SARS-CoV-2 variant infectivity predicted via comparative molecular dynamics simulation

abstract: Widespread human transmission of SARS-CoV-2 highlights the substantial public health, economic, and societal consequences of virus spillover from wildlife and also presents a repeated risk of reverse spillovers back to naive wildlife populations. We employ comparative statistical analyses of a large set of short-term molecular dynamic (MD) simulations to investigate potential human to bat (Rhinolophus macrotis) cross-species infectivity allowed by the binding of SARS-CoV-2 receptor-binding domain (RBD) to angiotensin-converting enzyme 2 (ACE2) across the bat progenitor strain and emerging human strain variants of concern (VOC). We statistically compare the dampening of atom motion during binding across protein sites upon the formation of the RBD/ACE2 binding interface using bat vs. human target receptors (i.e. bACE2 and hACE2). We report that while the bat progenitor viral strain RaTG13 shows some pre-adaption to binding hACE2, it also exhibits stronger overall affinity to bACE2. However, while the early emergent human strains and later VOCs exhibit robust binding to both hACE2 and bACE2, the delta and omicron variants exhibit evolutionary adaption of binding to hACE2. However, we conclude there is a still significant risk of mammalian cross-species infectivity of human VOCs during upcoming waves of infection as COVID-19 transitions from a pandemic to endemic status.


title: Read2Tree: scalable and accurate phylogenetic trees from raw reads

abstract: The inference of phylogenetic trees from raw sequencing reads is foundational to biology. However, state-of-the-art phylogenomics requires running complex pipelines, at significant computational and labour costs, with additional constraints in sequencing coverage, assembly and annotation quality. To overcome these challenges, we present Read2tree, which directly processes raw sequencing reads into groups of corresponding genes. In a benchmark encompassing a broad variety of datasets, our assembly-free approach was 10-100x faster than conventional approaches, and in most cases more accurate--the exception being when sequencing coverage was high and reference species very distant. To illustrate the broad applicability of the tool, we reconstructed a yeast tree of life of 435 species spanning 590 million years of evolution. Applied to Coronaviridae samples, Read2Tree accurately classified highly diverse animal samples and near-identical SARS-CoV-2 sequences on a single tree--thereby exhibiting remarkable breadth and depth. The speed, accuracy, and versatility of Read2Tree enables comparative genomics at scale.


title: An approximate line attractor in the hypothalamus that encodes an aggressive internal state

abstract: The hypothalamus plays a key role in regulating innate behaviors. It is widely believed to function as a system of  labeled lines, containing behavior-specific neurons with characteristic transcriptomic and connectomic profiles. This view however fails to explain why, although activation of estrogen receptor-1 (Esr1) expressing neurons in the ventromedial hypothalamus (VMHvl) promotes aggression, few VMHvl neurons are tuned to attack. To address this paradox, we adopted an unsupervised dynamical systems framework to analyze population activity among VMHvlEsr1 neurons during aggression. We discovered that this activity contains an "integration" dimension exhibiting slow-ramping dynamics and persistent activity that correlates with escalating aggressiveness. These dynamics are implemented as an approximate line attractor in state space. Our analysis suggests a function for VMHvl to encode the intensity of behavior-relevant motive states using line attractors. This view reconciles observational and perturbational studies of VMHvl, and reveals a new mode of neural computation in the hypothalamus.


title: An intranasal nanoparticle STING agonist has broad protective immunity against respiratory viruses and variants

abstract: Respiratory viral infections, especially Influenza (endemic) or SARS-CoV-2 (pandemic since 2020), cause morbidity and mortality worldwide. Despite remarkable progress in the development and deployment of vaccines, they are clearly impacted by the rapid emergence of viral variants. The development of an off-the-shelf, effective, safe, and low-cost drug for prophylaxis against respiratory viral infections is a major unmet medical need. Here, we developed NanoSTING, a liposomally encapsulated formulation of the endogenous STING agonist, 2-3 cGAMP, to function as an immunoantiviral. NanoSTING rapidly activates the bodys innate immune system to facilitate a broad-spectrum antiviral response against SARS-CoV-2 and influenza variants in hamsters and mice. We demonstrate that a single intranasal dose of NanoSTING can: (1) treat infections throughout the respiratory system and minimize clinical symptoms, (2) protect against highly pathogenic strains of SARS-CoV-2 (alpha and delta), (3) provide durable protection against reinfection from the same strains without the need for retreatment, (4) prevent transmission of the highly infectious SARS-CoV-2 Omicron strain, and (5) provide protection against both oseltamivir-sensitive and resistant strains of influenza. Mechanistically, administration of NanoSTING rapidly upregulated interferon-stimulated and antiviral pathways in both the nasal turbinates and lung. Our results support using NanoSTING as a thermostable, immunoantiviral with broad-spectrum antiviral properties making it appealing as a therapeutic for prophylactic or early post-exposure treatment.


title: BCG vaccination of Diversity Outbred mice induces cross-reactive antibodies to SARS-CoV-2 spike protein

abstract:
The Bacillus Calmette-Guerin (BCG) vaccine, the only vaccine against tuberculosis, induces cross-protection against pathogens unrelated to Mycobacterium, including viruses. Epidemiological studies have identified potential benefits of BCG vaccination against SARS-CoV-2 infection. While BCGs heterologous effects have been widely attributable to trained immunity, we hypothesized BCG vaccination could induce cross-reactive antibodies against the spike protein of SARS-CoV-2 Wuhan-Hu-1. The concentration of IgG reactive to SARS-CoV-2 spike protein from the sera of BCG-vaccinated, Diversity Outbred (DO) mice and C57BL/6J inbred mice was measured using ELISA. Sera from 10/15 BCG-vaccinated DO mice possessed more IgG reactive to recombinant spike protein than sera from BCG-vaccinated C57BL/6J mice and unvaccinated DO mice. Amino acid sequences common to BCG cell wall/membrane proteins and SARS-CoV-2 spike protein were identified as potential antigen candidates for future study. These results imply a humoral mechanism, influenced by genotype, by which BCG vaccination could confer immunity to SARS-CoV-2.

Graphic Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=80 SRC="FIGDIR/small/488640v1_ufig1.gif" ALT="Figure 1">
View larger version (23K):
org.highwire.dtl.DTLVardef@1ff9d16org.highwire.dtl.DTLVardef@a2302dorg.highwire.dtl.DTLVardef@8ee2borg.highwire.dtl.DTLVardef@4c7b55_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: The Spike protein of SARS-CoV-2 impairs lipid metabolism and increases susceptibility to lipotoxicity: implication for a role of Nrf2

abstract:
Background/objectivesCoronavirus disease 2019 (COVID-19) patients exhibit lipid metabolic alterations, but the mechanism remains unknown. In this study, we aimed to investigate whether the Spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) impairs lipid metabolism in host cells.

MethodsA Spike cell line in HEK293 was generated using the pcDNA vector carrying the Spike gene expression cassette. A control cell line was generated using the empty pcDNA vector. Gene expression profiles related to lipid metabolic, autophagic, and ferroptotic pathways were investigated. Palmitic acid (PA)-overload was used to assess lipotoxicity-induced necrosis.

ResultsAs compared with controls, the Spike cells showed a significant increase in lipid depositions on cell membranes as well as dysregulation of expression of a panel of molecules involved lipid metabolism, autophagy, and ferroptosis. The Spike cells showed an upregulation of nuclear factor erythroid 2-related factor 2 (Nrf2), a multifunctional transcriptional factor, in response to PA. Furthermore, the Spike cells exhibited increased necrosis in response to PA-induced lipotoxicity compared to control cells in a time- and dose-dependent manner via ferroptosis, which could be attenuated by the Nrf2 inhibitor trigonelline.

ConclusionsThe Spike protein impairs lipid metabolic and autophagic pathways in host cells, leading to increased susceptibility to lipotoxicity via ferroptosis which can be suppressed by a Nrf2 inhibitor. This data also suggests a central role of Nrf2 in Spike-induced lipid metabolic impairments.

HighlightsO_LIThe Spike protein increases lipid deposition in host cell membranes
C_LIO_LIThe Spike protein impairs lipid metabolic and autophagic pathways
C_LIO_LIThe Spike protein exaggerates PA-induced lipotoxicity in host cells via ferroptosis
C_LIO_LINrf2 inhibitor Trigonelline can mitigate the Spike protein-induced necrosis
C_LI

category: cell biology


title: Phylogeny and Metadata Network Database for Epidemiologic Surveillance

abstract: AO_SCPLOWBSTRACTC_SCPLOWThe ongoing SARS-CoV-2 pandemic has highlighted the difficulty in integrating disparate data sources for epidemiologic surveillance. To address this challenge, we have created a graph database to integrate phylogenetic trees, associated metadata, and community surveillance data for phylodynamic inference. As an example use case, we divided 22,713 SARS-CoV-2 samples into 5 groups, generated maximum likelihood trees, and inferred a potential transmission network from a forest of minimum spanning trees built on patristic distances between samples. We then used Cytoscape to visualize the resultant graphs.

category: evolutionary biology


title: Tracking the circulating SARS-CoV-2 variants in Turkey: Complete genome sequencing and molecular characterization of 1000 SARS-CoV-2 samples

abstract:
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly transmissible coronavirus and has caused a pandemic of acute respiratory disease, named  coronavirus disease 2019 (COVID-19). COVID-19 has a deep impact on public health as one of the most serious pandemics in the last century. Tracking SARS-CoV-2 is important for monitoring and assessing its evolution. This is only possible by detecting all mutations in the viral genome through genomic sequencing. Moreover, accurate detection of SARS-CoV-2 and tracking its mutations is also required for its correct diagnosis. Potential effects of mutations on the prognosis of the disease can be observed. Assignment of epidemiological lineages in an emerging pandemic requires efforts. To address this, we collected 1000 SARS-CoV-2 samples from different geographical regions in Turkey and analyze their genome comprehensively. To track the virus across Turkey we focus on 10 distinct cities in different geographic regions. Each SARS-CoV-2 genome was analyzed and named according to the nomenclature system of Nextclade and Pangolin Lineage. Furthermore, the frequency of the variations observed in 10 months was also determined by region. In this way, we have observed how the virus mutations and what kind of transmission mechanism it has. The effects of age and disease severity on lineage distribution were other considered parameters. The temporal rates of SARS-CoV-2 variants by time in Turkey were close to the global trend. This study is one of the most comprehensive whole genome analyses of SARS-CoV-2 that represents a general picture of the distribution of SARS-CoV-2 variations in Turkey in 2021.

Author SummarySince the outbreak of the COVID-19 pandemic in 2019, the viral genome of SARS-CoV-2 was analysed intensively all over the world both to detect its zoonotic origin and the emerging variants worldwide together with the variants effect on the prognosis and treatment, respectively, of the infection. Remarkable COVID-19 studies were also made in Turkey as it was in the rest of the world. To date, indeed, almost all studies on COVID-19 in Turkey either sequenced only a small number of the viral genome or analysed the viral genome which was obtained from online databases. In respect thereof, our study constitutes a milestone regarding both the huge sample size consisting of 1000 viral genomes and the widespread geographic origin of the viral genome samples. Our study provides new insights both into the SARS-CoV-2 landscape of Turkey and the transmission of the emerging viral pathogen and its interaction with its vertebrate host.


title: COVID-19 Infection and Transmission Includes Complex Sequence Diversity

abstract:
SARS-CoV-2 whole genome sequencing has played an important role in documenting the emergence of polymorphisms in the viral genome and its continuing evolution during the COVID-19 pandemic. Here we present data from over 360 patients to characterize the complex sequence diversity of individual infections identified during multiple variant surges (e.g., Alpha and Delta; requiring [≥] 80% genome coverage and [≥]100X read depth). Across our survey, we observed significantly increasing SARS-CoV-2 sequence diversity during the pandemic and frequent occurrence of multiple biallelic sequence polymorphisms in all infections. This sequence polymorphism shows that SARS-CoV-2 infections are heterogeneous mixtures. Convention for reporting microbial pathogens guides investigators to report a majority consensus sequence. In our study, we found that this approach would under-report at least 79% of the observed sequence variation. As we find that this sequence heterogeneity is efficiently transmitted from donors to recipients, our findings illustrate that infection complexity must be monitored and reported more completely to understand SARS-CoV-2 infection and transmission dynamics involving both immunocompetent and immunocompromised patients. Many of the nucleotide changes that would not be reported in a majority consensus sequence have now been observed as lineage defining SNPs in Omicron BA.1 and/or BA.2 variants. This suggests that minority alleles in earlier SARS-CoV-2 infections may play an important role in the continuing evolution of new variants of concern.

AUTHOR SUMMARYEvolution of the virus causing COVID-19 (SARS-CoV-2) has been associated with significant transmission surges. With evolution of SARS-CoV-2, evidence has accumulated regarding increased transmissibility of lineages, varying severity of illness, evasion of vaccines and diagnostic tests. Continuous tracking of SARS-CoV-2 lineage evolution distills very large and complex viral sequence data sets down to consensus sequences that report the majority nucleotide at each of over 29,000 positions in the SARS-CoV-2 genome. We observe that this eliminates considerable sequence variation and leads to a significant underestimation of SARS-CoV-2 infection diversity and transmission complexity. Additionally, concentration on the majority consensus sequence diverts attention from genetic variation that may contribute significantly to the continuing evolution of the COVID-19 pandemic.


title: Cryo-EM structures of Na+-pumping NADH-ubiquinone oxidoreductase from Vibrio cholerae

abstract: The Na+-pumping NADH-ubiquinone oxidoreductase (Na+-NQR) couples electron transfer from NADH to ubiquinone with Na+-pumping, generating an electrochemical Na+ gradient that is essential for energy-consuming reactions in bacteria. Since Na+-NQR is exclusively found in prokaryotes, it is a promising target for highly selective antibiotics. However, the molecular mechanism of inhibition is not well-understood for lack of the atomic structural information about an inhibitor-bound state. Here we present cryo-electron microscopy structures of Na+- NQR from Vibrio cholerae with or without a bound inhibitor at 2.5- to 3.1-[A] resolution. The structures reveal the arrangement of all six redox cofactors including riboflavin, whose position has been under debate, and a newly assigned 2Fe-2SNqrD/E cluster located between the membrane embedded NqrD and NqrE subunits. A large part of the hydrophilic NqrF near the cytoplasmic membrane surface is barely visible in the density map, suggesting a high degree of flexibility. This flexibility may be responsible to reducing the long distance between the 2Fe- 2S centers in NqrF and NqrD/E, consistent with physiologically relevant electron transfer. Two different types of specific inhibitors (korormicin A and aurachin D-42) bind to the N-terminal region of NqrB, which is disordered in the absence of inhibitors. The current inhibitor-bound structures reasonably explain our previous biochemical findings obtained by different chemistry-based experiments. This study provides a definite foundation for understanding the function of Na+-NQR and the molecular mechanism of its specific inhibitors to support molecular design of new antibiotics targeting the enzyme.


title: Predicting gene expression from cell morphology in human induced pluripotent stem cells

abstract: Purification is essential before differentiating human induced pluripotent stem cells (hiPSCs) into cells that fully express particular differentiation marker genes. High-quality iPSC clones are typically purified through gene expression profiling or visual inspection of the cell morphology; however, the relationship between the two methods remains unclear. We investigated the relationship between gene expression levels and morphology by analyzing live-cell phase-contrast images and mRNA profiles collected during the purification process. We employed this data and an unsupervised image feature extraction method to build a model that predicts gene expression levels from morphology. As a benchmark, we confirmed that the method can predict the gene expression levels from tissue images for cancer genes, performing as well as state-of-the-art methods. We then applied the method to iPSCs and identified two genes that are well-predicted from cell morphology. Although strong batch effects resulting from the reprogramming process preclude the ability to use the same model to predict across batches, prediction within a reprogramming batch is sufficiently robust to provide a practical approach for estimating expression levels of a few genes and monitoring the purification process.


title: Nucleosome sliding by the Chd1 chromatin remodeler relies on theintegrity of the DNA duplex

abstract: Chromatin remodelers use a helicase-type ATPase motor to shift DNA around the histone core. Although not directly reading out the DNA sequence, some chromatin remodelers are biased by DNA sequences, suggesting that they may be sensitive to properties of the DNA duplex. Here, we present a high-throughput method for determining nucleosome positioning in vitro using site-specific DNA cleavage coupled with next-generation sequencing. This method allowed us to systematically test how the introduction of poly(dA:dT) tracts and other perturbations affected the distribution of nucleosomes remodeled by the Chd1 remodeler. We found that Chd1 is sensitive to poly(dA:dT) tracts as short as 3 bp, and that its nucleosome sliding activity is severely perturbed by DNA mismatches and single-nucleotide insertions. These results suggest that remodelers rely on the integrity of duplex DNA for nucleosome sliding. We also discovered that DNA on the nucleosome can shift in the absence of a remodeler when multiple mismatches are placed at superhelix location 2 (SHL2). This DNA movement in response to a disruption of the double helix may explain why SHL2 is the preferred site of engagement by most chromatin remodelers.


title: Structures of the T cell potassium channel Kv1.3 with immunoglobulin modulators

abstract: The Kv1.3 potassium channel is expressed abundantly on activated T cells and mediates the cellular immune response. This role has made the channel a target for therapeutic immunomodulation to block its activity and suppress T cell activation. We determined structures of human Kv1.3 alone, with a nanobody inhibitor, and with an antibody-toxin fusion blocker. Rather than block the channel directly, four copies of the nanobody bind the tetramers voltage sensing domains and the pore domain to induce an inactive pore conformation. In contrast, the antibody-toxin fusion docks its toxin domain at the extracellular mouth of the channel to insert a critical lysine into the pore. The lysine induces an active conformation of the pore yet blocks ion permeation. This study visualizes Kv1.3 pore dynamics, defines two distinct mechanisms to suppress Kv1.3 channel activity with exogenous inhibitors, and provides a framework to aid development of emerging T cell immunotherapies.


title: Microbial mutualism generates multistable and oscillatory growth dynamics

abstract: Microbial communities typically comprise multiple different species with an intricate network of interactions, ranging from competitive to cooperative, between them. How does the nature of these inter-species interactions impact overall community behavior? While the influence of purely competitive interactions is well-studied, the opposite case of mutualistic interactions--which are also prevalent in many naturally-occurring communities--is poorly understood. Here, we address this gap in knowledge by mathematically modeling a well-mixed two-species community of aerobes and anaerobes having mutualistic metabolic interactions between them. Despite the simplicity of the model, we find that it reproduces three characteristic experimental findings. In particular, in response to changes in the fluxes of exogenously-supplied carbon and oxygen, the community adopts two distinct stable states with differing fractions of aerobes and anaerobes. These states are bistable, capable of arising under identical environmental conditions; transitions between the two are therefore history-dependent and can give rise to oscillations in the bacterial and chemical concentrations. Moreover, using the model, we establish biophysical principles describing how oxygen depletion and nutrient sharing jointly dictate the characteristics of the different states as well as the transitions between them. Altogether, this work thus helps disentangle and highlight the pivotal role of mutualism in governing the overall stability and functioning of microbial communities. Moreover, our model provides a foundation for future studies of more complex communities that play important roles in agriculture, environment, industry, and medicine.

category: systems biology