Christopher M. Hollenbeck

dDocent: a RADseq, variant-calling pipeline designed for population genomics of non-model organisms

Restriction-site associated DNA sequencing (RADseq) has become a powerful and useful approach for population genomics. Currently, no software exists that utilizes both paired-end reads from RADseq data to efficiently produce population-informative variant calls, especially for non-model organisms with large effective population sizes and high levels of genetic polymorphism. dDocent is an analysis pipeline with a user-friendly, command-line interface designed to process individually barcoded RADseq data (with double cut sites) into informative SNPs/Indels for population-level analyses. The pipeline, written in BASH, uses data reduction techniques and other stand-alone software packages to perform quality trimming and adapter removal, de novo assembly of RAD loci, read mapping, SNP and Indel calling, and baseline data filtering. Double-digest RAD data from population pairings of three different marine fishes were used to compare dDocent with Stacks, the first generally available, widely used pipeline for analysis of RADseq data. dDocent consistently identified more SNPs shared across greater numbers of individuals and with higher levels of coverage. This is due to the fact that dDocent quality trims instead of filtering, incorporates both forward and reverse reads (including reads with INDEL polymorphisms) in assembly, mapping, and SNP calling. The pipeline and a comprehensive user guide can be found at http://dDocent.wordpress.com.

Selection and sex-biased dispersal in a coastal shark: the influence of philopatry on adaptive variation

Sex‐biased dispersal is expected to homogenize nuclear genetic variation relative to variation in genetic material inherited through the philopatric sex. When site fidelity occurs across a heterogeneous environment, local selective regimes may alter this pattern. We assessed spatial patterns of variation in nuclear‐encoded, single nucleotide polymorphisms (SNPs) and sequences of the mitochondrial control region in bonnethead sharks (Sphyrna tiburo), a species thought to exhibit female philopatry, collected from summer habitats used for gestation. Geographic patterns of mtDNA haplotypes and putatively neutral SNPs confirmed female philopatry and male‐mediated gene flow along the northeastern coast of the Gulf of Mexico. A total of 30 outlier SNP loci were identified; alleles at over half of these loci exhibited signatures of latitude‐associated selection. Our results indicate that in species with sex‐biased dispersal, philopatry can facilitate sorting of locally adaptive variation, with the dispersing sex facilitating movement of potentially adaptive variation among locations and environments.

Contemporary population structure and post-glacial genetic demography in a migratory marine species, the blacknose shark, Carcharhinus acronotus

Patterns of population structure and historical genetic demography of blacknose sharks in the western North Atlantic Ocean were assessed using variation in nuclear‐encoded microsatellites and sequences of mitochondrial (mt)DNA. Significant heterogeneity and/or inferred barriers to gene flow, based on microsatellites and/or mtDNA, revealed the occurrence of five genetic populations localized to five geographic regions: the southeastern U.S Atlantic coast, the eastern Gulf of Mexico, the western Gulf of Mexico, Bay of Campeche in the southern Gulf of Mexico and the Bahamas. Pairwise estimates of genetic divergence between sharks in the Bahamas and those in all other localities were more than an order of magnitude higher than between pairwise comparisons involving the other localities. Demographic modelling indicated that sharks in all five regions diverged after the last glacial maximum and, except for the Bahamas, experienced post‐glacial, population expansion. The patterns of genetic variation also suggest that the southern Gulf of Mexico may have served as a glacial refuge and source for the expansion. Results of the study demonstrate that barriers to gene flow and historical genetic demography contributed to contemporary patterns of population structure in a coastal migratory species living in an otherwise continuous marine habitat. The results also indicate that for many marine species, failure to properly characterize barriers in terms of levels of contemporary gene flow could in part be due to inferences based solely on equilibrium assumptions. This could lead to erroneous conclusions regarding levels of connectivity in species of conservation concern.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 December 2011 - 31 January 2012: PERMANENT GENETIC RESOURCES NOTE

