Rapid range expansion increases genetic differentiation while causing limited reduction in genetic diversity in a damselfly

Many ectothermic species are currently expanding their geographic range due to global warming. This can modify the population genetic diversity and structure of these species because of genetic drift during the colonization of new areas. Although the genetic signatures of historical range expansions have been investigated in an array of species, the genetic consequences of natural, contemporary range expansions have received little attention, with the only studies available focusing on range expansions along a narrow front. We investigate the genetic consequences of a natural range expansion in the Mediterranean damselfly Coenagrion scitulum, which is currently rapidly expanding along a broad front in different directions. We assessed genetic diversity and genetic structure using 12 microsatellite markers in five centrally located populations and five recently established populations at the edge of the geographic distribution. Our results suggest that, although a marginal significant decrease in the allelic richness was found in the edge populations, genetic diversity has been preserved during the range expansion of this species. Nevertheless, edge populations were genetically more differentiated compared with core populations, suggesting genetic drift during the range expansion. The smaller effective population sizes of the edge populations compared with central populations also suggest a contribution of genetic drift after colonization. We argue and document that range expansion along multiple axes of a broad expansion front generates little reduction in genetic diversity, yet stronger differentiation of the edge populations.     

FAN1’s protection against CGG repeat expansion requires its nuclease activity and is FANCD2-independent

The Repeat Expansion Diseases, a large group of human diseases that includes the fragile X-related disorders (FXDs) and Huntington''s disease (HD), all result from expansion of a disease-specific microsatellite via a mechanism that is not fully understood. We have previously shown that mismatch repair (MMR) proteins are required for expansion in a mouse model of the FXDs, but that the FANCD2 and FANCI associated nuclease 1 (FAN1), a component of the Fanconi anemia (FA) DNA repair pathway, is protective. FAN1’s nuclease activity has been reported to be dispensable for protection against expansion in an HD cell model. However, we show here that in a FXD mouse model a point mutation in the nuclease domain of FAN1 has the same effect on expansion as a null mutation. Furthermore, we show that FAN1 and another nuclease, EXO1, have an additive effect in protecting against MSH3-dependent expansions. Lastly, we show that the loss of FANCD2, a vital component of the Fanconi anemia DNA repair pathway, has no effect on expansions. Thus, FAN1 protects against MSH3-dependent expansions without diverting the expansion intermediates into the canonical FA pathway and this protection depends on FAN1 having an intact nuclease domain.     

Species range expansion by beneficial mutations

A species' range can be limited when there is no genetic variation for a trait that allows for adaptation to more extreme environments. We study how range expansion occurs by the establishment of a new mutation that affects a quantitative trait in a spatially continuous population. The optimal phenotype for the trait varies linearly in space. The survival probabilities of new mutations affecting the trait are found by simulation. Shallow environmental gradients favour mutations that arise nearer to the range margin and that have smaller phenotypic effects than do steep gradients. Mutations that become established in shallow environmental gradients typically result in proportionally larger range expansions than those that establish in steep gradients. Mutations that become established in populations with high maximum growth rates tend to originate nearer to the range edge and to cause relatively smaller range expansion than mutations that establish in populations with low maximum growth rates. Under plausible parameter values, mutations that allow for range expansion tend to have large phenotypic effects (more than one phenotypic standard deviation) and cause substantial range expansions (15% or more). Sexual reproduction allows for larger range expansions and adaptation to more extreme environments than asexual reproduction.     

Evidence of rapid change in genetic structure and diversity during range expansion in a recovering large terrestrial carnivore

Recovery of natural populations occurs often with simultaneous or subsequent range expansions. According to population genetic theory, genetic structuring emerges at the expansion front together with decreasing genetic diversity, owing to multiple founder events. Thereupon, as the expansion proceeds and connectivity among populations is established, homogenization and a resurgence of genetic diversity are to be expected. Few studies have used a fine temporal scale combined with genetic sampling to track range expansions as they proceed in wild animal populations. As a natural experiment, the historical eradication of large terrestrial carnivores followed by their recovery and recolonization may facilitate empirical tests of these ideas. Here, using brown bear (Ursus arctos) as model species, we tested predictions from genetic theory of range expansion. Individuals from all over Finland were genotyped for every year between 1996 and 2010 using 12 validated autosomal microsatellite markers. A latitudinal shift of about 110 km was observed in the distribution and delineation of genetic clusters during this period. As the range expansion proceeded, we found, as theory predicts, that the degree of genetic structure decreased, and that both genetic variation and admixture increased. The genetic consequences of range expansions may first be detected after multiple generations, but we found major changes in genetic composition after just 1.5 generations, accompanied by population growth and increased migration. These rapid genetic changes suggest an ongoing concerted action of geographical and demographic expansion combined with substantial immigration of bears from Russia during the recovery of brown bears within the large ecosystem of northern Europe.     

