title: Multisensory correlation computations in the human brain uncovered by a time-resolved encoding model

abstract: Neural mechanisms that arbitrate between integrating and segregating multisensory information are essential for complex scene analysis and for the resolution of the multisensory correspondence problem. However, these mechanisms and their dynamics remain largely unknown, partly because classical models of multisensory integration are static. Here, we used the Multisensory Correlation Detector, a model that provides a good explanatory power for human behavior while incorporating dynamic computations. Participants judged whether sequences of auditory and visual signals originated from the same source (causal inference) or whether one modality was leading the other (temporal order), while being recorded with magnetoencephalography. To test the match between the Multisensory Correlation Detector dynamics and the magnetoencephalographic recordings, we developed a novel dynamic encoding-model approach of electrophysiological activity, which relied on temporal response functions. First, we confirm that the Multisensory Correlation Detector explains causal inference and temporal order patterns well. Second, we found strong fits of brain activity to the two outputs of the Multisensory Correlation Detector in temporo-parietal cortices, a region with known multisensory integrative properties. Finally, we report an asymmetry in the goodness of the fits, which were more reliable during the causal inference than during the temporal order judgment task. Overall, our results suggest the plausible existence of multisensory correlation detectors in the human brain, which explain why and how causal inference is strongly driven by the temporal correlation of multisensory signals.


title: H3K4 Methylation in β-cells Prevents Transcriptional Downregulation and Variance Associated with Type 2 Diabetes

abstract: Histone 3 lysine 4 trimethylation (H3K4me3) is associated with promoters of actively expressed genes, with genes important for cell identity frequently having exceptionally broad H3K4me3-enriched domains at their TSS. While H3K4 methylation is implicated in contributing to transcription, maintaining transcriptional stability, facilitating enhancer-promoter interactions, and preventing irreversible silencing, some studies suggest it has little functional impact. Therefore, the function of H3K4 methylation is not resolved. Insufficient insulin release by {beta}-cells is the primary etiology in type 2 diabetes (T2D) and is associated with the loss of expression of genes essential to normal {beta}-cell function. We find that H3K4me3 is reduced in islets from mouse models of diabetes and from human donors with T2D. Using a genetic mouse model to impair H3K4 methyltransferase activity of TrxG complexes, we find that reduction of H3K4 methylation significantly reduces insulin production and glucose-responsiveness and increases transcriptional entropy, indicative of a loss of {beta}-cell maturity. Genes that are downregulated by reduction to H3K4 methylation are concordantly downregulated in T2D. Loss of H3K4 methylation causes global dilution of epigenetic complexity but does not generally reduce gene expression - instead, genes related to {beta}-cell function and/or in particular chromatin environments are specifically affected. While neither H3K4me3 nor H3K4me1 are strictly required for the expression of many genes, the expression of genes with critical roles in {beta}-cell function becomes destabilized, with increased variance and decreased overall expression. Our data further suggests that, in absence of H3K4me3, promoter-associated H3K4me1 is sufficient to maintain expression. Together, these data implicate H3K4 methylation dysregulation as destabilizing {beta}-cell gene expression and contributing to {beta}-cell dysfunction in T2D.

category: cell biology


title: Dynamic variation in a socially exchanged fluid indicates a role in the regulation of ant colony development and maturation

abstract: Socially exchanged fluids, like seminal fluid and milk, allow individuals to molecularly influence conspecifics. Many social insects have a social circulatory system, where food and endogenously produced molecules are transferred mouth-to-mouth (stomodeal trophallaxis), connecting all the individuals in the society. To understand how these endogenous molecules relate to colony life, we used quantitative proteomics to investigate the trophallactic fluid within colonies of the carpenter ant Camponotus floridanus. We show that different stages of the colony life cycle circulate different types of proteins: young colonies prioritize direct carbohydrate processing; mature colonies prioritize accumulation and transmission of stored resources. Further, colonies circulate proteins implicated in oxidative stress, ageing, and social insect caste determination, potentially acting as superorganismal hormones. Brood-caring individuals that are also closer to the queen in the social network (nurses) showed higher abundance of oxidative stress-related proteins. Thus, trophallaxis behavior provides a mechanism for distributed metabolism in social insect societies.

category: evolutionary biology


title: Collective regulation of metabolism, development and longevity through a socially exchanged fluid


title: Coevolution between primates and venomous snakes revealed by α-neurotoxin susceptibility

abstract:
Evidence suggests venomous snakes and primates have evolved certain traits in response to a coevolutionary arms-race. In both clades, evolved traits include an increase in brain size and enhanced vision. Lineage specific traits include in primates an inherent fear of snakes, while cobras have evolved defensive toxins, hooding, aposematism and venom spitting. To strengthen the claims of coevolution between venomous snakes and primates, more evidence of coevolved traits is needed to highlight the importance of this arms-race. We report a significantly reduced susceptibility of snake venom -neurotoxins toward the -1 nicotinic acetylcholine receptor orthosteric site within the catarrhine primates. This trait is particularly amplified within the clade Homininae. This relationship is supported by post-synaptic neurotoxic symptoms of envenoming relative to prey species being much lower humans due to weak binding of -neurotoxins to human nicotinic receptors. Catarrhines are sympatric with many species of large, diurnal, neurotoxically venomous snakes and as such are likely to have had a long history of interaction with them. Conversely, the Lemuriformes and Platyrrhini are highly susceptible to binding of -neurotoxins, which is consistent with them occupying geographical locations either devoid of venomous snakes or areas with neurotoxic snakes that are small, fossorial, and nocturnal. These data are consistent with the snake detection theory in that they follow a similar pattern of evolved traits within specific primate clades that are sympatric with venomous snakes. These results add new strong evidence in support of snakes and primates coevolving through arms-races that shaped selection pressures for both lineages.

Significance StatementWe have discovered a pattern of primate susceptibility towards -neurotoxins that supports the theory of a coevolutionary arms-race between venomous snakes and primates.


title: Sex-dependent gene regulation of human atherosclerotic plaques by DNA methylation and transcriptome integration points to smooth muscle cell involvement in women.

abstract: Sex differences are evident in the clinical presentation and underlying histology of atherosclerotic disease with women developing more stable atherosclerotic lesions than men. It is unknown whether this is explained by sex differences in gene regulation in cellular compartments of atherosclerotic plaques. To study sex differences in gene regulation we performed genome-wide DNA methylation and transcriptomics analysis on plaques of 485 carotid endarterectomy patients (31% female). Sex-differential DNA methylation at 4,848 sites in the autosome was enriched for cell-fate commitment and developmental processes, and its deconvolution predicted more smooth muscle cells in females, as compared to more immune cells in males. RNA-sequencing of the same plaques corroborated the sex differences in DNA methylation predicted cell-types, in which genes that were higher expressed in females were enriched for TGF-beta signaling and extracellular matrix biology. In addition, female-biased genes were enriched for targeting by regulatory loci based on sex differential methylation. Lastly, by using single-cell RNA sequencing we showed that these female-biased genes are mostly expressed in smooth muscle cells, and higher expressed in smooth muscle cells from female (predominantly stable) plaques as compared to male (relatively unstable) plaques. Our approach identified female-biased genes in smooth muscle cells in fibrous atherosclerotic plaques. This points towards new mechanisms in smooth muscle cell biology of stable atherosclerotic plaques and offers new directions for research to develop new sex-specific therapeutics for atherosclerotic disease.

category: genomics


title: Structure of Ycf1p reveals the transmembrane domain TMD0 and the regulatory R region of ABCC transporters

abstract:
ATP binding cassette (ABC) proteins typically function in active transport of solutes across membranes. The ABC core structure is comprised of two transmembrane domains (TMD1 and TMD2) and two cytosolic nucleotide binding domains (NBD1 and NBD2). Some members of the C-subfamily of ABC (ABCC) proteins, including human multidrug resistance proteins (MRPs), also possess an N-terminal transmembrane domain (TMD0) that contains five transmembrane -helices and is connected to the ABC core by the L0 linker. While TMD0 was resolved in SUR1, the atypical ABCC protein that is part of the hetero-octameric ATP-sensitive K+ channel, little is known about the structure of TMD0 in monomeric ABC transporters. Here, we present the structure of yeast cadmium factor 1 protein (Ycf1p), a homologue of human MRP1, determined by electron cryomicroscopy (cryo-EM). Comparison of Ycf1p, SUR1, and a structure of MRP1 that showed TMD0 at low resolution demonstrates that TMD0 can adopt different orientations relative to the ABC core, including a 145{degrees} rotation between Ycf1p and SUR1. The cryo-EM map also reveals that segments of the regulatory (R) region, which links NBD1 to TMD2 and was poorly resolved in earlier ABCC structures, interacts with the L0 linker, NBD1, and TMD2. These interactions, combined with fluorescence quenching experiments of isolated NBD1 with and without the R region, suggests how post-translational modifications of the R region modulate ABC protein activity. Mapping known mutations from MRP2 and MRP6 onto the Ycf1p structure explains how mutations involving TMD0 and the R region of these proteins lead to disease.

Statement of SignificanceThe Ycf1p structure provides an atomic model for the TMD0 domain of ABCC transporters and for two segments of the regulatory (R) region that links NBD1 to TMD2. The orientation of TMD0 in Ycf1p differs from that seen in SUR1, the regulatory ABCC protein in KATP channels, demonstrating flexibility in TMD0/ABC core contacts. The structure suggests how post-translational modifications of the R region modulate ABC protein activity and provides a mechanistic understanding of several diseases that occur due to mutation of human homologues of Ycf1p.

category: biochemistry


title: Selectivity of Protein Interactions along the Aggregation Pathway of a-Synuclein

abstract:
The aggregation of -SYN follows a cascade of oligomeric, prefibrillar and fibrillar forms, culminating in the formation of Lewy Bodies (LB), the pathological hallmarks of Parkinsons Disease in neurons. Whilst -synuclein is a major contributor to LB, these dense accumulations of protein aggregates and tangles of fibrils contain over 70 different proteins. However, the potential for interactions between these proteins and the different aggregated species of -SYN is largely unknown. We hypothesized that the proteins present in the Lewy Bodies are trapped or pulled into the aggregates in a hierarchical manner, by binding at specific stages of the aggregation of -SYN.

