Key questions in the genetics and genomics of eco-evolutionary dynamics

Increasing acceptance that evolution can be ‘rapid'' (or ‘contemporary'') has generated growing interest in the consequences for ecology. The genetics and genomics of these ‘eco-evolutionary dynamics'' will be—to a large extent—the genetics and genomics of organismal phenotypes. In the hope of stimulating research in this area, I review empirical data from natural populations and draw the following conclusions. (1) Considerable additive genetic variance is present for most traits in most populations. (2) Trait correlations do not consistently oppose selection. (3) Adaptive differences between populations often involve dominance and epistasis. (4) Most adaptation is the result of genes of small-to-modest effect, although (5) some genes certainly have larger effects than the others. (6) Adaptation by independent lineages to similar environments is mostly driven by different alleles/genes. (7) Adaptation to new environments is mostly driven by standing genetic variation, although new mutations can be important in some instances. (8) Adaptation is driven by both structural and regulatory genetic variation, with recent studies emphasizing the latter. (9) The ecological effects of organisms, considered as extended phenotypes, are often heritable. Overall, the study of eco-evolutionary dynamics will benefit from perspectives and approaches that emphasize standing genetic variation in many genes of small-to-modest effect acting across multiple traits and that analyze overall adaptation or ‘fitness''. In addition, increasing attention should be paid to dominance, epistasis and regulatory variation.     

Phage and bacteria support mutual diversity in a narrowing staircase of coexistence

The competitive exclusion principle states that phage diversity M should not exceed bacterial diversity N. By analyzing the steady-state solutions of multistrain equations, we find a new constraint: the diversity N of bacteria living on the same resources is constrained to be M or M+1 in terms of the diversity of their phage predators. We quantify how the parameter space of coexistence exponentially decreases with diversity. For diversity to grow, an open or evolving ecosystem needs to climb a narrowing ‘diversity staircase'' by alternatingly adding new bacteria and phages. The unfolding coevolutionary arms race will typically favor high growth rate, but a phage that infects two bacterial strains differently can occasionally eliminate the fastest growing bacteria. This context-dependent fitness allows abrupt resetting of the ‘Red-Queen''s race'' and constrains the local diversity.     

Conjugative plasmids: vessels of the communal gene pool

Comparative whole-genome analyses have demonstrated that horizontal gene transfer (HGT) provides a significant contribution to prokaryotic genome innovation. The evolution of specific prokaryotes is therefore tightly linked to the environment in which they live and the communal pool of genes available within that environment. Here we use the term supergenome to describe the set of all genes that a prokaryotic ‘individual’ can draw on within a particular environmental setting. Conjugative plasmids can be considered particularly successful entities within the communal pool, which have enabled HGT over large taxonomic distances. These plasmids are collections of discrete regions of genes that function as ‘backbone modules’ to undertake different aspects of overall plasmid maintenance and propagation. Conjugative plasmids often carry suites of ‘accessory elements’ that contribute adaptive traits to the hosts and, potentially, other resident prokaryotes within specific environmental niches. Insight into the evolution of plasmid modules therefore contributes to our knowledge of gene dissemination and evolution within prokaryotic communities. This communal pool provides the prokaryotes with an important mechanistic framework for obtaining adaptability and functional diversity that alleviates the need for large genomes of specialized ‘private genes’.     

Short-term observations of marine bacterial and viral communities: patterns,connections and resilience

