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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Identification and molecular cloning of functional chicken IL-12.

By a combination of large-scale sequencing, bioinformatics, and traditional molecular biology, we identified the long-searched-for cDNA sequences encoding the homologues of the chicken IL-12p35 and IL-12p40 chains. These molecules are the first discovered nonmammalian IL-12 subunits. The homologies of the chicken IL-12p35 and IL-12p40 proteins to the corresponding known subunits of various species, i.e., humans, sheep, horse, cat, bovine, mouse, and woodchuck, ranged between 21 and 42%, respectively. The expression of IL-12 subunits was observed in lymphoid cells and proved to be dependent on the cell type and stimulus, while expression was not detected in stimulated primary chicken embryo fibroblast cells. Following transient expression of both molecules in COS-7 cells, we confirmed the necessity of heterodimerization into IL-12p70 to yield bioactivity as was also shown for its mammalian counterparts. The chicken IL-12p70 molecule, generated either by transient coexpression of monomeric IL-12p35 and monomeric IL-12p40 or as a fusion protein (as in a fusion linker construct), induced IFN-gamma synthesis and proliferative activity of freshly exposed chicken splenocytes. The high degree of functional similarity between chicken IL-12 and IL-12 of higher mammalian vertebrates, despite their poor sequence homology, illustrates the conservation and vital importance of the IL-12 molecule since the evolutionary dichotomy of birds and mammals >300 million years ago. In this article, we describe the first nonmammalian IL-12 molecule and show that this chicken IL-12 molecule is bioactive.

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.

Microarray analysis of Lyn-deficient B cells reveals germinal center-associated nuclear protein and other genes associated with the lymphoid germinal center.

Lyn is the only member of the Src family expressed in DT40 B cells, which provide a unique model to study the singular contribution of this protein tyrosine kinase (PTK) family to cell signaling. In these cells, gene ablation of Lyn leads to defective B cell receptor signaling. Complementary DNA array analysis of Lyn-deficient DT40 cells shows that the absence of Lyn leads to down-regulation of numerous genes encoding proteins involved in B cell receptor signaling, proliferation, control of transcription, immunity/inflammation response, and cytoskeletal organization. Most of these expression changes have not been previously associated with Lyn PTK signaling. They include alterations in mRNA levels of germinal center-associated nuclear protein (germinal center-associated DNA primase) (GANP), CD74, CD22, NF-kappaB, elongation factor 1alpha, CD79b, octamer binding factor 1, Ig H chain, stathmin, and gamma-actin. Changes in GANP expression were also confirmed in Lyn-deficient mice, suggesting that Lyn PTK has a unique function not compensated for by other Src kinases. Because Lyn-deficient mice have impaired development of germinal centers in spleen, the decreased expression of GANP in the Lyn-deficient DT40 cell line and Lyn-deficient mice suggests that Lyn controls the formation and proliferation of germinal centers via GANP. GANP promoter activity was higher in wild-type vs Lyn-deficient cells. Mutation of the PU.1 binding site reduced activity in wild-type cells and had no effect in Lyn-deficient cells. The presence of Lyn enhanced PU.1 expression in a Northern blot. Thus, the following new signaling pathway has been described: Lyn-->PU.1-->GANP.

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.

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.

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.