Genome Reconstruction Research Articles

Targeted enrichment of whole-genome SNPs from highly burned skeletal remains.

Genetic assessment of highly incinerated and/or degraded human skeletal material is a persistent challenge in forensic DNA analysis, including identifying victims of mass disasters. Few studies have investigated the impact of thermal degradation on whole-genome single-nucleotide polymorphism (SNP) quality and quantity using next-generation sequencing (NGS). We present whole-genome SNP data obtained from the bones and teeth of 27 fire victims using two DNA extraction techniques. Extracts were converted to double-stranded DNA libraries then enriched for whole-genome SNPs using unpublished biotinylated RNA baits and sequenced on an Illumina NextSeq 550 platform. Raw reads were processed using the EAGER (Efficient Ancient Genome Reconstruction) pipeline, and the SNPs filtered and called using FreeBayes and GATK (v. 3.8). Mixed-effects modeling of the data suggest that SNP variability and preservation is predominantly determined by skeletal element and burn category, and not by extraction type. Whole-genome SNP data suggest that selecting long bones, hand and foot bones, and teeth subjected to temperatures <350°C are the most likely sources for higher genomic DNA yields. Furthermore, we observed an inverse correlation between the number of captured SNPs and the extent to which samples were burned, as well as a significant decrease in the total number of SNPs measured for samples subjected to temperatures >350°C. Our data complement previous analyses of burned human remains that compare extraction methods for downstream forensic applications and support the idea of adopting a modified Dabney extraction technique when traditional forensic methods fail to produce DNA yields sufficient for genetic identification.

3D Genome Reconstruction from Partially Phased Hi-C Data

The 3-dimensional (3D) structure of the genome is of significant importance for many cellular processes. In this paper, we study the problem of reconstructing the 3D structure of chromosomes from Hi-C data of diploid organisms, which poses additional challenges compared to the better-studied haploid setting. With the help of techniques from algebraic geometry, we prove that a small amount of phased data is sufficient to ensure finite identifiability, both for noiseless and noisy data. In the light of these results, we propose a new 3D reconstruction method based on semidefinite programming, paired with numerical algebraic geometry and local optimization. The performance of this method is tested on several simulated datasets under different noise levels and with different amounts of phased data. We also apply it to a real dataset from mouse X chromosomes, and we are then able to recover previously known structural features.

ContScout: sensitive detection and removal of contamination from annotated genomes

Contamination of genomes is an increasingly recognized problem affecting several downstream applications, from comparative evolutionary genomics to metagenomics. Here we introduce ContScout, a precise tool for eliminating foreign sequences from annotated genomes. It achieves high specificity and sensitivity on synthetic benchmark data even when the contaminant is a closely related species, outperforms competing tools, and can distinguish horizontal gene transfer from contamination. A screen of 844 eukaryotic genomes for contamination identified bacteria as the most common source, followed by fungi and plants. Furthermore, we show that contaminants in ancestral genome reconstructions lead to erroneous early origins of genes and inflate gene loss rates, leading to a false notion of complex ancestral genomes. Taken together, we offer here a tool for sensitive removal of foreign proteins, identify and remove contaminants from diverse eukaryotic genomes and evaluate their impact on phylogenomic analyses.

Chronic shedding of a SARS-CoV-2 Alpha variant in wastewater

BackgroundCentral Michigan University (CMU) participated in a state-wide SARS-CoV-2 wastewater monitoring program since 2021. Wastewater samples were collected from on-campus sites and nine off-campus wastewater treatment plants servicing small metropolitan and rural communities. SARS-CoV-2 genome copies were quantified using droplet digital PCR and results were reported to the health department.ResultsOne rural, off-campus site consistently produced higher concentrations of SARS-CoV-2 genome copies. Samples from this site were sequenced and contained predominately a derivative of Alpha variant lineage B.1.1.7, detected from fall 2021 through summer 2023. Mutational analysis of reconstructed genes revealed divergence from the Alpha variant lineage sequence over time, including numerous mutations in the Spike RBD and NTD.ConclusionsWe discuss the possibility that a chronic SARS-CoV-2 infection accumulated adaptive mutations that promoted long-term infection. This study reveals that small wastewater treatment plants can enhance resolution of rare events and facilitate reconstruction of viral genomes due to the relative lack of contaminating sequences.

The Natural and Clinical History of Plague: From the Ancient Pandemics to Modern Insights.