Observation of short-term temporal variation in bacterial and viral communities is important for understanding patterns of aquatic microbial diversity. We collected surface seawater once daily for 38 consecutive days with seven more samples interspersed over 40 more days at one location ∼2 km from Santa Catalina Island, California. Bacterial communities were analyzed by automated ribosomal intergenic spacer analysis (ARISA) and viral communities were analyzed by terminal restriction fragment length polymorphism (TRFLP) of the conserved T4-like myoviral gene encoding the major capsid protein (g23). Common bacterial and viral taxa were consistently dominant, and relatively few displayed dramatic increases/decreases or ‘boom/bust'' patterns that might be expected from dynamic predator-prey interactions. Association network analysis showed most significant covariations (associations) occurred among bacterial taxa or among viral taxa and there were several modular (highly-interconnected) associations (P⩽0.005). Associations observed between bacteria and viruses (P⩽0.005) occurred with a median time lag of 2 days. Regression of all pairwise Bray-Curtis similarities between samples indicated a rate of bacterial community change that slows from 2.1%–0.18% per day over a week to 2 months; the rate stays around 0.4% per day for viruses. Our interpretation is that, over the scale of days, individual bacterial and viral OTUs can be dynamic and patterned; resulting in statistical associations regarded as potential ecological interactions. However, over the scale of weeks, average bacterial community variation is slower, suggesting that there is strong community-level ecological resilience, that is, a tendency to converge towards a ‘mean'' microbial community set by longer-term controlling factors.     

Huanglongbing alters the structure and functional diversity of microbial communities associated with citrus rhizosphere

The diversity and stability of bacterial communities present in the rhizosphere heavily influence soil and plant quality and ecosystem sustainability. The goal of this study is to understand how ‘Candidatus Liberibacter asiaticus'' (known to cause Huanglongbing, HLB) influences the structure and functional potential of microbial communities associated with the citrus rhizosphere. Clone library sequencing and taxon/group-specific quantitative real-time PCR results showed that ‘Ca. L. asiaticus'' infection restructured the native microbial community associated with citrus rhizosphere. Within the bacterial community, phylum Proteobacteria with various genera typically known as successful rhizosphere colonizers were significantly greater in clone libraries from healthy samples, whereas phylum Acidobacteria, Actinobacteria and Firmicutes, typically more dominant in the bulk soil were higher in ‘Ca. L. asiaticus''-infected samples. A comprehensive functional microarray GeoChip 3.0 was used to determine the effects of ‘Ca. L. asiaticus'' infection on the functional diversity of rhizosphere microbial communities. GeoChip analysis showed that HLB disease has significant effects on various functional guilds of bacteria. Many genes involved in key ecological processes such as nitrogen cycling, carbon fixation, phosphorus utilization, metal homeostasis and resistance were significantly greater in healthy than in the ‘Ca. L. asiaticus''-infected citrus rhizosphere. Our results showed that the microbial community of the ‘Ca. L. asiaticus''-infected citrus rhizosphere has shifted away from using more easily degraded sources of carbon to the more recalcitrant forms. Overall, our study provides evidence that the change in plant physiology mediated by ‘Ca. L. asiaticus'' infection could elicit shifts in the composition and functional potential of rhizosphere microbial communities. In the long term, these fluctuations might have important implications for the productivity and sustainability of citrus-producing agro-ecosystems.     

Evolution of symbiotic organs and endosymbionts in lygaeid stinkbugs

We investigated seed bugs of the genus Nysius (Insecta: Hemiptera: Lygaeidae) for their symbiotic bacteria. From all the samples representing 4 species, 18 populations and 281 individuals, specific bacterial 16S rRNA gene sequences were consistently identified, which formed a distinct clade in the Gammaproteobacteria. In situ hybridization showed that the bacterium was endocellularly localized in a pair of large bacteriomes that were amorphous in shape, deep red in color, and in association with gonads. In the ovary of adult females, the endosymbiont was also localized in the ‘infection zone'' in the middle of each germarium and in the ‘symbiont ball'' at the anterior pole of each oocyte, indicating vertical transmission of the endosymbiont through the ovarial passage. Phylogenetic analyses based on bacterial 16S rRNA, groEL and gyrB genes consistently supported a coherent monophyly of the Nysius endosymbionts. The possibility of a sister relationship to ‘Candidatus Kleidoceria schneideri'', the bacteriome-associated endosymbiont of a lygaeid bug Kleidocerys resedae, was statistically rejected, indicating independent evolutionary origins of the endosymbionts in the Lygaeidae. The endosymbiont genes consistently exhibited AT-biased nucleotide compositions and accelerated rates of molecular evolution, and the endosymbiont genome was only 0.6 Mb in size. The endosymbiont phylogeny was congruent with the host insect phylogeny, suggesting strict vertical transmission and host–symbiont co-speciation over evolutionary time. Based on these results, we discuss the evolution of bacteriomes and endosymbionts in the Heteroptera, most members of which are associated with gut symbiotic bacteria. The designation ‘Candidatus Schneideria nysicola'' is proposed for the endosymbiont clade.     

