Molecular Inversion Probes for targeted resequencing in non-model organisms

Applications that require resequencing of hundreds or thousands of predefined genomic regions in numerous samples are common in studies of non-model organisms. However few approaches at the scale intermediate between multiplex PCR and sequence capture methods are available. Here we explored the utility of Molecular Inversion Probes (MIPs) for the medium-scale targeted resequencing in a non-model system. Markers targeting 112 bp of exonic sequence were designed from transcriptome of Lissotriton newts. We assessed performance of 248 MIP markers in a sample of 85 individuals. Among the 234 (94.4%) successfully amplified markers 80% had median coverage within one order of magnitude, indicating relatively uniform performance; coverage uniformity across individuals was also high. In the analysis of polymorphism and segregation within family, 77% of 248 tested MIPs were confirmed as single copy Mendelian markers. Genotyping concordance assessed using replicate samples exceeded 99%. MIP markers for targeted resequencing have a number of advantages: high specificity, high multiplexing level, low sample requirement, straightforward laboratory protocol, no need for preparation of genomic libraries and no ascertainment bias. We conclude that MIP markers provide an effective solution for resequencing targets of tens or hundreds of kb in any organism and in a large number of samples

Linkage map of Lissotriton newts provides insight into the genetic basis of reproductive isolation

Linkage maps are widely used to investigate structure, function, and evolution of genomes. In speciation research, maps facilitate the study of the genetic architecture of reproductive isolation by allowing identification of genomic regions underlying reduced fitness of hybrids. Here we present a linkage map for European newts of the Lissotriton vulgaris species complex, constructed using two families of F2 L. montandoni × L. vulgaris hybrids. The map consists of 1146 protein-coding genes on 12 linkage groups, equal to the haploid chromosome number, with a total length of 1484 cM (1.29 cM per marker). It is notably shorter than two other maps available for salamanders, but the differences in map length are consistent with cytogenetic estimates of the number of chiasmata per chromosomal arm. Thus, large salamander genomes do not necessarily translate into long linkage maps, as previously suggested. Consequently, salamanders are an excellent model to study evolutionary consequences of recombination rate variation in taxa with large genomes and a similar number of chromosomes. A complex pattern of transmission ratio distortion (TRD) was detected: TRD occurred mostly in one family, in one breeding season, and was clustered in two genomic segments. This is consistent with environment-dependent mortality of individuals carrying L. montandoni alleles in these two segments and suggests a role of TRD blocks in reproductive isolation. The reported linkage map will empower studies on the genomic architecture of divergence and interactions between the genomes of hybridizing newts

Constraint and adaptation in newt Toll-like receptor genes

Acute die-offs of amphibian populations worldwide have been linked to the emergence of viral and fungal diseases. Inter and intraspecific immunogenetic differences may influence the outcome of infection. Toll-like receptors (TLRs) are an essential component of innate immunity and also prime acquired defenses. We report the first comprehensive assessment of TLR gene variation for urodele amphibians. The Lissotriton newt TLR repertoire includes representatives of 13 families and is compositionally most similar to that of the anuran Xenopus. Both ancient and recent gene duplications have occurred in urodeles, bringing the total number of TLR genes to at least 21. Purifying selection has predominated the evolution of newt TLRs in both long (∼70 Ma) and medium (∼18 Ma) timescales. However, we find evidence for both purifying and positive selection acting on TLRs in two recently diverged (2–5 Ma) allopatric evolutionary lineages (Lissotriton montandoni and L. vulgaris graecus). Overall, both forms of selection have been stronger in L. v. graecus, while constraint on most TLR genes in L. montandoni appears relaxed. The differences in selection regimes are unlikely to be biased by demographic effects because these were controlled by means of a historical demographic model derived from an independent data set of 62 loci. We infer that TLR genes undergo distinct trajectories of adaptive evolution in closely related amphibian lineages, highlight the potential of TLRs to capture the signatures of different assemblages of pathogenic microorganisms, and suggest differences between lineages in the relative roles of innate and acquired immunity

Limited diversity associated with duplicated class II MHC-DRB genes in the red squirrel population in the United Kingdom compared with continental Europe

The red squirrel (Sciurus vulgaris) population in the United Kingdom has declined over the last century and is now on the UK endangered species list. This is the result of competition from the eastern grey squirrel (S. carolinensis) which was introduced in the 19th century. However, recent evidence suggests that the rate of population decline is enhanced by squirrelpox disease, caused by a viral infection carried asymptomatically by grey squirrels but to which red squirrels are highly susceptible. Population genetic diversity provides some resilience to rapidly evolving or exotic pathogens. There is currently no data on genetic diversity of extant UK squirrel populations with respect to genes involved in disease resistance. Diversity is highest at loci involved in the immune response including genes clustered within the major histocompatibility complex (MHC). Using the class II DRB locus as a marker for diversity across the MHC region we genotyped 110 red squirrels from locations in the UK and continentalEurope. Twenty four Scvu-DRB alleles at two functional loci; Scvu-DRB1 and Scvu- DRB2, were identified. High levels of diversity were identified at both loci in the continental populations. In contrast, no diversity was observed at the Scvu-DRB2 locus in the mainland UK population while a high level of homozygosity was observed at the Scvu-DRB1 locus. The red squirrel population in the UK appears to lack the extensive MHC diversity associated with continental populations, a feature which may have contributed to their rapid decline

