Analysis of the swine tracheobronchial lymph node transcriptomic response to infection with a Chinese highly pathogenic strain of porcine reproductive and respiratory syndrome virus
© Miller et al.; licensee BioMed Central Ltd. 2012
Received: 31 July 2012
Accepted: 15 October 2012
Published: 30 October 2012
Porcine reproductive and respiratory syndrome virus (PRRSV) is a major pathogen of swine worldwide. Emergence in 2006 of a novel highly pathogenic PRRSV (HP-PRRSV) isolate in China necessitated a comparative investigation into the host transcriptome response in tracheobronchial lymph nodes (TBLN) 13 days post-infection with HP-PRRSV rJXwn06, PRRSV strain VR-2332 or sham inocula. RNA from each was prepared for next-generation sequencing. Amplified library constructs were directly sequenced and a list of sequence transcripts and counts was generated using an RNAseq analysis pipeline to determine differential gene expression. Transcripts were annotated and relative abundance was calculated based upon the number of times a given transcript was represented in the library.
Major changes in transcript abundance occurred in response to infection with either PRRSV strain, each with over 630 differentially expressed transcripts. The largest increase in transcript level for either virus versus sham-inoculated controls were three serum amyloid A2 acute-phase isoforms. However, the degree of up or down-regulation of transcripts following infection with HP-PRRSV rJXwn06 was greater than transcript changes observed with US PRRSV VR-2332. Also, of 632 significantly altered transcripts within the HP-PRRSV rJXwn06 library 55 were up-regulated and 69 were down-regulated more than 3-fold, whilst in the US PRRSV VR-2332 library only 4 transcripts were up-regulated and 116 were down-regulated more than 3-fold.
The magnitude of differentially expressed gene profiles detected in HP-PRRSV rJXwn06 infected pigs as compared to VR-2332 infected pigs was consistent with the increased pathogenicity of the HP-PRRSV in vivo.
Porcine reproductive and respiratory syndrome virus (PRRSV), the causative agent of PRRS in swine, is a member of the Arteriviridae family in the order Nidovirales. PRRSV causes highly significant economic losses to the swine industry worldwide  as a result of both reproductive failure (late-term abortions and stillbirths) in pregnant sows and respiratory disease (pneumonia) in nursery and grower/finishing pigs . Infection with PRRSV also predisposes pigs to infection by bacterial pathogens as well as other viral pathogens [3–7], as such, PRRSV is a key etiological agent of the porcine respiratory disease complex (PRDC). Clinical disease caused by PRRSV is highly variable, ranging from mild, subclinical infection to acute death of adult animals . Differences in virulence have been attributed to numerous factors including host genetics, management practices, and virus strain heterogeneity [9–16]. Relatively little is known about the interactions of PRRSV and host cells. The lymph node is an anatomic site where the innate immune response and adaptive immune system interface. Tracheobronchial lymph nodes (TBLN) in swine drain the lung field and provide the focal structure that can reproducibly be identified. Although the TBLN contains a number of cell types, sampling this tissue allows study of direct and indirect effects of an infectious agent on the lung and cells within the lymph node.
In 2006 a unique syndrome with high morbidity and mortality was recognized in growing pigs in China that was originally known as porcine high fever disease (PHFD) due to its uncertain etiology . Experimental infection of pigs in China with these novel viral isolates reproduced the clinical disease providing strong evidence for the role of PRRSV as the causal agent of PHFD. However, there was still a question as to whether there was some unknown agent in the PRRSV preparations that increased the severity of the clinical disease over what was expected for a “routine” PRRSV infection. This question was resolved when PHFD was reproduced in China with virus derived from an infectious clone of the JX143 PRRSV isolate  demonstrating that PRRSV isolates with a common genetic motif had a causal role in PHFD leading to this lineage of virus being called highly pathogenic PRRSV (HP-PRRSV). We imported a plasmid containing a full-length clone of the 2006 JXwn06 HP-PRRSV isolate  from which infectious virus (rJXwn06) was rescued. An animal study was conducted comparing the pathogenicity of HP-PRRSV isolate rJXwn06 with the North American prototype strain VR-2332 PRRSV . The objective of this report was to investigate gene expression profiles in porcine tracheobronchial lymph node (TBLN) during viral infection with HP-PRRSV rJXwn06 strain alongside of US PRRSV strain VR-2332 at a snapshot of 13 days post-infection using bioinformatics.
Results and discussion
Mapping short RNA-seq reads and estimating transcript expression levels
Genomic Short-read Nucleotide Alignment Program (GSNAP) was used for alignment and genome construction, and Cufflinks to determine if differential expression and changes in transcript abundance were statistically significant [21, 22]. The RNASeq yielded 55,527,464 reads for the control, 43,263,207 reads for the HP-PRRSV, and 34,555,783 for VR-2332 libraries after quality trimming and excluding any reads less than 25 bp. Cufflinks was used to measure transcript abundances in fragments per kilobase of exon per million fragments mapped (FPKM). The Cuffdiff output contained normalized FPKM for comparison between libraries (Additional file 1). These values were used to calculate the fold change (log2 transformed) in expression between the experimental unit and the control.
Table of transcript counts
Up-regulated genes >3-fold
Down-regulated genes >3-fold
HP-PRRSV rJXwn06 vs. Control
US PRRSV VR-2332 vs. Control
This derived catalog of expressed genes represents the first comparative analysis of the HP-PRRSV rJXwn06 and VR-2332-infected TBLN transcript abundance profiles and provides a database that informs us of genes involved in normal TBLN physiology, as well as genes whose abundance is altered by PRRSV infection.
Annotation of significantly differentially expressed genes
Gene annotation of all significant hits (Additional files 1 and 2) was then carried out using a MySQL database matching the Ensmbl (Sscrofa9.56) chromosome location of aligned transcripts to gene names. Gene IDs and log2 fold-change expression values for significant hits, that had FPKM values in both the control and the infected differential expression testing for transcripts (Cuffdiff output files), were then analyzed using the Ingenuity Pathway Analysis software. When comparing the TBLN transcriptome from sham-inoculated controls vs. the HP-PRRSV rJXwn06-infected pigs, 568 of the 632 gene IDs mapped to the Ingenuity Knowledge Base and 165 were up-regulated while 148 were down-regulated. In the TBLN of control vs. VR-2332-infected pigs, 528 of the 633 gene IDs mapped to the Ingenuity Knowledge Base and only 8 were up-regulated while 235 were down-regulated.
Top Ten HGNC named genes with their fold change obtained from the RNAseq data
US PRRSV VR-2332
Fold change (log2)
Fold change (log2)
This study produced transcriptional profiles of TBLNs from non-infected, HP-PRRSV rJXwn06 and US PRRSV VR-2332-infected pigs that provides insight into immune dysregulation elicited by the virus on host transcript abundance levels necessary for a effective immune response.
This RNA-Seq compendium extends the analyses of previous gene expression atlases performed using Affymetrix GeneChip technology and provides an example of new methods to accommodate the increase in transcriptome data obtained from next generation sequencing [29, 31, 32].
It is well established that many pathogens cause changes in expression of specific genes that act to protect the host and clear the infection. PRRSV strains differ in their dysregulation of the immune response to infection and delay in development of a protective immune response in vaccinated pigs [33, 34]. A higher number of significantly differentially expressed gene instances were detected in HP-PRRSV rJXwn06 than VR-2332 when normalized to control samples at a snapshot of 13 days post inoculation (dpi). As anticipated, some of the genes (e.g., resistin) and pathways identified would be expected to be involved in the host response to a severe disease. In the case of resistin, it would be expected that adipose tissue stores are being mobilized as part of the host response to infection, which includes a high fever typical of infection with HP-PRRSV. There are specific cellular proteins that regulate a protective immune response, for example the pro-inflammatory genes that were up-regulated to a greater extent in HP-PRRSV rJXwn06 than VR-2332 when normalized to control samples as observed when comparing the pathogenicity of HP-PRRSV isolate rJXwn06 with the North American prototype strain VR-2332 PRRSV. At 13 dpi HP-PRRSV rJXwn06 inoculated pigs had an interstitial pneumonia that was significantly more severe than theVR-2332 inoculated group which appeared to be convalescing . Future studies of these differentially expressed genes, their transcript abundance, protein level, and protein function will enhance our understanding of the interaction of PRRSV with the host. Identification of new virulence mechanisms of PRRSV may improve the prospects for rational design of more effective vaccines to limit viral replication and shedding.
Cells and virus
MARC-145 cells were cultured in minimum essential medium (EMEM, SAFC 56416C) with 10% fetal bovine serum at 37°C, 5% CO2. Wild-type (wt) Type 2 PRRSV strain VR-2332 (GenBank U87392), passage 6 on MARC-145 cells, was titrated and used for the swine study. Virus (rescued JXwn06; rJXwn06) was rescued from a cloned cDNA of Chinese highly pathogenic Type 2 PRRSV strain JXwn06 [pWSK-JXwn; GenBank EF641008, ] and passaged 3 times on MARC-145 cells for use in the swine study.
The animal use protocol was reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) of the National Animal Disease Center-USDA-Agricultural Research Service. Thirty-two 10-week-old cross-bred pigs were obtained from a U.S. high-health herd and were found to be free of PRRSV and influenza virus antibodies using commercially available enzyme-linked immunosorbent assay (ELISA) kits (HerdChek PRRS 2XR; IDEXX Laboratories, Westbrook, Maine) and NP ELISA (MultiS ELISA, IDEXX, Westbrook, Maine), respectively. Pigs were also confirmed negative for porcine circovirus type 2 by quantitative real-time PCR . One day prior to starting the experiment, pigs were bled, weighed and randomly assigned to one of four groups. Group 1 (N = 8) consisted of negative control pigs, which received an intranasal 2 ml sham inoculum of minimum essential media (MEM) on 0 dpi. Group 2 pigs (N = 12) were challenged intranasally with 2 ml of 1 × 106 50% tissue culture infective dose (TCID50)/ml of Chinese PRRSV strain rJXwn06 in Animal Biosafety Level-3-Agriculture (ABSL-3-Ag) housing, where they remained for the duration of the experiment. Group 3 consisted of naïve pigs (N = 4) that were placed in contact with Group 2 swine on 2 dpi. Group 4 pigs (N = 8) were challenged intranasally with 2 ml of 1 × 106 TCID50/ml of Type 2 prototype strain VR-2332. Groups 1 and 2 were housed in separate isolation rooms in an ABSL2 facility. Animal care and euthanasia were conducted in accordance with the Report of the AVMA Panel on Euthansia and under the supervision of IACUC of NADC. Serum and bronchoalveolar lung lavage fluid (BALF) were tested for infectious virus as described previously . Lungs were scored for gross lesions  and sections fixed for histopathology. Swabs were collected from BALF, and various sites for bacterial isolation .
Following humane euthanasia, tracheobronchial lymph nodes (TBLN) from in vivo HP-PRRSV rJXwn06 (N = 8), US PRRSV VR-2332 (N = 7), or sham-infected pigs (N = 8) were harvested at 13 days post-infection and total cellular RNA was prepared as follows. One gram of TBLN from each pig was collected immediately upon necropsy, minced and stored in RNAlater (Life Technologies, Grand Island, NY) at −80°C until homogenized for extraction of total RNA with MagMAX™-96 for Microarrays Total RNA Isolation Kit (Applied Biosystems, Carlsbad, CA) using the manufacturer's protocol. The integrity of the RNA was confirmed with a 2100 Bioanalyzer and RNA 6000 Nano-chip (Agilent, Santa Clara, CA). The samples used had an average RNA Integrity Number (RIN) value of 7.8 and 28S:18S rRNA ratio of 1.9.
cDNA library construction
cDNA libraries were constructed from pooled total cellular RNA from the TBLN in each treatment group using TruSeq Sample Prep Kits (Illumina Inc., San Diego, CA) and sequenced by 2 × 100 paired-end sequencing on an Illumina HiSeq 2000 instrument.
RNA-Seq (Whole Transcriptome Shotgun Sequencing) pipeline
Ingenuity pathway analysis (IPA)
Datasets representing genes with altered expression profile derived from RNAseq analyses were imported into the Ingenuity Pathway Analysis Tool (IPA Tool; Ingenuity® Systems, Redwood City, CA, USA; http://www.ingenuity.com). In IPA, differentially expressed genes were mapped to genetic networks available in the Ingenuity database and then ranked by score.
The basis of the IPA program consists of the Ingenuity Pathway Knowledge Base (IPKB) that is derived from known functions and interactions of genes published in the literature. Thus, the IPA Tool allows the identification of biological networks, global functions within the host and functional pathways of a particular dataset. The program also gives the significance value of the differentially expressed genes, the other genes with which it interacts, and how the products of the genes directly or indirectly act on each other, including those not involved in the microarray analysis. The networks created are ranked depending on the number of significantly expressed genes they contain and also list diseases that were most significant (Figure 1).
We thank the following members of the Virus and Prion Research Unit at the National Animal Disease Center: J. Huegel, J. Crabtree, A. Burow, D. Adolphson, S. Anderson, M. Kappes and A. Vorwald for technical assistance. We also gratefully acknowledge D. Alt of the Genomics Unit at the National Animal Disease Center, and A. Severin of Iowa State University NGS Bioinformatics, for assistance in data analysis. Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture.
- Neumann EJ, Kliebenstein JB, Johnson CD, Mabry JW, Bush EJ, Seitzinger AH, Green AL, Zimmerman JJ: Assessment of the economic impact of porcine reproductive and respiratory syndrome on swine production in the United States. J Am Vet Med Assoc. 2005, 227 (3): 385-392. 10.2460/javma.2005.227.385.View ArticlePubMed
- Zimmerman JJ, Yoon KJ, Wills RW, Swenson SL: General overview of PRRSV: a perspective from the United States. Vet Microbiol. 1997, 55 (1–4): 187-196.View ArticlePubMed
- Done SH, Paton DJ: Porcine reproductive and respiratory syndrome: clinical disease, pathology and immunosuppression. Vet Rec. 1995, 136 (2): 32-35. 10.1136/vr.136.2.32.View ArticlePubMed
- Feng W, Laster SM, Tompkins M, Brown T, Xu JS, Altier C, Gomez W, Benfield D, McCaw MB: In utero infection by porcine reproductive and respiratory syndrome virus is sufficient to increase susceptibility of piglets to challenge by Streptococcus suis type II. J Virol. 2001, 75 (10): 4889-4895. 10.1128/JVI.75.10.4889-4895.2001.PubMed CentralView ArticlePubMed
- Rossow KD, Collins JE, Goyal SM, Nelson EA, Christopher-Hennings J, Benfield DA: Pathogenesis of porcine reproductive and respiratory syndrome virus infection in gnotobiotic pigs. Vet Pathol. 1995, 32 (4): 361-373. 10.1177/030098589503200404.View ArticlePubMed
- Wills RW, Gray JT, Fedorka-Cray PJ, Yoon KJ, Ladely S, Zimmerman JJ: Synergism between porcine reproductive and respiratory syndrome virus (PRRSV) and Salmonella choleraesuis in swine [In Process Citation]. Vet Microbiol. 2000, 71 (3–4): 177-192.View ArticlePubMed
- Zeman D, Neiger R, Yaeger M, Nelson E, Benfield D, Leslie-Steen P, Thomson J, Miskimins D, Daly R, Minehart M: Laboratory investigation of PRRS virus infection in three swine herds. J Vet Diagn Invest. 1993, 5 (4): 522-528. 10.1177/104063879300500404.View ArticlePubMed
- Mengeling WL, Lager KM: A brief review of procedures and potential problems associated with the diagnosis of porcine reproductive and respiratory syndrome [In Process Citation]. Vet Res. 2000, 31 (1): 61-69.PubMed
- Goldberg TL, Hahn EC, Weigel RM, Scherba G: Genetic, geographical and temporal variation of porcine reproductive and respiratory syndrome virus in Illinois. J Gen Virol. 2000, 81 (Pt 1): 171-179.View ArticlePubMed
- Goldberg TL, Weigel RM, Hahn EC, Scherba G: Associations between genetics, farm characteristics and clinical disease in field outbreaks of porcine reproductive and respiratory syndrome virus [In Process Citation]. Prev Vet Med. 2000, 43 (4): 293-302. 10.1016/S0167-5877(99)00104-X.View ArticlePubMed
- Halbur PG, Paul PS, Frey ML, Landgraf J, Eernisse K, Meng XJ, Lum MA, Andrews JJ, Rathje JA: Comparison of the pathogenicity of two US porcine reproductive and respiratory syndrome virus isolates with that of the Lelystad virus. Vet Pathol. 1995, 32 (6): 648-660. 10.1177/030098589503200606.View ArticlePubMed
- Halbur PG, Paul PS, Meng XJ, Lum MA, Andrews JJ, Rathje JA: Comparative pathogenicity of nine US porcine reproductive and respiratory syndrome virus (PRRSV) isolates in a five-week-old cesarean- derived, colostrum-deprived pig model. J Vet Diagn Invest. 1996, 8 (1): 11-20. 10.1177/104063879600800103.View ArticlePubMed
- Keffaber KK: Reproductive failure of unknown etiology. Am Assoc Swine Pract News. 1989, 1: 1-9.
- Lewis CR, Ait-Ali T, Clapperton M, Archibald AL, Bishop S: Genetic perspectives on host responses to porcine reproductive and respiratory syndrome (PRRS). Viral Immunol. 2007, 20 (3): 343-358. 10.1089/vim.2007.0024.View ArticlePubMed
- Meng XJ: Heterogeneity of porcine reproductive and respiratory syndrome virus: implications for current vaccine efficacy and future vaccine development. Vet Microbiol. 2000, 74 (4): 309-329. 10.1016/S0378-1135(00)00196-6.View ArticlePubMed
- Wensvoort G: Lelystad virus and the porcine epidemic abortion and respiratory syndrome. Vet Res. 1993, 24 (2): 117-124.PubMed
- Tian K, Yu X, Zhao T, Feng Y, Cao Z, Wang C, Hu Y, Chen X, Hu D, Tian X, et al: Emergence of fatal PRRSV variants: unparalleled outbreaks of atypical PRRS in China and molecular dissection of the unique hallmark. PLoS One. 2007, 2 (6): e526-10.1371/journal.pone.0000526.PubMed CentralView ArticlePubMed
- Lv J, Zhang J, Sun Z, Liu W, Yuan S: An infectious cDNA clone of a highly pathogenic porcine reproductive and respiratory syndrome virus variant associated with porcine high fever syndrome. J Gen Virol. 2008, 89 (Pt 9): 2075-2079.View ArticlePubMed
- Zhou L, Zhang J, Zeng J, Yin S, Li Y, Zheng L, Guo X, Ge X, Yang H: The 30-amino-acid deletion in the Nsp2 of highly pathogenic porcine reproductive and respiratory syndrome virus emerging in China is not related to its virulence. J Virol. 2009, 83 (10): 5156-5167. 10.1128/JVI.02678-08.PubMed CentralView ArticlePubMed
- Guo B, Lager KM, Henningson JN, Miller LC, Schlink SN, Kappes MA, Kehrli ME, Brockmeier SL, Nicholson TL, Yang H-C, et al: Experimental infection of united states swine with a chinese highly pathogenic strain of porcine reproductive and respiratory syndrome virus. Virology. 2012, http://dx.doi.org/10.1016/j.virol.2012.09.013.
- Robertson G, Schein J, Chiu R, Corbett R, Field M, Jackman SD, Mungall K, Lee S, Okada HM, Qian JQ, et al: De novo assembly and analysis of RNA-seq data. Nat Methods. 2010, 7 (11): 909-912. 10.1038/nmeth.1517.View ArticlePubMed
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28 (5): 511-515. 10.1038/nbt.1621.PubMed CentralView ArticlePubMed
- Seal RL, Gordon SM, Lush MJ, Wright MW, Bruford EA: Genenames.org: the HGNC resources in 2011. Nucleic Acids Res. 2011, 39 (Database issue): D514-519.PubMed CentralView ArticlePubMed
- Dai MH, Xia T, Chen XD, Gan L, Feng SQ, Qiu H, Peng Y, Yang ZQ: Cloning and characterization of porcine resistin gene. Domest Anim Endocrinol. 2006, 30 (2): 88-97. 10.1016/j.domaniend.2005.06.003.View ArticlePubMed
- Chen H, Lunney JK, Cheng L, Li X, Cao J, Zhu M, Zhao S: Porcine S100A8 and S100A9: molecular characterizations and crucial functions in response to Haemophilus parasuis infection. Dev Comp Immunol. 2011, 35 (4): 490-500. 10.1016/j.dci.2010.11.017.View ArticlePubMed
- Yang W, Block ER: Effect of hypoxia and reoxygenation on the formation and release of reactive oxygen species by porcine pulmonary artery endothelial cells. J Cell Physiol. 1995, 164 (2): 414-423. 10.1002/jcp.1041640222.View ArticlePubMed
- Suzuki A, Yamada R, Yamamoto K: Citrullination by peptidylarginine deiminase in rheumatoid arthritis. Ann N Y Acad Sci. 2007, 1108: 323-339. 10.1196/annals.1422.034.View ArticlePubMed
- Cases S, Smith SJ, Zheng YW, Myers HM, Lear SR, Sande E, Novak S, Collins C, Welch CB, Lusis AJ, et al: Identification of a gene encoding an acyl CoA:diacylglycerol acyltransferase, a key enzyme in triacylglycerol synthesis. Proc Natl Acad Sci USA. 1998, 95 (22): 13018-13023. 10.1073/pnas.95.22.13018.PubMed CentralView ArticlePubMed
- Xiao S, Mo D, Wang Q, Jia J, Qin L, Yu X, Niu Y, Zhao X, Liu X, Chen Y: Aberrant host immune response induced by highly virulent PRRSV identified by digital gene expression tag profiling. BMC Genomics. 2010, 11: 544-10.1186/1471-2164-11-544.PubMed CentralView ArticlePubMed
- Lager KM, Faaberg KS, Guo B, Brockmeier SL, Henningson JN, Schlink SN, Miller LC, Kappes MA, Kehrli ME, Loving CL, et al: HP-PRRSV challenge of 4 and 10-week-old pigs. 5th PRRSpective Symposium: 2011; Bangkok, Thailand. Germany: Boehringer Ingelheim Animal Health GmbH; 2011;
- Miller LC, Neill JD, Harhay GP, Lager KM, Laegreid WW, Kehrli ME: In-depth global analysis of transcript abundance levels in porcine alveolar macrophages following infection with porcine reproductive and respiratory syndrome virus. Adv Virol. 2010, 2010: 864181.PubMed CentralView ArticlePubMed
- Zhou P, Zhai S, Zhou X, Lin P, Jiang T, Hu X, Jiang Y, Wu B, Zhang Q, Xu X, et al: Molecular characterization of transcriptome-wide interactions between highly pathogenic porcine reproductive and respiratory syndrome virus and porcine alveolar macrophages in vivo. Int J Biol Sci. 2011, 7 (7): 947-959.PubMed CentralView ArticlePubMed
- Lemke CD, Haynes JS, Spaete R, Adolphson D, Vorwald A, Lager K, Butler JE: Lymphoid hyperplasia resulting in immune dysregulation is caused by porcine reproductive and respiratory syndrome virus infection in neonatal pigs. J Immunol. 2004, 172 (3): 1916-1925.View ArticlePubMed
- Wesley RD, Lager KM, Kehrli ME: Infection with Porcine reproductive and respiratory syndrome virus stimulates an early gamma interferon response in the serum of pigs. Can J Vet Res. 2006, 70 (3): 176-182.PubMed CentralPubMed
- Opriessnig T, Yu S, Gallup JM, Evans RB, Fenaux M, Pallares F, Thacker EL, Brockus CW, Ackermann MR, Thomas P, et al: Effect of vaccination with selective bacterins on conventional pigs infected with type 2 porcine circovirus. Vet Pathol. 2003, 40 (5): 521-529. 10.1354/vp.40-5-521.View ArticlePubMed
- Faaberg KS, Kehrli ME, Lager KM, Guo B, Han J: In vivo growth of porcine reproductive and respiratory syndrome virus engineered nsp2 deletion mutants. Virus Res. 2010, 154 (1–2): 77-85.View ArticlePubMed
- Brockmeier SL, Loving CL, Nicholson TL, Palmer MV: Coinfection of pigs with porcine respiratory coronavirus and Bordetella bronchiseptica. Vet Microbiol. 2008, 128 (1–2): 36-47.View ArticlePubMed
- Costa V, Angelini C, De Feis I, Ciccodicola A: Uncovering the complexity of transcriptomes with RNA-Seq. J Biomed Biotechnol. 2010, 2010: 853916.PubMed CentralView ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.