Although the plastome has long been considered the workhorse of phylogenetic inference in plants, reliance on chloroplast data alone may limit the ability to identify inheritance patterns of polyploids, as well as introgression and hybridization events (Godden et al., 2012; Twyford and Ennos, 2012). The reliable estimation of an underlying species tree also depends upon the acquisition of multiple, unlinked loci—especially for recent and rapid radiations—shifting the focus toward the development of single- or low-copy nuclear gene regions for phylogenetic analyses. The process of identifying and developing these nuclear gene regions using traditional methods can be time consuming and costly, but the increasing availability of high-throughput sequencing data, as well as new bioinformatic approaches, allows for the efficient and cost-effective exploration of the nuclear genome.
Here we focus on developing a suite of putatively single-copy nuclear gene regions in Castilleja L. (Orobanchaceae; “the paintbrushes”), a clade rich with polyploid and hybrid taxa, and the product of an ongoing rapid radiation (Tank and Olmstead, 2008, 2009). Previous studies using the nuclear ribosomal ITS and ETS regions, the low-copy nuclear gene waxy, and the plastid trnL-F and rps16 intron regions hinted at cytonuclear discordance in some taxa, and most relationships among closely related taxa were unresolved (Tank and Olmstead, 2008, 2009). We recently developed primer combinations targeting the most variable regions of the plastome in Castilleja (Latvis et al., 2017a), and now present a companion set of nuclear primers with the goal of obtaining a resolved species tree for this challenging clade, as well as to aid in the detection of introgression and hybrid speciation. Primers for microsatellite markers have also been developed by Fant et al. (2013) for population-level investigations. In part, we follow the approach outlined by Blischak et al. (2014) to develop nuclear primer combinations from genome-skimming data, while following specifications for the Fluidigm Access Array microfluidic PCR system (Fluidigm, South San Francisco, California, USA) (see Latvis et al., 2017a). Thus, all primer combinations use the same annealing temperature of 60°C and may be amplified in parallel prior to high-throughput sequencing or traditional Sanger sequencing. We specifically target putatively single-copy genes from the conserved ortholog set (COSII) and pentatricopeptide repeat (PPR) domains, both of which have been highlighted for their phylogenetic utility in plants (COSII: Wu et al., 2006; PPR domain: Li et al., 2008; Yuan et al., 2009, 2010). We also test these primers for their broader applicability across Orobanchaceae following the approach outlined in Latvis et al. (2017a) with the goal of finding a subset of nuclear gene regions that would amplify across this more inclusive clade. Previous phylogenetic studies within Orobanchaceae employed the nuclear phytochrome genes PHYA and PHYB, the nuclear ribosomal ITS region, and the plastid matK and rps2 genes. Orobanchaceae is the largest clade of parasitic angiosperms, and plastome reduction and accelerated rates of molecular evolution in retained plastid genes have been documented (see discussion in Bennett and Mathews, 2006). Additionally, phytochrome genes regulate responses to light and can be significantly modified in parasites (Bennett and Mathews, 2006). Therefore, the development of additional single-copy nuclear regions would provide a much-needed source of phylogenetic information in Orobanchaceae, and may provide a more reliable estimate of branch lengths for further studies of diversification and character evolution.
Nuclear primer pairs designed for Castilleja (locus and region amplified), amplicon lengths, and validation results for Orobanchaceae and outgroup taxon Paulownia. All pairs were designed for an annealing temperature of 60°C (±1°C). Boldfaced rows correspond to core Orobanchaceae primers, defined by successful amplification in two or more major clades in Orobanchaceae (see Fig. 1).
METHODS AND RESULTS
We assembled contigs from raw reads of three low-coverage Castilleja genomes, C. cusickii Greenm., C.foliolosa Hook. & Arn., and C. tenuis (A. Heller) T. I. Chuang & Heckard (Latvis et al., 2017a; National Center for Biotechnology Information [NCBI] Sequence Read Archive [SRA] accession SRP100222) using CAP3 (Huang and Madan, 1999) with the default settings. The accessions were sequenced as 100-bp single-end reads on an Illumina HiSeq 2000 (Illumina, San Diego, California, USA), yielding ~12.5 million reads per taxon (Uribe-Convers et al., 2014) and an average depth of coverage of ~0.8×. These taxa include both annual and perennial lineages of Castilleja and span the phylogenetic breadth of the clade (Tank and Olmstead, 2008, 2009). These assemblies were then culled to include only contigs of 1 Kb or larger using a custom R script. The culled assemblies and script are available from the Dryad Digital Repository ( http://doi.org/10.5061/dryad.52v62; Latvis et al., 2017b).
To search for hits among our contigs, available COS sequences were obtained from Sol genomics ( https://solgenomics.net), and PPR loci were mined from the Mimulus L. genome on Phytozome (Hellsten et al., 2013; https://phytozome.jgi.doe.gov) using the 127 PPR orthologs identified in Arabidopsis Heynh. by Yuan et al. (2009) as references. Both gene sets may be found in Uribe-Convers et al. (2016) and were used to construct local BLAST databases for the search (-makeblastdb). We used TBLASTX to search each Castilleja CAP3 assembly (with contigs of 1 KB or greater) against both the COS and PPR databases, indicating tab-delimited output (-outfmt 6). Output files were filtered for alignment length >200 and a maximum E-value of 1e-10, and were culled to include only unique hits.
Hits shared between C. cusickii, C. foliolosa, and C. tenuis were placed together into individual FASTA files (data available from the Dryad repository: http://doi.org/10.5061/dryad.52v62; Latvis et al., 2017b), imported into Geneious R7 version 7.0.6 (Kearse et al., 2012), and aligned with MAFFT version 7.017b under the default settings (Katoh and Standley, 2013). We designed primer pairs using Primer3 (Untergasser et al., 2012) using the same specifications for the Fluidigm Access Array system as Latvis et al. (2017a), but with a size range between 400–525 bp and an optimal size of 500 bp. We designed 10–30 primer pairs for each identified locus and prioritized them based on desired size and the presence of multiple G or C bases at the 3′ end of the primers (GC clamp). This also allowed us to design overlapping sets of primers with the potential to combine them after sequencing to produce longer contigs for downstream phylogenetic analyses. Suitable primer pairs were validated for Castilleja with PCR following the same amplification protocol and using the same Castilleja accessions as in Latvis et al. (2017a) and visualized on an agarose gel. We present 87 nuclear primer combinations specifically designed and validated for Castilleja (Table 1).
To investigate whether any of these primer combinations would amplify successfully across Orobanchaceae, we searched for our selected Castilleja primers against an assembled low-coverage genome for Lamourouxia multifida Kunth using BLAST (Altschul et al., 1990). Lamourouxia multifida was sequenced on an Illumina HiSeq 2000 at the University of Oregon as 100-bp paired-end reads, and contigs were assembled using SPAdes (Bankevich et al., 2012) under the default settings. BLAST search parameters, assessment of suitable hits, and subsequent PCR validation with Lamourouxia virgata Kunth, Physocalyx major Mart., and Neobartsia filiformis (Wedd.) Uribe-Convers & Tank are described in Latvis et al. (2017a). Primer combinations with successful amplification in Lamourouxia Kunth and at least one other taxon were selected for further PCR testing with other major lineages in Orobanchaceae (sensu McNeal et al., 2013; Fig. 1). This second round of PCR validation follows Latvis et al. (2017a), except that two of the accessions used for testing were changed. As in Latvis et al. (2017a), we also included a negative control and conserved sequence-tagged “universal” primers for the trnL-F region as a positive control for all primer pairs. We used Neobartsia peruviana (Walp.) Uribe-Convers & Tank instead of N. filiformis, and Paulownia fortunei (Seem.) Hemsl. instead of P. elongata Siebold & Zucc. (Appendix 1). Of the 87 nuclear primer combinations specifically designed for Castilleja, we identified 27 with broader applicability in Orobanchaceae, chosen if they successfully amplified in Pedicularideae (Clade IV; including Castilleja, Lamourouxia, and Pedicularis L.) and at least one of the other major clades highlighted in Fig. 1. Validation results are presented in Table 1 with these “core Orobanchaceae” combinations boldfaced.
We present 87 nuclear primer pairs specifically designed for Castilleja that target COSII and PPR loci. Although we target the same putative single-copy nuclear domains as previous studies (Wu et al., 2006; Li et al., 2008; Yuan et al., 2010; Blischak et al., 2014; Uribe-Convers et al., 2016), we developed primers for different loci and present unique primer combinations in this study. As with our chloroplast primers (Latvis et al., 2017a), all combinations were designed with the Fluidigm microfluidic PCR system in mind, allowing for parallelization of amplification for downstream high-throughput sequencing platforms. Of these, we identify a set of 27 primer combinations with broader utility across Orobanchaceae. The development of primers for putatively single-copy nuclear loci will greatly enhance efforts to understand evolutionary history at multiple taxonomic scales, both for Castilleja and across Orobanchaceae.
This research was supported by resources at the Institute for Bioinformatics and Evolutionary Studies (IBEST; NIH/NCRR P20RR16448 and P20RR016454) and by the following awards from the National Science Foundation: DEB-1253463 (awarded to D.C.T.), DEB-1502061 (awarded to D.C.T. for S.J.J.), and DEB-1455399 (support for P.D.B.).
- Altschul, S. F., W. Gish, W. Miller, E. W. Meyer, and D. J. Lipman. 1990. Basic local alignment search tool. Journal of Molecular Biology 215: 403–410. Google Scholar
- Bankevich, A., S. Nurk, D. Antipov, A. Gurevich, M. Dvorkin, A. S. Kulikov, V. Lesin, et al. 2012. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. Journal of Computational Biology 19: 455–477. Google Scholar
- Bennett, J. R., and S. Mathews. 2006. Phylogeny of the parasitic plant family Orobanchaceae inferred from phytochrome A. American Journal of Botany 93: 1039–1051. Google Scholar
- Blischak, P. D., A. J. Wenzel, and A. D. Wolfe. 2014. Gene prediction and annotation in Penstemon (Plantaginaceae): A workflow for marker development from extremely low-coverage genome sequencing. Applications in Plant Sciences 2: 1400044. Google Scholar
- Fant, J. B., H. Wolf-Weinberg, D. C. Tank, and K. A. Skogen. 2013. Characterization of microsatellite loci in Castilleja sessiliflora and transferability to 24 Castilleja species (Orobanchaceae). Applications in Plant Sciences 1: 1200564. Google Scholar
- Godden, G. T., I. E. Jordan-Thaden, S. Chamala, A. A. Growl, N. Garcia, C. C. Germain-Aubrey, J. M. Heaney, et al. 2012. Making nextgeneration sequencing work for you: Approaches and practical considerations for marker development and phylogenetics. Plant Ecology and Diversity 5: 427–450. Google Scholar
- Hellsten, U., K. M. Wright, J. Jenkins, S. Shu, Y. Yuan, S. R. Wessler, J. Schmutz, et al. 2013. Fine-scale variation in meiotic recombination in Mimulus inferred from population shotgun sequencing. Proceedings of the National Academy of Sciences USA. Google Scholar
- Huang, X., and A. Madan. 1999. CAP3: A DNA sequence assembly program. Genome Research 9: 868–877. Google Scholar
- Katoh, K., and D. M. Standley. 2013. MAFFT Multiple Sequence Alignment Software version 7: Improvements in performance and usability. Molecular Biology and Evolution 30: 772–780. Google Scholar
- Kearse, M., R. Moir, A. Wilson, S. Stones-Havas, M. Cheung, S. Sturrock, S. Buxton, et al. 2012. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics (Oxford, England) 28: 1647–1649. Google Scholar
- Latvis, M., S. M. E. Mortimer, D. F. Morales-Briones, S. Torpey, S. Uribe-Convers, S. J. Jacobs, S. Mathews, and D. C. Tank. 2017a. Primers for Castilleja and their utility across Orobanchaceae: I. Chloroplast primers. Applications in Plant Sciences 5: 1700020. Google Scholar
- Li, M., J. Wunder, G. Bissoli, E. Scarponi, S. Gazzani, E. Barbaro, H. Saedler, and C. Varotto. 2008. Development of COS genes as universally applicable markers for phylogenetic reconstructions of closely related plant species. Cladistics 24: 727–745. Google Scholar
- McNeal, J. R., J. R. Bennett, A. D. Wolfe, and S. Mathews. 2013. Phylogeny and origins of holoparasitism in Orobanchaceae. American Journal of Botany 100: 971–983. Google Scholar
- Tank, D. C., and R. G. Olmstead. 2008. From annuals to perennials: Phylogeny of subtribe Castillejinae (Orobanchaceae). American Journal of Botany 95: 608–625. Google Scholar
- Tank, D. C., and R. G. Olmstead. 2009. The evolutionary origin of a second radiation of annual Castillejo (Orobanchaceae) species in South America: The role of long distance dispersal and allopolyploidy. American Journal of Botany 96: 1907–1921. Google Scholar
- Twyford, A. D., and R. A. Ennos. 2012. Next-generation hybridization and introgression. Heredity 108: 179–189. Google Scholar
- Untergasser, A., I. Cutcutache, T. Koressaar, J. Ye, B. C. Faircloth, M. Remm, and S. G. Rozen. 2012. Primer3–New capabilities and interfaces. Nucleic Acids Research 40: e115. Google Scholar
- Uribe-Convers, S., J. R. Duke, M. J. Moore, and D. C. Tank. 2014. A long-PCR based method for chloroplast genome enrichment and phylogenomics in angiosperms. Applications in Plant Sciences 2: 1300063. Google Scholar
- Uribe-Convers, S., M. L. Settles, and D. C. Tank. 2016. A phylogenomic approach based on PCR target enrichment and high throughput sequencing: Resolving the diversity within the South American species of Bartsia L. (Orobanchaceae). PLoS ONE 11: e0148203. Google Scholar
- Wu, F., L. A. Mueller, D. Crouzillat, V. Pétiard, and S. D. Tanksley. 2006. Combining bioinformatics and phylogenetics to identify large sets of single-copy orthologous genes (COSII) for comparative, evolutionary and systematic studies: A test case in the Euasterid plant clade. Genetics 174: 1407–1420. Google Scholar
- Yuan, Y.-W., C. Liu, H. E. Marx, and R. G. Olmstead. 2009. The pentatricopeptide repeat (PPR) gene family, a tremendous resource for plant phylogenetic studies. New Phytologist 182: 272–283. Google Scholar
- Yuan, Y., C. Liu, H. E. Marx, and R. G. Olmstead. 2010. An empirical demonstration of using pentatricopeptide repeat (PPR) genes as plant phylogenetic tools: Phylogeny of Verbenaceae and the Verbena complex. Molecular Phylogenetics and Evolution 54: 23–35. Google Scholar