Roles of avian herpesvirus microRNAs in infection, latency, and oncogenesis.

MicroRNAs have been reported for the avian herpesviruses Marek's disease virus 1 (MDV1; oncogenic), Marek's disease virus 2 (MDV2; non-oncogenic), herpesvirus of turkeys (HVT), and infectious laryngotracheitis virus (ILTV). No obvious phylogenetic relationships exist among the avian herpesvirus microRNAs, but the general genomic locations of microRNA clusters are conserved, with microRNAs being located in the repeat regions of the genomes. In some cases, microRNAs are antisense to open reading frames. Among MDV1 field isolates with different virulence properties, microRNAs are highly conserved, and variations that have been observed lie in putative promoter regions. One cluster of MDV1 microRNAs lies upstream of the meq gene, and this cluster is more highly expressed in tumors caused by an extremely virulent MDV1 isolate compared to tumors caused by a less virulent isolate. Several of the avian herpesvirus microRNAs are orthologs of microRNAs in other species. For example, mdv1-miR-M4 shares a seed sequence with gga-miR-155 (also shared with Kaposi sarcoma herpesvirus (KSHV) kshv-miR-K12), mdv2-miR-M21 shares a seed with miR-29b, and hvt-miR-H14 shares a seed sequence with miR-221. Functional analyses of avian herpesvirus microRNAs include a variety of in vitro assays to demonstrate potential function as well as the use of mutants that can exploit the ability to assess phenotypes experimentally in the natural host. This article is part of a Special Issue entitled:MicroRNA's in viral gene regulation.

Sequence conservation and differential expression of Marek's disease virus microRNAs.

Marek's disease virus (MDV), a herpesvirus that causes a lymphoproliferative disorder in chickens, encodes a number of microRNAs derived primarily from two locations in the MDV genome. One cluster of microRNA genes flanks the meq oncogene, and a second cluster is found within the latency-associated transcript (LAT) region. The sequences of MDV microRNAs from a collection of field and reference strains with various levels of virulence were compared and found to be highly conserved. However, microRNAs from the meq cluster were detected at higher levels in lymphomas caused by a form of the virus designated very virulent plus (vv+; strain 615K, also known as T. King) than in those caused by a less virulent (very virulent [vv]) form (RB1B). For example, levels of mdv1-miR-M4, which shares a seed sequence with miR-155, a microRNA implicated in B-cell lymphoma, were threefold higher and levels of mdv1-miR-M2*/3p were more than sixfold higher in vv+ MDV-induced tumors than in vv MDV-induced tumors. In contrast, levels of the microRNAs from the LAT cluster were equivalent in tumors produced by vv and vv+ strains. Additionally, mdv1-miR-M4 is the MDV microRNA most highly expressed in tumors, where it accounts for 72% of all MDV microRNAs, as determined by deep sequencing. These data suggest that the meq cluster microRNAs play an important role in the pathogenicity of MDV.

Deep sequencing of chicken microRNAs.

BACKGROUND: The use of new, deep sequencing technologies has greatly accelerated microRNA discovery. We have applied this approach to the identification of chicken microRNAs and to the comparison of microRNAs in chicken embryo fibroblasts (CEF) infected with Marek's disease virus (MDV) to those present in uninfected CEF. RESULTS: We obtained 125,463 high quality reads that showed an exact match to the chicken genome. The majority of the reads corresponded to previously annotated chicken microRNAs; however, the sequences of many potential novel microsRNAs were obtained. A comparison of the reads obtained in MDV-infected and uninfected CEF indicates that infection does not significantly perturb the expression profile of microRNAs. Frequently sequenced microRNAs include miR-221/222, which are thought to play a role in growth and proliferation. A number of microRNAs (e.g., let-7, miR-199a-1, 26a) are expressed at lower levels in MDV-induced tumors, highlighting the potential importance of this class of molecules in tumorigenesis. CONCLUSION: Deep sequencing technology is highly suited for small RNA discovery. This approach is independent of comparative sequence analysis, which has been the primary method used to identify chicken microRNAs. Our results have confirmed the expression of many microRNAs identified by sequence similarity and identified a pool of candidate novel microRNAs.

