Mateusz Konczal

Accuracy of allele frequency estimation using pooled RNA-Seq

For nonmodel organisms, genome‐wide information that describes functionally relevant variation may be obtained by RNA‐Seq following de novo transcriptome assembly. While sequencing has become relatively inexpensive, the preparation of a large number of sequencing libraries remains prohibitively expensive for population genetic analyses of nonmodel species. Pooling samples may be then an attractive alternative. To test whether pooled RNA‐Seq accurately predicts true allele frequencies, we analysed the liver transcriptomes of 10 bank voles. Each sample was sequenced both as an individually barcoded library and as a part of a pool. Equal amounts of total RNA from each vole were pooled prior to mRNA selection and library construction. Reads were mapped onto the de novo assembled reference transcriptome. High‐quality genotypes for individual voles, determined for 23 682 SNPs, provided information on ‘true’ allele frequencies; allele frequencies estimated from the pool were then compared with these values. ‘True’ frequencies and those estimated from the pool were highly correlated. Mean relative estimation error was 21% and did not depend on expression level. However, we also observed a minor effect of interindividual variation in gene expression and allele‐specific gene expression influencing allele frequency estimation accuracy. Moreover, we observed strong negative relationship between minor allele frequency and relative estimation error. Our results indicate that pooled RNA‐Seq exhibits accuracy comparable with pooled genome resequencing, but variation in expression level between individuals should be assessed and accounted for. This should help in taking account the difference in accuracy between conservatively expressed transcripts and these which are variable in expression level.

Development, validation and high-throughput analysis of sequence markers in nonmodel species

DNA sequences derived from multiple regions of the nuclear genome are essential for historical inferences in the fields of phylogeography and phylogenetics. The appropriate markers should be single‐copy, variable, easy to amplify from multiple samples and easy to sequence using high‐throughput technologies. This may be difficult to achieve for species lacking sequenced genomes and particularly challenging for species possessing large genomes, which consist mostly of repetitive sequences. Here, we present a cost‐effective, broadly applicable framework for designing, validating and high‐throughput sequencing of multiple markers in nonmodel species without sequenced genomes. We demonstrate its utility in two closely related species of newts, representatives of urodeles, a vertebrate group characterized by large genomes. We show that over 80 markers, c. 600 bp each, developed mainly from 3′ untranslated transcript regions (3′UTR) may be effectively multiplexed and sequenced. Data are further processed using standard, freely available bioinformatic tools, producing phase‐resolved sequences. The approach does not require barcoded PCR primers, and the cost of library preparation is independent of the number of markers investigated. We hope that this approach will be of broad interest for researchers working at the interface of population genetics and phylogenetics, exploring deep intraspecific genetic structure, species boundaries and phylogeographies of closely related species.

Initial Molecular-Level Response to Artificial Selection for Increased Aerobic Metabolism Occurs Primarily through Changes in Gene Expression

Experimental evolution combined with genome or transcriptome resequencing (Evolve and Resequence) represents a promising approach for advancing our understanding of the genetic basis of adaptation. Here, we applied this strategy to investigate the effect of selection on a complex trait in lines derived from a natural population of a small mammal. We analyzed the liver and heart transcriptomes of bank voles (Myodes [=Clethrionomys] glareolus) that had been selected for increased aerobic metabolism. The organs were sampled from 13th generation voles; at that point, the voles from four replicate selected lines had 48% higher maximum rates of oxygen consumption than those from four control lines. At the molecular level, the response to selection was primarily observed in gene expression: Over 300 genes were found to be differentially expressed between the selected and control lines and the transcriptome-wide pattern of expression distinguished selected lines from controls. No evidence for selection-driven changes of allele frequencies at coding sites was found: No single nucleotide polymorphism (SNP) changed frequency more than expected under drift alone and frequency changes aggregated over all SNPs did not separate selected and control lines. Nevertheless, among genes which showed highest differentiation in allele frequencies between selected and control lines we identified, using information about gene functions and the biology of the selected phenotype, plausible targets of selection; these genes, together with those identified in expression analysis, have been prioritized for further studies. Because our selection lines were derived from a natural population, the amount and the spectrum of variation available for selection probably closely approximated that typically found in populations of small mammals. Therefore, our results are relevant to the understanding of the molecular basis of complex adaptations occurring in natural vertebrate populations.

Genomic Response to Selection for Predatory Behavior in a Mammalian Model of Adaptive Radiation.

