title: Polymerase theta-helicase promotes end joining by stripping single-stranded DNA-binding proteins and bridging DNA ends

abstract: Homologous recombination-deficient cancers rely on DNA polymerase Theta (Pol{theta})-Mediated End Joining (TMEJ), an alternative double-strand break repair pathway. Pol{theta} is the only vertebrate polymerase that encodes an N-terminal superfamily 2 (SF2) helicase domain, but the role of this helicase domain in TMEJ remains unclear. Using single-molecule imaging, we demonstrate that Pol{theta}-helicase (Pol{theta}-h) is a highly processive single-stranded DNA (ssDNA) motor protein that can efficiently strip Replication Protein A (RPA) from ssDNA. Pol{theta}-h also has a limited capacity for disassembling RAD51 filaments but is not processive on doublestranded DNA. Pol{theta}-h can bridge two non-complementary DNA strands in trans. PARylation of Pol{theta}-h by PARP-1 resolves these DNA bridges. We conclude that Pol{theta}-h removes RPA and RAD51 filaments and mediates bridging of DNA overhangs to aid in polymerization by the Pol{theta} polymerase domain.

category: biophysics


title: AR-V7 exhibits non-canonical mechanisms of nuclear import and chromatin engagement in Castrate-Resistant Prostate Cancer

abstract: Expression of the AR splice variant, AR-V7, in prostate cancer is correlated with poor patient survival and resistance to AR targeted therapies and taxanes. Currently, there is no specific inhibitor of AR-V7, while the molecular mechanisms regulating its biological function are not well elucidated. Here we report that AR-V7 has unique biological features that functionally differentiate it from canonical AR-fl or from the second most prevalent variant, AR-v567. First, AR-V7 exhibits fast nuclear import kinetics via a pathway distinct from the nuclear localization signal dependent importin- /{beta} pathway used by AR-fl and AR-v567. We also show that the dimerization box domain, known to mediate AR dimerization and transactivation, is required for AR-V7 nuclear import but not for AR-fl. Once in the nucleus, AR-V7 is transcriptionally active, yet exhibits unusually high intranuclear mobility and transient chromatin interactions, unlike the stable chromatin association of liganded AR-fl. The high intranuclear mobility of AR-V7 together with its high transcriptional output, suggest a Hit-and-Run mode of transcription. Our findings reveal unique mechanisms regulating AR-V7 activity, offering the opportunity to develop selective therapeutic interventions.

category: cancer biology


title: New fossils of Jurassic ophiurid brittle stars (Ophiuroidea; Ophiurida) provide evidence for early clade evolution in the deep sea

abstract: 1.Understanding of the evolutionary history of the ophiuroids, or brittle stars, is hampered by a patchy knowledge of the fossil record. Especially the stem members of the living clades are poorly known, resulting in blurry concepts of the early clade evolution and imprecise estimates of divergence ages. Here, we describe new ophiuroid fossil from the Lower Jurassic of France, Luxembourg, and Austria and introduce the new taxa Ophiogojira labadiei gen. et sp. nov. from lower Pliensbachian shallow sublittoral deposits, Ophiogojira andreui gen. et sp. nov. from lower Toarcian shallow sublittoral deposits, and Ophioduplantiera noctiluca gen. et sp. nov. from late Sinemurian to lower Pliensbachian bathyal deposits. A Bayesian morphological phylogenetic analysis shows that Ophiogojira holds a basal position within the order Ophiurida, whereas Ophioduplantiera has a more crownward position within the ophiurid family Ophiuridae. The position of Ophioduplantiera in the evolutionary tree suggests that family-level divergences within the Ophiurida must have occurred before the late Sinemurian, and that ancient slope environments played an important role in fostering early clade evolution.

category: paleontology


title: The Copper Responsive Transcription Factor SPL7 Represses Key Abscisic Acid Biosynthetic Genes to Balance Growth and Drought Tolerance

abstract:
Plants adapt to adverse environments by turning on defense against abiotic stresses, which is mainly orchestrated by the phytohormone abscisic acid (ABA). But how ABA homeostasis is modulated to balance growth and stress responses is still largely unknown. Here we report that prior treatment of Arabidopsis seedling with high copper retardates growth but enhances draught tolerance at later stages by modulating ABA accumulation. Subsequent genetic, physiological, transcriptomic, and molecular investigations revealed that the copper responsive transcription factor SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 7 (SPL7) is a strong regulator of ABA accumulation. We showed that SPL7 is destabilized by high copper and consistently suppresses genes encoding three key oxygenases in the ABA biosynthetic pathway of land plants via binding to the GTAC copper response motifs in their promoters. These results revealed a new mechanism whereby copper availability, inversely reflected by SPL7 abundance, modulates de novo ABA biosynthesis to balance growth and drought tolerance.

One-sentence summaryHigh copper availability represses SPL7, releasing its suppression on key ABA biosynthetic genes and leading to increased ABA accumulation that inhibits growth but enhances drought tolerance.

category: plant biology


title: Dual control of MAPK activities by AP2C1 and MKP1 MAPK phosphatases regulates defence responses in Arabidopsis

abstract:
Mitogen-activated protein kinase (MAPK) cascades transmit environmental signals and induce stress and defence responses in plants. These signalling cascades are negatively controlled by specific phosphatases of the type 2C Ser/Thr protein phosphatase (PP2C) and dual-specificity phosphatase (DSP) families that inactivate stress-induced MAPKs; however, the interplay between phosphatases of these different types has remained unknown. Our work reveals that different Arabidopsis MAPK phosphatases, the PP2C-type AP2C1 and the DSP-type MKP1, exhibit both specific and overlapping functions in plant stress responses. Each single mutant and ap2c1 mkp1 double mutant displayed enhanced wound-induced activation of MAPKs MPK3, MPK4, and MPK6, as well as induction of a set of transcription factors. Moreover, ap2c1 mkp1 double mutants show an autoimmune-like response, associated with elevated levels the stress hormones salicylic acid and ethylene, and of the phytoalexin camalexin. Interestingly, this phenotype is reduced in ap2c1 mkp1 mpk6 triple mutants, suggesting that the autoimmune-like response is due to MPK6 misregulation. We conclude that the evolutionarily distant MAPK phosphatases AP2C1 and MKP1 contribute crucially to the tight control of MPK6 activity, ensuring appropriately balanced stress signalling and suppression of autoimmune-like responses during plant growth and development.

HighlightDouble MAPK phosphatase mutant plants ap2c1 mkp1 exhibit constitutive, autoimmune-like stress responses, dependent on their substrate MAPK MPK6.


title: Large-scale characterisation of the pregnancy vaginal microbiome and sialidase activity in a low-risk Chinese population

abstract:
Vaginal microbiota-host interactions are linked to preterm birth (PTB), which continues to be the primary cause of global childhood mortality. Despite the majority of PTB occuring in Asia, studies of the pregnancy vaginal microbiota are largely limited to Northern American and European populations. Here, we characterised the vaginal microbiome of 2689 pregnant Chinese women using metataxonomics and in a subset (n=823), the relationship between vaginal microbiota composition, sialidase activity and leukocyte presence and pregnancy outcomes. Vaginal microbiota were most frequently dominated by Lactobacillus crispatus or L. iners, with the latter associated with vaginal leukocyte presence. Women with high sialidase activity were enriched for bacterial vaginosis-associated genera including Gardnerella, Atopobium and Prevotella. Vaginal microbiota composition, high sialidase activity and/or leukocyte presence was not associated with PTB risk suggesting underlying differences in the vaginal microbiota and/or host immune responses of Chinese women, possibly accounting for low PTB rates in this population.

ImportanceSpecific vaginal microorganisms or  vaginal microbiota, are associated with preterm birth, which is the primary cause of death in children under 5yrs of age worldwide. Despite most preterm births occuring in Asia, almost all studies of the pregnancy vaginal microbiota have been limited to Northern American and European women. Here, we studied the vaginal microbiota in a large cohort of 2689 pregnant Chinese women and showed that it was most frequently dominated by Lactobacillus crispatus or L. iners. The latter was associated with leukocyte infiltration of vaginal secretions. Women with high activity of the enzyme sialidase, were frequently colonised by species associated with the common condition bacterial vaginosis, including Gardnerella, Atopobium and Prevotella species. Vaginal microbiota, high sialidase activity and/or leukocyte presence was not associated with preterm birth risk indicating differences in the microbe-host immune responses of Chinese women, possibly explaining why preterm birth rates in this population are low.

category: systems biology


title: NLoed: A Python package for nonlinear optimal experimental design in systems biology

abstract:
MotivationModelling in systems and synthetic biology relies on accurate parameter estimates and predictions. Accurate model calibration relies, in turn, on data, and on how well-suited the available data is to a particular modelling task. Optimal experimental design (OED) techniques can be used to identify experiments and data collection procedures that will most efficiently contribute to a given modelling objective. However, implementation of OED is limited by currently available software tools that are not well-suited for the diversity of nonlinear models and non-normal data commonly encountered in biological research. Moreover, existing OED tools do not make use of the state-of-the-art numerical tools, resulting in inefficient computation.

ResultsHere we present the NLoed software package. NLoed is an open-source Python library providing convenient access to OED methods, with particular emphasis on experimental design for systems biology research. NLoed supports a wide variety of nonlinear, multi-input/output, and dynamic models, and facilitates modelling and design of experiments over a wide variety of data types. To support OED investigations, the NLoed package implements maximum likelihood fitting and diagnostic tools, providing a comprehensive modelling workflow. NLoed offers an accessible, modular, and flexible OED tool-set suited to the wide variety of experimental scenarios encountered in systems biology research. We demonstrate NLOEDs capabilities by applying it to experimental design for characterization of a bacterial optogenetic system.

AvailabilityNLoed is available via pip from the PyPi repository; https://pypi.org/project/nloed/. Source code, documentation and examples can be found on Github at https://github.com/ingallslab/NLoed.

Contactbingalls@uwaterloo.ca

Supplementary informationSupplementary materials are available online.


title: A representative reference for MRI-based human axon radius assessment using light microscopy

