Positive-strand RNA virus genome replication organelles: structure, assembly, control.

Positive-strand RNA [(+)RNA] viruses include pandemic SARS-CoV-2, tumor-inducing hepatitis C virus, debilitating chikungunya virus (CHIKV), lethal encephalitis viruses, and many other major pathogens. (+)RNA viruses replicate their RNA genomes in virus-induced replication organelles (ROs) that also evolve new viral species and variants by recombination and mutation and are crucial virus control targets. Recent cryo-electron microscopy (cryo-EM) reveals that viral RNA replication proteins form striking ringed 'crowns' at RO vesicle junctions with the cytosol. These crowns direct RO vesicle formation, viral (-)RNA and (+)RNA synthesis and capping, innate immune escape, and transfer of progeny (+)RNA genomes into translation and encapsidation. Ongoing studies are illuminating crown assembly, sequential functions, host factor interactions, etc., with significant implications for control and beneficial uses of viruses.

Nodavirus RNA replication crown architecture reveals proto-crown precursor and viral protein A conformational switching.

Positive-strand RNA viruses replicate their genomes in virus-induced membrane vesicles, and the resulting RNA replication complexes are a major target for virus control. Nodavirus studies first revealed viral RNA replication proteins forming a 12-fold symmetric "crown" at the vesicle opening to the cytosol, an arrangement recently confirmed to extend to distantly related alphaviruses. Using cryoelectron microscopy (cryo-EM), we show that mature nodavirus crowns comprise two stacked 12-mer rings of multidomain viral RNA replication protein A. Each ring contains an ~19 nm circle of C-proximal polymerase domains, differentiated by strikingly diverged positions of N-proximal RNA capping/membrane binding domains. The lower ring is a "proto-crown" precursor that assembles prior to RNA template recruitment, RNA synthesis, and replication vesicle formation. In this proto-crown, the N-proximal segments interact to form a toroidal central floor, whose 3.1 Å resolution structure reveals many mechanistic details of the RNA capping/membrane binding domains. In the upper ring, cryo-EM fitting indicates that the N-proximal domains extend radially outside the polymerases, forming separated, membrane-binding "legs." The polymerase and N-proximal domains are connected by a long linker accommodating the conformational switch between the two rings and possibly also polymerase movements associated with RNA synthesis and nonsymmetric electron density in the lower center of mature crowns. The results reveal remarkable viral protein multifunctionality, conformational flexibility, and evolutionary plasticity and insights into (+)RNA virus replication and control.

Subdomain cryo-EM structure of nodaviral replication protein A crown complex provides mechanistic insights into RNA genome replication.

For positive-strand RNA [(+)RNA] viruses, the major target for antiviral therapies is genomic RNA replication, which occurs at poorly understood membrane-bound viral RNA replication complexes. Recent cryoelectron microscopy (cryo-EM) of nodavirus RNA replication complexes revealed that the viral double-stranded RNA replication template is coiled inside a 30- to 90-nm invagination of the outer mitochondrial membrane, whose necked aperture to the cytoplasm is gated by a 12-fold symmetric, 35-nm diameter "crown" complex that contains multifunctional viral RNA replication protein A. Here we report optimizing cryo-EM tomography and image processing to improve crown resolution from 33 to 8.5 Å. This resolves the crown into 12 distinct vertical segments, each with 3 major subdomains: A membrane-connected basal lobe and an apical lobe that together comprise the ∼19-nm-diameter central turret, and a leg emerging from the basal lobe that connects to the membrane at ∼35-nm diameter. Despite widely varying replication vesicle diameters, the resulting two rings of membrane interaction sites constrain the vesicle neck to a highly uniform shape. Labeling protein A with a His-tag that binds 5-nm Ni-nanogold allowed cryo-EM tomography mapping of the C terminus of protein A to the apical lobe, which correlates well with the predicted structure of the C-proximal polymerase domain of protein A. These and other results indicate that the crown contains 12 copies of protein A arranged basally to apically in an N-to-C orientation. Moreover, the apical polymerase localization has significant mechanistic implications for template RNA recruitment and (-) and (+)RNA synthesis.

Human papillomavirus oncogenes reprogram the cervical cancer microenvironment independently of and synergistically with estrogen.