In this study we uncover a map of interactions of a total of 65 proteins, against different species formed by -SYN. We measured binding to monomeric -SYN using AlphaScreen, a sensitive nano-bead assay for detection of protein-protein interactions. To access different oligomeric species, we made use of the pathological mutants of -SYN (A30P, G51D and A53T), which form oligomeric species with distinct properties. Finally, we used bacterially expressed recombinant -SYN to generate amyloid fibrils and measure interactions with a pool of GFP-tagged potential partners. Binding to oligomers and fibrils was measured by two-color coincidence detection (TCCD) on a single molecule spectroscopy setup. Overall, we demonstrate that LB components are selectively recruited to specific steps in the formation of the LB, explaining their presence in the inclusions. Only a few proteins were found to interact with -SYN monomers at detectable levels, and only a subset recognizes the oligomeric -SYN including autophagosomal proteins. We therefore propose a new model for the formation of Lewy Bodies, where selectivity of protein partners at different steps drives the arrangement of these structures, uncovering new ways to modulate aggregation.

Significance StatementThe molecular complexity of the Lewy Bodies has been a major hindrance to a bottom-up reconstruction of these inclusions, protein by protein. This work presents an extensive dataset of protein-protein interactions, showing that despite its small size and absence of structure, -SYN binds to specific partners in the LB, and that there is a clear selectivity of interactions between the different -SYN species along the self-assembly pathway. We use single-molecule methods to deconvolute number and size of the co-aggregates, to gain detailed information about the mechanisms of interaction. These observations constitute the basis for the elaboration of a global interactome of -SYN.


title: Reaching to inhibit a prepotent response: a wearable 3-axis accelerometer kinematic analysis

abstract: The present work explores the distinctive contribution of motor planning and control to human reaching movements. In particular, the movements were triggered by the selection of a prepotent response (Dominant) or, instead, by the inhibition of the prepotent response, that required the selection of an alternative one (Non-dominant). To this aim, we adapted a Go/No-Go task to investigate both the dominant and non-dominant movements of a cohort of 19 adults, utilizing kinematic measures to discriminate between the planning and control components of the two actions. To sample such measures, a low-cost, easy to use, 3-axis wrist worn accelerometer was put to good use to obtain raw acceleration data and to compute and break down its velocity components. The values obtained with such task indicate that with the inhibition of a prepotent response, the selection and execution of the alternative one yields both a longer reaction time and movement duration. Moreover, the peak velocity occurred later in time with respect to the dominant response, revealing that participants tended to indulge more in motor planning rather than in adjusting their movement along the way. Finally, comparing such results to the findings obtained by other means in literature, we discuss the feasibility of an accelerometer-based analysis to disentangle distinctive cognitive mechanisms of human movements.

category: bioinformatics


title: Allosteric Regulation of the EphA2 Receptor Intracellular Region by Serine/Threonine Kinases

abstract: Eph receptor tyrosine kinases play a key role in cell-cell communication. However, lack of structural information on the entire multi-domain intracellular region of any Eph receptor has hindered detailed understanding of their signaling mechanisms. Here, we use an integrative structural biology approach combining X-ray crystallography, small-angle X-ray scattering and hydrogen-deuterium exchange mass spectrometry, to gain the first insights into the structure and dynamics of the entire EphA2 intracellular region. EphA2 promotes cancer malignancy through a poorly understood non-canonical form of signaling that depends on serine/threonine phosphorylation of the linker connecting the EphA2 kinase and SAM domains. We uncovered two distinct molecular mechanisms that may function in concert to mediate the effects of linker phosphorylation through an orchestrated allosteric regulatory network. The first involves a shift in the equilibrium between a "closed" configuration of the EphA2 intracellular region and an "open" more extended configuration induced by the accumulation of phosphorylation sites in the linker. This implies that cooperation of multiple serine/threonine kinase signaling networks is necessary to promote robust EphA2 non-canonical signaling. The second involves allosteric rearrangements in the kinase domain and juxtamembrane segment induced by phosphorylation of some linker residues, suggesting a link between EphA2 non-canonical signaling and canonical signaling through tyrosine phosphorylation. Given the key role of EphA2 in cancer malignancy, this new knowledge can inform therapeutic strategies.

category: biophysics


title: Generating Agent-Based Multiscale Multicellular Spatiotemporal Models from Ordinary Differential Equations of Biological Systems, with Applications in Viral Infection

abstract:
The biophysics of an organism span scales from subcellular to organismal and include spatial processes like diffusion of molecules, cell migration, and flow of intravenous fluids. Mathematical biology seeks to explain biophysical processes in mathematical terms at, and across, all relevant spatial and temporal scales. While non-spatial, ordinary differential equation (ODE) models are often used and readily calibrated to experimental data, they do not explicitly represent spatial and stochastic features of a biological system, limiting their insights and applications. Spatial models describe biological systems with spatial information but are mathematically complex and computationally expensive, which limits the ability to calibrate and deploy them. In this work we develop a formal method for deriving cell-based, spatial, multicellular models from ODE models of population dynamics in biological systems, and vice-versa. We provide examples of generating spatiotemporal, multicellular models from ODE models of viral infection and immune response. In these models the determinants of agreement of spatial and non-spatial models are the degree of spatial heterogeneity in viral production and rates of extracellular viral diffusion and decay. We show how ODE model parameters can implicitly represent spatial parameters, and cell-based spatial models can generate uncertain predictions through sensitivity to stochastic cellular events, which is not a feature of ODE models. Using our method, we can test ODE models in a multicellular, spatial context and translate information to and from non-spatial and spatial models, which help to employ spatiotemporal multicellular models using calibrated ODE model parameters, investigate objects and processes implicitly represented by ODE model terms and parameters, and improve the reproducibility of spatial, stochastic models. We hope to employ our method to generate new ODE model terms from spatiotemporal, multicellular models, recast popular ODE models on a cellular basis, and generate better models for critical applications where spatial and stochastic features affect outcomes.

Statement of SignificanceOrdinary differential equations (ODEs) are widely used to model and efficiently simulate multicellular systems without explicit spatial information, while spatial models permit explicit spatiotemporal modeling but are mathematically complicated and computationally expensive. In this work we develop a method to generate stochastic, agent-based, multiscale models of multicellular systems with spatial resolution at the cellular level according to non-spatial ODE models. We demonstrate how to directly translate model terms and parameters between ODE and spatial models and apply non-spatial model terms to boundary conditions using examples of viral infection modeling, and show how spatial models can interrogate implicitly represented biophysical mechanisms in non-spatial models. We discuss strategies for co-developing spatial and non-spatial models and reconciling disagreements between them.


title: Enhanced diffusion by reversible binding to active polymers

abstract:
Cells are known to use reversible binding to active biopolymer networks to allow diffusive transport of particles in an otherwise impenetrable mesh. We here determine the motion of a particle that experiences random forces during binding and unbinding events while being constrained by attached polymers. Using Monte-Carlo simulations and a statistical mechanics model, we find that enhanced diffusion is possible with active polymers. However, this is possible only under optimum conditions that has to do with the relative length of the chains to that of the plate. For example, in systems where the plate is shorter than the chains, diffusion is maximum when many chains have the potential to bind but few remain bound at any one time. Interestingly, if the chains are shorter than the plate, we find that diffusion is maximized when more active chains remain transiently bound. The model provides insight into these findings by elucidating the mechanisms for binding-mediated diffusion in biology and design rules for macromolecular transport in transient synthetic polymers.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=92 SRC="FIGDIR/small/428690v1_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@1f09992org.highwire.dtl.DTLVardef@f26b4borg.highwire.dtl.DTLVardef@162368corg.highwire.dtl.DTLVardef@eb2b45_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Exploring Energy Landscapes of Intrinsically Disordered Proteins: Insights into Functional Mechanisms

abstract: Intrinsically disordered proteins (IDPs) lack a rigid 3D structure and populate a polymorphic ensemble of conformations. Because of the lack of a reference conformation, their energy landscape representation in terms of reaction coordinates presents a daunting challenge. Here, our newly developed Energy Landscape Visualization Method (ELViM), a reaction coordinate-free approach, shows its prime application to explore frustrated energy landscapes of an intrinsically disordered protein, Prostate-Associated Gene 4 (PAGE4). PAGE4 is a transcriptional coactivator that potentiates the oncogene c-Jun. Two kinases, namely HIPK1 and CLK2, phosphorylate PAGE4 generating variants phosphorylated at different serine/threonine residues (HIPK1-PAGE4 and CLK2-PAGE4, respectively) with opposing functions. While HIPK1-PAGE4 predominantly phosphorylates Thr51 and potentiates c-Jun, CLK2-PAGE4 hyper-phosphorylates PAGE4 and attenuates transactivation. To understand the underlying mechanisms of conformational diversity among different phosphoforms, we have analyzed their atomistic trajectories simulated using AWSEM forcefield and the energy landscapes were elucidated using ELViM. This method allows us to identify and compare the population distributions of different conformational ensembles of PAGE4 phosphoforms using the same effective phase space. The results reveal a predominant conformational ensemble with an extended C-terminal segment of WT PAGE4, which exposes a functional residue Thr51, implying its potential of undertaking a fly-casting mechanism while binding to its cognate partner. In contrast, for HIPK1-PAGE4, a compact conformational ensemble enhances its population sequestering phosphorylated-Thr51. This clearly explains the experimentally observed weaker affinity of HIPK1-PAGE4 for c-Jun. ELViM appears as a powerful tool especially to analyze the highly-frustrated energy landscape representation of IDPs where appropriate reaction coordinates are hard to apprehend.


title: Biomembranes undergo complex, non-axisymmetric deformations governed by Kirchhoff-Love kinematics and revealed by a three dimensional computational framework

