Thomas L. Parchman

Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery

BackgroundMassively parallel sequencing of cDNA is now an efficient route for generating enormous sequence collections that represent expressed genes. This approach provides a valuable starting point for characterizing functional genetic variation in non-model organisms, especially where whole genome sequencing efforts are currently cost and time prohibitive. The large and complex genomes of pines (Pinus spp.) have hindered the development of genomic resources, despite the ecological and economical importance of the group. While most genomic studies have focused on a single species (P. taeda), genomic level resources for other pines are insufficiently developed to facilitate ecological genomic research. Lodgepole pine (P. contorta) is an ecologically important foundation species of montane forest ecosystems and exhibits substantial adaptive variation across its range in western North America. Here we describe a sequencing study of expressed genes from P. contorta, including their assembly and annotation, and their potential for molecular marker development to support population and association genetic studies.ResultsWe obtained 586,732 sequencing reads from a 454 GS XLR70 Titanium pyrosequencer (mean length: 306 base pairs). A combination of reference-based and de novo assemblies yielded 63,657 contigs, with 239,793 reads remaining as singletons. Based on sequence similarity with known proteins, these sequences represent approximately 17,000 unique genes, many of which are well covered by contig sequences. This sequence collection also included a surprisingly large number of retrotransposon sequences, suggesting that they are highly transcriptionally active in the tissues we sampled. We located and characterized thousands of simple sequence repeats and single nucleotide polymorphisms as potential molecular markers in our assembled and annotated sequences. High quality PCR primers were designed for a substantial number of the SSR loci, and a large number of these were amplified successfully in initial screening.ConclusionsThis sequence collection represents a major genomic resource for P. contorta, and the large number of genetic markers characterized should contribute to future research in this and other pines. Our results illustrate the utility of next generation sequencing as a basis for marker development and population genomics in non-model species.

Stick insect genomes reveal natural selection's role in parallel speciation.

Stick to the Bush Can the underlying genetic changes driving the divergence of populations into new species be predicted or repeated? Soria-Carrasco et al. (p. 738) investigated the genetic changes observed after one generation when stick insect (Timema cristinae) populations were transplanted from their preferred host plants to alternative hosts. Diverged genetic regions were relatively small, with most loci showing divergence in a single population pair. However, the number of loci showing parallel divergence was greater than expected by chance. Thus, selection can drive parallel phenotypic evolution via parallel genetic changes. Parallel speciation in insects shows both convergent and divergent selection after one generation. Natural selection can drive the repeated evolution of reproductive isolation, but the genomic basis of parallel speciation remains poorly understood. We analyzed whole-genome divergence between replicate pairs of stick insect populations that are adapted to different host plants and undergoing parallel speciation. We found thousands of modest-sized genomic regions of accentuated divergence between populations, most of which are unique to individual population pairs. We also detected parallel genomic divergence across population pairs involving an excess of coding genes with specific molecular functions. Regions of parallel genomic divergence in nature exhibited exceptional allele frequency changes between hosts in a field transplant experiment. The results advance understanding of biological diversification by providing convergent observational and experimental evidence for selection’s role in driving repeatable genomic divergence.

Genome-wide association genetics of an adaptive trait in lodgepole pine

Pine cones that remain closed and retain seeds until fire causes the cones to open (cone serotiny) represent a key adaptive trait in a variety of pine species. In lodgepole pine, there is substantial geographical variation in serotiny across the Rocky Mountain region. This variation in serotiny has evolved as a result of geographically divergent selection, with consequences that extend to forest communities and ecosystems. An understanding of the genetic architecture of this trait is of interest owing to the wide‐reaching ecological consequences of serotiny and also because of the repeated evolution of the trait across the genus. Here, we present and utilize an inexpensive and time‐effective method for generating population genomic data. The method uses restriction enzymes and PCR amplification to generate a library of fragments that can be sequenced with a high level of multiplexing. We obtained data for more than 95 000 single nucleotide polymorphisms across 98 serotinous and nonserotinous lodgepole pines from three populations. We used a Bayesian generalized linear model (GLM) to test for an association between genotypic variation at these loci and serotiny. The probability of serotiny varied by genotype at 11 loci, and the association between genotype and serotiny at these loci was consistent in each of the three populations of pines. Genetic variation across these 11 loci explained 50% of the phenotypic variation in serotiny. Our results provide a first genome‐wide association map of serotiny in pines and demonstrate an inexpensive and efficient method for generating population genomic data.

Reciprocal selection causes a coevolutionary arms race between crossbills and lodgepole pine.