abstract: Non-invasive assessment of axon radii via MRI bears great potential for clinical and neuroscience research as it is a main determinant of the neuronal conduction velocity. However, there is a lack of representative histological reference data on the scale of the cross-section of MRI voxels for validating the MRI-visible, effective radius (reff). Because the current gold standard stems from neuroanatomical studies designed to estimate the frequency-weighted arithmetic mean radius (rarith) on small ensembles of axons, it is unsuited to estimate the tail-weighted reff. We propose CNN-based segmentation on high-resolution, large-scale light microscopy (lsLM) data to generate a representative reference for reff. In a human corpus callosum, we assessed estimation accuracy and bias of rarith and reff. Furthermore, we investigated whether mapping anatomy-related variation of rarith and reff is confounded by low-frequency variation of the image intensity, e.g., due to staining heterogeneity. Finally, we analyzed the potential error due to outstandingly large axons in reff. Compared to rarith, reff was estimated with higher accuracy (normalized-root-mean-square-error of reff: 7.2 %; rarith: 21.5 %) and lower bias (normalized-mean-bias-error of reff: -1.7 %; rarith: 16 %). While rarith was confounded by variation of the image intensity, variation of reff seemed anatomy-related. The largest axons contributed between 0.9 % and 3 % to reff. In conclusion, the proposed method accurately estimates reff at MRI voxel resolution across a human corpus callosum sample. Further investigations are required to assess generalization to brain areas with different axon radii ensembles.

category: neuroscience


title: Transsaccadic integration operates independently in different feature dimensions

abstract: Our knowledge about objects in our environment reflects an integration of current visual input with information from preceding gaze fixations. Such a mechanism may reduce uncertainty, but requires the visual system to determine which information obtained in different fixations should be combined or kept separate. To investigate the basis of this decision, we conducted three experiments. Participants viewed a stimulus in their peripheral vision, then made a saccade that shifted the object into the opposite hemifield. During the saccade, the object underwent changes of varying magnitude in two feature dimensions (Experiment 1: color and location, Experiments 2 and 3: color and orientation). Participants reported whether they detected any change and estimated one of the post-saccadic features. Integration of pre-saccadic with post-saccadic input was observed as a bias in estimates towards the pre-saccadic feature value. In all experiments, pre-saccadic bias weakened as the magnitude of the transsaccadic change in the estimated feature increased. Changes in the other feature, despite having a similar probability of detection, had no effect on integration. Results were quantitatively captured by an observer model where the decision whether to integrate information from sequential fixations is made independently for each feature and coupled to awareness of a feature change.


title: Scalable Bayesian GPFA with automatic relevance determination and discrete noise models

abstract: Latent variable models are ubiquitous in the exploratory analysis of neural population recordings, where they allow researchers to summarize the activity of large populations of neurons in lower dimensional  latent spaces. Existing methods can generally be categorized into (i) Bayesian methods that facilitate flexible incorporation of prior knowledge and uncertainty estimation, but which typically do not scale to large datasets; and (ii) highly parameterized methods without explicit priors that scale better but often struggle in the low-data regime. Here, we bridge this gap by developing a fully Bayesian yet scalable version of Gaussian process factor analysis (bGPFA), which models neural data as arising from a set of inferred latent processes with a prior that encourages smoothness over time. Additionally, bGPFA uses automatic relevance determination to infer the dimensionality of neural activity directly from the training data during optimization. To enable the analysis of continuous recordings without trial structure, we introduce a novel variational inference strategy that scales near-linearly in time and also allows for non-Gaussian noise models appropriate for electrophysiological recordings. We apply bGPFA to continuous recordings spanning 30 minutes with over 14 million data points from primate motor and somatosensory cortices during a self-paced reaching task. We show that neural activity progresses from an initial state at target onset to a reach-specific preparatory state well before movement onset. The distance between these initial and preparatory latent states is predictive of reaction times across reaches, suggesting that such preparatory dynamics have behavioral relevance despite the lack of externally imposed delay periods. Additionally, bGPFA discovers latent processes that evolve over slow timescales on the order of several seconds and contain complementary information about reaction time. These timescales are longer than those revealed by methods which focus on individual movement epochs and may reflect fluctuations in e.g. task engagement.


title: Functional Differences in the Neural Substrates of Auditory Cognition as a Consequence of Music Training

abstract:
Previous studies have demonstrated that musical deviants (syntactically irregular chords) elicit event related potentials/fields with negative polarity; specifically, the early right anterior negativity and the right anterior temporal negativity responses with peak latencies at ~200 ms and ~350 ms, respectively, post stimulus onset. Here, we investigated differences in the neural dynamics of the auditory perceptual system of individuals with music training compared to those with no music training. Magnetoencephalography was used to examine the neural response to a deviant sound when the auditory system was primed using stimulus entrainment to evoke an auditory gamma-band response between 31 Hz and 39 Hz, in 2 Hz steps. Participants responded to the harmonic relationship between the entrainment stimulus and the subsequent target stimulus. Gamma frequencies carry stimulus information; thus, the paradigm primed the auditory system with a known gamma frequency and evaluated any improvement in the brains response to a deviant stimulus. The entrainment stimuli did not elicit an early right anterior negativity response. Furthermore, the source location of the event-related field difference during the later right anterior temporal negativity response time-window varied depending on group and entrainment condition. In support of previous findings from research using this, and a functionally similar visual-priming paradigm, a 7 Hz phase modulation of gamma amplitude was found for non-musicians following 33 Hz stimulus entrainment. Overall, significant effects of gamma entrainment were found more frequently in the non-music brain. By contrast, musicians demonstrated a greater range of interactions with slower brain rhythms, indicative of increased top-down control.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=130 SRC="FIGDIR/small/446932v1_ufig1.gif" ALT="Figure 1">
View larger version (49K):
org.highwire.dtl.DTLVardef@1097a58org.highwire.dtl.DTLVardef@13f3c74org.highwire.dtl.DTLVardef@99fae1org.highwire.dtl.DTLVardef@125d351_HPS_FORMAT_FIGEXP  M_FIG C_FIG ConclusionO_LIMusicians auditory cognition relies more on top-down processes, while non-musicians rely on bottom-up processing, and therefore, their auditory cognition is facilitated by the entrained gamma-band response.
C_LIO_LIFirst neural evidence in support of faster reaction-time responses due to an interaction in phase of an entrained gamma-band response of 33 Hz and a slower endogenous theta rhythm (~7 Hz) - consistently reported in previous research.
C_LI


title: Verbally mediated cooperation is consistently associated with inter-brain synchrony in frontal and temporoparietal areas: A mini-review and meta-analysis

abstract: Single-brain neuroimaging studies have shown that human cooperation is associated with neural activity in frontal and temporoparietal regions. However, it remains unclear whether single-brain studies are informative about cooperation in real life, where people interact dynamically. Such dynamic interactions have become the focus of inter-brain studies. An advantageous technique in this regard is functional near-infrared spectroscopy (fNIRS) because it is less susceptible to movement artifacts than more conventional techniques like EEG or fMRI. We conducted a systematic review and the first quantitative meta-analysis of fNIRS hyperscanning of cooperation, based on thirteen studies with 890 participants. Overall, the meta-analysis revealed evidence of statistically significant inter-brain synchrony while people were cooperating, with large overall effect sizes in both frontal and temporoparietal areas. All thirteen studies observed significant inter-brain synchrony in the prefrontal cortex (PFC), suggesting that this region is particularly relevant for cooperative behavior. The consistency in these findings is unlikely to be due to task-related activations, given that the relevant studies used diverse cooperation tasks. Together, the present findings support the importance of inter-brain synchronization of frontal and temporoparietal regions in interpersonal cooperation. Moreover, the present article highlights the usefulness of meta-analyses as a tool for discerning patterns in inter-brain dynamics.


title: Cooperative behavior evokes inter-brain synchrony in the prefrontal and temporoparietal cortex: A systematic review and meta-analysis of fNIRS hyperscanning studies.


title: The timing of differentiation and potency of CD8 effector function is set by RNA binding proteins

abstract: CD8+ T cell differentiation into effector cells is initiated early after antigen encounter by signals from the T cell antigen receptor and costimulatory molecules. The molecular mechanisms that determine the timing and rate of differentiation however are not defined. Here we show that the RNA binding proteins (RBP) ZFP36 and ZFP36L1 limit the rate of differentiation of activated naive CD8+ T cells and the potency of the resulting cytotoxic lymphocytes. The RBP act in an early and short temporal window to enforce dependency on costimulation via CD28 for full T cell activation and effector differentiation by directly binding mRNA of NF-B, IRF8 and NOTCH1 transcription factors and cytokines, including IL2. Their absence in T cells, or the adoptive transfer of a small numbers of CD8+ T cells lacking the RBP, promotes resilience to influenza A virus infection without immunopathology. These findings highlight ZFP36 and ZFP36L1 as nodes for the integration of the early T cell activation signals determining the speed and quality of the CD8 response.

category: immunology


title: Modelling the effects of bioclimatic characteristics and climate change on the potential distribution of a monospecific species Colophospermum mopane (Benth.) Leonard in southern Africa.

abstract:
Global climate change is gradually changing species distribution and their patterns of diversity. Yet, factors that influence the local distribution and habitat preferences for southern African species remain largely unexplored. Here, we computed the suitable habitats in the southern African region for Colophospermum mopane (Benth.) using the maximum entropy (Maxent) modeling approach. We utilized one Global Circulation Model (GCM) and three Representative concentration pathways (RCPs) to determine the current and future distribution of C. mopane. The results showed that the distribution of C. mopane was mainly influenced by solar radiation, annual temperature range, and annual precipitation. According to the species response curves, this species preferred habitats with annual precipitation of 130-200 mm, an annual temperature range of 28{degrees} C, and elevations of about 500 m above sea level. The results highlight that the geographic range of C. mopane is likely to expand along the borders of Zambia and Zimbabwe in the future, particularly in the miombo plains. Conversely, suitable habitat areas reduce significantly in the eastern area of the southern African region, while the western areas expand. Overall, the appropriate habitat areas will likely decline in the 2050s under both RCPs and expand in the 2070s under the two scenarios. This knowledge is important for landscape planners and rangeland managers working to safeguard biodiversity from extinction.

