Understanding marine reserve function in a seascape genetics context: Nucella lapillus in Strangford Lough (Northern Ireland) as an example

A conservation priority in the marine environment is the establishment of ecologically coherent reserve networks. Since these networks will integrate existent reserves, an understanding of spatial genetic diversity and genetic connectivities between areas is necessary. Using Strangford Lough marine nature reserve (MNR) as a model, spatial genetic analyses were employed to evaluate the function of the lough. Samples of the marine gastropod Nucella lapillus (L.) from 7 locations in the reserve and adjacent areas were screened at 6 microsatellites. Genetic variation was temporally stable. Significant genetic structuring (F-ST = 0.133) was observed among samples. Genetic divergence and isolation by distance indicated reduced gene flow between the marine reserve and coastal samples relative to that between adjacent coastal samples. Partitioning of genetic variation between the reserve and coast was significant (AMOVA, 7.45%, p

Landscape effects on extremely fragmented populations of a rare solitary bee, Colletes floralis

Globally there is concern over the decline of bees, an ecologically important group of pollinating insects. Genetic studies provide insights into population structure that are crucial for conservation management but that would be impossible to obtain by conventional ecological methods. Yet conservation genetic studies of bees have primarily focussed on social species rather than the more species-rich solitary bees. Here we investigate the population structure of Colletes floralis, a rare and threatened solitary mining bee, in Ireland and Scotland using nine microsatellite loci. Genetic diversity was surprisingly as high in Scottish (Hebridean island) populations at the extreme northwestern edge of the species range as in mainland Irish populations further south. Extremely high genetic differentiation among populations was detected; multilocus FST was up to 0.53, and G’ST and Dest were even higher (maximum: 0.85 and 1.00 respectively). A pattern of isolation by distance was evident for sites separated by land. Water appears to act as a substantial barrier to gene flow yet sites separated by sea did not exhibit isolation by distance. Colletes floralis populations are extremely isolated and probably not in regional migration-drift equilibrium. GIS-based landscape genetic analysis reveals urban areas as a potential and substantial barrier to gene flow. Our results highlight the need for urgent site-specific management action to halt the decline of this and potentially other rare solitary bees.

Genetic differentiation across the social transition in a socially polymorphic sweat bee, Halictus rubicundus

Eusociality is widely considered a major evolutionary transition. The socially polymorphic sweat bee Halictus rubicundus, solitary in cooler regions of its holarctic range and eusocial in warmer parts, is an excellent model organism to address this transition, and specifically the question of whether sociality is associated with a strong barrier to gene flow between phenotypically divergent populations. Mitochondrial DNA (COI) from specimens collected across the British Isles, where both solitary and social phenotypes are represented, displayed limited variation, but placed all specimens in the same European lineage; haplotype network analysis failed to differentiate solitary and social lineages. Microsatellite genetic variability was high and enabled us to quantify genetic differentiation among populations and social phenotypes across Great Britain and Ireland. Results from conceptually different analyses consistently showed greater genetic differentiation between geographically distant populations, independently of their social phenotype, suggesting that the two social forms are not reproductively isolated. A landscape genetic approach revealed significant isolation by distance (Mantel test r = 0.622, p

Levels and patterns of population genetic diversity in the red seaweed Chondrus crispus (Florideophyceae): a direct comparison of single nucleotide polymorphisms and microsatellites

Single nucleotide polymorphisms (SNPs) are predicted to supersede microsatellites as the marker of choice for population genetic studies in the near future. To date, however, very few studies have directly compared both marker systems in natural populations, particularly in non-model organisms. In the present study, we compared the utility of SNPs and microsatellites for population genetic analysis of the red seaweed Chondrus crispus (Florideophyceae). Six SNP loci yielded very different patterns of intrapopulation genetic diversity compared to those obtained using seven moderately (mean 5.2 alleles) polymorphic microsatellite loci, although Bayesian clustering analysis gave largely congruent results between the two marker classes. A weak but significant pattern of isolation-by-distance was observed across scales from a few hundred metres to approximately 200?km using the combined SNP and microsatellite data set of 13 loci. Over larger scales, however, there was little correlation between genetic divergence and geographical distance. Our findings suggest that even a moderate number of SNPs is sufficient to determine patterns of genetic diversity across natural populations, and also highlight the fact that patterns of genetic variation in seaweeds arise through a complex interplay of short- and long-term natural processes, as well as anthropogenic influence.

Complex pattern of genetic structuring in the Atlantic salmon (Salmo salar L.) of the River Foyle system in northwest Ireland: disentangling the evolutionary signal from population stochasticity

Little is known about the microevolutionary processes shaping within river population genetic structure of aquatic organisms characterized by high levels of homing and spawning site fidelity. Using a microsatellite panel, we observed complex and highly significant levels of intrariver population genetic substructure and Isolation-by-Distance, in the Atlantic salmon stock of a large river system. Two evolutionary models have been considered explaining mechanisms promoting genetic substructuring in Atlantic salmon, the member-vagrant and metapopulation models. We show that both models can be simultaneously used to explain patterns and levels of population structuring within the Foyle system. We show that anthropogenic factors have had a large influence on contemporary population structure observed. In an analytical development, we found that the frequently used estimator of genetic differentiation, F-ST, routinely underestimated genetic differentiation by a factor three to four compared to the equivalent statistic Jost's D-est (Jost 2008). These statistics also showed a near-perfect correlation. Despite ongoing discussions regarding the usefulness of "adjusted" F-ST statistics, we argue that these could be useful to identify and quantify qualitative differences between populations, which are important from management and conservation perspectives as an indicator of existence of biologically significant variation among tributary populations or a warning of critical environmental damage.

Understanding macroalgal dispersal in a complex hydrodynamic environment: a combined population genetic and physical modelling approach

Gene flow in macroalgal populations can be strongly influenced by spore or gamete dispersal. This, in turn, is influenced by a convolution of the effects of current flow and specific plant reproductive strategies. Although several studies have demonstrated genetic variability in macroalgal populations over a wide range of spatial scales, the associated current data have generally been poorly resolved spatially and temporally. In this study, we used a combination of population genetic analyses and high-resolution hydrodynamic modelling to investigate potential connectivity between populations of the kelp Laminaria digitata in the Strangford Narrows, a narrow channel characterized by strong currents linking the large semi-enclosed sea lough, Strangford Lough, to the Irish Sea. Levels of genetic structuring based on six microsatellite markers were very low, indicating high levels of gene flow and a pattern of isolation-by-distance, where populations are more likely to exchange migrants with geographically proximal populations, but with occasional long-distance dispersal. This was confirmed by the particle tracking model, which showed that, while the majority of spores settle near the release site, there is potential for dispersal over several kilometres. This combined population genetic and modelling approach suggests that the complex hydrodynamic environment at the entrance to Strangford Lough can facilitate dispersal on a scale exceeding that proposed for L. digitata in particular, and the majority of macroalgae in general. The study demonstrates the potential of integrated physical–biological approaches for the prediction of ecological changes resulting from factors such as anthropogenically induced coastal zone changes.

Use of bacterial whole-genome sequencing to investigate local persistence and spread in bovine tuberculosis

Mycobacterium bovis is the causal agent of bovine tuberculosis, one of the most important diseases currently facing the UK cattle industry. Here, we use high-density whole genome sequencing (WGS) in a defined sub-population of M. bovis in 145 cattle across 66 herd breakdowns to gain insights into local spread and persistence. We show that despite low divergence among isolates, WGS can in principle expose contributions of under-sampled host populations to M. bovis transmission. However, we demonstrate that in our data such a signal is due to molecular type switching, which had been previously undocumented for M. bovis. Isolates from farms with a known history of direct cattle movement between them did not show a statistical signal of higher genetic similarity. Despite an overall signal of genetic isolation by distance, genetic distances also showed no apparent relationship with spatial distance among affected farms over distances <5 km. Using simulations, we find that even over the brief evolutionary timescale covered by our data, Bayesian phylogeographic approaches are feasible. Applying such approaches showed that M. bovis dispersal in this system is heterogeneous but slow overall, averaging 2 km/year. These results confirm that widespread application of WGS to M. bovis will bring novel and important insights into the dynamics of M. bovis spread and persistence, but that the current questions most pertinent to control will be best addressed using approaches that more directly integrate WGS with additional epidemiological data.

Population and landscape genetics of the stone marten and red fox in Portugal : implications for conservation management of common carnivores

Tese de doutoramento, Biologia (Biologia da Conservação), Universidade de Lisboa, Faculdade de Ciências, 2015

Kin structure and queen execution in the Argentine ant Linepithema humile

Every spring, workers of the Argentine Ant Linepithema humile kill a large proportion of queens within their nests, Although this behaviour inflicts a high energetic cost oil the colonies, its biological significance has remained elusive so far. An earlier study showed that the probability of a queen being executed is not related to her weight, fecundity, or age. Here we test the hypothesis that workers collectively eliminate queens to which they are less related, thereby increasing their inclusive fitness. We found no evidence for this hypothesis. Workers of a nest were on average not significantly less related to executed queens than to surviving ones. Moreover, a population genetic analysis revealed that workers were not genetically differentiated between nests. This means that workers of a given nest are equally related to any queen in the population and that there can be no increase in average worker-queen relatedness by selective elimination of queens. Finally, our genetic analyses also showed that, in contrast to workers, queens were significantly genetically differentiated between nests and that there was significant isolation by distance for queens.

Genetic and phenotypic population divergence on a microgeographic scale in brown trout.

Salmonid populations of many rivers are rapidly declining. One possible explanation is that habitat fragmentation increases genetic drift and reduces the populations' potential to adapt to changing environmental conditions. We measured the genetic and eco-morphological diversity of brown trout (Salmo trutta) in a Swiss stream system, using multivariate statistics and Bayesian clustering. We found large genetic and phenotypic variation within only 40 km of stream length. Eighty-eight percent of all pairwise F(ST) comparisons and 50% of the population comparisons in body shape were significant. High success rates of population assignment tests confirmed the distinctiveness of populations in both genotype and phenotype. Spatial analysis revealed that divergence increased with waterway distance, the number of weirs, and stretches of poor habitat between sampling locations, but effects of isolation-by-distance and habitat fragmentation could not be fully disentangled. Stocking intensity varied between streams but did not appear to erode genetic diversity within populations. A lack of association between phenotypic and genetic divergence points to a role of local adaptation or phenotypically plastic responses to habitat heterogeneity. Indeed, body shape could be largely explained by topographic stream slope, and variation in overall phenotype matched the flow regimes of the respective habitats.

Phylogeographie und Populationsgenetik der Arten der Gattung Fosterella (Bromeliaceae) in den südamerikanischen Anden

Die tropischen Anden sind eines der artenreichsten Gebiete der Erde. Fast die Hälfte der 45.000 in diesem Gebiet vorkommenden Gefäßpflanzenarten sind in den Anden endemisch (Myers et al. 2000). Die Gattung Fosterella (Bromeliaceae) ist eine den Anden zugeordnete Pflanzengruppe, denn die meisten ihrer 31 Arten kommen in den Anden vor. Achtzehn Arten sind kleinräumige Endemiten. Fosterella hat damit Modellcharakter für diese Region. In der vorliegenden Arbeit wurde die Evolution der Gattung in Raum und Zeit mithilfe der vergleichenden Sequenzierung von sechs plastidären Loci (atpB-rbcL, matK, psbB-psbH, rpl32-trnL, rps16-trnK, rps16-Intron) und einem nukleären Marker (PHYC) untersucht. Es wurden über 90 Akzessionen von 24 Fosterella-Arten untersucht. Mit 5,6 % informativer Merkmale innerhalb der Gattung war rpl32-trnL der informativste Chloroplastenmarker. Es wurden mit den kombinierten Sequenzdaten eine Maximum Parsimony-, eine Maximum Likelihood- und eine Bayes´sche Analyse berechnet. Weiterhin wurden biogeographische und ultrametrische Untersuchungen durchgeführt. Die 6-Locus-Phylogenie zeigt eine Aufteilung der monophyletischen Gattung Fosterella in sechs Gruppen, von denen vier – die penduliflora-, weddelliana-, weberbaueri- und micrantha-Gruppe - klar monophyletisch und gut gestützt sind. Die albicans- und die rusbyi-Gruppe bilden hingegen einen Komplex. Ultrametrische Analysen legen ein Alter der Gattung von ca. 9,6 Mio. Jahren nahe. Der geographische Ursprung von Fosterella befindet sich nach den vorliegenden biogeographischen Analysen in den Anden und nach der Biom-Analyse zu gleicher Wahrscheinlichkeit entweder in andinen Trockenwäldern (seasonally dry tropical forests, SDTFs) oder in azonalen Standorten des amazonischen Tieflands östlich der Anden. Es gab mehrere Ausbreitungsereignisse, von denen die beiden Fernausbreitungsereignisse nach Mittelamerika (F. micrantha) und in das zentrale Amazonasgebiet (F. batistana) die auffälligsten sind. Die feuchten Bergregenwälder (Yungas) der Anden wurden offenbar mehrfach unabhängig von Fosterella-Arten besiedelt. Insgesamt wurden elf nukleäre Marker (XDH, GS, RPB2, MS, ADH, MS, GLO/PI, CHS, FLO/LFY, NIAi3 und PHYC) auf ihre Anwendbarkeit für molekularsystematische Studien in Fosterella getestet. Davon konnten acht Marker erfolgreich mithilfe einer PCR amplifiziert werden. Die Fragmentgrößen lagen zwischen 350 bp und 1.500 bp. Nur für drei Loci (FLO/LFY, NIAi3 und PHYC) konnten lesbare DNA-Sequenzen in Fosterella erzeugt werden. FLO/LFY zeigte nur 1,5 % Variabilität innerhalb der Gattung. Der NIA-Locus erzeugte bei der Amplifikation mehrere Fragmente, die separat voneinander sequenziert wurden. Der Locus PHYC konnte hingegen aufgrund der guten Amplifizier- und Sequenzierbarkeit für das gesamte Probenset sequenziert werden. Dieser Marker zeigte eine Variabilität innerhalb der Gattung von 10,2 %, davon waren 6,8 % informativ. In der Phylogenie basierend auf PHYC ist Fosterella klar monophyletisch, innerhalb der Gattung zeigt sich jedoch an der Basis eine unaufgelöste Polytomie. Es lassen sich neun mehr oder weniger gut gestützte Artengruppen definieren – rusbyi-, villosula-, albicans-, weddelliana-, penduliflora-, weberbaueri-, micrantha-, robertreadii- und spectabilis-Gruppe - die sich in ihrer Zusammensetzung mit Ausnahme der weddelliana-Gruppe von den nach Chloroplastendaten definierten Gruppen unterscheiden. Viele Arten sind para- oder polyphyletisch, so z. B. F. albicans, F. penduliflora und F. rusbyi. Bei den beiden erstgenannten Arten weisen die unterschiedlichen Stellungen in Chloroplasten- und Kernphylogenie auf Hybridisierungsereignisse hin.

Estudio de las variaciones espaciales y temporales de la diversidad genética de la trucha común, Salmo trutta, en ríos de la Península Ibérica

Se ha analizado las causas de la distribución espacial de la variabilidad genética del ADN mitocondrial en poblaciones de trucha común de la cuenca del Duero y de los Pirineos Orientales. En total se han analizado de novo 49 localidades, 13 en la cuenca del río Duero y 36 en los principales ríos del Pirineo oriental. Además se analizaron las fluctuaciones temporales en 14 de las localidades del Pirineo Oriental. Estudios previos indican un marcado contraste de los patrones de diversidad entre ambos territorios. En la cuenca del río Duero los análisis confirmaron la presencia de los dos linajes matriarcales descritos previamente, el linaje Atlántico (AT) y el linaje Duero (DU). Los análisis de la varianza molecular (AMOVA) siguiendo una jerarquía hidrográfica sugirieron una alta estructuración de las poblaciones coincidente con los patrones ictiológicos observados en la cuenca. El linaje DU parece haber estado presente permanentemente en la cuenca interior del Duero, mientras que las zonas más próximas a la desembocadura han padecido diversas colonizaciones de trucha del linaje AT, que reflejarían los cambios climáticos ocurridos en el Cuaternario. Se ha detectado una discrepancia en el límite entre ambos grupos definidos por genes nucleares (alozimas) y el ADN mitocondrial. Estas discrepancias pueden ser debidas a un efecto más severo de la deriva genética en el ADN mitocondrial que en los marcadores nucleares. Sin embargo, en este trabajo se han observado evidencias a favor de selección en el ADN mitocondrial del linaje DU que también explicaría estas discrepancias. El análisis más exhaustivo en las cuencas de los Pirineos orientales, permitió detectar nuevos haplotipos mitocondriales de los linajes Adriático (AD) y Mediterráneo (ME). En esta región, los AMOVAs confirmaron que las diferencias entre poblaciones dentro de río son más importantes que las diferencias entre ríos. No obstante se observó un patrón de aislamiento por distancia en toda la zona, reflejo de la estructuración de las poblaciones en la cuenca del río Ebro. Además, aunque los AMOVAs mostraron que el componente temporal de la variación es inferior al espacial, las fluctuaciones temporales en la comparación matriarcal de las poblaciones resultaron estadísticamente significativas. Estas fluctuaciones están asociadas tanto a la deriva genética como a procesos de flujo génico entre poblaciones próximas. Dentro de las cuencas, los componentes de diferenciación entre afluentes son, en general, superiores a los obtenidos dentro de cada afluente, patrón que parece estar extendido en la trucha común. Los estudios a escala microgeográfica en la Noguera Vallferrera y Noguera Cardós (afluentes del Noguera Pallaresa) reprodujeron este patrón de diferenciación. Los tamaños efectivos y la tasa de migración entre ambos ríos fueron similares a los descritos en poblaciones noratlánticas. Los tamaños efectivos de las hembras (Nef), calculados a partir del ADN mitocondrial fueron menos de la mitad del tamaño efectivo total tanto en la Noguera Vallferrera como en el resto de localidades pirenaicas estudiadas. Estos bajos tamaños efectivos de las hembras serían también responsables de las fluctuaciones temporales observadas. Los ejemplares repoblados parecen hibridar poco con los nativos, pero su presencia podría intensificar indirectamente los procesos de deriva genética y complicar la conservación de los patrimonios genéticos nativos. Con la salvedad de la existencia de selección que favorece a los haplotipos del linaje DU, los procesos poblacionales que regulan la distribución de la variabilidad genética en la cuenca del Duero y en los Pirineos Orientales podrían ser parecidos y caracterizados por la existencia de múltiples demes interconectados a lo largo del curso fluvial.

Phylogeographic analysis of the chloroplast DNA variation in wild common bean (Phaseolus vulgaris L.) in the Americas

The wild common bean (Phaseolus vulgaris) is widely but discontinuously distributed from northern Mexico to northern Argentina on both sides of the Isthmus of Panama. Little is known on how the species has reached its current disjunct distribution. In this research, chloroplast DNA polymorphisms in seven non-coding regions were used to study the history of migration of wild P. vulgaris between Mesoamerica and South America. A penalized likelihood analysis was applied to previously published Leguminosae ITS data to estimate divergence times between P. vulgaris and its sister taxa from Mesoamerica, and divergence times of populations within P. vulgaris. Fourteen chloroplast haplotypes were identified by PCR-RFLP and their geographical associations were studied by means of a Nested Clade Analysis and Mantel Tests. The results suggest that the haplotypes are not randomly distributed but occupy discrete parts of the geographic range of the species. The current distribution of haplotypes may be explained by isolation by distance and by at least two migration events between Mesoamerica and South America: one from Mesoamerica to South America and another one from northern South America to Mesoamerica. Age estimates place the divergence of P. vulgaris from its sister taxa from Mesoamerica at or before 1.3 Ma, and divergence of populations from Ecuador-northern Peru at or before 0.6 Ma. As these ages are taken as minimum divergence times, the influence of past events, such as the closure of the Isthmus of Panama and the final uplift of the Andes, on the migration history and population structure of this species cannot be disregarded.

The automation and evaluation of nested clade phylogeographic analysis

Nested clade phylogeographic analysis (NCPA) is a popular method for reconstructing the demographic history of spatially distributed populations from genetic data. Although some parts of the analysis are automated, there is no unique and widely followed algorithm for doing this in its entirety, beginning with the data, and ending with the inferences drawn from the data. This article describes a method that automates NCPA, thereby providing a framework for replicating analyses in an objective way. To do so, a number of decisions need to be made so that the automated implementation is representative of previous analyses. We review how the NCPA procedure has evolved since its inception and conclude that there is scope for some variability in the manual application of NCPA. We apply the automated software to three published datasets previously analyzed manually and replicate many details of the manual analyses, suggesting that the current algorithm is representative of how a typical user will perform NCPA. We simulate a large number of replicate datasets for geographically distributed, but entirely random-mating, populations. These are then analyzed using the automated NCPA algorithm. Results indicate that NCPA tends to give a high frequency of false positives. In our simulations we observe that 14% of the clades give a conclusive inference that a demographic event has occurred, and that 75% of the datasets have at least one clade that gives such an inference. This is mainly due to the generation of multiple statistics per clade, of which only one is required to be significant to apply the inference key. We survey the inferences that have been made in recent publications and show that the most commonly inferred processes (restricted gene flow with isolation by distance and contiguous range expansion) are those that are commonly inferred in our simulations. However, published datasets typically yield a richer set of inferences with NCPA than obtained in our random-mating simulations, and further testing of NCPA with models of structured populations is necessary to examine its accuracy.

Population subdivision and molecular sequence variation: Theory and analysis of Drosophila ananassae data

Population subdivision complicates analysis of molecular variation. Even if neutrality is assumed, three evolutionary forces need to be considered: migration, mutation, and drift. Simplification can be achieved by assuming that the process of migration among and drift within subpopulations is occurring fast compared to Mutation and drift in the entire population. This allows a two-step approach in the analysis: (i) analysis of population subdivision and (ii) analysis of molecular variation in the migrant pool. We model population subdivision using an infinite island model, where we allow the migration/drift parameter Theta to vary among populations. Thus, central and peripheral populations can be differentiated. For inference of Theta, we use a coalescence approach, implemented via a Markov chain Monte Carlo (MCMC) integration method that allows estimation of allele frequencies in the migrant pool. The second step of this approach (analysis of molecular variation in the migrant pool) uses the estimated allele frequencies in the migrant pool for the study of molecular variation. We apply this method to a Drosophila ananassae sequence data set. We find little indication of isolation by distance, but large differences in the migration parameter among populations. The population as a whole seems to be expanding. A population from Bogor (Java, Indonesia) shows the highest variation and seems closest to the species center.