title: Comparative Genomics and Integrated Network Approach Unveiled Undirected Phylogeny Patterns, Co-mutational Hotspots, Functional Crosstalk and Regulatory Interactions in SARS-CoV-2

abstract:
SARS-CoV-2 pandemic resulted in 92 million cases in a span of one year. The study focuses on understanding population specific variations attributing its high rate of infections in specific geographical regions particularly in USA. Rigorous phylogenomic network analysis of complete SARS-CoV-2 genomes (245) inferred five central clades named a (ancestral), b, c, d and e (subtype e1 & e2). The clade d & e2 were found exclusively comprising of USA. Clades were distinguished by 10 co-mutational combinations in Nsp3, ORF8, Nsp13, S, Nsp12, Nsp2 and Nsp6. Our analysis revealed that only 67.46% of SNP mutations were at amino acid level. T1103P mutation in Nsp3 was predicted to increase protein stability in 238 strains except 6 strains which were marked as ancestral type; whereas co-mutation (P409L & Y446C) in Nsp13 were found in 64 genomes from USA highlighting its 100% co-occurrence. Docking highlighted mutation (D614G) caused reduction in binding of Spike proteins with ACE2, but it also showed better interaction with TMPRSS2 receptor contributing to high transmissibility among USA strains. We also found host proteins, MYO5A, MYO5B, MYO5C had maximum interaction with viral proteins (N, S, M). Thus, blocking the internalization pathway by inhibiting MYO5 proteins which could be an effective target for COVID-19 treatment. The functional annotations of the HPI network were found to be closely associated with hypoxia and thrombotic conditions confirming the vulnerability and severity of infection. We also screened CpG islands in Nsp1 & N conferring ability of SARS-CoV-2 to enter and trigger ZAP activity inside host cell.

ImportanceIn the current study we presented a global view of mutational pattern observed in SARS-CoV-2 virus transmission. This provided a who-infect-whom geographical model since the early pandemic. This is hitherto the most comprehensive comparative genomics analysis of full-length genomes for co-mutations at different geographical regions specially in USA strains. Compositional structural biology results suggested that mutations have balance of contrary forces effect on pathogenicity suggesting only few mutations to effective at translation level but not all. Novel HPI analysis and CpG predictions elucidates the proof of concept of hypoxia and thrombotic conditions in several patients. Thus, the current study focuses the understanding of population specific variations attributing high rate of SARS-CoV-2 infections in specific geographical regions which may eventually be vital for the most severely affected countries and regions for sharp development of custom-made vindication strategies.


title: Immune responses in pancreatic cancer may be restricted byprevalence of activated regulatory T-cells, dysfunctionalCD8+ T-cells, and senescent T-cells

abstract:
Pancreatic cancer has the worst prognosis of any human malignancy and leukocyte infiltration is a major prognostic marker of the disease. As current immunotherapies confer negligible survival benefits, there is a need to better characterise leukocytes in pancreatic cancer to identify better therapeutic strategies.

In this study, we analysed 32 human pancreatic cancer patients from two independent cohorts. A multi-parameter mass-cytometry analysis was performed on 32,000 T-cells from eight patients. Single-cell RNA sequencing dataset analysis was performed on a cohort of 24 patients. Multiplex immunohistochemistry imaging and spatial analysis were performed to map immune infiltration into the tumour microenvironment.

Regulatory T-cell populations demonstrated highly immunosuppressive states with high TIGIT, ICOS and CD39 expression. CD8+ T-cells were found to be either in senescence or an exhausted state. The exhausted CD8 T-cells had low PD-1 expression but high TIGIT and CD39 expression. These findings were corroborated in an independent pancreatic cancer single-cell RNA dataset from additional 24 patients.

These data suggest that T-cells are major players in the suppressive microenvironment of pancreatic cancer. Our work identifies novel therapeutic targets that should form the basis for rational design of a new generation of clinical trials in pancreatic ductal adenocarcinoma.

category: immunology


title: Preference for animate domain sounds in the fusiform gyrus of blind individuals is modulated by shape-action mapping

abstract: In high-level visual shape areas in the human brain, preference for inanimate objects is observed regardless of stimulation modality (visual/auditory/tactile) and subjects visual experience (sighted/blind individuals), whereas preference for animate entities seems robust only in the visual modality. Here, we test a hypothesis explaining this effect: visual shape representations can be reliably activated through different sensory modalities when they are systematically related to action system representations. We studied fMRI activations in congenitally blind and sighted subjects listening to animal, object, and human sounds. We found that, in blind individuals, the typical anatomical location of the fusiform face area responds to human facial expression sounds, with a clear mapping between the facial motor action and the resulting face shape, but not to speech or animal sounds. Using face areas activation in the blind subjects we could distinguish between specific facial expressions used in the study, but not between specific speech sounds. We conclude that auditory stimulation can reliably activate visual representations of those stimuli - inanimate or animate - for which shape and action computations are transparently related. Our study suggests that visual experience is not necessary for the development of functional preference for face-related information in the fusiform gyrus.

category: neuroscience


title: Simultaneously monitoring aquatic and riparian biodiversity using riverine water eDNA

abstract: Environmental DNA (eDNA) metabarcoding for biodiversity monitoring is a critical technical advance. Both aquatic and terrestrial biodiversity information can be detected in riverine water eDNA. However, it remains unverified whether riverine water eDNA can be used to simultaneously monitor aquatic and terrestrial biodiversity. Our specific objective was to assess the effectiveness of monitoring aquatic and riparian biodiversity using riverine water eDNA. We proposed that the monitoring effectiveness (the proportion of aquatic and terrestrial biodiversity information detected by riverine water eDNA samples) could be approximated by the transportation effectiveness of land-to-river and upstream-to-downstream biodiversity information flow. We conducted a case study in a watershed on the Qinghai-Tibet Plateau and estimated the effectiveness of using riverine water eDNA to monitor aquatic and riparian biodiversity based on comparing the operational taxonomic units (OTUs) and species assemblages of three taxonomic communities detected in riverine water eDNA samples and riparian soil eDNA samples in spring, summer, and autumn. The aquatic and riparian biodiversity of a watershed on the Qinghai-Tibet Plateau could be simultaneously effectively monitored using riverine water eDNA on summer or autumn rainy days. Monitoring bacterial communities was more efficient than monitoring eukaryotic communities. On summer rainy days, 43%-76% of riparian species could be detected in water eDNA samples, 92%-99% of upstream species could be detected in a 1-km downstream eDNA sample, and more than 50% of dead bioinformation (i.e., the bioinformation labeling the biological material without life activity and fertility) could be monitored 4-6 km downstream for eukaryotes and 13-19 km for bacteria. We encourage more studies on the monitoring effectiveness for each taxonomic community in other watersheds with different environmental conditions. We believe that in future ecological research, conservation and management, we could efficiently monitor and assess the aquatic and terrestrial biodiversity by simply using riverine water eDNA samples.

category: ecology


title: Low competitive status elicits aggression in healthy young men: behavioral and neural evidence

abstract: Winners are commonly assumed to compete more aggressively than losers. Here, we find overwhelming evidence for the opposite. We first demonstrate that low-ranking teams commit more fouls than they receive in top-tier soccer, ice hockey, and basketball mens leagues. We replicate this effect in the laboratory, showing that male participants deliver louder sound blasts to a rival when placed in a low-status position. Using neuroimaging, we characterize brain activity patterns that encode competitive status as well as those that facilitate status-dependent aggression in healthy young men. These analyses reveal three key findings. First, anterior hippocampus and striatum contain multivariate representations of competitive status. Second, interindividual differences in status-dependent aggression are linked with a sharper status differentiation in the striatum and with greater reactivity to status-enhancing victories in the dorsal anterior cingulate cortex. Third, activity in ventromedial, ventrolateral, and dorsolateral prefrontal cortex is associated with trial-wise increases in status-dependent aggressive behavior. Taken together, our results run counter to narratives glorifying aggression in competitive situations. Rather, we show that those in the lower ranks of skill-based hierarchies are more likely to behave aggressively and identify the potential neural basis of this phenomenon.


title: Lift observation conveys object weight distribution but partly enhances predictive lift planning

abstract:
Observation of object lifting allows updating of internal object representations for object weight, in turn enabling accurate scaling of fingertip forces when lifting the same object. Here, we investigated whether lift observation also enables updating of internal representations for an objects weight distribution. We asked participants to lift an inverted T-shaped manipulandum, of which the weight distribution could be changed, in turns with an actor. Participants were required to minimize object roll (i.e.,  lift performance) during lifting and were allowed to place their fingertips at self-chosen locations. The center of mass changed unpredictably every third to sixth trial performed by the actor and participants were informed that they would always lift the same weight distribution as the actor. Participants observed either erroneous (i.e., object rolling towards its heavy side) or skilled (i.e., minimized object roll) lifts. Lifting performance after observation was compared to lifts without prior observation and to lifts after active lifting, which provided haptic feedback about the weight distribution. Our results show that observing both skilled and erroneous lifts convey an objects weight distribution similar to active lifting, resulting in altered digit positioning strategies. However, minimizing object roll on novel weight distributions was only improved after observing error lifts and not after observing skilled lifts. In sum, these findings suggest that although observing motor errors and skilled motor performance enables updating of digit positioning strategy, only observing error lifts enables changes in predictive motor control when lifting objects with unexpected weight distributions.

New and noteworthyIndividuals are able to extract an objects size and weight by observing interactions with objects and subsequently integrate this information in their own motor repertoire. Here, we show that this ability extrapolates to weight distributions. Specifically, we highlighted that individuals can perceive an objects weight distribution during lift observation but can only partially embody this information when planning their own actions.


title: Patients with ventromedial prefrontal lesions show an implicit approach bias to angry faces

