Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Marine foundation species shifting ranges and genetic baselines

This thesis revealed the most importance factors shaping the distribution, abundance and genetic diversity of four marine foundation species. Environmental conditions, particularly sea temperatures, nutrient availability and ocean waves, played a primary role in shaping the spatial distribution and abundance of populations, acting on scales varying from tens of meters to hundreds of kilometres. Furthermore, the use of Species Distribution Models (SDMs) with biological records of occurrence and high-resolution oceanographic data, allowed predicting species distributions across time. This approach highlighted the role of climate change, particularly when extreme temperatures prevailed during glacial and interglacial periods. These results, when combined with mtDNA and microsatellite genetic variation of populations allowed inferring for the influence of past range dynamics in the genetic diversity and structure of populations. For instance, the Last Glacial Maximum produced important shifts in species ranges, leaving obvious signatures of higher genetic diversities in regions where populations persisted (i.e., refugia). However, it was found that a species’ genetic pool is shaped by regions of persistence, adjacent to others experiencing expansions and contractions. Contradicting expectations, refugia seem to play a minor role on the re(colonization) process of previously eroded populations. In addition, the available habitat area for expanding populations and the inherent mechanisms of species dispersal in occupying available habitats were also found to be fundamental in shaping the distributions of genetic diversity. However, results suggest that the high levels of genetic diversity in some populations do not rule out that they may have experienced strong genetic erosion in the past, a process here named shifting genetic baselines. Furthermore, this thesis predicted an ongoing retraction at the rear edges and extinctions of unique genetic lineages, which will impoverish the global gene pool, strongly shifting the genetic baselines in the future.

Phylogéographie comparée de la souris à pattes blanches et de la souris sylvestre, deux vecteurs de la maladie de Lyme au Québec

Mon étude vise à évaluer la propagation d’une zoonose en émergence au Québec, la maladie de Lyme, en conséquence du réchauffement climatique. Le pathogène responsable de cette infection, Borrelia burgdorferi, est transmis par l’intermédiaire d’une tique parasite, Ixodes scapularis, de plus en plus commune au Québec en raison de l’augmentation de la température moyenne du climat depuis les dernières décennies. Puisque la tique a une capacité de déplacement très restreinte, on s'attend à ce que sa dispersion soit liée à celle de son hôte primaire, soit la souris à pattes blanches (Peromyscus leucopus). Je décrirai donc d’abord les espèces impliquées, leur écologie et leur rôle dans ce système à trois niveaux (hôte/pathogène/vecteur). Puis, à l’aide de séquences d’ADN mitochondrial, je comparerai la phylogéographie des deux principales espèces de souris au Québec, la souris à pattes blanches et la souris sylvestre (P. maniculatus). Des analyses d’arbres et de réseaux d’haplotypes ont révélé des différences significatives dans la structure génétique et ainsi montré que les populations de P. leucopus seraient en expansion dans le sud du Québec. Cette étude nous a finalement permis d’émettre des hypothèses sur le patron d’établissement de la maladie de Lyme au Québec.

Development of Nicarbazin as a Reproductive Inhibitor for Resident Canada Geese

Expanding populations of resident Canada geese that remain in suburban and urban areas year-round often result in increased conflicts with humans. Non-lethal and humane means are needed for managing the size of Canada goose flocks residing near or on airports, golf courses, industrial parks, government sites, and city parks. A side effect of nicarbazin, a veterinary drug used to control coccidiosis in chickens, is decreased egg production and hatching. Exploiting this side effect, studies of nicarbazin for reducing the hatchability of eggs from Canada geese were conducted. An initial study in Coturnix quail verified reduction in hatchability in a species other than chickens. Because plasma nicarbazin was not routinely measured, a study in chickens was conducted to determine the relationship between plasma and egg nicarbazin. A comparative study in chickens, mallards, and Canada geese showed that nicarbazin absorption was lowest in geese. Studies in both penned and wild Canada geese showed that reduction in hatchability was possible but neither study used bait suitable for general field application. Bait development led to the OvoControl-G® (Innolytics LLC) bait, which resulted in reduction in hatchability of 51% at treated sites compared to control sites in the field. Previous studies showed that nicarbazin is practically non-toxic and is environmentally friendly; timing and management of baiting will minimize non-target hazards. OvoControl-G® 2500 ppm nicarbazin bait is recommended for incorporation into a comprehensive management plan as a reproductive inhibitor for use in controlling resident Canada goose flock sizes.

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.

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.

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.

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.