Variation in local carrying capacity and the individual fate of bacterial colonizers in the phyllosphere

Using a phyllosphere model system, we demonstrated that the term ‘carrying capacity'', as it is commonly used in microbial ecology, needs to be understood as the sum of many ‘local carrying capacities'' in order to better explain and predict the course and outcome of bacterial colonization of an environment. Using a green fluorescent protein-based bioreporter system for the quantification of reproductive success (RS) in individual Erwinia herbicola cells, we were able to reconstruct the contribution of individual immigrants to bacterial population sizes on leaves. Our analysis revealed that plant foliage represents to bacteria an environment where individual fate is determined by the local carrying capacity of the site where an immigrant cell lands. With increasing inoculation densities, the RS of most immigrants declined, suggesting that local carrying capacity under the tested conditions was linked to local nutrient availability. Fitting the observed experimental data to an adapted model of phyllosphere colonization indicated that there might exist three types of sites on leaves, which differ in their frequency of occurrence and local carrying capacity. Specifically, our data were consistent with a leaf environment that is characterized by few sites where individual immigrants can produce high numbers of offspring, whereas the remainder of the leaf offered an equal number of sites with low and medium RS. Our findings contribute to a bottom–up understanding of bacterial colonization of leaf surfaces, which includes a quantifiable role of chance in the experience at the individual level and in the outcome at the population level.     

The Association of Genotype-Based Inbreeding Coefficient with a Range of Physical and Psychological Human Traits

Across animal species, offspring of closely related mates exhibit lower fitness, a phenomenon called inbreeding depression. Inbreeding depression in humans is less well understood because mating between close relatives is generally rare and stigmatised, confounding investigation of its effect on fitness-relevant traits. Recently, the availability of high-density genotype data has enabled quantification of variation in distant inbreeding in ‘outbred’ human populations, but the low variance of inbreeding detected from genetic data in most outbred populations means large samples are required to test effects, and only a few traits have yet been studied. However, it is likely that isolated populations, or those with a small effective population size, have higher variation in inbreeding and therefore require smaller sample sizes to detect inbreeding effects. With a small effective population size and low immigration, Northern Finland is such a population. We make use of a sample of ∼5,500 ‘unrelated’ individuals in the Northern Finnish Birth Cohort 1966 with known genotypes and measured phenotypes across a range of fitness-relevant physical and psychological traits, including birth length and adult height, body mass index (BMI), waist-to-hip ratio, blood pressure, heart rate, grip strength, educational attainment, income, marital status, handedness, health, and schizotypal features. We find significant associations in the predicted direction between individuals'' inbreeding coefficient (measured by proportion of the genome in runs of homozygosity) and eight of the 18 traits investigated, significantly more than the one or two expected by chance. These results are consistent with inbreeding depression effects on a range of human traits, but further research is needed to replicate and test alternative explanations for these effects.     

Plasmid carriage can limit bacteria–phage coevolution

Coevolution with bacteriophages is a major selective force shaping bacterial populations and communities. A variety of both environmental and genetic factors has been shown to influence the mode and tempo of bacteria–phage coevolution. Here, we test the effects that carriage of a large conjugative plasmid, pQBR103, had on antagonistic coevolution between the bacterium Pseudomonas fluorescens and its phage, SBW25ϕ2. Plasmid carriage limited bacteria–phage coevolution; bacteria evolved lower phage-resistance and phages evolved lower infectivity in plasmid-carrying compared with plasmid-free populations. These differences were not explained by effects of plasmid carriage on the costs of phage resistance mutations. Surprisingly, in the presence of phages, plasmid carriage resulted in the evolution of high frequencies of mucoid bacterial colonies. Mucoidy can provide weak partial resistance against SBW25ϕ2, which may have limited selection for qualitative resistance mutations in our experiments. Taken together, our results suggest that plasmids can have evolutionary consequences for bacteria that go beyond the direct phenotypic effects of their accessory gene cargo.     