This article documents the addition of 473 microsatellite marker loci and 71 pairs of single‐nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Barteria fistulosa, Bombus morio, Galaxias platei, Hematodinium perezi, Macrocentrus cingulum Brischke (a.k.a. M. abdominalis Fab., M. grandii Goidanich or M. gifuensis Ashmead), Micropogonias furnieri, Nerita melanotragus, Nilaparvata lugens Stål, Sciaenops ocellatus, Scomber scombrus, Spodoptera frugiperda and Turdus lherminieri. These loci were cross‐tested on the following species: Barteria dewevrei, Barteria nigritana, Barteria solida, Cynoscion acoupa, Cynoscion jamaicensis, Cynoscion leiarchus, Cynoscion nebulosus, Cynoscion striatus, Cynoscion virescens, Macrodon ancylodon, Menticirrhus americanus, Nilaparvata muiri and Umbrina canosai. This article also documents the addition of 116 sequencing primer pairs for Dicentrarchus labrax.

Haplotyping RAD loci: an efficient method to filter paralogs and account for physical linkage.

Next‐generation sequencing of reduced‐representation genomic libraries provides a powerful methodology for genotyping thousands of single‐nucleotide polymorphisms (SNPs) among individuals of nonmodel species. Utilizing genotype data in the absence of a reference genome, however, presents a number of challenges. One major challenge is the trade‐off between splitting alleles at a single locus into separate clusters (loci), creating inflated homozygosity, and lumping multiple loci into a single contig (locus), creating artefacts and inflated heterozygosity. This issue has been addressed primarily through the use of similarity cut‐offs in sequence clustering. Here, two commonly employed, postclustering filtering methods (read depth and excess heterozygosity) used to identify incorrectly assembled loci are compared with haplotyping, another postclustering filtering approach. Simulated and empirical data sets were used to demonstrate that each of the three methods separately identified incorrectly assembled loci; more optimal results were achieved when the three methods were applied in combination. The results confirmed that including incorrectly assembled loci in population‐genetic data sets inflates estimates of heterozygosity and deflates estimates of population divergence. Additionally, at low levels of population divergence, physical linkage between SNPs within a locus created artificial clustering in analyses that assume markers are independent. Haplotyping SNPs within a locus effectively neutralized the physical linkage issue without having to thin data to a single SNP per locus. We introduce a Perl script that haplotypes polymorphisms, using data from single or paired‐end reads, and identifies potentially problematic loci.

Rapid genetic and morphologic divergence between captive and wild populations of the endangered Leon Springs pupfish, Cyprinodon bovinus

The Leon Springs pupfish (Cyprinodon bovinus) is an endangered species currently restricted to a single desert spring and a separate captive habitat in southwestern North America. Following establishment of the captive population from wild stock in 1976, the wild population has undergone natural population size fluctuations, intentional culling to purge genetic contamination from an invasive congener (Cyprinodon variegatus) and augmentation/replacement of wild fish from the captive stock. A severe population decline following the most recent introduction of captive fish prompted us to examine whether the captive and wild populations have differentiated during the short time they have been isolated from one another. If so, the development of divergent genetic and/or morphologic traits between populations could contribute to a diminished ability of fish from one location to thrive in the other. Examination of genomewide single nucleotide polymorphisms and morphologic variation revealed no evidence of residual C. variegatus characteristics in contemporary C. bovinus samples. However, significant genetic and morphologic differentiation was detected between the wild and captive populations, some of which might reflect local adaptation. Our results indicate that genetic and physical characteristics can diverge rapidly between isolated subdivisions of managed populations, potentially compromising the value of captive stock for future supplementation efforts. In the case of C. bovinus, our findings underscore the need to periodically inoculate the captive population with wild genetic material to help mitigate genetic, and potentially morphologic, divergence between them and also highlight the utility of parallel morphologic and genomic evaluation to inform conservation management planning.

Phylogenetic relationships of tropical eastern Pacific snappers (Lutjanidae) inferred from mitochondrial DNA sequences

