Forward-in-Time, Spatially Explicit Modeling Software to Simulate Genetic Lineages Under Selection

SELECTOR is a software package for studying the evolution of multiallelic genes under balancing or positive selection while simulating complex evolutionary scenarios that integrate demographic growth and migration in a spatially explicit population framework. Parameters can be varied both in space and time to account for geographical, environmental, and cultural heterogeneity. SELECTOR can be used within an approximate Bayesian computation estimation framework. We first describe the principles of SELECTOR and validate the algorithms by comparing its outputs for simple models with theoretical expectations. Then, we show how it can be used to investigate genetic differentiation of loci under balancing selection in interconnected demes with spatially heterogeneous gene flow. We identify situations in which balancing selection reduces genetic differentiation between population groups compared with neutrality and explain conflicting outcomes observed for human leukocyte antigen loci. These results and three previously published applications demonstrate that SELECTOR is efficient and robust for building insight into human settlement history and evolution.

Comparison of different methods to detect genetic barriers in a small mammal population

Habitat fragmentation and the consequently the loss of connectivity between populations can reduce the individuals interchange and gene flow, increasing the chances of inbreeding, and the increase the risk of local extinction. Landscape genetics is providing more and better tools to identify genetic barriers.. To our knowledge, no comparison of methods in terms of consistency has been made with observed data and species with low dispersal ability. The aim of this study is to examine the consistency of the results of five methods to detect barriers to gene flow in a Mediterranean pine vole population Microtus duodecimcostatus: F-statistics estimations, Non-Bayesian clustering, Bayesian clustering, Boundary detection and Simple/Partial Mantel tests. All methods were consistent in detecting the stream as a non-genetic barrier. However, no consistency in results among the methods were found regarding the role of the highway as a genetic barrier. Fst, Bayesian clustering assignment test and Partial Mantel test identifyed the highway as a filter to individual interchange. The Mantel tests were the most sensitive method. Boundary detection method (Monmonier’s Algorithm) and Non-Bayesian approaches did not detect any genetic differentiation of the pine vole due to the highway. Based on our findings we recommend that the genetic barrier detection in low dispersal ability populations should be analyzed with multiple methods such as Mantel tests, Bayesian clustering approaches because they show more sensibility in those scenarios and with boundary detection methods by having the aim of detect drastic changes in a variable of interest between the closest individuals. Although simulation studies highlight the weaknesses and the strengths of each method and the factors that promote some results, tests with real data are needed to increase the effectiveness of genetic barrier detection.

Genetic diversity and connectivity remain high in Holothuria polii (Delle Chiaje 1823) across a coastal lagoon-open sea environmental gradient

