Genomic investigations of unexplained acute hepatitis in children.

Since its first identification in Scotland, over 1,000 cases of unexplained paediatric hepatitis in children have been reported worldwide, including 278 cases in the UK1. Here we report an investigation of 38 cases, 66 age-matched immunocompetent controls and 21 immunocompromised comparator participants, using a combination of genomic, transcriptomic, proteomic and immunohistochemical methods. We detected high levels of adeno-associated virus 2 (AAV2) DNA in the liver, blood, plasma or stool from 27 of 28 cases. We found low levels of adenovirus (HAdV) and human herpesvirus 6B (HHV-6B) in 23 of 31 and 16 of 23, respectively, of the cases tested. By contrast, AAV2 was infrequently detected and at low titre in the blood or the liver from control children with HAdV, even when profoundly immunosuppressed. AAV2, HAdV and HHV-6 phylogeny excluded the emergence of novel strains in cases. Histological analyses of explanted livers showed enrichment for T cells and B lineage cells. Proteomic comparison of liver tissue from cases and healthy controls identified increased expression of HLA class 2, immunoglobulin variable regions and complement proteins. HAdV and AAV2 proteins were not detected in the livers. Instead, we identified AAV2 DNA complexes reflecting both HAdV-mediated and HHV-6B-mediated replication. We hypothesize that high levels of abnormal AAV2 replication products aided by HAdV and, in severe cases, HHV-6B may have triggered immune-mediated hepatic disease in genetically and immunologically predisposed children.

Genomic and geographical structure of human cytomegalovirus.

Human cytomegalovirus (CMV) has infected humans since the origin of our species and currently infects most of the world's population. Variability between CMV genomes is the highest of any human herpesvirus, yet large portions of the genome are conserved. Here, we show that the genome encodes 74 regions of relatively high variability each with 2 to 8 alleles. We then identified two patterns in the CMV genome. Conserved parts of the genome and a minority (32) of variable regions show geographic population structure with evidence for African or European clustering, although hybrid strains are present. We find no evidence that geographic segregation has been driven by host immune pressure affecting known antigenic sites. Forty-two variable regions show no geographical structure, with similar allele distributions across different continental populations. These "nongeographical" regions are significantly enriched for genes encoding immunomodulatory functions suggesting a core functional importance. We hypothesize that at least two CMV founder populations account for the geographical differences that are largely seen in the conserved portions of the genome, although the timing of separation and direction of spread between the two are not clear. In contrast, the similar allele frequencies among 42 variable regions of the genome, irrespective of geographical origin, are indicative of a second evolutionary process, namely balancing selection that may preserve properties critical to CMV biological function. Given that genetic differences between CMVs are postulated to alter immunogenicity and potentially function, understanding these two evolutionary processes could contribute important information for the development of globally effective vaccines and the identification of novel drug targets.

Persistent low-level variants in a subset of viral genes are highly predictive of poor outcome in immunocompromised patients with cytomegalovirus infection.

BACKGROUND: Human cytomegalovirus is the most common and serious opportunistic infection after solid organ and haematopoietic stem cell transplantation. In this study, we used whole-genome cytomegalovirus data to investigate viral factors associated with the clinical outcome. METHODS: We sequenced cytomegalovirus samples from 16 immunocompromised paediatric patients with persistent viraemia. 8/16 patients died of complications due to cytomegalovirus infection. We also sequenced samples from 35 infected solid organ adult recipients of whom one died with cytomegalovirus infection. RESULTS: We showed that samples from both groups have fixed variants at resistance sites and mixed infections. NGS sequencing also revealed non-fixed variants at resistance sites in most of the patients who died (6/9). A machine learning approach identified 10 genes with non-fixed variants in these patients. These genes formed a viral signature which discriminated patients with cytomegalovirus infection who died from those that survived with high accuracy (AUC=0.96). Lymphocyte numbers for a subset of patients showed no recovery post-transplant in the patients who died. CONCLUSIONS: We hypothesise that the viral signature identified in this study may be a useful biomarker for poor response to antiviral drug treatment and indirectly for poor T cell function, potentially identifying early, those patients requiring non-pharmacological interventions.

Evolutionary History of Endogenous Human Herpesvirus 6 Reflects Human Migration out of Africa.

