Unraveling Genomic Complexity at a Quantitative Disease Resistance Locus in Maize

Multiple disease resistance has important implications for plant fitness, given the selection pressure that many pathogens exert directly on natural plant populations and indirectly via crop improvement programs. Evidence of a locus conditioning resistance to multiple pathogens was found in bin 1.06 of the maize genome with the allele from inbred line “Tx303” conditioning quantitative resistance to northern leaf blight (NLB) and qualitative resistance to Stewart’s wilt. To dissect the genetic basis of resistance in this region and to refine candidate gene hypotheses, we mapped resistance to the two diseases. Both resistance phenotypes were localized to overlapping regions, with the Stewart’s wilt interval refined to a 95.9-kb segment containing three genes and the NLB interval to a 3.60-Mb segment containing 117 genes. Regions of the introgression showed little to no recombination, suggesting structural differences between the inbred lines Tx303 and “B73,” the parents of the fine-mapping population. We examined copy number variation across the region using next-generation sequencing data, and found large variation in read depth in Tx303 across the region relative to the reference genome of B73. In the fine-mapping region, association mapping for NLB implicated candidate genes, including a putative zinc finger and pan1. We tested mutant alleles and found that pan1 is a susceptibility gene for NLB and Stewart’s wilt. Our data strongly suggest that structural variation plays an important role in resistance conditioned by this region, and pan1, a gene conditioning susceptibility for NLB, may underlie the QTL.     

Growth-dependent bacterial susceptibility to ribosome-targeting antibiotics

Bacterial growth environment strongly influences the efficacy of antibiotic treatment, with slow growth often being associated with decreased susceptibility. Yet in many cases, the connection between antibiotic susceptibility and pathogen physiology remains unclear. We show that for ribosome-targeting antibiotics acting on Escherichia coli, a complex interplay exists between physiology and antibiotic action; for some antibiotics within this class, faster growth indeed increases susceptibility, but for other antibiotics, the opposite is true. Remarkably, these observations can be explained by a simple mathematical model that combines drug transport and binding with physiological constraints. Our model reveals that growth-dependent susceptibility is controlled by a single parameter characterizing the ‘reversibility’ of ribosome-targeting antibiotic transport and binding. This parameter provides a spectrum classification of antibiotic growth-dependent efficacy that appears to correspond at its extremes to existing binary classification schemes. In these limits, the model predicts universal, parameter-free limiting forms for growth inhibition curves. The model also leads to non-trivial predictions for the drug susceptibility of a translation mutant strain of E. coli, which we verify experimentally. Drug action and bacterial metabolism are mechanistically complex; nevertheless, this study illustrates how coarse-grained models can be used to integrate pathogen physiology into drug design and treatment strategies.     

Enhancements to the Rosetta Energy Function Enable Improved Identification of Small Molecules that Inhibit Protein-Protein Interactions

Protein-protein interactions are among today’s most exciting and promising targets for therapeutic intervention. To date, identifying small-molecules that selectively disrupt these interactions has proven particularly challenging for virtual screening tools, since these have typically been optimized to perform well on more “traditional” drug discovery targets. Here, we test the performance of the Rosetta energy function for identifying compounds that inhibit protein interactions, when these active compounds have been hidden amongst pools of “decoys.” Through this virtual screening benchmark, we gauge the effect of two recent enhancements to the functional form of the Rosetta energy function: the new “Talaris” update and the “pwSHO” solvation model. Finally, we conclude by developing and validating a new weight set that maximizes Rosetta’s ability to pick out the active compounds in this test set. Looking collectively over the course of these enhancements, we find a marked improvement in Rosetta’s ability to identify small-molecule inhibitors of protein-protein interactions.     

A Population Genetic Model for the Initial Spread of Partially Resistant Malaria Parasites under Anti-Malarial Combination Therapy and Weak Intrahost Competition

