Evolutionary forces shaping genomic islands of population differentiation in humans

Background Levels of differentiation among populations depend both on demographic and selective factors: genetic drift and local adaptation increase population differentiation, which is eroded by gene flow and balancing selection. We describe here the genomic distribution and the properties of genomic regions with unusually high and low levels of population differentiation in humans to assess the influence of selective and neutral processes on human genetic structure. Methods Individual SNPs of the Human Genome Diversity Panel (HGDP) showing significantly high or low levels of population differentiation were detected under a hierarchical-island model (HIM). A Hidden Markov Model allowed us to detect genomic regions or islands of high or low population differentiation. Results Under the HIM, only 1.5% of all SNPs are significant at the 1% level, but their genomic spatial distribution is significantly non-random. We find evidence that local adaptation shaped high-differentiation islands, as they are enriched for non-synonymous SNPs and overlap with previously identified candidate regions for positive selection. Moreover there is a negative relationship between the size of islands and recombination rate, which is stronger for islands overlapping with genes. Gene ontology analysis supports the role of diet as a major selective pressure in those highly differentiated islands. Low-differentiation islands are also enriched for non-synonymous SNPs, and contain an overly high proportion of genes belonging to the 'Oncogenesis' biological process. Conclusions Even though selection seems to be acting in shaping islands of high population differentiation, neutral demographic processes might have promoted the appearance of some genomic islands since i) as much as 20% of islands are in non-genic regions ii) these non-genic islands are on average two times shorter than genic islands, suggesting a more rapid erosion by recombination, and iii) most loci are strongly differentiated between Africans and non-Africans, a result consistent with known human demographic history.

Invasion and fixation of sex-reversal genes

We simulated a meta-population with random dispersal among demes but local mating within demes to investigate conditions under which a dominant female-determining gene W, with no individual selection advantage, can invade and become fixed in females, changing the population from male to female heterogamety. Starting with one mutant W in a single deme, the interaction of sex ratio selection and random genetic drift causes W to be fixed among females more often than a comparable neutral mutation with no influence on sex determination, even when YY males have slightly reduced viability. Meta-population structure and interdeme selection can also favour the fixation of W. The reverse transition from female to male heterogamety can also occur with higher probability than for a comparable neutral mutation. These results help to explain the involvement of sex-determining genes in the evolution of sex chromosomes and in sexual selection and speciation.

DEVELOPMENT AND APPLICATIONS OF GENOMIC RESOURCES FOR TWO NORTH AMERICAN HARDWOOD TAXA