Large C9orf72 Hexanucleotide Repeat Expansions Are Seen in Multiple Neurodegenerative Syndromes and Are More Frequent Than Expected in the UK Population

Hexanucleotide repeat expansions in C9orf72 are a major cause of frontotemporal lobar degeneration (FTLD) and amyotrophic lateral sclerosis (ALS). Understanding the disease mechanisms and a method for clinical diagnostic genotyping have been hindered because of the difficulty in estimating the expansion size. We found 96 repeat-primed PCR expansions: 85/2,974 in six neurodegenerative diseases cohorts (FTLD, ALS, Alzheimer disease, sporadic Creutzfeldt-Jakob disease, Huntington disease-like syndrome, and other nonspecific neurodegenerative disease syndromes) and 11/7,579 (0.15%) in UK 1958 birth cohort (58BC) controls. With the use of a modified Southern blot method, the estimated expansion range (smear maxima) in cases was 800–4,400. Similarly, large expansions were detected in the population controls. Differences in expansion size and morphology were detected between DNA samples from tissue and cell lines. Of those in whom repeat-primed PCR detected expansions, 68/69 were confirmed by blotting, which was specific for greater than 275 repeats. We found that morphology in the expansion smear varied among different individuals and among different brain regions in the same individual. Expansion size correlated with age at clinical onset but did not differ between diagnostic groups. Evidence of instability of repeat size in control families, as well as neighboring SNP and microsatellite analyses, support multiple expansion events on the same haplotype background. Our method of estimating the size of large expansions has potential clinical utility. C9orf72-related disease might mimic several neurodegenerative disorders and, with potentially 90,000 carriers in the United Kingdom, is more common than previously realized.     

On the accumulation of deleterious mutations during range expansions

We investigate the effect of spatial range expansions on the evolution of fitness when beneficial and deleterious mutations cosegregate. We perform individual‐based simulations of 1D and 2D range expansions and complement them with analytical approximations for the evolution of mean fitness at the edge of the expansion. We find that deleterious mutations accumulate steadily on the wave front during range expansions, thus creating an expansion load. Reduced fitness due to the expansion load is not restricted to the wave front, but occurs over a large proportion of newly colonized habitats. The expansion load can persist and represent a major fraction of the total mutation load for thousands of generations after the expansion. The phenomenon of expansion load may explain growing evidence that populations that have recently expanded, including humans, show an excess of deleterious mutations. To test the predictions of our model, we analyse functional genetic diversity in humans and find patterns that are consistent with our model.     

Features of trinucleotide repeat instability in vivo

Unstable repeats are associated with various types of cancer and have been implicated in more than 40 neurode-generative disorders. Trinucleotide repeats are located in non-coding and coding regions of the genome. Studies of bacteria, yeast, mice and man have helped to unravel some features of the mechanism of trinucleotide expansion. Looped DNA structures comprising trinucleotide repeats are processed during replication and/or repair to generate deletions or expansions. Most in vivo data are consistent with a model in which expansion and deletion occur by different mechanisms. In mammals, microsatellite instability is complex and appears to be influenced by genetic, epigenetic and developmental factors.     

Chloroplast microsatellites reveal colonization and metapopulation dynamics in the Canary Island pine

Chloroplast microsatellites are becoming increasingly popular markers for population genetic studies in plants, but there has been little focus on their potential for demographic inference. In this work the utility of chloroplast microsatellites for the study of population expansions was explored. First, we investigated the power of mismatch distribution analysis and the F(S) test with coalescent simulations of different demographic scenarios. We then applied these methods to empirical data obtained for the Canary Island pine (Pinus canariensis). The results of the simulations showed that chloroplast microsatellites are sensitive to sudden population growth. The power of the F(S) test and accuracy of demographic parameter estimates, such as the time of expansion, were reduced proportionally to the level of homoplasy within the data. The analysis of Canary Island pine chloroplast microsatellite data indicated population expansions for almost all sample localities. Demographic expansions at the island level can be explained by the colonization of the archipelago by the pine, while population expansions of different ages in different localities within an island could be the result of local extinctions and recolonization dynamics. Comparable mitochondrial DNA sequence data from a parasite of P. canariensis, the weevil Brachyderes rugatus, supports this scenario, suggesting a key role for volcanism in the evolution of pine forest communities in the Canary Islands.     