Human herpesvirus 6A and 6B (HHV-6) can integrate into the germline, and as a result, ∼70 million people harbor the genome of one of these viruses in every cell of their body. Until now, it has been largely unknown if 1) these integrations are ancient, 2) if they still occur, and 3) whether circulating virus strains differ from integrated ones. Here, we used next-generation sequencing and mining of public human genome data sets to generate the largest and most diverse collection of circulating and integrated HHV-6 genomes studied to date. In genomes of geographically dispersed, only distantly related people, we identified clades of integrated viruses that originated from a single ancestral event, confirming this with fluorescent in situ hybridization to directly observe the integration locus. In contrast to HHV-6B, circulating and integrated HHV-6A sequences form distinct clades, arguing against ongoing integration of circulating HHV-6A or "reactivation" of integrated HHV-6A. Taken together, our study provides the first comprehensive picture of the evolution of HHV-6, and reveals that integration of heritable HHV-6 has occurred since the time of, if not before, human migrations out of Africa.

Epstein-Barr virus (EBV) deletions as biomarkers of response to treatment of chronic active EBV.

Chronic active Epstein-Barr virus (CAEBV) disease is a rare condition characterised by persistent EBV infection in previously healthy individuals. Defective EBV genomes were found in East Asian patients with CAEBV. In the present study, we sequenced 14 blood EBV samples from three UK patients with CAEBV, comparing the results with saliva CAEBV samples and other conditions. We observed EBV deletions in blood, some of which may disrupt viral replication, but not saliva in CAEBV. Deletions were lost overtime after successful treatment. These findings are compatible with CAEBV being associated with the evolution and persistence of EBV+ haematological clones that are lost on successful treatment.

Blood Transcriptomic Stratification of Short-term Risk in Contacts of Tuberculosis.

BACKGROUND: The highest risk of tuberculosis arises in the first few months after exposure. We reasoned that this risk reflects incipient disease among tuberculosis contacts. Blood transcriptional biomarkers of tuberculosis may predate clinical diagnosis, suggesting they offer improved sensitivity to detect subclinical incipient disease. Therefore, we sought to test the hypothesis that refined blood transcriptional biomarkers of active tuberculosis will improve stratification of short-term disease risk in tuberculosis contacts. METHODS: We combined analysis of previously published blood transcriptomic data with new data from a prospective human immunodeficiency virus (HIV)-negative UK cohort of 333 tuberculosis contacts. We used stability selection as an alternative computational approach to identify an optimal signature for short-term risk of active tuberculosis and evaluated its predictive value in independent cohorts. RESULTS: In a previously published HIV-negative South African case-control study of patients with asymptomatic Mycobacterium tuberculosis infection, a novel 3-gene transcriptional signature comprising BATF2, GBP5, and SCARF1 achieved a positive predictive value (PPV) of 23% for progression to active tuberculosis within 90 days. In a new UK cohort of 333 HIV-negative tuberculosis contacts with a median follow-up of 346 days, this signature achieved a PPV of 50% (95% confidence interval [CI], 15.7-84.3) and negative predictive value of 99.3% (95% CI, 97.5-99.9). By comparison, peripheral blood interferon gamma release assays in the same cohort achieved a PPV of 5.6% (95% CI, 2.1-11.8). CONCLUSIONS: This blood transcriptional signature provides unprecedented opportunities to target therapy among tuberculosis contacts with greatest risk of incident disease.

Blood transcriptomic discrimination of bacterial and viral infections in the emergency department: a multi-cohort observational validation study.

BACKGROUND: There is an urgent need to develop biomarkers that stratify risk of bacterial infection in order to support antimicrobial stewardship in emergency hospital admissions. METHODS: We used computational machine learning to derive a rule-out blood transcriptomic signature of bacterial infection (SeptiCyte™ TRIAGE) from eight published case-control studies. We then validated this signature by itself in independent case-control data from more than 1500 samples in total, and in combination with our previously published signature for viral infections (SeptiCyte™ VIRUS) using pooled data from a further 1088 samples. Finally, we tested the performance of these signatures in a prospective observational cohort of emergency department (ED) patients with fever, and we used the combined SeptiCyte™ signature in a mixture modelling approach to estimate the prevalence of bacterial and viral infections in febrile ED patients without microbiological diagnoses. RESULTS: The combination of SeptiCyte™ TRIAGE with our published signature for viral infections (SeptiCyte™ VIRUS) discriminated bacterial and viral infections in febrile ED patients, with a receiver operating characteristic area under the curve of 0.95 (95% confidence interval 0.90-1), compared to 0.79 (0.68-0.91) for WCC and 0.73 (0.61-0.86) for CRP. At pre-test probabilities 0.35 and 0.72, the combined SeptiCyte™ score achieved a negative predictive value for bacterial infection of 0.97 (0.90-0.99) and 0.86 (0.64-0.96), compared to 0.90 (0.80-0.94) and 0.66 (0.48-0.79) for WCC and 0.88 (0.69-0.95) and 0.60 (0.31-0.72) for CRP. In a mixture modelling approach, the combined SeptiCyte™ score estimated that 24% of febrile ED cases receiving antibacterials without a microbiological diagnosis were due to viral infections. Our analysis also suggested that a proportion of patients with bacterial infection recovered without antibacterials. CONCLUSIONS: Blood transcriptional biomarkers offer exciting opportunities to support precision antibacterial prescribing in ED and improve diagnostic classification of patients without microbiologically confirmed infections.

