<title>Introduction</title> <p>Despite the genome sequence of the human and several other mammalian species having been determined ##REF##11237011##[1]##–##REF##12466850##[2]##, it has become evident that the transcriptional output of the genome is much more complex than the predicted protein coding genes. Although a mammalian genome contains about 20,000 protein coding genes, there are an additional 23,000 transcriptional units (TUs) comprised of RNA without protein coding potential ##REF##16141072##[3]##, which originate from more than 181,000 different transcripts in the mouse. As detected by whole genome tiling arrays, the human genome also shows comparable complexity ##REF##15790807##[4]##; reviewed in ##REF##16859803##[5]##). EST analysis has shown that RNA expression is frequently restricted to specific tissues. For instance, more than 110,000 3′ end cDNA clusters, out of a total of 171,000 3′ end EST clusters prepared during the Mouse cDNA Encyclopedia Project, were identified from only one library ##REF##12819125##[6]##, largely from the brain ##REF##12819138##[7]##. These data suggest that to comprehensively identify the RNAs involved in specific biological processes, or at least the tissue-specific RNA/mRNA isoforms, it is essential to prepare novel EST collections, preferably using normalization and subtraction to identify rarely expressed mRNAs.</p> <p>Apart from the importance for expression profiling studies, full-length cDNAs are essential for the understanding of the structure of the protein encoded in specific tissues. Full-length cDNAs also constitute an invaluable resource for future functional studies such as ectopic expression with lentiviral vectors ##REF##17018643##[8]##. Additionally, they provide specific information for gene start and termination sites ##REF##16507138##[9]##. In particular, 5′ -ends help to identify tissue-specific promoter elements to connect them to TUs, and at the same time to novel datasets that allow for mRNA expression profiling based on measuring transcriptional activity at the start site (or mRNA cap). Accordingly, we have previously developed the cap-analysis gene expression (CAGE) method (##REF##16489339##[10]##–##REF##14663149##[11]##), which is based upon the production of short (20–21 nt) tags corresponding to the 5′ end of capped RNA transcripts. After high-throughput profiling, these are aligned onto the genome to identify their promoter elements, and the specific activity at a certain start site is measured as the frequency of CAGE tags. Full-length cDNA collections, and 5′ ESTs are very beneficial to map CAGE tags to TUs, which exhibit alternative promoter usage in at least 50% of the cases (##REF##16645617##[12]##).</p> <p>Here, we have developed a full-length cDNA and microarray resource for studies with particular emphasis on the mouse visual cortex. Due to tissue restriction of mRNA expression ##REF##12819125##[6]##, we tested the resource for differential gene expression across different cell types by comparing the transcriptome of parvalbumin (<italic>pvalb</italic>) cells versus the rest of cortex. The <italic>pvalb</italic>-cell network has recently been identified as an important trigger for critical period plasticity in the postnatal visual cortex ##REF##16261181##[13]##. Sensory deprivation typically between postnatal day P20–P40 in mice leads to a functional and structural rewiring of visual cortical circuits that underlies an enduring loss of visual acuity (amblyopia). Direct manipulation of postnatal GABA function can delay or accelerate plasticity onset ##REF##16261181##[13]##, carrying potentially broad implications for therapies in adulthood as well as a deeper mechanistic understanding of cognitive developmental disorders.</p>
<title>Materials and Methods</title> <title>Sequencing and clustering</title> <p>76,999 5′ ESTs from the K2, K3, K4 and K5 libraries were sequenced with RISA sequencers ##REF##11076861##[28]## and deposited in DDBJ under the accession numbers BY237236–BY298532 (all the range), CJ142844–CJ159349 (a subset of the range). 367 additional ESTs were used but not deposited in DDBJ because they were superseded by higher quality reads for the same clone. Their sequence is available upon request. The 5′ ESTs were clustered using stackPACK ##REF##10568754##[16]## with the following external programs: d2_cluster ##REF##10568753##[29]## (word_size = 6, similarity_cutoff = 0.98, minimum_sequence_size = 50, window_size = 200, reverse_comparison = 1), phrap (P. Green, unpublished) (old_ace = 1, vector_bound = 0, forcelevel = 5, trim_score = 20, penalty = −2, gap_init = −4, gap_ext = −3, ins_gap_ext = −3, del_gap_ext = −3, maxgap = 30, flags = -retain_duplicates -node_space 8) and craw ##REF##9521931##[30]## (sig = 0.5, window_size = 100, ignore_first = 50). In addition, 75,757 3′ ESTs were deposited in DDBJ under the accession numbers BY599749–BY660055 and CJ279459–CJ294908.</p> <title>Annotation</title> <p>The annotation of full-length cDNAs which were sequenced after completion of the FANTOM3 project was made in a pipeline containing five filters (##SUPPL##0##Figures S1## and ##SUPPL##1##S2##). First, the FANTOM3 representative transcript set ##REF##16141072##[3]## and refSeq ##REF##17130148##[31]## were searched for identical DNA sequences (&gt; = 98% similarity over &gt; = 100 bp) using the fasta34 program from the fasta3 package ##REF##9514730##[32]## with the parameters -Q -H -n -d0 -m9 -E100. Clones with no identical match were searched for open reading frames identical or similar to the entries of Swiss-Prot ##REF##12520024##[33]##, TrEMBL ##REF##12520024##[33]## or refSeq using the fasty34 program from the fasta3 package, with the parameters -Q -H -n -d0 -m9 -E100. cDNA clones with no match were scanned for potential open reading frames (ESTs were skipped) using the DECODER ##REF##11242592##[34]##, geneid 1.2 ##REF##1619647##[35]## and orfind (Tatusov, unpublished) programs. ESTs were compared to the UniGene database ##REF##12519941##[36]## using fasta34 with the same parameters. Clones which were not caught by any of these filters were mapped on the mouse genome (release UCSC mm5 ##REF##12466850##[2]##) using the sim4 ##REF##9750195##[37]## program.</p> <title>CAGE analysis</title> <p>Cage libraries were prepared as in ##REF##16489339##[10]##, using the visual cortex of mice at age P21, P26, P54 and P71. The libraries were sequenced and analyzed as in ##REF##16645617##[12]##, and were deposited in DDBJ under the accession numbers AAAAN0000001–AAAAN0074253, AAAAI0000001–AAAAI0019356, AAAAM0000001–AAAAM0054286 and AAAAK0000001–AAAAK0041805 respectively. The cerebellar CAGE library is available under the accession numbers AAAAA0000001–AAAA0240780.</p> <title>Rearray and clone amplification</title> <p>cDNA clones were rearrayed with a Q-bot instrument (Genetix) from 20% glycerol/LB ampicillin stocks that were kept at −80°C. A growth check was performed on all glycerol stocks. For the growth check, 384-well plates were translated into four 96-well plates per 384-well plate and one bacterial culture per glycerol stock was inoculated. All plates were analyzed after 18 h of growth at 37°C by visual inspection. We performed one Plasmid DNA Preparation per Glycerol Stock using a Montage Purification Kit (Millipore) and our robotic systems. One quality control was performed on three clones per plate, and the yield was 2.66 µg DNA±1.73. For the preparation of PCR templates, we diluted the plasmid DNA preparation by 1∶500. 5 µl of template were amplified in 50 µl reactions in 96-well plates with 1.25 U Expand High Fidelity Taq/Tgo polymerase blend (Roche), 15 pmol of forward primer (<named-content content-type="gene">
<sc>tgtaaaacgacggccagt</sc>
</named-content>), 15 pmol of reverse primer (<named-content content-type="gene">
<sc>agcggataacaatttcacacagga</sc>
</named-content>) and 200 µM of each dNTP (TaKaRa) in 1× Expand High Fidelity buffer (Roche). The program of the thermocycler was 2 min at 94°C followed by 30 cycles of 15 s at 94°C, 30 s at 60°C, and 2 min 30 s at 68°C, plus a final step of 7 min at 72°C. The amplification products were then purified with Multiscreen PCR plates (Millipore) following manufacturer instruction, and dried. With this protocol, we could achieve an average yield of &gt;8 µg for 92.7% of the wells.</p> <title>Animals and tissue preparation</title> <p>We used C57Bl/6 BAC transgenic mice expressing the enhanced green fluorescent protein (EGFP) under the control of the <italic>parvalbumin</italic> regulatory sequences ##REF##12177202##[21]##. Animals were maintained on a 12 h light/dark cycle with access to food and water <italic>ad libitum</italic>. All procedures involving animals and their care were carried out in accordance with the directives of RIKEN's Institutional Animal Care and Use Committee. All experimental groups were sacrificed at a similar time of day to avoid possible circadian effects. Primary visual cortex (area V1) from anesthetized mice was dissected and homogenized by sonication.</p> <title>Microarray analysis</title> <p>The PCR-amplified cDNAs were resuspended in 50% DMSO (Nacalai Tesque, Kyoto, Japan), and microarrays were printed on UltraGAPS slides (Corning) using an OmniGrid arrayer (Genomic Solutions) with Microarray Stealth Spotting pins (SMP3B, TeleChem). Detailed information about clone names, coordinates and accession numbers have been deposited in NCBIs Gene Expression Omnibus (GEO, <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/geo/">http://www.ncbi.nlm.nih.gov/geo/</ext-link>) and are accessible through GEO Platform accession number GPL5681.</p> <p>Total RNA was extracted from 15,000–20,000 FACS-sorted cells (Wagatsuma <italic>et al,</italic> in prep.) using the RNEasy Mini kit from Qiagen according to the manufacturer's instructions, and stored at −80°C until use.</p> <p>The following protocol is for one slide, and was scaled accordingly if more than one slide was to be hybridized. Total RNA samples were mixed with 2 µl of SPIKE Mix (GE Healthcare Lucidea universal scorecard), and probes were prepared following the method of ##REF##12089562##[22]##, except that 3 µg of cDNA were used per slide, and that the probes were dried and resuspended in 6 µl of water. The probes were then denatured at 95°C for 2 min and then mixed with 7.5 µl 4× Microarray hybridization solution version 2 (Amersham, Cat. No RPK0325), 15 µl formamide, and 145 µl of hybridization buffer (2.5 ng/µl oligo(dT) 12–18 primers (Invitrogen), 500 ng/µl mouse Cot-1 DNA (Invitrogen), 0.5× Microarray hybridization buffer version 2 (Amersham, Cat. No RPK0325, stock concentration: 4×), 25% formamide (Nacalai Tesque), 50% ULTRAhyb buffer (Ambion)), at 50°C.</p> <p>The mixture of Cy3- and Cy5-labeled probes was applied at 75°C onto arrays, and hybridized for 3 h at 65°C, 3 h at 55°C and 12 h at 50°C in the dark, using an automatic hybridizer, GeneMachines HybStation (Genomic Solutions). After hybridization, slides were washed once in 2× saline-sodium citrate (SSC) and 0.2% SDS at 53°C, in 1× SSC and 0.2% SDS at 53°C, in 2× SSC at 24°C, and in 0.2× SSC for 5 min at 24°C using one cycle of 10 s flow and 1 s hold, and 15 (first wash) or 13 (other washes) cycles of 1 s flow and 10 s hold, and then dried by spinning for 1 min.</p> <p>The scanning was done using an Axon GenPix 4000B scanner (Axon Instruments, Union City, CA) at 10 µm resolution. PMT voltage settings were varied to obtain maximum signal intensities with &lt;5% probe saturation. TIFF images were captured and analyzed with GenePix 6.0 (Axon Ins.) software and Acuity 4.0 (Axon Ins.). We used the Autoflag function to exclude spots for which the signal intensity was not higher than the local background. In addition, we also flagged the spots for which the scanning area was lower than 100 pixels. We call “standard” the experiments in which the probes prepared from the fluorescent cells were labelled with Cy5, and “swapped” the ones where the label was Cy3.</p> <title>Statistical analysis</title> <p>Lowess (locally weighted scatter plot for smoother) analysis was used for normalizing microarray data. Within-print-tip-group normalization was conducted by using the loess.smooth function in the statistical software package R ##UREF##0##[38]##. In short, this lowess normalization approach is to use the dependence of <italic>M = log<sub>2</sub> (Cy5/Cy3)</italic> on <italic>A = log<sub>2</sub> (Cy5×Cy3)</italic> for estimating its normalizing curve, which is a function of <italic>A</italic>
##UREF##1##[39]##–##REF##11842121##[40]##. Due to its estimation procedure, its normalization curve is considered as estimated mostly based on unexpressed clones (which should be the majority) in each small segment of <italic>A</italic>
##REF##11842121##[40]## (we used 50 segments). We tried to further ensure this point in our lowess normalization; for each segment of <italic>A</italic> we defined a window that included 40% of the points that are nearest. Inside this window, the standard deviation of <italic>M</italic> values was calculated and lowess normalization curve was estimated by using only data points whose <italic>M</italic> values fell into a range within three times the standard deviation in each <italic>A</italic> segment. Similar results could be obtained using smaller cutoff values, e.g. 1 s.d., suggesting that this three s.d. is stringent enough (data not shown). Once the normalization curve was estimated, data points out of the range were also normalized.</p> <p>To ensure data quality, we then excluded some clones from later statistical analysis by using the following two criteria, called global filtering and individual filtering. For global filtering, overall quality of each sample was examined using 3 repeated spottings (called replicates and denoted by rep1, rep2, rep3 below). For each standard or swap dataset (of each sample), the correlation of <italic>M</italic> values was calculated between replicates (e.g. the correlation between rep1 and rep2) and averaged over three cases (i.e. rep1 vs. rep2, rep2 vs. rep3, and rep3 vs. rep1). In all standard and swap datasets except the sample 2 standard dataset, the averaged correlation was larger than 0.899, whereas it was only 0.663 in the sample 2 standard dataset. We inspected the sample 2 standard dataset by creating the color image of <italic>M</italic> values in the format of the original microarray plate, and found blurred traces largely stretched over the plate. Given these, we decided to discard the sample 2 standard dataset from later statistical analysis.</p> <p>Next, for individual filtering, we examined the clone-wise consistency between standard and swap datasets (within each sample). We used the lowess-adjusted <italic>M</italic> values for this purpose and first obtained the mean of the adjusted <italic>M</italic> values over replicates for standard and swap datasets, denoted by  and , respectively, where <italic>i</italic> and <italic>j</italic> are the indices of clones and samples, respectively, <italic>i.e. i = 1,…4541</italic>, <italic>j = 1,…5</italic>. Second, for each sample (each fixed <italic>j</italic>), ideally we expect a consistency such that  is nearly equal to , where the multiplier −1 is necessary because of the swapped dye labeling. Hence, for each sample (each fixed <italic>j</italic>), we calculated the standard deviation of  over clones (<italic>i.e.</italic> over the index <italic>i</italic>). Clones were flagged, if they were outside the range of two times the standard deviation; these flagged clones were excluded from later statistical analysis.</p> <p>We used two statistics, <italic>p</italic>-value for identifying differentially expressed genes and <italic>q</italic>-value for dealing with issues of multiple hypotheses testing, <italic>i.e.</italic> controlling false discovery rate (see below). In calculating these, dye-averaged value was used for each clone, given by multiplying −1 to  as . Note that due to the filtering analysis described above, sample 2-standard dataset and also clones flagged by individual filtering were excluded from the calculation. In cases when the number of samples for each clone became less than 3 due to the above filtering or some other bad conditions, the clone was excluded from calculation to secure further robustness. <italic>p</italic>-value was obtained by calculating <italic>t</italic>-statistics  with <italic>n−1</italic> degrees of freedom, where <italic>n</italic> was the number of samples taken into account for calculating the mean of each clone. <italic>q</italic>-value can be used to assess false discovery rate (FDR; more precisely, pFDR, see ##UREF##2##[41]## for details). For this calculation, we used QVALUE software at <ext-link ext-link-type="uri" xlink:href="http://genomine.org/qvalue/">http://genomine.org/qvalue/</ext-link>. In short, for a specified threshold <italic>r</italic> (<italic>0&lt;r≤1</italic>), we can write , where <italic>m</italic> indicates the number of <italic>p<sub>i</sub></italic>, and  is estimated by  where <italic>λ</italic> is a tuning parameter to adjust the degree of conservative estimate of FDR and <italic>q</italic> value ##REF##12883005##[23]##. We used cubic spline at <italic>λ = 1</italic> and the <italic>q</italic> value is estimated by .</p>
<title>Results and Discussion</title> <p>We previously prepared four full-length cDNA libraries from whole extracts of wild-type mouse visual cortex at P18, P24, P28 and P55 using the CAP trapper method ##REF##8938445##[14]##. The libraries were subjected to subtraction and normalization ##REF##14513556##[15]##. A total of 76,999 5′ ends and 75,757 3′ ends were sequenced, and a cluster analysis showed that the resulting libraries had an average redundancy of 1.7 ##REF##14513556##[15]##. These libraries were used as a starting point for the preparation of a non-redundant clone collection representing the transcriptome of the visual cortex across its postnatal development.</p> <p>We pooled the 5′ ESTs from all four libraries and grouped them into 6,818 clusters and 9,391 singletons using the stackPACK software ##REF##10568754##[16]##. Thus, the redundancy increased to 4.75 after <italic>in silico</italic> pooling. We counted the number of 5′ ends per cluster, and identified the genes with the most imbalanced expression patterns, in order to detect potential genes of interest for the study of the maturation of visual cortex. We used the method developed by Stekel <italic>et al.</italic>
##REF##11116099##[17]## for comparing digital expression patterns, and calculated a <italic>R</italic> statistic for each cluster, using an in-house program available upon request. This method did not allow a quantitative prediction of false discovery rate. Hence, the upper interval of <italic>R v</italic>alues which cannot be fitted to an exponentially decreasing distribution were kept as potentially significant, as suggested in Stekel <italic>et al.</italic>
##REF##11116099##[17]##. We then arbitrarily selected 11 as a cut-off value (##FIG##0##Figure 1##), which highlighted fourteen protein-coding genes (##TAB##0##Table 1##).</p> <p>One representative clone was selected for each cluster. Together with the singletons, they constituted 16,209 different clones that were selected for rearray and PCR amplification in order to set up a microarray platform. 400 cDNAs were full-length sequenced and annotated as FANTOM 3 clones ##REF##16141072##[3]##. Of these, 179 have been annotated as non-coding by the FANTOM 3 consortium. This suggests that our whole collection is well enriched in non-coding cDNAs. To annotate the clones for which only partial sequence was available, we compared their 5′ ends to Transcription Units (TUs) of FANTOM3's Representative Transcript Set (RTS) that includes ESTs (<ext-link ext-link-type="uri" xlink:href="http://ftp://fantom.gsc.riken.jp/RTPS/fantom3_mouse/primary_est_rtps">ftp://fantom.gsc.riken.jp/RTPS/fantom3_mouse/primary_est_rtps</ext-link>). TUs are clusters of full-length transcripts that contain a common core of genetic information ##REF##12466851##[18]##. 13,189 clones whose 5′ EST was co-located to a TU inherited their annotation and identification number. The full-length cDNAs that define the TUs originate mostly from libraries made totally or partially from brain tissues (brain, whole head, whole embryo), from libraries related to the immune system (immune cells are present in most tissues and therefore contribute to their transcriptome), and from a testis library. The 3,020 clones whose 5′ end did not match a FANTOM3 TU are those whose sequence satisfied our quality criteria for clustering, but were discarded by the more stringent filters of the FANTOM3 RTS. There are 580 FANTOM 3 ESTs that are unique to our visual cortex libraries according to the FANTOM3 representative transcript set, of which 292 have been included in our spotted clone collection.</p> <p>To confirm the expression of our clones in the visual cortex, we prepared CAGE libraries ##REF##16489339##[10]##–##REF##14663149##[11]## from the visual cortex of mice at four developmental stages (P21, P26, P54 and P71) and compared the overlap of the cDNA and CAGE libraries for gene detection. We mapped the CAGE tags to the mouse genome in order to assign them to TUs ##REF##16645617##[12]## and pooled them in one virtual library of 180,984 tags. 16,194 TUs were identified by this technique as being expressed in the visual cortex. Of these, 6,874 had a counterpart in our clone collection. Interestingly, 7,005 TUs of our collection also had a counterpart in a CAGE library of similar sequencing depth (211,541 tags) made from cerebellar tissue ##REF##16645617##[12]##. This suggests that our collection can also be useful in the analysis of non-cortical brain tissues (##FIG##1##Figure 2##).</p> <p>Cerebral cortex is composed of many cellular and neuronal types. Our spotted clone collection represents their combined transcriptomes. We sought to investigate the complexity of a homogeneous subset of cortical cells. Parvalbumin (<italic>pvalb</italic>) is a calcium-binding protein which is used as a molecular signature of some GABAergic interneuron cell types ##REF##9498163##[19]##–##REF##10223310##[20]##. Using a BAC transgenic line expressing the green fluorescent protein (EGFP) under the control of the <italic>pvalb</italic> regulatory sequences ##REF##12177202##[21]##, we dissociated and sorted fluorescent and non-fluorescent cells from mouse cerebral cortex at the peak of the visual critical period (P28–29) and compared their transcriptional profiles on microarrays made with the first 4,512 clones of our collection, each spotted three times. Cy5- and Cy3-labeled cDNAs were prepared from the total RNA of 15,000–20,000 cells using the method of ##REF##12089562##[22]##, and compared five samples using a dye-swap strategy. The data discussed in this publication have been deposited in NCBIs Gene Expression Omnibus (GEO, <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/geo/">http://www.ncbi.nlm.nih.gov/geo/</ext-link>) and are accessible through GEO Series accession number GSE8968.</p> <p>The slides were scanned with the photomultiplier set to maximize signal intensity while avoiding excess saturation. The fluorescence intensity values were normalized using LOWESS normalization. To ensure data quality, we first investigated the overall quality of each slide and decided to discard the data from one slide. Second, we examined data consistency using dye-swap information within each sample and excluded clones that did not meet criterion from later statistical analysis, lowering the total number of included clones to 4,250. We ranked them by <italic>p</italic>-values based on a <italic>t</italic>-test. To address multiple comparisons, we chose <italic>q</italic>-values as a false discovery rate-based measure of significance ##REF##12883005##[23]##. By setting the significance criterion at <italic>p</italic>&lt;0.05, 1,400 of the 4,250 clones were differentially expressed (<italic>p</italic>&lt;0.05) (##TAB##1##Table 2##, ##FIG##2##Figure 3A##). The <italic>q</italic>-value corresponding to criterion was 0.066, so the number of false positives was estimated to be at most 92 (92 = <italic>q×n</italic>).</p> <p>As a control for the quality of cell sorting, we checked for genes known to be absent from <italic>pvalb</italic>-expressing neurons. Both GABA-synthetic enzymes, <italic>GAD65</italic> and <italic>GAD67</italic>, were included as control genes, however, they are not ideal, as <italic>pvalb</italic>-cells constitute only a subset of GABAergic neuron ##REF##16261181##[13]##. Thus, the non-<italic>pvalb</italic> fraction of our homogenates would include a substantial amount of <italic>GAD</italic> gene expression. Instead, we used <italic>pvalb</italic> itself as a control. Moreover, the glutamate receptor <italic>Gria2</italic> (<italic>GluR2</italic>) was significantly underrepresented in the sample prepared from fluorescent cells (<italic>p</italic> = 0.00048, <italic>M</italic> = −0.77), as expected for <italic>pvalb</italic>-cells that express little or no GluR2 subunits ##REF##8011338##[24]##. Instead, <italic>Gria1</italic>, another receptor from the same family that is expressed in <italic>pvalb</italic>-expressing cells, was over-represented in the fluorescent samples (<italic>p</italic> = 0.0042, <italic>M</italic> = 1.2).</p> <p>We inspected the significantly over- or under-represented genes exhibiting the strongest expression differences between EGFP-positive and negative samples. ##TAB##2##Table 3## lists the 24 strongest variations in both directions. Prompted by the co-occurrence in ##TAB##2##Table 3## of genes implicated in ATP metabolism, such as subunits of the mitochondrial ATPase or components of the oxidative chain, we examined all other genes from these families, and found that they are particularly abundant among over-represented genes (##FIG##2##Figure 3B##, ##TAB##3##Table 4##). This may reflect the intense metabolic activity of the fast-spiking <italic>pvalb</italic>-expressing cells ##REF##18054855##[25]##.</p> <p>To confirm that the fluorescent cells consisted only of neurons, we checked for the expression of <italic>Myelin Basic Protein</italic> (<italic>Mbp</italic>). Our slides contained two different clones coding for the <italic>Mbp</italic> gene, both showing a strong underrepresentation (<italic>M</italic>≈−2.1 each). Despite this strong value, it was not prominently ranked since the <italic>p-</italic>value was slightly lower than 5%. These two <italic>Mbp</italic> clones exemplify a class of genes which display a strong variability in their differential expression across the analyzed samples (##FIG##2##Figure 3C##). Among these, we could also find another myelin-related gene, <italic>plp1</italic>, and a component of the MAP kinase pathway, <italic>Map2k6</italic>. In the case of the <italic>Mbp</italic> gene, we hypothesize that the source of the noise affecting the <italic>p-</italic>value was the quantity of oligodendrocytes included amongst the non-fluorescent cells, which can vary significantly with each dissection.</p> <p>We then sought to investigate the differential expression of non-coding transcripts. Since we do not have the full-length sequence of most of these clones, we relied on their annotation to identify those that are potentially non-coding. We reasoned that most protein-coding genes are given an original gene symbol, whereas the clones whose function was unknown keep their transient MGI (Mouse Genome Informatics) name for a long time. In our collection, these clones typically have a symbol constructed by adding “Rik” to the clone identifier. Using the FANTOM 3 genome browser, we inspected the mapping of the the top 100 clones (ranked by their <italic>p</italic>-value) which inherited a gene symbol including “Rik”. Only seven of these were likely to be non-coding (##FIG##3##Figure 4##).</p> <p>This low occurrence of non-coding clones in the differentially regulated set can be explained as follows. First, non-coding RNAs are often weakly expressed, which makes the detection of variation in expression harder. Second, it has been reported that more than one third of all transcribed sequences exist in polyadenylated and non-poly-adenylated forms ##REF##15790807##[4]##. Since we used oligo dT primers during the reverse-transcription, this aggravates low expression levels. Importantly, because the function of the non-coding transcriptome is still poorly characterized, it remains unknown whether technical reasons alone or biological considerations explain the low occurrence of differentially expressed non-coding transcripts. From this perspective, our results suggest that the non-coding transcriptome is homogeneous at the level of different mature cortical populations.</p> <p>We created a collection of 16,209 RIKEN clones by subtracting, normalizing and clustering cDNA libraries from the mouse visual cortex. This clone collection is representative of the brain transcriptome as confirmed by an independent transcriptional profiling using the CAGE technology. Comparison of a homogeneous interneuron population to a whole-tissue control with our cDNA microarray platform identified more than 30% of differentially expressed RNAs even taking a conservative approach. Besides full-length cDNA arrays, and the investment in equipment and resources that these technologies require, other alternatives strategies can be conceived. Possible choices are the oligonucleotide arrays. However, due to prevalence of alternative promoters, exons and polyadenylation site, one would ideally have to prepare multiple oligonucleotides for each coding and non-coding gene. Our clone collection aims at maximizing the possibility of implicating novel non-coding cDNAs in a neural function, and its advantage is that when a candidate is found, a full-length cDNA is available. For analysis focused on transcriptional control, we would rather recommend CAGE when the technology to make CAGE libraries from limited amount of tissue will be available.</p> <p>The large heterogeneity of gene expression between cells strongly suggests that future studies will require the separation of single neuron populations. This is particularly true for cells that are rare in a tissue, such as the <italic>pvalb</italic>-expressing neurons, whose specific transcriptome variations during development and brain plasticity would be masked by the large majority of unrelated cells. In general, these data emphasize once more the importance of dissecting neuronal sub-types when expression profiling to study specific brain regions. In summary, these findings show that our clone collection, assembled with an unbiased and hypothesis-free strategy, is suitable for probing the transcriptome of purified cell populations and is broad enough to provide control genes. Our publicly available cDNAs are thus a proven tool for the profiling and comparison of homogeneous cell types within the brain. This will help progress toward better defining these cell types at a transcriptomic level ##REF##16301166##[26]##–##REF##16962313##[27]##.</p>
<title>Results and Discussion</title> <p>We previously prepared four full-length cDNA libraries from whole extracts of wild-type mouse visual cortex at P18, P24, P28 and P55 using the CAP trapper method ##REF##8938445##[14]##. The libraries were subjected to subtraction and normalization ##REF##14513556##[15]##. A total of 76,999 5′ ends and 75,757 3′ ends were sequenced, and a cluster analysis showed that the resulting libraries had an average redundancy of 1.7 ##REF##14513556##[15]##. These libraries were used as a starting point for the preparation of a non-redundant clone collection representing the transcriptome of the visual cortex across its postnatal development.</p> <p>We pooled the 5′ ESTs from all four libraries and grouped them into 6,818 clusters and 9,391 singletons using the stackPACK software ##REF##10568754##[16]##. Thus, the redundancy increased to 4.75 after <italic>in silico</italic> pooling. We counted the number of 5′ ends per cluster, and identified the genes with the most imbalanced expression patterns, in order to detect potential genes of interest for the study of the maturation of visual cortex. We used the method developed by Stekel <italic>et al.</italic>
##REF##11116099##[17]## for comparing digital expression patterns, and calculated a <italic>R</italic> statistic for each cluster, using an in-house program available upon request. This method did not allow a quantitative prediction of false discovery rate. Hence, the upper interval of <italic>R v</italic>alues which cannot be fitted to an exponentially decreasing distribution were kept as potentially significant, as suggested in Stekel <italic>et al.</italic>
##REF##11116099##[17]##. We then arbitrarily selected 11 as a cut-off value (##FIG##0##Figure 1##), which highlighted fourteen protein-coding genes (##TAB##0##Table 1##).</p> <p>One representative clone was selected for each cluster. Together with the singletons, they constituted 16,209 different clones that were selected for rearray and PCR amplification in order to set up a microarray platform. 400 cDNAs were full-length sequenced and annotated as FANTOM 3 clones ##REF##16141072##[3]##. Of these, 179 have been annotated as non-coding by the FANTOM 3 consortium. This suggests that our whole collection is well enriched in non-coding cDNAs. To annotate the clones for which only partial sequence was available, we compared their 5′ ends to Transcription Units (TUs) of FANTOM3's Representative Transcript Set (RTS) that includes ESTs (<ext-link ext-link-type="uri" xlink:href="http://ftp://fantom.gsc.riken.jp/RTPS/fantom3_mouse/primary_est_rtps">ftp://fantom.gsc.riken.jp/RTPS/fantom3_mouse/primary_est_rtps</ext-link>). TUs are clusters of full-length transcripts that contain a common core of genetic information ##REF##12466851##[18]##. 13,189 clones whose 5′ EST was co-located to a TU inherited their annotation and identification number. The full-length cDNAs that define the TUs originate mostly from libraries made totally or partially from brain tissues (brain, whole head, whole embryo), from libraries related to the immune system (immune cells are present in most tissues and therefore contribute to their transcriptome), and from a testis library. The 3,020 clones whose 5′ end did not match a FANTOM3 TU are those whose sequence satisfied our quality criteria for clustering, but were discarded by the more stringent filters of the FANTOM3 RTS. There are 580 FANTOM 3 ESTs that are unique to our visual cortex libraries according to the FANTOM3 representative transcript set, of which 292 have been included in our spotted clone collection.</p> <p>To confirm the expression of our clones in the visual cortex, we prepared CAGE libraries ##REF##16489339##[10]##–##REF##14663149##[11]## from the visual cortex of mice at four developmental stages (P21, P26, P54 and P71) and compared the overlap of the cDNA and CAGE libraries for gene detection. We mapped the CAGE tags to the mouse genome in order to assign them to TUs ##REF##16645617##[12]## and pooled them in one virtual library of 180,984 tags. 16,194 TUs were identified by this technique as being expressed in the visual cortex. Of these, 6,874 had a counterpart in our clone collection. Interestingly, 7,005 TUs of our collection also had a counterpart in a CAGE library of similar sequencing depth (211,541 tags) made from cerebellar tissue ##REF##16645617##[12]##. This suggests that our collection can also be useful in the analysis of non-cortical brain tissues (##FIG##1##Figure 2##).</p> <p>Cerebral cortex is composed of many cellular and neuronal types. Our spotted clone collection represents their combined transcriptomes. We sought to investigate the complexity of a homogeneous subset of cortical cells. Parvalbumin (<italic>pvalb</italic>) is a calcium-binding protein which is used as a molecular signature of some GABAergic interneuron cell types ##REF##9498163##[19]##–##REF##10223310##[20]##. Using a BAC transgenic line expressing the green fluorescent protein (EGFP) under the control of the <italic>pvalb</italic> regulatory sequences ##REF##12177202##[21]##, we dissociated and sorted fluorescent and non-fluorescent cells from mouse cerebral cortex at the peak of the visual critical period (P28–29) and compared their transcriptional profiles on microarrays made with the first 4,512 clones of our collection, each spotted three times. Cy5- and Cy3-labeled cDNAs were prepared from the total RNA of 15,000–20,000 cells using the method of ##REF##12089562##[22]##, and compared five samples using a dye-swap strategy. The data discussed in this publication have been deposited in NCBIs Gene Expression Omnibus (GEO, <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/geo/">http://www.ncbi.nlm.nih.gov/geo/</ext-link>) and are accessible through GEO Series accession number GSE8968.</p> <p>The slides were scanned with the photomultiplier set to maximize signal intensity while avoiding excess saturation. The fluorescence intensity values were normalized using LOWESS normalization. To ensure data quality, we first investigated the overall quality of each slide and decided to discard the data from one slide. Second, we examined data consistency using dye-swap information within each sample and excluded clones that did not meet criterion from later statistical analysis, lowering the total number of included clones to 4,250. We ranked them by <italic>p</italic>-values based on a <italic>t</italic>-test. To address multiple comparisons, we chose <italic>q</italic>-values as a false discovery rate-based measure of significance ##REF##12883005##[23]##. By setting the significance criterion at <italic>p</italic>&lt;0.05, 1,400 of the 4,250 clones were differentially expressed (<italic>p</italic>&lt;0.05) (##TAB##1##Table 2##, ##FIG##2##Figure 3A##). The <italic>q</italic>-value corresponding to criterion was 0.066, so the number of false positives was estimated to be at most 92 (92 = <italic>q×n</italic>).</p> <p>As a control for the quality of cell sorting, we checked for genes known to be absent from <italic>pvalb</italic>-expressing neurons. Both GABA-synthetic enzymes, <italic>GAD65</italic> and <italic>GAD67</italic>, were included as control genes, however, they are not ideal, as <italic>pvalb</italic>-cells constitute only a subset of GABAergic neuron ##REF##16261181##[13]##. Thus, the non-<italic>pvalb</italic> fraction of our homogenates would include a substantial amount of <italic>GAD</italic> gene expression. Instead, we used <italic>pvalb</italic> itself as a control. Moreover, the glutamate receptor <italic>Gria2</italic> (<italic>GluR2</italic>) was significantly underrepresented in the sample prepared from fluorescent cells (<italic>p</italic> = 0.00048, <italic>M</italic> = −0.77), as expected for <italic>pvalb</italic>-cells that express little or no GluR2 subunits ##REF##8011338##[24]##. Instead, <italic>Gria1</italic>, another receptor from the same family that is expressed in <italic>pvalb</italic>-expressing cells, was over-represented in the fluorescent samples (<italic>p</italic> = 0.0042, <italic>M</italic> = 1.2).</p> <p>We inspected the significantly over- or under-represented genes exhibiting the strongest expression differences between EGFP-positive and negative samples. ##TAB##2##Table 3## lists the 24 strongest variations in both directions. Prompted by the co-occurrence in ##TAB##2##Table 3## of genes implicated in ATP metabolism, such as subunits of the mitochondrial ATPase or components of the oxidative chain, we examined all other genes from these families, and found that they are particularly abundant among over-represented genes (##FIG##2##Figure 3B##, ##TAB##3##Table 4##). This may reflect the intense metabolic activity of the fast-spiking <italic>pvalb</italic>-expressing cells ##REF##18054855##[25]##.</p> <p>To confirm that the fluorescent cells consisted only of neurons, we checked for the expression of <italic>Myelin Basic Protein</italic> (<italic>Mbp</italic>). Our slides contained two different clones coding for the <italic>Mbp</italic> gene, both showing a strong underrepresentation (<italic>M</italic>≈−2.1 each). Despite this strong value, it was not prominently ranked since the <italic>p-</italic>value was slightly lower than 5%. These two <italic>Mbp</italic> clones exemplify a class of genes which display a strong variability in their differential expression across the analyzed samples (##FIG##2##Figure 3C##). Among these, we could also find another myelin-related gene, <italic>plp1</italic>, and a component of the MAP kinase pathway, <italic>Map2k6</italic>. In the case of the <italic>Mbp</italic> gene, we hypothesize that the source of the noise affecting the <italic>p-</italic>value was the quantity of oligodendrocytes included amongst the non-fluorescent cells, which can vary significantly with each dissection.</p> <p>We then sought to investigate the differential expression of non-coding transcripts. Since we do not have the full-length sequence of most of these clones, we relied on their annotation to identify those that are potentially non-coding. We reasoned that most protein-coding genes are given an original gene symbol, whereas the clones whose function was unknown keep their transient MGI (Mouse Genome Informatics) name for a long time. In our collection, these clones typically have a symbol constructed by adding “Rik” to the clone identifier. Using the FANTOM 3 genome browser, we inspected the mapping of the the top 100 clones (ranked by their <italic>p</italic>-value) which inherited a gene symbol including “Rik”. Only seven of these were likely to be non-coding (##FIG##3##Figure 4##).</p> <p>This low occurrence of non-coding clones in the differentially regulated set can be explained as follows. First, non-coding RNAs are often weakly expressed, which makes the detection of variation in expression harder. Second, it has been reported that more than one third of all transcribed sequences exist in polyadenylated and non-poly-adenylated forms ##REF##15790807##[4]##. Since we used oligo dT primers during the reverse-transcription, this aggravates low expression levels. Importantly, because the function of the non-coding transcriptome is still poorly characterized, it remains unknown whether technical reasons alone or biological considerations explain the low occurrence of differentially expressed non-coding transcripts. From this perspective, our results suggest that the non-coding transcriptome is homogeneous at the level of different mature cortical populations.</p> <p>We created a collection of 16,209 RIKEN clones by subtracting, normalizing and clustering cDNA libraries from the mouse visual cortex. This clone collection is representative of the brain transcriptome as confirmed by an independent transcriptional profiling using the CAGE technology. Comparison of a homogeneous interneuron population to a whole-tissue control with our cDNA microarray platform identified more than 30% of differentially expressed RNAs even taking a conservative approach. Besides full-length cDNA arrays, and the investment in equipment and resources that these technologies require, other alternatives strategies can be conceived. Possible choices are the oligonucleotide arrays. However, due to prevalence of alternative promoters, exons and polyadenylation site, one would ideally have to prepare multiple oligonucleotides for each coding and non-coding gene. Our clone collection aims at maximizing the possibility of implicating novel non-coding cDNAs in a neural function, and its advantage is that when a candidate is found, a full-length cDNA is available. For analysis focused on transcriptional control, we would rather recommend CAGE when the technology to make CAGE libraries from limited amount of tissue will be available.</p> <p>The large heterogeneity of gene expression between cells strongly suggests that future studies will require the separation of single neuron populations. This is particularly true for cells that are rare in a tissue, such as the <italic>pvalb</italic>-expressing neurons, whose specific transcriptome variations during development and brain plasticity would be masked by the large majority of unrelated cells. In general, these data emphasize once more the importance of dissecting neuronal sub-types when expression profiling to study specific brain regions. In summary, these findings show that our clone collection, assembled with an unbiased and hypothesis-free strategy, is suitable for probing the transcriptome of purified cell populations and is broad enough to provide control genes. Our publicly available cDNAs are thus a proven tool for the profiling and comparison of homogeneous cell types within the brain. This will help progress toward better defining these cell types at a transcriptomic level ##REF##16301166##[26]##–##REF##16962313##[27]##.</p>
<p>Conceived and designed the experiments: Mf TH PC. Performed the experiments: Mf AW TK AAO YK FS. Analyzed the data: CP Mf NH TK HN. Contributed reagents/materials/analysis tools: KW HN. Wrote the paper: CP Mf AW NH TK AAO YK KW HN TH PC.</p> <title>Background</title> <p>The transcriptome of the cerebral cortex is remarkably homogeneous, with variations being stronger between individuals than between areas. It is thought that due to the presence of many distinct cell types, differences within one cell population will be averaged with the noise from others. Studies of sorted cells expressing the same transgene have shown that cell populations can be distinguished according to their transcriptional profile.</p> <title>Methodology</title> <p>We have prepared a low-redundancy set of 16,209 full-length cDNA clones which represents the transcriptome of the mouse visual cortex in its coding and non-coding aspects. Using an independent tag-based approach, CAGE, we confirmed the cortical expression of 72% of the clones. Clones were amplified by PCR and spotted on glass slides, and we interrogated the microarrays with RNA from flow-sorted fluorescent cells from the cerebral cortex of parvalbumin-egfp transgenic mice.</p> <title>Conclusions</title> <p>We provide an annotated cDNA clone collection which is particularly suitable for transcriptomic analysis in the mouse brain. Spotting it on microarrays, we compared the transcriptome of EGFP positive and negative cells in a parvalbumin-egfp transgenic background and showed that more than 30% of clones are differentially expressed. Our clone collection will be a useful resource for the study of the transcriptome of single cell types in the cerebral cortex.</p>
<title>Supporting Information</title>
<p>We would like to thank K. Nakano and J. Horiuchi at NTT Software for supporting database construction, M. Kato and A. Itoh for assistance in programming, F. Hori for managing EST sequencing data, H. Monyer (Heidelberg Univ) for kindly providing <italic>pvalb-EGFP</italic> mouse breeding pairs, H. Uruma for secretarial assistance, and My Persson for English corrections. This paper is dedicated to the memory of Sayaka Fujishima (b.1975-d.2006).</p>
<fig id="pone-0003012-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003012.g001</object-id><label>Figure 1</label><caption><title>Counts in discrete intervals of R-values (natural logarithmic scale).</title><p>The line was fit to the part of the curve which can be approximated by a logarithmic function. Its intersection with the x axis is suggested as a rough threshold estimate above which R-values are considered significant.</p></caption></fig> <fig id="pone-0003012-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003012.g002</object-id><label>Figure 2</label><caption><title>Venn diagram representing the overlap between transcription units (TUs) of our clone collection and two sets of CAGE tags originating from the cerebellum or visual cortex.</title><p>Only 15% of TUs do not have a counterpart in the CAGE libraries, suggesting that our clone collection is well suited to interrogate brain samples.</p></caption></fig> <fig id="pone-0003012-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003012.g003</object-id><label>Figure 3</label><caption><title>Statistical analysis of microarray results.</title><p>A: q and p values ranked by p-values. B: Vulcano plot showing the p-values as a function of expression changes. Blue diamonds: genes from the Nduf family. Red triangles: genes from the Atp5 family. Green dots: genes from the Cox6 family. C: Absolute value of M ranked by increasing values of p. The red bars represent the standard deviation of M. Some genes have a low rank despite a strong M value, often because of their higher variability, as indicated by the length of their vertical bars.</p></caption></fig> <fig id="pone-0003012-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003012.g004</object-id><label>Figure 4</label><caption><title>Non-coding clones which show over-representation (M&gt;0) or under-representation (M&lt;0) in the parvalbumin-egfp cells.</title><p>Pictures extracted from the FANTOM 3 genome browser. The RIKEN clone ID is indicated, followed by the p and M values. Black rulers: chromosomal coordinates. Red arrows: known promoters. Stronger promoters are represented by bigger arrows with a stronger color. Plain arrows: known full-length cDNAs (brown boxing: plus strand; magenta boxing: reverse strand). Pink areas in the cDNAs: predicted coding sequences. The 5′ ends of the selected clones (and the 3′ end, in the case of K23022K18) are represented by plain arrows in a yellow background. K330007E03: clone overlapping the promoter of a coding gene. K230014I16: variant of a spliced non-coding gene. K230322K18: intronic clone overlapping the last exon of a coding gene. K230002G07: unspliced clone sharing the promoter of a coding gene in a head-to-head arrangement. K230002P11: intergenic clone. The coding sequence detected in one of the cDNAs is unlikely to be productive as the ATG is starting the transcript. K230013M23: spliced clone matching cDNAs with no known coding sequence. K230008M20: clone forming a sense-antisense chain with two genes oriented head-to-head.</p></caption></fig>
<table-wrap id="pone-0003012-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003012.t001</object-id><label>Table 1</label><caption><title>EST clusters with a strongly imbalanced developmental expression pattern.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Clone ID</td><td align="left" rowspan="1" colspan="1">
<italic>R</italic>-value</td><td align="left" rowspan="1" colspan="1">EST count</td><td align="left" rowspan="1" colspan="1">Gene symbol</td><td align="left" rowspan="1" colspan="1">Full official name</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">K230007N19</td><td align="left" rowspan="1" colspan="1">15.5</td><td align="left" rowspan="1" colspan="1">43,9,23,8</td><td align="left" rowspan="1" colspan="1">Exdl2</td><td align="left" rowspan="1" colspan="1">exonuclease 3′-5′ domain-like 2</td></tr><tr><td align="left" rowspan="1" colspan="1">K230008D05</td><td align="left" rowspan="1" colspan="1">14.3</td><td align="left" rowspan="1" colspan="1">4,0,17,0</td><td align="left" rowspan="1" colspan="1">Kctd10</td><td align="left" rowspan="1" colspan="1">potassium channel tetramerisation domain containing 10</td></tr><tr><td align="left" rowspan="1" colspan="1">K330007G06</td><td align="left" rowspan="1" colspan="1">14.3</td><td align="left" rowspan="1" colspan="1">0,11,0,2</td><td align="left" rowspan="1" colspan="1">Aldh1b1</td><td align="left" rowspan="1" colspan="1">aldehyde dehydrogenase 1 family, member B1</td></tr><tr><td align="left" rowspan="1" colspan="1">K230004E20</td><td align="left" rowspan="1" colspan="1">13.8</td><td align="left" rowspan="1" colspan="1">26,5,16,2</td><td align="left" rowspan="1" colspan="1">Fbxw7</td><td align="left" rowspan="1" colspan="1">F-box and WD-40 domain protein 7, archipelago homolog Drosophila)</td></tr><tr><td align="left" rowspan="1" colspan="1">K430007K11</td><td align="left" rowspan="1" colspan="1">13.5</td><td align="left" rowspan="1" colspan="1">0,0,13,2</td><td align="left" rowspan="1" colspan="1">Kcnv1</td><td align="left" rowspan="1" colspan="1">potassium channel, subfamily V, member 1</td></tr><tr><td align="left" rowspan="1" colspan="1">K230005P19</td><td align="left" rowspan="1" colspan="1">13.5</td><td align="left" rowspan="1" colspan="1">22,12,1,20</td><td align="left" rowspan="1" colspan="1">Decr2</td><td align="left" rowspan="1" colspan="1">2-4-dienoyl-Coenzyme A reductase 2, peroxisomal</td></tr><tr><td align="left" rowspan="1" colspan="1">K330313B11</td><td align="left" rowspan="1" colspan="1">13.2</td><td align="left" rowspan="1" colspan="1">0,2,1,14</td><td align="left" rowspan="1" colspan="1">1500002O20Rik</td><td align="left" rowspan="1" colspan="1">RIKEN cDNA 1500002O20 gene</td></tr><tr><td align="left" rowspan="1" colspan="1">K230309E11</td><td align="left" rowspan="1" colspan="1">12.9</td><td align="left" rowspan="1" colspan="1">1,11,1,0</td><td align="left" rowspan="1" colspan="1">Ev1</td><td align="left" rowspan="1" colspan="1">Ena-vasodilator stimulated phosphoprotein</td></tr><tr><td align="left" rowspan="1" colspan="1">K230034K22</td><td align="left" rowspan="1" colspan="1">12.2</td><td align="left" rowspan="1" colspan="1">2,13,1,1</td><td align="left" rowspan="1" colspan="1">4632419K20Rik</td><td align="left" rowspan="1" colspan="1">RIKEN cDNA 4632419K20 gene</td></tr><tr><td align="left" rowspan="1" colspan="1">K230035O10</td><td align="left" rowspan="1" colspan="1">12.2</td><td align="left" rowspan="1" colspan="1">2,0,15,2</td><td align="left" rowspan="1" colspan="1">Srgap2</td><td align="left" rowspan="1" colspan="1">SLIT-ROBO Rho GTPase activating protein 3</td></tr><tr><td align="left" rowspan="1" colspan="1">K230025H24</td><td align="left" rowspan="1" colspan="1">12.1</td><td align="left" rowspan="1" colspan="1">2,0,16,10</td><td align="left" rowspan="1" colspan="1">Mfn2</td><td align="left" rowspan="1" colspan="1">mitofusin 2</td></tr><tr><td align="left" rowspan="1" colspan="1">K530003G22</td><td align="left" rowspan="1" colspan="1">12.0</td><td align="left" rowspan="1" colspan="1">0,0,0,9</td><td align="left" rowspan="1" colspan="1">Arid1b</td><td align="left" rowspan="1" colspan="1">AT rich interactive domain 1B (Swi1 like)</td></tr><tr><td align="left" rowspan="1" colspan="1">K230006M16</td><td align="left" rowspan="1" colspan="1">11.7</td><td align="left" rowspan="1" colspan="1">3,1,20,5</td><td align="left" rowspan="1" colspan="1">Rbm39</td><td align="left" rowspan="1" colspan="1">RNA binding motif 39</td></tr><tr><td align="left" rowspan="1" colspan="1">K430007G04</td><td align="left" rowspan="1" colspan="1">11.6</td><td align="left" rowspan="1" colspan="1">0,0,9,0</td><td align="left" rowspan="1" colspan="1">6330500D04Rik</td><td align="left" rowspan="1" colspan="1">RIKEN cDNA 6330500D04 gene</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003012-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003012.t002</object-id><label>Table 2</label><caption><title>Lowest <italic>p</italic>- and <italic>q</italic>-values for given number of ranked clones, and number of clones reaching given <italic>p</italic>- or <italic>q</italic>-values.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">top <italic>n</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>p</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>q</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>n×q</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>n</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>p</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>q</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>n×q</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">2000</td><td align="left" rowspan="1" colspan="1">0.12</td><td align="left" rowspan="1" colspan="1">0.11</td><td align="left" rowspan="1" colspan="1">220</td><td align="left" rowspan="1" colspan="1">1400</td><td align="left" rowspan="1" colspan="1">&lt;0.05</td><td align="left" rowspan="1" colspan="1">0.066</td><td align="left" rowspan="1" colspan="1">92</td></tr><tr><td align="left" rowspan="1" colspan="1">1000</td><td align="left" rowspan="1" colspan="1">0.025</td><td align="left" rowspan="1" colspan="1">0.046</td><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">1062</td><td align="left" rowspan="1" colspan="1">0.029</td><td align="left" rowspan="1" colspan="1">&lt;0.05</td><td align="left" rowspan="1" colspan="1">53</td></tr><tr><td align="left" rowspan="1" colspan="1">500</td><td align="left" rowspan="1" colspan="1">0.0073</td><td align="left" rowspan="1" colspan="1">0.027</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">617</td><td align="left" rowspan="1" colspan="1">&lt;0.01</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">18</td></tr><tr><td align="left" rowspan="1" colspan="1">200</td><td align="left" rowspan="1" colspan="1">0.0019</td><td align="left" rowspan="1" colspan="1">0.017</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">0.00025</td><td align="left" rowspan="1" colspan="1">&lt;0.01</td><td align="left" rowspan="1" colspan="1">0</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003012-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003012.t003</object-id><label>Table 3</label><caption><title>Clones showing the strongest statistically significant differential expression.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td colspan="4" align="left" rowspan="1">
<italic>under-represented</italic>
</td><td colspan="4" align="left" rowspan="1">
<italic>over-represented</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">SYMBOL</td><td align="left" rowspan="1" colspan="1">CLONE_ID</td><td align="left" rowspan="1" colspan="1">p</td><td align="left" rowspan="1" colspan="1">M</td><td align="left" rowspan="1" colspan="1">SYMBOL</td><td align="left" rowspan="1" colspan="1">CLONE_ID</td><td align="left" rowspan="1" colspan="1">p</td><td align="left" rowspan="1" colspan="1">M</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Crip2</td><td align="left" rowspan="1" colspan="1">K230313M18</td><td align="left" rowspan="1" colspan="1">7.58E-3</td><td align="left" rowspan="1" colspan="1">−2.69</td><td align="left" rowspan="1" colspan="1">Cplx1</td><td align="left" rowspan="1" colspan="1">K230003P17</td><td align="left" rowspan="1" colspan="1">4.44E-3</td><td align="left" rowspan="1" colspan="1">2.58</td></tr><tr><td align="left" rowspan="1" colspan="1">C1qb</td><td align="left" rowspan="1" colspan="1">K230028J21</td><td align="left" rowspan="1" colspan="1">1.00E-2</td><td align="left" rowspan="1" colspan="1">−2.68</td><td align="left" rowspan="1" colspan="1">Fgf9</td><td align="left" rowspan="1" colspan="1">K230018I07</td><td align="left" rowspan="1" colspan="1">9.07E-5</td><td align="left" rowspan="1" colspan="1">2.32</td></tr><tr><td align="left" rowspan="1" colspan="1">Syt5</td><td align="left" rowspan="1" colspan="1">K230336L04</td><td align="left" rowspan="1" colspan="1">2.31E-3</td><td align="left" rowspan="1" colspan="1">−2.40</td><td align="left" rowspan="1" colspan="1">Gtl2</td><td align="left" rowspan="1" colspan="1">K230042M22</td><td align="left" rowspan="1" colspan="1">3.93E-3</td><td align="left" rowspan="1" colspan="1">1.88</td></tr><tr><td align="left" rowspan="1" colspan="1">Ctxn</td><td align="left" rowspan="1" colspan="1">K230013F20</td><td align="left" rowspan="1" colspan="1">1.93E-3</td><td align="left" rowspan="1" colspan="1">−2.33</td><td align="left" rowspan="1" colspan="1">9130213B05Rik</td><td align="left" rowspan="1" colspan="1">K230034N12</td><td align="left" rowspan="1" colspan="1">9.41E-4</td><td align="left" rowspan="1" colspan="1">1.69</td></tr><tr><td align="left" rowspan="1" colspan="1">Sepp1</td><td align="left" rowspan="1" colspan="1">K230009C19</td><td align="left" rowspan="1" colspan="1">4.85E-4</td><td align="left" rowspan="1" colspan="1">−2.23</td><td align="left" rowspan="1" colspan="1">6330512M04Rik</td><td align="left" rowspan="1" colspan="1">K230318F21</td><td align="left" rowspan="1" colspan="1">5.29E-3</td><td align="left" rowspan="1" colspan="1">1.62</td></tr><tr><td align="left" rowspan="1" colspan="1">Apoe</td><td align="left" rowspan="1" colspan="1">K230004K09</td><td align="left" rowspan="1" colspan="1">1.66E-2</td><td align="left" rowspan="1" colspan="1">−2.23</td><td align="left" rowspan="1" colspan="1">Ndufa11</td><td align="left" rowspan="1" colspan="1">K230031N16</td><td align="left" rowspan="1" colspan="1">7.69E-3</td><td align="left" rowspan="1" colspan="1">1.61</td></tr><tr><td align="left" rowspan="1" colspan="1">Nrn1</td><td align="left" rowspan="1" colspan="1">K230014L09</td><td align="left" rowspan="1" colspan="1">1.19E-2</td><td align="left" rowspan="1" colspan="1">−2.13</td><td align="left" rowspan="1" colspan="1">Evl</td><td align="left" rowspan="1" colspan="1">K230309E11</td><td align="left" rowspan="1" colspan="1">1.06E-2</td><td align="left" rowspan="1" colspan="1">1.55</td></tr><tr><td align="left" rowspan="1" colspan="1">Snca</td><td align="left" rowspan="1" colspan="1">K230001F08</td><td align="left" rowspan="1" colspan="1">1.07E-3</td><td align="left" rowspan="1" colspan="1">−2.02</td><td align="left" rowspan="1" colspan="1">Vip</td><td align="left" rowspan="1" colspan="1">K230041M20</td><td align="left" rowspan="1" colspan="1">1.25E-3</td><td align="left" rowspan="1" colspan="1">1.54</td></tr><tr><td align="left" rowspan="1" colspan="1">6330527O06Rik</td><td align="left" rowspan="1" colspan="1">K230003M18</td><td align="left" rowspan="1" colspan="1">2.02E-3</td><td align="left" rowspan="1" colspan="1">−2.00</td><td align="left" rowspan="1" colspan="1">Btbd14a</td><td align="left" rowspan="1" colspan="1">K330002G16</td><td align="left" rowspan="1" colspan="1">1.81E-2</td><td align="left" rowspan="1" colspan="1">1.48</td></tr><tr><td align="left" rowspan="1" colspan="1">6330403K07Rik</td><td align="left" rowspan="1" colspan="1">K230304J17</td><td align="left" rowspan="1" colspan="1">4.34E-4</td><td align="left" rowspan="1" colspan="1">−1.97</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">K330002G06</td><td align="left" rowspan="1" colspan="1">1.96E-2</td><td align="left" rowspan="1" colspan="1">1.48</td></tr><tr><td align="left" rowspan="1" colspan="1">D11Ertd707e</td><td align="left" rowspan="1" colspan="1">K230009N15</td><td align="left" rowspan="1" colspan="1">6.26E-5</td><td align="left" rowspan="1" colspan="1">−1.94</td><td align="left" rowspan="1" colspan="1">Cox6c</td><td align="left" rowspan="1" colspan="1">K230002F10</td><td align="left" rowspan="1" colspan="1">4.77E-3</td><td align="left" rowspan="1" colspan="1">1.46</td></tr><tr><td align="left" rowspan="1" colspan="1">Ptn</td><td align="left" rowspan="1" colspan="1">K230334M04</td><td align="left" rowspan="1" colspan="1">6.56E-4</td><td align="left" rowspan="1" colspan="1">−1.91</td><td align="left" rowspan="1" colspan="1">Pcdhac2</td><td align="left" rowspan="1" colspan="1">K230029N10</td><td align="left" rowspan="1" colspan="1">2.58E-3</td><td align="left" rowspan="1" colspan="1">1.45</td></tr><tr><td align="left" rowspan="1" colspan="1">Schip1</td><td align="left" rowspan="1" colspan="1">K230006C03</td><td align="left" rowspan="1" colspan="1">1.35E-3</td><td align="left" rowspan="1" colspan="1">−1.83</td><td align="left" rowspan="1" colspan="1">Atp5g3</td><td align="left" rowspan="1" colspan="1">K330002E13</td><td align="left" rowspan="1" colspan="1">3.25E-4</td><td align="left" rowspan="1" colspan="1">1.44</td></tr><tr><td align="left" rowspan="1" colspan="1">Baiap2</td><td align="left" rowspan="1" colspan="1">K230319I19</td><td align="left" rowspan="1" colspan="1">1.76E-2</td><td align="left" rowspan="1" colspan="1">−1.83</td><td align="left" rowspan="1" colspan="1">2010107H07Rik</td><td align="left" rowspan="1" colspan="1">K330010C07</td><td align="left" rowspan="1" colspan="1">1.20E-3</td><td align="left" rowspan="1" colspan="1">1.44</td></tr><tr><td align="left" rowspan="1" colspan="1">3110035E14Rik</td><td align="left" rowspan="1" colspan="1">K230010N04</td><td align="left" rowspan="1" colspan="1">2.01E-2</td><td align="left" rowspan="1" colspan="1">−1.82</td><td align="left" rowspan="1" colspan="1">Vamp1</td><td align="left" rowspan="1" colspan="1">K230001K20</td><td align="left" rowspan="1" colspan="1">2.35E-3</td><td align="left" rowspan="1" colspan="1">1.43</td></tr><tr><td align="left" rowspan="1" colspan="1">H2-DMa</td><td align="left" rowspan="1" colspan="1">K230007M04</td><td align="left" rowspan="1" colspan="1">3.35E-4</td><td align="left" rowspan="1" colspan="1">−1.78</td><td align="left" rowspan="1" colspan="1">9530058B02Rik</td><td align="left" rowspan="1" colspan="1">K230010H11</td><td align="left" rowspan="1" colspan="1">6.03E-3</td><td align="left" rowspan="1" colspan="1">1.40</td></tr><tr><td align="left" rowspan="1" colspan="1">Hn1</td><td align="left" rowspan="1" colspan="1">K230005H22</td><td align="left" rowspan="1" colspan="1">1.68E-4</td><td align="left" rowspan="1" colspan="1">−1.62</td><td align="left" rowspan="1" colspan="1">Cox6c</td><td align="left" rowspan="1" colspan="1">K230006O16</td><td align="left" rowspan="1" colspan="1">8.83E-3</td><td align="left" rowspan="1" colspan="1">1.37</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">K230011F01</td><td align="left" rowspan="1" colspan="1">1.23E-2</td><td align="left" rowspan="1" colspan="1">−1.61</td><td align="left" rowspan="1" colspan="1">1810012P15Rik</td><td align="left" rowspan="1" colspan="1">K230027M24</td><td align="left" rowspan="1" colspan="1">1.26E-2</td><td align="left" rowspan="1" colspan="1">1.37</td></tr><tr><td align="left" rowspan="1" colspan="1">Ngef</td><td align="left" rowspan="1" colspan="1">K330017A06</td><td align="left" rowspan="1" colspan="1">3.67E-3</td><td align="left" rowspan="1" colspan="1">−1.51</td><td align="left" rowspan="1" colspan="1">Stac2</td><td align="left" rowspan="1" colspan="1">K230027M01</td><td align="left" rowspan="1" colspan="1">1.92E-3</td><td align="left" rowspan="1" colspan="1">1.34</td></tr><tr><td align="left" rowspan="1" colspan="1">Gap43</td><td align="left" rowspan="1" colspan="1">K230009L16</td><td align="left" rowspan="1" colspan="1">3.03E-2</td><td align="left" rowspan="1" colspan="1">−1.44</td><td align="left" rowspan="1" colspan="1">AI836003</td><td align="left" rowspan="1" colspan="1">K230306I16</td><td align="left" rowspan="1" colspan="1">2.78E-6</td><td align="left" rowspan="1" colspan="1">1.33</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">K230302E22</td><td align="left" rowspan="1" colspan="1">2.48E-4</td><td align="left" rowspan="1" colspan="1">−1.42</td><td align="left" rowspan="1" colspan="1">Ccne1</td><td align="left" rowspan="1" colspan="1">K230044G09</td><td align="left" rowspan="1" colspan="1">1.36E-3</td><td align="left" rowspan="1" colspan="1">1.31</td></tr><tr><td align="left" rowspan="1" colspan="1">Fbxl2</td><td align="left" rowspan="1" colspan="1">K230024L20</td><td align="left" rowspan="1" colspan="1">5.40E-3</td><td align="left" rowspan="1" colspan="1">−1.42</td><td align="left" rowspan="1" colspan="1">Pdlim3</td><td align="left" rowspan="1" colspan="1">K230306A09</td><td align="left" rowspan="1" colspan="1">6.48E-3</td><td align="left" rowspan="1" colspan="1">1.31</td></tr><tr><td align="left" rowspan="1" colspan="1">Pcsk2</td><td align="left" rowspan="1" colspan="1">K230002M17</td><td align="left" rowspan="1" colspan="1">1.59E-3</td><td align="left" rowspan="1" colspan="1">−1.37</td><td align="left" rowspan="1" colspan="1">Atp1b1</td><td align="left" rowspan="1" colspan="1">K230006G24</td><td align="left" rowspan="1" colspan="1">1.56E-2</td><td align="left" rowspan="1" colspan="1">1.29</td></tr><tr><td align="left" rowspan="1" colspan="1">Ctsb</td><td align="left" rowspan="1" colspan="1">K230008F16</td><td align="left" rowspan="1" colspan="1">1.35E-2</td><td align="left" rowspan="1" colspan="1">−1.37</td><td align="left" rowspan="1" colspan="1">2610205H19Rik</td><td align="left" rowspan="1" colspan="1">K230306A08</td><td align="left" rowspan="1" colspan="1">7.83E-6</td><td align="left" rowspan="1" colspan="1">1.26</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003012-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003012.t004</object-id><label>Table 4</label><caption><title>Differential expression of the Atp5, Cox6 or Nduf family genes.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">CLONE_ID</td><td align="left" rowspan="1" colspan="1">SYMBOL</td><td align="left" rowspan="1" colspan="1">p</td><td align="left" rowspan="1" colspan="1">M</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">K230031N16</td><td align="left" rowspan="1" colspan="1">Ndufa11</td><td align="left" rowspan="1" colspan="1">7.69E-3</td><td align="left" rowspan="1" colspan="1">1.61</td></tr><tr><td align="left" rowspan="1" colspan="1">K230002F10</td><td align="left" rowspan="1" colspan="1">Cox6c</td><td align="left" rowspan="1" colspan="1">4.77E-3</td><td align="left" rowspan="1" colspan="1">1.46</td></tr><tr><td align="left" rowspan="1" colspan="1">K330002E13</td><td align="left" rowspan="1" colspan="1">Atp5g3</td><td align="left" rowspan="1" colspan="1">3.25E-4</td><td align="left" rowspan="1" colspan="1">1.44</td></tr><tr><td align="left" rowspan="1" colspan="1">K230006O16</td><td align="left" rowspan="1" colspan="1">Cox6c</td><td align="left" rowspan="1" colspan="1">8.83E-3</td><td align="left" rowspan="1" colspan="1">1.37</td></tr><tr><td align="left" rowspan="1" colspan="1">K230309I10</td><td align="left" rowspan="1" colspan="1">Cox6b1</td><td align="left" rowspan="1" colspan="1">4.73E-4</td><td align="left" rowspan="1" colspan="1">1.22</td></tr><tr><td align="left" rowspan="1" colspan="1">K230001M19</td><td align="left" rowspan="1" colspan="1">Atp5h</td><td align="left" rowspan="1" colspan="1">1.99E-6</td><td align="left" rowspan="1" colspan="1">1.22</td></tr><tr><td align="left" rowspan="1" colspan="1">K230340G22</td><td align="left" rowspan="1" colspan="1">Ndufb9</td><td align="left" rowspan="1" colspan="1">2.84E-3</td><td align="left" rowspan="1" colspan="1">1.16</td></tr><tr><td align="left" rowspan="1" colspan="1">K230002E05</td><td align="left" rowspan="1" colspan="1">Ndufb8</td><td align="left" rowspan="1" colspan="1">2.84E-4</td><td align="left" rowspan="1" colspan="1">1.15</td></tr><tr><td align="left" rowspan="1" colspan="1">K230006I15</td><td align="left" rowspan="1" colspan="1">Ndufa8</td><td align="left" rowspan="1" colspan="1">5.80E-3</td><td align="left" rowspan="1" colspan="1">1.14</td></tr><tr><td align="left" rowspan="1" colspan="1">K230009H06</td><td align="left" rowspan="1" colspan="1">Ndufs5</td><td align="left" rowspan="1" colspan="1">3.36E-3</td><td align="left" rowspan="1" colspan="1">1.13</td></tr><tr><td align="left" rowspan="1" colspan="1">K230012O10</td><td align="left" rowspan="1" colspan="1">Ndufa2</td><td align="left" rowspan="1" colspan="1">2.13E-2</td><td align="left" rowspan="1" colspan="1">1.13</td></tr><tr><td align="left" rowspan="1" colspan="1">K230310N23</td><td align="left" rowspan="1" colspan="1">Ndufb7</td><td align="left" rowspan="1" colspan="1">9.34E-4</td><td align="left" rowspan="1" colspan="1">1.11</td></tr><tr><td align="left" rowspan="1" colspan="1">K330022I20</td><td align="left" rowspan="1" colspan="1">Ndufa6</td><td align="left" rowspan="1" colspan="1">4.48E-4</td><td align="left" rowspan="1" colspan="1">1.11</td></tr><tr><td align="left" rowspan="1" colspan="1">K230010M09</td><td align="left" rowspan="1" colspan="1">Atp5h</td><td align="left" rowspan="1" colspan="1">4.84E-5</td><td align="left" rowspan="1" colspan="1">1.09</td></tr><tr><td align="left" rowspan="1" colspan="1">K230310M07</td><td align="left" rowspan="1" colspan="1">Ndufa5</td><td align="left" rowspan="1" colspan="1">1.50E-2</td><td align="left" rowspan="1" colspan="1">1.09</td></tr><tr><td align="left" rowspan="1" colspan="1">K230044G12</td><td align="left" rowspan="1" colspan="1">Ndufc2</td><td align="left" rowspan="1" colspan="1">1.06E-3</td><td align="left" rowspan="1" colspan="1">1.06</td></tr><tr><td align="left" rowspan="1" colspan="1">K230319K23</td><td align="left" rowspan="1" colspan="1">Atp5g1</td><td align="left" rowspan="1" colspan="1">3.37E-3</td><td align="left" rowspan="1" colspan="1">1.04</td></tr><tr><td align="left" rowspan="1" colspan="1">K230012P07</td><td align="left" rowspan="1" colspan="1">Atp5j</td><td align="left" rowspan="1" colspan="1">3.80E-4</td><td align="left" rowspan="1" colspan="1">1.02</td></tr><tr><td align="left" rowspan="1" colspan="1">K230304A02</td><td align="left" rowspan="1" colspan="1">Ndufa5</td><td align="left" rowspan="1" colspan="1">8.35E-3</td><td align="left" rowspan="1" colspan="1">0.96</td></tr><tr><td align="left" rowspan="1" colspan="1">K230019O17</td><td align="left" rowspan="1" colspan="1">Ndufs2</td><td align="left" rowspan="1" colspan="1">2.80E-4</td><td align="left" rowspan="1" colspan="1">0.95</td></tr><tr><td align="left" rowspan="1" colspan="1">K230304I19</td><td align="left" rowspan="1" colspan="1">Ndufb10</td><td align="left" rowspan="1" colspan="1">7.66E-4</td><td align="left" rowspan="1" colspan="1">0.95</td></tr><tr><td align="left" rowspan="1" colspan="1">K330018J12</td><td align="left" rowspan="1" colspan="1">Atp5e</td><td align="left" rowspan="1" colspan="1">1.94E-3</td><td align="left" rowspan="1" colspan="1">0.92</td></tr><tr><td align="left" rowspan="1" colspan="1">K230039K18</td><td align="left" rowspan="1" colspan="1">Atp5c1</td><td align="left" rowspan="1" colspan="1">2.06E-4</td><td align="left" rowspan="1" colspan="1">0.89</td></tr><tr><td align="left" rowspan="1" colspan="1">K230050F20</td><td align="left" rowspan="1" colspan="1">Ndufb2</td><td align="left" rowspan="1" colspan="1">2.15E-3</td><td align="left" rowspan="1" colspan="1">0.85</td></tr><tr><td align="left" rowspan="1" colspan="1">K230053K02</td><td align="left" rowspan="1" colspan="1">Atp5o</td><td align="left" rowspan="1" colspan="1">1.98E-4</td><td align="left" rowspan="1" colspan="1">0.83</td></tr><tr><td align="left" rowspan="1" colspan="1">K230014B16</td><td align="left" rowspan="1" colspan="1">Atp5c1</td><td align="left" rowspan="1" colspan="1">6.17E-4</td><td align="left" rowspan="1" colspan="1">0.81</td></tr><tr><td align="left" rowspan="1" colspan="1">K230003F11</td><td align="left" rowspan="1" colspan="1">Ndufs8</td><td align="left" rowspan="1" colspan="1">4.23E-3</td><td align="left" rowspan="1" colspan="1">0.80</td></tr><tr><td align="left" rowspan="1" colspan="1">K230304N08</td><td align="left" rowspan="1" colspan="1">Ndufv2</td><td align="left" rowspan="1" colspan="1">3.12E-4</td><td align="left" rowspan="1" colspan="1">0.79</td></tr><tr><td align="left" rowspan="1" colspan="1">K330022A12</td><td align="left" rowspan="1" colspan="1">Ndufv1</td><td align="left" rowspan="1" colspan="1">7.91E-3</td><td align="left" rowspan="1" colspan="1">0.78</td></tr><tr><td align="left" rowspan="1" colspan="1">K230008K02</td><td align="left" rowspan="1" colspan="1">Ndufa10</td><td align="left" rowspan="1" colspan="1">1.25E-2</td><td align="left" rowspan="1" colspan="1">0.76</td></tr><tr><td align="left" rowspan="1" colspan="1">K230001C13</td><td align="left" rowspan="1" colspan="1">Cox6a1</td><td align="left" rowspan="1" colspan="1">7.63E-3</td><td align="left" rowspan="1" colspan="1">0.66</td></tr><tr><td align="left" rowspan="1" colspan="1">K330009A17</td><td align="left" rowspan="1" colspan="1">Atp5d</td><td align="left" rowspan="1" colspan="1">1.36E-2</td><td align="left" rowspan="1" colspan="1">0.60</td></tr><tr><td align="left" rowspan="1" colspan="1">K230053L04</td><td align="left" rowspan="1" colspan="1">Ndufs3</td><td align="left" rowspan="1" colspan="1">5.61E-3</td><td align="left" rowspan="1" colspan="1">0.49</td></tr><tr><td align="left" rowspan="1" colspan="1">K230005J19</td><td align="left" rowspan="1" colspan="1">Ndufaf1</td><td align="left" rowspan="1" colspan="1">5.43E-3</td><td align="left" rowspan="1" colspan="1">0.32</td></tr><tr><td align="left" rowspan="1" colspan="1">K230034B16</td><td align="left" rowspan="1" colspan="1">Atp5g2</td><td align="left" rowspan="1" colspan="1">1.27E-1</td><td align="left" rowspan="1" colspan="1">0.27</td></tr><tr><td align="left" rowspan="1" colspan="1">K230003B02</td><td align="left" rowspan="1" colspan="1">Ndufs1</td><td align="left" rowspan="1" colspan="1">8.05E-2</td><td align="left" rowspan="1" colspan="1">0.22</td></tr><tr><td align="left" rowspan="1" colspan="1">K230003P04</td><td align="left" rowspan="1" colspan="1">Atp5f1</td><td align="left" rowspan="1" colspan="1">4.42E-1</td><td align="left" rowspan="1" colspan="1">0.16</td></tr><tr><td align="left" rowspan="1" colspan="1">K230018C04</td><td align="left" rowspan="1" colspan="1">Atp5b</td><td align="left" rowspan="1" colspan="1">6.76E-1</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">K230001F17</td><td align="left" rowspan="1" colspan="1">Atp5a1</td><td align="left" rowspan="1" colspan="1">6.43E-1</td><td align="left" rowspan="1" colspan="1">0.07</td></tr><tr><td align="left" rowspan="1" colspan="1">K230302C16</td><td align="left" rowspan="1" colspan="1">Atp5s</td><td align="left" rowspan="1" colspan="1">7.52E-1</td><td align="left" rowspan="1" colspan="1">0.04</td></tr><tr><td align="left" rowspan="1" colspan="1">K230001C22</td><td align="left" rowspan="1" colspan="1">Atp5b</td><td align="left" rowspan="1" colspan="1">5.26E-1</td><td align="left" rowspan="1" colspan="1">−0.21</td></tr></tbody></table></alternatives></table-wrap>
<inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula>
<supplementary-material content-type="local-data" id="pone.0003012.s001"><label>Figure S1</label><caption><p>Annotation pipeline (databases).</p><p>(0.06 MB PNG)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003012.s002"><label>Figure S2</label><caption><p>Annotation pipeline (scores).</p><p>(0.09 MB PNG)</p></caption></supplementary-material>
<fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>This work was supported by a grant of the 6th Framework of the EU commission to the NFG consortium (to PC), a RIKEN presidential grant for collaborative research (to PC and TKH), a Grant-in-Aid for Scientific Research on Priority Areas - Integrative Brain Research - from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan (to TKH), a Japanese Society for the Promotion of Science (JSPS) grant. CP was supported by a fellowship from the JSPS (number P05880). HN is supported by Grant-in-Aids on Priority Area Research (C) from MEXT (number 18020034), and for Young Scientists (A) from JSPS, Japan (number 17680029).</p></fn></fn-group>
<graphic xlink:href="pone.0003012.g001"/> <graphic id="pone-0003012-t001-1" xlink:href="pone.0003012.t001"/> <graphic xlink:href="pone.0003012.g002"/> <graphic xlink:href="pone.0003012.g003"/> <graphic id="pone-0003012-t002-2" xlink:href="pone.0003012.t002"/> <graphic id="pone-0003012-t003-3" xlink:href="pone.0003012.t003"/> <graphic id="pone-0003012-t004-4" xlink:href="pone.0003012.t004"/> <graphic xlink:href="pone.0003012.g004"/> <inline-graphic xlink:href="pone.0003012.e001.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003012.e002.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003012.e003.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003012.e004.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003012.e005.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003012.e006.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003012.e007.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003012.e008.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003012.e009.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003012.e010.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003012.e011.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003012.e012.jpg" mimetype="image"/>
<media xlink:href="pone.0003012.s001.png"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003012.s002.png"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["38"], "element-citation": ["\n"], "collab": ["R Development Core Team"], "year": ["2007"], "article-title": ["R: A Language and Environment for Statistical Computing, Vienna, Austria."], "comment": ["ISBN 3-900051-07-0, URL. "], "ext-link": ["http://www.R-project.org"]}, {"label": ["39"], "element-citation": ["\n"], "surname": ["Cleverland", "Devlin"], "given-names": ["W", "S"], "year": ["1988"], "article-title": ["Locally weighted regression: An approach to regression analysis by local fitting."], "source": ["J Amer Statist Assoc"], "volume": ["83"], "fpage": ["596"], "lpage": ["610"]}, {"label": ["41"], "element-citation": ["\n"], "surname": ["Storey"], "given-names": ["JD"], "year": ["2003"], "article-title": ["The positive false discovery rate: a Bayesian interpretation and the q-value."], "source": ["Ann Statist"], "volume": ["31"], "fpage": ["2013"], "lpage": ["2035"]}]
PLoS One. 2008 Aug 20; 3(8):e3012