DNA hydroxymethylation differences underlie phenotypic divergence of somatic growth in Nile tilapia reared in common garden.

Phenotypic plasticity of metabolism and growth are essential for adaptation to new environmental conditions, such as those experienced during domestication. Epigenetic regulation plays a key role in this process but the underlying mechanisms are poorly understood, especially in the case of hydroxymethylation. Using reduced representation 5-hydroxymethylcytosine profiling, we compared the liver hydroxymethylomes in full-sib Nile tilapia with distinct growth rates (3.8-fold difference) and demonstrated that DNA hydroxymethylation is strongly associated with phenotypic divergence of somatic growth during the early stages of domestication. The 2677 differentially hydroxymethylated cytosines between fast- and slow-growing fish were enriched within gene bodies (79%), indicating a pertinent role in transcriptional regulation. Moreover, they were found in genes involved in biological processes related to skeletal system and muscle structure development, and there was a positive association between somatic growth and 5hmC levels in genes coding for growth factors, kinases and receptors linked to myogenesis. Single nucleotide polymorphism analysis revealed no genetic differentiation between fast- and slow-growing fish. In addition to unveiling a new link between DNA hydroxymethylation and epigenetic regulation of growth in fish during the initial stages of domestication, this study suggests that epimarkers may be applied in selective breeding programmes for superior phenotypes.

Hybridization of Atlantic puffins in the Arctic coincides with 20th-century climate change.

The Arctic is experiencing the fastest rates of global warming, leading to shifts in the distribution of its biota and increasing the potential for hybridization. However, genomic evidence of recent hybridization events in the Arctic remains unexpectedly rare. Here, we use whole-genome sequencing of contemporary and 122-year-old historical specimens to investigate the origin of an Arctic hybrid population of Atlantic puffins (Fratercula arctica) on Bjørnøya, Norway. We show that the hybridization between the High Arctic, large-bodied subspecies F. a. naumanni and the temperate, smaller-sized subspecies F. a. arctica began as recently as six generations ago due to an unexpected southward range expansion of F. a. naumanni. Moreover, we find a significant temporal loss of genetic diversity across Arctic and temperate puffin populations. Our observations provide compelling genomic evidence of the impacts of recent distributional shifts and loss of diversity in Arctic communities during the 20th century.

Teleost genomic repeat landscapes in light of diversification rates and ecology.

Repetitive DNA make up a considerable fraction of most eukaryotic genomes. In fish, transposable element (TE) activity has coincided with rapid species diversification. Here, we annotated the repetitive content in 100 genome assemblies, covering the major branches of the diverse lineage of teleost fish. We investigated if TE content correlates with family level net diversification rates and found support for a weak negative correlation. Further, we demonstrated that TE proportion correlates with genome size, but not to the proportion of short tandem repeats (STRs), which implies independent evolutionary paths. Marine and freshwater fish had large differences in STR content, with the most extreme propagation detected in the genomes of codfish species and Atlantic herring. Such a high density of STRs is likely to increase the mutational load, which we propose could be counterbalanced by high fecundity as seen in codfishes and herring.

Adaptation to the High-Arctic island environment despite long-term reduced genetic variation in Svalbard reindeer.

Typically much smaller in number than their mainland counterparts, island populations are ideal systems to investigate genetic threats to small populations. The Svalbard reindeer (Rangifer tarandus platyrhynchus) is an endemic subspecies that colonized the Svalbard archipelago ca. 6,000-8,000 years ago and now shows numerous physiological and morphological adaptations to its arctic habitat. Here, we report a de-novo chromosome-level assembly for Svalbard reindeer and analyze 133 reindeer genomes spanning Svalbard and most of the species' Holarctic range, to examine the genomic consequences of long-term isolation and small population size in this insular subspecies. Empirical data, demographic reconstructions, and forward simulations show that long-term isolation and high inbreeding levels may have facilitated the reduction of highly deleterious-and to a lesser extent, moderately deleterious-variation. Our study indicates that long-term reduced genetic diversity did not preclude local adaptation to the High Arctic, suggesting that even severely bottlenecked populations can retain evolutionary potential.

Adaptive protein evolution through length variation of short tandem repeats in Arabidopsis.

