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.

To address the role of insect herbivores in adaptation of plant populations and the persistence of selection through succession, we manipulated herbivory in a long-term field experiment. We suppressed insects in half of 16 plots over nine years and examined the genotypic structure and chemical defense of common dandelion (Taraxacum officinale), a naturally colonizing perennial apomictic plant. Insect suppression doubled dandelion abundance in the first few years, but had negligible effects thereafter. Using microsatellite DNA markers, we genotyped >2500 plants and demonstrate that insect suppression altered the genotypic composition of plots in both sampling years. Phenotypic and genotypic estimates of defensive terpenes and phenolics from the field plots allowed us to infer phenotypic plasticity and the response of dandelion populations to insect-mediated natural selection. The effects of insect suppression on plant chemistry were, indeed, driven both by plasticity and plant genotypic identity. In particular, di-phenolic inositol esters were more abundant in plots exposed to herbivory (due to the genotypic composition of the plots) and were also induced in response to herbivory. This field experiment thus demonstrates evolutionary sorting of plant genotypes in response to insect herbivores that was in same direction as the plastic defensive response within genotypes.

We developed nine polymorphic microsatellite loci for evening primrose (Oenothera biennis). These loci have two to 18 alleles per locus and observed heterozygosities ranging from 0 to 0.879 in a sample of 34 individuals. In a pattern consistent with the functionally asexual reproductive system of this species, 17/36 pairs of loci revealed significant linkage disequilibrium and three loci showed significant deviations from Hardy-Weinberg equilibrium. The loci will be informative in identifying genotypes in multigenerational field studies to assess changes in genotype frequencies.

BackgroundComparative genomic and/or transcriptomic analyses involving elasmobranchs remain limited, with genome level comparisons of the elasmobranch immune system to that of higher vertebrates, non-existent. This paper reports a comparative RNA-seq analysis of heart tissue from seven species, including four elasmobranchs and three teleosts, focusing on immunity, but concomitantly seeking to identify genetic similarities shared by the two lamnid sharks and the single billfish in our study, which could be linked to convergent evolution of regional endothermy.ResultsAcross seven species, we identified an average of 10,877 Swiss-Prot annotated genes from an average of 32,474 open reading frames within each species’ heart transcriptome. About half of these genes were shared between all species while the remainder included functional differences between our groups of interest (elasmobranch vs. teleost and endotherms vs. ectotherms) as revealed by Gene Ontology (GO) and selection analyses. A repeatedly represented functional category, in both the uniquely expressed elasmobranch genes (total of 259) and the elasmobranch GO enrichment results, involved antibody-mediated immunity, either in the recruitment of immune cells (Fc receptors) or in antigen presentation, including such terms as “antigen processing and presentation of exogenous peptide antigen via MHC class II”, and such genes as MHC class II, HLA-DPB1. Molecular adaptation analyses identified three genes in elasmobranchs with a history of positive selection, including legumain (LGMN), a gene with roles in both innate and adaptive immunity including producing antigens for presentation by MHC class II. Comparisons between the endothermic and ectothermic species revealed an enrichment of GO terms associated with cardiac muscle contraction in endotherms, with 19 genes expressed solely in endotherms, several of which have significant roles in lipid and fat metabolism.ConclusionsThis collective comparative evidence provides the first multi-taxa transcriptomic-based perspective on differences between elasmobranchs and teleosts, and suggests various unique features associated with the adaptive immune system of elasmobranchs, pointing in particular to the potential importance of MHC Class II. This in turn suggests that expanded comparative work involving additional tissues, as well as genome sequencing of multiple elasmobranch species would be productive in elucidating the regulatory and genome architectural hallmarks of elasmobranchs.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3411-x) contains supplementary material, which is available to authorized users.

Table S1. Number of individuals collected and sequenced per site. The number of individuals sequenced indicates the number of individuals with sequences retained after all quality filtering. Site ID # Samples Collected # Samples sequenced