Carcinogenesis Advance Access originally published online on August 27, 2007
Carcinogenesis 2007 28(12):2484-2490; doi:10.1093/carcin/bgm194
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Genome-wide detection of testis- and testicular cancer-specific alternative splicing
Department of Genetics and Center for Developmental Biology, College of Life Sciences, Wuhan University, Wuhan 430072, P. R. China
* To whom correspondence should be addressed. Fax: +86 27 68756253; Email: rjzhou{at}whu.edu.cn Correspondence may also be addressed to Hanhua Cheng. Tel: +86 27 87889294; Fax: +86 27 68756253;Email: hhcheng{at}whu.edu.cn
| Abstract |
|---|
|
|
|---|
Alternative pre-messenger RNA (mRNA) splicing is a key molecular event that allows for protein diversity and plays important roles in development and disease. Alternative pre-mRNA splicing regulations during spermatogenesis and alternative pre-mRNA splicing etiology in testicular tumorigenesis are yet to be characterized. By genome-wide analysis, here we describe alternative splicing features that distinguish distinctive patterns of alternative pre-mRNA splicing among human testis, testicular cancer and mouse testis. Through computationally subtractive analysis, we detected 80 testis-specific transcript candidates in human testis, 175 in human testicular cancer and 262 in mouse testis, which were integrated into a database. Reverse transcription–polymerase chain reaction confirmed that most of these transcript candidates from mouse testis were testis specific. Around 40% of the transcripts were from unknown/hypothetical genes, which were useful for further functional analysis. These transcripts were not overlapped, indicating lack of evolutionary conservation. Further chromosome mapping showed distinct chromosomal preference of alternative pre-mRNA splicing events. Comparison analysis indicated that alternative pre-mRNA splicing in human testicular tumor shared some characters/trends with those in mouse testis. Moreover, human testicular tumor tended to use rare splice sites and there were also distinct sequences adjacent dominant splice sites between normal testis and testicular tumor. These special features of alternative pre-mRNA splicing in human testicular tumor suggested that testicular tumorigenesis was involved in multiple steps/levels of alternative splicing events. Using alternative splicing as a potential source for new clinical diagnostic, prognostic and therapeutic strategies for treatment of testicular tumors seems to have a bright prospect.
Abbreviations: EST, expressed sequence tag; mRNA, messenger RNA; PCR, polymerase chain reaction; RT, reverse transcription
| Introduction |
|---|
|
|
|---|
Since the discovery that the coding information of metazoan genes is interrupted by introns (1,2), it still amazes how both intron removal and exon rearrangement are regulated precisely to produce correct proteomes in cell type- or developmental stage-specific manner. Alternative splicing, the process by which the exons of primary transcripts [pre-messenger RNAs (mRNAs)] from genes can be spliced in different arrangements to produce structurally and functionally distinct mRNA and protein variants, is the most widely used mechanism to enhance protein diversity of higher eukaryotic organisms. It was estimated that
35–70% of all human genes appear to undergo alternative splicing (3–8). Comparable amounts of alternative splicing were also detected in the fly, nematode and other vertebrates (8,9), emphasizing its major role in the biology of all animals. In certain pathological conditions, aberrantly spliced pre-mRNAs can be generated that go unnoticed by the nonsense-mediated mRNA decay pathway, and are, therefore, translated into aberrant proteins, which are involved in several human diseases (10–12). There have been also multifaceted studies into the relevance of alternative splicing to cancer at the level of individual genes and individual cancers (13). Currently, the analysis of cancer-specific alternative splicing is a promising step forward in potential source for new clinical diagnostic, prognostic and therapeutic strategies.
However, despite the growing interest on the impact of alternative splicing in various aspects of the biological processes, our understanding of alternative splicing is still very primitive and its regulatory mechanisms are mostly unknown (14,15). The unusually high levels of alternative splicing were often observed in the human brain and testis (16,17). Why is alternative splicing so abundant in the testis? Alternative splicing is probably very important throughout germ cell development. In order to advance our understanding of the biological significance of alternative splicing in humans testis and its tumorigenesis, it is essential to identify tissue-specific splicing events. Although tissue-specific alternative splicing patterns was analyzed by computationally determined splice (8,15,16), comparisons between normal and tumor, especially in human testis, combined with experimentally determined analysis, will provide more information on our understanding of human testicular tumorigenesis.
Here, we performed a genome-wide analysis of alternative splicing, combining with computationally subtractive screening and experimental verification, to compare tissue-specific alternative splicing patterns among human normal testis and testicular cancer and mouse testis. Our results distinguish distinctive patterns of alternative splicing and identify pronounced differences in the chromosomal distribution, splice types and splice-/flanking-splice sites among the three kinds of tissues, suggesting multiple and even complex role of alternative splicing in human testicular tumorigenesis.
| Materials and methods |
|---|
|
|
|---|
Data sources and filtration
Expressed sequence tag (EST) data of human and mouse were drawn from Cancer Genome Anatomy Project (CGAP) project (http://cgap.nci.nih.gov/Tissues/LibraryFinder). CGAP project collects EST libraries from all over the world and gives good tissue information. We downloaded all available EST libraries of human and mouse from CGAP libraries, Mammalian Gene Collection libraries and Open Reading frame EST Sequencing libraries. We sought to avoid mixing of multiple tissues. Among these libraries that signed pooled, normalized or subtracted were excluded because these procedures will affect tissue classification. All EST data in different tissues used were listed in Table I. Human testis libraries were mainly from adults (eight), five were unknown age and one from embryo. Testicular cancers were mainly seminoma and embryonal carcinoma. Mouse testis libraries were from adult, 15 day old and unknown age. All other libraries of the tissues were used regardless of age.
|
All collection data were then dealt with three procedures: repeat sequence masking (program, repeatmasker; repeat database, repbase; girnst server: www.girinst.org) to remove simple repeat in the dataset, vector and contamination masking (program, crossmatch; vector database, UniVec_Core; National Center for Biotechnology Information ftp server: ftp://ftp.ncbi.nih.gov/) to clean the vector sequences and a final cleaning of short and rubbish sequences (program seqclean from egassembler server: http://egassembler.hgc.jp). Alu repeats were included in following analysis. The filtered ESTs were available for the following analysis.
Computational procedures to identify alternative splicing
A local Basic Local Alignment Search Tool database was constructed for ESTs of each tissue. Testis-specific transcripts were identified based on a subtractive model showed in Figure 1a. Briefly, there are following three steps: (i) Testis EST dataset was blasted to itself. The e-value was set as 1e-30. Gaps (insertion and deletion) of ESTs were identified after alignments. Parameters to identify alternative splicing: the gap length, >10 bp; nucleotide identity, >95%. (ii) Testis ESTs were blasted to the ESTs from the other tissues. Parameters were the same as step (i). (iii) Subtractive ESTs were identified as testis-specific ESTs by insertion/deletion comparisons after the Basic Local Alignment Search Tool. Computer programs were written using Perl language.
|
EST/genomic sequence alignments, chromosome mapping and splice site analysis
To decrease errors in EST alignments and determine the chromosomal loci of each gene, we localized ESTs to genomic sequences using BLAST-Like Alignment Tool tools (http://genome.ucsc.edu). We used the default parameters and selected the best score results. The exon positions of each transcript on the chromosome were recorded and used to determine splice sites and gene structure. Splice sites of both 5' and 3' exon/intron boundaries were aligned online via http://weblogo.berkeley.edu/logo.cgi. We allowed an error of 10 bp in exon/intron boundary. Based on comparisons of EST/genomic alignment, two possible errors can be checked out: (i) candidate EST in the same gene was not on the same chromosome; (ii) candidate EST in the same gene was not in the same locus on the chromosome. The reasons that caused these errors mainly included EST sequencing error, pseudogene and multiple copy genes. The two cases were excluded as false positive in final database.
Reverse transcription–polymerase chain reaction
Reverse transcription–polymerase chain reaction (RT–PCR) was used to amplify transcripts. Total RNAs from mouse and human tissues were prepared by the RNeasy Mini kit (Qiagen, Shanghai, China) according to the manufacturer's instructions. Human testicular cancer sample (seminoma) was provided by the Clinic of Wuhan University Zhongnan Hospital. All the RNAs were digested by RNase-free DNase I and purified. About 3 µg RNAs were used as template for RT using 0.5 µg poly(T)20 primer and 200 U Moloney Murine Leukemia Virus reverse transcriptase (Promega, Madison, WI, USA). PCR was performed in a 20 µl reaction mix containing 10 mM Tris–HCl, pH 8.3, 1.5 mM MgCl2, 50 mM KCl, 150 mM deoxynucleoside triphosphate, 0.2 µM each primer and 1 U Taq DNA polymerase. Amplification conditions are 94°C, 30 s; 60–68°C (supplementary Table S1 is available at Carcinogenesis Online), 30 s and 72°C, 30 s for 35 cycles. All PCR products were resolved on 1.5% agarose gels run at 90 V/20 cm for 0.5 h in Tris-Borate-EDTA buffer.
| Results |
|---|
|
|
|---|
Detection of testis-specific alternative splicing
For analyzing tissue specificity, we have carefully classified all available EST libraries of human and mouse from CGAP, Mammalian Gene Collection and Open Reading frame EST Sequencing libraries into 19 distinct tissue classes to avoid mixing of multiple tissues. Our final classification consisted of human normal and tumor samples (4 353 301 ESTs) and mouse normal tissues (2 070 996 ESTs; Table I). Through computationally subtractive analysis, we detected 80 testis-specific transcript candidates in human testis, 175 in human testicular cancer and 262 in mouse testis. These transcripts belonged to 72 genes in human testis, 145 in human testicular cancer and 212 in mouse testis. Around 87% of these genes contained at least one testis-specific transcript candidate,
10% of them had two and a few of them had three to six transcripts (Figure 1b). We also analyzed overlapping among detected genes/transcripts in the three samples and found that there were only a few genes overlapped. To see if there was chromosome preference of alternative splicing events, we further mapped these transcripts onto chromosomes. Based on chromosome locations of transcripts, we observed that alternative splicing events in human normal testis were rich on chromosome 1, 12, 17 and 20; however, in testis cancer, they displayed strong preference to chromosome 2, 11, 12, 19 and X (Figure 1c). In mouse testis, alternative splicing events occurred mainly on chromosome 2, 7 and 17, although distributed on all chromosomes (Figure 1d).
To facilitate future study and reference of alternatively spliced genes of mammals, we constructed a database based on our analysis, which contained 80 tissue-specific transcripts in human testis, 175 in human testicular cancer and 262 in mouse testis, of which, 40% were unknown/hypothetical genes. It gave an overview of testis-specific alternative splicing in human and mouse and was especially useful for further study on testis oncogenesis. The main information of this database includes gene ID, gene alternative splicing patterns, gene structures, locations in chromosomes, sequences, DNA and protein sequences linked with the National Center for Biotechnology Information website. An example page showed the details of gene AI507611 in mouse (Figure 2). The database can be accessed at http://202.114.72.2:8080/input.jsp (password: 123454321).
|
To validate alternative splicing detected by the computational approaches, we designed RT–PCR primers to amplify testis-specific transcripts of mouse whose function information was not clear in public databases. Of 62 genes amplified, 26 transcripts were testis specific (Figure 3a) and 23 were dominantly expressed in testis (Figure 3b). The two groups accounted for 79% of all cases. Another eight genes (16%) were ubiquitously expressed in the tested tissues (Figure 3c). The remaining five genes cannot be amplified in our validation. Reasons for this may be related to EST sequencing errors, splicing errors or improper primers. Predicted transcripts specific for testicular cancer were verified by RT–PCR, which showed that some were testis cancer specific (Figure 3e1) and some were expressed at equal or lower level in testis cancer (Figure 3e2).
|
Patterns of the identified alternative splicing in testis
After comparisons between EST and genomic sequences on chromosomes, we detected five types of alternative splicing patterns: cassette alternative exon, alternative 5' splice site, alternative 3' splice site, intron retention and mutually exclusive alternative exons (Figure 4a). It was intriguing that the cassette alternative exon variants formed the most frequent category (50% of all the types). For example, by subtractive analysis, we detected that loss of the exon III by cassette alternative exon of mouse Tcte3 (t-complex-associated testis expressed 3) gene resulted in a testis-specific transcript 2 which deleted Tctex-1 domain (Figure 4c). Similarly, deletion of the exon VI by cassette alternative exon of human TMBIM4 (trans-membrane BAX inhibitor motif containing 4) gene generated a testis-specific short form of transcript 2, which lost the 3' part of the domain UPF0005 (Figure 4d). Alternative splicing may cause testis cancers. In support of this, we detected a testis cancer-specific transcript 1 of human PRDX3 (peroxiredoxin 3) gene, which was produced by the alternative 5' and 3' splice sites (Figure 4e). Indeed, we detected >110 of human testis cancer-specific alternatively splicing transcripts.
|
Comparisons of ratios of these types of alternative splicing patterns among human testis, testicular cancer and mouse testis showed that both alternative 3' splice site and intron retention had a higher proportion in cancer, however ratio of both cassette alternative exon and alternative 5' splice site decreased in cancer. In addition, in mouse testis, the type alternative 5' splice site was highest, whereas the cassette alternative exon was lowest among the three samples (Figure 4b). These results showed that there was significant preference in use of these types of alternative splicing in testicular cancer and mouse testis.
Diversification of the identified alternative splice sites
To see preference/diversification of alternative splice sites in these testis-expressed genes, we analyzed 1747 splice sites in human testis, 5159 in human testicular cancer and 7002 in mouse testis by comparison of each EST with its genomic sequence and mapping it onto chromosome. We detected four basic donor–acceptor splice sites: GT–AG, GC–AG, GG–AG and GT–GG, of which GT–AG was the most dominant site (90–93%; Figure 5a). Furthermore, several interesting points we found were as follows: (i) GG–AG was detected in both normal (0.38–0.47%) and cancer testis (0.74%). (ii) GT–GG was only detected in cancer (0.45%) and mouse testis (0.19%), not in normal human testis. (iii) Besides the four basic donor–acceptor splice sites, there were quite a number of rare sites especially in testicular cancer (6.31%). (iv) In all, 27% of the rare sites were shared by human testicular cancer and mouse testis, 21% of them were human testis cancer specific, whereas only a few of them were normal human testis specific (4%; Figure 5b). These results suggested that the rare splice sites may be involved in testis tumorigenesis; moreover, most of these rare splice sites were mouse derived.
|
We further compared sequences adjacent dominant splice sites among human testis cancer, testicular cancer and mouse testis. In 5' splice site region, downstream dominant splice site GT, though a consensus AAGT at position 3–6, less A at the position 3 (60%) and more AG at the position 4 (74%) and 5 (78%) were observed in testicular cancer compared with relevant positions in normal human testis (70%, 65% and 70%, respectively; Figure 6a). In the upstream of the 3' splice site AG, no preference of nucleotide at the –2 position and a C/T preference at the –1 position were observed in the three kinds of samples (Figure 6b). However, obvious T preference from position –3 to –9 was detected in testicular cancer, especially at the position –9 (64%), while in normal human testis the positions with the T preference was wider upstream position –3, especially less T (57%) at the position –9 and more T at the position 10 (64%; cancer, 42%) and 11 (63%; cancer, 50%), and less obvious T preference in mouse testis was observed.
|
| Discussion |
|---|
|
|
|---|
We have developed a method for precisely delineating alternative splice patterns by alignment EST/cDNA sequences first, and comparing the ESTs/cDNAs with their genomic sequences, and then mapping them onto chromosomes. These procedures eliminated many EST errors, pseudogene and multiple copy/repeat gene problems when data were from diverse EST databases. Considering the fact that genome-wide studies found transposable elements (mainly Alus) in nearly 4% of human protein-coding regions (18), the Alu repeats were not excluded in our subtractive analysis. The alternatively spliced transcripts were subject to subtractive screening of testis versus all the other tissues, and these analyses finally yielded testis-specific transcripts. Interestingly,
40% of the transcripts were from unknown/hypothetical genes, which were useful for further functional analysis. RT–PCR analysis further confirmed most of our computationally determined results. Nevertheless, there were a small number of transcripts from computationally subtractive results, which were ubiquitously expressed in the tested tissues. This may be caused by EST coverage, their fragmentary character and sequencing errors escaped from our analysis. In any case, we believe our database of 80 testis-specific transcript candidates in human testis, 175 in human testicular cancer and 262 in mouse testis can be a rich source of discovery for our studying spermatogenesis and especially testicular tumorigenesis. However, its major value is providing previously unknown patterns of alternative splicing in human testicular tumor. (i) Alternative splicing events in testicular tumor have chromosomal preference, a strong preference to chromosome 2, 11, 12, 19 and X. (ii) Alternative splicing in human testicular tumor shares some characters/trends with those in mouse testis. For example, in the usage of types of alternative splicing, both human testicular cancer and mouse testis tend to use types of both alternative 3' splice site and intron retention. Moreover, a considerable number (27%) of the rare splice sites are shared by human testicular cancer and mouse testis, and genes involved in synthesis/degradation are frequently spliced in these two groups. (iii) Human testicular tumor tends to use rare splice sites. (iv) There are distinct sequences adjacent dominant splice sites in human testicular tumor. These special features of alternative splicing in human testicular tumor suggest that testicular tumorigenesis is involved in multiple steps/levels of alternative splicing events. This view is consistent with results observed in testicular germ cell tumorigenesis, which is involved in multiple aspects including somatic genetic changes, allelic imbalance, multiple gene mutations and methylation patterns (19). Testicular cancer seems to need multiple genetic events to occur. Taken together, these data suggest a complex etiopathology of testicular tumorigenesis. Our data include discovery of many novel splice forms of testicular cancer-associated genes, alternative splicing patterns of testicular cancer and suggest a significant new direction for testicular cancer research.
Interestingly, these tissue-specific alternative splicing variants are not evolutionarily conserved between human and mouse. Whether would the biological significance of their existence questionable? The same results were obtained to show that vast majority of alternative splicing events represented in the human EST database are not conserved in mouse (20). Main explanation for this issue may be because most of these splicing events may represent tissue-, disease- or allele-specific splicing, or events for which phylogenetic distribution is highly restricted. And alternative splicing could be also the easiest road to diversification of the species. We believe that species-specific alternative splicing represents a fascinating aspect of it that has undoubtedly had a major impact on genome evolution (21). Moreover, analyses of substitution rates in the genomes of different mammalian species indicate that the exons of species-specific splice variants evolve more rapidly than constitutive or conserved alternatively spliced exons (22). Indeed, it will be important in future studies to decipher the interplay among species-specific alternative splicing, protein diversity and pathogenic alternative splicing in cancer.
Main interest of this paper is data mining for testis-specific transcripts by the genome-wide identification and subtractive analysis of alternative splicing variants. Some of the testis cancer-specific splice forms we detected may be medically important. However, this should be further verified in detail by experimental approaches such micro-array analysis and functional analysis, so as to determine which are initial triggers of tumorigenesis or just reflections of transformed state, because at present time, some factors such as grade of aggressiveness, aging and stress were not considered owing to limitation of data sources. Indeed these factors affected the alternative splicing events (23–26). Computational analyses have gone some way to providing a global view of tissue-specific splicing, but these studies are limited by inadequate transcript coverage, by difficulties in assigning ESTs to tissue or cell types or sorting different categories.
In conclusion, evidences are accumulating that just as alternative splicing is important for spermatogenesis, so aberrations of alternative splicing patterns are important for testicular cancer. Testicular tumorigenesis needs multiple genetic events to occur; therefore, processes at the top of genetic hierarchies such as alternative splicing could be the cause, and the prospect that we might start to use alternative splicing as a potential source for new clinical diagnostic, prognostic and therapeutic strategies. For example, novel therapeutic strategies could be designed to target the splicing machinery, which can be defective in some testicular tumors. Further efforts into alternative splicing etiology in testicular cancer seem bright.
| Supplemetary material |
|---|
|
|
|---|
Supplementary Table S1 can be found at http://carcin.oxfordjournals.org/
| Funding |
|---|
|
|
|---|
National Natural Science Foundation of China; the National Key Basic Research project (2006CB102103); the Program for New Century Excellent Talents in University; 111 project (B06018).
| Acknowledgments |
|---|
The authors thank Dr J.Feng for help in writing software.
Conflict of Interest Statement: None declared.
| References |
|---|
|
|
|---|
- Berget SM, et al. Spliced segments at the 5' terminus of adenovirus 2 late mRNA. Proc. Natl Acad. Sci. USA (1977) 74:3171–3175.
[Abstract/Free Full Text] - Chow LT, et al. An amazing sequence arrangement at the 5' ends of adenovirus 2 messenger RNA. Cell (1977) 12:1–8.[CrossRef][Web of Science][Medline]
- Mironov AA, et al. Frequent alternative splicing of human genes. Genome Res. (1999) 9:1288–1293.
[Abstract/Free Full Text] - Brett D, et al. EST comparison indicates 38% of human mRNAs contain possible alternative splice forms. FEBS Lett. (2000) 474:83–86.[CrossRef][Web of Science][Medline]
- Kan Z, et al. Selecting for functional alternative splices in ESTs. Genome Res. (2002) 12:1837–1845.
[Abstract/Free Full Text] - Modrek B, et al. A genomic view of alternative splicing. Nat. Genet. (2002) 30:13–19.[CrossRef][Web of Science][Medline]
- Johnson JM, et al. Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays. Science (2003) 302:2141–2144.
[Abstract/Free Full Text] - Kim N, et al. The ASAP II database: analysis and comparative genomics of alternative splicing in 15 animal species. Nucleic Acids Res. (2007) 35:D93–D98.
[Abstract/Free Full Text] - Brett D, et al. Alternative splicing and genome complexity. Nat. Genet. (2002) 30:29–30.[CrossRef][Web of Science][Medline]
- Caceres JF, et al. Alternative splicing: multiple control mechanisms and involvement in human disease. Trends Genet. (2002) 18:186–193.[CrossRef][Web of Science][Medline]
- Faustino NA, et al. Pre-mRNA splicing and human disease. Genes Dev. (2003) 17:419–437.
[Free Full Text] - Garcia-Blanco MA, et al. Alternative splicing in disease and therapy. Nat. Biotechnol. (2004) 22:535–546.[CrossRef][Web of Science][Medline]
- Venables JP. Unbalanced alternative splicing and its significance in cancer. Bioessays (2006) 28:378–386.[CrossRef][Web of Science][Medline]
- Ladd AN, et al. Finding signals that regulate alternative splicing in the post-genomic era. Genome Biol. (2002) 3:reviews0008.[Medline]
- Takeda J, et al. Large-scale identification and characterization of alternative splicing variants of human gene transcripts using 56,419 completely sequenced and manually annotated full-length cDNAs. Nucleic Acids Res. (2006) 34:3917–3928.
[Abstract/Free Full Text] - Xu Q, et al. Genome-wide detection of tissue-specific alternative splicing in the human transcriptome. Nucleic Acids Res. (2002) 30:3754–3766.
[Abstract/Free Full Text] - Yeo G, et al. Variation in alternative splicing across human tissues. Genome Biol. (2004) 5:R74.[CrossRef][Medline]
- Nekrutenko A, et al. Transposable elements are found in a large number of human protein-coding genes. Trends Genet. (2001) 17:619–621.[CrossRef][Web of Science][Medline]
- Horwich A, et al. Testicular germ-cell cancer. Lancet (2006) 367:754–765.[CrossRef][Web of Science][Medline]
- Yeo GW, et al. Identification and analysis of alternative splicing events conserved in human and mouse. Proc. Natl Acad. Sci. USA (2005) 102:2850–2855.
[Abstract/Free Full Text] - Blencowe BJ. Alternative splicing: new insights from global analyses. Cell (2006) 126:37–47.[CrossRef][Web of Science][Medline]
- Xing Y, et al. Evidence of functional selection pressure for alternative splicing events that accelerate evolution of protein subsequences. Proc. Natl Acad. Sci. USA (2005) 102:13526–13531.
[Abstract/Free Full Text] - Dias CS, et al. Regulation of hdm2 by stress-induced hdm2alt1 in tumor and nontumorigenic cell lines correlating with p53 stability. Cancer Res. (2006) 66:9467–9473.
[Abstract/Free Full Text] - Chandler DS, et al. Genotoxic stress induces coordinately regulated alternative splicing of the p53 modulators MDM2 and MDM4. Cancer Res. (2006) 66:9502–9508.
[Abstract/Free Full Text] - Rossi EL. Stress-induced alternative gene splicing in mind-body medicine. Adv. Mind Body Med. (2004) 20:12–19.[Medline]
- Chang SY, et al. Age and gender-dependent alternative splicing of P/Q-type calcium channel EF-hand. Neuroscience (2007) 145:1026–1036.[CrossRef][Web of Science][Medline]
![]()
CiteULike
Connotea
Del.icio.us What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||





