The study of microbial diversity over time and space is fundamental to the understanding of their ecology and evolution. The underlying processes driving these patterns are not fully resolved but can be studied using population genetic approaches. Here we investigated the population genetic structure of Gonyostomum semen, a bloom-forming phytoplankton species, across two continents. The species appears to be expanding in Europe, whereas similar trends are not observed in the USA. Our aim was to investigate if populations of Gonyostomum semen in Europe and in the USA are genetically differentiated, if there is population genetic structure within the continents, and what the potential drivers of differentiation are. To this end, we used a novel method based on single-amplified genomes combined with Restriction-site Associated DNA sequencing that allows de novo genotyping of natural single-cell isolates without the need for culturing. We amplified over 900 single-cell genomes from 25 lake populations across Europe and the USA and identified two distinct population clusters, one in Europe and another in the USA. Low genetic diversity in European populations supports the hypothesized recent expansion of Gonyostomum semen on this continent. Geographic population structure within each continent was associated with differences in environmental variables that may have led to ecological divergence of population clusters. Overall, our results show that single-amplified genomes combined with Restriction-site Associated DNA sequencing can be used to analyze microalgal population structure and differentiation based on single-cell isolates from natural, uncultured samples.

Dispersal can modify how bacterial community composition (BCC) changes in response to environmental perturbations, yet knowledge about the functional consequences of dispersal is limited. Here we hypothesized that changes in bacterial community production in response to a salinity disturbance depend on the possibility to recruit cells from different dispersal sources. To investigate this, we conducted an in situ mesocosm experiment where bacterial communities of an oligotrophic lake were exposed to different salinities (0, 18, 36 psu) for 2 weeks and subjected to dispersal of cells originating from sediments, air (mesocosms open to air deposition), both or none. BCC was determined using 454 pyrosequencing of the 16S rRNA gene and bacterial production was measured by (3) H leucine uptake. Bacterial production differed significantly among salinity treatments and dispersal treatments, being highest at high salinity. These changes were associated with changes in BCC and it was found that the identity of the main functional contributors differed at different salinities. Our results further showed that after a salinity perturbation, the response of bacterial communities depended on the recruitment of taxa, including marine representatives (e.g., Alphaproteobacteria Loktanella, Erythrobacter and the Gammaproteobacterium Rheiheimera) from dispersal sources, in which atmospheric deposition appeared to play a major role.

Ectotherm development rates often show adaptive divergence along climatic gradients, but the genetic basis for this variation is rarely studied. Here, we investigated the genetic basis for phenotypic variation in larval development in the moor frog Rana arvalis from five regions along a latitudinal gradient from Germany to northern Sweden. We focused on the C/EBP-1 gene, a transcription factor associated with larval development time. Allele frequencies at C/EBP-1 varied strongly among geographical regions. Overall, the distribution of alleles along the gradient was in concordance with the dual post-glacial colonization routes into Scandinavia, with a large number of alleles exclusively present along the southern colonization route. Only three of 38 alleles were shared between the routes. Analysis of contemporary selection on C/EBP-1 showed divergent selection among the regions, probably reflecting adaptation to the local environmental conditions, although this was especially strong between southern and northern regions coinciding also with lineages from different colonization routes. Overall, the C/EBP-1 gene has historically been under purifying selection, but two specific amino acid positions showed significant signals of positive selection. These positions showed divergence between southern and northern regions, and we suggest that they are functionally involved in the climatic adaptation of larval development. Using phenotypic data from a common garden experiment, we found evidence for specific C/EBP-1 alleles being correlated with larval development time, suggesting a functional role in adaptation of larval development to large-scale climatic variation.

We analysed diverse strains of Francisella tularensis subsp. holarctica to assess if its division into biovars I and II is associated with specific mutations previously linked to erythromycin resistance and to determine the distribution of this resistance trait across this subspecies. Three-hundred and fourteen F. tularensis subsp. holarctica strains were tested for erythromycin susceptibility and whole-genome sequences for these strains were examined for SNPs in genes previously associated with erythromycin resistance. Each strain was assigned to a global phylogenetic framework using genome-wide canonical SNPs. The contribution of a specific SNP to erythromycin resistance was examined using allelic exchange. The geographical distribution of erythromycin-resistant F. tularensis strains was further investigated by literature search. There was a perfect correlation between biovar II strains (erythromycin resistance) and the phylogenetic group B.12. Only B.12 strains had an A → C SNP at position 2059 in the three copies of the rrl gene. Introducing 2059C into an rrl gene of an erythromycin-susceptible F. tularensis strain resulted in resistance. An additional 1144 erythromycin-resistant strains were identified from the scientific literature, all of them from Eurasia. Erythromycin resistance in F. tularensis is caused by an A2059C rrl gene mutation, which exhibits a strictly clonal inheritance pattern found only in phylogenetic group B.12. This group is an extremely successful clone, representing the most common type of F. tularensis throughout Eurasia.

Canine degenerative myelopathy (DM) is a naturally occurring neurodegenerative disease with similarities to some forms of amyotrophic lateral sclerosis (ALS). Most dogs that develop DM are homozygous for a common superoxide dismutase 1 gene (SOD1) mutation. However, not all dogs homozygous for this mutation develop disease. We performed a genome-wide association analysis in the Pembroke Welsh Corgi (PWC) breed comparing DM-affected and -unaffected dogs homozygous for the SOD1 mutation. The analysis revealed a modifier locus on canine chromosome 25. A haplotype within the SP110 nuclear body protein (SP110) was present in 40% of affected compared with 4% of unaffected dogs (P = 1.5 × 10(-5)), and was associated with increased probability of developing DM (P = 4.8 × 10(-6)) and earlier onset of disease (P = 1.7 × 10(-5)). SP110 is a nuclear body protein involved in the regulation of gene transcription. Our findings suggest that variations in SP110-mediated gene transcription may underlie, at least in part, the variability in risk for developing DM among PWCs that are homozygous for the disease-related SOD1 mutation. Further studies are warranted to clarify the effect of this modifier across dog breeds.

Human population history in the Holocene was profoundly impacted by changes in lifestyle following the invention and adoption of food-production practices. These changes triggered significant increases in population sizes and expansions over large distances. Here we investigate the population history of the Fulani, a pastoral population extending throughout the African Sahel/Savannah belt. Based on genome-wide analyses we propose that ancestors of the Fulani population experienced admixture between a West African group and a group carrying both European and North African ancestries. This admixture was likely coupled with newly adopted herding practices, as it resulted in signatures of genetic adaptation in contemporary Fulani genomes, including the control element of the LCT gene enabling carriers to digest lactose throughout their lives. The lactase persistence (LP) trait in the Fulani is conferred by the presence of the allele T-13910, which is also present at high frequencies in Europe. We establish that the T-13910 LP allele in Fulani individuals analysed in this study lies on a European haplotype background thus excluding parallel convergent evolution. We furthermore directly link the T-13910 haplotype with the Lactase Persistence phenotype through a Genome Wide Association study (GWAS) and identify another genomic region in the vicinity of the SPRY2 gene associated with glycaemic measurements after lactose intake. Our findings suggest that Eurasian admixture and the European LP allele was introduced into the Fulani through contact with a North African population/s. We furthermore confirm the link between the lactose digestion phenotype in the Fulani to the MCM6/LCT locus by reporting the first GWAS of the lactase persistence trait. We also explored other signals of recent adaptation in the Fulani and identified additional candidates for selection to adapt to herding life-styles.

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear. Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes. Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.