Search

Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%–40%. 1 It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis). 2 Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17th–19th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large F ST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics. • The second plague pandemic homogenized ancestry in Trondheim • Gaelic ancestry is sharply reduced in post-pandemic Trondheim • Pervasive reference bias taints frequency differences observed between populations Gopalakrishnan et al. investigate the genomic signatures of the second plague pandemic on the residents of Trondheim in Norway. They find that the pandemic resulted in a sharp reduction in Gaelic ancestry and also find evidence of differential reference bias among their ancient samples, which reduces the reliability of selection analyses. [ABSTRACT FROM AUTHOR]

Current mitochondrial DNA (mtDNA) haplogroup classification tools map reads to a single reference genome and perform inference based on the detected mutations to this reference. This approach biases haplogroup assignments towards the reference and prohibits accurate calculations of the uncertainty in assignment. We present HaploCart, a probabilistic mtDNA haplogroup classifier which uses a pangenomic reference graph framework together with principles of Bayesian inference. We demonstrate that our approach significantly outperforms available tools by being more robust to lower coverage or incomplete consensus sequences and producing phylogenetically-aware confidence scores that are unbiased towards any haplogroup. HaploCart is available both as a command-line tool and through a user-friendly web interface. The C++ program accepts as input consensus FASTA, FASTQ, or GAM files, and outputs a text file with the haplogroup assignments of the samples along with the level of confidence in the assignments. Our work considerably reduces the amount of data required to obtain a confident mitochondrial haplogroup assignment. Author summary: Pangenome graphs are powerful and relatively nascent data structures for representing an entire collection of genomic sequences and their homology. Here we present HaploCart, a tool which leverages the power of pangenomics, in conjunction with maximum-likelihood estimation, to improve human mtDNA haplotype inference on single-source samples (i.e. the sample is not a mixture of multiple contributors, be they human or contaminant). In this context, mapping to many reference genomes at once vastly reduces the Eurocentric bias inherent in contemporary methods, and also improves haplotyping performance at low coverage depths. We show that HaploCart is far more accurate than competing programs on simulated and empirical datasets, and reports clade-level posterior probabilities that accurately reflect confidence in our phylogenetic assignments. Our work can easily be generalized to other haploid markers and suggests that pangenome-based approaches combined with Bayesian methods show promise for improving inference and mitigating ethnicity-related bias in a large class of bioinformatics problems involving sequencing data. [ABSTRACT FROM AUTHOR]

ABSTRACT The recently extirpated Zanzibar leopard was the only known African leopard (Panthera pardus spp.) population restricted exclusively to a major island habitat. Although its demise was driven through habitat change and conflict with humans, given its role as a keystone species for the Zanzibar Archipelago, its successful potential reintroduction might offer a means for helping preserve the natural biodiversity of its former habitat. Whether this is feasible, however, would be contingent on both whether closely related source populations can be identified on mainland Africa, and whether the Zanzibar form exhibited any special adaptations that might need to be considered when choosing such a source. In light of these questions, we genomically profiled two of the six known historic specimens, to explore whether they represent a realistic candidate for de‐extirpation through reintroduction. Our analyses indicate that despite its geographical separation, the Zanzibar leopard shared a close genetic relationship with mainland East African individuals. Furthermore, although its uniqueness as an island population was emphasised by genomic signatures of high inbreeding and increased mutation load, the latter similar to the level of the critically endangered Amur leopard (P. p. orientalis), we find no evidence of functionally significant genetic diversity unique to Zanzibar. We therefore conclude that should attempts to restore leopards to Zanzibar be considered, then mainland East African leopards would provide a suitable gene pool. [ABSTRACT FROM AUTHOR]