abstract: Biomembranes play a central role in various phenomena like locomotion of cells, cell-cell interactions, packaging and transport of nutrients, transmission of nerve impulses, and in maintaining organelle morphology and functionality. During these processes, the membranes undergo significant morphological changes through deformation, scission, and fusion. Modeling the underlying mechanics of such morphological changes has traditionally relied on reduced order axisymmetric representations of membrane geometry and deformation. Axisymmetric representations, while robust and extensively deployed, suffer from their inability to model symmetry breaking deformations and structural bifurcations. To address this limitation, a three-dimensional computational mechanics framework for high fidelity modeling of biomembrane deformation is presented. The proposed framework brings together Kirchhoff-Love thin-shell kinematics, Helfrich-energy based mechanics, and state-of-the-art numerical techniques for modeling deformation of surface geometries. Lipid bilayers are represented as spline-based surface discretizations immersed in a three-dimensional space; this enables modeling of a wide spectrum of membrane geometries, boundary conditions, and deformations that are physically admissible in a 3D space. The mathematical basis of the framework and its numerical machinery are presented, and their utility is demonstrated by modeling three classical, yet non-trivial, membrane deformation problems: formation of tubular shapes and their lateral constriction, Piezo1-induced membrane footprint generation and gating response, and the budding of membranes by protein coats during endocytosis. For each problem, the full three dimensional membrane deformation is captured, potential symmetry-breaking deformation paths identified, and various case studies of boundary and load conditions are presented. Using the endocytic vesicle budding as a case study, we also present a "phase diagram" for its symmetric and broken-symmetry states.


title: Quantification of in-plane flexoelectricity in lipid bilayers

abstract: Lipid bilayers behave as 2D dielectric materials that undergo polarization and deformation in the presence of an electric field. This effect has been previously modeled by continuum theories which assume a polarization field oriented normal to the membrane surface. However, the molecular architecture of the lipids reveals that the heqadgroup dipoles are primarily oriented tangential to the membrane surface. Here, we perform atomistic and coarse-grained molecular dynamics simulations to quantify the in-plane polarization undergone by a flat bilayer and a spherical vesicle in the presence of an applied electric field. We use these predictions to compute an effective in-plane flexoelectric coefficient for four different lipid types. Our findings provide the first molecular proof of the in-plane polarization undergone by lipid bilayers and furnish the material parameter required to quantify membrane-electric field interactions.


title: Glucose inhibits neutrophil migration and haemostasis, and accelerates diet-induced hyperlipidaemia in zebrafish embryos

abstract: Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals.

category: pathology


title: Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae


title: Macaque V1 responses to 2nd-order contrast-modulated stimuli and the possible subcortical and cortical contributions

abstract: V1 neurons as linear filters supposedly only respond to 1st-order luminance-modulated (LM) stimuli, but not 2nd-order contrast-modulated (CM) ones. To solve this difficulty, filter-rectify-filter models are proposed, in which first-stage filters respond to CM stimulus elements, and the nonlinear-rectified outputs are summed by a second-stage filter for CM stimulus representation. Correspondingly, neurophysiological evidence shows V1/A17 neurons less responsive to CM stimuli than V2/A18 neurons. Here we used two-photon calcium imaging to demonstrate substantial V1 responses to CM gratings with unimodally distributed LM/CM preferences. Moreover, LM responses are suppressed by LM and CM adaptations regardless of orientation, but CM responses are more suppressed by same-orientation LM and CM adaptations than by orthogonal ones. While LM adaptation results agree with the Hubel-Wiesel view of LGN contributions to V1 orientation responses, CM adaptation results, which include both orientation-unspecific and specific components, may suggest similar subcortical contributions plus additional refinement by recurrent intracortical interactions.

category: neuroscience


title: Pervasive false brain connectivity from electrophysiological signals

abstract: 1.Signals of brain electric neuronal activity, either invasively measured or non-invasively estimated, are commonly used for connectivity inference. One popular methodology assumes that the neural dynamics follow a multivariate autoregression, where the autoregressive coefficients represent the couplings among regions. If observation noise is present and ignored, as is common in practice, the estimated couplings are biased, affecting all forms of Granger-causality inference, both in time and in frequency domains. Significant nonsense coupling, i.e., nonsense connectivity, can appear when in reality there is none, since there is always observation noise in two possible forms: measurement noise, and activity from other brain regions due to volume conduction and low spatial resolution. This problem is critical, and is currently not being addressed, calling into question the validity of many Granger-causality reports in the literature. An estimation method that accounts for noise is based on an overdetermined system of high-order multivariate Yule-Walker equations, which give reduced variance estimators for the coupling coefficients of the unobserved signals. Simulation-based comparisons to other published methods are presented, demonstrating its adequate performance. In addition, simulation results are presented for a zero connectivity case with noisy observations, where the new method correctly reports no connectivity while classical analyses (as found in most software packages) report nonsense connectivity. For the sake of reproducible research, the supplementary material includes, in human readable format, all the time series data used here.


title: Impaired refinement of kinematic variability in Huntington disease mice on an automated home-cage forelimb motor task

abstract:
The effective development of novel therapies in mouse models of neurological disorders relies on behavioural assessments that provide accurate read-outs of neuronal dysfunction and/or degeneration. We designed an automated behavioural testing system ( PiPaw) which integrates an operant lever-pulling task directly into the mouse home-cage. This task is accessible to group-housed mice 24-hours per day, enabling high-throughput longitudinal analysis of forelimb motor learning. Moreover, this design eliminates the need for exposure to novel environments and minimizes experimenter interaction, significantly reducing two of the largest stressors associated with animal behaviour. Mice improved their performance of this task over one week of testing by reducing inter-trial variability of reward-related kinematic parameters (pull amplitude or peak velocity). In addition, mice displayed short-term improvements in reward rate, and a concomitant decrease in movement variability, over the course of brief (<10 minutes) bouts of task engagement. We used this system to assess motor learning in mouse models of the inherited neurodegenerative disorder, Huntington disease (HD). Despite having no baseline differences in task performance, Q175-FDN HD mice were unable to modulate the variability of their movements in order to increase reward on either short or long timescales. Task training was associated with a decrease in the amplitude of spontaneous excitatory activity recorded from striatal medium spiny neurons in the hemisphere contralateral to the trained forelimb in wildtype mice; however, no such changes were observed in Q175-FDN mice. This behavioural screening platform should prove useful for preclinical drug trials towards improved treatments in HD and other neurological disorders.

Significance StatementIn order to develop effective therapies for neurological disorders such as Huntington disease (HD), its important to be able to accurately and reliably assess the behaviour of mouse models of these conditions. Moreover, these behavioural assessments should provide an accurate readout of underlying neuronal dysfunction and/or degeneration. In this paper, we employed an automated behavioural testing system to assess motor learning in mice within their home-cage. Using this system, we were able to study motor abnormalities in HD mice with an unprecedented level of detail, and identified a specific behavioural deficit associated with an underlying impairment in striatal neuronal plasticity. These results validate the usefulness of this system for assessing behaviour in mouse models of HD and other neurological disorders.


title: Validation of the Munich Actimetry Sleep Detection Algorithm for estimating sleep-wake patterns from activity recordings

abstract:
Periods of sleep and wakefulness can be estimated from wrist-locomotor activity recordings via algorithms that identify periods of relative activity and inactivity. Here, we evaluated the performance of our Munich Actimetry Sleep Detection Algorithm (MASDA). MASDA uses a moving 24-hour-threshold and correlation procedure estimating relatively consolidated periods of sleep and wake.

MASDA was validated against sleep logs and polysomnography. Sleep-log validation was performed on 2 field samples collected over 54 and 34 days (median) in 34 adolescents and 28 young adults. Polysomnographic validation was performed on a clinical sample of 23 individuals undergoing 1 night of polysomnography. Epoch-by-epoch analyses were conducted and comparisons of sleep measures via Bland-Altman plots and correlations.

Compared with sleep logs, MASDA classified sleep with a median sensitivity of 80% (IQR = 75-86%) and specificity of 91% (87-92%). Mean onset and offset times were highly correlated (r = 0.86-0.91). Compared with polysomnography, MASDA reached a median sensitivity of 92% (85-100%), but low specificity of 33% (10-98%), owing to the low frequency of wake episodes in the nighttime polysomnographic recordings. MASDA overestimated sleep onset (~21 min) and underestimated wake after sleep onset (~26 min), while not performing systematically different from polysomnography in other sleep parameters.

These results demonstrate the validity of MASDA to faithfully estimate sleep-wake patterns in field studies. With its good performance across day- and nighttime, it enables analyses of sleep-wake patterns in long recordings performed to assess circadian and sleep regularity and is therefore an excellent objective alternative to sleep logs in field settings.


title: Lasting effects of a single psilocybin dose on resting-state functional connectivity in healthy individuals

abstract:
BackgroundPsilocybin is a psychedelic drug that has shown lasting positive effects on clinical symptoms and self-reported well-being following a single dose. There has been little research into the long-term effects of psilocybin on brain connectivity in humans.

AimsEvaluate changes in resting-state functional connectivity (RSFC) at one-week and three-months after one psilocybin dose in 10 healthy psychedelic-naive volunteers and explore associations between change in RSFC and related measures.

MethodsParticipants received 0.2-0.3 mg/kg psilocybin in a controlled setting. Participants completed resting-state fMRI scans at baseline, one-week and three-months post-administration and [11C]Cimbi-36 PET scans at baseline and one-week. We examined changes in within-network, between-network and region-to-region RSFC. We explored associations between changes in RSFC and psilocybin-induced phenomenology as well as changes in psychological measures and neocortex serotonin 2A receptor binding.

ResultsPsilocybin was well tolerated and produced positive changes in well-being. At one-week only, executive control network (ECN) RSFC was significantly decreased (Cohens d=-1.73, pFWE=0.010). We observed no other significant changes in RSFC at one-week or three-months, nor changes in region-to-region RSFC. Exploratory analyses indicated that decreased ECN RSFC at one-week predicted increased mindfulness at three-months (r =-0.65).

