The impact of non-lethal synergists on the population and evolutionary dynamics of host-pathogen interactions

Multiple pathogenic infections can influence disease transmission and virulence, and have important consequences for understanding the community ecology and epidemiology of host-pathogen interactions. Here the population and evolutionary dynamics of a host-pathogen interaction with free-living stages are explored in the presence of a non-lethal synergist that hosts must tolerate. Through the coupled effects on pathogen transmission, host mass gain and allometry it is shown how investing in tolerance to a non-lethal synergist can lead to a broad range of different population dynamics. The effects of the synergist on pathogen fitness are explored through a series of life-history trait trade-offs. Coupling trade-offs between pathogen yield and pathogen speed of kill and the presence of a synergist favour parasites that have faster speeds of kill. This evolutionary change in pathogen characteristics is predicted to lead to stable population dynamics. Evolutionary analysis of tolerance of the synergist (strength of synergy) and lethal pathogen yield show that decreasing tolerance allows alternative pathogen strategies to invade and replace extant strategies. This evolutionary change is likely to destabilise the host-pathogen interaction leading to population cycles. Correlated trait effects between speed of kill and tolerance (strength of synergy) show how these traits can interact to affect the potential for the coexistence of multiple pathogen strategies. Understanding the consequences of these evolutionary relationships is important for the both the evolutionary and population dynamics of host-pathogen interactions.     

PARTIAL RESISTANCE IN THE LINUM‐MELAMPSORA HOST–PATHOGEN SYSTEM: DOES PARTIAL RESISTANCE MAKE THE RED QUEEN RUN SLOWER?

Five levels of disease expression were scored in a cross-inoculation study of 120 host and 60 pathogen lines of wild flax Linum marginale and its rust fungus Melampsora lini sampled from six natural populations. Patterns of partial resistance showed clear evidence of gene-for-gene interactions, with particular levels of partial resistance occurring in specific host-pathogen combinations. Sympatric and putatively more highly coevolved host-pathogen combinations had a lower frequency of partial resistance types relative to allopatric combinations. Sympatric host-pathogen combinations also showed a lower diversity of resistance responses, but there was a trend toward a greater fraction of this variance being determined by pathogen-genotype × host-genotype interactions. In this system, there was no evidence that partial resistances slow host-pathogen coevolution. The analyses show that if variation is generated by among population host or pathogen dispersal, then coevolution occurs largely by pathogens overcoming the partial resistances that are generated.     

Pathogen interactions, population cycles, and phase shifts

Interspecific pathogen interactions can profoundly affect pathogen population dynamics and the efficacy of control strategies. However, many pathogens exhibit cyclic abundance patterns (e.g., seasonality), and temporal asynchrony between interacting pathogens could reduce the impact of those interactions. Here we use an extension of our previously published model to investigate the effects of cycles on pathogen interaction. We demonstrate that host immune memory can maintain the impact of an interaction, even when the effector pathogen abundance is low or the pathogen is absent. Paradoxically, immune memory can result in pathogens interacting more strongly when temporally out of phase. We find that interactions between species can result in changes to the temporal pattern of the affected species. We further demonstrate that this may be observed in a natural host-pathogen system. Given the continuing debate regarding the relevance of pathogen interactions in natural systems and increasing concern about treatment strategies for coinfections, both the discovery of a shift in cycle in empirical data and the mechanism by which we identified it are important. Finally, because the model structure used here is analogous to models of a simple predator-prey system, we also consider the consequences of these findings in the context of that system.     

Host‐plant genotypic diversity and community genetic interactions mediate aphid spatial distribution

