Dissecting neural pathways for forgetting in Drosophila olfactory aversive memory

Recent studies have identified molecular pathways driving forgetting and supported the notion that forgetting is a biologically active process. The circuit mechanisms of forgetting, however, remain largely unknown. Here we report two sets of Drosophila neurons that account for the rapid forgetting of early olfactory aversive memory. We show that inactivating these neurons inhibits memory decay without altering learning, whereas activating them promotes forgetting. These neurons, including a cluster of dopaminergic neurons (PAM-β′1) and a pair of glutamatergic neurons (MBON-γ4>γ1γ2), terminate in distinct subdomains in the mushroom body and represent parallel neural pathways for regulating forgetting. Interestingly, although activity of these neurons is required for memory decay over time, they are not required for acute forgetting during reversal learning. Our results thus not only establish the presence of multiple neural pathways for forgetting in Drosophila but also suggest the existence of diverse circuit mechanisms of forgetting in different contexts.Although forgetting commonly has a negative connotation, it is a functional process that shapes memory and cognition (1–4). Recent studies, including work in relatively simple invertebrate models, have started to reveal basic biological mechanisms underlying forgetting (5–15). In Drosophila, single-session Pavlovian conditioning by pairing an odor (conditioned stimulus, CS) with electric shock (unconditioned stimulus, US) induces aversive memories that are short-lasting (16). The memory performance of fruit flies is observed to drop to a negligible level within 24 h, decaying rapidly early after training and slowing down thereafter (17). Memory decay or forgetting requires the activation of the small G protein Rac, a signaling protein involved in actin remodeling, in the mushroom body (MB) intrinsic neurons (6). These so-called Kenyon cells (KCs) are the neurons that integrate CS–US information (18, 19) and support aversive memory formation and retrieval (20–22). In addition to Rac, forgetting also requires the DAMB dopamine receptor (7), which has highly enriched expression in the MB (23). Evidence suggests that the dopamine-mediated forgetting signal is conveyed to the MB by dopamine neurons (DANs) in the protocerebral posterior lateral 1 (PPL1) cluster (7, 24). Therefore, forgetting of olfactory aversive memory in Drosophila depends on a particular set of intracellular molecular pathways within KCs, involving Rac, DAMB, and possibly others (25), and also receives modulation from extrinsic neurons. Although important cellular evidence supporting the hypothesis that memory traces are erased under these circumstances is still lacking, these findings lend support to the notion that forgetting is an active, biologically regulated process (17, 26).Although existing studies point to the MB circuit as essential for forgetting, several questions remain to be answered. First, whereas the molecular pathways for learning and forgetting of olfactory aversive memory are distinct and separable (6, 7), the neural circuits seem to overlap. Rac-mediated forgetting has been localized to a large population of KCs (6), including the γ-subset, which is also critical for initial memory formation (21, 27). The site of action of DAMB for forgetting has yet to be established; however, the subgroups of PPL1-DANs implicated in forgetting are the same as those that signal aversive reinforcement and are required for learning (28–30). It leaves open the question of whether the brain circuitry underlying forgetting and learning is dissociable, or whether forgetting and learning share the same circuit but are driven by distinct activity patterns and molecular machinery (26). Second, shock reinforcement elicits multiple memory traces through at least three dopamine pathways to different subdomains in the MB lobes (28, 29). Functional imaging studies have also revealed Ca²⁺-based memory traces in different KC populations (31). It is poorly understood how forgetting of these memory traces differs, and it remains unknown whether there are multiple regulatory neural pathways. Notably, when PPL1-DANs are inactivated, forgetting still occurs, albeit at a lower rate (7). This incomplete block suggests the existence of an additional pathway(s) that conveys forgetting signals to the MB. Third, other than memory decay over time, forgetting is also observed through interference (32, 33), when new learning or reversal learning is introduced after training (6, 34, 35). Time-based and interference-based forgetting shares a similar dependence on Rac and DAMB (6, 7). However, it is not known whether distinct circuits underlie forgetting in these different contexts.In the current study, we focus on the diverse set of MB extrinsic neurons (MBENs) that interconnect the MB lobes with other brain regions, which include 34 MB output neurons (MBONs) of 21 types and ∼130 dopaminergic neurons of 20 types in the PPL1 and protocerebral anterior medial (PAM) clusters (36, 37). These neurons have been intensively studied in olfactory memory formation, consolidation, and retrieval in recent years (e.g., 24, 28–30, 38–48); however, their roles in forgetting have not been characterized except for the aforementioned PPL1-DANs. In a functional screen, we unexpectedly found that several Gal4 driver lines of MBENs showed significantly better 3-h memory retention when the Gal4-expressing cells were inactivated. The screen has thus led us to identify two types of MBENs that are not involved in initial learning but play important and additive roles in mediating memory decay. Furthermore, neither of these MBEN types is required for reversal learning, supporting the notion that there is a diversity of neural circuits that drive different forms of forgetting.     

Distinct dopamine neurons mediate reward signals for short- and long-term memories

Drosophila melanogaster can acquire a stable appetitive olfactory memory when the presentation of a sugar reward and an odor are paired. However, the neuronal mechanisms by which a single training induces long-term memory are poorly understood. Here we show that two distinct subsets of dopamine neurons in the fly brain signal reward for short-term (STM) and long-term memories (LTM). One subset induces memory that decays within several hours, whereas the other induces memory that gradually develops after training. They convey reward signals to spatially segregated synaptic domains of the mushroom body (MB), a potential site for convergence. Furthermore, we identified a single type of dopamine neuron that conveys the reward signal to restricted subdomains of the mushroom body lobes and induces long-term memory. Constant appetitive memory retention after a single training session thus comprises two memory components triggered by distinct dopamine neurons.Memory of a momentous event persists for a long time. Whereas some forms of long-term memory (LTM) require repetitive training (1–3), a highly relevant stimulus such as food or poison is sufficient to induce LTM in a single training session (4–7). Recent studies have revealed aspects of the molecular and cellular mechanisms of LTM formation induced by repetitive training (8–11), but how a single training induces a stable LTM is poorly understood (12).Appetitive olfactory learning in fruit flies is suited to address the question, as a presentation of a sugar reward paired with odor induces robust short-term memory (STM) and LTM (6, 7). Odor is represented by a sparse ensemble of the 2,000 intrinsic neurons, the Kenyon cells (13). A current working model suggests that concomitant reward signals from sugar ingestion cause associative plasticity in Kenyon cells that might underlie memory formation (14–20). A single activation session of a specific cluster of dopamine neurons (PAM neurons) by sugar ingestion can induce appetitive memory that is stable over 24 h (19), underscoring the importance of sugar reward to the fly.The mushroom body (MB) is composed of the three different cell types, α/β, α′/β′, and γ, which have distinct roles in different phases of appetitive memories (11, 21–25). Similar to midbrain dopamine neurons in mammals (26, 27), the structure and function of PAM cluster neurons are heterogeneous, and distinct dopamine neurons intersect unique segments of the MB lobes (19, 28–34). Further circuit dissection is thus crucial to identify candidate synapses that undergo associative modulation.By activating distinct subsets of PAM neurons for reward signaling, we found that short- and long-term memories are independently formed by two complementary subsets of PAM cluster dopamine neurons. Conditioning flies with nutritious and nonnutritious sugars revealed that the two subsets could represent different reinforcing properties: sweet taste and nutritional value of sugar. Constant appetitive memory retention after a single training session thus comprises two memory components triggered by distinct reward signals.     