High-risk human papillomaviruses (HPVs) infect epithelial cells and are causally associated with cervical cancer, but HPV infection is not sufficient for carcinogenesis. Previously, we reported that estrogen signaling in the stromal tumor microenvironment is associated with cervical cancer maintenance and progression. We have now determined how HPV oncogenes and estrogen treatment affect genome-wide host gene expression in laser-captured regions of the cervical epithelium and stroma of untreated or estrogen-treated nontransgenic and HPV-transgenic mice. HPV oncogene expression in the cervical epithelium elicited significant gene-expression changes in the proximal stromal compartment, and estrogen treatment uniquely affected gene expression in the cervical microenvironment of HPV-transgenic mice compared with nontransgenic mice. Several potential estrogen-induced paracrine-acting factors were identified in the expression profile of the cervical tumor microenvironment. The microenvironment of estrogen-treated HPV-transgenic mice was significantly enriched for chemokine/cytokine activity and inflammatory and immune functions associated with carcinogenesis. This inflammatory signature included several proangiogenic CXCR2 receptor ligands. A subset of the same CXCR2 ligands was likewise increased in cocultures of early-passage cells from human cervical samples, with levels highest in cocultures of cervical fibroblasts and cancer-derived epithelial cells. Our studies demonstrate that high-risk HPV oncogenes profoundly reprogram the tumor microenvironment independently of and synergistically with estrogen. These observations illuminate important means by which HPVs can cause cancer through alterations in the tumor microenvironment.

Molecular transitions from papillomavirus infection to cervical precancer and cancer: Role of stromal estrogen receptor signaling.

To study the multistep process of cervical cancer development, we analyzed 128 frozen cervical samples spanning normalcy, increasingly severe cervical intraepithelial neoplasia (CIN1- CIN3), and cervical cancer (CxCa) from multiple perspectives, revealing a cascade of progressive changes. Compared with normal tissue, expression of many DNA replication/repair and cell proliferation genes was increased in CIN1/CIN2 lesions and further sustained in CIN3, consistent with high-risk human papillomavirus (HPV)-induced tumor suppressor inactivation. The CIN3-to-CxCa transition showed metabolic shifts, including decreased expression of mitochondrial electron transport complex components and ribosomal protein genes. Significantly, despite clinical, epidemiological, and animal model results linking estrogen and estrogen receptor alpha (ERα) to CxCa, ERα expression declined >15-fold from normalcy to cancer, showing the strongest inverse correlation of any gene with the increasing expression of p16, a marker for HPV-linked cancers. This drop in ERα in CIN and tumor cells was confirmed at the protein level. However, ERα expression in stromal cells continued throughout CxCa development. Our further studies localized stromal ERα to FSP1+, CD34+, SMA- precursor fibrocytes adjacent to normal and precancerous CIN epithelium, and FSP1-, CD34-, SMA+ activated fibroblasts in CxCas. Moreover, rank correlations with ERα mRNA identified IL-8, CXCL12, CXCL14, their receptors, and other angiogenesis and immune cell infiltration and inflammatory factors as candidates for ERα-induced stroma-tumor signaling pathways. The results indicate that estrogen signaling in cervical cancer has dramatic differences from ERα+ breast cancers, and imply that estrogen signaling increasingly proceeds indirectly through ERα in tumor-associated stromal fibroblasts.

Organelle-like membrane compartmentalization of positive-strand RNA virus replication factories.

Positive-strand RNA virus genome replication is invariably associated with extensively rearranged intracellular membranes. Recent biochemical and electron microscopy analyses, including three-dimensional electron microscope tomographic imaging, have fundamentally advanced our understanding of the ultrastructure and function of organelle-like RNA replication factories. Notably, for a range of positive-strand RNA viruses embodying many major differences, independent studies have revealed multiple common principles. These principles include that RNA replication often occurs inside numerous virus-induced vesicles invaginated or otherwise elaborated from a continuous, often endoplasmic reticulum-derived membrane network. Where analyzed, each such vesicle typically contains only one or a few genome replication intermediates in conjunction with many copies of viral nonstructural proteins. In addition, these genome replication compartments often are closely associated with sites of virion assembly and budding. Our understanding of these complexes is growing, providing substantial new insights into the organization, coordination, and potential control of crucial processes in virus replication.

Key aspects of papillomavirus infection influence the host cervicovaginal microbiome in a preclinical murine papillomavirus (MmuPV1) infection model.