A novel interferon regulatory factor (IRF), IRF-10, has a unique role in immune defense and is induced by the v-Rel oncoprotein.

The cloning and functional characterization of a novel interferon regulatory factor (IRF), IRF-10, are described. IRF-10 is most closely related to IRF-4 but differs in both its constitutive and inducible expression. The expression of IRF-10 is inducible by interferons (IFNs) and by concanavalin A. In contrast to that of other IRFs, the inducible expression of IRF-10 is characterized by delayed kinetics and requires protein synthesis, suggesting a unique role in the later stages of an antiviral defense. Accordingly, IRF-10 is involved in the upregulation of two primary IFN-gamma target genes (major histocompatibility complex [MHC] class I and guanylate-binding protein) and interferes with the induction of the type I IFN target gene for 2',5'-oligo(A) synthetase. IRF-10 binds the interferon-stimulated response element site of the MHC class I promoter. In contrast to that of IRF-1, which has some of the same functional characteristics, the expression of IRF-10 is not cytotoxic for fibroblasts or B cells. The expression of IRF-10 is induced by the oncogene v-rel, the proto-oncogene c-rel, and IRF-4 in a tissue-specific manner. Moreover, v-Rel and IRF-4 synergistically cooperate in the induction of IRF-10 in fibroblasts. The level of IRF-10 induction in lymphoid cell lines by Rel proteins correlates with Rel transformation potential. These results suggest that IRF-10 plays a role in the late stages of an immune defense by regulating the expression some of the IFN-gamma target genes in the absence of a cytotoxic effect. Furthermore, IRF-10 expression is regulated, at least in part, by members of the Rel/NF-kappa B and IRF families.

Analysis of gene expression, copy number and palindrome formation with a Dt40 enriched cDNA microarray.

DT40 presents a unique opportunity to exploit newly available tools for chicken genomic analysis. A 13K chicken cDNA microarray representing 11447 non-overlapping ESTs has been developed. This array detects expression of 7086 DT40 genes of which_644 are over-expressed 3-fold or greater and 1585 are under-expressed 3-fold or greater relative to normal post-hatch bursal cell populations. Changes in RNA expression due to single gene alterations can be detected by expression profiling. For example, by this method, over expression of the oncogenic micro RNA bic up-regulates expression of VBP, a known regulator of Avian Leukosis Virus LTR- driven transcription with very little additional expression change, A degree of cytogenetic abnormality and instability of DT40 cells has been observed, which is characterized at the fine structure level using microarray-based comparative genome hybridization (array-CGH). The relationship between gene copy number and RNA expression levels can be assessed in the same tissue samples using the same microarray. A newly introduced technique for genome-wide analysis of palindrome formation (GAPF) detects long inverted repeats, or palindromes, which are early events in gene amplification and possibly other DNA structural change. Since both array CGH-detected copy number changes and GAPF-detected palindromes are abundant in DT40, these techniques, coupled with targeted gene deletion and replacement, may provide a powerful tool for analysis of genomic instability and its underlying genetic mechanisms.

Chicken genomics resource: sequencing and annotation of 35,407 ESTs from single and multiple tissue cDNA libraries and CAP3 assembly of a chicken gene index.

Its accessibility, unique evolutionary position, and recently assembled genome sequence have advanced the chicken to the forefront of comparative genomics and developmental biology research as a model organism. Several chicken expressed sequence tag (EST) projects have placed the chicken in 10th place for accrued ESTs among all organisms in GenBank. We have completed the single-pass 5'-end sequencing of 37,557 chicken cDNA clones from several single and multiple tissue cDNA libraries and have entered 35,407 EST sequences into GenBank. Our chicken EST sequences and those found in public databases (on July 1, 2004) provided a total of 517,727 public chicken ESTs and mRNAs. These sequences were used in the CAP3 assembly of a chicken gene index composed of 40,850 contigs and 79,192 unassembled singlets. The CAP3 contigs show a 96.7% match to the chicken genome sequence. The University of Delaware (UD) EST collection (43,928 clones) was assembled into 19,237 nonredundant sequences (13,495 contigs and 5,742 unassembled singlets). The UD collection contains 6,223 unique sequences that are not found in other public EST collections but show a 76% match to the chicken genome sequence. Our chicken contig and singlet sequences were annotated according to the highest BlastX and/or BlastN hits. The UD CAP3 contig assemblies and singlets are searchable by nucleotide sequence or key word (http://cogburn.dbi.udel.edu), and the cDNA clones are readily available for distribution from the chick EST website and clone repository (http://www.chickest.udel.edu). The present paper describes the construction and normalization of single and multiple tissue chicken cDNA libraries, high-throughput EST sequencing from these libraries, the CAP3 assembly of a chicken gene index from all public ESTs, and the identification of several nonredundant chicken gene sets for production of custom DNA microarrays.