From the Cover: PNAS Plus: Pan-neuronal imaging in roaming Caenorhabditis elegans

We present an imaging system for pan-neuronal recording in crawling Caenorhabditis elegans. A spinning disk confocal microscope, modified for automated tracking of the C. elegans head ganglia, simultaneously records the activity and position of ∼80 neurons that coexpress cytoplasmic calcium indicator GCaMP6s and nuclear localized red fluorescent protein at 10 volumes per second. We developed a behavioral analysis algorithm that maps the movements of the head ganglia to the animal’s posture and locomotion. Image registration and analysis software automatically assigns an index to each nucleus and calculates the corresponding calcium signal. Neurons with highly stereotyped positions can be associated with unique indexes and subsequently identified using an atlas of the worm nervous system. To test our system, we analyzed the brainwide activity patterns of moving worms subjected to thermosensory inputs. We demonstrate that our setup is able to uncover representations of sensory input and motor output of individual neurons from brainwide dynamics. Our imaging setup and analysis pipeline should facilitate mapping circuits for sensory to motor transformation in transparent behaving animals such as C. elegans and Drosophila larva.Understanding how brain dynamics creates behaviors requires quantifying the flow and transformation of sensory information to motor output in behaving animals. Optical imaging using genetically encoded calcium or voltage fluorescent probes offers a minimally invasive method to record neural activity in intact animals. The nematode Caenorhabditis elegans is particularly ideal for optical neurophysiology owing to its small size, optical transparency, compact nervous system, and ease of genetic manipulation. Imaging systems for tracking the activity of small numbers of neurons have been effective in determining their role during nematode locomotion and navigational behaviors like chemotaxis, thermotaxis, and the escape response (1–6). Recordings from large numbers of interconnected neurons are required to understand how neuronal ensembles carry out the systematic transformations of sensory input into motor patterns that build behavioral decisions.Several methods for fast 3D imaging of neural activity in a fixed imaging volume have been developed for different model organisms (7–14). High-speed light sheet microscopy, light field microscopy, multifocus microscopy, and two-photon structured illumination microscopy have proved effective for rapidly recording large numbers of neurons in immobilized, intact, transparent animals like larval zebrafish and nematodes (15–19). However, these methods are problematic when attempting to track many neurons within the bending and moving body of a behaving animal. Panneuronal recording in moving animals poses higher demands on spatial and temporal resolution. Furthermore, extracting neuronal signals from recordings in a behaving animal requires an effective analysis pipeline to segment image volumes into the activity patterns of discrete and identifiable neurons.Here, we use high-speed spinning disk confocal microscopy—modified for automated tracking using real-time image analysis and motion control software—to volumetrically image the head ganglia of behaving C. elegans adults at single-cell resolution. Our setup can simultaneously track ∼80 neurons with 0.45 × 0.45 × 2-μm resolution at 10 Hz. Activity was reported by the ultrasensitive calcium indicator GCaMP6s expressed throughout the cytosol under the control of the pan-neuronal rgef-1 promoter (a gift from D. Pilgrim, University of Alberta, Edmonton, Alberta, Canada) (20). To facilitate segmentation into individual identifiable neurons, nuclei were tracked using calcium-insensitive, nuclear-bound red fluorescent protein (RFP), TagRFP, under the control of another pan-neuronal rab-3 promoter (a gift from O. Hobert, Columbia University, New York) (21). We developed an image analysis pipeline that converts the gross movements of the head into the time-varying position and posture of the crawling worm, and converts fluorescence measurements into near simultaneous activity patterns of all imaged neurons.A similar approach to brainwide imaging in moving C. elegans using the same transgenic strain has recently been reported (22). Although both setups use customized spinning disk confocal microscopes, the strategies for tracking the moving neurons and analyzing behavioral and neural activity patterns are different. Nguyen et al. (22) use a low power objective to track the posture of the animal and a high power objective to locate and image the nerve ring. The advantage of our single objective setup is that it affords the flexibility, for example, to deliver thermosensory inputs using an opaque temperature controlled stage below the animal. The advantage of low-magnification imaging is that it provides a direct measurement of animal posture, which we must infer. These new technologies for pan-neuronal imaging in roaming animals now enables correlating brainwide dynamics to sensory inputs and motor outputs in transparent behaving animals like C. elegans and Drosophila larvae.     

Hyperpolarization-activated,cyclic nucleotide-gated cation channels in Aplysia: Contribution to classical conditioning

Hyperpolarization-activated, cyclic nucleotide-gated cation (HCN) channels are critical regulators of neuronal excitability, but less is known about their possible roles in synaptic plasticity and memory circuits. Here, we characterized the HCN gene organization, channel properties, distribution, and involvement in associative and nonassociative forms of learning in Aplysia californica. Aplysia has only one HCN gene, which codes for a channel that has many similarities to the mammalian HCN channel. The cloned acHCN gene was expressed in Xenopus oocytes, which displayed a hyperpolarization-induced inward current that was enhanced by cGMP as well as cAMP. Similarly to its homologs in other animals, acHCN is permeable to K⁺ and Na⁺ ions, and is selectively blocked by Cs⁺ and ZD7288. We found that acHCN is predominantly expressed in inter- and motor neurons, including LFS siphon motor neurons, and therefore tested whether HCN channels are involved in simple forms of learning of the siphon-withdrawal reflex in a semiintact preparation. ZD7288 (100 μM) significantly reduced an associative form of learning (classical conditioning) but had no effect on two nonassociative forms of learning (intermediate-term sensitization and unpaired training) or baseline responses. The HCN current is enhanced by nitric oxide (NO), which may explain the postsynaptic role of NO during conditioning. HCN current in turn enhances the NMDA-like current in the motor neurons, suggesting that HCN channels contribute to conditioning through this pathway.Hyperpolarization-activated, cyclic nucleotide-gated (HCN), cation nonselective ion channels generate hyperpolarization-activated inward currents (I_h) and thus tend to stabilize membrane potential (1–3). In addition, binding of cyclic nucleotides (cAMP and cGMP) to the C-terminal cyclic nucleotide binding domain (CNBD) enhances I_h and thus couples membrane excitability with intracellular signaling pathways (2, 4). HCN channels are widely important for numerous systemic functions such as hormonal regulation, heart contractility, epilepsy, pain, central pattern generation, sensory perception (4–15), and learning and memory (16–24).However, in previous studies it has been difficult to relate the cellular effects of HCN channels directly to their behavioral effects, because of the immense complexity of the mammalian brain. We have therefore investigated the role of HCN channels in Aplysia, which has a numerically simpler nervous system (25). We first identified and characterized an HCN gene in Aplysia, and showed that it codes for a channel that has many similarities to the mammalian HCN channel. We found that the Aplysia HCN channel is predominantly expressed in motor neurons including LFS neurons in the siphon withdrawal reflex circuit (26, 27). We therefore investigated simple forms of learning of that reflex in a semiintact preparation (28–30) and found that HCN current is involved in classical conditioning and enhances the NMDA-like current in the motor neurons. These results provide a direct connection between HCN channels and behavioral learning and suggest a postsynaptic mechanism of that effect. HCN current in turn is enhanced by nitric oxide (NO), a transmitter of facilitatory interneurons, and thus may contribute to the postsynaptic role of NO during conditioning.     