abstract: Damage to the ventromedial prefrontal cortex (VMPFC) can cause maladaptive social behavior, but the cognitive processes underlying these behavioral changes are still uncertain. Here, we tested whether patients with acquired VMPFC lesions show altered approach-avoidance tendencies to emotional facial expressions. Thirteen patients with focal VMPFC lesions and 31 age- and gender-matched healthy controls performed an implicit approach-avoidance task in which they either pushed or pulled a joystick depending on stimulus color. While controls avoided angry faces, VMPFC patients displayed an incongruent response pattern characterized by both increased approach and reduced avoidance of angry facial expressions. The approach bias was stronger in patients with higher self-reported impulsivity and disinhibition, and in those with larger lesions. We further used linear ballistic accumulator modelling to investigate latent parameters underlying approach-avoidance decisions. Controls displayed negative drift rates when approaching angry faces, whereas VMPFC lesions abolished this pattern. In addition, VMPFC patients had weaker response drifts than controls during avoidance. Finally, patients showed reduced drift rate variability and shorter non-decision times, indicating impulsive and rigid decision-making. Our findings thus suggest that VMPFC damage alters the pace of evidence accumulation in response to social signals, eliminating a default, protective avoidant bias and facilitating dysfunctional approach behavior.


title: Vegetative phase change in Populus tremula x alba

abstract:
O_LIPlants transition through juvenile and adult phases of vegetative development in a process known as vegetative phase change (VPC). In poplars (genus Populus) the differences between these stages are subtle, making it difficult to determine when this transition occurs. Previous studies of VPC in poplars have relied on plants propagated in vitro, leaving the natural progression of this process unknown.
C_LIO_LIWe examined developmental morphology of seed-grown and in vitro derived Populus tremula x alba (clone 717-1B4), and compared the phenotype of these, to transgenics with manipulated miR156 expression, the master regulator of VPC.
C_LIO_LIIn seed-grown plants, most traits changed from node-to-node during the first 3 months of development but remained constant after node 25. Many traits remained unchanged in clones over-expressing miR156, or were enhanced when miR156 was lowered, demonstrating their natural progression is regulated by the miR156/SPL pathway. The characteristic leaf fluttering of Populus is one of these miR156-regulated traits.
C_LIO_LIVegetative development in plants grown from culture mirrored that of seed-grown plants, allowing direct comparison between plants often used in research and those found in nature. These results provide a foundation for further research on the role of VPC in the ecology and evolution of this economically important genus.
C_LI

category: plant biology


title: Universal features of epidemic models under social distancing guidelines

abstract: Social distancing as a form of nonpharmaceutical intervention has been enacted in many countries as a form of mitigating the spread of COVID-19. There has been a large interest in mathematical modeling to aid in the prediction of both the total infected population and virus-related deaths, as well as to aid government agencies in decision making. As the virus continues to spread, there are both economic and sociological incentives to minimize time spent with strict distancing mandates enforced, and/or to adopt periodically relaxed distancing protocols, which allow for scheduled economic activity. The main objective of this study is to reduce the disease burden in a population, here measured as the peak of the infected population, while simultaneously minimizing the length of time the population is socially distanced, utilizing both a single period of social distancing as well as periodic relaxation. We derive a linear relationship among the optimal start time and duration of a single interval of social distancing from an approximation of the classic epidemic SIR model. Furthermore, we see a sharp phase transition region in start times for a single pulse of distancing, where the peak of the infected population changes rapidly; notably, this transition occurs well before one would intuitively expect. By numerical investigation of more sophisticated epidemiological models designed specifically to describe the COVID-19 pandemic, we see that all share remarkably similar dynamic characteristics when contact rates are subject to periodic or one-shot changes, and hence lead us to conclude that these features are universal in epidemic models. On the other hand, the nonlinearity of epidemic models leads to non-monotone behavior of the peak of infected population under periodic relaxation of social distancing policies. This observation led us to hypothesize that an additional single interval social distancing at a proper time can significantly decrease the infected peak of periodic policies, and we verified this improvement numerically. While synchronous quarantine and social distancing mandates across populations effectively minimize the spread of an epidemic over the world, relaxation decisions should not be enacted at the same time for different populations.

category: systems biology


title: Genetic influences on hub connectivity of the human connectome

abstract: Brain network hubs are both highly connected and highly inter-connected, forming a critical communication backbone for coherent neural dynamics. The mechanisms driving this organization are poorly understood. Using diffusion-weighted imaging in twins, we identify a major role for genes, showing that they preferentially influence connectivity strength between network hubs of the human connectome. Using transcriptomic atlas data, we show that connected hubs demonstrate tight coupling of transcriptional activity related to metabolic and cytoarchitectonic similarity. Finally, comparing over thirteen generative models of network growth, we show that purely stochastic processes cannot explain the precise wiring patterns of hubs, and that model performance can be improved by incorporating genetic constraints. Our findings indicate that genes play a strong and preferential role in shaping the functionally valuable, metabolically costly connections between connectome hubs.


title: miRbiom: Machine-Learning on Bayesian Causal Nets of RBP-miRNA interactions successfully predicts miRNA profiles

abstract: Formation of mature miRNAs and their expression is a highly controlled process. It is very much dependent upon the post-transcriptional regulatory events. Recent findings suggest that several RNA binding proteins beyond Drosha/Dicer are involved in the processing of miRNAs. Deciphering of conditional networks for these RBP-miRNA interactions may help to reason the spatio-temporal nature of miRNAs which can also be used to predict miRNA profiles. In this direction, >25TB of data from different platforms were studied (CLIP-seq/RNA-seq/miRNA-seq) to develop Bayesian causal networks capable of reasoning miRNA biogenesis. The networks ably explained the miRNA formation when tested across a large number of conditions and experimentally validated data. The networks were modeled into an XGBoost machine learning system where expression information of the network components was found capable to quantitatively explain the miRNAs formation levels and their profiles. The models were developed for 1,204 human miRNAs whose accurate expression level could be detected directly from the RNA-seq data alone without any need of doing separate miRNA profiling experiments like miRNA-seq or arrays. A first of its kind, miRbiom performed consistently well with high average accuracy (91%) when tested across a large number of experimentally established data from several conditions. It has been implemented as an interactive open access web-server where besides finding the profiles of miRNAs, their downstream functional analysis can also be done. miRbiom will help to get an accurate prediction of human miRNAs profiles in the absence of profiling experiments and will be an asset for regulatory research areas. The study also shows the importance of having RBP interaction information in better understanding the miRNAs and their functional projectiles where it also lays the foundation of such studies and software in future.

category: bioinformatics


title: COPII Sec23 proteins form isoform-specific ER exit sites with differential effects on polarized growth

abstract: COPII, a coat of proteins that form vesicles on the ER, mediates vesicle traffic from the ER to the Golgi. In contrast to metazoans that have few genes encoding each COPII component, plants have expanded these gene families leading to the hypothesis that plant COPII has functionally diversified. Here, we analyzed the gene families encoding for the Sec23/24 heterodimer in the moss Physcomitrium (Physcomitrella) patens. In P. patens, Sec23 and Sec24 gene families are each comprised of seven genes. Silencing the Sec23/24 genes revealed isoform specific contributions to polarized growth, with the closely related Sec23D/E and Sec24C/D essential for protonemal development. Focusing on the Sec23 gene family, we discovered that loss of Sec23D alters ER morphology, increases ER stress, inhibits trafficking to the Golgi and to the plasma membrane in tip growing protonemata. In contrast, the remaining five Sec23 genes are dispensable for tip growth. While Sec23A/B/C/F/G do not quantitatively affect ER to Golgi trafficking in protonemata, they do contribute to secretion to the plasma membrane. Of the three highly expressed Sec23 isoforms in protonemata, Sec23G forms ER exit sites that are larger than Sec23B and Sec23D and do not overlap with Sec23D. Furthermore, ER exit sites labeled by Sec23B or Sec23G form in the absence of Sec23D. These data suggest that Sec23D/E form unique ER exit sites contributing to secretion that is essential for tip growing protonemata.

category: cell biology


title: Maternal RND3/RhoE deficiency impairs placental mitochondrial function in preeclampsia by modulating PPARγ-UCP2 cascade

abstract: Preeclampsia (PE) is a life-threatening disease of pregnant women associated with severe hypertension, proteinuria, or multi-organ injuries. Mitochondrial-mediated placental oxidative stress plays a key role in the pathogenesis of PE. However, the underlying mechanism remains to be revealed. Here, we identify Rnd3, a small Rho GTPase, regulating placental mitochondrial reactive oxygen species (ROS). We showed that Rnd3 is down-regulated in primary trophoblasts isolated from PE patients. Loss of Rnd3 in trophoblasts resulted in excessive ROS generation, cell apoptosis, mitochondrial injury, and proton leakage from the respiratory chain. Moreover, Rnd3 overexpression partially rescues the mitochondrial defects and oxidative stress in human PE primary trophoblasts. Rnd3 physically interacts with the peroxisome proliferators-activated receptor {gamma} (PPAR{gamma}) and promotes the PPAR{gamma}-mitochondrial uncoupling protein 2 (UCP2) cascade. Forced expression of PPAR{gamma} rescues deficiency of Rnd3-mediated mitochondrial dysfunction. We conclude that Rnd3 acts as a novel protective factor in placental mitochondria through PPAR{gamma}-UCP2 signaling and highlight that downregulation of Rnd3 is a potential factor involved in PE pathogenesis.

category: pathology


title: Organization of DNA replication origin firing in Xenopus egg extracts : the role of intra-S checkpoint

abstract:
During cell division, the duplication of the genome starts at multiple positions called replication origins. Origin firing requires the interaction of rate-limiting factors with potential origins during the S(ynthesis)-phase of the cell cycle. Origins fire as synchronous clusters which is proposed to be regulated by the intra-S checkpoint. By modelling the unchallenged, the checkpoint-inhibited and the checkpoint protein Chk1 over-expressed replication pattern of single DNA molecules from Xenopus sperm chromatin replicated in egg extracts, we demonstrate that the quantitative modelling of data requires: 1) a segmentation of the genome into regions of low and high probability of origin firing; 2) that regions with high probability of origin firing escape intra-S checkpoint regulation and 3) the variability of the rate of DNA synthesis close to replication forks is a necessary ingredient that should be taken in to account in order to describe the dynamic of replication origin firing. This model implies that the observed origin clustering emerges from the apparent synchrony of origin firing in regions with high probability of origin firing and challenge the assumption that the intra-S checkpoint is the main regulator of origin clustering.

