title: High expression of VRT2 during wheat spikelet initiation increases the number of rudimentary basal spikelets

abstract:
Spikelets are the fundamental building blocks of Poaceae inflorescences and their development and branching patterns determine the various inflorescence architectures and grain yield of grasses. In wheat, the central spikelets produce the most and largest grains, while spikelet size gradually decreases acro- and basipetally, giving rise to the characteristic lanceolate shape of wheat spikes. The acropetal gradient correlates with the developmental age of spikelets, however the basal spikelets are developed first and the cause of their small size and rudimentary development is unclear. Here, we adapted G&T-seq, a low-input transcriptomics approach, to characterise gene expression profiles within spatial sections of individual spikes before and after the establishment of the lanceolate shape. We observed larger differences in gene expression profiles between the apical, central and basal sections of a single spike than between any section belonging to consecutive developmental timepoints. We found that SVP MADS-box transcription factors, including VRT-A2, are expressed highest in the basal section of the wheat spike and display the opposite expression gradient to flowering E-class SEP1 genes. Based on multi-year field trials and transgenic lines, we show that higher expression of VRT-A2 in the basal sections of the spike is associated with increased numbers of rudimentary basal spikelets. Our results, supported by computational modelling, suggest that the delayed transition of basal spikelets from vegetative to floral developmental programmes results in the lanceolate shape of wheat spikes. This study highlights the value of spatially resolved transcriptomics to gain new insights into developmental genetics pathways of grass inflorescences.

One sentence summaryLarge transcriptional gradients exist within a wheat spike and are associated with rudimentary basal spikelet development, resulting in the characteristic lanceolate shape of wheat spikes.


title: High expression of VRT2 increases the number of rudimentary basal spikelets in wheat

category: plant biology


title: Mnemonic content and hippocampal patterns shape judgments of time

abstract:
Our experience of time can feel dilated or compressed, rather than reflecting true "clock time." Although many contextual factors influence the subjective perception of time, it is unclear how memory accessibility plays a role in constructing our experience of and memory for time. Here, we used a combination of behavioral and fMRI measures to ask the question of how memory is incorporated into temporal duration judgments. Behaviorally, we found that event boundaries, which have been shown to disrupt ongoing memory integration processes, result in the temporal compression of duration judgments. Additionally, using a multivoxel pattern similarity analysis of fMRI data, we found that greater temporal pattern change in the left hippocampus within individual trials was associated with longer duration judgments. Together, these data suggest that mnemonic processes play a role in constructing representations of time.

Statement of RelevanceOur everyday experiences convey a powerful truth: That our perception of time often diverges from the reality of time. When enjoying an active vacation with family, time moves quickly: hours go by in minutes. When sitting through an unnecessary meeting, time moves slowly: minutes go by in hours. What is the origin of these phenomenologically compelling illusions of time perception? Past research has examined how a range of specific factors, from emotions to blinking, contribute to the distortion of time. Here, in contrast, we evaluate how the content and accessibility of our memories shapes time perception. We show that context shifts, known to disrupt memory processing, also lead to robust contractions of perceived time. We discuss how both effects -- memory disruptions and time distortions -- may be linked via the hippocampus.

"An hour, once it lodges in the queer element of the human spirit, may be stretched to fifty or a hundred times its clock length; on the other hand, an hour may be accurately represented on the timepiece of the mind by one second." -- Virginia Woolf

category: neuroscience


title: Transcranial magnetic stimulation of the occipital cortex interferes with foot movements in blind individuals

abstract: Research in blind individuals has shown that after visual loss, the occipital cortex can be reorganized and repurposed for nonvisual perception and cognitive functions. However, no studies have directly examined the involvement of the visual cortex in motor function. Here, we show that a rhythmic foot movement performed by blind individuals can be disrupted by transcranial magnetic stimulation (TMS) to their primary and secondary visual cortex (V1/V2). This disruptive effect of TMS was absent for sighted participants. Our result suggests that the visual cortex of blind individuals is involved in sensorimotor control. This is the first experimental evidence that functional repurposing of the human visual cortex is not be restricted to perception and cognitive functions, but also extends to motor function.


title: Diverse but unique astrocytic phenotypes during embryonic stem cell differentiation, culturing and aging

abstract: Astrocytes are resident glia cells of the central nervous system (CNS) that play complex and heterogeneous roles in brain development, homeostasis and disease. Since their vast involvement in health and disease is becoming increasingly recognized, suitable and reliable tools for studying these cells in vivo and in vitro are of utmost importance. One of the key challenges hereby is to adequately mimic their context-dependent in vivo phenotypes and functions in vitro. To better understand the spectrum of astrocytic variations in defined settings we performed a side-by-side-comparison of embryonic stem cell (ESC)-derived astrocytes as well as primary neonatal and adult astrocytes, revealing major differences on a functional and transcriptomic level, specifically on proliferation, migration, calcium signaling and cilium activity. Our results highlight the need to carefully consider the choice of astrocyte origin and phenotype with respect to age, isolation and culture protocols based on the respective biological question.


title: Predictive framework of macrophage activation

abstract:
Macrophages populate every organ during homeostasis and disease, displaying features of tissue imprinting and heterogeneous activation. The disjointed picture of macrophage biology that emerged from these observations are a barrier for integration across models or with in vitro macrophage activation paradigms. We set out to contextualize macrophage heterogeneity across mouse tissues and inflammatory conditions, specifically aiming to define a common framework of macrophage activation. We built a predictive model with which we mapped the activation of macrophages across 12 tissues and 25 biological conditions, finding a striking commonality and finite number of transcriptional profiles, which we modelled as defined stages along four conserved activation paths. We verified this model with adoptive cell transfer experiments and identified transient RELM[a] expression as a feature of macrophage tissue engraftment. We propose that this integrative approach of macrophage classification allows the establishment of a common predictive framework of macrophage activation in inflammation and homeostasis.

One Sentence SummaryWe propose an integrative approach of macrophage classification that allows the establishment of a common framework of macrophage activation in inflammation and homeostasis.

category: immunology


title: A role of cryptochrome for magnetic field-dependent improvement of sleep quality, lifespan and motor function in Drosophila

abstract:
Understanding the molecular genetic basis of animal magnet reception has been one of the big challenges in biology. Recently it was discovered that the magnetic sense of Drosophila melanogaster is mediated by the ultraviolet (UV)-A/blue light photoreceptor cryptochrome (Cry).

Here, using the fruit fly as a magnet receptive model organism, we show that exposure to a specific AC magnetic field during night time affected the health of the fly. AC magnetic field exposure showed lifespan extension under starvation, sleep improvement and prevention of decreased motor function. In contrast, all the health improvement effect was not observed in cryptochrome mutant flies (cryb). We showed that AC magnetic field exposure prevented motor dysfunction in Gauchers disease model Drosophila. The data suggests that magnetic field-dependent improvement of sleep quality, lifespan and motor function is mediated through a cry-dependent pathway in animals.

category: animal behavior and cognition


title: Bacteria.guru: comparative transcriptomics and co-expression database for bacterial pathogens

abstract:
SummaryThe bacterial kingdom comprises unicellular prokaryotes able to establish symbioses from mutualism to parasitism. To combat bacterial pathogenicity, we need an enhanced understanding of gene function and regulation, which will mediate the development of novel antimicrobials. Gene expression can predict gene function, but there lacks a database enabling expansive inter- and intraspecific exploration of gene expression profiles and co-expression networks for bacteria. To address this, we integrated the genomic and transcriptomic data of the 17 most notorious and studied bacterial pathogens, creating bacteria.guru, an interactive database that can identify, visualize, and compare gene expression profiles, co-expression networks, functionally enriched clusters, and gene families across species. Through illustrating antibiotic resistance mechanisms in P. aeruginosa, we demonstrate that bacteria.guru could potentially aid the discovery of multi-faceted antibiotic targets. Hence, we believe bacteria.guru will facilitate future bacterial research.

AvailabilityThe database and co-expression networks are freely available from https://bacteria.guru/. The sample annotations are found in the supplemental data.

category: bioinformatics


title: Automatic segmentation of the rat brain hippocampus in MRI after traumatic brain injury

abstract: AO_SCPLOWBSTRACTC_SCPLOWRegistration-based methods are commonly used in the anatomical segmentation of magnetic resonance (MR) brain images. However, they are sensitive to the presence of deforming brain pathologies that may interfere with the alignment of the atlas image with the target image. Our goal was to develop an algorithm for automated segmentation of the normal and injured rat hippocampus. We implemented automated segmentation using a U-Net-like Convolutional Neural Network (CNN). of sham-operated experimental controls and rats with lateral-fluid-percussion induced traumatic brain injury (TBI) on MR images and trained ensembles of CNNs. Their performance was compared to three registration-based methods: single-atlas, multi-atlas based on majority voting and Similarity and Truth Estimation for Propagated Segmentations (STEPS). Then, the automatic segmentations were quantitatively evaluated using six metrics: Dice score, Hausdorff distance, precision, recall, volume similarity and compactness using cross-validation. Our CNN and multi-atlas -based segmentations provided excellent results (Dice scores > 0.90) despite the presence of brain lesions, atrophy and ventricular enlargement. In contrast, the performance of singe-atlas registration was poor (Dice scores < 0.85). Unlike registration-based methods, which performed better in segmenting the contralateral than the ipsilateral hippocampus, our CNN-based method performed equally well bilaterally. Finally, we assessed the progression of hippocampal damage after TBI by applying our automated segmentation tool. Our data show that the presence of TBI, time after TBI, and whether the location of the hippocampus was ipsilateral or contralateral to the injury explained hippocampal volume (p = 0.029, p < 0.001, and p < 0.001 respectively).


title: Automated, high-dimensional evaluation of physiological aging and resilience in outbred mice