Neural substrate for higher-order learning in an insect: Mushroom bodies are necessary for configural discriminations

Learning theories distinguish elemental from configural learning based on their different complexity. Although the former relies on simple and unambiguous links between the learned events, the latter deals with ambiguous discriminations in which conjunctive representations of events are learned as being different from their elements. In mammals, configural learning is mediated by brain areas that are either dispensable or partially involved in elemental learning. We studied whether the insect brain follows the same principles and addressed this question in the honey bee, the only insect in which configural learning has been demonstrated. We used a combination of conditioning protocols, disruption of neural activity, and optophysiological recording of olfactory circuits in the bee brain to determine whether mushroom bodies (MBs), brain structures that are essential for memory storage and retrieval, are equally necessary for configural and elemental olfactory learning. We show that bees with anesthetized MBs distinguish odors and learn elemental olfactory discriminations but not configural ones, such as positive and negative patterning. Inhibition of GABAergic signaling in the MB calyces, but not in the lobes, impairs patterning discrimination, thus suggesting a requirement of GABAergic feedback neurons from the lobes to the calyces for nonelemental learning. These results uncover a previously unidentified role for MBs besides memory storage and retrieval: namely, their implication in the acquisition of ambiguous discrimination problems. Thus, in insects as in mammals, specific brain regions are recruited when the ambiguity of learning tasks increases, a fact that reveals similarities in the neural processes underlying the elucidation of ambiguous tasks across species.Learning can be categorized into two levels of complexity termed elemental and configural (nonelemental) (1–3). Simple and unambiguous links between events characterize elemental learning (4). By contrast, ambiguity and nonlinearity characterize configural learning, where associations involve conjunctions of elemental stimuli, which may have different, contradictory outcomes. As a consequence, solving configural tasks typically requires treating stimulus conjunctions as being different from the simple sum of their elemental components (5–8). For example, in a negative patterning task (9–11), subjects have to discriminate a nonreinforced conjunction of two elements A and B from its reinforced elements (i.e., AB– vs. A+ and B+), which requires treating AB as being different from the simple sum of A and B (12, 13). The ambiguity of the task lies in the fact that each element (A and B) is as often reinforced (when presented alone) as nonreinforced (when presented as a compound). In mammals, different brain structures have been associated with these two learning forms: Whereas the hippocampus seems to be dispensable for learning elemental associations (6, 8), it is required for fast formation of conjunctive representations during learning tasks, such as spatial learning or contextual fear conditioning (6, 8, 10, 14–19). Moreover, the cortical system is necessary to form configural representations over extended training, thus supporting the learning of nonlinear discriminations,Here, we ask whether the specialization of different brain centers for learning tasks of different complexity is a property that can be extended to an insect brain. Insects offer the possibility of studying sophisticated behaviors and simultaneously accessing the neural bases of these behaviors (20). Several studies have shown that insects, in particular the honey bee Apis mellifera, possess higher-order cognitive abilities (5, 21), which raises the question of which neural mechanisms support these capacities in a brain whose size is only 1 mm³ (22).The mushroom bodies (MBs) are paired structures in the insect brain that have been historically associated with olfactory learning and memory. Their function has been extensively studied in a variety of elemental learning protocols, mainly in the honey bee and the fruit fly Drosophila melanogaster (23–29). In both species, MBs play a fundamental role for the encoding, storing, and retrieval of appetitive and aversive elemental memories, but no study has clearly established their role for nonelemental learning and memory (30). In fruit flies, this missing information may be due to the incapacity of these insects to solve nonelemental problems, such as negative patterning (31). By contrast, honey bees exhibit elaborated nonelemental learning abilities (32–36), which have been suggested to require intact MB function (5).Here, we used a combination of nonelemental conditioning protocols, disruption of MB function, and optophysiological recordings of neural activity to determine whether MBs are necessary for nonelemental forms of learning. Our results show that acquisition of olfactory patterning discriminations is impaired in bees in which neural activity in the MBs was blocked by procaine injection (37, 38), but not in control animals injected with saline solution. By contrast, MB blockade by procaine affected neither olfactory processing upstream of the MBs nor elemental olfactory discriminations. To uncover the neural mechanisms underlying the necessity of MBs for patterning discriminations, we focused on GABAergic feedback neurons (39), which provide inhibitory feedback to the MBs of the bee (40–43). We blocked GABAergic signaling by locally injecting picrotoxin (PTX), a GABA antagonist, into the MB calyces or into the MB lobes. We show that GABAergic feedback to the calyces—but not to the lobes—is required for patterning discriminations. These results uncover a previously unidentified role for MBs: namely, the disambiguation between elemental and conjunctive odor representations, thus supporting the learning of nonlinear discriminations.     

From the Cover: PNAS Plus: Single-molecule motions and interactions in live cells reveal target search dynamics in mismatch repair