Genetic variation in original and colonizing Drosophila buzzatii populations analysed by microsatellite loci isolated with a new PCR screening method

A new polymerase chain reaction-based screening method for microsatellites is presented. Using this method, we isolated 12 microsatellite loci from Drosophila buzzatii, two of which were X-linked. We applied the other 10 microsatellite loci to the analysis of genetic variation in five natural populations of D. buzzatii. Two populations were from the species' original distribution in Argentina, whereas the other three were from Europe (two) and Australia that were colonized 200 and 65 years ago, respectively. Allelic variation was much larger in the original populations than in the colonizing ones and there was a tendency to decreased heterozygosity in the colonizing populations. We used three different statistical procedures for detecting population bottlenecks. All procedures suggested that the low variability in the populations in the Old World was not the result of the recent population decline, but was due to a founder effect followed by a population expansion. In fact, one procedure which detects population expansions and declines based on the genealogical history of microsatellite data suggested that an expansion had taken place in all the colonized populations.     

Interruptions in the expanded ATTCT repeat of spinocerebellar ataxia type 10: repeat purity as a disease modifier?

Spinocerebellar ataxia type 10 (SCA10) is one of numerous genetic disorders that result from simple repeat expansions. SCA10 is caused by expansion of an intronic ATTCT pentanucleotide repeat tract. It is clinically characterized by progressive ataxia, seizures, and anticipation, which can vary within and between families. We report two SCA10 families showing distinct frequencies of seizures and correlations of repeat length with age at onset. One family displayed uninterrupted ATTCT expansions, whereas the other showed multiple interruptions of the repeat by nonconsensus repeat units, which differed both in the length and/or sequence of the repeat unit. Disease-causing microsatellite expansions have been assumed to be composed of uninterrupted pure repeats. Our findings for SCA10 challenge this convention and suggest that the purity of the expanded repeat element may be a disease modifier.     

Characterization of Demographic Expansions From Pairwise Comparisons of Linked Microsatellite Haplotypes

This work extends the methods of demographic inference based on the distribution of pairwise genetic differences between individuals (mismatch distribution) to the case of linked microsatellite data. Population genetics theory describes the distribution of mutations among a sample of genes under different demographic scenarios. However, the actual number of mutations can rarely be deduced from DNA polymorphisms. The inclusion of mutation models in theoretical predictions can improve the performance of statistical methods. We have developed a maximum-pseudolikelihood estimator for the parameters that characterize a demographic expansion for a series of linked loci evolving under a stepwise mutation model. Those loci would correspond to DNA polymorphisms of linked microsatellites (such as those found on the Y chromosome or the chloroplast genome). The proposed method was evaluated with simulated data sets and with a data set of chloroplast microsatellites that showed signal for demographic expansion in a previous study. The results show that inclusion of a mutational model in the analysis improves the estimates of the age of expansion in the case of older expansions.     

On the genealogy of a duplicated microsatellite

When a microsatellite locus is duplicated in a diploid organism, a single pair of PCR primers may amplify as many as four distinct alleles. To study the evolution of a duplicated microsatellite, we consider a coalescent model with symmetric stepwise mutation. Conditional on the time of duplication and a mutation rate, both in a model of completely unlinked loci and in a model of completely linked loci, we compute the probabilities for a sampled diploid individual to amplify one, two, three, or four distinct alleles with one pair of microsatellite PCR primers. These probabilities are then studied to examine the nature of their dependence on the duplication time and the mutation rate. The mutation rate is observed to have a stronger effect than the duplication time on the four probabilities, and the unlinked and linked cases are seen to behave similarly. Our results can be useful for helping to interpret genetic variation at microsatellite loci in species with a very recent history of gene and genome duplication.     

Demographic expansion of parasitic nematodes of livestock based on mitochondrial DNA regions that conflict with the infinite-sites model

Mitochondrial ND4 sequences of populations of four species of parasitic nematodes of livestock were subjected to demographic analyses. Deviation from selective neutrality was detectable using the frequency spectrum of segregating sites and highly negative neutrality statistics. However, the mitochondrial data sets do not comply with the infinite-sites model that underlies these tests, and as a consequence, it was not established whether these features are solely a result of population expansion, or whether aspects of the molecular evolution of these mitochondrial regions are also involved. Coalescent analyses based on Fu's Fs neutrality test, which incorporated estimates of rate heterogeneity, the transition-transversion ratio and nucleotide bias, as well as analyses that are fairly robust to deviations from the infinite-sites model supported population expansion. Also analyses that do not depend on the infinite-sites model suggested historical population expansion of these nematodes. The very similar time since expansion, the absence of signatures of positive selection in ND4 and the logical association with human demography imply that selective sweeps of mitochondrial variants are less probable, and that expansion is the most likely scenario for the parasitic nematodes of livestock. The methods used to characterize the expansion have different assumptions and emphasize different aspects of expansions. The resulting restrictions on the interpretation of expansions are outlined.     

