Human Oncogenic Epstein–Barr Virus in Water and Human Blood Infection of Communities in Phayao Province, Thailand

: Water can contain pathogenic viruses. Many studies on RNA virus sources have shown that water can transmit them. However, there are few reports on pathogenic DNA virus transmission through water, such as adenovirus, which pose a widespread public health risk. Therefore, this study aimed to show waterborne viral transmission by detecting viruses in pooled human whole blood samples, tap water, and natural water from Mueang District, Phayao Province, Thailand, using a metagenomic approach. Viral prevalence in whole blood samples was measured by polymerase chain reaction (PCR) and quantitative PCR (qPCR), and environmental factors that affect viral infection were assessed. Metagenomics results showed that Epstein–Barr virus (EBV) members were among the prominent cancer-associated oncogenic DNA viruses detected in human blood and all water types similar to the EBV reference sequence (NC_007605). There were 59 out of 813 (7.26%) human whole blood samples that were positive for EBV DNA based on PCR and qPCR for the EBNA-1 and EBNA-2 genes. Water-and blood-borne human oncogenic EBV should be a concern in tap water treatment and blood transfusion in patients, respectively. Therefore, the detection of EBV in water suggests that transmission via water is possible and should be investigated further.


Introduction
Rural residents have high disparities, even in access to clean water, which is a fundamental factor in life [1][2][3].Poverty, poor edification, low health knowledge, lifestyle, insufficient sanitation, poor nutrition, exercise, and environmental factors (such as clean water) have been identified as primary factors leading to the high infectious disease incidence among women in developing countries [4,5].Similarly, antibiotic use was larger in rural (23.6%) than in urban (20.2%) areas compared [3].The International Agency for Research on Cancer (IARC) estimates that 15-20% of cancers are associated with infectious agents [6].
Water 2023, 15, 323 2 of 18 Several DNA and RNA viruses from various cancer-associated viral families are found in water [7].A leading problem is an increasing trend in oral squamous cell carcinoma (OSCC) in patients aged under 45 years (20-34 years) in the Netherlands between 1989 and 2018, with an annual percentage change of 2.4% (95% confidence interval (CI): 1.1-3.8%)[8,9], associated with oncogenic DNA viruses.Ninety-eight studies found that colorectal, breast, kidney, pancreas, and uterine cancers are increasing in younger age groups; although, the mechanisms remain unclear [10].This study is concerned with the quality of water that might be a source of carcinogens, such as viruses, which can cause human disease and many cancers.
Plant, human, and animal viruses can be transmitted through water.Viruses isolated from trees, soil, and water from the abovementioned areas may originate from infected plant debris.Human and animal waste has also been identified as a possible source of enteric viruses in water [11][12][13].About 25% of the global population consumes fecally contaminated water, which may contain bacteria, protozoa, and viruses that can cause many diseases in humans [14], and pose a public health risk [15].Tap-water-borne viruses such as norovirus [16] and torque teno virus (TTV) are of particular concern since many cause disease [17].Human adenovirus, rotavirus, and hepatitis A virus (HAV) have been found in treated wastewater treatment [15,18].Adenovirus, astrovirus, norovirus, hepatitis E virus (HEV), HAV, coronavirus, poliovirus, coxsackievirus, echovirus, and rotavirus have been found in drinking water [19].The conventional methods used to detect waterborne viruses include polymerase chain reaction (PCR), quantitative PCR (qPCR), and nested PCR [15].
The oral route is the primary route of EBV transmission.EBV is known to be transmitted via saliva but rarely spreads through semen, blood, or organ transplantation [23,24].EBV can be spread by using objects, such as a toothbrush, drinking glass, or eating utensils, which an infected person recently used or by sharing their food [32].Nevertheless, organ transplantation and blood transfusion can cause EBV spread, and blood transfusion recipients should be considered at particularly high risk [32].Half of all mothers and one-week-old newborns were EBV positive based on nested PCR and 60.6% based on immunoglobulin G (IgG) [33,34].EBV infection was not observed before 3 months of age but slowly emerged thereafter, reaching a cumulative rate of 1.7%, 11.6%, and 21.5% at 6, 12, and 14 months, respectively [35].The transmission routes leading to EBV infection between three months and adulthood other than human-to-human transmission remain unclear, and it is unknown whether it can be transmitted through environmental sources such as water.
EBV has a double-stranded linear DNA genome containing approximately 100 genes surrounded by a protein capsid.Latent EBV nuclear antigen (EBNA) genes include EBNA-1, EBNA-2, EBNA-3A, EBNA-3B, EBNA-3C, EBNA-LP, LMP-1, and LMP-2A/2B [36,37].EBV infection and the associated immune response occur mainly during the first years of life but may also occur during adolescence or adulthood [38,39].EBV-positive antibodies are present in the serum of >90% of children and adults worldwide.Another study in Thailand found that the prevalence of anti-EBV IgG was 34.9-95% [40,41].However, serum antibody IgG levels were unrelated to the disease [40,42].The prevalence of EBV infection in Thailand using specific immunoglobulin M (IgM) antibodies was 11-13%, representing acute or primary infection associated with EBV DNA.An analysis of anti-EBV IgM and EBV DNA in 80 serum samples from children aged <2 years found that three (3.8%) were positive for IgM and none were positive for EBV DNA using PCR [40,41].
EBNA-1 DNA was detected in 72.5% of whole blood samples from Indonesian NPC patients using qPCR [30].EBNA-1, EBNA-2, and LMP-1 genes were detected in 56% of whole blood samples from Chinese HIV patients tested for EBV DNA using nested PCR [43].The LMP-1 gene was detected in 61% of healthy blood donors in Qatar [44].In addition, the EBNA-2 gene was detected in 52.6% of Qatari samples tested using qPCR [45].Moreover, 9.5% of whole blood samples were positive for EBV DNA in diffuse large B cell lymphoma [46].The prevalence of EBV-positive diffuse large B-cell lymphoma in East Asians (8.7-11.4%)was higher than in Western countries (5%) [32,33].A study of blood donors in Ouagadougou, Burkina Faso, found a lower EBV level of 5.1% using qPCR [47].The PCR-based prevalence of EBV in normal oral exfoliated cells (the most commonly found EBV infection) was 19% in Thailand, lower than in other countries [48].In situ hybridization, Southern blotting, and dot blotting have all been used to diagnose and monitor primary EBV infection.Measuring EBV DNA in blood samples by qPCR is important for monitoring EBV-associated disease [32,[49][50][51].Molecular techniques such as qPCR can be used to detect EBV in various samples, such as serum, whole blood, unfractionated blood, tissue biopsy, and peripheral blood mononuclear cells (PBMCs) [49].A good correlation between EBV DNA loads in whole blood and PBMCs was found by qPCR [52,53].Few studies have investigated EBV viremia in healthy individuals [47].
Environmental factors, including sunlight/vitamin D, smoking, and body mass index (BMI), affect EBV infection [38].EBV infection is not recognized as a national problem in Thailand due to insufficient knowledge.Information on EBV pathogenesis can be found only in textbooks and among a small group of researchers in Thailand.However, the increasing incidence of many cancers has led to appreciable data directly associating cancer with EBV infection.The presence of EBV in water may be associated with human blood transmission, which poses a widespread public health risk.
Phayao Province, Thailand, is a rural small province with a population of 472,356 in 2019.Mueang District, in the center of Phayao Province.Mae Ka and Wiang Subdistrict in the center of the Mueang District.Wiang Subdistrict residents use tap water from the Kwan Phayao lake (19 • 10 32.5 N 99 • 52 04.0 E), while those of the Mae Ka Subdistrict (No. 1) use tap water from Huai Na Poi Brook (19 • 00 59.2 N 99 • 53 01.9 E) and also Kwan Phayao lake.Water from Kwan Phayao is contaminated by sewage from communities, livestock, and hospitals.Nevertheless, the Provincial Waterworks Authority sources tap water from Kwan Phayao.In contrast, tap water sourced from Huai Na Poi, originates from a deep and wide forest basin, with water flowing in from the mountains and is uncontaminated by humans.However, it is relatively cloudy and smells due to humus deposition.Because tap water from Huai Na Poi is directly produced by the villages, it has a lower quality, evident from its turbidity.
Therefore, this study aimed to demonstrate that water is a source of viral transmission by detecting viruses in pooled human whole blood, tap water (Mae Ka and Wiang Subdistricts), and natural water (Huai Na Poi and Kwan Phayao) in the Mueang District, Phayao Province, Thailand, using metagenomics.Tap water from the Mae Ka and Wiang Subdistricts was collected for quality evaluation.The prevalence of oncogenic DNA (EBV) viruses in whole blood samples was determined using PCR and qPCR.Environmental factors that affect viral infection were also evaluated.

Water Sampling Methods
Tap water from Mae Ka Subdistrict (No. 1) and Wiang Subdistrict was collected from Mueang Phayao District, Phayao Province (5 L for chemical analysis and 1 L for microbial analysis).Water quality was tested in duplicate by chemical analysis to detect turbidity, color (Pt.Co), pH, total dissolved solids, total hardness (calcium carbonate [CaCO 3 ]), chloride (Cl − ), fluoride (F − ), nitrate (NO 3− ), nitrite (NO 2− ), arsenic (As), cadmium (Cd), chromium (Cr), copper (Cu), iron (Fe), lead (Pb), manganese (Mn), mercury (Hg), sulfate (SO 4 2− ), and zinc (Zn).Microbial analysis was performed for Escherichia coli (E.coli) and total coliform bacteria.The Provincial Waterworks Authority's standard value was used as a reference.We used an in-house testing method (TE-CH-037) based on standard methods for examining water and wastewater: APHA, AWWA WEF, 23rd Five L each of tap water (Mae Ka Subdistrict (No. 1) and Wiang Subdistrict) and natural water (Huai Na Poi and Kwan Phayao) were collected from Mueang Phayao District, Phayao Province, for water filtration and metagenomic analysis.Tap water was collected by swabbing the faucet with alcohol wipes, turning on the water, letting it run for about 1 min, and then rinsing the container a few times before collecting the sample.Natural water samples were collected from under the surface in sterile glass containers and stored on ice for future metagenomic analysis.

Water Filtration
A total of 5 L of each water sample was filtered using a Rocker 300 vacuum pump (Waller Chemical Thailand Co., Ltd., Bangkok, Thailand) using MF-Millipore pore size 0.025 µm MCE filter nitrocellulose membranes (MF-MilliporeTM, 47 mm; Merck Millipore Ltd., Tullagreen, Carrigtwohill, Co. Cork, Ireland) and on top with 0.45 µm membrane filter.Next, DNA was extracted from all microorganisms, including viruses, on the filter membranes using DNA extraction from filter membranes method.It should be noted that the smallest viruses are about 0.02-0.4µm in diameter (EBV is approximately 0.12-0.18µm in diameter) [54][55][56][57].Microfiltration is generally considered appropriate for quantifying microbes and viruses 0.025-10.0µm [58,59] in water samples [60,61].

Specimens
Eight hundred thirteen human whole blood samples were collected from donors living around the natural water sources in the Mae Ka and Wiang Subdistricts, Phayao Province.Samples were collected from individuals aged 3-90 years, along with information on their sex, health status, and life history, including congenital disease, BMI, family cancer history, exercise, alcohol consumption, smoking, secondhand smoke, drinking water source, cleaning of water for consumption, and type of water used for brushing teeth, serving spoons, and eating fresh fruit and vegetables.The sample size was calculated as N = Z 2  1−a P(1 − P)/d 2 for an EBV prevalence (P) of about 2% [62]: Z = 1.96 for a 95% confidence level and d = 0.01.This study was approved by the Committee on Human Research Ethics in Health Sciences and Science and Technology, University of Phayao (Mae Ka, Thailand; 1.3/023/63).All procedures involving human participants performed in the study were in accordance with the ethical standards of the Declaration of Helsinki, the Belmont Report, the Council for International Organizations of Medical Sciences guidelines, and the International Conference on Harmonization in Good Clinical Practice.Informed consent was obtained from all subjects and/or their legal guardian(s).

DNA Extraction 2.4.1. DNA Extraction from Filter Membranes
Filter membranes were used for total DNA extraction and metagenomic analysis.Cell lysis used a cell lysis solution 4 mL soak and mix by vortex (10 mM Tris-hydrochloric acid (pH 7.8), 5 mM ethylenediaminetetraacetic acid (EDTA), 0.5% sodium dodecyl sulfate, and 1 mL lysozyme) and proteinase K 0.5 mL (stock 20 mg/mL) incubated at 37 • C and 55 • C for 30 min and 60 min, respectively, and allowed to cool at room temperature (RT) for 20 min.Potassium acetate (5M) 4 mL was added to precipitate the protein, and samples were centrifuged at 12,000× g at 4 • C for 10 min.DNA was purified using (1:1) phenol: chloroform: isoamyl alcohol (25:24:1).Isopropanol was used for DNA precipitation, and Tris-EDTA (TE) buffer was used to dissolve the DNA total volume 30 mL/membranes set (0.45 µm + 0.025 µm).The DNA was stored at −20 • C for future metagenomic analysis.Total extracted DNA was confirmed by agarose gel electrophoresis.

DNA Extraction from Whole Blood Samples
DNA was extracted from whole blood samples (300 µL) using a Genomic DNA Isolation Kit (PDC11-0100; BIO-HELIX Co., Xindian Dist, New Taipei City, Taiwan).Briefly, 900 µL of buffer CR was added and mixed by inversion.The tubes were incubated at RT for 10 min and then centrifuged for 5 min at 4000× g.Next, 300 µL of buffer CC was added to the resuspended cells and incubated at 60 • C for 10 min until the sample lysates were clear.Then, 400 µL of buffer CB was added and shaken vigorously, followed by centrifugation at 12,000× g for 1 min.Column CC was placed in a 2 mL collection tube.The clear supernatant from the previous step was transferred completely to the column and centrifuged at 14,000× g for 30 s.The flow-through was discarded, 400 µL of the buffer W1 was added to the column CC, and samples were centrifuged at 14,000× g for 30 s.The flow-through was discarded, and 600 µL of buffer W2 was added to column CC and centrifuged at 14,000× g for 30 s.The flow-through was discarded, and the column was centrifuged at 14,000× g for 2 min.Finally, buffer TE (50-200 µL) was added and centrifuged for 2 min at 14,000× g.The isolated DNA was then stored at −20 • C.

Metagenomics
The DNA extracted from the tap water, natural water and the pooled human whole blood samples were analyzed by shotgun metagenomic sequencing.The extracted DNA represented all viruses in the water and whole blood samples.Fraction steps were checked using Agilent 2100 and qPCR to confirm that sufficient target DNA was present in the end repairing range.For tail construction, sequencing adapters were ligated, and PCR enrichment steps were used to create the library.Sequencing was performed using the Illumina platform after library clustering with paired-end reads.

Geographic Information System (GIS)
A GIS was used to connect the data on oncogenic virus prevalence in the Mueang District and show their distribution.GIS connected oncogenic virus data to a map and integrated the location data of 15 subdistricts, including Wiang, Mae Tam, Mae Na Ruea, Ban Tun, Ban Tam, Ban Tom, Mae Puem, Mae Ka, Ban Mai, Cham Pa Wai, Tha Wang Thong, Mae Sai, Ban Sang, Tha Champi, and San Pa Muang.

Statistical Analysis
The data were analyzed using the IBM SPSS software version 16, and p < 0.05 was considered statistically significant.Pearson's Chi-squared test was used to compare categorical variables between groups.An independent Student's t-test was used to compare separate mean ± standard deviation (SD) sets.

Water Quality Analysis
Tap water sources from the Mae Ka (No. 1) and Wiang Subdistricts were collected from Mueang Phayao District, Phayao Province.Chemical and microbial analysis indicated that tap water from Mae Ka (No. 1) contained 0.330 mg/L Fe, 0.135 mg/L Mn, turbidity of 17.235 mg/L, and E. coli and total coliform bacteria numbers above the 2017 standard value from the Thai Provincial Waterworks Authority [65].Chemical and microbial analysis of tap water from the Wiang Subdistrict showed that it also contained total coliform bacteria numbers above the standard value [65].The results are shown in Table 1.
Chemical and microbial analysis of tap water from Mae Ka (No. 1) showed it contained higher than standard values, which are harmful to consumers.Its turbidity was visible to the naked eye, while tap water from the Wiang Subdistrict was clear.Tap water from both sources was contaminated with total coliform bacteria, which are harmful to the health of individuals in the community.

Metagenomics Analysis
Extracted DNA from tap water, natural water, and pooled whole blood samples were analyzed by shotgun metagenomic sequencing.Surprisingly, the metagenomic analysis indicated the presence of pathogenic bacteria such as E. coli, Salmonella enterica, Salmonella spp., Staphylococcus aureus, Legionella spp., Pseudomonas aeruginosa, Vibrio cholerae, Serratia marcescens, Klebsiella pneumoniae, Shigella sonnei, Shigella flexneri, and Burkholderia pseudomallei in all the water samples, including tap water.The metagenomics analysis also revealed Cryptosporidium and Giardia spp. in all water samples, including tap water.Additionally, unidentified (virus groups) comprised 0.2% of all microbes.
Herpesvirales was found in human blood and all water samples (Figure 1).Equid gammaherpesvirus 5 was found in all water samples.Herpesviridae was the predominant oncogenic DNA virus, such as oncogenic HHV-8, in both natural water sources.Herpesviridae was the predominant disease-causing DNA virus found in human blood and water samples.In pooled human blood samples, EBV comprised 2% of total Herpesviridae.In addition, EBV comprised 54% of all Herpesviridae in Mae Ka (No. 1) tap water, 89% in Wiang subdistrict tap water, 43% in the Huai Na Poi water, and 37% in the Kwan Phayao water.Tap water from Mae Ka (No. 1) and Wiang Subdistrict, natural water from Huai Na Poi and Kwan Phayao, and pooled human blood samples contained 6.42, 100.12, 6.36, 7.66 EBV copies/L, and 4.53 EBV copies/100 blood samples.This study shows that EBV is found in tap and natural water sources and human blood samples.Human betaherpesvirus 6B was found in the natural water from Kwan Phayao, Human betaherpesvirus 6A was found in Mae Ka (No. 1) tap water and human blood sample.
ples.In pooled human blood samples, EBV comprised 2% of total Herpesviridae.In addition, EBV comprised 54% of all Herpesviridae in Mae Ka (No. 1) tap water, 89% in Wiang subdistrict tap water, 43% in the Huai Na Poi water, and 37% in the Kwan Phayao water.Tap water from Mae Ka (No. 1) and Wiang Subdistrict, natural water from Huai Na Poi and Kwan Phayao, and pooled human blood samples contained 6.42, 100.12, 6.36, 7.66 EBV copies/L, and 4.53 EBV copies/100 blood samples.This study shows that EBV is found in tap and natural water sources and human blood samples.Human betaherpesvirus 6B was found in the natural water from Kwan Phayao, Human betaherpesvirus 6A was found in Mae Ka (No. 1) tap water and human blood sample.The sequences of EBV from natural and tap water were similar to EBV reference strain (RefSeq) NC_007605 (171,823 bp DNA circular, Human gammaherpesvirus 4, complete genome, version NC_007605.1),mostly at location about 800-102 kb.The sequence of EBV from human blood was similar to EBV reference strain (RefSeq) NC_007605 (171,823 bp DNA circular, Human gammaherpesvirus 4, complete genome, version NC_007605.1)at location about 800-102 kb but not the same location with water.Therefore, in all samples (human blood, tap water, and natural water) were found Human gammaherpesvirus 4 (such as at the location of EBNA-1 gene) with mostly percentage identity up to 100% when using cutoff >50 bp.However, when using a <50 bp cutoff, the result found that human blood samples and all types of water were similar with EBV reference strain (RefSeq) NC_007605 at the same location with percentage identity mostly up to 100% (Figure 2, Supplementary Figures S1 and S2).The sequences of EBV from natural and tap water were similar to EBV reference strain (RefSeq) NC_007605 (171,823 bp DNA circular, Human gammaherpesvirus 4, complete genome, version NC_007605.1),mostly at location about 800-102 kb.The sequence of EBV from human blood was similar to EBV reference strain (RefSeq) NC_007605 (171,823 bp DNA circular, Human gammaherpesvirus 4, complete genome, version NC_007605.1)at location about 800-102 kb but not the same location with water.Therefore, in all samples (human blood, tap water, and natural water) were found Human gammaherpesvirus 4 (such as at the location of EBNA-1 gene) with mostly percentage identity up to 100% when using cutoff >50 bp.However, when using a <50 bp cutoff, the result found that human blood samples and all types of water were similar with EBV reference strain (RefSeq) NC_007605 at the same location with percentage identity mostly up to 100% (Figure 2, Supplementary Figures S1 and S2).

EBV DNA Detection by PCR and qPCR and Environment Factors
Human blood samples were collected from 813 blood donors in Phayao Province, comprising 231 men (28.41%) and 582 women (71.59%) aged 3-90 years (mean 52.01 ± 18.65 years).PCR identified 7/813 (0.86%) EBV-positive samples (EBNA-2 gene), while qPCR identified 6/813 (0.74%) and 53/813 (6.52%) samples based on EBNA-2 and EBNA-1, respectively.EBV prevalence based on all three methods was 59/813 (7.26%;Table 2).EBV prevalence differed significantly by method (p > 0.05).EBNA-1 gene detection by qPCR showed higher sensitivity in blood than EBNA-2 gene detection.EBNA-2 gene detection was also inconsistent between PCR and qPCR.Age, congenital disease, family cancer history, alcohol consumption, secondhand smoke, drinking water source (bottled water), and tap water used for brushing teeth and eating fresh fruit were associated with EBV DNA positivity in blood, depending on the detection method used.BMI, exercise, smoking, drinking water source (boiled or filtered), cleaning of water for consumption, and water used with vegetables were not significantly associated with EBV DNA detection, depending on the detection method used.The results are shown in Table 2 and Figure 3.
eating fresh fruit were associated with EBV DNA positivity in blood, depending on the detection method used.BMI, exercise, smoking, drinking water source (boiled or filtered), cleaning of water for consumption, and water used with vegetables were not significantly associated with EBV DNA detection, depending on the detection method used.The results are shown in Table 2 and Figure 3. Life history (yes/no) and risk factors for EBV infection (%) in blood based on the three methods.The color represents the percentage of exposure to environmental risk factors: 0%-5%, blue; 5%-10%, green; 10%-15%, yellow; and 15%-20%, red.Male, congenital disease, family cancer history, exercise, and boiled/filtered drinking water were not risk factors for EBV DNA detection in blood.Not using tap water for brushing teeth, alcohol consumption, secondhand smoke exposure, and not eating fresh fruit weekly were risk factor factors for EBV DNA detection in blood.Note: not using tap water for brushing teeth reflects the use of other water sources, such as rain/artesian water.Key: *, statistically significant (p < 0.05).Individuals with congenital disease had higher PCR-based EBNA-2 gene positivity (6/386 (1.55%)) than those with no congenital disease (1/427 (0.23%); p = 0.042).Individuals with a family cancer history had higher qPCR-based EBNA-2 gene positivity (4/184 (2.17%)) than those with no family cancer history (2/629 (0.32%); p = 0.010; odds ratio = 6.967; 95% CI: 1.266-38.344).Individuals who did not use bottled water had a higher qPCR-based EBNA-1 gene positivity (8/63 (12.69%)) than those who used bottled drinking water (45/750 (6.00%); p = 0.039; odds ratio = 2.279; 95% CI: 1.023-5.074).

GIS
To assess the epidemic spread of EBV in various sub-districts in Phayao province, Thailand, GIS was used to connect and integrate EBV prevalence data from human blood samples across 15 subdistricts in the Mueang District and show their distribution as a map (Figure 4).

GIS
To assess the epidemic spread of EBV in various sub-districts in Phayao province, Thailand, GIS was used to connect and integrate EBV prevalence data from human blood samples across 15 subdistricts in the Mueang District and show their distribution as a map (Figure 4).

Discussion
The 2017 report by the Thai Provincial Waterworks Authority [65] indicated that total coliform bacteria, S. aureus, Salmonella spp., Clostridium perfringens, E. coli, and Shigella spp.were not detected.However, this study detected all these pathogens in both tap water sources.Therefore, tap water quality is not passing the standards of the Provincial Waterworks Authority, affecting the health of the residents in the district.
The Thai Metropolitan Waterworks Authority has increased the testing requirements for viruses to determine the safety of tap water [66].This study did not address RNA viruses.However, it detected pathogens in both tap water sources that are prohibited by these standards, including total coliform bacteria, S. aureus, Salmonella spp., E. coli, Shigella spp., Legionella spp., P. aeruginosa, V. cholerae, C. perfringens, Cryptosporidium spp., and Giardia spp.Producing tap water from dirty water sources can cause many problems, and its safety should be taken seriously.Even though the tap water sources from Phayao are unsuitable for drinking, most of the population uses tap water to brush their teeth, wash vegetables and fruit, or wash dishes, which may be a channel for bacterial and viral diseases.In particular, the oncogenic viruses found in this study (EBV) were present in large quantities in water sources.
This study found that the sequence of oncogenic EBV from natural water, tap water, and human blood sample was similar to the EBV reference strain, Human gammaherpesvirus 4, complete genome (version NC_007605.1).Possibly, there are associations between oncogenic EBV in water and humans and further study should be conducted on the risk of oncogenic EBV in water being transmitted to humans.
Herpesviruses (equine herpesvirus type-1) have been found to persist in various water environments for up to three weeks in controlled experiments [67].This study found equid gammaherpesvirus 5 in all water sources, which causes chronic lung disease and equine multinodular pulmonary fibrosis in adult horses [68].The herpes simplex virus (HSV) survived 4 h in tap water, 24 h in distilled water, and 4.5 h on plastic surfaces, suggesting that fomites such as these may be nonvenereal routes for HSV transmission in spa water [69].Therefore, oncogenic EBV in water suggests that transmission to humans is possible and should be investigated further.
This study is the first to discuss the safety of natural water, tap water and community health in Thailand in the EBV context.Currently, non-communicable diseases account for 71% of deaths in Thailand, followed by cancer (17%) [70].Many factors, including viruses, cause these diseases.Water is frequently a vehicle for viral transmission.Virus survival is higher in sterile water [11][12][13].This study provides an overview of viral communities that exist in natural water sources in the large Phayao District using metagenomics.Therefore, it can serve as a future guideline for public health agencies and the Provincial Waterworks Authority.
This study found that individuals aged 11-50 were most likely to be infected due to the increased frequency of close social contact with infected individuals.Additionally, EBV infections might play a role in individual immune responses, which are more severe in adolescence or adulthood [38].EBV prevalence differed based on the detection method used.Therefore, several methods or genes are required for EBV detection.This study found the low prevalence of blood-borne EBV infection.Thailand is in a tropical zone and has year-round sunshine.Sunlight and vitamin D are believed to protect against autoimmunity by increasing the number of CD8+ T cells available to control low-level EBV infections [71].
EBV transits between epithelial and B cells.The first innate barrier to infection is the oropharyngeal epithelium.EBV replicates mainly in epithelial cells in the oropharynx.Therefore, EBV is believed to undergo replication to be efficiently transmitted and spread to other hosts via saliva [72].Contact through saliva, such as kissing, sharing glasses, and not using a serving spoon or utensils, may cause EBV infection in epithelial cells but not in blood.Therefore, this study did not find an association between these factors and blood-borne EBV infection.EBV spreading from epithelial cells to the blood might be associated with other factors, such as host immune responses.EBV DNA in whole blood is a superior prognostic and monitoring factor in many diseases, such as diffuse large B-cell lymphoma [46].Therefore, the presence of blood-borne EBV may be clinically significant [38] and for blood transfusion.
This study found fewer blood-borne EBV infections in individuals who used tap water for brushing their teeth compared to rain or artesian water, which is still used in rural areas.The Mueang District is the center of Phayao Province, which is a small province, and the authority does not provide tap water.In some areas, villages still produce tap water, and rain and artesian water are still used.Therefore, water quality might not be appropriately controlled, and contact with infectious agents via water may be a concern based on the water source used.
The EBNA-1 gene is important for cell immortalization [73], and the EBNA-2 gene is essential for establishing a latent EBV infection and B-cell immortalization [74,75].Individuals not using bottled water had a higher EBNA-1 gene positivity than those who used bottled water.While bottled water is generally considered safe, its sterilization process remains of concern [76,77].
This study found that environmental factors such as age, alcohol consumption, secondhand smoke exposure, and not eating fresh fruit were risk factors for blood-borne EBV infection.Smoking was associated with EBV infection [78] and increased risk of EBV-positive Hodgkin lymphoma (odds ratio = 1.4; 95% CI: 1.1-1.9)[79], gastric cancer [80], and NPC [81].However, no mention of secondhand smoke exposure was found in previous studies.In addition, they found no significant association between EBV and alcohol consumption in NPC [80,81] or gastric cancer [80].However, this study found that blood-borne EBV infection was associated with alcohol consumption.Nevertheless, several studies have found an association between habitual alcohol consumption and NPC and OSCC risk [82][83][84][85].
This study found that eating fresh fruit may be protective against NPC caused by blood-borne EBV infection [86,87].This study also found that BMI was not associated with blood-borne EBV infection, different from Keegan et al. [88].While obesity has been shown to impact the host immune response to infections that can induce inflammation, its association remains unclear [38].A family cancer history has been previously associated with NPC [89], according to this study.

Conclusions
The findings of this study support the conclusion that water in developing countries continues to be unsafe.Human oncogenic EBV was present in natural and tap water.Finally, basic knowledge about tap water sterilization processes and viral control should be promoted.Blood-borne oncogenic EBV should be a concern for blood transfusion in patients.Therefore, the detection of oncogenic EBV in water suggests that transmission via water is possible and should be investigated further.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/10.3390/w15020323/s1,Supplementary File S1: metagenomic analysis using a cutoff of more than 50 bp and Supplementary File S2: metagenomic analysis using a cutoff of less than 50 bp.Supplementary Figures S1 (The sequence location of EBV from blood and water on EBV) and S2 (The alignment sequence of EBV from water compare with EBV) in Supplementary File S3.

Figure 1 .
Figure 1.Herpesvirales in human blood and all water samples.Percentage (%) values represent the percentage of the virus among all viruses.Key: HB, pooled human blood samples; MTW, tap water from Mae Ka (No. 1); PTW, tap water from Wiang Subdistrict; KNW, natural water from Kwan Phayao; NNW, natural water from Huai Na Poi.

Figure 1 .
Figure 1.Herpesvirales in human blood and all water samples.Percentage (%) values represent the percentage of the virus among all viruses.Key: HB, pooled human blood samples; MTW, tap water from Mae Ka (No. 1); PTW, tap water from Wiang Subdistrict; KNW, natural water from Kwan Phayao; NNW, natural water from Huai Na Poi.

Figure 2 .
Figure 2. The examples sequence location of EBV from blood and water on EBV RefSeq NC_007605 (171,823 bp complete genome) when using cutoff >50 bp.EBV were mostly found at the location about 800kb-102kb.Key: HB, pooled human blood samples; MTW, tap water from Mae Ka (No. 1); PTW, tap water from Wiang Subdistrict; KNW, natural water from Kwan Phayao; NNW, natural water from Huai Na Poi.

Figure 2 .
Figure 2. The examples sequence location of EBV from blood and water on EBV RefSeq NC_007605 (171,823 bp complete genome) when using cutoff >50 bp.EBV were mostly found at the location about 800 kb-102 kb.Key: HB, pooled human blood samples; MTW, tap water from Mae Ka (No. 1); PTW, tap water from Wiang Subdistrict; KNW, natural water from Kwan Phayao; NNW, natural water from Huai Na Poi.

Figure 3 .
Figure3.Life history (yes/no) and risk factors for EBV infection (%) in blood based on the three methods.The color represents the percentage of exposure to environmental risk factors: 0%-5%, blue; 5%-10%, green; 10%-15%, yellow; and 15%-20%, red.Male, congenital disease, family cancer history, exercise, and boiled/filtered drinking water were not risk factors for EBV DNA detection in blood.Not using tap water for brushing teeth, alcohol consumption, secondhand smoke exposure, and not eating fresh fruit weekly were risk factor factors for EBV DNA detection in blood.Note: not using tap water for brushing teeth reflects the use of other water sources, such as rain/artesian water.Key: *, statistically significant (p < 0.05).

Figure 3 .
Figure3.Life history (yes/no) and risk factors for EBV infection (%) in blood based on the three methods.The color represents the percentage of exposure to environmental risk factors: 0-5%, blue; 5-10%, green; 10-15%, yellow; and 15-20%, red.Male, congenital disease, family cancer history, exercise, and boiled/filtered drinking water were not risk factors for EBV DNA detection in blood.Not using tap water for brushing teeth, alcohol consumption, secondhand smoke exposure, and not eating fresh fruit weekly were risk factor factors for EBV DNA detection in blood.Note: not using tap water for brushing teeth reflects the use of other water sources, such as rain/artesian water.Key: *, statistically significant (p < 0.05).

Figure 4 .
Figure 4. GIS-connected oncogenic EBV prevalence data from human blood samples in the 15 Mueang subdistricts highlights the EBV distribution.EBV-positive means EBV-positive in human blood.

Table 1 .
The water turbidity, color, pH, and trace element analysis.