The Seychelles magpie‐robin's (SMR) five island populations exhibit some of the lowest recorded levels of genetic diversity among endangered birds, and high levels of inbreeding. These populations collapsed during the 20th century, and the species was listed as Critically Endangered in the IUCN Red List in 1994. An assisted translocation‐for‐recovery program initiated in the 1990s increased the number of mature individuals, resulting in its downlisting to Endangered in 2005. Here, we explore the temporal genomic erosion of the SMR based on a dataset of 201 re‐sequenced whole genomes that span the past ~150 years. Our sample set includes individuals that predate the bottleneck by up to 100 years, as well as individuals from contemporary populations established during the species recovery program. Despite the SMR's recent demographic recovery, our data reveal a marked increase in both the genetic load and realized load in the extant populations when compared to the historical samples. Conservation management may have reduced the intensity of selection by increasing juvenile survival and relaxing intraspecific competition between individuals, resulting in the accumulation of loss‐of‐function mutations (i.e. severely deleterious variants) in the rapidly recovering population. In addition, we found a 3‐fold decrease in genetic diversity between temporal samples. While the low genetic diversity in modern populations may limit the species' adaptability to future environmental changes, future conservation efforts (including IUCN assessments) may also need to assess the threats posed by their high genetic load. Our computer simulations highlight the value of translocations for genetic rescue and show how this could halt genomic erosion in threatened species such as the SMR. [ABSTRACT FROM AUTHOR]

Ancestral Coast Salish societies in the Pacific Northwest kept long-haired “woolly dogs” that were bred and cared for over millennia. However, the dog wool–weaving tradition declined during the 19th century, and the population was lost. In this study, we analyzed genomic and isotopic data from a preserved woolly dog pelt from “Mutton,” collected in 1859. Mutton is the only known example of an Indigenous North American dog with dominant precolonial ancestry postdating the onset of settler colonialism. We identified candidate genetic variants potentially linked with their distinct woolly phenotype. We integrated these data with interviews from Coast Salish Elders, Knowledge Keepers, and weavers about shared traditional knowledge and memories surrounding woolly dogs, their importance within Coast Salish societies, and how colonial policies led directly to their disappearance. [ABSTRACT FROM AUTHOR]

The island of St Helena played a crucial role in the suppression of the transatlantic slave trade. Strategically located in the middle of the South Atlantic, it served as a staging post for the Royal Navy and reception point for enslaved Africans who had been "liberated" from slave ships intercepted by the British. In total, St Helena received approximately 27,000 liberated Africans between 1840 and 1867. Written sources suggest that the majority of these individuals came from West Central Africa, but their precise origins are unknown. Here, we report the results of ancient DNA analyses that we conducted as part of a wider effort to commemorate St Helena's liberated Africans and to restore knowledge of their lives and experiences. We generated partial genomes (0.1–0.5×) for 20 individuals whose remains had been recovered during archaeological excavations on the island. We compared their genomes with genotype data for over 3,000 present-day individuals from 90 populations across sub-Saharan Africa and conclude that the individuals most likely originated from different source populations within the general area between northern Angola and Gabon. We also find that the majority (17/20) of the individuals were male, supporting a well-documented sex bias in the latter phase of the transatlantic slave trade. The study expands our understanding of St Helena's liberated African community and illustrates how ancient DNA analyses can be used to investigate the origins and identities of individuals whose lives were bound up in the story of slavery and its abolition. During the mid-19th century, around 27,000 "liberated Africans" were brought to St Helena island as part of the Royal Navy's efforts to suppress the transatlantic trade in enslaved Africans. Ancient DNA analyses, carried out on behalf of the community on St Helena, now shed light on their origins and identities. [ABSTRACT FROM AUTHOR]

Motivation Machine learning (ML) methods are motivated by the need to automate information extraction from large datasets in order to support human users in data-driven tasks. This is an attractive approach for integrative joint analysis of vast amounts of omics data produced in next generation sequencing and other -omics assays. A systematic assessment of the current literature can help to identify key trends and potential gaps in methodology and applications. We surveyed the literature on ML multi-omic data integration and quantitatively explored the goals, techniques and data involved in this field. We were particularly interested in examining how researchers use ML to deal with the volume and complexity of these datasets. Results Our main finding is that the methods used are those that address the challenges of datasets with few samples and many features. Dimensionality reduction methods are used to reduce the feature count alongside models that can also appropriately handle relatively few samples. Popular techniques include autoencoders, random forests and support vector machines. We also found that the field is heavily influenced by the use of The Cancer Genome Atlas dataset, which is accessible and contains many diverse experiments. Availability and implementation All data and processing scripts are available at this GitLab repository: https://gitlab.com/polavieja%5flab/ml%5fmulti-omics%5freview/ or in Zenodo: https://doi.org/10.5281/zenodo.7361807. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

