Impact of range expansions on current human genomic diversity

The patterns of population genetic diversity depend to a large extent on past demographic history. Most human populations are known to have gone recently through a series of range expansions within and out of Africa, but these spatial expansions are rarely taken into account when interpreting observed genomic diversity, possibly because they are difficult to model. Here we review available evidence in favour of range expansions out of Africa, and we discuss several of their consequences on neutral and selected diversity, including some recent observations on an excess of rare neutral and selected variants in large samples. We further show that in spatially subdivided populations, the sampling strategy can severely impact the resulting genetic diversity and be confounded by past demography. We conclude that ignoring the spatial structure of human population can lead to some misinterpretations of extant genetic diversity.

Expansion load: recessive mutations and the role of standing genetic variation

Expanding populations incur a mutation burden – the so-called expansion load. Previous studies of expansion load have focused on codominant mutations. An important consequence of this assumption is that expansion load stems exclusively from the accumulation of new mutations occurring in individuals living at the wave front. Using individual-based simulations, we study here the dynamics of standing genetic variation at the front of expansions, and its consequences on mean fitness if mutations are recessive. We find that deleterious genetic diversity is quickly lost at the front of the expansion, but the loss of deleterious mutations at some loci is compensated by an increase of their frequencies at other loci. The frequency of deleterious homozygotes therefore increases along the expansion axis, whereas the average number of deleterious mutations per individual remains nearly constant across the species range. This reveals two important differences to codominant models: (i) mean fitness at the front of the expansion drops much faster if mutations are recessive, and (ii) mutation load can increase during the expansion even if the total number of deleterious mutations per individual remains constant. We use our model to make predictions about the shape of the site frequency spectrum at the front of range expansion, and about correlations between heterozygosity and fitness in different parts of the species range. Importantly, these predictions provide opportunities to empirically validate our theoretical results. We discuss our findings in the light of recent results on the distribution of deleterious genetic variation across human populations and link them to empirical results on the correlation of heterozygosity and fitness found in many natural range expansions.

Changes in morphology and physiology of an East Mediterranean sponge in different habitats

The sponge Tetilla sp. (Tetractinomorpha: Tetillidae) is a common species in the eastern Mediterranean. This sponge inhabits four different habitat types differing in wave impact and irradiance levels. Two of these habitats (a shallow cave and deep water) are characterized by relatively calm water, whereas the other two (shallow exposed site and tide pools) are in turbulent water with high energy flow. The present study examined the influence of physical (depth, illumination and water motion) and biotic factors on morphology, skeletal plasticity and reproductive traits among the four spatially separated populations. Sponges from tidal pools had significantly larger body volume than sponges from deep water and from shallow caves (ANOVA: tidal-deep P< 0.0001; tidal-shallow caves P< 0.05). Sponges from exposed habitats were significantly larger than deep-water sponges (ANOVA: P=0.01). In addition, individuals from tide pools and from the exposed habitat had a significantly higher proportion of structural silica than sponges from the calmer deep water and from the cave sites. Oxea spicules in sponges from the calm habitats were significantly shorter than in those from the tidal pools and the exposed habitats. The percentage of spicules out of a sponge's dry weight in individuals transplanted from deep (calm) to shallow (turbulent) water significantly increased by 21.9&PLUSMN; 12.9%. The new spicule percentage did not differ significantly from that of sponges originally from shallow water. Oocyte diameter differed significantly between habitats. The maximal size of mature eggs was found in deep-water sponges in June (97&PLUSMN; 5 μ m). In the shallow habitats, a smaller maximal oocyte diameter was found in the cave, in May (56.5&PLUSMN; 3 μ m). Furthermore, oocyte density in shallow-water sponges was highest in May and decreased in June (with 88.2&PLUSMN; 9 and 19.3&PLUSMN; 9 oocytes mm(-2), respectively). At the same time (June), oocyte density of deep-water sponges had just reached its maximum (155&PLUSMN; 33.7 oocytes mm(-2)). The difference in oocyte size and density between deep- and shallow-water individuals indicates an earlier gamete release in the shallow sponge population. The results suggest that plasticity in skeletal design of this sponge indicates a trade off between spicule production and investment in reproduction.

Does colonization asymmetry matter in metapopulations?