Genetic variation in plants can influence the community structure of associated species, through both direct and indirect interactions. Herbivorous insects are known to feed on a restricted range of plants, and herbivore preference and performance can vary among host plants within a species due to genetically based traits of the plant (e.g., defensive compounds). In a natural system, we expect to find genetic variation within both plant and herbivore communities and we expect this variation to influence species interactions. Using a three‐species plant‐aphid model system, we investigated the effect of genetic diversity on genetic interactions among the community members. Our system involved a host plant (Hordeum vulgare) that was shared by an aphid (Sitobion avenae) and a hemi‐parasitic plant (Rhinanthus minor). We showed that aphids cluster more tightly in a genetically diverse host‐plant community than in a genetic monoculture, with host‐plant genetic diversity explaining up to 24% of the variation in aphid distribution. This is driven by differing preferences of the aphids to the different plant genotypes and their resulting performance on these plants. Within the two host‐plant diversity levels, aphid spatial distribution was influenced by an interaction among the aphid's own genotype, the genotype of a competing aphid, the origin of the parasitic plant population, and the host‐plant genotype. Thus, the overall outcome involves both direct (i.e., host plant to aphid) and indirect (i.e., parasitic plant to aphid) interactions across all these species. These results show that a complex genetic environment influences the distribution of herbivores among host plants. Thus, in genetically diverse systems, interspecific genetic interactions between the host plant and herbivore can influence the population dynamics of the system and could also structure local communities. We suggest that direct and indirect genotypic interactions among species can influence community structure and processes.     

Evolution of Pathogen Specialisation in a Host Metapopulation: Joint Effects of Host and Pathogen Dispersal

Metapopulation processes are important determinants of epidemiological and evolutionary dynamics in host-pathogen systems, and are therefore central to explaining observed patterns of disease or genetic diversity. In particular, the spatial scale of interactions between pathogens and their hosts is of primary importance because migration rates of one species can affect both spatial and temporal heterogeneity of selection on the other. In this study we developed a stochastic and discrete time simulation model to specifically examine the joint effects of host and pathogen dispersal on the evolution of pathogen specialisation in a spatially explicit metapopulation. We consider a plant-pathogen system in which the host metapopulation is composed of two plant genotypes. The pathogen is dispersed by air-borne spores on the host metapopulation. The pathogen population is characterised by a single life-history trait under selection, the infection efficacy. We found that restricted host dispersal can lead to high amount of pathogen diversity and that the extent of pathogen specialisation varied according to the spatial scale of host-pathogen dispersal. We also discuss the role of population asynchrony in determining pathogen evolutionary outcomes.     

Host pathogen protein interactions predicted by comparative modeling

Pathogens have evolved numerous strategies to infect their hosts, while hosts have evolved immune responses and other defenses to these foreign challenges. The vast majority of host-pathogen interactions involve protein-protein recognition, yet our current understanding of these interactions is limited. Here, we present and apply a computational whole-genome protocol that generates testable predictions of host-pathogen protein interactions. The protocol first scans the host and pathogen genomes for proteins with similarity to known protein complexes, then assesses these putative interactions, using structure if available, and, finally, filters the remaining interactions using biological context, such as the stage-specific expression of pathogen proteins and tissue expression of host proteins. The technique was applied to 10 pathogens, including species of Mycobacterium, apicomplexa, and kinetoplastida, responsible for "neglected" human diseases. The method was assessed by (1) comparison to a set of known host-pathogen interactions, (2) comparison to gene expression and essentiality data describing host and pathogen genes involved in infection, and (3) analysis of the functional properties of the human proteins predicted to interact with pathogen proteins, demonstrating an enrichment for functionally relevant host-pathogen interactions. We present several specific predictions that warrant experimental follow-up, including interactions from previously characterized mechanisms, such as cytoadhesion and protease inhibition, as well as suspected interactions in hypothesized networks, such as apoptotic pathways. Our computational method provides a means to mine whole-genome data and is complementary to experimental efforts in elucidating networks of host-pathogen protein interactions.     

VARIATION IN PATHOGENICITY AMONG AND WITHIN POPULATIONS OF THE FUNGUS PHOMOPSIS SUBORDINARIA INFECTING PLANTAGO LANCEOLATA

Pathogenicity of isolates of the fungus Phomopsis subordinaria, sampled in three scarcely or heavily infected populations of Plantago lanceolata, was investigated on three different host genotypes. The expression of the pathogen appeared to be quantitative rather than qualitative in character, which suggests polygenic inheritance of host susceptibility. Significant statistical interaction between pathogen and host pointed to some degree of physiological specialization between them. None of the individual host-pathogen combinations was found to contribute significantly to the interaction. Differences in mean pathogenicity between the pathogen populations could not explain the different intensities of disease observed in the field. As the variation in susceptibility between populations of the host at the same three locations also cannot account for the differences in intensity of disease in the field, it can be concluded that environmental factors (in particular weevils that spread the disease) are important for the development of the disease. In one of the populations, the spatial scale at which variation within the pathogen occurs was determined. It appeared that the pathogen varied in pathogenicity in the field, even among scapes within an individual host plant. The consequences of this scale of variation in the pathogen are discussed for the dynamics and evolution of the pathosystem.     

