Rapid sympatric ecological differentiation of crater lake cichlid fishes within historic times

Background: After a volcano erupts, a lake may form in the cooled crater and become an isolated aquatic ecosystem. This makes fishes in crater lakes informative for understanding sympatric evolution and ecological diversification in barren environments. From a geological and limnological perspective, such research offers insight about the process of crater lake ecosystem establishment and speciation. In the present study we use genetic and coalescence approaches to infer the colonization history of Midas cichlid fishes (Amphilophus cf. citrinellus) that inhabit a very young crater lake in Nicaragua-the ca. 1800 year-old Lake Apoyeque. This lake holds two sympatric, endemic morphs of Midas cichlid: one with large, hypertrophied lips (~20% of the total population) and another with thin lips. Here we test the associated ecological, morphological and genetic diversification of these two morphs and their potential to represent incipient speciation.
Results: Gene coalescence analyses [11 microsatellite loci and mitochondrial DNA (mtDNA) sequences] suggest that crater lake Apoyeque was colonized in a single event from the large neighbouring great lake Managua only about 100 years ago. This founding in historic times is also reflected in the extremely low nuclear and mitochondrial genetic diversity in Apoyeque. We found that sympatric adult thin- and thick-lipped fishes occupy distinct ecological trophic niches. Diet, body shape, head width, pharyngeal jaw size and shape and stable isotope values all differ significantly between the two lip-morphs. The eco-morphological features pharyngeal jaw shape, body shape, stomach contents and stable isotopes (d15N) all show a bimodal distribution of traits, which is compatible with the expectations of an initial stage of ecological speciation under disruptive selection. Genetic differentiation between the thin- and thick-lipped population is weak at mtDNA sequence (FST = 0.018) and absent at nuclear microsatellite loci (FST < 0.001).
Conclusions: This study provides empirical evidence of eco-morphological differentiation occurring very quickly after the colonization of a new and vacant habitat. Exceptionally low levels of neutral genetic diversity and inference from coalescence indicates that the Midas cichlid population in Apoyeque is much younger (ca. 100 years or generations old) than the crater itself (ca. 1 800 years old). This suggests either that the crater remained empty for many hundreds of years after its formation or that remnant volcanic activity prevented the establishment of a stable fish population during the early life of the crater lake. Based on our findings of eco-morphological variation in the Apoyeque Midas cichlids, and known patterns of adaptation in Midas cichlids in general, we suggest that this population may be in a very early stage of speciation (incipient species), promoted by disruptive selection and ecological diversification.

DiveRsity: An R package for the estimation and exploration of population genetics parameters and their associated errors