Despite the considerable evidence showing that dispersal between habitat patches is often asymmetric, most of the metapopulation models assume symmetric dispersal. In this paper, we develop a Monte Carlo simulation model to quantify the effect of asymmetric dispersal on metapopulation persistence. Our results suggest that metapopulation extinctions are more likely when dispersal is asymmetric. Metapopulation viability in systems with symmetric dispersal mirrors results from a mean field approximation, where the system persists if the expected per patch colonization probability exceeds the expected per patch local extinction rate. For asymmetric cases, the mean field approximation underestimates the number of patches necessary for maintaining population persistence. If we use a model assuming symmetric dispersal when dispersal is actually asymmetric, the estimation of metapopulation persistence is wrong in more than 50% of the cases. Metapopulation viability depends on patch connectivity in symmetric systems, whereas in the asymmetric case the number of patches is more important. These results have important implications for managing spatially structured populations, when asymmetric dispersal may occur. Future metapopulation models should account for asymmetric dispersal, while empirical work is needed to quantify the patterns and the consequences of asymmetric dispersal in natural metapopulations.

Dynamique d’une population faunique en expansion sous différents scénarios climatiques : le dindon sauvage (Meleagris gallopavo)

Les changements climatiques récents ont mené à l’expansion de la répartition de plusieurs espèces méridionales, mais ont aussi causé l’extinction locale d’espèces se retrouvant à la limite de leur tolérance environnementale. Ces populations en expansion peuvent favoriser différentes stratégies d’histoire de vie en répondant à différents facteurs limitants. Dans cette thèse, je vise à déterminer et quantifier l’effet du climat et des évènements extrêmes sur le cycle de vie complet d’une espèce en expansion (le dindon sauvage) pour comprendre les changements au niveau populationnel ainsi que les mécanismes impliqués dans l’expansion de la distribution d’une espèce. J’ai défini les évènements extrêmes de pluie, d’épaisseur de neige au sol et de température, comme un évènement dont la fréquence est plus rare que le 10e et 90e percentile. En utilisant l’approche « Measure-Understand-Predict » (MUP), j’ai tout d’abord suivi trois populations le long d’un gradient latitudinal de sévérité hivernale pour mesurer l’effet de variables météorologiques sur la dynamique des populations. La survie des dindons sauvages diminuait drastiquement lorsque l’accumulation de neige au sol dépassait 30 cm pour une période de 10 jours et diminuait également avec la température. Au printemps, la persistance de la neige affectait négativement le taux d’initiation de la nidification et l’augmentation de la pluie diminuait la survie des nids. Dans une deuxième étape, j’ai examiné l’impact des évènements climatiques extrêmes et des processus démographiques impliqués dans l’expansion du dindon, liés à la théorie des histoires de vie pour comprendre la relation entre la dynamique de ces populations en expansions avec le climat. J’ai démontré que la fréquence des évènements extrêmes hivernaux et, d’une façon moins importante, les évènements extrêmes estivaux limitaient l’expansion nordique des dindons sauvages. J’ai appuyé, à l’aide de données empiriques et de modélisation, les hypothèses de la théorie classique des invasions biologiques en montrant que les populations en établissement priorisaient les paramètres reproducteurs tandis que la survie adulte était le paramètre démographique affectant le plus la dynamique des populations bien établies. De plus, les populations les plus au nord étaient composées d’individus plus jeunes ayant une espérance de vie plus faible, mais avaient un potentiel d’accroissement plus élevé que les populations établies, comme le suggère cette théorie. Finalement, j’ai projeté l’impact de la récolte sur la dynamique des populations de même que le taux de croissance de cette espèce en utilisant les conditions climatiques futures projetées par les modèles de l’IPCC. Les populations en établissement avaient un taux de récolte potentiel plus élevé, mais la proportion de mâles adultes, possédant des caractéristiques recherchées par les chasseurs, diminuait plus rapidement que dans les populations établies. Dans le futur, la fréquence des évènements extrêmes de pluie devrait augmenter tandis que la fréquence des évènements extrêmes de température hivernale et d’accumulation de neige au sol devraient diminuer après 2060, limitant probablement l’expansion nordique du dindon sauvage jusqu’en 2100. Cette thèse améliore notre compréhension des effets météorologiques et du climat sur l’expansion de la répartition des espèces ainsi que les mécanismes démographiques impliqués, et nous a permis de prédire la probabilité de l’expansion nordique de la répartition du dindon sauvage en réponse aux changements climatiques.

An R package for simulating metapopulation dynamics and range expansion under environmental change