Sequence comparison of prefrontal cortical brain transcriptome from a tame and an aggressive silver fox (Vulpes vulpes)

<p>Abstract</p> <p>Background</p> <p>Two strains of the silver fox (<it>Vulpes vulpes</it>), with markedly different behavioral phenotypes, have been developed by long-term selection for behavior. Foxes from the tame strain exhibit friendly behavior towards humans, paralleling the sociability of canine puppies, whereas foxes from the aggressive strain are defensive and exhibit aggression to humans. To understand the genetic differences underlying these behavioral phenotypes fox-specific genomic resources are needed.</p> <p>Results</p> <p>cDNA from mRNA from pre-frontal cortex of a tame and an aggressive fox was sequenced using the Roche 454 FLX Titanium platform (> 2.5 million reads & 0.9 Gbase of tame fox sequence; >3.3 million reads & 1.2 Gbase of aggressive fox sequence). Over 80% of the fox reads were assembled into contigs. Mapping fox reads against the fox transcriptome assembly and the dog genome identified over 30,000 high confidence fox-specific SNPs. Fox transcripts for approximately 14,000 genes were identified using SwissProt and the dog RefSeq databases. An at least 2-fold expression difference between the two samples (p < 0.05) was observed for 335 genes, fewer than 3% of the total number of genes identified in the fox transcriptome.</p> <p>Conclusions</p> <p>Transcriptome sequencing significantly expanded genomic resources available for the fox, a species without a sequenced genome. In a very cost efficient manner this yielded a large number of fox-specific SNP markers for genetic studies and provided significant insights into the gene expression profile of the fox pre-frontal cortex; expression differences between the two fox samples; and a catalogue of potentially important gene-specific sequence variants. This result demonstrates the utility of this approach for developing genomic resources in species with limited genomic information.</p

Selective landscapes in newt immune genes inferred from patterns of nucleotide variation

Host–pathogen interactions may result in either directional selection or in pressure for the maintenance of polymorphism at the molecular level. Hence signatures of both positive and balancing selection are expected in immune genes. Because both overall selective pressure and specific targets may differ between species, large-scale population genomic studies are useful in detecting functionally important immune genes and comparing selective landscapes between taxa. Such studies are of particular interest in amphibians, a group threatened worldwide by emerging infectious diseases. Here, we present an analysis of polymorphism and divergence of 634 immune genes in two lineages of Lissotriton newts: L. montandoni and L. vulgaris graecus. Variation in newt immune genes has been shaped predominantly by widespread purifying selection and strong evolutionary constraint, implying long-term importance of these genes for functioning of the immune system. The two evolutionary lineages differ in the overall strength of purifying selection which can partially be explained by demographic history but may also signal differences in long-term pathogen pressure. The prevalent constraint notwithstanding, 23 putative targets of positive selection and 11 putative targets of balancing selection were identified. The latter were detected by composite tests involving the demographic model and further validated in independent population samples. Putative targets of balancing selection encode proteins which may interact closely with pathogens but include also regulators of immune response. The identified candidates will be useful for testing whether genes affected by balancing selection are more prone to interspecific introgression than other genes in the genome

Location learning in virtual environments: The effect of saliency of landmarks and boundaries

In the present paper, we investigated the learning of object location in a virtual reality environment. Previous studies suggested that learning is essentially different when participants rely on landmarks (i.e., single stationary objects in the space) or boundaries (i.e., terrain features), the former requiring associative reinforcement learning, while the latter relying on incidental learning. However, in previous studies the saliency of the spatial cues and the position of objects relative to spatial cues were not carefully controlled, which could account for the specific learning effect found. In our study, we controlled these effects. We presented a virtual arena on a computer screen, featuring a landmark and a boundary cue, and four everyday objects, the location of which had to be learned and recalled. In the learning phase, either one or both cues were present, while in the test phase only one cue was present. Contrary to previous results, we did not obtain the overshadowing effect, as there was no difference in the location errors between the four conditions. At the same time, we found that location errors were influenced by the relative spatial position of the object and cue. This suggests that it is the distance from the cue and not the cue type, which is the main factor behind localization errors

Genotypes for BS11 genes in Lissotriton vulgaris vulgaris

Genotypes called with GATK in Lissotriton vulgaris vulgari

SLiM simulation script moderate NFD M0

Forward simulation script for the scenario: moderate negative frequency-dependent selection, migration M0 = 0. Script run in SLiM v2.5. Example command to run the script: "./slim -d simnum=00000 -d u=1.17e-7 -d locus=0 example_simulation_6_pops_00000_0.txt". example_simulation_6_pops_00000_0.txt - name of the SLiM script with one simulation specification; simnum - simulation number; u - locus mutation rate; locus - locus number