Coastal lagoons represent habitats with widely heterogeneous environmental conditions, particularly as regards salinity and temperature,which fluctuate in both space and time. These characteristics suggest that physical and ecological factors could contribute to the genetic divergence among populations occurring in coastal lagoon and opencoast environments. This study investigates the genetic structure of Holothuria polii at a micro-geographic scale across theMar Menor coastal lagoon and nearbymarine areas, estimating the mitochondrial DNA variation in two gene fragments, cytochrome oxidase I (COI) and 16S rRNA (16S). Dataset of mitochondrial sequences was also used to test the influence of environmental differences between coastal lagoon andmarine waters on population genetic structure. All sampled locations exhibited high levels of haplotype diversity and low values of nucleotide diversity. Both genes showed contrasting signals of genetic differentiation (non-significant differences using COI and slight differences using 16S, which could due to different mutation rates or to differential number of exclusive haplotypes. We detected an excess of recent mutations and exclusive haplotypes, which can be generated as a result of population growth. However, selective processes can be also acting on the gene markers used; highly significant generalized additive models have been obtained considering genetic data from16S gene and independent variables such as temperature and salinity.

Spatial patterns of genetic diversity in Mediterranean eagle owl Bubo bubo populations

Little information is available on the patterns of genetic connectivity in owls. We studied the genetic structure of the eagle owl Bubo bubo (Linnaeus, 1758) in southeastern Spain at two different spatial scales. Seven microsatellites previously described for this species were used, although only six loci amplified correctly. The observed low genetic variation could be explained by the short dispersal distance, high mortality rate and high degree of monogamy shown by this large nocturnal predator. As expected, the highest genetic isolation was detected in the geographically most isolated population. Significant genetic differentiation was found among study units separated by less than 50 km. The territorial analysis showed interesting connectivity patterns related with the gene flow and turnover rate of the breeding individuals. The lowest genetic diversity was found in the region with the largest population, which could imply incipient inbreeding.

Small-scale genetic structure of Cerastoderma glaucum in a lagoonal environment: potential significance of habitat discontinuity and unstable population dynamics

Environmental heterogeneity in coastal lagoons is expected to facilitate local adaptation in response to different ecological conditions, causing significant genetic structuring within lagoon populations at a small scale and also differentiation between lagoons. However, these patterns and processes of genetic structuring are still poorly understood. The aims of our study were (1) to seek genetic structure at a small scale in Cerastoderma glaucum inside the Mar Menor coastal lagoon using a mitochondrial DNA marker (COI) that has previously detected genetic differentiation inside the lagoon in other species and (2) to evaluate the influence of extreme environmental conditions and habitat discontinuity on its genetic composition. The results indicate high levels of haplotype diversity and low values of nucleotide diversity. COI data provide evidence of significant population differentiation among some localities within the lagoon. Limited gene flow and unstable population dynamics (i.e. fluctuations in population size caused by local extinction and recolonization), probably due to the high environmental heterogeneity, could generate the small-scale genetic divergence detected between populations within the lagoon.

Genetic diversity across geographical scales in marine coastal ecosystems: Holothuria arguinensis a model species

Coastal lagoons are considered one of the most productive areas of our planet harboring a large variety of habitats. Their transitional character, between terrestrial and marine environments, creates a very particular ecosystem with important variations of its environmental conditions. The organisms that are able to survive on these ecosystems frequently experience strong selective pressures and constrictions to gene flowwith marine populations, which could contribute to genetic divergence among populations inhabiting coastal lagoon and marine environments. Therefore, the main aims of this study are to asses the genetic diversity and population structure of Holothuria arguinensis across geographical ranges, to test the hypothesis of coastal lagoons as hotspots of genetic diversity in the Ria Formosa lagoon, and to determine the role of exporting standing genetic variation from the lagoon to open sea and their implications to recent geographical expansion events. To reach these objectives, we investigate the genetic structure of H. arguinensis using two mitochondrial DNA markers (COI and 16S) at different spatial scales: i) small, inside Ria Formosa coastal lagoon, South Portugal; 2) large, including most of the geographical distribution of this species (South and Western Portuguese coast and Canary islands); these results will allow us to compare the genetic diversity of lagoonal and marine populations of H. arguinensis. On this framework, its recent geographical expansion events, recorded by Rodrigues (2012) and González-Wangüemert and Borrero-Pérez (2012), will be analyzed considering the potential contribution from lagoonal genetic pool. Non-significant genetic structure and high haplotypic diversity were found inside the Ria Formosa coastal lagoon. Both genes were unable to detect significant genetic differentiation among lagoonal and marine localities, suggesting a high rate of gene flow. The results supported our hypotheses that coastal lagoons are not only acting as hotspots of genetic diversity, but also contributing for the genetic variability of the species, working as a source of new haplotypes and enhancing adaptation to the high variable conditions. Different genetic patterns of colonization were found on H. arguinensis, but they must be studied more deeply.

Genetic considerations on the introduction of farmed fish in marine protected areas: The case of study of white seabream restocking in the Gulf of Castellammare (Southern Tyrrhenian Sea)

Human exploitation has drastically reduced the abundance and distribution of several marine fish and invertebrate populations through overfishing and habitat destruction. Restocking can potentially mitigate these impacts and help to reconstitute depleted stocks but genetic repercussions must be considered. In the present study, the degree of genetic similarity between white seabream (Diplodus sargus Linnaeus 1758) individuals reared for restocking purposes and the receiving population in the Gulf of Castellammare fishery reserve (Sicily, Italy) was assessed using microsatellites. We also inferred the spatial pattern of the genetic structure of D. sargus and connectivity along Sicilian coasts. The farmed population showed significant heterozygosity deficiency in 6 loci and an important reduction in the number of alleles, which could indicate an incipient inbreeding. Both the farmed population and the target one for restocking (Castellammare fishery reserve), showed high and significant values of genetic differentiation due to different allele frequencies, number of privative alleles and total number of alleles. These findings indicate a low degree of genetic similarity between both populations, therefore this restocking initiative is not advisable. The genetic connectivity pattern, highly consistent with oceanographic currents, identified two distinct metapopulations of white seabream around Sicily. Thus it is recommended to utilize broods from the same metapopulation for restocking purposes to provide a better genetic match to the wild populations.