Few studies have shown both reciprocal selection and reciprocal adaptations for a coevolving system in the wild. The goal of our study was to determine whether the patterns of selection on Rocky Mountain lodgepole pine (Pinus contorta spp. latifolia) and red crossbills (Loxia curvirostra complex) were concordant with earlier published evidence of reciprocal adaptations in lodgepole pine and crossbills on isolated mountain ranges in the absence of red squirrels (Tamiasciurus hudsonicus). We found that selection (directional) by crossbills on lodgepole pine where Tamiasciurus are absent was divergent from the selection (directional) exerted by Tamiasciurus on lodgepole pine. This resulted in divergent selection between areas with and without Tamiasciurus that was congruent with the geographic patterns of cone variation. In the South Hills, Idaho, where Tamiasciurus are absent and red crossbills are thought to be coevolving with lodgepole pine, crossbills experienced stabilizing selection on bill size, with cone structure as the agent of selection. These results show that crossbills and lodgepole pine exhibit reciprocal adaptations in response to reciprocal selection, and they provide insight into the traits mediating and responding to selection in a coevolutionary arms race.

DIVERSIFYING COEVOLUTION BETWEEN CROSSBILLS AND BLACK SPRUCE ON NEWFOUNDLAND

Abstract Coevolution is increasingly recognized as an important process structuring geographic variation in the form of selection for many populations. Here we consider the importance of a geographic mosaic of coevolution to patterns of crossbill (Loxia) diversity in the northern boreal forests of North America. We examine the relationships between geographic variation in cone morphology, bill morphology, and feeding performance to test the hypothesis that, in the absence of red squirrels (Tamiasciurus hudsonicus), black spruce (Picea mariana) has lost seed defenses directed at Tamiasciurus and that red crossbills (L. curvirostra) and black spruce have coevolved in an evolutionary arms race. Comparisons of cone morphology and several indirect lines of evidence suggest that black spruce has evolved defenses in response to Tamiasciurus on mainland North America but has lost these defenses on Newfoundland. Cone traits that deter crossbills, including thicker scales that require larger forces to separate, are elevated in black spruce on Newfoundland, and larger billed crossbills have higher feeding performances than smaller billed crossbills on black spruce cones from Newfoundland. These results imply that the large bill of the Newfoundland crossbill (L. c. percna) evolved as an adaptation to the elevated cone defenses on Newfoundland and that crossbills and black spruce coevolved in an evolutionary arms race on Newfoundland during the last 9000 years since glaciers retreated. On the mainland where black spruce is not as well defended against crossbills, the small‐billed white‐winged crossbill (L. leucoptera leucoptera) is more efficient and specializes on seeds in the partially closed cones. Finally, reciprocal adaptations between crossbills and conifers are replicated in black spruce and Rocky Mountain lodgepole pine (Pinus contorta ssp. latifolia), with coevolution most pronounced in isolated populations where Tamiasciurus are absent as a competitor. This study further supports the role of Tamiasciurus in determining the selection mosaic for crossbills and suggests that a geographic mosaic of coevolution has been a prominent factor underlying the diversification of North American crossbills.

The genomic consequences of adaptive divergence and reproductive isolation between species of manakins

The processes of adaptation and speciation are expected to shape genomic variation within and between diverging species. Here we analyze genomic heterogeneity of genetic differentiation and introgression in a hybrid zone between two bird species (Manacus candei and M. vitellinus) using 59 100 SNPs, a whole genome assembly, and Bayesian models. Measures of genetic differentiation ( FST ) and introgression (genomic cline center [α] and rate [β]) were highly heterogeneous among loci. We identified thousands of loci with elevated parameter estimates, some of which are likely to be associated with variation in fitness in Manacus populations. To analyze the genomic organization of differentiation and introgression, we mapped SNPs onto a draft assembly of the M. vitellinus genome. Estimates of FST , α, and β were autocorrelated at very short physical distances (< 100 bp), but much less so beyond this. In addition, average statistical associations (linkage disequilibrium) between SNPs were generally low and were not higher in admixed populations than in populations of the parental species. Although they did not occur with a constant probability across the genome, loci with elevated FST , α, and β were not strongly co‐localized in the genome. Contrary to verbal models that predict clustering of loci involved in adaptation and isolation in discrete genomic regions, these results are consistent with the hypothesis that genetic regions involved in adaptive divergence and reproductive isolation are scattered throughout the genome. We also found that many loci were characterized by both exceptional genetic differentiation and introgression, consistent with the hypothesis that loci involved in isolation are also often characterized by a history of divergent selection. However, the concordance between isolation and differentiation was only partial, indicating a complex architecture and history of loci involved in isolation.