Author summaryDNA replication is one of the fundamental cell functions. The genome of eukaryotic organisms is duplicated from multiple positions named replication origins. Single molecule experiments allow to visualise the dynamics of spatio-temporal patterns created during replication process. The dynamic of replication process is regulated by checkpoints. However, the influence and the role of checkpoint regulation in the dynamics of spatio-temporal patterns of DNA replication is not understood. In this work we build a minimal, process-based and data rooted numerical model that allows to decipher the impact of checkpoint regulation on the dynamics of spatio-temporal pattern of DNA replication.


title: The Biomechanical Basis of Biased Epithelial Tube Elongation in Lung and Kidney Development

abstract:
During lung development, epithelial branches expand preferentially in longitudinal direction. This bias in outgrowth has been linked to a bias in cell shape and in the cell division plane. How this bias arises is unknown. Here, we show that biased epithelial outgrowth occurs independent of the surrounding mesenchyme, of preferential turn-over of the extracellular matrix at the bud tips, and of FGF signalling. There is also no evidence for actin-rich filopodia at the bud tips. Rather, we find epithelial tubes to be collapsed during early lung and kidney development, and we observe fluid flow in the narrow tubes. By simulating the measured fluid flow inside segmented narrow epithelial tubes, we show that the shear stress levels on the apical surface are sufficient to explain the reported bias in cell shape and outgrowth. We use a cell-based vertex model to confirm that apical shear forces, unlike constricting forces, can give rise to both the observed bias in cell shapes and tube elongation. We conclude that shear stress may be a more general driver of biased tube elongation beyond its established role in angiogenesis.

Summary StatementWe systematically analysed the requirements for biased elongating outgrowth of epithelial tubes during lung and kidney development, and provide evidence that fluid-flow induced shear stress drives their biased elongation.

category: developmental biology


title: Movement signaling in ventral pallidum and dopaminergic midbrain is gated by behavioral state in singing birds

abstract:
Movement-related neuronal discharge in ventral tegmental area (VTA) and ventral pallidum (VP) is inconsistently observed across studies. One possibility is that some neurons are movement-related and others are not. Another possibility is that the precise behavioral conditions matter - that a single neuron can be movement related under certain behavioral states but not others. We recorded single VTA and VP neurons in birds transitioning between singing and non-singing states, while monitoring body movement with microdrive-mounted accelerometers. Many VP and VTA neurons exhibited body movement-locked activity exclusively when the bird was not singing. During singing, VP and VTA neurons could switch off their tuning to body movement and become instead precisely time-locked to specific song syllables. These changes in neuronal tuning occurred rapidly at state boundaries. Our findings show that movement-related activity in limbic circuits can be gated by behavioral context.

Significance statementNeural signals in the limbic system have long been known to represent body movements as well as reward. Here we show that single neurons dramatically change their tuning to movements when a bird starts to sing.


title: Tracking the contribution of inductive bias to individualized internal models

abstract: Internal models capture the regularities of the environment and are central to understanding how humans adapt to environmental statistics. In general, the correct internal model is unknown to observers, instead approximate and transient ones are recruited. However, experimenters assume an ideal observer model, which captures stimulus structure but ignores the diverging hypotheses that humans form during learning. We combine non-parametric Bayesian methods and probabilistic programming to infer rich and dynamic individualised internal models from response times in an implicit visuomotor sequence learning task. We identify two contributors to the internal model: the ideal observer model and a Markov model capturing only immediate temporal dependencies between observations. Individual learning curves revealed internal models initially dominated by the Markov model, which was later traded-off with the ideal observer model. Thus, our results reveal a structured inductive bias that varies across individuals both in strength and persistence but is consistent in overall structure.


title: Phase response analyses support a relaxation oscillator model of locomotor rhythm generation in Caenorhabditis elegans

abstract: Neural circuits coordinate with muscles and sensory feedback to generate motor behaviors appropriate to an animals environment. In C. elegans, the mechanisms by which the motor circuit generates undulations and modulates them based on the environment are largely unclear. We quantitatively analyzed C. elegans locomotion during free movement and during transient optogenetic muscle inhibition. Undulatory movements were highly asymmetrical with respect to the duration of bending and unbending during each cycle. Phase response curves induced by brief optogenetic inhibition of head muscles showed gradual increases and rapid decreases as a function of phase at which the perturbation was applied. A relaxation oscillator model based on proprioceptive thresholds that switch the active muscle moment was developed and is shown to quantitatively agree with data from free movement, phase responses, and previous results for gait adaptation to mechanical loadings. Our results suggest a neuromuscular mechanism underlying C. elegans motor pattern generation within a compact circuit.


title: Conformational entropy limits the transition from nucleation to elongation in amyloid aggregation

abstract:
The formation of {beta}-sheet rich amyloid fibrils in Alzheimers disease and other neurodegenerative disorders is limited by a slow nucleation event. To understand the initial formation of {beta}-sheets from disordered peptides, we used all-atom simulations to parameterize a lattice model that treats each amino acid as a binary variable with {beta} and non-{beta} states. We show that translational and conformational entropy give the nascent {beta}-sheet an anisotropic surface tension which can be used to describe the nucleus with two-dimensional Classical Nucleation Theory. Since translational entropy depends on concentration, the aspect ratio of the critical {beta}-sheet changes with protein concentration. Our model explains the transition from the nucleation phase to elongation as the point where the {beta}-sheet core becomes large enough to overcome the conformational entropy cost to straighten the terminal molecule. At this point the {beta}-strands in the nucleus spontaneously elongate, which results in a larger binding surface to capture new molecules. These results suggest that nucleation is relatively insensitive to sequence differences in co-aggregation experiments because the nucleus only involves a small portion of the peptide.

SIGNIFICANCEThe conversion of soluble proteins to amyloid aggregates is associated with many neurodegenerative diseases. Experiments have shown that this conversion occurs by a slow nucleation step followed by rapid growth. This work identifies the principle contributions to the free energy barrier that separates these two stages. It also shows how factors like protein concentration, sidechain interactions, and interactions with the environment can modify the barrier and affect nucleation times.

category: biophysics


title: Probing machine-learning classifiers using noise, bubbles, and reverse correlation

abstract:
BackgroundMany scientific fields now use machine-learning tools to assist with complex classification tasks. In neuroscience, automatic classifiers may be useful to diagnose medical images, monitor electrophysiological signals, or decode perceptual and cognitive states from neural signals. However, such tools often remain black-boxes: they lack interpretability. A lack of interpretability has obvious ethical implications for clinical applications, but it also limits the usefulness of these tools to formulate new theoretical hypotheses.

New methodWe propose a simple and versatile method to help characterize the information used by a classifier to perform its task. Specifically, noisy versions of training samples or, when the training set is unavailable, custom-generated noisy samples, are fed to the classifier. Multiplicative noise, so-called "bubbles", or additive noise are applied to the input representation. Reverse correlation techniques are then adapted to extract either the discriminative information, defined as the parts of the input dataset that have the most weight in the classification decision, and represented information, which correspond to the input features most representative of each category.

ResultsThe method is illustrated for the classification of written numbers by a convolutional deep neural network; for the classification of speech versus music by a support vector machine; and for the classification of sleep stages from neurophysiological recordings by a random forest classifier. In all cases, the features extracted are readily interpretable.

Comparison with Existing MethodsQuantitative comparisons show that the present method can match state-of-the art interpretation methods for convolutional neural networks. Moreover, our method uses an intuitive and well-established framework in neuroscience, reverse correlation. It is also generic: it can be applied to any kind of classifier and any kind of input data.

ConclusionsWe suggest that the method could provide an intuitive and versatile interface between neuroscientists and machine-learning tools.

HighlightsO_LIThe heuristics of black-box classifiers can be probed with noisy inputs
C_LIO_LIThe relevant features can be visualised in the input representation space
C_LIO_LIThe method applies to any kind of data such as 2D images or 1D time series
C_LIO_LIIt applies to any classifier such as deep neural networks, support vector machines, random forests
C_LI


title: An expanding manifold of transmodal regions characterizes adolescent reconfigurations of structural connectome organization

abstract:
AO_SCPLOWBSTRACTC_SCPLOWAdolescence is a critical time for the continued maturation of brain networks. Here, we assessed structural connectome development in a large longitudinal sample ranging from childhood to young adulthood. By projecting high-dimensional connectomes into compact manifold spaces, we identified a marked expansion of structural connectomes with the strongest effects in transmodal regions during adolescence. Findings reflected increased within-module connectivity together with increased segregation, indicating increasing differentiation of higher-order association networks from the rest of the brain. Projection of subcortico-cortical connectivity patterns into these manifolds showed parallel alterations in pathways centered on the caudate and thalamus. Connectome findings were contextualized via spatial transcriptome association analysis, highlighting genes enriched in cortex, thalamus, and striatum. Statistical learning of cortical and subcortical manifold features at baseline and their maturational change predicted measures of intelligence at follow-up. Our findings demonstrate that connectome manifold learning can bridge the conceptual and empirical gaps between macroscale network reconfigurations, microscale processes, and cognitive outcomes in adolescent development.

IO_SCPLOWMPACTC_SCPLOWManifold learning of longitudinal brain network data provides novel insights into adolescent structural connectome maturation, and how multiple scales of cortical and subcortical organization interact in typical neurodevelopment.


title: EMC2: A versatile algorithm for robust tracking of calcium dynamics from individual neurons in behaving animals

