Predicting coding potential from genome sequence: application to betaherpesviruses infecting rats and mice

Decoding murine cytomegalovirus

The genomes of both human cytomegalovirus (HCMV) and murine cytomegalovirus (MCMV) were first sequenced over 20 years ago. Similar to HCMV, the MCMV genome had initially been proposed to harbor ≈170 open reading frames (ORFs). More recently, omics approaches revealed HCMV gene expression to be substantially more complex comprising several hundred viral ORFs. Here, we provide a state-of-the art reannotation of lytic MCMV gene expression based on integrative analysis of a large set of omics data. Our data reveal 365 viral transcription start sites (TiSS) that give rise to 380 and 454 viral transcripts and ORFs, respectively. The latter include >200 small ORFs, some of which represented the most highly expressed viral gene products. By combining TiSS profiling with metabolic RNA labelling and chemical nucleotide conversion sequencing (dSLAM-seq), we provide a detailed picture of the expression kinetics of viral transcription. This not only resulted in the identification of a novel MCMV immediate early transcript encoding the m166.5 ORF, which we termed ie4, but also revealed a group of well-expressed viral transcripts that are induced later than canonical true late genes and contain an initiator element (Inr) but no TATA- or TATT-box in their core promoters. We show that viral upstream ORFs (uORFs) tune gene expression of longer viral ORFs expressed in cis at translational level. Finally, we identify a truncated isoform of the viral NK-cell immune evasin m145 arising from a viral TiSS downstream of the canonical m145 mRNA. Despite being ≈5-fold more abundantly expressed than the canonical m145 protein it was not required for downregulating the NK cell ligand, MULT-I. In summary, our work will pave the way for future mechanistic studies on previously unknown cytomegalovirus gene products in an important virus animal model.

COVID-19: Attacks Immune Cells and Interferences With Antigen Presentation Through MHC-Like Decoy System

The high mortality of coronavirus disease 2019 is related to poor antigen presentation and lymphopenia. Cytomegalovirus and the herpes family encode a series of major histocompatibility complex (MHC)-like molecules required for targeted immune responses to achieve immune escape. In this present study, domain search results showed that many proteins of the severe acute respiratory syndrome coronavirus 2 virus had MHC-like domains, which were similar to decoys for the human immune system. MHC-like structures could bind to MHC receptors of immune cells (such as CD4 + T-cell, CD8 + T-cell, and natural killer-cell), interfering with antigen presentation. Then the oxygen free radicals generated by E protein destroyed immune cells after MHC-like of S protein could bind to them. Mutations in the MHC-like region of the viral proteins such as S promoted weaker immune resistance and more robust transmission. S 127-194 were the primary reason for the robust transmission of delta variants. The S 144-162 regulated the formation of S trimer. The mutations of RdRP: G671S and N: D63G of delta variant caused high viral load. S 62-80 of alpha, beta, lambda variants were the important factor for fast-spreading. S 616-676 and 1014-1114 were causes of high mortality for gamma variants infections. These sites were in the MHC-like structure regions.

Membraneless Compartmentalization of Nuclear Assembly Sites during Murine Cytomegalovirus Infection

Extensive reorganization of infected cells and the formation of large structures known as the nuclear replication compartment (RC) and cytoplasmic assembly compartment (AC) is a hallmark of beta-herpesvirus infection. These restructurings rely on extensive compartmentalization of the processes that make up the virus manufacturing chain. Compartmentalization of the nuclear processes during murine cytomegalovirus (MCMV) infection is not well described. In this study, we visualized five viral proteins (pIE1, pE1, pM25, pm48.2, and pM57) and replicated viral DNA to reveal the nuclear events during MCMV infection. As expected, these events can be matched with those described for other beta and alpha herpesviruses and contribute to the overall picture of herpesvirus assembly. Imaging showed that four viral proteins (pE1, pM25, pm48.2, and pM57) and replicated viral DNA condense in the nucleus into membraneless assemblies (MLAs) that undergo a maturation sequence to form the RC. One of these proteins (pM25), which is also expressed in a cytoplasmic form (pM25l), showed similar MLAs in the AC. Bioinformatics tools for predicting biomolecular condensates showed that four of the five proteins had a high propensity for liquid–liquid phase separation (LLPS), suggesting that LLPS may be a mechanism for compartmentalization within RC and AC. Examination of the physical properties of MLAs formed during the early phase of infection by 1,6-hexanediol treatment in vivo revealed liquid-like properties of pE1 MLAs and more solid-like properties of pM25 MLAs, indicating heterogeneity of mechanisms in the formation of virus-induced MLAs. Analysis of the five viral proteins and replicated viral DNA shows that the maturation sequence of RC and AC is not completed in many cells, suggesting that virus production and release is carried out by a rather limited number of cells. This study thus lays the groundwork for further investigation of the replication cycle of beta-herpesviruses, and the results should be incorporated into plans for high-throughput and single-cell analytic approaches.

The m15 Locus of Murine Cytomegalovirus Modulates Natural Killer Cell Responses to Promote Dissemination to the Salivary Glands and Viral Shedding

As the largest herpesviruses, the 230 kb genomes of cytomegaloviruses (CMVs) have increased our understanding of host immunity and viral escape mechanisms, although many of the annotated genes remain as yet uncharacterised. Here we identify the m15 locus of murine CMV (MCMV) as a viral modulator of natural killer (NK) cell immunity. We show that, rather than discrete transcripts from the m14, m15 and m16 genes as annotated, there are five 3′-coterminal transcripts expressed over this region, all utilising a consensus polyA tail at the end of the m16 gene. Functional inactivation of any one of these genes had no measurable impact on viral replication. However, disruption of all five transcripts led to significantly attenuated dissemination to, and replication in, the salivary glands of multiple strains of mice, but normal growth during acute infection. Disruption of the m15 locus was associated with heightened NK cell responses, including enhanced proliferation and IFNγ production. Depletion of NK cells, but not T cells, rescued salivary gland replication and viral shedding. These data demonstrate the identification of multiple transcripts expressed by a single locus which modulate, perhaps in a concerted fashion, the function of anti-viral NK cells.

Helixer: cross-species gene annotation of large eukaryotic genomes using deep learning

Motivation

Current state-of-the-art tools for the de novo annotation of genes in eukaryotic genomes have to be specifically fitted for each species and still often produce annotations that can be improved much further. The fundamental algorithmic architecture for these tools has remained largely unchanged for about two decades, limiting learning capabilities. Here, we set out to improve the cross-species annotation of genes from DNA sequence alone with the help of deep learning. The goal is to eliminate the dependency on a closely related gene model while also improving the predictive quality in general with a fundamentally new architecture.

Results

We present Helixer, a framework for the development and usage of a cross-species deep learning model that improves significantly on performance and generalizability when compared to more traditional methods. We evaluate our approach by building a single vertebrate model for the base-wise annotation of 186 animal genomes and a separate land plant model for 51 plant genomes. Our predictions are shown to be much less sensitive to the length of the genome than those of a current state-of-the-art tool. We also present two novel post-processing techniques that each worked to further strengthen our annotations and show in-depth results of an RNA-Seq based comparison of our predictions. Our method does not yet produce comprehensive gene models but rather outputs base pair wise probabilities.

Availability and implementation

The source code of this work is available at https://github.com/weberlab-hhu/Helixer under the GNU General Public License v3.0. The trained models are available at https://doi.org/10.5281/zenodo.3974409

Supplementary information

Supplementary data are available at Bioinformatics online.

Rapid adaptation to human protein kinase R by a unique genomic rearrangement in rhesus cytomegalovirus

