- Open Access
Development of a rapid and inexpensive method to reveal natural antisense transcripts
© Collani and Barcaccia; licensee BioMed Central Ltd. 2012
Received: 10 May 2012
Accepted: 31 August 2012
Published: 12 September 2012
Natural antisense transcripts (NATs) are a group of RNAs encoded within a cell that have transcript complementarity to other RNA transcripts. NATs have been identified in multiple eukaryotes, including humans, mice, yeast and several plants, and are known to play crucial roles in gene regulation and modification via RNA interference, alternative splicing and genomic imprinting. NATs are also involved in several human diseases.
We describe a novel method to detect the occurrence of target NATs in specific plant tissues. This method differs from the others currently used in molecular biology laboratories for a number of reasons, particularly the simplicity and versatility of application, low cost and lower material requirement. We demonstrate that NATs can be detected by using diluted cDNA, avoiding the need for a large amount of RNA, thus differing from basic techniques, such as Northern blot hybridisation and reverse-transcription PCR amplification. Furthermore, our method also allows the precise detection of long NATs and their cloning into plasmid vectors for downstream applications. We also reported the first case of a tissue-specific NAT occurring in Oleaceae family and, the antisense orientation of this transcript, allows the splicing of two introns otherwise impossible in the sense orientation.
This method is the first that combines the polymerisation and cleavage activity of DNA polymerase and exonuclease enzymes, respectively, to discover NATs in living organisms. It may simplify the discovery of NATs in plants providing a new strategy for an easy identification and characterization of this group of RNA molecules. Furthermore, since NATs are found in multiple eukaryotes, our method can be easily applied to a wide range of organisms, including human, mice and yeast.
Natural antisense transcripts (NATs) are RNAs with sequences that have complementarity to other endogenous RNAs . NATs are classified into two main categories: cis-NATs, transcribed in cis from opposing DNA strands at the same genomic locus, and trans-NATs, transcribed in trans at separate loci . The cis-NATs, which are transcribed from the reverse-complementary strand of an annotated gene and, hence, are fully or partially overlapping with their respective mRNAs, are also known as antisense RNAs (asRNAs) . Several methods were developed to predict asRNAs in different species, and it was estimated that asRNAs occur for a proportion of annotated genes equal to 22-26% of the human genome [3–5], 15-29% of the mouse genome [5–7], and approximately 9% of the Arabidopsis thaliana genome [8, 9]. Databases that include thousands of predicted NATs for both the animal and plant kingdoms were also developed and are now available on-line [5, 10].
In recent years, NATs were found to be implicated in a wide range of aspects of gene expression in eukaryotes, including genomic imprinting, RNA interference, translational regulation, alternative splicing, X-inactivation and RNA editing [1, 11–15]. It is also known that NATs may be involved in several human diseases and in several responses to stress in plants [14, 16].
In the last decade, the massive application of NGS technologies allowed the implementation of the large repertoires of genomic and transcribed sequence data for the in silico prediction of NATs [5, 17–23]. Nevertheless, an in vitro validation of their presence in specific cell populations and the timing of expression are crucial to determine the role played by a given NAT. To date, the detection of NATs within tissues or organs is achieved by reverse transcription PCR-mediated amplification (RT-PCR) and Northern blot hybridisation techniques, which are based on the use of a large amount of total RNA and require RNase-free reagents and disposables to avoid degradation of the RNA.
Moreover, Northern blot hybridisations are time consuming and require labelled probes, and certain RNAs produce secondary structures during the reverse transcription step of RT-PCR that can act as primers for retro-transcriptase, hence determining primer-independent cDNA synthesis. It has been shown that retro-transcriptases that also have RNase H activity [specially designed for reverse transcription with different amount of RNA and any additional RNase H digestion step, such as Omniscript® Reverse Transcription (Qiagen)] can reduce the amount of primer-independent cDNA synthesis , but the cost of this enzyme remains relatively high, and the method needs to be optimised for each case. Furthermore, the presence of RNase H activity during first-strand synthesis may increase the degradation of the template mRNA, resulting in decreased full-length and first-strand cDNA [see technical note reported in protocols of the most common reverse-transcriptase enzymes, such as SuperScript® III (Invitrogen)]. Our method overcomes these problems as it is based on the use of a cDNA template. However, it is also possible that single-strand cDNA might produce secondary structures and, thus, prime DNA-dependent DNA-polymerase, producing primer-independent DNA synthesis. Using tagged primers that contain a known sequence at the 5’ end and specific primers designed on the known sequence (TAG) during the next amplification is sufficient to overcome this technical issue .
Here, we report a simple method to detect NATs within a population of cDNAs that avoids the direct use of RNA, as it is based on diluted cDNA templates (up to 1:25) retro-transcribed from small amounts of total RNA (up to 100 ng). Theoretically, if a pair of primers is able to amplify a gene in a diluted cDNA sample, then the method can be efficiently applied for the discovery of NATs. It is worth noting that, during RNA retro-transcription using gene-specific primers (GSPs), some GSPs fail to prime cDNA synthesis, even though they work well in PCR using DNA templates (see technical note reported in protocols of the most common reverse-transcriptase enzymes). Our method is more reliable because it is based on the same cDNA template as PCR.
Our method was developed and tested using an anther-specific NAT found in olive (Olea europaea L.) that is related to SLG of Brassica, OeSLG, a gene encoding a protein that functions as an enhancer of the self-incompatible response in some Brassica lines, even though the presence of the protein is not crucial for the rejection of self-pollen. Our overall results demonstrated the usefulness of the method to reveal natural antisense transcripts, and its distinctive features are critically discussed.
Detection of tissue specific transcripts
The amplicons were cloned and sequenced. Intriguingly, three forms of OeSLG were isolated from anther tissue: one sense OeSLG, as in the other tissues, and two forms differing from the sense OeSLG by the lack of one or two internal regions. In detail, the three transcripts consisted of one long form (1,269 bp), equivalent to that expected for OeSLG, one intermediate form (1,181 bp), lacking an 88 nt region, and one short form (780 bp), which was characterised by the lack of the same 88 nt region and an additional 401 nt region (Additional file 1). When we compared the antisense sequences of the intermediate and the short forms to that expected from the OeSLG coding sequence, we found that both missing regions were flanked by the splice donor (GU) and the splice acceptor (AG) sequence sites. We speculated that these regions were introns excised from the antisense sequences under the control of an anther-specific promoter. To demonstrate our hypothesis, we tested our new method to reveal the putative NATs.
Detection of NAT using different tissues
The method was applied to three different tissues, pistils, anthers and leaves, using one 5’-tagged forward primer designed to anneal to the start codon of the OeSLG gene (OeSLG_full-length_5’-tag_FOR) to detect the sense transcripts and one 5’-tagged reverse primer designed to anneal to the stop codon of the OeSLG gene (OeSLG_full-length_5’-tag_REV) to detect the antisense transcripts. One 5’-tagged forward primer (OeEF_5’-tag_FOR) was also added to each reaction to detect the sense transcripts of the OeEF gene used as housekeeping gene. The final amplifications of the OeSLG gene were performed using one primer designed to anneal to the TAG (primer_TAG) and GSPs designed to anneal to the stop codon (OeSLG_full-length_REV) and start codon (OeSLG_full-length_FOR) to amplify the full-length sense and antisense transcripts, respectively. OeEF was amplified using the TAG primer and one internal reverse GSP (OeEF_REV). The results showed that the sense transcripts of OeSLG were present in the pistils, anthers and leaves, whereas the antisense transcripts of OeSLG were present only in the anthers (Figure 2, panel B).
As expected, two different antisense forms were detected in the anthers because the two different splicing variants were found to be specifically expressed in this organ: the intermediate (1,181 bp) and short (780 bp) forms. To validate the first result, we also tested internal primers for OeSLG and OeEF. As in the above experiments, the sense transcripts of the OeSLG gene were detected using the forward 5’-tagged primer (OeSLG_full-length_5’-tag_FOR) and its antisense transcripts using the reverse 5’-tagged primer (OeSLG_full-length_5’-tag_REV); the sense transcript of the OeEF gene was detected using the forward 5’-tagged primer OeEF_5’-tag_FOR. A further reaction was performed without adding primers to test the synthesis of the complementary strands in a primer-independent manner. The final amplifications were performed using the TAG primer in combination with the internal reverse GSP OeSLG_internal_5’_REV to detect the sense transcripts, the internal forward GSP OeSLG_internal_3’_FOR to detect the antisense transcripts, and the reverse GSP OeEF_REV (Figure 2, panel C). The results showed the expression of the sense OeSLG transcripts in the pistils, anthers and leaves and also confirmed the specificity of the antisense transcripts only in the anthers. The reactions performed without adding primers showed no signal, meaning that the amplicons obtained were reliable.
Our method proved to be suitable to demonstrate the presence of antisense transcripts specifically expressed in anthers, and this orientation allowed the splicing of introns otherwise not possible in the sense orientation due to the lack of the splice donor and acceptor sites.
Interaction between sense and antisense OeSLG transcripts in anthers
The results showed the expected bands, meaning that interaction between the two forms, the entire sense and spliced antisense OeSLG, is possible.
RT-PCR analysis with gene-specific primers
Our method provided several advantages compared to the methods normally used to detect NATs, such as Northern blot hybridisation and RT-PCR amplification. First, the present method avoids the use of large amounts of RNA and all of the problems associated with the use of RNA, such as the risk of RNA degradation and requirement of RNase-free conditions. It is also occasionally difficult to collect a large amount of RNA from some tissues or organs, therefore working with diluted cDNA allows for larger number of experiments using the same amount of RNA. Moreover, some fresh tissues are not available throughout the year, such as woody plants, which cannot be grown in greenhouses and that flower only once per year. Thus, the first-strand cDNA synthesised from limited amount of RNA can be diluted and used as a template, assuring that experiments can be conducted without waiting for the next flowering period.
It is known that during RNA retro-transcription using GSPs, a common method used to detect NATs in specific tissues or organs, some GSPs fail to prime cDNA synthesis, even though they work well in PCR using DNA templates (technical note reported in protocols of the most common reverse-transcriptase enzymes). Because our method is based on cDNA, if a given pair of primers is able to amplify a gene using a DNA template, then it is more reliable to use this pair to prime the synthesis of the complementary strand of the target cDNA compared to the retro-transcription of RNA.
Our method is less expensive in comparison to the commercial kits or reagents required for RT-PCR amplification and for Northern blot hybridisation and also reduces the risks associated with the use of toxic reagents involved in the detection procedures.
During recent years, RNA-dependent DNA polymerase enzymes were constantly modified, and some of them now can also function at a high temperature, such as 50°C, increasing the stringency and the specificity of those reactions using GSPs . Yet, after a short step of annealing the primer, our method is performed at an even higher temperature, at 72°C, allowing more specificity when compared to RT-PCR.
To the best of our knowledge, this method is the first that combines the more amenable condition of using DNA instead of RNA with the polymerisation and cleavage activity of DNA polymerase and exonuclease enzymes, respectively, allowing for the precise detection of NATs in living organisms and providing several advantages compared to the conventional methods currently used. We also characterized the first NATs in Oleaceae family and we demonstrated its expression restricted to the anthers. Further studies should be aimed at elucidating the biological function of these NATs in olive anthers. Furthermore, since NATs are found in multiple eukaryotes, our method can be easily applied to a wide range of organisms, including human, mice and yeast.
Plant material, RNA extraction and cDNA synthesis
Total RNA was obtained from the pistils, anthers, leaves, branches and roots of the olive cultivar Leccino using the RNeasy Plant Mini Kit (Qiagen) and was treated by DNase (Qiagen) according to the manufacturer’s protocols. The RNA samples were analysed by agarose gel electrophoresis to test the integrity and were quantified using a nano-drop spectrophotometer.
The method consists of four main steps: first-strand cDNA synthesis by RNA retro-transcription, synthesis of the complementary strand of target cDNAs, degradation of the ssDNAs in solution, amplification by PCR of those dsDNAs that were not degraded (Figure 1).
First-strand cDNAs were produced from 250 ng of total RNA by RT-PCR and SuperScript III (Invitrogen) using Oligo(dT) primers. The samples were then treated with RNaseH (Invitrogen) according to the manufacturer’s protocol.
The synthesis of the complementary strand of the target cDNA was performed using 5 μl of dilute cDNA (1:25) as the template in a final volume of 25 μl of a normal PCR reaction (2.5 μl 10 X Reaction Buffer, 1.2 50 mM MgCl2, 1 μl 10 mM dNTPs, 1 U Taq and 13.1 μl MilliQ water). Then, 1 μl of 10 μM 5’-tagged primer based on the strand to be amplified was added: a forward primer to detect the sense transcript and a reverse primer to detect the putative antisense transcripts, as the cDNA is a ssDNA orientated in the 3’→5’ direction of the mRNA (5’→3’). Along with the gene-specific primer, 1 μl of 10 μM 5’-tagged forward primer annealing to Elongation Factor α1, a housekeeping gene, was added to each reaction for a positive control. An initial denaturation step at 94°C for 5 minutes, subsequent annealing at 60°C (the melting temperature of the primer) for 30 seconds and an elongation phase at 72°C for 30 minutes allowed a conventional DNA-dependent DNA polymerase to synthesise the complementary strand of the target cDNA. At the end of this step, those cDNAs having homology with the primer will be dsDNA, whereas all of the other cDNAs will remain as ssDNA.
The degradation of the ssDNA was performed by a treatment with exonuclease, an enzyme that recognises ssDNAs, releasing deoxyribonucleoside 5’-monophosphates in a stepwise manner. Exonuclease I (2 U, Fermentas) was added to each reaction and incubated at 37°C for 1 hour, followed by 15 minutes at 80°C to inactivate the enzyme. This step degraded all of the ssDNAs; in other words, those cDNAs that are not preserved by the complementary strand formed during the previous step will be degraded.
Name and sequence (5’ – 3’) of primers used for the experiments
One reaction was also performed without primers to test for possible cDNA secondary structures, triggering the primer-independent synthesis of complementary strands. This reaction was also used to test whether the method can be applied to detect possible interactions between the sense and antisense transcripts.
Test of interaction between sense and antisense OeSLG
To test the interaction between the sense and antisense OeSLG, the full-length sequences were cloned and amplified using primers designed to include the start and stop codons. The products were extracted with phenol-chloroform-isoamyl alcohol and precipitated with ethanol, separated using agarose gel electrophoresis to evaluate the integrity and quantified using a nano-drop spectrophotometer. To employ a stoichiometric balance between the antisense (780 bp) and the sense (1,269 bp) transcripts, we used a 1:1.62 ratio of the two transcripts. Three different reactions were performed. The first included 1,620 ng of sense OeSLG, the second consisted of 810 ng and 500 ng of sense and antisense OeSLG, respectively, and the third consisted of 1,000 ng of antisense OeSLG. Annealing buffer (10 mM Tris [pH 8,0], 50 mM NaCl, and 1 mM EDTA [pH 8,0]) was added to the DNA to a final volume of 50 μl, and the reactions were heated at 94°C for five minutes and cooled at room temperature for three hours. An aliquot was assayed by 1% agarose gel electrophoresis.
RT-PCR analysis using gene-specific primers
In order to get a further validation of our method, we performed RT-PCRs using tagged gene-specific primers (GSPs) in RNAs isolated from pistils and anthers. In particular, we used a reverse primer (i.e. OeSLG_full-length_5’-tag_REV) designed to specifically detect the sense transcripts and a forward primer (i.e. OeSLG_full-length_5’-tag_FOR) designed to specifically detect the antisense transcripts. Furthermore we added to each reaction OeEF_REV to detect the Elengation Factor α1 used as housekeeping gene. The subsequent amplifications were carried out using primer_TAG in combination with OeSLG_internal_FOR to detect sense transcripts and in combination with OeSLG_internal_REV to detect antisense transcripts. OeEF was amplified by using OeEF_FOR and OeEF_REV (Table 1).
This research was supported by Progetto Strategico MiPAAF “OLEA - Genomica e Miglioramento genetico dell’olivo”, D.M. 27011/7643/10. Thanks are due to Fiammetta Alagna and Luciana Baldoni, Institute of Plant Genetics – CNR (Italy) to have kindly supplied olive plant materials.
- Lavorgna G, Dahary D, Lehner B, Sorek R, Sanderson C, Casari G: In search of antisense. Trends Biochem Sci. 2004, 29 (2): 88-94. 10.1016/j.tibs.2003.12.002.View ArticlePubMedGoogle Scholar
- Georg J, Honsel A, Voss B, Rennenberg H, Hess W: A long antisense RNA in plant chloroplasts. New Phytol. 2010, 186 (3): 615-622. 10.1111/j.1469-8137.2010.03203.x.View ArticlePubMedGoogle Scholar
- Yelin R, Dahary D, Sorek R, Levanon E, Goldstein O, Shoshan A, Diber A, Biton S, Tamir Y, Khosravi R, Nemzer S, Pinner E, Walach S, Bernstein J, Savitsky K, Rotman G: Widespread occurrence of antisense transcription in the human genome. Nat Biotechnol. 2003, 21 (4): 379-386. 10.1038/nbt808.View ArticlePubMedGoogle Scholar
- Chen JJ, Sun M, Kent W, Huang X, Xie H, Wang W, Zhou G, Shi R, Rowley J: Over 20% of human transcripts might form sense-antisense pairs. Nucleic Acids Res. 2004, 32 (16): 4812-4820. 10.1093/nar/gkh818.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Y, Liu XS, Liu QR, Wei L: Genome-wide in silico identification and analysis of cis natural antisense transcripts (cis-NATs) in ten species. Nucleic Acids Res. 2006, 34 (12): 3465-3475. 10.1093/nar/gkl473.PubMed CentralView ArticlePubMedGoogle Scholar
- Okazaki Y, Furuno M, Kasukawa T, Adachi J, Bono H, Kondo S, Nikaido I, Osato N, Saito R, Suzuki H, Yamanaka I, Kiyosawa H, Yagi K, Tomaru Y, Hasegawa Y, Nogami A, Schonbach C, Gojobori T, Baldarelli R, Hill D, Bult C, Hume D, Quackenbush J, Schriml L, Kanapin A, Matsuda H, Batalov S, Beisel K, Blake J, Bradt D: Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs. Nature. 2002, 420 (6915): 563-573. 10.1038/nature01266.View ArticlePubMedGoogle Scholar
- Kiyosawa H, Yamanaka I, Osato N, Kondo S, Hayashizaki Y: Antisense transcripts with FANTOM2 clone set and their implications for gene regulation. Genome Res. 2003, 13 (6): 1324-1334. 10.1101/gr.982903.PubMed CentralView ArticlePubMedGoogle Scholar
- Jen CH, Michalopoulos I, Westhead D, Meyer P: Natural antisense transcripts with coding capacity in Arabidopsis may have a regulatory role that is not linked to double-stranded RNA degradation. Genome Biol. 2005, 6 (6): R51-10.1186/gb-2005-6-6-r51.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang XJ, Gaasterland T, Chua N: Genome-wide prediction and identification of cis-natural antisense transcripts in Arabidopsis thaliana. Genome Biol. 2005, 6 (4): R30-10.1186/gb-2005-6-4-r30.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen DJ, Yuan C, Zhang J, Zhang Z, Bai L, Meng Y, Chen L, Chen M: PlantNATsDB: a comprehensive database of plant natural antisense transcripts. Nucleic Acids Res. 2012, 40: D1187-D1193. 10.1093/nar/gkr823.PubMed CentralView ArticlePubMedGoogle Scholar
- Kiyosawa H, Abe K: Speculations on the role of natural antisense transcripts in mammalian X chromosome evolution. Cytogenet Genome Res. 2002, 99 (1–4): 151-156.View ArticlePubMedGoogle Scholar
- Werner A, Carlile M, Swan D: What do natural antisense transcripts regulate?. RNA Biol. 2009, 6 (1): 43-48. 10.4161/rna.6.1.7568.View ArticlePubMedGoogle Scholar
- Tufarelli C, Stanley J, Garrick D, Sharpe J, Ayyub H, Wood W, Higgs D: Transcription of antisense RNA leading to gene silencing and methylation as a novel cause of human genetic disease. Nat Genet. 2003, 34 (2): 157-165. 10.1038/ng1157.View ArticlePubMedGoogle Scholar
- Yu WQ, Gius D, Onyango P, Muldoon-Jacobs K, Karp J, Feinberg A, Cui H: Epigenetic silencing of tumour suppressor gene p15 by its antisense RNA. Nature. 2008, 451 (7175): 202-206. 10.1038/nature06468.PubMed CentralView ArticlePubMedGoogle Scholar
- Ohhata T, Hoki Y, Sasaki H, Sado T: Crucial role of antisense transcription across the Xist promoter in Tsix-mediated Xist chromatin modification. Development. 2008, 135 (2): 227-235.View ArticlePubMedGoogle Scholar
- Mercer TR, Dinger M, Mattick J: Long non-coding RNAs: insights into functions. Nat Rev Genet. 2009, 10 (3): 155-159. 10.1038/nrg2521.View ArticlePubMedGoogle Scholar
- Katayama S, Tomaru Y, Kasukawa T, Waki K, Nakanishi M, Nakamura M, Nishida H, Yap CC, Suzuki M, Kawai J, Suzuki H, Carninci P, Hayashizaki Y, Wells C, Frith M, Ravasi T, Pang KC, Hallinan J, Mattick J, Hume DA, Lipovich L, Batalov S, Engström PG, Mizuno Y, Faghihi MA, Sandelin A, Chalk AM, Mottagui-Tabar S, Liang Z, Lenhard B: Antisense Transcription in the Mammalian Transcriptome. Science. 2005, 309 (5740): 1564-1566.View ArticlePubMedGoogle Scholar
- Cheng J, Kapranov P, Drenkow J, Dike S, Brubaker S, Patel S, Long J, Stern D, Tammana H, Helt G, Sementchenko V, Piccolboni A, Bekiranov S, Bailey D, Ganesh M, Ghosh S, Bell I, Gerhard D, Gingeras T: Transcriptional maps of 10 human chromosomes at 5-nucleotide resolution. Science. 2005, 308 (5725): 1149-1154. 10.1126/science.1108625.View ArticlePubMedGoogle Scholar
- Røsok Ø, Sioud M: Systematic identification of sense-antisense transcripts in mammalian cells. Nat Biotechnol. 2004, 22 (1): 104-108. 10.1038/nbt925.View ArticlePubMedGoogle Scholar
- Ge XJ, Wu Q, Jung Y, Chen J, Wang S: A large quantity of novel human antisense transcripts detected by LongSAGE. Bioinformatics. 2006, 22 (20): 2475-2479. 10.1093/bioinformatics/btl429.View ArticlePubMedGoogle Scholar
- Ge XJ, Rubinstein W, Jung Y, Wu Q: Genome-wide analysis of antisense transcription with Affymetrix exon array. BMC Genomics. 2008, 9: 27-10.1186/1471-2164-9-27.PubMed CentralView ArticlePubMedGoogle Scholar
- He YP, Vogelstein B, Velculescu V, Papadopoulos N, Kinzler K: The Antisense Transcriptomes of Human Cells. Science. 2008, 322 (5909): 1855-1857. 10.1126/science.1163853.PubMed CentralView ArticlePubMedGoogle Scholar
- Core LJ, Waterfall J, Lis J: Nascent RNA Sequencing Reveals Widespread Pausing and Divergent Initiation at Human Promoters. Science. 2008, 322 (5909): 1845-1848. 10.1126/science.1162228.PubMed CentralView ArticlePubMedGoogle Scholar
- Haddad F: Potential pitfalls in the accuracy of analysis of natural sense-antisense RNA pairs by reverse transcription-PCR. BMC Biotechnol. 2007, 7: 21-10.1186/1472-6750-7-21.PubMed CentralView ArticlePubMedGoogle Scholar
- Purcell MK, Hart S, Kurath G, Winton J: Strand- specific, real-time RT-PCR assays for quantification of genomic and positive-sense RNAs of the fish rhabdovirus, Infectious hematopoietic necrosis virus. J Virol Methods. 2006, 132 (1–2): 18-24.View ArticlePubMedGoogle Scholar
- Collani S, Alagna F, Colao C, Galla G, Ramina A, Baldoni L, Perrotta G, Muleo R, Barcaccia G: Self-incompatibility in olive: a new hypothesis on the S-locus genes controlling pollen-pistil interaction. Acta Horticulturae. 2012, in pressGoogle Scholar
- Faghihi MA, Wahlestedt C: Regulatory roles of natural antisense transcripts. Nat Rev Mol Cell Biol. 2009, 10 (9): 637-643. 10.1038/nrm2738.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.