abstract: Measuring the activity of neuronal populations with calcium imaging can capture emergent functional properties of neuronal circuits with single cell resolution. However, the motion of freely behaving animals, together with the intermittent detectability of calcium sensors, can hinder automatic long-term monitoring of neuronal activity and the subsequent functional characterization of neural circuits. We report the development and open-source implementation of a multi-step cellular tracking algorithm (Elastic Motion Correction and Concatenation or EMC2) that compensates for the intermittent disappearance of moving neurons by integrating local deformation information from detectable neurons. We demonstrate the accuracy and versatility of our algorithm using calcium imaging data from two-photon volumetric microscopy in visual cortex of awake mice, and from confocal microscopy in behaving Hydra, which experiences major body deformation during its contractions. We quantify the performance of our algorithm using ground truth manual tracking of neurons, along with synthetic time-lapse sequences, covering a wide range of particle motions and detectability parameters. By monitoring the calcium activity of the same neuronal populations along several days in layer 2/3 of visual cortex of mice, we report important turn-over within the active neurons across days, with only few neurons that remained active across days. Then, by combining automatic tracking of single neuron activity over long time-lapse sequences with statistical clustering, we characterize and map neuronal ensembles in behaving Hydra. We document the existence three major non-overlapping ensembles of neurons (CB, RP1 and RP2) whose activity correlates with contractions and elongations. Our results prove that the EMC2 algorithm can be used as a robust platform for neuronal tracking in behaving animals.


title: Evolution and spread of SARS-CoV-2 likely to be affected by climate

abstract:
COVID-19 pandemic has been extensively studied by many researchers. However, it is still unclear why it was restricted to higher latitudes during the initial days and later cascaded in the tropics. Here, we analyzed 176 SARS-CoV-2 genomes across different latitudes and climate (Koppens climate) that provided insights about within species virus evolution and its relation to abiotic factors. Two genetically variant groups, named as G1 and G2 were identified, well defined by four mutations. The G1 group (ancestor), is mainly restricted to warm and moist, temperate climate (Koppens C climate) while its descendent G2 group surpasses the climatic restrictions of G1, initially cascading into neighboring cold climate (D) of higher latitudes and later into hot climate of the tropics (A). It appears that the gradation of temperate climate (Cfa-Cfb) to cold climate (Dfa-Dfb) drives the evolution of G1 into G2 variant group which later adapted to tropical climate (A) as well. It seems this virus followed inverse latitudinal gradient in the beginning due to its preference towards temperate (C) and cold climate (D). Our work elucidates virus evolutionary studies combined with climatic studies can provide crucial information about the pathogenesis and natural spreading pathways in such outbreaks which is hard to achieve through individual studies. Mutational insights gained may help design an efficacious vaccine.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=107 SRC="FIGDIR/small/147074v3_ufig1.gif" ALT="Figure 1">
View larger version (28K):
org.highwire.dtl.DTLVardef@f39a3org.highwire.dtl.DTLVardef@1910598org.highwire.dtl.DTLVardef@6516corg.highwire.dtl.DTLVardef@cbb48c_HPS_FORMAT_FIGEXP  M_FIG C_FIG In BriefThe authors elucidate adaptation of SARS-CoV-2 to different climates by studying phylogenetics and the distribution of strains on Koppens climate map.

HighlightsO_LIPhylogenetic network divides SARS-CoV-2 strains into two variant groups, G1 and G2.
C_LIO_LIG1 strains is restricted to Koppens "temperate" climate (mainly Cfa-Cfb).
C_LIO_LIG2 strains has evolved from G1 to sustain in other climates mainly "humid-continental" (Dfa-Dfb) and "tropical-savannah" (Aw) climate.
C_LI

category: evolutionary biology


title: CD209L/L-SIGN and CD209/DC-SIGN act as receptors for SARS-CoV-2

abstract: As the COVID-19 pandemic continues to spread, investigating the processes underlying the interactions between SARS-CoV-2 and its hosts is of high importance. Here, we report the identification of CD209L/L-SIGN and the related protein CD209/DC-SIGN as receptors capable of mediating SARS-CoV-2 entry into human cells. Immunofluorescence staining of human tissues revealed prominent expression of CD209L in the lung and kidney epithelium and endothelium. Multiple biochemical assays using a purified recombinant SARS-CoV-2 spike receptor binding domain (S-RBD) or S1 encompassing both NTB and RBD and ectopically expressed CD209L and CD209 revealed that CD209L and CD209 interact with S-RBD. CD209L contains two N-glycosylation sequons, at sites N92 and N361, but we determined that only site N92 is occupied. Removal of the N-glycosylation at this site enhances the binding of S-RBD with CD209L. CD209L also interacts with ACE2, suggesting a role for heterodimerization of CD209L and ACE2 in SARS-CoV-2 entry and infection in cell types where both are present. Furthermore, we demonstrate that human endothelial cells are permissive to SARS-CoV-2 infection and interference with CD209L activity by knockdown strategy or with soluble CD209L inhibits virus entry. Our observations demonstrate that CD209L and CD209 serve as alternative receptors for SARS-CoV-2 in disease-relevant cell types, including the vascular system. This property is particularly important in tissues where ACE2 has low expression or is absent, and may have implications for antiviral drug development.


title: Epstein-Barr virus tegument protein BGLF2 in exosomes released from virus-producer cells assists de novo infection

abstract: Viruses must adapt to the environment of their host cells to establish infection and persist. Diverse mammalian cells, including virus-infected cells, secrete extracellular vesicles such as exosomes containing proteins and miRNAs, and use these vesicles to mediate intercellular communications. However, the roles of exosomes in viral infection remain unclear. Here we screened viral proteins to identify those responsible for the exosome-mediated upregulation of Epstein-Barr virus (EBV) infection. We found BGLF2 protein encapsulated in exosomes, which were released from EBV-infected cells. BGLF2 protein is a tegument protein that exists the space between the envelope and the nucleocapsid, and it is released into the cytoplasm shortly after infection. BGLF2 protein-containing exosomes enhanced viral gene expression and repressed innate immunity, thereby assisting the EBV infection. In summary, the EBV tegument protein BGLF2 is encapsulated not only encapsulated in viral particles, but also in exosomes secreted from infected cells. Therefore, BGLF2 may play a crucial role in establishing EBV latent infection.

category: microbiology


title: Single-cell analysis reveals the pan-cancer invasiveness-associated transition of adipose-derived stromal cells into COL11A1-expressing cancer-associated fibroblasts

abstract:
During the last ten years, many research results have been referring to a particular type of cancer-associated fibroblasts associated with poor prognosis, invasiveness, metastasis and resistance to therapy in multiple cancer types, characterized by a gene expression signature with prominent presence of genes COL11A1, THBS2 and INHBA. Identifying the underlying biological mechanisms responsible for their creation may facilitate the discovery of targets for potential pan-cancer therapeutics. Using a novel computational approach for single-cell gene expression data analysis identifying the dominant cell populations in a sequence of samples from patients at various stages, we conclude that these fibroblasts are produced by a pan-cancer cellular transition originating from a particular type of adipose-derived stromal cells naturally present in the stromal vascular fraction of normal adipose tissue, having a characteristic gene expression signature. Focusing on a rich pancreatic cancer dataset, we provide a detailed description of the continuous modification of the gene expression profiles of cells as they transition from APOD-expressing adipose-derived stromal cells to COL11A1-expressing cancer-associated fibroblasts, identifying the key genes that participate in this transition. These results also provide an explanation to the well-known fact that the adipose microenvironment contributes to cancer progression.

Author summaryComputational analysis of rich gene expression data at the single-cell level from cancer biopsies can lead to biological discoveries about the nature of the disease. Using a computational methodology that identifies the gene expression profile of the dominant cell population for a particular cell type in the microenvironment of tumors, we observed that there is a remarkably continuous modification of this profile among patients, corresponding to a cellular transition. Specifically, we found that the starting point of this transition has a unique characteristic signature corresponding to cells that are naturally residing in normal adipose tissue. We also found that the endpoint of the transition has another characteristic signature corresponding to a particular type of cancer-associated fibroblasts with prominent expression of gene COL11A1, which has been found strongly associated with invasiveness, metastasis and resistance to therapy in multiple cancer types. Our results provide an explanation to the well-known fact that the adipose tissue contributes to cancer progression, shedding light on the biological mechanism by which tumor cells interact with the adipose microenvironment. We provide a detailed description of the changing profile during the transition, identifying associated genes as potential targets for pan-cancer therapeutics inhibiting the underlying mechanism.

category: cancer biology


title: Spatiotemporal whole-brain dynamics of auditory patterns recognition

abstract: Music is a non-verbal human language, built on logical structures and articulated in balanced hierarchies between sounds, offering excellent opportunities to explore how the brain creates meaning for complex spatiotemporal auditory patterns. Using the high temporal resolution of magnetoencephalography in 70 participants, we investigated their unfolding brain dynamics during the recognition of previously memorized J.S. Bachs musical patterns from prelude in C minor BWV 847 compared to novel patterns matched in terms of entropy and information content. Remarkably, the recognition of the memorized music ignited a widespread brain network comprising primary auditory cortex, superior temporal gyrus, insula, frontal operculum, cingulate gyrus, orbitofrontal cortex, basal ganglia, thalamus and hippocampus. Furthermore, measures of both brain activity and functional connectivity presented an overall increase over time, following the evolution and unfolding of the memorized musical patterns. Specifically, while the auditory cortex responded mainly to the first tones of the patterns, the activity and synchronization of higher-order brain areas such as cingulate, frontal operculum, hippocampus and orbitofrontal cortex largely increased over time, arguably representing the key whole-brain mechanisms for conscious recognition of auditory patterns as predicted by the global neuronal workspace hypothesis. In conclusion, our study described the fine-grained whole-brain activity and functional connectivity dynamics responsible for processing and recognition of previously memorized music. Further, the study highlights how the use of musical patterns in combination with a wide array of analytical tools and neuroscientific measures spanning from decoding to fast neural phase synchronization can shed new light on meaningful, complex cognitive processes.


title: Perturbed DNA methylation by sustained overexpression of Gadd45b induces chromatin disorganization, DNA strand breaks and dopaminergic neuron death in mice

