Exploring the unmapped DNA and RNA reads in a songbird genome.

Background: A widely used approach in next-generation sequencing projects is the alignment of reads to a reference genome. Despite methodological and hardware improvements which have enhanced the efficiency and accuracy of alignments, a significant percentage of reads frequently remain unmapped. Usually, unmapped reads are discarded from the analysis process, but significant biological information and insights can be uncovered from these data. We explored the unmapped DNA (normal and bisulfite treated) and RNA sequence reads of the great tit (Parus major) reference genome individual. From the unmapped reads we generated de novo assemblies, after which the generated sequence contigs were aligned to the NCBI non-redundant nucleotide database using BLAST, identifying the closest known matching sequence. Results: Many of the aligned contigs showed sequence similarity to different bird species and genes that were absent in the great tit reference assembly. Furthermore, there were also contigs that represented known P. major pathogenic species. Most interesting were several species of blood parasites such as Plasmodium and Trypanosoma. Conclusions: Our analyses revealed that meaningful biological information can be found when further exploring unmapped reads. For instance, it is possible to discover sequences that are either absent or misassembled in the reference genome, and sequences that indicate infection or sample contamination. In this study we also propose strategies to aid the capture and interpretation of this information from unmapped reads

Within-season changes in habitat use of forest-dwelling boreal bats

Bats utilize forests as roosting sites and feeding areas. However, it has not been documented how bats utilize these habitats in the boreal zone with methods afforded by recent technological advances. Forest structure and management practices can create a variety of three-dimensional habitats for organisms capable of flight, such as bats. Here, we study the presence of boreal bats in a forest forming a mosaic of different age classes, dominant tree species, canopy cover, soil fertility, and other environmental variables, throughout their active season in the summer using passive ultrasound detectors. Our results indicate a preference for mature forest by Eptesicus nilssonii and a pooled set of Myotis bats. Both groups of bats also showed temporal changes in their habitat use regarding forest age. In June and July, both groups occurred more often in mature than young forests, but from August onwards, the difference in occurrence became less evident in Myotis and disappeared completely in E. nilssonii. In addition, E. nilssonii was more often present in forests with low canopy cover, and its occurrence shifted from coniferous forests to deciduous forests during the season. The results reflect the within-season dynamics of bat communities and their ability to utilize different types of forest as environmental conditions change. Yet, the results most importantly emphasize the importance of mature forests to bat diversity and the need to conserve such environments in the boreal zone.Peer reviewe

On the hadronic contribution to sterile neutrino production

Sterile neutrinos with masses in the keV range are considered to be a viable candidate for warm dark matter. The rate of their production through active-sterile neutrino transitions peaks, however, at temperatures of the order of the QCD scale, which makes it difficult to estimate their relic abundance quantitatively, even if the mass of the sterile neutrino and its mixing angle were known. We derive here a relation, valid to all orders in the strong coupling constant, which expresses the production rate in terms of the spectral function associated with active neutrinos. The latter can in turn be expressed as a certain convolution of the spectral functions related to various mesonic current-current correlation functions, which are being actively studied in other physics contexts. In the naive weak coupling limit, the appropriate Boltzmann equations can be derived from our general formulae.Comment: 28 pages. v2: small clarifications added, published versio

Recent natural selection causes adaptive evolution of an avian polygenic trait

We used extensive data from a long-term study of great tits (Parus major) in the United Kingdom and Netherlands to better understand how genetic signatures of selection translate into variation in fitness and phenotypes. We found that genomic regions under differential selection contained candidate genes for bill morphology and used genetic architecture analyses to confirm that these genes, especially the collagen gene COL4A5, explained variation in bill length. COL4A5 variation was associated with reproductive success, which, combined with spatiotemporal patterns of bill length, suggested ongoing selection for longer bills in the United Kingdom. Last, bill length and COL4A5 variation were associated with usage of feeders, suggesting that longer bills may have evolved in the United Kingdom as a response to supplementary feeding

Evolutionary signals of selection on cognition from the great tit genome and methylome

For over 50 years, the great tit (Parus major) has been a model species for research in evolutionary, ecological and behavioural research; in particular, learning and cognition have been intensively studied. Here, to provide further insight into the molecular mechanisms behind these important traits, we de novo assemble a great tit reference genome and whole-genome re-sequence another 29 individuals from across Europe. We show an overrepresentation of genes related to neuronal functions, learning and cognition in regions under positive selection, as well as increased CpG methylation in these regions. In addition, great tit neuronal non-CpG methylation patterns are very similar to those observed in mammals, suggesting a universal role in neuronal epigenetic regulation which can affect learning-, memory- and experience-induced plasticity. The high-quality great tit genome assembly will play an instrumental role in furthering the integration of ecological, evolutionary, behavioural and genomic approaches in this model species.</p