Summary: We present a new R package, diveRsity, for the calculation of various diversity statistics, including common diversity partitioning statistics (?, G) and population differentiation statistics (D, GST ', ? test for population heterogeneity), among others. The package calculates these estimators along with their respective bootstrapped confidence intervals for loci, sample population pairwise and global levels. Various plotting tools are also provided for a visual evaluation of estimated values, allowing users to critically assess the validity and significance of statistical tests from a biological perspective. diveRsity has a set of unique features, which facilitate the use of an informed framework for assessing the validity of the use of traditional F-statistics for the inference of demography, with reference to specific marker types, particularly focusing on highly polymorphic microsatellite loci. However, the package can be readily used for other co-dominant marker types (e.g. allozymes, SNPs). Detailed examples of usage and descriptions of package capabilities are provided. The examples demonstrate useful strategies for the exploration of data and interpretation of results generated by diveRsity. Additional online resources for the package are also described, including a GUI web app version intended for those with more limited experience using R for statistical analysis. © 2013 British Ecological Society.

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

The evolution of jumping performance in anurans: morphological correlates and ecological implications

We investigated the evolution of anuran locomotor performance and its morphological correlates as a function of habitat use and lifestyles. We reanalysed a subset of the data reported by Zug (Smithson. Contrib. Zool. 1978; 276: 1-31) employing phylogenetically explicit statistical methods (n = 56 species), and assembled morphological data on the ratio between hind-limb length and snout-vent length (SVL) from the literature and museum specimens for a large subgroup of the species from the original paper (n = 43 species). Analyses using independent contrasts revealed that classifying anurans into terrestrial, semi-aquatic, and arboreal categories cannot distinguish between the effects of phylogeny and ecological diversification in anuran locomotor performance. However, a more refined classification subdividing terrestrial species into `fossorials` and `non-fossorials`, and arboreal species into `open canopy`, `low canopy` and `high canopy`, suggests that part of the variation in locomotor performance and in hind-limb morphology can be attributed to ecological diversification. In particular, fossorial species had significantly lower jumping performances and shorter hind limbs than other species after controlling for SVL, illustrating how the trade-off between burrowing efficiency and jumping performance has resulted in morphological specialization in this group.

Spatial patterns, landscape genetics and post virus recovery of blacklip abalone, Haliotis rubra (Leach), in the western commercial fishing zone of Victoria

Historically, collecting nearshore habitat information has been problematic. Existing methods, such as aerial and satellite image interpretation are limited due to the attenuation of light in the water column obscuring the seabed structure. The advent of airborne bathymetric LiDAR (Light Detection and Ranging) systems (laser scanning of the seabed) now provides high-resolution seabed ‘images’ in areas that were previously difficult to survey. LiDAR imagery is available for the entire coastline of Victoria, Australia to depths of around 25 m, after being initially collected for climate change modelling by the Future Coasts Program (http://www.climatechange.vic.gov.au/adapting-to-climate-change/future-coasts). This dataset has provided the opportunity to test its applicability to inform fisheries management. Detailed geophysical information combined with spatially explicit AbTrack GPS located fisheries records and targeted genetic sampling is used in this study to provide a better understanding of the extent of available fishing grounds, direction of fishing effort and stock population structure within the Victorian western zone abalone fishery.
The species distribution modelling technique MaxEnt was used to produce a potential habitat suitability map for abalone in an attempt to capture the effective footprint of the  fishery. Also, by interrogating the spatially defined effort localities, we demonstrate an approach that may be used to identify areas where fishing effort is concentrated, and how this parameter changes temporally.
Despite barriers to adult dispersal (soft sediment barriers between reef patches), the genetic study indicates that larval movement is able to homogenize the gene pool over  large geographic distances. The western, central and eastern zone abalone stocks in Victoria were found to be a single large panmictic unit. This indicates high levels of stock connectivity and no obvious impacts of Abalone Viral Ganglioneuritis (AVG) on the genetic health of western zone stocks. We used detailed seafloor structure information interpreted from LiDAR to inform a replicated hierarchical fine scale genetic sampling design. We demonstrated that there may be extensive migration among abalone stocks across the Victorian abalone fishery.
This is contrary to previous studies that suggest recruitment is highly localised. In combination, these findings provide a valuable insight into the biology of H. rubra and immediate benefits for fisheries management. We discuss these results in the context of predicting resilience and adaptive potential of H. rubra stocks to environmental pressures and the spread of heritable diseases.
Adoption pathways are also provided to benefit future stock augmentation activities to catalyse the recovery of AVG affected reef codes. As larval dispersal is likely to be spatially and temporally variable, some AVG affected stocks are likely to recover through natural recruitment, while others will benefit from augmentation activities to ‘kick-start’ stock recovery. Evidence of neutral genetic homogeneity across Victorian reef codes suggests that the relocation of animals is unlikely to have significant genetic risks; however the potential for locally adaptive genetic differences may exist, and should be taken into consideration in future stock augmentation planning.
When combined, the spatial and genetic analyses provide valuable insights into stock productivity within the western zone fishery. Reefs appear to be expansive and support much available habitat, and the movement of larvae among reef structures is likely to be extensive in this region. Consequently, we propose that colonisation success and productivity is likely to be driven by ecological factors such as resources and/or competition, or physical factors such as wave exposure.

Non invasive methods for genetic analysis applied to ecological and behavioral studies in Latino-America

Documenting the presence and abundance of the neotropical mammals is the first step for understanding their population ecology, behavior and genetic dynamics in designing conservation plans. The combination of field research with molecular genetics techniques are new tools that provide valuable biological information avoiding the disturbance in the ecosystems, trying to minimize the human impact in the process to gather biological information. The objective of this paper is to review the available non invasive sampling techniques that have been used in Neotropical mammal studies to apply to determine the presence and abundance, population structure, sex ratio, taxonomic diagnostic using mitochondrial markers, and assessing genetic variability using nuclear markers. There are a wide range of non invasive sampling techniques used to determine the species identification that inhabit an area such as searching for tracks, feces, and carcasses. Other useful equipment is the camera traps that can generate an image bank that can be valuable to assess species presence and abundance by morphology. With recent advances in molecular biology, it is now possible to use the trace amounts of DNA in feces and amplify it to analyze the species diversity in an area, and the genetic variability at intraspecific level. This is particularly helpful in cases of sympatric and cryptic species in which morphology failed to diagnose the taxonomic status of several species of brocket deer of the genus Mazama.

The influence of landscape characteristics and home-range size on the quantification of landscape-genetics relationships

A common approach used to estimate landscape resistance involves comparing correlations of ecological and genetic distances calculated among individuals of a species. However, the location of sampled individuals may contain some degree of spatial uncertainty due to the natural variation of animals moving through their home range or measurement error in plant or animal locations. In this study, we evaluate the ways that spatial uncertainty, landscape characteristics, and genetic stochasticity interact to influence the strength and variability of conclusions about landscape-genetics relationships. We used a neutral landscape model to generate 45 landscapes composed of habitat and non-habitat, varying in percent habitat, aggregation, and structural connectivity (patch cohesion). We created true and alternate locations for 500 individuals, calculated ecological distances (least-cost paths), and simulated genetic distances among individuals. We compared correlations between ecological distances for true and alternate locations. We then simulated genotypes at 15 neutral loci and investigated whether the same influences could be detected in simple Mantel tests and while controlling for the effects of isolation-by distance using the partial Mantel test. Spatial uncertainty interacted with the percentage of habitat in the landscape, but led to only small reductions in correlations. Furthermore, the strongest correlations occurred with low percent habitat, high aggregation, and low to intermediate levels of cohesion. Overall genetic stochasticity was relatively low and was influenced by landscape characteristics.

Simulating Population Genetics of Pathogen Vectors in Changing Landscapes: Guidelines and Application with Triatoma brasiliensis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Conservation genetics of European wildcat (Felis silvestris silvestris): a wide and integrating analysis protocol for admixture inferences and population structure

Introgression of domestic cat genes into European wildcat (Felis silvestris silvestris) populations and reduction of wildcats’ range in Europe, leaded by habitat loss and fragmentation, are considered two of the main conservation problems for this endangered feline. This thesis addressed the questions related with the artificial hybridization and populations’ fragmentation, using a conservation genetics perspective. We combined the use of highly polymorphic loci, Bayesian statistical inferences and landscape analyses tools to investigate the origin of the geographic-genetic substructure of European wildcats (Felis silvestris silvestris) in Italy and Europe. The genetic variability of microsatellites evidenced that European wildcat populations currently distributed in Italy differentiated in, and expanded from two distinct glacial refuges during the Last Glacial Maximum. The genetic and geographic substructure detected between the eastern and western sides of the Apennine ridge, resulted by adaptation to specific ecological conditions of the Mediterranean habitats. European wildcat populations in Europe are strongly structured into 5 geographic-genetic macro clusters corresponding to: the Italian peninsular & Sicily; Balkans & north-eastern Italy; Germany eastern; central Europe; and Iberian Peninsula. Central European population might have differentiated in the extra-Mediterranean Würm ice age refuge areas (Northern Alps, Carpathians, and the Bulgarian mountain systems), while the divergence among and within the southern European populations might have resulted by the Pleistocene bio geographical framework of Europe, with three southern refugia localized in the Balkans, Italian Peninsula and Iberia Peninsula. We further combined the use of most informative autosomal SNPs with uniparental markers (mtDNA and Y-linked) for accurately detecting parental genotypes and levels of introgressive hybridization between European wild and domestic cats. A total of 11 hybrids were identified. The presence of domestic mitochondrial haplotypes shared with some wild individuals led us to hypnotize the possibility that ancient introgressive events might have occurred and that further investigation should be recommended.

Is hatchery stocking a help or harm? Evidence, limitations and future directions in ecological and genetic surveys

Hatchery fish stocking for stock enhancement has been operated at a massive and global scale. However, the use of hatchery fish as a means of stock enhancement is highly controversial, and little is known about its effects on wild stock and consequences for stock enhancement. Here we review the scientific literature on this subject in order to address a fundamental - question is hatchery stocking a help or harm for wild stock and stock enhancement? We summarized 266 peer-reviewed papers that were published in the last 50 years, which describe empirical case studies on ecology and genetics of hatchery stocks and their effects on stock enhancement. Specifically, we asked whether hatchery stock and wild stock differed in fitness and the level of genetic variation, and whether stocking affected population abundance. Seventy studies contained comparisons between hatchery and wild stocks, out of which 23 studies showed significantly negative effects of hatchery rearing on the fitness of stocked fish, and 28 studies showed reduced genetic variation in hatchery populations. None of these studies suggested a positive genetic effect on the fitness of hatchery-reared individuals after release. These results suggest that negative effects of hatchery rearing are not just a concern but undeniably present in many aquaculture species. In a few cases, however, no obvious effect of hatchery rearing was observed, and a positive contribution of hatchery stock to the abundance of fish populations was indicated. These examples suggest that there is a chance to improve hatchery practices and mitigate the negative effects on wild stocks, although scientific data supporting the positive effect on stock enhancement are largely missing at this moment. Technically, microsatellite-based parentage assignments have been proven as a useful tool for the evaluation of reproductive fitness in natural settings, which is a key for stock enhancement by hatchery-based stocking. We discuss implications of these results, as well as their limitations and future directions. (C) 2010 Elsevier B.V. All rights reserved.

Divergence along a steep ecological gradient in lake whitefish (Coregonus sp.)

To understand mechanisms structuring diversity in young adaptive radiations, quantitative and unbiased information about genetic and phenotypic diversity is much needed. Here, we present the first in-depth investigation of whitefish diversity in a Swiss lake, with continuous spawning habitat sampling in both time and space. Our results show a clear cline like pattern in genetics and morphology of populations sampled along an ecological depth gradient in Lake Neuchâtel. Divergent natural selection appears to be involved in shaping this cline given that trait specific P(ST)-values are significantly higher than F(ST)-values when comparing populations caught at different depths. These differences also tend to increase with increasing differences in depth, indicating adaptive divergence along a depth gradient, which persists despite considerable gene flow between adjacent demes. It however remains unclear, whether the observed pattern is a result of currently stable selection-gene flow balance, incipient speciation, or reverse speciation due to anthropogenic habitat alteration causing two formerly divergent species to collapse into a single gene pool.