The human pathogen Yersinia pestis is responsible for bubonic, septicemic, and pneumonic plague. A deeply comprehensive overview of its historical context, bacteriological characteristics, genomic analysis based on ancient DNA (aDNA) and modern strains, and its impact on historical and actual human populations, is explored. The results from multiple studies have been synthesized to investigate the origins of plague, its transmission, and effects on different populations. Additionally, molecular interactions of Y. pestis, from its evolutionary origins to its adaptation to flea-born transmission, and its impact on human and wild populations are considered. The characteristic combinations of aDNA patterns, which plays a decisive role in the reconstruction and analysis of ancient genomes, are reviewed. Bioinformatics is fundamental in identifying specific Y. pestis lineages, and automated pipelines are among the valuable tools in implementing such studies. Plague, which remains among human history's most lethal infectious diseases, but also other zoonotic diseases, requires the continuous investigation of plague topics. This can be achieved by improving molecular and genetic screening of animal populations, identifying ecological and social determinants of outbreaks, increasing interdisciplinary collaborations among scientists and public healthcare providers, and continued research into the characterization, diagnosis, and treatment of these diseases.

Spatiotemporal dynamics of giant viruses within a deep freshwater lake reveal a distinct dark-water community.

Giant viruses (GVs) significantly regulate the ecological dynamics of diverse ecosystems. Although metagenomics has expanded our understanding of their diversity and ecological roles played in marine environments, little is known about GVs of freshwater ecosystems. Most previous studies have employed short-read sequencing and therefore resulted in fragmented genomes, hampering accurate assessment of genetic diversity. We sought to bridge this knowledge gap and overcome previous technical limitations. We subjected spatiotemporal (2 depths × 12months) samples from Lake Biwa to metagenome-assembled genome reconstruction enhanced by long-read metagenomics. This yielded 293 GV metagenome-assembled genomes. Of these, 285 included previously unknown species in five orders of nucleocytoviruses and the first representatives of freshwater mirusviruses, which exhibited marked divergence from marine-derived lineages. The good performance of our long-read metagenomic assembly was demonstrated by the detection of 42 (14.3%) genomes composed of single contigs with completeness values >90%. GVs were partitioned across water depths, with most species specific to either the sunlit epilimnion or the dark hypolimnion. Epilimnion-specific members tended to be transient and exhibit short and intense abundance peaks, in line with the fact that they regulate the surface algal blooms. During the spring bloom, mirusviruses and members of three nucleocytovirus families were among the most abundant viruses. In contrast, hypolimnion-specific ones, including a mirusvirus genome, were typically more persistent in the hypolimnion throughout the water-stratified period, suggesting that they infect hosts specific to the hypolimnion and play previously unexplored ecological roles in dark water microbial ecosystems.

Temporal colonization and metabolic regulation of the gut microbiome in neonatal oxen at single nucleotide resolution.

The colonization of microbes in the gut is key to establishing a healthy host-microbiome symbiosis for newborns. We longitudinally profiled the gut microbiome in a model consisting of 36 neonatal oxen from birth up to 2months postpartum and carried out microbial transplantation to reshape their gut microbiome. Genomic reconstruction of deeply sequenced fecal samples resulted in a total of 3931 metagenomic-assembled genomes from 472 representative species, of which 184 were identified as new species when compared with existing databases of oxen. Single nucleotide level metagenomic profiling shows a rapid influx of microbes after birth, followed by dynamic shifts during the first few weeks of life. Microbial transplantation was found to reshape the genetic makeup of 33 metagenomic-assembled genomes (FDR < 0.05), mainly from Prevotella and Bacteroides species. We further linked over 20 million microbial single nucleotide variations to 736 plasma metabolites, which enabled us to characterize 24 study-wide significant associations (P < 4.4 × 10-9) that identify the potential microbial genetic regulation of host immune and neuro-related metabolites, including glutathione and L-dopa. Our integration analyses further revealed that microbial genetic variations may influence the health status and growth performance by modulating metabolites via structural regulation of their encoded proteins. For instance, we found that the albumin levels and total antioxidant capacity were correlated with L-dopa, which was determined by single nucleotide variations via structural regulations of metabolic enzymes. The current results indicate that temporal colonization and transplantation-driven strain replacement are crucial for newborn gut development, offering insights for enhancing newborn health and growth.

CoCoPyE: feature engineering for learning and prediction of genome quality indices.

