Quantitative Analysis of Mechanisms That Govern Red Blood Cell Age Structure and Dynamics during Anaemia

Mathematical modelling has proven an important tool in elucidating and quantifying mechanisms that govern the age structure and population dynamics of red blood cells (RBCs). Here we synthesise ideas from previous experimental data and the mathematical modelling literature with new data in order to test hypotheses and generate new predictions about these mechanisms. The result is a set of competing hypotheses about three intrinsic mechanisms: the feedback from circulating RBC concentration to production rate of immature RBCs (reticulocytes) in bone marrow, the release of reticulocytes from bone marrow into the circulation, and their subsequent ageing and clearance. In addition we examine two mechanisms specific to our experimental system: the effect of phenylhydrazine (PHZ) and blood sampling on RBC dynamics. We performed a set of experiments to quantify the dynamics of reticulocyte proportion, RBC concentration, and erythropoietin concentration in PHZ-induced anaemic mice. By quantifying experimental error we are able to fit and assess each hypothesis against our data and recover parameter estimates using Markov chain Monte Carlo based Bayesian inference. We find that, under normal conditions, about 3% of reticulocytes are released early from bone marrow and upon maturation all cells are released immediately. In the circulation, RBCs undergo random clearance but have a maximum lifespan of about 50 days. Under anaemic conditions reticulocyte production rate is linearly correlated with the difference between normal and anaemic RBC concentrations, and their release rate is exponentially correlated with the same. PHZ appears to age rather than kill RBCs, and younger RBCs are affected more than older RBCs. Blood sampling caused short aperiodic spikes in the proportion of reticulocytes which appear to have a different developmental pathway than normal reticulocytes. We also provide evidence of large diurnal oscillations in serum erythropoietin levels during anaemia.     

Antiviral antibodies are necessary to prevent cytotoxic T-lymphocyte escape in mice infected with a coronavirus

Mutation within virus-derived CD8 T-cell epitopes can effectively abrogate cytotoxic T-lymphocyte (CTL) recognition and impede virus clearance in infected hosts. These so-called “CTL escape variant viruses” are commonly selected during persistent infections and are associated with rapid disease progression and increased disease severity. Herein, we tested whether antiviral antibody-mediated suppression of virus replication and subsequent virus clearance were necessary for preventing CTL escape in coronavirus-infected mice. We found that compared to wild-type mice, B-cell-deficient mice did not efficiently clear infectious virus, uniformly developed clinical disease, and harbored CTL escape variant viruses. These data directly demonstrate a critical role for antiviral antibody in protecting from the selective outgrowth of CTL escape variant viruses.     

Ecological rules for the assembly of microbiome communities

Humans and many other hosts establish a diverse community of beneficial microbes anew each generation. The order and identity of incoming symbionts is critical for health, but what determines the success of the assembly process remains poorly understood. Here we develop ecological theory to identify factors important for microbial community assembly. Our method maps out all feasible pathways for the assembly of a given microbiome—with analogies to the mutational maps underlying fitness landscapes in evolutionary biology. Building these “assembly maps” reveals a tradeoff at the heart of the assembly process. Ecological dependencies between members of the microbiota make assembly predictable—and can provide metabolic benefits to the host—but these dependencies may also create barriers to assembly. This effect occurs because interdependent species can fail to establish when each relies on the other to colonize first. We support our predictions with published data from the assembly of the preterm infant microbiota, where we find that ecological dependence is associated with a predictable order of arrival. Our models also suggest that hosts can overcome barriers to assembly via mechanisms that either promote the uptake of multiple symbiont species in one step or feed early colonizers. This predicted importance of host feeding is supported by published data on the impacts of breast milk in the assembly of the human microbiome. We conclude that both microbe to microbe and host to microbe interactions are important for the trajectory of microbiome assembly.

Humans and many other hosts establish a diverse community of beneficial microbes anew each generation, but what determines the success of the assembly process remains poorly understood. This study develops ecological theory that reveals the rules underlying the assembly of such host-associated microbiota.     

Competing for blood: the ecology of parasite resource competition in human malaria–helminth co‐infections