abstract: Behavior and physiology are essential readouts in many studies but have not benefited from the high-dimensional data revolution that has transformed molecular and cellular phenotyping. To address this, we developed an approach that combines commercially available automated phenotyping hardware with a systems biology analysis pipeline to generate a high-dimensional readout of mouse behavior/physiology, as well as intuitive and health-relevant summary statistics (resilience and biological age). We used this platform to longitudinally evaluate aging in hundreds of outbred mice across an age range from 6 months to 3.4 years. In contrast to the assumption that aging can only be measured at the limits of animal ability via challenge-based tasks, we observed widespread physiological and behavioral aging starting in early life. Using network connectivity analysis, we found that organism-level resilience exhibited an accelerating decline with age that was distinct from the trajectory of individual phenotypes. We developed a method, Combined Aging and Survival Prediction of Aging Rate (CASPAR), for jointly predicting chronological age and survival time and showed that the resulting model is able to predict both variables simultaneously, a behavior that is not captured by separate age and mortality prediction models. This study provides a uniquely high-resolution view of physiological aging in mice and demonstrates that systems-level analysis of physiology provides insights not captured by individual phenotypes. The approach described here allows aging, and other processes that affect behavior and physiology, to be studied with sophistication and rigor.

category: systems biology


title: Asynchronous and coherent dynamics in balanced excitatory-inhibitory populations

abstract: Dynamic excitatory-inhibitory (E-I) balance is a paradigmatic mechanism invoked to explain the irregular low firing activity observed in the cortex. However, we will show that the E-I balance can be at the origin of other regimes observable in the brain. The analysis is performed by combining extensive simulations of sparse E-I networks composed of N spiking neurons with analytical investigations of low dimensional neural mass models. The bifurcation diagrams, derived for the neural mass model, allow to classify the possible asynchronous and coherent behaviours emerging in balanced E-I networks with structural heterogeneity for any finite in-degree K. In the limit N >> K >> 1 both supra and sub-threshold balanced asynchronous regimes can be observed in our system. Due to the heterogeneity the asynchronous states are characterized by the splitting of the neurons in three groups: silent, fluctuation and mean driven. These features are consistent with experimental observations reported for heterogeneous neural circuits. The coherent rhythms observed in our system can range from periodic and quasi-periodic collective oscillations (COs) to coherent chaos. These rhythms are characterized by regular or irregular temporal fluctuations joined to spatial coherence somehow similar to coherent fluctuations observed in the cortex over multiple spatial scales. The COs can emerge due to two different mechanisms. A first mechanism similar to the pyramidal-interneuron gamma (PING) one, usually invoked for the emergence of{gamma} -oscillations. The second mechanism is intimately related to the presence of current fluctuations, which sustain COs characterized by an essentially simultaneous bursting of the two populations. We observe period-doubling cascades involving the PING-like COs finally leading to the appearance of coherent chaos. Fluctuation driven COs are usually observable in our system as quasi-periodic collective motions characterized by two incommensurate frequencies. However, for sufficiently strong current fluctuations we report a novel mechanism of frequency locking among collective rhythms promoted by these intrinsic fluctuations. Our analysis suggest that despite PING-like or fluctuation driven COS are observable for any finite in-degree K, in the limit N >> K >> 1 these solutions finally result in two coexisting balanced regimes: an asynchronous and a fully synchronized one.


title: A Master Autoantigen-ome Links Alternative Splicing, Female Predilection, and COVID-19 to Autoimmune Diseases

abstract: Chronic and debilitating autoimmune sequelae pose a grave concern for the post-COVID-19 pandemic era. Based on our discovery that the glycosaminoglycan dermatan sulfate (DS) displays peculiar affinity to apoptotic cells and autoantigens (autoAgs) and that DS-autoAg complexes cooperatively stimulate autoreactive B1 cell responses, we compiled a database of 751 candidate autoAgs from six human cell types. At least 657 of these have been found to be affected by SARS-CoV-2 infection based on currently available multi-omic COVID data, and at least 400 are confirmed targets of autoantibodies in a wide array of autoimmune diseases and cancer. The autoantigen-ome is significantly associated with various processes in viral infections, such as translation, protein processing, and vesicle transport. Interestingly, the coding genes of autoAgs predominantly contain multiple exons with many possible alternative splicing variants, short transcripts, and short UTR lengths. These observations and the finding that numerous autoAgs involved in RNA-splicing showed altered expression in viral infections suggest that viruses exploit alternative splicing to reprogram host cell machinery to ensure viral replication and survival. While each cell type gives rise to a unique pool of autoAgs, 39 common autoAgs associated with cell stress and apoptosis were identified from all six cell types, with several being known markers of systemic autoimmune diseases. In particular, the common autoAg UBA1 that catalyzes the first step in ubiquitination is encoded by an X-chromosome escape gene. Given its essential function in apoptotic cell clearance and that X-inactivation escape tends to increase with aging, UBA1 dysfunction can therefore predispose aging women to autoimmune disorders. In summary, we propose a model of how viral infections lead to extensive molecular alterations and host cell death, autoimmune responses facilitated by autoAg-DS complexes, and ultimately autoimmune diseases. Overall, this master autoantigen-ome provides a molecular guide for investigating the myriad of autoimmune sequalae to COVID-19 and clues to the rare but reported adverse effects of the currently available COVID vaccines.


title: Two-dimensional multiplexed assay for rapid and deep SARS-CoV-2 serology profiling and for machine learning prediction of neutralization capacity

abstract: Antibody responses serve as the primary protection against SARS-CoV-2 infection through neutralization of viral entry into cells. We have developed a two-dimensional multiplex bead binding assay (2D-MBBA) that quantifies multiple antibody isotypes against multiple antigens from a single measurement. Here, we applied our assay to profile IgG, IgM and IgA levels against the spike antigen, its receptor-binding domain and natural and designed mutants. Machine learning algorithms trained on the 2D-MBBA data substantially improve the prediction of neutralization capacity against the authentic SARS-CoV-2 virus of serum samples of convalescent patients. The algorithms also helped identify a set of antibody isotype-antigen datasets that contributed to the prediction, which included those targeting regions outside the receptor-binding interface of the spike protein. We applied the assay to profile samples from vaccinated, immune-compromised patients, which revealed differences in the antibody profiles between convalescent and vaccinated samples. Our approach can rapidly provide deep antibody profiles and neutralization prediction from essentially a drop of blood without the need of BSL-3 access and provides insights into the nature of neutralizing antibodies. It may be further developed for evaluating neutralizing capacity for new variants and future pathogens.


title: Liposomal Dexamethasone Reduces A/H1N1 Influenza-Associated Morbidity in Mice

abstract:
Re-emerging viral threats have continued to challenge the medical and public health systems. It has become clear that a significant number of severe viral infection cases are due to an overreaction of the immune system, which leads to hyperinflammation. In this study, we aimed to demonstrate the therapeutic efficacy of the dexamethasone nanomedicine in controlling the symptoms of influenza infection. We found that the A/Wisconsin/WSLH34939/2009 (H1N1) infection induced severe pneumonia in mice with a death rate of 80%, accompanied by significant epithelial cell damage, infiltration of immune cells, and accumulation of pro-inflammatory cytokines in the airway space. Moreover, the intranasal delivery of liposomal dexamethasone during disease progression reduced the death rate by 20%. It also significantly reduced the protein level of tumor necrosis factor-alpha (TNF), interleukin1{beta} (IL-1{beta}), IL-6, and the C-X-C motif chemokine ligand 2 (CXCL2) as well as the number of infiltrated immune cells in the bronchoalveolar lavage fluids as compared to the free drug. It was found that the liposomal dexamethasone was mainly distributed into the monocyte/macrophages in the lungs, suggesting its mode of action via the specific delivery of the drug into myeloid cells as a major cell population for inducing the cytokine storm. Taken together, the intranasal delivery of liposomal dexamethasone may serve as a novel promising therapeutic strategy for the treatment of influenza A-induced pneumonia.

ImportanceInfluenza A virus causes annual epidemics and sporadic pandemics of respiratory diseases. An excessive immune response, called a cytokine storm, caused by the influenza virus is the most common complication of influenza infection that is associated with high levels of morbidity and mortality. Here, we report that the liposomal dexamethasone targeting the monocytes/macrophages during disease progression reduced the death rate as well as the protein level of inflammatory cytokines and the number of immune cells. Therefore, the findings of this study may be applied to address the limitations of the current strategies for the treatment of influenza.


title: Single-sequence protein structure prediction using language models from deep learning

abstract: AlphaFold2 and related systems use deep learning to predict protein structure from co-evolutionary relationships encoded in multiple sequence alignments (MSAs). Despite dramatic, recent increases in accuracy, three challenges remain: (i) prediction of orphan and rapidly evolving proteins for which an MSA cannot be generated, (ii) rapid exploration of designed structures, and (iii) understanding the rules governing spontaneous polypeptide folding in solution. Here we report development of an end-to-end differentiable recurrent geometric network (RGN) able to predict protein structure from single protein sequences without use of MSAs. This deep learning system has two novel elements: a protein language model (AminoBERT) that uses a Transformer to learn latent structural information from millions of unaligned proteins and a geometric module that compactly represents C backbone geometry. RGN2 outperforms AlphaFold2 and RoseTTAFold (as well as trRosetta) on orphan proteins and is competitive with designed sequences, while achieving up to a 106-fold reduction in compute time. These findings demonstrate the practical and theoretical strengths of protein language models relative to MSAs in structure prediction.


title: Cooled SPAD array detector for low light-dose fluorescence laser scanning microscopy

