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.

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.
The field of environmental DNA (eDNA) research has been rapidly expanding in recent years, resulting in unprecedented advances in a range of biological monitoring applications. Environmental DNA research provides a non-invasive and cost-effective approach for the study and management of wild populations and invasive species, by using a forensics approach to the extraction and identification of DNA fragments released as organisms travel through and interact with the environment [1][2][3][4][5][6][7] . Environmental DNA analysis is also being applied to issues of human and animal health-for example, in pathogen, parasite and pollen monitoring 1,[8][9][10] . This includes the rapidly emerging field of human eDNA-based pathogen detection from human wastewater. Such approaches developed quickly during the early stages of the  Article https://doi.org/10.1038/s41559-023-02056-2 sea turtle eDNA than human eDNA present (turtle RPTM/human RPTM ratio, 63.50 and 18.50, respectively). The high rate of human eDNA recovered from wild samples (Fig. 1b and Supplementary Table 1) is particularly relevant as our study sites were not directly adjacent to areas of dense human habitation (that is, cities or towns). The relative paucity of human-aligning eDNA reads in the water negative field control and the rehabilitation tank samples (Supplementary Table 1) confirms that the vast majority of human-aligning reads did not originate from contamination during sample processing. In total, 1.8 million paired-end human-aligning reads (300 base pairs (bp) per read pair) were recovered as by-catch from this eDNA study ( Fig. 1b and Supplementary Table 1).
The human Y chromosome is a fast-evolving reduced chromosome, which is divergent in structure and gene content from even humans' closest extant relative, the chimpanzee 25 . Y chromosomes are not shared among all vertebrate species; they are specific to therian mammals. Sea turtles (our target study animal) do not possess Y chromosomes, instead having temperature-dependent sex determination. The human Y chromosome is therefore a useful genome region for confirming the presence of genuine human reads in these complex metagenomic eDNA sequencing samples. Human Y chromosome genomic bycatch reads were detected in all samples (Extended Data Fig. 1b), despite the Y chromosome's small size and the default exclusion of reads originating from human females. Estuarine samples tended to have higher human Y chromosome RPTM values than oceanic samples (Extended Data Fig. 1b).
Having determined that HGB can occur from wildlife eDNA sampling, we next investigated the feasibility of intentional human eDNA recovery from sampling sites of suspected high and low human eDNA release. Species-specific qPCR assays for the quantification of human eDNA revealed that human DNA was readily recoverable from water sites located close to towns (Figs. 1c and 2 and Extended Data Figs. 1c, 2 and 3). These findings were replicated in both subtropical Florida and temperate Ireland. Minimal human eDNA was detected in a mountain tributary (Goldmines River) of the Avoca River, Co. Wicklow, Ireland, close to its source and above the line of human habitation (Fig. 1c, Extended Data Fig. 1c and Supplementary Table 2). In contrast, human eDNA was detected once the Avoca River enters Arklow town, with human eDNA levels in the river increasing as it flows through the town, before being diluted when the river enters the Irish Sea ( Fig. 1c and Extended Data Figs. 1c and 2a,b). More human eDNA was present in the samples within the town even though less than half of the volume of water could be filtered compared with the non-town samples, due to increased filter clogging induced by more turbid water (Supplementary  Table 2). Similarly, no human eDNA was detected from a private well, above the human habitation line ( Fig. 1c and Extended Data Fig. 1c). Non-human eukaryotic eDNA was recovered from all Irish samples, demonstrating successful eDNA extraction (Extended Data Fig. 3). Similar to the Irish samples, human eDNA was readily detectible from water samples taken near the city of St. Augustine, Florida, while no human eDNA was detectible by qPCR from ocean water collected on an incoming tide from the ocean beyond the Matanzas Inlet, Florida ( Fig. 2a and Extended Data Fig. 4).
We next compared human eDNA recovery from beach sand footprints with that from an island with restricted human access. No human eDNA was detected by qPCR in sand samples from a remote area of Rattlesnake Island, Florida, part of the Fort Matanzas National Monument ( Fig. 2b and Extended Data Fig. 5). This part of the island is inaccessible to the general public; our access was facilitated by the US National Park Service. By contrast, human eDNA was readily detectible from beach sand samples recovered from human footprints (Fig. 2b).
Then, we compared the recovery of human eDNA from air in rooms where humans were present with that from air in rooms where humans were absent. No human eDNA was detected by qPCR in negative field controls (eDNA filters left in rooms with humans present pandemic and have already been repurposed for other pathogens such as monkeypox, poliovirus and tuberculosis 1,[11][12][13][14][15] . Environmental DNA has been successfully obtained from a range of sample types including air, soil, terrestrial and aquatic sediments, water (marine, freshwater and wastewater), permafrost, snow and ice cores 10,16,17 .
Environmental DNA research has traditionally relied primarily on targeted methodologies, such as quantitative PCR (qPCR) and metabarcoding-based next-generation sequencing, and early applications focused on bacterial communities 18 . However, continued improvements in deep sequencing technology and novel bioinformatics refinements mean that untargeted shotgun-sequencing-based approaches are becoming feasible (Extended Data Fig. 1a), which more fully capture the true extent of genetic diversity within a sample 1,8,19 . Shotgun sequencing is set to become more labour-and cost-effective than qPCR or metabarcoding in the near future while providing the least biased biodiversity assessments, thus providing the broadest possible presence and abundance information across all taxa. We have recently shown that untargeted shotgun deep sequencing (the direct sequencing of total eDNA with no prior enrichment or selection) can provide both host and pathogen sequence data 8,17 , while also simultaneously capturing all other biodiversity within an environmental sample. Similar to biodiversity assessments, shotgun sequencing of wastewater samples could be applied to monitor all human pathogens simultaneously but would also probably capture a large volume of human genomic data.
While there is a plethora of beneficial applications of eDNA, we postulate that an unintended negative consequence of eDNA approaches might be the capture of human genomic information (human genetic bycatch (HGB); Fig. 1a). Beneficial applications of human-focused eDNA sampling can also be envisaged. Currently, human DNA is rarely (if ever) the intended target of eDNA studies, leaving the field with a lack of specific human-related regulatory guidelines or ethical approvals. Current targeted qPCR and metabarcoding-based eDNA approaches do not recover any substantial human genomic information. However, as eDNA shifts towards shotgun sequencing, potentially large volumes of human eDNA will be retrieved, including sufficient data to identify and phenotype human individuals. Obtaining genetic data from identifiable persons requires informed consent 20 . Legal and ethical frameworks are common in studies involving humans and studies that generate patient data, albeit with continued debate regarding whether such policies are sufficiently rigorous in relation to informed consent, data ownership and data protection [20][21][22][23][24] .
To ascertain whether human genomic DNA could be harvested from eDNA data, we aligned the sequencing data previously generated 8,17 as part of our wildlife and pathogen eDNA projects against the human reference genome. Having demonstrated the occurrence of HGB, we next applied species-specific qPCR to quantify the level of human eDNA in environmental water samples from sites distant from and close to human habitation, from human footprints in beach sand and from occupied and unoccupied room air ( Supplementary Fig. 1). Finally, we applied long-read shotgun sequencing and short-read sequencing human exome enrichment to obtain human-aligning sequences to reconstruct informative human haplotypes (genetic ancestry and mutations) from eDNA.

Results
Human-aligning reads were detected in all samples ( Fig. 1b and Supplementary Table 1) of untargeted shotgun deep sequencing from water and sand eDNA generated for wildlife and pathogen monitoring 8,17 . Furthermore, in some wild (non-rehabilitation tank) water samples, human-aligning reads were detected at levels almost as high as that of our main study species, the green sea turtle 17 (green sea turtle reads per ten million total reads (RPTM)/human RPTM per wild water sample average ratio, 1.39; minimum ratio, 1.03; maximum ratio, 2.54). As expected, the rehabilitation sand and tank water samples had more     Article https://doi.org/10.1038/s41559-023-02056-2 during sampling, but not connected to a vacuum pump) or from vacuum-pumped air from rooms in which no humans were present (Fig. 2c). Human eDNA was recovered from rooms with humans present as participants went about regular working activities (Fig. 2c). Human eDNA was recovered even though air samples from rooms with humans present were collected from a sterile veterinary hospital environment. Differences in the level of human eDNA detected were broadly in close agreement between the two independent species-specific eDNA assays (LILRB2 and ZNF285). This indicates that both are suitable for qPCR-based human eDNA applications (Figs. 1c and 2a,b and Extended Data Fig. 1c).
Having quantified the level of human eDNA in each sample by qPCR, we selected seven samples for Oxford Nanopore MinION sequencing (Supplementary Table 3) to confirm that human genomic information could be derived from intentional human eDNA sampling. MinION shotgun sequencing was conducted with six intentional human eDNA samples (collected in 2022) and our longest-travelled water negative field control sample (collected in 2018, travelled over 904 km). Oxford Nanopore sequencing was selected as it generates longer read lengths (Extended Data Fig. 6a,b) than Illumina sequencing. Without any enrichment (shotgun sequencing), the samples of human eDNA in water, human footprint sand and room air returned thousands of human-aligning reads, while negative field control water and no-human-site sand eDNA samples (from a snake trail) had only 2 to 26 human-aligning reads ( Fig. 3a and Supplementary Table 3). For all human-positive substrates, coverages were relatively even across the human reference genome, with reads from all chromosomes being detected (Fig. 3b). Nanopore sequencing revealed that eDNA is not necessarily fragmented DNA, as even without employing high-molecular-weight extraction methodology or long-read library preparation protocols (short-read nanopore buffer was used), we recovered human eDNA single reads up to 148,969 bp long (air eDNA; 120,998 bp for water and 39,229 bp for sand) ( Fig. 3b and Supplementary Table 4). The longest human-mitochondrial-aligning single read was 16,535 bp, only 34 bp shorter than the full-length mitochondrial reference genome (Supplementary Table 4). The average read length across all human-positive nanopore samples was 1,514 bp (Supplementary Table 4).
To demonstrate that human eDNA could be used for applications beyond mere quantification, we examined known human genetic variants in the nanopore eDNA data to determine whether eDNA-based ancestry and disease susceptibility applications may be feasible. Deletions annotated in gnomAD (v.2.1) 26 could be detected in all three shotgun human-positive eDNA substrates, mainly in water but also to a lesser extent in sand and air samples ( Fig. 4a and Supplementary Table 5). The longest deletion detected within a single read was 40,738 bp, a common copy-number polymorphism in European and Latino populations, located on human chromosome 2 (Fig. 4a). Genes associated with the deletions (deletions within or adjacent to the gene) from the Moultrie Creek B water sample have a range of functions including neuron differentiation, the regulation of double-strand break repair and the regulation of proteolysis. Deletions in or near prominent cancer-associated genes (for example, ALK, LIN28B, PDGFD and WNT7B) were also detected (Extended Data Fig. 7a and Supplementary Table 5). Other structural variant types (insertions and duplications) were also detectable, when analysed with the EPI2ME Sniffles-based structural variant caller (Extended Data Fig. 7b). Human mitochondrial reads from the shotgun nanopore sequencing were also assessed for confirmed mitochondrial pathogenicity-associated alleles (MitoTip 27 and ClinGen 28 ); we detected seven mitochondrial mutations (six from water eDNA and one from air eDNA) associated with a range of diseases, including autism, diabetes, eye diseases and cardiac diseases (Supplementary Table 6).
We next assessed the feasibility of human exome enrichment on water and sand eDNA samples. Five eDNA samples were human exome enriched and then sequenced on an Illumina NovaSeq6000. Three of these samples had also been subjected to long-read sequencing (Supplementary Table 4). Despite exome enrichment, the negative field control samples (NFC sand eDNA Rattlesnake Island site 2 sample 1 and NFC Moultrie Creek) generated very few human-aligning reads (Extended Data Fig. 7c and Supplementary Table 4). For all three human-positive eDNA samples sequenced, exome enrichment increased the proportion of human-aligning reads ( Fig. 4b and Extended Data Fig. 7c). Post-exome enrichment, the water eDNA sample had 48% of sequenced reads aligning to the human genome (20.72 billion human-aligning bases), representing 473× coverage of the targeted human exome regions ( Fig. 4b and Supplementary Table 4).
We then examined the human population-genomic-level resolution present within our eDNA samples. The detected deletions themselves (Fig. 4a) occur at varying allele frequencies in different human populations (Supplementary Table 5). We also conducted haplogroup and haplotype analysis on all sequenced intentional human-eDNA-positive samples: four nanopore samples, three exome-enriched samples and the wildlife-focused sand and water Illumina sequencing with the highest number of human-aligning reads. Every eDNA sample assessed enabled population genomic analysis, with proportions of specific haplogroups and haplotypes varying between samples (Fig. 4c). Generally, haplotypes from European-Indian ancestry, particularly H2a2a1, were the most widely recovered from each substrate type (Fig. 4c). For samples with known participants (footprint and room air), the haplotyping matched the participant profile, while for the water eDNA samples, the results broadly matched the demographics of the area from which the sample was taken (https:// www.census.gov/quickfacts/fact/table/staugustinesouthcdpflorida, staugustinecityflorida/PST045221). Even with no prior enrichment, reads containing information on genetic ancestry and parental origin have been recovered ( Fig. 4c and Supplementary Table 5).
Each sample also contained a complex mix of microbial reads, demonstrating the utility of eDNA (including air eDNA) to simultaneously detect humans, animals and microbes (Extended Data Fig. 8a,b). This included the room-air sample (from a sea turtle hospital), which had airborne pathogenic (tumour-implicated) sea turtle herpesvirus and papillomavirus DNA as well as sea turtle DNA ( Fig. 4d and Extended Data Figs. 8c and 9a).

Discussion
The rapid expansion of the field of eDNA coupled with increasingly cost-effective deep-sequencing technologies promises an exponential increase in the generation of shotgun data from environmental samples over the coming years 1,16,19,29,30 . This heightens the likelihood of unintended large-scale HGB from the general public in eDNA studies. In the present study, we report that HGB was found in all field eDNA samples. These samples had been collected primarily for the detection of non-human species, marine turtles, animal pathogens and metagenomics 8,17 . With no human enrichment prior to shotgun sequencing and with sampling having been conducted in areas of relatively low human habitation densities, we nevertheless inadvertently captured a substantial amount of human genomic data. There is likely to be a continuum across which HGB may be beneficial, neutral or exploitative-for example, samples collected from a popular tourist beach are likely to be less ethically concerning than wastewater monitoring of small, defined, stable populations. Similarly, areas with greater human population numbers are more likely to return more human data. When the sites were intentionally tested, far more human eDNA was indeed recovered from aquatic sampling sites in which water had passed through populated areas. Furthermore, targeted sequencing of human eDNA from complex metagenomic samples was readily achievable with existing exome enrichment approaches. Having demonstrated the propensity of eDNA approaches to capture human eDNA equally well as that of other species, we will now consider the potential beneficial applications and ethical implications of these findings (Box 1). Article https://doi.org/10.1038/s41559-023-02056-2 At a minimum, HGB may complicate approvals for eDNA-based deep-sequencing projects, potentially requiring them to be reviewed by human-focused ethical review boards, even when humans are not the intended target. As eDNA approaches are non-invasive, limited approvals are currently needed, even when investigating endangered species. Added layers of regulation may therefore complicate the implementation and widespread adoption of eDNA approaches, hampering conservation and non-human genome research applications.
Whole human genome (hg38) aligning reads from intentional high-human and low-human eDNA samples (no-enrichment Oxford Nanopore shotgun long-read sequencing) a b      Table 5. b, Percentage of sequenced reads from eDNA metagenomic samples that aligned (minimap2) to the human genome, compared with all reads generated for each sample (regardless of species of origin). Nanopore shotgun sequencing is compared with Illumina exome enrichment. Note that the same male footprint sand eDNA sample and the same high-human water eDNA sample was used for both shotgun sequencing (nanopore) and exome enrichment (Illumina). c, Haplogroup and haplotype analysis of human-mitochondrial-aligning reads for each sequenced (shotgun or exome-enriched) intentional human sample and for the water and sand bycatch Illumina shotgun-sequenced samples (sea turtle focused) with the highest number of human-aligning reads. Pie charts within the same box are from sequencing data for which the same original sample was used (that is, for shotgun and exome enrichment). d, Quantification of green sea turtle and sea turtle herpesvirus (ChHV5) and papillomavirus (CmPV1) aligning reads (minimap2) from the same room-air eDNA nanopore shotgun sequencing data as were used for identifying human-aligning reads (Fig. 3a).
Article https://doi.org/10.1038/s41559-023-02056-2 The capture of non-genetic human information such as speech (eco-acoustic data) or human images (camera traps) has posed similar conundrums for conservationists as HGB. In those cases, codes of conduct, enhanced education and automated data filtering have been suggested to alleviate privacy and consent concerns [31][32][33][34] .
The deposition of sequencing data in a public repository is a prerequisite for scientific publication and open data initiatives. However, if eDNA sequencing datasets also include human genomic information, this may produce a potential conflict. An open ethical question is therefore whether any such data would have to be pre-filtered to remove human sequencing data prior to deposition, similar to the de-identification of patient data 35 . This in turn raises issues relating to transparency and the correct filtering criteria. Furthermore, with whom should the responsibility of policing such a solution lie: journals, data repositories, ethical committees or investigators themselves? Another alternative is attempting to actively block the sequencing of human DNA from eDNA samples. However, a recent study reported that even when human blockers were applied, human reads were still detected 36 .
Current deep-sequencing approaches suffer from a somewhat similar issue as HGB, whereby genomic information from personnel who process samples may be inadvertently sequenced 37,38 . However, in terms of ethical considerations, this inadvertent genome capture from contamination occurs from experienced personnel who are well versed in the technology. Environmental-DNA-based inadvertent capturing of genomic information from the general public is a more complex ethical conundrum, as the genomic information captured comes from individuals mostly unaware of the technology use and oblivious to the fact that their genetic information has been inadvertently obtained. Complicating matters, researchers are already exploring the potential of eDNA for individual-level identification and tracking in wildlife populations [39][40][41] . HGB in eDNA metagenomics studies therefore raises the intriguing question of whether individual humans could also be identifiable from eDNA data; this is likely, given the unfragmented nature of eDNA reported here. This could require eDNA studies to obtain prior informed consent, a near-impossible task. Given the large geographical range over which eDNA can travel 1,6,42,43 (particularly for aquatic samples [44][45][46], it is impractical to imagine that prior informed consent could realistically be obtained. Even with no targeting or enrichment, our MinION shotgun sequencing was able to identify the genetic ancestry within pooled human populations and to identify variants associated with disease susceptibility. The long eDNA reads sequenced here suggest that with targeted enrichment of informative genomic locations, one could achieve individual-level identification even from pooled samples. Given the advanced targeting approaches already in existence for genomic regions of interest (for example, from disease and population genomics research, nanopore adaptive sequencing [47][48][49][50][51] or the exome enrichment used in this study), targeting specific regions or variants from human eDNA samples is already completely feasible.
While our data demonstrate the possibility of human genomic bycatch in wildlife-oriented eDNA studies, they also highlight the feasibility of targeted human eDNA applications. Particularly sensitive or informative regions of the human genome could be intentionally targeted. Human privacy issues are particularly relevant to the increased use of pathogen detection from human wastewater during the COVID-19 pandemic 1,11,12 . Even in relation to pathogens, the legal and ethical implications of wastewater monitoring have not been adequately considered. Wastewater is currently utilized to detect illicit drugs, antidepressants, stress markers and alcohol consumption 52 . Our findings show that such ethical dilemmas need to be considered not only for wastewater sampling but for the field of eDNA as a whole, with HGB feasible from air, aquatic and substrate environmental sampling.
Given the demonstrated feasibility of human eDNA analysis, including disease-associated loci calling, it is likely that beneficial applications may exist for such human-orientated eDNA approaches (Box 2). For example, one could utilize wastewater or air eDNA-based sampling to correlate the level of pathogens with the abundance of susceptibility loci in a given population 1,13,15,53,54 . This may be especially feasible given the ever-increasing portability of DNA sequencing/ surveillance technology 1,29 . It is also immediately feasible to consider intentional human and pathogen eDNA sampling from specific sources, such as substrate footprints or filtered-air eDNA (as we report here) from specific locations-for example, rooms in the home, public indoor or outdoor spaces, hospital wards, or even special diagnostic rooms designed to filter the host and pathogen eDNA shed by a single patient over a short period. Here we report that room-air sampling in a clinical setting (a veterinary hospital) recovered host (animal patient), human (staff) and vertebrate viral pathogens, suggesting future novel monitoring approaches in human and animal medical settings and beyond. Alternatively, water eDNA collected before wastewater leaves the home could be utilized, linking to continual health monitoring and chronic disease management initiatives 55 . This is particularly pertinent to routine continual health monitoring for somatic mutations, which arise later in life and can be drivers of life-threatening diseases such as cancer. Human eDNA approaches can be combined with advances in the ongoing genomic medicine revolution 24 to enable novel applications, as evidenced here by the ability to detect specific human disease-associated mutations from eDNA and to enrich for human genomic regions of interest (for example, exome) and simultaneously detect microbes, including pathogenic viruses.
Haplotyping and phylogenetic analysis of individuals (whale sharks, sea turtles and kakapo) from eDNA samples have already been achieved, although to date these have tended to be samples dominated by a single individual 17 this list and highlights how haplotyping can be achieved from complex multi-individual (metagenomic) samples with diverse population structures. Given its less biased nature, human eDNA may also help redress the balance in relation to the lack of sequenced genomes from diverse populations and underrepresented rare human alleles in genomic databases. Traditionally, these databases have been skewed towards populations of European ancestry 57 . Separately, genetic analysis is already being applied to combat the illegal wildlife trade by identifying the source animal populations 58,59 . Similarly, it may be possible to adopt such approaches in targeting human eDNA recovered from illegal wildlife parts (poacher/trader handling) to broadly identify the geographic human populations through which the illegal parts transited. Human eDNA could also be employed to identify sites of archaeological importance; such possibilities include sampling of remote stagnant waterbodies (for example, bogs) to identify undiscovered sacrificial sites, or eDNA sampling to help in the recovery of more recent human remains. Even search-and-rescue missions could use eDNA techniques if currently in-trial air eDNA sampling by drones is successful. Recently, it has been postulated that human eDNA extracted from the air could also have novel forensic and criminal investigation applications [60][61][62] . Our human population genomics assessment of airborne DNA is an early step in this direction, providing proof-of-concept of the recovery of shotgun-sequencing-quality human DNA from room air. Both ancient and contemporary DNA have their own established ethical procedures, although in those disciplines too, the pace of advancement has outpaced the dialogue about research ethics [63][64][65][66] . Vertebrate air eDNA is a recent field of study, and improved capture devices is an active area of research from which human airborne DNA applications could benefit 36 . Water-based eDNA approaches have also been postulated to be beneficial for forensic investigations 67,68 . Forensic science has well-established ethics procedures, although issues pertaining to racial discrimination persit 69,70 . Taken together, forensic applications can be envisaged for air, soil and water eDNA approaches, especially in light of the intact non-fragmented nature of eDNA reported here.
While benefits of human eDNA analysis may exist, there could conceivably be more worrying applications of such technology. The accumulation of population-level genomic databases is currently highly desirable, being a valuable research and commercial commodity 71,72 . While some projects maintain the genomic data and subsequent findings in public ownership, the lucrative economic potential of such databases means that a host of private companies have also been rapidly accumulating them (for example, DNA ancestry companies and genomic medicine companies) and selling access 72,73 . Countries and corporate entities are racing to create ever-larger pan-genomic patient/population datasets. Examples include the 100,000 Genomes Project and the new Genome UK Strategy (https://www.gov.uk/government/publications/ genome-uk-the-future-of-healthcare); the Biobank project, which aims to sequence half a million genomes (https://www.sanger.ac.uk/collaboration/uk-biobank-whole-genome-sequencing-project/); and the National Institutes of Health All of Us Research Program (https://allofus. nih.gov/), a project to sequence the genomes of one million US citizens. Such pan-genomic activities have the potential for great public good, including advanced medical and pharmaceutical applications, but they have been an ethical minefield regarding ownership, data protection, insurance coverage and privacy issues 24,72,[74][75][76] . Examples from digital databases also clearly show that wherever valuable databases exist, there is an ever-present temptation to exploit the data for financial gain, whether legally or otherwise. The Cambridge Analytica scandal, which involved complex harvesting of data from Facebook profiles 77 , is one such example, and Google's harvesting of medical records with individual identifying information still intact is another 78 .
Human eDNA applications present similar concerns: they could be employed for the surveillance of individuals, minority groups (genetic ancestry) or genetically driven disabilities, or to obtain genomic information from local populations without their knowledge or consent, including from 'valuable' genetically diverse indigenous populations (uncontacted tribes, particular ethnic groups and so on). Such scenarios expand the current issues regarding the commodification of genomic information, which is already a particularly acute concern in relation to indigenous peoples 72,79-81 . These novel human eDNA applications may also enable unscrupulous (though possibly not technically illegal) entities to perform covert mass capture of genomic data from populations, to further expand population-level genomic databases. Furthermore, such resources can then be utilized for enhanced monitoring of individuals for less savoury purposes. Genetic/genomic surveillance is a serious ethical concern, with documented human rights abuses having already occurred whereby national DNA databases were used with other surveillance data to monitor minority populations 70,82 . The application of human eDNA approaches could further undermine genetic consent, limiting the ability of threatened minorities to withhold their genetic information. In the future, human eDNA could also be utilized to determine whether members of a genetically distinct group were present in a given population-for example, through wastewater monitoring or air filtering at checkpoints, in urban areas or in private dwellings. Such potential is particularly chilling given the propensity of humans to carry out ethnic persecution and genocide throughout our history. Many of the ethical issues raised here are more likely to arise when sampling sites are closer to human urban areas. Similarly, issues pertaining to individuals or ethnic groups are more likely to be pertinent in areas with small, stable populations of humans contributing to the recovered eDNA-for example, a specific watershed or portion of a sewage system. It is clear that in the short term the convergence of massively powerful deep-sequencing technologies with improvements in ecological eDNA sample capture approaches raises serious questions about human genomic bycatch, requiring discussion and regulatory consideration from scientific, ethical, regulatory and social perspectives. Regulators, researchers, funders and other stakeholders should develop responses to the ethical implications of HGB and intentional human eDNA applications. Such planning should be initiated immediately, pre-empting the technology becoming even more widespread, affordable and entrenched, be it for beneficial or exploitative applications. Such ex ante planning is crucial for ensuring that laws and ethics stay ahead of emerging technology. Conversely, these same eDNA approaches can open up novel beneficial applications in areas from human health to criminal forensics.

Sample collection, DNA extraction and sequencing
All laboratory procedures from sampling through final analysis (sequencing or qPCR) were conducted in a way that minimized any human DNA contamination from investigators. Although the Illumina HiSeq eDNA samples had also been sampled and utilized for sea turtle and pathogen research, the original study design for these samples included assessment of the human-aligning reads from these shotgun data. Therefore, all appropriate precautions to avoid contamination were taken. This included no contact of investigators or samplers with the substrates being sampled (water or sand), new nitrile gloves being used for each sample collection, frequent glove changes throughout sample processing, cleaning of equipment and benchtops with bleach prior to use, and negative field controls being treated identically to genuine field samples throughout all processes from extraction through qPCR/sequencing (see below). All eDNA samples were extracted in labs that do not process human samples: a sea turtle lab in Florida and a chick/mouse developmental biology lab in Ireland. The Irish qPCRs were set up in a laminar flow cabinet after a protracted 1 h exposure of UV of the cabinet and equipment. It was originally postulated that the sea turtle rehabilitation tank water sample (2017, HiSeq) would contain more human DNA than the field water samples, due to hospitals staffs' interactions with the tank water and sea turtle patients. For all subsequent eDNA sampling, we had already established that the 2017 tank sample contained less human eDNA than the 2017 field samples. We therefore continued to observe strict procedures to avoid investigator contamination in all subsequent sampling in order to be able to investigate HGB. While we intended in advance to assess the level of human eDNA recovered from the 2017 and 2020 samples, no human factors (such as population density) were considered when selecting sampling sites. Site selection was purely based on the wildlife eDNA considerations of these sites-that is, site selection was not biased towards areas with high human densities. Conversely, site selection for all 2022 samples was intentionally directed towards sites with high or low human perturbations (primarily based on human population density).
Previously extracted 8 green and loggerhead sea turtle DNA from tissue samples was used for species specificity qPCR tests. DNA was extracted from tissue using a Qiagen DNeasy Blood and Tissue kit (Qiagen, cat. no. 69504) according to the manufacturer's instructions. SH-SY5Y cells (ATCC, cat. no. CRL-2266) were gifted by the Loesgen lab (University of Florida), and no human cells were processed in the same lab as the eDNA samples. DNA was extracted using a Qiagen DNeasy Blood and Tissue kit according to the manufacturer's instructions. SH-SY5Y genomic DNA was used for generating standard curves. To avoid any potential contamination, the standard curves were run only after all eDNA qPCRs had been completed.
HGB samples. Seawater (rehabilitation tank and oceanic) and beach sand eDNA samples were collected between 2017 and 2021 and sequenced as part of our sea turtle and sea turtle pathogen research 8,17 . All samples originated from Florida, US 8,17 . The tank sample is a pooled 8,17 sample containing eDNA extracted from seawater from five rehabilitation tanks at the University of Florida's Whitney Laboratory for Marine Bioscience and Sea Turtle Hospital. The four main tanks housed juvenile green sea turtles and were 240 cm in diameter and had a full volume of 2,270 l, and a smaller loggerhead post-hatchling tank had a volume of 480.75 l. The separately extracted eDNA from each tank was pooled in equal volumes prior to library preparation.
Intentional human samples. River, estuarine, seawater (oceanic) and beach sand samples were collected between May and July 2022 (Supplementary Table 2). Negative field control samples of water and sand were also collected and processed as per the study samples. For negative field control water sampling, 1 l of MilliQ water (Florida) or 1 l of Qiagen Nuclease-free water (Ireland, cat. no. 129117) was transported from the laboratory to the rehabilitation or wild sampling locations and stored in a cool box with the environmental samples to monitor for potential contamination during sampling, transportation and processing. For sea turtle negative field control sand sampling, 50 ml of beach sand was collected away from suspected turtle presence (that is, away from sea turtle tracks or obvious human activity) on each sampling trip. For human negative field control sand sampling, 50 ml of beach sand was collected away from suspected human activity (such as footprints) on each sampling trip and from a restricted-access location on Rattlesnake Island, part of the Fort Matanzas National Monument managed by the US National Park Service. The water and sand negative field controls were filtered and extracted alongside the other collected sand and water samples from each sampling trip and subjected to the same next-generation sequencing conditions and qPCR conditions (intentional human samples). The standard volume of seawater filtered (0.22 μm pore Millipore Sterivex-GP Pressure Filter Units (Merk Millipore, cat. no. SVGPL10RC)) was 500 ml for each DNA sample 8,17 . Samples of less than 500 ml (Supplementary Tables 1 and 3) were a lower volume due to debris clogging the filter and preventing larger volumes being filtered. One sample ('tidal pool') had 1 l of seawater filtered. For sand eDNA, a 50 ml tube was filled with sand from each sampling event, with 10 ml of this sand used per individual eDNA extraction 17 . Human-present air samples were collected from a 280 ft 2 room while the participants went about their daily work activities (that is, they could enter and exit the room throughout the sampling period), with a maximum of the same six participants using the room for a portion of the sampling period. The room was air-conditioned (outside air) and had an external door that was opened and closed, and occasionally left open for certain work procedures. Negative field control samples of air were also collected and processed as per the study samples. For air eDNA, two types of negative field controls were collected: (1) a filter kept in the room being sampled for the duration of sampling, but with no air being pulled through it (that is, no pump), and (2) air filtered (with a pump) from a room with no humans present at the time of filtering or during the previous 24 h. Human-related sampling was conducted with University of Florida Institutional Review Board (IRB-01) ethical approval under project number IRB202201336, with all participants providing informed consent. Four participants (three Article https://doi.org/10.1038/s41559-023-02056-2 female and one male) provided sand footprint samples, and six participants (five female and one male) provided room-air samples (pooled room air). Participation was on a voluntary basis, and the participants received no compensation.
Prior to filtration and between every sample, the laboratory surfaces and filtration equipment (all standard laboratory equipment involved in the washing and filtration process, as well as the filtration pump itself) were disinfected with 70% ethanol, and the sampling equipment (collection bottles) was disinfected (washed thoroughly) with 10% bleach and rinsed thoroughly with deionized water.
Sand. 1X TE-IDTE pH 8.0 1X TE Solution (Integrated DNA Technologies, cat. no. 11-05-01-09) was added to each individual sand sample in individual 50 ml Falcon conical centrifuge tubes, at approximately two times the volume of sand (10 ml of sand and 20 ml of 1X TE). The samples were shaken gently by hand and then set on a rocking platform for 1 h at room temperature, with additional gentle shaking by hand every 15 min. The samples were rested until sand had sunk to the bottom of each tube; then, the supernatant was immediately pipetted into a 60 ml sterile BD luer lock syringe (Fisher Scientific, cat. no. 136898). The samples were then hand-filtered using 60 ml BD luer lock syringes through 0.22 μm Sterivex-GP Pressure Filter Units (Millipore, cat. no. SVGPL10RC) and capped with B.Braun luer lock caps (Medline, cat. no. BMGTMR2000B). Finally, 740 μl of Buffer ATL and 60 μl of Proteinase K from a Qiagen DNeasy Blood and Tissue Kit (Qiagen, cat. no. 69504) were added to each sample, and they were placed in 50 ml Falcon conical centrifuge tubes in a rolling incubator overnight (24-26 h) at 56 °C.
Water. The water samples were pumped (by hand (in Ireland) or electronically (in Florida)) through 0.22 μm Sterivex-GP Pressure Filter Units (Millipore, cat. no. SVGPL10RC) and capped with B.Braun luer lock caps (Medline, cat. no. BMGTMR2000B). Hand-pumping was with sterile 60 ml BD luer lock syringes. Electronic pumping was done using a GeoTech Peristaltic Pump Series II. Then, 740 μl of Buffer ATL and 60 μl of Proteinase K from a Qiagen DNeasy Blood and Tissue Kit were added to each sample, and they were placed in 50 ml Falcon conical centrifuge tubes in a rolling incubator overnight (24-26 h) at 56 °C.
Air. For air eDNA sampling, room air was passed through 0.22 μm pore Millipore Sterivex-GP Pressure Filter Units (Merk Millipore, cat. no. SVGPL10RC), using a Welch vacuum pump (2019LD-4112) or a GeoTech Peristaltic Pump Series II. Environmental DNA was extracted 8,17 as for the water and sand samples, except that only 20 μl of Proteinase K was added per filter and the 56 °C ATL Buffer and Proteinase K incubation was conducted for 1 h.
DNA extraction and sequencing. For all three sample types (sand, water and air), after the 56 °C incubation, the solutions of Buffer ATL (after water or air filtration, or after sand washing and filtration), Proteinase K and eDNA were transferred from the Sterivex-GP Pressure Filter Units to 2 ml microcentrifuge tubes using 10 ml BD Slip Tip Sterile Syringes (Fisher Scientific, cat. no. 14823434). DNA was then isolated using a modified Qiagen DNeasy Blood and Tissue Kit protocol 8,17,83 . Following incubation, equal volumes of 800 μl (sand and water) of AL Buffer and 800 μl (sand and water) of ice-cold ethanol were added to each sample, and they were vortexed vigorously and microcentrifuged after each addition. For the air samples, 400-500 μl of each solution was used (as a lower volume of ATL solution is recovered from the initially dry Sterivex filters). Each sample was loaded into a DNeasy spin column and centrifuged at 6,000g for 1 min, with flow-through being discarded after each spin (repeated until the entire contents of each sample were spun through the spin column). Then, 500 μl of Buffer AW1 was added, and the samples were centrifuged at 6,000g for 1 min. Next, 500 μl of Buffer AW2 was added, the samples were centrifuged at 16,000g for 3 min, flow-through was removed and they were spun for an additional 1 min at 16,000g. DNA was eluted with 70 μl of AE Buffer (incubated at 70 °C before being added to the spin column), incubated on the column at room temperature for 7 min and centrifuged into a 1.5 ml microcentrifuge tube at 6,000g for 1 min. DNA concentration was measured on a ThermoScientific Nanodrop 2000 Spectrophotometer (Fisher Scientific), and the samples were stored at −20 °C until qPCR or shotgun sequencing.
Library preparation and Illumina shotgun sequencing were conducted at the University of Florida's Interdisciplinary Center for Biotechnology Research Core Facilities. Four water eDNA samples (2017) were sequenced on an Illumina HiSeq 3000, and all subsequent Illumina samples (seven water and ten sand eDNA) were sequenced on an Illumina NovaSeq 6000 (Supplementary Table 1) 8,17 . All Oxford Nanopore samples were sequenced on a MinION in the Duffy lab at the University of Florida's Whitney Laboratory for Marine Bioscience. This device was not previously used for any human samples. We used a personal Min-ION sequencer instead of the previous core facility high-throughput Illumina sequencer as there were some human-aligning reads in our Illumina negative field control sample, although these were 13 to 38 times fewer human reads than were recovered from environmentally derived eDNA samples (Supplementary Table 1 45 Mb of exonic content) exome capture analysis, according to the manufacturer's instructions. Briefly, eDNA was fragmentized, adapters ligated and amplified for nine cycles and purified with Ampure XP beads (Beckman Coulter, cat. no. A63881). The five samples were pooled with equal volume for one exome hybridization according to the user guide, and the probe hybridization was performed at 95 °C for 5 min, one cycle of 1 min each, starting at 94 °C for the first cycle, then decreasing 2 °C per cycle for 18 cycles, and a hold for 90 min at the final temperature. After the exome hybridization, the probe capture, wash and elute was performed, followed by library enrichment with ten cycles of PCR amplification. Exome sequencing was conducted on an Illumina NovaSeq 6000 S4 flow cell for 2× 150 bp cycles aiming for 50 million reads per sample (with negative field controls expected to return fewer reads due to a lack of human eDNA).

Bioinformatic analysis
All bioinformatic tools were utilized using the default parameters, unless otherwise stated. The Galaxy platform (https://usegalaxy.eu/) was used for bioinformatic analysis 84,85 , with NanoGalaxy also used for nanopore sequenced data 86 7 (refs. 104,105) was used to detect structural variants in each sample, using the specific settings 'non-germline' and 'minsupport = 1' to increase sensitivity to single split-read resolution. The results of individual runs were used as inputs for a second, combined rearrangement calling by Sniffles' multi-sample input support. To screen for known human deletions, outputs were processed using the vcfR (v.1.12.0) library in R 109 . Deletion call regions were matched against the gnomAD v.2.1 structural variants database 26 , using a lenient distance threshold of up to 5% of the total event length around both the 5′ and 3′ deletion break-end positions. We manually verified the resulting copy number variants' ONT sequencing support using IGV 110 . Deletion hits were visualized on GRCh37 chromosome maps using the chromoMap (v.4.1.1) library in R 111 . Functional clustering of deletion gene hits was performed using HumanBase (https://humanbase.net/) functional gene network analysis 112 , accessed 21 October 2022.

Quantitative PCR
Two human Applied Biosystem pre-validated Taqman Gene Expression qPCR assays directed against the LILRB2 gene and the ZNF285 gene (assay IDs Hs01629548_s1 and Hs00603276_s1, respectively) were selected for use as species-specific human assays, on the basis of having no cross reactivity with over 27 other species from mice to plants (https://www.thermofisher.com/order/genome-database/; mouse, rat, Arabidopsis, C. elegans, fruit fly, bovine, dog, Chinese hamster, goat, white-tufted-ear marmoset, guinea pig, zebrafish, horse, chicken, soybean, cynomolgus monkey, sheep, rabbit, rice, rhesus monkey, baker's yeast, fission yeast, pig, bread wheat, wine grape, western clawed frog and maize) and having both primers and probe within a single exon (that is, detect DNA). We also showed that these assays did not cross-react with green sea turtle or loggerhead sea turtle DNA (Extended Data Fig. 9b). A pan-eukaryotic 18S rRNA gene (Applied Biosystem, 4352930E) pre-validated Taqman Gene Expression assay, which also has both primers and probe within a single exon (that is, detects DNA), was used to quantify the total level of pan-eukaryotic DNA in each of the Irish samples. Green sea turtle species-specific qPCR from air eDNA was conducted as previously developed and validated assays targeting the 16S rRNA gene, including the use of a 499 bp synthetic gene fragment for standard curve generation 8,17 (forward primer, TGCAAAAGCGGGAATAACAC; reverse primer, TCGCCCCAAC-CAAAAATATAG; FAM labelled (ZEN and Iowa Black double quenched) probe, CAACTATCTATACCCACTCACTCTAAGGACCTATAA (synthesized by Integrated DNA Technologies)). For the green sea turtle 16S rRNA synthetic gene fragment sequence, see Supplementary Table 7.
The qPCR reaction mixtures were performed on 384-well plates in a total volume of 10 μl per well: 5 μl of TaqMan Fast Advanced Master Mix (Fisher Scientific, cat. no. 4444557), 3.5 μl of nuclease-free water (Fisher Scientific), 0.5 μl of the respective assay (primer/probe mix, manufacturer-supplied concentration) and 1 μl DNA template (or, for no-template controls, an additional 1 μl of nuclease-free water per well). Depending on the sample type and the sample volume available, each biological sample was run in three to six technical replicates (three technical replicates for tissue samples and four to six technical replicates for eDNA samples). Negative field controls had the same number of technical replicates run as their corresponding eDNA samples. No template controls were run in triplicate on every qPCR plate. qPCR reactions were performed on an Applied Biosystems QuantStudio 6 Pro (Florida samples) or an Applied Biosystems QuantStudio 7 Flex (Irish samples) with the following cycling parameters: 95 °C for 20 s for one cycle, followed by 45 cycles of 95 °C for 1 s and 60 °C for 20 s. The qPCR results were plotted with BoxPlotR 113 (http://shiny.chemgrid.org/ boxplotr/) with every datapoint displayed. Tukey whiskers (extending to data points that are less than 1.5× the interquartile range away from the first/third quartile) were utilized for every box plot. One box is graphed per single sample, consisting of all qPCR technical replicate wells for that sample. Biological replicates are denoted by the letters A-D at the end of the sample name. Biological replicates are not pooled on any box plots; each sample is denoted by its own box.  Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/ licenses/by/4.0/. Oxford Nanopore eDNA, N50 values for each total library (no enrichment Oxford Nanopore shotgun long read sequencing)

N50 length, kb
Extended Data Fig. 6 | Oxford Nanopore long read shotgun sequencing study of eDNA samples post human DNA qPCR-based quantification. 6a) Longest individual human reference genome (T2T) aligning read (in kb) generated from each nanopore sequenced sample. 6b) N50 length (in kb) from each nanopore sequenced sample, for total sequenced reads, that is pre any reference genome alignment. Note, negative field controls had reduced overall numbers of reads (not reflected in N50 values) reflecting a lack of DNA in the sample (Supplemental Table 3).
Each qPCR reaction is a 10 μl reaction containing 1 μl of extracted genomic DNA template. Note that both assays were species-specific with no cross-reactivity (no amplification) in either sea turtle tissue sample. Loggerhead and green sea turtles can be added to the list of 27 other species (from mice to plants, see methods section) experimentally validated with which these human assays do not cross-react.