abstract: Heterochromatin disorganization is a key hallmark of aging and DNA methylation state is currently the main molecular predictor of chronological age. The most frequent neurodegenerative diseases like Parkinson disease and Alzheimers disease are age-related but how the aging process and chromatin alterations are linked to neurodegeneration is unknown. Here, we investigated the consequences of viral overexpression of Gadd45b, a multifactorial protein involved in active DNA demethylation, in the midbrain of wild-type mice. Gadd45b overexpression induces global and stable changes in DNA methylation, particularly on gene bodies of genes related to neuronal functions. DNA methylation changes were accompanied by perturbed H3K9me3-marked heterochromatin and increased DNA damage. Prolonged Gadd45b expression resulted in dopaminergic neuron degeneration accompanied by altered expression of candidate genes related to heterochromatin maintenance, DNA methylation or Parkinson disease. Gadd45b overexpression rendered midbrain dopaminergic neurons more vulnerable to acute oxidative stress. Heterochromatin disorganization and DNA demethylation resulted in derepression of mostly young LINE-1 transposable elements, a potential source of DNA damage, prior to Gadd45b-induced neurodegeneration. Our data implicate that alterations in DNA methylation and heterochromatin organization, LINE-1 derepression and DNA damage can represent important contributors in the pathogenic mechanisms of dopaminergic neuron degeneration with potential implications for Parkinson disease.


title: A rapid, highly sensitive and open-access SARS-CoV-2 detection assay for laboratory and home testing

abstract: Global efforts to combat the Covid-19 pandemic caused by SARS-CoV-2 still heavily rely on RT-qPCR-based diagnostic tests. However, their high cost, moderate throughput and reliance on sophisticated equipment limit widespread implementation. Loop-mediated isothermal amplification after reverse transcription (RT-LAMP) is an alternative detection method that has the potential to overcome these limitations. We present a rapid, robust, sensitive and versatile RT-LAMP based SARS-CoV-2 detection assay. Our forty-minute procedure bypasses a dedicated RNA isolation step, is insensitive to carry-over contamination, and uses a hydroxynaphthol blue (HNB)-based colorimetric readout, which allows robust SARS-CoV-2 detection from various sample types. Based on this assay, we have substantially increased sensitivity and scalability by a simple nucleic acid enrichment step (bead-LAMP), established a pipette-free version for home testing (HomeDip-LAMP), and developed open source enzymes that can be produced in any molecular biology setting. Our advanced, universally applicable RT-LAMP assay is a major step towards population-scale SARS-CoV-2 testing.

category: molecular biology


title: A whole-cortex probabilistic diffusion tractography connectome

abstract:
The WU-Minn Human Connectome Project (HCP) is a publicly-available dataset containing state-of-art structural, functional, and diffusion-MRI for over a thousand healthy subjects. While the planned scope of the HCP included an anatomical connectome, resting-state functional-MRI forms the bulk of the HCPs current connectomic output. We address this by presenting a full-cortex connectome derived from probabilistic diffusion tractography and organized into the HCP-MMP1.0 atlas. Probabilistic methods and large sample sizes are preferable for whole-connectome mapping as they increase the fidelity of traced low-probability connections. We find that overall, connection strengths are lognormally distributed and decay exponentially with tract length, that connectivity reasonably matches macaque histological tracing in homologous areas, that contralateral homologs and left-lateralized language areas are hyperconnected, and that hierarchical similarity influences connectivity. We compare the diffusion-MRI connectome to existing resting-state fMRI and cortico-cortico evoked potential connectivity matrices and find that it is more similar to the latter. This work helps fulfill the promise of the HCP and will make possible comparisons between the underlying structural connectome and functional connectomes of various modalities, brain states, and clinical conditions.

Significance StatementThe tracts between cortical parcels can be estimated from diffusion MRI, but most studies concentrate on only the largest connections. Here we present an atlas, the largest and most detailed of its kind, showing connections among all cortical parcels. Connectivity is relatively enhanced between frontotemporal language areas and homologous contralateral locations. We find that connectivity decays with fiber tract distance more slowly than predicted by brain volume and that structural and stimulation-derived connectivity are more similar to each other than to resting-state functional MRI correlations. The connectome presented is publicly available and organized into a commonly used scheme for defining brain areas in order to enable ready comparison to other brain imaging datasets of various modalities.


title: The crystal structure of the Ca2+-ATPase 1 from Listeria monocytogenes reveals a pump primed for dephosphorylation

abstract: Many bacteria export intracellular calcium using active transporters homologous to the sarco/endoplasmic reticulum Ca2+-ATPase (SERCA). Here we present three crystal structures of Ca2+-ATPase 1 from Listeria monocytogenes (LMCA1). Structures with BeF3- mimicking a phosphoenzyme state reveal a closed state, which is intermediate between the outward-open E2P and the proton-occluded E2-P* conformations known for SERCA. It suggests that LMCA1 in the E2P state is pre-organized for dephosphorylation upon Ca2+ release, consistent with the rapid dephosphorylation observed in single-molecule studies. An arginine side-chain occupies the position equivalent to calcium binding site I in SERCA, leaving a single Ca2+-binding site in LMCA1, corresponding to SERCA site II. Observing no putative transport pathways dedicated to protons, we infer a direct proton counter transport through the Ca2+ exchange pathways. The LMCA1 structures provide insight into the evolutionary divergence and conserved features of this important class of ion transporters.

category: biochemistry


title: De novo deletions and duplications at recombination hotspots in mouse germlines

abstract: Numerous DNA double-strand breaks (DSBs) arise during meiosis to initiate homologous recombination. These DSBs are usually repaired faithfully, but here we uncover a new type of mutational event in which deletions form via joining of ends from two closely spaced DSBs (double cuts) within a single hotspot or at adjacent hotspots on the same or different chromatids. Deletions occur in normal meiosis but are much more frequent when DSB formation is dysregulated in the absence of the ATM kinase. Events between chromosome homologs point to multi-chromatid damage and aborted gap repair. Some deletions contain DNA from other hotspots, indicating that double cutting at distant sites creates substrates for insertional mutagenesis. End joining at double cuts can also yield tandem duplications or extrachromosomal circles. Our findings highlight the importance of DSB regulation and reveal a previously hidden potential for meiotic mutagenesis that is likely to affect human health and genome evolution.


title: Transcriptome and metabolome analysis of crGART, a novel cell model of de novo purine synthesis deficiency: Alterations in CD36 expression and activity

abstract: In humans, GART [phosphoribosylglycinamide formyltransferase (EC 2.1.2.2) / phosphoribosylglycinamide synthetase (EC 6.3.4.13) / phosphoribosylaminoimidazole synthetase (EC 6.3.3.1)] is a trifunctional protein which catalyzes the second, third, and fifth reactions of the ten step de novo purine synthesis (DNPS) pathway. The second step of DNPS is conversion of phosphoribosylamine (5-PRA) to glycineamide ribonucleotide (GAR). 5-PRA is extremely unstable under physiological conditions and is unlikely to accumulate in the absence of GART activity. Recently, a HeLa cell line null mutant for GART was constructed via CRISPR-Cas9 mutagenesis. This cell line, crGART, is an important cellular model of DNPS inactivation that does not accumulate DNPS pathway intermediates. In the current study, we characterized the crGART versus HeLa transcriptomes in purine-supplemented and purine-depleted growth conditions. We observed multiple transcriptome changes and discuss pathways and ontologies particularly relevant to Alzheimer disease and Down syndrome. We also report initial analysis of the Cluster of Differentiation (CD36) gene, which is highly expressed in crGART versus HeLa.


title: Decoding the microstructural properties of white matter using realistic models

abstract:
Multi-echo gradient echo (ME-GRE) magnetic resonance signal evolution in white matter has a strong dependence on the orientation of myelinated axons with respect to the main static field. Although analytical solutions have been able to predict some of the white matter (WM) signal behaviour of the hollow cylinder model, it has been shown that realistic models of WM offer a better description of the signal behaviour observed.

In this work, we present a pipeline to (i) generate realistic 2D WM models with their microstructure based on real axon morphology with adjustable fiber volume fraction (FVF) and g-ratio. We (ii) simulate their interaction with the static magnetic field to be able to simulate their MR signal. For the first time, we (iii) demonstrate that realistic 2D WM models can be used to simulate a MR signal that provides a good approximation of the signal obtained from a real 3D WM model derived from electron microscopy. We then (iv) demonstrate in silico that 2D WM models can be used to predict microstructural parameters in a robust way if ME-GRE multi-orientation data is available and the main fiber orientation in each pixel is known using DTI. A deep learning network was trained and characterized in its ability to recover the desired microstructural parameters such as FVF, g-ratio, free and bound water transverse relaxation and magnetic susceptibility. Finally, the network was trained to recover these micro-structural parameters from an ex vivo dataset acquired in 9 orientations with respect to the magnetic field and 12 echo times. We demonstrate that this is an overdetermined problem and that as few as 3 orientations can already provide comparable results for some of the decoded metrics.

[Highlights] - A pipeline to generate realistic white models of arbitrary fiber volume fraction and g-ratio is presented; - We present a methodology to simulated the gradient echo signal from segmented 2D and 3D models of white matter, which takes into account the interaction of the static magnetic field with the anisotropic susceptibility of the myelin phospholipids; - Deep Learning Networks can be used to decode microstructural white matter parameters from the signal of multi-echo multi-orientation data;


title: Hypomyelination leads to alterations in inhibitory function and parvalbumin-positive neuron density in auditory cortex

abstract: For a long time, myelin was thought to be restricted to excitatory neurons, and studies on dysmyelination focused primarily on excitatory cells. Recent evidence showed that axons of inhibitory neurons in the neocortex are also myelinated, but the role of myelin on inhibitory circuits remains unknown. Here we studied the impact of mild hypomyelination on both excitatory and inhibitory connectivity in the primary auditory cortex (A1) with well-characterized mouse models of hypomyelination due to loss of oligodendrocyte ErbB receptor signaling. Using laser-scanning photostimulation, we found that mice with mild hypomyelination have reduced functional inhibitory connections to A1 L2/3 neurons without changes in excitatory connections, resulting in altered excitatory/inhibitory balance. These effects are not associated with altered expression of GABAergic and glutamatergic synaptic components, but with reduced density of parvalbumin-positive (PV+) neurons, which reflects reduced PV expression by interneurons rather than PV+ neuronal loss. While immunostaining shows that hypomyelination occurs in both PV+ and PV- axons, there is a strong correlation between MBP and PV expression suggesting that myelination influences PV expression. Together, the results demonstrate that mild hypomyelination impacts A1 neuronal networks, reducing inhibitory activity, and shifting networks towards excitation.


