Synteny of Prunus and other model plant species
© Jung et al. 2009
Received: 03 September 2008
Accepted: 10 February 2009
Published: 10 February 2009
Skip to main content
© Jung et al. 2009
Received: 03 September 2008
Accepted: 10 February 2009
Published: 10 February 2009
Fragmentary conservation of synteny has been reported between map-anchored Prunus sequences and Arabidopsis. With the availability of genome sequence for fellow rosid I members Populus and Medicago, we analyzed the synteny between Prunus and the three model genomes. Eight Prunus BAC sequences and map-anchored Prunus sequences were used in the comparison.
We found a well conserved synteny across the Prunus species - peach, plum, and apricot - and Populus using a set of homologous Prunus BACs. Conversely, we could not detect any synteny with Arabidopsis in this region. Other peach BACs also showed extensive synteny with Populus. The syntenic regions detected were up to 477 kb in Populus. Two syntenic regions between Arabidopsis and these BACs were much shorter, around 10 kb. We also found syntenic regions that are conserved between the Prunus BACs and Medicago. The array of synteny corresponded with the proposed whole genome duplication events in Populus and Medicago. Using map-anchored Prunus sequences, we detected many syntenic blocks with several gene pairs between Prunus and Populus or Arabidopsis. We observed a more complex network of synteny between Prunus-Arabidopsis, indicative of multiple genome duplication and subsequence gene loss in Arabidopsis.
Our result shows the striking microsynteny between the Prunus BACs and the genome of Populus and Medicago. In macrosynteny analysis, more distinct Prunus regions were syntenic to Populus than to Arabidopsis.
Prunus belongs to Rosaceae, the third most economically important plant family in the United States and other temperate regions of the world [1, 2]. Within Rosaceae, Prunus contains the most diverse array of crops: fruits such as peach, apricot, plum, sweet cherry, and sour cherry, nuts like almond, and lumber trees like black cherry. Other important fruit producing crops in Rosaceae include apple (Malus), pear (Pyrus), raspberry/blackberry (Rubus) and strawberry (Fragaria). In addition, Rosaceae contains a wide variety of ornamental plants including roses, flowering cherry, crabapple and quince.
Significant conservation of the genomes among the Prunus member species has been shown by comparative mapping studies . For example, comparisons of the anchor marker positions on the Prunus reference map with those on 13 other maps showed that the genomes of the diploid Prunus species are essentially collinear . Large collinear blocks among different genera in Rosaceae, such as Prunus and Malus, were also detected . Similar conserved collinearity of shared markers has also been observed among other closely related species within the grasses , legumes [5, 6] and between potato and tomato in the Solanaceae .
With the increase in available genomic sequence data, more laboratories have looked at collinearity at the genome sequence level. The synteny between rice and other cereals are shown to be high , and substantial micro-collinearity among legumes was reported [9–11]. These comparisons have revealed that the conserved syntenic regions, detected by marker collinearity, were often interrupted by small genome rearrangements such as insertions, deletions, inversions and translocations [12–14]. Large scale sequence data has also enabled the detection of micro-collinear regions in less closely related species without apparent macrosynteny . Small scale microsynteny was observed when a 276 kb region of the model legume Medicago (rosid I) was compared with the genomic sequence of Arabidopsis (rosid II) . Conversely, no macrosynteny was found when the two genomes were compared using 172 mapped markers of Medicago . Significant levels of local synteny were also detected covering segments of 1 Mb of Arabidopsis and regions of <5 cM in lettuce and sunflower, members of the asterids, even though the macrosyntenic patterns covering larger segments were not evident . A recent study detected microsynteny between coffee, another asterid species, and the Arabidopsis genome using putative orthologous sequences and BAC ends . We have previously detected conserved syntenic regions between Prunus and Arabidopsis using 475 peach ESTs anchored to Prunus genetic maps and 1097 peach ESTs anchored to BAC contigs . The conserved syntenic regions were short and fragmentary, and often the Prunus regions matched to more than one Arabidopsis block. This complex network of microsynteny, often including non-collinear regions, between Arabidopsis and the distantly related species suggests multiple genome duplications followed by differential gene loss [15–21]. Large segmental or whole genome duplication followed by gene loss seems to be prevalent in the evolution of all flowering plant genomes, and the degree of genome instability seems to be higher in some species of Brassicaceae, such as Arabidopsis, and in some species of Poaceae as well [22–24].
Synteny analysis across species provides insight on the evolutionary relationships between different lineages and also the opportunity to study the relationship between genome structure and function of organisms. The micro-collinear regions conserved between model organisms like Arabidopsis and lesser-studied organisms can also facilitate marker saturation and candidate gene searches. For example, micro-collinearity data between Arabidopsis and rice has been utilized in the improvement of phylogenetic resolution of the expansin gene family . Multiple rounds of polyploidization or large segmental duplication followed by gene loss, however, can greatly obscure the synteny, impeding the transfer of genomic knowledge from model species genomes to those of less well characterized species.
The recent availability of the whole genome sequence of Populus trichocarpa provides another resource in the detection of conserved synteny among plant genomes . Also available are the partially sequenced Medicago genome . Populus and Medicago belong to rosid I which also includes Prunus; hence Prunus is evolutionarily closer to them than Arabidopsis which belongs to rosid II. Populus is of particular interest to compare with Prunus since Populus is completely sequenced, is a fellow tree species and it appears to be more stable compared to the Arabidopsis genome; nucleotide substitution, tandem gene duplication and gross chromosomal rearrangement proceed more slowly in Populus then in Arabidopsis . The detection of melon BAC regions that are more conserved in Populus than in Medicago or Arabidopsis  and the report of papaya BACs that are more collinear in Populus than in Arabidopsis  also suggest that the Populus genome may be useful in the exploration of distantly related species genomes.
With near completion of the peach physical map  and the development of EST-derived genetic markers, the numbers of sequences that are anchored to Prunus genetic maps or the peach physical map have doubled since our previous synteny analysis between Prunus and Arabidopsis. Also available to us are four homologous Prunus BAC sequences and four additional peach BAC sequences. Using these data, we analyzed the degree and the pattern of the conserved syntenic region between Prunus and Populus or Arabidopsis. We also compared the Prunus BAC sequences with the partially sequenced Medicago genome to gain further insight on the genome evolution of the related plant species.
Sequences from four Prunus BACs, and four additional peach BACs, 028F08, 082I18, pPn31C7, and PpN089G02, were used to assess the degree of conserved synteny with the Populus, Medicago, and Arabidopsis genomes. The four Prunus BACs were selected using the same probe to study disease-resistant genes in Prunus species, so they potentially represent homologous genomic regions. These include one peach BAC with two contigs (058P54-C23 and 058P54-C24), one apricot BAC (AprC27), and two plum BACs (Plum045O02 and Plum080O24). Gene prediction was conducted with FGENESH . The two peach BACs, 028F08 and 082I18, overlap by 2 kb and two genes, predicted by FGENESH program (see Methods). The combined sequences of the two BACs are 77.7 kb long and contain a total of 38 predicted genes (028F08-082I18_1 to 028F08-082I18_38). The peach BACs pPn31C7 and PpN089G02 were 48.8 kb and 132.2 kb long and had nine (pPn31C7_1 to pPn31C7_9) and 42 (PpN089G02_1 to PpN089G02_42) predicted genes, respectively. The two contigs of a peach BAC and the three Prunus BACs that contain putative disease-resistant genes, 058P54-C23, 058P54-C24, AprC27, Plum045O02 and Plum080O24, had 19 (058P54-C23_1 to058P54-C23_19), 9 (058P54-C24_1 to 058P54-C24_9), 11 (AprC27_1 to AprC27_11), 13 (045O02_1 to 045O02_13) and 4 (080O24_1 to 080O24_4) predicted genes, respectively. The sizes of the peach, apricot and plum BAC sequences were 90 kb, 41 kb and 70.1 kb, respectively. Cumulatively, the total length of peach BACs used in our analysis was 350 kb.
The predicted gene sequences in these Prunus BACs and their Populus, Medicago and Arabidopsis homologs were used to find conserved syntenic regions. Syntenic groups were identified when the distance between the two adjacent matches was less than 200 kb and when the syntenic regions contain at least four gene pairs. We first detected collinear syntenic regions using the DAGchainer program  and then merged the overlapping syntenic regions. This method gave us a more comprehensive view of the syntenic regions which include sections where gene contents are conserved but not the gene order due to small-scale genome rearrangements.
Synteny quality within the individual homologus block between the Prunus BACs and Populus or Medicago.
No. genes in BAC
No. genes in Chr
*synteny quality (%)
The E values of most of gene matches in the syntenic regions were considerably lower than our cut-off value, 1E-6 (see Methods). 96.3% and 97.1% of the gene matches in Prunus-Populus and Prunus-Medicago, respectively, had an E value of less than 1E-10. The median E values were 1E-93 and 1E-109 in Prunus-Populus and Prunus-Medicago, respectively.
To assess the degree of macrosynteny conservation, we used two sets of map-anchored Prunus sequences in the analysis. One set comprised of 1093 sequences that are anchored to the TxE Prunus reference map . The majority of these sequences are peach fruit ESTs but also include sequences from markers that were directly used in mapping. Some ESTs are anchored to multiple positions and the number of anchored sequences in each linkage group is as follows: 281 in G1, 364 in G2, 268 in G3, 424 in G4, 221 in G5, 209 in G6, 65 in G7 and 310 in G8. Another set comprised of peach EST sequences that are anchored to the peach physical map. This data comprised of 2140 EST sequences that are anchored to 1500 BAC contigs and their Populus and Arabidopsis homologs.
These map-anchored sequences and their Populus and Arabidopsis homologs were used to find conserved syntenic regions. The syntenic groups were selected when the distance between the two adjacent matches were less than the maximum distance and when the syntenic regions contain at least four gene pairs (See Methods).
Characteristics of syntenic blocks conserved between Prunus map anchored sequences and Arabidopsis or Populus.
# of syntenic regions
# collinear groups
# gene pairs
(Populus or Arabidopsis)
before the merge
after the merge
424 kb - 3.26 Mb
5.1 - 15.7 cM2
(2.8 - 8.6 Mb)
197 kb - 1 Mb
4.6 - 15.7 cM2
(2.5 - 8.6 Mb)
158 kb - 1.1 Mb
(95 - 432.5 kb)
(81.3 - 533.8 kb)
Distribution of syntenic blocks in Prunus TxE map (A) and peach physical map (B).
Synteny with Populus
Synteny with Arabidopsis
# syntenic groups
region in each LG contating syntenic groups (cM)
# syntenic groups
Region in each LG contating syntenic groups (cM)
4 (I, XII, XIV)
25.6 - 42.75
7 (Chr1, 2, 3, 5)
25.6 - 41.3
3 (II, XII, XIV)
14 - 19.7
7 (Chr1, 2, 3, 4, 5)
7.4 - 19.7
26.1 - 35.5
3 (Chr1, 3)
26.1 - 39.1
Synteny with Populus
Synteny with Arabidopsis
# syntenic groups
region in each contig contating syntenic groups
# syntenic groups
region in each contig contating syntenic groups*
2 (Chr2, Chr3)
5 (II, VI, IX, XIII, XIV)
14 (Chr1, Chr2, Chr3, Chr4, Chr5)
2 (II, V)
10 groups matching to 5 contigs
17 groups matching to 3 contigs
The position of 28 major trait loci affecting agronomic characters found in various Prunus species have been previously established in the TxE map using the data from different linkage maps anchored with the TxE reference map . We were interested in finding any syntenic regions conserved between the TxE Prunus map and the Arabidopsis or Populus genomes covering these trait loci. Several syntenic regions matched to trait-loci containing regions in G3 of the Prunus TxE map (Figure 6). The trait loci contained in these syntenic regions are Ag for anther color of peach and almond (yellow/anthocyanic), PcP for the polycarpel trait of peach , and Fc for flower color (pink/pale pink) in peach . All these syntenic regions span 8.3 cM to 15 cM in the TxE map and contained either two or three trait loci (Figure 6). Among the four Populus blocks that matched to these regions, two blocks in LG_XII and LG_XIV seem to represent duplicated blocks, since two genes matched to the same peach sequences and two genes matched to the sequences that are anchored to the same position in the TxE map. The other two Populus blocks in LG_I and LG_XII may also represent duplicated regions since one of the genes matched to the same peach sequence and the rest matched to the sequences that are anchored to the same position in the TxE map. There were seven Arabidopsis blocks that matched to these regions. Some blocks share one to two genes that match to the same peach sequences and other blocks share two to three genes that match to the sequences that are anchored to the same position in the TxE map. The complex network of synteny in Arabidopsis in this region is also indicative of multiple independent genome duplications and rearrangements in the Arabidopsis genome.
Number of syntenic groups between Populus or Arabidopsis genomes and Prunus genetic map or peach physical map that are detected at various significance thresholds.
Total # of groups
We compared the level of conserved synteny between the Prunus genome and the whole genome sequences of two model organisms, Arabidopsis and Populus, and the partially sequenced Medicago genome. Since the whole genome sequence of peach or other Prunus sequences is not currently available, we employed available Prunus BAC sequences to assess the level of conserved microsynteny and map-anchored Prunus sequences to assess the pattern of synteny throughout the genome. In the Populus genome, we found well-conserved syntenic regions in all the BACs studied, and the syntenic regions cover almost the entire length of the BACs in some cases. All the BACs showed conserved synteny with two regions in two different Populus chromosomes, supporting the proposed event of the whole genome duplication in Populus . The conserved gene content in the two different Populus genomes showed that different sets of genes were lost in the duplicated Populus regions. The order of genes in these BACs and both of the Populus regions were very well conserved, without signs of genome translocations. The synteny quality within each block was also considerable. On the contrary, we did not detect that level of microsynteny between Prunus BACs and the Arabidopsis genome. The two microsyntenic regions between these contained only four and five gene matches. We also found well-conserved syntenic regions between Prunus BACs and the partially sequenced Medicago genome. Interestingly, the homologous BACs of plum, peach, and apricot detected three syntenic regions in three different chromosomes of Medicago. These three regions may be indicative of large-scale or segmental duplication events in the evolutionary history of Medicago, in addition to the whole genome duplication event proposed to have occurred within the rosid I clade. The other two peach BACs detected only one syntenic region in Medicago. This may be explained by the fact that the Medicago Genome Release 1.0 constitutes only 38–47% of the entire genome. Another possibility is that these peach BACs do have only one conserved syntenic region in Medicago genome. It has shown that the level of internal synteny within both Medicago and Lotus is much lower than the intergenomic synteny between the two, presumably through significant gene loss and rearrangement after the whole genome duplication event .
Our study analyzed a fine-level microsynteny between Prunus and the three model dicot organisms. Our results substantiate our previous findings that the level of gene-order conservation between peach and Arabidopsis is very fragmentary [17, 36]. The microsynteny between the Prunus and Populus genome, however, is much better conserved, promising the utility of the Populus genome in the study of the Prunus genome and vice versa when the peach genome is sequenced (underway by the Joint Genome Institute). Our results also suggest that the Prunus genome may be closer to Populus than to Medicago, which is evolutionarily closer to Prunus than Populus. Similar results have been reported; two melon BACs showed more conserved synteny to Populus than to the evolutionarily closer Medicago . Lai et al.  reported that ordered papaya BAC end sequences showed a higher level of synteny with poplar than with the more closely related Arabidopsis which both belong to the Brassicales.
We also employed Prunus sequences that are anchored to the Prunus genetic map or the peach physical map to study the level of conserved synteny at the whole genome level. Only the completely sequenced genomes of Populus and Arabidopsis were used in this analysis. A number of syntenic regions that contain several gene pairs were detected between Prunus map anchored sequences and both Populus and Arabidopsis genomes. The number of the syntenic regions was actually higher in Arabidopsis than in Populus. However, in most cases, multiple Arabidopsis regions were syntenic to the overlapping regions in the Prunus TxE map or peach BAC contigs, resulting in a complex network of synteny. In contrast, the syntenic regions in Populus detected more distinct regions in Prunus. We found blocks in both the Prunus genetic map and peach physical map that were syntenic to two different homologous regions in Populus genome. However, we did not find any Prunus regions that detected homologous regions in the Arabidopsis genome despite the fact that multiple Arabidopsis genomic regions were syntenic to overlapping Prunus regions. Our findings are consistent with previous reports that the Arabidopsis genome went through multiple large genome duplication events followed by frequent gene loss [37, 38]. The multiple Arabidopsis genome regions with synteny to overlapping Prunus regions, which we observed, may represent the duplicated regions that subsequently went through gene loss. Simillion et al.  have reported that high frequency of gene loss after duplication in Arabidopsis reduces collinearity and that, at extreme, the duplicated regions no longer share homologous genes. They have shown that they could detect these highly degenerated duplicated blocks by indirect comparison with other segments. The synteny obscured by severe subsequent gene loss can impede the transfer of genomic knowledge from one species to another, suggesting the potential problem of using Arabidopsis in the study of other species such as those in Prunus. In contrast, the Populus genome is reported to be a more stable and homologous genome; blocks that have arisen from the most recent genome-wide duplication event can be clearly detected . Our results corroborate these previous findings and suggest that Populus is better model to study synteny with Prunus genomes.
We report the evidence of well-conserved microsynteny between Prunus BACs and two plant model species, Populus and Medicago. The observed network of synteny supported the whole genome duplication events, proposed to have occurred in Populus and Medicago. Interestingly, the level of synteny conservation seemed higher between Prunus-Populus than in Prunus-Medicago, even though Medicago, rather than Populus, is evolutionarily closer to Prunus. In comparison, the level of conserved synteny between the Prunus BAC and Arabidopsis was insignificant. We also analyzed the extent of conserved synteny between Prunus map-anchored sequences and the completely sequenced genomes of Populus and Arabidopsis. We detected a number of syntenic regions that contain several gene pairs between Prunus map anchored sequences and both Populus and Arabidopsis genomes. In comparison to the synteny between Populus and Prunus, we observed a more complex network of synteny between Arabidopsis and Prunus in which multiple Arabidopsis regions shows synteny to overlapping regions in Prunus.
Three sets of Prunus sequences were used in the analysis of conserved synteny between Prunus and the whole genome sequences of Arabidopsis and Populus: Prunus BAC sequences, Prunus sequences anchored to the Prunus reference map [1, 39], and the peach EST sequences anchored to the peach physical map . Prunus BAC sequences were also compared with the partially sequenced Medicago genome .
The sequences of four peach BACs, 028F08 (AC154900), 082I18 (AC154901), pPn31C7 (AF467900), and PpN089G02 (DQ863257), were downloaded from Genbank. The sizes of the sequences were 66.2 kb, 73.3 kb, 48.8 kb, and 132.3 kb, respectively. We also downloaded the sequences of four Prunus BACs that contain a putative disease resistant gene from the Genome Database for Rosaceae (GDR) . These include one peach BAC with two contigs (058P54-C23 and 058P54-C24), one apricot BAC (AprC27), and two plum BACs (Plum045O02 and Plum080O24). The two peach BACs, 028F08 (AC154900), 082I18 (AC154901) are a contiguous genomic segment. Gene prediction was done using the FGENESH program .
Prunus sequences that are anchored to the TxE Prunus reference map were downloaded from GDR. The majority of the 1093 anchored sequences were derived from peach fruit ESTs . Some ESTs were directly anchored to the map by co-hybridization to BACs with genetic markers, and others were indirectly anchored through hybridization to BACs that belong to the same BAC contig that contain anchored genetic markers. The positions (cM) of the genetic markers were used as the positions for the genetically anchored ESTs on the map. The set also included sequences of markers that are directly used in mapping. The third set, also downloaded from GDR, comprised of 2140 EST sequences that were anchored to 1500 BAC contigs of the peach physical map. The position of the individual BACs in the BAC contigs were used as the positions ESTs in the peach physical map. For the ESTs that are anchored to multiple overlapping BACs in a BAC contig, the innermost left and right positions were assigned.
The sequence data (TAIR6_pep_20060907) and the chromosome coordinate data (sv_gene.data) of the Arabidopsis translated proteins were downloaded from the Arabidopsis Information Resources (TAIR) database . The sequence data (Poptr1_FilteredModels_proteins.fa) and the chromosome coordinate data (Poptr1_FilteredModels.gff) of the Populus translated proteins were downloaded from the JGI web site http://genome.jgi-psf.org/Poptr1/Poptr1.download.ftp.html. Medicago genome release version 1.0 data was downloaded from http://medicago.org. International Medicago Genome Annotation Group (IMGAG) annotated protein sequences (MTpep1.txt) were used for our analysis.
The predicted protein sequences of the Prunus BACs were compared with the Arabidopsis, Populus, and Medicago proteins by pairwise comparison using the BLASTP program. The top ten matches with an E value less than 10-6 were used for further analysis. Syntenic groups were selected when the distance between the two adjacent matches were less than 200 kb and when the syntenic regions contain at least four gene pairs.
Mapped Prunus sequences that are homologous to the Arabidopsis or Populus proteins were determined by pairwise comparison using the BLASTX program. The top ten matches with E value less than 10-6 were used for further analysis. The syntenic groups were selected when the distance between the two adjacent matches were less than the maximum distance and when the syntenic regions contain at least four gene pairs. For the synteny analysis between the TxE genetic map and Populus or Arabidopsis, the maximum distance was set as the approximate genomic distance that covers about 100 genes - 1 Mb for Populus and 0.5 Mb for Arabidopsis. This was approximated by the genome size, 550 Mb for Populus  and 157 Mb for Arabidopsis , and the predicted number of genes, 45,555 for Populus (annotation v1.1) and 32,041 (TAIR7 Genome release) for Arabidopsis. For sequences anchored in the TxE Prunus genetic map, we selected a less strict parameter, 5 cM, as the maximum distance since most of the sequences are not directly mapped, but are indirectly anchored to the position by the association with other mapped markers (see above). 5 cM is 0.95% of the TxE map (524 cM)  and 0.5 Mb that we used for Arabidopsis is about 0.35% of the genome (157 Mb). We used a smaller maximum distance - 0.5 Mb for Populus and 0.25 Mb for Arabidopsis - for the synteny analysis between the peach physical map and Populus or Arabidopsis, since the sequences are positioned to BAC contigs, which are much smaller than chromosomes or linkage groups. The maximum distance in peach BAC contig was set to 200 unit lengths (unit given in the FPC output), which corresponds to 0.25 Mb. One unit length corresponds approximately to 1.25 kb, estimated from the total physical length of the peach physical map, 303 Mb,  and the total unit length of the BAC contigs, 242772.
We employed DAGchainer  to detect collinear segments conserved in Prunus and other genomes of Arabidopsis, Medicago and Populus. When we used map-anchored Prunus sequences, we found many cases in which multiple collinear segments, detected by the software, reside in the same or overlapping genome regions. This occurs when the two genomes in comparison contain syntenic regions in which gene contents are conserved but the gene order is not well conserved, which is caused by differential small genome rearrangements such as conversions, translocations and duplications in the evolutionary history of two genomes. In these regions, DAGchainer can detect more than one collinear syntenic region in the same genomic block by selecting the gene pairs with the conserved order. To solve this problem, we merged the overlapping syntenic regions that are generated by DAGchainer. The syntenic groups were merged when they overlap or when they are separated by less than the maximum distance between the two adjacent matches in both genomes. The merged syntenic groups represent syntenic regions in which the gene contents are conserved but not the gene order. This procedure provides us with a more accurate assessment of the degree of conserved synteny between the two genomes. The Prunus blocks that match to more than one site in the Populus or Arabidopsis genome were detected by selecting syntenic regions that share more than three Prunus sequences. All the in-house developed scripts are available to the public when requested.
To determine whether the syntenic groups we report were detected by chance, we tested the statistical significance for each group. Both of the Populus and Arabidopsis genomes were randomized by leaving the locations the same but permuting the gene names. We analyzed 1000 simulated Populus and Arabidopsis genomes for the occurrence of each conserved syntenic group and calculated the probability of the match occurring by chance.
Genome Database of Rosaceae
The Arabidopsis Information Resources
International Medicago Genome Annotation Group.
The work was funded by National Science Foundation Plant Genome Program (#0320544 to D.M.); United States Department of Agriculture Cooperative State Research, Education and Extension Service-National Research Initiative-Plant Genome Program (#2005-35300-15452 to A.A.). The authors would like to thank Dr Laura Georgi from Clemson University for providing the plum BAC library, Dr Richard McCombie from the Lita Annenberg Hazen Genome Center at Cold Spring Harbor Laboratory for assistance in finishing of peach BACs, 028F08 (AC154900), 082I18 (AC154901) and Dr Chris Saski from the Clemson University Genomics Institute for assistance in finishing the four Prunus BACs; 058P54, AprC27, Plum045O02 and Plum080O24.
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.