Environmental variation across the range of wild species can lead to local adaptations. The Baltic Sea was formed when the Fenno-Scandian ice sheet retreated around 12 thousand years ago, creating a new brackish water habitat colonised by both marine and freshwater fish species. The European perch (Perca fluviatilis) is a predatory freshwater fish with a large geographical distribution across Eurasia, where it inhabits a wide range of environmental niches. In the Baltic Sea region it has even developed a specialised brackish water perch variant that can tolerate environmental salinity levels, which are lethal to freshwater perch. However, very little is known about the underlying mechanisms facilitating the colonisation and adaptation of perch to the Baltic Sea. Here, we use Genotyping-By-Sequencing data from six freshwater and six brackish water localities to disclose the evolutionary relationship between freshwater and brackish water perch. Our results show that the brackish water perch occurs in multiple distinct genetic clusters. We find that gene flow between brackish water perch with full access to the sea likely led to lower levels of differentiation and higher diversity than in freshwater perch. Selection analyses suggest that genomic adaptation played a role in the colonisation of the Baltic Sea and that the top three regions under selection harbour salinity tolerance genes. We complete by discussing the implications of our findings for management of brackish water perch in the western Baltic sea. [Display omitted] • GBS data from 12 perch populations, six with brackish and six with freshwater origin. • Brackish water populations were found in all three major genetic clades. • Top three regions under selection harboured salinity tolerance genes. [ABSTRACT FROM AUTHOR]

Preserving genetic diversity and adaptive potential while avoiding inbreeding depression is crucial for the long-term conservation of natural populations. Despite demographic increases, traces of past bottleneck events at the genomic level should be carefully considered for population management. From this perspective, the peninsular Italian wolf is a paradigmatic case. After being on the brink of extinction in the late 1960s, peninsular Italian wolves rebounded and recolonized most of the peninsula aided by conservation measures, including habitat and legal protection. Notwithstanding their demographic recovery, a comprehensive understanding of the genomic consequences of the historical bottleneck in Italian wolves is still lacking. To fill this gap, we sequenced whole genomes of thirteen individuals sampled in the core historical range of the species in Central Italy to conduct population genomic analyses, including a comparison with wolves from two highly-inbred wolf populations (i.e., Scandinavia and Isle Royale). We found that peninsular Italian wolves, despite their recent recovery, still exhibit relatively low genetic diversity, a small effective population size, signatures of inbreeding, and a non-negligible genetic load. Our findings indicate that the peninsular Italian wolf population is still susceptible to bottleneck legacies, which could lead to local inbreeding depression in case of population reduction or fragmentations. This study emphasizes the importance of considering key genetic parameters to design appropriate long-term conservation management plans., (© The American Genetic Association. 2024.)

The multi-millenia long history between dogs and humans has placed them at the forefront of archeological and genomic research. Despite ongoing efforts including the analysis of ancient dog and wolf genomes, many questions remain regarding their geographic and temporal origins, and the microevolutionary processes that led to the diversity of breeds today. Although ancient genomes provide valuable information, their use is hindered by low depth of coverage and post-mortem damage, which inhibits confident genotype calling. In the present study, we assess how genotype imputation of ancient dog and wolf genomes, utilising a large reference panel, can improve the resolution provided by ancient datasets. Imputation accuracy was evaluated by down-sampling high coverage dog and wolf genomes to 0.05-2x coverage and comparing concordance between imputed and high coverage genotypes. We measured the impact of imputation on principal component analyses and runs of homozygosity. Our findings show high (R 2 >0.9) imputation accuracy for dogs with coverage as low as 0.5x and for wolves as low as 1.0x. We then imputed a dataset of 90 ancient dog and wolf genomes, to assess changes in inbreeding during the last 10,000 years of dog evolution. Ancient dog and wolf populations generally exhibited lower inbreeding levels than present-day individuals. Interestingly, regions with low ROH density maintained across ancient and present-day samples were significantly associated with genes related to olfaction and immune response. Our study indicates that imputing ancient canine genomes is a viable strategy that allows for the use of analytical methods previously limited to high-quality genetic data.