The Genetics of Postglacial Colonization and Peripherality in Northwestern Populations of the Spring Peeper, Pseudacris crucifer

Glaciation over the Pleistocene induced dramatic range fluctuations for species across North America such that postglacial recolonization by southern refugial lineages has characterized the genetic structure of northern North American species. Based on the leading edge model of postglacial range expansion, dispersal and rapid population growth in these northern taxa is expected to produce vast areas of genetic homogeneity. Previous work on the widely distributed spring peeper (Pseudacris crucifer) revealed six distinct mitochondrial lineages that diverged between 3-11 mya, expanding and contracting with glacial cycles. Beginning 16,000 yBP, receding glaciers permitted Eastern lineage refugia residing in the southern Appalachians to migrate northward into the St. Lawrence Valley then westward through most of central Canada. Peripheral populations at the northwestern range limit of P. crucifer in central Manitoba are likely descended from this westward expanding Eastern lineage. According to the central-marginal hypothesis, founder effects from colonization as well as limited gene flow is expected to reveal genetic differentiation and lower genetic diversity in peripheral populations. The goal of my study is to further our understanding of peripheral range dynamics in peripheral Manitoba populations of P. crucifer by determining their genetic affinity and diversity relative to more central populations in Ontario and Minnesota. In this study I amplified and aligned cytochrome b sequences from sample sites across central Manitoba to reconstruct a Bayesian phylogeny for P. crucifer; additionally, microsatellite loci were genotyped to estimate genetic diversity. Results from this study affirmed Eastern lineage descent for peripheral Manitoba sites by aligning with Ontario. Initial colonization by the Interior lineage between glacial retreat and the appearance of arid vicariance events may explain the apparent introgression of non-Eastern lineages in Manitoba. However, genetic diversity measured in expected heterozygosity (H¬e) was not found to be significantly different in Manitoba genotypes. Greater isolation by distance and inbreeding relative to Ontario and Minnesota is likely the primary driver of genetic variation in these sites. Further sampling is necessary to generate a more complete genetic population structure for P. crucifer.

Multiple origins of European populations of giant liver fluke Fascioloides magna (Trematoda: Fasciolidae), an invasive liver parasite of ruminants

The giant liver fluke, Fascioloides magna, liver parasite of free-living and domestic ruminants of Europe and North America, was analysed in order to determine the origin of European populations and to reveal the biogeography of this originally North American parasite on the European continent. The previously selected variable fragments of the mitochondrial cytochrome c oxidase subunit I (cox1; 384 bp) and nicotinamid dehydrogenase subunit I (nad1; 405 bp) were applied as a tool. The phylogenetic trees and haplotype networks were constructed and the level of genetic structuring was evaluated using population genetic tools. In F. magna individuals originating from all European natural foci (Italy, Czech Republic, Danube floodplain forests) and from four of five major North American enzootic areas, 16 cox1 and 18 nad1 haplotypes were determined. The concatenated sequence set produced 22 distinct haplotypes. The European fluke populations were less diverse than those from North America in that they contained proportionately fewer haplotypes (8), while more substantial level of genetic diversity and higher number of haplotypes (15) were recorded in North America. Only one haplotype was shared between the European (Italy) and North American (USA/Oregon and Canada/Alberta) flukes supporting a western North American origin of the Italian F. magna population. Haplotypes found in Italy were distinct from those determined in the remaining European localities what indicates that introduction of F. magna onto the European continent is a result of more than one event. In Czech focus, a south-eastern US origin of giant liver fluke was revealed. Identical haplotypes, common for parasites from Czech Republic and from expanding focus of Danube floodplain forests, implies introduction of F. magna to the Danube region from an already established Czech focus.