ConclusionsThese findings in a small cohort indicate that psilocybin affects ECN function within the psychedelic "afterglow" period. Our findings implicate ECN modulation as mediating psilocybin-induced, long-lasting increases in mindfulness. Although our findings implicate a neural pathway mediating lasting psilocybin effects, it is notable that changes in neuroimaging measures at three-months, when personality changes are observed, remain to be identified.


title: The SARS-CoV-2 Y453F mink variant displays a striking increase in ACE-2 affinity but does not challenge antibody neutralization

abstract: SARS-CoV-2 transmission from humans to animals has been reported for many domesticated species, including cats, dogs and minks. Identification of novel spike gene mutations appearing in minks has raised major concerns about potential immune evasion and challenges for the global vaccine strategy. The genetic variant, known as "cluster-five", arose among farmed minks in Denmark and resulted in a complete shutdown of the worlds largest mink production. However, the functional properties of this new variant are not established. Here we present functional data on the Y453F cluster-five receptor-binding domain (RBD) and show that it does not decrease established humoral immunity or affect the neutralizing response in a vaccine model based on wild-type RBD or spike. However, it binds the human ACE-2 receptor with a four-fold higher affinity suggesting an enhanced transmission capacity and a possible challenge for viral control.

category: immunology


title: Heterologous vaccination regimens with self-amplifying RNA and Adenoviral COVID vaccines induce superior immune responses than single dose vaccine regimens in mice

abstract: Several vaccines have demonstrated efficacy against SARS-CoV-2 mediated disease, yet there is limited data on the immune response induced by heterologous vaccination regimens using alternate vaccine modalities. Here, we present a detailed description of the immune response, in mice, following vaccination with a self-amplifying RNA (saRNA) vaccine and an adenoviral vectored vaccine (ChAdOx1 nCoV-19/AZD1222) against SARS-CoV-2. We demonstrate that antibody responses are higher in two dose heterologous vaccination regimens than single dose regimens. Neutralising titres after heterologous prime-boost were at least comparable or higher than the titres measured after homologous prime boost vaccination with viral vectors. Importantly, the cellular immune response after a heterologous regimen is dominated by cytotoxic T cells and Th1+ CD4 T cells which is superior to the response induced in homologous vaccination regimens in mice. These results underpin the need for clinical trials to investigate the immunogenicity of heterologous regimens with alternate vaccine technologies.


title: Heterologous vaccination regimens with self-amplifying RNA and Adenoviral COVID vaccines induce robust immune responses in mice


title: RNAi and Ino80 complex control rate limiting translocation step that moves rDNA to eroding telomeres

abstract: The discovery of HAATIrDNA, a mode of telomerase-negative survival in which canonical telomeres are replaced with ribosomal DNA (rDNA) repeats that acquire chromosome end-protection capability, raised crucial questions as to how rDNA tracts  jump to eroding, nonhomologous chromosome ends. Here we show that HAATIrDNA formation is initiated and limited by a single translocation that juxtaposes rDNA from Chromosome (Chr) III onto subtelomeric elements (STE) on Chr I or II; this rare reaction requires the RNAi pathway and the Ino80 nucleosome remodeling complex (Ino80C), thus defining an unforeseen relationship between these two machineries. The unique STE-rDNA junction created by this initial translocation is efficiently copied to the remaining STE chromosome ends, without the need for RNAi or Ino80C, forming HAATIrDNA. Intriguingly, both the RNAi and Ino80C machineries contain a component that plays dual roles in HAATI subtype choice. Dcr1 of the RNAi pathway and Iec1 of the Ino80C both promote HAATIrDNA formation as part of their respective canonical machineries, but both also inhibit formation of the exceedingly rare HAATISTE (in which STE sequences mobilize throughout the genome and assume chromosome end protection capacity) in non-canonical, pathway-independent manners. This work provides a glimpse into a previously unrecognized crosstalk between RNAi and Ino80C in controlling unusual translocation reactions that establish telomere-free linear chromosome ends.

category: molecular biology


title: Does insect herbivory suppress ecosystem productivity? Evidence from a temperate woodland

abstract: Our current understanding of the relationship between insect herbivory and ecosystem productivity is limited. Previous studies have typically quantified only leaf area loss, or have been conducted during outbreak years. These set-ups often ignore the physiological changes taking place in the remaining plant tissue after insect attack, or may not represent typical, non-outbreak herbivore densities. Here, we estimate the amount of carbon lost to insect herbivory in a temperate deciduous woodland both through leaf area loss and, notably, through changes in leaf gas exchange in non-consumed leaves under non-outbreak densities of insects. We calculate how net primary productivity changes with decreasing and increasing levels of herbivory, and estimate what proportion of the carbon involved in the leaf area loss is transferred further in the food web. We estimate that the net primary productivity of an oak stand under ambient levels of herbivory is 54 - 69% lower than that of a completely intact stand. The effect of herbivory quantified only as leaf area loss (0.1 Mg C ha-1 yr-1) is considerably smaller than when the effects of herbivory on leaf physiology are included (8.5 Mg C ha-1 yr-1). We propose that the effect of herbivory on primary productivity is non-linear and mainly determined by changes in leaf gas exchange. We call for replicated studies in other systems to validate the relationship between insect herbivory and ecosystem productivity described here.

category: ecology


title: The evolutionary making of SARS-CoV-2

abstract: A mechanistic understanding of how SARS-CoV-2 (sarbecovirus, betacoronavirus) infects human cells is emerging, but the evolutionary trajectory that gave rise to this pathogen is poorly understood. Here we scan SARS-CoV-2 protein sequences in-silico for innovations along the evolutionary lineage starting with the last common ancestor of coronaviruses. SARS-CoV-2 substantially differs from viruses outside sarbecovirus both in its set of encoded proteins and in their domain architectures, indicating divergent functional demands. Within sarbecoviruses, sub-domain level profiling using predicted linear epitopes reveals how the primary interface between host cell and virus, the spike, was gradually reshaped. The only epitope that is private to SARS-CoV-2 overlaps with the furin cleavage site, a "switch" that modulates spikes conformational landscape in response to host-cell interaction. This cleavage site has fundamental relevance for both immune evasion and cell infection, and the apparently ongoing evolutionary fine-tuning of its use by SARS-CoV-2 should be monitored.


title: Estimating repeat spectra and genome length from low-coverage genome skims with RESPECT

abstract: The cost of sequencing the genome is dropping at a much faster rate compared to assembling and finishing the genome. The use of lightly sampled genomes (genome-skims) could be transformative for genomic ecology, and results using k-mers have shown the advantage of this approach in identification and phylogenetic placement of eukaryotic species. Here, we revisit the basic question of estimating genomic parameters such as genome length, coverage, and repeat structure, focusing specifically on estimating the k-mer repeat spectrum. We show using a mix of theoretical and empirical analysis that there are fundamental limitations to estimating the k-mer spectra due to ill-conditioned systems, and that has implications for other genomic parameters. We get around this problem using a novel constrained optimization approach (Spline Linear Programming), where the constraints are learned empirically. On reads simulated at 1X coverage from 66 genomes, our method, REPeat SPECTra Estimation (RESPECT), had < 1.5% error in length estimation compared to 34% error previously achieved. In shotgun sequenced read samples with contaminants, RESPECT length estimates had median error 4%, in contrast to other methods that had median error 80%. Together, the results suggest that low-pass genomic sequencing can yield reliable estimates of the length and repeat content of the genome. The RESPECT software will be publicly available at https://github.com/shahab-sarmashghi/RESPECT.git


title: Multiomics single-cell analysis of human pancreatic islets reveals novel cellular states in health and type 1 diabetes

abstract: Type 1 diabetes (T1D) is an autoimmune disease of only partially defined etiology in which immune cells destroy insulin-producing beta cells. Using single-cell transcriptomics and an advanced analytical strategy to assess pancreatic islets of T1D, autoantibody-positive, and non-diabetic organ donors, we identified both canonical cell types and rare insulin-expressing cells with a hybrid mixture of endocrine and exocrine gene signatures within all donors. We further found elevated expression of MHC Class II pathway genes in exocrine ductal cells of T1D donors, which we confirmed through CyTOF, in situ imaging mass cytometry, and immunofluorescence analysis. Taken together, our multimodal analyses identify novel cell types and processes that may contribute to T1D immunopathogenesis and provide new cellular and molecular insights into human pancreas function.


title: The regulatory landscape of cells in the developing mouse cerebellum

abstract: Organ development is orchestrated by cell- and time-specific gene regulatory networks. Here we investigated the regulatory basis of mouse cerebellum development from early neurogenesis to adulthood. By acquiring snATAC-seq profiles for ~90,000 cells spanning eleven stages, we mapped all major cerebellar cell types and identified candidate cis-regulatory elements (CREs). We detected extensive spatiotemporal heterogeneity among progenitor cells and characterized the regulatory programs underlying the differentiation of cerebellar neurons. Although CRE activity is predominantly cell type- and time-specific, periods of greater regulatory change are shared across cell types. There is a universal decrease in CRE conservation and pleiotropy during development and differentiation, but the degree of evolutionary constraint differs between cerebellar cell types. Our work delineates the developmental and evolutionary dynamics of gene regulation in cerebellar cells and provides general insights into mammalian organ development.


title: NUCLEIC TRANSFORMER: DEEP LEARNING ONNUCLEIC ACIDS WITH SELF-ATTENTION ANDCONVOLUTIONS

abstract: AO_SCPLOWBSTRACTC_SCPLOWMuch work has been done to apply machine learning and deep learning to genomics tasks, but these applications usually require extensive domain knowledge and the resulting models provide very limited interpretability. Here we present the Nucleic Transformer, a conceptually simple but effective and interpretable model architecture that excels in a variety of DNA/RNA tasks. The Nucleic Transformer processes nucleic acid sequences with self-attention and convolutions, two deep learning techniques that have proved dominant in the fields of computer vision and natural language processing. We demonstrate that the Nucleic Transformer can be trained in both supervised and unsupervised fashion without much domain knowledge to achieve high performance with limited amounts of data in Escherichia coli promoter classification, viral genome identification, and degradation properties of COVID-19 mRNA vaccine candidates. Additionally, we showcase extraction of promoter motifs from learned attention and how direct visualization of self-attention maps assists informed decision making using deep learning models.