Correction to: Blood transcriptomic discrimination of bacterial and viral infections in the emergency department: a multi-cohort observational validation study.

An amendment to this paper has been published and can be accessed via the original article.

Interleukin-31 promotes pathogenic mechanisms underlying skin and lung fibrosis in scleroderma.

OBJECTIVES: Cytokines released by infiltrating T cells may promote mechanisms leading to fibrosis in scleroderma. The aim of this study was to investigate the role of the Th2 cytokine IL-31, and its receptor IL-31RA, in scleroderma skin and lung fibrosis. METHODS: IL-31 was measured by ELISA of plasma, and by immunochemistry of fibrotic skin and lung tissue of scleroderma patients. The receptor, IL-31RA, was assayed by qPCR of tissue resident cells. Next-generation sequencing was used to profile the responses of normal skin fibroblasts to IL-31. In wild-type Balb/c mice, IL-31 was administered by subcutaneous mini pump, with or without additional TGFß, and the fibrotic reaction measured by histology and ELISA of plasma. RESULTS: IL-31 was present at high levels in plasma and fibrotic skin and lung lesions in a subset of scleroderma patients, and the receptor overexpressed by downstream cells relevant to the disease process, including skin and lung fibroblasts, through loss of epigenetic regulation by miR326. In skin fibroblasts, IL-31 induced next generation sequencing profiles associated with cellular growth and proliferation, anaerobic metabolism and mineralization, and negatively associated with angiogenesis and vascular repair, as well as promoting phenotype changes including migration and collagen protein release via pSTAT3, resembling the activation state in the disease. In mice, IL-31 induced skin and lung fibrosis. No synergy was seen with TGFß, which supressed IL-31RA. CONCLUSION: IL-31/IL-31RA is confirmed as a candidate pro-fibrotic pathway, which may contribute to skin and lung fibrosis in a subset of scleroderma patients.

New insights into the genetics of primary open-angle glaucoma based on meta-analyses of intraocular pressure and optic disc characteristics.

Primary open-angle glaucoma (POAG), the most common optic neuropathy, is a heritable disease. Siblings of POAG cases have a ten-fold increased risk of developing the disease. Intraocular pressure (IOP) and optic nerve head characteristics are used clinically to predict POAG risk. We conducted a genome-wide association meta-analysis of IOP and optic disc parameters and validated our findings in multiple sets of POAG cases and controls. Using imputation to the 1000 genomes (1000G) reference set, we identified 9 new genomic regions associated with vertical cup-disc ratio (VCDR) and 1 new region associated with IOP. Additionally, we found 5 novel loci for optic nerve cup area and 6 for disc area. Previously it was assumed that genetic variation influenced POAG either through IOP or via changes to the optic nerve head; here we present evidence that some genomic regions affect both IOP and the disc parameters. We characterized the effect of the novel loci through pathway analysis and found that pathways involved are not entirely distinct as assumed so far. Further, we identified a novel association between CDKN1A and POAG. Using a zebrafish model we show that six6b (associated with POAG and optic nerve head variation) alters the expression of cdkn1a. In summary, we have identified several novel genes influencing the major clinical risk predictors of POAG and showed that genetic variation in CDKN1A is important in POAG risk.

Sequence Variation of Epstein-Barr Virus: Viral Types, Geography, Codon Usage, and Diseases.

