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.

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.

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.

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.

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.

The Alpha variant was not associated with excess nosocomial SARS-CoV-2 infection in a multi-centre UK hospital study.

OBJECTIVES: Recently emerging SARS-CoV-2 variants have been associated with an increased rate of transmission within the community. We sought to determine whether this also resulted in increased transmission within hospitals. METHODS: We collected viral sequences and epidemiological data of patients with community and healthcare associated SARS-CoV-2 infections, sampled from 16th November 2020 to 10th January 2021, from nine hospitals participating in the COG-UK HOCI study. Outbreaks were identified using ward information, lineage and pairwise genetic differences between viral sequences. RESULTS: Mixed effects logistic regression analysis of 4184 sequences showed healthcare-acquired infections were no more likely to be identified as the Alpha variant than community acquired infections. Nosocomial outbreaks were investigated based on overlapping ward stay and SARS-CoV-2 genome sequence similarity. There was no significant difference in the number of patients involved in outbreaks caused by the Alpha variant compared to outbreaks caused by other lineages. CONCLUSIONS: We find no evidence to support it causing more nosocomial transmission than previous lineages. This suggests that the stringent infection prevention measures already in place in UK hospitals contained the spread of the Alpha variant as effectively as other less transmissible lineages, providing reassurance of their efficacy against emerging variants of concern.

SARS-CoV-2 lineage B.1.1.7 is associated with greater disease severity among hospitalised women but not men: multicentre cohort study.

BACKGROUND: SARS-CoV-2 lineage B.1.1.7 has been associated with an increased rate of transmission and disease severity among subjects testing positive in the community. Its impact on hospitalised patients is less well documented. METHODS: We collected viral sequences and clinical data of patients admitted with SARS-CoV-2 and hospital-onset COVID-19 infections (HOCIs), sampled 16 November 2020 to 10 January 2021, from eight hospitals participating in the COG-UK-HOCI study. Associations between the variant and the outcomes of all-cause mortality and intensive therapy unit (ITU) admission were evaluated using mixed effects Cox models adjusted by age, sex, comorbidities, care home residence, pregnancy and ethnicity. FINDINGS: Sequences were obtained from 2341 inpatients (HOCI cases=786) and analysis of clinical outcomes was carried out in 2147 inpatients with all data available. The HR for mortality of B.1.1.7 compared with other lineages was 1.01 (95% CI 0.79 to 1.28, p=0.94) and for ITU admission was 1.01 (95% CI 0.75 to 1.37, p=0.96). Analysis of sex-specific effects of B.1.1.7 identified increased risk of mortality (HR 1.30, 95% CI 0.95 to 1.78, p=0.096) and ITU admission (HR 1.82, 95% CI 1.15 to 2.90, p=0.011) in females infected with the variant but not males (mortality HR 0.82, 95% CI 0.61 to 1.10, p=0.177; ITU HR 0.74, 95% CI 0.52 to 1.04, p=0.086). INTERPRETATION: In common with smaller studies of patients hospitalised with SARS-CoV-2, we did not find an overall increase in mortality or ITU admission associated with B.1.1.7 compared with other lineages. However, women with B.1.1.7 may be at an increased risk of admission to intensive care and at modestly increased risk of mortality.

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.

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.

Exaggerated IL-17A activity in human in vivo recall responses discriminates active tuberculosis from latent infection and cured disease.

Host immune responses at the site of Mycobacterium tuberculosis infection can mediate pathogenesis of tuberculosis (TB) and onward transmission of infection. We hypothesized that pathological immune responses would be enriched at the site of host-pathogen interactions modeled by a standardized tuberculin skin test (TST) challenge in patients with active TB compared to those without disease, and interrogated immune responses by genome-wide transcriptional profiling. We show exaggerated interleukin-17A (IL-17A) and T helper 17 (TH17) responses among 48 individuals with active TB compared to 191 with latent TB infection, associated with increased neutrophil recruitment and matrix metalloproteinase-1 expression, both involved in TB pathogenesis. Curative antimicrobial treatment reversed these observed changes. Increased IL-1ß and IL-6 responses to mycobacterial stimulation were evident both in circulating monocytes and in molecular changes at the site of TST in individuals with active TB, supporting a model in which monocyte-derived IL-1ß and IL-6 promote TH17 differentiation within tissues. Modulation of these cytokine pathways may provide a rational strategy for host-directed therapy in active TB.

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.

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.

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.