Instability of the fragile X syndrome repeat in mice: the effect of age, diet and mutations in genes that affect DNA replication, recombination and repair proficiency

Repeat expansion diseases such as fragile X syndrome (FXS) result from increases in the size of a specific tandem repeat array. In addition to large expansions, small changes in repeat number and deletions are frequently seen in FXS pedigrees. No mouse model accurately recapitulates all aspects of this instability, particularly the occurrence of large expansions. This may be due to differences between mice and humans in CIS and/or TRANS-acting factors that affect repeat stability. The identification of such factors may help reveal the expansion mechanism and allow the development of suitable animal models for these disorders. We have examined the effect of age, dietary folate, and mutations in the Werner's syndrome helicase (WRN) and TRP53 genes on FXS repeat instability in mice. WRN facilitates replication of the FXS repeat and enhances Okazaki fragment processing, thereby reducing the incidence of processes that have been suggested to lead to expansion. p53 is a protein involved in DNA damage surveillance and repair. We find two types of repeat instability in these mice, small changes in repeat number that are seen at frequencies approaching 100%, and large deletions which occur at a frequency of about 10%. The frequency of these events was independent of WRN, p53, parental age, or folate levels. The large deletions occur at the same frequency in mice homozygous and heterozygous for the repeat suggesting that they are not the result of an interallelic recombination event. In addition, no evidence of large expansions was seen. Our data thus show that the absence of repeat expansions in mice is not due to a more efficient WRN protein or p53-mediated error correction mechanism, and suggest that these proteins, or the pathways in which they are active, may not be involved in expansion in humans either. Moreover, the fact that contractions occur in the absence of expansions suggests that these processes occur by different mechanisms.     

Niche expansion leads to small-scale adaptive divergence along an elevation gradient in a medium-sized passerine bird

Niche expansion can lead to adaptive differentiation and speciation, but there are few examples from contemporary niche expansions about how this process is initiated. We assess the consequences of a niche expansion by Mexican jays (Aphelocoma ultramarina) along an elevation gradient. We predicted that jays at high elevation would have straighter bills adapted to feeding on pine cones, whereas jays at low elevation would have hooked bills adapted to feeding on acorns. We measured morphological and genetic variation of 95 adult jays and found significant differences in hook length between elevations in accordance with predictions, a pattern corroborated by analysis at the regional scale. Genetic results from microsatellite and mtDNA variation support phenotypic differentiation in the presence of gene flow coupled with weak, but detectable genetic differentiation between high- and low-elevation populations. These results demonstrate that niche expansion can lead to adaptive divergence despite gene flow between parapatric populations along an elevation gradient, providing information on a key precursor to ecological speciation.     

ATR protects the genome against CGG.CCG-repeat expansion in Fragile X premutation mice

Fragile X mental retardation syndrome is a repeat expansion disease caused by expansion of a CGG·CCG-repeat tract in the 5′ UTR of the FMR1 gene. In humans, small expansions occur more frequently on paternal transmission while large expansions are exclusively maternal in origin. It has been suggested that expansion is the result of aberrant DNA replication, repair or recombination. To distinguish amongst these possibilities we crossed mice containing 120 CGG·CCG-repeats in the 5′ UTR of the mouse Fmr1 gene to mice with mutations in ATR, a protein important in the cellular response to stalled replication forks and bulky DNA lesions. We show here that ATR heterozygosity results in increased expansion rates of maternally, but not paternally, transmitted alleles. In addition, age-related somatic expansions occurred in mice of both genders that were not seen in ATR wild-type animals. Some ATR-sensitive expansion occurs in postmitotic cells including haploid gametes suggesting that aberrant DNA repair is responsible. Our data suggest that two mechanisms of repeat expansion exist that may explain the small and large expansions seen in humans. In addition, our data provide an explanation for the maternal bias of large expansions in humans and the lower incidence of these expansions in mice.     

Genetic footprints reveal geographic patterns of expansion in Fennoscandian red foxes