Development of a cDNA array for chicken gene expression analysis.

BACKGROUND: The application of microarray technology to functional genomic analysis in the chicken has been limited by the lack of arrays containing large numbers of genes. RESULTS: We have produced cDNA arrays using chicken EST collections generated by BBSRC, University of Delaware and the Fred Hutchinson Cancer Research Center. From a total of 363,838 chicken ESTs representing 24 different adult or embryonic tissues, a set of 11,447 non-redundant ESTs were selected and added to an existing collection of clones (4,162) from immune tissues and a chicken bursal cell line (DT40). Quality control analysis indicates there are 13,007 useable features on the array, including 160 control spots. The array provides broad coverage of mRNAs expressed in many tissues; in addition, clones with expression unique to various tissues can be detected. CONCLUSIONS: A chicken multi-tissue cDNA microarray with 13,007 features is now available to academic researchers from http://genomics@fhcrc.org. Sequence information for all features on the array is in GenBank, and clones can be readily obtained. Targeted users include researchers in comparative and developmental biology, immunology, vaccine and agricultural technology. These arrays will be an important resource for the entire research community using the chicken as a model.

Conservation of biological properties of the CD40 ligand, CD154 in a non-mammalian vertebrate.

Signals delivered by the CD40 ligand, CD154, have crucial roles in immune responses in mammals, being required for development of germinal centres, maturation of T-dependent antibody responses, and generation of B-cell memory. To determine whether these functions were conserved in a non-mammalian species, a putative chicken CD 154 cDNA was used to make an oligomeric fusion protein, and to raise monoclonal antibodies. The antibodies detected surface expression on activated T-cells. The fusion protein detected expression of a receptor on B-cells, thrombocytes and macrophages. Biological effects of the fusion protein included induction of NO synthesis in a macrophage cell line, enhancement of splenic B-cell survival, and induction of apoptosis in a bursal lymphoma cell line. These observations demonstrated substantial functional equivalence with mammalian CD 154 and thus provided evidence for the early evolutionary emergence of the set of functions associated with this molecule, and its central role in the vertebrate immune system.

Jak3-regulated genes: DNA array analysis of concanavalin a-interleukin-2-activated chicken T cells treated with a specific jak3 inhibitor.

Janus kinase 3 (Jak3) is important in the activation and proliferation of lymphoid cells and binds to the common gamma subunit of several cytokine receptors, including the interleukin-2 (IL-2) receptor (IL-2R). DNA arrays were used to measure mRNA levels of a large number of genes regulated by signaling through the Jak3 tyrosine kinase pathway by blocking concanavalin A (ConA)-IL-2-activated chicken splenic T cells with a specific Jak3 inhibitor (WHI-P154). Of the 635 genes detected by arrays containing about 1200 cDNAs, 12 were upregulated in control cells compared with inhibitor-treated cells, and 6 were expressed at higher levels in the inhibitor-treated group. By identifying genes that are directly or indirectly regulated by Jak3, we can gain insight into the roles of this key intermediate in avian T cell activation and further our understanding of intracellular signaling networks in the immune response.

Emerging roles of chicken and viral microRNAs in avian disease.

