Worldwide human mitochondrial haplogroup distribution from urban sewage

Time-series sewage metagenomics distinguishes seasonal, human-derived and environmental microbial communities potentially allowing source-attributed surveillance

Sewage metagenomics has risen to prominence in urban population surveillance of pathogens and antimicrobial resistance (AMR). Unknown species with similarity to known genomes cause database bias in reference-based metagenomics. To improve surveillance, we seek to recover sewage genomes and develop a quantification and correlation workflow for these genomes and AMR over time. We use longitudinal sewage sampling in seven treatment plants from five major European cities to explore the utility of catch-all sequencing of these population-level samples. Using metagenomic assembly methods, we recover 2332 metagenome-assembled genomes (MAGs) from prokaryotic species, 1334 of which were previously undescribed. These genomes account for ~69% of sequenced DNA and provide insight into sewage microbial dynamics. Rotterdam (Netherlands) and Copenhagen (Denmark) show strong seasonal microbial community shifts, while Bologna, Rome, (Italy) and Budapest (Hungary) have occasional blooms of Pseudomonas-dominated communities, accounting for up to ~95% of sample DNA. Seasonal shifts and blooms present challenges for effective sewage surveillance. We find that bacteria of known shared origin, like human gut microbiota, form communities, suggesting the potential for source-attributing novel species and their ARGs through network community analysis. This could significantly improve AMR tracking in urban environments.

mtDNA "nomenclutter" and its consequences on the interpretation of genetic data

Population-based studies of human mitochondrial genetic diversity often require the classification of mitochondrial DNA (mtDNA) haplotypes into more than 5400 described haplogroups, and further grouping those into hierarchically higher haplogroups. Such secondary haplogroup groupings (e.g., "macro-haplogroups") vary across studies, as they depend on the sample quality, technical factors of haplogroup calling, the aims of the study, and the researchers' understanding of the mtDNA haplogroup nomenclature. Retention of historical nomenclature coupled with a growing number of newly described mtDNA lineages results in increasingly complex and inconsistent nomenclature that does not reflect phylogeny well. This "clutter" leaves room for grouping errors and inconsistencies across scientific publications, especially when the haplogroup names are used as a proxy for secondary groupings, and represents a source for scientific misinterpretation. Here we explore the effects of phylogenetically insensitive secondary mtDNA haplogroup groupings, and the lack of standardized secondary haplogroup groupings on downstream analyses and interpretation of genetic data. We demonstrate that frequency-based analyses produce inconsistent results when different secondary mtDNA groupings are applied, and thus allow for vastly different interpretations of the same genetic data. The lack of guidelines and recommendations on how to choose appropriate secondary haplogroup groupings presents an issue for the interpretation of results, as well as their comparison and reproducibility across studies. To reduce biases originating from arbitrarily defined secondary nomenclature-based groupings, we suggest that future updates of mtDNA phylogenies aimed for the use in mtDNA haplogroup nomenclature should also provide well-defined and standardized sets of phylogenetically meaningful algorithm-based secondary haplogroup groupings such as "macro-haplogroups", "meso-haplogroups", and "micro-haplogroups". Ideally, each of the secondary haplogroup grouping levels should be informative about different human population history events. Those phylogenetically informative levels of haplogroup groupings can be easily defined using TreeCluster, and then implemented into haplogroup callers such as HaploGrep3. This would foster reproducibility across studies, provide a grouping standard for population-based studies, and reduce errors associated with haplogroup nomenclatures in future studies.

Balancing public health and group privacy: Ethics, rights, and obligations for wastewater surveillance systems

As the threat of COVID-19 recedes, wastewater surveillance - unlike other pandemic-era public health surveillance methods - seems here to stay. Concerns have been raised, however, about the potential risks that wastewater surveillance might pose towards group privacy. Existing scholarship has focused upon using ethics- or human rights-based frameworks as a means of balancing the public health objectives of wastewater surveillance and the potential risks it might pose to group privacy. However, such frameworks greatly lack enforceability. In order to further the strong foundation laid by such frameworks - while addressing their lack of enforceability - this paper proposes the idea of the 'obligation' as an alternative way to regulate wastewater surveillance systems. The legal codification of said obligations provides a method of ensuring that wastewater surveillance systems can be deployed effectively and equitably. Our paper proposes that legal obligations for wastewater surveillance can be created and enforced through transparent and purposeful legislation (which would include limits on power and grant institutions substantial oversight) as well as paying heed to non-legislative legal means of enforcement, such as through courts or contracts. Introducing legal obligations for wastewater surveillance could therefore be highly useful to researchers, policymakers, corporate technologists, and government agencies working in this field.