title: Hypomyelination reduces parvalbumin interneuron density and auditory cortex inhibitory function


title: A fast image processing toolbox for all-optical closed-loop control of neuronal activities

abstract: The development of in vivo imaging and optogenetic tools makes it possible to control neural circuit activities in an all-optical, closed-loop manner [1], but such applications are limited by the lack of software for online analysis of neuronal imaging data. We developed an analysis software ORCA (Online Real-time activity and offline Cross-session Analysis), which performs image registration, neuron segmentation, and activity extraction at over 100 frames per second, fast enough to support real-time detection and readout of neural activity. Our active neuron detection algorithm is purely statistical, achieving much higher speed than previous methods. We demonstrated closed-loop control of neurons that were identified on the fly, without prior recording or image processing. ORCA also includes a cross-session alignment module that efficiently tracks neurons across multiple sessions. In summary, ORCA is a powerful toolbox for fast imaging data analysis, and provides a solution for all-optical closed-loop control of neuronal activity.


title: Real-time image processing toolbox for all-optical closed-loop control of neuronal activities


title: Wandering minds, sleepy brains: lapses of attention and slow waves in wakefulness.

abstract: Attentional lapses are ubiquitous and can negatively impact performance. They correlate with mind wandering, or thoughts that are unrelated to ongoing tasks and environmental demands. In other cases, the stream of consciousness itself comes to a halt and the mind goes blank. What is happening in the brain that leads to these mental states? To understand the neural mechanisms underlying attentional lapses, we cross-analysed the behaviour, subjective experience and neural activity of healthy participants performing a task. Random interruptions prompted participants to indicate their mental states as task-focused, mind-wandering or mind-blanking. High-density electroencephalography revealed the occurrence of spatially and temporally localized slow waves, a pattern of neural activity characteristic of the transition toward sleep. These slow waves accompanied behavioural markers of lapses and preceded reports of mind wandering and mind blanking. Furthermore, the location of slow waves distinguished sluggish versus impulsive behaviours, mind wandering versus mind blanking. Our results suggest attentional lapses share a common physiological origin: the emergence of local sleep-like activity within the awake brain.


title: Predicting lapses of attention with sleep-likeslow waves.


title: Integration of genetic, transcriptomic, and clinical data provides insight into 16p11.2 and 22q11.2 CNV genes

abstract: Deletions and duplications of the multigenic 16p11.2 and 22q11.2 copy number variants (CNVs) are associated with brain-related disorders including schizophrenia, intellectual disability, obesity, bipolar disorder, and autism spectrum disorder (ASD). The contribution of individual CNV genes to each of these phenotypes is unknown, as is the contribution of CNV genes to subtler health impacts. Hypothesizing that DNA copy number acts via RNA expression, we attempted a novel in silico fine-mapping approach in non-carriers using both GWAS and biobank data. We first asked whether expression level of a CNV gene impacts risk for a known brain-related phenotype(s). Using transcriptomic imputation, we tested for association within GWAS for schizophrenia, IQ, BMI, bipolar disorder, and ASD. We found individual genes in 16p11.2 associated with schizophrenia, BMI, and IQ (SPN), using conditional analysis to identify INO80E as the driver of schizophrenia, and SPN and INO80E as drivers of BMI. Second, we used a biobank containing electronic health data to compare the medical phenome of CNV carriers to controls within 700,000 individuals to investigate a spectrum of health effects, identifying novel and previously observed traits. Third, we used genotypes for over 48,000 biobank individuals to perform phenome-wide association studies between imputed expressions of 16p11.2 and 22q11.2 genes and over 1,500 health traits, finding seventeen significant gene-trait pairs, including psychosis (NPIPB11, SLX1B) and mood disorders (SCARF2), and overall enrichment of mental traits. Our results demonstrate how integration of genetic and clinical data aids in understanding CNV gene function, and implicate pleiotropy and multigenicity in CNV biology.

category: genetics


title: Heterogeneity of murine periosteum progenitors involved in fracture healing

abstract: The periosteum is the major source of cells involved in fracture healing. We sought to characterize progenitor cells and their contribution to bone fracture healing. The periosteum is highly enriched for progenitor cells, including Sca1+ cells, CFU-F and label-retaining cells compared to the endosteum and bone marrow. Using lineage tracing, we demonstrate that SMA identifies long-term, slow-cycling, self-renewing osteochondroprogenitors in the adult periosteum that are functionally important for bone formation during fracture healing. In addition, Col2.3CreER-labeled osteoblast cells contribute around 10% of osteoblasts, but no chondrocytes in fracture calluses. Most periosteal osteochondroprogenitors following fracture, can be targeted by SMACreER. Previously identified skeletal stem cell populations were common in periosteum, but contained high proportions of mature osteoblasts. We have demonstrated that the periosteum is highly enriched for skeletal progenitor cells and there is heterogeneity in the populations of cells that contribute to mature lineages during periosteal fracture healing.


title: Only two out of five articles by New Zealand researchers are free-to-access: a multiple API study of access, citations, cost of Article Processing Charges (APC), and the potential to increase the proportion of open access

abstract: We studied journal articles published by researchers at all eight of New Zealand universities in 2017 to determine how many were freely accessible on the web. We wrote software code to harvest data from multiple sources, code that we now share to enable others to reproduce our work on their own sample set. In May 2019, we ran our code to determine which of the 2017 articles were open at that time and by what method; where those articles would have incurred an Article Processing Charge (APC) we calculated the cost if those charges had been paid. Where articles were not freely available we determined whether the policies of publishers in each case would have allowed deposit in a non-commercial repository (Green open access). We also examined citation rates for different types of access. We found that, of our 2017 sample set, about two out of every five articles were freely accessible without payment or subscription (41%). Where research was explicitly said to be funded by New Zealands major research funding agencies, the proportion was slightly higher at 45%. Where open articles would have incurred an APC we estimated an average cost per article of USD1,682 (for publications where all articles require an APC, that is, Gold open access) and USD2,558 (where APC payment is optional, Hybrid open access) at a total estimated cost of USD1.45m. Of the paid options, Gold is by far more common for New Zealand researchers (82% Gold, 18% Hybrid). In terms of citations, our analysis aligned with previous studies that suggest a correlation between publications being freely accessible and, on balance, slightly higher rates of citation. This is not seen across all types of open access, however, with Diamond OA achieving the lowest rates. Where articles were not freely accessible we found that a very large majority of them (88% or 3089 publications) could have been legally deposited in an institutional repository. Similarly, only in a very small number of cases had a version deposited in the repository of a New Zealand university made the difference between the publication being freely accessible or not (125 publications). Given that most New Zealand researchers support research being open, there is clearly a large gap between belief and practice in New Zealands research ecosystem.

category: scientific communication and education


title: Virus-derived peptide inhibitors of the herpes simplex virus type 1 nuclear egress complex

abstract: Herpesviruses infect a majority of the human population, establishing lifelong latent infections for which there is no cure. Periodic viral reactivation spreads infection to new hosts while causing various disease states particularly detrimental in the immunocompromised. Efficient viral replication, and ultimately the spread of infection, is dependent on the nuclear egress complex (NEC), a conserved viral heterodimer that helps translocate mature viral capsids from the nucleus to the cytoplasm where they mature into infectious virions. Here, we have identified peptides capable of inhibiting the membrane-budding activity of the NEC in vitro. Biophysical characterization of these peptides using circular dichroism suggests that secondary structure, rather than amino acid sequence, influences the extent of inhibition. Current therapeutics that target viral DNA replication machinery are rendered ineffective by drug resistance due to viral mutations. Our results establish a basis for the development of an alternative class of inhibitors against nuclear egress, an essential step in herpesvirus replication, potentially expanding the current repertoire of available therapeutics.


title: Olfactory-Trigeminal Integration in the Primary Olfactory Cortex

abstract:
Humans naturally integrate signals from the olfactory and intranasal trigeminal systems. A tight interplay has been demonstrated between these two systems, and yet the underlying neural circuitry that mediates olfactory-trigeminal integration remains poorly understood. Using functional magnetic resonance imaging (fMRI), combined with psychophysics, this study investigated the neural mechanisms underlying olfactory-trigeminal integration. Fifteen participants with normal olfactory function performed a localization task with air-puff stimuli, phenylethyl alcohol (PEA; rose odor), or a combination thereof while being scanned. The ability to localize PEA to either nostril was at chance. Yet, its presence significantly improved the localization accuracy of weak, but not strong, air-puffs, relative to air-puff localization without concomitant PEA when both stimuli were delivered concurrently to the same nostril, but not when different nostrils received the two stimuli. This enhancement in localization accuracy, exemplifying the principles of spatial coincidence and inverse effectiveness in multisensory integration, was associated with multisensory integrative activity in the primary olfactory (POC), orbitofrontal (OFC), superior temporal (STC), inferior parietal (IPC) and cingulate cortices, and in the cerebellum. Multisensory enhancement in most of these regions, except the OFC, correlated with behavioral multisensory enhancement, as did increases in connectivity between some of these regions. We interpret these findings as indicating that the POC is part of a distributed brain network mediating integration between the olfactory and trigeminal systems.

HIGHLIGHTSO_LIPsychophysical and neuroimaging study of olfactory-mechanosensory (OM) integration
C_LIO_LIBehavior, cortical activity and network connectivity show OM integration
C_LIO_LIOM integration obeys principles of inverse effectiveness and spatial coincidence
C_LIO_LIBehavioral and neural measures of OM integration are correlated
C_LI