Human papillomaviruses (HPVs) are the most common sexually transmitted infection in the United States and are a major etiological agent of cancers in the anogenital tract and oral cavity. Growing evidence suggests changes in the host microbiome are associated with the natural history and ultimate outcome of HPV infection. We sought to define changes in the host cervicovaginal microbiome during papillomavirus infection, persistence, and pathogenesis using the murine papillomavirus (MmuPV1) cervicovaginal infection model. Cervicovaginal lavages were performed over a time course of MmuPV1 infection in immunocompetent female FVB/N mice and extracted DNA was analyzed by qPCR to track MmuPV1 viral copy number. 16S ribosomal RNA (rRNA) gene sequencing was used to determine the composition and diversity of microbial communities throughout this time course. We also sought to determine whether specific microbial communities exist across the spectrum of MmuPV1-induced neoplastic disease. We, therefore, performed laser-capture microdissection to isolate regions of disease representing all stages of neoplastic disease progression (normal, low- and high-grade dysplasia, and cancer) from female reproductive tract tissue sections from MmuPV1-infected mice and performed 16S rRNA sequencing. Consistent with other studies, we found that the natural murine cervicovaginal microbiome is highly variable across different experiments. Despite these differences in initial microbiome composition between experiments, we observed that MmuPV1 persistence, viral load, and severity of disease influenced the composition of the cervicovaginal microbiome. These studies demonstrate that papillomavirus infection can alter the cervicovaginal microbiome.IMPORTANCEHuman papillomaviruses (HPVs) are the most common sexually transmitted infection in the United States. A subset of HPVs that infect the anogenital tract (cervix, vagina, anus) and oral cavity cause at least 5% of cancers worldwide. Recent evidence indicates that the community of microbial organisms present in the human cervix and vagina, known as the cervicovaginal microbiome, plays a role in HPV-induced cervical cancer. However, the mechanisms underlying this interplay are not well-defined. In this study, we infected the female reproductive tract of mice with a murine papillomavirus (MmuPV1) and found that key aspects of papillomavirus infection and disease influence the host cervicovaginal microbiome. This is the first study to define changes in the host microbiome associated with MmuPV1 infection in a preclinical animal model of HPV-induced cervical cancer. These results pave the way for using MmuPV1 infection models to further investigate the interactions between papillomaviruses and the host microbiome.

Intersection of the multivesicular body pathway and lipid homeostasis in RNA replication by a positive-strand RNA virus.

Like many positive-strand RNA viruses, brome mosaic virus (BMV) RNA replication occurs in membrane-invaginated vesicular compartments. BMV RNA replication compartments show parallels with membrane-enveloped, budding retrovirus virions, whose release depends on the cellular multivesicular body (MVB) sorting pathway. BMV RNA replication compartments are not released from their parent membranes, but might depend on MVB functions for membrane invagination. Prior results show that BMV RNA replication is severely inhibited by deletion of the crucial MVB gene DOA4 or BRO1. We report here that involvement of DOA4 and BRO1 in BMV RNA replication is not dependent on the MVB pathway's membrane-shaping functions but rather is due to their roles in recycling ubiquitin from MVB cargos. We show that deleting DOA4 or BRO1 inhibits the ubiquitination- and proteasome-dependent activation of homologous transcription factors Mga2p and Spt23p, which regulate many lipid metabolism genes, including the fatty acid desaturase gene OLE1, which is essential for BMV RNA replication. However, Mga2p processing and BMV RNA replication are restored by supplementing free ubiquitin, which is depleted in doa4Δ and bro1Δ cells. The results identify Mga2p and Spt23p processing and lipid regulation as sensitive targets of ubiquitin depletion and correctly predict multiple effects of modulating additional host genes RFU1, UBP6, and UFD3. Our results also show that BMV RNA replication depends on additional Mga2p-regulated genes likely involved in lipid metabolism beyond OLE1. Among other points, these findings show the potential for blocking viral RNA replication by modulating lipid synthesis at multiple levels.

Crowning Touches in Positive-Strand RNA Virus Genome Replication Complex Structure and Function.

Positive-strand RNA viruses, the largest genetic class of eukaryotic viruses, include coronaviruses and many other established and emerging pathogens. A major target for understanding and controlling these viruses is their genome replication, which occurs in virus-induced membrane vesicles that organize replication steps and protect double-stranded RNA intermediates from innate immune recognition. The structure of these complexes has been greatly illuminated by recent cryo-electron microscope tomography studies with several viruses. One key finding in diverse systems is the organization of crucial viral RNA replication factors in multimeric rings or crowns that among other functions serve as exit channels gating release of progeny genomes to the cytosol for translation and encapsidation. Emerging results suggest that these crowns serve additional important purposes in replication complex assembly, function, and interaction with downstream processes such as encapsidation. The findings provide insights into viral function and evolution and new bases for understanding, controlling, and engineering positive-strand RNA viruses.