Using Mahalanobis distance to compare genomic signatures between bacterial plasmids and chromosomes

Plasmids are ubiquitous mobile elements that serve as a pool of many host beneficial traits such as antibiotic resistance in bacterial communities. To understand the importance of plasmids in horizontal gene transfer, we need to gain insight into the ‘evolutionary history’ of these plasmids, i.e. the range of hosts in which they have evolved. Since extensive data support the proposal that foreign DNA acquires the host's nucleotide composition during long-term residence, comparison of nucleotide composition of plasmids and chromosomes could shed light on a plasmid's evolutionary history. The average absolute dinucleotide relative abundance difference, termed δ-distance, has been commonly used to measure differences in dinucleotide composition, or ‘genomic signature’, between bacterial chromosomes and plasmids. Here, we introduce the Mahalanobis distance, which takes into account the variance–covariance structure of the chromosome signatures. We demonstrate that the Mahalanobis distance is better than the δ-distance at measuring genomic signature differences between plasmids and chromosomes of potential hosts. We illustrate the usefulness of this metric for proposing candidate long-term hosts for plasmids, focusing on the virulence plasmids pXO1 from Bacillus anthracis, and pO157 from Escherichia coli O157:H7, as well as the broad host range multi-drug resistance plasmid pB10 from an unknown host.     

Anaerobic ammonia oxidation in a fertilized paddy soil

Evidence for anaerobic ammonium oxidation in a paddy field was obtained in Southern China using an isotope-pairing technique, quantitative PCR assays and 16S rRNA gene clone libraries, along with nutrient profiles of soil cores. A paddy field with a high load of slurry manure as fertilizer was selected for this study and was shown to contain a high amount of ammonium (6.2–178.8 mg kg⁻¹). The anaerobic oxidation of ammonium (anammox) rates in this paddy soil ranged between 0.5 and 2.9 nmolN per gram of soil per hour in different depths of the soil core, and the specific cellular anammox activity observed in batch tests ranged from 2.9 to 21 fmol per cell per day. Anammox contributed 4–37% to soil N2 production, the remainder being due to denitrification. The 16S rRNA gene sequences of surface soil were closely related to the anammox bacteria ‘Kuenenia'', ‘Anammoxoglobus'' and ‘Jettenia''. Most of the anammox 16S rRNA genes retrieved from the deeper soil were affiliated to ‘Brocadia''. The retrieval of mainly bacterial amoA sequences in the upper part of the paddy soil indicated that nitrifying bacteria may be the major source of nitrite for anammox bacteria in the cultivated horizon. In the deeper oxygen-limited parts, only archaeal amoA sequences were found, indicating that archaea may produce nitrite in this part of the soil. It is estimated that a total loss of 76 g N m⁻² per year is linked to anammox in the paddy field.     

Pervasive sign epistasis between conjugative plasmids and drug-resistance chromosomal mutations

Multidrug-resistant bacteria arise mostly by the accumulation of plasmids and chromosomal mutations. Typically, these resistant determinants are costly to the bacterial cell. Yet, recently, it has been found that, in Escherichia coli bacterial cells, a mutation conferring resistance to an antibiotic can be advantageous to the bacterial cell if another antibiotic-resistance mutation is already present, a phenomenon called sign epistasis. Here we study the interaction between antibiotic-resistance chromosomal mutations and conjugative (i.e., self-transmissible) plasmids and find many cases of sign epistasis (40%)--including one of reciprocal sign epistasis where the strain carrying both resistance determinants is fitter than the two strains carrying only one of the determinants. This implies that the acquisition of an additional resistance plasmid or of a resistance mutation often increases the fitness of a bacterial strain already resistant to antibiotics. We further show that there is an overall antagonistic interaction between mutations and plasmids (52%). These results further complicate expectations of resistance reversal by interdiction of antibiotic use.     

