Estimating microhaplotype allele frequencies from low-coverage or pooled sequencing data.

BACKGROUND: Microhaplotypes have the potential to be more cost-effective than SNPs for applications that require genetic panels of highly variable loci. However, development of microhaplotype panels is hindered by a lack of methods for estimating microhaplotype allele frequency from low-coverage whole genome sequencing or pooled sequencing (pool-seq) data. RESULTS: We developed new methods for estimating microhaplotype allele frequency from low-coverage whole genome sequence and pool-seq data. We validated these methods using datasets from three non-model organisms. These methods allowed estimation of allele frequency and expected heterozygosity at depths routinely achieved from pooled sequencing. CONCLUSIONS: These new methods will allow microhaplotype panels to be designed using low-coverage WGS and pool-seq data to discover and evaluate candidate loci. The python script implementing the two methods and documentation are available at https://www.github.com/delomast/mhFromLowDepSeq .

Interrelationships and biogeography of the New World pufferfish genus Sphoeroides (Tetraodontiformes: Tetraodontidae) inferred using ultra-conserved DNA elements.

Colonization of the New World by marine taxa has been hypothesized to have occurred through the Tethys Sea or by crossing the East Pacific Barrier. To better understand patterns and timing of diversification, geological events can be coupled with time calibrated phylogenetic hypotheses to infer major drivers of diversification. Phylogenetic relationships among members of Sphoeroides, a genus of four toothed pufferfishes (Tetraodontiformes: Tetraodontidae) which are found nearly exclusively in the New World (eastern Pacific and western Atlantic), were reconstructed using sequences from ultra-conserved DNA elements, nuclear markers with clear homology among many vertebrate taxa. Hypotheses derived from concatenated maximum-likelihood and species tree summary methods support a paraphyletic Sphoeroides, with Colomesus deeply nested within the genus. Analyses also revealed S. pachygaster, a pelagic species with a cosmopolitan distribution, as the sister taxon to the remainder of Sphoeroides and recovered distinct lineages within S. pachygaster, indicating that this cosmopolitan species may represent a species complex. Ancestral range reconstruction may suggest the genus colonized the New World through the eastern Pacific before diversifying in the western Atlantic, though date estimates for these events are uncertain due to the lack of reliable fossil record for the genus.

Implications of Large-Effect Loci for Conservation: A Review and Case Study with Pacific Salmon.

The increasing feasibility of assembling large genomic datasets for non-model species presents both opportunities and challenges for applied conservation and management. A popular theme in recent studies is the search for large-effect loci that explain substantial portions of phenotypic variance for a key trait(s). If such loci can be linked to adaptations, 2 important questions arise: 1) Should information from these loci be used to reconfigure conservation units (CUs), even if this conflicts with overall patterns of genetic differentiation? 2) How should this information be used in viability assessments of populations and larger CUs? In this review, we address these questions in the context of recent studies of Chinook salmon and steelhead (anadromous form of rainbow trout) that show strong associations between adult migration timing and specific alleles in one small genomic region. Based on the polygenic paradigm (most traits are controlled by many genes of small effect) and genetic data available at the time showing that early-migrating populations are most closely related to nearby late-migrating populations, adult migration differences in Pacific salmon and steelhead were considered to reflect diversity within CUs rather than separate CUs. Recent data, however, suggest that specific alleles are required for early migration, and that these alleles are lost in populations where conditions do not support early-migrating phenotypes. Contrasting determinations under the US Endangered Species Act and the State of California's equivalent legislation illustrate the complexities of incorporating genomics data into CU configuration decisions. Regardless how CUs are defined, viability assessments should consider that 1) early-migrating phenotypes experience disproportionate risks across large geographic areas, so it becomes important to identify early-migrating populations that can serve as reliable sources for these valuable genetic resources; and 2) genetic architecture, especially the existence of large-effect loci, can affect evolutionary potential and adaptability.

These aren't the loci you'e looking for: Principles of effective SNP filtering for molecular ecologists.

