UBC Research Data

Data from: The population genomics of sunflowers and genomic determinants of protein evolution revealed by RNAseq Renaut, Sébastien; Grassa, Christopher J.; Moyers, Brook T.; Kane, Nolan C.; Rieseberg, Loren H.; Grassa, Christopher; Moyers, Brook; Kane, Nolan; Rieseberg, Loren


Few studies have investigated the causes of evolutionary rate variation among plant nuclear genes, especially in recently diverged species still capable of hybridizing in the wild. The recent advent of Next Generation Sequencing (NGS) permits investigation of genome wide rates of protein evolution and the role of selection in generating and maintaining divergence. Here, we use individual whole-transcriptome sequencing (RNAseq) to refine our understanding of the population genomics of a wild species of sunflowers (Helianthus spp.) and the factors that affect rates of protein evolution. We aligned 35 GB of transcriptome sequencing data and identified 433,257 polymorphic sites (SNPs) in a reference transcriptome comprising 16,312 genes. Using SNP markers, we identified strong population clustering largely corresponding to the three species analyzed here (Helianthus annuus, H. petiolaris, H. debilis), with one distinct early generation hybrid. Then, we calculated the proportions of adaptive substitution fixed by selection (alpha) and identified gene ontology categories with elevated values of alpha. The “response to biotic stimulus” category had the highest mean alpha across the three interspecific comparisons, implying that natural selection imposed by other organisms plays an important role in driving protein evolution in wild sunflowers. Finally, we examined the relationship between protein evolution (dN/dS ratio) and several genomic factors predicted to co-vary with protein evolution (gene expression level, divergence and specificity, genetic divergence [FST], and nucleotide diversity [pi]). We find that variation in rates of protein divergence was correlated with gene expression level and specificity, consistent with results from a broad range of taxa and timescales. This would in turn imply that these factors govern protein evolution both at a microevolutionary and macroevolutionary timescale. Our results contribute to a general understanding of the determinants of rates of protein evolution and the impact of selection on patterns of polymorphism and divergence.; Usage notes
reference transcriptomefasta file containing reference transcriptome (16312 contigs) described in the method section.TranscriptomeRef_mar30.fa
snp_tableTable containing genotypes for all 433,257 SNPs. First column contains contig name were the SNP is found, second column the reference allele, third column the position of the SNP within the contig, last three columns Fst values for the three comparisons (annuus-debilis, annuus-petiolaris, debilis-petiolaris). Species identification of individuals can be found in table 4. File was compressed (.zip).raw read count per geneNumber of reads aligned, per contig, per individuals. Column one contains contig name, column two length of the contig.raw_read_counts

Item Media

Item Citations and Data

Usage Statistics