HighlightsO_LIHigh reliability of models in habitat suitability modelling for C. mopane
C_LIO_LISolar radiation is the most significant variable for the current distribution of mopane.
C_LIO_LIClimate change is and will reduce the habitat suitability of our target species.
C_LI

category: ecology


title: Western Amazon was a center of Neotropical fish dispersal, as evidenced by the continental-wide time-stratified biogeographic analysis of the hyper-diverse Hypostomus catfish

abstract: The Amazon is probably the most diverse realm on Earth, and is considered to be the primary source of diversity and a center of dispersal for Neotropical terrestrial organisms. Yet, the assumption that the Amazon basin is a primordial place of fish species origination and dispersal into other drainages still need to be tested. We addressed this issue by inferring a time-stratified biogeographic history and reconstructing the ancestral habitat preference of Hypostomus, a continentally widespread and species-rich Neotropical genus. We found that Hypostomus emerged in the Western Amazon ([~]14.7 Ma), when the Western Amazon River was flowing northwards and disconnected from the Eastern Amazon. We show that dispersal events in the first half of Hypostomus evolution occurred from the Western Amazon into adjacent basins, initiating its Neotropical radiation. The ancestral preferred habitat consisted in small rivers with running waters, a predominant habitat in river headwaters. Because of strong niche conservatism in the early evolution of Hypostomus, we suggest that most of the out-of-Western-Amazon dispersal occurred via headwater captures. The radiation of Hypostomus was further promoted by major reconfigurations of river basins, which opened dispersal opportunities into new drainages. Diversification in habitat preference coincided with colonization of basins already occupied by congenerics, indicative of niche shifts triggered by inter-specific competition and species coexistence. By analyzing the evolutionary history of Hypostomus, we show that Western Amazon was the main center of fish dispersal in the Neotropical Region from Middle Miocene to the present, supporting the cradle hypothesis of fish origination and dispersal.

category: evolutionary biology


title: Identification of a motif in Trm732 required for 2'-O-methylation of the tRNA anticodon loop by Trm7

abstract: Posttranscriptional tRNA modifications are essential for proper gene expression, and defects in the enzymes that perform tRNA modifications are associated with numerous human disorders. Throughout eukaryotes, 2-O-methylation of residues 32 and 34 of the anticodon loop of tRNA is important for proper translation, and in humans, lack of these modifications results in non-syndromic X-linked intellectual disability. In yeast, the methyltransferase Trm7 forms a complex with Trm732 to 2-O-methylate tRNA residue 32 and with Trm734 to 2-O-methylate tRNA residue 34. Trm732 and Trm734 are required for the methylation activity of Trm7, but the role of these auxiliary proteins is not clear. Additionally, Trm732 and Trm734 homologs are implicated in biological processes not directly related to translation, suggesting that these proteins may have additional cellular functions. To identify critical amino acids in Trm732, we generated variants and tested their ability to function in yeast cells. We identified a conserved RRSAGLP motif in the conserved DUF2428 domain of Trm732 that is required for tRNA modification activity by both yeast Trm732 and its human homolog THADA. The identification of Trm732 variants that lack tRNA modification activity will help to determine if other biological functions ascribed to Trm732 and THADA are directly due to tRNA modification, or to secondary effects due to other functions of these proteins.

category: genetics


title: Genome-wide approaches delineate the additive, epistatic, and pleiotropic nature of variants controlling fatty acid composition in peanut (Arachis hypogaea L.).

abstract:
The fatty acid composition of seed oil is a major determinant of the flavor, shelf-life, and nutritional quality of peanuts. Major QTLs controlling high oil content, high oleic content, and low linoleic content have been characterized in several seed oil crop species. Here we employ genome-wide association approaches on a recently genotyped collection of 787 plant introduction accessions in the USDA peanut core collection, plus selected improved cultivars, to discover markers associated with the natural variation in fatty acid composition, and to explain the genetic control of fatty acid composition in seed oils.

Overall, 251 single nucleotide polymorphisms (SNPs) had significant trait associations with the measured fatty acid components. Twelve SNPs were associated with two or three different traits. Of these loci with apparent pleiotropic effects, 10 were associated with both oleic (C18:1) and linoleic acid (C18:2) content at different positions in the genome. In all 10 cases, the favorable allele had an opposite effect - increasing and lowering the concentration, respectively, of oleic and linoleic acid. The other traits with pleiotropic variant control were palmitic (C16:0), behenic (C22:0), lignoceric (C24:0), gadoleic (C20:1), total saturated, and total unsaturated fatty acid content. One hundred (100) of the significantly associated SNPs were located within 1000 kbp of 55 genes with fatty acid biosynthesis functional annotations. These genes encoded, among others: ACCase carboxyl transferase subunits, and several fatty acid synthase II enzymes.

With the exception of gadoleic (C20:1) and lignoceric (C24:0) acid content, which occur at relatively low abundance in cultivated peanut, all traits had significant SNP interactions exceeding a stringent Bonferroni threshold ( = 1%). We detected 7,682 pairwise SNP interactions affecting the relative abundance of fatty acid components in the seed oil. Of these, 627 SNP pairs had at least one SNP within 1000 kbp of a gene with fatty acid biosynthesis functional annotation. We evaluated 168 candidate genes underlying these SNP interactions. Functional enrichment and protein-to-protein interactions supported significant interactions (p- value < 1.0E-16) among the genes evaluated. These results show the complex nature of the biology and genes underlying the variation in seed oil fatty acid composition and contribute to an improved genotype-to-phenotype map for fatty acid variation in peanut seed oil.

Key phrasesSNP Genotyping, Genome-wide Association Study (GWAS), GWAS of interacting SNPs (GWASi), Pleiotropy, Seed fatty acid composition, Oleic-Linoleic acid ratio.

category: genomics


title: Single-Cell Sequencing Reveals Lineage-Specific Dynamic Genetic Regulation of Gene Expression During Human Cardiomyocyte Differentiation

abstract: Dynamic and temporally specific gene regulatory changes may underlie unexplained genetic associations with complex disease. During a dynamic process such as cellular differentiation, the overall cell type composition of a tissue (or an in vitro culture) and the gene regulatory profile of each cell can both experience significant changes over time. To identify these dynamic effects in high resolution, we collected single-cell RNA-sequencing data over a differentiation time course from induced pluripotent stem cells to cardiomyocytes, sampled at 7 unique time points in 19 human cell lines. We employed a flexible approach to map dynamic eQTLs whose effects vary significantly over the course of bifurcating differentiation trajectories, including many whose effects are specific to one of these two lineages. Our study design allowed us to distinguish true dynamic eQTLs affecting a specific cell lineage from expression changes driven by potentially non-genetic differences between cell lines such as cell composition. Additionally, we used the cell type profiles learned from single-cell data to deconvolve and re-analyze data from matched bulk RNA-seq samples. Using this approach, we were able to identify a large number of novel dynamic eQTLs in single cell data while also attributing dynamic effects in bulk to a particular lineage. Overall, we found that using single cell data to uncover dynamic eQTLs can provide new insight into the gene regulatory changes that occur among heterogeneous cell types during cardiomyocyte differentiation.


title: Structural basis for sterol sensing by Scap and Insig

abstract: The sterol regulatory element-binding protein (SREBP) pathway senses the cellular cholesterol level through sterol regulated association between Scap and Insig. Despite the recent structural determination of the transmembrane domains of human Scap and Insig-2 bound to 25-hydroxycholesterol (25HC), the structure and regulatory mechanism of the luminal domains of Scap by cholesterol remains elusive. Here, combining cryo-EM analysis and artificial intelligence-facilitated structural prediction, we report the structure of the human Scap/Insig-2 complex in the presence of digitonin instead of 25HC. Despite the lack of sequence similarity, the structure of the luminal domain Loop 1 and a co-folded segment in Loop 7 of Scap resembles that of the luminal/extracellular domain in NPC1 and related proteins. Comparison of the sterol-loaded structures of these proteins provides clues of the regulation of Loop 1/7 interaction by cholesterol. We also show that the structure of Scap(D428A), which suppresses SREBP activation under sterol depletion, is identical to WT when complexed with Insig-2, although the gain of function may also involve a later step in protein trafficking.

category: biochemistry


title: Molecular mechanism of quorum sensing inhibition in Streptococcus by the phage protein paratox

abstract: Streptococcus pyogenes, or Group A Streptococcus, is a Gram-positive bacterium that can be both a human commensal and pathogen. Central to this dichotomy are temperate bacteriophages that incorporate into the bacterial genome as a prophage. These genetic elements encode both the phage proteins as well as toxins harmful to the human host. One such conserved phage protein paratox (Prx) is always found encoded adjacent to the toxin genes and this linkage is preserved during transduction. Within Streptococcus pyogenes, Prx functions to inhibit the quorum-sensing ComRS receptor-signal pair that is the master regulator of natural competence, or the ability to uptake endogenous DNA. Specifically, Prx directly binds and inhibits the receptor ComR by unknown mechanism. To understand how Prx inhibits ComR at the molecular level we pursued an X-ray crystal structure of Prx bound to ComR. The structural data supported by solution X-ray scattering data demonstrate that Prx induces a conformational change in ComR to directly access the DNA binding domain. Furthermore, electromobility shift assays and competition binding assays reveal that Prx effectively uncouples the inter-domain conformational change that is required for activation of ComR by the signaling molecule XIP. Although to our knowledge the molecular mechanism of quorum-sensing inhibition by Prx is unique, it is analogous to the mechanism employed by the phage protein Aqs1 in Pseudomonas aeruginosa. Together, this demonstrates an example of convergent evolution between Gram-positive and Gram-negative phages to inhibit quorum-sensing, and highlights the versatility of small phage proteins.


title: TDFragMapper: a visualization tool for evaluating experimental parameters in top-down proteomics

abstract:
MotivationWe present a new software-tool allowing an easy visualization of fragment ions and thus a rapid evaluation of key experimental parameters on the sequence coverage obtained for the MS/MS analysis of intact proteins. Our tool can deal with multiple fragmentation methods.