Ecological theory suggests that co‐infecting parasite species can interact within hosts directly, via host immunity and/or via resource competition. In mice, competition for red blood cells (RBCs) between malaria and bloodsucking helminths can regulate malaria population dynamics, but the importance of RBC competition in human hosts was unknown. We analysed infection density (i.e. the concentration of parasites in infected hosts), from a 2‐year deworming study of over 4000 human subjects. After accounting for resource‐use differences among parasites, we find evidence of resource competition, priority effects and a competitive hierarchy within co‐infected individuals. For example reducing competition via deworming increased Plasmodium vivax densities 2.8‐fold, and this effect is limited to bloodsucking hookworms. Our ecological, resource‐based perspective sheds new light into decades of conflicting outcomes of malaria–helminth co‐infection studies with significant health and transmission consequences. Beyond blood, investigating within‐human resource competition may bring new insights for improving human health.     

CpG Dinucleotide Frequencies Reveal the Role of Host Methylation Capabilities in Parvovirus Evolution

Parvoviruses are rapidly evolving viruses that infect a wide range of hosts, including vertebrates and invertebrates. Extensive methylation of the parvovirus genome has been recently demonstrated. A global pattern of methylation of CpG dinucleotides is seen in vertebrate genomes, compared to “fractional” methylation patterns in invertebrate genomes. It remains unknown if the loss of CpG dinucleotides occurs in all viruses of a given DNA virus family that infect host species spanning across vertebrates and invertebrates. We investigated the link between the extent of CpG dinucleotide depletion among autonomous parvoviruses and the evolutionary lineage of the infected host. We demonstrate major differences in the relative abundance of CpG dinucleotides among autonomous parvoviruses which share similar genome organization and common ancestry, depending on the infected host species. Parvoviruses infecting vertebrate hosts had significantly lower relative abundance of CpG dinucleotides than parvoviruses infecting invertebrate hosts. The strong correlation of CpG dinucleotide depletion with the gain in TpG/CpA dinucleotides and the loss of TpA dinucleotides among parvoviruses suggests a major role for CpG methylation in the evolution of parvoviruses. Our data present evidence that links the relative abundance of CpG dinucleotides in parvoviruses to the methylation capabilities of the infected host. In sum, our findings support a novel perspective of host-driven evolution among autonomous parvoviruses.     

Diverse Phage-Encoded Toxins in a Protective Insect Endosymbiont

The lysogenic bacteriophage APSE infects “Candidatus Hamiltonella defensa,” a facultative endosymbiont of aphids and other sap-feeding insects. This endosymbiont has established a beneficial association with aphids, increasing survivorship following attack by parasitoid wasps. Although APSE and “Ca. Hamiltonella defensa” are effectively maternally transmitted between aphid generations, they can also be horizontally transferred among insect hosts, which results in genetically distinct “Ca. Hamiltonella defensa” strains infecting the same aphid species and sporadic distributions of both APSE and “Ca. Hamiltonella defensa” among hosts. Aphids infected only with “Ca. Hamiltonella defensa” have significantly less protection than those infected with both “Ca. Hamiltonella defensa” and APSE. This protection has been proposed to be connected to eukaryote-targeted toxins previously discovered in the genomes of two characterized APSE strains. In this study, we have sequenced partial genomes from seven additional APSE strains to address the evolution and extent of toxin variation in this phage. The APSE lysis region has been a hot spot for nonhomologous recombination of novel virulence cassettes. We identified four new toxins from three protein families, Shiga-like toxin, cytolethal distending toxin, and YD-repeat toxins. These recombination events have also resulted in reassortment of the downstream lysozyme and holin genes. Analysis of the conserved APSE genes flanking the variable toxin cassettes reveals a close phylogenetic association with phage sequences from two other facultative endosymbionts of insects. Thus, phage may act as a conduit for ongoing gene exchange among heritable endosymbionts.     

Division of labor between the pore-1 loops of the D1 and D2 AAA+ rings coordinates substrate selectivity of the ClpAP protease

