Identifying adaptive genetic divergence among populations from genome scans

The identification of signatures of natural selection in genomic surveys has become an area of intense research, stimulated by the increasing ease with which genetic markers can be typed. Loci identified as subject to selection may be functionally important, and hence (weak) candidates for involvement in disease causation. They can also be useful in determining the adaptive differentiation of populations, and exploring hypotheses about speciation. Adaptive differentiation has traditionally been identified from differences in allele frequencies among different populations, summarised by an estimate of F-ST. Low outliers relative to an appropriate neutral population-genetics model indicate loci subject to balancing selection, whereas high outliers suggest adaptive (directional) selection. However, the problem of identifying statistically significant departures from neutrality is complicated by confounding effects on the distribution of F-ST estimates, and current methods have not yet been tested in large-scale simulation experiments. Here, we simulate data from a structured population at many unlinked, diallelic loci that are predominantly neutral but with some loci subject to adaptive or balancing selection. We develop a hierarchical-Bayesian method, implemented via Markov chain Monte Carlo (MCMC), and assess its performance in distinguishing the loci simulated under selection from the neutral loci. We also compare this performance with that of a frequentist method, based on moment-based estimates of F-ST. We find that both methods can identify loci subject to adaptive selection when the selection coefficient is at least five times the migration rate. Neither method could reliably distinguish loci under balancing selection in our simulations, even when the selection coefficient is twenty times the migration rate.

Quantitative genetic divergence and standing genetic (co)variance in thermal reaction norms along latitude

Although the potential to adapt to warmer climate is constrained by genetic trade-offs, our understanding of how selection and mutation shape genetic (co)variances in thermal reaction norms is poor. Using 71 isofemale lines of the fly Sepsis punctum, originating from northern, central, and southern European climates, we tested for divergence in juvenile development rate across latitude at five experimental temperatures. To investigate effects of evolutionary history in different climates on standing genetic variation in reaction norms, we further compared genetic (co)variances between regions. Flies were reared on either high or low food resources to explore the role of energy acquisition in determining genetic trade-offs between different temperatures. Although the latter had only weak effects on the strength and sign of genetic correlations, genetic architecture differed significantly between climatic regions, implying that evolution of reaction norms proceeds via different trajectories at high latitude versus low latitude in this system. Accordingly, regional genetic architecture was correlated to region-specific differentiation. Moreover, hot development temperatures were associated with low genetic variance and stronger genetic correlations compared to cooler temperatures. We discuss the evolutionary potential of thermal reaction norms in light of their underlying genetic architectures, evolutionary histories, and the materialization of trade-offs in natural environments.

Speciation as an active force in promoting genetic evolution

There is a growing appreciation among evolutionary biologists that the rate and tempo of molecular evolution might often be altered at or near the time of speciation, i.e. that speciation is in some way a special time for genes. Molecular phylogenies frequently reveal increased rates of genetic evolution associated with speciation and other lines of investigation suggest that various types of abrupt genomic disruption can play an important role in promoting speciation via reproductive isolation. These phenomena are in conflict with the gradual view of molecular evolution that is implicit in much of our thinking about speciation and in the tools of modern biology. This raises the prospect of studying the molecular evolutionary consequences of speciation per se and studying the footprint of speciation as an active force in promoting genetic divergence. Here we discuss the reasons to believe that speciation can play such a role and elaborate on possible mechanisms for accelerated rates of evolution following speciation. We provide an example of how it is possible detect whether accelerated bursts of evolution occur in neutral and/or adaptive regions of genes and discuss the implications of rapid episodes of change for conventional models of molecular evolution. Speciation might often owe more to ephemeral and essentially arbitrary events that cause reproductive isolation than to the gradual and regular tug of natural selection that draws a species into a new niche.

Molecular genetic differentiation in earthworms inhabiting a heterogeneous Pb-polluted landscape

A Pb-mine site situated on acidic soil, but comprising of Ca-enriched islands around derelict buildings was used to study the spatial pattern of genetic diversity in Lumbricus rubellus. Two distinct genetic lineages ('A' and 'B'), differentiated at both the mitochondrial (mtDNA COII) and nuclear level (AFLPs) were revealed with a mean inter-lineage mtDNA sequence divergence of approximately 13%, indicative of a cryptic species complex. AFLP analysis indicates that lineage A individuals within one central 'ecological island' site are uniquely clustered, with little genetic overlap with lineage A individuals at the two peripheral sites. FTIR microspectroscopy of Pb-sequestering chloragocytes revealed different phosphate profiles in residents of adjacent acidic and calcareous islands. Bioinformatics found over-representation of Ca pathway genes in ESTPb libraries. Subsequent sequencing of a Ca-transport gene, SERCA, revealed mutations in the protein's cytosolic domain. We recommend the mandatory genotyping of all individuals prior to field-based ecotoxicological assays, particularly those using discriminating genomic technologies.

Genetic structure and systematic relationships within the Ophrys fuciflora aggregate (Orchidaceae: Orchidinae): high diversity in Kent and a wind-induced discontinuity bisecting the Adriatic