A total of 2,206 base pairs of coding sequences of mitochondrial DNA from nine of 10 lutjanine (Lutjanidae) species endemic to the eastern Pacific Ocean were used in a phylogenetic analysis to infer species relationships among lutjanines on either side of the Isthmus of Panama. Monotypic Hoplopagrus (guentherii) appears to be at the base of all lutjanines assayed to date and represents the first of five, eastward dispersal events of lutjanines from the Indo-Pacific to the Americas. Diversification of lutjanines in the Americas appears to have occurred primarily in the western Atlantic Ocean (Caribbean Sea) and been aided by an early Miocene eastward-flowing current and by more optimal lutjanine habitat in the warmer Caribbean waters. Five geminate species pairs in five, strongly supported clades were identified. In three, the trans-isthmian geminate species in the eastern Pacific belonged to clades in which basal and other members of the clade are endemic to the western Atlantic, suggesting that the ancestor to each species geminate pair dispersed westward into the eastern Pacific prior to closure of connections between the two oceans; the opposite occurred in one of the clades. The fifth clade contained only two species, one on either side of the Isthmus. Diversification in life-history and other traits in each geminate species pair appears to have been slight. Five lutjanine species possessing horizontal blue stripes (the blue-lined complex) and distributed from the western Indian Ocean to the eastern Pacific formed a strongly supported monophyletic clade. The lone blue-lined species in the eastern Pacific (Lutjanus viridis) could be conspecific with its highly dispersive sister species (Lutjanus kasmira) from the Indo-Pacific. The chronology of lutjanine diversification and closure of seawater exchange between the two oceans is briefly discussed.

A comparison of neutral and immune genetic variation in Atlantic salmon, Salmo salar L. in Chilean aquaculture facilities.

Genetic diversity was assessed in samples of cultured Atlantic salmon, Salmo salar L., obtained from facilities in Chile between 2005 and 2010, a period of time during which the infectious pathogens Infectious Salmon Anemia (ISA) virus, Caligus rogercresseyi (sea lice), and Piscirickettsia salmonis (salmon rickettsial syndrome) were common. Two panels of microsatellite markers were utilized: one with microsatellites with no known gene associations (neutral) and one featuring microsatellites linked to putative immune-related genes (immune-related). Allelic richness and gene diversity across samples were significantly greater in neutral loci as compared to immune-related loci. Both diversity measures were homogeneous among samples for immune-related loci and heterogeneous among samples for neutral loci. Immune-related loci were identified as FST outliers in pairwise comparisons of samples at a 10-fold higher frequency than neutral loci. These results indicate that neutral and immune-related portions of the Atlantic salmon genome may have differed in response to the gauntlet of pathogens and that monitoring of specific, well characterized immune-related loci as well as neutral loci in cultured species could be useful when disease control and prevention is a goal.

Linkage Mapping and Comparative Genomics of Red Drum (Sciaenops ocellatus) Using Next-Generation Sequencing

Developments in next-generation sequencing allow genotyping of thousands of genetic markers across hundreds of individuals in a cost-effective manner. Because of this, it is now possible to rapidly produce dense genetic linkage maps for nonmodel species. Here, we report a dense genetic linkage map for red drum, a marine fish species of considerable economic importance in the southeastern United States and elsewhere. We used a prior microsatellite-based linkage map as a framework and incorporated 1794 haplotyped contigs derived from high-throughput, reduced representation DNA sequencing to produce a linkage map containing 1794 haplotyped restriction-site associated DNA (RAD) contigs, 437 anonymous microsatellites, and 44 expressed sequence-tag-linked microsatellites (EST-SSRs). A total of 274 candidate genes, identified from transcripts from a preliminary hydrocarbon exposure study, were localized to specific chromosomes, using a shared synteny approach. The linkage map will be a useful resource for red drum commercial and restoration aquaculture, and for better understanding and managing populations of red drum in the wild.