ClpAP, an ATP-dependent protease consisting of ClpA, a double-ring hexameric unfoldase of the ATPases associated with diverse cellular activities superfamily, and the ClpP peptidase, degrades damaged and unneeded proteins to support cellular proteostasis. ClpA recognizes many protein substrates directly, but it can also be regulated by an adapter, ClpS, that modifies ClpA’s substrate profile toward N-degron substrates. Conserved tyrosines in the 12 pore-1 loops lining the central channel of the stacked D1 and D2 rings of ClpA are critical for degradation, but the roles of these residues in individual steps during direct or adapter-mediated degradation are poorly understood. Using engineered ClpA hexamers with zero, three, or six pore-1 loop mutations in each ATPases associated with diverse cellular activities superfamily ring, we found that active D1 pore loops initiate productive engagement of substrates, whereas active D2 pore loops are most important for mediating the robust unfolding of stable native substrates. In complex with ClpS, active D1 pore loops are required to form a high affinity ClpA•ClpS•substrate complex, but D2 pore loops are needed to “tug on” and remodel ClpS to transfer the N-degron substrate to ClpA. Overall, we find that the pore-1 loop tyrosines in D1 are critical for direct substrate engagement, whereas ClpS-mediated substrate delivery requires unique contributions from both the D1 and D2 pore loops. In conclusion, our study illustrates how pore loop engagement, substrate capture, and powering of the unfolding/translocation steps are distributed between the two rings of ClpA, illuminating new mechanistic features that may be common to double-ring protein unfolding machines.     

Re-Examination of Globally Flat Space-Time

In the following, we offer a novel approach to modeling the observed effects currently attributed to the theoretical concepts of “dark energy,” “dark matter,” and “dark flow.” Instead of assuming the existence of these theoretical concepts, we take an alternative route and choose to redefine what we consider to be inertial motion as well as what constitutes an inertial frame of reference in flat space-time. We adopt none of the features of our current cosmological models except for the requirement that special and general relativity be local approximations within our revised definition of inertial systems. Implicit in our ideas is the assumption that at “large enough” scales one can treat objects within these inertial systems as point-particles having an insignificant effect on the curvature of space-time. We then proceed under the assumption that time and space are fundamentally intertwined such that time- and spatial-translational invariance are not inherent symmetries of flat space-time (i.e., observable clock rates depend upon both relative velocity and spatial position within these inertial systems) and take the geodesics of this theory in the radial Rindler chart as the proper characterization of inertial motion. With this commitment, we are able to model solely with inertial motion the observed effects expected to be the result of “dark energy,” “dark matter,” and “dark flow.” In addition, we examine the potential observable implications of our theory in a gravitational system located within a confined region of an inertial reference frame, subsequently interpreting the Pioneer anomaly as support for our redefinition of inertial motion. As well, we extend our analysis into quantum mechanics by quantizing for a real scalar field and find a possible explanation for the asymmetry between matter and antimatter within the framework of these redefined inertial systems.     

A Mapping of Drug Space from the Viewpoint of Small Molecule Metabolism

Small molecule drugs target many core metabolic enzymes in humans and pathogens, often mimicking endogenous ligands. The effects may be therapeutic or toxic, but are frequently unexpected. A large-scale mapping of the intersection between drugs and metabolism is needed to better guide drug discovery. To map the intersection between drugs and metabolism, we have grouped drugs and metabolites by their associated targets and enzymes using ligand-based set signatures created to quantify their degree of similarity in chemical space. The results reveal the chemical space that has been explored for metabolic targets, where successful drugs have been found, and what novel territory remains. To aid other researchers in their drug discovery efforts, we have created an online resource of interactive maps linking drugs to metabolism. These maps predict the “effect space” comprising likely target enzymes for each of the 246 MDDR drug classes in humans. The online resource also provides species-specific interactive drug-metabolism maps for each of the 385 model organisms and pathogens in the BioCyc database collection. Chemical similarity links between drugs and metabolites predict potential toxicity, suggest routes of metabolism, and reveal drug polypharmacology. The metabolic maps enable interactive navigation of the vast biological data on potential metabolic drug targets and the drug chemistry currently available to prosecute those targets. Thus, this work provides a large-scale approach to ligand-based prediction of drug action in small molecule metabolism.     

The mismatched nucleotides encoded in vaccinia virus flip-and-flop hairpin telomeres serve an essential role in virion maturation