Transmission rates and adaptive evolution of pathogens in sympatric heterogeneous plant populations

Diversification in agricultural cropping patterns is widely practised to delay the build-up of virulent races that can overcome host resistance in pathogen populations. This can lead to balanced polymorphism, but the long-term consequences of this strategy for the evolution of crop pathogen populations are still unclear. The widespread occurrence of sibling species and reproductively isolated sub-species among fungal and oomycete plant pathogens suggests that evolutionary divergence is common. This paper develops a mathematical model of host-pathogen interactions using a simple framework of two hosts to analyse the influences of sympatric host heterogeneity on the long-term evolutionary behaviour of plant pathogens. Using adaptive dynamics, which assumes that sequential mutations induce small changes in pathogen fitness, we show that evolutionary outcomes strongly depend on the shape of the trade-off curve between pathogen transmission on sympatric hosts. In particular, we determine the conditions under which the evolutionary branching of a monomorphic into a dimorphic population occurs, as well as the conditions that lead to the evolution of specialist (single host range) or generalist (multiple host range) pathogen populations.     

Demographic and genetic invasion history of a 9-year-old roadside population of Bunias orientalis L. (Brassicaceae)

The population history of a 9-year-old roadside population of the invasive plant Bunias orientalis was reconstructed by demographic analysis including size, position, age (determined by herbchronology) and RAPD-PCR patterns of individual plants. We evaluated emerging patterns of population growth and genetic structure during a full period of population development under typical site conditions (anthropogenic disturbance) and their possible consequences for the invasion potential of the species. The population has grown rapidly and continuously (though with slowing geometric population increase) during the 9 years since its foundation, filling the space available in the study area. Genetic variation (RAPD markers) was already high in the founder cohorts and remained at the same level throughout population development (variance fluctuations <15%). Both results may be related to the mowing management at the site which seems to promote population growth of B. orientalis relative to other co-occuring species and to prevent the genetic drift and the development of spatial genetic structure that would be expected under isolation-by-distance models. Large founder plants had comparatively low genetic variance and were more closely related to younger cohorts than were small founder plants, indicating that selection acted during population development. Overall, the current anthropogenic disturbance regimes may contribute to high genetic variability by artificially increasing gene flow and thereby promoting the adaptability of invasive species to the often unpredictable conditions at disturbed sites. Our approach using retrospective demographic investigation allows the detection of spatio-temporal microscale patterns in genetic and phenotypic variation. Thus it allows a thorough understanding of local invasions of perennial herbaceous plants. Received: 23 November 1998 / Accepted: 14 April 1999     

LOCAL MALADAPTATION IN THE ANTHER-SMUT FUNGUS MICROBOTRYUM VIOLACEUM TO ITS HOST PLANT SILENE LATIFOLIA: EVIDENCE FROM A CROSS-INOCULATION EXPERIMENT

Conventional wisdom holds that parasites evolve more rapidly than their hosts and are therefore locally adapted, that is, better at exploiting sympatric than allopatric hosts. We studied local adaptation in the insect-transmitted fungal pathogen Microbotryum violaceum and its host plant Silene latifolia. Infection success was tested in sympatric (local) and allopatric (foreign) combinations of pathogen and host from 14 natural populations from a metapopulation. Seedlings from up to 10 seed families from each population were exposed to sporidial suspensions from each of four fungal strains derived from the same population, from a near-by population (< 10 km distance), and from two populations at an intermediate (< 30 km) and remote (< 170 km) distance, respectively. We obtained significant pathogen X plant interactions in infection success (proportion of diseased plants) at both fungal population and strain level. There was an overall pattern of local maladaptation of this pathogen: average fungal infection success was significantly lower on sympatric hosts (mean proportion of diseased plants = 0.32 ± 0.03 SE) than on allopatric hosts (0.40 ± 0.02). Five of the 14 fungal populations showed no strong reduction in infection success on sympatric hosts, and three even tended to perform better on sympatric hosts. This pattern is consistent with models of time-lagged cycles predicting patterns of local adaptation in host-parasite systems to emerge only on average. Several factors may restrict the evolutionary potential of this pathogen relative to that of its host. First, a predominantly selfing breeding system may limit its ability to generate new virulence types by sexual recombination, whereas the obligately outcrossing host 5. latifolia may profit from rearrangement of resistance alleles by random mating. Second, populations often harbor only a few infected individuals, so virulence variation may be further reduced by drift. Third, migration rates among host plant populations are much higher than among pathogen populations, possibly because pollinators prefer healthy over diseased plants. Migration among partly isolated populations may therefore introduce novel host plant resistance variants more often than novel parasite virulence variants. That migration contributes to the coevolutionary dynamics in this system is supported by the geographic pattern of infectivity. Infection success increased over the first 10–km range of host-pathogen population distances, which is likely the natural range of gene exchange.     