Multifunctional Protein A Is the Only Viral Protein Required for Nodavirus RNA Replication Crown Formation.

Positive-strand RNA virus RNA genome replication occurs in membrane-associated RNA replication complexes (RCs). Nodavirus RCs are outer mitochondrial membrane invaginations whose necked openings to the cytosol are "crowned" by a 12-fold symmetrical proteinaceous ring that functions as the main engine of RNA replication. Similar protein crowns recently visualized at the openings of alphavirus and coronavirus RCs highlight their broad conservation and functional importance. Using cryo-EM tomography, we earlier showed that the major nodavirus crown constituent is viral protein A, whose polymerase, RNA capping, membrane interaction and multimerization domains drive RC formation and function. Other viral proteins are strong candidates for unassigned EM density in the crown. RNA-binding RNAi inhibitor protein B2 co-immunoprecipitates with protein A and could form crown subdomains that protect nascent viral RNA and dsRNA templates. Capsid protein may interact with the crown since nodavirus virion assembly has spatial and other links to RNA replication. Using cryoelectron tomography and complementary approaches, we show that, even when formed in mammalian cells, nodavirus RC crowns generated without B2 and capsid proteins are functional and structurally indistinguishable from mature crowns in infected Drosophila cells expressing all viral proteins. Thus, the only nodaviral factors essential to form functional RCs and crowns are RNA replication protein A and an RNA template. We also resolve apparent conflicts in prior results on B2 localization in infected cells, revealing at least two distinguishable pools of B2. The results have significant implications for crown structure, assembly, function and control as an antiviral target.

An amphipathic alpha-helix controls multiple roles of brome mosaic virus protein 1a in RNA replication complex assembly and function.

Brome mosaic virus (BMV) protein 1a has multiple key roles in viral RNA replication. 1a localizes to perinuclear endoplasmic reticulum (ER) membranes as a peripheral membrane protein, induces ER membrane invaginations in which RNA replication complexes form, and recruits and stabilizes BMV 2a polymerase (2a(Pol)) and RNA replication templates at these sites to establish active replication complexes. During replication, 1a provides RNA capping, NTPase and possibly RNA helicase functions. Here we identify in BMV 1a an amphipathic alpha-helix, helix A, and use NMR analysis to define its structure and propensity to insert in hydrophobic membrane-mimicking micelles. We show that helix A is essential for efficient 1a-ER membrane association and normal perinuclear ER localization, and that deletion or mutation of helix A abolishes RNA replication. Strikingly, mutations in helix A give rise to two dramatically opposite 1a function phenotypes, implying that helix A acts as a molecular switch regulating the intricate balance between separable 1a functions. One class of helix A deletions and amino acid substitutions markedly inhibits 1a-membrane association and abolishes ER membrane invagination, viral RNA template recruitment, and replication, but doubles the 1a-mediated increase in 2a(Pol) accumulation. The second class of helix A mutations not only maintains efficient 1a-membrane association but also amplifies the number of 1a-induced membrane invaginations 5- to 8-fold and enhances viral RNA template recruitment, while failing to stimulate 2a(Pol) accumulation. The results provide new insights into the pathways of RNA replication complex assembly and show that helix A is critical for assembly and function of the viral RNA replication complex, including its central role in targeting replication components and controlling modes of 1a action.

MicroRNA 29c is down-regulated in nasopharyngeal carcinomas, up-regulating mRNAs encoding extracellular matrix proteins.

Using highly sensitive microarray-based procedures, we identified eight microRNAs (miRNAs) showing robust differential expression between 31 laser-capture-microdissected nasopharyngeal carcinomas (NPCs) and 10 normal healthy nasopharyngeal epithelial samples. In particular, miRNA mir-29c was expressed at one-fifth the levels in tumors as in normal epithelium. In NPC tumors, the lower mir-29c levels correlated with higher levels of multiple mRNAs whose 3' UTRs can bind mir-29c at target sequences conserved across many vertebrates. In cultured cells, introduction of mir-29c down-regulated these genes at the level of mRNA and inhibited expression of luciferase encoded by vectors having the 3' UTRs of these genes. Moreover, for each of several genes tested, mutating the mir-29c target sites in the 3' UTR abrogated mir-29c-induced inhibition of luciferase expression. Most of the mir-29c-targeted genes identified encode extracellular matrix proteins, including multiple collagens and laminin gamma1, that are associated with tumor cell invasiveness and metastatic potential, prominent characteristics of NPC. Thus, we identify eight miRNAs differentially expressed in NPC and demonstrate the involvement of one in regulating genes involved in metastasis.

