De novo assembly of the common marmoset transcriptome from NextGen mRNA sequences
© Maudhoo et al.; licensee BioMed Central Ltd. 2014
Received: 15 August 2014
Accepted: 4 September 2014
Published: 19 September 2014
Nonhuman primates are important for both biomedical studies and understanding human evolution. Although research in these areas has mostly focused on Old World primates, such as the rhesus macaque, the common marmoset (Callithrix jacchus), a New World primate, offers important advantages in comparison to other primates, such as an accelerated lifespan. To conduct Next Generation expression studies or to study primate evolution, a high quality annotation of the marmoset genome is required. The availability of marmoset transcriptome data from five tissues, including both raw sequences and assembled transcripts, will aid in the annotation of the newly released marmoset assembly.
RNA was extracted from five tissues: skeletal muscle, bladder and hippocampus from a male common marmoset, and cerebral cortex and cerebellum from a female common marmoset. All five RNA samples were sequenced on the Illumina HiSeq 2000 platform. Sequences were deposited in the NCBI Sequence Read Archive. Transcripts were assembled, annotated and deposited in the NCBI Transcriptome Shotgun Assembly database.
We have provided a high quality annotation of 51,163 transcripts with full-length coding sequence. This set represented a total of 10,833 unique genes. In addition to providing empirical support for the existence of these 10,833 genes, we also provide sequence information for 2,422 genes that were not previously identified in the Ensembl annotation of the marmoset genome.
KeywordsCallithrix jacchus Common marmoset Transcriptome mRNA-seq Assembly
Nonhuman primates have, and will continue to be, important translational biomedical models for normative biological and behavioral phenomena, and for exploring atypical and/or pathological states relevant to human health and wellbeing . In the past, Old World primates (e.g., rhesus macaques) and hominoid primates (e.g., chimpanzees) have been widely utilized by the biomedical community as the primary nonhuman primate models, justified primarily on the basis of the close phylogenetic relatedness of these species to human beings, and hence high homology in biological and behavioral traits. However, Old World and hominoid primates are long-lived and costly to maintain, making the use of these primates in longitudinal research problematic from both scientific and financial perspectives.
The New World marmoset monkey (Order Primates, Parvorder Platyrrhini, Family Cebidae, Subfamily Callitrichinae) constitutes an important alternative model to Old World and hominoid primates. From a management and safety perspective, marmosets are easier and less expensive to maintain and pose substantially lower potential for zoonoses than Old World and hominoid primates . These small-bodied monkeys have accelerated life histories relative to other primates, reaching reproductive age at approximately 12–18 months of age, and an average captive lifespan of seven years, thus making developmental and longitudinal studies more practicable . The marmoset has proven utility as a model in a wide range of biomedical research, including immunology, reproduction, neuroscience, and developmental biology (reviewed in [2–4]). Furthermore, genomic research with this species has already yielded important insights [5, 6].
A draft version of the marmoset (Callithrix jacchus) assembly has recently been published . To conduct Next Generation gene expression studies and to study primate evolution, a high quality annotation of the marmoset genome is required. Here, we present marmoset transcriptome data from five tissues, including 30.2 Gb of raw sequences and 209,966 assembled transcripts representing the full length CDS of 10,833 unique genes. These resources should be helpful in the annotation of the newly-released marmoset assembly.
Two 2-year-old common marmosets (Callithrix jacchus), (animal number 2218, female and animal number 2219, male) were euthanized due to gastrointestinal illness and associated weight loss. Marmosets were humanely euthanized with appropriate veterinary supervision according to procedures approved by the UNO/UNMC Institutional Animal Care and Use Committee (#12-099-12). The animals were initially anesthetized with gaseous isoflurane, and then received a lethal overdose of sodium pentobarbital (i.v., 150 mg/kg). Immediately upon deep sedation, tissues were collected, placed in cryovials and immediately frozen in liquid nitrogen. The cryovials were stored at -80°C. Tissues included: skeletal muscle, bladder and cerebral hippocampus (from animal number 2218) and cerebellum and cerebral cortex (from animal number 2219).
Frozen tissues were thawed and total RNA was extracted with TRIzol and purified with an RNeasy Minikit (Qiagen). Extracted RNA was stored at -80°C until analysis. RNA extraction quality was assessed with an Agilent 2100 Bioanalyzer using the RNA 6000 Nano Kit (Agilent). All samples used had RNA Integrity Numbers greater than 8.0.
Sequence and assembled transcript accessions
Reads were filtered to remove genomic contamination. To accomplish this, the reads were aligned with human RefSeq (version 57) mRNA sequences using BLAST (BLAST + v2.2.25) . All reads with a reported alignment length of less than 87 and greater than 115 nucleotides were removed from the input file. We also removed a read if its mate was removed.
Submitted contig statistics
This strategy was designed to avoid false positives. However, because we used human transcripts to identify orthologs, we will likely miss annotating potential marmoset-specific genes.
We used our assembled marmoset transcripts as queries with BLASTx  (BLAST + v2.2.25) to identify orthologous human RefSeq proteins. We used the full-length coding sequence (CDS) ranges obtained from the BLASTx output files to derive putative marmoset protein sequences. Only transcripts with both start and stop codons were annotated. In addition, only transcripts with conceptually derived proteins with alignment gap lengths of less than or equal to 15, protein length differences of less than or equal to 15, and protein identities of greater than or equal to 70% with respect to their human orthologs were annotated. Both annotated and unannotated transcripts were deposited in the TSA database (Table 1).
A total of 209,966 transcripts were submitted to NCBI. From this group, we annotated 51,163 transcripts with full-length CDS. This set represented a total of 10,833 unique genes. Conceptual translations (derived from BLASTx) had an average identity and similarity of 96.56% and 97.92% to their human orthologs, respectively.
Comparison with Ensembl annotations
Availability of supporting data
The datasets supporting this article are available at the National Center for Biotechnology Information (NCBI) accessioned under Bioproject ID: 203643. Illumina sequences were deposited in the SRA under accessions [SRX285593, SRX285538, SRX285591, SRX285592, SRX285594]. Assembled transcripts were deposited in the Transcriptome Shotgun Assembly (TSA) database under accessions [GAMQ01000001-GAMQ01033528, AMT01000001-GAMT01041104, GAMS01000001-GAMS01044498, GAMR01000001-GAMR01044112, GAMP01000001-GAMP01046724]. Supporting BLASTx results, and ISA-TAB metadata are also available from the GigaScience database .
Basic local alignment search tool
Sequence read archive
Transcriptome shotgun assembly.
We thank Alok Dhar at the UNMC DNA Sequencing Core facility for library preparation and Illumina sequencing and Sanjit Pandey for bioinformatics support. This research was supported, in part, by funds from the Nebraska Research Initiative awarded to RBN. Partial support for the marmoset colony housed at the Callitrichid Research Center at UNO derives from NIH HD042882 awarded to JAF.
- Phillips KA, Bales KL, Capitanio JP, Conley A, Czoty PW, ‘t Hart BA, Hopkins WD, Hu SL, Miller LA, Nader MA, Nathanielsz PW, Rogers J, Shively CA, Voytko ML: Why primate models matter. Am J Primatol. 2014, 76: 801-827. 10.1002/ajp.22281.View ArticlePubMedPubMed CentralGoogle Scholar
- Abbott DH, Barnett DK, Colman RJ, Yamamoto ME, Schultz-Darken NJ: Aspects of common marmoset basic biology and life history important for biomedical research. Comp Med. 2003, 53: 339-350.PubMedGoogle Scholar
- Tardif SD, Araujo A, Arruda MF, French JA, Sousa MB, Yamamoto ME: Reproduction and aging in marmosets and tamarins. Interdiscip Top Gerontol. 2008, 36: 29-48.View ArticlePubMedGoogle Scholar
- Mansfield K: Marmoset models commonly used in biomedical research. Comp Med. 2003, 53: 383-392.PubMedGoogle Scholar
- Ward JM, Vallender EJ: The resurgence and genetic implications of New World primates in biomedical research. Trends Genet. 2012, 28: 586-591. 10.1016/j.tig.2012.09.003.View ArticlePubMedPubMed CentralGoogle Scholar
- The Marmoset Genome Sequencing and Analysis Consortium: The common marmoset genome provides insight into primate biology and evolution. Nature Genet. 2014, 46: 850-857. 10.1038/ng.3042.View ArticlePubMed CentralGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1016/S0022-2836(05)80360-2.View ArticlePubMedGoogle Scholar
- Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18: 821-829. 10.1101/gr.074492.107.View ArticlePubMedPubMed CentralGoogle Scholar
- Schulz MH, Zerbino DR, Vingron M, Birney E: Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 2012, 28: 1086-1092. 10.1093/bioinformatics/bts094.View ArticlePubMedPubMed CentralGoogle Scholar
- Gish W, States DJ: Identification of protein coding regions by database similarity search. Nature Genet. 1993, 3: 266-272. 10.1038/ng0393-266.View ArticlePubMedGoogle Scholar
- Maudhoo MD, Ren D, Gradnigo JS, Gibbs RM, Lubker AC, Moriyama EN, French JA, Norgren RB: Supporting materials for the “De novo assembly of the common marmoset transcriptome from NextGen mRNA sequences”. GigaScience Database. GigaScience Database. 2014,http://dx.doi.org/10.5524/100103,Google Scholar
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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.