Poxvirus genomes consist of a linear duplex DNA that ends in short inverted and complementary hairpin structures. These elements also encode loops and mismatches that likely serve a role in genome packaging and perhaps replication. We constructed mutant vaccinia viruses (VACV) where the native hairpins were replaced by altered forms and tested effects on replication, assembly, and virulence. Our studies showed that structure, not sequence, likely determines function as one can replace an Orthopoxvirus (VACV) hairpin with one copied from a Leporipoxvirus with no effect on growth. Some loops can be deleted from VACV hairpins with little effect, but VACV bearing too few mismatches grew poorly and we couldn’t recover viruses lacking all mismatches. Further studies were conducted using a mutant bearing only one of six mismatches found in wild-type hairpins (SΔ1Δ3–6). This virus grew to ~20-fold lower titers, but neither DNA synthesis nor telomere resolution was affected. However, the mutant exhibited a particle-to-PFU ratio 10-20-fold higher than wild-type viruses and p4b/4b core protein processing was compromised, indicating an assembly defect. Electron microscopy showed that SΔ1Δ3–6 mutant development was blocked at the immature virus (IV) stage, which phenocopies known effects of I1L mutants. Competitive DNA binding assays showed that recombinant I1 protein had less affinity for the SΔ1Δ3–6 hairpin than the wild-type hairpin. The SΔ1Δ3–6 mutant was also attenuated when administered to SCID-NCR mice by tail scarification. Mice inoculated with viruses bearing wild-type hairpins exhibited a median survival of 30–37 days, while mice infected with SΔ1Δ3–6 virus survived >70 days. Persistent infections favor genetic reversion and genome sequencing detected one example where a small duplication near the hairpin tip likely created a new loop. These observations show that mismatches serve a critical role in genome packaging and provide new insights into how VACV “flip and flop” telomeres are arranged.     

Adaptive search space pruning in complex strategic problems

People have limited computational resources, yet they make complex strategic decisions over enormous spaces of possibilities. How do people efficiently search spaces with combinatorially branching paths? Here, we study players’ search strategies for a winning move in a “k-in-a-row” game. We find that players use scoring strategies to prune the search space and augment this pruning by a “shutter” heuristic that focuses the search on the paths emanating from their previous move. This strong pruning has its costs—both computational simulations and behavioral data indicate that the shutter size is correlated with players’ blindness to their opponent’s winning moves. However, simulations of the search while varying the shutter size, complexity levels, noise levels, branching factor, and computational limitations indicate that despite its costs, a narrow shutter strategy is the dominant strategy for most of the parameter space. Finally, we show that in the presence of computational limitations, the shutter heuristic enhances the performance of deep learning networks in these end-game scenarios. Together, our findings suggest a novel adaptive heuristic that benefits search in a vast space of possibilities of a strategic game.     

Bracoviruses recruit host integrases for their integration into caterpillar’s genome

Some DNA viruses infect host animals usually by integrating their DNAs into the host genome. However, the mechanisms for integration remain largely unknown. Here, we find that Cotesia vestalis bracovirus (CvBV), a polydnavirus of the parasitic wasp C. vestalis (Haliday), integrates its DNA circles into host Plutella xylostella (L.) genome by two distinct strategies, conservatively and randomly, through high-throughput sequencing analysis. We confirmed that the conservatively integrating circles contain an essential “8+5” nucleotides motif which is required for integration. Then we find CvBV circles are integrated into the caterpillar’s genome in three temporal patterns, the early, mid and late stage-integration. We further identify that three CvBV-encoded integrases are responsible for some, but not all of the virus circle integrations, indeed they mainly participate in the processes of early stage-integration. Strikingly, we find two P. xylostella retroviral integrases (PxIN1 and PxIN2) are highly induced upon wasp parasitism, and PxIN1 is crucial for integration of some other early-integrated CvBV circles, such as CvBV_04, CvBV_12 and CvBV_24, while PxIN2 is important for integration of a late-integrated CvBV circle, CvBV_21. Our data uncover a novel mechanism in which CvBV integrates into the infected host genome, not only by utilizing its own integrases, but also by recruiting host enzymes. These findings will strongly deepen our understanding of how bracoviruses regulate and integrate into their hosts.     

Non-Specific Manipulation of Gammarid Behaviour by P. minutus Parasite Enhances Their Predation by Definitive Bird Hosts