One hundred thirty-eight new Epstein-Barr virus (EBV) genome sequences have been determined. One hundred twenty-five of these and 116 from previous reports were combined to produce a multiple-sequence alignment of 241 EBV genomes, which we have used to analyze variation within the viral genome. The type 1/type 2 classification of EBV remains the major form of variation and is defined mostly by EBNA2 and EBNA3, but the type 2 single-nucleotide polymorphisms (SNPs) at the EBNA3 locus extend into the adjacent gp350 and gp42 genes, whose products mediate infection of B cells by EBV. A small insertion within the BART microRNA region of the genome was present in 21 EBV strains. EBV from saliva of U.S. patients with chronic active EBV infection aligned with the wild-type EBV genome with no evidence of WZhet rearrangements. The V3 polymorphism in the Zp promoter for BZLF1 was found to be frequent in nasopharyngeal carcinoma cases from both Hong Kong and Indonesia. Codon usage was found to differ between latent and lytic cycle EBV genes, and the main forms of variation of the EBNA1 protein have been identified.IMPORTANCE Epstein-Barr virus causes most cases of infectious mononucleosis and posttransplant lymphoproliferative disease. It contributes to several types of cancer, including Hodgkin's lymphoma, Burkitt's lymphoma, diffuse large B cell lymphoma, nasopharyngeal carcinoma, and gastric carcinoma. EBV genome variation is important because some of the diseases associated with EBV have very different incidences in different populations and geographic regions, and differences in the EBV genome might contribute to these diseases. Some specific EBV genome alterations that appear to be significant in EBV-associated cancers are already known, and current efforts to make an EBV vaccine and antiviral drugs should also take account of sequence differences in the proteins used as targets.

A meta-analysis of 120 246 individuals identifies 18 new loci for fibrinogen concentration.

Genome-wide association studies have previously identified 23 genetic loci associated with circulating fibrinogen concentration. These studies used HapMap imputation and did not examine the X-chromosome. 1000 Genomes imputation provides better coverage of uncommon variants, and includes indels. We conducted a genome-wide association analysis of 34 studies imputed to the 1000 Genomes Project reference panel and including ∼120 000 participants of European ancestry (95 806 participants with data on the X-chromosome). Approximately 10.7 million single-nucleotide polymorphisms and 1.2 million indels were examined. We identified 41 genome-wide significant fibrinogen loci; of which, 18 were newly identified. There were no genome-wide significant signals on the X-chromosome. The lead variants of five significant loci were indels. We further identified six additional independent signals, including three rare variants, at two previously characterized loci: FGB and IRF1. Together the 41 loci explain 3% of the variance in plasma fibrinogen concentration.

A genome-wide association study of intra-ocular pressure suggests a novel association in the gene FAM125B in the TwinsUK cohort.

Glaucoma is a major cause of blindness in the world. To date, common genetic variants associated with glaucoma only explain a small proportion of its heritability. We performed a genome-wide association study of intra-ocular pressure (IOP), an underlying endophenotype for glaucoma. The discovery phase of the study was carried out in the TwinsUK cohort (N = 2774) analyzing association between IOP and single nucleotide polymorphisms (SNPs) imputed to HapMap2. The results were validated in 12 independent replication cohorts of European ancestry (combined N = 22 789) that were a part of the International Glaucoma Genetics Consortium. Expression quantitative trait locus (eQTL) analyses of the significantly associated SNPs were performed using data from the Multiple Tissue Human Expression Resource (MuTHER) Study. In the TwinsUK cohort, IOP was significantly associated with a number of SNPs at 9q33.3 (P = 3.48 × 10(-8) for rs2286885, the most significantly associated SNP at this locus), within the genomic sequence of the FAM125B gene. Independent replication in a composite panel of 12 cohorts revealed consistent direction of effect and significant association (P = 0.003, for fixed-effect meta-analysis). Suggestive evidence for an eQTL effect of rs2286885 was observed for one of the probes targeting the coding region of the FAM125B gene. This gene codes for a component of a membrane complex involved in vesicular trafficking process, a function similar to that of the Caveolin genes (CAV1 and CAV2) which have previously been associated with primary open-angle glaucoma. This study suggests a novel association between SNPs in FAM125B and IOP in the TwinsUK cohort, though further studies to elucidate the functional role of this gene in glaucoma are necessary.

Digenic inheritance involving a muscle-specific protein kinase and the giant titin protein causes a skeletal muscle myopathy.