title: Nucleic Transformer: Deep Learning on Nucleic Acids with Self-attention and Convolutions


title: SMILE: Mutual Information Learning for Integration of Single Cell Omics Data

abstract: Deep learning approaches have empowered single-cell omics data analysis in many ways, generating new insights from complex cellular systems. As there is an increasing need for single cell omics data to be integrated across sources, types, and features of data, the challenges of integrating single-cell omics data are rising. Here, we present a deep clustering algorithm that learns discriminative representation for single-cell data via maximizing mutual information, SMILE (Single-cell Mutual Information Learning). Using a unique cell-pairing design, SMILE successfully integrates multi-source single-cell transcriptome data, removing batch effects and projecting similar cell types, even from different tissues, into the same representation space. SMILE can also integrate data from two or more modalities, such as joint profiling technologies using singlecell ATAC-seq, RNA-seq, DNA methylation, Hi-C, and ChIP data. SMILE works well even when feature types are unmatched, such as genes for RNA-seq and genome wide peaks for ATAC-seq.


title: Strobemers: an alternative to k-mers for sequence comparison

abstract:
K-mer-based methods are widely used in bioinformatics for various types of sequence comparison. However, a single mutation will mutate k consecutive k-mers and makes most k-mer based applications for sequence comparison sensitive to variable mutation rates. Many techniques have been studied to overcome this sensitivity, e.g., spaced k-mers and k-mer permutation techniques, but these techniques do not handle indels well. For indels, pairs or groups of small k-mers are commonly used, but these methods first produce k-mer matches, and only in a second step, a pairing or grouping of k-mers is performed. Such techniques produce many redundant k-mer matches due to the size of k.

Here, we propose strobemers as an alternative to k-mers for sequence comparison. Intuitively, strobemers consist of linked minimizers. We use simulated data to show that strobemers provide more evenly distributed sequence matches and are less sensitive to different mutation rates than k-mers and spaced k-mers. Strobemers also produce a higher match coverage across sequences. We further implement a proof-of-concept sequence matching tool StrobeMap, and use synthetic and biological Oxford Nanopore sequencing data to show the utility of using strobemers for sequence comparison in different contexts such as sequence clustering and alignment scenarios. A reference implementation of our tool StrobeMap together with code for analyses is available at https://github.com/ksahlin/strobemers.


title: SUFI: An automated approach to spectral unmixing of fluorescent biological images