Sequencing reduced-representation libraries of restriction site-associated DNA (RADseq) to identify single nucleotide polymorphisms (SNPs) is quickly becoming a standard methodology for molecular ecologists. Because of the scale of RADseq data sets, putative loci cannot be assessed individually, making the process of filtering noise and correctly identifying biologically meaningful signal more difficult. Artefacts introduced during library preparation and/or bioinformatic processing of SNP data can create patterns that are incorrectly interpreted as indicative of population structure or natural selection. Therefore, it is crucial to carefully consider types of errors that may be introduced during laboratory work and data processing, and how to minimize, detect and remove these errors. Here, we discuss issues inherent to RADseq methodologies that can result in artefacts during library preparation and locus reconstruction resulting in erroneous SNP calls and, ultimately, genotyping error. Further, we describe steps that can be implemented to create a rigorously filtered data set consisting of markers accurately representing independent loci and compare the effect of different combinations of filters on four RAD data sets. At last, we stress the importance of publishing raw sequence data along with final filtered data sets in addition to detailed documentation of filtering steps and quality control measures.

Molecular phylogenetics of New World searobins (Triglidae; Prionotinae).

Phylogenetic relationships among members of the New World searobin genera Bellator and Prionotus (Family Triglidae, Subfamily Prionotinae) and among other searobins in the families Triglidae and Peristediidae were investigated using both mitochondrial and nuclear DNA sequences. Phylogenetic hypotheses derived from maximum likelihood and Bayesian methodologies supported a monophyletic Prionotinae that included four well resolved clades of uncertain relationship; three contained species in the genus Prionotus and one contained species in the genus Bellator. Bellator was always recovered within the genus Prionotus, a result supported by post hoc model testing. Two nominal species of Prionotus (P. alatus and P. paralatus) were not recovered as exclusive lineages, suggesting the two may comprise a single species. Phylogenetic hypotheses also supported a monophyletic Triglidae but only if armored searobins (Family Peristediidae) were included. A robust morphological assessment is needed to further characterize relationships and suggest classification of clades within Prionotinae; for the time being we recommend that Bellator be considered a synonym of Prionotus. Relationships between armored searobins (Family Peristediidae) and searobins (Family Triglidae) and relationships within Triglidae also warrant further study.

Contrasting patterns of sequence variation in steelhead populations reflect distinct evolutionary processes.

Multiple evolutionary processes influence genome-wide allele frequencies and quantifying effects of genetic drift, and multiple forms of selection remain challenging in natural populations. Here, we investigate variation at major effect loci in contrast to patterns of neutral drift across a wide collection of steelhead (Oncorhynchus mykiss) populations that have declined in abundance due to anthropogenic impacts. Whole-genome resequencing of 74 populations of steelhead revealed genome-wide patterns (~8 million SNPs) consistent with expected neutral population structure. However, allelic variation at major effect loci associated with adult migration timing (chromosome 28: GREB1L/ROCK1) and age at maturity (chromosome 25: SIX6) reflected how selection has acted on phenotypic variation in contrast with neutral structure. Variation at major effect loci was influenced by evolutionary processes with differing signals between the strongly divergent Coastal and Inland lineages, while allele frequencies within and among populations within the Inland lineage have been driven by local natural selection as well as recent anthropogenic influences. Recent anthropogenic effects appeared to have influenced the frequency of major effect alleles including artificial selection for specific traits in hatchery stocks with subsequent gene flow into natural populations. Selection from environmental factors at various scales has also likely influenced variation for major effect alleles. These results reveal evolutionary mechanisms that influence allele frequencies at major effect loci that are critical for conservation of phenotypic traits and life history variation of this protected species.

Genetic variation associated with adult migration timing in lineages of Steelhead and Chinook Salmon in the Columbia River.