The exploration of the microbial world has been greatly advanced by the reconstruction of genomes from metagenomic sequence data. However, the rapidly increasing number of metagenome-assembled genomes has also resulted in a wide variation in data quality. It is therefore essential to quantify the achieved completeness and possible contamination of a reconstructed genome before it is used in subsequent analyses. The classical approach for the estimation of quality indices solely relies on a relatively small number of universal single-copy genes. Recent tools try to extend the genomic coverage of estimates for an increased accuracy. We developed CoCoPyE, a fast tool based on a novel 2-stage feature extraction and transformation scheme. First, it identifies genomic markers and then refines the marker-based estimates with a machine learning approach. In our simulation studies, CoCoPyE showed a more accurate prediction of quality indices than the existing tools. While the CoCoPyE web server offers an easy way to try out the tool, the freely available Python implementation enables integration into existing genome reconstruction pipelines. CoCoPyE provides a new approach to assess the quality of genome data. It complements and improves existing tools and may help researchers to better distinguish between low-quality draft and high-quality genome assemblies in metagenome sequencing projects.

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing.

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

Reconstruction of private genomes through reference-based genotype imputation

BackgroundGenotype imputation is an essential step in genetic studies to improve data quality and statistical power. Public imputation servers are widely used by researchers to impute their data using otherwise access-controlled reference panels of high-fidelity genomes held by these servers.ResultsWe report evidence against the prevailing assumption that providing access to panels only indirectly via imputation servers poses a negligible privacy risk to individuals in the panels. To this end, we present algorithmic strategies for adaptively constructing artificial input samples and interpreting their imputation results that lead to the accurate reconstruction of reference panel haplotypes. We illustrate this possibility on three reference panels of real genomes for a range of imputation tools and output settings. Moreover, we demonstrate that reconstructed haplotypes from the same individual could be linked via their genetic relatives using our Bayesian linking algorithm, which allows a substantial portion of the individual’s diploid genome to be reassembled. We also provide population genetic estimates of the proportion of a panel that could be linked when an adversary holds a varying number of genomes from the same population.ConclusionsOur results show that genomes in imputation server reference panels can be vulnerable to reconstruction, implying that additional safeguards may need to be considered. We suggest possible mitigation measures based on our findings. Our work illustrates the value of adversarial algorithms in uncovering new privacy risks to help inform the genomics community towards secure data sharing practices.

Reconstruction of Ancestral Genomes as a Key to Understanding the Early Evolution of Vertebrate Genotype

Reconstruction of Ancestral Genomes as a Key to Understanding the Early Evolution of Vertebrate Genotype

Dancing the Nanopore limbo – Nanopore metagenomics from small DNA quantities for bacterial genome reconstruction

BackgroundWhile genome-resolved metagenomics has revolutionized our understanding of microbial and genetic diversity in environmental samples, assemblies of short-reads often result in incomplete and/or highly fragmented metagenome-assembled genomes (MAGs), hampering in-depth genomics. Although Nanopore sequencing has increasingly been used in microbial metagenomics as long reads greatly improve the assembly quality of MAGs, the recommended DNA quantity usually exceeds the recoverable amount of DNA of environmental samples. Here, we evaluated lower-than-recommended DNA quantities for Nanopore library preparation by determining sequencing quality, community composition, assembly quality and recovery of MAGs.ResultsWe generated 27 Nanopore metagenomes using the commercially available ZYMO mock community and varied the amount of input DNA from 1000 ng (the recommended minimum) down to 1 ng in eight steps. The quality of the generated reads remained stable across all input levels. The read mapping accuracy, which reflects how well the reads match a known reference genome, was consistently high across all libraries. The relative abundance of the species in the metagenomes was stable down to input levels of 50 ng. High-quality MAGs (> 95% completeness, ≤ 5% contamination) could be recovered from metagenomes down to 35 ng of input material. When combined with publicly available Illumina reads for the mock community, Nanopore reads from input quantities as low as 1 ng improved the quality of hybrid assemblies.ConclusionOur results show that the recommended DNA amount for Nanopore library preparation can be substantially reduced without any adverse effects to genome recovery and still bolster hybrid assemblies when combined with short-read data. We posit that the results presented herein will enable studies to improve genome recovery from low-biomass environments, enhancing microbiome understanding.

Genomics of Yoonia sp. Isolates (Family Roseobacteraceae) from Lake Zhangnai on the Tibetan Plateau.