The genomics of adaptation to climate in European great tit (Parus major) populations

The recognition that climate change is occurring at an unprecedented rate means that there is increased urgency in understanding how organisms can adapt to a changing environment. Wild great tit (Parus major) populations represent an attractive ecological model system to understand the genomics of climate adaptation. They are widely distributed across Eurasia and they have been documented to respond to climate change. We performed a Bayesian genome-environment analysis, by combining local climate data with single nucleotide polymorphisms genotype data from 20 European populations (broadly spanning the species’ continental range). We found 36 genes putatively linked to adaptation to climate. Following an enrichment analysis of biological process Gene Ontology (GO) terms, we identified over-represented terms and pathways among the candidate genes. Because many different genes and GO terms are associated with climate variables, it seems likely that climate adaptation is polygenic and genetically complex. Our findings also suggest that geographical climate adaptation has been occurring since great tits left their Southern European refugia at the end of the last ice age. Finally, we show that substantial climate-associated genetic variation remains, which will be essential for adaptation to future changes

Dynamical inference of hidden biological populations.

Population fluctuations in a predator-prey system are analyzed for the case where the number of prey could be determined, subject to measurement noise, but the number of predators was unknown. The problem of how to infer the unmeasured predator dynamics, as well as the model parameters, is addressed. Two solutions are suggested. In the first of these, measurement noise and the dynamical noise in the equation for predator population are neglected; the problem is reduced to a one-dimensional case, and a Bayesian dynamical inference algorithm is employed to reconstruct the model parameters. In the second solution a full-scale Markov Chain Monte Carlo simulation is used to infer both the unknown predator trajectory, and also the model parameters, using the one-dimensional solution as an initial guess

The great tit HapMap project: a continental‐scale analysis of genomic variation in a songbird

A major aim of evolutionary biology is to understand why patterns of genomic diversity vary within taxa and space. Large-scale genomic studies of widespread species are useful for studying how environment and demography shape patterns of genomic divergence. Here, we describe one of the most geographically comprehensive surveys of genomic variation in a wild vertebrate to date; the great tit (Parus major) HapMap project. We screened ca 500,000 SNP markers across 647 individuals from 29 populations, spanning ~30 degrees of latitude and 40 degrees of longitude – almost the entire geographical range of the European subspecies. Genome-wide variation was consistent with a recent colonisation across Europe from a South-East European refugium, with bottlenecks and reduced genetic diversity in island populations. Differentiation across the genome was highly heterogeneous, with clear ‘islands of differentiation’, even among populations with very low levels of genome-wide differentiation. Low local recombination rates were a strong predictor of high local genomic differentiation (FST), especially in island and peripheral mainland populations, suggesting that the interplay between genetic drift and recombination causes highly heterogeneous differentiation landscapes. We also detected genomic outlier regions that were confined to one or more peripheral great tit populations, probably as a result of recent directional selection at the species' range edges. Haplotype-based measures of selection were related to recombination rate, albeit less strongly, and highlighted population-specific sweeps that likely resulted from positive selection. Our study highlights how comprehensive screens of genomic variation in wild organisms can provide unique insights into spatio-temporal evolutionary dynamics

Parus major Genome sequencing and assembly

Seasonal timing is a key life-history trait with major fitness consequences. The small songbird Parus major (great tit) for decades has been a model to study fitness traits like e.g. avian timing of reproduction. The research is closely linked to the impact of global climate change on timing and its consequences. Linking quantitative genetic variation in life-history traits with polymorphisms in specific genes is essential for better understanding the causes and consequences of the diversity in these traits. Genetic variation in life-history traits in wild songbirds has been demonstrated in many, often long-term, studies throughout the world. Linking this variation to genomic information requires the development of the necessary genomics tools specifically aimed at these non-model species. The assembly and annotation of the genome of the great tit will greatly enhance the further use of the great tit as a model species in this research field

Parus major Genome sequencing and assembly

Seasonal timing is a key life-history trait with major fitness consequences. The small songbird Parus major (great tit) for decades has been a model to study fitness traits like e.g. avian timing of reproduction. The research is closely linked to the impact of global climate change on timing and its consequences. Linking quantitative genetic variation in life-history traits with polymorphisms in specific genes is essential for better understanding the causes and consequences of the diversity in these traits. Genetic variation in life-history traits in wild songbirds has been demonstrated in many, often long-term, studies throughout the world. Linking this variation to genomic information requires the development of the necessary genomics tools specifically aimed at these non-model species. The assembly and annotation of the genome of the great tit will greatly enhance the further use of the great tit as a model species in this research field