Hardwoods comprise about half of the biomass of forestlands in North America and present many uses including economic, ecological and aesthetic functions. Forest trees rely on the genetic variation within tree populations to overcome the many biotic, abiotic, anthropogenic factors which are further worsened by climate change, that threaten their continued survival and functionality. To harness these inherent genetic variations of tree populations, informed knowledge of the genomic resources and techniques, which are currently lacking or very limited, are imperative for forest managers. The current study therefore aimed to develop genomic microsatellite markers for the leguminous tree species, honey locust, Gleditsia triacanthos L. and test their applicability in assessing genetic variation, estimation of gene flow patterns and identification of a full-sib mapping population. We also aimed to test the usefulness of already developed nuclear and gene-based microsatellite markers in delineation of species and taxonomic relationships between four of the taxonomically difficult Section Lobatae species (Quercus coccinea, Q. ellipsoidalis, Q. rubra and Q. velutina. We recorded 100% amplification of G. triacanthos genomic microsatellites developed using Illumina sequencing techniques in a panel of seven unrelated individuals with 14 of these showing high polymorphism and reproducibility. When characterized in 36 natural population samples, we recorded 20 alleles per locus with no indication for null alleles at 13 of the 14 microsatellites. This is the first report of genomic microsatellites for this species. Honey locust trees occur in fragmented populations of abandoned farmlands and pastures and is described as essentially dioecious. Pollen dispersal if the main source of gene flow within and between populations with the ability to offset the effects of random genetic drift. Factors known to influence gene include fragmentation and degree of isolation, which make the patterns gene flow in fragmented populations of honey locust a necessity for their sustainable management. In this follow-up study, we used a subset of nine of the 14 developed gSSRs to estimate gene flow and identify a full-sib mapping population in two isolated fragments of honey locust. Our analyses indicated that the majority of the seedlings (65-100% - at both strict and relaxed assignment thresholds) were sired by pollen from outside the two fragment populations. Only one selfing event was recorded confirming the functional dioeciousness of honey locust and that the seed parents are almost completely outcrossed. From the Butternut Valley, TN population, pollen donor genotypes were reconstructed and used in paternity assignment analyses to identify a relatively large full-sib family comprised of 149 individuals, proving the usefulness of isolated forest fragments in identification of full-sib families. In the Ames Plantation stand, contemporary pollen dispersal followed a fat-tailed exponential-power distribution, an indication of effective gene flow. Our estimate of δ was 4,282.28 m, suggesting that insect pollinators of honey locust disperse pollen over very long distances. The high proportion of pollen influx into our sampled population implies that our fragment population forms part of a large effectively reproducing population. The high tendency of oak species to hybridize while still maintaining their species identity make it difficult to resolve their taxonomic relationships. Oaks of the section Lobatae are famous in this regard and remain unresolved at both morphological and genetic markers. We applied 28 microsatellite markers including outlier loci with potential roles in reproductive isolation and adaptive divergence between species to natural populations of four known interfertile red oaks, Q. coccinea, Q. ellpsoidalis, Q. rubra and Q. velutina. To better resolve the taxonomic relationships in this difficult clade, we assigned individual samples to species, identified hybrids and introgressive forms and reconstructed phylogenetic relationships among the four species after exclusion of genetically intermediate individuals. Genetic assignment analyses identified four distinct species clusters, with Q. rubra most differentiated from the three other species, but also with a comparatively large number of misclassified individuals (7.14%), hybrids (7.14%) and introgressive forms (18.83%) between Q. ellipsoidalis and Q. velutina. After the exclusion of genetically intermediate individuals, Q. ellipsoidalis grouped as sister species to the largely parapatric Q. coccinea with high bootstrap support (91 %). Genetically intermediate forms in a mixed species stand were located proximate to both potential parental species, which supports recent hybridization of Q. velutina with both Q. ellipsoidalis and Q. rubra. Analyses of genome-wide patterns of interspecific differentiation can provide a better understanding of speciation processes and taxonomic relationships in this taxonomically difficult group of red oak species.

Adaptive rather than non-adaptive evolution of Mimulus guttatus in its invasive range

Adaptive and non-adaptive evolutionary processes are likely to play important roles in biological invasions but their relative importance has hardly ever been quantified. Moreover, although genetic differences between populations in their native versus invasive ranges may simply reflect different positions along a genetic latitudinal cline, this has rarely been controlled for. To study non-adaptive evolutionary processes in invasion of Mimulus guttatus, we used allozyme analyses on offspring of seven native populations from western North America, and three and four invasive populations from Scotland and New Zealand, respectively. To study quantitative genetic differentiation, we grew 2474 plants representing 17 native populations and the seven invasive populations in a common greenhouse environment under temporarily and permanently wet soil conditions. The absence of allozyme differentiation between the invasive and native range indicates that multiple genotypes had been introduced to Scotland and New Zealand, and suggests that founder effects and genetic drift played small, if any, roles in shaping genetic structure of invasive M. guttatus populations. Plants from the invasive and native range did not differ in phenology, floral traits and sexual and vegetative reproduction, and also not in plastic responses to the watering treatments. However, plants from the invasive range produced twice as many flower-bearing upright side branches than the ones from the native populations. Further, with increasing latitude of collection, vegetative reproduction of our experimental plants increased while sexual reproduction decreased. Plants from the invasive and native range shared these latitudinal clines. Because allozymes showed that the relatedness between native and invasive populations did not depend on latitude, this suggests that plants in the invasive regions have adapted to the local latitude. Overall, our study indicates that quantitative genetic variation of M. guttatus in its two invasive regions is shaped by adaptive evolutionary processes rather than by non-adaptive ones. (C) 2007 Gesellschaft fur Okologie. Published by Elsevier GmbH. All rights reserved.

Divergent evolutionary processes associated with colonization of offshore islands

Oceanic islands have been a test ground for evolutionary theory, but here, we focus on the possibilities for evolutionary study created by offshore islands. These can be colonized through various means and by a wide range of species, including those with low dispersal capabilities. We use morphology, modern and ancient sequences of cytochrome b (cytb) and microsatellite genotypes to examine colonization history and evolutionary change associated with occupation of the Orkney archipelago by the common vole (Microtus arvalis), a species found in continental Europe but not in Britain. Among possible colonization scenarios, our results are most consistent with human introduction at least 5100 bp (confirmed by radiocarbon dating). We used approximate Bayesian computation of population history to infer the coast of Belgium as the possible source and estimated the evolutionary timescale using a Bayesian coalescent approach. We showed substantial morphological divergence of the island populations, including a size increase presumably driven by selection and reduced microsatellite variation likely reflecting founder events and genetic drift. More surprisingly, our results suggest that a recent and widespread cytb replacement event in the continental source area purged cytb variation there, whereas the ancestral diversity is largely retained in the colonized islands as a genetic ‘ark’. The replacement event in the continental M. arvalis was probably triggered by anthropogenic causes (land-use change). Our studies illustrate that small offshore islands can act as field laboratories for studying various evolutionary processes over relatively short timescales, informing about the mainland source area as well as the island.

Adaptive evolution during an ongoing range expansion: the invasive bank vole (Myodes glareolus) in Ireland

Range expansions are extremely common, but have only recently begun to attract attention in terms of their genetic consequences. As populations expand, demes at the wave front experience strong genetic drift, which is expected to reduce genetic diversity and potentially cause ‘allele surfing’, where alleles may become fixed over a wide geographical area even if their effects are deleterious. Previous simulation models show that range expansions can generate very strong selective gradients on dispersal, reproduction, competition and immunity. To investigate the effects of range expansion on genetic diversity and adaptation, we studied the population genomics of the bank vole (Myodes glareolus) in Ireland. The bank vole was likely introduced in the late 1920s and is expanding its range at a rate of ~2.5 km/year. Using genotyping-by-sequencing, we genotyped 281 bank voles at 5979 SNP loci. Fourteen sample sites were arranged in three transects running from the introduction site to the wave front of the expansion. We found significant declines in genetic diversity along all three transects. However, there was no evidence that sites at the wave front had accumulated more deleterious mutations. We looked for outlier loci with strong correlations between allele frequency and distance from the introduction site, where the direction of correlation was the same in all three transects. Amongst these outliers, we found significant enrichment for genic SNPs, suggesting the action of selection. Candidates for selection included several genes with immunological functions and several genes that could influence behaviour.

Genomics of Divergence along a Continuum of Parapatric Population Differentiation

The patterns of genomic divergence during ecological speciation are shaped by a combination of evolutionary forces. Processes such as genetic drift, local reduction of gene flow around genes causing reproductive isolation, hitchhiking around selected variants, variation in recombination and mutation rates are all factors that can contribute to the heterogeneity of genomic divergence. On the basis of 60 fully sequenced three-spined stickleback genomes, we explore these different mechanisms explaining the heterogeneity of genomic divergence across five parapatric lake and river population pairs varying in their degree of genetic differentiation. We find that divergent regions of the genome are mostly specific for each population pair, while their size and abundance are not correlated with the extent of genome-wide population differentiation. In each pair-wise comparison, an analysis of allele frequency spectra reveals that 25–55% of the divergent regions are consistent with a local restriction of gene flow. Another large proportion of divergent regions (38–75%) appears to be mainly shaped by hitchhiking effects around positively selected variants. We provide empirical evidence that alternative mechanisms determining the evolution of genomic patterns of divergence are not mutually exclusive, but rather act in concert to shape the genome during population differentiation, a first necessary step towards ecological speciation.

Spatiotemporal dynamics of Puumala hantavirus associated with its rodent host, Myodes glareolus

Many viruses significantly impact human and animal health. Understanding the population dynamics of these viruses and their hosts can provide important insights for epidemiology and virus evolution. Puumala virus (PUUV) is a European hantavirus that may cause regional outbreaks of hemorrhagic fever with renal syndrome in humans. Here, we analyzed the spatiotemporal dynamics of PUUV circulating in local populations of its rodent reservoir host, the bank vole (Myodes glareolus) during eight years. Phylogenetic and population genetic analyses of all three genome segments of PUUV showed strong geographical structuring at a very local scale. There was a high temporal turnover of virus strains in the local bank vole populations, but several virus strains persisted through multiple years. Phylodynamic analyses showed no significant changes in the local effective population sizes of PUUV, although vole numbers and virus prevalence fluctuated widely. Microsatellite data demonstrated also a temporally persisting subdivision between local vole populations, but these groups did not correspond to the subdivision in the virus strains. We conclude that restricted transmission between vole populations and genetic drift play important roles in shaping the genetic structure and temporal dynamics of PUUV in its natural host which has several implications for zoonotic risks of the human population.

Abies – Circum-Mediterranean firs in Europe: distribution, habitat, usage and threats

Most European firs occur predominantly in small to medium-sized populations in the Mediterranean region, sometimes with fragmented and limited distributions, except for silver fir (Abies alba). They all are genetically closely related and can easily hybridise, perhaps as a consequence of late speciation during the late Quaternary. Circum-Mediterranean firs occur principally in mountain areas with medium to high precipitations rates which are mostly concentrated during the winter period. The species are able to tolerate long droughts in summer and tend to form pure stands when in optimal habitats. In the past firs have been extensively logged for construction and fire wood and their stands were replaced by other more disturbance adapted species or converted into rural areas. Nowadays with the exception of silver fir and Caucasian fir (Abies nordmanniana), circum-Mediterranean firs do not have a wide commercial interest. In Turkey they are still exploited for timber wood, while other firs have an ornamental use in gardening. Great importance is given to their preservation, especially to those populations which have very limited areas and specimens, with the creation of protected reserves and conservation programmes. Wild fires, livestock grazing and genetic drift represent actually their main threats.

MATHEMATICAL STUDIES ON THE EVOLUTIONARY CHANGE OF DNA SEQUENCES

With the aim of understanding the mechanism of molecular evolution, mathematical problems on the evolutionary change of DNA sequences are studied. The problems studied and the results obtained are as follows: (1) Estimation of evolutionary distance between nucleotide sequences. Studying the pattern of nucleotide substitution for the case of unequal substitution rates, a new mathematical formula for estimating the average number of nucleotide substitutions per site between two homologous DNA sequences is developed. It is shown that this formula has a wider applicability than currently available formulae. A statistical method for estimating the number of nucleotide changes due to deletion and insertion is also developed. (2) Biases of the estimates of nucleotide substitutions obtained by the restriction enzyme method. The deviation of the estimate of nucleotide substitutions obtained by the restriction enzyme method from the true value is investigated theoretically. It is shown that the amount of the deviation depends on the nucleotides in the recognition sequence of the restriction enzyme used, unequal rates of substitution among different nucleotides, and nucleotide frequences, but the primary factor is the unequal rates of nucleotide substitution. When many different kinds of enzymes are used, however, the amount of average deviation is generally small. (3) Distribution of restriction fragment lengths. To see the effect of undetectable restriction fragments and fragment differences on the estimate of nucleotide differences, the theoretical distribution of fragment lengths is studied. This distribution depends on the type of restriction enzymes used as well as on the relative frequencies of four nucleotides. It is shown that undetectability of small fragments or fragment differences gives a serious underestimate of nucleotide substitutions when the length-difference method of estimation is used, but the extent of underestimation is small when the site-difference method is used. (4) Evolutionary relationships of DNA sequences in finite populations. A mathematical theory on the expected evolutionary relationships among DNA sequences (nucleons) randomly chosen from the same or different populations is developed under the assumption that the evolutionary change of nucleons is determined solely by mutation and random genetic drift. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author). UMI ^