MutS is responsible for initiating the correction of DNA replication errors. To understand how MutS searches for and identifies rare base-pair mismatches, we characterized the dynamic movement of MutS and the replisome in real time using superresolution microscopy and single-molecule tracking in living cells. We report that MutS dynamics are heterogeneous in cells, with one MutS population exploring the nucleoid rapidly, while another MutS population moves to and transiently dwells at the replisome region, even in the absence of appreciable mismatch formation. Analysis of MutS motion shows that the speed of MutS is correlated with its separation distance from the replisome and that MutS motion slows when it enters the replisome region. We also show that mismatch detection increases MutS speed, supporting the model for MutS sliding clamp formation after mismatch recognition. Using variants of MutS and the replication processivity clamp to impair mismatch repair, we find that MutS dynamically moves to and from the replisome before mismatch binding to scan for errors. Furthermore, a block to DNA synthesis shows that MutS is only capable of binding mismatches near the replisome. It is well-established that MutS engages in an ATPase cycle, which is necessary for signaling downstream events. We show that a variant of MutS with a nucleotide binding defect is no longer capable of dynamic movement to and from the replisome, showing that proper nucleotide binding is critical for MutS to localize to the replisome in vivo. Our results provide mechanistic insight into the trafficking and movement of MutS in live cells as it searches for mismatches.DNA mismatch repair (MMR) is the highly conserved process responsible for correcting DNA replication errors (1). Although replication errors occur infrequently in bacteria (∼1 error per 31 million bp) (2), the consequences of MMR failure on human health are severe (3). MutS is the first protein involved in the MMR pathway, and it is responsible for detecting rare base-pairing errors. In Bacillus subtilis, MutS then recruits MutL, an endonuclease in most bacteria and eukaryotic organisms, to incise the DNA (4, 5). After MutL incision, the error-containing strand is removed, and the DNA is resynthesized to complete error correction (6).The mechanism by which MutS homologs locate a single mismatch among millions of correctly paired nucleotides has been studied extensively by bulk biochemistry, in vitro atomic force microscopy and single-molecule imaging, and visualizing MutS using in vivo cell biology approaches (7–20). In vitro single-molecule studies largely indicate that MutS operates as a searching clamp diffusing along DNA in a 1D search process (8, 10, 13). In this model, after mismatch recognition, MutS can dwell at the mismatch before exchanging ADP for ATP, converting into a stable ATP-bound sliding clamp, and subsequently, diffusing away from the mismatch at a faster rate in search of MutL and possible strand discrimination signals (8, 10, 15). Studying the mechanistic steps of the search process with DNA curtains has provided evidence that MutS may identify errors through a combination of 1D sliding and a 3D pathway (15). These results show that MutS predominantly searches DNA by 1D diffusion but is also capable of mismatch recognition through a 3D pathway, in which MutS binds the mismatch without engaging in a prior 1D search. The 3D pathway for mismatch recognition has been proposed to help MutS circumvent protein barriers that would exist in vivo. Additional analysis of MutS on DNA curtains suggests that DNA binding proteins largely prevent 1D diffusion by MutS homologs (14). These results indicate that a search mechanism limited to 1D diffusion would present a significant challenge for identifying replication errors in vivo. Therefore, it is still not clear how MutS searches for errors within the supercoiled, compacted, and protein-laden nucleoid in the crowded 3D environment of a living bacterial cell.Each of these in vitro single-molecule studies has provided important insight into the target search dynamics of MutS, and these studies and other bulk studies have elucidated many of the mechanistic steps leading to mismatch identification on naked DNA in vitro. Of course, most of these studies have been limited to analysis of MutS in isolation on protein-free DNA; they often do not incorporate replication proteins or other in vivo obstacles that are likely to impact the in vivo search process. The search dynamics of MutS have not been investigated at single-molecule resolution for any MutS homolog in live cells.The need to understand how the biochemistry of MutS homologs translates into mismatch identification in live cells has led to several studies across many organisms using a variety of approaches, including fluorescent protein fusions and bulk fluorescence imaging. These studies have shown that MutS homologs form foci that colocalize with the DNA replication machinery (replisome) and that interaction of MutS with replication processivity clamps [proliferating cell nuclear antigen (PCNA) or β-clamp] is important for focus formation (16, 17, 19–22). In Saccharomyces cerevisiae, MMR is also temporally coupled to DNA replication (23). However, the colocalization between bacterial MutS and the replisome is far from absolute. For example, in the absence of an exogenous mutagen, bulk microscopy detects fluorescent MutS foci in only ∼9% of B. subtilis cells. Furthermore, only about one-half of these 9% of cells with a MutS focus show colocalization of MutS to the replisome (20). Therefore, only about 4.5% of cells under normal growth conditions in B. subtilis show colocalization between MutS and the replisome (20). Obviously, the behavior of MutS in the remaining 91% of cells without a MutS focus has been impossible to determine using ensemble fluorescence techniques.Most of the MutS pool cannot be visualized by bulk microscopy, because the large fluorescent foci that are visible by ensemble fluorescence imaging account for only 10% of cellular MutS (that is, ∼8 of 80 dimers present in a cell) (24). Therefore, prior studies have been unable to image and quantify the behavior of MutS outside of large static foci visible during ensemble fluorescence microscopy. Therefore, to understand MutS location and dynamics in vivo, only a more sensitive method with higher spatial and temporal resolutions can unambiguously determine the extent to which MutS is enriched near the replisome and answer two open questions. Is MutS enrichment at the replisome constitutive in bacteria? Does MutS search for mismatches away from the replisome in vivo?To further understand the process of MMR inside cells and gain insight into the movement and location of MutS molecules in vivo throughout MMR, we probed the spatial distribution, dynamics, and genomic distribution of MutS in live bacteria with superresolution fluorescence imaging (25–27), single-molecule tracking (28–30), and ChIP-sequencing (ChIP-seq) approaches. We investigated the effect of replisome association, active DNA synthesis, mismatch binding, nucleotide binding, and presence of MutL on the dynamics of MutS movement in vivo. This study applies single-molecule imaging to a dedicated DNA repair pathway and captures the trafficking behavior of MutS, a protein that is conserved from bacteria to humans. A high degree of spatial and temporal resolution has allowed us to view the dynamic movement and search process of single MutS molecules during MMR in vivo.     

Convergent evolution toward an improved growth rate and a reduced resistance range in Prochlorococcus strains resistant to phage