The domestic pigeon's exceptional phenotypic diversity was key in developing Darwin's Theory of Evolution and establishing the concept of artificial selection. However, unlike its domestic counterpart, its wild progenitor, the rock dove Columba livia has received considerably less attention. Therefore, questions regarding its domestication, evolution, taxonomy, and conservation status remain unresolved. We generated whole-genome sequencing data from 65 historical rock doves that represent all currently recognized subspecies and span the species' original geographic distribution. Our dataset includes 3 specimens from Darwin's collection, and the type specimens of 5 different taxa. We characterized their population structure, genomic diversity, and gene-flow patterns. Our results show the West African subspecies C. l. gymnocyclus is basal to rock doves and domestic pigeons, and suggests gene-flow between the rock dove's sister species C. rupestris, and the ancestor of rock doves after its split from West African populations. These genomes allowed us to propose a model for the evolution of the rock dove in light of the refugia theory. We propose that rock dove genetic diversity and introgression patterns derive from a history of allopatric cycles and dispersion waves during the Quaternary glacial and interglacial periods. To explore the rock dove domestication history, we combined our new dataset with available genomes from domestic pigeons. Our results point to at least 1 domestication event in the Levant that gave rise to all domestic breeds analysed in this study. Finally, we propose a species-level taxonomic arrangement to reflect the evolutionary history of the West African rock dove populations., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The tiger (Panthera tigris) is a charismatic megafauna species that originated and diversified in Asia and probably experienced population contraction and expansion during the Pleistocene, resulting in low genetic diversity of modern tigers. However, little is known about patterns of genomic diversity in ancient populations. Here we generated whole-genome sequences from ancient or historical (100-10,000 yr old) specimens collected across mainland Asia, including a 10,600-yr-old Russian Far East specimen (RUSA21, 8× coverage) plus six ancient mitogenomes, 14 South China tigers (0.1-12×) and three Caspian tigers (4-8×). Admixture analysis showed that RUSA21 clustered within modern Northeast Asian phylogroups and partially derived from an extinct Late Pleistocene lineage. While some of the 8,000-10,000-yr-old Russian Far East mitogenomes are basal to all tigers, one 2,000-yr-old specimen resembles present Amur tigers. Phylogenomic analyses suggested that the Caspian tiger probably dispersed from an ancestral Northeast Asian population and experienced gene flow from southern Bengal tigers. Lastly, genome-wide monophyly supported the South China tiger as a distinct subspecies, albeit with mitochondrial paraphyly, hence resolving its longstanding taxonomic controversy. The distribution of mitochondrial haplogroups corroborated by biogeographical modelling suggested that Southwest China was a Late Pleistocene refugium for a relic basal lineage. As suitable habitat returned, admixture between divergent lineages of South China tigers took place in Eastern China, promoting the evolution of other northern subspecies. Altogether, our analysis of ancient genomes sheds light on the evolutionary history of tigers and supports the existence of nine modern subspecies., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

The taxonomic status of the now likely extirpated Korean Peninsula wolf has been extensively debated, with some arguing it represents an independent wolf lineage, Canis coreanus . To investigate the Korean wolf's genetic affiliations and taxonomic status, we sequenced and analysed the genomes of a Korean wolf dated to the beginning of the 20th century, and a captive wolf originally from the Pyongyang Central Zoo. Our results indicated that the Korean wolf bears similar genetic ancestry to other regional East Asian populations, therefore suggesting it is not a distinct taxonomic lineage. We identified regional patterns of wolf population structure and admixture in East Asia with potential conservation consequences in the Korean Peninsula and on a regional scale. We find that the Korean wolf has similar genomic diversity and inbreeding to other East Asian wolves. Finally, we show that, in contrast to the historical sample, the captive wolf is genetically more similar to wolves from the Tibetan Plateau; hence, Korean wolf conservation programmes might not benefit from the inclusion of this specimen., (© 2023 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

The Sicilian wolf remained isolated in Sicily from the end of the Pleistocene until its extermination in the 1930s-1960s. Given its long-term isolation on the island and distinctive morphology, the genetic origin of the Sicilian wolf remains debated. We sequenced four nuclear genomes and five mitogenomes from the seven existing museum specimens to investigate the Sicilian wolf ancestry, relationships with extant and extinct wolves and dogs, and diversity. Our results show that the Sicilian wolf is most closely related to the Italian wolf but carries ancestry from a lineage related to European Eneolithic and Bronze Age dogs. The average nucleotide diversity of the Sicilian wolf was half of the Italian wolf, with 37-50% of its genome contained in runs of homozygosity. Overall, we show that, by the time it went extinct, the Sicilian wolf had high inbreeding and low-genetic diversity, consistent with a population in an insular environment., Competing Interests: The authors declare no competing interests., (© 2023 The Authors.)