ResultsWe demonstrate that TDFragMapper can rapidly highlight the experimental fragmentation parameters that are critical to the characterization of intact proteins of various size using top-down proteomics.

AvailabilityTDFragMapper, a demonstration video and user tutorial are freely available at https://msbio.pas-teur.fr/tdfragmapper, for academic use; all data are thus available from the ProteomeXchange consortium (identifier PXD024643).

Contactdiogobor@gmail.com or julia.chamot-rooke@pasteur.fr


title: Unique Features of Different Classes of G-Protein-Coupled Receptors Revealed from Sequence Coevolutionary and Structural Analysis

abstract:
G-protein-coupled receptors (GPCRs) are the largest family of human membrane proteins and represent the primary targets of about one third of currently marketed drugs. Despite the critical importance, experimental structures have been determined for only a limited portion of GPCRs and functional mechanisms of GPCRs remain poorly understood. Here, we have constructed novel sequence coevolutionary models of the A and B classes of GPCRs and compared them with residue contact frequency maps generated with available experimental structures. Significant portions of structural residue contacts were successfully detected in the sequence-based covariational models. "Exception" residue contacts predicted from sequence coevolutionary models but not available structures added missing links that were important for GPCR activation and allosteric modulation. Moreover, we identified distinct residue contacts involving different sets of functional motifs for GPCR activation, such as the Na+ pocket, CWxP, DRY, PIF and NPxxY motifs in the class A and the HETx and PxxG motifs in the class B. Finally, we systematically uncovered critical residue contacts tuned by allosteric modulation in the two classes of GPCRs, including those from the activation motifs and particularly the extracellular and intracellular loops in class A GPCRs. These findings provide a promising framework for rational design of ligands to regulate GPCR activation and allosteric modulation.

SignificanceG-protein-coupled receptors (GPCRs) play key roles in cellular signaling and serve as the primary targets of [~]1/3 of currently marketed drugs. In this work, we have presented the first analysis of both residue sequence coevolution and structural contact maps in different classes of GPCRs. We have inferred pathways for GPCR signal transduction that could not be determined from structural analysis alone. Distinct residue contacts have been identified in the signaling pathways of class A and B GPCRs. Our combined sequence coevolutionary and structural contact analysis has thus revealed important insights into the mechanism of GPCR signal transduction, which is expected to facilitate rational drug design of the GPCRs.

category: bioinformatics


title: Investigating the concept of accessibility for predicting novel RNA-RNA interactions

abstract:
State-of-the-art methods for predicting novel trans RNA-RNA interactions use the so-called accessibility as key concept. It estimates whether a region in a given RNA sequence is accessible for forming trans interactions, using a thermodynamic model which quantifies its secondary structure features. RNA-RNA interactions are then predicted by finding the minimum free energy base-pairing between the two transcripts, taking into account the accessibility as energy penalty.

We investigated the underlying assumptions of this approach using the two methods RNAO_SCPLOWPLEXC_SCPLOW and IO_SCPLOWNTAC_SCPLOWRNA on two datasets, containing sRNA-mRNA and snoRNA-rRNA interactions, respectively.

We find that (1) known trans RNA-RNA interactions frequently overlap regions containing RNA structure features, (2) the estimated accessibility reflects sRNA structures fairly well, but often disagrees with structure annotations of longer transcripts, (3) the prediction performance of RNA-RNA interaction prediction methods is independent of the quality of the estimated accessibility profiles, and (4) one important overall effect of accessibility profiles is to prevent the thermodynamic model from predicting too long interactions.

Based on our findings, we conclude that the accessibility concept to the minimum free energy approach to predicting novel RNA-RNA interactions has conceptual limitations and discuss potential ways of improving the field in the future.


title: Deep Learning Approach to Genomic Breakage Study from Primary Sequence

abstract:
Identifying the molecular mechanisms related to genomic breakage is an important goal of cancer mechanism studies. Among the diverse location of the breakpoints of structural variants, the fusion genes, which have the breakpoints in the gene bodies and typically identified from RNA-seq data, can provide a highlighted structural variant resource for studying the genomic breakages with expression and potential pathogenic impacts. In this study, we developed FusionAI which utilizes deep learning to predict gene fusion breakpoints based on primary sequences and let us identify fusion breakage code and genomic context. FusionAI leverages the known fusion breakpoints to provide a prediction model of the fusion genes from the primary genomic sequences via deep learning, thereby helping researchers a more accurate selection of fusion genes and better understand genomic breakage.

HighlightsO_LIFusionAI, a 9-layer deep neural network, predicts fusion gene breakpoints from a DNA sequence
C_LIO_LIFusonAI reduce the cost and effort for validating fusion genes by decreasing specificity
C_LIO_LIHigh feature importance scored regions were apart 100nt on average from the exon junction breakpoints
C_LIO_LIHigh feature importance scored regions overlapped with 44 different human genomic features
C_LIO_LITranscription factor fusion genes are targeted by the GC-rich motif TFs
C_LIO_LIFusionAI gives less scores to the non-disease derived breakpoints
C_LI


title: Protein Interaction Network-based Deep Learning Framework for Identifying Disease-Associated Human Proteins

abstract: Infectious diseases in humans appear to be one of the most primary public health issues. Identification of novel disease-associated proteins will furnish an efficient recognition of the novel therapeutic targets. Here, we develop a Graph Convolutional Network (GCN)-based model called PINDeL to identify the disease-associated host proteins by integrating the human Protein Locality Graph and its corresponding topological features. Because of the amalgamation of GCN with the protein interaction network, PINDeL achieves the highest accuracy of 83.45% while AUROC and AUPRC values are 0.90 and 0.88, respectively. With high accuracy, recall, F1-score, specificity, AUROC, and AUPRC, PINDeL outperforms other existing machine-learning and deep-learning techniques for disease gene/protein identification in humans. Application of PINDeL on an independent dataset of 24320 proteins, which are not used for training, validation, or testing purposes, predicts 6448 new disease-protein associations of which we verify 3196 disease-proteins through experimental evidence like disease ontology, Gene Ontology, and KEGG pathway enrichment analyses. Our investigation informs that experimentally-verified 748 proteins are indeed responsible for pathogen-host protein interactions of which 22 disease-proteins share their association with multiple diseases such as cancer, aging, chem-dependency, pharmacogenomics, normal variation, infection, and immune-related diseases. This unique Graph Convolution Network-based prediction model is of utmost use in large-scale disease-protein association prediction and hence, will provide crucial insights on disease pathogenesis and will further aid in developing novel therapeutics.


title: Reverse-Complement Equivariant Networks for DNA Sequences

abstract: As DNA sequencing technologies keep improving in scale and cost, there is a growing need to develop machine learning models to analyze DNA sequences, e.g., to decipher regulatory signals from DNA fragments bound by a particular protein of interest. As a double helix made of two complementary strands, a DNA fragment can be sequenced as two equivalent, so-called Reverse Complement (RC) sequences of nucleotides. To take into account this inherent symmetry of the data in machine learning models can facilitate learning. In this sense, several authors have recently proposed particular RC-equivariant convolutional neural networks (CNNs). However, it remains unknown whether other RC-equivariant architectures exist, which could potentially increase the set of basic models adapted to DNA sequences for practitioners. Here, we close this gap by characterizing the set of all linear RC-equivariant layers, and show in particular that new architectures exist beyond the ones already explored. We further discuss RC-equivariant pointwise nonlinearities adapted to different architectures, as well as RC-equivariant embeddings of k-mers as an alternative to one-hot encoding of nucleotides. We show experimentally that the new architectures can outperform existing ones.


title: A novel tubulin binding molecule drives differentiation of acute myeloid leukaemia cells

abstract: Acute Myeloid Leukaemia (AML) continues to have a poor prognosis, especially in the elderly. One reason for this is that many treatment regimens are not well tolerated by elderly patients. Much current focus is on the development of therapies that can target specific vulnerabilities of AML while having fewer toxic side effects. However, despite much recent progress in developing better drugs, many patients with AML still die within a year of diagnosis, partly due to the fact that it is difficult to identify therapeutic targets that are effective across multiple AML subtypes. One common factor across AML subtypes is the presence of a block in differentiation. Thus screening for compounds that can overcome this block in genetically diverse AML models should allow for the identification of agents that are not dependent on a specific mutation for their efficacy. Here, we used a phenotypic screen to identify novel compounds that stimulate differentiation in several AML cell lines. Lead compounds were shown to decrease tumour burden and to increase survival in vivo. Using multiple complementary target deconvolution approaches, these compounds were revealed to be anti-mitotic tubulin disruptors that cause differentiation by inducing a G2-M mitotic arrest. Together, these results reveal a novel function for tubulin disruptors in causing differentiation of AML cells.


title: Novel temporal and spatial patterns of metastatic colonization from rapid-autopsy tumor biopsies

abstract:
BackgroundMetastatic breast cancer is a deadly disease with a low 5-year survival rate. Tracking metastatic spread in living patients is difficult, and thus poorly understood.

ResultsVia rapid autopsy, we have collected 30 tumor samples over 3 timepoints and across 8 organs from a triple-negative metastatic breast cancer patient. The large number of sites sampled, together with deep whole genome sequencing and advanced computational analysis, allowed us to comprehensively reconstruct the tumors evolution at subclonal resolution. The most unique, previously not reported aspect of the tumors evolution we observed in this patient was the presence of "subclone incubators", i.e. already metastatic sites where substantial tumor evolution occurred before colonization of additional sites and organs by subclones that evolved at the incubator site. Overall, we identified four discrete waves of metastatic expansions, each of which resulted in a number of new, genetically similar metastasis sites that also enriched for particular organs (e.g. abdominal vs bone and brain). The lung played a critical role in facilitating metastatic spread in this patient: the lung was the first site of metastatic escape from the primary breast lesion; subclones at this site were the source of all four subsequent metastatic waves; and multiple sites in the lung acted as subclone incubators. Finally, functional annotation revealed that many known driver or metastasis-promoting tumor mutations in this patient were shared by some, but not all metastatic sites, highlighting the need for more comprehensive surveys of a patients metastases for effective clinical intervention.

