Haploid, diploid, and pooled exome capture recapitulate features of biology and paralogy in two non-model tree species
Lind, Brandon et al. (2021), Haploid, diploid, and pooled exome capture recapitulate features of biology and paralogy in two non-model tree species, Dryad, Dataset, https://doi.org/10.5061/dryad.k0p2ngf7w
Despite their suitability for studying evolution, many conifer species have large and repetitive giga-genomes (16-31Gbp) that create hurdles to producing high coverage SNP datasets that capture diversity from across the entirety of the genome. Due in part to multiple ancient whole genome duplication events, gene family expansion and subsequent evolution within Pinaceae, false diversity from the misalignment of paralog copies creates further challenges in accurately and reproducibly inferring evolutionary history from sequence data. Here, we leverage the cost-saving benefits of pool-seq and exome-capture to discover SNPs in two conifer species, Douglas-fir (Pseudotsuga menziesii var. menziesii (Mirb.) Franco, Pinaceae) and jack pine (Pinus banksiana Lamb., Pinaceae). We show, using minimal baseline filtering, that allele frequencies estimated from pooled individuals show a strong positive correlation with those estimated by sequencing the same population as individuals (r > 0.948), on par with such comparisons made in model organisms. Further, we highlight the utility of haploid megagametophyte tissue for identifying sites that are likely due to misaligned paralogs. Together with additional minor filtering, we show that it is possible to remove many of the loci with large frequency estimate discrepancies between individual and pooled sequencing approaches, improving the correlation further (r > 0.973). Our work addresses bioinformatic challenges in non-model organisms with large and complex genomes, highlights the use of megagametophyte tissue for the identification of paralog sites, and suggests the combination of pool-seq and exome capture to be robust for further evolutionary hypothesis testing in these systems.
This data was collected from natural populations. Exome-capture pool-seq (20 diploid individuals), individually sequenced (the same diploid individuals in pools), and haploid megagametophyte data from a single individual was sequenced on an Illumina HiSeq4000 instrument at Centre d'expertise et de services Génome Québec, Montréal, Canada. Sequence data was processed according to bioinformatic best practices.
All code to analyze these files is also attached.
Each file of code is saved as jupyter notebook format (.ipynb) and as .html. HTML can be used to view the notebook without launching a jupyter kernel.