With the discovery of a major effect region (GREB1L, ROCK1) for adult migration timing in genomes of both Chinook Salmon and Steelhead, several subsequent studies have investigated the effect size and distribution of early and late migration alleles among populations in the Columbia River. Here, we synthesize the results of these studies for the major lineages of Chinook Salmon and Steelhead that include highly distinct groups in the interior Columbia River that exhibit atypical life histories from most coastal lineage populations of these two species. Whole-genome studies with high marker density have provided extensive insight into SNPs most associated with adult migration timing, and suites of markers for each species have been genotyped in large numbers of individuals to further validate phenotypic effects. For Steelhead, the largest phenotypic effect sizes have been observed in the coastal lineage (36% of variation for passage timing at Bonneville Dam; 43% of variation for tributary arrival timing) compared to the inland lineage (7.5% of variation for passage timing at Bonneville Dam; 8.4% of variation for tributary arrival timing) that overwinter in freshwater prior to spawning. For Chinook Salmon, large effect sizes have been observed in all three lineages for multiple adult migration phenotypes (Coastal lineage: percentage of variation of 27.9% for passage timing at Bonneville Dam, 28.7% for arrival timing for spawning; Interior ocean type: percentage of variation of 47.6% for passage timing at Bonneville Dam, 39.6% for tributary arrival timing, 77.9% for arrival timing for spawning; Interior stream type: percentage of variation of 35.3% for passage at Bonneville Dam, 9.8% for tributary arrival timing, 4.7% for arrival timing for spawning). Together, these results have extended our understanding of genetic variation associated with life history diversity in distinct populations of the Columbia River, however, much research remains necessary to determine the causal mechanism for this major effect region on migration timing in these species.

A genomic region associated with iteroparous spawning phenology is linked with age-at-maturity in female steelhead trout.

Age-at-maturity and iteroparity are two life history variations of steelhead trout (Oncorhynchus mykiss) that are believed to increase population resilience and stability. While repeat-spawning individuals are thought to have historically made up a substantial portion of the reproductive population in the Columbia River and the majority of females still attempt outmigration as kelts, return rates of repeat-spawner are low throughout the basin and below 1% for the furthest migrating stocks. Notably, outmigrating adults exhibit variation in rematuration phenology, displaying either "consecutive" (reproduce immediately the following season) or "skip" (delay spawning for future seasons) spawning patterns. Here, we use low coverage whole genome sequencing of consecutive versus skip spawning female Columbia River steelhead from two populations to test for genomic differences between these two iteroparous phenotypes. We identified genomic regions on several chromosomes which were associated with the phenology of iteroparity, including a region on chromosome 25 containing two genes, estradiol receptor beta (ERß) and glycoprotein hormone beta-5 (GPHB5), which, in mammals, are estrogen-sensitive and expressed in reproductive tissues. Allele frequencies in this ERß/GPHB5 region differed among female steelhead of different age at maturity, but not males. These genes also shared an island of linkage disequilibrium with the SIX6 gene, 600Kbp away on the same chromosome, a region of known association with age-at-maturity. These observations contribute to growing evidence that age-at-maturity and the phenology of iteroparity are determined by overlapping physiological processes and genetic pathways.

Efficient species identification for Pacific salmon genetic monitoring programs.

Genetic monitoring of Pacific salmon in the Columbia River basin provides crucial information to fisheries managers that is otherwise challenging to obtain using traditional methods. Monitoring programs such as genetic stock identification (GSI) and parentage-based tagging (PBT) involve genotyping tens of thousands of individuals annually. Although rare, these large sample collections inevitably include misidentified species, which exhibit low genotyping success on species-specific Genotyping-in-Thousands by sequencing (GT-seq) panels. For laboratories involved in large-scale genotyping efforts, diagnosing non-target species and reassigning them to the appropriate monitoring program can be costly and time-consuming. To address this problem, we identified 19 primer pairs that exhibit consistent cross-species amplification among salmonids and contain 51 species informative variants. These genetic markers reliably discriminate among 11 salmonid species and two subspecies of Cutthroat Trout and have been included in species-specific GT-seq panels for Chinook Salmon, Coho Salmon, Sockeye Salmon, and Rainbow Trout commonly used for Pacific salmon genetic monitoring. The majority of species-informative amplicons (16) were newly identified from the four existing GT-seq panels, thus demonstrating a low-cost approach to species identification when using targeted sequencing methods. A species-calling script was developed that is tailored for routine GT-seq genotyping pipelines and automates the identification of non-target species. Following extensive testing with empirical and simulated data, we demonstrated that the genetic markers and accompanying script accurately identified species and are robust to missing genotypic data and low-frequency, shared polymorphisms among species. Finally, we used these tools to identify Coho Salmon incidentally caught in the Columbia River Chinook Salmon sport fishery and used PBT to determine their hatchery of origin. These molecular and computing resources provide a valuable tool for Pacific salmon conservation in the Columbia River basin and demonstrate a cost-effective approach to species identification for genetic monitoring programs.