abstract:
The single-photon timing and sensitivity performance and the imaging ability of asynchronous-readout single-photon avalanche diode (SPAD) array detectors have opened up enormous perspectives in fluorescence (lifetime) laser scanning microscopy (FLSM), such as super-resolution image scanning microscopy and high-information content fluorescence fluctuation spectroscopy (FFS). However, the strengths of these FLSM techniques depend on the many different characteristics of the detector, such as dark-noise, photon-detection efficiency, after-pulsing probability, and optical-cross talk, whose overall optimization is typically a trade-off between these characteristics. To mitigate this trade-off, we present a novel SPAD array detector with an active cooling system, which substantially reduces the dark-noise without significantly deteriorating any other detector characteristics. In particular, we show that lowering the temperature of the sensor to -15{degrees}C significantly improves the signal-to-noise ratio due to a 10-fold decrease in the dark-count rate compared to room temperature. As a result, for imaging, the laser power can be decreased by more than a factor of three, which is particularly beneficial for live-cell super-resolution imaging, as demonstrated in fixed and living cells expressing GFP-tagged proteins. For FFS, together with the benefit of the reduced laser power, we show that cooling the detector is necessary to remove artifacts in the correlation function, such as spurious negative correlations observed in the hot elements of the detector, i.e., elements whose dark-noise is substantially higher than the median value. Overall, this detector represents a further step towards the integration of SPAD array detectors in any FLSM system.

SIGNIFICANCESingle-photon avalanche diode (SPAD) array detectors are revolutionizing fluorescence laser-scanning microscopy (FLSM). Thanks to their single-photon timing and sensitivity ability and their imaging faculty, a SPAD array detector transforms any FLSM into a super-resolution microscope, and opens a whole range of possibilities for the study of sample dynamics by means of fluorescence fluctuation spectroscopy (FFS). However, dark-noise can be a severe problem for both imaging and FFS. For imaging, the signal overcomes noise only for a relatively high illumination intensity, which can be detrimental for live-cell experiments. For FFS, the noise leads to artifacts in the correlation curves, potentially leading to wrong conclusions about the sample. We show that lowering the temperature of the detector to -15{degrees}C solves both problems

category: biophysics


title: Redox state and cellular uptake of copper is regulated by the N-terminus of human Copper Transporter-1

abstract: Copper(I) is essential for all life forms. Though Cu(II) is most abundant state in environment, its reduction to Cu(I) is prerequisite for bio-utilization, by a mechanism that is uncharacterized. We show that in human Copper Transporter-1, two amino-terminal methionine-histidine clusters and neighbouring aspartates distinctly binds Cu(II) and Cu(I) preceding its import. The endocytosis of hCTR1 from basolateral membrane of polarized epithelia to Common-Recycling-Endosomes is dependent on copper reduction and Cu(I) coordination by methionines. The transient binding of both Cu(II) and Cu(I) during the reduction process facilitated by aspartates acts as another crucial determinant of hCTR1 endocytosis. Mutating 7Met-Gly-Met9 and Asp13 abrogates copper uptake and endocytosis that is correctable by reduced and non-reoxidizable Cu(I). Histidines clusters are crucial for hCTR1 functioning at limiting copper. Finally, we show that two N-terminal His-Met-Asp clusters exhibit functional complementarity in regulating Cu(I)-induced hCTR1 endocytosis. We propose a mechanistic model where His-Met-Asp residues of amino-terminal hCTR1 coordinates copper and maintains its reduced state crucial for uptake.

category: cell biology


title: Three-axis classification of mouse lung mesenchymal cells reveals two populations of myofibroblasts

abstract: The mesenchyme consists of heterogeneous cell populations that support neighboring structures and are integral to intercellular signaling. Despite such importance, mesenchymal cell types are poorly defined morphologically and molecularly, lagging behind their counterparts in the epithelial, endothelial, and immune lineages. Leveraging single-cell RNA-seq, three-dimensional imaging, and lineage tracing, we classify the mouse lung mesenchyme into three proximal-distal axes that are associated with the endothelium, epithelium, and interstitium, respectively. From proximal to distal, (1) the vascular axis includes vascular smooth muscle cells and pericytes that transition as arterioles and venules ramify into capillaries; (2) the epithelial axis includes airway smooth muscle cells and two populations of myofibroblasts: ductal myofibroblasts, surrounding alveolar ducts and marked by CDH4, HHIP, and Lgr6, which persist post-alveologenesis, and alveolar myofibroblasts, surrounding alveoli and marked by high expression of PDGFRA, which undergo developmental apoptosis; (3) the interstitial axis, residing between the epithelial and vascular trees and sharing a newly-identified marker MEOX2, includes fibroblasts in the bronchovascular bundle and the alveolar interstitium that are marked by IL33/DNER/PI16 and Wnt2, respectively. Single-cell imaging reveals distinct morphology of each mesenchymal cell population. This classification provides a conceptual and experimental framework applicable to other organs.

category: developmental biology


title: Evolutionary flexibility in routes to mat formation by Pseudomonas

abstract: Many bacteria form mats at the air-liquid interface of static microcosms. These structures typically involve the secretion of exopolysaccharide(s), the production of which is often controlled by the secondary messenger c-di-GMP. Mechanisms of mat formation have been particularly well characterized in Pseudomonas fluorescens SBW25; mutations that lead to an increase in c-di-GMP production by diguanylate cyclases (WspR, AwsR, or MwsR) result in the secretion of cellulose, and mat formation. Here, we characterize and compare mat formation in two close relatives of SBW25: Pseudomonas simiae PICF7 and Pseudomonas fluorescens A506. We find that PICF7 - the strain more closely related to SBW25 - can form mats through mutations affecting the activity of the same three diguanylate cyclases as SBW25. However, instead of cellulose, these mutations activate the production of the Pel exopolysaccharide. We also provide evidence for at least two further - as yet uncharacterized - routes to PICF7 mat formation. P. fluorescens A506, while retaining the same mutational routes to mat formation as SBW25 and PICF7, forms mats by a semi-heritable mechanism that likely culminates in Pga and/or Psl production. Overall, our results demonstrate a high level of evolutionary flexibility in the molecular and structural routes to mat formation, even among close relatives.

category: evolutionary biology


title: Evolution of functional diversity in the holozoan tyrosine kinome

abstract: The emergence of multicellularity is strongly correlated with the expansion of tyrosine kinases, a conserved family of signaling enzymes that regulates pathways essential for cell-to-cell communication. Although tyrosine kinases have been classified from several model organisms, a molecular-level understanding of tyrosine kinase evolution across all holozoans is currently lacking. Using a hierarchical sequence constraint-based classification of diverse holozoan tyrosine kinases, we construct a new phylogenetic tree that identifies two ancient clades of cytoplasmic and receptor tyrosine kinases separated by the presence of an extended insert segment in the kinase domain connecting the D and E helices. Present in nearly all receptor tyrosine kinases, this fast-evolving insertion imparts diverse functionalities such as post-translational modification sites and regulatory interactions. Eph and EGFR receptor tyrosine kinases are two exceptions which lack this insert, each forming an independent lineage characterized by unique functional features. We also identify common constraints shared across multiple tyrosine kinase families which warrant the designation of three new subgroups: Src Module (SrcM), Insulin Receptor Kinase-Like (IRKL), and Fibroblast, Platelet-derived, Vascular, and growth factor Receptors (FPVR). Subgroup-specific constraints reflect shared autoinhibitory interactions involved in kinase conformational regulation. Conservation analyses describe how diverse tyrosine kinase signaling functions arose through the addition of family-specific motifs upon subgroup-specific features and co-evolving protein domains. We propose the oldest tyrosine kinases, IRKL, SrcM, and Csk, originated from unicellular pre-metazoans and were co-opted for complex multicellular functions. The increased frequency of oncogenic variants in more recent tyrosine kinases suggests that lineage-specific functionalities are selectively altered in human cancers.


title: A general hidden Markov model for multi-parent populations

abstract: A common step in the analysis of multi-parent populations is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of multi-parent populations. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly-available data on Diversity Outbred and Collaborative Cross mice

category: genetics


title: A generic hidden Markov model for multi-parent populations


title: A distinct core regulatory module enforces oncogene expression in KMT2A-rearranged leukemia

abstract: A small set of lineage-restricted transcription factors (TFs), termed core regulatory circuitry (CRC), control cell identity and malignant transformation. Here, we integrated gene dependency, chromatin architecture and TF perturbation datasets to characterize 31 core TFs in acute myeloid leukemia (AML). Contrary to a widely accepted model, we detected a modular CRC structure with hierarchically organized, partially redundant and only sparsely interconnected modules of core TFs controlling distinct genetic programs. Rapid TF degradation followed by measurement of genome-wide transcription rates revealed that core TFs directly regulate dramatically fewer genes than previously assumed. Leukemias carrying KMT2A (MLL) rearrangements depend on the IRF8/MEF2 axis to directly enforce expression of the key oncogenes MYC, HOXA9 and BCL2. Our datasets provide an evolving model of CRC organization in human cells, and a resource for further inquiries into and therapeutic targeting of aberrant transcriptional circuits in cancer.

category: genomics


title: IMPDH1 retinal variants control filament architecture to tune allosteric regulation

abstract: IMP dehydrogenase (IMPDH), a key regulatory enzyme in purine nucleotide biosynthesis, dynamically assembles filaments in response to changes in metabolic demand. Humans have two isoforms: IMPDH2 filaments reduce sensitivity to feedback inhibition by the downstream product GTP, while IMPDH1 assembly remains uncharacterized. IMPDH1 plays a unique role in retinal metabolism, and point mutants cause blindness and disrupt GTP regulation. Here, in a series of cryo-EM structures we show that IMPDH1 assembles polymorphic filaments with different assembly interfaces in active and inhibited states. Retina-specific splice variants introduce structural elements that reduce sensitivity to GTP inhibition, including stabilization of the active filament form. Finally, we show that IMPDH1 disease mutations fall into two classes: one disrupts GTP regulation and the other has no in vitro phenotype. These findings provide a foundation for understanding the role of IMPDH1 in retinal function and disease and demonstrate the diverse mechanisms by which metabolic enzyme filaments are allosterically regulated.

category: biochemistry


title: Increasing seed thiamin content impacts stored carbon partitioning and subsequent seedling stress tolerance