In digenic inheritance, pathogenic variants in two genes must be inherited together to cause disease. Only very few examples of digenic inheritance have been described in the neuromuscular disease field. Here we show that predicted deleterious variants in SRPK3, encoding the X-linked serine/argenine protein kinase 3, lead to a progressive early onset skeletal muscle myopathy only when in combination with heterozygous variants in the TTN gene. The co-occurrence of predicted deleterious SRPK3/TTN variants was not seen among 76,702 healthy male individuals, and statistical modeling strongly supported digenic inheritance as the best-fitting model. Furthermore, double-mutant zebrafish (srpk3-/-; ttn.1+/-) replicated the myopathic phenotype and showed myofibrillar disorganization. Transcriptome data suggest that the interaction of srpk3 and ttn.1 in zebrafish occurs at a post-transcriptional level. We propose that digenic inheritance of deleterious changes impacting both the protein kinase SRPK3 and the giant muscle protein titin causes a skeletal myopathy and might serve as a model for other genetic diseases.

Candidate gene study of macular response to supplemental lutein and zeaxanthin.

Supplementation with carotenoids is proposed to protect against age-related macular degeneration. There is, however, considerable variability in retinal macular pigment response, which may be due to underlying genetic variation. The purpose of this study was to determine whether genetic factors, which have been previously associated with cross-sectional macular pigment levels in the retina or serum lutein, also influence response to supplementation. To this end we conducted an association study in 310 subjects from the TwinsUK cohort between variants in 8 candidate genes and serum lutein and retinal macular pigment optical density (MPOD) levels before and after supplementation. Four variants were associated with MPOD response to supplementation (p < 0.05): rs11057841 (SCARB1), rs4926339 (RPE65), rs1929841 (ABCA1) and rs174534 (FADS1). We also confirmed previous associations between rs6564851 near BMCO1 (p < 0.001) and rs11057841 within SCARB1 (p = 0.01) and baseline measures of serum lutein; while the latter was also associated with MPOD response, none of the BMCO1 variants were. Finally, there was evidence for association between variants near RPE65 and ELOVL2 and changes in lutein concentration after supplementation. This study is the first to show association between genetic variants and response to carotenoids supplementation. Our findings suggest an important link between MP response and the biological processes of carotenoids transport and fatty acid metabolism.

Age of myopia onset in a British population-based twin cohort.

PURPOSE: School-age myopia is becoming more common in Asia and North America; data from the United Kingdom has suggested a significant amount of myopia develops after the age of 17 years. Age of spectacle wear has been used as a proxy of myopia severity in a recent large genome-wide association study. The purpose of this study was to examine the age of onset of spectacle wear in a large British twin cohort, to examine the reliability and reproducibility of self-reported age of onset as a proxy measure of myopia severity, and to see if there is evidence in the UK of a rising prevalence of myopia. METHODS: Non-cycloplegic autorefraction was performed on over 6000 subjects from the Twins UK cohort, a large, well-characterized volunteer cohort of British, predominantly Caucasian female twins, between 1998 and 2010. Questionnaires asking age of first spectacle wear were conducted in 2003 and 2008. Myopia was defined as worse than or equal to -1.00 Dioptres, and adult onset myopia as occurring on or after the age of 17 years. RESULTS: Autorefractive data was available on 6097 participants at a mean age of 53 years. The mean S.E. was -0.36 D (S.D. 2.67, range -25.13 to +9.38). 1705 subjects (28%) were myopic with a mean refractive error of -3.54 (S.D. 2.51, range -25.13 to -1.00) and the median age of first glasses wear was 15 years (mean 18.4 years, S.D. 12.24, range 0-74). Of those who provided an age at which they first wore glasses in both questionnaire sources (n = 628), there was median difference in response of 0 years (S.D. 7.18, mean 0.7, maximum 53). A statistically significant cohort effect for increased myopia prevalence across a range of age groups between 1998-1999 and 2008-2010 was identified, with myopia prevalence increasing from 27% to 34% in those aged 50-54 and from 16% to 32% in those aged 55-59. CONCLUSIONS: Almost half the myopes in this UK-based population wore glasses after the age of 17; further research into adult-onset myopia is required. Although self-reported age of glasses is reproducible and reflects severity, it only explains approximately 15% of the variance of spherical equivalent, so is a rough proxy of refractive error, but still may be useful in large-scale population studies without access to refraction. We have demonstrated a significant cohort effect for increased myopia prevalence in the UK population over a 10-year period.

Haplotype assignment of longitudinal viral deep sequencing data using covariation of variant frequencies.

