Paul D. Blischak

Accounting for genotype uncertainty in the estimation of allele frequencies in autopolyploids

Despite the increasing opportunity to collect large‐scale data sets for population genomic analyses, the use of high‐throughput sequencing to study populations of polyploids has seen little application. This is due in large part to problems associated with determining allele copy number in the genotypes of polyploid individuals (allelic dosage uncertainty–ADU), which complicates the calculation of important quantities such as allele frequencies. Here, we describe a statistical model to estimate biallelic SNP frequencies in a population of autopolyploids using high‐throughput sequencing data in the form of read counts. We bridge the gap from data collection (using restriction enzyme based techniques [e.g. GBS, RADseq]) to allele frequency estimation in a unified inferential framework using a hierarchical Bayesian model to sum over genotype uncertainty. Simulated data sets were generated under various conditions for tetraploid, hexaploid and octoploid populations to evaluate the models performance and to help guide the collection of empirical data. We also provide an implementation of our model in the R package polyfreqs and demonstrate its use with two example analyses that investigate (i) levels of expected and observed heterozygosity and (ii) model adequacy. Our simulations show that the number of individuals sampled from a population has a greater impact on estimation error than sequencing coverage. The example analyses also show that our model and software can be used to make inferences beyond the estimation of allele frequencies for autopolyploids by providing assessments of model adequacy and estimates of heterozygosity.

Gene Prediction and Annotation in Penstemon (Plantaginaceae): A Workflow for Marker Development from Extremely Low-Coverage Genome Sequencing

Premise of the study: Penstemon (Plantaginaceae) is a large and diverse genus endemic to North America. However, determining the phylogenetic relationships among its 280 species has been difficult due to its recent evolutionary radiation. The development of a large, multilocus data set can help to resolve this challenge. Methods: Using both previously sequenced genomic libraries and our own low-coverage whole-genome shotgun sequencing libraries, we used the MAKER2 Annotation Pipeline to identify gene regions for the development of sequencing loci from six extremely low-coverage Penstemon genomes (∼0.005×−0.007×). We also compared this approach to BLAST searches, and conducted analyses to characterize sequence divergence across the species sequenced. Results: Annotations and gene predictions were successfully added to more than 10,000 contigs for potential use in downstream primer design. Primers were then designed for chloroplast, mitochondrial, and nuclear loci from these annotated sequences. MAKER2 identified longer gene regions in all six Penstemon genomes when compared with BLASTN and BLASTX searches. The average level of sequence divergence among the six species was 7.14%. Discussion: Combining bioinformatics tools into a workflow that produces annotations can be useful for creating potential phylogenetic markers from thousands of sequences even when genome coverage is extremely low and reference data are only available from distant relatives. Furthermore, the output from MAKER2 contains information about important gene features, such as exon boundaries, and can be easily integrated with visualization tools to facilitate the process of marker development.

HyDe: A Python Package for Genome-Scale Hybridization Detection