Prochlorococcus is an abundant marine cyanobacterium that grows rapidly in the environment and contributes significantly to global primary production. This cyanobacterium coexists with many cyanophages in the oceans, likely aided by resistance to numerous co-occurring phages. Spontaneous resistance occurs frequently in Prochlorococcus and is often accompanied by a pleiotropic fitness cost manifested as either a reduced growth rate or enhanced infection by other phages. Here, we assessed the fate of a number of phage-resistant Prochlorococcus strains, focusing on those with a high fitness cost. We found that phage-resistant strains continued evolving toward an improved growth rate and a narrower resistance range, resulting in lineages with phenotypes intermediate between those of ancestral susceptible wild-type and initial resistant substrains. Changes in growth rate and resistance range often occurred in independent events, leading to a decoupling of the selection pressures acting on these phenotypes. These changes were largely the result of additional, compensatory mutations in noncore genes located in genomic islands, although genetic reversions were also observed. Additionally, a mutator strain was identified. The similarity of the evolutionary pathway followed by multiple independent resistant cultures and clones suggests they undergo a predictable evolutionary pathway. This process serves to increase both genetic diversity and infection permutations in Prochlorococcus populations, further augmenting the complexity of the interaction network between Prochlorococcus and its phages in nature. Last, our findings provide an explanation for the apparent paradox of a multitude of resistant Prochlorococcus cells in nature that are growing close to their maximal intrinsic growth rates.Large bacterial populations are present in the oceans, playing important roles in primary production and the biogeochemical cycling of matter. These bacterial communities are highly diverse (1–4) yet form stable and reproducible bacterial assemblages under similar environmental conditions (5–7).These bacteria are present together with high abundances of viruses (phages) that have the potential to infect and kill them (8–11). Although studied only rarely in marine organisms (12–16), this coexistence is likely to be the result of millions of years of coevolution between these antagonistic interacting partners, as has been well documented for other systems (17–20). From the perspective of the bacteria, survival entails the selection of cells that are resistant to infection, preventing viral production and enabling the continuation of the cell lineage. Resistance mechanisms include passively acquired spontaneous mutations in cell surface molecules that prevent phage entry into the cell and other mechanisms that actively terminate phage infection intracellularly, such as restriction–modification systems and acquired resistance by CRISPR-Cas systems (21, 22). Mutations in the phage can also occur that circumvent these host defenses and enable the phage to infect the recently emerged resistant bacterium (23).Acquisition of resistance by bacteria is often associated with a fitness cost. This cost is frequently, but not always, manifested as a reduction in growth rate (24–27). Recently, an additional type of cost of resistance was identified, that of enhanced infection whereby resistance to one phage leads to greater susceptibility to other phages (14, 15, 28).Over the years, a number of models have been developed to explain coexistence in terms of the above coevolutionary processes and their costs (16, 29–32). In the arms race model, repeated cycles of host mutation and virus countermutation occur, leading to increasing breadths of host resistance and viral infectivity. However, experimental evidence generally indicates that such directional arms race dynamics do not continue indefinitely (25, 33, 34). Therefore, models of negative density-dependent fluctuations due to selective trade-offs, such as kill-the-winner, are often invoked (20, 33, 35, 36). In these models, fluctuations are generally considered to occur between rapidly growing competition specialists that are susceptible to infection and more slowly growing resistant strains that are considered defense specialists. Such negative density-dependent fluctuations are also likely to occur between strains that have differences in viral susceptibility ranges, such as those that would result from enhanced infection (30).The above coevolutionary processes are considered to be among the major mechanisms that have led to and maintain diversity within bacterial communities (32, 35, 37–39). These processes also influence genetic microdiversity within populations of closely related bacteria. This is especially the case for cell surface-related genes that are often localized to genomic islands (14, 40, 41), regions of high gene content, and gene sequence variability among members of a population. As such, populations in nature display an enormous degree of microdiversity in phage susceptibility regions, potentially leading to an assortment of subpopulations with different ranges of susceptibility to coexisting phages (4, 14, 30, 40).Prochlorococcus is a unicellular cyanobacterium that is the numerically dominant photosynthetic organism in vast oligotrophic expanses of the open oceans, where it contributes significantly to primary production (42, 43). Prochlorococcus consists of a number of distinct ecotypes (44–46) that form stable and reproducible population structures (7). These populations coexist in the oceans with tailed double-stranded DNA phage populations that infect them (47–49).Previously, we found that resistance to phage infection occurs frequently in two high-light–adapted Prochlorococcus ecotypes through spontaneous mutations in cell surface-related genes (14). These genes are primarily localized to genomic island 4 (ISL4) that displays a high degree of genetic diversity in environmental populations (14, 40). Although about a third of Prochlorococcus-resistant strains had no detectable associated cost, the others came with a cost manifested as either a slower growth rate or enhanced infection by other phages (14). In nature, Prochlorococcus seems to be growing close to its intrinsic maximal growth rate (50–52). This raises the question as to the fate of emergent resistant Prochlorococcus lineages in the environment, especially when resistance is accompanied with a high growth rate fitness cost.To begin addressing this question, we investigated the phenotype of Prochlorococcus strains with time after the acquisition of resistance. We found that resistant strains evolved toward an improved growth rate and a reduced resistance range. Whole-genome sequencing and PCR screening of many of these strains revealed that these phenotypic changes were largely due to additional, compensatory mutations, leading to increased genetic diversity. These findings suggest that the oceans are populated with rapidly growing Prochlorococcus cells with varying degrees of resistance and provide an explanation for how a multitude of presumably resistant Prochlorococcus cells are growing close to their maximal known growth rate in nature.     

Major diversification of voltage-gated K+ channels occurred in ancestral parahoxozoans

