Search

Ancient environmental DNA (aeDNA) is a crucial source of information for past environmental reconstruction. However, the computational analysis of aeDNA involves the inherited challenges of ancient DNA (aDNA) and the typical difficulties of eDNA samples, such as taxonomic identification and abundance estimation of identified taxonomic groups. Current methods for aeDNA fall into those that only perform mapping followed by taxonomic identification and those that purport to do abundance estimation. The former leaves abundance estimates to users, while methods for the latter are not designed for large metagenomic datasets and are often imprecise and challenging to use. Here, we introduce euka, a tool designed for rapid and accurate characterisation of aeDNA samples. We use a taxonomy-based pangenome graph of reference genomes for robustly assigning DNA sequences and use a maximum-likelihood framework for abundance estimation. At the present time, our database is restricted to mitochondrial genomes of tetrapods and arthropods but can be expanded in future versions. We find euka to outperform current taxonomic profiling tools and their abundance estimates. Crucially, we show that regardless of the filtering threshold set by existing methods, euka demonstrates higher accuracy. Furthermore, our approach is robust to sparse data, which is idiosyncratic of aeDNA, detecting a taxon with an average of 50 reads aligning. We also show that euka is consistent with competing tools on empirical samples. euka's features are fine-tuned to deal with the challenges of aeDNA, making it a simple-to-use, all-in-one tool. It is available on GitHub: https://github.com/grenaud/vgan. euka enables researchers to quickly assess and characterise their sample, thus allowing it to be used as a routine screening tool for aeDNA.

The application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale and heterogeneous nature of multi-modal data makes integration and inference a non-trivial task. We developed a deep-learning-based framework, multi-omics variational autoencoders (MOVE), to integrate such data and applied it to a cohort of 789 people with newly diagnosed type 2 diabetes with deep multi-omics phenotyping from the DIRECT consortium. Using in silico perturbations, we identified drug-omics associations across the multi-modal datasets for the 20 most prevalent drugs given to people with type 2 diabetes with substantially higher sensitivity than univariate statistical tests. From these, we among others, identified novel associations between metformin and the gut microbiota as well as opposite molecular responses for the two statins, simvastatin and atorvastatin. We used the associations to quantify drug-drug similarities, assess the degree of polypharmacy and conclude that drug effects are distributed across the multi-omics modalities.

Since the influenza pandemic in 2009, there has been an increased focus on swine influenza A virus (swIAV) surveillance. This paper describes the results of the surveillance of swIAV in Danish swine from 2011 to 2018. In total, 3800 submissions were received with a steady increase in swIAV positive submissions, reaching 56% in 2018. Full genome sequences were obtained from 129 swIAV positive samples. Altogether, 17 different circulating genotypes were identified including six novel reassortants harboring human seasonal IAV gene segments. The phylogenetic analysis revealed substantial genetic drift and also evidence of positive selection occurring mainly in antigenic sites of the hemagglutinin protein and confirmed the presence of a swine divergent cluster among the H lpdm09Nx (clade 1 A.3.3.2) viruses. The results provide essential data for the control of swIAV in pigs and emphasize the importance of contemporary surveillance for discovering novel swIAV strains posing a potential threat to the human population.

Beginning in late 2017, highly pathogenic avian influenza (HPAI) H5N6 viruses caused outbreaks in wild birds and poultry in several European countries. H5N6 viruses were detected in 43 wild birds found dead throughout Denmark. Most of the Danish virus-positive dead birds were found in the period from February to April 2018. However, unlike the rest of Europe, sporadic HPAI H5N6-positive dead wild birds were detected in Denmark in July, August, September, and December 2018, with the last positive bird being found in January 2019. HPAI viruses were not detected in active surveillance of apparently healthy wild birds. In this study, we use full genome sequencing and phylogenetic analysis to investigate the wild bird HPAI H5N6 viruses found in Denmark. The Danish viruses were found to be closely related to those of contemporary HPAI H5N6 viruses detected in Europe. Their sequences formed two clusters indicating that at least two or more introductions of H5N6 into Denmark occurred. Notably, all viruses detected in the latter half of 2018 and in 2019 grouped into the same cluster. The H5N6 viruses appeared to have been maintained undetected in the autumn 2018.

