High degree of cryptic population differentiation in the Baltic Sea herring Clupea harengus

Population genomics reveals seahorses (Hippocampus erectus) of the western mid-Atlantic coast to be residents rather than vagrants

Understanding population structure and areas of demographic persistence and transients is critical for effective species management. However, direct observational evidence to address the geographic scale and delineation of ephemeral or persistent populations for many marine fishes is limited. The Lined seahorse (Hippocampus erectus) can be commonly found in three western Atlantic zoogeographic provinces, though inhabitants of the temperate northern Virginia Province are often considered tropical vagrants that only arrive during warm seasons from the southern provinces and perish as temperatures decline. Although genetics can locate regions of historical population persistence and isolation, previous evidence of Virginia Province persistence is only provisional due to limited genetic sampling (i.e., mitochondrial DNA and five nuclear loci). To test alternative hypotheses of historical persistence versus the ephemerality of a northern Virginia Province population we used a RADseq generated dataset consisting of 11,708 single nucleotide polymorphisms (SNP) sampled from individuals collected from the eastern Gulf of Mexico to Long Island, NY. Concordant results from genomic analyses all infer three genetically divergent subpopulations, and strongly support Virginia Province inhabitants as a genetically diverged and a historically persistent ancestral gene pool. These results suggest that individuals that emerge in coastal areas during the warm season can be considered "local" and supports offshore migration during the colder months. This research demonstrates how a large number of genes sampled across a geographical range can capture the diversity of coalescent histories (across loci) while inferring population history. Moreover, these results clearly demonstrate the utility of population genomic data to infer peripheral subpopulation persistence in difficult-to-observe species.

Population genomics of marine fishes: next-generation prospects and challenges

Over the past few years, technological advances have facilitated giant leaps forward in our ability to generate genome-wide molecular data, offering exciting opportunities for gaining new insights into the ecology and evolution of species where genomic information is still limited. Marine fishes are valuable organisms for advancing our understanding of evolution on historical and contemporary time scales, and here we highlight areas in which research on these species is likely to be particularly important in the near future. These include possibilities for gaining insights into processes on ecological time scales, identifying genomic signatures associated with population divergence under gene flow, and determining the genetic basis of phenotypic traits. We also consider future challenges pertaining to the implementation of genome-wide coverage through next-generation sequencing and genotyping methods in marine fishes. Complications associated with fast decay of linkage disequilibrium, as expected for species with large effective population sizes, and the possibility that adaptation is associated with both soft selective sweeps and polygenic selection, leaving complex genomic signatures in natural populations, are likely to challenge future studies. However, the combination of high genome coverage and new statistical developments offers promising solutions. Thus, the next generation of studies is likely to truly facilitate the transition from population genetics to population genomics in marine fishes. This transition will advance our understanding of basic evolutionary processes and will offer new possibilities for conservation and management of valuable marine resources.

New dimensions of tropical diversity: an inordinate fondness for insect molecules, taxa, and trophic interactions

Most known insect species are involved in chemically mediated plant-insect multi-trophic interactions, and recent syntheses point to a substantial gap in our understanding of trophic interaction diversity, especially in the tropics. One approach to filling this gap is to examine relationships between genomic, metabolomic, taxonomic, and trophic interaction diversity via quantifying and comparing these dimensions of biodiversity at multiple scales. Innovative approaches to research on the origins and maintenance of tropical insect diversity should merge traditional approaches to natural history and taxonomy with modern measures of interaction diversity, genetic variation, and phytochemical diversity. These approaches will elucidate relationships between plant chemistry, specialization, climate, and different dimensions of biodiversity.

Seascape continuity plays an important role in determining patterns of spatial genetic structure in a coral reef fish