Spatial and temporal aspects of synaptic plasticity

The notion that changes in synaptic efficacy underlie learning and memory processes is now widely accepted even if definitive proof of the synaptic plasticity and memory hypothesis is still lacking. When learning occurs, patterns of neural activity representing the occurrence of events cause changes in the strength of synaptic connections within the brain. Reactivation of these altered connections constitutes the experience of memory for these events and for other events with which they may be associated. These statements summarize a long-standing theory of memory formation that we refer to as the synaptic plasticity and memory hypothesis. Since activity-dependent synaptic plasticity is induced at appropriate synapses during memory formation, and is both necessary and sufficient for the information storage, we can speculate that a methodological study of the synapse will help us understand the mechanism of learning. Random events underlie a wide range of biological processes as diverse as genetic drift and molecular diffusion, regulation of gene expression and neural network function. Additionally spatial variability may be important especially in systems with nonlinear behavior. Since synapse is a complex biological system we expect that stochasticity as well as spatial gradients of different enzymes may be significant for induction of plasticity. ^ In that study we address the question "how important spatial and temporal aspects of synaptic plasticity may be". We developed methods to justify our basic assumptions and examined the main sources of variability of calcium dynamics. Among them, a physiological method to estimate the number of postsynaptic receptors as well as a hybrid algorithm for simulating postsynaptic calcium dynamics. Additionally we studied how synaptic geometry may enhance any possible spatial gradient of calcium dynamics and how that spatial variability affect plasticity curves. Finally, we explored the potential of structural synaptic plasticity to provide a metaplasticity mechanism specific for the synapse. ^