Multi-locus species tree for the Amazonian peacock basses (Cichlidae: Cichla): emergent phylogenetic signal despite limited nuclear variation.

The inference of phylogenies of closely related species is obstructed by phenomena such as porous species boundaries and deep coalescence, and is often exacerbated by low levels of nucleotide variation among most loci surveyed in phylogenetic studies. We investigated the utility of twenty-one nuclear loci that had a range of 5-40 (median of 14) variable sites per locus to estimate the phylogeny of the genus Cichla, a group of 15 Neotropical cichlid fishes that began to diverge in the early to mid Miocene. We found that under a concatenated approach, the least variable loci, while contributing less to the overall phylogenetic signal (posterior node support), nevertheless provided information that increased support for the final tree. Moreover, this was not a result of misdirection by mutational noise, as the inference from all data was far superior to those from reduced datasets (those with more variable loci) in terms of the relative precision of posterior tree space. Phylogenetic methods that allowed each locus to have a separate genealogy, including Bayesian concordance analysis and a multispecies coalescent model, provided phylogenies that were also compatible with the concatenated tree in terms of the eight recently delimited species of Cichla, albeit with somewhat diminished support for some branches. In contrast, described species that still regularly exchange genes showed unstable relationships among analyses: not a surprising result from analyses that assume that gene tree heterogeneity results from incomplete lineage sorting and not gene flow. Importantly, we also observed that the confidence intervals for node ages in the coalescent analyses were quite wide, and likely susceptible to influence of the prior on node density (e.g. birth-death).

PoolParty2: An integrated pipeline for analysing pooled or indexed low-coverage whole-genome sequencing data to discover the genetic basis of diversity.

Whole-genome sequencing data allow survey of variation from across the genome, reducing the constraint of balancing genome sub-sampling with estimating recombination rates and linkage between sampled markers and target loci. As sequencing costs decrease, low-coverage whole-genome sequencing of pooled or indexed-individual samples is commonly utilized to identify loci associated with phenotypes or environmental axes in non-model organisms. There are, however, relatively few publicly available bioinformatic pipelines designed explicitly to analyse these types of data, and fewer still that process the raw sequencing data, provide useful metrics of quality control and then execute analyses. Here, we present an updated version of a bioinformatics pipeline called PoolParty2 that can effectively handle either pooled or indexed DNA samples and includes new features to improve computational efficiency. Using simulated data, we demonstrate the ability of our pipeline to recover segregating variants, estimate their allele frequencies accurately, and identify genomic regions harbouring loci under selection. Based on the simulated data set, we benchmark the efficacy of our pipeline with another bioinformatic suite, angsd, and illustrate the compatibility and complementarity of these suites using angsd to generate genotype likelihoods as input for identifying linkage outlier regions using alignment files and variants provided by PoolParty2. Finally, we apply our updated pipeline to an empirical dataset of low-coverage whole genomic data from population samples of Columbia River steelhead trout (Oncorhynchus mykiss), results from which demonstrate the genomic impacts of decades of artificial selection in a prominent hatchery stock. Thus, we not only demonstrate the utility of PoolParty2 for genomic studies that combine sequencing data from multiple individuals, but also illustrate how it compliments other bioinformatics resources such as angsd.

Simultaneous delimitation of species and quantification of interspecific hybridization in Amazonian peacock cichlids (genus cichla) using multi-locus data.