Genetic Studies of Leptin Concentrations Implicate Leptin in the Regulation of Early Adiposity.

Leptin influences food intake by informing the brain about the status of body fat stores. Rare LEP mutations associated with congenital leptin deficiency cause severe early-onset obesity that can be mitigated by administering leptin. However, the role of genetic regulation of leptin in polygenic obesity remains poorly understood. We performed an exome-based analysis in up to 57,232 individuals of diverse ancestries to identify genetic variants that influence adiposity-adjusted leptin concentrations. We identify five novel variants, including four missense variants, in LEP, ZNF800, KLHL31, and ACTL9, and one intergenic variant near KLF14. The missense variant Val94Met (rs17151919) in LEP was common in individuals of African ancestry only, and its association with lower leptin concentrations was specific to this ancestry (P = 2 × 10-16, n = 3,901). Using in vitro analyses, we show that the Met94 allele decreases leptin secretion. We also show that the Met94 allele is associated with higher BMI in young African-ancestry children but not in adults, suggesting that leptin regulates early adiposity.

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.

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.

Fetal inheritance of chromosomally integrated human herpesvirus 6 predisposes the mother to pre-eclampsia.

Pre-eclampsia (typically characterized by new-onset hypertension and proteinuria in the second half of pregnancy) represents a major determinant of the global burden of disease1,2. Its pathophysiology involves placental dysfunction, but the mechanism is unclear. Viral infection can cause organ dysfunction, but its role in placentally related disorders of human pregnancy is unknown3. We addressed this using RNA sequencing metagenomics4-6 of placental samples from normal and complicated pregnancies. Here, we show that human herpesvirus 6 (HHV-6, A or B) RNA was detected in 6.1% of cases of pre-eclampsia and 2.2% of other pregnancies. Fetal genotyping demonstrated that 70% of samples with HHV-6 RNA in the placenta exhibited inherited, chromosomally integrated HHV-6 (iciHHV-6). We genotyped 467 pre-eclampsia cases and 3,854 controls and found an excess of iciHHV-6 in the cases (odds ratio of 2.8, 95% confidence intervals of 1.4-5.6, P = 0.008). We validated this finding by comparing iciHHV-6 in a further 740 cases with controls from large-scale population studies (odds ratio of 2.5, 95% confidence intervals of 1.4-4.4, P = 0.0013). We conclude that iciHHV-6 results in the transcription of viral RNA in the human placenta and predisposes the mother to pre-eclampsia.

Pathogenic Activation of Mesenchymal Stem Cells Is Induced by the Disease Microenvironment in Systemic Sclerosis.

OBJECTIVE: In systemic sclerosis (SSc), a persistent tissue repair process leads to progressive fibrosis of the skin and internal organs. The role of mesenchymal stem cells (MSCs), which characteristically initiate and regulate tissue repair, has not been fully evaluated. We undertook this study to investigate whether dividing metakaryotic MSCs are present in SSc skin and to examine whether exposure to the disease microenvironment activates MSCs and leads to transdifferentiation. METHODS: Skin biopsy material from patients with recent-onset diffuse SSc was examined by collagenase spread of 1-mm-thick surface-parallel sections, in order to identify dividing metakaryotic stem cells in each tissue plane. Adipose-derived MSCs from healthy controls were treated with dermal blister fluid (BF) from patients with diffuse SSc and profiled by next-generation sequencing, or they were evaluated for phenotypic changes relevant to SSc. Differential responses of dermal fibroblasts were studied in parallel. RESULTS: MSC-like cells undergoing active metakaryotic division were identified in SSc sections (but not control sections) most prominently in the deep dermis and adjacent to damaged microvessels, in both clinically involved and uninvolved skin. Furthermore, exposure to SSc BF caused selective MSC activation, inducing a myofibroblast signature, while reducing signatures of vascular repair and adipogenesis and enhancing migration and contractility. Microenvironmental factors implicated in inducing transdifferentiation included the profibrotic transforming growth factor ß, the presence of lactate, and mechanosensing, while the microenvironment Th2 cytokine, interleukin-31, enhanced osteogenic commitment (calcinosis). CONCLUSION: Dividing MSC-like cells are present in the SSc disease microenvironment where multiple factors, likely acting in concert, promote transdifferentiation and lead to a complex and resistant disease state.