Spatial and temporal variation in host–parasitoid interactions: leafcutter ant hosts and their phorid parasitoids

1. Parasitoid–host interactions are important components of ecological communities. Although parasitoid–host interactions are strongly shaped by evolutionary history, the abundance of both the parasitoid and the host may have a role in determining the nature of the interaction once phylogenetic relationships are considered. 2. Leafcutter ants are hosts of phorid parasitoids and represent a well‐defined and specialised module within a larger network of ant–symbiont interactions. A low specificity host taxa and a positive association between host abundance and parasitoid interaction frequency were expected due to the close phylogenetic relatedness of the hosts. 3. The interactions among all species of leafcutter ants and their parasitoids were quantified in two localities with different species richness. This study also characterised the spatial‐temporal variability of these interactions, determined the patterns of parasitoid specificity and host selection, and tested for an association between host abundance and parasitoid interaction frequency. 4. Contrary to expectation, most parasitoid species were highly specialised and interaction frequency for parasitoid species was not related to host abundance. All host ant species were attacked by more than one phorid species. Some phorid species used more than one host species and showed preference for the same species over space and time, suggesting that there are physiological and/or behavioural restrictions on host use. 5. These results show that there is a tendency for specialisation even when hosts are highly similar in their ecology. From a biological control perspective, these parasitoids may be effective candidates, due to the high specificity of some species and little host‐use variation through time.     

Population dynamics and sex ratio of a parasitoid altered by fungal-infected diet of host butterfly

Variation of host quality affects population dynamics of parasitoids, even at the landscape scale. What causes host quality to vary and the subsequent mechanisms by which parasitoid population dynamics are affected can be complex. Here, we examine the indirect interaction of a plant pathogen with a parasitoid wasp. Under laboratory conditions, parasitoids from hosts fed fungus-infected plants weighed less than those from hosts fed uninfected plants, indicating that the fungus causes the hosts to be of poor quality. However, parasitoids reared from hosts fed fungal-infected diet also tended to be female, a characteristic associated with high host quality. The pathogen, herbivore and parasitoid persist regionally as metapopulations in a shared landscape in Aland, Finland. In an analysis of the metapopulation dynamics of the parasitoid over 6 years, the probability of colonization of a host population increased by more than twofold in patches occupied by the plant pathogen. While we cannot determine that the relationship is causal, a compelling explanation is that the plant pathogen facilitates the establishment by the parasitoid by increasing the fraction of female offspring. This is a novel mechanism of spatial multi-trophic level interactions.     

Differential trends in the codon usage patterns in HIV-1 genes

Host-pathogen interactions underlie one of the most complex evolutionary phenomena resulting in continual adaptive genetic changes, where pathogens exploit the host's molecular resources for growth and survival, while hosts try to eliminate the pathogen. Deciphering the molecular basis of host-pathogen interactions is useful in understanding the factors governing pathogen evolution and disease propagation. In host-pathogen context, a balance between mutation, selection, and genetic drift is known to maintain codon bias in both organisms. Studies revealing determinants of the bias and its dynamics are central to the understanding of host-pathogen evolution. We considered the Human Immunodeficiency Virus (HIV) type 1 and its human host to search for evolutionary signatures in the viral genome. Positive selection is known to dominate intra-host evolution of HIV-1, whereas high genetic variability underlies the belief that neutral processes drive inter-host differences. In this study, we analyze the codon usage patterns of HIV-1 genomes across all subtypes and clades sequenced over a period of 23 years. We show presence of unique temporal correlations in the codon bias of three HIV-1 genes illustrating differential adaptation of the HIV-1 genes towards the host preferred codons. Our results point towards gene-specific translational selection to be an important force driving the evolution of HIV-1 at the population level.     

