NetCooperate: a network-based tool for inferring host-microbe and microbe-microbe cooperation

Background

Host-microbe and microbe-microbe interactions are often governed by the complex exchange of metabolites. Such interactions play a key role in determining the way pathogenic and commensal species impact their host and in the assembly of complex microbial communities. Recently, several studies have demonstrated how such interactions are reflected in the organization of the metabolic networks of the interacting species, and introduced various graph theory-based methods to predict host-microbe and microbe-microbe interactions directly from network topology. Using these methods, such studies have revealed evolutionary and ecological processes that shape species interactions and community assembly, highlighting the potential of this reverse-ecology research paradigm.

Results

NetCooperate is a web-based tool and a software package for determining host-microbe and microbe-microbe cooperative potential. It specifically calculates two previously developed and validated metrics for species interaction: the Biosynthetic Support Score which quantifies the ability of a host species to supply the nutritional requirements of a parasitic or a commensal species, and the Metabolic Complementarity Index which quantifies the complementarity of a pair of microbial organisms’ niches. NetCooperate takes as input a pair of metabolic networks, and returns the pairwise metrics as well as a list of potential syntrophic metabolic compounds.

Conclusions

The Biosynthetic Support Score and Metabolic Complementarity Index provide insight into host-microbe and microbe-microbe metabolic interactions. NetCooperate determines these interaction indices from metabolic network topology, and can be used for small- or large-scale analyses. NetCooperate is provided as both a web-based tool and an open-source Python module; both are freely available online at http://elbo.gs.washington.edu/software_netcooperate.html.     

A critical three-way junction is conserved in budding yeast and vertebrate telomerase RNAs

The telomerase ribonucleoprotein copies a short template within its integral RNA moiety onto eukaryotic chromosome ends, compensating for incomplete replication and degradation. Non-template regions of telomerase RNA (TER) are also crucial for telomerase function, yet they are highly divergent in sequence among species and their roles are largely unclear. Using both phylogenetic and mutational analyses, we predicted secondary structures for TERs from Kluyveromyces budding yeast species. A comparison of these secondary structure models with the published model for the Saccharomyces cerevisiae TER reveals a common arrangement into three long arms, a templating domain in the center and several conserved elements in the same positions within the structure. One of them, a three-way junction element, is highly conserved in budding yeast TERs. This element also shows sequence and structure similarity to the critical CR4-CR5 activating domain of vertebrate TERs. Mutational analysis in Kluyveromyces lactis confirmed that this element, and in particular the residues conserved across yeast and vertebrates, is critical for telomerase action both in vivo and in vitro. These findings demonstrate that despite the extreme divergence of TER sequences from different organisms, they do share conserved elements, which presumably carry out common roles in telomerase function.     

Gene loss rate: a probabilistic measure for the conservation of eukaryotic genes

The rate of conservation of a gene in evolution is believed to be correlated with its biological importance. Recent studies have devised various conservation measures for genes and have shown that they are correlated with several biological characteristics of functional importance. Specifically, the state-of-the-art propensity for gene loss (PGL) measure was shown to be strongly correlated with gene essentiality and its number of protein–protein interactions (PPIs). The observed correlation between conservation and functional importance varies however between conservation measures, underscoring the need for accurate and general measures for the rate of gene conservation. Here we develop a novel maximum-likelihood approach to computing the rate in which a gene is lost in evolution, motivated by the same principles as those underlying PGL. However, in difference to PGL which considers only the most parsimonious ancestral states of the internal nodes of the phylogenetic tree relating the species, our approach weighs in a probabilistic manner all possible ancestral states, and includes the branch length information as part of the probabilistic model. In application to data of 16 eukaryotic genomes, our approach shows higher correlations with experimental data than PGL, including data on gene lethality, level of connectivity in a PPI network and coherence within functionally related genes.     

An end to endless forms: epistasis, phenotype distribution bias, and nonuniform evolution