BACKGROUND: MicroRNAs are short RNAs (~22 nt) expressed by plants, animals and viruses that regulate gene expression post-transcriptionally, and their importance is highlighted by distinct patterns of expression in many physiological processes, including development, hematopoeisis, stress resistance, and disease. Our group has characterized the microRNAs encoded by the avian herpesviruses; namely, oncogenic Marek's disease (MD) virus (MDV1), non-oncogenic MDV (MDV2) herpesvirus of turkeys (HVT), and infectious laryngotracheitis virus (ILTV). METHODS: MicroRNAs encoded by the avian herpesviruses were identified using next generation sequencing technologies (454, Illumina). RESULTS: The microRNAs of each the avian herpesviruses have unique sequences, but the genomic locations are similar, in that the microRNAs tend to be clustered in the rapidly evolving repeat regions of the viral genomes. For a given viral species the microRNA sequence is highly conserved in different strains with the exception of a virulence-associated polymorphism in the putative promoter of the MDV1 microRNAs upstream of the meq oncogene. These microRNAs are relatively highly expressed in tumors produced by very virulent MDV1 isolates compared to tumors produced by less virulent strains. MDV1 and HVT encode homologs of the host microRNA, miR-221, which targets a gene important in cell cycle regulation. MDV1 encodes a microRNA (mdv1-miR-M4) that shares a seed sequence with miR-155, a microRNA important in immune function. Mdv-miR-M4 is highly expressed in MDV induced tumors, while miR-155 is present at very low levels. CONCLUSIONS: MicroRNAs are highly conserved among different field strains of MDV1, and they are expressed in lytic and latent infections and in MDV1-derived tumors. This suggests that these small molecules are very important to the virus, and roles in immune evasion, anti-apoptosis, or proliferation are likely.

A microRNA of infectious laryngotracheitis virus can downregulate and direct cleavage of ICP4 mRNA.

Viral microRNAs regulate gene expression using either translational repression or mRNA cleavage and decay. Two microRNAs from infectious laryngotracheitis virus (ILTV), iltv-miR-I5 and iltv-miR-I6, map antisense to the ICP4 gene. Post-transcriptional repression by these microRNAs was tested against a portion of the ICP4 coding sequence cloned downstream of firefly luciferase. Luciferase activity was downregulated by approximately 60% with the iltv-miR-I5 mimic. Addition of an iltv-miR-I5 antagomiR or mutagenesis of the target seed sequence alleviated this effect. The iltv-miR-I5 mimic, when co-transfected with a plasmid expressing ICP4, reduced ICP4 transcript levels by approximately 50%, and inhibition was relieved by an iltv-miR-I5 antagomiR. In infected cells, iltv-miR-I5 mediated cleavage at the canonical site, as indicated by modified RACE analysis. Thus, in this system, iltv-miR-I5 decreased ILTV ICP4 mRNA levels via transcript cleavage and degradation. Downregulation of ICP4 could impact the balance between the lytic and latent states of the virus in vivo.

MicroRNAs of Gallid and Meleagrid herpesviruses show generally conserved genomic locations and are virus-specific.

Many herpesviruses, including Marek's disease viruses (MDV1 and MDV2), encode microRNAs. In this study, we report microRNAs of two related herpesviruses, infectious laryngotracheitis virus (ILTV) and herpesvirus of turkeys (HVT), as well as additional MDV2 microRNAs. The genome locations, but not microRNA sequences, are conserved among all four of these avian herpesviruses. Most are clustered in the repeats flanking the unique long region (I/TR(L)), except in ILTV which lacks these repeats. Two abundant ILTV microRNAs are antisense to the immediate early gene ICP4. A homologue of host microRNA, gga-miR-221, was found among the HVT microRNAs. Additionally, a cluster of HVT microRNAs was found in a region containing two locally duplicated segments, resulting in paralogous HVT microRNAs with 96-100% identity. The prevalence of microRNAs in the genomic repeat regions as well as in local repeats suggests the importance of genetic plasticity in herpesviruses for microRNA evolution and preservation of function.

Marek's disease virus up-regulates major histocompatibility complex class II cell surface expression in infected cells.