To develop public-health policies that extend the lifespan of affordable anti-malarial drugs as effective treatment options, it is necessary to understand the evolutionary processes leading to the origin and spread of mutations conferring drug resistance in malarial parasites. We built a population-genetic model for the emergence of resistance under combination drug therapy. Reproductive cycles of parasites are specified by their absolute fitness determined by clinical parameters, thus coupling the evolutionary-genetic with population-dynamic processes. Initial mutations confer only partial drug-resistance. Therefore, mutant parasites rarely survive combination therapy and within-host competition is very weak among parasites. The model focuses on the early phase of such unsuccessful recurrent mutations. This ends in the rare event of mutants enriching in an infected individual from which the successful spread of resistance over the entire population is initiated. By computer simulations, the waiting time until the establishment of resistant parasites is analysed. Resistance spreads quickly following the first appearance of a host infected predominantly by mutant parasites. This occurs either through a rare transmission of a resistant parasite to an uninfected host or through a rare failure of drugs in removing “transient” mutant alleles. The emergence of resistance is delayed with lower mutation rate, earlier treatment, higher metabolic cost of resistance, longer duration of high drug dose, and higher drug efficacy causing a stronger reduction in the sensitive and resistant parasites’ fitnesses. Overall, contrary to other studies’ proposition, the current model based on absolute fitness suggests that aggressive drug treatment delays the emergence of drug resistance.     

Evolutionary and Topological Properties of Genes and Community Structures in Human Gene Regulatory Networks

The diverse, specialized genes present in today’s lifeforms evolved from a common core of ancient, elementary genes. However, these genes did not evolve individually: gene expression is controlled by a complex network of interactions, and alterations in one gene may drive reciprocal changes in its proteins’ binding partners. Like many complex networks, these gene regulatory networks (GRNs) are composed of communities, or clusters of genes with relatively high connectivity. A deep understanding of the relationship between the evolutionary history of single genes and the topological properties of the underlying GRN is integral to evolutionary genetics. Here, we show that the topological properties of an acute myeloid leukemia GRN and a general human GRN are strongly coupled with its genes’ evolutionary properties. Slowly evolving (“cold”), old genes tend to interact with each other, as do rapidly evolving (“hot”), young genes. This naturally causes genes to segregate into community structures with relatively homogeneous evolutionary histories. We argue that gene duplication placed old, cold genes and communities at the center of the networks, and young, hot genes and communities at the periphery. We demonstrate this with single-node centrality measures and two new measures of efficiency, the set efficiency and the interset efficiency. We conclude that these methods for studying the relationships between a GRN’s community structures and its genes’ evolutionary properties provide new perspectives for understanding evolutionary genetics.     

Construction of disease-specific cytokine profiles by associating disease genes with immune responses

The pathogenesis of many inflammatory diseases is a coordinated process involving metabolic dysfunctions and immune response—usually modulated by the production of cytokines and associated inflammatory molecules. In this work, we seek to understand how genes involved in pathogenesis which are often not associated with the immune system in an obvious way communicate with the immune system. We have embedded a network of human protein-protein interactions (PPI) from the STRING database with 14,707 human genes using feature learning that captures high confidence edges. We have found that our predicted Association Scores derived from the features extracted from STRING’s high confidence edges are useful for predicting novel connections between genes, thus enabling the construction of a full map of predicted associations for all possible pairs between 14,707 human genes. In particular, we analyzed the pattern of associations for 126 cytokines and found that the six patterns of cytokine interaction with human genes are consistent with their functional classifications. To define the disease-specific roles of cytokines we have collected gene sets for 11,944 diseases from DisGeNET. We used these gene sets to predict disease-specific gene associations with cytokines by calculating the normalized average Association Scores between disease-associated gene sets and the 126 cytokines; this creates a unique profile of inflammatory genes (both known and predicted) for each disease. We validated our predicted cytokine associations by comparing them to known associations for 171 diseases. The predicted cytokine profiles correlate (p-value<0.0003) with the known ones in 95 diseases. We further characterized the profiles of each disease by calculating an “Inflammation Score” that summarizes different modes of immune responses. Finally, by analyzing subnetworks formed between disease-specific pathogenesis genes, hormones, receptors, and cytokines, we identified the key genes responsible for interactions between pathogenesis and inflammatory responses. These genes and the corresponding cytokines used by different immune disorders suggest unique targets for drug discovery.     

Exploring a Local Genetic Interaction Network Using Evolutionary Replay Experiments