Studies of the evolution of development characterize the way in which gene regulatory dynamics during ontogeny constructs and channels phenotypic variation. These studies have identified a number of evolutionary regularities: (1) phenotypes occupy only a small subspace of possible phenotypes, (2) the influence of mutation is not uniform and is often canalized, and (3) a great deal of morphological variation evolved early in the history of multicellular life. An important implication of these studies is that diversity is largely the outcome of the evolution of gene regulation rather than the emergence of new, structural genes. Using a simple model that considers a generic property of developmental maps-the interaction between multiple genetic elements and the nonlinearity of gene interaction in shaping phenotypic traits-we are able to recover many of these empirical regularities. We show that visible phenotypes represent only a small fraction of possibilities. Epistasis ensures that phenotypes are highly clustered in morphospace and that the most frequent phenotypes are the most similar. We perform phylogenetic analyses on an evolving, developmental model and find that species become more alike through time, whereas higher-level grades have a tendency to diverge. Ancestral phenotypes, produced by early developmental programs with a low level of gene interaction, are found to span a significantly greater volume of the total phenotypic space than derived taxa. We suggest that early and late evolution have a different character that we classify into micro- and macroevolutionary configurations. These findings complement the view of development as a key component in the production of endless forms and highlight the crucial role of development in constraining biotic diversity and evolutionary trajectories.     

Evolution of learning in fluctuating environments: when selection favors both social and exploratory individual learning

Cumulative cultural change requires organisms that are capable of both exploratory individual learning and faithful social learning. In our model, an organism's phenotype is initially determined innately (by its genotypic value) or by social learning (copying a phenotype from the parental generation), and then may or may not be modified by individual learning (exploration around the initial phenotype). The environment alternates periodically between two states, each defined as a certain range of phenotypes that can survive. These states may overlap, in which case the same phenotype can survive in both states, or they may not. We find that a joint social and exploratory individual learning strategy-the strategy that supports cumulative culture-is likely to spread when the environmental states do not overlap. In particular, when the environmental states are contiguous and mutation is allowed among the genotypic values, this strategy will spread in either moderately or highly stable environments, depending on the exact nature of the individual learning applied. On the other hand, natural selection often favors a social learning strategy without exploration when the environmental states overlap. We find only partial support for the "consensus" view, which holds that individual learning, social learning, and innate determination of behavior will evolve at short, intermediate, and long environmental periodicities, respectively.     

Ex vivo rapamycin generates apoptosis-resistant donor Th2 cells that persist in vivo and prevent hemopoietic stem cell graft rejection

Because ex vivo rapamycin generates murine Th2 cells that prevent Graft-versus-host disease more potently than control Th2 cells, we hypothesized that rapamycin would generate Th2/Tc2 cells (Th2/Tc2.R cells) that abrogate fully MHC-disparate hemopoietic stem cell rejection more effectively than control Th2/Tc2 cells. In a B6-into-BALB/c graft rejection model, donor Th2/Tc2.R cells were indeed enriched in their capacity to prevent rejection; importantly, highly purified CD4+ Th2.R cells were also highly efficacious for preventing rejection. Rapamycin-generated Th2/Tc2 cells were less likely to die after adoptive transfer, accumulated in vivo at advanced proliferative cycles, and were present in 10-fold higher numbers than control Th2/Tc2 cells. Th2.R cells had a multifaceted, apoptosis-resistant phenotype, including: 1) reduced apoptosis after staurosporine addition, serum starvation, or CD3/CD28 costimulation; 2) reduced activation of caspases 3 and 9; and 3) increased anti-apoptotic Bcl-xL expression and reduced proapoptotic Bim and Bid expression. Using host-versus-graft reactivity as an immune correlate of graft rejection, we found that the in vivo efficacy of Th2/Tc2.R cells 1) did not require Th2/Tc2.R cell expression of IL-4, IL-10, perforin, or Fas ligand; 2) could not be reversed by IL-2, IL-7, or IL-15 posttransplant therapy; and 3) was intact after therapy with Th2.R cells relatively devoid of Foxp3 expression. We conclude that ex vivo rapamycin generates Th2 cells that are resistant to apoptosis, persist in vivo, and effectively prevent rejection by a mechanism that may be distinct from previously described graft-facilitating T cells.     

The DR1 and DR6 First Exons of Human Herpesvirus 6A Are Not Required for Virus Replication in Culture and Are Deleted in Virus Stocks That Replicate Well in T-Cell Lines