title: Neural Flip-Flops II: Short-Term Memory and Electroencephalography

abstract:
By enabling many brain structures state changes, the explicit cascaded oscillators proposed here can generate the rhythmic neural activity found in EEGs. The function of such synchronization in information processing systems is timing error avoidance. The narrow requirement for the oscillator input pulse duration suggests a possible relationship to the abnormal electrical activity characteristic of epileptic seizures. Together, flip-flops and synchronization by oscillators suggest a resolution to the longstanding controversy of whether short-term memory depends on neurons firing persistently or in brief, coordinated bursts.

The proposed cascade of oscillators consists of a ring oscillator and four toggle flip-flops connected in sequence. The novel oscillator and toggle are composed of three and six neurons, respectively. Their operation depends only on minimal properties of excitatory and inhibitory inputs.

The hypothesis that cascaded oscillators produce EEG phenomena implies that the distribution of EEG frequencies is determined by just two parameters, the mean (d) and standard deviation ({sigma}d) of the delay times of neurons that make up the initial oscillators in the cascades. For example, if d and {sigma}d are measured in milliseconds, the boundary separating the alpha and beta frequency bands is O_FD O_INLINEFIG[Formula]C_INLINEFIGC_FD

With 4 and 1.5 ms being the best available estimates for d and {sigma}d, respectively, this predicted boundary value is 14.9 Hz, which is within the range of commonly cited estimates obtained empirically from EEGs. Four and 1.5 ms also accurately predict the peaks and other boundaries of the five major EEG frequency bands.

Significance statementThe neuronal model proposed here implies several major aspects of electroencephalography. The matched periods of neural activity found in EEGs, and their wide distribution across the brain and across the frequency spectrum, follow from selective pressure for a biologically useful function: timing error avoidance for diverse brain functions in the trade-off between speed and accuracy. This activity can be achieved with a simple organization of synaptic connections and minimal neuron capabilities of excitation and inhibition. The multimodal distribution of EEG frequencies is an explicit function of the mean and variance of neuron delay times. The model suggests a relationship to epileptic seizures and a resolution to a short-term memory controversy. Two EEG characteristics make other models implausible.


title: Neural Flip-Flops II: The Role of Cascaded Oscillators in Short-Term Memory, EEGs, and Epilepsy


title: CoronaHiT: High throughput sequencing of SARS-CoV-2 genomes

abstract: The COVID-19 pandemic has spread to almost every country in the world since it started in China in late 2019. Controlling the pandemic requires a multifaceted approach including whole genome sequencing to support public health interventions at local and national levels. One of the most widely used methods for sequencing is the ARTIC protocol, a tiling PCR approach followed by Oxford Nanopore sequencing (ONT) of up to 96 samples at a time. There is a need, however, for a flexible, platform agnostic, method that can provide multiple throughput options depending on changing requirements as the pandemic peaks and troughs. Here we present CoronaHiT, a method capable of multiplexing up to 96 small genomes on a single MinION flowcell or >384 genomes on Illumina NextSeq, using transposase mediated addition of adapters and PCR based addition of barcodes to ARTIC PCR products. We demonstrate the method by sequencing 95 and 59 SARS-CoV-2 genomes for routine and rapid outbreak response runs, respectively, on Nanopore and Illumina platforms and compare to the standard ARTIC LoCost nanopore method. Of the 154 samples sequenced using the three approaches, genomes with [≥] 90% coverage (GISAID criteria) were generated for 64.3% of samples for ARTIC LoCost, 71.4% for CoronaHiT-ONT, and 76.6% for CoronaHiT-Illumina and have almost identical clustering on a maximum likelihood tree. In conclusion, we demonstrate that CoronaHiT can multiplex up to 96 SARS-CoV-2 genomes per MinION flowcell and that Illumina sequencing can be performed on the same libraries, which will allow significantly higher throughput. CoronaHiT provides increased coverage for higher Ct samples, thereby increasing the number of high quality genomes that pass the GISAID QC threshold. This protocol will aid the rapid expansion of SARS-CoV-2 genome sequencing globally, to help control the pandemic.

category: genomics


title: Synthetic Heparan Sulfate Mimetic Pixatimod (PG545) Potently Inhibits SARS-CoV-2 By Disrupting The Spike-ACE2 interaction

abstract: Heparan sulfate (HS) is a cell surface polysaccharide recently identified as a co-receptor with the ACE2 protein for recognition of the S1 spike protein on SARS-CoV-2 virus, providing a tractable new target for therapeutic intervention. Clinically-used heparins demonstrate inhibitory activity, but world supplies are limited, necessitating alternative solutions. Synthetic HS mimetic pixatimod is a drug candidate for cancer with immunomodulatory and heparanase-inhibiting properties. Here we show that pixatimod binds to and destabilizes the SARS-CoV-2 spike protein receptor binding domain (S1-RBD), and directly inhibits its binding to human ACE2, consistent with molecular modelling identification of multiple molecular contacts and overlapping pixatimod and ACE2 binding sites. Assays with multiple clinical isolates of live SARS-CoV-2 virus show that pixatimod potently inhibits infection of monkey Vero E6 and human bronchial epithelial cells at concentrations within its safe therapeutic dose range. Furthermore, in a K18-hACE2 mouse model pixatimod demonstrates that pixatimod markedly attenuates SARS-CoV-2 viral titer and COVID-19-like symptoms. This demonstration of potent anti-SARS-CoV-2 activity establishes proof-of-concept for targeting the HS-Spike protein-ACE2 axis with synthetic HS mimetics. Together with other known activities of pixatimod our data provides a strong rationale for its clinical investigation as a potential multimodal therapeutic to address the COVID-19 pandemic.


title: Systems genetics analyses in Diversity Outbred mice inform human bone mineral density GWAS and identify Qsox1 as a novel determinant of bone strength

abstract: Genome-wide association studies (GWASs) for osteoporotic traits have identified over 1000 associations; however, their impact has been limited by the difficulties of causal gene identification and a strict focus on bone mineral density (BMD). Here, we used Diversity Outbred (DO) mice to directly address these limitations by performing the first systems genetics analysis of 55 complex skeletal phenotypes. We applied a network approach to cortical bone RNA-seq data to discover 72 genes likely to be causal for human BMD GWAS associations, including the novel genes SERTAD4 and GLT8D2. We also performed GWAS in the DO for a wide-range of bone traits and identified Qsox1 as a novel gene influencing cortical bone accrual and bone strength. Our results provide a new perspective on the genetics of osteoporosis and highlight the ability of the mouse to inform human genetics.


title: Wolbachia-conferred antiviral protection is determined by developmental temperature

abstract:
Wolbachia is a maternally transmitted bacterium widespread in arthropods and filarial nematodes, and confers strong antiviral protection in Drosophila melanogaster and other insects. Wolbachia-transinfected Aedes aegypti are currently being deployed to fight transmission of dengue and Zika viruses. However, the mechanism of antiviral protection and factors influencing it are still not fully understood. Here we show that temperature modulates Wolbachia-conferred protection in Drosophila melanogaster. Temperature after infection directly impacts Drosophila C virus replication and modulates Wolbachia protection. At higher temperatures virus proliferates more and is more lethal, while Wolbachia confers lower protection. Strikingly, host developmental temperature is a determinant of Wolbachia-conferred antiviral protection. While there is a strong protection when flies are raised from egg to adult at 25{degrees}C, the protection is highly reduced or completely abolished when flies develop at 18{degrees}C. However, Wolbachia-induced changes during development are not sufficient to limit virus-induced mortality, as Wolbachia is still required to be present in adults at the time of infection. This developmental effect is general, since it was present in different host genotypes, Wolbachia variants and upon infection with different viruses. Overall, we show that Wolbachia-conferred antiviral protection is temperature dependent, being present or absent depending on the environmental conditions. This interaction likely impacts Wolbachia-host interactions in nature and, as a result, frequencies of host and symbionts in different climates. Dependence of Wolbachia-mediated pathogen blocking on developmental temperature could be used to dissect the mechanistic bases of protection and should be considered by programmes deploying Wolbachia as an antiviral agent in the field.

Significance StatementInsects are often infected with beneficial intracellular bacteria. The bacterium Wolbachia can protect insects from pathogenic viruses. This effect can be used to prevent transmission of dengue and Zika viruses by Wolbachia-infected mosquitoes. To deploy Wolbachia in the field successfully and understand the biology of insects in the wild we need to discover which factors affect Wolbachia-conferred antiviral protection. Here we show that the temperature in which insects develop from eggs to adults can determine presence or absence of antiviral protection. The environment, therefore, influences this insect-bacterium interaction. Our work may help to provide insights into the mechanism of viral blocking by Wolbachia and inform programs using Wolbachia in mosquito-borne disease control.


title: JIP3 links lysosome transport to regulation of multiple components of the axonal cytoskeleton

abstract: Lysosome axonal transport is important for the clearance of cargoes sequestered by the endocytic and autophagic pathways. Building on observations that mutations in the JIP3 (MAPK8IP3) gene result in lysosome-filled axonal swellings, we analyzed the impact of JIP3 depletion on the cytoskeleton of human neurons. Dynamic focal lysosome accumulations were accompanied by disruption of the axonal periodic scaffold (spectrin, F-actin and myosin II) throughout each affected axon. Additionally, axonal microtubule organization was locally disrupted at each lysosome-filled swelling. This local axonal microtubule disorganization was accompanied by accumulations of both F-actin and myosin II. These results indicate that transport of axonal lysosomes is functionally interconnected with mechanisms that control the organization and maintenance of the axonal cytoskeleton. They have potential relevance to human neurological disease arising from JIP3 mutations as well as for neurodegenerative diseases associated with the focal accumulations of lysosomes within axonal swellings such as Alzheimers disease.


title: Oxytocin excites BNST interneurons and inhibits BNST output neurons to the central amygdala