The bacterial and mitochondrial ribosomal A-site molecular switches possess different conformational substates

The A site of the small ribosomal subunit participates in the fidelity of decoding by switching between two states, a resting ‘off’ state and an active decoding ‘on’ state. Eight crystal structures of RNA duplexes containing two minimal decoding A sites of the Homo sapiens mitochondrial wild-type, the A1555G mutant or bacteria have been solved. The resting ‘off’ state of the mitochondrial wild-type A site is surprisingly different from that of the bacterial A site. The mitochondrial A1555G mutant has two types of the ‘off’ states; one is similar to the mitochondrial wild-type ‘off’ state and the other is similar to the bacterial ‘off’ state. Our present results indicate that the dynamics of the A site in bacteria and mitochondria are different, a property probably related to the small number of tRNAs used for decoding in mitochondria. Based on these structures, we propose a hypothesis for the molecular mechanism of non-syndromic hearing loss due to the mitochondrial A1555G mutation.     

Evolution of body shape in sympatric versus non-sympatric Tropheus populations of Lake Tanganyika

Allopatric speciation often yields ecologically equivalent sister species, so that their secondary admixis enforces competition. The shores of Lake Tanganyika harbor about 120 distinct populations of the cichlid genus Tropheus, but only some are sympatric. When alone, Tropheus occupies a relatively broad depth zone, but in sympatry, fish segregate by depth. To assess the effects of competition, we studied the partial co-occurrence of Tropheus moorii ‘Kaiser'' and ‘Kirschfleck'' with Tropheus polli. A previous study demonstrated via standardized breeding experiments that some observed differences between Tropheus ‘Kaiser'' living alone and in sympatry with T. polli have a genetic basis despite large-scale phenotypic plasticity. Using geometric morphometrics and neutral genetic markers, we now investigated whether sympatric populations differ consistently in body shape from populations living alone and if the differences are adaptive. We found significant differences in mean shape between non-sympatric and sympatric populations, whereas all sympatric populations of both color morphs clustered together in shape space. Sympatric populations had a relatively smaller head, smaller eyes and a more anterior insertion of the pectoral fin than non-sympatric populations. Genetically, however, non-sympatric and sympatric ‘Kaiser'' populations clustered together to the exclusion of ‘Kirschfleck''. Genetic distances, but not morphological distances, were correlated with geographic distances. Within- and between-population covariance matrices for T. moorii populations deviated from proportionality. It is thus likely that natural selection acts on both phenotypic plasticity and heritable traits and that both factors contribute to the observed shape differences. The consistency of the pattern in five populations suggests ecological character displacement.     

When Phase Contrast Fails: ChainTracer and NucTracer,Two ImageJ Methods for Semi-Automated Single Cell Analysis Using Membrane or DNA Staining

Within bacterial populations, genetically identical cells often behave differently. Single-cell measurement methods are required to observe this heterogeneity. Flow cytometry and fluorescence light microscopy are the primary methods to do this. However, flow cytometry requires reasonably strong fluorescence signals and is impractical when bacteria grow in cell chains. Therefore fluorescence light microscopy is often used to measure population heterogeneity in bacteria. Automatic microscopy image analysis programs typically use phase contrast images to identify cells. However, many bacteria divide by forming a cross-wall that is not detectable by phase contrast. We have developed ‘ChainTracer’, a method based on the ImageJ plugin ObjectJ. It can automatically identify individual cells stained by fluorescent membrane dyes, and measure fluorescence intensity, chain length, cell length, and cell diameter. As a complementary analysis method we developed ''NucTracer'', which uses DAPI stained nucleoids as a proxy for single cells. The latter method is especially useful when dealing with crowded images. The methods were tested with Bacillus subtilis and Lactococcus lactis cells expressing a GFP-reporter. In conclusion, ChainTracer and NucTracer are useful single cell measurement methods when bacterial cells are difficult to distinguish with phase contrast.     

Single Cell Analysis of a Bacterial Sender-Receiver System