The metapopulation paradigm is central in ecology and conservation biology to understand the dynamics of spatially-structured populations in fragmented landscapes. Metapopulations are often studied using simulation modelling, and there is an increasing demand of user-friendly software tools to simulate metapopulation responses to environmental change. Here we describe the MetaLandSim R package, mwhich integrates ideas from metapopulation and graph theories to simulate the dynamics of real and virtual metapopulations. The package offers tools to (i) estimate metapopulation parameters from empirical data, (ii) to predict variation in patch occupancy over time in static and dynamic landscapes, either real or virtual, and (iii) to quantify the patterns and speed of metapopulation expansion into empty landscapes. MetaLandSim thus provides detailed information on metapopulation processes, which can be easily combined with land use and climate change scenarios to predict metapopulation dynamics and range expansion for a variety of taxa and ecological systems.

Evaluating Effects of Everglades Restoration on American Crocodile Populations in South Florida Using a Spatially-Explicit, Stage-Based Population Model

The distribution and abundance of the American crocodile (Crocodylus acutus) in the Florida Everglades is dependent on the timing, amount, and location of freshwater flow. One of the goals of the Comprehensive Everglades Restoration Plan (CERP) is to restore historic freshwater flows to American crocodile habitat throughout the Everglades. To predict the impacts on the crocodile population from planned restoration activities, we created a stage-based spatially explicit crocodile population model that incorporated regional hydrology models and American crocodile research and monitoring data. Growth and survival were influenced by salinity, water depth, and density-dependent interactions. A stage-structured spatial model was used with discrete spatial convolution to direct crocodiles toward attractive sources where conditions were favorable. The model predicted that CERP would have both positive and negative impacts on American crocodile growth, survival, and distribution. Overall, crocodile populations across south Florida were predicted to decrease approximately 3 % with the implementation of CERP compared to future conditions without restoration, but local increases up to 30 % occurred in the Joe Bay area near Taylor Slough, and local decreases up to 30 % occurred in the vicinity of Buttonwood Canal due to changes in salinity and freshwater flows.

Corrected mean-field models for spatially-dependent advection-diffusion-reaction phenomena

In the exclusion-process literature, mean-field models are often derived by assuming that the occupancy status of lattice sites is independent. Although this assumption is questionable, it is the foundation of many mean-field models. In this work we develop methods to relax the independence assumption for a range of discrete exclusion process-based mechanisms motivated by applications from cell biology. Previous investigations that focussed on relaxing the independence assumption have been limited to studying initially-uniform populations and ignored any spatial variations. By ignoring spatial variations these previous studies were greatly simplified due to translational invariance of the lattice. These previous corrected mean-field models could not be applied to many important problems in cell biology such as invasion waves of cells that are characterised by moving fronts. Here we propose generalised methods that relax the independence assumption for spatially inhomogeneous problems, leading to corrected mean-field descriptions of a range of exclusion process-based models that incorporate (i) unbiased motility, (ii) biased motility, and (iii) unbiased motility with agent birth and death processes. The corrected mean-field models derived here are applicable to spatially variable processes including invasion wave type problems. We show that there can be large deviations between simulation data and traditional mean-field models based on invoking the independence assumption. Furthermore, we show that the corrected mean-field models give an improved match to the simulation data in all cases considered.

Anthropogenic landscape change promotes asymmetric dispersal and limits regional patch occupancy in a spatially structured bird population

1. Local extinctions in habitat patches and asymmetric dispersal between patches are key processes structuring animal populations in heterogeneous environments. Effective landscape conservation requires an understanding of how habitat loss and fragmentation influence demographic processes within populations and movement between populations. 2. We used patch occupancy surveys and molecular data for a rainforest bird, the logrunner (Orthonyx temminckii), to determine (i) the effects of landscape change and patch structure on local extinction; (ii) the asymmetry of emigration and immigration rates; (iii) the relative influence of local and between-population landscapes on asymmetric emigration and immigration; and (iv) the relative contributions of habitat loss and habitat fragmentation to asymmetric emigration and immigration. 3. Whether or not a patch was occupied by logrunners was primarily determined by the isolation of that patch. After controlling for patch isolation, patch occupancy declined in landscapes experiencing high levels of rainforest loss over the last 100 years. Habitat loss and fragmentation over the last century was more important than the current pattern of patch isolation alone, which suggested that immigration from neighbouring patches was unable to prevent local extinction in highly modified landscapes. 4. We discovered that dispersal between logrunner populations is highly asymmetric. Emigration rates were 39% lower when local landscapes were fragmented, but emigration was not limited by the structure of the between-population landscapes. In contrast, immigration was 37% greater when local landscapes were fragmented and was lower when the between-population landscapes were fragmented. Rainforest fragmentation influenced asymmetric dispersal to a greater extent than did rainforest loss, and a 60% reduction in mean patch area was capable of switching a population from being a net exporter to a net importer of dispersing logrunners. 5. The synergistic effects of landscape change on species occurrence and asymmetric dispersal have important implications for conservation. Conservation measures that maintain large patch sizes in the landscape may promote asymmetric dispersal from intact to fragmented landscapes and allow rainforest bird populations to persist in fragmented and degraded landscapes. These sink populations could form the kernel of source populations given sufficient habitat restoration. However, the success of this rescue effect will depend on the quality of the between-population landscapes.