abstract: Multispectral fluorescence imaging coupled with linear unmixing is a form of image data collection and analysis that uses multiple fluorescent dyes - each measuring a specific biological signal - that are simultaneously measured and subsequently "unmixed" to provide a read-out for each individual signal. This strategy allows for measuring multiple signals in a single data capture session - for example, multiple proteins or RNAs in tissue slices or cultured cells, but can often result in mixed signals and bleed-through problems across dyes. Existing spectral unmixing algorithms are not optimized for challenging biological specimens such as postmortem human brain tissue, and often require manual intervention to extract spectral signatures. We therefore developed an intuitive, automated, and flexible package called SUFI: spectral unmixing of fluorescent images (https://github.com/LieberInstitute/SUFI). This package unmixes multispectral fluorescence images by automating the extraction of spectral signatures using Vertex Component Analysis, and then performs one of three unmixing algorithms derived from remote sensing. We demonstrate these remote sensing algorithms performance on four unique biological datasets and compare the results to unmixing results obtained using ZEN Black software (Zeiss). We lastly integrate our unmixing pipeline into the computational tool dotdotdot that is used to quantify individual RNA transcripts at single cell resolution in intact tissues and perform differential expression analysis of smFISH data, and thereby provide a one-stop solution for multispectral fluorescence image analysis and quantification. In summary, we provide a robust, automated pipeline to assist biologists with improved spectral unmixing of multispectral fluorescence images.


title: SUFI: An automated approach to spectral unmixing of fluorescent multiplex images captured in mouse and postmortem human brain tissues


title: Development of spike receptor-binding domain nanoparticle as a vaccine candidate against SARS-CoV-2 infection in ferrets

abstract: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a causative agent of COVID-19 pandemic, enters host cells via the interaction of its Receptor-Binding Domain (RBD) of Spike protein with host Angiotensin-Converting Enzyme 2 (ACE2). Therefore, RBD is a promising vaccine target to induce protective immunity against SARS-CoV-2 infection. In this study, we report the development of RBD protein-based vaccine candidate against SARS-CoV-2 using self-assembling H. pylori-bullfrog ferritin nanoparticles as an antigen delivery. RBD-ferritin protein purified from mammalian cells efficiently assembled into 24-mer nanoparticles. 16-20 months-old ferrets were vaccinated with RBD-ferritin nanoparticles (RBD-nanoparticles) by intramuscular or intranasal inoculation. All vaccinated ferrets with RBD-nanoparticles produced potent neutralizing antibodies against SARS-CoV-2. Strikingly, vaccinated ferrets demonstrated efficient protection from SARS-CoV-2 challenge, showing no fever, body weight loss and clinical symptoms. Furthermore, vaccinated ferrets showed rapid clearance of infectious viruses in nasal washes and lungs as well as viral RNA in respiratory organs. This study demonstrates the Spike RBD-nanoparticle as an effective protein vaccine candidate against SARS-CoV-2.


title: Pseudomonas can prevent the parasitic fungus, while keeping the crop fungus unaffected, in the gardens of Odontotermes obesus.

abstract: Insects that farm monocultures of fungi are canonical examples of nutritional symbiosis as well as independent evolution of agriculture in non-human animals. But just like in human agriculture, these fungal crops face constant threat of invasion by weeds which, if unchecked, takes over the crop fungus. In fungus-growing termites, the crop fungus (Termitomyces) faces such challenges from the parasitic fungus Pseudoxylaria. The mechanism by which Pseudoxylaria is suppressed is not known. However, evidence suggests that some bacterial secondary symbionts can serve as defensive mutualists by preventing the growth of Pseudoxylaria. However, such secondary symbionts must possess the dual, yet contrasting, capabilities of suppressing the weedy fungus while keeping the growth of the crop fungus unaffected. This study describes the isolation, identification and culture-dependent estimation of the roles of several such putative defensive mutualists from the colonies of the wide-spread fungus-growing termite from India, Odontotermes obesus. From the 38 bacterial cultures tested, a strain of Pseudomonas showed significantly greater suppression of the weedy fungus than the crop fungus. Moreover, a 16S rRNA pan-microbiome survey, using the Nanopore platform, revealed Pseudomonas to be a part of the core microbiota of Odontotermes obesus. A meta-analysis of microbiota composition across different species of Odontotermes also confirms the wide-spread prevalence of Pseudomonas within this termite. These evidence indicate that Pseudomonas could be playing the role of defensive mutualist within Odontotermes.


title: Long-read assembly and comparative evidence-based reanalysis of Cryptosporidium genome sequences reveal new biological insights

abstract: Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community still relies on a fragmented reference genome sequence from 2004. Incomplete reference sequences hamper experimental design and interpretation. We have generated a new C. parvum IOWA genome assembly supported by PacBio and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species C. parvum, C. hominis and C. tyzzeri. The new C. parvum IOWA reference genome assembly is larger, gap free and lacks ambiguous bases. This chromosomal assembly recovers 13 of 16 possible telomeres and raises a new hypothesis for the remaining telomeres and associated subtelomeric regions. Comparative annotation revealed that most "missing" orthologs are found suggesting that species differences result primarily from structural rearrangements, gene copy number variation and SNVs in C. parvum, C. hominis and C. tyzzeri. We made >1,500 C. parvum annotation updates based on experimental evidence. They included new transporters, ncRNAs, introns and altered gene structures. The new assembly and annotation revealed a complete DNA methylase Dnmt2 ortholog. 190 genes under positive selection including many new candidates were identified using the new assembly and annotation as reference. Finally, possible subtelomeric amplification and variation events in C. parvum are detected that reveal a new level of genome plasticity that will both inform and impact future research.


title: The effects of maternal social connectivity and integration on offspring survival in a marmot

abstract:
In social species, maternal social relationships, in addition to direct care, impact offspring survival but much of what we know about these effects comes from studies of obligately social and cooperatively breeding species. Yellow-bellied marmots (Marmota flaviventer) are a facultatively social species whose social groups vary in composition, size, and cohesiveness. This natural variation in sociality and cooperative breeding behavior makes yellow-bellied marmots an ideal species within which to study the effects of maternal affiliative and agonistic social behavior on offspring. We used social network analysis to investigate the relationship between maternal social connectivity and integration on offspring summer and yearly survival, with the hypothesis that offspring with more affiliative mothers are more likely to survive than the offspring of more agonistic mothers. However, we found the inverse to be true: pups born to mothers who received more affiliative interactions were less likely to survive while the offspring of mothers who were more highly integrated into agonistic networks had enhanced survival. Overall, maternal social network measures were positively and negatively correlated with offspring survival, indicating that pups are influenced by their mothers social world, often in contradictory ways. Relative predation risk and colony location also mediated the effects of social relationships on pup survival. This study contributes to a small but growing body of work that demonstrates that specific attributes of sociality have specific consequences and that by adopting an attribute-focused view of sociality we are better able to understand how environmental conditions mediate the costs and benefits of sociality.

Lay SummaryMaternal social relationships can impact offspring survival but much of what we know about these effects comes from studies of obligately social species. In faculatively social yellow-bellied marmots we found that pups born to mothers who received more affiliative interactions were less likely to survive while the offspring of mothers who were more highly integrated into agonistic networks had enhanced survival. Overall, pups are influenced by their mothers social world, often in contradictory ways.

category: animal behavior and cognition


title: Gene network module changes associated with the vertebrate fin to limb transition

abstract: Evolutionary phenotypic transitions, such as the fin to limb transition in vertebrate evolution, result from changes in associated genes and their interactions, often in response to changing environment. Identifying the associated changes in gene networks is vital to achieve a better understanding of these transitions. Previous experimental studies have been typically limited to manipulating a small number of genes. To expand the number of analyzed genes and hence, biological knowledge, we computationally isolated and compared the gene modules for paired fins (pectoral fin, pelvic fin) of fishes (zebrafish) to those of the paired limbs (forelimb, hindlimb) of mammals (mouse) using quality-enhanced gene networks from zebrafish and mouse. We ranked module genes according to their weighted-degrees and identified the highest-ranking hub genes, which were important for the module stability. Further, we identified genes conserved during the fin to limb transition and investigated the fates of zebrafish-specific and mouse-specific module genes in relation to their involvements in newly emerged or lost anatomical structures during the aquatic to terrestrial vertebrate transition. This paper presents the results of our investigations and demonstrates a general network-based computational workflow to study evolutionary phenotypic transitions involving diverse model organisms and anatomical entities.


title: Recombination and low-diversity confound homoplasy-based methods to detect the effect of SARS-CoV-2 mutations on viral transmissibility

abstract: The SARS-CoV-2 variant carrying the Spike protein mutation G614 was first detected in late January 2020 and within a few months became the dominant form globally. In the months that followed, many studies, both in vitro and in animal models, showed that variants carrying this mutation were more infectious and more readily transmitted than the ancestral Wuhan form. Here we investigate why a recently published study by van Dorp et al. failed to detect such higher transmissibility of the G614 variant using homoplasy-based methods. We show that both low diversity and recombination confound the methods utilized by van Dorp et al. and significantly decrease their sensitivity. Furthermore, though they claim no evidence of recombination in their dataset, we and several other studies identify a subset of the sequences as recombinants, possibly enough to affect their statistic adversely.


title: A genetic approach to dissect the role of prefoldins in Arabidopsis

abstract: The prefoldin complex (PFDc) was identified in humans as co-chaperone of the cytosolic chaperonin TRiC/CCT. It is conserved in eukaryotes and is composed of subunits PFD1 to 6. PFDc-TRiC/CCT operates folding actin and tubulins. In addition to this function, PFDs participate in a wide range of cellular processes, both in the cytoplasm and in the nucleus, and their malfunction cause developmental alterations and disease in animals, and altered growth and environmental responses in yeast and plants. Genetic analyses in yeast indicate that not all functions performed by PFDs require the participation of the canonical complex. The lack of systematic genetic analyses in higher eukaryotes makes it difficult to discern whether PFDs participate in a particular process as canonical complex or in alternative configurations, i.e. as individual subunits or in other complexes. To tackle this question, and on the premise that the canonical complex cannot be formed if one subunit is missing, we have prepared an Arabidopsis mutant deficient in the six prefoldins, and compared various growth and environmental responses with those of the individual pfd. In this way, we demonstrate that the PFDc is required to delay flowering, for seed germination, or to respond to high salt stress, whereas two or more PFDs redundantly attenuate the response to osmotic stress. A coexpression analysis of differentially expressed genes in the sextuple mutant has identified several transcription factors, such as ABI5 or PIF4, acting downstream of PFDs. Furthermore, it has made possible to assign novel roles for PFDs, for instance, in the response to warm temperature.

category: plant biology


title: The composition and function of Enterococcus faecalismembrane vesicles

abstract: Membrane vesicles (MVs) contribute to various biological processes in bacteria, including virulence factor delivery, host immune evasion, and cross-species communication. MVs are frequently being discharged from the surface of both Gram-negative and Gram-positive bacteria during growth. In some Gram-positive bacteria, genes affecting MV biogenesis have been identified, but the mechanism of MV formation is unknown. In Enterococcus faecalis, a causative agent of life-threatening bacteraemia and endocarditis, neither mechanisms of MV formation nor their role in virulence has been examined. Since MVs of many bacterial species are implicated in host-pathogen interactions, biofilm formation, horizontal gene transfer, and virulence factor secretion in other species, we sought to identify, describe, and functionally characterize MVs from E. faecalis. Here we show that E. faecalis releases MVs that possess unique lipid and protein profiles, distinct from the intact cell membrane, and are enriched in lipoproteins. MVs of E. faecalis are specifically enriched in unsaturated lipids that might provide membrane flexibility to enable MV formation, providing the first insights into the mechanism of MV formation in this Gram-positive organism.

category: microbiology


title: Comparative Genomics and Environmental Distribution of Large dsDNA viruses in the family Asfarviridae

abstract: The Asfarviridae is a family of Nucleo-Cytoplasmic Large DNA Viruses (NCLDV) of which African swine fever virus (ASFV) is the most well-characterized. Recently the discovery of several Asfarviridae members other than ASFV has suggested that this family represents a diverse and cosmopolitan group of viruses, but the genomics and distribution of this family have not been studied in detail. To this end we analyzed five complete genomes and 35 metagenome-assembled genomes (MAGs) of viruses from this family to shed light on their evolutionary relationships and environmental distribution. The Asfarvirus MAGs derive from diverse marine, freshwater, and terrestrial habitats, underscoring the broad environmental distribution of this family. We present phylogenetic analyses using conserved marker genes and whole-genome comparison of pairwise average amino acid identity values, revealing a high level of genomic divergence across disparate Asfarviruses. Further, we found that Asfarviridae genomes encode genes with diverse predicted metabolic roles and detectable sequence homology to proteins in bacteria, archaea, and different eukaryotes, highlighting the genomic chimerism that is a salient feature of NCLDV. Our read mapping from Tara oceans metagenomic data also revealed that three Asfarviridae MAGs were present in multiple marine samples, indicating that they are widespread in the ocean. In one of these MAGs we identified four marker genes with >95% amino acid identity to genes sequenced from a virus that infects the dinoflagellate Heterocapsa circularisquama (HcDNAV). This suggests a potential host for this MAG, which would thereby represent a near-complete genome of a dinoflagellate-infecting giant virus. Together, these results show that Asfarviridae are ubiquitous, comprise similar sequence divergence as other NCLDV families, and include several members that are widespread in the ocean and potentially infect ecologically important protists.


title: Acidovorax pan-genome reveals specific functional traits for plant beneficial and pathogenic plant-associations

abstract: Beta-proteobacteria belonging to the genus Acidovorax have been described from various environments. Many strains can interact with a range of hosts, including humans and plants, forming neutral, beneficial or detrimental associations. In the frame of this study we investigated genomic properties of 52 bacterial strains of the genus Acidovorax, isolated from healthy roots of Lotus japonicus, with the intent of identifying traits important for effective plant-growth promotion. Based on single-strain inoculation bioassays with L. japonicus, performed in a gnotobiotic system, we could distinguish 7 robust plant-growth promoting strains from strains with no significant effects on plant-growth. We could show that the genomes of the two groups differed prominently in protein families linked to sensing and transport of organic acids, production of phytohormones, as well as resistance and production of compounds with antimicrobial properties. In a second step, we compared the genomes of the tested isolates with those of plant pathogens and free-living strains of the genus Acidovorax sourced from public repositories. Our pan-genomics comparison revealed features correlated with commensal and pathogenic lifestyle. We could show that commensals and pathogens differ mostly in their ability to use plant-derived lipids and in the type of secretion-systems being present. Most free-living Acidovorax strains did not harbour any secretion-systems. Overall, our data indicated that Acidovorax strains undergo extensive adaptations to their particular lifestyle by horizontal uptake of novel genetic information and loss of unnecessary genes.


title: From heterogeneity to homogeneity: coordination of siderophore gene expression among clonal cells of the bacterium Pseudomonas aeruginosa

abstract: There has been great progress in understanding how bacterial groups coordinate social actions, such as biofilm formation, swarming and public-goods secretion. Less clear, however, is whether the seemingly coordinated responses observed at the group level actually mirror what individual cells do. Here, we use a microscopy approach to simultaneously quantify the investment of individual cells of the bacterium Pseudomonas aeruginosa into two public goods, the siderophores pyochelin and pyoverdine. Using gene expression as a proxy for investment, we initially observed no coordination but high heterogeneity and bimodality in siderophore gene expression across cells. With increasing cell density, gene expression became homogenized across cells, accompanied by a shift from pyochelin to pyoverdine expression. We found positive correlations in the expression of pyochelin and pyoverdine genes across cells, and show that cell-to-cell variation is driven by differences in cellular metabolic states. We propose a model explaining how variation in internal iron stocks can spur initial erratic gene expression, while siderophore-mediated signalling and intra-cellular feedbacks later on can induce highly coordinated gene expression and synchronized shifts from pyochelin to pyoverdine. Our work provides new insights into bacterial collective decision-making processes and reveals a three-phase chronobiological siderophore investment cycle in P. aeruginosa.


title: Hyperpiliation reduces Pseudomonas aeruginosa pathogenicity by impeding type III secretion

abstract:
Type IVa pili (T4aP) are important virulence factors for many bacterial pathogens. Previous studies suggested that the retraction ATPase, PilT, modulates pathogenicity due to its critical role in pilus dynamics and twitching motility. Here we use a Caenorhabditis elegans slow killing model to show that hyperpiliation, not loss of pilus retraction, reduces virulence of Pseudomonas aeruginosa strains PAK and PA14 by interfering with function of the contact-dependent type III secretion system (T3SS). Hyperactivating point mutations in the P. aeruginosa PilSR two-component system that controls transcription of the major pilin gene, pilA, increased levels of surface pili to the same extent as deleting pilT, without impairing twitching motility. These functionally hyperpiliated PilSR mutants had significant defects in pathogenicity that were rescued by deleting pilA or by increasing the length of T3SS needles via deletion of the needle-length regulator, PscP. Hyperpiliated pilT deletion or pilO point mutants showed similar PilA-dependent impairments in virulence, validating the phenotype. Together, our data support a model where a surfeit of pili prevents effective engagement of contact-dependent virulence factors. These findings suggest that the role of T4aP retraction in virulence should be revised.

SIGNIFICANCEPseudomonas aeruginosa is a major contributor to hospital-acquired infections and particularly problematic due to its intrinsic resistance to many front-line antibiotics. Strategies to combat this and other important pathogens include development of anti-virulence therapeutics. We show that the pathogenicity of P. aeruginosa is impaired when the amount of type IVa pili (T4aP) expressed on the cell surface increases, independent of the bacterias ability to twitch. We propose that having excess T4aP on the cell surface can physically interfere with productive engagement of the contact-dependent type III secretion toxin delivery system. A better understanding of how T4aP modulate interaction of bacteria with target cells will improve the design of therapeutics targeting components involved in regulation of T4aP expression and function, to reduce the clinical burden of P. aeruginosa and other T4aP-expressing bacteria.


title: Circadian dynamics of the teleost skin immune-microbiome interface

abstract: Circadian rhythms of host immune activity and their microbiomes are likely pivotal to health and disease resistance. The integration of chronotherapeutic approaches to disease mitigation in managed animals, however, is yet to be realised. In aquaculture, light manipulation is commonly used to enhance growth and control reproduction but may have unknown negative consequences for animal health. Infectious diseases are a major barrier to sustainable aquaculture and understanding the circadian dynamics of fish immunity and crosstalk with the microbiome is urgently needed. We demonstrate daily rhythms in fish skin immune expression and microbiomes, that are modulated by photoperiod and parasitic infection. We identify putative associations of host clock and immune gene profiles with microbial composition. Our results suggest circadian perturbation that shifts the magnitude and timing of immune and microbiota activity, is detrimental to fish health. This study represents a valuable foundation for investigating the utility of chronotherapies in aquaculture, and more broadly contributes to our understanding of circadian health in vertebrates.


title: The DEAD-box RNA helicases RhlE2 is a global regulator of Pseudomonas aeruginosa lifestyle and pathogenesis

abstract: The RhlE DEAD-box RNA helicase protein family is widespread among Proteobacteria, but it is the least understood due to the lack of a clear biological function. Here, we study the two RhlE homologs present in the opportunistic pathogen Pseudomonas aeruginosa. RhlE1 and RhlE2 diverged during P. aeruginosa evolution; our data indicate that this resulted in a non-redundant biological role, a distinct molecular action and an enzymatic activity differentially stimulated by RNA. Whereas RhlE1 is specifically necessary for bacteria growth in cold, we show that RhlE2 acts as global post-transcriptional regulator, affecting the level of hundreds of cellular transcripts and multiple functionalities indispensable not only for P. aeruginosa environmental adaptation, but also for its virulence. The global action of RhlE2 relies on a unique C-terminal extension, which establishes an RNA-dependent interaction with the RNase E endonuclease and the cellular RNA degradation machinery.


title: Effect of concentration and hydraulic reaction time on the removal of pharmaceutical compounds in a membrane bioreactor inoculated with activated sludge

abstract: Pharmaceuticals are often not fully removed in wastewater treatment plants (WWTPs) and are thus being detected at trace levels in water bodies all over the world posing a risk to numerous organisms. These organic micropollutants (OMPs) reach WWTPs at concentrations sometimes too low to serve as growth substrate for microorganisms, thus co-metabolism is thought to be the main conversion mechanism. In this study, the microbial removal of six pharmaceuticals was investigated in a membrane bioreactor at increasing concentrations (4-800 nM) of the compounds and using three different hydraulic retention times (HRT; 1, 3.5, 5 days). The bioreactor was inoculated with activated sludge from a Dutch WWTP and fed with ammonium, acetate, and methanol as main growth substrates to stimulate and mimic co-metabolism in a WWTP. Each pharmaceutical compound had a different average removal efficiency: acetaminophen (100%) > fluoxetine (50%) > metoprolol (25%) > diclofenac (20%) > metformin (15%) > carbamazepine (10%). Higher pharmaceutical influent concentrations proportionally increased the removal rate of each compound, but surprisingly not the removal percentage. Furthermore, only metformin removal improved to 80-100% when HRT or biomass concentration was increased in the reactor. Microbial community changes were followed with 16S rRNA gene amplicon sequencing in response to the increment of supplied pharmaceutical concentration: it was found that Nitrospirae and Planctomycetes 16S rRNA relative gene abundance decreased, whereas Acidobacteria and Bacteroidetes increased. Remarkably, the Dokdonella genus, previously implicated in acetaminophen metabolism, showed a 30-fold increase in abundance at the highest (800 nM) concentration of pharmaceuticals applied. Taken together, these results suggest that the incomplete removal of most pharmaceutical compounds in WWTPs is neither dependent on concentration nor HRT. Accordingly, we propose a chemical equilibrium or a growth substrate limitation as the responsible mechanisms of the incomplete removal. Finally, Dokdonella could be the main acetaminophen degrader under activated sludge conditions, and non-antimicrobial pharmaceuticals might still be toxic to relevant WWTP bacteria.


title: Resurrection of a viral internal ribosome entry site from a 700 year old ancient Northwest Territories cripavirus

abstract:
The dicistrovirus intergenic region internal ribosome entry site (IGR IRES) uses an unprecedented streamlined mechanism whereby the IRES adopts a triple-pseudoknot (PK) structure to directly bind to the conserved core of the ribosome and drive translation from a non-AUG codon. The origin of this IRES mechanism is not known. Previously, a partial fragment of a divergent dicistrovirus RNA genome, named ancient Northwest territories cripavirus (aNCV), was extracted from 700-year-old caribou feces trapped in a subarctic ice patch. Structural prediction of the aNCV IGR sequence generated a secondary structure similar to contemporary IGR IRES structures. There are, however, subtle differences including 105 nucleotides upstream of the IRES of unknown function. Using filter binding assays, we showed that the aNCV IGR IRES could bind to purified salt-washed human ribosomes and compete with a prototypical IGR IRES for ribosomes. Toeprinting analysis using primer extension pinpointed the putative start site of the aNCV IGR at a GCU alanine codon adjacent to PKI. Using a bicistronic reporter RNA, the aNCV IGR IRES can direct internal ribosome entry in vitro in a manner dependent on the integrity of the PKI domain. Lastly, we generated a chimeric virus clone by swapping the aNCV IRES into the cricket paralysis virus infectious clone. The chimeric infectious clone with an aNCV IGR IRES supported translation and virus infection. The characterization and resurrection of a functional IGR IRES from a divergent 700-year-old virus provides a historical framework in the importance of this viral translational mechanism.

IMPORTANCEInternal ribosome entry sites are RNA structures that are used by some positive-sense monopartite RNA viruses to drive viral protein synthesis. The origin of internal ribosome entry sites is not known. Using biochemical approaches, we demonstrate that an RNA structure from an ancient viral genome that was discovered from a 700-year-old caribou feces trapped in subarctic ice is functionally similar to modern internal ribosome entry sites. We resurrect this ancient RNA mechanism by demonstrating that it can support virus infection in a contemporary virus clone, thus providing insights into the origin and evolution of this viral strategy.


title: Spatial organization of transcript elongation and splicing kinetics

abstract: The organization of the genome in three-dimensional space has been shown to play an important role in gene expression. Specifically, facets of genomic interaction such as topologically associated domains (TADs) have been shown to regulate transcription by bringing regulatory elements into close proximity1. mRNA production is an intricate process with multiple control points including regulation of Pol II elongation and the removal of non-coding sequences via pre-mRNA splicing2. The connection between genomic compartments and the kinetics of RNA biogenesis and processing has been largely unexplored. Here, we measure Pol II elongation and splicing kinetics genome-wide using a novel technique that couples nascent RNA-seq with a mathematical model of transcription and co-transcriptional RNA processing. We uncovered multiple layers of spatial organization of these rates: the rate of splicing is coordinated across introns within individual genes, and both elongation and splicing rates are coordinated within TADs, as are alternative splicing outcomes. Overall, our work establishes that the kinetics of transcription and splicing are coordinated by the spatial organization of the genome and suggests that TADs are a major platform for coordination of alternative splicing.


title: Generation of a novel Nkx6-1 Venus fusion reporter mouse line

abstract: Nkx6-1 is a member of the Nkx family of homeodomain transcription factors (TF) that regulates motor neuron development, neuron specification and pancreatic endocrine and {beta}-cell differentiation. To facilitate the isolation and tracking of Nkx6-1-expressing cells, we have generated a novel Nkx6-1 Venus Fusion (Nkx6-1-VF) reporter allele. The Nkx6-1-VF knock-in reporter is regulated by endogenous cis-regulatory elements of Nkx6-1 and the fluorescent protein fusion does not interfere with the TF function, as homozygous mice are viable and fertile. In addition, the nuclear localization of Nkx6-1-VF protein reflects the endogenous Nkx6-1 protein distribution. During embryonic pancreas development, the reporter protein marks the pancreatic ductal progenitors and the endocrine lineage, but is absent in the exocrine compartment. Moreover, the levels of Nkx6-1-VF reporter is upregulated upon {beta}-cell differentiation during the major wave of endocrinogenesis. In the adult islets of Langerhans, the reporter protein is exclusively found in insulin-secreting {beta}-cells. Importantly, the Venus reporter activities allows successful tracking of {beta}-cells in live-cell imaging and their specific isolation by flow sorting. In summary, the generation of Nkx6-1-VF reporter line provides a unique tool to study the spatio-temporal expression pattern of this TF during organ development and isolate and track Nkx6-1-expressing cells such as pancreatic {beta}-cells, but also neurons and motor neurons in health and disease.

category: developmental biology


title: Biodiversity Image Quality Metadata Augments Convolutional Neural Network Classification of Fish Species

abstract: Biodiversity image repositories are crucial sources of training data for machine learning approaches to biological research. Metadata, specifically metadata about object quality, is putatively an important prerequisite to selecting sample subsets for these experiments. This study demonstrates the importance of image quality metadata to a species classification experiment involving a corpus of 1935 fish specimen images which were annotated with 22 metadata quality properties. A small subset of high quality images produced an F1 accuracy of 0.41 compared to 0.35 for a taxonomically matched subset of low quality images when used by a convolutional neural network approach to species identification. Using the full corpus of images revealed that image quality differed between correctly classified and misclassified images. We found the visibility of all anatomical features was the most important quality feature for classification accuracy. We suggest biodiversity image repositories consider adopting a minimal set of image quality metadata to support future machine learning projects.


title: Reliance on polyfunctional tissue leads to a reproduction-immunity tradeoff due to inherent constraint

abstract: The use of one tissue for multiple purposes can result in constraints, impaired function, and tradeoffs. The insect fat body performs remarkably diverse functions including metabolic control, reproductive provisioning, and systemic immune responses. Immunity and reproduction are observed to trade off in many organisms, although the mechanistic basis for the tradeoff is generally unknown. More generally, how do polyfunctional tissues simultaneously execute multiple distinct physiological functions? Using single-nucleus sequencing, we determined the Drosophila melanogaster fat body executes diverse basal functions with heterogenous cellular subpopulations. However, as an emergency function, the immune response engages the entire tissue. We found that reproductively active females exhibit impaired capacity to produce new protein in response to infection, resulting in the reproduction-immunity tradeoff. We suggest that such inherent internal limitations may provide a general explanation for the wide prevalence of physiological and evolutionary tradeoffs.


title: The genomes of precision edited cloned calves show no evidence for off-target events or increased de novo mutagenesis

abstract: Animal health and welfare are at the forefront of public concern and the agricultural sector is responding by prioritising the selection of welfare-relevant traits in their breeding schemes. In some cases, welfare-enhancing traits such as horn-status (i.e., polled) or diluted coat colour, which could enhance heat tolerance, may not segregate in breeds of primary interest, highlighting gene-editing tools such as the CRISPR-Cas9 technology as an approach to rapidly introduce variation into these populations. A major limitation preventing the acceptance of CRISPR-Cas9 mediated gene-editing, however, is the potential for off-target mutagenesis, which has raised concerns about the safety and ultimate applicability of this technology. Here, we present a clone-based study design that has allowed a detailed investigation of off-target and de novo mutagenesis in a cattle line bearing edits in the PMEL gene for diluted coat-colour. No off-target events were detected from high depth whole genome sequencing performed in precursor cell-lines and resultant calves cloned from those edited and non-edited cell lines. Long molecule sequencing at the edited site and plasmid-specific PCRs did not reveal structural variations and/or plasmid integration events in edited samples. Furthermore, an in-depth analysis of de novo mutations across samples revealed that the mutation frequency and spectra were unaffected by editing status. Cells in culture, however, had a distinct mutation signature where de novo mutations were predominantly C>A mutations, and in cloned calves they were predominantly T>G mutations, deviating from the expected excess of C>T mutations. We conclude that the gene-edited cells and calves in this study did not present a higher mutation load than unedited controls. Cell culture and somatic cell nuclear transfer cloning processes contributed the major source of contrast in mutational profile between samples.


title: Resequencing 250 soybean accessions: new insights into genes associated with agronomic traits and genetic networks

abstract: Limited knowledge on genomic diversity and the functional genes associated with soybean variety traits has resulted in slow breeding progress. We sequenced the genome of 250 soybean landraces and cultivars from China, America and Europe, and investigated their population structure, genetic diversity and architecture and selective sweep regions of accessions. We identified five novel agronomically important genes and studied the effects of functional mutations in respective genes. We found candidate genes GSTT1, GL3 and GSTL3 associated with isoflavone content, CKX3 associated with yield traits, and CYP85A2 associated with both architecture and yield traits. Our phenotype-gene network analysis revealed that hub nodes play a role in complex phenotypic associations. In this work, we describe novel agronomic trait associated genes and a complex genetic network, providing a valuable resource for future soybean molecular breeding.


title: Genomic and phenotypic biology of novel strains of Dickeya zeae isolated from pineapple and taro in Hawaii: insights into genome plasticity, pathogenicity, and virulence determinants

abstract: Dickeya zeae, a bacterial plant pathogen in the family Pectobacteriaceae, is responsible for a wide range of diseases on potato, maize, rice, banana, pineapple, taro and ornamentals and significantly reduces crop production; D. zeae causes soft rot of taro (Colocasia esculenta) and heart rot of pineapple (Ananas comosus). In this study, we used Pacific Biosciences to sequence two high quality complete genomes of novel strains of D. zeae: PL65 (size - 4.74997 MB; depth - 701; GC - 53.3%) and A5410 (size - 4.7792 MB; depth - 558; GC - 53.6%) isolated from economically important Hawaiian crops, taro and pineapple, respectively. Additional complete genomes of D. zeae representing two additional hosts (rice and banana), and other species for taxonomic comparison, were retrieved from the NCBI GenBank genome database. The genomic analyses indicated truncated type III and IV secretion systems (T3SS and T4SS) in the taro strain, which only harbors 1 and 2 genes of T3SS and T4SS, respectively, and showed high heterogeneity in the type VI secretion system. Unlike the EC1 strain, neither the PL65 nor the A5410 genome harbors the zeamine biosynthesis gene cluster, which plays a key role in bacterial virulence. The ANI and dDDH percentages between the two genomes were 94.47 and 57.00, respectively. In this study, we compared major virulence factors (plant cell wall-degrading extracellular enzymes and protease) produced by D. zeae strains and virulence ability on taro corms and pineapple. Both strains produced protease, pectate lyases and cellulases but no significant quantitative differences were observed (p>0.05) among the strains. All the strains produced symptoms on taro corms and pineapple leaves. Strain PL65 developed symptoms faster than the others. Our study highlights genetic constituents of pathogenicity determinants and genomic heterogeneity that will help understand the virulence mechanisms and aggressiveness of this plant pathogen.


title: Comparative genomics reveals factors associated with phenotypic expression of Wolbachia

abstract: Wolbachia is a widespread, vertically transmitted bacterial endosymbiont known for manipulating arthropod reproduction. Its most common form of reproductive manipulation is Cytoplasmic incompatibility (CI), observed when a modification in the male sperm leads to embryonic lethality unless a compatible rescue factor is present in the female egg. CI attracts scientific attention due to its implications for host speciation and in the use of Wolbachia for controlling vector-borne diseases. However, our understanding of CI is complicated by the complexity of the phenotype, whose expression depends on both symbiont and host factors. In the present study, we perform a comparative analysis of nine complete Wolbachia genomes with known CI properties in the same genetic host background, Drosophila simulans STC. We describe genetic differences between closely related strains and uncover evidence that phages and other mobile elements contribute to rapid evolution of both genomes and phenotypes of Wolbachia. Additionally, we identify both known and novel genes associated with the modification and rescue functions of CI. We combine our observations with published phenotypic information and discuss how variability in cif genes, novel CI-associated genes and Wolbachia titer might contribute to poorly understood aspects of CI such as strength and bidirectional incompatibility. We speculate that high titer CI strains could be better at invading new hosts already infected with a CI Wolbachia, due to a higher rescue potential, and suggest that titer might thus be a relevant parameter to consider for future strategies using CI Wolbachia in biological control.


title: Redox potential linked to water loss from proteins in evolution and development

abstract: Reactions involving water and oxygen are basic features of geological and biological processes. To understand how life interacts with its environment requires monitoring interactions with H2O and O2 at timescales relevant to not only organismal growth but also billions of years of geobiological evolution. Here, chemical transformations intrinsic to evolution and development were characterized by analyzing data from recent phylostratigraphic and proteomic studies. This two-stage analysis involves obtaining chemical metrics (carbon oxidation state and stoichiometric hydration state) from the elemental compositions of proteins followed by modeling the relative stabilities of target proteins against a proteomic background to infer thermodynamic parameters (oxygen fugacity, water activity, and virtual redox potential (Eh)). The main results of this study are a rise in carbon oxidation state of proteins spanning the time of the Great Oxidation Event, a rise in virtual redox potential that coincides with the likely emergence of aerobic metabolism, a rise in carbon oxidation state of proteins inferred from the transcriptome in late stages of Bacillus subtilis biofilm growth, and a drop in stoichiometric hydration state of the fruit fly developmental proteome at the same time as a drop in organismal water content. Stoichiometric hydration state also decreases for proteins with more recent gene ages and through stages of biofilm development, leading to predictions of higher hydration potentials at earlier time points. By building chemical representations of protein evolution and developmental proteomes, exciting new types of geobiochemical proxies can be developed.


title: A thermodynamic model for water activity and redox potential in evolution and development


title: TAK1: a potential therapeutic target for retinal neovascularization

abstract: Retinal neovascularization, or pathological angiogenesis in the retina, is a leading cause of blindness in developed countries. Transforming growth factor-{beta}-activated kinase 1 (TAK1) is a mitogen-activated protein kinase kinase kinase (MAPKKK) activated by TGF-{beta}1 and other pro-inflammatory cytokines. TAK1 is also a key mediator of inflammation, innate immune responses, apoptosis and tissue homeostasis and plays an important role in physiological angiogenesis. Its role in pathological angiogenesis, particularly in retinal neovascularization, remains unclear. We investigated the regulatory role of TAK1 in pathological angiogenesis in the retina. Transcriptome analysis of human retina featuring retinal neovascularization revealed enrichment of known TAK1-mediated signaling pathways. Selective inhibition of TAK1 activation by 5Z-7-oxozeaenol attenuated aberrant retinal angiogenesis in rats following oxygen-induced retinopathy. Transcriptome profiling revealed that TAK1 activation in human microvascular endothelial cells under TNF stimulation led to increase the gene expression related to cytokines and leukocyte-endothelial interaction, mainly through nuclear factor kappa B (NF{kappa}B) signaling pathways. These results reveal that inhibition of TAK1 signaling may have therapeutic value for the treatment of pathological angiogenesis in the retina.


title: TAK1 blockade as a therapy for retinal neovascularization