Understanding how genes interact is a central challenge in biology. Experimental evolution provides a useful, but underutilized, tool for identifying genetic interactions, particularly those that involve non-loss-of-function mutations or mutations in essential genes. We previously identified a strong positive genetic interaction between specific mutations in KEL1 (P344T) and HSL7 (A695fs) that arose in an experimentally evolved Saccharomyces cerevisiae population. Because this genetic interaction is not phenocopied by gene deletion, it was previously unknown. Using “evolutionary replay” experiments, we identified additional mutations that have positive genetic interactions with the kel1-P344T mutation. We replayed the evolution of this population 672 times from six timepoints. We identified 30 populations where the kel1-P344T mutation reached high frequency. We performed whole-genome sequencing on these populations to identify genes in which mutations arose specifically in the kel1-P344T background. We reconstructed mutations in the ancestral and kel1-P344T backgrounds to validate positive genetic interactions. We identify several genetic interactors with KEL1, we validate these interactions by reconstruction experiments, and we show these interactions are not recapitulated by loss-of-function mutations. Our results demonstrate the power of experimental evolution to identify genetic interactions that are positive, allele specific, and not readily detected by other methods, shedding light on an underexplored region of the yeast genetic interaction network.     

The Interaction between Osmotic- and Pressure-induced Water Flow in Plant Roots

This paper presents a general model for coupled solute and water flow through plant roots based on the thermodynamics of irreversible processes. The model explains in a straight-forward manner such experimentally observed phenomena as changes in root resistance, increased solute flux, and apparent negative resistance, which have been reported for root systems under the influence of a hydrostatic pressure gradient. These apparent anomalies are explained on the basis of the interaction between the osmotic and hydrostatic driving forces and the well known “sweeping away” or dilution effect. We show that with a constant hydraulic conductivity the only features necessary to explain these phenomena are some type of membrane or membranelike structure and a mechanism for actively accumulating solutes.     

Adaptive Landscape by Environment Interactions Dictate Evolutionary Dynamics in Models of Drug Resistance

The adaptive landscape analogy has found practical use in recent years, as many have explored how their understanding can inform therapeutic strategies that subvert the evolution of drug resistance. A major barrier to applications of these concepts is a lack of detail concerning how the environment affects adaptive landscape topography, and consequently, the outcome of drug treatment. Here we combine empirical data, evolutionary theory, and computer simulations towards dissecting adaptive landscape by environment interactions for the evolution of drug resistance in two dimensions—drug concentration and drug type. We do so by studying the resistance mediated by Plasmodium falciparum dihydrofolate reductase (DHFR) to two related inhibitors—pyrimethamine and cycloguanil—across a breadth of drug concentrations. We first examine whether the adaptive landscapes for the two drugs are consistent with common definitions of cross-resistance. We then reconstruct all accessible pathways across the landscape, observing how their structure changes with drug environment. We offer a mechanism for non-linearity in the topography of accessible pathways by calculating of the interaction between mutation effects and drug environment, which reveals rampant patterns of epistasis. We then simulate evolution in several different drug environments to observe how these individual mutation effects (and patterns of epistasis) influence paths taken at evolutionary “forks in the road” that dictate adaptive dynamics in silico. In doing so, we reveal how classic metrics like the IC₅₀ and minimal inhibitory concentration (MIC) are dubious proxies for understanding how evolution will occur across drug environments. We also consider how the findings reveal ambiguities in the cross-resistance concept, as subtle differences in adaptive landscape topography between otherwise equivalent drugs can drive drastically different evolutionary outcomes. Summarizing, we discuss the results with regards to their basic contribution to the study of empirical adaptive landscapes, and in terms of how they inform new models for the evolution of drug resistance.     

The Antibiotic Dosage of Fastest Resistance Evolution: Gene Amplifications Underpinning the Inverted-U

To determine the dosage at which antibiotic resistance evolution is most rapid, we treated Escherichia coli in vitro, deploying the antibiotic erythromycin at dosages ranging from zero to high. Adaptation was fastest just below erythromycin’s minimal inhibitory concentration (MIC) and genotype-phenotype correlations determined from whole genome sequencing revealed the molecular basis: simultaneous selection for copy number variation in three resistance mechanisms which exhibited an “inverted-U” pattern of dose-dependence, as did several insertion sequences and an integron. Many genes did not conform to this pattern, however, reflecting changes in selection as dose increased: putative media adaptation polymorphisms at zero antibiotic dosage gave way to drug target (ribosomal RNA operon) amplification at mid dosages whereas prophage-mediated drug efflux amplifications dominated at the highest dosages. All treatments exhibited E. coli increases in the copy number of efflux operons acrAB and emrE at rates that correlated with increases in population density. For strains where the inverted-U was no longer observed following the genetic manipulation of acrAB, it could be recovered by prolonging the antibiotic treatment at subMIC dosages.     