A recent phylogenetic study based on multiple datasets is used as the framework for a more detailed examination of one of the ten molecularly circumscribed groups identified, the Ophrys fuciflora aggregate. The group is highly morphologically variable, prone to phenotypic convergence, shows low levels of sequence divergence and contains an unusually large proportion of threatened taxa, including the rarest Ophrys species in the UK. The aims of this study were to (a) circumscribe minimum resolvable genetically distinct entities within the O. fuciflora aggregate, and (b) assess the likelihood of gene flow between genetically and geographically distinct entities at the species and population levels. Fifty-five accessions sampled in Europe and Asia Minor from the O. fuciflora aggregate were studied using the AFLP genetic fingerprinting technique to evaluate levels of infraspecific and interspecific genetic variation and to assess genetic relationships between UK populations of O. fuciflora s.s. in Kent and in their continental European and Mediterranean counterparts. The two genetically and geographically distinct groups recovered, one located in England and central Europe and one in south-eastern Europe, are incongruent with current species delimitation within the aggregate as a whole and also within O. fuciflora s.s. Genetic diversity is higher in Kent than in the rest of western and central Europe. Gene flow is more likely to occur between populations in closer geographical proximity than those that are morphologically more similar. Little if any gene flow occurs between populations located in the south-eastern Mediterranean and those dispersed throughout the remainder of the distribution, revealing a genetic discontinuity that runs north-south through the Adriatic. This discontinuity is also evident in other clades of Ophrys and is tentatively attributed to the long-term influence of prevailing winds on the long-distance distribution of pollinia and especially seeds. A cline of gene flow connects populations from Kent and central and southern Europe; these individuals should therefore be considered part of an extensive meta-population. Gene flow is also evident among populations from Kent, which appear to constitute a single metapopulation. They show some evidence of hybridization, and possibly also introgression, with O. apifera.

Genomic and genetic analyses of diversity and plant interactions of Pseudomonas fluorescens

Background: Pseudomonas fluorescens are common soil bacteria that can improve plant health through nutrient cycling, pathogen antagonism and induction of plant defenses. The genome sequences of strains SBW25 and Pf0-1 were determined and compared to each other and with P. fluorescens Pf-5. A functional genomic in vivo expression technology (IVET) screen provided insight into genes used by P. fluorescens in its natural environment and an improved understanding of the ecological significance of diversity within this species. Results: Comparisons of three P. fluorescens genomes (SBW25, Pf0-1, Pf-5) revealed considerable divergence: 61% of genes are shared, the majority located near the replication origin. Phylogenetic and average amino acid identity analyses showed a low overall relationship. A functional screen of SBW25 defined 125 plant-induced genes including a range of functions specific to the plant environment. Orthologues of 83 of these exist in Pf0-1 and Pf-5, with 73 shared by both strains. The P. fluorescens genomes carry numerous complex repetitive DNA sequences, some resembling Miniature Inverted-repeat Transposable Elements (MITEs). In SBW25, repeat density and distribution revealed 'repeat deserts' lacking repeats, covering approximately 40% of the genome. Conclusions: P. fluorescens genomes are highly diverse. Strain-specific regions around the replication terminus suggest genome compartmentalization. The genomic heterogeneity among the three strains is reminiscent of a species complex rather than a single species. That 42% of plant-inducible genes were not shared by all strains reinforces this conclusion and shows that ecological success requires specialized and core functions. The diversity also indicates the significant size of genetic information within the Pseudomonas pan genome.

Adaptive divergence and speciation among sexual and pseudoviviparous populations of Festuca

Pseudovivipary is an environmentally induced flowering abnormality in which vegetative shoots replace seminiferous (sexual) inflorescences. Pseudovivipary is usually retained in transplantation experiments, indicating that the trait is not solely induced by the growing environment. Pseudovivipary is the defining characteristic of Festuca vivipara, and arguably the only feature separating this species from its closest seminiferous relative, Festuca ovina. We performed phylogenetic and population genetic analysis on sympatric F. ovina and F. vivipara samples to establish whether pseudovivipary is an adaptive trait that accurately defines the separation of genetically distinct Festuca species. Chloroplast and nuclear marker-based analyses revealed that variation at a geographical level can exceed that between F. vivipara and F. ovina. We deduced that F. vivipara is a recent species that frequently arises independently within F. ovina populations and has not accumulated significant genetic differentiation from its progenitor. We inferred local gene flow between the species. We identified one amplified fragment length polymorphism marker that may be linked to a pseudovivipary-related region of the genome, and several other markers provide evidence of regional local adaptation in Festuca populations. We conclude that F. vivipara can only be appropriately recognized as a morphologically and ecologically distinct species; it lacks genetic differentiation from its relatives. This is the first report of a ‘failure in normal flowering development’ that repeatedly appears to be adaptive, such that the trait responsible for species recognition constantly reappears on a local basis.

Progressive erosion of genetic and epigenetic variation in callus-derived cocoa (Theobroma cacao) plants

Relatively little is known about the timing of genetic and epigenetic forms of somaclonal variation arising from callus growth. We surveyed for both types of change in cocoa (Theobroma cacao) plants regenerated from calli of various ages, and also between tissues from the source trees. For genetic change, we used 15 single sequence repeat (SSR) markers from four source trees and from 233 regenerated plants. For epigenetic change, we used 386 methylation-sensitive amplified polymorphism (MSAP) markers on leaf and explant (staminode) DNA from two source trees and on leaf DNA from 114 regenerants. Genetic variation within source trees was limited to one slippage mutation in one leaf. Regenerants were far more variable, with 35% exhibiting at least one mutation. Genetic variation initially accumulated with culture age but subsequently declined. MSAP (epigenetic) profiles diverged between leaf and staminode samples from source trees. Multivariate analysis revealed that leaves from regenerants occupied intermediate eigenspace between leaves and staminodes of source plants but became progressively more similar to source tree leaves with culture age. Statistical analysis confirmed this rather counterintuitive finding that leaves of ‘late regenerants’ exhibited significantly less genetic and epigenetic divergence from source leaves than those exposed to short periods of callus growth.