abstract: Thiamin and thiamin pyrophosphate (TPP) are essential components for the function of enzymes involved in the metabolism of carbohydrates and amino acids in living organisms. In addition to its role as a cofactor, thiamin plays a key role in resistance against biotic and abiotic stresses in plants. Most of the studies used exogenous thiamin to enhance stress tolerance in plants. In this study, we achieved this objective by genetically engineering Arabidopsis thaliana and Camelina sativa for the seed-specific co-overexpression of the Arabidopsis thiamin biosynthetic genes Thi4, ThiC, and ThiE. Elevated thiamin content in the seeds of transgenic plants was accompanied by the enhanced expression levels of transcripts encoding thiamin cofactor-dependent enzymes. Furthermore, seed germination and root growth in thiamin over-producing lines were more tolerant to oxidative stress caused by salt and paraquat treatments. The transgenic seeds also accumulated more oil (up to16.4% in Arabidopsis and17.9% in C. sativa) and carbohydrate but less protein than the control seeds. The same results were also observed in TPP over-producing Arabidopsis plants generated by the seed-specific overexpression of TPK1. Together, our findings suggest that thiamin and TPP over-production in transgenic lines confer a boosted abiotic stress tolerance and alter the seed carbon partitioning as well.

category: bioengineering


title: Nanofibrous PEDOT-Carbon Composite on Flexible Probes for Soft Neural Interfacing

abstract: In this study, we report a flexible implantable 4-channel microelectrode probe coated with highly porous and robust nanocomposite of poly(3,4-ethylenedioxythiophene) (PEDOT) and carbon nanofiber (CNF) as a solid doping template for high-performance in vivo neuronal recording and stimulation. A simple yet well-controlled deposition strategy was developed via in situ electrochemical polymerization technique to create a porous network of PEDOT and CNFs on a flexible 4-channel gold microelectrode probe. Different morphological and electrochemical characterizations showed that they exhibit remarkable and superior electrochemical properties, yielding microelectrodes combining high surface area, low impedance (16.8 {+/-} 2 M{Omega}.m2 at 1 kHz) and elevated charge injection capabilities (7.6 {+/-} 1.3 mC/cm2) that exceed those of pure and composite PEDOT layers. In addition, the PEDOT-CNF composite electrode exhibited extended biphasic charge cycle endurance, resulting in a negligible physical delamination or degradation for long periods of electrical stimulation. In vitro testing on mouse brain slices showed that they can record spontaneous oscillatory field potentials as well as single-unit action potentials and allow to safely deliver electrical stimulation for evoking field potentials. The combined superior electrical properties, durability and 3D microstructure topology of the PEDOT-CNF composite electrodes demonstrate outstanding potential for developing future neural surface interfacing applications.


title: Identifying robust functional modules using three-body correlations in Escherichia coli

abstract: Understanding the underlying structure of a gene regulatory network is crucial to understand the biological functions of genes or groups of genes. A common strategy to investigate it is to find community structure of these networks. However, methods of finding these communities are often sensitive to noise in the gene expression data and the inherent stochasticity of the community detection algorithms. Here we introduce an approach for identifying functional groups and their hierarchical organization in gene co-expression networks from expression data. A network describing the relatedness in the expression profiles of genes is first inferred using an information theoretic approach. Community structure within the inferred network is found by using modularity maximization. This community structure is further refined using three-body structural correlations to robustly identify important functional gene communities. We apply this approach to the expression data of Escherichia coli genes and identify 25 robust groups, many of which show key associations with important biological functions as demonstrated by Gene Ontology (GO) term enrichment analysis. Thus, our approach makes specific and novel predictions about the function of these genes.


title: The amyloid concentric β-barrel hypothesis: Models of amyloid beta 42 oligomers and annular protofibrils

abstract: Amyloid beta (A{beta}) peptides, a major contributor to Alzheimers disease, occur in differing lengths, each of which forms a multitude of assembly types. The most toxic soluble oligomers are formed by A{beta}42; some of which have antiparallel {beta}-sheets. Previously, our group proposed molecular models of A{beta}42 hexamers in which the C-terminus third of the peptide (S3) forms an antiparallel 6-stranded {beta}-barrel that is surrounded by an antiparallel barrel formed by the more polar N-terminus (S1) and middle (S2) portions. These hexamers were proposed to act as seeds from which dodecamers, octadecamers, both smooth and beaded annular protofibrils, and transmembrane channels form. Since then, numerous aspects of our models have been supported by experimental findings. Recently, NMR-based structures have been proposed for A{beta}42 tetramers and octamers, and NMR studies have been reported for oligomers composed of ~ 32 monomers. Here we propose a range of concentric {beta}-barrel models and compare their dimensions to image-averaged electron micrographs of both beaded annular protofibrils (bAPFs) and smooth annular protofibrils (sAPFs) of A{beta}42. The smaller oligomers have 6, 8, 12, 16, and 18 monomers. These beads string together to form necklace-like bAPFs. These gradually morph into sAPFs in which a S3 {beta}-barrel is shielded on one or both sides by {beta}-barrels formed from S1 and S2 segments.


title: The role of actin and myosin II in the cell cortex of adhered and suspended cells

abstract: Adhesion induces dramatic morphological and mechanical changes to cells, which are reflected by changes to the actin cortex. Among the many different proteins involved in this sub-membranous layer, motor proteins (e.g., nonmuscle myosin II [NMII]) and actin nucleators (e.g., Arp2/3, formins) are known to have significant influences on its dynamics and structure. The different roles of NMII, Arp2/3, and formins in the dynamics, structure, and mechanics of the actin cortex depend on the adhesion state of the cell. In this study, we unravel the interplay between the dynamics, structure, and mechanics of the actin cortex in adhered cells and in cells in suspension. We show that treatments with extrinsic cellular perturbants lead to alterations of all three properties that are correlated. However, intrinsic actin cortex variations between different cell adhesion states lead to unexpected correlations. Surprisingly, we find that NMII minifilaments have a minor influence on the actin cortex. Using new microscopy techniques, we show that NMII minifilaments are not localized within the actin cortex, as previously thought, but concentrated in a layer beneath it. Our treatments affecting Arp2/3 and formin reveal correlations between the actin cortex characteristics. Our data build towards a comprehensive understanding of the actin cortex. This understanding allows the prediction and control of cortical changes, which is essential for the study of general cellular processes, such as cell migration, metastasis, and differentiation.


title: Stop Bickering! Reconciling Signaling Pathway Databases with Network Topologies

abstract: A major goal of molecular systems biology is to understand the coordinated function of genes or proteins in response to cellular signals and to understand these dynamics in the context of disease. Signaling pathway databases such as KEGG, NetPath, NCI-PID, and Panther describe the molecular interactions involved in different cellular responses. While the same pathway may be present in different databases, prior work has shown that the particular proteins and interactions differ across database annotations. However, to our knowledge no one has attempted to quantify their structural differences. It is important to characterize artifacts or other biases within pathway databases, which can provide a more informed interpretation for downstream analyses. In this work, we consider signaling pathways as graphs and we use topological measures to study their structure. We find that topological characterization using graphlets (small, connected subgraphs) distinguishes signaling pathways from appropriate null models of interaction networks. Next, we quantify topological similarity across pathway databases. Our analysis reveals that the pathways harbor database-specific characteristics implying that even though these databases describe the same pathways, they tend to be systematically different from one another. We show that pathway-specific topology can be uncovered after accounting for database-specific structure. This work present the first step towards elucidating common pathway structure beyond their specific database annotations.


title: Using single-cell models to predict the functionality of synthetic circuits at the population scale

abstract: Mathematical modeling has become a major tool to guide the characterization and synthetic construction of cellular processes. However, models typically lose their capacity to explain or predict experimental outcomes as soon as any, even minor, modification of the studied system or its operating conditions is implemented. This limits our capacity to fully comprehend the functioning of natural biological processes and is a major roadblock for the de novo design of complex synthetic circuits. Here, using a specifically constructed yeast optogenetic differentiation system as an example, we show that a simple deterministic model can explain system dynamics in given conditions but loses validity when modifications to the system are made. On the other hand, deploying theory from stochastic chemical kinetics and developing models of the systems components that simultaneously track single-cell and population processes allows us to quantitatively predict emerging dynamics of the system without any adjustment of model parameters. We conclude that carefully characterizing the dynamics of cell-to-cell variability using appropriate modeling theory may allow one to unravel the complex interplay of stochastic single-cell and population processes and to predict the functionality of composed synthetic circuits in growing populations before the circuit is constructed.


title: Learning directed acyclic graphs for ligands and receptors based on spatially resolved transcriptomic analysis of ovarian cancer

abstract: A critical bottleneck step towards understanding the immune activation and suppression mechanisms in tumor samples is to identify transcriptional signals underlying the cell-cell communication between tumor cells and immune/stromal cells in tumor microenvironments (TME). Cell-to-cell communication extensively relies on interactions between secreted ligands and cell-surface receptors, which create a highly connected signaling network through many ligand-receptor paths. The latest advance in in situ omics analyses, such as spatial transcriptomic (ST) analysis, provide unique opportunities to directly characterize ligand-receptor signaling networks that powers cell-cell communication, which has not been feasible based on either bulk or single-cell omics data. In this paper, we focus on high grade serous ovarian cancer (HGSC), and propose a novel statistical method, DAGBagST, to characterize the ligand-receptor interaction networks between adjacent tumor and stroma cells in ovarian tumors based on spatial transcriptomic data. DAGBagST utilizes a directed acyclic graph (DAG) model with a novel approach to handle the zero-inflated distribution observed in the ST data. It also leverages existing ligand-receptor regulation databases as prior information, and employs a bootstrap aggregation strategy to achieve robust network estimation. We applied DAGBagST to ST datasets of tumor samples from four HGSC patients, and identified common and distinct ligand-receptor regulations between adjacent tumor and stromal grids across multiple tumors. These results cast light on biological processes relating to the communication between tumor and TME cells in these ovarian tumors. An open-source R package for performing DAGBagST will be shared upon publication.