Towards a mechanistic understanding of reciprocal drug–microbiome interactions

Broad‐spectrum antibiotics target multiple gram‐positive and gram‐negative bacteria, and can collaterally damage the gut microbiota. Yet, our knowledge of the extent of damage, the antibiotic activity spectra, and the resistance mechanisms of gut microbes is sparse. This limits our ability to mitigate microbiome‐facilitated spread of antibiotic resistance. In addition to antibiotics, non‐antibiotic drugs affect the human microbiome, as shown by metagenomics as well as in vitro studies. Microbiome–drug interactions are bidirectional, as microbes can also modulate drugs. Chemical modifications of antibiotics mostly function as antimicrobial resistance mechanisms, while metabolism of non‐antibiotics can also change the drugs’ pharmacodynamic, pharmacokinetic, and toxic properties. Recent studies have started to unravel the extensive capacity of gut microbes to metabolize drugs, the mechanisms, and the relevance of such events for drug treatment. These findings raise the question whether and to which degree these reciprocal drug–microbiome interactions will differ across individuals, and how to take them into account in drug discovery and precision medicine. This review describes recent developments in the field and discusses future study areas that will benefit from systems biology approaches to better understand the mechanistic role of the human gut microbiota in drug actions.     

Modeling the ecology of parasitic plasmids

Plasmids are autonomous genetic elements that can be exchanged between microorganisms via horizontal gene transfer (HGT). Despite the central role they play in antibiotic resistance and modern biotechnology, our understanding of plasmids’ natural ecology is limited. Recent experiments have shown that plasmids can spread even when they are a burden to the cell, suggesting that natural plasmids may exist as parasites. Here, we use mathematical modeling to explore the ecology of such parasitic plasmids. We first develop models of single plasmids and find that a plasmid’s population dynamics and optimal infection strategy are strongly determined by the plasmid’s HGT mechanism. We then analyze models of co-infecting plasmids and show that parasitic plasmids are prone to a “tragedy of the commons” in which runaway plasmid invasion severely reduces host fitness. We propose that this tragedy of the commons is averted by selection between competing populations and demonstrate this effect in a metapopulation model. We derive predicted distributions of unique plasmid types in genomes—comparison to the distribution of plasmids in a collection of 17,725 genomes supports a model of parasitic plasmids with positive plasmid–plasmid interactions that ameliorate plasmid fitness costs or promote the invasion of new plasmids.Subject terms: Theoretical ecology, Microbial ecology     

Global Quantitative Modeling of Chromatin Factor Interactions

Chromatin is the driver of gene regulation, yet understanding the molecular interactions underlying chromatin factor combinatorial patterns (or the “chromatin codes”) remains a fundamental challenge in chromatin biology. Here we developed a global modeling framework that leverages chromatin profiling data to produce a systems-level view of the macromolecular complex of chromatin. Our model ultilizes maximum entropy modeling with regularization-based structure learning to statistically dissect dependencies between chromatin factors and produce an accurate probability distribution of chromatin code. Our unsupervised quantitative model, trained on genome-wide chromatin profiles of 73 histone marks and chromatin proteins from modENCODE, enabled making various data-driven inferences about chromatin profiles and interactions. We provided a highly accurate predictor of chromatin factor pairwise interactions validated by known experimental evidence, and for the first time enabled higher-order interaction prediction. Our predictions can thus help guide future experimental studies. The model can also serve as an inference engine for predicting unknown chromatin profiles — we demonstrated that with this approach we can leverage data from well-characterized cell types to help understand less-studied cell type or conditions.     

A One-Degree-of-Freedom Test for Supra-Multiplicativity of SNP Effects