Detecting patterns of spatial genetic structure (SGS) can help identify intrinsic and extrinsic barriers to gene flow within metapopulations. For marine organisms such as coral reef fishes, identifying these barriers is critical to predicting evolutionary dynamics and demarcating evolutionarily significant units for conservation. In this study, we adopted an alternative hypothesis-testing framework to identify the patterns and predictors of SGS in the Caribbean reef fish Elacatinus lori. First, genetic structure was estimated using nuclear microsatellites and mitochondrial cytochrome b sequences. Next, clustering and network analyses were applied to visualize patterns of SGS. Finally, logistic regressions and linear mixed models were used to identify the predictors of SGS. Both sets of markers revealed low global structure: mitochondrial ΦST=0.12, microsatellite FST=0.0056. However, there was high variability among pairwise estimates, ranging from no differentiation between sites on contiguous reef (ΦST=0) to strong differentiation between sites separated by ocean expanses≥20 km (maximum ΦST=0.65). Genetic clustering and statistical analyses provided additional support for the hypothesis that seascape discontinuity, represented by oceanic breaks between patches of reef habitat, is a key predictor of SGS in E. lori. Notably, the estimated patterns and predictors of SGS were consistent between both sets of markers. Combined with previous studies of dispersal in E. lori, these results suggest that the interaction between seascape continuity and the dispersal kernel plays an important role in determining genetic connectivity within metapopulations.

Next-generation sequencing for molecular ecology: a caveat regarding pooled samples

We develop a model based on the Dirichlet-compound multinomial distribution (CMD) and Ewens sampling formula to predict the fraction of SNP loci that will appear fixed for alternate alleles between two pooled samples drawn from the same underlying population. We apply this model to next-generation sequencing (NGS) data from Baltic Sea herring recently published by (Corander et al., 2013, Molecular Ecology, 2931-2940), and show that there are many more fixed loci than expected in the absence of genetic structure. However, we show through coalescent simulations that the degree of population structure required to explain the fraction of alternatively fixed SNPs is extraordinarily high and that the surplus of fixed loci is more likely a consequence of limited representation of individual gene copies in the pooled samples, than it is of population structure. Our analysis signals that the use of NGS on pooled samples to identify divergent SNPs warrants caution. With pooled samples, it is hard to diagnose when an NGS experiment has gone awry; especially when NGS data on pooled samples are of low read depth with a limited number of individuals, it may be worthwhile to temper claims of unexpected population differentiation from pooled samples, pending verification with more reliable methods or stricter adherence to recommended sampling designs for pooled sequencing e.g. Futschik & Schlötterer 2010, Genetics, 186, 207; Gautier et al., 2013a, Molecular Ecology, 3766-3779). Analysis of the data and diagnosis of problems is easier and more reliable (and can be less costly) with individually barcoded samples. Consequently, for some scenarios, individual barcoding may be preferable to pooling of samples.

Genetic diversity and population structure of Mongolian domestic Bactrian camels (Camelus bactrianus)

The tradition of animal husbandry in the context of a nomadic lifestyle has been of great significance in the Mongolian society. Both Bactrian camels and horses have been invaluable for the survival and development of human activities in the harsh arid environment of the Mongolian steppe. As camels offer unique and sustainable opportunities for livestock production in marginal agro-ecological zones, we investigated the current genetic diversity of three local Mongolian camel breeds and compared their levels of variation with common native Mongolian camels distributed throughout the country. Based on mitochondrial and nuclear markers, we found levels of genetic diversity in Mongolian populations similar to that reported for Chinese Bactrian camels and for dromedaries. Little differentiation was detected between single breeds, except for a small group originating from the northwestern Mongolian Altai. We found neither high inbreeding levels in the different breeds nor evidence for a population decline. Although the Mongolian camel census size has severely declined over the past 20 years, our analyses suggest that there still exists a stable population with adequate genetic variation for continued sustainable utilization.

Populations of a cyprinid fish are self-sustaining despite widespread feminization of males

BACKGROUND

Treated effluents from wastewater treatment works can comprise a large proportion of the flow of rivers in the developed world. Exposure to these effluents, or the steroidal estrogens they contain, feminizes wild male fish and can reduce their reproductive fitness. Long-term experimental exposures have resulted in skewed sex ratios, reproductive failures in breeding colonies, and population collapse. This suggests that environmental estrogens could threaten the sustainability of wild fish populations.