Understanding the genomic differentiation between marine and non-marine aquatic microbes remains a compelling question in ecology. While previous research has identified several lacustrine lineages within the predominantly marine Roseobacteraceae family, limited genomic data have constrained our understanding of their ecological adaptation mechanisms. In this study, we isolated four novel Yoonia strains from a brackish lake on the Tibetan Plateau. These strains have diverged from their marine counterparts within the same genus, indicating a recent habitat transition event from marine to non-marine environments. Metabolic comparisons and ancestral genomic reconstructions in a phylogenetic framework reveal metabolic shifts in salinity adaptation, compound transport, aromatics degradation, DNA repair, and restriction systems. These findings not only corroborate the metabolic changes commonly observed in other non-marine Roseobacters but also unveil unique adaptations, likely reflecting the localized metabolic changes in responses to Tibetan Plateau environments. Collectively, our study expands the known genomic diversity of non-marine Roseobacteraceae lineages and enhances our understanding of microbial adaptations to lacustrine ecosystems.

Complete genome reconstruction of the global and European regional dispersal history of the lumpy skin disease virus.

Lumpy skin disease virus (LSDV) has a complex epidemiology involving multiple strains, recombination, and vaccination. Its DNA genome provides limited genetic variation to trace outbreaks in space and time. Sequencing of LSDV whole genomes has also been patchy at global and regional scales. Here, we provide the first fine-grained whole genome sequence sampling of a constrained LSDV outbreak (southeastern Europe, 2015-2017), which we analyze along with global publicly available genomes. We formally evaluate the past occurrence of recombination events as well as the temporal signal that is required for calibrating molecular clock models and subsequently conduct a time-calibrated spatially explicit phylogeographic reconstruction. Our study further illustrates the importance of accounting for recombination events before reconstructing global and regional dynamics of DNA viruses. More LSDV whole genomes from endemic areas are needed to obtain a comprehensive understanding of global LSDV dispersal dynamics.

Emergence of Elizabethkingia miricola in the Brazilian Amazon region

Infectious diseases are not only a threat to humans and animals but also to an entire ecosystem if the three axis of One Health are considered. Elizabethkingia miricola is a water-borne and opportunistic pathogen frequently evolving to high morbidity and mortality outcomes. More recently, reports on E. miricola infections causing death in different anuran species are increasing. Here, we reported the emergence, in a clinical setting in Brazil (Amazon region), of an E. miricola (EM15) with an almost untreatable antibiotic resistance phenotype that was closely related to the lineage involved with anuran infections worldwide. A genomic reconstruction revealed that EM15 belonged to a cluster represented by genomes from human and animal (anurans) sources recovered from China, Europe, and the Americas from 2002 to 2021, and from the Space Station Mir condensation water. Their resistome and virulome were also in agreement with their genetic relationship since they harboured the same set of resistance and putative virulence genes. Therefore, the emergence of EM15 in a clinical setting in the Amazon region, a hotspot area of anuran diversity and home to Amerindian groups that live in close relationship with the wildlife, may represent a warning to the public and environmental health.

Extension of the Segatella copri complex to 13 species with distinct large extrachromosomal elements and associations with host conditions

The Segatella copri (formerly Prevotella copri) complex (ScC) comprises taxa that are key members of the human gut microbiome. It was previously described to contain four distinct phylogenetic clades. Combining targeted isolation with large-scale metagenomic analysis, we defined 13 distinct Segatella copri-related species, expanding the ScC complex beyond four clades. Complete genome reconstruction of thirteen strains from seven species unveiled the presence of genetically diverse large circular extrachromosomal elements. These elements are consistently present in most ScC species, contributing to intra- and inter-species diversities. The nine species-level clades present in humans display striking differences in prevalence and intra-species genetic makeup across human populations. Based on a meta-analysis, we found reproducible associations between members of ScC and the male sex and positive correlations with lower visceral fat and favorable markers of cardiometabolic health. Our work uncovers genomic diversity within ScC, facilitating a better characterization of the human microbiome.

Comparative genomics of Mycoplasma feriruminatoris, a fast-growing pathogen of wild Caprinae.