Cryo-electron microscopy of nodavirus RNA replication organelles illuminates positive-strand RNA virus genome replication.

The nodavirus flock house virus recently provided a well-characterized model for the first cryo-electron microscope tomography of membrane-bound, positive-strand RNA ((+)RNA) virus genome replication complexes (RCs). The resulting first views of RC organization and complementary biochemical results showed that the viral RNA replication vesicle is tightly packed with the dsRNA genomic RNA replication intermediate, and that (+)ssRNA replication products are released through the vesicle neck to the cytosol through a 12-fold symmetric ring or crown of multi-functional viral RNA replication proteins, which likely also contribute to viral RNA synthesis. Subsequent studies identified similar crown-like RNA replication protein complexes in alphavirus and coronavirus RCs, indicating related mechanisms across highly divergent (+)RNA viruses. As outlined in this review, these results have significant implications for viral function, evolution and control.

Fundamental differences in cell cycle deregulation in human papillomavirus-positive and human papillomavirus-negative head/neck and cervical cancers.

Human papillomaviruses (HPV) are associated with nearly all cervical cancers, 20% to 30% of head and neck cancers (HNC), and other cancers. Because HNCs also arise in HPV-negative patients, this type of cancer provides unique opportunities to define similarities and differences of HPV-positive versus HPV-negative cancers arising in the same tissue. Here, we describe genome-wide expression profiling of 84 HNCs, cervical cancers, and site-matched normal epithelial samples in which we used laser capture microdissection to enrich samples for tumor-derived versus normal epithelial cells. This analysis revealed that HPV(+) HNCs and cervical cancers differed in their patterns of gene expression yet shared many changes compared with HPV(-) HNCs. Some of these shared changes were predicted, but many others were not. Notably, HPV(+) HNCs and cervical cancers were found to be up-regulated in their expression of a distinct and larger subset of cell cycle genes than that observed in HPV(-) HNC. Moreover, HPV(+) cancers overexpressed testis-specific genes that are normally expressed only in meiotic cells. Many, although not all, of the hallmark differences between HPV(+) HNC and HPV(-) HNC were a direct consequence of HPV and in particular the viral E6 and E7 oncogenes. This included a novel association of HPV oncogenes with testis-specific gene expression. These findings in primary human tumors provide novel biomarkers for early detection of HPV(+) and HPV(-) cancers, and emphasize the potential value of targeting E6 and E7 function, alone or combined with radiation and/or traditional chemotherapy, in the treatment of HPV(+) cancers.

Cryo-electron tomography reveals novel features of a viral RNA replication compartment.

Positive-strand RNA viruses, the largest genetic class of viruses, include numerous important pathogens such as Zika virus. These viruses replicate their RNA genomes in novel, membrane-bounded mini-organelles, but the organization of viral proteins and RNAs in these compartments has been largely unknown. We used cryo-electron tomography to reveal many previously unrecognized features of Flock house nodavirus (FHV) RNA replication compartments. These spherular invaginations of outer mitochondrial membranes are packed with electron-dense RNA fibrils and their volumes are closely correlated with RNA replication template length. Each spherule's necked aperture is crowned by a striking cupped ring structure containing multifunctional FHV RNA replication protein A. Subtomogram averaging of these crowns revealed twelve-fold symmetry, concentric flanking protrusions, and a central electron density. Many crowns were associated with long cytoplasmic fibrils, likely to be exported progeny RNA. These results provide new mechanistic insights into positive-strand RNA virus replication compartment structure, assembly, function and control.

Genes involved in DNA repair and nitrosamine metabolism and those located on chromosome 14q32 are dysregulated in nasopharyngeal carcinoma.