The degree of antigenic drift in swine influenza A viruses (swIAV) has historically been regarded as minimal compared to that of human influenza A virus strains. However, as surveillance activities on swIAV have increased, more isolates have been characterized, revealing a high level of genetic and antigenic differences even within the same swIAV lineage. The objective of this study was to investigate the level of genetic drift in one enzootically infected swine herd over one year. Nasal swabs were collected monthly from sows (n = 4) and piglets (n = 40) in the farrowing unit, and from weaners (n = 20) in the nursery. Virus from 1-4 animals were sequenced per month. Analyses of the sequences revealed that the hemagglutinin (HA) gene was the main target for genetic drift with a substitution rate of 7.6 × 10-3 substitutions/site/year and evidence of positive selection. The majority of the mutations occurred in the globular head of the HA protein and in antigenic sites. The phylogenetic tree of the HA sequences displayed a pectinate typology, where only a single lineage persists and forms the ancestor for subsequent lineages. This was most likely caused by repeated selection of a single immune-escape variant, which subsequently became the founder of the next wave of infections.

Vaccines against classical swine fever have proven very effective in protecting pigs from this deadly disease. However, little is known about how vaccination impacts the selective pressures acting on the classical swine fever virus (CSFV). Here we use high-throughput sequencing of viral genomes to investigate evolutionary changes in virus populations following the challenge of naïve and vaccinated pigs with the highly virulent CSFV strain "Koslov". The challenge inoculum contained an ensemble of closely related viral sequences, with three major haplotypes being present, termed A, B, and C. After the challenge, the viral haplotype A was preferentially located within the tonsils of naïve animals but was highly prevalent in the sera of all vaccinated animals. We find that the viral population structure in naïve pigs after infection is very similar to that in the original inoculum. In contrast, the viral population in vaccinated pigs, which only underwent transient low-level viremia, displayed several distinct changes including the emergence of 16 unique non-synonymous single nucleotide polymorphisms (SNPs) that were not detectable in the challenge inoculum. Further analysis showed a significant loss of heterogeneity and an increasing positive selection acting on the virus populations in the vaccinated pigs. We conclude that vaccination imposes a strong selective pressure on viruses that subsequently replicate within the vaccinated animal.

Influenza A virus (IAV) is a highly contagious pathogen in pigs. Swine IAV (swIAV) infection causes respiratory disease and is thereby a challenge for animal health, animal welfare and the production economy. In Europe, the most widespread strategy for controlling swIAV is implementation of sow vaccination programs, to secure delivery of protective maternally derived antibodies (MDAs) to the newborn piglets. In this study we report a unique case, where a persistently swIAV (A/sw/Denmark/P5U4/2016(H1N1)) infected herd experienced an acute outbreak with a new swIAV subtype (A/sw/Denmark/HB4280U1/2017(H1N2)) and subsequently decided to implement a mass sow vaccination program. Clinical registrations, nasal swabs and blood samples were collected from four different batches of pigs before and after vaccination. Virus isolation, sequencing of the virus strain and hemagglutinin inhibition (HI) tests were performed on samples collected before and during the outbreak and after implementation of mass sow vaccination. After implementation of the sow mass vaccination, the time of infection was delayed and the viral load significantly decreased. An increased number of pigs, however, tested positive at two consecutive sampling times indicating prolonged shedding. In addition, a significantly smaller proportion of the 10–12 weeks old pigs were seropositive by the end of the study, indicating an impaired induction of antibodies against swIAV in the presence of MDAs. Sequencing of the herd strains revealed major differences in the hemagglutinin gene of the strain isolated before- and during the acute outbreak despite that, the two strains belonged to the same HA lineage. The HI tests confirmed a limited degree of cross-reaction between the two strains. Furthermore, the sequencing results of the hemagglutinin gene obtained before and after implementation of mass sow vaccination revealed an increased substitution rate and an increase in positively selected sites in the globular hea

