Urinary Virome Perturbations in Kidney Transplantation
Mar 18, 2022
Contact: Audrey Hu Whatsapp/hp: 0086 13880143964 Email: audrey.hu@wecistanche.com
Tara K. Sigdel et al
The human microbiome is important for health and plays a role in essential metabolic functions and protection from certain pathogens. Conversely, dysbiosis of the microbiome is seen in the context of various diseases. Recent studies have highlighted that a complex microbial community containing hundreds of bacteria colonizes the healthy urinary tract, but little is known about the human urinary viruses in health and disease. To evaluate the human urinary virome in the context of kidney transplantation (tx), variations in the composition of the urinary virome were evaluated in urine samples from normal healthy volunteers as well as patients with kidney disease after they had undergone kidney tx. Liquid chromatography-mass spectrometry/mass spectrometry analysis was undertaken on a selected cohort of 142 kidneys tx patients and normal healthy controls, from a larger biobank of 770 kidney biopsy matched urine samples. In addition to analysis of normal healthy control urine, the cohort of kidney tx patients had biopsy-confirmed phenotype classification, coincident with the urine sample analyzed, of stable grafts (STA), acute rejection, BK virus nephritis, and chronic allograft nephropathy. We identified 37 unique viruses, 29 of which are being identified for the first time in human urine samples. The composition of the human urinary virome differs in health and kidney injury, and the distribution of viral proteins in the urinary tract may be further impacted by IS exposure, diet, and environmental, dietary, or cutaneous exposure to various insecticides and pesticides.
Keywords: kidney transplantation, virome, urine, proteomics, biomarkers
Cistanche tubulosa prevents kidney disease, click here to get the sample
INTRODUCTION
The human microbiome has been studied extensively in various biofluids, such as blood (1), urine (2–4), saliva (5, 6), cerebrospinal fluid (7), and bronchoalveolar lavage (8, 9), for its influence on human health and disease (10). Several reports on the assessment of microbiomes in different regions of the human body have been reported including lung and gut (11–15). However, only little is known about the urinary microbiome and its changes in the context of kidney injury after kidney transplantation (tx) (16, 17). Most human microbiome studies have mapped 16S rRNA for bacterial profiling (18), and there is only a small body of data examining the urine virome (19, 20), and a single published study has used next-generation sequencing to map viral genomic components in the urine of kidney transplant patients (19).
Dysbiosis of the microbiome is associated with multiple diseases such as inflammatory bowel disease, colon cancer, obesity, and pulmonary disease (21–25), but little is known about alterations of the virome in human body fluids, such as saliva, bronchoalveolar lavage, and urine. It is recognized that various components of the microbiome perform essential functions including biosynthesis of cofactors and vitamins, metabolism of essential compounds, and barrier protection from pathogens (26–28). Similar dysbiotic and protective functions may also be ascribed to the urine virome; hence, it is important to ascertain the composition of the human virome in health and disease. For purposes of this study, we chose to conduct these studies in urine in health and in the context of IS exposure after kidney tx, during stable kidney function, and in the context of kidney tx injury. The focus of the study was to evaluate the repertoire of viruses as determined by the peptides present in human urine and their perturbations after IS exposure and various causes of acute, chronic, and infectious kidney injury.
The microbiome changes over time and correlates with organism diversity. Microbiome analysis can be segregated into both structural analysis based on operational taxonomic units based on sequence phylogeny, and functional analysis based on metagenomics sequencing and proteomics (29, 30). Identified microbiota in different diseases can be further studied by cultivation, functional metagenomics, and multiplexed immunofluorescence and in situ hybridization. The NIH human microbiome project has published the human microbiome in 15 body sites from 300 individuals (31).
TREATMENT OF KIDNEY DISEASE: CISTANCHE
MATERIALS AND METHODS
A total of 142 unique samples were evaluated from a biorepository containing 2016 collected by IRB approved informed consent from adult and pediatric samples from the kidney tx programs at Stanford University and University of California San Francisco, between urine samples of which 770 were accompanied with matched kidney tx box with centralized pathology histology reads and compartment scores using the standardized Banff schema (32) for scoring kidney tx by injury (Figure 1). The study was approved by The Human Research Protection Program of the University of California, San Francisco. The urine samples were phenotyped based on the matched kidney bx pathology into five groups: healthy control (HC; n = 9), Stable graft (STA; n = 40), acute rejection (AR; n = 37), chronic allograft nephropathy (CAN; n = 39), and BK virus nephritis (BKVN; n = 17). Urine was centrifuged 2,000 ×g at 4°C for 20 min to get rid of urine sediments. The supernatant was passed through a filter membrane of 10 kDa to remove native peptides from intact proteins larger than 10 kDa in size. The total protein was then trypsin digested and the resulting tryptic peptides were analyzed by the LC-MS platform (Orbitrap Velos MS). The detailed methods of protein preparation and analyses are reported elsewhere (33).
The MSGF plus customized algorithm generated by our group (https://omics.pnl.gov/software/ms-gf), was used to search MS/ MS spectra against the combined human protein sequence database and the NCBI viral database. Peptides were initially identified from database searching applying the following criteria: MSGF spectrum E-value (a probability value of the peptide to MS/MS spectrum match with the lower value the higher probability to be a correct match) to be <10-10, Peptide level Q-value (false discovery rate estimated by targeted-decoy database search) to be <0.01, and mass measurement error <10 ppm (±5 ppm). The decoy database searching methodology was used to confirm the final false discovery rate at the unique peptide level to be <1%. Due to the anticipated higher false discovery rate for peptides from viral proteins, more stringent filtering criteria with MSGF spectrum E value to be <1E-13 was applied. The false discovery rate was estimated to be nearly 0% based on the well-accepted target decoy searching strategy because no decoy hits were observed following this stringent cutoff. Data are shown as percentages and mean ± SD. Comparisons of different categories are done using ANOVA and p values of <0.05 are considered significant.
RESULTS
Our group has previously published a detailed analysis of biologically relevant human proteins in these urine samples collected from kidney transplant recipients with different graft injury phenotypes, as confirmed by matched kidney transplant histopathology on the biopsy, collected at the same time as the urine sample; this data has been deposited in the proteomic MassIVE repository (accession MSV000079262) and in the ProteomeXchange repository (accession PXD002761) (33). In this study, we only focused on the identification and analysis of viral proteins in the same cohort of kidney transplant patients, with the inclusion also of age- and gender-matched healthy control human urine samples to evaluate viral proteins in both health and kidney injury. It is important to note that of the total analyzed human and viral proteins in urine, viral proteins alone constitute <0.2% of the total identified proteins, highlighting the very low abundance of rather rare viral proteins in human urine, irrespective of the type of kidney injury, and irrespective of baseline immunosuppression usage.
The results presented in this paper come from viral peptide mapping, unlike genomic sequencing of viral components as reported by previous publications (19). As an initial analysis, we focus on the prevalence of urinary viral proteins specific to each kidney transplant phenotype of STA, AR, CAN, and BKVN, and the variations noted over the healthy control urine virome. Urinary viral protein data for each sample was evaluated for sample abundance relative to the mean level of that virus in the entire sampled population. We found that on an average, 22% (range 4–67%) of kidney transplant patients had different viral proteins detected in their urine, excluding patients with BKVN, where detection of the BK virus in urine is to be expected, as these patients have an infection with BK virus in their urinary tract and in the kidney transplant. The distribution of different viral proteins in different tx injury states and normal health is shown in the prevalence heat map (Figure 2). Overall, a total of 57 viral proteins from 37 unique viruses were identified, many of these being unexpected and previously not described as “commensals” in human urine. Of these unique viruses, 8 viruses were single-stranded RNA viruses, 1 virus was a single-stranded RNA-RT virus, and 28 viruses were double-stranded DNA viruses (Table 1). Out of the 142 patients, all patients had at least one viral protein in their urine, with an average of 10.23 ± 4.57 viral proteins/sample.
Prevalence of Human Urine Viral Proteins
The healthy control group examined had 20 unique viral proteins with a group average of 13.67 ± 1.32 viral proteins. Kidney tx patients with maintenance IS and stable kidney tx function had a statistically significantly lower number of unique viruses in the group (8.08 ± 2.78; p = 7e-10) when compared to the healthy cohort (normal kidney function and a normal immune system), without any IS exposure. Nine new viral proteins are instead detected in the STA cohort compared to healthy controls, belonging to the following viruses: Junin virus, Cotesia congregata bracovirus, Agrotis segetum nucleopolyhedrovirus, Psittacid herpesvirus 1, Pseudocowpow virus, Spodoptera frugiperda multiple nucleopolyhedroviruses, Japanese yam mosaic virus, Cowpea mottle virus, and Helicoverpa zea single nucleopolyhedrovirus. Kidney injury after tx (despite continued IS exposure), results in an overall increase in the number of unique viral proteins over the STA tx cohort, with a significant increase in urine viral proteins in the CAN patient group (11.46 ± 3.52; p = 8e-10), and the BKVN patient group (15.76 ± 5.65; p = 13.7e-10). The repertoire of urine viral proteins appears to be quite distinct in different tx injury categories (Figure 3). The prevalence of BKV viral proteins in urine increases to 60–70% in AR, 70–80% in patients in CAN, and 100% in patients with BKVN (Figure 2) highlighting that increasing abundance of the BK virus in the urine may result from augmentation of IS, as seen in the AR category, and with greater time post-tx, as seen in the CAN category. BKVN urine samples show the maximal divergence of urinary viral proteins, as expected. Four viral proteins that are consistently present in all other samples, inclusive of normal healthy controls (Canarypox virus, Tacaribe virus, Simian virus 12, Acidianus filamentous virus 8), are no longer observed in the BKVN cohort. As BKV is a DNA virus, we examined if the changes in the urinary virome are largely related to the emergence of new DNA viruses in BKVN. Five new viral proteins are noted only in the BKVN cohort (Ectropis obliqua nucleopolyhedrovirus, Fowl adenovirus D, Taura syndrome virus, Invertebrate iridescent virus 6, Gallid herpesvirus 2). 4/5 of BKVN cohort viruses were double-stranded DNA viruses (Table 1).
The abundance of Human Urine Viral Proteins
We next evaluated the abundance of different viral proteins across all tx patients and HC urine samples. Among the 37 unique viruses examined, the following viruses were the most abundant in STA IS patients (Acidianus filamentous virus 8, Agrotis segetum nucleopolyhedrovirus, BK polyomavirus, Canarypox virus, Cotesia congregata bracovirus, Cowpea mottle virus, Glossina pallidipes salivary gland, hypertrophy virus, Simian virus 12, Tacaribe virus), and increased abundance of the Marsellavirus is seen in both CAN and BKVN. Despite intrarenal infection with the BKV virus in BKVN patients, we observe that low levels of the BK virus are detected also in all examined tx categories and healthy controls, likely as the BK virus is normal human urine commensal, with its persistence in the human uroepithelium at low levels in almost all examined HC and kidney tx patients is not unexpected.
TREATMENT OF KIDNEY DISEASE: CISTANCHE
DISCUSSION
A study the role of the microbiome in organ transplantation (13, 17, 34) takes into account the impact of food intake, digestion, metabolism, and modulation; however, dysbiosis of microbiota due to the tx and IS medicines is a contributing factor that decreases in the baseline predominant microbes and also results in a loss of overall diversity, and the emergence of few new dominant microbial populations (35). A very interesting finding from this study is that only 8 out of the 37 viruses identified in this dataset have been previously described in humans, and many of these have also been described to be pathogenic. Human infections have been described with some of the identified viruses, suggesting that their presence in urine may be related to the pathogenesis of general systemic or underlying renal injury. Junin virus (an arenavirus) infection can result in clinical human disease inclusive of fever, as well as an entity known as Argentine hemorrhagic fever (36). Tacaribe virus is an arenavirus that can also cause human fevers and hemorrhagic disease (36). The HHV-6 virus is a highly prevalent virus in children and causes fever, diarrhea, and rashes (37). Rhinoviruses are single-stranded RNA viruses that are the most infectious agent in humans and the most likely culprit for the common cold (38). Acanthocystis surface chlorella virus 1 has been in 44% of healthy humans (39) and has only been described in the oropharynx and is not known to be pathogenic to humans. The Marseille virus is a nucleocytoplasmic large DNA virus that has been described in the blood and stool of patients with nonfebrile lymphadenopathy (40).
The SV40 virus is a DNA polyomavirus that was likely introduced into the human population through contaminated vaccines (41) and infection with this virus is frequently co-infected with BK virus infection, this being an almost invariant finding in the tx population. Recent studies have also highlighted the role of chronic SV40 infection with human cancers (42). Antibodies against the T antigen of SV40 cross-reacts with the T antigens of BK and JC viruses, which are all in the polyomavirus family, and the SV40 stain is used for diagnosis of BK virus infection in BKVN. BK virus exposure, in the kidney, is seen in 90% of the normal population, and our data suggest that the urinary virome has traces of BK viral proteins in all the healthy controls sampled.
BK virus replication is seen in 10–60% of kidney tx patients who have been described to shed BK virus in their urine, which is confirmed by our data (Figure 2). BK virus is latent in renal tubular epithelial cells and its prevalence is known to increase immune dysfunction and immunosuppression exposure (43). In this study, we could not identify any urinary virome-specific differences that could distinguish the AR and BKVN groups, despite these conditions being immunologically distinct and requiring very different treatment approaches, i.e., minimization of immunosuppression in BKVN and increased immunosuppression in AR. The AR category was be confirmed and there was a negative SV40 stain in the box tissue examined. Nevertheless, this data suggests that the detection of BKVN may be under-reported in kidney tx patients, and histological changes may be patchy resulting in BKVN disease underdiagnosis.
Of the remaining 19 viruses that have not been previously described in humans, the following have been noted to be used as pesticides/insecticides (E. obliqua nucleopolyhedrovirus, Euproctis pseudoconspersa nucleopolyhedrovirus, Helicoverpa armigera nucleopolyhedrovirus, H. zea single nucleopolyhedrovirus). Pesticides in public health use are intended to limit the potential for the disease but they have been known to make their way into human tissues/bodily fluids if not properly handled and that organic diets limit the number of pesticides in children's urine (44, 45). The finding of these viral proteins in human urine samples raises a question of environmental exposure to these previously undiscovered, possible contaminants.
The finding of 100% prevalence of some viruses only in the healthy control group suggests that there are viral commensals that may exist either without harm or may even be hypothesized to assist in immune adaptation. Phylogenetic tree analysis of these viruses may be interesting, as some viruses such as bacteriophages play an important role in health by removing pathogenic bacteria and helping boost innate immunity (16). A reduction in the protective viruses or phages may result in increased bacterial proliferation which is seen in the context of IS after kidney tx. Metagenomic sequencing of the new viral species found in tx patients under IS may identify if IS exposure has changed the resistance patterns or resistome of some of the viruses, especially those that are more prone to drive inflammation and immunity (16). Silver nanoparticle treatment has been shown to reduce gut viral and bacterial populations that are pro-inflammatory (46). Mucinophilic microbiomes may be more prone to provide non-host protective immunity, and the loss of any eosinophilic viruses/microbes with IS may also be a driver of the increased risk of cystitis and urinary tract infections seen in patients after tx.
The microbiome intersects with states of health and disease. The novel detection of a large number of viruses used as pesticides/insecticides is a surprising finding. Identification of viruses that are not expected to be present in the urine points to a possibility of their origin to enter bodily systems through ingestion (contamination of food, water, or other drinks such as milk) or cutaneous absorption (e.g., improper handwashing). In fact, according to federal law, a small residue of pesticides/insecticide contamination of human food is recognized and accepted (41, 47), even though no one to date has examined their presence in urine. It is also possible that these viruses colonize or invade the urinary tract by ascending infection, especially in the context of IS in a kidney tx patient, where there is lowered host immune defense.
Further validation of these results by evaluation of the human virome repertoire in different geographic and demographic cohorts of patients with different causes of kidney injury will better identify if the particular repertoire of viral proteins observed in this study is specific to a patient, their geography, their demography, or the kidney disease subtype. Additional validation studies by urine viral PCR assays can provide for rapid assessment of the clinical impact of these viruses in other kidney injury cohorts. Interestingly, in the repertoire of identified viral peptides, we did not observe any phages. This is attributed to the methodology used for peptide preparation from urine supernatant, which does not capture bacteria (which get pelleted out), the host of phage viruses (48). The virome report based on DNA sequence analysis also did not identify phages in the published report (19). The urine virome study that used purified bacteriophage isolation using cesium chloride density gradient has reported the identification of phages (20). These articles all support that the methodology of sample prep influences the eventual findings.
Colonization, resistance, and microbial ecology has been described and well studied in the context of microbial infections (49), and it has been shown that antibiotic treatment (in the context of treating infections in the animals that form food products or as part of the antiviral prophylaxis for tx patients in the first 3–6 months post-tx) eliminates many commensal bacterial and viral species from the gut and other body fluids and reduces antimicrobial defense (50). Colonization by microbiota in early life shapes the immune system and creates a window of opportunity for a homeostatic state, which if disrupted, in the context of tx and IS, can contribute to inflammation. Precision medicine approaches in the future may want to customize therapy from a “microbiota centered” and “host centered” approach to precision medicine. Until recently, we have not paid much attention to the food microbiome and the understanding that many emerging infections may be vectored by foods. Current efforts to develop a metagenome tracker tool that will evaluate the impact of different infections—bacterial, viral, and fungal and their fingerprints in the context of the food that we eat, that may be influenced by ambient temperature, animal cohort poor survival causation, and agricultural practices inclusive of soil treatment (methyl bromide treatment of soil as a pesticide result in a huge increase in Bacillus species), water sources [different microbial loads in pond vs. well water (51)], genetics of the food plant cultivar that may result in altered plant defense to microbes, plants, and pesticides (51). Viral and bacterial lipoproteins may be responsible for immune-modulatory properties, and recent studies also suggest a role for the interaction of innate immunity for regulating microbiome diversity and controlling infection, by specific cell lineages such as mucosal invariant T cells (MAIT cells). Hence, additional studies that focus on the causation of microbial diversity of human tissues and its impact on inflammation and immune responses will be paramount to conduct as part of future research.
TREATMENT OF KIDNEY DISEASE: CISTANCHE
ETHICS STATEMENT
All the study samples were collected from pediatric and young adult recipients transplanted between the years 2000 and 2011 at Lucile Packard Children’s Hospital of Stanford University. The study was approved by the ethics committees of Stanford University Medical School and UCSF Medical Center. All adult patients and parents/guardians of non-adult patients provided written informed consent to participate in the research, in full adherence to the Declaration of Helsinki.
AUTHOR CONTRIBUTIONS
MS and TS participated in the design of the study; SN, TS, NM, and MS participated in the writing of the article and the interpretation of the results; NM performed the statistical analyses and evaluation of results; CN, KB-J, and W-JQ generated and processed the raw data.
FUNDING
We acknowledge support from the University of California San Francisco. The authors acknowledge the funding support from NIDDK R01DK083447 (to MS), DP3 DK110844 (W-JQ), and P41 GM103493 (W-JQ).
TREATMENT OF KIDNEY DISEASE: CISTANCHE
REFERENCE
1. Potgieter M, Bester J, Kell DB, Pretorius E. The dormant blood microbiome in chronic, inflammatory diseases. FEMS Microbiol Rev (2015) 39(4):567–91. doi:10.1093/female/fuv013
2. Zhang Y, Zhao F, Deng Y, Zhao Y, Ren H. Metagenomic and metabolomic analysis of the toxic effects of the tri-chloroacetamide-induced gut microbiome and urine metabolome perturbations in mice. J Proteome Res (2015) 14(4):1752–61. doi:10.1021/pr5011263
3. Lewis DA, Brown R, Williams J, White P, Jacobson SK, Marchesi JR, et al. The human urinary microbiome; bacterial DNA in voided urine of asymptomatic adults. Front Cell Infect Microbiol (2013) 3:41. doi:10.3389/fcimb.2013.00041
4. Fouts DE, Pieper R, Szpakowski S, Pohl H, Knoblach S, Suh MJ, et al. Integrated next-generation sequencing of 16S rDNA and metaproteomics differentiate the healthy urine microbiome from asymptomatic bacteriuria in neuropathic bladder associated with spinal cord injury. J Transl Med (2012) 10:174. doi:10.1186/1479-5876-10-174
5. Lim Y, Totsika M, Morrison M, Punyadeera C. The saliva microbiome profiles are minimally affected by collection methods or DNA extraction protocols. Sci Rep (2017) 7(1):8523. doi:10.1038/s41598-017-07885-3
6. Nasidze I, Quinque D, Li J, Li M, Tang K, Stoneking M. Comparative analysis of human saliva microbiome diversity by barcoded pyrosequencing and cloning approaches. Anal Biochem (2009) 391(1):64–8. doi:10.1016/j.ab.2009.04.034
7. Wang Y, Kasper LH. The role of the microbiome in central nervous system disorders. Brain Behav Immun (2014) 38:1–12. doi:10.1016/j.bbi.2013.12.015
8. Lee SH, Sung JY, Yong D, Chun J, Kim SY, Song JH, et al. Characterization of the microbiome in bronchoalveolar lavage fluid of patients with lung cancer compared with benign mass-like lesions. Lung Cancer (2016) 102:89–95. doi:10.1016/j.lungcan.2016.10.016
9. Twigg HL III, Morris A, Ghedin E, Curtis JL, Huffnagle GB, Crothers K, et al. Use of bronchoalveolar lavage to assess the respiratory microbiome: signal in the noise. Lancet Respir Med (2013) 1(5):354–6. doi:10.1016/S2213-2600 (13)70117-6
10. Human Microbiome Project Consortium. A framework for human microbiome research. Nature (2012) 486(7402):215–21. doi:10.1038/nature11209
11. Kroemer A, Elsabbagh AM, Matsumoto CS, Zasloff M, Fishbein TM. The microbiome and its implications in intestinal transplantation. Curr Opin Organ Transplant (2016) 21(2):135–9. doi:10.1097/MOT.0000000000000278
12. Weber D, Oefner PJ, Hiergeist A, Koestler J, Gessner A, Weber M, et al. Low urinary indoxyl sulfate levels early after transplantation reflect a disrupted microbiome and are associated with poor outcome. Blood (2015) 126(14):1723–8. doi:10.1182/blood-2015-04-638858
13. Vindigni SM, Surawicz CM. The gut microbiome: a clinically significant player in transplantation? Expert Rev Clin Immunol (2015) 11(7):781–3. doi:10.158 6/1744666X.2015.1043894
14. Becker J, Poroyko V, Bhorade S. The lung microbiome after lung transplantation. Expert Rev Respir Med (2014) 8(2):221–31. doi:10.1586/17476348.20 14.890518
15. Hartman AL, Lough DM, Barupal DK, Fiehn O, Fishbein T, Zasloff M, et al. Human gut microbiome adopts an alternative state following small bowel transplantation. Proc Natl Acad Sci U S A (2009) 106(40):17187–92. doi:10.1073/pnas.0904847106
16. Rani A, Ranjan R, McGee HS, Andropolis KE, Panchal DV, Hajjiri Z, et al.Urinary microbiome of kidney transplant patients reveals dysbiosis with potential for antibiotic resistance. Transl Res (2017) 181:59–70. doi:10.1016/j. trsl.2016.08.008
17. Fricke WF, Maddox C, Song Y, Bromberg JS. Human microbiota characterization in the course of renal transplantation. Am J Transplant (2014) 14(2):416–27. doi:10.1111/ajt.12588
18. Lee JR, Muthukumar T, Dadhania D, Toussaint NC, Ling L, Pamer E, et al. Gut microbial community structure and complications after kidney transplantation: a pilot study. Transplantation (2014) 98(7):697–705. doi:10.1097/ TP.0000000000000370
19. Rani A, Ranjan R, McGee HS, Metwally A, Hajjiri Z, Brennan DC, et al. A diverse virome in kidney transplant patients contains multiple viral subtypes with distinct polymorphisms. Sci Rep (2016) 6:33327. doi:10.1038/srep33327
20. Santiago-Rodriguez TM, Ly M, Bonilla N, Pride DT. The human urine virome is associated with urinary tract infections. Front Microbiol (2015) 6:14. doi:10.3389/fmicb.2015.00014