Human herpesvirus 6A (HHV-6A) and HHV-6B are lymphotropic viruses which replicate in cultured activated cord blood mononuclear cells (CBMCs) and in T-cell lines. Viral genomes are composed of 143-kb unique (U) sequences flanked by ∼8- to 10-kb left and right direct repeats, DR_L and DR_R. We have recently cloned HHV-6A (U1102) into bacterial artificial chromosome (BAC) vectors, employing DNA replicative intermediates. Surprisingly, HHV-6A BACs and their parental DNAs were found to contain short ∼2.7-kb DRs. To test whether DR shortening occurred during passaging in CBMCs or in the SupT1 T-cell line, we compared packaged DNAs from various passages. Restriction enzymes, PCR, and sequencing analyses have shown the following. (i) Early (1992) viral preparations from CBMCs contained ∼8-kb DRs. (ii) Viruses currently propagated in SupT1 cells contained ∼2.7-kb DRs. (iii) The deletion spans positions 60 to 5545 in DR_L, including genes encoded by DR1 through the first exon of DR6. The pac-2-pac-1 packaging signals, the DR7 open reading frame (ORF), and the DR6 second exon were not deleted. (iv) The DR_R sequence was similarly shortened by 5.4 kb. (v) The DR1 through DR6 first exon sequences were deleted from the entire HHV-6A BACs, revealing that they were not translocated into other genome locations. (vi) When virus initially cultured in CBMCs was passaged in SupT1 cells no DR shortening occurred. (vii) Viral stocks possessing short DRs replicated efficiently, revealing the plasticity of herpesvirus genomes. We conclude that the DR deletion occurred once, producing virus with advantageous growth “conquering” the population. The DR1 gene and the first DR6 exon are not required for propagation in culture.Human herpesvirus 6 (HHV-6) is a member of the Betaherpesvirus subfamily, as recently reviewed (46). The virus can enter hematopoietic cells, including T cells, B cells, natural killer (NK) cells, monocytes, and dendritic cells (DCs), as well as nonhematopoietic cells, as reviewed in references 8, 17, and 46. In culture, the virus replicates in activated peripheral blood lymphocytes (PBLs), cord blood mononuclear cells (CBMCs), and in T-cell lines (1, 17, 46). HHV-6 isolates fall into two distinct classes designated as HHV-6A and HHV-6B variants. The two variants can be distinguished by their restriction enzyme patterns, antigenicity, DNA sequences, and disease association (1, 36, 46). HHV-6B is the causative agent of roseola infantum, a prevalent children''s disease characterized by high fever and skin rash (47). In rare cases, the virus exhibits neurotropism and has been found in children experiencing convulsions up to lethal encephalitis (1, 21, 46, 48).HHV-6B reactivation from latency was found to occur in patients receiving immunosuppressive treatment in bone marrow and other transplantations. This was associated with febrile illness, delayed transplant engraftment, and neurological involvement, up to lethal encephalitis (5, 13, 34, 46). HHV-6A has thus far no clear disease association, although several studies have suggested central nervous system (CNS) tropism, including aggravation of symptoms in patients with multiple sclerosis (MS) (6, 14, 33, 41).HHV-6A and HHV-6B share general genomic architecture. The unit-length DNA molecules are approximately 160 kb, composed of a 143-kb unique (U) segment flanked by left and right direct repeats (DR_L and DR_R, respectively) (19, 24, 27, 46). The DRs are of sizes 8 to 10 kb in different viral isolates (2, 19, 24, 46). In both the HHV-6A and HHV-6B genomes, the herpesvirus conserved cleavage/packaging signals pac-1 and pac-2 (9, 15, 17) are located at the left and the right termini of the DRs (17, 19, 46). The PubMed sequence for the U1102 strain (accession no. {"type":"entrez-nucleotide","attrs":{"text":"NC_001664","term_id":"1344462938","term_text":"NC_001664"}}NC_001664) starts with the pac-1 signal at positions 1 to 56, followed by multiple copies of perfect and imperfect telomere-like sequences, up to position 418. It was suggested that the telomeric repeats may have originated from host cell chromosomal telomeres (43). Additionally, the DR encodes several open reading frames (ORFs), four of which are dealt with in our paper: (i) the spliced DR1 at positions 501 to 759 and 843 to 2653; (ii) DR5 at positions 3738 to 4164; (iii) the spliced DR6 at positions 4725 to 5028 and 5837 to 6720; and (iv) an ORF of DR7, at positions 5629 to 6720, partially overlapping the DR6 gene (20). Hollsberg and coworkers (37) have recently found that the homologous gene in HHV-6B encodes a nuclear protein that forms a complex with viral DNA processivity factor p41. Gompels and coworkers have also shown that DR1 and DR6 are partly homologous to the human cytomegalovirus (HCMV) US22 gene family. Both have a CXC motif: DR1 with homology to the HCMV US26 gene and DR6 with homology to the HCMV US22 gene (20). The map continues with reiterated perfect hexanucleotide telomeric sequences (GGGTAA)_n at positions 7655 to 8008 (19, 43). The number of telomeric repeats was found to vary in different viral strains (2, 43). The DR terminates with the pac-2 signal.We have recently cloned the intact HHV-6A genome into bacterial artificial chromosomes (BACs), by direct cloning of unit-length DNA produced from circular or head-to-tail replication intermediates into modified BAC vectors containing the green fluorescent protein (GFP) marker and ampicillin-puromycin (Amp-Puro) selection cassette (3). Surprisingly, the HHV-6A BAC clones as well as the parental HHV-6A (U1102) propagated in our laboratory in SupT1 cells were found to contain DRs of ∼2.7 kb instead of the expected ∼8- to 10-kb DRs, as in the early publications (19, 24, 27, 46) and in the PubMed sequence. This has raised the following questions. When did the deleted DRs arise? What was the detailed structure of deleted DRs?HHV-6 was discovered by Gallo and colleagues in 1986 (35), and viral isolates were obtained in multiple laboratories from AIDS patients, patients with lymphoproliferative disorders, and patients with roseola infantum (12, 26, 42, 45, 47). The isolates were propagated initially in activated PBLs and CBMCs and then in continuous T-cell lines, including HSB-2, J-JHAN, SupT1, Molt-3, and MT-4 (11, 45). The U1102 strain isolated by Downing and colleagues (12) was contributed to our laboratory by Robert Honess and was propagated first in activated PBLs and CBMCs (11, 18, 36, 45) and then in J-JHAN and SupT1 T cells (4, 30). To answer the question with regard to the origin of the short DRs and their structure, we have compared earlier viral HHV-6A passages with the currently propagated virus and the HHV-6A BAC clones. We describe here the detailed structure of the DR_L and DR_R of the “new” virus, containing the short ∼2.7-kb DR. We show that the deletion contained the left multiple repeats of telomere-like sequences and the ORFs from DR1 up to the DR6 first exon. Review of viral passaging since 1992 indicated that the deletion occurred spontaneously. The deleted viruses were stably and efficiently propagated in SupT1 T cells, indicating that the DR1 and DR6 first exons are not essential for virus in vitro replication.     