ConclusionsOur analysis revealed the presence of substantial tumor evolution at metastatic incubator sites, with potentially important clinical implications. Our study demonstrated that sampling of a large number of metastatic sites affords unprecedented detail for studying metastatic evolution.


title: Alpha-to-beta- and gamma-band activity reflect predictive coding in affective visual processing

abstract: Processing of negative affective pictures typically leads to desynchronization of alpha-to-beta frequencies (ERD) and synchronization of gamma frequencies (ERS). Given that in predictive coding higher frequencies have been associated with prediction errors, while lower frequencies have been linked to expectations, we tested the hypothesis that alpha-to-beta ERD and gamma ERS induced by aversive pictures are associated with expectations and prediction errors, respectively. We recorded EEG while volunteers were involved in a probabilistically cued affective picture task using three different negative valences to produce expectations and prediction errors. Our data show that alpha-to-beta band activity was related to the expected valence of the stimulus as predicted by a cue. The absolute mismatch of the expected and actual valence, which denotes an absolute prediction error was related to gamma band activity. This demonstrates that top-down predictions and bottom-up prediction errors are represented in specific spectral patterns associated with affective picture processing.


title: Lack of p62 impairs glycogen aggregation and exacerbates pathology in a mouse model of myoclonic epilepsy of Lafora

abstract:
BackgroundLafora disease (LD) is a fatal childhood-onset dementia characterized by the extensive accumulation of glycogen aggregates--the so-called Lafora Bodies (LBs)--in several organs. The accumulation of LBs in the brain underlies the neurological phenotype of the disease. LBs are composed of abnormal glycogen and various associated proteins, including p62, an autophagy adaptor that participates in the aggregation and clearance of misfolded proteins.

MethodsTo study the role of p62 in the formation of LBs and its participation in the pathology of LD, we generated a mouse model of the disease (malinKO) lacking p62.

ResultsDeletion of p62 prevented LB accumulation in skeletal muscle and cardiac tissue. In the brain, the absence of p62 altered LB morphology and increased susceptibility to epilepsy.

ConclusionsThese results demonstrate that p62 participates in the formation of LBs and suggest that the sequestration of abnormal glycogen into LBs is a protective mechanism through which to reduce the deleterious consequences of its accumulation in the brain.


title: Lung epithelial signaling mediates early vaccine-induced CD4+ T cell activation and Mtb control

abstract:
Tuberculosis (TB) is one of the leading causes of death due to a single infectious agent. The development of a TB vaccine that induces durable and effective immunity to Mycobacterium tuberculosis (Mtb) infection is urgently needed. Early and superior Mtb control can be induced in M. bovis Bacillus Calmette-Guerin (BCG) vaccinated hosts when the innate immune response is targeted to generate effective vaccine-induced immunity. In the present study, we show that innate activation of DCs is critical for mucosal localization of clonally activated vaccine-induced CD4+ T cells in the lung, and superior early Mtb control. In addition, our study reveals that Th1/Th17 cytokine axis play an important role in superior vaccine induced immunity. Our studies also show that activation of nuclear factor kappa-light-chain-enhancer of activated B cells (NF{kappa}{beta}) pathway in lung epithelial cells is critical for the mucosal localization of activated vaccine-induced CD4+ T cells for rapid Mtb control. Thus, our study provides novel insights into the immune mechanisms that can overcome TB vaccine bottlenecks and provide early rapid Mtb control.

ImportanceTuberculosis is a leading cause of death due to single infectious agent accounting 1.4 million deaths each year. The only licensed vaccine BCG is not effective due to variable efficacy. In our study, we determined the early immune events necessary for achieving complete protection in BCG vaccinated host. Our study reveals that innate activation of DCs can mediate superior and early Mtb control in BCG vaccinated host through lung epithelial cell signaling and localization of clonal activated, Mtb antigen specific, cytokine producing CD4+ T cells within the lung parenchyma and airways. Thus, our study provides novel insights into the immune mechanisms that can overcome TB vaccine bottlenecks and provide early rapid Mtb control.


title: Single cell clonal analysis identifies an AID-dependent pathway of plasma cell differentiation

abstract: Germinal centers (GC) are microstructures where B cells that have been activated by antigen can improve the affinity of their B cell receptors and differentiate into memory B cells (MBCs) or antibody secreting plasma cells. Activation Induced Deaminase (AID) initiates antibody diversification in GCs by somatic hypermutation and class switch recombination. Here we have addressed the role of AID in the terminal differentiation of GC B cells by combining single cell transcriptome and immunoglobulin clonal analysis in a mouse model that traces AID-experienced cells. We identified 8 transcriptional clusters that include dark zone and light zone GC subsets, plasmablasts/plasma cells (PB), 4 subsets of MBCs and a novel prePB subset, which shares the strongest clonal relationships with PBs. Mice lacking AID have various alterations in the size and expression profiles of these transcriptional clusters. We find that AID deficiency leads to a reduced proportion of prePB cells and severely impairs transitions between the prePB and the PB subsets. Thus, AID shapes the differentiation fate of GC B cells by enabling PB generation from a prePB state.


title: Indexable signal amplification for multiparametric imaging

abstract: Multiparametric imaging allows researchers to measure the expression of many biomarkers simultaneously, allowing detailed characterization of cell microenvironments. One such technique, CODEX, allows fluorescence imaging of >30 proteins in a single tissue section. In the commercial CODEX system, primary antibodies are conjugated to DNA barcodes. This modification can result in antibody dysfunction, and development of a custom antibody panel can be very costly and time consuming as trial and error of modified antibodies proceeds. To address these challenges, we developed novel tyramide-conjugated DNA barcodes that can be used with primary antibodies via peroxidase-conjugated secondary antibodies. This approach results in signal amplification and imaging without the need to conjugate primary antibodies. When combined with commercially available barcode-conjugated primary antibodies, we can very quickly develop working antibody panels. We also present methods to perform antibody staining using a commercially available automated tissue stainer and in situ hybridization imaging on a CODEX platform.


title: Fluctuating environment can negate cheater success due to speed-agility trade-off

abstract: Stability of microbial cooperation through common goods is susceptible to cheating. Evidence suggests that cheating plays a less prominent role in many natural systems than hitherto predicted by models of eco-evolutionary dynamics and evolutionary game theory. While several cheater negating factors such as spatial segregation have been identified, most consider single-nutrient regimes. Here we propose a cheater-suppressing mechanism based on previous experimental observations regarding the biochemical trade-off between growth speed and delay in switching to alternative nutrients. As changing the nutrient source requires redistribution of enzymatic resources to different metabolic pathways, the advantage in speed is offset by lower agility due to longer time required for resource re-allocation. Using an in silico model system of sucrose utilisation by Saccharomyces cerevisiae, we find that a tradeoff between growth rate and diauxic lag duration can supress cheaters under fluctuating nutrient availability and thereby stabilise cooperation. The resulting temporal dynamics constrain cheaters despite their competitive benefit for the growth on the primary nutrient via avoided public goods synthesis costs. We further show that this speed-agility trade-off can function in synergy with spatial segregation to avoid the collapse of the community due to the cheaters. Taken together, the growth-agility trade-off may contribute to cheater suppression in microbial ecosystems experiencing fluctuating environments, such as plant root microbiota and gut microbiota.


title: A large chromosomal inversion shapes gene expression in seaweed flies (Coelopa frigida)

abstract: Inversions often underlie complex adaptive traits, but the genic targets inside them are largely unknown. Gene expression profiling provides a powerful way to link inversions with their phenotypic consequences. We examined the effects of the Cf-Inv(1) inversion in the seaweed fly Coelopa frigida on gene expression variation across sexes and life stages. Our analyses revealed that Cf-Inv(1) shapes global expression patterns but the extent of this effect is variable with much stronger effects in adults than larvae. Furthermore, within adults, both common as well as sex specific patterns were found. The vast majority of these differentially expressed genes mapped to Cf-Inv(1). However, genes that were differentially expressed in a single context (i.e. in males, females or larvae) were more likely to be located outside of Cf-Inv(1). By combining our findings with genomic scans for environmentally associated SNPs, we were able to pinpoint candidate variants in the inversion that may underlie mechanistic pathways that determine phenotypes. Together the results in this study, combined with previous findings, support the notion that the polymorphic Cf-Inv(1) inversion in this species is a major factor shaping both coding and regulatory variation resulting in highly complex adaptive effects.


title: A fast, general synteny detection engine

abstract: The increasingly widespread availability of genomic data has created a growing need for fast, sensitive and scalable comparative analysis methods. A key aspect of comparative genomic analysis is the study of synteny, co-localized gene clusters shared among genomes due to descent from common ancestors. Synteny can provide unique insight into the origin, function, and evolution of genome architectures, but methods to identify syntenic patterns in genomic datasets are often inflexible and slow, and use diverse definitions of what counts as likely synteny. Moreover, the reliable identification of putatively syntenic regions (i.e., whether they are truly indicative of homology) with different lengths and signal to noise ratios can be difficult to quantify. Here, we present Mology, a fast, flexible, alignment-free, nonparametric method to detect regions of syntenic elements among genomes or other datasets. The core algorithm operates on consecutive, rank-ordered elements, which could be genes, operons, motifs, sequence fragments, or any other orderable element. It is agnostic to the physical distance between distinct elements and also to directionality and order within syntenic regions, although such considerations can be addressed post hoc. We describe the underlying statistical theory behind our analysis method, and employ a Monte Carlo approach to estimate the false positive rate and positive predictive values for putative syntenic regions. We also evaluate how varying amounts of noise affect recovery of true syntenic regions among Saccharomycetaceae yeast genomes with up to ~100 million years of divergence. We discuss different strategies for recursive application of our method on syntenic regions with sparser signal than considered here, as well as the general applicability of the core algorithm.


title: PHARP: A pig haplotype reference panel for genotype imputation