Intrinsically disordered protein regions are of high importance for biotic and abiotic stress responses in plants. Tracts of identical amino acids accumulate in these regions and can vary in length over generations because of expansions and retractions of short tandem repeats at the genomic level. However, little attention has been paid to what extent length variation is shaped by natural selection. By environmental association analysis on 2514 length variable tracts in 770 whole-genome sequenced Arabidopsis thaliana, we show that length variation in glutamine and asparagine amino acid homopolymers, as well as in interaction hotspots, correlate with local bioclimatic habitat. We determined experimentally that the promoter activity of a light-stress gene depended on polyglutamine length variants in a disordered transcription factor. Our results show that length variations affect protein function and are likely adaptive. Length variants modulating protein function at a global genomic scale has implications for understanding protein evolution and eco-evolutionary biology.

Atlantic cod (Gadus morhua) MHC I localizes to endolysosomal compartments independently of cytosolic sorting signals.

Major histocompatibility complex (MHC) class I and II are crucial for the adaptive immune system because they are involved in peptide presentation to T cells. Until recently, it was believed that MHC genes and their associated immune components had been conserved since their evolutionary emergence in jawed fish. However, sequencing of the Atlantic cod (Gadus morhua) genome revealed a loss of MHC class II genes, and an extreme expansion of MHC class I genes. These findings lead to the hypothesis that a loss of the MHC class II pathway coincided with a more versatile use of MHC class I, but so far there is no direct experimental evidence in support of this. To gain a deeper understanding of the function of the expanded MHC class I, we selected five MHC class I gene variants representing five of the six clades identified in previous studies and investigated their intracellular localization in human and Atlantic cod larval cells. Intriguingly, we uncovered that all selected MHC class I variants localize to endolysosomal compartments in Atlantic cod cells. Additionally, by introducing point mutations or deletions in the cytosolic tail, we found that hypothetical sorting signals in the MHC class I cytosolic tail do not influence MHC class I trafficking. Moreover, we demonstrated that in Atlantic cod, tapasin and MHC class I colocalize on endolysosomes suggesting that peptide-loading assistance and stabilization of MHC class I occurs outside the endoplasmic reticulum. Altogether, our results demonstrate that MHC class I from Atlantic cod is sorted to the endolysosomal system, which may indicate that it interacts with exogenous peptides for potential cross presentation.

Ancient DNA evidence for the ecological globalization of cod fishing in medieval and post-medieval Europe.

Understanding the historical emergence and growth of long-range fisheries can provide fundamental insights into the timing of ecological impacts and the development of coastal communities during the last millennium. Whole-genome sequencing approaches can improve such understanding by determining the origin of archaeological fish specimens that may have been obtained from historic trade or distant water. Here, we used genome-wide data to individually infer the biological source of 37 ancient Atlantic cod specimens (ca 1050-1950 CE) from England and Spain. Our findings provide novel genetic evidence that eleventh- to twelfth-century specimens from London were predominantly obtained from nearby populations, while thirteenth- to fourteenth-century specimens were derived from distant sources. Our results further suggest that Icelandic cod was indeed exported to London earlier than previously reported. Our observations confirm the chronology and geography of the trans-Atlantic cod trade from Newfoundland to Spain starting by the early sixteenth century. Our findings demonstrate the utility of whole-genome sequencing and ancient DNA approaches to describe the globalization of marine fisheries and increase our understanding regarding the extent of the North Atlantic fish trade and long-range fisheries in medieval and early modern times.

Salinity-induced transcriptome profiles in marine and freshwater threespine stickleback after an abrupt 6-hour exposure.

Saltwater and freshwater environments have opposing physiological challenges, yet, there are fish species that are able to enter both habitats during short time spans, and as individuals they must therefore adjust quickly to osmoregulatory contrasts. In this study, we conducted an experiment to test for plastic responses to abrupt salinity changes in two populations of threespine stickleback, Gasterosteus aculeatus, representing two ecotypes (freshwater and ancestral saltwater). We exposed both ecotypes to abrupt native (control treatment) and non-native salinities (0‰ and 30‰) and sampled gill tissue for transcriptomic analyses after 6 h of exposure. To investigate genomic responses to salinity, we analyzed four different comparisons; one for each ecotype (in their control and exposure salinity; (1) and (2), one between ecotypes in their control salinity (3), and the fourth comparison included all transcripts identified in (3) that did not show any expressional changes within ecotype in either the control or the exposed salinity (4)). Abrupt salinity transfer affected the expression of 10 and 1530 transcripts for the saltwater and freshwater ecotype, respectively, and 1314 were differentially expressed between the controls, including 502 that were not affected by salinity within ecotype (fixed expression). In total, these results indicate that factors other than genomic expressional plasticity are important for osmoregulation in stickleback, due to the need for opposite physiological pathways to survive the abrupt change in salinity.