We investigate a 2,000-year genetic transect through Scandinavia spanning the Iron Age to the present, based on 48 new and 249 published ancient genomes and genotypes from 16,638 modern individuals. We find regional variation in the timing and magnitude of gene flow from three sources: the eastern Baltic, the British-Irish Isles, and southern Europe. British-Irish ancestry was widespread in Scandinavia from the Viking period, whereas eastern Baltic ancestry is more localized to Gotland and central Sweden. In some regions, a drop in current levels of external ancestry suggests that ancient immigrants contributed proportionately less to the modern Scandinavian gene pool than indicated by the ancestry of genomes from the Viking and Medieval periods. Finally, we show that a north-south genetic cline that characterizes modern Scandinavians is mainly due to the differential levels of Uralic ancestry and that this cline existed in the Viking Age and possibly earlier., Competing Interests: Declaration of interests A.H., K.H.S.M., K.S., and S.S.E. are employees of deCODE genetics., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

Identifying the specific human leukocyte antigen (HLA) allele combination of an individual is crucial in organ donation, risk assessment of autoimmune and infectious diseases and cancer immunotherapy. However, due to the high genetic polymorphism in this region, HLA typing requires specialized methods. We investigated the performance of five next-generation sequencing (NGS) based HLA typing tools with a non-restricted license namely HLA*LA, Optitype, HISAT-genotype, Kourami and STC-Seq. This evaluation was done for the five HLA loci, HLA-A, -B, -C, -DRB1 and -DQB1 using whole-exome sequencing (WES) samples from 829 individuals. The robustness of the tools to lower depth of coverage (DOC) was evaluated by subsampling and HLA typing 230 WES samples at DOC ranging from 1X to 100X. The HLA typing accuracy was measured across four typing resolutions. Among these, we present two clinically-relevant typing resolutions (P group and pseudo-sequence), which specifically focus on the peptide binding region. On average, across the five HLA loci examined, HLA*LA was found to have the highest typing accuracy. For the individual loci, HLA-A, -B and -C, Optitype's typing accuracy was the highest and HLA*LA had the highest typing accuracy for HLA-DRB1 and -DQB1. The tools' robustness to lower DOC data varied widely and further depended on the specific HLA locus. For all Class I loci, Optitype had a typing accuracy above 95% (according to the modification of the amino acids in the functionally relevant portion of the HLA molecule) at 50X, but increasing the DOC beyond even 100X could still improve the typing accuracy of HISAT-genotype, Kourami, and STC-seq across all five HLA loci as well as HLA*LA's typing accuracy for HLA-DQB1. HLA typing is also used in studies of ancient DNA (aDNA), which is often based on sequencing data with lower quality and DOC. Interestingly, we found that Optitype's typing accuracy is not notably impaired by short read length or by DNA damage, which is typical of aDNA, as long as the DOC is sufficiently high., Competing Interests: NT, MK and TT are employed by Evaxion Biotech and have a financial stake in the company. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Thuesen, Klausen, Gopalakrishnan, Trolle and Renaud.)

The grey wolf (Canis lupus) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage (Canis familiaris) lived 1-8 . Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000-30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located., (© 2022. The Author(s).)

Retracing pathways of historical species introductions is fundamental to understanding the factors involved in the successful colonization and spread, centuries after a species' establishment in an introduced range. Numerous plants have been introduced to regions outside their native ranges both intentionally and accidentally by European voyagers and early colonists making transoceanic journeys; however, records are scarce to document this. We use genotyping-by-sequencing and genotype-likelihood methods on the selfing, global weed, Plantago major , collected from 50 populations worldwide to investigate how patterns of genomic diversity are distributed among populations of this global weed. Although genomic differentiation among populations is found to be low, we identify six unique genotype groups showing very little sign of admixture and low degree of outcrossing among them. We show that genotype groups are latitudinally restricted, and that more than one successful genotype colonized and spread into the introduced ranges. With the exception of New Zealand, only one genotype group is present in the Southern Hemisphere. Three of the most prevalent genotypes present in the native Eurasian range gave rise to introduced populations in the Americas, Africa, Australia, and New Zealand, which could lend support to the hypothesis that P. major was unknowlingly dispersed by early European colonists. Dispersal of multiple successful genotypes is a likely reason for success. Genomic signatures and phylogeographic methods can provide new perspectives on the drivers behind the historic introductions and the successful colonization of introduced species, contributing to our understanding of the role of genomic variation for successful establishment of introduced taxa., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Iwanycki Ahlstrand, Gopalakrishnan, Vieira, Bieker, Meudt, Dunbar-Co, Rothfels, Martinez-Swatson, Maldonado, Hassemer, Shipunov, Bowers, Gardner, Xu, Ghorbani, Amano, Grace, Pringle, Bishop, Manzanilla, Cotrim, Blaney, Zubov, Choi, Yesil, Bennett, Vimolmangkang, El-Seedi, Staub, Li, Boldbaatar, Hislop, Caddy, Muasya, Saslis-Lagoudakis, Gilbert, Zerega and Rønsted.)