Trophically-transmitted parasites often change the phenotype of their intermediate hosts in ways that increase their vulnerability to definitive hosts, hence favouring transmission. As a “collateral damage”, manipulated hosts can also become easy prey for non-host predators that are dead ends for the parasite, and which are supposed to play no role in transmission strategies. Interestingly, infection with the acanthocephalan parasite Polymorphus minutus has been shown to reduce the vulnerability of its gammarid intermediate hosts to non-host predators, whose presence triggered the behavioural alterations expected to favour trophic transmission to bird definitive hosts. Whilst the behavioural response of infected gammarids to the presence of definitive hosts remains to be investigated, this suggests that trophic transmission might be promoted by non-host predation risk. We conducted microcosm experiments to test whether the behaviour of P. minutus-infected gammarids was specific to the type of predator (i.e. mallard as definitive host and fish as non-host), and mesocosm experiments to test whether trophic transmission to bird hosts was influenced by non-host predation risk. Based on the behaviours we investigated (predator avoidance, activity, geotaxis, conspecific attraction), we found no evidence for a specific fine-tuned response in infected gammarids, which behaved similarly whatever the type of predator (mallard or fish). During predation tests, fish predation risk did not influence the differential predation of mallards that over-consumed infected gammarids compared to uninfected individuals. Overall, our results bring support for a less sophisticated scenario of manipulation than previously expected, combining chronic behavioural alterations with phasic behavioural alterations triggered by the chemical and physical cues coming from any type of predator. Given the wide dispersal range of waterbirds (the definitive hosts of P. minutus), such a manipulation whose efficiency does not depend on the biotic context is likely to facilitate its trophic transmission in a wide range of aquatic environments.     

Contact Trees: Network Visualization beyond Nodes and Edges

Node-Link diagrams make it possible to take a quick glance at how nodes (or actors) in a network are connected by edges (or ties). A conventional network diagram of a “contact tree” maps out a root and branches that represent the structure of nodes and edges, often without further specifying leaves or fruits that would have grown from small branches. By furnishing such a network structure with leaves and fruits, we reveal details about “contacts” in our ContactTrees upon which ties and relationships are constructed. Our elegant design employs a bottom-up approach that resembles a recent attempt to understand subjective well-being by means of a series of emotions. Such a bottom-up approach to social-network studies decomposes each tie into a series of interactions or contacts, which can help deepen our understanding of the complexity embedded in a network structure. Unlike previous network visualizations, ContactTrees highlight how relationships form and change based upon interactions among actors, as well as how relationships and networks vary by contact attributes. Based on a botanical tree metaphor, the design is easy to construct and the resulting tree-like visualization can display many properties at both tie and contact levels, thus recapturing a key ingredient missing from conventional techniques of network visualization. We demonstrate ContactTrees using data sets consisting of up to three waves of 3-month contact diaries over the 2004-2012 period, and discuss how this design can be applied to other types of datasets.     

Confidence Intervals for fMRI Activation Maps

Neuroimaging activation maps typically color voxels to indicate whether the blood oxygen level-dependent (BOLD) signals measured among two or more experimental conditions differ significantly at that location. This data presentation, however, omits information critical for interpretation of experimental results. First, no information is represented about trends at voxels that do not pass the statistical test. Second, no information is given about the range of probable effect sizes at voxels that do pass the statistical test. This leads to a fundamental error in interpreting activation maps by naïve viewers, where it is assumed that colored, “active” voxels are reliably different from uncolored “inactive” voxels. In other domains, confidence intervals have been added to data graphics to reduce such errors. Here, we first document the prevalence of the fundamental error of interpretation, and then present a method for solving it by depicting confidence intervals in fMRI activation maps. Presenting images where the bounds of confidence intervals at each voxel are coded as color allows readers to visually test for differences between “active” and “inactive” voxels, and permits for more proper interpretation of neuroimaging data. Our specific graphical methods are intended as initial proposals to spur broader discussion of how to present confidence intervals for fMRI data.     

Novel method for evaluating the health condition of mice in space through a video downlink