Invading and expanding: range dynamics and ecological consequences of the greater white-toothed shrew (Crocidura russula) invasion in Ireland

Establishing how invasive species impact upon pre-existing species is a fundamental question in ecology and conservation biology. The greater white-toothed shrew (Crocidura russula) is an invasive species in Ireland that was first recorded in 2007 and which, according to initial data, may be limiting the abundance/distribution of the pygmy shrew (Sorex minutus), previously Ireland’s only shrew species. Because of these concerns, we undertook an intensive live-trapping survey (and used other data from live-trapping, sightings and bird of prey pellets/nest inspections collected between 2006 and 2013) to model the distribution and expansion of C. russula in Ireland and its impacts on Ireland’s small mammal community. The main distribution range of C. russula was found to be approximately 7,600 km2 in 2013, with established outlier populations suggesting that the species is dispersing with human assistance within the island. The species is expanding rapidly for a small mammal, with a radial expansion rate of 5.5 km/yr overall (2008–2013), and independent estimates from live-trapping in 2012–2013 showing rates of 2.4–14.1 km/yr, 0.5–7.1 km/yr and 0–5.6 km/yr depending on the landscape features present. S. minutus is negatively associated with C. russula. S. minutus is completely absent at sites where C. russula is established and is only present at sites at the edge of and beyond the invasion range of C. russula. The speed of this invasion and the homogenous nature of the Irish landscape may mean that S. minutus has not had sufficient time to adapt to the sudden appearance of C. russula. This may mean the continued decline/disappearance of S. minutus as C. russula spreads throughout the island.

Implications of climate change for expanding cities worldwide

This paper analyses the trends of the changing environmental effects within growing megacities as their diameters exceed 50–100 km and their populations rise beyond 30 million people. The authors consider how these effects are influenced by climate change, to which urban areas themselves contribute, caused by their increasing greenhouse gas emissions associated with rapidly expanding energy use. Other environmental and social factors are assessed, quantitatively and qualitatively, using detailed modelling of urban mesoscale meteorology, which shows how these factors can lead to large conurbations becoming more vulnerable to climatic and environmental hazards. The paper discusses the likely changes in meteorological and hydrological hazards in urban areas, both as the climate changes and the sizes of urban areas grow. Examples are given of how these risks are being reduced through innovations in warning and response systems, planning and infrastructure design, which should include refuges against extreme natural disasters. Policies are shown to be more effective when they are integrated and based on substantial community involvement. Some conclusions are drawn regarding how policies for the natural and artificial environment and for reducing many kinds of climate and hazard risk are related to future designs and planning of infrastructure and open spaces.

Towards global drought early warning capability: expanding international cooperation for the development of a framework for global drought monitoring and forecasting

Drought is a global problem that has far-reaching impacts and especially 47 on vulnerable populations in developing regions. This paper highlights the need for a Global Drought Early Warning System (GDEWS), the elements that constitute its underlying framework (GDEWF) and the recent progress made towards its development. Many countries lack drought monitoring systems, as well as the capacity to respond via appropriate political, institutional and technological frameworks, and these have inhibited the development of integrated drought management plans or early warning systems. The GDEWS will provide a source of drought tools and products via the GDEWF for countries and regions to develop tailored drought early warning systems for their own users. A key goal of a GDEWS is to maximize the lead time for early warning, allowing drought managers and disaster coordinators more time to put mitigation measures in place to reduce the vulnerability to drought. To address this, the GDEWF will take both a top-down approach to provide global real-time drought monitoring and seasonal forecasting, and a bottom-up approach that builds upon existing national and regional systems to provide continental to global coverage. A number of challenges must be overcome, however, before a GDEWS can become a reality, including the lack of in-situ measurement networks and modest seasonal forecast skill in many regions, and the lack of infrastructure to translate data into useable information. A set of international partners, through a series of recent workshops and evolving collaborations, has made progress towards meeting these challenges and developing a global system.