Individual and population-level responses to ocean acidification

Ocean acidification is predicted to have detrimental effects on many marine organisms and ecological processes. Despite growing evidence for direct impacts on specific species, few studies have simultaneously considered the effects of ocean acidification on individuals (e.g. consequences for energy budgets and resource partitioning) and population level demographic processes. Here we show that ocean acidification increases energetic demands on gastropods resulting in altered energy allocation, i.e. reduced shell size but increased body mass. When scaled up to the population level, long-term exposure to ocean acidification altered population demography, with evidence of a reduction in the proportion of females in the population and genetic signatures of increased variance in reproductive success among individuals. Such increased variance enhances levels of short-term genetic drift which is predicted to inhibit adaptation. Our study indicates that even against a background of high gene flow, ocean acidification is driving individual- and population-level changes that will impact eco-evolutionary trajectories.

Reconstructing the Timing and Dispersion Routes of HIV-1 Subtype B Epidemics in The Caribbean and Central America: A Phylogenetic Story.

The Caribbean and Central America are among the regions with highest HIV-1B prevalence worldwide. Despite of this high virus burden, little is known about the timing and the migration patterns of HIV-1B in these regions. Migration is one of the major processes shaping the genetic structure of virus populations. Thus, reconstruction of epidemiological network may contribute to understand HIV-1B evolution and reduce virus prevalence. We have investigated the spatio-temporal dynamics of the HIV-1B epidemic in The Caribbean and Central America using 1,610 HIV-1B partial pol sequences from 13 Caribbean and 5 Central American countries. Timing of HIV-1B introduction and virus evolutionary rates, as well as the spatial genetic structure of the HIV-1B populations and the virus migration patterns were inferred. Results revealed that in The Caribbean and Central America most of the HIV-1B variability was generated since the 80 s. At odds with previous data suggesting that Haiti was the origin of the epidemic in The Caribbean, our reconstruction indicated that the virus could have been disseminated from Puerto Rico and Antigua. These two countries connected two distinguishable migration areas corresponding to the (mainly Spanish-colonized) Easter and (mainly British-colonized) Western islands, which indicates that virus migration patterns are determined by geographical barriers and by the movement of human populations among culturally related countries. Similar factors shaped the migration of HIV-1B in Central America. The HIV-1B population was significantly structured according to the country of origin, and the genetic diversity in each country was associated with the virus prevalence in both regions, which suggests that virus populations evolve mainly through genetic drift. Thus, our work contributes to the understanding of HIV-1B evolution and dispersion pattern in the Americas, and its relationship with the geography of the area and the movements of human populations.