The wMel Wolbachia strain blocks dengue and invades caged Aedes aegypti populations

Dengue fever is the most important mosquito-borne viral disease of humans with more than 50 million cases estimated annually in more than 100 countries. Disturbingly, the geographic range of dengue is currently expanding and the severity of outbreaks is increasing. Control options for dengue are very limited and currently focus on reducing population abundance of the major mosquito vector, Aedes aegypti. These strategies are failing to reduce dengue incidence in tropical communities and there is an urgent need for effective alternatives. It has been proposed that endosymbiotic bacterial Wolbachia infections of insects might be used in novel strategies for dengue control. For example, the wMelPop-CLA Wolbachia strain reduces the lifespan of adult A. aegypti mosquitoes in stably transinfected lines. This life-shortening phenotype was predicted to reduce the potential for dengue transmission. The recent discovery that several Wolbachia infections, including wMelPop-CLA, can also directly influence the susceptibility of insects to infection with a range of insect and human pathogens has markedly changed the potential for Wolbachia infections to control human diseases. Here we describe the successful transinfection of A. aegypti with the avirulent wMel strain of Wolbachia, which induces the reproductive phenotype cytoplasmic incompatibility with minimal apparent fitness costs and high maternal transmission, providing optimal phenotypic effects for invasion. Under semi-field conditions, the wMel strain increased from an initial starting frequency of 0.65 to near fixation within a few generations, invading A. aegypti populations at an accelerated rate relative to trials with the wMelPop-CLA strain. We also show that wMel and wMelPop-CLA strains block transmission of dengue serotype 2 (DENV-2) in A. aegypti, forming the basis of a practical approach to dengue suppression.

A spatially explicit multi-competitor coexistence model of penaeid (shrimp) distribution on the Australian Great Barrier Reef

A spatially explicit multi-competitor coexistence model was developed for meta-populations of prawns (shrimp) occupying habitat patches across the Great Barrier Reef, where dispersal was localised and dispersal rates varied between species. Prawns were modelled as individuals moving to and from patches or cells according to pre-set decision rules. The landscape was simulated as a matrix of cells with each cell having a spatially explicit survival index for each species. Mixed species prawn assemblages moved over this simplified spatially explicit landscape. A low level of chronic random environmental disturbance was assumed (cyclone and tropical storm damage) with additional acute spatially confined disturbance due to commercial trawling, modelled as an increase in mortality affecting inter-specific competition. The general form of the results was for increased disturbance to favour good-colonising "generalist" species at the expense of good-competitor "specialists". Increasing fishing mortality (local patch extinctions) combined with poor colonising ability resulted in low equilibrium abundance for even the best competitor, while in the same circumstances the poorest competitor but best coloniser could have the highest equilibrium abundance. This mimics the switch from high-value prawn species to lower-value prawn species as trawl effort increases, reflected in historic catch and effort logbook data and reported anecdotaly from the north Queensland trawl fleet. To match the observed distribution and behaviour of prawn assemblages, a combination inter-species competition, a spatially explicit landscape, and a defined pattern of disturbance (trawling) was required. Modelling this combination could simulate not only general trends in spatial distribution of each of prawn species but also localised concentrations observed in the survey data

Effects of fishing on tropical reef associated shark populations on the Great Barrier Reef