abstract: Pigs not only function as a major meat source worldwide but also are commonly used as an animal model for studying human complex traits. A large haplotype reference panel has been used to facilitate efficient phasing and imputation of relatively sparse genome-wide microarray chips and low-coverage sequencing data. Using the imputed genotypes in the downstream analysis, such as GWASs, TWASs, eQTL mapping and genomic prediction (GS), is beneficial for obtaining novel findings. However, currently, there is still a lack of publicly available and high-quality pig reference panels with large sample sizes and high diversity, which greatly limits the application of genotype imputation in pigs. In response, we built the pig Haplotype Reference Panel (PHARP) database. PHARP provides a reference panel of 2,012 pig haplotypes at 34 million SNPs constructed using whole-genome sequence data from more than 49 studies of 71 pig breeds. It also provides Web-based analytical tools that allow researchers to carry out phasing and imputation consistently and efficiently. PHARP is freely accessible at http://alphaindex.zju.edu.cn/PHARP/index.php. We demonstrate its applicability for pig commercial 50K SNP arrays, by accurately imputing 2.6 billion genotypes at a concordance rate value of 0.971 in 81 Large White pigs (~ 17x sequencing coverage). We also applied our reference panel to impute the low-density SNP chip into the high-density data for three GWASs and found novel significantly associated SNPs that might be casual variants.


title: How do packet losses affect measures of averaged neural signals?

abstract: Recent advances in implanted device development have enabled chronic streaming of neural data to external devices allowing for long timescale, naturalistic recordings. However, characteristic data losses occur during wireless transmission. Estimates for the duration of these losses are typically uncertain reducing signal quality and impeding analyses. To characterize the effect of these losses on recovery of averaged neural signals, we simulated neural time series data for a typical event-related potential (ERP) experiment. We investigated how the signal duration and the degree of timing uncertainty affected the offset of the ERP, its duration in time, its amplitude, and the ability to resolve small differences corresponding to different task conditions. Simulations showed that long timescale signals were generally robust to the effects of packet losses apart from timing offsets while short timescale signals were significantly delocalized and attenuated. These results provide clarity on the types of signals that can be resolved using these datasets and provide clarity on the restrictions imposed by data losses on typical analyses.

category: bioengineering


title: Piezo1 ion channels inherently function as independent mechanotransducers

abstract: Piezo1 is a mechanically activated ion channel involved in sensing forces in various cell types and tissues. Cryo-electron microscopy has revealed that the Piezo1 structure is bowl-shaped and capable of inducing membrane curvature via its extended footprint, which indirectly suggest that Piezo1 ion channels may bias each others spatial distribution and interact functionally. Here, we use cell-attached patch-clamp electrophysiology and pressure-clamp stimulation to functionally examine large numbers of membrane patches from cells expressing Piezo1 endogenously at low levels and cells overexpressing Piezo1 at high levels. Our data, together with stochastic simulations of Piezo1 spatial distributions, show that both at endogenous densities (1-2 channels/m2), and at non-physiological densities (10-20 channels/m2) predicted to cause substantial footprint overlap, Piezo1 density has no substantial effect on its pressure sensitivity or open probability in the nominal absence of membrane tension. The results suggest that Piezo channels, at densities likely to be physiologically relevant, inherently behave as independent mechanotransducers. We propose this property is essential for cells to transduce forces homogeneously across the entire cell membrane.


title: Mechanistic studies of non-canonical amino acid mutagenesis

abstract: Over the past decade, harnessing the cellular protein synthesis machinery to incorporate non-canonical amino acids (ncAAs) into tailor-made peptides has significantly advanced many aspects of molecular science. More recently, groundbreaking progress in our ability to engineer this machinery for improved ncAA incorporation has led to significant enhancements of this powerful tool for biology and chemistry. By revealing the molecular basis for the poor or improved incorporation of ncAAs, mechanistic studies of ncAA incorporation by the protein synthesis machinery have tremendous potential for informing and directing such engineering efforts. In this chapter, we describe a set of complementary biochemical and single-molecule fluorescence assays that we have adapted for mechanistic studies of ncAA incorporation. Collectively, these assays provide data that can guide engineering of the protein synthesis machinery to expand the range of ncAAs that can be incorporated into peptides and increase the efficiency with which they can be incorporated, thereby enabling the full potential of ncAA mutagenesis technology to be realized.


title: Enteric coronavirus infection and treatment modeled with an immunocompetent human intestine-on-a-chip

abstract: Many patients infected with coronaviruses, such as SARS-CoV-2 and NL63 that use ACE2 receptors to infect cells, exhibit gastrointestinal symptoms and viral proteins are found in the human gastrointestinal tract, yet little is known about the inflammatory and pathological effects of coronavirus infection on the human intestine. Here, we used a human intestine-on-a-chip (Intestine Chip) microfluidic culture device lined by patient organoid-derived intestinal epithelium interfaced with human vascular endothelium to study host cellular and inflammatory responses to infection with NL63 coronavirus. These organoid-derived intestinal epithelial cells dramatically increased their ACE2 protein levels when cultured under flow in the presence of peristalsis-like mechanical deformations in the Intestine Chips compared to when cultured statically as organoids or in Transwell inserts. Infection of the intestinal epithelium with NL63 on-chip led to inflammation of the endothelium as demonstrated by loss of barrier function, increased cytokine production, and recruitment of circulating peripheral blood mononuclear cells (PMBCs). Treatment of NL63 infected chips with the approved protease inhibitor drug, nafamostat, inhibited viral entry and resulted in a reduction in both viral load and cytokine secretion, whereas remdesivir, one of the few drugs approved for COVID19 patients, was not found to be effective and it also was toxic to the endothelium. This model of intestinal infection was also used to test the effects of other drugs that have been proposed for potential repurposing against SARS-CoV-2. Taken together, these data suggest that the human Intestine Chip might be useful as a human preclinical model for studying coronavirus related pathology as well as for testing of potential anti-viral or anti-inflammatory therapeutics.


title: TUSX: an accessible toolbox for transcranial ultrasound simulation

abstract: AO_SCPLOWBSTRACTC_SCPLOWNormally, the complicated nature of acoustic simulation makes it infeasible for most research groups doing individual transcranial ultrasound studies, hindering interpretation of results and complicating the determination of safety limits. We present here an open-source MATLAB toolbox to perform acoustic simulations using subject-specific medical images for transcranial ultrasound experiments. This toolbox, Transcranial Ultrasound Simulation Toolbox (TUSX), consists of an integrated processing pipeline that takes in structural MR or CT images, processes them for accurate simulation, and runs the simulations using k-Wave, an existing open-source acoustics package. We describe here the processing TUSX performs, along with its reasoning. We also validate its output using real-world pressure measurements in a water tank.


title: Pathogenic neutrophilia drives acute respiratory distress syndrome in severe COVID-19 patients

abstract:
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the ensuing COVID-19 pandemic have caused [~]40 million cases and over 648,000 deaths in the United States alone. Troubling disparities in COVID-19-associated mortality emerged early, with nearly 70% of deaths confined to Black/African-American (AA) patients in some areas, yet targeted studies within this demographic are scant. Multi-omics single-cell analyses of immune profiles from airways and matching blood samples of Black/AA patients revealed low viral load, yet pronounced and persistent pulmonary neutrophilia with advanced features of cytokine release syndrome and acute respiratory distress syndrome (ARDS), including exacerbated production of IL-8, IL-1{beta}, IL-6, and CCL3/4 along with elevated levels of neutrophil elastase and myeloperoxidase. Circulating S100A12+/IFITM2+ mature neutrophils are recruited via the IL-8/CXCR2 axis, which emerges as a potential therapeutic target to reduce pathogenic neutrophilia and constrain ARDS in severe COVID-19.

Graphical AbstractThe lung pathology due to severe COVID-19 is marked by a perpetual pathogenic neutrophilia, leading to acute respiratory distress syndrome (ARDS) even in the absence of viral burden. Circulating mature neutrophils are recruited to the airways via IL-8 (CXCL8)/CXCR2 chemotaxis. Recently migrated neutrophils further differentiate into a transcriptionally active and hyperinflammatory state, with an exacerbated expression of IL-8 (CXCL8), IL-1{beta} (IL1B), CCL3, CCL4, neutrophil elastase (NE), and myeloperoxidase (MPO) activity. Airway neutrophils and recruited inflammatory monocytes further increase their production of IL-8 (CXCL8), perpetuating lung neutrophilia in a feedforward loop. MdCs and T cells produce IL-1{beta} and TNF, driving neutrophils reprogramming and survival.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=142 SRC="FIGDIR/small/446468v2_ufig1.gif" ALT="Figure 1">
View larger version (43K):
org.highwire.dtl.DTLVardef@81fd3aorg.highwire.dtl.DTLVardef@181e63org.highwire.dtl.DTLVardef@172fedcorg.highwire.dtl.DTLVardef@ba55a7_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: CNS antigen-specific control of early B-lineage development in the meninges

abstract: The random V(D)J recombination of immunoglobulins (Ig) loci often creates autoreactive B cell progenitors expressing self-recognized B-cell receptors (BCRs)1, which are eliminated or inactivated through an autoantigen-dependent central tolerance checkpoint to prevent autoimmune reactions2,3, a process thought to be restricted to the bone marrow (BM) in the adult mammals4. Here we report that early developing B cells are also present in the meninges of mice at all ages. Single cell RNA-sequencing (scRNA-seq) analysis revealed a consecutive trajectory of meningeal developing B cells in mice and non-human primates (NHPs). Parabiosis together with lineage tracing of hematopoietic stem cells (HSCs) showed that meningeal developing B cells are continuously replenished from the HSC-derived progenitors via a circulationindependent route. Importantly, autoreactive immature B cells which recognize myelin oligodendrocyte glycoprotein (MOG)5, a central nervous system (CNS)-specific antigen, are eliminated from the meninges but not BM. Furthermore, genetic deletion of MOG restored the self-reactive B cells in the meninges. Thus, we propose that meninges function as a unique reservoir for B cell development, allowing in situ negative selection of CNS-antigen-autoreactive B cells to ensure a local non-self-reactive immune repertoire.


title: The mutational landscape of SARS-CoV-2 variants diversifies T cell targets in an HLA supertype-dependent manner