Polymorphisms in nitrosamine metabolism, DNA repair, and immune response genes have been associated with nasopharyngeal carcinoma (NPC). Studies have suggested chromosomal regions involved in NPC. To shed light on NPC etiology, we evaluated host gene expression patterns in 31 NPC and 10 normal nasopharyngeal tissue specimens using the Affymetrix Human Genome U133 Plus 2.0 Array. We focused on genes in five a priori biological pathways and chromosomal locations. Rates of differential expression within these prespecified lists and overall were tested using a bootstrap method. Differential expression was observed for 7.6% of probe sets overall. Elevations in rate of differential expression were observed within the DNA repair (13.7%; P = 0.01) and nitrosamine metabolism (17.5%; P = 0.04) pathways. Differentially expressed probe sets within the DNA repair pathway were consistently overexpressed (93%), with strong effects observed for PRKDC, PCNA, and CHEK1. Differentially expressed probe sets within the nitrosamine metabolism pathway were consistently underexpressed (100%), with strong effects observed for NQ01, CYP2B6, and CYP2E1. No significant evidence of increases in rate of differential expression was seen within the immune/inflammatory pathway. A significant elevation in rate of differential expression was noted for chromosome 4p15.1-4q12 (13.0%; P = 0.04); both overexpression and underexpression were evident (38% and 62%, respectively). An elevation in the rate of differential expression on chromosome 14q32 was observed (11.3%; P = 0.06) with a consistent pattern of gene underexpression (100%; P < 0.0001). These effects were similar when excluding late-stage tumors. Our results suggest that nitrosamine activation and DNA repair are important in NPC. The consistent down-regulation of expression on chromosome 14q32 suggests loss of heterozygosity in this region.

Lack of association between EBV and breast carcinoma.

Multiple conflicting findings have been presented which indicate that EBV may be found in anywhere from 0% to 51% of breast carcinomas. When EBV has been found causally associated with other human cancers, its DNA and one or more of its viral products have been detected in most tumor cells of a given biopsy. To test whether EBV has such an association with breast cancer, we measured the number of viral DNA molecules per cell in matched normal and tumor biopsies from 45 patients using real-time quantitative PCR. In no case could EBV DNA consistently be detected, with either of two different probes, at levels above 0.1 molecules per cell in two sections of the tumor samples. These levels of detection match those detected in EBV-negative cell lines and therefore likely represent noise in the assays. Equally importantly, the distribution of these low signals was the same between tumors and their matched normal controls. We conclude that EBV does not contribute to the development of breast cancers as it does to epithelial cancers such as nasopharyngeal and gastric carcinomas or to Burkitt's and Hodgkin's lymphomas.

Cytoplasmic viral replication complexes.

Many viruses that replicate in the cytoplasm compartmentalize their genome replication and transcription in organelle-like structures that enhance replication efficiency and protection from host defenses. In particular, recent studies with diverse positive-strand RNA viruses have further elucidated the ultrastructure of membrane-bound RNA replication complexes and how these complexes function in close coordination with virion assembly and budding. The structure, function, and assembly of some positive-strand RNA virus replication complexes have parallels and potential evolutionary links with the replicative cores of double-strand RNA virus and retrovirus virions and more general similarities with the replication factories of cytoplasmic DNA viruses.

Genome-wide expression profiling reveals EBV-associated inhibition of MHC class I expression in nasopharyngeal carcinoma.

To identify the molecular mechanisms by which EBV-associated epithelial cancers are maintained, we measured the expression of essentially all human genes and all latent EBV genes in a collection of 31 laser-captured, microdissected nasopharyngeal carcinoma (NPC) tissue samples and 10 normal nasopharyngeal tissues. Global gene expression profiles clearly distinguished tumors from normal healthy epithelium. Expression levels of six viral genes (EBNA1, EBNA2, EBNA3A, EBNA3B, LMP1, and LMP2A) were correlated among themselves and strongly inversely correlated with the expression of a large subset of host genes. Among the human genes whose inhibition was most strongly correlated with increased EBV gene expression were multiple MHC class I HLA genes involved in regulating immune response via antigen presentation. The association between EBV gene expression and inhibition of MHC class I HLA expression implies that antigen display is either directly inhibited by EBV, facilitating immune evasion by tumor cells, and/or that tumor cells with inhibited presentation are selected for their ability to sustain higher levels of EBV to take maximum advantage of EBV oncogene-mediated tumor-promoting actions. Our data clearly reflect such tumor promotion, showing that deregulation of key proteins involved in apoptosis (BCL2-related protein A1 and Fas apoptotic inhibitory molecule), cell cycle checkpoints (AKIP, SCYL1, and NIN), and metastasis (matrix metalloproteinase 1) is closely correlated with the levels of EBV gene expression in NPC.