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.
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. 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. 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.
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.
Display omitted
•Sphoeroides (Tetraodontidae) is recovered as paraphyletic.•Colomesus should be considered a junior synonym of Sphoeroides.•Sphoeroides pachygaster may harbor cryptic ...diversity.•Analyses suggest colonization of the New World through the eastern Pacific.
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.
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).
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.
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.
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.
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.
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.
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.
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 (R2 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 R2 = 0.139–0.170; females R2 = 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.