Abstract.— The analysis of hybridization and gene flow among closely related taxa is a common goal for researchers studying speciation and phylogeography. Many methods for hybridization detection use simple site pattern frequencies from observed genomic data and compare them to null models that predict an absence of gene flow. The theory underlying the detection of hybridization using these site pattern probabilities exploits the relationship between the coalescent process for gene trees within population trees and the process of mutation along the branches of the gene trees. For certain models, site patterns are predicted to occur in equal frequency (i.e., their difference is 0), producing a set of functions called phylogenetic invariants. In this article, we introduce HyDe, a software package for detecting hybridization using phylogenetic invariants arising under the coalescent model with hybridization. HyDe is written in Python and can be used interactively or through the command line using pre‐packaged scripts. We demonstrate the use of HyDe on simulated data, as well as on two empirical data sets from the literature. We focus in particular on identifying individual hybrids within population samples and on distinguishing between hybrid speciation and gene flow. HyDe is freely available as an open source Python package under the GNU GPL v3 on both GitHub (https://github.com/pblischak/HyDe) and the Python Package Index (PyPI: https://pypi.python.org/pypi/phyde).

SNP genotyping and parameter estimation in polyploids using low-coverage sequencing data

Motivation Genotyping and parameter estimation using high throughput sequencing data are everyday tasks for population geneticists, but methods developed for diploids are typically not applicable to polyploid taxa. This is due to their duplicated chromosomes, as well as the complex patterns of allelic exchange that often accompany whole genome duplication (WGD) events. For WGDs within a single lineage (autopolyploids), inbreeding can result from mixed mating and/or double reduction. For WGDs that involve hybridization (allopolyploids), alleles are typically inherited through independently segregating subgenomes. Results We present two new models for estimating genotypes and population genetic parameters from genotype likelihoods for auto‐ and allopolyploids. We then use simulations to compare these models to existing approaches at varying depths of sequencing coverage and ploidy levels. These simulations show that our models typically have lower levels of estimation error for genotype and parameter estimates, especially when sequencing coverage is low. Finally, we also apply these models to two empirical datasets from the literature. Overall, we show that the use of genotype likelihoods to model non‐standard inheritance patterns is a promising approach for conducting population genomic inferences in polyploids. Availability and implementation A C ++ program, EBG, is provided to perform inference using the models we describe. It is available under the GNU GPLv3 on GitHub: https://github.com/pblischak/polyploid‐genotyping.

Population genetics of Penstemon albomarginatus (Plantaginaceae), a rare Mojave Desert species of conservation concern

Penstemon albomarginatus is a psammophytic endemic of the Mojave Desert, found only in deep sand and dune habitats of San Bernardino County, California, Mohave County, Arizona, and Clark and Nye Counties in Nevada. We used six microsatellite loci to assess genetic differentiation and diversity for 228 individuals across the 12 known populations of this rare species. A slight heterozygote deficiency was found in two populations, but most populations show no signs of inbreeding. Results show a geographic pattern of northern populations being more closely related to one another compared to all other geographic regions. Genetic diversity was greatest in the southern populations, with decreasing amounts of diversity observed with latitude. In general, the geographic pattern of genetic diversity among all populations suggests a post-glacial dispersal from south-to-north. Our results are discussed in the framework of anthropogenic pressures on deep sand habitats of the Mojave Desert.

Primers for Castilleja and their Utility Across Orobanchaceae: II. Single-Copy Nuclear Loci

Premise of the study: We developed primers targeting nuclear loci in Castilleja with the goal of reconstructing the evolutionary history of this challenging clade. These primers were tested across other major clades in Orobanchaceae to assess their broader utility. Methods and Results: We assembled low-coverage genomes for three taxa in Castilleja and developed primer combinations for the single-copy conserved ortholog set (COSII) and the pentatricopeptide repeat (PPR) gene family. These primer combinations were designed to take advantage of the Fluidigm microfluidic PCR platform and are well suited for high-throughput sequencing applications. Eighty-seven primers were designed for Castilleja, and 27 were found to have broader utility in Orobanchaceae. Conclusions: These results demonstrate the utility of these primers, not only across Castilleja, but for other lineages within Orobanchaceae as well. This expanded molecular toolkit will be an asset to future phylogenetic studies in Castilleja and throughout Orobanchaceae.

Genomic signatures of sympatric speciation with historical and contemporary gene flow in a tropical anthozoan

Sympatric diversification is increasingly thought to have played an important role in the evolution of biodiversity around the globe. However, an in situ sympatric origin for co-distributed taxa is difficult to demonstrate empirically because different evolutionary processes can lead to similar biogeographic outcomes-especially in ecosystems with few hard barriers to dispersal that can facilitate allopatric speciation followed by secondary contact (e.g. marine habitats). Here we use a genomic (ddRADseq), model-based approach to delimit a cryptic species complex of tropical sea anemones that are co-distributed on coral reefs throughout the Tropical Western Atlantic. We use coalescent simulations in fastsimcoal2 to test competing diversification scenarios that span the allopatric-sympatric continuum. We recover support that the corkscrew sea anemone Bartholomea annulata (Le Sueur, 1817) is a cryptic species complex, co-distributed throughout its range. Simulation and model selection analyses suggest these lineages arose in the face of historical and contemporary gene flow, supporting a sympatric origin, but an alternative secondary contact model also receives appreciable model support. Leveraging the genome of Exaiptasia pallida we identify five loci under divergent selection between cryptic B. annulata lineages that fall within mRNA transcripts or CDS regions. Our study provides a rare empirical, genomic example of sympatric speciation in a tropical anthozoan-a group that includes reef-building corals. Finally, these data represent the first range-wide molecular study of any tropical sea anemone, underscoring that anemone diversity is under described in the tropics, and highlighting the need for additional systematic studies into these ecologically and economically important species.

Fluidigm2PURC: Automated processing and haplotype inference for double-barcoded PCR amplicons

Premise of the Study Targeted enrichment strategies for phylogenomic inference are a time‐ and cost‐efficient way to collect DNA sequence data for large numbers of individuals at multiple, independent loci. Automated and reproducible processing of these data is a crucial step for researchers conducting phylogenetic studies. Methods and Results We present Fluidigm2PURC, an open source Python utility for processing paired‐end Illumina data from double‐barcoded PCR amplicons. In combination with the program PURC (Pipeline for Untangling Reticulate Complexes), our scripts process raw FASTQ files for analysis with PURC and use its output to infer haplotypes for diploids, polyploids, and samples with unknown ploidy. We demonstrate the use of the pipeline with an example data set from the genus Thalictrum (Ranunculaceae). Conclusions Fluidigm2PURC is freely available for Unix‐like operating systems on GitHub (https://github.com/pblischak/fluidigm2purc) and for all operating systems through Docker (https://hub.docker.com/r/pblischak/fluidigm2purc).