Comparison of the optimal and suboptimal quantity of mitotype libraries using next-generation sequencing

Optimizing analysis parameters and sample input is crucial in forensic genetics methods to generate reliable results, and even more so when working with muti-copy mitochondrial DNA (mtDNA) and low-quality samples. This study compared mitotypes based on next-generation sequencing (NGS) results derived from the same samples at two different sequencing library concentrations-30 pM and 0.3 pM. Thirty femur samples from the Second World War were used as a model for poorly preserved DNA. Quantitative PCR (qPCR) method targeting 113 bp long fragment was employed to assess the quantity of mitogenomes. HID Ion Chef™ Instrument with Precision ID mtDNA Control Region Panel was used for library preparation and templating. Sequencing was performed with Ion GeneStudio™ S5 System. Reference haplotypes were determined from sequencing samples at 30 pM library input. Haplotypes were compared between optimal (30 pM) and suboptimal (0.3 pM) library inputs. Often the difference in haplotypes was length heteroplasmy, which in line with other studies shows that this type of variant is not reliable for interpretation in forensics. Excluding length variants at positions 573, 309, and 16,193, 56.7% of the samples matched, and in two samples, no sequence was obtained at suboptimal library input. The rest of the samples differed between optimal and suboptimal library input. To conclude, genotyping and analyzing low-quantity libraries derived from low-quality aged skeletonized human remains therefore must be done with caution in forensic genetics casework.

Measure quantity of mitochondrial DNA in aged bones or calculate it from nuclear DNA quantitative PCR results?

Mitochondrial DNA (mtDNA) is of great value in forensics to procure information about a person when a next of kin, personal belongings, or other sources of nuclear DNA (nDNA) are unavailable, or nDNA is lacking in quality and quantity. The quality and reliability of the results depend greatly on ensuring optimal conditions for the given method, for instance, the optimal input of the copy number (CN) in next-generation sequencing (NGS) methods. The unavailability of commercial quantitative PCR (qPCR) methods to determine mtDNA CN creates the necessity to rely on recommendations to infer mtDNA CN from nDNA yield. Because nDNA yield varies between individuals, tissues, parts of the same tissue, and because mtDNA CN varies between tissues, such assumptions must be examined for a specific context, rather than be generalized. This study compares mtDNA CN calculated from nDNA yield and qPCR measured mtDNA CN. Seventy-five femurs from the Second World War victims were used as samples; they were cut below the greater trochanter, surface contaminants were removed by mechanical and chemical cleaning, samples were fully demineralized, and DNA was isolated. PowerQuant® Kit (Promega) was used to analyze DNA yield. An in-house method was used to determine mtDNA CN. Comparison of mtDNA CN from nDNA derived calculations and measured mtDNA CN highlighted vast differences. The results emphasize the need to perform qPCR to assess mtDNA CN before NGS analyses of aged bones' mitogenomes rather than estimating mtDNA CN from nDNA yield to ensure the quality and reliability of the results of NGS analysis.

Inadvertent human genomic bycatch and intentional capture raise beneficial applications and ethical concerns with environmental DNA

The field of environmental DNA (eDNA) is advancing rapidly, yet human eDNA applications remain underutilized and underconsidered. Broader adoption of eDNA analysis will produce many well-recognized benefits for pathogen surveillance, biodiversity monitoring, endangered and invasive species detection, and population genetics. Here we show that deep-sequencing-based eDNA approaches capture genomic information from humans (Homo sapiens) just as readily as that from the intended target species. We term this phenomenon human genetic bycatch (HGB). Additionally, high-quality human eDNA could be intentionally recovered from environmental substrates (water, sand and air), holding promise for beneficial medical, forensic and environmental applications. However, this also raises ethical dilemmas, from consent, privacy and surveillance to data ownership, requiring further consideration and potentially novel regulation. We present evidence that human eDNA is readily detectable from 'wildlife' environmental samples as human genetic bycatch, demonstrate that identifiable human DNA can be intentionally recovered from human-focused environmental sampling and discuss the translational and ethical implications of such findings.

Global Distribution and Diversity of Prevalent Sewage Water Plasmidomes