BACKGROUND: Introgression likely plays a significant role in evolution, but understanding the extent and consequences of this process requires a clear identification of species boundaries in each focal group. The delimitation of species, however, is a contentious endeavor. This is true not only because of the inadequacy of current tools to identify species lineages, but also because of the inherent ambiguity between natural populations and species paradigms. The result has been a debate about the supremacy of various species concepts and criteria. Here, we utilized multiple separate sources of molecular data, mtDNA, nuclear sequences, and microsatellites, to delimit species under a polytypic species concept (PTSC) and estimate the frequency and genomic extent of introgression in a Neotropical genus of cichlid fishes (Cichla). We compared our inferences of species boundaries and introgression under this paradigm to those when species are identified under a diagnostic species concept (DSC). RESULTS: We find that, based on extensive molecular data and an inclusive species concept, 8 separate biological entities should be recognized rather than the 15 described species of Cichla. Under the PTSC, fewer individuals are expected to exhibit hybrid ancestry than under the DSC (~2% vs. ~12%), but a similar number of the species exhibit introgression from at least one other species (75% vs. 60%). Under either species concept, the phylogenetic breadth of introgression in this group is notable, with both sister species and species from different major mtDNA clades exhibiting introgression. CONCLUSIONS: Introgression was observed to be a widespread phenomenon for delimited species in this group. While several instances of introgressive hybridization were observed in anthropogenically altered habitats, most were found in undisturbed natural habitats, suggesting that introgression is a natural but ephemeral part of the evolution of many tropical species. Nevertheless, even transient introgression may facilitate an increase in genetic diversity or transfer of adaptive mutations that have important consequences in the evolution of tropical biodiversity.

Species-level phylogeny of 'Satan's perches' based on discordant gene trees (Teleostei: Cichlidae: Satanoperca Günther 1862).

Neotropical rivers are home to the largest assemblage of freshwater fishes, but little is known about the phylogeny of these fishes at the species level using multi-locus molecular markers. Here, we present a phylogeny for all known species of the genus Satanoperca, a widespread group of Neotropical cichlid fishes, based on analysis of six unlinked genetic loci. To test nominal and proposed species limits for this group, we surveyed mtDNA sequence variation among 320 individuals representing all know species. Most nominal species were supported by this approach but we determined that populations in the Xingu, Tapajós, and Araguaia+Paraná Rivers are likely undescribed species, while S. jurupari and S. mapiritensis did not show clear genetic distinction. To infer a phylogeny of these putative species, we conducted maximum likelihood and Bayesian non-clock and relaxed clock analyses of concatenated data from three genes (one mitochondrial, two nuclear). We also used a multi-species coalescent model to estimate a species tree from six unlinked loci (one mitochondrial, five nuclear). The topologies obtained were congruent with other results, but showed only minimal to moderate support for some nodes, suggesting that more loci will be needed to satisfactorily estimate the distribution of coalescent histories within Satanoperca. We determined that this variation results from topological discordance among separate gene trees, likely due to differential sorting of ancestral polymorphisms.

Genomic region associated with run timing has similar haplotypes and phenotypic effects across three lineages of Chinook salmon.

Conserving life-history variation is a stated goal of many management programs, but the most effective means by which to accomplish this are often far from clear. Early- and late-migrating forms of Chinook salmon (Oncorhynchus tshawytscha) face unequal pressure from natural and anthropogenic forces that may alter the impacts of genetic variation underlying heritable migration timing. Genomic regions of chromosome 28 are known to be strongly associated with migration variation in adult Chinook salmon, but it remains unclear whether there is consistent association among diverse lineages and populations in large basins such as the Columbia River. With high-throughput genotyping (GT-seq) and phenotyping methods, we examined the association of genetic variation in 28 markers (spanning GREB1L to ROCK1 of chromosome 28) with individual adult migration timing characteristics gleaned from passive integrated transponder recordings of over 5000 Chinook salmon from the three major phylogeographic lineages that inhabit the Columbia River Basin. Despite the strong genetic differences among them in putatively neutral genomic regions, each of the three lineages exhibited very similar genetic variants in the chromosome 28 region that were significantly associated with adult migration timing phenotypes. This is particularly notable for the interior stream-type lineage, which exhibits an earlier and more constrained freshwater entry than the other lineages. In both interior stream-type and interior ocean-type lineages of Chinook salmon, heterozygotes of the most strongly associated linkage groups had largely intermediate migration timing relative to homozygotes, and results indicate codominance or possibly marginal partial dominance of the allele associated with early migration. Our results lend support to utilization of chromosome 28 variation in tracking and predicting run timing in these lineages of Chinook salmon in the Columbia River.