High interannual variability in connectivity and genetic pool of a temperate clingfish matches oceanographic transport predictions

Adults of most marine benthic and demersal fish are site-attached, with the dispersal of their larval stages ensuring connectivity among populations. In this study we aimed to infer spatial and temporal variation in population connectivity and dispersal of a marine fish species, using genetic tools and comparing these with oceanographic transport. We focused on an intertidal rocky reef fish species, the shore clingfish Lepadogaster lepadogaster, along the southwest Iberian Peninsula, in 2011 and 2012. We predicted high levels of self-recruitment and distinct populations, due to short pelagic larval duration and because all its developmental stages have previously been found near adult habitats. Genetic analyses based on microsatellites countered our prediction and a biophysical dispersal model showed that oceanographic transport was a good explanation for the patterns observed. Adult sub-populations separated by up to 300 km of coastline displayed no genetic differentiation, revealing a single connected population with larvae potentially dispersing long distances over hundreds of km. Despite this, parentage analysis performed on recruits from one focal site within the Marine Park of Arrábida (Portugal), revealed self-recruitment levels of 2.5% and 7.7% in 2011 and 2012, respectively, suggesting that both long- and short-distance dispersal play an important role in the replenishment of these populations. Population differentiation and patterns of dispersal, which were highly variable between years, could be linked to the variability inherent in local oceanographic processes. Overall, our measures of connectivity based on genetic and oceanographic data highlight the relevance of long-distance dispersal in determining the degree of connectivity, even in species with short pelagic larval durations.

Genetic signature of the northward expansion of the Egyptian mongoose Herpestes ichneumon (Herpestidae) in the Iberian Peninsula

In the last three decades, the range of the Egyptian mongoose (Herpestes ichneumon) has increased in the Iberian Peninsula. A panel of microsatellites was used to confront the patterns of genetic diversity of the species with the scenario of its recent northward expansion in its Iberian range. Evidence of substructure and significant genetic differentiation within the studied population were recorded, with a central-northern subpopulation (CNorth) and a southern subpopulation (S). Northward range expansion was supported by the observed allelic frequencies, diversity parameters, and observed heterozygosity of the studied loci, with S showing a higher allelic diversity and a higher number of private alleles than CNorth. Patterns of isolation-by-distance and isolation-by-barrier as a result of the Tagus River were demonstrated, suggesting that the river acted as a semi-permeable barrier, possibly leading to genetic differentiation of the studied population. The observed individuals from CNorth in southern locations and individuals from S in central/northern areas might comprise evidence for long range dispersals across the studied range. A bottleneck event after population expansion was supported by a significant heterozygosity deficiency in CNorth, which is in agreement with a scenario of founder events occurring in recently colonized areas after the crossing of the Tagus River.

Genetic structure and mating type analysis of the Pyricularia oryzae population causing widespread epidemics in southern Brazil.

Blast is a major disease of rice in Brazil, the largest rice-producing country outside Asia. This study aimed to assess the genetic structure and mating-type frequency in a contemporary Pyricularia oryzae population, which caused widespread epidemics during the 2012/13 season in the Brazilian lowland subtropical region. Symptomatic leaves and panicles were sampled at flooded rice fields in the states of Rio Grande do Sul (RS, 34 fields) and Santa Catarina (SC, 21 fields). The polymorphism at ten simple sequence repeats (SSR or microsatellite) loci and the presence of MAT1-1 or MAT1-2 idiomorphs were assessed in a population comprised of 187 isolates. Only the MAT1-2 idiomorph was found and 162 genotypes were identified by the SSR analysis. A discriminant analysis of principal components (DAPC) of SSR data resolved four genetic groups, which were strongly associated with the cultivar of origin of the isolates. There was high level of genotypic diversity and moderate level of gene diversity regardless whether isolates were grouped in subpopulations based on geographic region, cultivar host or cultivar within region. While regional subpopulations were weakly differentiated, high genetic differentiation was found among subpopulations comprised of isolates from different cultivars. The data suggest that the rice blast pathogen population in southern Brazil is comprised of clonal lineages that are adapting to specific cultivar hosts. Farmers should avoid the use of susceptible cultivars over large areas and breeders should focus at enlarging the genetic basis of new cultivars.

