Evaluating assembly and variant calling software for strain-resolved analysis of large DNA viruses.

Infection with human cytomegalovirus (HCMV) can cause severe complications in immunocompromised individuals and congenitally infected children. Characterizing heterogeneous viral populations and their evolution by high-throughput sequencing of clinical specimens requires the accurate assembly of individual strains or sequence variants and suitable variant calling methods. However, the performance of most methods has not been assessed for populations composed of low divergent viral strains with large genomes, such as HCMV. In an extensive benchmarking study, we evaluated 15 assemblers and 6 variant callers on 10 lab-generated benchmark data sets created with two different library preparation protocols, to identify best practices and challenges for analyzing such data. Most assemblers, especially metaSPAdes and IVA, performed well across a range of metrics in recovering abundant strains. However, only one, Savage, recovered low abundant strains and in a highly fragmented manner. Two variant callers, LoFreq and VarScan2, excelled across all strain abundances. Both shared a large fraction of false positive variant calls, which were strongly enriched in T to G changes in a 'G.G' context. The magnitude of this context-dependent systematic error is linked to the experimental protocol. We provide all benchmarking data, results and the entire benchmarking workflow named QuasiModo, Quasispecies Metric determination on omics, under the GNU General Public License v3.0 (https://github.com/hzi-bifo/Quasimodo), to enable full reproducibility and further benchmarking on these and other data.

Microbial exposure drives polyclonal expansion of innate γδ T cells immediately after birth.

Starting at birth, the immune system of newborns and children encounters and is influenced by environmental challenges. It is still not completely understood how γδ T cells emerge and adapt during early life. Studying the composition of T cell receptors (TCRs) using next-generation sequencing (NGS) in neonates, infants, and children can provide valuable insights into the adaptation of T cell subsets. To investigate how neonatal γδ T cell repertoires are shaped by microbial exposure after birth, we monitored the γ-chain (TRG) and δ-chain (TRD) repertoires of peripheral blood T cells in newborns, infants, and young children from Europe and sub-Saharan Africa. We identified a set of TRG and TRD sequences that were shared by all children from Europe and Africa. These were primarily public clones, characterized by simple rearrangements of Vγ9 and Vδ2 chains with low junctional diversity and usage of non-TRDJ1 gene segments, reminiscent of early ontogenetic subsets of γδ T cells. Further profiling revealed that these innate, public Vγ9Vδ2+ T cells underwent an immediate TCR-driven polyclonal proliferation within the first 4 wk of life. In contrast, γδ T cells using Vδ1+ and Vδ3+TRD rearrangements did not significantly expand after birth. However, different environmental cues may lead to the observed increase of Vδ1+ and Vδ3+TRD sequences in the majority of African children. In summary, we show how dynamic γδ TCR repertoires develop directly after birth and present important differences among γδ T cell subsets.

Diagnostic Yield and Therapeutic Consequences of Targeted Next-Generation Sequencing in Sporadic Primary Immunodeficiency.

INTRODUCTION: Primary immunodeficiencies (PIDs) are a heterogeneous group of disorders characterized by increased susceptibility to infections, immune dysregulation, and/or malignancy. Genetic studies, especially during the last decade, led to a better understanding of the pathogenesis of primary immunodeficiencies and contributed to their classification into distinct monogenic disorders falling under one of the >430 currently known inborn errors of immunity (IEI). The growing availability of molecular genetic testing resulted in the increasing identification of patients with IEI. Here, we evaluated the diagnostic yield and the clinical consequences of targeted next-generation sequencing (tNGS) in a cohort of 294 primary immunodeficiency patients, primarily consisting of cases with sporadic primary antibody deficiency. METHOD: We have custom designed a tNGS panel to sequence a cohort of PID patients. Agilent's HaloPlex Target Enrichment System for Illumina was used for DNA target enrichment. RESULTS: tNGS identified a definite or predicted pathogenic variant in 15.3% of patients. The highest diagnostic rate was observed among patients with combined immunodeficiency or immune dysregulation, for whom genetic diagnosis may affect therapeutic decision-making. CONCLUSION: Next-generation sequencing has changed diagnostic assignment and paved the way for targeted therapeutic intervention with agents directed at reverting the disease-causing molecular abnormality or its pathophysiological consequences. Therefore, such targeted therapies and identifying the genetic basis of PID can be essential for patients with manifested immune dysregulation as conventional immunomodulatory regimens may exert an immunosuppressive effect, aggravating their immunodeficiency or may only inadequately control autoimmune or lymphoproliferative manifestations.

Murine Cytomegalovirus M25 Proteins Sequester the Tumor Suppressor Protein p53 in Nuclear Accumulations.

To ensure productive infection, herpesviruses utilize tegument proteins and nonstructural regulatory proteins to counteract cellular defense mechanisms and to reprogram cellular pathways. The M25 proteins of mouse cytomegalovirus (MCMV) belong to the betaherpesvirus UL25 gene family that encodes viral proteins implicated with regulatory functions. Through affinity purification and mass spectrometric analysis, we discovered the tumor suppressor protein p53 as a host factor interacting with the M25 proteins. M25-p53 interaction in infected and transfected cells was confirmed by coimmunoprecipitation. Moreover, the proteins colocalized in nuclear dot-like structures upon both infection and inducible expression of the two M25 isoforms. p53 accumulated in wild-type MCMV-infected cells, while this did not occur upon infection with a mutant lacking the M25 gene. Both M25 proteins were able to mediate the effect, identifying them as the first CMV proteins responsible for p53 accumulation during infection. Interaction with M25 proteins led to substantial prolongation of the half-life of p53. In contrast to the higher abundance of the p53 protein in wild-type MCMV-infected cells, the transcript levels of the prominent p53 target genes Cdkn1a and Mdm2 were diminished compared to cells infected with the ΔM25 mutant, and this was associated with reduced binding of p53 to responsive elements within the respective promoters. Notably, the productivity of the M25 deletion mutant was partially rescued on p53-negative fibroblasts. We propose that the MCMV M25 proteins sequester p53 molecules in the nucleus of infected cells, reducing their availability for activating a subset of p53-regulated genes, thereby dampening the antiviral role of p53.IMPORTANCE Host cells use a number of factors to defend against viral infection. Viruses are, however, in an arms race with their host cells to overcome these defense mechanisms. The tumor suppressor protein p53 is an important sensor of cell stress induced by oncogenic insults or viral infections, which upon activation induces various pathways to ensure the integrity of cells. Viruses have to counteract many functions of p53, but complex DNA viruses such as cytomegaloviruses may also utilize some p53 functions for their own benefit. In this study, we discovered that the M25 proteins of mouse cytomegalovirus interact with p53 and mediate its accumulation during infection. Interaction with the M25 proteins sequesters p53 molecules in nuclear dot-like structures, limiting their availability for activation of a subset of p53-regulated target genes. Understanding the interaction between viral proteins and p53 may allow to develop new therapeutic strategies against cytomegalovirus and other viruses.

The Essential Human Cytomegalovirus Proteins pUL77 and pUL93 Are Structural Components Necessary for Viral Genome Encapsidation.

UNLABELLED: Several essential viral proteins are proposed to participate in genome encapsidation of human cytomegalovirus (HCMV), among them pUL77 and pUL93, which remain largely uncharacterized. To gain insight into their properties, we generated an HCMV mutant expressing a pUL77-monomeric enhanced green fluorescent protein (mGFP) fusion protein and a pUL93-specific antibody. Immunoblotting demonstrated that both proteins are incorporated into capsids and virions. Conversely to data suggesting internal translation initiation sites within the UL93 open reading frame (ORF), we provide evidence that pUL93 synthesis commences at the first start codon. In infected cells, pUL77-mGFP was found in nuclear replication compartments and dot-like structures, colocalizing with capsid proteins. Immunogold labeling of nuclear capsids revealed that pUL77 is present on A, B, and C capsids. Pulldown of pUL77-mGFP revealed copurification of pUL93, indicating interaction between these proteins, which still occurred when capsid formation was prevented. Correct subnuclear distribution of pUL77-mGFP required pUL93 as well as the major capsid protein (and thus probably the presence of capsids), but not the tegument protein pp150 or the encapsidation protein pUL52, demonstrating that pUL77 nuclear targeting occurs independently of the formation of DNA-filled capsids. When pUL77 or pUL93 was missing, generation of unit-length genomes was not observed, and only empty B capsids were produced. Taken together, these results show that pUL77 and pUL93 are capsid constituents needed for HCMV genome encapsidation. Therefore, the task of pUL77 seems to differ from that of its alphaherpesvirus orthologue pUL25, which exerts its function subsequent to genome cleavage-packaging. IMPORTANCE: The essential HCMV proteins pUL77 and pUL93 were suggested to be involved in viral genome cleavage-packaging but are poorly characterized both biochemically and functionally. By producing a monoclonal antibody against pUL93 and generating an HCMV mutant in which pUL77 is fused to a fluorescent protein, we show that pUL77 and pUL93 are capsid constituents, with pUL77 being similarly abundant on all capsid types. Each protein is required for genome encapsidation, as the absence of either pUL77 or pUL93 results in a genome packaging defect with the formation of empty capsids only. This distinguishes pUL77 from its alphaherpesvirus orthologue pUL25, which is enriched on DNA-filled capsids and exerts its function after the viral DNA is packaged. Our data for the first time describe an HCMV mutant with a fluorescent capsid and provide insight into the roles of pUL77 and pUL93, thus contributing to a better understanding of the HCMV encapsidation network.

Expression of the human cytomegalovirus UL11 glycoprotein in viral infection and evaluation of its effect on virus-specific CD8 T cells.

UNLABELLED: The human cytomegalovirus (CMV) UL11 open reading frame (ORF) encodes a putative type I transmembrane glycoprotein which displays remarkable amino acid sequence variability among different CMV isolates, suggesting that it represents an important virulence factor. In a previous study, we have shown that UL11 can interact with the cellular receptor tyrosine phosphatase CD45, which has a central role for signal transduction in T cells, and treatment of T cells with large amounts of a soluble UL11 protein inhibited their proliferation. In order to analyze UL11 expression in CMV-infected cells, we constructed CMV recombinants whose genomes either encode tagged UL11 versions or carry a stop mutation in the UL11 ORF. Moreover, we examined whether UL11 affects the function of virus-specific cytotoxic T lymphocytes (CTLs). We found that the UL11 ORF gives rise to several proteins due to both posttranslational modification and alternative translation initiation sites. Biotin labeling of surface proteins on infected cells indicated that only highly glycosylated UL11 forms are present at the plasma membrane, whereas less glycosylated UL11 forms were found in the endoplasmic reticulum. We did not find evidence of UL11 cleavage or secretion of a soluble UL11 version. Cocultivation of CTLs recognizing different CMV epitopes with fibroblasts infected with a UL11 deletion mutant or the parental strain revealed that under the conditions applied UL11 did not influence the activation of CMV-specific CD8 T cells. For further studies, we propose to investigate the interaction of UL11 with CD45 and the functional consequences in other immune cells expressing CD45. IMPORTANCE: Human cytomegalovirus (CMV) belongs to those viruses that extensively interfere with the host immune response, yet the precise function of many putative immunomodulatory CMV proteins remains elusive. Previously, we have shown that the CMV UL11 protein interacts with the leukocyte common antigen CD45, a cellular receptor tyrosine phosphatase with a central role for signal transduction in T cells. Here, we examined the proteins expressed by the UL11 gene in CMV-infected cells and found that at least one form of UL11 is present at the cell surface, enabling it to interact with CD45 on immune cells. Surprisingly, CMV-expressed UL11 did not affect the activity of virus-specific CD8 T cells. This finding warrants investigation of the impact of UL11 on CD45 functions in other leukocyte subpopulations.

Viral modulation of type II interferon increases T cell adhesion and virus spread.

During primary varicella zoster virus (VZV) infection, infected lymphocytes drive primary viremia, causing systemic dissemination throughout the host, including the skin. This results in cytokine expression, including interferons (IFNs), which partly limit infection. VZV also spreads from skin keratinocytes to lymphocytes prior to secondary viremia. It is not clear how VZV achieves this while evading the cytokine response. Here, we show that VZV glycoprotein C (gC) binds IFN-γ and modifies its activity, increasing the expression of a subset of IFN-stimulated genes (ISGs), including intercellular adhesion molecule 1 (ICAM1), chemokines and immunomodulatory genes. The higher ICAM1 protein level at the plasma membrane of keratinocytes facilitates lymphocyte function-associated antigen 1-dependent T cell adhesion and expression of gC during infection increases VZV spread to peripheral blood mononuclear cells. This constitutes the discovery of a strategy to modulate IFN-γ activity, upregulating a subset of ISGs, promoting enhanced lymphocyte adhesion and virus spread.

Detecting patches of protein sites of influenza A viruses under positive selection.

Influenza A viruses are single-stranded RNA viruses capable of evolving rapidly to adapt to environmental conditions. Examples include the establishment of a virus in a novel host or an adaptation to increasing immunity within the host population due to prior infection or vaccination against a circulating strain. Knowledge of the viral protein regions under positive selection is therefore crucial for surveillance. We have developed a method for detecting positively selected patches of sites on the surface of viral proteins, which we assume to be relevant for adaptive evolution. We measure positive selection based on dN/dS ratios of genetic changes inferred by considering the phylogenetic structure of the data and suggest a graph-cut algorithm to identify such regions. Our algorithm searches for dense and spatially distinct clusters of sites under positive selection on the protein surface. For the hemagglutinin protein of human influenza A viruses of the subtypes H3N2 and H1N1, our predicted sites significantly overlap with known antigenic and receptor-binding sites. From the structure and sequence data of the 2009 swine-origin influenza A/H1N1 hemagglutinin and PB2 protein, we identified regions that provide evidence of evolution under positive selection since introduction of the virus into the human population. The changes in PB2 overlap with sites reported to be associated with mammalian adaptation of the influenza A virus. Application of our technique to the protein structures of viruses of yet unknown adaptive behavior could identify further candidate regions that are important for host-virus interaction.

The human cytomegalovirus UL11 protein interacts with the receptor tyrosine phosphatase CD45, resulting in functional paralysis of T cells.

Human cytomegalovirus (CMV) exerts diverse and complex effects on the immune system, not all of which have been attributed to viral genes. Acute CMV infection results in transient restrictions in T cell proliferative ability, which can impair the control of the virus and increase the risk of secondary infections in patients with weakened or immature immune systems. In a search for new immunomodulatory proteins, we investigated the UL11 protein, a member of the CMV RL11 family. This protein family is defined by the RL11 domain, which has homology to immunoglobulin domains and adenoviral immunomodulatory proteins. We show that pUL11 is expressed on the cell surface and induces intercellular interactions with leukocytes. This was demonstrated to be due to the interaction of pUL11 with the receptor tyrosine phosphatase CD45, identified by mass spectrometry analysis of pUL11-associated proteins. CD45 expression is sufficient to mediate the interaction with pUL11 and is required for pUL11 binding to T cells, indicating that pUL11 is a specific CD45 ligand. CD45 has a pivotal function regulating T cell signaling thresholds; in its absence, the Src family kinase Lck is inactive and signaling through the T cell receptor (TCR) is therefore shut off. In the presence of pUL11, several CD45-mediated functions were inhibited. The induction of tyrosine phosphorylation of multiple signaling proteins upon TCR stimulation was reduced and T cell proliferation was impaired. We therefore conclude that pUL11 has immunosuppressive properties, and that disruption of T cell function via inhibition of CD45 is a previously unknown immunomodulatory strategy of CMV.

Inference of genotype-phenotype relationships in the antigenic evolution of human influenza A (H3N2) viruses.

Distinguishing mutations that determine an organism's phenotype from (near-) neutral 'hitchhikers' is a fundamental challenge in genome research, and is relevant for numerous medical and biotechnological applications. For human influenza viruses, recognizing changes in the antigenic phenotype and a strains' capability to evade pre-existing host immunity is important for the production of efficient vaccines. We have developed a method for inferring 'antigenic trees' for the major viral surface protein hemagglutinin. In the antigenic tree, antigenic weights are assigned to all tree branches, which allows us to resolve the antigenic impact of the associated amino acid changes. Our technique predicted antigenic distances with comparable accuracy to antigenic cartography. Additionally, it identified both known and novel sites, and amino acid changes with antigenic impact in the evolution of influenza A (H3N2) viruses from 1968 to 2003. The technique can also be applied for inference of 'phenotype trees' and genotype-phenotype relationships from other types of pairwise phenotype distances.