abstract: The rapid, global dispersion of SARS-CoV-2 since its initial identification in December 2019 has led to the emergence of a diverse range of variants. The initial concerns regarding the virus were quickly compounded with concerns relating to the impact of its mutated forms on viral infectivity, pathogenicity and immunogenicity. To address the latter, we seek to understand how the mutational landscape of SARS-CoV-2 has shaped HLA-restricted T cell immunity at the population level during the first year of the pandemic, before mass vaccination. We analyzed a total of 330,246 high quality SARS-CoV-2 genome assemblies sampled across 143 countries and all major continents. Strikingly, we found that specific mutational patterns in SARS-CoV-2 diversify T cell epitopes in an HLA supertype-dependent manner. In fact, we observed that proline residues are preferentially removed from the proteome of prevalent mutants, leading to a predicted global loss of SARS-CoV-2 T cell epitopes in individuals expressing HLA-B alleles of the B7 supertype family. In addition, we show that this predicted global loss of epitopes is largely driven by a dominant C-to-U mutation type at the RNA level. These results indicate that B7 supertype-associated epitopes, including the most immunodominant ones, were more likely to escape CD8+ T cell immunosurveillance during the first year of the pandemic. Together, our study lays the foundation to help understand how SARS-CoV-2 mutants shape the repertoire of T cell targets and T cell immunity across human populations. The proposed theoretical framework has implications in viral evolution, disease severity, vaccine resistance and herd immunity.


title: CD4+ T cell lymphopenia and dysfunction in severe COVID-19 disease is autocrine TNF-α/TNFRI-dependent

abstract:
Lymphopenia is common in severe COVID-19 disease, yet the mechanisms are poorly understood. In 148 patients with severe COVID-19, we found lymphopenia was associated with worse survival. CD4+ lymphopenia predominated, with lower CD4+/CD8+ ratios in severe COVID-19 compared to recovered, mild disease (p<0.0001). In severe disease, immunodominant CD4+ T cell responses to Spike-1(S1) produced increased in vitro TNF-, but impaired proliferation and increased susceptibility to activation-induced cell death (AICD). CD4+TNF-+ T cell responses inversely correlated with absolute CD4+ counts from severe COVID-19 patients (n=76; R=-0.744, P<0.0001). TNF- blockade including infliximab or anti-TNFRI antibodies strikingly rescued S1-specific CD4+ proliferation and abrogated S1-AICD in severe COVID-19 patients (P<0.001). Single-cell RNAseq demonstrated downregulation of Type-1 cytokines and NF{kappa}B signaling in S1-stimulated CD4+ cells with infliximab treatment. Lung CD4+ T cells in severe COVID-19 were reduced and produced higher TNF- versus PBMC. Together, our findings show COVID-19-associated CD4+ lymphopenia and dysfunction is autocrine TNF-/TNFRI-dependent and therapies targeting TNF- may be beneficial in severe COVID-19.

One Sentence SummaryAutocrine TNF-/TNFRI regulates CD4+ T cell lymphopenia and dysfunction in severe COVID-19 disease.


title: KIF18B is a cell-type specific regulator of spindle orientation in the epidermis

abstract: Proper spindle orientation is required for asymmetric cell division and the establishment of complex tissue architecture. In the developing epidermis, spindle orientation requires a conserved cortical protein complex of LGN/NuMA/dynein-dynactin. However, how microtubule dynamics are regulated to interact with this machinery and properly position the mitotic spindle is not fully understood. Furthermore, our understanding of the processes that link spindle orientation during asymmetric cell division to cell fate specification in distinct tissue contexts remains incomplete. We report a role for the microtubule catastrophe factor KIF18B in regulating microtubule dynamics to promote spindle orientation in keratinocytes. During mitosis, KIF18B accumulates at the cell cortex, colocalizing with the conserved spindle orientation machinery. In vivo we find that KIF18B is required for oriented cell divisions within the hair placode, the first stage of hair follicle morphogenesis, but is not essential in the interfollicular epidermis. Disrupting spindle orientation in the placode, using mutations in either KIF18B or NuMA, results in aberrant cell fate marker expression of hair follicle progenitor cells. These data functionally link spindle orientation to cell fate decisions during hair follicle morphogenesis. Taken together, our data demonstrate a role for regulated microtubule dynamics in spindle orientation in epidermal cells. This work also highlights the importance of spindle orientation during asymmetric cell division to dictate cell fate specification.

category: cell biology


title: Volume growth in animal cells is cell cycle dependent and shows additive fluctuations

abstract: The way proliferating animal cells coordinate the growth of their mass, volume, and other relevant size parameters is a long-standing question in biology. Studies focusing on cell mass have identified patterns of mass growth as a function of time and cell cycle phase, but little is known about volume growth. To address this question, we improved our fluorescence exclusion method of volume measurement (FXm) and obtained 1700 single-cell volume growth trajectories of HeLa cells. We find that, during most of the cell cycle, volume growth is close to exponential and proceeds at a higher rate in S-G2 than in G1. Comparing the data with a mathematical model, we establish that the cell-to-cell variability in volume growth arises from constant-amplitude fluctuations in volume steps, rather than fluctuations of the underlying specific growth rate. We hypothesize that such "additive noise" could emerge from the processes that regulate volume adaptation to biophysical cues, such as tension or osmotic pressure.


title: Nitrogen fixation in mesoscale eddies of the North Pacific Subtropical Gyre: patterns and mechanisms

abstract: 1.Mesoscale eddies have been shown to support elevated dinitrogen (N2) fixation rates (NFRs) and abundances of N2-fixing microorganisms (diazotrophs), but the mechanisms underlying these observations are not well understood. We explored relationships among NFRs and cyanobacterial diazotroph abundances in eddy pairs of opposite polarity sampled in the North Pacific Subtropical Gyre and compared our observations to seasonal trends from the Hawaii Ocean Time-series (HOT) program. Consistent with previous reports, we found that NFRs were anomalously high for this region (up to 3.7-fold above previous monthly HOT observations) in the centers of the sampled anticyclones, coinciding with elevated abundances of Crocosphaera in the summertime. We then coupled our field-based observations, together with transcriptomic analyses of nutrient stress marker genes and ecological models, to evaluate potential mechanisms controlling diazotroph abundance and activity at the mesoscale. Specifically, we evaluated the role of biological (via estimates of growth and grazing rates) and physical controls on populations of Crocosphaera, Trichodesmium, and diatom symbionts. Our results suggest that increased Crocosphaera abundances associated with summertime anticyclones resulted from the alleviation of phosphate limitation, allowing cells to grow at rates exceeding grazing losses. In contrast, distributions of larger, buoyant taxa (Trichodesmium and diatom symbionts) appeared less affected by eddy-driven biological controls. Instead, they appeared driven by physical dynamics along frontal boundaries that separate cyclonic and anticyclonic eddies. Together, the interplay of eddy-specific changes in bottom-up control, top-down control, and the physical accumulation of cells likely explains the elevated diazotroph abundances and NFRs associated with anticyclones and eddy fronts.


title: What should a poor mother do? Influence of host plant quality on oviposition strategy and behavior in a polyphagous moth

abstract: To maximise fitness, individuals may apply different reproductive strategies. Such strategies could be phenotypically plastic and vary depending on the environment. For example, when resources are limited females often face a trade-off between investing in offspring quantity and quality, and how she balances this trade-off may depend on the environment. For phytophagous insects, and especially generalist insects, variation in host plant quality could have large effects on mating, reproduction and offspring performance. Here, we study if the polyphagous moth Spodoptera littoralis, which selects host plants through experience-based preference induction, also has a flexible allocation between egg weight and egg number as well as in temporal egg-laying behavior depending on larval host plant species. We found that S. littoralis has a canalized egg size and that an increased reproductive investment is made in egg quantity rather than egg quality. This increased investment depends on larval host plant species, probably reflecting parental condition. The constant egg weight may be due to physiological limitations or to limited possibilities to increase fitness through larger offspring size. We furthermore found that differences in onset of egg-laying is mainly due to differences in mating propensity between individuals raised on different host plant species. Thus, females do not seem to make a strategic reproductive investment in challenging environments. Instead, the low-quality host plant induces less and later reproduction, which could have consequences for population dynamics in the field.


title: Sorting droplets into many outlets

abstract: Droplet microfluidics is a commercially successful technology, widely used in single cell sequencing and droplet PCR. Combining droplet making with droplet sorting has also been demonstrated, but so far found limited use, partly due to difficulties in scaling manufacture with injection molded plastics. We introduce a droplet sorting system with several new elements, including: 1) an electrode design combining metallic and ionic liquid parts, 2) a modular, multi-sorting fluidic design with features for keeping inter-droplet distances constant, 3) using timing parameters calculated from fluorescence or scatter signal triggers to precisely actuate dozens of sorting electrodes, 4) droplet collection techniques, including ability to collect a single droplet, and 5) a new emulsion breaking method to collect aqueous samples for downstream analysis. We use these technologies to build a fluorescence based cell sorter that sorts with high purity. We also show that these microfluidic designs can be translated into injection molded thermoplastic, suitable for industrial production. Finally, we tally the advantages and limitations of these devices.


title: ProMoCell and ProModb: web services for analyzing interaction-based functionally localized protein modules in a cell

abstract: The modular organization of a cell which can be determined by its interaction network allows us to understand a mesh of cooperation among the functional modules. Therefore, cellular level identification of functional modules aids in understanding the functional and structural characteristics of the biological network of a cell and also assists in determining or comprehending the evolutionary signal. We develop ProMoCell that performs real-time web scraping for generating clusters of the cellular level functional units of an organism. ProMoCell constructs the Protein Locality Graphs and clusters the cellular level functional units of an organism by utilizing experimentally verified data from various online sources. Also, we develop ProModb, a database service that houses precomputed whole-cell protein-protein interaction network-based functional modules of an organism using ProMoCell. Our web service is entirely synchronized with the KEGG pathway database and allows users to generate spatially localized protein modules for any organism belonging to the KEGG genome using its real-time web scraping characteristics. Hence, the server will host as many organisms as is maintained by the KEGG database. Our web services provide the users a comprehensive and integrated tool for an efficient browsing and extraction of the spatial locality-based protein locality graph and the functional modules constructed by gathering experimental data from several interaction databases and pathway maps. We believe that our web services will be beneficial in pharmacological research, where a novel research domain called modular pharmacology has initiated the study on the diagnosis, prevention, and treatment of deadly diseases using functional modules.