If genetic architectures of various quantitative traits are similar, as studies on model organisms suggest, comparable selection pressures should produce similar molecular patterns for various traits. To test this prediction, we used a laboratory model of vertebrate adaptive radiation to investigate the genetic basis of the response to selection for predatory behavior and compare it with evolution of aerobic capacity reported in an earlier work. After 13 generations of selection, the proportion of bank voles (Myodes [=Clethrionomys] glareolus) showing predatory behavior was five times higher in selected lines than in controls. We analyzed the hippocampus and liver transcriptomes and found repeatable changes in allele frequencies and gene expression. Genes with the largest differences between predatory and control lines are associated with hunger, aggression, biological rhythms, and functioning of the nervous system. Evolution of predatory behavior could be meaningfully compared with evolution of high aerobic capacity, because the experiments and analyses were performed in the same methodological framework. The number of genes that changed expression was much smaller in predatory lines, and allele frequencies changed repeatably in predatory but not in aerobic lines. This suggests that more variants of smaller effects underlie variation in aerobic performance, whereas fewer variants of larger effects underlie variation in predatory behavior. Our results thus contradict the view that comparable selection pressures for different quantitative traits produce similar molecular patterns. Therefore, to gain knowledge about molecular-level response to selection for complex traits, we need to investigate not only multiple replicate populations but also multiple quantitative traits.

De novo transcriptome assembly facilitates characterisation of fast-evolving gene families, MHC class I in the bank vole ( Myodes glareolus )

The major histocompatibility complex (MHC) plays a central role in the adaptive immune response and is the most polymorphic gene family in vertebrates. Although high-throughput sequencing has increasingly been used for genotyping families of co-amplifying MHC genes, its potential to facilitate early steps in the characterisation of MHC variation in nonmodel organism has not been fully explored. In this study we evaluated the usefulness of de novo transcriptome assembly in characterisation of MHC sequence diversity. We found that although de novo transcriptome assembly of MHC I genes does not reconstruct sequences of individual alleles, it does allow the identification of conserved regions for PCR primer design. Using the newly designed primers, we characterised MHC I sequences in the bank vole. Phylogenetic analysis of the partial MHC I coding sequence (2–4 exons) of the bank vole revealed a lack of orthology to MHC I of other Cricetidae, consistent with the high gene turnover of this region. The diversity of expressed alleles was characterised using ultra-deep sequencing of the third exon that codes for the peptide-binding region of the MHC molecule. High allelic diversity was demonstrated, with 72 alleles found in 29 individuals. Interindividual variation in the number of expressed loci was found, with the number of alleles per individual ranging from 5 to 14. Strong signatures of positive selection were found for 8 amino acid sites, most of which are inferred to bind antigens in human MHC, indicating conservation of structure despite rapid sequence evolution.

Transcriptomics of Intralocus Sexual Conflict: Gene Expression Patterns in Females Change in Response to Selection on a Male Secondary Sexual Trait in the Bulb Mite

Intralocus sexual conflict (IASC) prevents males and females from reaching their disparate phenotypic optima and is widespread, but little is known about its genetic underpinnings. In Rhizoglyphus robini, a mite species with alternative male morphs, elevated sexual dimorphism of the armored fighter males (compared to more feminized scramblers males) was previously reported to be associated with increased IASC. Because IASC persists if gene expression patterns are correlated between sexes, we compared gene expression patterns of males and females from the replicate lines selected for increased proportion of fighter or scrambler males (F- and S-lines, respectively). Specifically, we tested the prediction that selection for fighter morph caused correlated changes in gene expression patterns in females. We identified 532 differentially expressed genes (FDR < 0.05) between the F-line and S-line males. Consistent with the prediction, expression levels of these genes also differed between females from respective lines. Thus, significant proportion of genes differentially expressed between sexually selected male phenotypes showed correlated expression levels in females, likely contributing to elevated IASC in F-lines reported in a previous study.

Genomics of end-Pleistocene population replacement in a small mammal

Current species distributions at high latitudes are the product of expansion from glacial refugia into previously uninhabitable areas at the end of the last glaciation. The traditional view of postglacial colonization is that southern populations expanded their ranges into unoccupied northern territories. Recent findings on mitochondrial DNA (mtDNA) of British small mammals have challenged this simple colonization scenario by demonstrating a more complex genetic turnover in Britain during the Pleistocene–Holocene transition where one mtDNA clade of each species was replaced by another mtDNA clade of the same species. Here, we provide evidence from one of those small mammals, the bank vole (Clethrionomys glareolus), that the replacement was genome-wide. Using more than 10 000 autosomal SNPs we found that similar to mtDNA, bank vole genomes in Britain form two (north and south) clusters which admix. Therefore, the genome of the original postglacial colonists (the northern cluster) was probably replaced by another wave of migration from a different continental European population (the southern cluster), and we gained support for this by modelling with approximate Bayesian computation. This finding emphasizes the importance of analysis of genome-wide diversity within species under changing climate in creating opportunities for sophisticated testing of population history scenarios.