Genomic evidence for the parallel evolution of coastal forms in the Senecio lautus complex

Instances of parallel ecotypic divergence where adaptation to similar conditions repeatedly cause similar phenotypic changes in closely related organisms are useful for studying the role of ecological selection in speciation. Here we used a combination of traditional and next generation genotyping techniques to test for the parallel divergence of plants from the Senecio lautus complex, a phenotypically variable groundsel that has adapted to disparate environments in the South Pacific. Phylogenetic analysis of a broad selection of Senecio species showed that members of the S. lautus complex form a distinct lineage that has diversified recently in Australasia. An inspection of thousands of polymorphisms in the genome of 27 natural populations from the S. lautus complex in Australia revealed a signal of strong genetic structure independent of habitat and phenotype. Additionally, genetic differentiation between populations was correlated with the geographical distance separating them, and the genetic diversity of populations strongly depended on geographical location. Importantly, coastal forms appeared in several independent phylogenetic clades, a pattern that is consistent with the parallel evolution of these forms. Analyses of the patterns of genomic differentiation between populations further revealed that adjacent populations displayed greater genomic heterogeneity than allopatric populations and are differentiated according to variation in soil composition. These results are consistent with a process of parallel ecotypic divergence in face of gene flow.

High connectivity in rastrelliger kanagurta: Influence of historical signatures and migratory behaviour inferred from mtDNA Cytochrome b

Phylogeographic patterns and population structure of the pelagic Indian mackerel, Rastrelliger kanagurta were examined in 23 populations collected from the Indonesian-Malaysian Archipelago (IMA) and the West Indian Ocean (WIO). Despite the vast expanse of the IMA and neighbouring seas, no evidence for geographical structure was evident. An indication that R. kanagurta populations across this region are essentially panmictic. This study also revealed that historical isolation was insufficient for R. kanagurta to attain migration drift equilibrium. Two distinct subpopulations were detected between the WIO and the IMA (and adjacent populations); interpopulation genetic variation was high. A plausible explanation for the genetic differentiation observed between the IMA and WIO regions suggest historical isolation as a result of fluctuations in sea levels during the late Pleistocene. This occurrence resulted in the evolution of a phylogeographic break for this species to the north of the Andaman Sea.

The spatiotemporal dynamics of Tribolium castaneum (Herbst): adult flight and gene flow.

Tribolium castaneum (Herbst) has been used as a model organism to develop and test important ecological and evolutionary concepts and is also a major pest of grain and grain products globally. This beetle species is assumed to be a good colonizer of grain storages through anthropogenic movement of grain, and active dispersal by flight is considered unlikely. Studies using T. castaneum have therefore used confined or walking insects. We combine an ecological study of dispersal with an analysis of gene flow using microsatellites to investigate the spatiotemporal dynamics and adult flight of T. castaneum in an ecological landscape in eastern Australia. Flying beetles were caught in traps at grain storages and in fields at least 1 km from the nearest stored grain at regular intervals for an entire year. Significantly more beetles were trapped at storages than in fields, and almost no beetles were caught in native vegetation reserves many kilometres from the nearest stored grain. Genetic differentiation between beetles caught at storages and in fields was low, indicating that although T. castaneum is predominantly aggregated around grain storages, active dispersal takes place to the extent that significant gene flow occurs between them, mitigating founder effects and genetic drift. By combining ecological and molecular techniques, we reveal much higher levels of active dispersal through adult flight in T. castaneum than previously thought. We conclude that the implications of adult flight to previous and future studies on this model organism warrant consideration.

Population structure of Pyrenophora teres, the causal agent of net blotch of barley