abstract:
The dorsolateral bed nucleus of the stria terminalis (BNSTDL) has high expression of oxytocin (OT) receptors (OTR), which were shown to facilitate cued fear. However, the role of OTR in the modulation of BNSTDL activity remains elusive. BNSTDL contains GABA-ergic neurons classified based on intrinsic membrane properties into three types. Using in vitro patch-clamp recordings in male rats, we demonstrate that OT selectively excites and increases spontaneous firing rate of Type I BNSTDL neurons. As a consequence, OT increases the frequency, but not amplitude, of spontaneous inhibitory post-synaptic currents (sIPSCs) selectively in Type II neurons, an effect abolished by OTR antagonist or tetrodotoxin, and reduces spontaneous firing rate in these neurons. These results suggest an indirect effect of OT in Type II neurons, which is mediated via OT-induced increase in firing of Type I interneurons. As Type II BNSTDL neurons were shown projecting to the central amygdala (CeA), we also recorded from retrogradelly labeled BNST[->]CeA neurons and we show that OT increases the frequency of sIPSC in these Type II BNST[->]CeA output neurons. In contrast, in Type III neurons, OT reduces the amplitude, but not frequency, of both sIPSCs and evoked IPSCs via a postsynaptic mechanism without changing their intrinsic excitability. We present a model of fine-tuned modulation of BNSTDL activity by OT, which selectively excites BNSTDL interneurons and inhibits Type II BNST[->]CeA output neurons. These results suggest that OTR in the BNST might facilitate cued fear by inhibiting the BNST[->]CeA neurons.

HighlightsO_LIOxytocin directly excites and increases spontaneous firing of Type I BNST interneurons
C_LIO_LIOxytocin indirectly inhibits Type II BNST neurons
C_LIO_LIOxytocin inhibits Type II BNST output neurons to the central amygdala (BNST[->]CeA)
C_LIO_LIOxytocin reduces GABA-ergic transmission in Type III BNST neurons
C_LIO_LIOxytocin might facilitate cued fear by inhibiting the Type II BNST[->]CeA neurons
C_LI

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=174 HEIGHT=200 SRC="FIGDIR/small/169466v3_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@134326dorg.highwire.dtl.DTLVardef@1dad252org.highwire.dtl.DTLVardef@1263dd6org.highwire.dtl.DTLVardef@1a5c9e9_HPS_FORMAT_FIGEXP  M_FIG C_FIG

category: neuroscience


title: Increased Expression of Chondroitin Sulfotransferases following AngII may Contribute to Pathophysiology Underlying Covid-19 Respiratory Failure: Impact may be Exacerbated by Decline in Arylsulfatase B Activity

abstract: The precise mechanisms by which Covid-19 infection leads to hypoxia and respiratory failure have not yet been elucidated. Interactions between sulfated glycosaminoglycans (GAGs) and the SARS-CoV-2 spike glycoprotein have been identified as participating in viral adherence and infectivity. The spike glycoprotein binds to respiratory epithelium through the angiotensin converting enzyme 2 (ACE2) receptor, which endogenously interacts with Angiotensin (Ang) II to yield Angiotensin 1-7. In this report, we show that stimulation of human vascular smooth muscle cells by Ang II leads to increased mRNA expression of two chondroitin sulfotransferases (CHST11 and CHST15), which are required for synthesis of chondroitin 4-sulfate (C4S) and chondroitin 4,6-disulfate (CSE), respectively. Also, increased total sulfated GAGs, increased sulfotransferase activity, and increased expression of the proteoglycans biglycan, syndecan, perlecan, and versican followed treatment by Ang II. Candesartan, an Angiotensin II receptor blocker (Arb), largely, but incompletely, inhibited these increases, and the differences from baseline remained significant. These results suggest that another effect of Ang II also contributes to the increased expression of chondroitin sulfotransferases, total sulfated GAGs, and proteoglycans. We hypothesize that activation of ACE2 may contribute to these increases and suggest that the SARS-CoV-2 spike glycoprotein interaction with ACE2 may also increase chondroitin sulfotransferases, sulfated GAGs, and proteoglycans and thereby contribute to viral adherence to bronchioalveolar cells and to respiratory compromise in SARS-CoV-2 infection.


title: Simple cloning of large natural product biosynthetic gene clusters from Streptomyces by an engineered CRISPR/Cas12a system

abstract: Directly cloning of biosynthetic gene clusters (BGCs) from microbial genomes has been revolutionizing the natural product-based drug discovery. However, it is still very challenging to efficiently clone, for example, large (> 80kb) and GC-rich (> 70%), streptomycete originating BGCs. In this study, we developed a simple, fast yet efficient and low-cost in vitro platform for direct cloning large BGCs from streptomycete genomic DNA, named as CAT-FISHING (CRISPR/Cas12a- and Agarose plug-based sysTem for Fast bIoSyntHetIc geNe cluster cloninG), by combining the advantages of CRISPR/Cas12a cleavage and bacterial artificial chromosome (BAC) library construction. CAT-FISHING was demonstrated by directly cloning large DNA fragments ranging from 47 to 139 kb with GC content of > 70% from the S. albus J1074 genome in a relatively efficient manner. Moreover, surugamides, encoded by a captured 87-kb BGC with GC content of 76%, was heterologously expressed in a Streptomyces chassis. These results indicate that CAT-FISHING is a powerful platform for BGCs batch cloning, which would be greatly beneficial to the natural products-based drug discovery. We believe that this system will lead a renaissance of interest in microorganisms as a source for drug development.

category: synthetic biology


title: Human papillomavirus type 16 sub-lineages and integration in cancer

abstract: Our lab has been intrigued by the fact that viral genomes often take on the role of mobile elements to perpetuate their existence in a complex organisms genome. Multiple DNA viruses such as Epstein-Barr virus, hepatitis B virus, and human papillomavirus (HPV) can invade their host genome, as "genomic parasites". We have been investigating HPV type 16 (HPV16), which is a prominent human tumour virus. In our recent in vitro work using 3D organoids, a common variant of HPV16s coding region elicited early integration into the host genome compared to the HPV16 prototype sequence. Next-generation sequencing (NGS) data confirmed a transcriptomic profile of increased proliferation and chromosomal instability--both hallmarks of cancer. Epidemiologically, this variant is associated with a high cervical cancer incidence. To take inquiries a step further, we investigated variant-specific integration across HPV16-related cancers using NGS data from population-derived clinical samples in The Cancer Genome Atlas (TCGA)-curated database. Data were analyzed for HPV16 positivity, sub-lineage, and viral-host integration using a bioinformatic pipeline of open-source tools, including HPVDetector. Here, we report the analysis of 120 cervical cancer cases comprising HPV16 positive and negative samples as well as their different sub-lineage and integration status. The integration signature between variant and prototype did not differ quantitatively but qualitatively: that of the variant being related to hypoxia/energetics (Warburg effect) and that of the prototype being much more varied to include host immune abrogation and cancer pathways activation. We conclude by discussing challenges and future directions for expanding these analyses.


title: Learning Mitigates Genetic Drift

abstract: Genetic drift is a basic evolutionary principle describing random changes in allelic frequencies, with far-reaching consequences in various topics ranging from species conservation efforts to speciation. The conventional approach assumes that genetic drift has the same effect on all populations undergoing the same change in size, regardless of different behaviors and history of the populations. However, here we reason that processes leading to a systematic increase of individuals chances of survival, such as learning or immunological memory, can mitigate the effects of genetic drift even if the overall mortality rate in the population does not change. We further test this notion in an agent-based model monitoring allele frequencies in a population of prey, either able or not able to learn. Importantly, both these populations start with the same effective size and have the same and constant overall mortality rates. Our results demonstrate that even under these conditions, learning can mitigate drift if generations overlap. Furthermore, this effect holds regardless if the population is haploid or diploid or whether it reproduces sexually or asexually. Therefore, our findings demonstrate that learning is an overlooked factor affecting the effective population size. These findings may be of importance not only for basic evolutionary theory but also for other fields using the concept of genetic drift.


title: Computational Studies of Substrate Transport and Specificity in a Phospholipid Flippase

abstract: Type 4 P-type ATPases are lipid flippases which help maintain asymmetric phospholipid distribution in eukaryotic membranes by driving unidirectional translocation of phospholipid substrates. Recent cryo-EM and crystal structures have provided a detailed view of flippases, and we here use molecular dynamics simulations to study the first steps of phospholipid transport and lipid substrate specificity. Our simulations and new cryo-EM structure shows phospholipid binding to a groove and subsequent movement towards the centre of the membrane, and reveal a preference for phosphatidylserine lipids. We find that only the lipid head group stays in the groove while the lipid tails remain in the membrane, thus visualizing how flippases have evolved to transport large substrates. The flippase also induces deformation and thinning of the outer leaflet facilitating lipid recruitment. Our simulations provide insight into substrate binding to flippases and suggest that multiple sites and steps in the functional cycle contribute to substrate selectivity.


title: Novel PE and APC Tandems: Additional Near-Infrared Fluorochromes for Use in Spectral Flow Cytometry

abstract: Recent advances in flow cytometry instrumentation and fluorochrome chemistries have greatly increased fluorescent conjugated antibody combinations that can be used reliably and easily in routine experiments. The Cytek Aurora flow cytometer has three excitation lasers (405nm, 488nm, and 640nm) and incorporates the latest Avalanche Photodiode (APD) technology, demonstrating significant improvement in sensitivity for fluorescent emission signals longer than 800nm. However, there are limited commercially available fluorochromes capable of excitation with peak emission signals beyond 800nm. To address this gap, we engineered 6 new fluorochromes: PE-750, PE-800, PE-830 for the 488nm laser and APC-750, APC-800, APC-830 for the 640nm laser. Utilizing the principal of fluorescence resonance energy transfer (FRET), these novel structures were created by covalently linking a protein donor dye with an organic small molecule acceptor dye. Additionally, each of these fluorochrome conjugates were shown to be compatible with Fixation/Permeabilization buffer reagents, and demonstrated acceptable brightness and stability when conjugated to antigen-specific monoclonal antibodies. These six novel fluorochrome reagents can increase the numbers of fluorochromes that can be used on a 3-laser spectral flow cytometer.

category: immunology