The acquisition of hypovirulence in host-pathogen systems with three trophic levels

A major focus of research on the dynamics of host-pathogen interactions has been the evolution of pathogen virulence, which is defined as the loss in host fitness due to infection. It is usually assumed that changes in pathogen virulence are the result of selection to increase pathogen fitness. However, in some cases, pathogens have acquired hypovirulence by themselves becoming infected with hyperparasites. For example, the chestnut blight fungus Cryphonectria parasitica has become hypovirulent in some areas by acquiring a double-stranded RNA hyperparasite that debilitates the pathogen, thereby reducing its virulence to the host. In this article, we develop and analyze a mathematical model of the dynamics of host-pathogen interactions with three trophic levels. The system may be dominated by either uninfected (virulent) or hyperparasitized (hypovirulent) pathogens, or by a mixture of the two. Hypovirulence may allow some recovery of the host population, but it can also harm the host population if the hyperparasite moves the transmission rate of the pathogen closer to its evolutionarily stable strategy. In the latter case, the hyperparasite is effectively a mutualist of the pathogen. Selection among hyperparasites will often minimize the deleterious effects, or maximize the beneficial effects, of the hyperparasite on the pathogen. Increasing the frequency of multiple infections of the same host individual promotes the acquisition of hypovirulence by increasing the opportunity for horizontal transmission of the hyperparasite. This effect opposes the usual theoretical expectation that multiple infections promote the evolution of more virulent pathogens via selection for rapid growth within hosts.     

Density‐related effects on the infectivity and aggressiveness of a sterilising smut in a wild population of Digitaria sanguinalis

Understanding host–pathogen evolutionary dynamics needs characterisation and quantification of processes occurring at many spatiotemporal scales. With this aim, the effects of smut on a naturally infected population of the summer annual Digitaria sanguinalis were followed for 4 years in an uncropped field. The main purpose of the study was to quantify the effects of within‐population density on the infectivity and the aggressiveness of the pathogen in a range of densities that occurred naturally. The infectivity‐related variable measured was the proportion of smutted plants at the end of each growing season; proportions were analysed using a generalised linear model with a binomial distribution considering the year, the density and their interaction as effects. The aggressiveness‐related variables chosen were the number of smutted inflorescences per plant and per area, obtained over the last 2 years; they were analysed by means of ancova considering disease status (seeded or smutted), year, density and all the interactions between them. Although the disease is monocyclic, results showed clearly that infectivity increased with plant density. The number of inflorescences per plant was 1.5 times higher in smutted plants than in healthy plants throughout the range of densities. This variable declined when density increased, but as the infectivity increased at a higher rate, the aggressiveness also increased with density. The surprising results on infectivity are discussed in the context of current knowledge of plant–pathogen interaction dynamics, as well as neighbour effects on pathogen aggressiveness. Moreover, the results could be useful to develop weed biological control strategies.     

Cellular versus viral microRNAs in host–virus interaction

MicroRNAs (miRNAs) mark a new paradigm of RNA-directed gene expression regulation in a wide spectrum of biological systems. These small non-coding RNAs can contribute to the repertoire of host-pathogen interactions during viral infection. This interplay has important consequences, both for the virus and the host. There have been reported evidences of host-cellular miRNAs modulating the expression of various viral genes, thereby playing a pivotal role in the host–pathogen interaction network. In the hide-and-seek game between the pathogens and the infected host, viruses have evolved highly sophisticated gene-silencing mechanisms to evade host-immune response. Recent reports indicate that virus too encode miRNAs that protect them against cellular antiviral response. Furthermore, they may exploit the cellular miRNA pathway to their own advantage. Nevertheless, our increasing knowledge of the host–virus interaction at the molecular level should lead us toward possible explanations to viral tropism, latency and oncogenesis along with the development of an effective, durable and nontoxic antiviral therapy. Here, we summarize the recent updates on miRNA-induced gene-silencing mechanism, modulating host–virus interactions with a glimpse of the miRNA-based antiviral therapy for near future.     

Evolutionary diversification through hybridization in a wild host-pathogen interaction