Background: Animal protein production is increasingly looking towards microbiome-associated services such as the design of new and better probiotic solutions to further improve gut health and production sustainability. Here, we investigate the functional effects of bacteria-based pro- and synbiotic feed additives on microbiome-associated functions in relation to growth performance in the commercially important rainbow trout (Oncorhynchus mykiss). We combine complementary insights from multiple omics datasets from gut content samples, including 16S bacterial profiling, whole metagenomes, and untargeted metabolomics, to investigate bacterial metagenome-assembled genomes (MAGs) and their molecular interactions with host metabolism., Results: Our findings reveal that (I) feed additives changed the microbiome and that rainbow trout reared with feed additives had a significantly reduced relative abundance of the salmonid related Candidatus Mycoplasma salmoninae in both the mid and distal gut content, (II) genome resolved metagenomics revealed that alterations of microbial arginine biosynthesis and terpenoid backbone synthesis pathways were directly associated with the presence of Candidatus Mycoplasma salmoninae, and (III) differences in the composition of intestinal microbiota among feed types were directly associated with significant changes of the metabolomic landscape, including lipids and lipid-like metabolites, amino acids, bile acids, and steroid-related metabolites., Conclusion: Our results demonstrate how the use of multi-omics to investigate complex host-microbiome interactions enable us to better evaluate the functional potential of probiotics compared to studies that only measure overall growth performance or that only characterise the microbial composition in intestinal environments. Video Abstract., (© 2022. The Author(s).)