Clarification of the criteria for managing animal health is essential to increase the reliability of experiments and ensure transparency in animal welfare. For experiments performed in space, there is no consensus on how to care for animals owing to technical issues, launch mass limitation, and human resources. Some biological processes in mammals, such as musculoskeletal or immune processes, are altered in the space environment, and mice in space can be used to simulate morbid states, such as senescence acceleration. Thus, there is a need to establish a novel evaluation method and evaluation criteria to monitor animal health. Here, we report a novel method to evaluate the health of mice in space through a video downlink in a series of space experiments using the Multiple Artificial-gravity Research System (MARS). This method was found to be more useful in evaluating animal health in space than observations and body weight changes of the same live mice following their return to Earth. We also developed criteria to evaluate health status via a video downlink. These criteria, with “Fur condition” and “Respiratory” as key items, provided information on the daily changes in the health status of mice and helped to identify malfunctions at an early stage. Our method and criteria led to the success of our missions, and they will help establish appropriate rules for space experiments in the future.     

Peptide Inhibitors of Dengue-Virus Entry Target a Late-Stage Fusion Intermediate

The mechanism of membrane fusion by “class II” viral fusion proteins follows a pathway that involves large-scale domain rearrangements of the envelope glycoprotein (E) and a transition from dimers to trimers. The rearrangement is believed to proceed by an outward rotation of the E ectodomain after loss of the dimer interface, followed by a reassociation into extended trimers. The ∼55-aa-residue, membrane proximal “stem” can then zip up along domain II, bringing together the transmembrane segments of the C-terminus and the fusion loops at the tip of domain II. We find that peptides derived from the stem of dengue-virus E bind stem-less E trimer, which models a conformational intermediate. In vitro assays demonstrate that these peptides specifically block viral fusion. The peptides inhibit infectivity with potency proportional to their affinity for the conformational intermediate, even when free peptide is removed from a preincubated inoculum before infecting cells. We conclude that peptides bind virions before attachment and are carried with virions into endosomes, the compartment in which acidification initiates fusion. Binding depends on particle dynamics, as there is no inhibition of infectivity if preincubation and separation are at 4°C rather than 37°C. We propose a two-step model for the mechanism of fusion inhibition. Targeting a viral entry pathway can be an effective way to block infection. Our data, which support and extend proposed mechanisms for how the E conformational change promotes membrane fusion, suggest strategies for inhibiting flavivirus entry.     

Social Distance Evaluation in Human Parietal Cortex

Across cultures, social relationships are often thought of, described, and acted out in terms of physical space (e.g. “close friends” “high lord”). Does this cognitive mapping of social concepts arise from shared brain resources for processing social and physical relationships? Using fMRI, we found that the tasks of evaluating social compatibility and of evaluating physical distances engage a common brain substrate in the parietal cortex. The present study shows the possibility of an analytic brain mechanism to process and represent complex networks of social relationships. Given parietal cortex''s known role in constructing egocentric maps of physical space, our present findings may help to explain the linguistic, psychological and behavioural links between social and physical space.     

Bayesian Analysis of Congruence of Core Genes in Prochlorococcus and Synechococcus and Implications on Horizontal Gene Transfer

It is often suggested that horizontal gene transfer is so ubiquitous in microbes that the concept of a phylogenetic tree representing the pattern of vertical inheritance is oversimplified or even positively misleading. “Universal proteins” have been used to infer the organismal phylogeny, but have been criticized as being only the “tree of one percent.” Currently, few options exist for those wishing to rigorously assess how well a universal protein phylogeny, based on a relative handful of well-conserved genes, represents the phylogenetic histories of hundreds of genes. Here, we address this problem by proposing a visualization method and a statistical test within a Bayesian framework. We use the genomes of marine cyanobacteria, a group thought to exhibit substantial amounts of HGT, as a test case. We take 379 orthologous gene families from 28 cyanobacteria genomes and estimate the Bayesian posterior distributions of trees – a “treecloud” – for each, as well as for a concatenated dataset based on putative “universal proteins.” We then calculate the average distance between trees within and between all treeclouds on various metrics and visualize this high-dimensional space with non-metric multidimensional scaling (NMMDS). We show that the tree space is strongly clustered and that the universal protein treecloud is statistically significantly closer to the center of this tree space than any individual gene treecloud. We apply several commonly-used tests for incongruence/HGT and show that they agree HGT is rare in this dataset, but make different choices about which genes were subject to HGT. Our results show that the question of the representativeness of the “tree of one percent” is a quantitative empirical question, and that the phylogenetic central tendency is a meaningful observation even if many individual genes disagree due to the various sources of incongruence.