Monitoring gene expression dynamics on the single cell level provides important information on cellular heterogeneity and stochasticity, and potentially allows for more accurate quantitation of gene expression processes. We here study bacterial senders and receivers genetically engineered with components of the quorum sensing system derived from Aliivibrio fischeri on the single cell level using microfluidics-based bacterial chemostats and fluorescence video microscopy. We track large numbers of bacteria over extended periods of time, which allows us to determine bacterial lineages and filter out subpopulations within a heterogeneous population. We quantitatively determine the dynamic gene expression response of receiver bacteria to varying amounts of the quorum sensing inducer N-3-oxo-C6-homoserine lactone (AHL). From this we construct AHL response curves and characterize gene expression dynamics of whole bacterial populations by investigating the statistical distribution of gene expression activity over time. The bacteria are found to display heterogeneous induction behavior within the population. We therefore also characterize gene expression in a homogeneous bacterial subpopulation by focusing on single cell trajectories derived only from bacteria with similar induction behavior. The response at the single cell level is found to be more cooperative than that obtained for the heterogeneous total population. For the analysis of systems containing both AHL senders and receiver cells, we utilize the receiver cells as ‘bacterial sensors’ for AHL. Based on a simple gene expression model and the response curves obtained in receiver-only experiments, the effective AHL concentration established by the senders and their ‘sending power’ is determined.     

Amelioration of the cost of conjugative plasmid carriage in Eschericha coli K12

Although plasmids can provide beneficial functions to their host bacteria, they might confer a physiological or energetic cost. This study examines how natural selection may reduce the cost of carrying conjugative plasmids with drug-resistance markers in the absence of antibiotic selection. We studied two plasmids, R1 and RP4, both of which carry multiple drug resistance genes and were shown to impose an initial fitness cost on Escherichia coli. To determine if and how the cost could be reduced, we subjected plasmid-containing bacteria to 1100 generations of evolution in batch cultures. Analysis of the evolved populations revealed that plasmid loss never occurred, but that the cost was reduced through genetic changes in both the plasmids and the bacteria. Changes in the plasmids were inferred by the demonstration that evolved plasmids no longer imposed a cost on their hosts when transferred to a plasmid-free clone of the ancestral E. coli. Changes in the bacteria were shown by the lowered cost when the ancestral plasmids were introduced into evolved bacteria that had been cured of their (evolved) plasmids. Additionally, changes in the bacteria were inferred because conjugative transfer rates of evolved R1 plasmids were lower in the evolved host than in the ancestral host. Our results suggest that once a conjugative bacterial plasmid has invaded a bacterial population it will remain even if the original selection is discontinued.     

Explaining Bacterial Dispersion on Leaf Surfaces with an Individual-Based Model (PHYLLOSIM)

We developed the individual-based model PHYLLOSIM to explain observed variation in the size of bacterial clusters on plant leaf surfaces (the phyllosphere). Specifically, we tested how different ‘waterscapes’ impacted the diffusion of nutrients from the leaf interior to the surface and the growth of individual bacteria on these nutrients. In the ‘null’ model or more complex ‘patchy’ models, the surface was covered with a continuous water film or with water drops of equal or different volumes, respectively. While these models predicted the growth of individual bacterial immigrants into clusters of variable sizes, they were unable to reproduce experimentally derived, previously published patterns of dispersion which were characterized by a much larger variation in cluster sizes and a disproportionate occurrence of clusters consisting of only one or two bacteria. The fit of model predictions to experimental data was about equally poor (<5%) regardless of whether the water films were continuous or patchy. Only by allowing individual bacteria to detach from developing clusters and re-attach elsewhere to start a new cluster, did PHYLLOSIM come much closer to reproducing experimental observations. The goodness of fit including detachment increased to about 70–80% for all waterscapes. Predictions of this ‘detachment’ model were further supported by the visualization and quantification of bacterial detachment and attachment events at an agarose-water interface. Thus, both model and experiment suggest that detachment of bacterial cells from clusters is an important mechanism underlying bacterial exploration of the phyllosphere.