Genotyping single nucleotide polymorphisms and inferring ploidy by amplicon sequencing for polyploid, ploidy-variable organisms.

Whole genome duplication is hypothesized to have played a critical role in the evolution of several major taxa, including vertebrates, and while many lineages have rediploidized, some retain polyploid genomes. Additionally, variation in ploidy can occur naturally or be artificially induced within select plant and animal species. Modern genetic techniques have not been widely applied to polyploid or ploidy-variable species, in part due to the difficulty of obtaining genotype data from polyploids. In this study, we demonstrate a strategy for developing an amplicon sequencing panel of single nucleotide polymorphisms for high-throughput genotyping of polyploid organisms. We then develop a method to infer ploidy of individuals from amplicon sequencing data that is generalized to apply to any ploidy and does not require prior identification of heterozygous genotypes. Combining these two techniques will allow researchers to both infer ploidy and generate ploidy-aware genotypes with the same amplicon sequencing panel. We demonstrate this approach with white sturgeon Acipenser transmontanus, a ploidy-variable (octoploid, decaploid and dodecaploid) imperiled species under conservation management in the Pacific Northwest and obtained a panel of 325 loci. These loci were validated by examining inheritance in known-cross families, and the ploidy inference method was validated with known ploidy samples. We provide scripts that adapt existing pipelines to genotype polyploids and an R package for application of the ploidy inference method. We expect that these techniques will empower studies of genetic variation and inheritance in polyploid organisms that vary in ploidy level, either naturally or as a result of artificial propagation practices.

Steelhead (Oncorhynchus mykiss) lineages and sexes show variable patterns of association of adult migration timing and age-at-maturity traits with two genomic regions.

As life history diversity plays a critical role in supporting the resilience of exploited populations, understanding the genetic basis of those life history variations is important for conservation management. However, effective application requires a robust understanding of the strength and universality of genetic associations. Here, we examine genetic variation of single nucleotide polymorphisms in genomic regions previously associated with migration phenology and age-at-maturity in steelhead (Oncorhynchus mykiss) from the Columbia River. We found chromosome 28 markers (GREB1L, ROCK1 genes) explained significant variance in migration timing in both coastal and inland steelhead. However, strength of association was much greater in coastal than inland steelhead (R 2 0.51 vs. 0.08), suggesting that genomic background and challenging inland migration pathways may act to moderate effects of this region. Further, we found that chromosome 25 candidate markers (SIX6 gene) were significantly associated with age and size at first return migration for inland steelhead, and this pattern was mediated by sex in a predictable pattern (males R 2 = 0.139-0.170; females R 2 = 0.096-0.111). While this encourages using these candidate regions in predicting life history characteristics, we suggest that stock-specific associations and haplotype frequencies will be useful in guiding implementation of genetic assays to inform management.

The largest fish in the world's biggest river: Genetic connectivity and conservation of Arapaima gigas in the Amazon and Araguaia-Tocantins drainages.

