Figures
Abstract
Background
Cephalopoda are a class of Mollusca species found in all the world's oceans. They are an important model organism in neurobiology. Unfortunately, the lack of neuronal molecular sequences, such as ESTs, transcriptomic or genomic information, has limited the development of molecular neurobiology research in this unique model organism.
Results
With high-throughput Illumina Solexa sequencing technology, we have generated 59,859 high quality sequences from 12,918,391 paired-end reads. Using BLASTx/BLASTn, 12,227 contigs have blast hits in the Swissprot, NR protein database and NT nucleotide database with E-value cutoff 1e−5. The comparison between the Octopus vulgaris central nervous system (CNS) library and the Aplysia californica/Lymnaea stagnalis CNS ESTs library yielded 5.93%/13.45% of O. vulgaris sequences with significant matches (1e−5) using BLASTn/tBLASTx. Meanwhile the hit percentage of the recently published Schistocerca gregaria, Tilapia or Hirudo medicinalis CNS library to the O. vulgaris CNS library is 21.03%–46.19%. We constructed the Phylogenetic tree using two genes related to CNS function, Synaptotagmin-7 and Synaptophysin. Lastly, we demonstrated that O. vulgaris may have a vertebrate-like Blood-Brain Barrier based on bioinformatic analysis.
Conclusion
This study provides a mass of molecular information that will contribute to further molecular biology research on O. vulgaris. In our presentation of the first CNS transcriptome analysis of O. vulgaris, we hope to accelerate the study of functional molecular neurobiology and comparative evolutionary biology.
Citation: Zhang X, Mao Y, Huang Z, Qu M, Chen J, Ding S, et al. (2012) Transcriptome Analysis of the Octopus vulgaris Central Nervous System. PLoS ONE 7(6): e40320. https://doi.org/10.1371/journal.pone.0040320
Editor: Christian Schönbach, Kyushu Institute of Technology, Japan
Received: December 27, 2011; Accepted: June 4, 2012; Published: June 29, 2012
Copyright: © 2012 Zhang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by the National fund for oceanography research in Public Interest (201005013) and the National Key Technology R&D Program (2011BAD13). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Cephalopods are a group of unique Mollusk species that evolved in the Late Cambrian. There are around 800 extant species in Cephalopoda, none of which live in fresh water. Cephalopods possess many primitive traits of invertebrates (the genital chamber is connected with the body cavity; there is no kidney in the early embryonic development period, and the genital duct derives from the coelomic duct), but also exhibit more advanced features such as the vertebrate-like eye, no larval phase in ontogenesis, and a highly centralized nervous system. Cephalopods' hybrid evolutionary characteristics suggest that the Cephalopoda are a special branch of Mollusca. The Cephalopods' nervous system is the most complex of the invertebrates. This has led to its use in research. For example, the squid giant axon was used to prove the mechanism of action potentials by John Carew Eccles, Alan Lloyd Hodgkin and Andrew Fielding Huxley, for which they won the 1963 Nobel Prize. Recently, the use of Cephalopods in neurobiology research has been limited. In contrast, molecular neurobiology research using high-throughput transcriptome sequencing technology in other Mollusk species has been increasing rapidly. Leonid L Moroz established a central nervous system (CNS) database of Aplysia californica for gene discovery, expression profiling, and characterization of signaling pathways in 2006 [1]. In 2009, Z.P. Feng and Z. Zhang provided a Lymnaea stagnalis CNS EST database and used it to illustrate the substantial genetic diversity within A. californica [2].
Through the evolutionary process, the structure and function of the central nervous system has become more complex. Some of these adaptations are important for physical protection of the CNS. The Blood-Brain Barrier is one example. It has been demonstrated to be an important structure for maintaining the basic stability of the internal environment. In 1885, Paul Ehrlich found that vital dyes can stain all mammalian organs except the brain. Edwin Goldmann demonstrated that if dye were injected into the cerebro-spinal fluid, then the central nervous system (CNS) would be stained, but no other organs. These results suggest that there exists a barrier between the CNS and the peripheral organs. Blood-Brain Barrier function is controlled by the brain's microvascular endothelial cells (BMVEC), basement membrane, and four different types of neighboring glial cells such as astrocytes, perivascular pericytes and microglia. Neurons are also important for the function of the CNS and BBB. The brain's microvascular endothelial cells (EC), and cellular connection between the BMVEC were demonstrated to play a key role in BBB function [3]. EC supply high electric impedance, specific receptors and transporters, while the specialized junctions eliminate gaps between EC and prevent free substances interchange between the blood and brain parenchymal space [4]–[7]. Research related to the BBB are not only in vertebrates; in invertebrates such as Drosophila there also exists a barrier system to protect the CNS, albeit a different one [8]. It seems that the barrier system was necessary in organisms with a more active lifestyle. Unlike other Mollusk species, Cephalopods live a more active life. Cephalopods also have a CNS structure and closed vascular system similar to vertebrates. As of yet, there is no direct evidence that Cephalopods have a vertebrate-like BBB. Discovery of a vertebrate-like BBB would accelerate the understanding of BBB structure, function and origin.
O. vulgaris is one of the commercially valuable benthic species among Cephalopods. It is widely distributed in the Pacific, Atlantic and Indian Oceans, as well as the Mediterranean Sea. It belongs to the phylum Mollusca (Cephalopoda, Dibranchiata, Octopoda, Octopodidae). The economic value and wide distribution make it easy to sample, and aid in its use as a model organism. In recent years, most O. vulgaris research has been focused on larval growth [9]–[12], physiology [13]–[16], ecological toxicology [17]–[22] and fisheries [23]–[26]. Due to a lack of molecular sequence information, the development of molecular biology research in O. vulgaris is limited.
In this study, we dissect a complete adult O. vulgaris for CNS tissue to prepare the first non-normalized cDNA library. We use high-throughput Illumina Solexa sequencing technology to acquire a mass of CNS molecular information about the Cephalopod O. vulgaris. After assembly, gene name annotation and GO/KEGG annotation, this transcriptome data was used to develop a better understanding of O. vulgaris, such as a molecular foundation of its biological characteristics and phylogenetic analysis of O. vulgaris versus other published CNS transcriptome data. It was also a great resource for mining information about genes related to the CNS and vertebrate-like Blood-Brain Barrier. We present the first CNS transcriptome analysis of Cephalopoda species O. vulgaris aimed to accelerate the study of functional molecular neurobiology and comparative evolutionary biology in Cephalopoda.
Results and Discussion
1. Illumina sequencing and sequence assembly
We dissected a complete adult O. vulgaris for CNS, including cerebral ganglion, visceral ganglion, pedal ganglion, and part of the peripheral nervous system (Figure 1). The non-normalized cDNA library is constructed following the manufacturer's instructions (Illumina). High-throughput SOLEXA paired-end sequencing yielded a total of 13,753,396 reads with length of 90 bp and Q20 percentage equal 96.38%. 834,465 sequences were removed because of low quality. These high-quality sequences were assembled with VELVET [27] and a summary table for assembly was deposited to Table S1. For best assembly results, four important parameters were tested, including K-mers, coverage, minimum contigs length, and numbers of contigs. 12,918,931 reads were assembled into 59,859 contigs with a maximum length of 8,970 bp and the N50 length of 450 bp. Then the CD-HIT program [28] was performed and the results showed that only 48 (<0.1%) sequences had significant similarities (>98% identity and >95% coverage) against other sequences within the dataset, indicating the assembly was accomplished and high quality. The length distribution for all contigs is presented in Figure 2. Of these, a total of 31,315 open reading frames (ORFs≥50 aa) have been detected (Figure S1).
2. BLASTX searches in Swissprot and NR protein database
Contig gene name annotation was archived through BLASTx and BLASTn searches against the Swissprot,NR (NCBI non-redundant protein database), and NT (NCBI nucleotide sequences database). Using Perl, the description of the most relevant hits with E-value less than 1e−5 were assigned to the query sequences. This search revealed that only 10,412 (17.39%) and 1,815 (3.03%) sequences have a significance blast hit (1e−5). The reason why most of the sequences (79.58%) didn't have a significant blast hit is probably caused by potentially novel genes and a lack of molecular information of closely related species in Cephalopods.
Table 1 is a list of the 20 sequences, sorted by read number, which revealed the most expressed proteins or peptides. All sequences with high read number have been annotated. This result may illustrate the vivid biological characteristics of the CNS transcriptome, including several neurohormones, neuro-consisted proteins and stress-response proteins. Some proteins may involve features specific to O. vulgaris, such as fatty acid-binding proteins and apolopophorins which reflect the strong lipid metabolism of O. vulgaris, and retinol dehydrogenase which convert retinol into retinal to maintain vision.
3. Functional annotation based on GO and KEGG analysis
In order to determine the functions of these sequences, 12,227 genes were selected to annotate with the GO database by Gominer [29]. Figure 3 shows the gene number of GO functional annotation analysis in level 2.
By selecting the defining GO term with keywords “neuro” and “nervous”, more than 50 GO terms are involved and their matching genes ranged from 1 to 338. The largest GO categorie related to CNS function is nervous system developments (GO index: 0007399). Table 2 presents the number of genes matched in 15 GO terms that are relevant to neuronal functions and have the largest numbers matching, demonstrating our transcript database has a great number of neural sequences that will provide a great resource for further research of the development, function and regulatory mechanisms of the O. vulgaris central nervous system.
KEGG (Kyoto Encyclopedia of Genes and Genomes) is a pathway-based categorization of orthologous genes that provides useful information for predicting functional profiles of genes [30]. To further demonstrate the biological pathways that are active in the CNS of O. vulgaris, 3,825 genes were mapped into the signaling pathways in the KEGG.
4. Comparisons with the previously published Mollusc species A. californica and L. stagnalis CNS database
The published CNS transcriptome database of A. californica [2] and L. stagnalis [1] give us an opportunity to understand the diversity between O. vulgaris and other Mollusk species. Tables 3 and 4 provides the data of BLASTn or tBLASTx hits between O. vulgaris and L. stagnalis (or A. californica). In the comparison between O. vulgaris and L. stagnalis ESTs (10,375 ESTs which assemble into 7,712 unique sequences) by BLASTn and tBLASTx, at E-value cutoff of 1e−5, about 7.83% (BLASTn) and 26.53% (tBLASTx) of L. stagnalis cDNAs have hits in O. vulgaris CNS library. In the mean time, approximately 1.03% (BLASTn) and 5.77% (tBLASTx) of the sequences in O. vulgaris CNS library have a hit in L. stagnalis cDNAs. Considering these two databases are constructed from CNS tissues, the main cause of low percentages of blast hits could be the difference in their genetic codes. In addition, these results may be explained by: 1) the sequences below 200 bp are unlikely to have a blast hit 2) it is possible that some genes exist in both species, but have expression levels below the detection limit of the transcriptome studies, and 3) differing expression profiles between different species or developmental stages.
The CNS database of A. californica contains two small EST databases: Normalized from Pedal-Pleural ganglia and normalized CNS library (juvenile 1), which provided 40607 sequences after assemble. 3.97% (BLASTn) and 18.00% (tBLASTx) A. californica ESTs have a hit in our transcriptome data, while 1.90% (BLASTn) and 12.85% (tBLASTx) of O. vulgaris sequences have a significant match in A. californica. The low similarity may be due to: 1) the huge diversity of the two species; 2) the CNS database of A. californica contains both adult and juvenile CNS tissue (the expression of genes may vary at different growth stages); and 3) though the two parts of A. californica database are normalized, some sequences may exist in both databases that significantly downgrade the hit percentage of O. vulgaris cDNAs to the A. californica ESTs database.
5. Comparison with CNS sequence databases of other model organisms
Meanwhile, we chose three different CNS databases belonging to a fish (Tilapia) [31], an arthropod (Schistocerca gregaria) [32], and an annelid (Hirudo medicinalis) [33], all of which were recently published in NCBI, to discuss the diversity among these model organisms in different phylums. We ran BLASTn/tBLASTx searches against the CNS database of fish, dessert locust, and leech. Results are shown in Table 5. With E-value cutoff 1e−5, 3544 cDNAs have a hit in Tilapia, the number of hit in other organisms are 5870(Locust) and 6568(Leech). Considering the scale of these CNS databases and the Mollusk species L. stagnalis EST database, these hit percentages suggest that O. vulgaris has similarity with all organisms except Annelid (H. medicinalis). This illustrates that O. vulgaris may have both original and evolution characters at the transcriptome level, and this comparison will be helpful for further research in homologous or non-homologous genes of O. vulgaris.
Figure 4 shows the distribution of sequence hits in the locust, fish and leech CNS database. 8,695 sequences match the condition and only 2,213 (25.45%) sequences exist in both 3 databases that are much lower than the same comparison in post studies, which indicate the low percentage of highly conserved genes. This comparison reflects O. vulgaris' large number of homologous genes in both vertebrates and invertebrates, indicating that it may be a good model organism for studying neurobiology.
6. Phylogenetic analysis of two CNS-related genes
We identified genes related to nervous system functions that existed in the four databases (Leech, Mollusks, Locust and Fish) by applying a Perl script to select the genes meeting the conditions of length (>1,000 bp), E-value (<1e−10), and coverage (>30). Then we built the corresponding phylogenetic trees. The tree in Figure 5 is generated by synaptotagmin-7; this protein is a member in the synaptotagmin family believed to be important in the docking and fusion of synaptic vesicles with the plasma membrane, such as neurotransmitter release [34], [35]. We found the sequence of O. vulgaris is most related to the Mollusk species A. californica. Synaptotagmin-7 widely exists in metazoan and the phylogenetic tree performs well in sorting different categories of animal, indicating that it is a good choice for explaining system evolution. Furthermore, the synaptotagmin-7 of Mollusks is more closely related to arthropods and fish that consistently use the results in comparison to the CNS database between O. vulgaris and Locust/Tilapia.
We further chose synaptophysin because they were found in A. californica and L. stagnalis and are correspond to the phylogenetic tree (Figure 6). Synaptophysin is a protein that acts as a marker for neuroendocrine tumors and quantification of synapses [36]. Even without an ortholog gene in A. californica or L. stagnalis, the result still reveals the closer relation of O. vulgaris to vertebrates.
7. A putative vertebrate-like Blood-Brain Barrier in O. vulgaris
Cephalopods have a well-developed nervous system, reflected in its precise internal structure and clear healing between ganglions covered with cerebral cortex. In addition, around the CNS there is a cartilaginous skull for protection. These characteristics are similar to vertebrates, indicating that the Cephalopods may be an evolutionary transition to the brain functions of vertebrates. By studying these transitions, we can clearly understand the process of its occurrence. The Blood-Brain Barrier is a kind of internal barrier system related to internal immunity that blocks pathogenic microorganisms and other macromolecules through the blood circulation into brain tissue to maintain the basic stability of the internal environment. It also has the important biological role of maintaining the normal physiological state of the central nervous system [3]. The vertebrate Blood-Brain Barrier has three histological bases: brain microvascular endothelial cells (BMVEC) and junctions between BMVEC, continuous basement membrane around the BMVEC and five different types of neighboring glial cells such as astrocytes, perivascular pericytes, microglia, and surrounding neurons. In addition, to ensure that central neurotransmitters are not able to pass the BBB and maintain stability of neurotransmitter concentrations, BMEVC has a unique enzyme system to inactivate the central neurotransmitters such as monoamine oxidase [37], AAAD, and COMT. We hope to interrogate our CNS dataset and compare it with three other invertebrate model organisms, to find whether Cephalopods have specific gene expression indicative of a vertebrate-like BBB.
The related genes involving junctions between EC, pathways across BBB, and the enzyme barrier system were selected for clarifying whether the molecular foundation of vertebrate-like BBB exists in four invertebrates (O. vulgaris, A. californica, S. gregaria, H. medicinalis). First, we focused on all of the important genes involved in tight junctions (TJ) and adherence junctions (AJ) structure such as Claudins [4], [38], Occludin, Junctional adhesion molecules(JAM) [39], Cytoplasmatic proteins [40]–[42] and Cadherins [43], Catenins [44]. At the beginning we searched the gene name annotation result of O. vulgaris, two TJ related genes and two AJ related genes were found when E-value<e−10. In addition, we carried out a tBLASTn search by using the amino sequences of the TJ and AJ related genes downloaded from NCBI against all CNS Datasets. These results were filtered with a Perl script including a number of strict conditions to ensure accuracy: E-value<e−10, the number of amino acids that align to the query sequences >100 and covered more than 80% of subject sequences. After two screening steps, O. vulgaris cDNA hit nine of twelve proteins. Considering the alignment features of BLAST and that some proteins have a low conservation between different species, these results indicate that the number of actual proteins will be more than its hit number. This result reveals that O. vulgaris has complete BMVEC junctions. After the same filter steps, A. californica cDNA hit only five proteins, S. gregaria and H. medicinalis hit four illustrated that the rest of organisms do not have a significant matching in proteins involving junctions between EC (Table 6).
Besides specific junctions, the elaborate systems for transporting macromolecules and elimination of neurotransmitter are important part of BBB functions [4], [45]. In O. vulgaris we found all of the transporters and enzymes that were well-studied exist after the filter steps, except the Caveolae [46], [47]. While in A. californica and S. gregaria, the miss of glucose transporter-1 [37], [48] reveals the absence of vertebrate-like BBB in these two species. It is impossible for the CNS to function with no energy supply when a physical barrier is established [41], [47]. Another result from exploring enzymes such as monoamine oxidase [37], COMT [49], and AAAD [50] also suggests that O. vulgaris has a thorough enzyme system for eliminating the cause of neurotransmitter feedback on the central nervous system. The other organisms are keeping their poor performance in searching the specific transporters and enzymes (Table 7).
Furthermore,we investigated the profiles of genes involved in TJ signaling pathway, which is considered to responsible for the barrier properties, between O. vulgaris and A. californica (Figure S2). In default parameters, the number of genes that mapped into the TJ signaling pathway were different (V:47, A:32). Meanwhile this comparative analysis presented three important unique genes in O. vulgaris: two transmembrane proteins (JAM and Claudins) and a cytoplasmic TJ accessory protein ZO-1. Claudins are considered to be responsible for permeability restriction in TJ [51] and JAM is involved in various of TJ function such as cell-to-cell adhesion, organizing structure, taking part in the formation of TJ as an integral membrane protein together with claudins [43]. ZO-1 existed as a carboxy-terminal region, which binds to actin and links the TJ to the cytoskeleton, acts as a central organizer of the TJ complex [52]. This result clearly implied that O. vulgaris have a more complex and integrated TJ functions than the model organism A. californica.
To verify the constitutive expression of the genes corresponding to the BBB structure and function, specific primers were designed based on assembled contigs and quantitative real-time PCR were performed. All of the genes showed a ubiquitous expression in all examined tissues, including brain, liver, heart, gill and muscle (Figure S3). The transcription pattern of these BBB relevant genes had been determined in vertebrates and different kinds of tissues indicated that these genes may not only be involved in BBB but also participated in other physiological functions. For example, ZO-1 mostly expressed in endothelial and epithelial cells forming the TJ assembly, but it still expressed in other tissues not forming TJ that may be involved in signal transduction at cell-cell junctions [43]. Glucose transporter 1 expressed in erythrocytes and also in the endothelial cells of barrier tissues, it also has been identified in muscle, fat and tissues with acute insulin-stimulated glucose transport [53]. Monoamine oxidases are found in neurons, astroglia and also found in the liver, gastrointestinal tract, and placenta that catalyze the oxidative deamination of monoamines [54]. The same results that are similar to above were displayed in other genes related to TJ (myosin [55]), AJ (α-Catenin and β-Catenin [44], Cadherins [56]), specific transporters (P-glyprptein [57], Multidrug resistance-associated protein [4], Organic anion transporter [58]), and enzymes (Catechol-O-methyltransferase [58], Aromatic L-amino acid decarboxylase [59]). These observations not only implied the accuracy of contig assembly, but also demonstrated that as vertebrates, all the target genes can be expressed in different tissues including CNS in O. vulgaris.
It is certainly to be noted that the CNS of O. vulgaris has a large number of proteins involved in specific junctions, transporters, and enzymes which are definitely indispensable to form an incredible system that may possesses most vertebrate BBB functions. The results in A. californica and S. gregaria indicate that the species with open vascular system may have a different strategy for protecting the basic stability of the internal environment like D. melanogaster. The low hit percentage in H. medicinalis demonstrated that although H. medicinalis have a similar circulatory system, the loose structure of central nervous system limited the development of vertebrate-like BBB system. Based on the results above, only O. vulgaris has the molecular basis of the vertebrate-like BBB, highlighting its use as a model organism for the in-depth study of phylogenetics, structure and function of the BBB.
Materials and Methods
1. Preparation of central nervous system samples
O. vulgaris (female) used in this study were obtained from the fish market in Xiamen, Fujian province. Before dissection they were maintained in aquaria. The CNS was dissected under Zeiss steRED Lumar V12 dissection system, including cerebral ganglion, visceral ganglion, pedal ganglion, and part of peripheral nerve. These were isolated and prepared for total RNA extraction.
2. RNA isolation, library preparation and sequencing
Using RNAiso Plus (TaKaRa), the total RNA was extracted in accordance with the manufacturer's instructions. The quality of the isolated RNA was checked by electrophoresis on a 1.5% agarose gel and by absorption spectroscopy. The library was constructed using the Mate Pair Library Preparation Kit. We enriched the mRNA by using the magnetic beads which contain Oligo (dT). While filled in the fragmentation buffer, we randomly broke mRNA into small pieces. Random hexamers were used for the synthesis of the first strand of cDNA. The complementary strand was synthesized. With several steps of purification, adaptor addition and cDNA length selection, this library was used for sequencing on an Illumina HiSeq 2000. A paired-end Solexa sequencing strategy was used for sequencing to better understand the gene expression profile and for convenient assembly of the entire transcriptome de novo.
3. Transcript assembly and annotation
Transcripts were assembled using VELVET software, 59,859 contigs were generated and 31,909 contigs with length ≥200 bp were submitted to the TSA (Transcriptome Shotgun Assembly Sequences Database) (Genbank TSA Acc. Number: JR435555–JR467463). Contigs with length <200 bp were deposited in table S2. To annotate these sequences we aligned them using the Swissprot, NR protein databases and NT nucleotide database. Sequences producing an E-value<1e−5 were considered a hit. Read number was determined based on the relationship between contig length and contig coverage. Functional GO annotations were performed using GOminer and KEGG pathway analysis was carried out using KASS (KEGG Automatic Annotation Server) with default parameters and 5 model organisms in vertebrate (Homo sapiens, Mus musculus, Gallus gallus, Xenopus laevis, Danio rerio) that were chosen as genes data set.
4. Comparison to post published CNS databases and phylogenetic analysis
BLASTn and tBLASTx were carried out for testing the similarity between O. vulgaris and other model organisms' CNS database (A. californica, L. stagnalis, S. gregaria, H. medicinalis, Tilapia) at both the nucleotide and amino acid level with a threshold E-value less than 1e−5. Phylogenetic analysis was carried out using the sequences that had been annotated and met the conditions of sequence length (>1,000 bp), E-value (<1e−10) and coverage (>30). Then the sequences were selected for building the phylogenetic tree using the neighbor-joining method.
5. Identification of genes involved in the vertebrate-like Blood-Brain Barrier
To identify the putative genes involved in the vertebrate-like Blood-Brain Barrier, the sequences from this study and published invertebrate CNS database were screened using confirmed protein sequences in NCBI, which had been demonstrated to have involved specific structures or functions of vertebrate BBB. tBLASTn was used for scanning and a Perl script for screening the blast results with E-value<1e−10. The number of amino acids that align to the query sequences>100 and covered more than 80% of subject sequence ensure that the alignment between sequences are of global significance.
6. Gene verification by Real-time PCR
The expression of 14 target genes related to Blood-Brain Barrier were Validated by quantitative Real-time PCR (qRT-PCR) using in ABI 7900HT Real-time Detection System (Applied Biosystems, USA) with SYBR® Premix DimerEraser (TaKaRa, Japan), and the house-keeping gene β-actin was used as internal standard gene. Primers for all target genes as well asβ-actin gene were listed in Table 8.
Total RNA was prepared from the tissues (brain, liver, heart, gill and muscle), and the first-strand cDNA synthesis was the same as described in section 1. PCR amplifications were performed in triplicate wells under following conditions: initial denaturation at 95°C for 30 s followed by 40 cycles of 95°C 5 s, 55°C for 30 s and 72°C 10 s. Dissociation analysis of amplification products was performed at the end of the PCR reaction. After the PCR program, data were analyzed with ABI 7900HT SDS software 2.3 (Applied Biosystems, USA). To maintain consistency, the baseline was set automatically by the software. The cycle threshold values were converted into the equivalent target amount using the external calibration curves (correlation coefficient >0.99), which were based on 15 standard targets of recombinant PMD19-T containing 14 target genes and β-actin gene fragment, respectively. All data were given in terms of relative mRNA expressed as means±SE.
Supporting Information
Figure S1.
ORF distribution of contigs obtained from O. vulgaris central nervous system (CNS) transcriptome library.
https://doi.org/10.1371/journal.pone.0040320.s001
(TIF)
Figure S2.
The tight junction pathway in O. vulgaris and A. californica .
https://doi.org/10.1371/journal.pone.0040320.s002
(RAR)
Figure S3.
Tissue distribution of 14 BBB-related genes.
https://doi.org/10.1371/journal.pone.0040320.s003
(RAR)
Table S2.
Contigs with length less than 200 bp.
https://doi.org/10.1371/journal.pone.0040320.s005
(TXT)
Acknowledgments
The authors gratefully thank Shuibo Chen for obtaining experimental animals; Ji Qing and Biyun Feng for assisting at the microdissection of O. vulgaris CNS tissues; Yunyun Zhuang for technical support. In addition all collaborators who were involved in O. vulgaris cultures are acknowledged.
Author Contributions
Conceived and designed the experiments: XZ YM JH SD. Performed the experiments: XZ MQ JH TS. Analyzed the data: XZ ZH MQ JC SD. Contributed reagents/materials/analysis tools: XZ YM ZH JC JH. Wrote the paper: XZ MQ JC SD. Wrote the PERL scripts used in analysis: XZ ZH JC.
References
- 1. Feng ZP, Zhang Z, van Kesteren RE, Straub VA, van Nierop P, et al. (2009) Transcriptome analysis of the central nervous system of the mollusc Lymnaea stagnalis. BMC Genomics 10: 451.
- 2. Moroz LL, Edwards JR, Puthanveettil SV, Kohn AB, Ha T, et al. (2006) Neuronal Transcriptome of Aplysia: Neuronal Compartments and Circuitry. Cell 127: 1453–1467.
- 3. Cardoso FL, Brites D, Brito MA (2010) Looking at the blood–brain barrier: Molecular anatomy and possible investigation approaches. Brain Research Reviews 64: 328–363.
- 4. Zlokovic BV (2008) The blood-brain barrier in health and chronic neurodegenerative disorders. Neuron 57: 178–201.
- 5. Zheng W, Aschner M, Ghersi-Egea J-F (2003) Brain barrier systems: a new frontier in metal neurotoxicological research. Toxicology and Applied Pharmacology 192: 1–11.
- 6.
Weksler BB (2005) Blood-brain barrier-specific properties of a human adult brain endothelial cell line. The FASEB Journal.
- 7. Hawkins RA, O'Kane RL, Simpson IA, Vina JR (2006) Structure of the blood-brain barrier and its role in the transport of amino acids. J Nutr 136: 218S–226S.
- 8. Schwabe T, Bainton RJ, Fetter RD, Heberlein U, Gaul U (2005) GPCR Signaling Is Required for Blood-Brain Barrier Formation in Drosophila. Cell 123: 133–144.
- 9. Quintana D, Domingues P, García S (2008) Effect of two artificial wet diets agglutinated with gelatin on feed and growth performance of common octopus (Octopus vulgaris) sub-adults. Aquaculture 280: 161–164.
- 10. Seixas P, Rey-Méndez M, Valente LMP, Otero A (2008) Producing juvenile Artemia as prey for Octopus vulgaris paralarvae with different microalgal species of controlled biochemical composition. Aquaculture 283: 83–91.
- 11. Villanueva R, Escudero JM, Deulofeu R, Bozzano A, Casoliva C (2009) Vitamin A and E content in early stages of cephalopods and their dietary effects in Octopus vulgaris paralarvae. Aquaculture 286: 277–282.
- 12. Seixas P, Rey-Méndez M, Valente LMP, Otero A (2010) High DHA content in Artemia is ineffective to improve Octopus vulgaris paralarvae rearing. Aquaculture 300: 156–162.
- 13. Hormiga JA, Almansa E, Sykes AV, Torres NV (2010) Model based optimization of feeding regimens in aquaculture: Application to the improvement of Octopus vulgaris viability in captivity. Journal of Biotechnology 149: 209–214.
- 14. Shomrat T, Feinstein N, Klein M, Hochner B (2010) Serotonin is a facilitatory neuromodulator of synaptic transmission and “reinforces” long-term potentiation induction in the vertical lobe of Octopus vulgaris. Neuroscience 169: 52–64.
- 15. Gutnick T, Byrne Ruth A, Hochner B, Kuba M (2011) Octopus vulgaris Uses Visual Information to Determine the Location of Its Arm. Current Biology 21: 460–462.
- 16. Prato E, Portacci G, Biandolino F (2010) Effect of diet on growth performance, feed efficiency and nutritional composition of Octopus vulgaris. Aquaculture 309: 203–211.
- 17. Raimundo J, Vale C, Canário J, Branco V, Moura I (2010) Relations between mercury, methyl-mercury and selenium in tissues of Octopus vulgaris from the Portuguese Coast. Environmental Pollution 158: 2094–2100.
- 18. Raimundo J, Vale C, Duarte R, Moura I (2010) Association of Zn, Cu, Cd and Pb with protein fractions and sub-cellular partitioning in the digestive gland of Octopus vulgaris living in habitats with different metal levels. Chemosphere 81: 1314–1319.
- 19. Raimundo J, Costa PM, Vale C, Costa MH, Moura I (2010) DNA damage and metal accumulation in four tissues of feral Octopus vulgaris from two coastal areas in Portugal. Ecotoxicology and Environmental Safety 73: 1543–1547.
- 20. Raimundo J, Costa PM, Vale C, Costa MH, Moura I (2010) Metallothioneins and trace elements in digestive gland, gills, kidney and gonads of Octopus vulgaris. Comparative Biochemistry and Physiology Part C: Toxicology & Pharmacology 152: 139–146.
- 21. Morita T, Otosaka S, Fujimoto K, Nishiuchi K, Kimoto K, et al. (2010) Detection and temporal variation of 60Co in the digestive glands of the common octopus, Octopus vulgaris, in the East China Sea. Marine Pollution Bulletin 60: 1193–1199.
- 22. Monteiro A, Costa PR (2011) Distribution and selective elimination of paralytic shellfish toxins in different tissues of Octopus vulgaris. Harmful Algae 10: 732–737.
- 23. Vargas-Yáñez M, Moya F, García-Martínez M, Rey J, González M, et al. (2009) Relationships between Octopus vulgaris landings and environmental factors in the northern Alboran Sea (Southwestern Mediterranean). Fisheries Research 99: 159–167.
- 24. Perales-Raya C, Bartolomé A, García-Santamaría MT, Pascual-Alayón P, Almansa E (2010) Age estimation obtained from analysis of octopus (Octopus vulgaris Cuvier, 1797) beaks: Improvements and comparisons. Fisheries Research 106: 171–176.
- 25. Sobrino I, Juarez A, Rey J, Romero Z, Baro J (2011) Description of the clay pot fishery in the Gulf of Cadiz (SW Spain) for Octopus vulgaris: Selectivity and exploitation pattern. Fisheries Research 108: 283–290.
- 26. Vázquez-Rowe I, Moreira MT, Feijoo G (2012) Environmental assessment of frozen common octopus (Octopus vulgaris) captured by Spanish fishing vessels in the Mauritanian EEZ. Marine Policy 36: 180–188.
- 27. Zerbino DR, Birney E (2008) Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res 18: 821–829.
- 28. Li W, Godzik A (2006) Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22: 1658–1659.
- 29. Zeeberg BR, Feng W, Wang G, Wang MD, Fojo AT, et al. (2003) GoMiner: a resource for biological interpretation of genomic and proteomic data. Genome Biol 4: R28.
- 30. Kanehisa M, Goto S (2000) KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 28: 27–30.
- 31. Lee B-Y, Howe AE, Conte MA, D'Cotta H, Pepey E, et al. (2010) An EST resource for tilapia based on 17 normalized libraries and assembly of 116,899 sequence tags. BMC Genomics 11: 278.
- 32. Zheng J, Badisco L, Huybrechts J, Simonet G, Verlinden H, et al. (2011) Transcriptome Analysis of the Desert Locust Central Nervous System: Production and Annotation of a Schistocerca gregaria EST Database. PLoS ONE 6: e17274.
- 33. Macagno ER, Gaasterland T, Edsall L, Bafna V, Soares MB, et al. (2010) Construction of a medicinal leech transcriptome database and its application to the identification of leech homologs of neural and innate immune genes. BMC Genomics 11: 407.
- 34. Sutton RB, Davletov BA, Berghuis AM, Sudhof TC, Sprang SR (1995) Structure of the first C2 domain of synaptotagmin I: a novel Ca2+/phospholipid-binding fold. Cell 80: 929–938.
- 35. Fukuda M, Aruga J, Niinobe M, Aimoto S, Mikoshiba K (1994) Inositol-1,3,4,5-tetrakisphosphate binding to C2B domain of IP4BP/synaptotagmin II. J Biol Chem 269: 29206–29211.
- 36. Calhoun ME, Jucker M, Martin LJ, Thinakaran G, Price DL, et al. (1996) Comparative evaluation of synaptophysin-based methods for quantification of synapses. J Neurocytol 25: 821–828.
- 37. Kalaria RN, Harik SI (1987) Blood-brain barrier monoamine oxidase: enzyme characterization in cerebral microvessels and other tissues from six mammalian species, including human. J Neurochem 49: 856–864.
- 38. Ueno M (2007) Molecular anatomy of the brain endothelial barrier: an overview of the distributional features. Curr Med Chem 14: 1199–1206.
- 39. Wolburg H, Lippoldt A (2002) Tight junctions of the blood-brain barrier: development, composition and regulation. Vascul Pharmacol 38: 323–337.
- 40. Cordenonsi M, D'Atri F, Hammar E, Parry DA, Kendrick-Jones J, et al. (1999) Cingulin contains globular and coiled-coil domains and interacts with ZO-1, ZO-2, ZO-3, and myosin. J Cell Biol 147: 1569–1582.
- 41. Persidsky Y, Ramirez SH, Haorah J, Kanmogne GD (2006) Blood-brain barrier: structural components and function under physiologic and pathologic conditions. J Neuroimmune Pharmacol 1: 223–236.
- 42. Sandoval KE, Witt KA (2008) Blood-brain barrier tight junction permeability and ischemic stroke. Neurobiol Dis 32: 200–219.
- 43. Vorbrodt AW, Dobrogowska DH (2003) Molecular anatomy of intercellular junctions in brain endothelial and epithelial barriers: electron microscopist's view. Brain Res Brain Res Rev 42: 221–242.
- 44. Vorbrodt AW, Li S, Brown WT, Ramakrishna N (2008) Increased expression of beta-catenin in brain microvessels of a segmentally trisomic (Ts65Dn) mouse model of Down syndrome. Brain Cell Biol 36: 203–211.
- 45. Choi YK, Kim KW (2008) Blood-neural barrier: its diversity and coordinated cell-to-cell communication. BMB Rep 41: 345–352.
- 46. Simionescu M, Gafencu A, Antohe F (2002) Transcytosis of plasma macromolecules in endothelial cells: A cell biological survey. Microscopy Research and Technique 57: 269–288.
- 47. Wolburg H, Noell S, Mack A, Wolburg-Buchholz K, Fallier-Becker P (2009) Brain endothelial cells and the glio-vascular complex. Cell Tissue Res 335: 75–96.
- 48. Vannucci SJ, Maher F, Simpson IA (1997) Glucose transporter proteins in brain: delivery of glucose to neurons and glia. Glia 21: 2–21.
- 49. Begley DJ (1996) The blood-brain barrier: principles for targeting peptides and drugs to the central nervous system. J Pharm Pharmacol 48: 136–146.
- 50. Hardebo JE, Owman C (1980) Barrier mechanisms for neurotransmitter monoamines and their precursors at the blood-brain interface. Ann Neurol 8: 1–31.
- 51. Tsukita S, Furuse M, Itoh M (2001) Multifunctional strands in tight junctions. Nat Rev Mol Cell Biol 2: 285–293.
- 52. Erickson KK, Sundstrom JM, Antonetti DA (2007) Vascular permeability in ocular disease and the role of tight junctions. Angiogenesis 10: 103–117.
- 53. Gould GW, Holman GD (1993) The glucose transporter family: structure, function and tissue-specific expression. Biochem J 295 (Pt 2): 329–341.
- 54. Shih JC, Chen K (2004) Regulation of MAO-A and MAO-B gene expression. Curr Med Chem 11: 1995–2005.
- 55. Wagner MC, Barylko B, Albanesi JP (1992) Tissue distribution and subcellular localization of mammalian myosin I. J Cell Biol 119: 163–170.
- 56. Navarro P, Ruco L, Dejana E (1998) Differential localization of VE- and N-cadherins in human endothelial cells: VE-cadherin competes with N-cadherin for junctional localization. J Cell Biol 140: 1475–1484.
- 57. Cordon-Cardo C, O'Brien JP, Boccia J, Casals D, Bertino JR, et al. (1990) Expression of the multidrug resistance gene product (P-glycoprotein) in human normal and tumor tissues. J Histochem Cytochem 38: 1277–1287.
- 58. Cha SH, Sekine T, Fukushima JI, Kanai Y, Kobayashi Y, et al. (2001) Identification and characterization of human organic anion transporter 3 expressing predominantly in the kidney. Mol Pharmacol 59: 1277–1286.
- 59. Albert VR, Allen JM, Joh TH (1987) A single gene codes for aromatic L-amino acid decarboxylase in both neuronal and non-neuronal tissues. J Biol Chem 262: 9404–9411.