We examined the origins and functional evolution of the Shaker and KCNQ families of voltage-gated K⁺ channels to better understand how neuronal excitability evolved. In bilaterians, the Shaker family consists of four functionally distinct gene families (Shaker, Shab, Shal, and Shaw) that share a subunit structure consisting of a voltage-gated K⁺ channel motif coupled to a cytoplasmic domain that mediates subfamily-exclusive assembly (T1). We traced the origin of this unique Shaker subunit structure to a common ancestor of ctenophores and parahoxozoans (cnidarians, bilaterians, and placozoans). Thus, the Shaker family is metazoan specific but is likely to have evolved in a basal metazoan. Phylogenetic analysis suggested that the Shaker subfamily could predate the divergence of ctenophores and parahoxozoans, but that the Shab, Shal, and Shaw subfamilies are parahoxozoan specific. In support of this, putative ctenophore Shaker subfamily channel subunits coassembled with cnidarian and mouse Shaker subunits, but not with cnidarian Shab, Shal, or Shaw subunits. The KCNQ family, which has a distinct subunit structure, also appears solely within the parahoxozoan lineage. Functional analysis indicated that the characteristic properties of Shaker, Shab, Shal, Shaw, and KCNQ currents evolved before the divergence of cnidarians and bilaterians. These results show that a major diversification of voltage-gated K⁺ channels occurred in ancestral parahoxozoans and imply that many fundamental mechanisms for the regulation of action potential propagation evolved at this time. Our results further suggest that there are likely to be substantial differences in the regulation of neuronal excitability between ctenophores and parahoxozoans.Voltage-gated K⁺ channels are highly conserved among bilaterian metazoans and play a central role in the regulation of excitation in neurons and muscle. Understanding the functional evolution of these channels may therefore provide important insights into how neuromuscular excitation evolved within the Metazoa. Three major gene families, Shaker, KCNQ, and Ether-a-go-go (EAG) encode all voltage-gated K⁺ channels in bilaterians (1, 2). In this study, we examine the functional evolution and origins of the Shaker and KCNQ gene families. Shaker family channels can be definitively identified by a unique subunit structure that includes both a voltage-gated K⁺ channel core and a family-specific cytoplasmic domain within the N terminus known as the T1 domain. T1 mediates assembly of Shaker family subunits into functional tetrameric channels (3, 4). KCNQ channels are also tetrameric but lack a T1 domain and use a distinct coiled-coil assembly domain in the C terminus (5, 6). KCNQ channels can be identified by the presence of this family-specific assembly motif and high amino acid conservation within the K⁺ channel core. Both channel families are found in cnidarians (1, 7) and thus predate the divergence of cnidarians and bilaterians, but their ultimate evolutionary origins have not yet been defined.Shaker family K⁺ channels serve diverse roles in the regulation of neuronal firing and can be divided into four gene subfamilies based on function and sequence homology: Shaker, Shab, Shal, and Shaw (8, 9). The T1 assembly domain is only compatible between subunits from the same gene subfamily (4, 10) and thus serves to keep the subfamilies functionally segregated. Shaker subfamily channels activate rapidly near action potential threshold and range from rapidly inactivating to noninactivating. Multiple roles for Shaker channels in neurons and muscles have been described, but their most unique and fundamental role may be that of axonal action potential repolarization. Shaker channels are clustered to the axon initial segment and nodes of Ranvier in vertebrate neurons (11–13) and underlie the delayed rectifier in squid giant axons (14). The Shaker subfamily is diverse in cnidarians (15, 16), and the starlet sea anemone Nematostella vectensis has functional orthologs of most identified Shaker current types observed in bilaterians (16).The Shab and Shal gene subfamilies encode somatodendritic delayed rectifiers and A currents, respectively (17–20). Shab channels are important for maintaining sustained firing (21, 22), whereas the Kv4-based A current modulates spike threshold and frequency (17). Shab and Shal channels are present in cnidarians, but cnidarian Shab channels have not been functionally characterized, and the only cnidarian Shal channels expressed to date display atypical voltage dependence and kinetics compared with bilaterian channels (23). Shaw channels are rapid, high-threshold channels specialized for sustaining fast firing in vertebrates (24, 25) but have a low activation threshold and may contribute to resting potential in Drosophila (19, 26, 27). A Caenorhabditis elegans Shaw has slow kinetics but a high activation threshold (28), and a single expressed cnidarian Shaw channel has the opposite: a low activation threshold but relatively fast kinetics (29). Thus, the ancestral properties and function of Shaw channels is not yet understood. Further functional characterization of cnidarian Shab, Shal, and Shaw channels would provide a better understanding of the evolutionary status of the Shaker family in early parahoxozoans.KCNQ family channels underlie the M current in vertebrate neurons (30) that regulates subthreshold excitability (31). The M current provides a fundamental mechanism for regulation of firing threshold through the Gq G-protein pathway because KCNQ channels require phosphatidylinositol 4,5-bisphosphate (PIP₂) for activation (32, 33). PIP₂ hydrolysis and subsequent KCNQ channel closure initiated by Gq-coupled receptors produces slow excitatory postsynaptic potentials, during which the probability of firing is greatly increased (32, 33). The key functional adaptations of KCNQ channels for this physiological role that can be observed in vitro are (i) a requirement for PIP₂ to couple voltage-sensor activation to pore opening (34, 35), and (ii) a hyperpolarized voltage–activation curve that allows channels to open below typical action potential thresholds. Both key features are found in vertebrate (30, 34, 36–38), Drosophila (39), and C. elegans (40) KCNQ channels, suggesting they may have been present in KCNQ channels in a bilaterian ancestor. Evolution of the M current likely represented a major advance in the ability to modulate the activity of neuronal circuits, but it is not yet clear when PIP₂-dependent KCNQ channels first evolved.Here, we examine the origins and functional evolution of the Shaker and KCNQ gene families. If we assume the evolution of neuronal signaling provided a major selective pressure for the functional diversification of voltage-gated K⁺ channels, then we can hypothesize that the appearance of these gene families might accompany the emergence of the first nervous systems or a major event in nervous system evolution. Recent phylogenies that place the divergence of ctenophores near the root of the metazoan tree suggest that the first nervous systems, or at least the capacity to make neurons, may have been present in a basal metazoan ancestor (41–43) (Fig. S1). One hypothesis then is that much of the diversity of metazoan voltage-gated channels should be shared between ctenophores and parahoxozoans [cnidarians, bilaterians, and placozoans (44)]. However, genome analysis indicates that many “typical” neuronal genes are missing in ctenophores and the sponges lack a nervous system, leading to the suggestion that extant nervous systems may have evolved independently in ctenophores and parahoxozoans (42, 45). Thus, a second hypothesis is that important steps in voltage-gated K⁺ channel evolution might have occurred separately in ctenophores and parahoxozoans. We tested these hypotheses by carefully examining the phylogenetic distribution and functional evolution of Shaker and KCNQ family K⁺ channels. Our results support a model in which major innovations in neuromuscular excitability occurred specifically within the parahoxozoan lineage.     

Mechanisms of NDV-3 vaccine efficacy in MRSA skin versus invasive infection