Influenza A virus (IAV) is a highly contagious pathogen in pigs. Swine IAV (swIAV) infection causes respiratory disease and is thereby a challenge for animal health, animal welfare and the production economy. In Europe, the most widespread strategy for controlling swIAV is implementation of sow vaccination programs, to secure delivery of protective maternally derived antibodies (MDAs) to the newborn piglets. In this study we report a unique case, where a persistently swIAV (A/sw/Denmark/P5U4/2016(H1N1)) infected herd experienced an acute outbreak with a new swIAV subtype (A/sw/Denmark/HB4280U1/2017(H1N2)) and subsequently decided to implement a mass sow vaccination program. Clinical registrations, nasal swabs and blood samples were collected from four different batches of pigs before and after vaccination. Virus isolation, sequencing of the virus strain and hemagglutinin inhibition (HI) tests were performed on samples collected before and during the outbreak and after implementation of mass sow vaccination. After implementation of the sow mass vaccination, the time of infection was delayed and the viral load significantly decreased. An increased number of pigs, however, tested positive at two consecutive sampling times indicating prolonged shedding. In addition, a significantly smaller proportion of the 10–12 weeks old pigs were seropositive by the end of the study, indicating an impaired induction of antibodies against swIAV in the presence of MDAs. Sequencing of the herd strains revealed major differences in the hemagglutinin gene of the strain isolated before- and during the acute outbreak despite that, the two strains belonged to the same HA lineage. The HI tests confirmed a limited degree of cross-reaction between the two strains. Furthermore, the sequencing results of the hemagglutinin gene obtained before and after implementation of mass sow vaccination revealed an increased substitution rate and an increase in positively selected sites in the globular head o

The description of comammox Nitrospira spp., performing complete ammonia-to-nitrate oxidation, and their co-occurrence with canonical β-proteobacterial ammonia oxidizing bacteria (β-AOB) in the environment, calls into question the metabolic potential of comammox Nitrospira and the evolutionary history of their ammonia oxidation pathway. We report four new comammox Nitrospira genomes, constituting two novel species, and the first comparative genomic analysis on comammox Nitrospira. Unlike canonical Nitrospira, comammox Nitrospira genomes lack genes for assimilatory nitrite reduction, suggesting that they have lost the potential to use external nitrite nitrogen sources. By contrast, compared to canonical Nitrospira, comammox Nitrospira harbor a higher diversity of urea transporters and copper homeostasis genes and lack cyanate hydratase genes. Additionally, the two comammox clades differ in their ammonium uptake systems. Contrary to β-AOB, comammox Nitrospira genomes have single copies of the two central ammonia oxidation pathway operons. Similar to ammonia oxidizing archaea and some oligotrophic AOB strains, they lack genes involved in nitric oxide reduction. Furthermore, comammox Nitrospira genomes encode genes that might allow efficient growth at low oxygen concentrations. Regarding the evolutionary history of comammox Nitrospira, our analyses indicate that several genes belonging to the ammonia oxidation pathway could have been laterally transferred from β-AOB to comammox Nitrospira. We postulate that the absence of comammox genes in other sublineage II Nitrospira genomes is the result of subsequent loss.

Inflammatory bowel disease (IBD) is a chronic intestinal disorder, with two main types: Crohn's disease (CD) and ulcerative colitis (UC), whose molecular pathology is not well understood. The majority of IBD-associated SNPs are located in non-coding regions and are hard to characterize since regulatory regions in IBD are not known. Here we profile transcription start sites (TSSs) and enhancers in the descending colon of 94 IBD patients and controls. IBD-upregulated promoters and enhancers are highly enriched for IBD-associated SNPs and are bound by the same transcription factors. IBD-specific TSSs are associated to genes with roles in both inflammatory cascades and gut epithelia while TSSs distinguishing UC and CD are associated to gut epithelia functions. We find that as few as 35 TSSs can distinguish active CD, UC, and controls with 85% accuracy in an independent cohort. Our data constitute a foundation for understanding the molecular pathology, gene regulation, and genetics of IBD.