Deviation from multiplicativity of genetic risk factors is biologically plausible and might explain why Genome-wide association studies (GWAS) so far could unravel only a portion of disease heritability. Still, evidence for SNP-SNP epistasis has rarely been reported, suggesting that 2-SNP models are overly simplistic. In this context, it was recently proposed that the genetic architecture of complex diseases could follow limiting pathway models. These models are defined by a critical risk allele load and imply multiple high-dimensional interactions. Here, we present a computationally efficient one-degree-of-freedom “supra-multiplicativity-test” (SMT) for SNP sets of size 2 to 500 that is designed to detect risk alleles whose joint effect is fortified when they occur together in the same individual. Via a simulation study we show that the SMT is powerful in the presence of threshold models, even when only about 30–45% of the model SNPs are available. In addition, we demonstrate that the SMT outperforms standard interaction analysis under recessive models involving just a few SNPs. We apply our test to 10 consensus Alzheimer’s disease (AD) susceptibility SNPs that were previously identified by GWAS and obtain evidence for supra-multiplicativity () that is not attributable to either two-way or three-way interaction.     

Modeling the Building Blocks of Biodiversity

Background

Networks of single interaction types, such as plant-pollinator mutualisms, are biodiversity’s “building blocks”. Yet, the structure of mutualistic and antagonistic networks differs, leaving no unified modeling framework across biodiversity’s component pieces.

Methods/Principal Findings

We use a one-dimensional “niche model” to predict antagonistic and mutualistic species interactions, finding that accuracy decreases with the size of the network. We show that properties of the modeled network structure closely approximate empirical properties even where individual interactions are poorly predicted. Further, some aspects of the structure of the niche space were consistently different between network classes.

Conclusions/Significance

These novel results reveal fundamental differences between the ability to predict ecologically important features of the overall structure of a network and the ability to predict pair-wise species interactions.     

Discovering sparse control strategies in neural activity

Biological circuits such as neural or gene regulation networks use internal states to map sensory input to an adaptive repertoire of behavior. Characterizing this mapping is a major challenge for systems biology. Though experiments that probe internal states are developing rapidly, organismal complexity presents a fundamental obstacle given the many possible ways internal states could map to behavior. Using C. elegans as an example, we propose a protocol for systematic perturbation of neural states that limits experimental complexity and could eventually help characterize collective aspects of the neural-behavioral map. We consider experimentally motivated small perturbations—ones that are most likely to preserve natural dynamics and are closer to internal control mechanisms—to neural states and their impact on collective neural activity. Then, we connect such perturbations to the local information geometry of collective statistics, which can be fully characterized using pairwise perturbations. Applying the protocol to a minimal model of C. elegans neural activity, we find that collective neural statistics are most sensitive to a few principal perturbative modes. Dominant eigenvalues decay initially as a power law, unveiling a hierarchy that arises from variation in individual neural activity and pairwise interactions. Highest-ranking modes tend to be dominated by a few, “pivotal” neurons that account for most of the system’s sensitivity, suggesting a sparse mechanism of collective control.     

In Silico Analysis of Antibiotic Resistance Genes in the Gut Microflora of Individuals from Diverse Geographies and Age-Groups

The spread of antibiotic resistance, originating from the rampant and unrestrictive use of antibiotics in humans and livestock over the past few decades has emerged as a global health problem. This problem has been further compounded by recent reports implicating the gut microbial communities to act as reservoirs of antibiotic resistance. We have profiled the presence of probable antibiotic resistance genes in the gut flora of 275 individuals from eight different nationalities. For this purpose, available metagenomic data sets corresponding to 275 gut microbiomes were analyzed. Sequence similarity searches of the genomic fragments constituting each of these metagenomes were performed against genes conferring resistance to around 240 antibiotics. Potential antibiotic resistance genes conferring resistance against 53 different antibiotics were detected in the human gut microflora analysed in this study. In addition to several geography/country-specific patterns, four distinct clusters of gut microbiomes, referred to as ‘Resistotypes’, exhibiting similarities in their antibiotic resistance profiles, were identified. Groups of antibiotics having similarities in their resistance patterns within each of these clusters were also detected. Apart from this, mobile multi-drug resistance gene operons were detected in certain gut microbiomes. The study highlighted an alarmingly high abundance of antibiotic resistance genes in two infant gut microbiomes. The results obtained in the present study presents a holistic ‘big picture’ on the spectra of antibiotic resistance within our gut microbiota across different geographies. Such insights may help in implementation of new regulations and stringency on the existing ones.     

Molecular interactions involved in proton-dependent gating in KcsA potassium channels