title: The control and training of single motor units in isometric tasks are constrained by a common synaptic input signal

abstract: Recent developments in neural interfaces enable the real-time and non-invasive tracking of motor neuron spiking activity. Such novel interfaces provide a promising basis for human motor augmentation by extracting potential high-dimensional control signals directly from the human nervous system. However, it is unclear how flexibly humans can control the activity of individual motor neurones to effectively increase the number of degrees-of-freedom available to coordinate multiple effectors simultaneously. Here, we provided human subjects (N=7) with real-time feedback on the discharge patterns of pairs of motor units (MUs) innervating a single muscle (tibialis anterior) and encouraged them to independently control the MUs by tracking targets in a 2D space. Subjects learned control strategies to achieve the target-tracking task for various combinations of MUs. These strategies rarely corresponded to a volitional control of independent input signals to individual MUs. Conversely, MU activation was consistent with a common input to the MU pair, while individual activation of the MUs in the pair was predominantly achieved by alterations in de-recruitment order that could be explained with history-dependent changes in motor neuron excitability. These results suggest that flexible MU control based on independent synaptic inputs to single MUs is not a simple to learn control strategy.


title: A Data-Dependent Acquisition Ladder for Ultrasensitive (Neuro)Proteomics

abstract: Measurement of broad types of proteins from a small number of cells to single cells would help to better understand the nervous system but requires significant leaps in high-resolution mass spectrometry (HRMS) sensitivity. Microanalytical capillary electrophoresis electrospray ionization (CE-ESI) offers a path to ultrasensitive proteomics by integrating scalability with sensitivity. We report here a data acquisition strategy that expands the detectable and quantifiable proteome in trace amounts of digests using CE-ESI-HRMS. Data-dependent acquisition (DDA) was programmed to progressively exclude high-intensity peptide signals during repeated measurements. These nested experiments formed rungs of our "DDA ladder." The method was tested for replicates analyzing ~500 pg of protein digest from cultured hippocampal (primary) neurons (mouse), which estimates to the total amount of protein from a single neuron. Analysis of net amounts approximating to ~10 neurons identified 428 nonredundant proteins (415 quantified), an ~35% increase over traditional DDA. The identified proteins were enriched in neuronal marker genes and molecular pathways of neurobiological importance. The DDA ladder deepens the detectable proteome from trace amounts of proteins, expanding the analytical toolbox of neuroscience.


title: Concentration-dependent change in hypothalamic neuronal transcriptome by the dietary fatty acids: oleic and palmitic acids

abstract: Prenatal high-fat diet exposure increases hypothalamic neurogenesis events in embryos and programs offspring to be obesity-prone. The molecular mechanism involved in these dietary effects of neurogenesis are unknown. This study investigated the effects of oleic and palmitic acids, which are abundant in a high-fat diet, on the hypothalamic neuronal transcriptome and how these changes impact neurogenesis events. The results show differential effects of low and high concentrations of oleic or palmitic acid treatment on differential gene transcription. Gene ontology analysis uncovered significant gene enrichment in several cellular pathways involved in gene regulation and protein production, particularly with proliferation, migration, and cell survival. The enriched signaling pathways include Wnt, integrin, PDGF, and apoptosis, in addition endocrine function signaling pathways CCKR and GnRH. Further examination of proliferation and migration show low concentrations of oleic acid to stimulate proliferation and high concentrations of both oleic and palmitic acid to stimulate apoptosis. Oleic acid also reduced hypothalamic neuronal migration, with little effects by palmitic acid. The results show direct impact of the two most abundant fatty acids in a high fat diet to directly impact hypothalamic neuronal proliferation and migration. The results also uncovered signaling pathways affected by oleic and palmitic acid and suggest a mechanism of prenatal high-fat diet induced neurogenesis events is through these two abundant fatty acids.


title: Altered low frequency brain rhythms precede changes in gamma power during tauopathy

abstract: Alzheimers disease and other dementias are associated with disruptions of electrophysiological brain activity, including low frequency and gamma rhythms. Many of these dementias are also associated with the malfunction of the membrane associated protein tau. Tauopathy disrupts neuronal function and the stability of synapses and is a key driver of neurodegeneration. Here we ask how brain rhythms are affected by tauopathy, at different stages of its progression. We performed local field potential recordings from visual cortex of rTg4510 and control animals at early stages of neurodegeneration (5 months) and at a more advanced stage where pathology is evident (8 months). We measured brain activity in the presence or absence of external visual stimulation, and while monitoring pupil diameter and locomotion to establish animal behavioural states. At 5 months, before substantial pathology, we found an increase in low frequency rhythms during resting state in tauopathic animals. This was because tauopathic animals entered intermittent periods of increased neural synchronisation, where activity across a wide band of low frequencies was strongly correlated. At 8 months, when the degeneration was more advanced, the increased synchronisation and low frequency power was accompanied by a reduction in power in the gamma range, with diverse effects across different components of the gamma rhythm. Our results indicate that slower rhythms are impaired earlier than gamma rhythms in tauopathy, suggesting that electrophysiological measurements can indicate both the presence and progression of tauopathic degeneration.


title: MapToCleave: high-throughput profiling of microRNA biogenesis in living cells

abstract: Previous large-scale studies have uncovered many features that determine the processing of microRNA (miRNA) precursors, however, they have been conducted in vitro. Here we introduce MapToCleave, a new method to simultaneously profile processing of thousands of distinct RNA structures in living cells. Our new in cell method captures essentially all the biogenesis features that have been discovered through near two decades of in vitro studies - providing support for both approaches. We find that miRNA precursors with a stable lower basal stem are more efficiently processed and also have higher expression in vivo in tissues from twenty animal species. We systematically compare the importance of known and novel sequence and structural features and test biogenesis of miRNA precursors from ten animal and plant species in human cells. Lastly, we provide evidence that the GHG motif better predicts processing when defined as a structure rather than sequence motif, consistent with recent cryo-EM studies. In summary, we apply a new screening assay in living cells to reveal the importance of lower basal stem stability for miRNA processing and in vivo expression.

category: molecular biology


title: Structures of full-length glycoprotein hormone receptor signaling complexes

abstract: Luteinizing hormone (LH) and chorionic gonadotropin (CG) are members of the glycoprotein hormone family essential to human reproduction and are important therapeutic drugs. They activate the same G protein-coupled receptor, LHCGR, by binding to the large extracellular domain (ECD). Here we report four cryo-EM structures of LHCGR, two wildtype receptor structures in the inactive and active states, and two constitutively active mutated receptor structures. The active structures are bound to CG and Gs heterotrimer, with one of the structure also containing the allosteric agonist, Org43553. The structures reveal a distinct  push and pull" mechanism of receptor activation, in which the ECD is pushed by the bound hormone and pulled by the extended hinge loop next to the transmembrane domain (TMD). A highly conserved 10-residue fragment (P10) from the hinge C-terminal loop at the ECD-TMD interface functions as a tethered agonist to induce conformational changes in TMD and G-protein coupling. Org43553 binds to a TMD pocket and interacts directly with P10 that further stabilizes the receptor in the active conformation. Together, these structures provide a common model for understanding glycoprotein hormone signal transduction and dysfunction, and inspire the search for clinically suitable small molecular compounds to treat endocrine diseases.


title: St18 specifies globus pallidus projection neuron identity in MGE lineage

abstract: The medial ganglionic eminence (MGE) is a progenitor domain in the subpallium that produces both locally-projecting interneurons which undergo tangential migration in structures such as the cortex as well as long-range projection neurons that occupy subcortical nuclei. Very little is known about the transcriptional mechanisms specifying the migratory behavior and axonal projection patterns of these two broad classes of MGE-derived neurons. In this study, we identify St18 as a novel transcriptional determinant specifying projection neuron fate in the MGE lineage. St18 is transiently expressed in the MGE subventricular zone (SVZ) and mantle, and we assessed its function using an ES cell-based model of MGE development. Induction of St18 is sufficient to direct ES-derived MGE neurons to adopt a projection neuron-like identity as defined by migration and morphology. Using genetic loss-of-function in mice, we find that St18 is required for the production of globus pallidus pars externa (GPe) prototypic projection neurons. Single cell RNA sequencing revealed that St18 regulates MGE output of specific neuronal populations: in the absence of St18, we observe a large expansion of cortical interneurons at the expense of putative GPe neurons. Through gene expression analysis we identified a downstream effector of St18, Cbx7, which is a component of Polycomb repressor complex 1. We find that Cbx7 is essential for projection neuron-like migration and is not involved in St18-mediated projection neuron-like morphology. Our results characterize a novel transcriptional determinant that directs GPe prototypic projection neuron identity. Further, we identified a downstream target of St18, Cbx7, which regulates only the migratory behavior of long-range projection neurons, suggesting that specific features of MGE projection neuron identity may be governed in a compartmentalized fashion by distinct transcriptional modules downstream of St18.


title: A rear-engine drives adherent tissue migration in vivo

abstract: During animal embryogenesis, homeostasis and disease, tissues push and pull on their surroundings to move forward. Although the force-generating machinery is known, it is unknown how tissues exert physical stresses on their substrate to generate motion in vivo. Here, we identify the force transmission machinery, the substrate, and the stresses that a tissue, the zebrafish posterior lateral line primordium, generates during its migration. We find that the primordium couples actin flow through integrins to the basement membrane for forward movement. Talin/integrin-mediated coupling is required for efficient migration and its loss is partly compensated for by increased actin flow. Using Embryogram, an approach to measure stresses in vivo, we show that the primordiums rear exerts high stresses, indicating that this tissue pushes itself forward with its back. This unexpected strategy likely also underlies the motion of other tissues in animals


title: Selection bias in mutation accumulation

