SSR markers in transcripts of genes linked to post-transcriptional and transcriptional regulatory functions during vegetative and reproductive development of Elaeis guineensis
© Tranbarger et al; licensee BioMed Central Ltd. 2011
Received: 12 July 2011
Accepted: 3 January 2012
Published: 3 January 2012
The oil palm (Elaeis guineensis Jacq.) is a perennial monocotyledonous tropical crop species that is now the world's number one source of edible vegetable oil, and the richest dietary source of provitamin A. While new elite genotypes from traditional breeding programs provide steady yield increases, the long selection cycle (10-12 years) and the large areas required to cultivate oil palm make genetic improvement slow and labor intensive. Molecular breeding programs have the potential to make significant impacts on the rate of genetic improvement but the limited molecular resources, in particular the lack of molecular markers for agronomic traits of interest, restrict the application of molecular breeding schemes for oil palm.
In the current study, 6,103 non-redundant ESTs derived from cDNA libraries of developing vegetative and reproductive tissues were annotated and searched for simple sequence repeats (SSRs). Primer pairs from sequences flanking 289 EST-SSRs were tested to detect polymorphisms in elite breeding parents and their crosses. 230 of these amplified PCR products, 88 of which were polymorphic within the breeding material tested. A detailed analysis and annotation of the EST-SSRs revealed the locations of the polymorphisms within the transcripts, and that the main functional category was related to transcription and post-transcriptional regulation. Indeed, SSR polymorphisms were found in sequences encoding AP2-like, bZIP, zinc finger, MADS-box, and NAC-like transcription factors in addition to other transcriptional regulatory proteins and several RNA interacting proteins.
The identification of new EST-SSRs that detect polymorphisms in elite breeding material provides tools for molecular breeding strategies. The identification of SSRs within transcripts, in particular those that encode proteins involved in transcriptional and post-transcriptional regulation, will allow insight into the functional roles of these proteins by studying the phenotypic traits that cosegregate with these markers. Finally, the oil palm EST-SSRs derived from vegetative and reproductive development will be useful for studies on the evolution of the functional diversity within the palm family.
Oil palm (Elaeis guineensis Jacq.), a perennial monocotyledonous tropical crop species that belongs to the family Arecaceae, is now the world's number one source of edible vegetable oil, and also the richest dietary source of provitamin A. While the worldwide demand for palm oil increases each year, new elite genotypes from traditional breeding programs provide a yield increase of only 1% per year and the long selection cycle (10-12 years) makes genetic improvement slow . Furthermore, to increase overall oil productivity without new expansion of oil palm cultivation in tropical forest regions with high biodiversity, there is a great need to develop molecular markers for molecular assisted breeding programs targeted to facilitate genetic improvement in yield, and as markers of other important agronomic characters of interest.
Microsatellite markers or simple sequence repeats (SSRs) are tandem DNA repeats from 1-6 bp that are found throughout the coding and non-coding regions of eukaryotic genomes. Non-coding SSRs are often highly polymorphic, co-dominant and simple to detect, and therefore easily adapted to use in high-throughput PCR-based genotyping. They have also been developed for a wide number of crop species and used for various important applications such as genome mapping, diversity studies, and QTL analysis . Due to their highly polymorphic nature, non-coding SSRs are especially useful for fingerprinting or varietal identification studies, but have limited use for studies with more distantly related species . While SSRs derived from non-coding genomic DNA are not transcribed, SSRs identified within transcript sequences can be associated to a function and linked more easily to a phenotypic trait of interest, making them useful for functional diversity studies [2, 3]. In addition, ESTs within genic regions are more transferable for use in diversity studies with more distantly related species. Furthermore, the presence of SSRs in transcribed regions can result in changes in function, transcription or translation. Indeed, SSRs in the coding regions that result in amino acid changes can cause either gain or loss of function, while the presence of SSRs in the 5'UTR can affect transcription or translation, and SSRs in the 3'UTR can affect splicing [3–5].
In the case of oil palm, previous studies have reported the identification of putative SSRs within available EST data [6–8]. However, very few EST-SSRs have been tested nor their usefulness been compared with SSRs identified from the non-coding parts of the genome [9, 10]. Indeed, the genetic maps available for oil palm are mainly based on anonymous non-coding SSRs, AFLPs or RAPDs [9, 11–14]. The oil palm EST-SSRs identified thus far that reveal polymorphisms are mainly from genes that lack similarity with known sequences or encode proteins with unknown function [7, 8]. In fact, in these two recently published articles from oil palm, only two EST-SSRs reported were similar to sequences with known functions. Furthermore, the RNAs used to produce the ESTs for those SSR searches were derived from a narrow range of tissue sources, mainly in vitro materials [6–8]. Despite the relatively low number of EST-SSR markers developed for oil palm, the few that have been tested for interspecies transferability indicate great promise for utilization for comparative genomic studies [7, 8]. Therefore, a strategy to identify EST-SSRs is not only important for diversity studies as a basis for molecular breeding strategies with oil palm, but in addition, markers in conserved coding regions allow easy transferability to other species within the Arecaceae family and provide tools for evolutionary and functional diversity studies [3–5].
The current study has the objective to identify SSRs in ESTs derived from oil palm and to evaluate their utility as molecular markers with plant material used in genetic improvement programs. In particular, we focused on the identification of SSRs in ESTs that originate from developing vegetative and reproductive tissues, and examine their potential for use in mapping and molecular breeding, in addition to functional diversity and genomics analyses within the Arecaceae family.
Results and discussion
Characteristics of SSRs derived from oil palm ESTs
Summary of intra-library analysis of ESTs derived from oil palm developmental related cDNA libraries used for the SSR analysis
cDNA Library Source Material
Shoot apex (normal)
Jouannic et al. 2005
Shoot apex (abnormal1)
Jouannic et al. 2005
Early somatic embryos2 (SSH) BAP3/no BAP
Early somatic embryos2 (SSH) no BAP/BAP
Jouannic et al. 2005
Jouannic et al. 2005
Male inflorescences SSH normal/abnormal1
Beule et al. 2011
Male inflorescences SSH abnormal1/normal
Beule et al. 2011
embryogenic suspension cells SSH +2,4-D4
Lin et al. 2009
embryogenic suspension cells SSH -2,4-D5
Lin et al. 2009
Zygotic Embryos un-normalized
Jouannic et al. 2005
Zygotic Embryos SSH-normalized
Total ESTs analyzed
Results compiled from cluster and SSR analysis of the 12 oil palm cDNA libraries
Steps of Analyses
SSRs found within unigene set
EST-SSRs with similarities found by BLASTX at NCBI
EST-SSRs similar to Unknown, predicted or hypothetical protein
EST-SSRs with no similarities in NCBI databases
GO biological process annotations
GO molecular function annotations
GO cellular component annotations
EST-SSRs with annotations and possible primer pairs identified in flanking sequences
EST-SSR primer pairs synthesized and tested for polymorphisms in a cross LM2TxDA10D
EST-SSR primer pairs that amplified PCR product
From earlier studies, a total of 544 genomic SSRs were identified from a total of 378 unisequences, or 243,943 bp [9, 10]. The distribution of the repeat number for SSR dinucleotide motifs found in the genomic SSRs was different between the EST-SSRs and the genomic SSRs (Figure 1c). Indeed, a peak of 6 repeats (the minimum repeat number cut-off parameter used to search for SSRs) was observed for the dinucleotide EST-SSRs, while the genomic dinucleotide SSRs had a distribution peak at 17-18 repeats. A higher quantity of low repeat numbers in the coding versus genomic SSRs may reflect the higher selective pressure of the coding portion compared to the noncoding portion of the genome.
SSRs that detect polymorphisms within LM2T and DA10D crosses
List of 40 EST-SSRs with similarities to known genes
GASA, gibberelic acid stimulated
Rubber elongation factor
fiber protein Fb2
UniGOS12 (GOLGI SNARE 12); SNARE binding
Protein Destination and Storage
26S proteasome regulatory particle non-ATPase sub-unit 12
zinc ion binding protein
20S proteasome alpha subunit E
60S ribosomal protein L24
proliferating cell nuclear proteinP120
rac GTPase activating protein 3
Transcription and Post-transcription
AP2 domain-containing transcription factor
PHD finger family protein
bZIP transcription factor
C3H-related transcription factor
MADS box transcription factor (AGL2/SEPALLATA) subfamily
NAM; No apical meristem (NAM) NAC-like protein
nucleolar phosphoprotein (RNA binding domain)
RNA binding (RRM/RBD/RNP motifs) family protein
RNA binding protein
BCAS2 protein (spliceosome associated protein)
putative splicing factor 3b, subunit 3 (RNA binding)
UniABC transporter family, cholesterol/phospholipid flippase
Unknown or Unclassified Proteins
An examination of the position of the SSRs within the transcripts revealed that the majority (18) was within the open reading frame (ORF) of the transcript, while eleven were within the 5' untranslated region (UTR), eight within the 3' UTR and 3 overlapped between the 5' UTR start codon and the ORF. Therefore, the majority of the SSRs identified affect the amino acid sequence of the gene product and thus may alter the gene function via a frameshift mutation, while the remainder of the SSRs found in the UTR regions could have an effect on transcription, translation or splicing of gene products .
An examination of the distribution of the EST-SSRs found from the vegetative and reproductive libraries used in the present study revealed five EST-SSRs were derived from the apex (A01 and A11), three from the female inflorescence (F01), fourteen from the male inflorescence (M01, M11 and M21), eighteen from the somatic embryo (E11, E21, S31 and S41) libraries, while no SSRs were found in the sequences from the zygotic libraries (Tables 1 and 3). EST-SSR markers that are associated with a given vegetative or reproductive phase may be useful for studies focused on the inter- and intra-specific functional diversity underlying these tissues in the Arecaceae family.
Transcriptional regulation, in particular through the activity of transcription factors, is known to play a central role in the plant growth and development, and during the evolution of plant form [23, 24]. A survey of the SSRs in the genome of rice and Arabidopsis indicated that amongst the most common GO categories were those related to the nucleus, transcription factor activity, nucleotide binding and DNA binding . Furthermore, a study with humans also found an enrichment in variable repeats in transcripts involved in transcriptional regulation and development . In the present study, SSR polymorphisms were detected in 6 transcripts that encode proteins similar to those that interact with RNA (spliceosome or RNA binding proteins), five similar to transcription factors (TF) including AP2-like, bZIP, zinc finger, MADS-box, and NAC-like TFs, and two transcriptional regulatory proteins including a PHD finger family protein and a retinoblastoma-binding protein. Interestingly, transcriptional regulators are not only central to the evolution of plant form, but also are associated with domestication of crop species [23, 26]. In particular, MADS-box genes were frequent targets of selection during maize domestication . The high number of polymorphisms in transcripts that encode proteins involved in transcriptional regulation in oil palm elite breeding material raises the question of a possible relation to the improvement gained from the reciprocal recurrent selection scheme developed for oil palm. However, the relatively low number of ESTs examined must be taken into account and a conclusive analysis awaits the availability of the genomic sequence of oil palm. Future objectives include the examination of the phenotypic consequences of these SSRs in different oil palm genetic material, and within the Arecaceae family as a whole to determine their relevance to the functional diversity observed.
SSRs in transcripts encoding proteins involved in transcriptional regulation and other functions found from the current study provide pertinent markers for applications such as mapping, molecular breeding and QTL analysis, in addition to the potential for uses in functional diversity studies within the oil palm and between other palm species. In particular, the identification of SSRs in transcripts related to transcriptional control will allow studies aimed at understanding the functional role of these genes in relation to the emerging domestication of the oil palm. However, it should be noted that due to the limited proportion of polymorphic SSRs present in the coding regions, it is important to develop the full range of other potential polymorphic markers in order to combine structural and functional genomics studies on a large genome-scale to allow marker-assisted selection in oil palm.
Plant material production
The preparation of embryogenic suspension cells and RNA extractions for the suppression subtractive hybridization (SSH) library constructions from the 30-day proliferation cycle and after 16 days of liquid pretreatment to initiate somatic embryogenesis was performed and described previously . In addition, a portion of pretreated embryogenic suspension cells initiated to undergo somatic embryogenesis was plated on solid agar plates containing the basal medium with or without 6-benzylaminopurine (synthetic cytokinin) for further somatic embryo development and collected after 7 days for RNA extractions and SSH library constructions. The material collected for the shoot apex, female and male inflorescences, and zygotic embryos for the unnormalized library constructions was described previously . The material for the normal and abnormal male inflorescences SSH libraries was described and performed previously . The zygotic embryos (3-5.5 months of development) were isolated from tenera palm seeds collected from trees (Deli x La Mé origin) cultivated at CRAPP Pobé Station, Benin.
cDNA library construction
The unnormalized libraries (A0, A1, F0, M0 and Z0) were constructed previously  and the SSH libraries were constructed as described previously . The zygotic embryo normalized SSH library was constructed using cDNA made from RNA extracted by RNAeasy lipid (Qiagen) from zygotic embryos (Table 1). The same cDNA was used for both the driver and tester library normalization.
EST generation, analysis, annotation and data mining to identify SSR markers
The ESTs originating from the SSH cDNA libraries (Table 1, libraries E1, E2, M1, M2, S3, S4 and Z1) were generated using standard high throughput sequencing by GATC Biotech AG, Germany. The DNA templates were subjected to single pass automated sequencing using the ABI3730 (Perkin Elmer, Foster City, CA, USA). The sequences were then subjected to an automated procedure to verify cleanse, store and analyze sequences as previously described . The automated analyses allowed the identification of potential unigenes (contigs plus singletons) through simultaneous cluster analysis. Finally, to assign putative functions to the ESTs, BLASTX http://www.ncbi.nlm.nih.gov/BLAST/ was used to compare sequences with the GenBank non-redundant protein sequence database as previously described . The ESTs were manually assigned to functional categories based on a previous catalogue system . In addition, Gene Ontology (GO, http://geneontology.org/)-based annotation was performed using Blast2GO to assign GO molecular function, biological process and cellular component terms . The sequences were analyzed using BLASTX against a GO-based plant uniprot database with an E-value cutoff of 10-10. To identify SSRs within the oil palm EST collection, the online SSR Analysis Tool (SAT; http://sat.cirad.fr/sat) with the default parameters for the SSRIT program was used . The complete list of Elaeis guineensis EST-SSR loci with their EST GenBank accession numbers, derived primer pairs, melting temperatures and predicted PCR product sizes are included in Additional file 1: Table S1. An annealing temperature of 52°C and an MgCl2 concentration of 0.6 mM was used for the PCR reactions performed as described previously . The ESTs from the libraries E1, E2 and Z1 were submitted to GenBank and were assigned the accession numbers JK668500-JK669437, JK669438-JK669618 and JK668122-JK668499 respectively.
Expressed Sequence Tag
Open Reading Frame
Simple Sequence Repeat
Suppression Subtractive Hybridization
We express our gratitude for the financial support from the French MAE, French Embassy in Thailand, the Thai Ministry of Education, the National Center for Genetic Engineering and Biotechnology (BIOTEC) Thailand who participated within the context of a French/Thailand bilateral PHC THAILANDE 2007-2008 project (code 16589YK) that allowed this study to be performed. Institutional funding was also provided by IRD and CIRAD. We thank Frederique Aberlenc-Bertossi for zygotic and Thierry Beulé for the somatic embryos respectively used to construct the cDNA libraries, the CRAPP Pobé oil palm research station, Benin for the genetic material and logistical support, and Stéphane Jouannic for technical assistance for the suppression subtractive hybridization cDNA library construction. We also thank Christine Tranchant and Virginie Pomiès for their excellent bioinformatics and technical assistance respectively.
- Corley RHV, Tinker PB: The Oil Palm.4edition.Oxford, UK: Blackwell Science;2003.View ArticleGoogle Scholar
- Varshney RK, Graner A, Sorrells ME: Genomics-assisted breeding for crop improvement. Trends Plant Sci. 2005, 10 (12): 621-630. 10.1016/j.tplants.2005.10.004.PubMedView ArticleGoogle Scholar
- Varshney RK, Graner A, Sorrells ME: Genic microsatellite markers in plants: features and applications. Trends Biotechnol. 2005, 23 (1): 48-55. 10.1016/j.tibtech.2004.11.005.PubMedView ArticleGoogle Scholar
- Lawson MJ, Zhang L: Housekeeping and tissue-specific genes differ in simple sequence repeats in the 5'-UTR region. Gene. 2008, 407 (1-2): 54-62. 10.1016/j.gene.2007.09.017.PubMedView ArticleGoogle Scholar
- Li YC, Korol AB, Fahima T, Nevo E: Microsatellites within genes: structure, function, and evolution. Mol Biol Evol. 2004, 21 (6): 991-1007. 10.1093/molbev/msh073.PubMedView ArticleGoogle Scholar
- Low E-TL, Alias H, Boon SH, Shariff EM, Tan CY, Ooi LC, Cheah SC, Raha AR, Wan KL, Singh R: Oil palm (Elaeis guineensis Jacq.) tissue culture ESTs: identifying genes associated with callogenesis and embryogenesis. BMC Plant Biol. 2008, 8: 62-10.1186/1471-2229-8-62.PubMedPubMed CentralView ArticleGoogle Scholar
- Singh R, Zaki NM, Ting N-C, Rosli R, Tan S-G, Low E-TL, Ithnin M, Cheah SC: Exploiting an oil palm EST database for the development of gene-derived SSR markers and their exploitation for assessment of genetic diversity. Biologia. 2008, 63 (2): 227-235. 10.2478/s11756-008-0041-z.View ArticleGoogle Scholar
- Ting N-C, Zaki NM, Rosli R, Low E-TL, Ithnin M, Cheah SC, Tan S-G, Singh R: SSR mining in oil palm EST database: application in oil palm germplasm diversity studies. J Genet. 2010, 89 (2): 135-145. 10.1007/s12041-010-0053-7.PubMedView ArticleGoogle Scholar
- Billotte N, Marseillac N, Risterucci AM, Adon B, Brottier P, Baurens FC, Singh R, Herran A, Asmady H, Billot C, et al: Microsatellite-based high density linkage map in oil palm (Elaeis guineensis Jacq.). Theor Appl Genet. 2005, 110 (4): 754-765. 10.1007/s00122-004-1901-8.PubMedView ArticleGoogle Scholar
- Billotte N, Risterucci AM, Barcelos E, Noyer JL, Amblard P, Baurens FC: Development, characterisation, and across-taxa utility of oil palm (Elaeis guineensis Jacq.) microsatellite markers. Genome. 2001, 44 (3): 413-425.PubMedView ArticleGoogle Scholar
- Cochard B, Adon B, Rekima S, Billotte N, de Desmier Chenon R, Koutou A, Nouy B, Omore A, Purba AR, Glazsmann JC, et al: Geographic and genetic structure of African oil palm diversity suggests new approaches to breeding. Tree Genetics and Genomes. 2009, 5: 493-504. 10.1007/s11295-009-0203-3.View ArticleGoogle Scholar
- Mayes S, Jack PL, Marshall DF, Corley RHV: Construction of a RFLP genetic linkage map for oil palm (Elaeis guineensis Jacq). Genome. 1997, 40 (1): 116-122. 10.1139/g97-016.PubMedView ArticleGoogle Scholar
- Moretzsohn MC, Nunes CDM, Ferreira ME, Grattapaglia D: RAPD linkage mapping of the shell thickness locus in oil palm (Elaeis guineensis Jacq.). Theor Appl Genet. 2000, 100 (1): 63-70. 10.1007/s001220050009.View ArticleGoogle Scholar
- Rance KA, Mayes S, Price Z, Jack PL, Corley RHV: Quantitative trait loci for yield components in oil palm (Elaeis guineensis Jacq.). Theor Appl Genet. 2001, 103 (8): 1302-1310. 10.1007/s122-001-8204-z.View ArticleGoogle Scholar
- Beulé T, Camps C, Debiesse S, Tranchant C, Dussert S, Sabau X, Jaligot E, Alwee SSRS, Tregear JW: Transcriptome analysis reveals differentially expressed genes associated with the mantled homeotic flowering abnormality in oil palm (Elaeis guineensis). Tree Genet Genomes. 2011, 7: 169-182. 10.1007/s11295-010-0323-9.View ArticleGoogle Scholar
- Jouannic S, Argout X, Lechauve F, Fizames C, Borgel A, Morcillo F, Aberlenc-Bertossi F, Duval Y, Tregear J: Analysis of expressed sequence tags from oil palm (Elaeis guineensis). FEBS Lett. 2005, 579 (12): 2709-2714. 10.1016/j.febslet.2005.03.093.PubMedView ArticleGoogle Scholar
- Lin HC, Morcillo F, Dussert S, Tranchant-Dubreuil C, Tregear JW, Tranbarger TJ: Transcriptome analysis during somatic embryogenesis of the tropical monocot Elaeis guineensis: evidence for conserved gene functions in early development. Plant Mol Biol. 2009, 70 (1-2): 173-192. 10.1007/s11103-009-9464-3.PubMedView ArticleGoogle Scholar
- Meunier J, Gascon JP: Le schéma général d'amélioration du palmier à huile à l' I.R.H.O. Oleagineux. 1972, 27: 1-12.Google Scholar
- Lehmann T, Hawley WA, Collins FH: An Evaluation of Evolutionary Constraints on Microsatellite Loci Using Null Alleles. Genetics. 1996, 144: 1155-1163.PubMedPubMed CentralGoogle Scholar
- Gupta PK, Rustgi S, Mir RR: Array-based high-throughput DNA markers for crop improvement. Heredity. 2008, 101 (1): 5-18. 10.1038/hdy.2008.35.PubMedView ArticleGoogle Scholar
- Jones N, Ougham H, Thomas H, Pasakinskiene I: Markers and mapping revisited: finding your gene. New Phytol. 2009, 183 (4): 935-966. 10.1111/j.1469-8137.2009.02933.x.PubMedView ArticleGoogle Scholar
- Bevan M, Bancroft I, Bent E, Love K, Goodman H, Dean C, Bergkamp R, Dirkse W, Van Staveren M, Stiekema W, et al: Analysis of 1.9 Mb of contiguous sequence from chromosome 4 of Arabidopsis thaliana. Nature. 1998, 391 (6666): 485-488. 10.1038/35140.PubMedView ArticleGoogle Scholar
- Doebley J, Lukens L: Transcriptional regulators and the evolution of plant form. Plant Cell. 1998, 10 (7): 1075-1082.PubMedPubMed CentralView ArticleGoogle Scholar
- Riechmann JL, Heard J, Martin G, Reuber L, Jiang C, Keddie J, Adam L, Pineda O, Ratcliffe OJ, Samaha RR, et al: Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000, 290 (5499): 2105-2110. 10.1126/science.290.5499.2105.PubMedView ArticleGoogle Scholar
- Legendre M, Pochet N, Pak T, Verstrepen KJ: Sequence-based estimation of minisatellite and microsatellite repeat variability. Genome Res. 2007, 17 (12): 1787-1796. 10.1101/gr.6554007.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhao Q, Thuillet AC, Uhlmann NK, Weber A, Rafalski JA, Allen SM, Tingey S, Doebley J: The role of regulatory genes during maize domestication: evidence from nucleotide polymorphism and gene expression. Genetics. 2008, 178 (4): 2133-2143. 10.1534/genetics.107.085506.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhao Q, Weber AL, McMullen MD, Guill K, Doebley J: MADS-box genes of maize: frequent targets of selection during domestication. Genet Res (Camb). 2011, 93 (1): 65-75. 10.1017/S0016672310000509.PubMed CentralView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talon M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36 (10): 3420-3435. 10.1093/nar/gkn176.PubMedPubMed CentralView ArticleGoogle Scholar
- Dereeper A, Argout X, Billot C, Rami JF, Ruiz M: SAT, a flexible and optimized Web application for SSR marker development. BMC Bioinformatics. 2007, 8: 465-10.1186/1471-2105-8-465.PubMedPubMed CentralView ArticleGoogle 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, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.