Mycoplasma feriruminatoris is a fast-growing Mycoplasma species isolated from wild Caprinae and first described in 2013. M. feriruminatoris isolates have been associated with arthritis, kerato conjunctivitis, pneumonia and septicemia, but were also recovered from apparently healthy animals. To better understand what defines this species, we performed a genomic survey on 14 strains collected from free-ranging or zoo-housed animals between 1987 and 2017, mostly in Europe. The average chromosome size of the M. feriruminatoris strains was 1,040±0,024 kbp, with 24 % G+C and 852±31 CDS. The core genome and pan-genome of the M. feriruminatoris species contained 628 and 1312 protein families, respectively. The M. feriruminatoris strains displayed a relatively closed pan-genome, with many features and putative virulence factors shared with species from the M. mycoides cluster, including the MIB-MIP Ig cleavage system, a repertoire of DUF285 surface proteins and a complete biosynthetic pathway for galactan. M. feriruminatoris genomes were found to be mostly syntenic, although repertoires of mobile genetic elements, including Mycoplasma Integrative and Conjugative Elements, insertion sequences, and a single plasmid varied. Phylogenetic- and gene content analyses confirmed that M. feriruminatoris was closer to the M. mycoides cluster than to the ruminant species M. yeatsii and M. putrefaciens. Ancestral genome reconstruction showed that the emergence of the M. feriruminatoris species was associated with the gain of 17 gene families, some of which encode defence enzymes and surface proteins, and the loss of 25 others, some of which are involved in sugar transport and metabolism. This comparative study suggests that the M. mycoides cluster could be extended to include M. feriruminatoris. We also find evidence that the specific organization and structure of the DnaA boxes around the oriC of M. feriruminatoris may contribute to drive the remarkable fast growth of this minimal bacterium.

Microbial degradation and assimilation of veratric acid in oxic and anoxic groundwaters.

Microbial communities are key players in groundwater ecosystems. In this dark environment, heterotrophic microbes rely on biomass produced by the activity of lithoautotrophs or on the degradation of organic matter seeping from the surface. Most studies on bacterial diversity in groundwater habitats are based on 16S gene sequencing and full genome reconstructions showing potential metabolic pathways used in these habitats. However, molecular-based studies do not allow for the assessment of population dynamics over time or the assimilation of specific compounds and their biochemical transformation by microbial communities. Therefore, in this study, we combined DNA-, phospholipid fatty acid-, and metabolomic-stable isotope probing to target and identify heterotrophic bacteria in the groundwater setting of the Hainich Critical Zone Exploratory (CZE), focusing on 2 aquifers with different physico-chemical conditions (oxic and anoxic). We incubated groundwater from 4 different wells using either 13C-labeled veratric acid (a lignin-derived compound) (single labeling) or a combination of 13CO2 and D-labeled veratric acid (dual labeling). Our results show that heterotrophic activities dominate all groundwater sites. We identified bacteria with the potential to break down veratric acid (Sphingobium or Microbacterium). We observed differences in heterotrophic activities between the oxic and anoxic aquifers, indicating local adaptations of bacterial populations. The dual labeling experiments suggested that the serine pathway is an important carbon assimilation pathway and that organic matter was an important source of hydrogen in the newly produced lipids. These experiments also yielded different labeled taxa compared to the single labeling experiments, showing that there exists a complex interaction network in the groundwater habitats.

A practical assembly guideline for genomes with various levels of heterozygosity.

Although current long-read sequencing technologies have a long-read length that facilitates assembly for genome reconstruction, they have high sequence errors. While various assemblers with different perspectives have been developed, no systematic evaluation of assemblers with long reads for diploid genomes with varying heterozygosity has been performed. Here, we evaluated a series of processes, including the estimation of genome characteristics such as genome size and heterozygosity, de novo assembly, polishing, and removal of allelic contigs, using six genomes with various heterozygosity levels. We evaluated five long-read-only assemblers (Canu, Flye, miniasm, NextDenovo and Redbean) and five hybrid assemblers that combine short and long reads (HASLR, MaSuRCA, Platanus-allee, SPAdes and WENGAN) and proposed a concrete guideline for the construction of haplotype representation according to the degree of heterozygosity, followed by polishing and purging haplotigs, using stable and high-performance assemblers: Redbean, Flye and MaSuRCA.

More than a decade of genetic research on the Denisovans.

Denisovans, a group of now extinct humans who lived in Eastern Eurasia in the Middle and Late Pleistocene, were first identified from DNA sequences just over a decade ago. Only ten fragmentary remains from two sites have been attributed to Denisovans based entirely on molecular information. Nevertheless, there has been great interest in using genetic data to understand Denisovans and their place in human history. From the reconstruction of a single high-quality genome, it has been possible to infer their population history, including events of admixture with other human groups. Additionally, the identification of Denisovan DNA in the genomes of present-day individuals has provided insights into the timing and routes of dispersal of ancient modern humans into Asia and Oceania, as well as the contributions of archaic DNA to the physiology of present-day people. In this Review, we synthesize more than a decade of research on Denisovans, reconcile controversies and summarize insights into their population history and phenotype. We also highlight how our growing knowledge about Denisovans has provided insights into our own evolutionary history.