Cytomegaloviruses (CMVs) are generally unable to cross species barriers, in part because prolonged coevolution with one host species limits their ability to evade restriction factors in other species. However, the limitation in host range is incomplete. For example, rhesus CMV (RhCMV) can replicate in human cells, albeit much less efficiently than in rhesus cells. Previously we reported that the protein kinase R (PKR) antagonist encoded by RhCMV, rTRS1, has limited activity against human PKR but is nonetheless necessary and sufficient to enable RhCMV replication in human fibroblasts (HF). We now show that knockout of PKR in human cells or treatment with the eIF2B agonist ISRIB, which overcomes the translational inhibition resulting from PKR activation, augments RhCMV replication in HF, indicating that human PKR contributes to the inefficiency of RhCMV replication in HF. Serial passage of RhCMV in HF reproducibly selected for viruses with improved ability to replicate in human cells. The evolved viruses contain an inverted duplication of the terminal 6.8 kb of the genome, including rTRS1. The duplication replaces ~11.8 kb just downstream of an internal sequence element, pac1-like, which is very similar to the pac1 cleavage and packaging signal found near the terminus of the genome. Plaque-purified evolved viruses produced at least twice as much rTRS1 as the parental RhCMV and blocked the PKR pathway more effectively in HF. Southern blots revealed that unlike the parental RhCMV, viruses with the inverted duplication isomerize in a manner similar to HCMV and other herpesviruses that have internal repeat sequences. The apparent ease with which this duplication event occurs raises the possibility that the pac1-like site, which is conserved in Old World monkey CMV genomes, may serve a function in facilitating rapid adaptation to evolutionary obstacles.

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.

Mouse Cytomegalovirus M34 Encodes a Non-essential, Nuclear, Early-Late Expressed Protein Required for Efficient Viral Replication

Human cytomegalovirus (HCMV) is a prototypic betaherpesvirus which causes severe manifestations in individuals with impaired or immature immunity. To investigate cytomegalovirus-induced pathogenesis and virus-specific immune responses, mouse cytomegalovirus (MCMV) infections in mice are employed as accepted small animal model. MCMV and HCMV share co-linear genomes and encode several homologous proteins. Due to the size and complexity of CMV genomes, the molecular functions of numerous cytomegaloviral gene products remain to be elucidated. While the essential nature of viral genes highlights their biological relevance, it renders functional studies particularly cumbersome by precluding experiments in the infection context. The HCMV-encoded protein pUL34 binds the HCMV genome and regulates viral gene expression (e.g., of US3). Several groups provided compelling evidence that UL34 is essential for HCMV replication. MCMV encodes the homologous protein pM34 (34% identical and 55% similar). Based on unsuccessful attempts to reconstitute M34-deficient virus from a bacterial artificial chromosome (BAC), M34 was previously classified as essential for MCMV replication. To characterize pM34 during viral infection, we engineered and analyzed an MCMV mutant expressing an HA-epitope-tagged pM34 which was expressed with early-late kinetics and localized in the nucleus. Additionally, we generated an M34-deficient (“ΔM34”) MCMV-BAC by replacing the entire M34 coding sequence by a kanamycin resistance cassette. The deletion of M34 was confirmed by Southern blot and PCR. Unexpectedly, we could reconstitute replicating ΔM34-MCMV upon transfection of the BAC DNA into mouse embryonic fibroblasts. The absence of M34 from the genome of the replicating ΔM34-MCMV was also confirmed. Accordingly, a ΔM34-MCMV, in which the kanamycin cassette was excised by frt/Flp-mediated recombination, was also replication competent. In order to corroborate the absence of pM34 protein, the M34 deletion was recapitulated on the background of M34HA, which yielded replicating virus devoid of detectable pM34HA protein. The replication of MCMVs lacking M34 was found to be 10- to 100-fold reduced as compared to wt-MCMV which might explain previous unsuccessful reconstitution attempts conducted by others. Taken together, our findings reveal that MCMV remains replication competent despite the absence of M34, enabling functional studies in the infection context.

Ribosome signatures aid bacterial translation initiation site identification

BACKGROUND

While methods for annotation of genes are increasingly reliable, the exact identification of translation initiation sites remains a challenging problem. Since the N-termini of proteins often contain regulatory and targeting information, developing a robust method for start site identification is crucial. Ribosome profiling reads show distinct patterns of read length distributions around translation initiation sites. These patterns are typically lost in standard ribosome profiling analysis pipelines, when reads from footprints are adjusted to determine the specific codon being translated.

RESULTS

Utilising these signatures in combination with nucleotide sequence information, we build a model capable of predicting translation initiation sites and demonstrate its high accuracy using N-terminal proteomics. Applying this to prokaryotic translatomes, we re-annotate translation initiation sites and provide evidence of N-terminal truncations and extensions of previously annotated coding sequences. These re-annotations are supported by the presence of structural and sequence-based features next to N-terminal peptide evidence. Finally, our model identifies 61 novel genes previously undiscovered in the Salmonella enterica genome.

CONCLUSIONS

Signatures within ribosome profiling read length distributions can be used in combination with nucleotide sequence information to provide accurate genome-wide identification of translation initiation sites.

Complete Genome Sequence of Elephant Endotheliotropic Herpesvirus 4, the First Example of a GC-Rich Branch Proboscivirus

A novel group of mammalian DNA viruses called elephant endotheliotropic herpesviruses (EEHVs) belonging to the Proboscivirus genus has been associated with nearly 100 cases of highly lethal acute hemorrhagic disease in young Asian elephants worldwide. The complete 180-kb genomes of prototype strains from three AT-rich branch viruses, EEHV1A, EEHV1B, and EEHV5, have been published. However, less than 6 kb of DNA sequence each from EEHV3, EEHV4, and EEHV7 showed them to be a hugely diverged second major branch with GC-rich characteristics. Here, we determined the complete 206-kb genome of EEHV4(Baylor) directly from trunk wash DNA by next-generation sequencing and de novo assembly procedures. Among a total of 119 genes with an overall colinear organization similar to those of the AT-rich EEHVs, major features of EEHV4 include a family of 26 paralogous 7xTM and vGPCR-like genes plus 25 novel or missing genes. The genome also contains an unusual distribution of tracts of 5 to 11 successive A or T nucleotides in intergenic domains between the mostly much higher GC content protein coding regions. Furthermore, an extremely high GC-rich bias in the third wobble position of codons clearly delineates the coding regions for many but not all proteins. There are also two novel captured cellular genes, including a C-type lectin (vECTL) and an O-linked acetylglucosamine transferase (vOGT), as well as an unusually large and complex Ori-Lyt dyad symmetry domain. Finally, 30 kb from a second strain proved to include three small chimeric domains, indicating the existence of distinct EEHV4A and EEHV4B subtypes. IMPORTANCE Multiple species of herpesviruses from three different lineages of the Proboscivirus genus (EEHV1/6, EEHV2/5, and EEHV3/4/7) infect both Asian and African elephants, but lethal hemorrhagic disease is largely confined to Asian elephant calves and is predominantly associated with EEHV1. Milder disease caused by EEHV5 or EEHV4 is being increasingly recognized as well, but little is known about the latter, which is estimated to have diverged at least 35 million years ago from the others within a distinctive GC-rich branch of the Proboscivirus genus. Here, we have determined the complete genomic DNA sequence of a strain of EEHV4 obtained from a trunk wash sample collected from a surviving Asian elephant calf undergoing asymptomatic shedding during convalescence after an acute hemorrhagic disease episode. This represents the first example from among the three known GC-rich branch Proboscivirus species to be assembled and fully annotated. Several distinctive features of EEHV4 compared to AT-rich branch genomes are described.

Comparison of the Gene Coding Contents and Other Unusual Features of the GC-Rich and AT-Rich Branch Probosciviruses