Increasing rates of life-threatening infections and decreasing susceptibility to antibiotics urge development of an effective vaccine targeting Staphylococcus aureus. This study evaluated the efficacy and immunologic mechanisms of a vaccine containing a recombinant glycoprotein antigen (NDV-3) in mouse skin and skin structure infection (SSSI) due to methicillin-resistant S. aureus (MRSA). Compared with adjuvant alone, NDV-3 reduced abscess progression, severity, and MRSA density in skin, as well as hematogenous dissemination to kidney. NDV-3 induced increases in CD3+ T-cell and neutrophil infiltration and IL-17A, IL-22, and host defense peptide expression in local settings of SSSI abscesses. Vaccine induction of IL-22 was necessary for protective mitigation of cutaneous infection. By comparison, protection against hematogenous dissemination required the induction of IL-17A and IL-22 by NDV-3. These findings demonstrate that NDV-3 protective efficacy against MRSA in SSSI involves a robust and complementary response integrating innate and adaptive immune mechanisms. These results support further evaluation of the NDV-3 vaccine to address disease due to S. aureus in humans.The bacterium Staphylococcus aureus is the leading cause of skin and skin structure infections (SSSIs), including cellulitis, furunculosis, and folliculitis (1–4), and a common etiologic agent of impetigo (5), erysipelas (6), and superinfection in atopic dermatitis (7). This bacterium is a significant cause of surgical or traumatic wound infections (8, 9), as well as decuibitus and diabetic skin lesions (10). Moreover, SSSI is an important risk factor for systemic infection. The skin is a key portal of entry for hematogenous dissemination, particularly in association with i.v. catheters. S. aureus is now the second most common bloodstream isolate in healthcare settings (11), and SSSI is a frequent source of invasive infections such as pneumonia or endocarditis (12, 13). Despite a recent modest decline in rates of methicillin-resistant S. aureus (MRSA) infection in some cohorts (13), infections due to S. aureus remain a significant problem (14, 15). Even with appropriate therapy, up to one-third of patients diagnosed with S. aureus bacteremia succumb—accounting for more attributable annual deaths than HIV, tuberculosis, and viral hepatitis combined (16).The empiric use of antibiotics in healthcare-associated and community-acquired settings has increased S. aureus exposure to these agents, accelerating selection of resistant strains. As a result, resistance to even the most recently developed agents is emerging at an alarming pace (17, 18). The impact of this trend is of special concern in light of high rates of mortality associated with invasive MRSA infection (e.g., 15–40% in bacteremia or endocarditis), even with the most recently developed antistaphylococcal therapeutics (19, 20). Moreover, patients who experience SSSI due to MRSA exhibit high 1-y recurrence rates, often prompting surgical debridement (21) and protracted antibiotic treatment.Infections due to MRSA are a special concern in immune-vulnerable populations, including hemodialysis (22), neutropenic (23, 24), transplantation (25), and otherwise immunosuppressed patients (26, 27), and in patients with inherited immune dysfunctions (28–31) or cystic fibrosis (32). Patients having deficient interleukin 17 (IL-17) or IL-22 responses (e.g., signal transduction mediators STAT3, DOCK8, or CARD9 deficiencies) exhibit chronic or “cold” abscesses, despite high densities of pathogens such as S. aureus (33, 34). For example, patients with Chronic Granulomatous Disease (CGD; deficient Th1 and oxidative burst response) have increased risk of disseminated S. aureus infection. In contrast, patients with Job’s Syndrome (deficient Th17 response) typically have increased risk to SSSI and lung infections, but less so for systemic S. aureus bacteremia (35, 36). This pattern contrasts that observed in neutropenic or CGD patients (37). These themes suggest efficacious host defenses against MRSA skin and invasive infections involve complementary but distinct molecular and cellular immune responses.From these perspectives, vaccines or immunotherapeutics that prevent or lessen severity of MRSA infections, or that enhance antibiotic efficacy, would be significant advances in patient care and public health. However, to date, there are no licensed prophylactic or therapeutic vaccine immunotherapies for S. aureus or MRSA infection. Unfortunately, efforts to develop vaccines targeting S. aureus capsular polysaccharide type 5 or 8 conjugates, or the iron-regulated surface determinant B protein, have not been successful thus far (38, 39). Likewise, passive immunization using monoclonal antibodies targeting the S. aureus adhesin clumping factor A (ClfA, tefibazumab) (40) or lipoteichoic acid (pagibaximab) (41) have not shown efficacy against invasive infections in human clinical studies to date. Moreover, the striking recurrence rates of SSSI due to MRSA imply that natural exposure does not induce optimal preventive immunity or durable anamnestic response to infection or reinfection. Thus, significant challenges exist in the development of an efficacious vaccine targeting diseases caused by S. aureus (42) that are perhaps not optimally addressed by conventional approaches.The NDV-3 vaccine reflects a new strategy to induce durable immunity targeting S. aureus. Its immunogen is engineered from the agglutinin-like sequence 3 (Als3) adhesin/invasin of Candida albicans, which we discovered to be a structural homolog of S. aureus adhesins (43). NDV-3 is believed to cross-protect against S. aureus and C. albicans due to sequence (T-cell) and conformational (B-cell) epitopes paralleled in both organisms (44). Our prior data have shown that NDV-3 is efficacious in murine models of hematogenous and mucosal candidiasis (45), as well as S. aureus bacteremia (46–48). Recently completed phase I clinical trials demonstrate the safety, tolerability, and immunogenicity of NDV-3 in humans (49).     

Novel serologic biomarkers provide accurate estimates of recent Plasmodium falciparum exposure for individuals and communities

Tools to reliably measure Plasmodium falciparum (Pf) exposure in individuals and communities are needed to guide and evaluate malaria control interventions. Serologic assays can potentially produce precise exposure estimates at low cost; however, current approaches based on responses to a few characterized antigens are not designed to estimate exposure in individuals. Pf-specific antibody responses differ by antigen, suggesting that selection of antigens with defined kinetic profiles will improve estimates of Pf exposure. To identify novel serologic biomarkers of malaria exposure, we evaluated responses to 856 Pf antigens by protein microarray in 186 Ugandan children, for whom detailed Pf exposure data were available. Using data-adaptive statistical methods, we identified combinations of antibody responses that maximized information on an individual’s recent exposure. Responses to three novel Pf antigens accurately classified whether an individual had been infected within the last 30, 90, or 365 d (cross-validated area under the curve = 0.86–0.93), whereas responses to six antigens accurately estimated an individual’s malaria incidence in the prior year. Cross-validated incidence predictions for individuals in different communities provided accurate stratification of exposure between populations and suggest that precise estimates of community exposure can be obtained from sampling a small subset of that community. In addition, serologic incidence predictions from cross-sectional samples characterized heterogeneity within a community similarly to 1 y of continuous passive surveillance. Development of simple ELISA-based assays derived from the successful selection strategy outlined here offers the potential to generate rich epidemiologic surveillance data that will be widely accessible to malaria control programs.Many countries have extensive programs to reduce the burden of Plasmodium falciparum (Pf), the parasite responsible for most malaria morbidity and mortality (1). Effectively using limited resources for malaria control or elimination and evaluating interventions require accurate measurements of the risk of being infected with Pf (2–15). To reflect the rate at which individuals are infected with Pf in a useful way, metrics used to estimate exposure in a community need to account for dynamic changes over space and time, especially in response to control interventions (16–18).A variety of metrics can be used to estimate Pf exposure, but tools that are more precise and low cost are needed for population surveillance. Existing metrics have varying intrinsic levels of precision and accuracy and are subject to a variety of extrinsic factors, such as cost, time, and availability of trained personnel (19). For example, entomological measurements provide information on mosquito to human transmission for a community but are expensive, require specially trained staff, and lack standardized procedures, all of which reduce precision and/or make interpretation difficult (19–22). Parasite prevalence can be measured by detecting parasites in the blood of individuals from a cross-sectional sample of a community and is, therefore, relatively simple and inexpensive to perform, but results may be imprecise, especially in areas of low transmission (19, 23), and biased by a number of factors, including immunity and access to antimalarial treatment (5, 6, 19, 23–25). The burden of symptomatic disease in a community can be estimated from routine health systems data; however, such data are frequently unreliable (5, 26–28) and generally underestimate the prevalence of Pf infection in areas of intense transmission. Precise and quantitative information about exposure at an individual level can be reliably obtained from cohort studies by measuring the incidence of asymptomatic and/or symptomatic Pf infection (i.e., by measuring the molecular force of infection) (29–35). Unfortunately, the expense of cohort studies limits their use to research settings. The end result is that most malaria-endemic regions lack reliable, timely data on Pf exposure, limiting the capabilities of malaria control programs to guide and evaluate interventions.Serologic assays offer the potential to provide incidence estimates for symptomatic and asymptomatic Pf infection, which are currently obtained from cohort studies, at the cost of cross-sectional studies (36–38). Although Pf infections are transient, a record of infection remains detectable in an individual’s antibody profile. Thus, appropriately chosen antibody measurements integrated with age can provide information about an individual’s exposure history. Antibodies can be measured by simple ELISAs and obtained from dried blood spots, which are easy to collect, transport, and store (39–41). Serologic responses to Pf antigens have been explored as potential epidemiological tools (42–45), and estimated rates of seroconversion to well-characterized Pf antigens accurately reflect stable rates of exposure in a community, whereas distinct changes in these rates are obtained from successful interventions (22, 39, 41, 46–53). However, current serologic assays are not designed to detect short-term or gradual changes in Pf exposure or measure exposure to infection at an individual level. The ability to calibrate antibody responses to estimates of exposure in individuals could allow for more flexible sampling of a population (e.g., not requiring age stratification), improve accuracy of exposure estimates from small sample sizes, and better characterize heterogeneity in exposure within a community.Different Pf antigens elicit antibody responses with different magnitudes and kinetics, providing a large and diverse set of potential biomarkers for exposure (38, 54–58). We hypothesized that new and more highly informative serologic biomarkers better able to characterize an individual’s recent exposure history could be identified by analyzing antibody responses to a large number of candidate Pf antigens in participants with well-characterized exposure histories. To test this hypothesis, we probed plasma from participants in two cohort studies in Uganda against a protein microarray containing 856 Pf antigens. The primary aim of this analysis was to identify responses to select antigens that were most informative of recent exposure using robust, data-adaptive statistical methods. Each participant’s responses to these selected antigens were used as predictors for two primary outcomes of their recent exposure to Pf: (i) days since last Pf infection and (ii) the incidence of symptomatic malaria in the last year. These individual-level estimates were then aggregated across a population to assess community-level malaria exposure. The selection strategy presented here identified accurate biomarkers of exposure for children living in areas of moderate to high Pf exposure and illustrates the utility of this flexible and broadly applicable approach.     