abstract: Mutation accumulation (MA) experiments, in which de novo mutations are sampled and subsequently characterized, are an essential tool in understanding the processes underlying evolution. In microbial populations, MA protocols typically involve a period of population growth between severe bottlenecks, such that a single individual can form a visible colony. While it has long been appreciated that the action of positive selection during this growth phase cannot be eliminated, it is typically assumed to be negligible. Here, we quantify the effect of both positive and negative selection in MA studies, demonstrating that selective effects can substantially bias the distribution of fitness effects (DFE) and mutation rates estimated from typical MA protocols in microbes. We then present a simple correction for this bias which applies to both beneficial and deleterious mutations, and can be used to correct the observed DFE in multiple environments. Finally, we use simulated MA experiments to illustrate the extent to which the MA-inferred DFE differs from the underlying true DFE, and demonstrate that the proposed correction accurately reconstructs the true DFE over a wide range of scenarios. These results highlight that positive selection during microbial MA experiments is in fact not negligible, but can be corrected to gain a more accurate understanding of fundamental evolutionary parameters.


title: TAS-Seq: a robust and sensitive amplification method for beads-based scRNA-seq

abstract: Single-cell RNA-sequencing (scRNA-seq) is valuable for analyzing cellular heterogeneity. Cell composition accuracy is critical for analyzing cell-cell interaction networks from scRNA-seq data. We developed terminator-assisted solid-phase cDNA amplification and sequencing (TAS-Seq) for scRNA-seq based on a terminator, terminal transferase, and nanowell/bead-based scRNA-seq platform; TAS-Seq showed high tolerance to variations in the terminal transferase reaction, which complicate the handling of existing terminal transferase-based scRNA-seq methods. In murine and human lung samples, TAS-Seq yielded scRNA-seq data that were highly correlated with flow-cytometric data, showing higher gene-detection sensitivity and more robust detection of important cell-cell interactions and expression of growth factors/interleukins in cell subsets than 10X Chromium v2 and Smart-seq2. Expanding TAS-Seq application will improve understanding and atlas construction of lung biology at the single-cell level.


title: TAS-Seq: a robust and sensitive amplification method for bead-based scRNA-seq


title: Ancestral contributions to contemporary European complex traits

abstract:
The contemporary European genetic makeup formed in the last 8000 years as the combination of three main genetic components: the local Western Hunter-Gatherers, the incoming Neolithic Farmers from Anatolia and the Bronze Age component from the Pontic Steppes. When meeting into the post-Neolithic European environment, the genetic variants accumulated during their three distinct evolutionary histories mixed and came into contact with new environmental challenges.

Here we investigate how this genetic legacy reflects on the complex trait landscape of contemporary European populations, using the Estonian Biobank as a case study.

For the first time we directly connect the phenotypic information available from biobank samples with the genetic similarity to these ancestral groups, both at a genome-wide level and focusing on genomic regions associated with each of the 27 complex traits we investigated. We also found SNPs connected to pigmentation, cholesterol, sleep, diastolic blood pressure, and body mass index (BMI) to show signals of selection following the post Neolithic admixture events. We recapitulate existing knowledge about pigmentation traits, corroborate the connection between Steppe ancestry and height and highlight novel associations. Among others, we report the contribution of Hunter Gatherer ancestry towards high BMI and low blood cholesterol levels.

Our results show that the ancient components that form the contemporary European genome were differentiated enough to contribute ancestry-specific signatures to the phenotypic variability displayed by contemporary individuals in at least 11 out of 27 of the complex traits investigated here.


title: The impact of eye closure on anticipatory alpha activity in a tactile discrimination task

abstract:
One of the very first observations made regarding alpha oscillations (8-14 Hz), is that they increase in power over posterior areas when awake participants close their eyes. Recent work, especially in the context of (spatial) attention, suggests that alpha activity reflects a mechanism of functional inhibition. However, it remains unclear how eye closure impacts anticipatory alpha modulation observed in attention paradigms, and how this affects subsequent behavioral performance. Here, we recorded magnetoencephalography (MEG) in 33 human participants performing a tactile discrimination task with their eyes open vs. closed. We replicated the hallmarks of previous somatosensory spatial attention studies: alpha lateralization across the somatosensory cortices as well as alpha increase over posterior regions. Furthermore, we found that eye closure leads to (i) reduced task performance, (ii) widespread increase in alpha power, and (iii) reduced anticipatory visual alpha modulation (iv) with no effect on somatosensory alpha lateralization. Regardless of whether participants had their eyes open or closed, increased posterior alpha power and somatosensory alpha lateralization improved their performance. Thus, we provide evidence that eye closure does not alter the impact of anticipatory alpha modulations on behavioral performance. We propose there is an optimal posterior alpha level for somatosensory task performance, which can be achieved through a combination of eye closure and top-down anticipatory attention.

Significance StatementAlpha oscillations are dominant when awake participants have their eyes closed. Furthermore, alpha is known to modulate with anticipatory attention, and has been ascribed a role of active functional inhibition. Surprisingly, the link between anticipatory alpha and eye closure remains unclear. Here we collected MEG data while human participants performed a tactile discrimination task either with their eyes open or closed. Eye closure led to a widespread increase in alpha power, and affected anticipatory visual alpha modulation but not somatosensory alpha lateralization. Importantly, eye closure did not affect the correlation between alpha and task performance. Our findings provide novel insights into how eye closure impacts anticipatory alpha modulation, and how optimal alpha levels for task performance can be achieved differently.


title: Contribution of arm movements to recovery after a trip in older adults.

abstract:
Falls are common in daily life, often caused by trips and slips and, particularly in older adults, with serious consequences. Although arm movements play an important role in balance control, there is limited research into the role of arm movements during balance recovery after tripping in older adults. We investigated how older adults use their arms to recover from a trip and the difference in the effects of arm movements between fallers (n=5) and non-fallers (n=11).

Sixteen older males and females (69.7{+/-}2.3 years) walked along a walkway and were occasionally tripped over suddenly appearing obstacles. We analysed the first trip using a biomechanical model based on full-body kinematics and force-plate data to calculate whole body orientation during the trip and recovery phase. With this model, we simulated the effects of arm movements at foot-obstacle impact and during trip recovery on body orientation.

Apart from an increase in sagittal plane forward body rotation at touchdown in fallers, we found no significant differences between fallers and non-fallers in the effects of arm movements on trip recovery. Like earlier studies in young adults, we found that arm movements during the recovery phase had most favourable effects in the transverse plane: by delaying the transfer of angular momentum of the arms to the body, older adults rotated the tripped side more forward thereby allowing for a larger recovery step. Older adults that are prone to falling might improve their balance recovery after tripping by learning to prolong ongoing arm movements.


title: Contribution of arm movements to balance recovery after tripping in older adults.


title: Bi-modal Variational Autoencoders for Metabolite Identification Using Tandem Mass Spectrometry

abstract: A grand challenge of analytical chemistry is the identification of unknown molecules based on tandem mass spectrometry (MS/MS) spectra. Current metabolite annotation approaches are often manual or partially automated, and commonly rely on a spectral database to search from or train machine learning classifiers on. Unfortunately, spectral databases are often instrument specific and incomplete due to the limited availability of compound standards or a molecular database, which limits the ability of methods utilizing them to predict novel molecule structures. We describe a generative modeling approach that can leverage the vast amount of unpaired and/or unlabeled molecule structures and MS/MS spectra to learn general rules for synthesizing molecule structures and MS/MS spectra. The approach is based on recent work using semi-supervised deep variational autoencoders to learn joint latent representations of multiple and complex modalities. We show that adding molecule structures with no spectra to the training set improves the prediction quality on spectra from a structure disjoint dataset of new molecules, which is not possible using bi-modal supervised approaches. The described methodology provides a demonstration and framework for how recent advances in semi-supervised machine learning can be applied to overcome bottlenecks in missing annotations and noisy data to tackle unaddressed problems in the life sciences where large volumes of data are available.


title: NanoNet: Rapid end-to-end nanobody modeling by deep learning at sub angstrom resolution

abstract: Antibodies are a rapidly growing class of therapeutics. Recently, single domain camelid VHH antibodies, and their recognition nanobody domain (Nb) appeared as a cost-effective highly stable alternative to full-length antibodies. There is a growing need for high-throughput epitope mapping based on accurate structural modeling of the variable domains that share a common fold and differ in the Complementarity Determining Regions (CDRs). We develop a deep learning end-to-end model, NanoNet, that given a sequence directly produces the 3D coordinates of the C[a] atoms of the entire VH domain. For the Nb test set, NanoNet achieves 1.7[A] overall average RMSD and 3.0[A] average RMSD for the most variable CDR3 loops. The accuracy for antibody VH domains is even higher: overall average RMSD < 1[A] and 2.2[A] RMSD for CDR3. NanoNet runtimes allow generation of ~1M nanobody structures in less than an hour on a standard CPU computer enabling high-throughput structure modeling.


title: ExhauFS: exhaustive search-based feature selection for classification and survival regression

abstract:
MotivationFeature selection is one of the main techniques used to prevent overfitting in machine learning applications. The most straightforward approach for feature selection is exhaustive search: one can go over all possible feature combinations and pick up the model with the highest accuracy. This method together with its optimizations were actively used in biomedical research, however, publicly available implementation is missing.

ResultsWe present ExhauFS - the user-friendly command-line implementation of the exhaustive search approach for classification and survival regression. Aside from tool description, we included three application examples in the manuscript to comprehensively review the implemented functionality. First, we executed ExhauFS on a toy cervical cancer dataset to illustrate basic concepts. Then, a multi-cohort microarray and RNA-seq breast cancer datasets were used to construct gene signatures for 5-year recurrence classification. Finally, Cox survival regression models were used to fit isomiR signatures for overall survival prediction for patients with colorectal cancer.

AvailabilitySource codes and documentation of ExhauFS are available on GitHub: https://github.com/s-a-nersisyan/ExhauFS.

Contactsnersisyan@hse.ru


title: Cross-organ transcriptomic comparison reveals universal factors during maturation