Arapaima, pirarucu or paiche (Arapaima gigas) is one of the largest freshwater fish in the world, and has a long history of commercial exploitation in the Amazon region. To estimate levels of genetic variability and historical and recent connectivity in Arapaima, we examined variation in eleven microsatellite DNA markers in individuals from 22 localities in Brazil, Colombia, and Peru. The results of analysis of molecular variance, Bayesian clustering and discriminant analysis of principal components showed that Arapaima in our samples represents two major populations, one in the Amazonas and one in the Araguaia-Tocantins River basins. The Amazonas population is further structured by isolation-by-distance with the hydrologically largely unconnected Amapá locality representing the eastern-most extreme of this continuum; gene flow predominates at distances of less than 1500 km with localities separated by over 2000 km dominated by genetic drift and effectively forming different populations. We saw no evidence of multiple species of Arapaima in the Amazonas basin, and analysis of pairwise genetic divergence (FST) with Mantel tests and correlograms indicated that this largest population exhibits a large-scale pattern of isolation-by-distance, with which results from MIGRATE-N agreed. The degree and significance of genetic divergence indicates that most sampled localities represent demographically independent sub-populations, although we did identify several recent migration events between both proximal and more distant localities. The levels of genetic diversity were heterogeneous across sites, including low genetic diversity, effective population sizes, and evidence of genetic bottlenecks in several places. On average the levels of gene diversity and rarefied allelic richness were higher for localities along the Amazonas mainstem than in the tributaries, despite these being the areas of highest fishing pressure, while the lowest values were found in tributary headwaters, where landscape modification is a significant threat. We recommend that managers consider the regional and local threats to these populations and tailor strategies accordingly, strategies which should ensure the ability of young A. gigas to disperse through floodplain corridors to maintain genetic diversity among otherwise sedentary adult sub-populations.

One species or four? Yes!...and, no. Or, arbitrary assignment of lineages to species obscures the diversification processes of Neotropical fishes.

Species are fundamental units in many biological disciplines, but there is continuing disagreement as to what species are, how to define them, and even whether the concept is useful. While some of this debate can be attributed to inadequate data and insufficient statistical frameworks in alpha taxonomy, an equal part results from the ambiguity over what species are expected to represent by the many who use them. Here, mtDNA data, microsatellite data, and sequence data from 17 nuclear loci are used in an integrated and quantitative manner to resolve the presence of evolutionary lineages, their contemporary and historical structure, and their correspondence to species, in a species complex of Amazonian peacock "bass" cichlids (Cichla pinima sensu lato). Results suggest that the historical narrative for these populations is more complex than can be portrayed by recognizing them as one, two, or four species: their history and contemporary dynamics cannot be unambiguously rendered as discrete units (taxa) at any level without both choosing the supremacy of one delimitation criterion and obscuring the very information that provides insight into the diversification process. This calls into question the utility of species as a rank, term, or concept, and suggests that while biologists may have a reasonable grasp of the structure of evolution, our methods of conveying these insights need updating. The lack of correspondence between evolutionary phenomena and discrete species should serve as a null hypothesis, and researchers should focus on quantifying the diversity in nature at whatever hierarchical level it occurs.

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.

Testing mitochondrial capture and deep coalescence in Amazonian cichlid fishes (Cichlidae: Cichla).

Hybridization and introgression have important consequences in evolution, such as increasing the genetic diversity and adaptive potential of a species. One of their most conspicuous footprints is discordance among gene trees or between genes and phenotypes. However, most studies that report introgression fail to disprove the null hypothesis that genetic incongruence may result from stochastic sorting of ancestral allelic polymorphisms. In the case of ancient introgression, these two processes may be especially difficult to distinguish topologically, but they make different predictions about the patterns of coalescence among loci. Here we apply three methods, molecular dating, multispecies coalescent models, and gene tree simulation under coalescence, to compare these two hypotheses that explain the polyphyletic mtDNA of the butterfly peacock bass, Cichla orinocensis. In comparison with a species tree based on 20 unlinked nuclear loci, we determined that mtDNA divergences were too recent to be explained by ancestral polymorphism. Similarly, coalescent species tree branches were significantly shorter when putative introgressed mtDNA was incorporated, and simulations showed the mtDNA topology to be unlikely under lineage sorting only. We conclude that introgression approximately 1.5 million years ago resulted in capture by C. orinocensis of an mtDNA lineage ancestral to the modern subspecies C. oc. monoculus.