Distribution of haplotypes from a chromosome 21 region distinguishes multiple prehistoric human migrations

Despite mounting genetic evidence implicating a recent origin of modern humans, the elucidation of early migratory gene-flow episodes remains incomplete. Geographic distribution of haplotypes may show traces of ancestral migrations. However, such evolutionary signatures can be erased easily by recombination and mutational perturbations. A 565-bp chromosome 21 region near the MX1 gene, which contains nine sites frequently polymorphic in human populations, has been found. It is unaffected by recombination and recurrent mutation and thus reflects only migratory history, genetic drift, and possibly selection. Geographic distribution of contemporary haplotypes implies distinctive prehistoric human migrations: one to Oceania, one to Asia and subsequently to America, and a third one predominantly to Europe. The findings with chromosome 21 are confirmed by independent evidence from a Y chromosome phylogeny. Loci of this type will help to decipher the evolutionary history of modern humans.

Population dynamics, demographic stochasticity, and the evolution of cooperation

A basic evolutionary problem posed by the Iterated Prisoner’s Dilemma game is to understand when the paradigmatic cooperative strategy Tit-for-Tat can invade a population of pure defectors. Deterministically, this is impossible. We consider the role of demographic stochasticity by embedding the Iterated Prisoner’s Dilemma into a population dynamic framework. Tit-for-Tat can invade a population of defectors when their dynamics exhibit short episodes of high population densities with subsequent crashes and long low density periods with strong genetic drift. Such dynamics tend to have reddened power spectra and temporal distributions of population size that are asymmetric and skewed toward low densities. The results indicate that ecological dynamics are important for evolutionary shifts between adaptive peaks.