Lymphocyte subsets in Atlantic cod (Gadus morhua) interrogated by single-cell sequencing.

Atlantic Cod (Gadus morhua) has lost the major histocompatibility complex class II presentation pathway. We recently identified CD8-positive T cells, B cells, and plasma cells in cod, but further characterisation of lymphocyte subsets is needed to elucidate immune adaptations triggered by the absence of CD4-positive T lymphocytes. Here, we use single-cell RNA sequencing to examine the lymphocyte heterogeneity in Atlantic cod spleen. We describe five T cell subsets and eight B cell subsets and propose a B cell trajectory of differentiation. Notably, we identify a subpopulation of T cells that are CD8-negative. Most of the CD8-negative T lymphocytes highly express the homologue of monocyte chemotactic protein 1b, and another subset of CD8-negative T lymphocytes express the homologue of the scavenger receptor m130. Uncovering the multiple lymphocyte cell sub-clusters reveals the different immune states present within the B and T cell populations, building a foundation for further work.

Ancient DNA reveals a southern presence of the Northeast Arctic cod during the Holocene.

Climate change has been implicated in an increased number of distributional shifts of marine species during the last century. Nonetheless, it is unclear whether earlier climatic fluctuations had similar impacts. We use ancient DNA to investigate the long-term spawning distribution of the Northeast Arctic cod (skrei) which performs yearly migrations from the Barents Sea towards spawning grounds along the Norwegian coast. The distribution of these spawning grounds has shifted northwards during the last century, which is thought to be associated with food availability and warming temperatures. We genetically identify skrei specimens from Ruskeneset in west Norway, an archaeological site located south of their current spawning range. Remarkably, 14C analyses date these specimens to the late Holocene, when temperatures were warmer than present-day conditions. Our results either suggest that temperature is not the only driver influencing the spawning distribution of Atlantic cod, or could be indicative of uncertainty in palaeoclimate reconstructions in this region. Regardless, our findings highlight the utility of aDNA to reconstruct the historical distribution of economically important fish populations and reveal the complexity of long-term ecological interactions in the marine environment.

An accurate assignment test for extremely low-coverage whole-genome sequence data.

Genomic assignment tests can provide important diagnostic biological characteristics, such as population of origin or ecotype. Yet, assignment tests often rely on moderate- to high-coverage sequence data that can be difficult to obtain for fields such as molecular ecology and ancient DNA. We have developed a novel approach that efficiently assigns biologically relevant information (i.e., population identity or structural variants such as inversions) in extremely low-coverage sequence data. First, we generate databases from existing reference data using a subset of diagnostic single nucleotide polymorphisms (SNPs) associated with a biological characteristic. Low-coverage alignment files are subsequently compared to these databases to ascertain allelic state, yielding a joint probability for each association. To assess the efficacy of this approach, we assigned haplotypes and population identity in Heliconius butterflies, Atlantic herring, and Atlantic cod using chromosomal inversion sites and whole-genome data. We scored both modern and ancient specimens, including the first whole-genome sequence data recovered from ancient Atlantic herring bones. The method accurately assigns biological characteristics, including population membership, using extremely low-coverage data (as low as 0.0001x) based on genome-wide SNPs. This approach will therefore increase the number of samples in evolutionary, ecological and archaeological research for which relevant biological information can be obtained.

Polygenic plague resistance in the great gerbil uncovered by population sequencing.

Pathogens can elicit high selective pressure on hosts, potentially altering genetic diversity over short evolutionary timescales. Intraspecific variation in immune response is observable as variable survivability from specific infections. The great gerbil (Rhombomys opimus) is a rodent plague host with a heterogenic but highly resistant phenotype. Here, we investigate the genomic basis for plague-resistant phenotypes by exposing wild-caught great gerbils to plague (Yersinia pestis). Whole genome sequencing of 10 survivors and 10 moribund individuals revealed a subset of genomic regions showing elevated differentiation. Gene ontology analysis of candidate genes in these regions demonstrated enrichment of genes directly involved in immune functions, cellular metabolism and the regulation of apoptosis as well as pathways involved in transcription, translation, and gene regulation. Transcriptomic analysis revealed that the early activated great gerbil immune response to plague consisted of classical components of the innate immune system. Our approach combining challenge experiments with transcriptomics and population level sequencing, provides new insight into the genetic background of plague-resistance and confirms its complex nature, most likely involving multiple genes and pathways of both the immune system and regulation of basic cellular functions.

A Chromosome-Level Genome Assembly of the Reed Warbler (Acrocephalus scirpaceus).