RESULTS

Here we tested this hypothesis by examining population genetic structures and effective population sizes (N(e)) of wild roach (Rutilus rutilus L.) living in English rivers contaminated with estrogenic effluents. N(e) was estimated from DNA microsatellite genotypes using approximate Bayesian computation and sibling assignment methods. We found no significant negative correlation between N(e) and the predicted estrogen exposure at 28 sample sites. Furthermore, examination of the population genetic structure of roach in the region showed that some populations have been confined to stretches of river with a high proportion of estrogenic effluent for multiple generations and have survived, apparently without reliance on immigration of fish from less polluted sites.

CONCLUSIONS

These results demonstrate that roach populations living in some effluent-contaminated river stretches, where feminization is widespread, are self-sustaining. Although we found no evidence to suggest that exposure to estrogenic effluents is a significant driving factor in determining the size of roach breeding populations, a reduction in N(e) of up to 65% is still possible for the most contaminated sites because of the wide confidence intervals associated with the statistical model.

High degree of genetic differentiation in marine three-spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three-spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene-based genome-scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection-as determined with several outlier detection methods-was low (F(ST) = 0.033 or 0.011, respectively), whereas average FST for directionally selected markers was significantly higher (F(ST) = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.

Oceanographic connectivity and environmental correlates of genetic structuring in Atlantic herring in the Baltic Sea

Marine fish often show little genetic structuring in neutral marker genes, and Atlantic herring (Clupea harengus) in the Baltic Sea are no exception; historically, very low levels of population differentiation (F ST ≍ 0.002) have been found, despite a high degree of interpopulation environmental heterogeneity in salinity and temperature. Recent exome sequencing and SNP studies have however shown that many loci are under selection in this system. Here, we combined population genetic analyses of a large number of transcriptome-derived microsatellite markers with oceanographic modelling to investigate genetic differentiation and connectivity in Atlantic herring at a relatively fine scale within the Baltic Sea. We found evidence for weak but robust and significant genetic structuring (F ST = 0.008) explainable by oceanographic connectivity. Genetic differentiation was also associated with site differences in temperature and salinity, with the result driven by the locus Her14 which appears to be under directional selection (F ST = 0.08). The results show that Baltic herring are genetically structured within the Baltic Sea, and highlight the role of oceanography and environmental factors in explaining this structuring. The results also have implications for the management of herring fisheries, the most economically important fishery in the Baltic Sea, suggesting that the current fisheries management units may be in need of revision.

RAD-Seq derived markers flank the shell colour and banding loci of the Cepaea nemoralis supergene

Studies on the classic shell colour and banding polymorphism of the land snail Cepaea played a crucial role in establishing the importance of natural selection in maintaining morphological variation. Cepaea is also a pre-eminent model for ecological genetics because the outward colour and banding phenotype is entirely genetically determined, primarily by a 'supergene' of at least five loci. Unfortunately, progress in understanding the evolution and maintenance of the Cepaea polymorphism stalled, partly because of a lack of genetic markers. With a view to re-establish Cepaea as a prominent model of molecular ecology, we made six laboratory crosses of Cepaea nemoralis, five of which segregated for shell ground colour (C) and the presence or absence of bands (B). First, scoring of colour and banding in 323 individuals found no recombination between the C and B loci of the supergene. Second, using restriction site-associated DNA sequencing (RAD-Seq) of two parents and 22 offspring, we identified 44 anonymous markers putatively linked to the colour (C) and banding (B) loci. The genotype of eleven of the most promising RAD-Seq markers was independently validated in the same 22 offspring, then up to a further 146 offspring were genotyped. The closest RAD-Seq markers scored are within ~0.6 centimorgan (cM) of the C-B supergene linkage group, with the combined loci together forming a 35.8 cM linkage map of markers that flank both sides of the Cepaea C-B supergene.