Avian wing length is an important trait that covaries with the ecology and migratory behaviour of a species and tends to change rapidly when the conditions are altered. Long-distance migrants typically have longer wings than short-distance migrants and sedentary species, and long-winged species also tend to be more dispersive. Although the substantial heritability of avian wing length is well established, the identification of causal genes has remained elusive. Based on large-scale genotyping of 1404 informative single nucleotide polymorphisms (SNP) in a captive population of 1067 zebra finches, we here show that the within-population variation of relative wing length (h(2) = 0.74 ± 0.05) is associated with standing genetic variation in at least six genomic regions (one genome-wide significant and five suggestive). The variance explained by these six quantitative trait loci (QTL) sums to 36.8% of the phenotypic variance (half of the additive genetic variance), although this likely is an overestimate attributable to the Beavis effect. As avian wing length is primarily determined by the length of the primary feathers, we then searched for candidate genes that are related to feather growth. Interestingly, all of the QTL signals co-locate with Wnt growth factors and closely interacting genes (Wnt3a, Wnt5a, Wnt6, Wnt7a, Wnt9a, RhoU and RhoV). Our findings therefore suggest that standing genetic variation in the Wnt genes might be linked to avian wing morphology, although there are many other genes that also fall within the confidence regions.

A multidrug-resistant lineage of Staphylococcus epidermidis named ST215 is a common cause of prosthetic joint infections and other deep surgical site infections in Northern Europe, but is not present elsewhere. The increasing resistance among S. epidermidis strains is a global concern. We used whole-genome sequencing to characterize ST215 from healthcare settings. We completed the genome of a ST215 isolate from a Swedish hospital using short and long reads, resulting in a circular 2,676,787 bp chromosome and a 2,326 bp plasmid. The new ST215 genome was placed in phylogenetic context using 1,361 finished public S. epidermidis reference genomes. We generated 10 additional short-read ST215 genomes and 11 short-read genomes of ST2, which is another common multidrug-resistant lineage at the same hospital. We studied recombination's role in the evolution of ST2 and ST215, and found multiple recombination events averaging 30-50 kb. By comparing the results of antimicrobial susceptibility testing for 31 antimicrobial drugs with the genome content encoding antimicrobial resistance in the ST215 and ST2 isolates, we found highly similar resistance traits between the isolates, with 22 resistance genes being shared between all the ST215 and ST2 genomes. The ST215 genome contained 29 genes that were historically identified as virulence genes of S. epidermidis ST2. We established that in the nucleotide sequence stretches identified as recombination events, virulence genes were overrepresented in ST215, while antibiotic resistance genes were overrepresented in ST2. This study features the extensive antibiotic resistance and virulence gene content in ST215 genomes. ST215 and ST2 lineages have similarly evolved, acquiring resistance and virulence through genomic recombination. The results highlight the threat of new multidrug-resistant S. epidermidis lineages emerging in healthcare settings.

Ancient remains found in permafrost represent a rare opportunity to study past ecosystems. Here, we present an exceptionally well-preserved ancient bird carcass found in the Siberian permafrost, along with a radiocarbon date and a reconstruction of its complete mitochondrial genome. The carcass was radiocarbon dated to approximately 44-49 ka BP, and was genetically identified as a female horned lark. This is a species that usually inhabits open habitat, such as the steppe environment that existed in Siberia at the time. This near-intact carcass highlights the potential of permafrost remains for evolutionary studies that combine both morphology and ancient nucleic acids.

To identify regions of aberrant DNA methylation in acute lymphoblastic leukemia (ALL) cells of different subtypes on a genome-wide scale. Whole-genome bisulfite sequencing (WGBS) was used to determine the DNA methylation levels in cells from four pediatric ALL patients of different subtypes. The findings were confirmed by 450k DNA methylation arrays in a large patient set. Compared with mature B or T cells WGBS detected on average 82,000 differentially methylated regions per patient. Differentially methylated regions are enriched to CpG poor regions, active enhancers and transcriptional start sites. We also identified approximately 8000 CpG islands with variable intermediate DNA methylation that seems to occur as a result of stochastic de novo methylation. WGBS provides an unbiased view and novel insights into the DNA methylome of ALL cells.

Most mutations in coding regions of the genome are deleterious, causing selection to favor mechanisms that minimize the mutational load over time [1-5]. DNA replication during cell division is a major source of new mutations. It is therefore important to limit the number of cell divisions between generations, particularly for large and long-lived organisms [6-9]. The germline cells of animals and the slowly dividing cells in plant meristems are adaptations to control the number of mutations that accumulate over generations [9-11]. Fungi lack a separated germline while harboring species with very large and long-lived individuals that appear to maintain highly stable genomes within their mycelia [8, 12, 13]. Here, we studied genomic mutation accumulation in the fairy-ring mushroom Marasmius oreades. We generated a chromosome-level genome assembly using a combination of cutting-edge DNA sequencing technologies and re-sequenced 40 samples originating from six individuals of this fungus. The low number of mutations recovered in the sequencing data suggests the presence of an unknown mechanism that works to maintain extraordinary genome integrity over vegetative growth in M. oreades. The highly structured growth pattern of M. oreades allowed us to estimate the number of cell divisions leading up to each sample [14, 15], and from this data, we infer an incredibly low per mitosis mutation rate (3.8 × 10 -12 mutations per site and cell division) as one of several possible explanations for the low number of identified mutations.

Transcontinental migration is a fascinating example of how animals can respond to climatic oscillation. Yet, quantitative data on fitness components are scarce, and the resulting population genetic consequences are poorly understood. Migratory divides, hybrid zones with a transition in migratory behaviour, provide a natural setting to investigate the micro-evolutionary dynamics induced by migration under sympatric conditions. Here, we studied the effects of migratory programme on survival, trait evolution and genome-wide patterns of population differentiation in a migratory divide of European barn swallows. We sampled a total of 824 individuals from both allopatric European populations wintering in central and southern Africa, respectively, along with two mixed populations from within the migratory divide. While most morphological characters varied by latitude consistent with Bergmann's rule, wing length co-varied with distance to wintering grounds. Survival data collected during a 5-year period provided strong evidence that this covariance is repeatedly generated by disruptive selection against intermediate phenotypes. Yet, selection-induced divergence did not translate into genome-wide genetic differentiation as assessed by microsatellites, mtDNA and >20 000 genome-wide SNP markers; nor did we find evidence of local genomic selection between migratory types. Among breeding populations, a single outlier locus mapped to the BUB1 gene with a role in mitotic and meiotic organization. Overall, this study provides evidence for an adaptive response to variation in migration behaviour continuously eroded by gene flow under current conditions of nonassortative mating. It supports the theoretical prediction that population differentiation is difficult to achieve under conditions of gene flow despite measurable disruptive selection.

We investigated associations between genetic variation in candidate genes and on a genome-wide scale with warfarin maintenance dose, time in therapeutic range (TTR), and risk of major bleeding. In total, 982 warfarin-treated patients from the RE-LY trial were studied. After adjusting for SNPs in VKORC1 and CYP2C9, SNPs in DDHD1 (rs17126068) and NEDD4 (rs2288344) were associated with dose. Adding these SNPs and CYP4F2 (rs2108622) to a base model increased R(2) by 2.9%. An SNP in ASPH (rs4379440) was associated with TTR (-6.8% per minor allele). VKORC1 was associated with time less than INR 2.0. VKORC1 and CYP2C9 were associated with time more than INR 3.0, but not with major bleeding. We identified two novel genes associated with warfarin maintenance dose and one gene associated with TTR. These genes need to be replicated in an independent cohort.