The reed warbler (Acrocephalus scirpaceus) is a long-distance migrant passerine with a wide distribution across Eurasia. This species has fascinated researchers for decades, especially its role as host of a brood parasite, and its capacity for rapid phenotypic change in the face of climate change. Currently, it is expanding its range northwards in Europe, and is altering its migratory behavior in certain areas. Thus, there is great potential to discover signs of recent evolution and its impact on the genomic composition of the reed warbler. Here, we present a high-quality reference genome for the reed warbler, based on PacBio, 10×, and Hi-C sequencing. The genome has an assembly size of 1,075,083,815 bp with a scaffold N50 of 74,438,198 bp and a contig N50 of 12,742,779 bp. BUSCO analysis using aves_odb10 as a model showed that 95.7% of BUSCO genes were complete. We found unequivocal evidence of two separate macrochromosomal fusions in the reed warbler genome, in addition to the previously identified fusion between chromosome Z and a part of chromosome 4A in the Sylvioidea superfamily. We annotated 14,645 protein-coding genes, and a BUSCO analysis of the protein sequences indicated 97.5% completeness. This reference genome will serve as an important resource, and will provide new insights into the genomic effects of evolutionary drivers such as coevolution, range expansion, and adaptations to climate change, as well as chromosomal rearrangements in birds.

Complex population structure of the Atlantic puffin revealed by whole genome analyses.

The factors underlying gene flow and genomic population structure in vagile seabirds are notoriously difficult to understand due to their complex ecology with diverse dispersal barriers and extensive periods at sea. Yet, such understanding is vital for conservation management of seabirds that are globally declining at alarming rates. Here, we elucidate the population structure of the Atlantic puffin (Fratercula arctica) by assembling its reference genome and analyzing genome-wide resequencing data of 72 individuals from 12 colonies. We identify four large, genetically distinct clusters, observe isolation-by-distance between colonies within these clusters, and obtain evidence for a secondary contact zone. These observations disagree with the current taxonomy, and show that a complex set of contemporary biotic factors impede gene flow over different spatial scales. Our results highlight the power of whole genome data to reveal unexpected population structure in vagile marine seabirds and its value for seabird taxonomy, evolution and conservation.

Length variation in short tandem repeats affects gene expression in natural populations of Arabidopsis thaliana.

The genetic basis for the fine-tuned regulation of gene expression is complex and ultimately influences the phenotype and thus the local adaptation of natural populations. Short tandem repeats (STRs) consisting of repetitive DNA motifs have been shown to regulate gene expression. STRs are variable in length within a population and serve as a heritable, but semi-reversible, reservoir of standing genetic variation. For sessile organisms, such as plants, STRs could be of major importance in fine-tuning gene expression as a response to a shifting local environment. Here, we used a transcriptome dataset from natural accessions of Arabidopsis thaliana to investigate population-wide gene expression patterns in light of genome-wide STR variation. We empirically modeled gene expression as a response to the STR length within and around the gene and demonstrated that an association between gene expression and STR length variation is unequivocally present in the sampled population. To support our model, we explored the promoter activity in a transcriptional regulator involved in root hair formation and provided experimentally determined causality between coding sequence length variation and promoter activity. Our results support a general link between gene expression variation and STR length variation in A. thaliana.

Genomic stability through time despite decades of exploitation in cod on both sides of the Atlantic.

The mode and extent of rapid evolution and genomic change in response to human harvesting are key conservation issues. Although experiments and models have shown a high potential for both genetic and phenotypic change in response to fishing, empirical examples of genetic responses in wild populations are rare. Here, we compare whole-genome sequence data of Atlantic cod (Gadus morhua) that were collected before (early 20th century) and after (early 21st century) periods of intensive exploitation and rapid decline in the age of maturation from two geographically distinct populations in Newfoundland, Canada, and the northeast Arctic, Norway. Our temporal, genome-wide analyses of 346,290 loci show no substantial loss of genetic diversity and high effective population sizes. Moreover, we do not find distinct signals of strong selective sweeps anywhere in the genome, although we cannot rule out the possibility of highly polygenic evolution. Our observations suggest that phenotypic change in these populations is not constrained by irreversible loss of genomic variation and thus imply that former traits could be reestablished with demographic recovery.

Innovation in Nucleotide-Binding Oligomerization-Like Receptor and Toll-Like Receptor Sensing Drives the Major Histocompatibility Complex-II Free Atlantic Cod Immune System.