abstract: Various cell types can be derived from stem cells. However, these cells are immature and do not match their adult counterparts in functional capabilities, limiting their use in disease modeling and cell therapies. Thus, it is crucial to understand the mechanisms of maturation in vivo. However, it is unknown if there are genes and pathways conserved across organs during maturation. To address this, we performed a time-series analysis of the transcriptome of the mouse heart, brain, liver, and kidney and analyzed their trajectories over time. In addition, gene regulatory networks were reconstructed to determine overlapping expression patterns. Based on these, we identified commonly upregulated and downregulated pathways across all four organs. Key upstream regulators were also predicted based on the temporal expression of downstream genes. These findings suggest the presence of universal regulators during organ maturation, which may help us develop a general strategy to mature stem cell-derived cells in vitro.


title: Learning the mutational landscape of the cancer genome

abstract: An ongoing challenge to better understand and treat cancer is to distinguish neutral mutations, which do not affect tumor fitness, from those that provide a proliferative advantage. However, the variability of mutation rates has limited our ability to model patterns of neutral mutations and therefore identify cancer driver mutations. Here, we predict cancer-specific mutation rates genome-wide by leveraging deep neural networks to learn mutation rates within kilobase-scale regions and then refining these estimates to test for evidence of selection on combinations of mutations by comparing observed to expected mutation counts. We mapped mutation rates for 37 cancer types and used these maps to identify new putative drivers in understudied regions of the genome including cryptic alternative-splice sites, 5 untranslated regions and infrequently mutated genes. These results, available for exploration via web interface, indicate the potential for high-resolution neutral mutation models to empower further driver discovery as cancer sequencing cohorts grow.

category: cancer biology


title: Benchmarking topological accuracy of bacterial phylogenomic workflows using in silico evolution

abstract:
Phylogenetic analyses are widely used in microbiological research, for example to trace the progression of bacterial outbreaks based on whole-genome sequencing data. In practice, multiple analysis steps such as de novo assembly, alignment and phylogenetic inference are combined to form phylogenetic workflows. Comprehensive benchmarking of the accuracy of complete phylogenetic workflows is lacking.

To benchmark different phylogenetic workflows, we simulated bacterial evolution under a wide range of evolutionary models, varying the relative rates of substitution, insertion, deletion, gene duplication, gene loss and lateral gene transfer events. The generated datasets corresponded to a genetic diversity usually observed within bacterial species ([≥]95% average nucleotide identity). We replicated each simulation three times to assess replicability. In total, we benchmarked seventeen distinct phylogenetic workflows using 8 different simulated datasets.

We found that recently developed k-mer alignment methods such as kSNP and SKA achieve similar accuracy as reference mapping. The high accuracy of k-mer alignment methods can be explained by the large fractions of genomes these methods can align, relative to other approaches. We also found that the choice of de novo assembly algorithm influences the accuracy of phylogenetic reconstruction, with workflows employing SPAdes or SKESA outperforming those employing Velvet. Finally, we found that the results of phylogenetic benchmarking are highly variable between replicates.

We conclude that for phylogenomic reconstruction k-mer alignment methods are relevant alternatives to reference mapping at species level, especially in the absence of suitable reference genomes. We show de novo genome assembly accuracy to be an underappreciated parameter required for accurate phylogenomic reconstruction.

Impact statementPhylogenetic analyses are crucial to understand the evolution and spread of microbes. Among their many applications is the reconstruction of transmission events which can provide information on the progression of pathogen outbreaks. For example, to investigate foodborne outbreaks such as the 2011 outbreak of Escherichia coli O104:H4 across Europe. As different microbes evolve differently, it is important to know which phylogenetic workflows are most accurate when working with diverse bacterial data. However, benchmarks usually consider only a limited dataset. We therefore employed a range of simulated evolutionary scenarios and benchmarked seventeen phylogenetic workflows on these simulated datasets. An advantage of our simulation approach is that we know a priori what the outcome of the analyses should be, allowing us to benchmark accuracy. We found significant differences between phylogenetic workflows and were able to dissect which factors contribute to phylogenetic analysis accuracy. Taken together, this new information will hopefully enable more accurate phylogenetic analysis of bacterial outbreaks.

Data summaryA Zenodo repository is available at https://doi.org/10.5281/zenodo.5036179 containing all simulated genomes, all alignments produced by phylogenetic workflows and .csv files summarising the topological accuracies of phylogenies produced based on these alignments. Code is available at https://github.com/niekh-13/phylogenetic_workflows.


title: GlycoSHIELD: a versatile pipeline to assess glycan impact on protein structures

abstract: More than 75% of surface and secreted proteins are modified by covalent addition of complex sugars through N- and O-glycosylation. Unlike proteins, glycans do not typically adopt specific secondary structures and remain very mobile, influencing protein dynamics and interactions with other molecules. Glycan conformational freedom impairs complete structural elucidation of glycoproteins. Computer simulations may be used to model glycan structure and dynamics. However, such simulations typically require thousands of computing hours on specialized supercomputers, thus limiting routine use. Here, we describe a reductionist method that can be implemented on personal computers to graft ensembles of realistic glycan conformers onto static protein structures in a matter of minutes. Using this open-source pipeline, we reconstructed the full glycan cover of SARS-CoV-2 Spike protein (S-protein) and a human GABAA receptor. Focusing on S-protein, we show that GlycoSHIELD recapitulates key features of extended simulations of the glycosylated protein, including epitope masking, and provides new mechanistic insights on N-glycan impact on protein structural dynamics.


title: Contribution of single mutations to selected SARS-CoV-2 emerging variants Spike antigenicity

abstract: Towards the end of 2020, multiple variants of concern (VOCs) and variants of interest (VOIs) have arisen from the original SARS-CoV-2 Wuhan-Hu-1 strain. Mutations in the Spike protein are highly scrutinized for their impact on transmissibility, pathogenesis and vaccine efficacy. Here, we contribute to the growing body of literature on emerging variants by evaluating the impact of single mutations on the overall antigenicity of selected variants and their binding to the ACE2 receptor. We observe a differential contribution of single mutants to the global variants phenotype related to ACE2 interaction and antigenicity. Using biolayer interferometry, we observe that enhanced ACE2 interaction is mostly modulated by a decrease in off-rate. Finally, we made the interesting observation that the Spikes from tested emerging variants bind better to ACE2 at 37{degrees}C compared to the D614G variant. Whether improved ACE2 binding at higher temperature facilitates emerging variants transmission remain to be demonstrated.

category: microbiology


title: ComF is a key mediator in single-stranded DNA transport and handling during natural transformation

abstract: Natural transformation plays a major role in the spreading of antibiotic resistances and virulence factors. Whilst bacterial species display specificities in the molecular machineries allowing transforming DNA capture and integration into their genome, the ComF(C) protein is essential for natural transformation in all Gram-positive and - negative species studied. Despite this, its role remains largely unknown. Here, we show that Helicobacter pylori ComF is not only involved in DNA transport through the cell membrane, but it also required for the handling of the ssDNA once it is delivered into the cytoplasm. ComF crystal structure revealed the presence of a zinc-finger motif and a putative phosphoribosyl transferase domain, both necessary for its in vivo activity. ComF is a membrane-associated protein with affinity for single-stranded DNA. Collectively, our results suggest that ComF provides the link between the transport of the transforming DNA into the cytoplasm and its handling by the recombination machinery.


title: Opposing roles of p38α phosphorylation and arginine methylation in driving TDP-43 proteinopathy.

abstract: Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder typically characterized by insoluble inclusions of hyperphosphorylated TDP-43. The mechanisms underlying toxic TDP-43 accumulation are not understood. Persistent activation of p38 mitogen-activated protein kinase (MAPK) is implicated in ALS. However, it is unclear how p38 MAPK affects TDP-43 proteinopathy. Here, we demonstrate that inhibition of p38 MAPK reduces pathological TDP-43 phosphorylation, aggregation, cytoplasmic mislocalization, and neurotoxicity. We establish that p38 MAPK phosphorylates TDP-43 at pathological serine 409/410 (S409/S410) and serine 292 (S292), which reduces TDP-43 liquid-liquid phase separation (LLPS) but allows pathological TDP-43 aggregation. Moreover, we show that protein arginine methyltransferase 1 methylates TDP-43 at R293. Importantly, S292 phosphorylation reduces R293 methylation, and R293 methylation reduces S409/S410 phosphorylation. R293 methylation permits TDP-43 LLPS and reduces pathological TDP-43 aggregation. Thus, strategies to reduce p38-mediated TDP-43 phosphorylation and promote R293 methylation could have therapeutic utility for ALS and related TDP-43 proteinopathies.


title: Are members of the Anopheles fluviatilis complex conspecific?

abstract: Anopheles fluviatilis sensu lato, a primary malaria vector in India, was identified to be comprised of four cryptic species, provisionally designated as species S, T, U and V. However, Kumar et al. (Mol Ecol Resour, 2013;13:354-61) considered all of the then known three members of this species complex (S, T and U) conspecific. The specific status of species S and T was refuted based on the lack of sufficient barcode gap in mitochondrial-CO1 and the perceived presence of heterozygotes in populations as detected through one of the two species-specific PCR assays employed for the cryptic species identification. The existence of species U was refuted claiming that earlier investigations have already refuted their existence. This conclusion is concerning because of the differential public health implications of members of the Fluviatilis Complex. Here we discuss problems associated with the CO1-based barcode approach for delimitation of cryptic species, the perceived heterozygosity between species S and T based on a species-specific PCR assay, and interpretation of published reports. We demonstrated that fixed differences do exist in the ITS2-rDNA sequence of species S and T with no evidence of heterozygotes in sympatric populations and, that the observed heterozygosity by Kumar et al. in the ITS2-based species diagnostic PCR is due to the high mispriming tendency of the T-specific primer with species S. We infer that mitochondrial DNA-based  barcoding gap, an arbitrary threshold recommended for species delimitation, alone, is inadequate to delimit the members of An. fluviatilis complex.