title: Deciphering biological evolution exploiting the topology of Protein Locality Graph

abstract: The conventional sequence comparison-based evolutionary studies ignore other evolutionary constraints like interaction among proteins, functions of proteins and genes etc. A lot of speculations exist in literature regarding the presence of species divergence at the level of the Protein Interaction Networks. Additionally, it has been conjectured that the intra-module connections stay conserved whereas the inter-module connections change during evolution. The most important components of the biological networks are the functional modules which are more conserved among the evolutionary closer species. Here, we demonstrate an alternative method to decipher biological evolution by exploiting the topology of a spatially localized Protein Interaction Network called Protein Locality Graph (PLG). Our lossless graph compression from PLG to a power graph called Protein Cluster Interaction Network (PCIN) results in a 90% size reduction and aids in improving computational time. Further, we exploit the topology of PCIN and demonstrate our capability of deriving the correct species tree by focusing on the cross-talk between the protein modules exclusively. Our results provide new evidence that traces of evolution are not only present at the level of the Protein-Protein Interactions, but are also very much present at the level of the inter-module interactions.


title: The Foundational data initiative for Parkinsons disease (FOUNDIN-PD): enabling efficient translation from genetic maps to mechanism

abstract: In the FOUNdational Data INitiative for Parkinsons Disease (FOUNDIN-PD) we sought to produce a multi-layered molecular dataset in a large cohort of 95 Induced pluripotent stem cells (iPSC) lines at multiple timepoints during differentiation to dopaminergic (DA) neurons, a major affected cell type in Parkinsons Disease (PD). The lines are derived from the Parkinsons Progression Markers Initiative (PPMI) study that includes both people with PD and unaffected individuals across a wide range of polygenic risk scores (PRS) with both risk variants identified by genome-wide association studies (GWAS), and monogenic causal alleles. We generated genetic, epigenetic, regulatory, transcriptomic, proteomic, and longitudinal cellular imaging data from iPSC-derived DA neurons to understand key molecular relationships between disease associated genetic variation and proximate molecular events in a PD relevant cell-type. Analyses of all data modalities collected in FOUNDIN-PD suggest that the differentiation to DA neurons, while not fully mature, was successful and robust. Interrogation of PD genetic risk in this relevant cellular context may elucidate the functional effects of some of these risk variants alone or in combination with other variants. These data reveal that DA neurons derived from human iPSC provide a valuable cellular context and foundational atlas for modeling PD-related genetic risk. In addition to making the data and analyses for this molecular atlas readily available, we have integrated these data into the browsable FOUNDIN-PD data portal (https://www.foundinpd.org) to be used as a resource for understanding the molecular pathogenesis of PD.


title: Opportunity or catastrophe? Effect of sea salt on host-parasite survival and reproduction

abstract:
Seawater intrusion caused by anthropogenic climate change may affect freshwater species and their parasites. While brackish water certainly impacts freshwater systems globally, its impact on disease transmission is largely unknown. This study examined the effect of artificial seawater on host-parasite interactions using a freshwater snail host, Biomphalaria alexandrina, and the human trematode parasite Schistosoma mansoni. Four components were analyzed to evaluate the impact of increasing salinity on disease transmission: snail survival, snail reproduction, infection prevalence, and the survival of the parasite infective stage (cercariae). We found a decrease in snail survival, snail egg mass production, and snail infection prevalence as salinity increases. However, cercarial survival peaked at an intermediate salinity value. Our results suggest that seawater intrusion into freshwaters has the potential to decrease schistosome transmission to humans.

Author SummaryClimate change will have numerous impacts on many systems, including host-parasite systems. One mechanism by which climate change with impact host-parasite interactions is by rising sea levels flooding coastal regions, increasing salinity in many freshwaters. Host-parasite interactions are a key component of freshwater ecosystems, but the effects of sea water intrusion on host-parasite dynamics are largely unknown. In this study, we quantify the effects of sea salt concentration on the model host-parasite system, Biomphalaria alexandrina and Schistosoma mansoni. We demonstrate a significant, negative relationship between sea salt concentration and host-parasite survival and reproduction. The increase in freshwater salinity associated with sea level rise has the potential to decrease parasite transmission and disease burden in humans and wildlife.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=120 SRC="FIGDIR/small/446887v1_ufig1.gif" ALT="Figure 1">
View larger version (23K):
org.highwire.dtl.DTLVardef@529bd5org.highwire.dtl.DTLVardef@174cc80org.highwire.dtl.DTLVardef@11e0912org.highwire.dtl.DTLVardef@e5ab7a_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Improving Timeliness in the Neglected Tropical Diseases Preventive Chemotherapy Donation Supply Chain through Information Sharing: A Retrospective Empirical Analysis

abstract:
BackgroundBillions of doses of medicines are donated for mass drug administrations in support of the World Health Organizations "Roadmap to Implementation," which aims to control, eliminate, and eradicate Neglected Tropical Diseases (NTDs). The supply chain to deliver these medicines is complex, with fragmented data systems and limited visibility on performance. This study empirically evaluates the impact of an online supply chain performance measurement system, "NTDeliver," providing understanding of the value of information sharing towards the success of global health programs.

MethodsRetrospective secondary data was extracted from NTDeliver, which included 1,484 shipments for four critical medicines ordered by over 100 countries between February 28, 2006 and December 31, 2018. We applied statistical regression models to analyze the impact on key performance metrics, comparing data before and after the system was implemented.

FindingsThe results suggest information sharing has a positive impact on three performance indicators: purchase order timeliness ({beta}=1.01, p<0.000), arrival timeliness ({beta}=0.53, p=0.09), and--most importantly--delivery timeliness ({beta}=0.73, p=0.03). Three variables indicated an increased positive impact when the data is publicly shared: shipment timeliness ({beta}=2.57, p=0.001), arrival timeliness ({beta}=2.88, p=0{middle dot}003), and delivery timeliness ({beta}=2.82, p=0.01).

ConclusionsOur findings suggest that information sharing between the NTD program partners can help drive improved performance in the supply chain, and even more so when data is shared publicly. Given the large volume of medicine and the significant number of people requiring these medicines, information sharing has the potential to provide improvements to global health programs affecting the health of tens to hundreds of millions of people

Author SummaryThe supply chain to deliver donated preventive chemotherapy medicines is complex due to the many stakeholders and partnerships participating, as well as challenging because the logistics are further complicated by delivery to remote destinations in developing countries. As MDA campaigns involve treating hundreds of thousands to millions of patients in endemic regions within entire countries over the course of days or weeks, close coordination and timing of medicine delivery is critical. Inefficiencies caused by fragmented data systems and limited transparency on supply chain performance further challenges the ability to identify shipment issues and explore the root cause of the issues. Prior to 2016, delivery was performing below standards, lagging as much as 40% below the WHO target of 80% on-time delivery. These delays result in wasted medicine donations, increased program costs, delayed MDAs, or sometimes even completely missed MDAs. In September 2016, an online supply chain performance measurement system, "NTDeliver," was launched by the NTD Supply Chain Forum (a public-private partnership focused on managing and improving the PC donation supply chain) to enhance supply chain performance and information transparency. Our findings suggest that information sharing through NTDeliver is positively associated with performance at several key stages in NTD supply chain and that information sharing has more substantial positive impact on performance when the information is made publicly accessible, focused towards country program managers. The study findings support investment in supply chain systems and commitment to data transparency, in the context of a growing focus on supply chain investment in NTD programs.

category: microbiology


title: Oil palm expansion increases the vectorial capacity of dengue vectors in Malaysian Borneo

abstract:
Changes in land-use and the associated shifts in environmental conditions can have large effects on the transmission and emergence of disease. Mosquito-borne disease are particularly sensitive to these changes because mosquito growth, reproduction, survival and susceptibility to infection are all thermally sensitive traits, and land use change dramatically alters local microclimate. Predicting disease transmission under environmental change is increasingly critical for targeting mosquito-borne disease control and for identifying hotspots of disease emergence. Mechanistic models offer a powerful tool for improving these predications. However, these approaches are limited by the quality and scale of temperature data and the thermal response curves that underlie predictions. Here, we used fine-scale temperature monitoring and a combination of empirical, laboratory and temperature-dependent estimates to estimate the vectorial capacity of Aedes albopictus mosquitoes across a tropical forest - oil palm plantation conversion gradient in Malaysian Borneo. We found that fine-scale differences in temperature between logged forest and oil palm plantation sites were not sufficient to produce differences in temperature-dependent trait estimates using published thermal performance curves. However, when measured under field conditions a key parameter, adult abundance, differed significantly between land-use types, resulting in estimates of vectorial capacity that were 1.5 times higher in plantations than in forests. The prediction that oil palm plantations would support mosquito populations with higher vectorial capacity was robust to uncertainties in our adult survival estimates. These results provide a mechanistic basis for understanding the effects of forest conversion on mosquito-borne disease risk, and a framework for interpreting emergent relationships between land-use and disease transmission. As rising demand for palm oil products drives continued expansion of plantations, these findings have important implications for conservation, land management and public health policy at the global scale.

AUTHOR SUMMARYThe large scale modification of landscapes by humans has contributed to the rise of emerging and re-emerging mosquito-borne diseases. While the association between athropogenic land-use change and these shifts in disease risk are frequently observed, our understanding of exactly how land-use change mechanistically alters disease risk remains unclear. Changes in local environmental conditions (e.g. temperature) may play an important role, due to the effects on mosquito life-history, but are rarely measured at scales relevant these small-bodied ectotherms. Here we measure the impact of tropical forest conversion to oil palm plantation on each of the components that determine the potential for mosquitoes to transmit pathogens (vectorial capacity). By combining fine-scale temperature data from the field, published temperature responses of mosquitoes and field data on adult mosquito populations, we show that land-use change from forest to plantation can be expected to increase vectorial capacity by 150%. Our results highlight the need to advance field research into fundamental mosquito ecology, and to more critically evaluate the increased risk of Aedes-borne disease in dynamic working landscapes against the benefits of economic development.