Coevolutionary outcomes between interacting species are predicted to vary across landscapes, as environmental conditions, gene flow, and the strength of selection vary among populations. Using a combination of molecular, experimental, and field approaches, we describe how broad-scale patterns of environmental heterogeneity, genetic divergence, and regional adaptation have the potential to influence coevolutionary processes in the Linum marginale-Melampsora lini plant-pathogen interaction. We show that two genetically and geographically divergent pathogen lineages dominate interactions with the host across Australia, and demonstrate a hybrid origin for one of the lineages. We further demonstrate that the geographic divergence of the two lineages of M. lini in Australia is related to variation among lineages in virulence, life-history characteristics, and response to environmental conditions. When correlated with data describing regional patterns of variation in host resistance diversity and mating system these observations highlight the potential for gene flow and geographic selection mosaics to generate and maintain coevolutionary diversification in long-standing host-pathogen interactions.     

Variation in host resistance and pathogen selective value in the interaction between Pinus sylvestris and the fungus Crumenulopsis sororia

There have been many studies of plant pathogen evolution in systems showing gene-for-gene control of host resistance. However little is known about situations, exemplified by Scots pine, Pinus sylvestris, and its fungal pathogen Crumenulopsis sororia, where variation in host resistance is quantitative. In a field experiment genetically marked isolates of C. sororia from three natural populations were reciprocally inoculated on 1- and 2-year-old branch tissue of P. sylvestris in the three sites from which they had been collected. Quantitative variation in host resistance was measured by comparing the performance of the same inocula on different host populations, individuals and tissues. The selective value of isolates derived from different populations was estimated by comparing the frequency of genotypes in lesion re-isolations with those in the initial inoculum mixtures. Host resistance varied significantly among populations, individuals within populations and between 1- and 2-year-old branch tissue of P. sylvestris. Large differences in the relative selective values of C. sororia isolates from different populations were detected. The selective value of pathogens was independent of the host population on which they were inoculated. However, their selective value did depend on the age of the tissue on which they grew. The implications of these results for modelling evolution in pathogen-host interactions that lack gene-for-gene determination of host resistance are discussed.     

Host-pathogen interactions: a proteomic view

Host-pathogen interactions reflect the balance of host defenses and pathogen virulence mechanisms. Advances in proteomic technologies now afford opportunities to compare protein content between complex biologic systems ranging from cells to animals and clinical samples. Thus, it is now possible to characterize host-pathogen interactions from a global proteomic view. Most reports to date focus on cataloging protein content of pathogens and identifying virulence-associated proteins or proteomic alterations in host response. A more in-depth understanding of host-pathogen interactions has the potential to improve our mechanistic understanding of pathogenicity and virulence, thereby defining novel therapeutic and vaccine targets. In addition, proteomic characterization of the host response can provide pathogen-specific host biomarkers for rapid pathogen detection and characterization, as well as for early and specific detection of infectious diseases. A review of host-pathogen interactions focusing on proteomic analyses of both pathogen and host will be presented. Relevant genomic studies and host model systems will be also be discussed.     

Genetic identity and genotype × genotype interactions between symbionts outweigh species level effects in an insect microbiome

Microbial symbionts often alter the phenotype of their host. Benefits and costs to hosts depend on many factors, including host genotype, symbiont species and genotype, and environmental conditions. Here, we present a study demonstrating genotype-by-genotype (G×G) interactions between multiple species of endosymbionts harboured by an insect, and the first to quantify the relative importance of G×G interactions compared with species interactions in such systems. In the most extensive study to date, we microinjected all possible combinations of five Hamiltonella defensa and five Fukatsuia symbiotica (X-type; PAXS) isolates into the pea aphid, Acyrthosiphon pisum. We applied several ecological challenges: a parasitoid wasp, a fungal pathogen, heat shock, and performance on different host plants. Surprisingly, genetic identity and genotype × genotype interactions explained far more of the phenotypic variation (on average 22% and 31% respectively) than species identity or species interactions (on average 12% and 0.4%, respectively). We determined the costs and benefits associated with co-infection, and how these compared to corresponding single infections. All phenotypes were highly reliant on individual isolates or interactions between isolates of the co-infecting partners. Our findings highlight the importance of exploring the eco-evolutionary consequences of these highly specific interactions in communities of co-inherited species.Subject terms: Microbial ecology, Community ecology, Symbiosis, Microbiome, Evolution