Genomic consequences of multiple speciation processes in a stick insect

Diverse geographical modes and mechanisms of speciation are known, and individual speciation genes have now been identified. Despite this progress, genome-wide outcomes of different evolutionary processes during speciation are less understood. Here, we integrate ecological and spatial information, mating trials, transplantation data and analysis of 86 130 single nucleotide polymorphisms (SNPs) in eight populations (28 pairwise comparisons) of Timema cristinae stick insects to test the effects of different factors on genomic divergence in a system undergoing ecological speciation. We find patterns consistent with effects of numerous factors, including geographical distance, gene flow, divergence in host plant use and climate, and selection against maladaptive hybridization (i.e. reinforcement). For example, the number of highly differentiated ‘outlier loci’, allele-frequency clines and the overall distribution of genomic differentiation were recognizably affected by these factors. Although host use has strong effects on phenotypic divergence and reproductive isolation, its effects on genomic divergence were subtler and other factors had pronounced effects. The results demonstrate how genomic data can provide new insights into speciation and how genomic divergence can be complex, yet predictable. Future work could adopt experimental, mapping and functional approaches to directly test which genetic regions are affected by selection and determine their physical location in the genome.

Genomics of isolation in hybrids

Hybrid zones are common in nature and can offer critical insights into the dynamics and components of reproductive isolation. Hybrids between diverged lineages are particularly informative about the genetic architecture of reproductive isolation, because introgression in an admixed population is a direct measure of isolation. In this paper, we combine simulations and a new statistical model to determine the extent to which different genetic architectures of isolation leave different signatures on genome-level patterns of introgression. We found that reproductive isolation caused by one or several loci of large effect caused greater heterogeneity in patterns of introgression than architectures involving many loci with small fitness effects, particularly when isolating factors were closely linked. The same conditions that led to heterogeneous introgression often resulted in a reasonable correspondence between outlier loci and the genetic loci that contributed to isolation. However, demographic conditions affected both of these results, highlighting potential limitations to the study of the speciation genomics. Further progress in understanding the genomics of speciation will require large-scale empirical studies of introgression in hybrid zones and model-based analyses, as well as more comprehensive modelling of the expected levels of isolation with different demographies and genetic architectures of isolation.

Experimental evidence for ecological selection on genome variation in the wild

Understanding natural selections effect on genetic variation is a major goal in biology, but the genome-scale consequences of contemporary selection are not well known. In a release and recapture field experiment we transplanted stick insects to native and novel host plants and directly measured allele frequency changes within a generation at 186 576 genetic loci. We observed substantial, genome-wide allele frequency changes during the experiment, most of which could be attributed to random mortality (genetic drift). However, we also documented that selection affected multiple genetic loci distributed across the genome, particularly in transplants to the novel host. Host-associated selection affecting the genome acted on both a known colour-pattern trait as well as other (unmeasured) phenotypes. We also found evidence that selection associated with elevation affected genome variation, although our experiment was not designed to test this. Our results illustrate how genomic data can identify previously underappreciated ecological sources and phenotypic targets of selection.

Patterns of genetic variation in the adaptive radiation of New World crossbills (Aves: Loxia )

Incipient species groups or young adaptive radiations such as crossbills (Aves: Loxia) present the opportunity to investigate directly the processes occurring during speciation. New World crossbills include white‐winged crossbills (Loxia leucoptera), Hispaniolan crossbills (Loxia megaplaga), and red crossbills (Loxia curvirostra complex), the last of which is comprised of at least nine morphologically and vocally differentiated forms (‘call types’) where divergent natural selection for specialization on different conifer resources has been strongly implicated as driving diversification. Here we use amplified fragment length polymorphism (AFLP) markers to investigate patterns of genetic variation across populations, call types, and species of New World crossbills. Tree‐based analyses using 440 AFLP loci reveal strongly supported clustering of the formally recognized species, but did not separate individuals from the eight call types in the red crossbill complex, consistent with recent divergence and ongoing gene flow. Analyses of genetic differentiation based on inferred allele frequency variation however, reveal subtle but significant levels of genetic differentiation among the different call types of the complex and indicate that between call‐type differentiation is greater than that found among different geographic locations within call types. Interpreted in light of evidence of divergent natural selection and strong premating reproductive isolation, the observed genetic differentiation suggests restricted gene flow among sympatric call types consistent with the early stages of ecological speciation.