Generating single-copy nuclear gene data for a recent adaptive radiation

Abstract

Recent adaptive radiations provide an exceptional opportunity to understand the processes of speciation and adaptation. However, reconstructing the phylogenetic history of recent and rapidly evolving clades often requires the use of multiple, independent gene genealogies. Nuclear introns are an obvious source of the necessary data but their use is often limited because degenerate primers can amplify paralogous loci. To identify PCR primers for a large number of loci in an especially rapid adaptive radiation, that of the flowering plant genus Aquilegia (Ranunculaceae), we developed an efficient method for amplifying multiple single-copy nuclear loci by sequencing a modest number of clones from a cDNA library and designing PCR primers; with one primer anchored in the 3′ untranslated region (3′-UTR) and one primer in the coding region of each gene. Variation between paralogous loci evolves more quickly in 3′-UTR regions compared to adjacent exons, and therefore we achieved high specificity for isolating orthologous loci. Furthermore, we were able to identify genes containing large introns by amplifying genes from genomic DNA and comparing the PCR product size to that predicted from their cDNA sequence. In Aquilegia eight out of eleven loci were isolated with this method and six of these loci had introns. Among four genes sequenced for samples spanning the phylogenetic breadth of the genus, we found sequence variation at levels similar to that observed in ITS, further supporting the recent and rapid radiation in Aquilegia. We assessed the orthology of amplification products by phylogenetic congruence among loci, the presence of two well established phylogenetic relationships, and similarity among loci for levels of sequence variation. Higher levels of variation among samples for one locus suggest possible paralogy. Overall, this method provides an efficient means of isolating predominantly single-copy loci from both low and high-copy gene families, providing ample nuclear variation for reconstructing species-level phylogenies in non-model taxa.