Spring barley is the most important crop in Finland based on cultivated land area. Net blotch, a disease caused by Pyrenophora teres Drech., is the most damaging disease of barley in Finland. The pressure to improve the economics and efficiency of agriculture has increased the need for more efficient plant protection methods. Development of durable host-plant resistance to net blotch is a promising possibility. However, deployment of disease resistant crops could initiate selection pressure on the pathogen (P. teres) population. The aim of this study was to understand the population biology of P. teres and to estimate the evolutionary potential of P. teres under selective pressure following deployment of resistance genes and application of fungicides. The study included mainly Finnish P. teres isolates. Population samples from Russia and Australia were also included. Using AFLP markers substantial genotypic variation in P. teres populations was identified. Differences among isolates were least within Finnish fields and significantly higher in Krasnodar, Russia. Genetic differentiation was identified among populations from northern Europe and from Australia, and between the two forms P. teres f. teres (PTT, net form of net blotch) and P. teres f. maculata (PTM, spot form of net blotch) in Australia. Differentiation among populations was also identified based on virulence between Finnish and Russian populations, and based on prochloraz (fungicide) tolerance in the Häme region in Finland. Surprisingly only PTT was recovered from Finland and Russia although both forms were earlier equally common in Finland. The reason for the shift in occurrence of forms in Finland remained uncertain. Both forms were found within several fields in Australia. Sexual reproduction of P. teres was supported by recover of both mating types in equal ratio in those areas although the prevalence of sexual mating seems to be less in Finland than in Australia. Population from Krasnodar was an exception since only one mating type was found in there. Based on the substantial high genotypic variation in Krasnodar it was suggested go represent an old P. teres population, whereas the Australian samples were suggested to represent newer populations. In conclusion, P. teres populations are differentiated at several levels. Human assistance in dispersal of P. teres on infected barley seed is obvious and decreases the differentiation among populations. This can increase the plant protection problems caused by this pathogen. P. teres is capable of sexual reproduction in several areas but the prevalence varies. Based on these findings it is apparent that P. teres has the potential to pose more serious problems in barley cultivation if plant protection is neglected. Therefore, good agricultural practices, including crop rotation and the use of healthy seed, are recommended.

Ecological impacts of Phlebiopsis gigantea biocontrol treatment against Heterobasidion spp. as revealed by fungal community profiling and population analyses

Wood decay fungi belonging to the species complex Heterobasidion annosum sensu lato are among the most common and economically important species causing root rot and stem decay in conifers of the northern temperate regions. New infections by these pathogens can be suppressed by tree stump treatments using chemical or biological control agents. In Finland, the corticiaceous fungus Phlebiopsis gigantea has been formulated into a commercial biocontrol agent called Rotstop (Verdera Ltd.). This thesis addresses the ecological impacts of Rotstop biocontrol treatment on the mycoflora of conifer stumps. Locally, fungal communities within Rotstop-treated and untreated stumps were analyzed using a novel method based on DGGE profiling of small subunit ribosomal DNA fragments amplified directly from wood samples. Population analyses for P. gigantea and H. annosum s.l. were conducted to evaluate possible risks associated with local and/or global distribution of the Rotstop strain. Based on molecular community profiling by DGGE, we detected a few individual wood-inhabiting fungal species (OTUs) that seemed to have suffered or benefited from the Rotstop biocontrol treatment. The DGGE analyses also revealed fungal diversity not retrieved by cultivation and some fungal sequence types untypical for decomposing conifer wood. However, statistical analysis of DGGE community profiles obtained from Rotstop-treated and untreated conifer stumps revealed that the Rotstop treatment had not caused a statistically significant reduction in the species diversity of wood-inhabiting fungi within our experimental forest plots. Locally, ISSR genotyping of cultured P. gigantea strains showed that the Rotstop biocontrol strain was capable of surviving up to six years within treated Norway spruce stumps, while in Scots pine stumps it was sooner replaced by successor fungal species. In addition, the spread of resident P. gigantea strains into Rotstop-treated forest stands seemed effective in preventing the formation of genetically monomorphic populations in the short run. On a global scale, we detected a considerable level of genetic differentiation between the interfertile European and North American populations of P. gigantea. These results strongly suggest that local biocontrol strains should be used in order to prevent global spread of P. gigantea and hybrid formation between geographically isolated populations. The population analysis for H. annosum s.l. revealed a collection of Chinese fungal strains that showed a high degree of laboratory fertility with three different allopatric H. annosum s.l. taxa. However, based on the molecular markers, the Chinese strains could be clearly affiliated with the H. parviporum taxonomical cluster, which thus appears to have a continuous distribution range from Europe through southern Siberia to northern China. Keywords: Rotstop, wood decay, DGGE, ISSR fingerprinting, ribosomal DNA