Three data sets were examined to define the level of interaction of reef associated sharks with the commercial Coral Reef Fin Fish Fishery within the Great Barrier Reef (GBR). Data were examined from fishery logbooks, an observer program within the fishery and a fishery-independent survey conducted as part of the Effects of Line Fishing (ELF) Experiment. The majority of the identified catch was comprised of grey reef (62-72%), whitetip reef (16-29%) and blacktip reef (6-13%) sharks. Logbook data revealed spatially and temporally variable landings of shark from the GBR. Catch per unit effort (CPUE) through time was stable for the period from 1989 to 2006 with no evidence of increase or decline. Data from observer and ELF data sets indicated no differences in CPUE among regions. The ELF data set demonstrated that CPUE was higher in Marine National Park zones (no fishing) when compared to General Use zones (open to fishing). The ongoing and consistent catches of reef sharks in the fishery and effectiveness of no-fishing zones suggest that management zones within the GBR Marine Park are effective at protecting a portion of the reef shark population from exploitation.

Habitat use and population dynamics of brown bears (Ursus arctos) in Scandinavia

Large carnivore populations are currently recovering from past extirpation efforts and expanding back into their original habitats. At the same time human activities have resulted in very few wilderness areas left with suitable habitats and size large enough to maintain populations of large carnivores without human contact. Consequently the long-term future of large carnivores depends on their successful integration into landscapes where humans live. Thus, understanding their behaviour and interaction with surrounding habitats is of utmost importance in the development of management strategies for large carnivores. This applies also to brown bears (Ursus arctos) that were almost exterminated from Scandinavia and Finland at the turn of the century, but are now expanding their range with the current population estimates being approximately 2600 bears in Scandinavia and 840 in Finland. This thesis focuses on the large-scale habitat use and population dynamics of brown bears in Scandinavia with the objective to develop modelling approaches that support the management of bear populations. Habitat analysis shows that bear home ranges occur mainly in forested areas with a low level of human influence relative to surrounding areas. Habitat modelling based on these findings allows identification and quantification of the potentially suitable areas for bears in Scandinavia. Additionally, this thesis presents novel improvements to home range estimation that enable realistic estimates of the effective area required for the bears to establish a home range. This is achieved through fitting to the radio-tracking data to establish the amount of temporal autocorrelation and the proportion of time spent in different habitat types. Together these form a basis for the landscape-level management of the expanding population. Successful management of bears requires also assessment of the consequences of harvest on the population viability. An individual-based simulation model, accounting for the sexually selected infanticide, was used to investigate the possibility of increasing the harvest using different hunting strategies, such as trophy harvest of males. The results indicated that the population can sustain twice the current harvest rate. However, harvest should be changed gradually while carefully monitoring the population growth as some effects of increased harvest may manifest themselves only after a time-delay. The results and methodological improvements in this thesis can be applied to the Finnish bear population and to other large carnivores. They provide grounds for the further development of spatially-realistic management-oriented models of brow bear dynamics that can make projections of the future distribution of bears while accounting for the development of human activities.

Influence of Water Depth on the Rate of Expansion of Giant Cutgrass Populations and Management Implications

Giant cutgrass ( Zizaniopsis miliacea (Michx.) Doell. & Asch.), a tall emergent grass native to the southeastern United States, was studied in Lake Seminole where it formed large expanding stands, and Lake Alice where it was confined to a stable narrow fringe.

Loss of allele diversity in introduced populations of the hermaphroditic bay scallop Argopecten irradians

The bay scallop Argopecten irradians is a hermaphroditic bivalve native to the Atlantic coast of the United States that was introduced to China for aquaculture production in 1982. It now supports a major aquaculture industry in China. Introduced species often start with limited genetic variability, which is problematic for the further selective breeding. Bay scallop aquaculture is exclusively hatchery based and as the initial introduction consisted of only 26 scallops, there have been concerns about inbreeding and inbreeding depression in cultured populations in China. In this study, eleven simple sequence repeat (SSR) markers were used to compare genetic variation in cultured populations from China with that in a natural population from the east coast of America. Although the difference in heterozygosity was small, the Chinese populations lost 9 of the 45 alleles (20%) found in the wild population. The reduced allele diversity suggests that the Chinese bay scallop populations experienced a bottleneck in genetic diversity that remains significant despite several recent introductions of new stocks aimed at expanding the gene pool. The loss of allele diversity may affect future efforts in selective breeding and domestication, and results of this study highlight the need for additional introductions, advanced breeding programs that minimize inbreeding and continued genetic monitoring. (c) 2007 Elsevier B.V. All rights reserved.