Many herpesviruses modulate major histocompatibility complex (MHC) expression on the cell surface as an immune evasion mechanism. We report here that Marek's disease virus (MDV), a lymphotrophic avian alphaherpesvirus, up-regulates MHC class II cell surface expression in infected cells, contrary to all other herpesviruses examined to date. This MDV-induced class II up-regulation was detected both in vitro and in vivo. This effect was not solely an indirect effect of interferon, which is a highly potent natural inducer of MHC class II expression, since MHC class II up-regulation in cultured primary fibroblast cells was confined to the infected cells only. MHC class II up-regulation was also observed in infected cells of the bursa of Fabricius during the lytic phase of MDV infection in birds and upon reactivation of MDV from latency in an MDV-transformed cell line. As MDV is a strictly cell-associated virus and requires activated T cells for its life cycle, this up-regulation of MHC class II in infected cells may contribute to virus spread within the infected host by increasing the chance of contact between productively infected cells and susceptible activated T cells.

Marek's disease virus encodes MicroRNAs that map to meq and the latency-associated transcript.

MicroRNAs (miRNAs) are a class of small (approximately 22-nucleotide) regulatory molecules that block translation or induce degradation of target mRNAs. These have been identified in a wide range of organisms, including viruses. In particular, the oncogenic gammaherpesviruses Kaposi's sarcoma herpesvirus and Epstein-Barr virus encode miRNAs that could potentially regulate either viral or host genes. To determine if Marek's disease virus (MDV), an oncogenic alphaherpesvirus of chickens, encodes miRNAs, we isolated small RNAs from MDV-infected chicken embryo fibroblasts (CEF) and used the 454 Life Sciences sequencing technology to obtain the sequences of 13,679 candidate host and viral small RNAs. Eight miRNAs were found, five of which flank the meq oncogene and three that map to the latency-associated transcript (LAT) region of the genome. The meq gene is unique to pathogenic serotypes of MDV and is transcriptionally active during latency and transformation, and the LAT region of the MDV genome is antisense to the immediate-early gene ICP4. Secondary structure analysis predicted that the regions flanking the miRNAs could form hairpin precursors. Northern blot analysis confirmed expression of all miRNAs in MDV-infected CEF, MDV-induced tumors, and MDV lymphoblastoid cell lines. We propose that the MDV miRNAs function to enable MDV pathogenesis and contribute to MDV-induced transformation of chicken T cells.

Transcriptome analysis for the chicken based on 19,626 finished cDNA sequences and 485,337 expressed sequence tags.

We present an analysis of the chicken (Gallus gallus) transcriptome based on the full insert sequences for 19,626 cDNAs, combined with 485,337 EST sequences. The cDNA data set has been functionally annotated and describes a minimum of 11,929 chicken coding genes, including the sequence for 2260 full-length cDNAs together with a collection of noncoding (nc) cDNAs that have been stringently filtered to remove untranslated regions of coding mRNAs. The combined collection of cDNAs and ESTs describe 62,546 clustered transcripts and provide transcriptional evidence for a total of 18,989 chicken genes, including 88% of the annotated Ensembl gene set. Analysis of the ncRNAs reveals a set that is highly conserved in chickens and mammals, including sequences for 14 pri-miRNAs encoding 23 different miRNAs. The data sets described here provide a transcriptome toolkit linked to physical clones for bioinformaticians and experimental biologists who wish to use chicken systems as a low-cost, accessible alternative to mammals for the analysis of vertebrate development, immunology, and cell biology.

Patterns of gene expression in the developing chick thymus.

Gene expression in thymic T cells during late embryogenesis and early growth in chicks was examined using cDNA microarrays. Gene expression patterns were profiled into nine clusters by using self-organizing maps (SOM) clustering analysis. The expression patterns for a set of genes confirmed current information on the development of immune response. Expression of cell surface markers (MHC class I alpha chain, MHC class II associated invariant chain, CD8 beta chain, CD18, and beta2-microglobulin), and genes involved in the innate immune response (NK lysin-like) increased with age, and these patterns were consistent with an increase in the immune responsiveness of the young chick. The expression of cytokine receptor common gamma chain (gammac), death receptor-3 (DR3), and TCR alpha chain increased up to 1 day of age and then decreased. DR3 could play a role in the apoptosis during T-cell maturation, while the differential expression of TCR genes could reflect regulation of the rearrangement of TCR genes and TCR-mediated signal transduction during T cell development. Three genes coding for previously uncharacterized proteins are included in the clusters. These gene expression profiling studies provide background information on the developing chick immune system and provide preliminary functional information on unknown proteins.

