Characterization of the transcripts of human cytomegalovirus UL144
© He et al; licensee BioMed Central Ltd. 2011
Received: 12 April 2011
Accepted: 14 June 2011
Published: 14 June 2011
The genome of human cytomegalovirus (HCMV) has been studied extensively, particularly in the UL/b' region. In this study, transcripts of one of the UL/b' genes, UL144, were identified in 3 HCMV isolates obtained from urine samples of congenitally infected infants.
Northern blot hybridization, cDNA library screening, and RACE-PCR were used.
We identified at least 4 differentially regulated 3'-coterminal transcripts of UL144 in infected cells of 1,300, 1,600, 1,700, and 3,500 nucleotides (nt). The 1600 nt transcript was the major form of UL144 mRNA. The largest transcript initiated from the region within the UL141 open reading frame (ORF) and included UL141, UL142, UL143, UL144, and UL145 ORFs.
These findings reveal the complex nature of the transcription of the UL144 gene in clinical isolates.
Keywordscytomegalovirus transcription UL144
Cytomegaloviruses (CMVs) are members of the subfamily Betaherpesvirinae of the family Herpesviridae. Human cytomegalovirus (HCMV) is a common cause of congenital viral infections and a frequent opportunistic pathogen in transplant recipients and AIDS patients [1–3]. The relative genomic complexity of HCMV is mirrored by its biological characteristics, as the most relevant cellular reservoirs of the latent virus and its sites for permissive replication have not been conclusively established and its pathogenesis is not well understood. Like other CMVs, HCMV has a very specific host range, but within a permissive host it enters and replicates in a wide variety of cell types [3, 4].
The HCMV genome consists of 230 to 235 kb of double stranded DNA and more than 160 predicted open reading frames (ORFs) [5–8]. The overall nucleotide sequence of strains isolated from unrelated sources is relatively conserved. However, the sequences of specific ORFs can be highly variable, including RL6, RL12, UL4, UL18, UL55 (gB), UL73 (gN), UL74 (gO), UL139, UL144, and UL146 [5, 9–18].
The HCMV UL144 ORF encodes a homologue of the herpesvirus entry mediator (HVEM) , which is a member of the tumor necrosis factor receptor (TNFR) superfamily and contributes to the ability of HCMV to escape immune clearance. Although the UL144 gene has been well characterized at the DNA and protein levels, no detailed analysis of its mRNA transcripts has been reported. In the current study, we characterized the transcripts of the UL144 gene in 3 clinical strains.
Cells and virus
Human embryonic lung fibroblast (HELF) cells were grown in Dulbecco's modified Eagle's medium supplemented with 10% (vol/vol) fetal calf serum. Three HCMV clinical isolates, named H, Xu, and Ch, were categorized as Group C, A, and B, according to the schema of Lurain et al. The isolates were derived from urine samples of three congenitally HCMV-infected children in the Pediatrics Department at the Affiliated Shengjing Hospital of China Medical University. All samples were collected with the permission of the infants' parents and the research had obtained the approval of the Hospital Ethical Committee. HELF cells were infected with the virus at a multiplicity of infection of approximately 10 plaque-forming units (PFU) per cell. Virus DNA of HCMV isolate H was extracted as previously described  and its UL/b' region was sequenced using a previously described shotgun sequencing method .
Total RNA was extracted from uninfected and HCMV-infected HELF cells using a standard guanidium isothyocyanate and phenol:chloroform method . The extracts were treated with DNase I (Ambion, USA). Immediate early (IE) RNA was prepared in the presence of 100 μ g/mL cycloheximide (Sigma, USA) at 24 hours post-infection (hpi), early (E) RNA was extracted in the presence of 100 μ g/mL phosphonoacetic acid (Sigma) at 48 hpi, late RNA was extracted in the absence of blocking agents at 96 hpi, and control RNA was extracted from uninfected cells .
Screening of the cDNA library
A cDNA library of the HCMV clinical strain H had been previously constructed using the pBluescript II SK vector . The primary library consisted of 1.12 × 106 recombinant clones/mL. Nearly 95% of the vectors contained inserts, as determined by blue-white plaque screening on NZY agar plates containing 5-bromo-4-chloro-3-indolyl-β-D-galactopyranoside (X-Gal) and isopropyl-β-D-thio-galactoside (IPTG). The average length of the inserts was 1.2 kb as determined by PCR with M13 primers and the Gel Image System.
Primers used for PCR amplification, Northern-blot and selection of cDNA clones.
Name of Primers
Position in UL144 ORF
Position in GQ981646
A complete list of all primers is provided in Table 1.
Northern blot hybridization
RNA was prepared from uninfected cells and cells infected with HCMV clinical isolates H, Ch, and Xu at IE, E, and L phases. RNA was analyzed by Northern blot to identify the transcripts of UL144.
Five micrograms of RNA was size-fractionated in a 1.2% denaturing agarose gel and transferred to a membrane (N-Hybond, Amersham Biosciences) by capillary blotting. A UL144 gene-specific probe was labeled with digoxigenin-dUTP, hybridized to the membrane, and detected by its NBT/BCIP chromogenic signal according to the manufacturer's instructions (DIG Northern starter kit, Roche Diagnostics, Mannheim, Germany). RNA representing the noncoding strand was generated using T7 RNA polymerase and labeled by incorporation of digoxigenin-dUTP. Three probes were produced according to the different sequences of isolates H, Ch, and Xu using the primers 144For2 and 144Prob (Table 1). The labeled probes were antisense to the 531 bp of the complete UL144 protein coding sequences (nucleotides 11584-12114, GQ981646 as reference) of the three isolates. Hybridization to the RNA blots was carried out at 55°C overnight. The blots were visualized by the addition of an anti-digoxigenin antibody labeled with alkaline phosphatase followed by the addition of the substrate of 5-bromo-4-chloro-3-indolyl phosphate (BCIP) and nitroblue tetrazolium (NBT). DIG-labeled RNA Molecular Weight Marker I 0.31-6.95 kb (Roche Diagnostics) was used to estimate the size of the bands by the logarithmic relationship between molecular weight and distance migrated.
5' RACE To define the 5' ends of transcripts in the UL144 intergenic region, RNA preparations of the isolates H, Ch, and Xu were analyzed by 5' RACE, which was carried out with the 5'-Full RACE kit (Takara, Japan) according to the manufacturer's instructions. Briefly, the RNA preparation was treated with calf intestine alkaline phosphatase (CIAP, 4.8 U/μg RNA) for 1 h at 50°C. Then the purified RNA was treated with tobacco acid pyrophosphorylase (TAP, 0.25 U/μg RNA) for 1 h at 37°C. The CIAP/TAP-treated RNA was ligated to the supplied 5' RACE adaptor (1 μL of 15 μM adaptor/2 μg RNA) with bacteriophage T4 RNA ligase (20 U/μg RNA) for 1 h at 16°C. Reverse transcription was carried out at 42°C for 1 h with 5 U M-MLV reverse transcriptase and 5 μM of random 9-mer primers. Control reactions were carried out in parallel using RNA treated without M-MLV [M-MLV(-)]or TAP [TAP(-)]. All reactions were performed in the presence of 10 U RNase Inhibitor. cDNA was amplified by nested PCR using 1.25 U LA Taq DNA polymerase (Takara) in a 50 μL reaction mixture including LA PCR buffer (50 mM KCl, 10 mM Tris-HCl, pH 8.3), 1.5 mM MgCl2, and 2.5 mM of each of the four deoxynucleoside triphosphates. The primers used in the primary PCR reaction were the 5' RACE outer primer and 144Rev1 (see Table 1 for primer coordinates within the H isolate sequence) and those in the second PCR reaction were the 5' RACE inner primer and 144Rev2. The PCR conditions included an initial denaturation step at 94°C for 3 min, 20-25 cycles of 94°C for 30 sec, 55°C for 30 sec, and 72°C for 3 min, followed by a final extension step at 72°C for 10 min.
3' RACE 3' RACE was carried out with the 3'-Full RACE Core Set Ver. 2.0 (Takara) following the manufacturer's instructions. Briefly, 500 ng of RNA was used as a template to synthesize cDNA with AMV Reverse transcriptase for 30 min at 50°C, 5 min at 99°C and 5 min at 5°C. The external reverse primer oligo dT-adaptor primer with a poly (T) tract was used to prime cDNA synthesis. cDNA of the 3 clinical isolates, including isolates H, Ch, and Xu, was then amplified by nested PCR with an adaptor primer and the external forward primer 144For2 and the inner forward primer 144For3.
RT-PCR RT-PCR was carried out using the RNA PCR Kit (AMV) Ver. 3.0 (Takara) according to the manufacturer's instructions. Reverse transcription was performed at 42°C for 1 h with 200 U of AMV Reverse Transcriptase (Takara) and 25 mM of oligo dT-adaptor primer. PCR was conducted using 10 μL of the synthesized cDNA and primers 144For1 and 144Rev2 (Table 1) with Ex Taq DNA polymerase. The PCR conditions for all the amplification reactions included an initial incubation at 94°C for 3 min, followed by 40 cycles of denaturation at 94°C for 40 sec, annealing at 50°C for 40 sec, and extension at 72°C for 2 min, with a final incubation at 72°C for10 min. Control reactions using RNA treated without AMV [AMV (-)] and water as templates were included in each set of PCR amplifications.
Cloning and sequencing
RACE and RT-PCR products were separated by electrophoresis on a 2% agarose gel and visualized with ethidium bromide. Bands were gel purified and the recovered fragments were ligated into a TA cloning vector, pCR (Invitrogen, USA). Inserts were sequenced by the dideoxynucleotide method with Sequenase (Invitrogen) and M13For and M13Rev primers using an Applied Biosystems automated sequencer.
Identification of UL144 transcripts by cDNA library screening
Twenty-two distinct clones were identified from the cDNA library and all were sequenced successfully. The 3' ends of the sequences occurred within a 3-nucleotide (nt) window extending from nt 12,844-12,846 (GQ981646 as reference). Five groups of UL144 transcripts were derived from the 22 cDNA sequences according to their different 5' ends. Among the 22 cDNA sequences, 15 (named as C1) transcripts were 1565-1612 nt in length with 5' ends at nt 11,276-11,278 (GQ981646 as reference) and different lengths of poly A tails, 4 (named as C3) were 1317 nt with their 5' ends at nt 11,528, and 1 (named as C2) was 1775 nt with its 5' end at nt 11,091. The other two transcripts were 1230 nt and 1161 nt in length (named as C4 and C5, respectively) and started at nt 11,635 and 11,704, respectively, indicating that both of the transcripts initiate within the UL144 ORF.
Mapping the 5' end of UL144 transcript using 5' RACE-PCR
Mapping the 3' end of UL144 transcript using 3' RACE-PCR
As shown in Figure 1B, 3' RACE-PCR yielded two products of 1.3 kb and 500 bp in length. Five UL144-specific clones with the 1.3 kb band from isolates H (1), Ch (1), and Xu (3) were sequenced. All the sequences terminated at nt 12,845, except for one from isolate Xu, which terminated at nt 12,844. Sequencing of the 500 bp band revealed it to be a nonspecific product.
Identification of UL144 transcripts by Northern blot
The UL144 gene is located in the UL/b' region of the HCMV genome, which contains approximately 20 kb of DNA encoding at least 19 ORFs but is dispensable for growth in vitro. The UL144 ORF encodes a structural homologue of the herpesvirus entry mediator (HVEM) [9, 26], a member of the TNFR superfamily, but it does not bind any known TNF-related ligands [9, 26, 27]. It negatively regulates Th1 cells by binding to the B- and T-cell lymphocyte activator BTLA , mimicking the inhibitory co-signaling function of HVEM. Like LMP-1 of the Epstein-Barr virus, UL144 activates TNFR-activated factor 6 (TRAF6), leading to NF-κB (nuclear factor kappa-light-chain-enhancer of activated B cells) activation in early infection . It induces expression of the macrophage-derived chemokine (MDC) [27, 29, 30], also termed CCL22 [chemokine (C-C motif) ligand 22], which attracts Th2 and regulatory T cells. Taken together, it is possible that upregulation of CCL22 by the UL144 protein may help HCMV evade immune surveillance, resulting in increased viral spread and dissemination. There is significant strain-specific sequence variability in the UL144 ORF [9, 18, 26, 31, 32] and these UL144 genotypes are associated with congenital infection outcome [31, 33]. However, no data has been reported thus far on the UL144 mRNA structure.
Transcripts of the HCMV UL144 gene demonstrated in this study.
The position of 5' end of the transcript
Isolates in which transcripts were identified
Amounts. of cDNA clones from the library
Strength of reaction in Northern blot
Position and forms of the TATA-box
Position and forms of the CAT-box
In this study, all 22 UL144-specific cDNA clones terminated within a 3 nt window located 100 bp downstream of the UL145 gene. The single poly(A) signal, which is downstream of the UL145 stop codon, supports the possibility of a common 3' terminal end for these transcripts. The 3' terminal ends of the transcripts were confirmed by 3' RACE using RNA from all the three isolates. The presence of 3' coterminal transcripts around the UL144 region is similar to that found in other genomic regions of HCMV, including the UL93-UL99 ORFs , UL146-UL132 ORFs , and M73 and M73.5 genes . It is possible that their different 5' non-coding sequences may have regulatory functions, warranting further investigation.
We obtained some small RNA fragments in our cDNA library and by 5' RACE that initiates inside the UL144 ORF. These RNA fragments are most likely products from incomplete transcription. Transcripts of approximately 3,000 nt in the RNA of isolates H and Ch and approximately 2,000 nt in the RNA of isolate Xu were present in our Northern blots but were not confirmed by cDNA library screening or 5' RACE-PCR, indicating that they may be due to non-specific hybridization.
In this study, we have demonstrated that the transcripts of UL144 have a complex transcriptional pattern of overlapping transcripts that differ at their 5' initiation sites in different UL144 genotype isolates. At least 4 transcripts were identified in the 3 clinical isolates studied. These findings reveal the complex nature of the transcription of the UL144 gene in clinical isolates.
List of abbreviations
herpesvirus entry mediator
open reading frame
tumor necrosis factor receptor
Human embryonic lung fibroblast
plaque forming units
TNFR-activated factor 6
nuclear factor kappa-light-chain-enhancer of activated B cells
chemokine (C-C motif) ligand 22.
This study was supported by the National Natural Science Fund of China (30672248, 30770109, 30700916, 30801254 and 30901625).
- Alford C, Stagno S, Pass R, Britt W: Congenital and perinatal cytomegalovirus infections. Rev Infect Dis 1990,12(Suppl 7):S745-753.View ArticlePubMed
- Pass RF: Congenital cytomegalovirus infection and hearing loss. Herpes 2005, 12: 50-55.PubMed
- Ho M: Cytomegalovirus biology and infection. 2nd edition. NEW YORK: plenum publishing corp; 1991.View Article
- Sinzger C, Grefte A, Plachter B, Gouw A, The T, Jahn G: Fibroblasts, epithelial cells, endothelial cells and smooth muscle cells are major targets of human cytomegalovirus infection in lung and gastrointestinal tissues. J Gen Virol 1995,76(Pt 4):741-750.View ArticlePubMed
- Dolan A, Cunningham C, Hector R, Hassan-Walker A, Lee L, Addison C, Dargan D, McGeoch D, Gatherer D, Emery V, Griffiths P, Sinzger C, McSharry B, Wilkinson G, Davison A: Genetic content of wild-type human cytomegalovirus. J Gen Virol 2004, 85: 1301-1312. 10.1099/vir.0.79888-0View ArticlePubMed
- Murphy E, Yu D, Grimwood J, Schmutz J, Dickson M, Jarvis M, Hahn G, Nelson J, Myers R, Shenk T: Coding potential of laboratory and clinical strains of human cytomegalovirus. Proc Natl Acad Sci USA 2003, 100: 14976-14981. 10.1073/pnas.2136652100PubMed CentralView ArticlePubMed
- Chee M, Bankier A, Beck S, Bohni R, Brown C, Cerny R, Horsnell T, Hutchison Cr, Kouzarides T, Martignetti J: Analysis of the protein-coding content of the sequence of human cytomegalovirus strain AD169. Curr Top Microbiol Immunol 1990, 154: 125-169.PubMed
- Davison AJ, Dolan A, Akter P, Addison C, Dargan DJ, Alcendor DJ, McGeoch DJ, Hayward GS: The human cytomegalovirus genome revisited: comparison with the chimpanzee cytomegalovirus genome. J Gen Virol 2003, 84: 17-28. 10.1099/vir.0.18606-0View ArticlePubMed
- Lurain NS, Kapell KS, Huang DD, Short JA, Paintsil J, Winkfield E, Benedict CA, Ware CF, Bremer JW: Human cytomegalovirus UL144 open reading frame: sequence hypervariability in low-passage clinical isolates. J Virol 1999, 73: 10040-10050.PubMed CentralPubMed
- Pignatelli S, Dal Monte P, Landini M: gpUL73 (gN) genomic variants of human cytomegalovirus isolates are clustered into four distinct genotypes. J Gen Virol 2001, 82: 2777-2784.View ArticlePubMed
- Prichard M, Penfold M, Duke G, Spaete R, Kemble G: A review of genetic differences between limited and extensively passaged human cytomegalovirus strains. Rev Med Virol 2001, 11: 191-200. 10.1002/rmv.315View ArticlePubMed
- Hassan-Walker A, Okwuadi S, Lee L, Griffiths P, Emery V: Sequence variability of the alpha-chemokine UL146 from clinical strains of human cytomegalovirus. J Med Virol 2004, 74: 573-579. 10.1002/jmv.20210View ArticlePubMed
- Rasmussen L, Geissler A, Winters M: Inter- and intragenic variations complicate the molecular epidemiology of human cytomegalovirus. J Infect Dis 2003, 187: 809-819. 10.1086/367900View ArticlePubMed
- Valés-Gómez M, Shiroishi M, Maenaka K, Reyburn H: Genetic variability of the major histocompatibility complex class I homologue encoded by human cytomegalovirus leads to differential binding to the inhibitory receptor ILT2. J Virol 2005, 79: 2251-2260. 10.1128/JVI.79.4.2251-2260.2005PubMed CentralView ArticlePubMed
- Bar M, Shannon-Lowe C, Geballe A: Differentiation of human cytomegalovirus genotypes in immunocompromised patients on the basis of UL4 gene polymorphisms. J Infect Dis 2001, 183: 218-225. 10.1086/317939View ArticlePubMed
- Davison A, Akter P, Cunningham C, Dolan A, Addison C, Dargan D, Hassan-Walker A, Emery V, Griffiths P, Wilkinson G: Homology between the human cytomegalovirus RL11 gene family and human adenovirus E3 genes. J Gen Virol 2003, 84: 657-663. 10.1099/vir.0.18856-0View ArticlePubMed
- Chou S, Dennison K: Analysis of interstrain variation in cytomegalovirus glycoprotein B sequences encoding neutralization-related epitopes. J Infect Dis 1991, 163: 1229-1234. 10.1093/infdis/163.6.1229View ArticlePubMed
- Bale JF, Petheram SJ, Robertson M, Murph JR, Demmler G: Human cytomegalovirus a sequence and UL144 variability in strains from infected children. J Med Virol 2001, 65: 90-96. 10.1002/jmv.2006View ArticlePubMed
- Montgomery R, Warner M, Lum B, Spear P: Herpes simplex virus-1 entry into cells mediated by a novel member of the TNF/NGF receptor family. Cell 1996, 87: 427-436. 10.1016/S0092-8674(00)81363-XView ArticlePubMed
- Alderete J, Jarrahian S, Geballe A: Translational effects of mutations and polymorphisms in a repressive upstream open reading frame of the human cytomegalovirus UL4 gene. J Virol 1999, 73: 8330-8337.PubMed CentralPubMed
- Zabarovsky E, Kashuba V, Pettersson B, Petrov N, Zakharyev V, Gizatullin R, Lebedeva T, Bannikov V, Pokrovskaya E, Zabarovska V: Shot-gun sequencing strategy for long-range genome mapping: a pilot study. Genomics 1994, 21: 495-500. 10.1006/geno.1994.1307View ArticlePubMed
- Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem 1987, 162: 156-159.View ArticlePubMed
- Scott G, Barrell B, Oram J, Rawlinson W: Characterisation of transcripts from the human cytomegalovirus genes TRL7, UL20a, UL36, UL65, UL94, US3 and US34. Virus Genes 2002, 24: 39-48. 10.1023/A:1014033920070View ArticlePubMed
- Ma YP, Ruan Q, Ji YH, Wang N, Li ML, Qi Y, He R, Sun ZR, Ren GW: Novel transcripts of human cytomegalovirus clinical strain found by cDNA library screening. Genetics and Molecular Research 2011, 10: 566-575. 10.4238/vol10-2gmr1059View ArticlePubMed
- Cha T, Tom E, Kemble G, Duke G, Mocarski E, Spaete R: Human cytomegalovirus clinical isolates carry at least 19 genes not found in laboratory strains. J Virol 1996, 70: 78-83.PubMed CentralPubMed
- Benedict CA, Butrovich KD, Lurain NS, Corbeil J, Rooney I, Schneider P, Tschopp J, Ware CF: Cutting edge: a novel viral TNF receptor superfamily member in virulent strains of human cytomegalovirus. J Immunol 1999, 162: 6967-6970.PubMed
- Poole E, King C, Sinclair J, Alcami A: The UL144 gene product of human cytomegalovirus activates NFkappaB via a TRAF6-dependent mechanism. EMBO J 2006, 25: 4390-4399. 10.1038/sj.emboj.7601287PubMed CentralView ArticlePubMed
- Cheung T, Humphreys I, Potter K, Norris P, Shumway H, Tran B, Patterson G, Jean-Jacques R, Yoon M, Spear P, Murphy K, Lurain N, Benedict C, Ware C: Evolutionarily divergent herpesviruses modulate T cell activation by targeting the herpesvirus entry mediator cosignaling pathway. Proc Natl Acad Sci USA 2005, 102: 13218-13223. 10.1073/pnas.0506172102PubMed CentralView ArticlePubMed
- Nakayama T, Hieshima K, Nagakubo D, Sato E, Nakayama M, Kawa K, Yoshie O: Selective induction of Th2-attracting chemokines CCL17 and CCL22 in human B cells by latent membrane protein 1 of Epstein-Barr virus. J Virol 2004, 78: 1665-1674. 10.1128/JVI.78.4.1665-1674.2004PubMed CentralView ArticlePubMed
- Schultheiss U, Püschner S, Kremmer E, Mak T, Engelmann H, Hammerschmidt W, Kieser A: TRAF6 is a critical mediator of signal transduction by the viral oncogene latent membrane protein 1. EMBO J 2001, 20: 5678-5691. 10.1093/emboj/20.20.5678PubMed CentralView ArticlePubMed
- Arav-Boger R, Battaglia C, Lazzarotto T, Gabrielli L, Zong J, Hayward G, Diener-West M, Landini M: Cytomegalovirus (CMV)-encoded UL144 (truncated tumor necrosis factor receptor) and outcome of congenital CMV infection. J Infect Dis 2006, 194: 464-473. 10.1086/505427View ArticlePubMed
- Picone O, Costa J, Chaix M, Ville Y, Rouzioux C, Leruez-Ville M: Human cytomegalovirus UL144 gene polymorphisms in congenital infections. J Clin Microbiol 2005, 43: 25-29. 10.1128/JCM.43.1.25-29.2005PubMed CentralView ArticlePubMed
- Arav-Boger R, Willoughby R, Pass R, Zong J, Jang W, Alcendor D, Hayward G: Polymorphisms of the cytomegalovirus (CMV)-encoded tumor necrosis factor-alpha and beta-chemokine receptors in congenital CMV disease. J Infect Dis 2002, 186: 1057-1064. 10.1086/344238View ArticlePubMed
- Welch AR, McNally LM, Gibson W: Cytomegalovirus assembly protein nested gene family: four 3'-coterminal transcripts encode four in-frame, overlapping proteins. J Virol 1991, 65: 4091-4100.PubMed CentralPubMed
- Lurain N, Fox A, Lichy H, Bhorade S, Ware C, Huang D, Kwan S, Garrity E, Chou S: Analysis of the human cytomegalovirus genomic region from UL146 through UL147A reveals sequence hypervariability, genotypic stability, and overlapping transcripts. Virol J 2006, 3: 4. 10.1186/1743-422X-3-4PubMed CentralView ArticlePubMed
- Scalzo AA, Forbes CA, Smith LM, Loh LC: Transcriptional analysis of human cytomegalovirus and rat cytomegalovirus homologues of the M73/M73.5 spliced gene family. Archives of Virology 2009, 154: 65-75. 10.1007/s00705-008-0274-8View ArticlePubMed
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.