Aleutian mink disease virus (AMDV) is a frequently encountered pathogen associated with mink farming. Previous phylogenetic analyses of AMDV have been based on shorter and more conserved parts of the genome, e.g. the partial NS1 gene. Such fragments are suitable for detection but are less useful for elucidating transmission pathways while sequencing entire viral genomes provides additional informative sites and often results in better-resolved phylogenies. We explore how whole-genome sequencing can benefit investigations of AMDV transmission by reconstructing the relationships between AMDV field samples from a Danish outbreak. We show that whole-genome phylogenies are much better resolved than those based on the partial NS1 gene sequences extracted from the same alignment. Well-resolved phylogenies contain more information about the underlying transmission trees and are useful for understanding the spread of a pathogen. In the main case investigated here, the transmission path suggested by the tree structure was supported by epidemiological data. The use of molecular clock models further improved tree resolution and provided time estimates for the viral ancestors consistent with the proposed direction of spread. It was however impossible to infer transmission pathways from the partial NS1 gene tree, since all samples from the case farms branched out from a single internal node. A sliding window analysis showed that there were no shorter genomic regions providing the same phylogenetic resolution as the entire genome. Altogether, these results suggest that phylogenetic analyses based on whole-genome sequencing taking into account sampling dates and epidemiological data is a promising set of tools for clarifying AMDV transmission.

Oncogenic human papillomaviruses (HPVs) are in most cases eliminated by intervention of T cells. As many other pathogens, these oncogenic HPVs belong to an ancient and diverse virus family. Therefore, we found it relevant to investigate the potential and limitations of inducing a broad response - either by inducing cross-reactive T cells or by administering a polyvalent vaccine. To test these strategies, we designed 3 ancestral and 2 circulating sequences based on the two domains of the E1 and E2 proteins of papillomaviruses (PVs) that exhibit the highest degree of conservation in comparison to the other PV proteins. The PV sequences were fused to a T cell adjuvant, the murine invariant chain and encoded in a recombinant adenoviral vector which was administered to naïve outbred mice. By measuring T cell responses induced by these different vaccines and towards peptide pools representing 3 circulating strains and a putative ancestor of oncogenic HPVs, we showed that the ancestral vaccine antigen has to be approximately 90% identical to the circulating PVs before a marked drop of ~90% mean CD8+ T cell responses ensues. Interestingly, the combination of two or three type-specific PV vaccines did not induce a significant decrease of the CD8+ T cell response to the individual targeted PV types. Polyvalent HPV vaccine based on the E1 and E2 proteins seem to be capable of triggering responses towards more than one type of PV while the cross-reactivity of ancestral vaccine seems insufficient in consideration of the sequence diversity between HPV types.

Nitrification, the biological oxidation of ammonium to nitrate, is a fundamental process in the nitrogen cycle and plays an important role in natural and engineered systems. Throughout the last century, nitrification was assumed to be a two-step process executed by two different functional groups, ammonia oxidizing prokaryotes (AOP) and nitrite oxidizing bacteria (NOB). Recently, several articles have shown the capability of a single microorganism, belonging to the genus Nitrospira, to carry out the complete oxidation of ammonia to nitrate (comammox). Nitrospira spp. are widespread in both natural and engineered ecosystems associated with nitrogen cycling and different species are frequently observed to coexist in the same environment. Besides recent discoveries pointing towards versatile metabolism in some Nitrospira species, little is known about the functional potential of the two comammox Nitrospira clades, and the factors involved in niche-partitioning between comammox and canonical Nitrospira. A comparative genomics analysis was conducted with five genomes recovered from a groundwater-fed rapid sand filter (including both comammox clades and a nitrite-oxidizing Nitrospira population genome) and high quality published Nitrospira genomes, to reveal distinct genomic features within Nitrospira. In addition, we investigated the evolution of the ammonia oxidation pathway in comammox Nitrospira. This analysis revealed distinct genetic capabilities of the different comammox clades and canonical Nitrospira which can help to explain the coexistence and niche partitioning of Nitrospira spp. These divergences range from the nitrogen source utilization capacity to the ability for electron donor versatility, and other characteristics such as stress response. With respect to the evolutionary history of comammox Nitrospira, our analysis indicates transfer events with betaproteobacterial ammonia oxidizers. In addition, transfer events between comammox clade A and clade B w