The effect of phenotypic plasticity on evolution in multipeaked fitness landscapes

When facing the challenge of developing an individual that best fits its environment, nature demonstrates an interesting combination of two fundamentally different adaptive mechanisms: genetic evolution and phenotypic plasticity. Following numerous computational models, it has become the accepted wisdom that lifetime acclimation (e.g. via learning) smooths the fitness landscape and consequently accelerates evolution. However, analytical studies, focusing on the effect of phenotypic plasticity on evolution in simple unimodal landscapes, have often found that learning hinders the evolutionary process rather than accelerating it. Here, we provide a general framework for studying the effect of plasticity on evolution in multipeaked landscapes and introduce a rigorous mathematical analysis of these dynamics. We show that the convergence rate of the evolutionary process in a given arbitrary one-dimensional fitness landscape is dominated by the largest descent (drawdown) in the landscape and provide numerical evidence to support an analogous dominance also in multidimensional landscapes. We consider several schemes of phenotypic plasticity and examine their effect on the landscape drawdown, identifying the conditions under which phenotypic plasticity is advantageous. The lack of such a drawdown in unimodal landscapes vs. its dominance in multipeaked landscapes accounts for the seemingly contradictory findings of previous studies.     

Mom1 Is a Semi-Dominant Modifier of Intestinal Adenoma Size and Multiplicity in Min/+ Mice

The intestinal tumor multiplicity in mice heterozygous for Apc(Min) is strongly modulated by genetic background. On the sensitive C57BL/6J (B6) background, mice develop large numbers of intestinal adenomas. The AKR/J (AKR) strain carries alleles that correlate with a strong reduction in tumor multiplicity. To study the effect of one of these modifiers, Mom1, we have generated a mouse line in which the AKR allele of Mom1 is carried on the sensitive B6 genetic background. This strain was produced by using a marker-assisted selection method to eliminate unlinked AKR alleles more rapidly. The application and efficiency of this method are discussed. We used this strain to determine that Mom1 affects both tumor multiplicity and tumor size in a semi-dominant fashion.