Multiple species of herpesviruses from three different lineages of the Proboscivirus genus (EEHV1/6, EEHV2/5, and EEHV3/4/7) infect either Asian or African elephants, but the highly lethal hemorrhagic disease is largely confined to Asian elephant calves and is predominantly associated with EEHV1. In the accompanying paper [P. D. Ling et al., mSphere 1(3):e00081-15, http://dx.doi.org/10.1128/mSphere.00081-15], we report the complete 206-kb genome of EEHV4, the third different species causing disease in Asian elephants and the first example of a GC-rich branch proboscivirus. To gain insights into the nature and differential properties of these two very anciently diverged lineages of elephant herpesviruses, we describe here several additional unusual features found in the complete GC-rich genome of EEHV4 with particular emphasis on patterns of divergence as well as common unique features that are distinct from those of all other herpesviruses, such as the enlarged AT-rich intergenic domains and gene families, including the large number of vGPCR-like proteins.

Nearly 100 cases of lethal acute hemorrhagic disease in young Asian elephants have been reported worldwide. All tested cases contained high levels of elephant endotheliotropic herpesvirus (EEHV) DNA in pathological blood or tissue samples. Seven known major types of EEHVs have been partially characterized and shown to all belong to the novel Proboscivirus genus. However, the recently determined 206-kb EEHV4 genome proved to represent the prototype of a GC-rich branch virus that is very distinct from the previously published 180-kb EEHV1A, EEHV1B, and EEHV5A genomes, which all fall within an alternative AT-rich branch. Although EEHV4 retains the large family of 7xTM and vGPCR-like genes, six are unique to either just one or the other branch. While both branches display a highly enriched distribution of A and T tracts in intergenic domains, they are generally much larger within the GC-rich branch. Both branches retain the vGCNT1 acetylglucosamine transferase and at least one vOX-2 gene, but the two branches differ by 25 genes overall, with the AT-rich branch encoding a fucosyl transferase (vFUT9) plus two or three more vOX2 proteins and an immunoglobulin-like gene family that are all absent from the GC-rich branch. Several envelope glycoproteins retain only 15 to 20% protein identity or less across the two branches. Finally, the two plausible predicted transcriptional regulatory proteins display no homology at all to those in the alpha-, beta-, or gammaherpesvirus subfamilies. These results reinforce our previous proposal that the probosciviruses should be designated a new subfamily of mammalian herpesviruses.

IMPORTANCE Multiple species of herpesviruses from three different lineages of the Proboscivirus genus (EEHV1/6, EEHV2/5, and EEHV3/4/7) infect either Asian or African elephants, but the highly lethal hemorrhagic disease is largely confined to Asian elephant calves and is predominantly associated with EEHV1. In the accompanying paper [P. D. Ling et al., mSphere 1(3):e00081-15, 10.1128/mSphere.00081-15], we report the complete 206-kb genome of EEHV4, the third different species causing disease in Asian elephants and the first example of a GC-rich branch proboscivirus. To gain insights into the nature and differential properties of these two very anciently diverged lineages of elephant herpesviruses, we describe here several additional unusual features found in the complete GC-rich genome of EEHV4 with particular emphasis on patterns of divergence as well as common unique features that are distinct from those of all other herpesviruses, such as the enlarged AT-rich intergenic domains and gene families, including the large number of vGPCR-like proteins.

Rat cytomegalovirus (RCMV) English isolate and a newly identified Berlin isolate share similarities with but are separate as an anciently diverged clade from Mouse CMV and the Maastricht isolate of RCMV

The genome of the rat cytomegalovirus (RCMV) English isolate (MuHV-8) differs significantly from the RCMV Maastricht isolate (MuHV-2) and other cytomegaloviruses (CMVs) in its size, base composition and genomic content. Analysis of the RCMV-Berlin isolate, MuHV-8, revealed that the two MuHV-8 isolates are highly similar in genome size and content, indicating that the smaller genome size (202 946 bp) compared to other known CMVs was not the result of an accidental deletion during passage in tissue culture. Surprisingly, the proteins encoded in MuHV-8 shared more overall similarity with their orthologues from mouse CMV (MuHV-1) compared to their orthologues in rat CMV (MuHV-2). Phylogenetic analyses of conserved viral genes showed that the two MuHV-8 isolates are from the same species and represent a unique clade that is distinct from other rodent CMVs.

Quantitative frame analysis and the annotation of GC-rich (and other) prokaryotic genomes. An application to Anaeromyxobacter dehalogenans

Motivation: Graphical representations of contrasts in GC usage among codon frame positions (frame analysis) provide evidence of genes missing from the annotations of prokaryotic genomes of high GC content but the qualitative approach of visual frame analysis prevents its applicability on a genomic scale.

Results: We developed two quantitative methods for the identification and statistical characterization in sequence regions of three-base periodicity (hits) associated with open reading frame structures. The methods were implemented in the N-Profile Analysis Computational Tool (NPACT), which highlights in graphical representations inconsistencies between newly identified ORFs and pre-existing annotations of coding-regions. We applied the NPACT procedures to two recently annotated strains of the deltaproteobacterium Anaeromyxobacter dehalogenans, identifying in both genomes numerous conserved ORFs not included in the published annotation of coding regions.

Availability and implementation: NPACT is available as a web-based service and for download at http://genome.ufl.edu/npact.

Contact:lucianob@ufl.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Tuberous Sclerosis Complex Protein 2-Independent Activation of mTORC1 by Human Cytomegalovirus pUL38

The mammalian target of rapamycin complex 1 (mTORC1) controls cell growth and anabolic metabolism and is a critical host factor activated by human cytomegalovirus (HCMV) for successful infection. The multifunctional HCMV protein pUL38 previously has been reported to activate mTORC1 by binding to and antagonizing tuberous sclerosis complex protein 2 (TSC2) (J. N. Moorman et al., Cell Host Microbe 3:253-262, 2008, http://dx.doi.org/10.1016/j.chom.2008.03.002). pUL38 also plays a role in blocking endoplasmic reticulum stress-induced cell death during HCMV infection. In this study, we showed that a mutant pUL38 lacking the N-terminal 24 amino acids (pHA-UL3825-331) was fully functional in suppressing cell death during infection. Interestingly, pHA-UL3825-331 lost the ability to interact with TSC2 but retained the ability to activate mTORC1, although to a lesser extent than full-length pHA-UL38. Recombinant virus expressing pHA-UL3825-331 replicated with ∼10-fold less efficiency than the wild-type virus at a low multiplicity of infection (MOI), but it grew similarly well at a high MOI, suggesting an MOI-dependent importance of pUL38-TSC2 interaction in supporting virus propagation. Site-directed mutational analysis identified a TQ motif at amino acid residues 23 and 24 as critical for pUL38 interaction with TSC2. Importantly, when expressed in isolation, the TQ/AA substitution mutant pHA-UL38 TQ/AA was capable of activating mTORC1 just like pHA-UL3825-331. We also created TSC2-null U373-MG cell lines by CRISPR genome editing and showed that pUL38 was capable of further increasing mTORC1 activity in TSC2-null cells. Therefore, this study identified the residues important for pUL38-TSC2 interaction and demonstrated that pUL38 can activate mTORC1 in both TSC2-dependent and -independent manners.

IMPORTANCE

HCMV, like other viruses, depends exclusively on its host cell to propagate. Therefore, it has developed methods to protect against host stress responses and to usurp cellular processes to complete its life cycle. mTORC1 is believed to be important for virus replication, and HCMV maintains high mTORC1 activity despite the stressful cellular environment associated with infection. mTORC1 inhibitors suppressed HCMV replication in vitro and reduced the incidence of HCMV reactivation in transplant recipients. We demonstrated that mTORC1 was activated by HCMV protein pUL38 in both TSC2-dependent and TSC2-independent manners. The pUL38-independent mode of mTORC1 activation also has been reported. These novel findings suggest the evolution of sophisticated approaches whereby HCMV activates mTORC1, indicating its importance in the biology and pathogenesis of HCMV.

Complexity and Diversity of the NKR-P1:Clr (Klrb1:Clec2) Recognition Systems

The NKR-P1 receptors were identified as prototypical natural killer (NK) cell surface antigens and later shown to be conserved from rodents to humans on NK cells and subsets of T cells. C-type lectin-like in nature, they were originally shown to be capable of activating NK cell function and to recognize ligands on tumor cells. However, certain family members have subsequently been shown to be capable of inhibiting NK cell activity, and to recognize proteins encoded by a family of genetically linked C-type lectin-related ligands. Some of these ligands are expressed by normal, healthy cells, and modulated during transformation, infection, and cellular stress, while other ligands are upregulated during the immune response and during pathological circumstances. Here, we discuss historical and recent developments in NKR-P1 biology that demonstrate this NK receptor–ligand system to be far more complex and diverse than originally anticipated.

Transcription of true late (γ2) cytomegalovirus genes requires UL92 function that is conserved among beta- and gammaherpesviruses

Human cytomegalovirus-encoded UL92 plays an essential role in viral replication that has not been resolved. We show here that this gene controls the accumulation of true late (γ2) viral transcripts, a property shared with several other recently evaluated genes (UL79, UL87, UL91, and UL95) conserved among beta- and gammaherpesviruses. When the UL92 mutant virus was evaluated, function was fully complemented by either the natural protein or the homologous Rh127 protein from rhesus cytomegalovirus. N-terminal epitope-tagged UL92 protein is functional, follows complex early-late expression kinetics, and localizes in the nucleus within viral replication compartments. UL92 severely impacts the late (72-h postinfection) expression of nine genes encoding virion proteins (UL32, UL55, UL73, UL75, UL80, UL86, UL99, and UL115), as well as UL91 and itself, but does not influence the levels of UL44, UL82, or UL83 accumulation. Although viral DNA is made at normal levels, viral capsid accumulation in the nucleus is severely compromised in UL92 mutant virus-infected cells, and mature virions are not observed in the cytoplasm. Taken together, UL92 is a key regulator of late viral gene expression, apparently functioning with four other beta- or gammaherpesvirus gene products in a pattern that appears reminiscent of gene regulation in T4 DNA bacteriophage.

Dual analysis of the murine cytomegalovirus and host cell transcriptomes reveal new aspects of the virus-host cell interface

Major gaps in our knowledge of pathogen genes and how these gene products interact with host gene products to cause disease represent a major obstacle to progress in vaccine and antiviral drug development for the herpesviruses. To begin to bridge these gaps, we conducted a dual analysis of Murine Cytomegalovirus (MCMV) and host cell transcriptomes during lytic infection. We analyzed the MCMV transcriptome during lytic infection using both classical cDNA cloning and sequencing of viral transcripts and next generation sequencing of transcripts (RNA-Seq). We also investigated the host transcriptome using RNA-Seq combined with differential gene expression analysis, biological pathway analysis, and gene ontology analysis. We identify numerous novel spliced and unspliced transcripts of MCMV. Unexpectedly, the most abundantly transcribed viral genes are of unknown function. We found that the most abundant viral transcript, recently identified as a noncoding RNA regulating cellular microRNAs, also codes for a novel protein. To our knowledge, this is the first viral transcript that functions both as a noncoding RNA and an mRNA. We also report that lytic infection elicits a profound cellular response in fibroblasts. Highly upregulated and induced host genes included those involved in inflammation and immunity, but also many unexpected transcription factors and host genes related to development and differentiation. Many top downregulated and repressed genes are associated with functions whose roles in infection are obscure, including host long intergenic noncoding RNAs, antisense RNAs or small nucleolar RNAs. Correspondingly, many differentially expressed genes cluster in biological pathways that may shed new light on cytomegalovirus pathogenesis. Together, these findings provide new insights into the molecular warfare at the virus-host interface and suggest new areas of research to advance the understanding and treatment of cytomegalovirus-associated diseases.

Cytomegalovirus UL91 is essential for transcription of viral true late (γ2) genes

Human cytomegalovirus-encoded UL91 is a betagamma gene that is essential for viral replication. Here we show that the 111-amino-acid (aa) UL91 protein controls accumulation of true-late (γ2) viral transcripts. The primate betaherpesvirus conserved N-terminal region from aa 1 to 71 is sufficient to fully reconstitute function. Evaluation of viral DNA, RNA, and antigen revealed that UL91 protein is expressed with leaky-late (γ1) kinetics, localizes in the nucleus without influencing viral DNA synthesis, and must be present from 48 h postinfection to support full expression of late viral transcripts and proteins. In the absence of UL91, viral capsid assembly in the nucleus of infected cells is significantly reduced, and mature, cytoplasmic virions fail to form. Taken together, the evidence shows that UL91 regulates late viral gene expression by a mechanism that is apparently conserved in betaherpesviruses and gammaherpesviruses.

Gene products of the embedded m41/m41.1 locus of murine cytomegalovirus differentially influence replication and pathogenesis

Cytomegaloviruses utilize overlapping and embedded reading frames as a way to efficiently package and express all genes necessary to carry out a complex lifecycle. Murine cytomegalovirus encodes a mitochondrial-localized inhibitor of Bak oligomerization (vIBO) from m41.1, a reading frame that is embedded within the m41 gene. The m41.1-encoded mitochondrial protein and m41-encoded Golgi-localized protein have both been implicated in cell death suppression; however, their contribution to viral infection within the host has not been investigated. Here, we report that mitochondrial-localized m41.1 (vIBO) is required for optimal viral replication in macrophages and has a modest impact on dissemination in infected mice. In contrast, Golgi-localized m41 protein is dispensable during acute infection and dissemination as well as for latency. All together, these data indicate that the primary evolutionary focus of this locus is to maintain mitochondrial function through inhibition of Bak-mediated death pathways in support of viral pathogenesis.

The genome of murine cytomegalovirus is shaped by purifying selection and extensive recombination

The herpesvirus lifestyle results in a long-term interaction between host and invading pathogen, resulting in exquisite adaptation of virus to host. We have sequenced the genomes of nine strains of murine cytomegalovirus (a betaherpesvirus), isolated from free-living mice trapped at locations separated geographically and temporally. Despite this separation these genomes were found to have low levels of nucleotide variation. Of the more than 160 open reading frames, almost 90% had a dN/dS ratio of amino acid substitutions of less than 0.6, indicating the level of purifying selection on the coding potential of MCMV. Examination of selection acting on individual genes at the codon level however indicates some level of positive selection, with 0.03% of codons showing strong evidence for positive selection. Conversely, 1.3% of codons show strong evidence of purifying selection. Alignments of both genome sequences and coding regions suggested that high levels of recombination have shaped the MCMV genome.

How the virus outsmarts the host: function and structure of cytomegalovirus MHC-I-like molecules in the evasion of natural killer cell surveillance

Natural killer (NK) cells provide an initial host immune response to infection by many viral pathogens. Consequently, the viruses have evolved mechanisms to attenuate the host response, leading to improved viral fitness. One mechanism employed by members of the β-herpesvirus family, which includes the cytomegaloviruses, is to modulate the expression of cell surface ligands recognized by NK cell activation molecules. A novel set of cytomegalovirus (CMV) genes, exemplified by the mouse m145 family, encode molecules that have structural and functional features similar to those of host major histocompatibility-encoded (MHC) class I molecules, some of which are known to contribute to immune evasion. In this review, we explore the function, structure, and evolution of MHC-I-like molecules of the CMVs and speculate on the dynamic development of novel immunoevasive functions based on the MHC-I protein fold.

Temporal profiling of the coding and noncoding murine cytomegalovirus transcriptomes

The global transcriptional program of murine cytomegalovirus (MCMV), involving coding, noncoding, and antisense transcription, remains unknown. Here we report an oligonucleotide custom microarray platform capable of measuring both coding and noncoding transcription on a genome-wide scale. By profiling MCMV wild-type and immediate-early mutant strains in fibroblasts, we found rapid activation of the transcriptome by 6.5 h postinfection, with absolute dependency on ie3, but not ie1 or ie2, for genomic programming of viral gene expression. Evidence is also presented to show, for the first time, genome-wide noncoding and bidirectional transcription at late stages of MCMV infection.

Re-evaluation of the genome sequence of guinea pig cytomegalovirus

Congenital infection by human cytomegalovirus (HCMV) is a major cause of birth defects and developmental abnormalities. Since guinea pig cytomegalovirus (GPCMV) crosses the placenta and causes infection in utero, GPCMV models are useful for studies of the mechanisms of transplacental transmission. During our characterization of a genomic locus required for GPCMV dissemination in animals, we found that the nucleotide sequence in and around the nearby immediate-early genes in our lineage of GPCMV strain 22122 [designated GPCMV (ATCC-P5)] showed clear differences from that reported previously for the same strain [designated GPCMV (UMN)] passaged extensively in vitro. Since in vitro passaging of HCMV is known to result in genetic alterations, especially in the UL128-UL131A locus, and loss of growth ability in particular cell types, in this study we determined the complete genome sequence of GPCMV (ATCC-P5), which grows efficiently in animals. A total of 359 differences were identified between the genome sequences of GPCMV (UMN) and GPCMV (ATCC-P5), and these resulted in structural differences in 29 protein-encoding regions. In addition, some genes predicted from our analysis but not from GPCMV (UMN) are well conserved among cytomegaloviruses. An additional 18 passages of GPCMV (ATCC-P5) in vitro generated no further marked alterations in these genes or in the locus corresponding to the HCMV UL128-UL131A. Our analyses indicate that the published sequence of GPCMV (UMN) contains a substantial number of sequencing errors and, possibly, some mutations resulting from a long history of passaging in vitro. Our re-evaluation of the genetic content of GPCMV will provide a solid foundation for future studies.

B7-mediated costimulation of CD4 T cells constrains cytomegalovirus persistence

Cytomegalovirus (CMV) utilizes multiple strategies to modulate immunity and promote lifelong, persistent/latent infection, including suppressing T cell activation pathways. Here we examined the role of B7 costimulatory ligands in establishing immune détente from both the host and virus perspectives. Mice lacking both B7.1 and B7.2 showed reduced early expansion of CMV-specific CD4 T cells, consequently allowing for enhanced levels of persistent virus replication. In turn, a CMV mutant lacking expression of the m138 and m147.5 gene products, which restrict B7.1 and B7.2 expression in infected antigen-presenting cells, induced a more robust CD4 T cell response and showed decreased persistence. Together, these data reveal a requirement for B7-mediated signaling in regulating the CMV-specific CD4 T cell response and establishing host-virus equilibrium.

Modulation of host innate and adaptive immune defenses by cytomegalovirus: timing is everything

Human cytomegalovirus (HCMV) (HHV-5, a beta-herpesvirus) causes the vast majority of infection-related congenital birth defects, and can trigger severe disease in immune suppressed individuals. The high prevalence of societal infection, the establishment of lifelong persistence and the growing number of immune-related diseases where HCMV is touted as a potential promoter is slowly heightening public awareness to this virus. The millions of years of co-evolution between CMV and the immune system of its host provides for a unique opportunity to study immune defense strategies, and pathogen counterstrategies. Dissecting the timing of the cellular and molecular processes that regulate innate and adaptive immunity to this persistent virus has revealed a complex defense network that is shaped by CMV immune modulation, resulting in a finely tuned host-pathogen relationship.

Multifaceted evasion of the interferon response by cytomegalovirus

Human cytomegalovirus (HCMV), which infects the majority of the population worldwide, causes few, if any, symptoms in otherwise healthy people but is responsible for considerable morbidity and mortality in immunocompromised patients and in congenitally infected newborns. The evolutionary success of HCMV depends in part on its ability to evade host defense systems. Here we review recent progress in elucidating the remarkable assortment of mechanisms employed by HCMV and the related beta-herpesviruses, murine cytomegaloviruses (MCMV) and rhesus cytomegaloviruses (RhCMV), for counteracting the host interferon (IFN) response. Very early after infection, cellular membrane sensors such as the lymphotoxin beta receptor initiate the production of antiviral cytokines including type I IFNs. However, virion factors, such as pp65 (ppUL83) and viral proteins made soon after infection including the immediate early gene 2 protein (pUL122), repress this response by interfering with steps in the activation of IFN regulatory factor 3 and NF-kappaB. CMVs then exert a multi-pronged attack on downstream IFN signaling. HCMV infection results in decreased accumulation and phosphorylation of the IFN signaling kinases Jak1 and Stat2, and the MCMV protein pM27 mediates Stat2 down-regulation, blocking both type I and type II IFN signaling. The HCMV immediate early gene 1 protein (pUL123) interacts with Stat2 and inhibits transcriptional activation of IFN-regulated genes. Infection also causes reduction in the abundance of p48/IRF9, a component of the ISGF3 transcription factor complex. Furthermore, CMVs have multiple genes involved in blocking the function of IFN-induced effectors. For example, viral double-stranded RNA-binding proteins are required to prevent the shutoff of protein synthesis by protein kinase R, further demonstrating the vital importance of evading the IFN response at multiple levels during infection.

The human cytomegalovirus UL36 gene controls caspase-dependent and -independent cell death programs activated by infection of monocytes differentiating to macrophages

The cellular protease caspase-8 activates extrinsic apoptosis and also functions to promote monocyte-to-macrophage differentiation. Differentiation-induced alterations to antiviral caspase-8-dependent cell death pathways are unclear. Here, we show THP-1 monocyte-to-macrophage differentiation alters the specific cell death pathways activated in response to human cytomegalovirus (HCMV) infection. Employing viruses with mutations in UL36, the gene that encodes the viral inhibitor of caspase-8 activation (vICA), our data indicate that both caspase-dependent and -independent death pathways are activated in response to infection. Activation of caspase-dependent and -independent cell death responses restricted growth of vICA-deficient viruses, and vICA/pUL36 inhibited either response. Thus, these studies also reveal that the UL36 gene controls a caspase-independent cell death pathway. The impact of caspases on control of antiviral responses differed at early and late stages of macrophage differentiation. Early in differentiation, vICA-deficient virus-induced cell death was dependent on caspases and inhibited by the pan-caspase inhibitor z-VAD(OMe)-fluoromethyl ketone. In contrast, virus-induced death at late times of differentiation was caspase independent. Additional unlabeled and fluorescent inhibitors indicated that caspase-8 promoted death from within infected cells at early but not late stages of differentiation. These data highlight the multifunctional role of vICA/pUL36 as HCMV encounters various antiviral responses during macrophage differentiation.

Stability of murine cytomegalovirus genome after in vitro and in vivo passage

While large DNA viruses are thought to have low mutation rates, only a small fraction of their genomes have been analyzed at the single-nucleotide level. Here, we defined the genetic stability of murine cytomegalovirus (MCMV) by whole-genome sequencing. Independently assembled sequences of three sister plaques showed only two single-base-pair substitutions after in vitro passage. In vivo-passaged MCMV likewise demonstrated low mutation rates, comparable to those after in vitro passage, indicating high genome stability of MCMV at the single-nucleotide level in the absence of obvious selection pressure.

Cytomegaloviruses inhibit Bak- and Bax-mediated apoptosis with two separate viral proteins

Apoptosis of infected cells can limit virus replication and serves as an innate defense mechanism against viral infections. Consequently, viruses delay apoptosis by expressing antiapoptotic proteins, many of which structurally resemble the cellular antiapoptotic protein Bcl-2. Like Bcl-2, the viral analogs inhibit apoptosis by preventing activation and/or oligomerization of the proapoptotic mitochondrial proteins Bax and Bak. Here we show that cytomegaloviruses (CMVs) have adopted a different strategy. They encode two separate mitochondrial proteins that lack obvious sequence similarities to Bcl-2-family proteins and specifically counteract either Bax or Bak. We identified a small mitochondrion-localized protein encoded by the murine CMV open reading frame (ORF) m41.1, which functions as a viral inhibitor of Bak oligomerization (vIBO). It blocks Bak-mediated cytochrome c release and Bak-dependent induction of apoptosis. It protects cells from cell death-inducing stimuli together with the previously identified Bax-specific inhibitor viral mitochondria-localized inhibitor of apoptosis (vMIA) (encoded by ORF m38.5). Similar vIBO proteins are encoded by CMVs of rats, and possibly by other CMVs as well. These results suggest a non-redundant function of Bax and Bak during viral infection, and a benefit for CMVs derived from the ability to inhibit Bak and Bax separately with two viral proteins.

Genome sequence divergences and functional variations in human cytomegalovirus strains

Genome sequences of numerous and wide-ranging species have been completed, but genome-wide sequence variation patterns linked to biological functions are just starting to be investigated. Here, by comparatively analyzing the genome variation patterns of human cytomegalovirus (HCMV) genomes, we revealed large sequence divergences and functional variations existing in HCMV genomes. They are divergent in genome-size, inversion, orientation and coding potential, even within conserved genes, including nucleotide polymorphism, DNA strand composition asymmetry, and evolutionary rate variation in conserved genes. These divergences in conserved genes are linked to HCMV biology. Codon usage variation of conserved genes located in the negative DNA strand is significantly different between HCMV strains, and this variation associates with virion production and virulence factor, suggesting that the negative DNA strand primarily contributes to virion production and virulence factor in HCMV. In addition, we also revealed that genes functioning for entry and egress are the most adaptable, and that those for transcription and replication are the most conserved in HCMV genomes. The conserved-transcription system is generally controlled by a genome-wide motif GCGC revealed in this study by Chaos map analysis. Our findings demonstrated that genome sequences of HCMV are generally divergent and these divergences directly reflect viral biology.

Analysis of the nucleotide sequence of the guinea pig cytomegalovirus (GPCMV) genome

In this report we describe the genomic sequence of guinea pig cytomegalovirus (GPCMV) assembled from a tissue culture-derived bacterial artificial chromosome clone, plasmid clones of viral restriction fragments, and direct PCR sequencing of viral DNA. The GPCMV genome is 232,678 bp, excluding the terminal repeats, and has a GC content of 55%. A total of 105 open reading frames (ORFs) of > 100 amino acids with sequence and/or positional homology to other CMV ORFs were annotated. Positional and sequence homologs of human cytomegalovirus open reading frames UL23 through UL122 were identified. Homology with other cytomegaloviruses was most prominent in the central ~60% of the genome, with divergence of sequence and lack of conserved homologs at the respective genomic termini. Of interest, the GPCMV genome was found in many cases to bear stronger phylogenetic similarity to primate CMVs than to rodent CMVs. The sequence of GPCMV should facilitate vaccine and pathogenesis studies in this model of congenital CMV infection.

Evolution of the G+C Content Frontier in the Rat Cytomegalovirus Genome

Within the 230138 bp of the rat cytomegalovirus (RCMV) genome, the G+C content changes abruptly at position 142644, constituting a G+C content frontier. To the left of this point, overall G+C content is 69.2%, and to the right it is only 47.6%. A region of extremely low G+C content (33.8%) is found in the 5 kb immediately to the right of the frontier, in which there are no predicted coding sequences. To the right of position 147501, the G+C content rises and predicted coding sequences reappear. However, these genes are much shorter (average 848 bp, 50% G+C) than those in the left two-thirds of the genome (average 1462 bp, 70% G+C). Whole genome alignment of several viruses indicates that the initial ultra-low G+C region appeared in the common ancestor of the genera Cytomegalovirus and Muromegalovirus, and that the lowering of G+C in the right third has been a subsequent process in the lineage leading to RCMV. The left two-thirds of RCMV has stop codon occurrences at 67.5% of their expected level, based on a modified Markov chain model of stop codon distribution, and the corresponding figure for the right third is 78%. Therefore, despite heavy mutation pressure, selective constraint has operated in the right third of the RCMV genome to maintain a degree of gene length unusual for such low G+C sequences.

Control of apoptosis by human cytomegalovirus

Caspase-dependent apoptosis has an important role in controlling viruses, and as a result, viruses often encode proteins that target this pathway. Caspase-dependent apoptosis can be activated from within the infected cell as an intrinsic response to replication-associated stresses or through death-inducing signals produced extrinsically by immune cells. Cytomegaloviruses (CMVs) encode a mitochondria-localized inhibitor of apoptosis, vMIA, and a viral inhibitor of caspase activation, vICA, the functional homologs of Bcl-2 related and c-FLIP proteins, respectively. Evidence from viral mutants deleting either vMIA or vICA suggests that each is necessary and sufficient to promote survival of infected cells undergoing caspase-dependent apoptosis. Additional proteins, including pUL38, IE1(491a), and IE2(579aa), can prevent apoptosis induced by various stimuli, while viruses with deletions of UL38, M45, or m41 undergo apoptosis. The viral RNA, beta2.7, binds mitochondrial respiratory complex I, maintains ATP production late in infection, and prevents death induced by a mitochondrial poison. Thus, CMV alters cell intrinsic defenses employing apoptosis, and multiple viral gene products together control death-inducing stimuli to promote survival.

Preliminary characterization of murine cytomegaloviruses with insertional and deletional mutations in the M34 open reading frame

A viable virus could not be recovered from a mutant murine cytomegalovirus (MCMV) BAC in which the M34 ORF had been deleted (BACDeltaM34). In contrast, an M34 mutant virus (RcM34), in which the M34 ORF was interrupted by transposon insertion at nt 44,827 of the Smith MCMV BAC, was attenuated in replication both in tissue culture and in SCID mice. Similarly, mutant virus Rc3'DeltaM34, in which the 3'-end was deleted from nt 44,724 to nt 45,647, produced similar replication kinetics in tissue culture to RcM34 while BAC5'DeltaM34, in which the 5'-end from nt 43,083 up to nt 44,896 was deleted, was non-viable like BACDeltaM34. A transcript analysis of wt and RcM34 virus-infected cells showed that a truncated transcript encoding a putative protein of 624 amino acids was produced by RcM34, of which the amino terminal 582 amino acids would be identical to the predicted wt 854 amino acids product. Recent, re-annotations of the MCMV genome have identified three putative M34 overlapping ORFs (m33.1, m34.1, and m34.2) that may be interrupted in the above mutants. All three were transcribed in RcM34 virus-infected cells confirming that the RcM34 virus phenotype was probably due to interruption of the M34 ORF. These results suggest that M34, like human CMV UL34, is an essential gene.

Murine cytomegalovirus m38.5 protein inhibits Bax-mediated cell death

Many viruses encode proteins that inhibit the induction of programmed cell death at the mitochondrial checkpoint. Murine cytomegalovirus (MCMV) encodes the m38.5 protein, which localizes to mitochondria and protects human HeLa cells and fibroblasts from apoptosis triggered by proteasome inhibitors but not from Fas-induced apoptosis. However, the ability of this protein to suppress the apoptosis of murine cells and its role during MCMV infection have not been investigated previously. Here we show that m38.5 is expressed at early time points during MCMV infection. Cells infected with MCMVs lacking m38.5 showed increased sensitivity to cell death induced by staurosporine, MG132, or the viral infection itself compared to the sensitivity of cells infected with wild-type MCMV. This defect was eliminated when an m38.5 or Bcl-X(L) gene was inserted into the genome of a deletion mutant. Using fibroblasts deficient in the proapoptotic Bcl-2 family proteins Bak and/or Bax, we further demonstrated that m38.5 protected from Bax- but not Bak-mediated apoptosis and interacted with Bax in infected cells. These results consolidate the role of m38.5 as a viral mitochondrion-localized inhibitor of apoptosis and its functional similarity to the human cytomegalovirus UL37x1 gene product. Although the m38.5 gene is not homologous to the UL37x1 gene at the sequence level, m38.5 is conserved among rodent cytomegaloviruses. Moreover, the fact that MCMV-infected cells are protected from both Bak- and Bax-mediated cell death suggests that MCMV possesses an additional, as-yet-unidentified mechanism to block Bak-mediated apoptosis.

Cytomegalovirus proteins vMIA and m38.5 link mitochondrial morphogenesis to Bcl-2 family proteins

Apoptosis is a host defense mechanism against viruses that can be subverted by viral gene products. Human cytomegalovirus encodes viral mitochondria-localized inhibitor of apoptosis (vMIA; also known as pUL37x1), which is targeted to mitochondria and functions as a potent cell death suppressor by binding to and inhibiting proapoptotic Bcl-2 family members Bax and Bak. vMIA expression also dramatically alters mitochondrial morphology, causing the fragmentation of these organelles. A potential ortholog of vMIA, m38.5, which was identified in murine cytomegalovirus, has been shown to localize to mitochondria and protect against chemically induced apoptosis by unknown mechanisms. Despite sharing negligible homology with vMIA and no region detectably corresponding to the vMIA Bax-binding domain, we find that m38.5, like vMIA, binds to Bax and recruits Bax to mitochondria. Interestingly, m38.5 and vMIA appear to block Bax downstream of translocation to mitochondria and after an initial stage of Bax conformational change. In contrast to vMIA, m38.5 neither binds to Bak nor causes mitochondrial fragmentation. Consistently with Bax-selective inactivation by m38.5, m38.5 fragments mitochondria in Bak knockout (KO) cells and protects Bak KO cells from apoptosis better than Bax KO cells. Thus, vMIA and m38.5 share some, but not all, features of apoptosis regulation through Bcl-2 family interaction and allow the dissection of Bax translocation into discrete steps.

Calculation of splicing potential from the Alternative Splicing Mutation Database

BACKGROUND

The Alternative Splicing Mutation Database (ASMD) presents a collection of all known mutations inside human exons which affect splicing enhancers and silencers and cause changes in the alternative splicing pattern of the corresponding genes.

FINDINGS

An algorithm was developed to derive a Splicing Potential (SP) table from the ASMD information. This table characterizes the influence of each oligonucleotide on the splicing effectiveness of the exon containing it. If the SP value for an oligonucleotide is positive, it promotes exon retention, while negative SP values mean the sequence favors exon skipping. The merit of the SP approach is the ability to separate splicing signals from a wide range of sequence motifs enriched in exonic sequences that are attributed to protein-coding properties and/or translation efficiency. Due to its direct derivation from observed splice site selection, SP has an advantage over other computational approaches for predicting alternative splicing.

CONCLUSION

We show that a vast majority of known exonic splicing enhancers have highly positive cumulative SP values, while known splicing silencers have core motifs with strongly negative cumulative SP values. Our approach allows for computation of the cumulative SP value of any sequence segment and, thus, gives researchers the ability to measure the possible contribution of any sequence to the pattern of splicing.

Horizontal gene transfer and nucleotide compositional anomaly in large DNA viruses

Background

DNA viruses have a wide range of genome sizes (5 kb up to 1.2 Mb, compared to 0.16 Mb to 1.5 Mb for obligate parasitic bacteria) that do not correlate with their virulence or the taxonomic distribution of their hosts. The reasons for such large variation are unclear. According to the traditional view of viruses as gifted "gene pickpockets", large viral genome sizes could originate from numerous gene acquisitions from their hosts. We investigated this hypothesis by studying 67 large DNA viruses with genome sizes larger than 150 kb, including the recently characterized giant mimivirus. Given that horizontally transferred DNA often have anomalous nucleotide compositions differing from the rest of the genome, we conducted a detailed analysis of the inter- and intra-genome compositional properties of these viruses. We then interpreted their compositional heterogeneity in terms of possible causes, including strand asymmetry, gene function/expression, and horizontal transfer.

Results

We first show that the global nucleotide composition and nucleotide word usage of viral genomes are species-specific and distinct from those of their hosts. Next, we identified compositionally anomalous (cA) genes in viral genomes, using a method based on Bayesian inference. The proportion of cA genes is highly variable across viruses and does not exhibit a significant correlation with genome size. The vast majority of the cA genes were of unknown function, lacking homologs in the databases. For genes with known homologs, we found a substantial enrichment of cA genes in specific functional classes for some of the viruses. No significant association was found between cA genes and compositional strand asymmetry. A possible exogenous origin for a small fraction of the cA genes could be confirmed by phylogenetic reconstruction.

Conclusion

At odds with the traditional dogma, our results argue against frequent genetic transfers to large DNA viruses from their modern hosts. The large genome sizes of these viruses are not simply explained by an increased propensity to acquire foreign genes. This study also confirms that the anomalous nucleotide compositions of the cA genes is sometimes linked to particular biological functions or expression patterns, possibly leading to an overestimation of recent horizontal gene transfers.

Murine cytomegalovirus open reading frame m29.1 augments virus replication both in vitro and in vivo

Murine cytomegalovirus mutant Rc29, with a premature stop codon mutation in the m29 open reading frame (ORF), produced no apparent phenotype in cell culture or following infection of BALB/c mice. In contrast, a similar mutant virus, Rc29.1, with a premature stop codon mutation in its m29.1 ORF, showed reduced virus yields (2-3 log(10) p.f.u. ml(-1)) in tissue culture. Mutant virus yields in BALB/c mice were delayed, reduced ( approximately 1 log(10) p.f.u. per tissue) and persisted less well in salivary glands compared with wild-type (wt) and revertant (Rv29.1) virus. In severe combined immunodeficiency mice, Rc29.1 virus showed delayed and reduced replication initially in all tissues (liver, spleen, kidneys, heart, lung and salivary glands). This delayed death until 31 days post-infection (p.i.) compared with wt (23 days p.i.) but at death virus yields were similar to wt. m29 gene transcription was initiated at early times post-infection, while production of a transcript from ORF m29.1 in the presence of cycloheximide indicated that it was an immediate-early gene. ORFs m29.1 and M28 are expressed from a bicistronic message, which is spliced infrequently. However, it is likely that each ORF expresses its own protein, as antiserum derived in rabbits to the m29.1 protein expressed in bacteria from the m29.1 ORF detected only one protein in Western blot analysis of the size predicted for the m29.1 protein. Our results suggest that neither ORF is essential for virus replication but m29.1 is important for optimal viral growth in vitro and in vivo.

Cytomegalovirus evasion of innate immunity by subversion of the NKR-P1B:Clr-b missing-self axis

Cytomegaloviruses are known to encode several gene products that function to subvert MHC-dependent immune recognition. Here we characterize a rat cytomegalovirus (RCMV) C-type lectin-like (RCTL) gene product with homology to the Clr ligands for the NKR-P1 receptors. RCMV infection rapidly extinguished host Clr-b expression, thereby sensitizing infected cells to killing by natural killer (NK) cells. However, the RCTL protein functioned as a decoy ligand to protect infected cells from NK killing via direct interaction with the NKR-P1B inhibitory receptor. In vivo, an RCTL mutant virus displayed diminished virulence in an NK-dependent and strain-specific manner, suggesting that host NKR-P1 polymorphisms have evolved to avert the viral decoy mechanism while maintaining Clr-b recognition to preserve self tolerance. These findings reveal a unique strategy adopted by cytomegaloviruses to evade MHC-independent self-nonself discrimination. The existence of lectin-like genes in several poxviruses suggests that this may represent a common theme for viral evasion of innate immunity.

Extensive sequence variation exists among isolates of murine cytomegalovirus within members of the m02 family of genes

Murine cytomegalovirus (MCMV) is a widely used model for human cytomegalovirus (HCMV) and has facilitated many important discoveries about the biology of CMVs. Most of these studies are conducted using the laboratory MCMV strains Smith and K181. However, wild-derived isolates of MCMV, like HCMV clinical isolates, exhibit genetic variation from laboratory strains, particularly at the ends of their genomes in areas containing known or putative immune-evasion and tropism genes. This study analysed the nucleotide sequence of the m02-m05 region, within the m02 gene family, of a number of laboratory and wild-derived MCMV isolates, and found a large degree of variation in both the sequence and arrangement of genes. A new open reading frame (ORF), designated m03.5, was found to be present in a number of wild isolates of MCMV in place of m03. Two distinct isolates, W8 and W8211, were found to possess both m03 and m03.5. Both m03 and m03.5 had early transcription kinetics and the encoded proteins could be detected on the cell surface, consistent with a possible role in immune evasion through binding to host-cell proteins. These data show that gene duplication and sequence variation occur within different isolates of MCMV found in the wild. As this variation among strains may alter the function of genes, these findings should be considered when analysing gene function or host-virus interactions in laboratory models.

Human cytomegalovirus UL38 protein blocks apoptosis

Apoptosis is an innate cellular defense response to viral infection. The slow-replicating human cytomegalovirus (HCMV) blocks premature death of host cells prior to completion of the infection cycle. In this study, we report that the HCMV UL38 gene encodes a cell death inhibitory protein. A mutant virus lacking the pUL38 coding sequence, ADdlUL38, grew poorly in human fibroblasts, failed to accumulate viral DNA to wild-type levels, and induced excessive death of infected cells. Cells expressing pUL38 were resistant to cell death upon infection and effectively supported the growth of ADdlUL38. Cells infected with the pUL38-deficient virus showed morphological changes characteristic of apoptosis, including cell shrinkage, membrane blebbing, vesicle release, and chromatin condensation and fragmentation. The proteolytic cleavage of two key enzymes involved in apoptosis, namely, caspase 3 and poly(ADP-ribose) polymerase, was activated upon ADdlUL38 infection, and the cleavage was blocked in cells expressing pUL38. The pan-caspase inhibitor Z-VAD-FMK largely restored the growth of ADdlUL38 in normal fibroblasts, indicating that the defective growth of the mutant virus mainly resulted from premature death of host cells. Furthermore, cells expressing pUL38 were resistant to cell death induced by a mutant adenovirus lacking the antiapoptotic E1B-19K protein or by thapsigargin, which disrupts calcium homeostasis in the endoplasmic reticulum. Taken together, these results indicate that the HCMV protein pUL38 suppresses apoptosis, blocking premature death of host cells to facilitate efficient virus replication.

Double-stranded RNA binding by a heterodimeric complex of murine cytomegalovirus m142 and m143 proteins

In response to viral infection, cells activate a variety of antiviral responses, including several that are triggered by double-stranded (ds) RNA. Among these are the protein kinase R and oligoadenylate synthetase/RNase L pathways, both of which result in the shutoff of protein synthesis. Many viruses, including human cytomegalovirus, encode dsRNA-binding proteins that prevent the activation of these pathways and thereby enable continued protein synthesis and viral replication. We have extended these analyses to another member of the beta subfamily of herpesviruses, murine cytomegalovirus (MCMV), and now report that products of the m142 and m143 genes together bind dsRNA. Coimmunoprecipitation experiments demonstrate that these two proteins interact in infected cells, consistent with their previously reported colocalization. Jointly, but not individually, the proteins rescue replication of a vaccinia virus mutant with a deletion of the dsRNA-binding protein gene E3L (VVDeltaE3L). Like the human cytomegalovirus dsRNA-binding protein genes TRS1 and IRS1, m142 and m143 are members of the US22 gene family. We also found that two other members of the MCMV US22 family, M23 and M24, encode dsRNA-binding proteins, but they do not rescue VVDeltaE3L replication. These results reveal that MCMV, like many other viruses, encodes dsRNA-binding proteins, at least two of which can inhibit dsRNA-activated antiviral pathways. However, unlike other well-studied examples, the MCMV proteins appear to act in a heterodimeric complex.

Genes of murine cytomegalovirus exist as a number of distinct genotypes

Murine cytomegaloviruses encode a number of genes which modulate polymorphic host immune responses. We suggest that these viral genes should themselves therefore exhibit sequence polymorphism. Additionally, clinical isolates of human cytomegalovirus (HCMV) have been shown to vary extensively from the common laboratory strains. Almost all research conducted on murine cytomegalovirus (MCMV) has used the laboratory strains Smith and K181, which have been extensively passaged in vitro and in vivo since isolation. Using the heteroduplex mobility assay (HMA) to determine levels of sequence variation 11 MCMV genes were examined from 26 isolates of MCMV from wild mice, as well as both laboratory strains. Both the HMA and sequencing of selected genes demonstrated that whilst certain genes (M33, mck-2, m147.5, m152) were highly conserved, others (m04, m06, M44, m138, m144, m145 and m155) contained significant sequence variation. Several of these genes (m06, m144 and m155) exist in wild MCMV strains as one of several distinct genotypes.

Experimental confirmation of global murine cytomegalovirus open reading frames by transcriptional detection and partial characterization of newly described gene products

Murine cytomegalovirus (MCMV) and human CMV (HCMV) share many features making the mouse system a potential small-animal model for HCMV. Although the genomic DNA sequence and the predicted open reading frames (ORFs) of MCMV have been determined, experimental evidence that the ORFs are actually transcribed has been lacking. We developed an MCMV global-DNA microarray that includes all previously predicted ORFs and 14 potential ones. A total of 172 ORFs were confirmed to be transcribed, including 7 newly discovered ORFs not previously predicted. No gene products from 10 previously predicted ORFs were detected by either DNA microarray analysis or reverse transcriptase PCR in MCMV-infected mouse fibroblasts, although 2 of those were expressed in a macrophage cell line, suggesting that potential gene products from these open reading frames are silenced in fibroblasts and required in macrophages. Immunohistochemical localization of the six newly described ORF products and three recently identified ones in cells transfected with the respective construct revealed four of the products in the nucleus and five in mitochondria. Analysis of two ORFs using site-directed mutagenesis showed that deletion of one of the mitochondrion-localized gene products led to significantly decreased replication in fibroblasts.

Predicting calcium-binding sites in proteins - a graph theory and geometry approach

Identifying calcium-binding sites in proteins is one of the first steps towards predicting and understanding the role of calcium in biological systems for protein structure and function studies. Due to the complexity and irregularity of calcium-binding sites, a fast and accurate method for predicting and identifying calcium-binding protein is needed. Here we report our development of a new fast algorithm (GG) to detect calcium-binding sites. The GG algorithm uses a graph theory algorithm to find oxygen clusters of the protein and a geometric algorithm to identify the center of these clusters. A cluster of four or more oxygen atoms has a high potential for calcium binding. High performance with about 90% site sensitivity and 80% site selectivity has been obtained for three datasets containing a total of 123 proteins. The results suggest that a sphere of a certain size with four or more oxygen atoms on the surface and without other atoms inside is necessary and sufficient for quickly identifying the majority of the calcium-binding sites with high accuracy. Our finding opens a new avenue to visualize and analyze calcium-binding sites in proteins facilitating the prediction of functions from structural genomic information.

Using hidden Markov models and observed evolution to annotate viral genomes

MOTIVATION

ssRNA (single stranded) viral genomes are generally constrained in length and utilize overlapping reading frames to maximally exploit the coding potential within the genome length restrictions. This overlapping coding phenomenon leads to complex evolutionary constraints operating on the genome. In regions which code for more than one protein, silent mutations in one reading frame generally have a protein coding effect in another. To maximize coding flexibility in all reading frames, overlapping regions are often compositionally biased towards amino acids which are 6-fold degenerate with respect to the 64 codon alphabet. Previous methodologies have used this fact in an ad hoc manner to look for overlapping genes by motif matching. In this paper differentiated nucleotide compositional patterns in overlapping regions are incorporated into a probabilistic hidden Markov model (HMM) framework which is used to annotate ssRNA viral genomes. This work focuses on single sequence annotation and applies an HMM framework to ssRNA viral annotation. A description of how the HMM is parameterized, whilst annotating within a missing data framework is given. A Phylogenetic HMM (Phylo-HMM) extension, as applied to 14 aligned HIV2 sequences is also presented. This evolutionary extension serves as an illustration of the potential of the Phylo-HMM framework for ssRNA viral genomic annotation.

RESULTS

The single sequence annotation procedure (SSA) is applied to 14 different strains of the HIV2 virus. Further results on alternative ssRNA viral genomes are presented to illustrate more generally the performance of the method. The results of the SSA method are encouraging however there is still room for improvement, and since there is overwhelming evidence to indicate that comparative methods can improve coding sequence (CDS) annotation, the SSA method is extended to a Phylo-HMM to incorporate evolutionary information. The Phylo-HMM extension is applied to the same set of 14 HIV2 sequences which are pre-aligned. The performance improvement that results from including the evolutionary information in the analysis is illustrated.

Contribution of GADD45 family members to cell death suppression by cellular Bcl-xL and cytomegalovirus vMIA

Mammalian cells and viruses encode inhibitors of programmed cell death that localize to mitochondria and suppress apoptosis initiated by a wide variety of inducers. Mutagenesis was used to probe the role of a predicted alpha-helical region within the hydrophobic antiapoptotic domain (AAD) of cytomegalovirus vMIA, the UL37x1 gene product. This region was found to be essential for cell death suppression activity. A screen for proteins that interacted with the AAD of functional vMIA but that failed to interact with mutants identified growth arrest and DNA damage 45 (GADD45alpha), a cell cycle regulatory protein activated by genotoxic stress, as a candidate cellular binding partner. GADD45alpha interaction required the AAD alpha-helical character that also dictated GADD45alpha-mediated enhancement of death suppression. vMIA mutants that failed to interact with GADD45alpha were completely nonfunctional in cell death suppression, and any of the three GADD45 family members (GADD45alpha, GADD45beta/MyD118, or GADD45gamma/OIG37/CR6/GRP17) was able to cooperate with vMIA; however, none influenced cell death when introduced into cells alone. GADD45alpha was found to increase vMIA protein levels comparably to treatment with protease inhibitors MG132 and ALLN. Targeted short interfering RNA knockdown of all three GADD45 family members maximally reduced vMIA activity, and this reduction was abrogated by additional GADD45alpha. Interestingly, GADD45 family members were also able to bind and enhance cell death suppression by Bcl-xL, a member of the Bcl-2 family of cell death suppressors, suggesting a direct cooperative link between apoptosis and the proteins that regulate the DNA damage response.