Combining samples for genetic association is standard practice in human genetic analysis of complex traits, but is rarely undertaken in rodent genetics. Here, using 23 phenotypes and genotypes from two independent laboratories, we obtained a sample size of 3076 commercially available outbred mice and identified 70 loci, more than double the number of loci identified in the component studies. Fine-mapping in the combined sample reduced the number of likely causal variants, with a median reduction in set size of 51%, and indicated novel gene associations, including Pnpo, Ttll6, and GM11545 with bone mineral density, and Psmb9 with weight. However, replication at a nominal threshold of 0.05 between the two component studies was low, with less than one-third of loci identified in one study replicated in the second. In addition to overestimates in the effect size in the discovery sample (Winner's Curse), we also found that heterogeneity between studies explained the poor replication, but the contribution of these two factors varied among traits. Leveraging these observations, we integrated information about replication rates, study-specific heterogeneity, and Winner's Curse corrected estimates of power to assign variants to one of four confidence levels. Our approach addresses concerns about reproducibility and demonstrates how to obtain robust results from mapping complex traits in any genome-wide association study., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

We present a genome assembly from an individual male Canis lupus orion (the grey wolf, subspecies: Greenland wolf; Chordata; Mammalia; Carnivora; Canidae). The genome sequence is 2,447 megabases in span. The majority of the assembly (98.91%) is scaffolded into 40 chromosomal pseudomolecules, with the X and Y sex chromosomes assembled., Competing Interests: Competing interests: Jonathan Threlfall was employed by F1000 Research Ltd until January 2021., (Copyright: © 2021 Sinding MHS et al.)

Dogs have been essential to life in the Siberian Arctic for over 9,500 y, and this tight link between people and dogs continues in Siberian communities. Although Arctic Siberian groups such as the Nenets received limited gene flow from neighboring groups, archaeological evidence suggests that metallurgy and new subsistence strategies emerged in Northwest Siberia around 2,000 y ago. It is unclear if the Siberian Arctic dog population was as continuous as the people of the region or if instead admixture occurred, possibly in relation to the influx of material culture from other parts of Eurasia. To address this question, we sequenced and analyzed the genomes of 20 ancient and historical Siberian and Eurasian Steppe dogs. Our analyses indicate that while Siberian dogs were genetically homogenous between 9,500 to 7,000 y ago, later introduction of dogs from the Eurasian Steppe and Europe led to substantial admixture. This is clearly the case in the Iamal-Nenets region (Northwestern Siberia) where dogs from the Iron Age period (∼2,000 y ago) possess substantially less ancestry related to European and Steppe dogs than dogs from the medieval period (∼1,000 y ago). Combined with findings of nonlocal materials recovered from these archaeological sites, including glass beads and metal items, these results indicate that Northwest Siberian communities were connected to a larger trade network through which they acquired genetically distinctive dogs from other regions. These exchanges were part of a series of major societal changes, including the rise of large-scale reindeer pastoralism ∼800 y ago., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

The Japanese or Honshū wolf was one the most distinct gray wolf subspecies due to its small stature and endemicity to the islands of Honshū, Shikoku, and Kyūshū. Long revered as a guardian of farmers and travellers, it was persecuted from the 17th century following a rabies epidemic, which led to its extinction in the early 20th century. To better understand its evolutionary history, we sequenced the nuclear genome of a 19th century Honshū wolf specimen to an average depth of coverage of 3.7✕. We find Honshū wolves were closely related to a lineage of Siberian wolves that were previously believed to have gone extinct in the Late Pleistocene, thereby extending the survival of this ancient lineage until the early 20th century. We also detected significant gene flow between Japanese dogs and the Honshū wolf, corroborating previous reports on Honshū wolf dog interbreeding., Competing Interests: The authors declare no conflict of interest., (© 2020 The Author(s).)

From ontogenesis to homeostasis, the phenotypes of complex organisms are shaped by the bidirectional interactions between the host organisms and their associated microbiota. Current technology can reveal many such interactions by combining multi-omic data from both hosts and microbes. However, exploring the full extent of these interactions requires careful consideration of study design for the efficient generation and optimal integration of data derived from (meta)genomics, (meta)transcriptomics, (meta)proteomics, and (meta)metabolomics. In this perspective, we introduce the holo-omic approach that incorporates multi-omic data from both host and microbiota domains to untangle the interplay between the two. We revisit the recent literature on biomolecular host-microbe interactions and discuss the implementation and current limitations of the holo-omic approach. We anticipate that the application of this approach can contribute to opening new research avenues and discoveries in biomedicine, biotechnology, agricultural and aquacultural sciences, nature conservation, as well as basic ecological and evolutionary research., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Lions are one of the world's most iconic megafauna, yet little is known about their temporal and spatial demographic history and population differentiation. We analyzed a genomic dataset of 20 specimens: two ca. 30,000-y-old cave lions ( Panthera leo spelaea ), 12 historic lions ( Panthera leo leo/Panthera leo melanochaita ) that lived between the 15th and 20th centuries outside the current geographic distribution of lions, and 6 present-day lions from Africa and India. We found that cave and modern lions shared an ancestor ca. 500,000 y ago and that the 2 lineages likely did not hybridize following their divergence. Within modern lions, we found 2 main lineages that diverged ca. 70,000 y ago, with clear evidence of subsequent gene flow. Our data also reveal a nearly complete absence of genetic diversity within Indian lions, probably due to well-documented extremely low effective population sizes in the recent past. Our results contribute toward the understanding of the evolutionary history of lions and complement conservation efforts to protect the diversity of this vulnerable species., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

The Arctic is warming at an unprecedented rate, with unknown consequences for endemic fauna. However, Earth has experienced severe climatic oscillations in the past, and understanding how species responded to them might provide insight into their resilience to near-future climatic predictions. Little is known about the responses of Arctic marine mammals to past climatic shifts, but narwhals ( Monodon monoceros ) are considered one of the endemic Arctic species most vulnerable to environmental change. Here, we analyse 121 complete mitochondrial genomes from narwhals sampled across their range and use them in combination with species distribution models to elucidate the influence of past and ongoing climatic shifts on their population structure and demographic history. We find low levels of genetic diversity and limited geographic structuring of genetic clades. We show that narwhals experienced a long-term low effective population size, which increased after the Last Glacial Maximum, when the amount of suitable habitat expanded. Similar post-glacial habitat release has been a key driver of population size expansion of other polar marine predators. Our analyses indicate that habitat availability has been critical to the success of narwhals, raising concerns for their fate in an increasingly warming Arctic.

The protozoan Plasmodium vivax is responsible for 42% of all cases of malaria outside Africa. The parasite is currently largely restricted to tropical and subtropical latitudes in Asia, Oceania, and the Americas. Though, it was historically present in most of Europe before being finally eradicated during the second half of the 20th century. The lack of genomic information on the extinct European lineage has prevented a clear understanding of historical population structuring and past migrations of P. vivax. We used medical microscope slides prepared in 1944 from malaria-affected patients from the Ebro Delta in Spain, one of the last footholds of malaria in Europe, to generate a genome of a European P. vivax strain. Population genetics and phylogenetic analyses placed this strain basal to a cluster including samples from the Americas. This genome allowed us to calibrate a genomic mutation rate for P. vivax, and to estimate the mean age of the last common ancestor between European and American strains to the 15th century. This date points to an introduction of the parasite during the European colonization of the Americas. In addition, we found that some known variants for resistance to antimalarial drugs, including Chloroquine and Sulfadoxine, were already present in this European strain, predating their use. Our results shed light on the evolution of an important human pathogen and illustrate the value of antique medical collections as a resource for retrieving genomic information on pathogens from the past., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Domestic dogs have been central to life in the North American Arctic for millennia. The ancestors of the Inuit were the first to introduce the widespread usage of dog sledge transportation technology to the Americas, but whether the Inuit adopted local Palaeo-Inuit dogs or introduced a new dog population to the region remains unknown. To test these hypotheses, we generated mitochondrial DNA and geometric morphometric data of skull and dental elements from a total of 922 North American Arctic dogs and wolves spanning over 4500 years. Our analyses revealed that dogs from Inuit sites dating from 2000 BP possess morphological and genetic signatures that distinguish them from earlier Palaeo-Inuit dogs, and identified a novel mitochondrial clade in eastern Siberia and Alaska. The genetic legacy of these Inuit dogs survives today in modern Arctic sledge dogs despite phenotypic differences between archaeological and modern Arctic dogs. Together, our data reveal that Inuit dogs derive from a secondary pre-contact migration of dogs distinct from Palaeo-Inuit dogs, and probably aided the Inuit expansion across the North American Arctic beginning around 1000 BP.

Natural selection in domestic dogs is of great interest in evolutionary biology since dogs have migrated to every inhabited continent of the world alongside humans, and adapted to diverse environments. Here, we explored their demographic history and genetic basis of adaptation to the tropical African environment using whole genome analyses of 19 African indigenous dogs from Nigeria. Demographic analysis suggests that the ancestors of these dogs migrated into Africa from Eurasia 14,000 years ago and underwent a severe founder effect before population expansion. Admixture analysis further reveals that African dog genomes contain about 1.88-3.50% introgression from African golden wolves (Canis anthus). Population genetic analysis identifies 50 positively selected genes linked with immunity, angiogenesis, ultraviolet protection, as well as insulin secretion and sensitivity that may contribute to adaptation to tropical conditions. One of the positively selected genes, adhesion G protein-coupled receptor E1 (ADGRE1), has also been found to be association with severe malaria resistance in African human populations. Functional assessments showed that ADGRE1 provides protective host defense against Plasmodium infections. This result, together with the fact that the inflammatory response to canine babesiosis is similar to complicated falciparum malaria in humans, support the dogs as a model for the study of malaria control and treatment., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

TagSNP selection is an important step in designing case control association studies. Among selection methods that have proliferated, the ones based on pairwise LD measurement are attractive for the purpose of designing association studies. The goal is to minimize the number of markers selected for genotyping in a particular platform and therefore reduce genotyping cost while simultaneously representing information provided by all other markers. Depending on the platform, it is also important to select sets that are robust against occasional genotyping failure. An array of methods has been proposed to effectively select these tagSNPs using various criteria. In this study, we extend the algorithms used in FESTA, a computer program we previously developed for picking tagSNPs using r2 criteria. We applied FESTA to the HapMap whole chromosome data in two different populations, and we also performed a power analysis for case-control association studies using simulated data. FESTA chooses 294322 tagSNPs in the autosomes in the CEPH samples. The YORUBA samples require 61.5% more tagSNPs than the CEPH samples. The power study showed that limiting ourselves to only tagSNPs, instead of choosing all SNPs in the interval for an association study, results in a power loss of only about 5-10%.