Population expansions of boreal species are among the most substantial ecological consequences of climate change, potentially transforming both structure and processes of northern ecosystems. Despite their importance, little is known about expansion dynamics of boreal species. Red foxes (Vulpes vulpes) are forecasted to become a keystone species in northern Europe, a process stemming from population expansions that began in the 19th century. To identify the relative roles of geographic and demographic factors and the sources of northern European red fox population expansion, we genotyped 21 microsatellite loci in modern and historical (1835–1941) Fennoscandian red foxes. Using Bayesian clustering and Bayesian inference of migration rates, we identified high connectivity and asymmetric migration rates across the region, consistent with source‐sink dynamics, whereby more recently colonized sampling regions received immigrants from multiple sources. There were no clear clines in allele frequency or genetic diversity as would be expected from a unidirectional range expansion from south to north. Instead, migration inferences, demographic models and comparison to historical red fox genotypes suggested that the population expansion of the red fox is a consequence of dispersal from multiple sources, as well as in situ demographic growth. Together, these findings provide a rare glimpse into the anatomy of a boreal range expansion and enable informed predictions about future changes in boreal communities.     

Dynamics of microsatellite divergence under stepwise mutation and proportional slippage/point mutation models

Recently Kruglyak, Durrett, Schug, and Aquadro showed that microsatellite equilibrium distributions can result from a balance between polymerase slippage and point mutations. Here, we introduce an elaboration of their model that keeps track of all parts of a perfect repeat and a simplification that ignores point mutations. We develop a detailed mathematical theory for these models that exhibits properties of microsatellite distributions, such as positive skewness of allele lengths, that are consistent with data but are inconsistent with the predictions of the stepwise mutation model. We use our theoretical results to analyze the successes and failures of the genetic distances (delta(mu))(2) and D(SW) when used to date four divergences: African vs. non-African human populations, humans vs. chimpanzees, Drosophila melanogaster vs. D. simulans, and sheep vs. cattle. The influence of point mutations explains some of the problems with the last two examples, as does the fact that these genetic distances have large stochastic variance. However, we find that these two features are not enough to explain the problems of dating the human-chimpanzee split. One possible explanation of this phenomenon is that long microsatellites have a mutational bias that favors contractions over expansions.     

Joint inference of microsatellite mutation models, population history and genealogies using transdimensional Markov Chain Monte Carlo

We provide a framework for Bayesian coalescent inference from microsatellite data that enables inference of population history parameters averaged over microsatellite mutation models. To achieve this we first implemented a rich family of microsatellite mutation models and related components in the software package BEAST. BEAST is a powerful tool that performs Bayesian MCMC analysis on molecular data to make coalescent and evolutionary inferences. Our implementation permits the application of existing nonparametric methods to microsatellite data. The implemented microsatellite models are based on the replication slippage mechanism and focus on three properties of microsatellite mutation: length dependency of mutation rate, mutational bias toward expansion or contraction, and number of repeat units changed in a single mutation event. We develop a new model that facilitates microsatellite model averaging and Bayesian model selection by transdimensional MCMC. With Bayesian model averaging, the posterior distributions of population history parameters are integrated across a set of microsatellite models and thus account for model uncertainty. Simulated data are used to evaluate our method in terms of accuracy and precision of estimation and also identification of the true mutation model. Finally we apply our method to a red colobus monkey data set as an example.     

Northern range expansion of European populations of the wasp spider Argiope bruennichi is associated with global warming–correlated genetic admixture and population‐specific temperature adaptations

Poleward range expansions are observed for an increasing number of species, which may be an effect of global warming during the past decades. However, it is still not clear in how far these expansions reflect simple geographical shifts of species ranges, or whether new genetic adaptations play a role as well. Here, we analyse the expansion of the wasp spider Argiope bruennichi into Northern Europe during the last century. We have used a range‐wide sampling of contemporary populations and historical specimens from museums to trace the phylogeography and genetic changes associated with the range shift. Based on the analysis of mitochondrial, microsatellite and SNP markers, we observe a higher level of genetic diversity in the expanding populations, apparently due to admixture of formerly isolated lineages. Using reciprocal transplant experiments for testing overwintering tolerance, as well as temperature preference and tolerance tests in the laboratory, we find that the invading spiders have possibly shifted their temperature niche. This may be a key adaptation for survival in Northern latitudes. The museum samples allow a reconstruction of the invasion's genetic history. A first, small‐scale range shift started around 1930, in parallel with the onset of global warming. A more massive invasion of Northern Europe associated with genetic admixture and morphological changes occurred in later decades. We suggest that the latter range expansion into far Northern latitudes may be a consequence of the admixture that provided the genetic material for adaptations to new environmental regimes. Hence, global warming could have facilitated the initial admixture of populations and this resulted in genetic lineages with new habitat preferences.