The absence of MHC class II antigen presentation and multiple pathogen recognition receptors in the Atlantic cod has not impaired its immune response however how underlying mechanisms have adapted remains largely unknown. In this study, ex vivo cod macrophages were challenged with various bacterial and viral microbe-associated molecular patterns (MAMP) to identify major response pathways. Cytosolic MAMP-PRR pathways based upon the NOD-like receptors (NLRs) and RIG-I-like receptors (RLRs) were identified as the critical response pathways. Our analyses suggest that internalization of exogenous ligands through scavenger receptors drives both pathways activating transcription factors like NF-kB (Nuclear factor-kappa B) and interferon regulatory factors (IRFs). Further, ligand-dependent differential expression of a unique TLR25 isoform and multiple NLR paralogues suggests (sub)neofunctionalization toward specific immune defensive strategies. Our results further demonstrate that the unique immune system of the Atlantic cod provides an unprecedented opportunity to explore the evolutionary history of PRR-based signaling in vertebrate immunity.

Single-Cell Transcriptome Profiling of Immune Cell Repertoire of the Atlantic Cod Which Naturally Lacks the Major Histocompatibility Class II System.

The Atlantic cod's unusual immune system, entirely lacking the Major Histocompatibility class II pathway, has prompted intriguing questions about what mechanisms are used to combat bacterial infections and how immunological memory is generated. By single-cell RNA sequencing we here report an in-depth characterisation of cell types found in immune tissues, the spleen and peripheral blood leukocytes of Atlantic cod. Unbiased transcriptional clustering revealed eleven distinct immune cell signatures. Resolution at the single cell level enabled characterisation of the major cell subsets including the cytotoxic T cells, B cells, erythrocytes, thrombocytes, neutrophils, and macrophages. Additionally, to our knowledge we are the first to uncover cell subsets in Atlantic cod which may represent dendritic cells, natural killer-like cells, and a population of cytotoxic cells expressing GATA-3, a master transcription factor of T helper 2 cells. We further identify putative gene markers for each cluster and describe the relative proportions of each cell type in the spleen and peripheral blood leukocytes. Of the major haematopoietic cell populations, the lymphocytes make up 55 and 68% of the spleen and peripheral blood leukocytes respectively, while the myeloid cells make up 45 and 32%. By single-cell analysis, this study provides the most detailed molecular and cellular characterisation of the immune system of the Atlantic cod so far.

The Genome of the Great Gerbil Reveals Species-Specific Duplication of an MHCII Gene.

The great gerbil (Rhombomys opimus) is a social rodent living in permanent, complex burrow systems distributed throughout Central Asia, where it serves as the main host of several important vector-borne infectious pathogens including the well-known plague bacterium (Yersinia pestis). Here, we present a continuous annotated genome assembly of the great gerbil, covering over 96% of the estimated 2.47-Gb genome. Taking advantage of the recent genome assemblies of the sand rat (Psammomys obesus) and the Mongolian gerbil (Meriones unguiculatus), comparative immunogenomic analyses reveal shared gene losses within TLR gene families (i.e., TLR8, TLR10, and the entire TLR11-subfamily) for Gerbillinae, accompanied with signs of diversifying selection of TLR7 and TLR9. Most notably, we find a great gerbil-specific duplication of the MHCII DRB locus. In silico analyses suggest that the duplicated gene provides high peptide binding affinity for Yersiniae epitopes as well as Leishmania and Leptospira epitopes, putatively leading to increased capability to withstand infections by these pathogens. Our study demonstrates the power of whole-genome sequencing combined with comparative genomic analyses to gain deeper insight into the immunogenomic landscape of the great gerbil and its close relatives.

Evolutionary selection of biofilm-mediated extended phenotypes in Yersinia pestis in response to a fluctuating environment.

Yersinia pestis is transmitted from fleas to rodents when the bacterium develops an extensive biofilm in the foregut of a flea, starving it into a feeding frenzy, or, alternatively, during a brief period directly after feeding on a bacteremic host. These two transmission modes are in a trade-off regulated by the amount of biofilm produced by the bacterium. Here by investigating 446 global isolated Y. pestis genomes, including 78 newly sequenced isolates sampled over 40 years from a plague focus in China, we provide evidence for strong selection pressures on the RNA polymerase ω-subunit encoding gene rpoZ. We demonstrate that rpoZ variants have an increased rate of biofilm production in vitro, and that they evolve in the ecosystem during colder and drier periods. Our results support the notion that the bacterium is constantly adapting-through extended phenotype changes in the fleas-in response to climate-driven changes in the niche.