Sewage water from around the world contains an abundance of short plasmids, several of which harbor antimicrobial resistance genes (ARGs). The global dynamics of plasmid-derived antimicrobial resistance and functions are only starting to be unveiled. Here, we utilized a previously created data set of 159,332 assumed small plasmids from 24 different global sewage samples. The detailed phylogeny, as well as the interplay between their protein domains, ARGs, and predicted bacterial host genera, were investigated to understand sewage plasmidome dynamics globally. A total of 58,429 circular elements carried genes encoding plasmid-related features, and MASH distance analyses showed a high degree of diversity. A single (yet diverse) cluster of 520 predicted Acinetobacter plasmids was predominant among the European sewage water. Our results suggested a prevalence of plasmid-backbone gene combinations over others. This could be related to selected bacterial genera that act as bacterial hosts. These combinations also mirrored the geographical locations of the sewage samples. Our functional domain network analysis identified three groups of plasmids. However, these backbone domains were not exclusive to any given group, and Acinetobacter was the dominant host genus among the theta-replicating plasmids, which contained a reservoir of the macrolide resistance gene pair msr(E) and mph(E). Macrolide resistance genes were the most common in the sewage plasmidomes and were found in the largest number of unique plasmids. While msr(E) and mph(E) were limited to Acinetobacter, erm(B) was disseminated among a range of Firmicutes plasmids, including Staphylococcus and Streptococcus, highlighting a potential reservoir of antibiotic resistance for these pathogens from around the globe. IMPORTANCE Antimicrobial resistance is a global threat to human health, as it inhibits our ability to treat infectious diseases. This study utilizes sewage water plasmidomes to identify plasmid-derived features and highlights antimicrobial resistance genes, particularly macrolide resistance genes, as abundant in sewage water plasmidomes in Firmicutes and Acinetobacter hosts. The emergence of macrolide resistance in these bacteria suggests that macrolide selective pressure exists in sewage water and that the resident bacteria can readily acquire macrolide resistance via small plasmids.

Extreme intraspecific divergence in mitochondrial haplotypes makes the threespine stickleback fish an emerging evolutionary mutant model for mito-nuclear interactions

Mitochondrial DNA is primarily maternally inherited in most animals and evolves about 10 times faster than biparentally inherited nuclear DNA. Mitochondrial dysfunction (mt-dys) arises when interactions between the co-evolving mitochondrial and nuclear genomes are perturbed in essential processes like oxidative phosphorylation (OXPHOS). Over time mt-dys can lead to mitochondrial diseases (mt-diseases), which are surprisingly prevalent and include common diseases such as Alzheimer's, Parkinson's, and diabetes. Unfortunately, the strong impact that intraspecific mitochondrial and nuclear genetic variation has on mt-disease complicates its study and the development of effective treatments. Animal models have advanced our understanding of mt-disease but their relevance to human conditions is often limited by their relatively low nuclear genetic diversity. Many traditional laboratory models also typically have a single mitochondrial haplotype (mitotype), in stark contrast to over 5,000 mitotypes in humans worldwide. The threespine stickleback fish has an evolutionary history that has made it a favorable evolutionary mutant model (EMM) for studying mito-nuclear interactions and possibly mt-diseases. EMMs are species with naturally evolved states that mimic maladaptive human diseases. In threespine stickleback, a period of isolation followed by introgression of the mitochondrial genome from a sister species resulted in the maintenance of two distinct mitochondrial haplotypes which continue to segregate within many populations of wild stickleback. The existence of two mitogenomes segregating in numerous genetically diverse populations provides a unique system for exploring complex mito-nuclear dynamics. Here we provide the first complete coding region analysis of the two threespine stickleback mitotypes, whose mitogenomic divergence exceeds that of other mammalian models for mitochondrial disease and even that between ancient and modern humans. We find that divergence is not uniform across the mitogenome, but primarily impacts protein coding genes, and significantly impacts proteins in Complex I of OXPHOS. The full characterization of these highly divergent intraspecific mitotypes provides a foundation for the development of threespine stickleback as an EMM for mito-nuclear interactions.

Pilot Screening of Cell-Free mtDNA in NIPT: Quality Control, Variant Calling, and Haplogroup Determination

Clinical tests based on whole-genome sequencing are generally focused on a single task approach, testing one or several parameters, although whole-genome sequencing (WGS) provides us with large data sets that can be used for many supportive analyses. In spite of low genome coverage, data of WGS-based non-invasive prenatal testing (NIPT) contain fully sequenced mitochondrial DNA (mtDNA). This mtDNA can be used for variant calling, ancestry analysis, population studies and other approaches that extend NIPT functionality. In this study, we analyse mtDNA pool from 645 cell-free DNA (cfDNA) samples of pregnant women from different regions of Russia, explore the effects of transportation and storing conditions on mtDNA content, analyse effects, frequency and location of mitochondrial variants called from samples and perform haplogroup analysis, revealing the most common mitochondrial superclades. We have shown that, despite the relatively low sequencing depth of unamplified mtDNA from cfDNA samples, the mtDNA analysis in these samples is still an informative instrument suitable for research and screening purposes.