Longitudinal deep sequencing of viruses can provide detailed information about intra-host evolutionary dynamics including how viruses interact with and transmit between hosts. Many analyses require haplotype reconstruction, identifying which variants are co-located on the same genomic element. Most current methods to perform this reconstruction are based on a high density of variants and cannot perform this reconstruction for slowly evolving viruses. We present a new approach, HaROLD (HAplotype Reconstruction Of Longitudinal Deep sequencing data), which performs this reconstruction based on identifying co-varying variant frequencies using a probabilistic framework. We illustrate HaROLD on both RNA and DNA viruses with synthetic Illumina paired read data created from mixed human cytomegalovirus (HCMV) and norovirus genomes, and clinical datasets of HCMV and norovirus samples, demonstrating high accuracy, especially when longitudinal samples are available.

Pericyte derived chemokines amplify neutrophil recruitment across the cerebrovascular endothelial barrier.

Excessive neutrophil extravasation can drive immunopathology, exemplified in pyogenic meningitis caused by Streptococcus pneumoniae infection. Insufficient knowledge of the mechanisms that amplify neutrophil extravasation has limited innovation in therapeutic targeting of neutrophil mediated pathology. Attention has focussed on neutrophil interactions with endothelia, but data from mouse models also point to a role for the underlying pericyte layer, as well as perivascular macrophages, the only other cell type found within the perivascular space in the cerebral microvasculature. We tested the hypothesis that human brain vascular pericytes (HBVP) contribute to neutrophil extravasation in a transwell model of the cerebral post-capillary venule. We show that pericytes augment endothelial barrier formation. In response to inflammatory cues, they significantly enhance neutrophil transmigration across the endothelial barrier, without increasing the permeability to small molecules. In our model, neither pericytes nor endothelia responded directly to bacterial stimulation. Instead, we show that paracrine signalling by multiple cytokines from monocyte derived macrophages drives transcriptional upregulation of multiple neutrophil chemokines by pericytes. Pericyte mediated amplification of neutrophil transmigration was independent of transcriptional responses by endothelia, but could be mediated by direct chemokine translocation across the endothelial barrier. Our data support a model in which microbial sensing by perivascular macrophages generates an inflammatory cascade where pericytes serve to amplify production of neutrophil chemokines that are translocated across the endothelial barrier to act directly on circulating neutrophils. In view of the striking redundancy in inflammatory cytokines that stimulate pericytes and in the neutrophil chemokines they produce, we propose that the mechanism of chemokine translocation may offer the most effective therapeutic target to reduce neutrophil mediated pathology in pyogenic meningitis.

GRACy: A tool for analysing human cytomegalovirus sequence data.

Modern DNA sequencing has instituted a new era in human cytomegalovirus (HCMV) genomics. A key development has been the ability to determine the genome sequences of HCMV strains directly from clinical material. This involves the application of complex and often non-standardized bioinformatics approaches to analysing data of variable quality in a process that requires substantial manual intervention. To relieve this bottleneck, we have developed GRACy (Genome Reconstruction and Annotation of Cytomegalovirus), an easy-to-use toolkit for analysing HCMV sequence data. GRACy automates and integrates modules for read filtering, genotyping, genome assembly, genome annotation, variant analysis, and data submission. These modules were tested extensively on simulated and experimental data and outperformed generic approaches. GRACy is written in Python and is embedded in a graphical user interface with all required dependencies installed by a single command. It runs on the Linux operating system and is designed to allow the future implementation of a cross-platform version. GRACy is distributed under a GPL 3.0 license and is freely available at https://bioinformatics.cvr.ac.uk/software/ with the manual and a test dataset.

Using Whole Genome Sequences to Investigate Adenovirus Outbreaks in a Hematopoietic Stem Cell Transplant Unit.

A recent surge in human mastadenovirus (HAdV) cases, including five deaths, amongst a haematopoietic stem cell transplant population led us to use whole genome sequencing (WGS) to investigate. We compared sequences from 37 patients collected over a 20-month period with sequences from GenBank and our own database of HAdVs. Maximum likelihood trees and pairwise differences were used to evaluate genotypic relationships, paired with the epidemiological data from routine infection prevention and control (IPC) records and hospital activity data. During this time period, two formal outbreaks had been declared by IPC, while WGS detected nine monophyletic clusters, seven were corroborated by epidemiological evidence and by comparison of single-nucleotide polymorphisms. One of the formal outbreaks was confirmed, and the other was not. Of the five HAdV-associated deaths, three were unlinked and the remaining two considered the source of transmission. Mixed infection was frequent (10%), providing a sentinel source of recombination and superinfection. Immunosuppressed patients harboring a high rate of HAdV positivity require comprehensive surveillance. As a consequence of these findings, HAdV WGS is being incorporated routinely into clinical practice to influence IPC policy contemporaneously.