title: Intermediary role of lung alveolar type 1 cells in epithelial repair upon Sendai virus infection

abstract: The lung epithelium forms the first barrier against respiratory pathogens and noxious chemicals; however, little is known about how >90% of this barrier - made of alveolar type 1 (AT1) cells - responds to injury, in contrast to our accumulating knowledge of epithelial progenitor and stem cells whose importance lies in their ability to restore the barrier. Using Sendai virus to model natural infection in mice, we combine 3D imaging, lineage-tracing, and single-cell genomics to show that AT1 cells have an intermediary role by persisting in areas depleted of alveolar type 2 (AT2) cells, mounting an interferon response, and receding from invading airway cells. Sendai virus infection mobilizes airway cells to form alveolar SOX2+ clusters without differentiating into AT1 or AT2 cells, as shown in influenza models. Intriguingly, large AT2-cell-depleted areas remain covered by AT1 cells, which we name "AT2-less regions", and are replaced by SOX2+ clusters spreading both basally and luminally around AT1 cell extensions. AT2 cell proliferation and differentiation are largely confined to topologically distal regions - the end of airspace that could be in the periphery or middle of the lung - and form de novo alveolar surface, with limited contribution to in situ repair of AT2-less regions. Time course single-cell RNA-seq and AT1-cell interactome analyses suggest enhanced recognition of AT1 cells by immune cells and altered growth signals. Our comprehensive spatiotemporal and genome-wide study highlights the hitherto unappreciated role of AT1 cells during Sendai virus infection and possibly other injury-repair processes.


title: Molecular basis of UHRF1 allosteric activation for synergistic histone modification binding by PI5P

abstract: Chromatin marks are recognized by distinct binding modules many of which are embedded in multidomain proteins or complexes. How the different protein functionalities of complex chromatin modulators are regulated is often unclear. Using a combination of biochemical, biophysical and structural approaches we delineated the regulation of the H3unmodified and H3K9me binding activities of the multidomain UHRF1 protein. The phosphoinositide PI5P interacts with two distant flexible linker regions of UHRF1 in a mode that is dependent on the polar head group and the acyl part of the phospholipid. The associated conformational rearrangements stably position the H3unmodified and H3K9me3 histone recognition modules of UHRF1 for multivalent and synergistic binding of the H3 tail. Our work highlights a novel molecular function for PI5P outside of the context of lipid mono- or bilayers and establishes a molecular paradigm for the allosteric regulation of complex, multidomain chromatin modulators by small cellular molecules.


title: Better together - Relative retention time plus spectral matching improves automated glycan characterization using PGC-nLC-IT-ESI-MS/MS

abstract: Porous Graphitized Carbon nano-liquid chromatography tandem mass spectrometry (PGC-nLC-MS/MS) is a glycomics technique with the unique capacity to differentiate isobaric glycans. The lack of suitable software tools integrating chromatography and MS-information delivered by PGC-nLC-MS/MS has been limiting fast and robust glycan identification and quantitation. We report a LC-system-independent strategy called GlycoRRT that combines relative retention time (RRT) and negative ion fragment spectra analyses for isobaric structure-specific glycomics of PGC-nLC-MS/MS data. The GlycoRRT toolset is fully customizable and easily adaptable enabling semi-automated high-throughput structural assignments. The current library contains over 200 entries and their individual meta-data (MS instrumentation, experimental conditions, retention times, fragmentation profiles and glycan structural diagnostic ion features) relevant for reliable data analyses. The GlycoRRT workflow was employed to map the N- and O-glycome in blood group matched human plasma and urine as well as decipher Immunoglobulin (IgG) glycosylation features from 13 different animal species. We have also developed visualization tools to enable a consistent, reliable, and reproducible analysis of large sets of multidimensional PGC-nLC-MS/MS glycomics data. This comprehensive glycan resource provides the glycan map of human and animal species, will serve as a reference in dissecting the role of glycans in host pathogen interaction and zoonotic disease transmission.


title: Analysis of 329,942 SARS-CoV-2 records retrieved from GISAID database

abstract:
BackgroundThe 31st of December 2019 was when the World Health Organization received a report about an outbreak of pneumonia of unknown etiology in the Chinese city of Wuhan. The outbreak was the result of the novel virus labeled as SARS-CoV-2, which spread to about 220 countries and caused approximately 3,311,780 deaths, infecting more than 159,319,384 people by May 12th, of 2021. The virus caused a worldwide pandemic leading to panic, quarantines, and lockdowns - although none of its predecessors from the coronavirus family have ever achieved such a scale. The key to understanding the global success of SARS-CoV-2 is hidden in its genome.

Materials and MethodsWe retrieved data for 329,942 SARS-CoV-2 records uploaded to the GISAID database from the beginning of the pandemic until the 8th of January 2021. To process the data, a Python variant detection script was developed, using pairwise2 from the BioPython library. Pandas, Matplotlib, and Seaborn, were applied to visualize the data. Genomic coordinates were obtained from the UCSC Genome Browser (https://genome.ucsc.edu/). Sequence alignments were performed for every gene separately. Genomes less than 26,000 nucleotides long were excluded from the research. Clustering was performed using HDBScan.

ResultsHere, we addressed the genetic variability of SARS-CoV-2 using 329,942 worldwide samples. The analysis yielded 155 genome variations (SNPs and deletions) in more than 0.3% of the sequences. Nine common SNPs were present in more than 20% of the samples. Clustering results suggested that a proportion of people (2.46%) were infected with a distinct subtype of the B.1.1.7 variant. The subtype may be characterized by four to six additional mutations, with four being a more frequent option (G28881A, G28882A, and G28883[C] in the N gene, A23403G in S, A28095T in ORF8, G25437T in ORF3a). Two clusters were formed by mutations in the samples uploaded predominantly by Denmark and Australia, which may indicate the emergence of "Danish" and "Australian" variants. Five clusters were linked to increased/decreased age, shifted gender ratio, or both. According to a correlation coefficient matrix, 69 mutations correlate with at least one other mutation (correlation coefficient greater than 0.7). We also addressed the completeness of the GISAID database, where between 77% and 93% of the fields were either left blank or filled incorrectly. Metadata mining analysis has led to a hypothesis about gender inequality in medical care in certain countries. Finally, we found ORF6 and E as the most conserved genes (96.15% and 94.66% of the sequences totally match the reference, respectively), making them potential targets for vaccines and treatment. Our results indicate areas of the SARS-CoV-2 genome that researchers can focus on for further structural and functional analysis.


title: Llamanade: an open-source computational pipeline for robust nanobody humanization

abstract:
Nanobodies (Nbs) have recently emerged as a promising class of antibody fragments for biomedical and therapeutic applications. Despite having marked physicochemical properties, Nbs are derived from camelids and may require "humanization" to improve translational potentials for clinical trials. Here we have systematically analyzed the sequence and structural properties of Nbs based on NGS (next-generation sequencing) databases and high-resolution structures. Our analysis reveals substantial framework diversities and underscores the key differences between Nbs and human Immunoglobulin G (IgG) antibodies. We identified conserved residues that may contribute to enhanced solubility, structural stability, and antigen-binding, providing insights into Nb humanization. Based on big data analysis, we developed "Llamanade, a user-friendly, open-source to facilitate rational humanization of Nbs. Using Nb sequence as input, Llamanade provides information on the sequence features, model structures, and optimizes solutions to humanize Nbs. The full analysis for a given Nb takes less than a minute on a local computer. To demonstrate the robustness of this tool, we applied it to successfully humanize a cohort of structurally diverse and highly potent SARS-CoV-2 neutralizing Nbs. Llamanade is freely available and will be easily accessible on a web server to support the development of a rapidly expanding repertoire of therapeutic Nbs into safe and effective trials.

Author SummaryCamelid Nbs are characterized by small size, excellent pharmacological properties and high flexibility in bioengineering for therapeutic development. However, Nbs are "xeno" antibodies, which require "humanization" to improve their translational potential. Currently, there is a lack of systematic investigation of Nbs to rationally guide humanization. No dedicated software has been developed for this purpose. Here, we report the development of Llamanade, an open-source computational pipeline and the first dedicated software to facilitate rational humanization of Nbs.

To subjectively evaluate Llamanade, we used it to humanize a cohort of structurally diverse and ultrapotent antiviral Nbs against SARS-CoV-2. Robust humanization by Llamanade significantly improved the humanness level of Nbs to closely resemble fully human IgGs. Importantly, these highly humanized antiviral Nbs remained excellent solubility and comparably high bioactivities to the non-humanized Nb precursors. We envision that Llamanade will help advance Nb research into therapeutic development.


title: Sequence dependence of transient Hoogsteen base-pairing in DNA

abstract: Hoogsteen (HG) base-pairing is characterized by a 180{degrees} rotation of the purine base with respect to the Watson-Crick-Franklin (WCF) motif. Recently, it has been found that both conformations coexist in a dynamical equilibrium and that several biological functions require HG pairs. This relevance has motivated experimental and computational investigations of the base-pairing transition. However, a systematic simulation of sequence variations has remained out of reach. Here, we employ advanced path-based methods to perform unprecedented free-energy calculations. Our methodology enables us to study the different mechanisms of purine rotation, either remaining inside or after flipping outside of the double helix. We study seven different sequences, which are neighbor variations of a well-studied A{middle dot}T pair in A6-DNA. We observe the known effect of A{middle dot}T steps favoring HG stability, and find evidence of triple-hydrogen-bonded neighbors hindering the inside transition. More importantly, we identify a dominant factor: the direction of the A rotation, with the 6-ring pointing either towards the longer or shorter segment of the chain, respectively relating to a lower or higher barrier. This highlights the role of DNAs relative flexibility as a modulator of the WCF/HG dynamic equilibrium. Additionally, we provide a robust methodology for future HG proclivity studies.