The pathogen Batrachochytrium dendrobatidis disturbs the frog skin microbiome during a natural epidemic and experimental infection

Symbiotic microbial communities may interact with infectious pathogens sharing a common host. The microbiome may limit pathogen infection or, conversely, an invading pathogen can disturb the microbiome. Documentation of such relationships during naturally occurring disease outbreaks is rare, and identifying causal links from field observations is difficult. This study documented the effects of an amphibian skin pathogen of global conservation concern [the chytrid fungus Batrachochytrium dendrobatidis (Bd)] on the skin-associated bacterial microbiome of the endangered frog, Rana sierrae, using a combination of population surveys and laboratory experiments. We examined covariation of pathogen infection and bacterial microbiome composition in wild frogs, demonstrating a strong and consistent correlation between Bd infection load and bacterial community composition in multiple R. sierrae populations. Despite the correlation between Bd infection load and bacterial community composition, we observed 100% mortality of postmetamorphic frogs during a Bd epizootic, suggesting that the relationship between Bd and bacterial communities was not linked to variation in resistance to mortal disease and that Bd infection altered bacterial communities. In a controlled experiment, Bd infection significantly altered the R. sierrae microbiome, demonstrating a causal relationship. The response of microbial communities to Bd infection was remarkably consistent: Several bacterial taxa showed the same response to Bd infection across multiple field populations and the laboratory experiment, indicating a somewhat predictable interaction between Bd and the microbiome. The laboratory experiment demonstrates that Bd infection causes changes to amphibian skin bacterial communities, whereas the laboratory and field results together strongly support Bd disturbance as a driver of bacterial community change during natural disease dynamics.Symbiotic interactions between microbes and multicellular organisms are ubiquitous. In recent years, research to understand the complex microbial communities living in or on multicellular organisms (termed the microbiome) has sparked fundamental changes in our understanding of the biology of metazoans (1–5). The microbiome can affect host health directly by influencing metabolism (6), development (7), inflammation (8), or behavior (9), but it may also influence host health indirectly through interactions with infectious pathogens. The microbiome may interact with pathogens through competition for resources, release of antimicrobial compounds, contact-dependent antagonism, or modulation of the host immune response (10), and an “imbalanced” microbiome may leave the host more susceptible to pathogen infection (11, 12). At the same time, an invading pathogen may disrupt the microbiome (10, 13–15). Thus, the microbiome may play a role in disease resistance, or may itself be disturbed or altered by invading pathogens. Although a wealth of recent research has described associations between microbiome composition and a variety of syndromes in both humans and animals (16–25), documentation of microbiome responses to natural epidemics of known infectious pathogens is rare.Chytridiomycosis is an emerging infectious disease of amphibians caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd). Bd is an aquatic fungus that infects the skin of amphibians and disrupts osmoregulation, a critical function of amphibian skin (26). Chytridiomycosis can be fatal, and the severity of disease symptoms has been linked to Bd load, which is a measure of the density of Bd cells infecting the host (27, 28). Bd has a broad host range spanning hundreds of amphibian species, and has been implicated in population extinctions and species declines worldwide (29–34). Efforts to understand and mitigate the effects of Bd have led to research examining the potential for symbiotic bacteria to increase resistance to infection by the pathogen (35, 36). Bacterial species isolated from the skin of amphibians have been shown to inhibit the growth of Bd and other fungal pathogens in culture (37–39), possibly by producing antifungal metabolites (40, 41). In a controlled laboratory experiment, inundation of Rana muscosa with the bacterium Janthinobacterium lividium protected frogs from subsequent Bd infection (42). These and other studies highlight the possible role of bacteria in resistance to chytridiomycosis, but critical questions remain. First, most research has focused on the ability of cultured bacteria to prevent Bd infection, whereas very little is known about whether Bd infection alters the diverse skin microbiome. Examining this latter concept is critical both to a basic understanding of how the microbiome interacts with pathogens and to conservation efforts because Bd-induced perturbations of the microbiome could undermine attempts to mitigate effects of Bd infection through augmentation with particular bacteria. A second knowledge gap is the paucity of comprehensive culture-independent assessments of the amphibian microbiome, which are important because the vast majority of environmental and symbiotic microbes are not readily cultured, and culture-based methods can lead to severe underestimates of diversity and biased assessment of community composition (43). Few studies have applied next-generation sequencing methods to characterize the microbial communities on amphibian skin (44–47), and, to our knowledge, none have done so in the context of Bd infection. A final challenge to understanding interactions between Bd and bacteria stems from the difficulties of drawing direct connections between laboratory and field studies. Laboratory studies are essential for definitive identification of cause and effect. However, complex natural microbiomes can be impossible to recreate in the laboratory, and field studies are needed to show whether processes identified in the laboratory are relevant in nature.We present paired laboratory and field studies using high-throughput 16S amplicon pyrosequencing both to document associations between Bd infection and the amphibian skin bacterial microbiome in nature and to deduce causal relationships in an experiment. Our work centers on the Sierra Nevada yellow-legged frog, Rana sierrae, which is severely threatened by, and has already suffered drastic declines due to, Bd (28, 48). We surveyed frogs from four distinct R. sierrae populations to test if differences in skin bacterial communities are associated with the intensity of pathogen infection. We then conducted a laboratory experiment to establish causal relationships underlying Bd-bacterial community associations. The data establish a strong effect of Bd infection on the composition of the amphibian skin bacterial microbiome that is consistent between the laboratory experiment and naturally occurring Bd dynamics in wild frog populations.