The bacterial potassium channel KcsA is gated open by the binding of protons to amino acids on the intracellular side of the channel. We have identified, via channel mutagenesis and x-ray crystallography, two pH-sensing amino acids and a set of nearby residues involved in molecular interactions that influence gating. We found that the minimal mutation of one histidine (H25) and one glutamate (E118) near the cytoplasmic gate completely abolished pH-dependent gating. Mutation of nearby residues either alone or in pairs altered the channel’s response to pH. In addition, mutations of certain pairs of residues dramatically increased the energy barriers between the closed and open states. We proposed a Monod–Wyman–Changeux model for proton binding and pH-dependent gating in KcsA, where H25 is a “strong” sensor displaying a large shift in pK_a between closed and open states, and E118 is a “weak” pH sensor. Modifying model parameters that are involved in either the intrinsic gating equilibrium or the pK_a values of the pH-sensing residues was sufficient to capture the effects of all mutations.     

Cycling Empirical Antibiotic Therapy in Hospitals: Meta-Analysis and Models

The rise of resistance together with the shortage of new broad-spectrum antibiotics underlines the urgency of optimizing the use of available drugs to minimize disease burden. Theoretical studies suggest that coordinating empirical usage of antibiotics in a hospital ward can contain the spread of resistance. However, theoretical and clinical studies came to different conclusions regarding the usefulness of rotating first-line therapy (cycling). Here, we performed a quantitative pathogen-specific meta-analysis of clinical studies comparing cycling to standard practice. We searched PubMed and Google Scholar and identified 46 clinical studies addressing the effect of cycling on nosocomial infections, of which 11 met our selection criteria. We employed a method for multivariate meta-analysis using incidence rates as endpoints and find that cycling reduced the incidence rate/1000 patient days of both total infections by 4.95 [9.43–0.48] and resistant infections by 7.2 [14.00–0.44]. This positive effect was observed in most pathogens despite a large variance between individual species. Our findings remain robust in uni- and multivariate metaregressions. We used theoretical models that reflect various infections and hospital settings to compare cycling to random assignment to different drugs (mixing). We make the realistic assumption that therapy is changed when first line treatment is ineffective, which we call “adjustable cycling/mixing”. In concordance with earlier theoretical studies, we find that in strict regimens, cycling is detrimental. However, in adjustable regimens single resistance is suppressed and cycling is successful in most settings. Both a meta-regression and our theoretical model indicate that “adjustable cycling” is especially useful to suppress emergence of multiple resistance. While our model predicts that cycling periods of one month perform well, we expect that too long cycling periods are detrimental. Our results suggest that “adjustable cycling” suppresses multiple resistance and warrants further investigations that allow comparing various diseases and hospital settings.     

Oncologic Trogocytosis of an Original Stromal Cells Induces Chemoresistance of Ovarian Tumours

Background

The microenvironment plays a major role in the onset and progression of metastasis. Epithelial ovarian cancer (EOC) tends to metastasize to the peritoneal cavity where interactions within the microenvironment might lead to chemoresistance. Mesothelial cells are important actors of the peritoneal homeostasis; we determined their role in the acquisition of chemoresistance of ovarian tumours.

Methodology/Principal Findings

We isolated an original type of stromal cells, referred to as “Hospicells” from ascitis of patients with ovarian carcinosis using limiting dilution. We studied their ability to confer chemoresistance through heterocellular interactions. These stromal cells displayed a new phenotype with positive immunostaining for CD9, CD10, CD29, CD146, CD166 and Multi drug resistance protein. They preferentially interacted with epithelial ovarian cancer cells. This interaction induced chemoresistance to platin and taxans with the implication of multi-drug resistance proteins. This contact enabled EOC cells to capture patches of the Hospicells membrane through oncologic trogocytosis, therefore acquiring their functional P-gp proteins and thus developing chemoresistance. Presence of Hospicells on ovarian cancer tissue micro-array from patients with neo-adjuvant chemotherapy was also significantly associated to chemoresistance.

Conclusions/Significance

This is the first report of trogocytosis occurring between a cancer cell and an original type of stromal cell. This interaction induced autonomous acquisition of chemoresistance. The presence of stromal cells within patient''s tumour might be predictive of chemoresistance. The specific interaction between cancer cells and stromal cells might be targeted during chemotherapy.