Five hundred sixty-five triples of chicken, human, and mouse candidate orthologs.

The Human Genome Project has provided abundant gene sequence information on human and important model organisms. The chicken is well positioned from an evolutionary standpoint to serve as a link between higher and lower organisms, particularly mammals, and amphibia and fish. In this study we used stringent criteria to select 565 triples of chicken, human, and mouse candidate orthologs. We analyze the sequences with respect to nucleotide and amino acid similarities. This analysis also allows measurement of evolutionary distances of different proteins. We found that chicken-human and chicken-mouse sequence identities are highly correlated; similarly for chicken-human and chicken-mouse evolutionary distances. With chicken as the out-group, we found that mouse has a higher substitution rate than human, supporting the generation-time effect hypothesis. We also described the transversion bias, which is the preference for some transversions than others in nucleotide substitutions. We demonstrated that there are statistically significant properties in the differences of orthologous sequences. The differential patterns, in combination with sequence similarity analysis, may lead to the identification of genes that are very divergent from the mammalian orthologs.

Herpesvirus of turkeys: microarray analysis of host gene responses to infection.

Herpesvirus of turkeys (HVT) provides an economically important live vaccine for prevention of Marek's disease (MD) of chickens. MD, characterized by both immunosuppression and T-cell lymphoma, is caused by another herpesvirus termed Marek's disease virus (MDV). Microarrays were used to investigate the response of chicken embryonic fibroblasts (CEF) to infection with HVT. Genes responding to HVT infection include several induced by interferon along with others modulating signal transduction, transcription, scaffolding proteins, and the cytoskeleton. Results are compared with earlier studies examining the responses of CEF cells to infection with MDV.

Differential cytokine expression in avian cells in response to invasion by Salmonella typhimurium, Salmonella enteritidis and Salmonella gallinarum.

Salmonella enterica is a facultative intracellular pathogen that is capable of causing disease in a range of hosts. Although human salmonellosis is frequently associated with consumption of contaminated poultry and eggs, and the serotypes Salmonella gallinarum and Salmonella pullorum are important world-wide pathogens of poultry, little is understood of the mechanisms of pathogenesis of Salmonella in the chicken. Type III secretion systems play a key role in host cell invasiveness and trigger the production of pro-inflammatory cytokines during invasion of mammalian hosts. This results in a polymorphonuclear cell influx that contributes to the resulting enteritis. In this study, a chicken primary cell culture model was used to investigate the cytokine responses to entry by the broad host range serotypes S. enteritidis and S. typhimurium, and the host specific serotype S. gallinarum, which rarely causes disease outside its main host, the chicken. The cytokines interleukin (IL)-1ss, IL-2, IL-6 and interferon (IFN)-gamma were measured by quantitative RT-PCR, and production of IL-6 and IFN-gamma was also determined through bioassays. All serotypes were invasive and had little effect on the production of IFN-gamma compared with non-infected cells; S. enteritidis invasion caused a slight down-regulation of IL-2 production. For IL-1ss production, infection with S. typhimurium had little effect, whilst infection with S. gallinarum or S. enteritidis caused a reduction in IL-1ss mRNA levels. Invasion of S. typhimurium and S. enteritidis caused an eight- to tenfold increase in production of the pro-inflammatory cytokine IL-6, whilst invasion by S. gallinarum caused no increase. These findings correlate with the pathogenesis of Salmonella in poultry. S. typhimurium and S. enteritidis invasion produces a strong inflammatory response, that may limit the spread of Salmonella largely to the gut, whilst S. gallinarum does not induce an inflammatory response and may not be limited by the immune system, leading to the severe systemic disease fowl typhoid.