| Human Molecular Genetics | Pages |
Large multi-chromosomal duplications encompass many members of the olfactory receptor gene family in the human genome
Introduction
Results
Many OR-containing clones hybridize to multiple genomic locations
Characterization of the extent and structure of the OR-containing duplications
Large-scale heteromorphism of these OR clusters is minimal in humans
Differences among primates in the genomic organization of OR duplications
Discussion
Large duplications in the human genome
Homology among >15 locations comprises both OR and non-OR sequences
The molecular structures of paralogous regions
Plasticity in the cytogenetic arrangement of OR duplications in primates
The structure and distribution of OR-containing paralogs suggest a hierarchical model for the duplication process
A subset of OR-clusters is significantly diverged from other OR-containing regions
OR-locations missing from this collection of genomic clones
Fraction of the human genome devoted to the OR family
Materials and Methods
Isolation of OR-containing genomic clones
FISH
PCR typing of monochromosomal somatic cell hybrid panel
Acknowledgements
References
Large multi-chromosomal duplications encompass many members of the olfactory receptor gene family in the human genome
INTRODUCTION
Humans can discriminate thousands of odors (1). This capability is due to the expression of a diverse repertoire of odorant receptors (OR) in the specialized sensory neurons in the olfactory neuroepithelium. The OR receptors are encoded in the genome by a large family of genes, whose coding regions are short (~1 kb) (2,3). ORs are members of the much larger family of G-protein-coupled receptors with seven transmembrane segments. Genes that are very closely related to ORs are expressed in the tongue and, surprisingly, in testes (4,5). A large subset of OR genes, ranging in similarity from 45 to 100%, can be amplified from the human genome using degenerate PCR primers in two highly conserved transmembrane regions (6). Remarkably, over two-thirds of the human OR sequences captured in this way are apparent pseudogenes (6).
Members of the OR gene family are distributed among many locations in the human genome (5-8). When a pool of OR-specific sequences was used as a probe for fluorescence in situ hybridization (FISH), signals were observed at >25 locations situated on all but a few chromosomes (6). PCR analyses of flow-sorted chromosomes confirmed the dispersed nature of the human OR family (6). Each OR location is likely to comprise multiple genes and/or pseudogenes (5,8,9; V. Brand-Arpon, S. Rouquier, H. Massa, P. de Jong, C. Ferraz, P.A. Ioannou, J.G. Demaille, B.J. Trask and D. Giorgi, submitted for publication).
Only a single OR is expressed in each neuron (9-13). Therefore, the multi-chromosomal distribution of OR genes presents a conundrum for understanding the transcriptional regulation of the OR gene family. The distributed nature of the human OR gene family is in sharp contrast to the organization of antigen-receptor gene families, where diversity is generated by joining elements from a single genomic cluster of possible components. For the OR family, transcriptional-control mechanisms must contend with the multiplicity of gene locations to ensure that only one OR is expressed in each neuron (and from a single allele) (10-13) and that each receptor is expressed in the appropriate zone within the neuroepithelium (9,11,12,14), while insuring that a diverse repertoire of receptors is expressed in the tissue as a whole. It is anticipated that a study of the genomic sequences surrounding transcribed OR genes will reveal how this set of demands is fulfilled at the molecular level.
The history of events leading to the multitude of OR-containing locations in the human genome and the plethora of pseudogenes is also insufficiently understood. The FISH results with OR-specific sequences (6) lead to two hypotheses about the evolutionary relationships among the many OR-containing locations. One possibility is that processed pseudogenes have been inserted into some sites by retrotransposition. In this case, homology would be limited to the transcribed portions of OR sequences. Alternatively, sites may be related as a consequence of interchromosomal duplications of large genomic segments. In this case, the homology is expected to extend to sequences that flank and/or lie between the OR genes. Given the many OR pseudogenes in the human genome, it is relevant to ask whether blocks of pseudogenes have been duplicated. Although several clusters of mouse OR genes are associated with large genomic duplications (15), and a subset of OR genes is repeated near multiple human telomeres within a much larger unit of duplication (7), the evolutionary relationships among most of the OR-containing sites in the genome are not yet known. The observation of extreme polymorphism among humans in the distribution of a subtelomeric block of OR genes (7) also raised the possibility that phenotypic variation could arise through variation in gene copy number or genomic context. It is therefore important to assess the plasticity and variability of the portion of the genome devoted to the OR family.
As a step towards understanding the evolution and transcriptional control of this gene family, we describe here the identification and cytogenetic characterization of genomic clones encompassing OR genes and pseudogenes. We show that the majority of OR clusters contains sequences that are shared by >15 sites in the genome. Only a few OR-containing regions, e.g. at 1q21-22, 1q44, 17p13 and 19p13.2, are sufficiently diverged from the others that clones derived from them hybridize to unique locations in the genome. We show that duplications of large genomic regions have accompanied the spread of OR genes to many locations. Finally, we demonstrate that the chromosomal distribution of OR duplicons has changed during the divergence of the great apes.
RESULTS
Many OR-containing clones hybridize to multiple genomic locations
Genomic clones containing OR-like sequences were identified from two clone libraries, a BAC library of ~80-100-kb inserts cloned from total human genomic DNA and a cosmid library constructed using DNA from flow-sorted chromosome 3. The initial screen was performed by hybridization with a complex pool of OR-specific sequences. This OR pool was PCR amplified from total human genomic DNA using degenerate primers recognizing evolutionarily conserved regions of the OR proteins (see Materials and Methods). Positives were confirmed by (i) PCR amplification of a product of the expected size using the degenerate OR primers, (ii) sequencing of cloned PCR products and/or (iii) Southern blot hybridization of EcoRI-digested DNA using the OR sequence pool as a probe (data not shown). Of the clearly positive clones, 52 BACs and 22 chromosome 3 cosmids were characterized cytogenetically by FISH.
Over half of the 74 clones produced FISH signals at more than one genomic location. Figure
Figure 1. A metaphase spread showing the hybridization signals produced by the OR-containing BAC 51D11. By analyzing 10 such metaphases, it was possible to discriminate recurrent sites of specific hybridization from sporadic background signals. Twenty sites of specific hybridization were detected (Fig. 2A) on the 13 indicated chromosomes. More than one region of hybridization is evident on several chromosomes (2, 3, 7, 11, 12 and 13). Note the marked difference in signal intensity among the locations, which is reflected in the average signal scores that are plotted in Figure 2A.
A

B

Figure 2. Summary of the hybridization signals produced by 74 OR-containing clones identified in a total genomic BAC library and a chromosome 3 specific cosmid library. (A) 32 clones produced signals consistently at 17 sites in the genome, and their signals are summarized in this panel. The black dots summarize the results of 14 OR-BACs, and the red dots summarize the results of 18 OR-containing chromosome 3 cosmids. Each row of dots corresponds to the signals observed by a particular clone. The size of each dot indicates the average signal intensity observed at a given location with a particular clone. The names of the clones are indicated in the two columns to the right. A small T is drawn at the 24 locations where signals were observed with more than four of the clones. Polymorphism (see text) may partially account for ambiguity in the assignment of clones in 2p. The stacks are positioned over the center of the hybridizing regions; the mapping precision is approximately half of a 400-level band. (B) Summary of the hybridization sites of BACs or cosmids that hybridized to three or fewer locations in the genome. Twenty-nine BACs hybridized to single sites; these sites are marked with white circles. Each circle denotes a separate clone, i.e. four different clones mapped to 1q21-22. Four chromosome 3 specific cosmids mapped to single sites, which are marked with black circles. Eight clones mapped to two sites, indicated with pairs of matching triangles. One clone mapped to the three sites indicated with a green square. The relative sizes of the triangles and square reflect the relative hybridization efficiency/intensity at the two or three sites, respectively. The names of the clones mapping to multiple locations are indicated in the figure. The single-locus BACs, listed in order of location, are: 1q21-22: 267C3, 460D10, 821D9, 980D1; 1q43-44: 176F8, 850H7, 992G12; 3q11.2-13.1: 748C6; 4q13-21: 984H6; 5q23-31: 995D3; 5q34: 303F10; 9q32-proximal q34: 17E12, 378E10, 855A10, 963F3, 987D11, 966G7, 996G8; 11p15: 626C11, 978C7; 13q21: 969B7 [*, 13q21 is the predominant site of 969B7, but this clone also produced dim signals infrequently at some of the same locations as the multi-site clones in (A)]; 14q32: 858F6; 17p12-13: 45F12, 284E5, 588A4; 19p13.3: 32E9; 19p13.1-13.2: 3F6, 272A4, 378D8. The single-locus chromosome 3 specific cosmids are 3p14-21: 13, 54; 3q13.3-21: 28; 3q28-29: 45. (The LL03NC01 designations are given in Materials and Methods.)
The multi-chromosomal distribution of BAC 51D11 is typical of 14 (26%) BACs and 18 (82%) chromosome 3 cosmids. Figure
The remaining 38 OR-BACs and four OR-containing chromosome 3 cosmids exhibited more restricted chromosomal distributions than the set typified by 51D11. The FISH results of these clones are summarized in Figure
Eight BAC clones produced signals at two locations. Two patterns are worth drawing attention to because each was observed with two BACs. An OR-related segment, represented by clones 365G5 and 148D6 (grey and dark purple triangles in Fig.
Finally, four (18%) of the OR chromosome 3 cosmids and 28 (54%) of the OR-BACs behaved as simple single-copy probes (black and white circles, respectively, in Fig.
Characterization of the extent and structure of the OR-containing duplications
We wished to determine if the duplications that resulted in homology among the 17+ sites, detected by the clones shown in Figures
(i) The intensity of the FISH signals at most of the cross-hybridizing locations was significantly greater with the genomic clones than with short OR-specific probes.
(ii) The addition of excess unlabeled OR-sequences to biotinylated cosmid 8 failed to attenuate the resulting FISH signals significantly at any of the locations (data not shown).
(iii) PCR assays designed from the sequences at the ends of several cosmids confirm the multiplicity of their cytogenetic locations. These primers were designed to avoid common interspersed repeats and OR homology. Each of the 24 human chromosomes, isolated in monochromosomal rodent somatic cell hybrid lines, was subjected to PCR amplification with these primer-pairs. In contrast to FISH, PCR assays chromosomes for a small region of a clone and demands that the sequence be duplicated with sufficient similarity that PCR-amplification occurs. As expected, each PCR assay yielded a product on chromosome 3, the origin of the cosmids. Primers designed from the T7 end of clone 32 [which produced a FISH signal on chromosomes 2-4, 7-14 and 21 (Fig.
Figure 3. PCR analyses of non-OR portions of clones that FISH-map to multiple locations. (A) Primer pair 32MF/32IR designed from the T7 end of cosmid 32, whose FISH distribution is summarized in Figure 2A. (B) Primer pair U45021 and L45210, designed from sequence 5[prime] of the cluster of OR sequences in the 3p13 contig (PCR assay D in Figure 5). (C) Primer pair U66952 and L67132, designed from the sequence lying between two OR sequences in the 3p13 contig (PCR assay E in Figure 5). The FISH results of cosmids 3, 48 and 81, which overlap one or both of the latter two assays, are summarized in Figures 2A and 4. The lanes are numbered to indicate the chromosomes contained in each hybrid cell line. M, mouse genomic DNA; C, Chinese hamster genomic DNA; H, human genomic DNA; H2O, water control. (iv) FISH analyses of clones derived from a 250 kb contig from 3p13 encompassing several OR genes (V. Brand-Arpon et al., submitted for publication) allowed us to characterize the extent and coarse structure of the homologous duplications. The FISH results of eight clones spanning this contig are summarized in Figure Figure 4. FISH results of eight clones derived from 3p13 and distributed across a 250 kb contig (V. Brand-Arpon et al., submitted for publication). The position of each clone is indicated in Figure 5A. The FISH results of each clone are summarized with a different color. Each stack of dots indicates the relative signal intensity observed in each of 20 chromosomes in 10 metaphases. The diameter of the dots is proportional to the signal intensity, which was scored on a scale of 0-4. A stack of 20 dots signifies that signals were observed on all analyzed chromosomes at that site; shorter stacks indicate hybridization efficiencies <100%. Each group of stacks is aligned with the midpoint of the hybridization range. Signals could typically be assigned with a precision of half of a 400-level band. Note that the black and white stacks for cosmids 96 and 97 are confined to two locations on chromosome 3, and PAC 169 and cosmid 88 (red and green stacks, respectively) show a different pattern than do cosmids 3, 26 and 48. (v) PCR assays confirm the multi-chromosomal distribution of sequences across the 3p13 contig (Fig. (vi) Sequence comparisons demonstrate that ~30 kb of the 3p13 sequence (DDBJ/EMBL/GenBank accession no. AF042089; V. Brand-Arpon et al., submitted for publication) is duplicated within a PAC clone (PAC pDJ392a17) that has been assigned to chromosome 11 (G.A. Evans et al., unpublished; GenBank accession no. AC000385) (Fig.
Large-scale heteromorphism of these OR clusters is minimal in humans
Because we had previously observed striking polymorphism in copy number and chromosomal location of one cluster of OR genes that map exclusively near telomeres (7), we analyzed the FISH patterns of 17 of the OR clones on two or more individuals. This set includes 10 OR-BACs that hybridized to a single location and five OR-BACs and three OR-cosmids that hybridized to multiple sites. With one exception, the patterns of hybridization were very similar on the different individuals. Slight differences could be ascribed to differences in hybridization efficiency and statistical sampling. We observed only a subtle polymorphism in the hybridization pattern on the p-arm of chromosome 2. Several of the multi-site clones produced signals at two locations within the p11.2-p13 interval on chromosomes of two individuals, while only one or the other site was labeled in other individuals. (For simplicity, signals observed at both locations are lumped into a single group in the clone summary in Fig.
Differences among primates in the genomic organization of OR duplications
In order to judge the plasticity in the cytogenetic organization of OR-containing duplications, we compared the FISH patterns of two of the multi-site OR clones in human, chimpanzee, gorilla, orangutan, gibbon and baboon. The results with cosmid 3 in the first four species are summarized in Figure
Figure 5. A simplified map of a 3p13 contig encompassing four OR sequences (V. Brand-Arpon et al., submitted for publication). The locations of the eight clones spanning the contig that were analyzed by FISH are indicated by the horizontal lines. In addition, we also FISH-mapped three clones that were identified by screening the chromosome 3 cosmid library with the T7 end of cosmid 88 (cosmids 111, 113 and 117) and four clones identified by screening the library with the T7 end of PAC 169 (cosmids 202, 204, 209, 215 and 219). These clones may derive from paralogous regions elsewhere on chromosome 3, and their overlap has not yet been established. These clones and cosmids 88, 96 and 97 lack detectable OR sequences. The thick line designates the 106 kb region that has been sequenced (V. Brand-Arpon et al., submitted for publication). The locations of the four OR sequences and an MLCK pseudogene are indicated by the small arrows. The OR sequence in PAC 169 was identified by PCR amplification of the PAC with the OR3B/OR5B primers. Its orientation and exact position is not known, but it appears to lie in the region not overlapped by cosmids 202-219, 88 or 111-117, because the OR3B/5B primers do not amplify OR sequences from these clones. (The LL03NC01 designations for these clones are given in Materials and Methods.) (B) PCR results using PCR-primer pairs designed from the 106 kb sequenced contig of 3p13 (assays C-F) and from the end-sequences of clones in the contig (A, B and G) (see Materials and Methods for details). The positions of the primer pairs are indicated in (A). All chromosomes in the Coriell monochromosomal panel were assayed; only the positive chromosomes are listed in the table. Gels showing assays D and E are shown in Figure 3B and C, respectively. Figure 6. A dot-matrix comparison of homologous regions in 106 kb sequence from 3p13 (GenBank accession no. AF042089) and the reverse-complement of the sequence of PAC 392a17 (G.A. Evans et al., unpublished data; GenBank accession no. AC000385), which has been assigned to chromosome 11. The matrix was generated using the ABI Inherit program, using the following parameter settings: window size, 20; offset, 10; match, 85%. Diagonal lines in the matrix indicate regions of homology between the two sequences. A 20-24 kb duplication unit composed of several distinct homology units can be discerned and is denoted by the colored segments above and to the right of the matrix. The 12 conspicuous blocks of sequence matches correspond to alignments of the three VNTR-like regions of the 106 kb sequence (at approximate nucleotide positions 48330-52320, 56950-60315 and 70530-72620 in AF042089) with the four VNTR regions in the PAC (designated by the grey blocks). These VNTRs are composed of imperfect repeats of a unit averaging 63 bp in length (the largest block in the 3p13 sequence contains ~40 copies of this repeat). A portion of the GAG region of an HERV-E endogenous retroviral element comprises the homologous segment indicated in purple. The regions corresponding to the OR sequences are indicated in red. The ORs are numbered from 1 to 6, with the red and black colors of the numbers indicating the two groups of three highly homologous ORs. The regions denoted in turquoise are ~90% identical, but their structure varies due to a variety of Alu and L1 insertions (darker tones of turquoise). The positions of common interspersed repeats detected with RepeatMasker (http://ftp.genome.washington.edu/RM/RepeatMasker.html ) are indicated along the bottom of the dot-matrix for the 3p13 sequence and along the right edge for the PAC. The duplication unit is not repeated in the regions of the 3p13 sequence 5[prime] and 3[prime] of the segment shown in the dot-matrix or elsewhere in the PAC. The PAC's sequence ends at position 165491. Figure 7. A comparison of the FISH results with OR-containing cosmid 3 (LL03NC01-4B17) on human, chimpanzee, gorilla and orangutan chromosomes. The chromosomes are numbered according to the human karyotype. Cosmid 3 is part of the 3p13 contig (Figure 5). For each species, 10 metaphase spreads were analyzed and signals at each location were scored on scale of 1-4. The size of the symbols is proportional to the average score in 20 chromosomes, and thus reflects combination signal intensity and hybridization efficiency. Differences in the FISH patterns that were observed with two or more tested clones are discussed in the text. We used cosmid 202 to analyze the duplications in gibbon (Hylobates lar). This clone produces signals at six locations, on gibbon chromosomes 1, 4, 7, 8, 11 and 12. These sites correspond to regions on human chromosomes 7, 3, 8, 16, 11 and 3, respectively (16; data not shown). In contrast, the brightest and/or most frequently labeled locations in human are two sites on 3, 4p16, 8p, 11p, 11q and 16p (plus dimmer and less frequent signals on 3, 7 and 12). Thus, like orangutan, gibbon lacks cross-hybridizing sequences at regions corresponding to HSA 4p16, one of the brightest locations in the human genome. The sacred baboon (Papio hamadryas) is more distantly related to human than is gibbon, yet chromosome-painting experiments suggest that fewer chromosomal rearrangements exchanges have occurred during its divergence from human than have occurred along the branches separating gibbon and human (17). Cosmid 202 produced signals at at least eight locations on baboon chromosomes. The brightest lie in regions corresponding to HSA 1 (two sites), 3, 7, 16 and 19. No signals were observed in regions thought to correspond to HSA 4, 8 or 11 (PHA 5, 8 or 14), the most prominent locations in human.
DISCUSSION
Large duplications in the human genome
It has long been recognized that genomes are shaped by evolutionary processes that include the duplication of chromosomal segments (18). As a consequence of these past duplications, the human genome is riddled with repeats, ranging from the small and ubiquitous Alu-elements to large low-copy repeats (7,9-30 and references therein). The scope and complexity of the genomic duplications involving the OR family reported here are unprecedented, however, and are perhaps matched only by the extensive paralogy of regions near telomeres (30), where OR genes have also been found (7). Of the 44 locations identified here with OR-containing genomic clones, all but 11 (24%) contain sequences that cross-hybridize to one or more other locations in the genome. [If we consider only the 32 locations identified with two or more OR-containing clones, all but 5 (16%) contain sequences related to at least one other location in the genome.] Most of these regions were not previously known to have close relatives elsewhere in the genome. Of the 74 OR-containing BACs and cosmids we characterized, 57% hybridize to two or more locations. The similarity of these interchromosomal duplications is sufficiently high (~90%) and extensive (>150 kb) that cross-hybridization is detected under conventional FISH conditions. Note that hybridization of Alu and L1 elements and other high-copy repeats to the chromosomes is blocked in these experiments. Our most striking finding is the sequence homology among a set of >20 locations distributed on 13 chromosomes (Fig.
Homology among >15 locations comprises both OR and non-OR sequences
Several lines of evidence lead us to conclude that the sequences common to these many locations comprise a combination of OR and non-OR sequences. Pools of small OR-specific probes produced significant accumulations of FISH signals (6) at or near 21 of the 24 locations detected here with four or more multi-site clones (within the error of FISH localizations carried out in separate laboratories). Thus, OR sequences lie in most, if not all, of these related locations. However, the FISH signals with the multi-site clones were significantly more robust than those produced by the OR-specific probes, and this pattern was not significantly altered by suppression with unlabeled OR sequences. Clones spanning ~150 kb of a 250 kb contig of 3p13 (V. Brand-Arpon et al., submitted for publication) cross-hybridize to 17 locations. One of these clones, cosmid 88, lacks OR sequences, yet its multi-site pattern is similar to those of its OR-containing neighbors. PCR assays designed from the ends of various multi-site clones as well as in OR-free regions of the 3p13-contig demonstrate that many of these paralogous duplications include non-OR sequences that are sufficiently similar to serve as a template for PCR.
The molecular structures of paralogous regions
The structures of the paralogous regions are intriguingly complex. Clones derived from a cosmid-walk across an OR-cluster on 3p13 allow us to divide this region coarsely into four zones, each with a different chromosomal distribution. The region 3[prime] of the cluster of OR sequences is duplicated on two sites on chromosome 3. The region encompassing the three OR genes appears to be duplicated at least in part on 18 sites. By walking 5[prime] of the OR cluster, a third zone is reached that is duplicated in several additional locations (16p, 19p, 19q and 21q), but is less homologous to some regions (2, 13 and 14). Clones identified by screening the chromosome 3 library with the 5[prime] most distal end of the 3p13 contig identify a fourth zone. Sequences in this zone hybridize to 3q28-29 in addition to all the sites identified by the third zone except those on 19 and 21. Although its structure and sequence variants distinguish the 3p13 region from its paralogs (V. Brand-Arpon et al., submitted for publication), our cosmid walks have so far failed to reach the boundaries of the duplicated region.
A comparison of the sequence of a PAC ascribed to chromosome 11 and 106 kb of 3p13 sequence allows us to describe the structure of the second zone of homology within the 3p13 contig in more detail (Fig.
We remarked that clones derived from 3p13 produced signals at many sites, such as 3q13-21, 4p16, 7p22, 8p23, 11p15, 11q12-13 and 16p13.3, that were as bright or brighter than the signals at the native 3p13 location. These results indicate that portions of the homologous region are likely to exist in multiple copies at these locations. In contrast, several sites, such as on chromosomes 2, 13 and 14, were labeled only dimly and infrequently with the multi-site clones. These sites also showed little to no cross-hybridization with clones lying 5[prime] of the cluster of 3 OR sequences in the 3p13 contig. These observations, combined with the positive results on these chromosomes for PCR assays that lie between the OR sequences within the 22 kb unit, suggest that homology at these sites is restricted to one or a few copies of this duplication unit.
Several regions appear to be a composite of OR-containing segments with many relatives in the genome and single-copy or more restricted OR-containing segments. For example, 11p15, 9q32-34 and 13q21 contain sequences that cross-hybridize to many multi-site clones, as well as segments that behave as single-copy probes (compare Fig.
Atypical intensities of several of the multi-site clones at particular locations belie their genomic origin. For example, clone 5C11 hybridizes particularly intensely to 2p12-13, which is dimly labeled by most of the other clones, 663A2 gives a very bright signal at 2q23, 989A4 at 13q21, and 922E9 at 21q22. It is likely that these clones contain sequences that overlap the duplicated region, but extend into relatively unique sequence at these locations.
Plasticity in the cytogenetic arrangement of OR duplications in primates
Previously, we detected marked polymorphism among humans using a subtelomeric block of OR genes, with copy-number ranging from 7 to 11 and wide diversity in chromosome locations (7). In contrast, we observe little large-scale polymorphism among humans using the subset of the OR family surveyed here.
However, our data show that much of the OR subgenome has been in flux at the cytological level during the divergence of primates. Only chimpanzee and human show a similar pattern of duplications with the multi-site clones. The most significant difference between these species is diminution of the signal on 4p in chimpanzee relative to human and gorilla.
Our results illustrate that many more molecular changes distinguish the genomes of primates than are evident from the conservation in banding patterns among the great apes (31) or the contiguity of regions labeled with chromosome-paints in baboon. A variety of changes in location and intensity of OR-containing segments have occurred during primate evolution. Of the six primates analyzed, only three (human, chimpanzee and gorilla, have sequences cross-hybridizing to the multi-site clones on 4p. Signal is lacking in these same species on 16p, where a large block of cross-hybridizing sequences is detected in orangutan, gibbon and baboon. The most parsimonious explanation for these two changes is that they are the result of the translocation of a block of sequence from 16p to 4p along the branch prior to human-chimpanzee-gorilla divergence. This translocation probably involved only a portion of the region devoted to OR sequences on 16p, since sequences derived from the 5[prime] end of the 3p13-contig map to 16p in the three species, and some humans carry a block of other OR sequences on 16p (7). Other significant changes that have occurred during primate evolution are the loss of cross-hybridizing sequences from chromosome 7 in gorilla, from chromosome 11 in orangutan and baboon, and from chromosome 8 in baboon. In addition, the sequence has been gained (or simply retained) on chromosome 19, and two sites on chromosome 1 in baboon.
Our findings, combined with recent observations for other subtelomeric (7,30), pericentromeric (19-21) and OR (32; V. Brand-Arpon et al., submitted for publication) sequences, indicate that large-scale differences among primate genomes have been grossly underestimated. The frequently quoted 2-10% sequence variation underestimates the degree of genomic variation among these species. Few of the sites detected with the multi-site clones have been stable during the divergence of primates. Our findings are consistent with the idea that copy number differences (amplification/attrition) and contextual changes (location changes) of sequences contribute to phenotypic differences among organisms and may be more significant for speciation events than subtle differences at the nucleotide level. It will be important to compare the collection of expressed OR sequences to learn how these gross structural changes affect the repertoire of olfactory receptors in the nasal epithelium or testes of these species.
The structure and distribution of OR-containing paralogs suggest a hierarchical model for the duplication process
A multi-step process of duplications and rearrangements is required to explain our observations. OR sequences and VNTR blocks appear to have duplicated to form the larger ~22 kb duplication unit, which duplicated both intra- and interchromosomally. We suspect that the large blocks of the ~63 bp repeats at the junctions of the 22 kb duplication facilitate the recombination events that lead to duplications (and deletions) of OR genes. The pGAG, Alu and Line insertions within the units may prove useful to date the duplication events and establish the relationships among some paralogs. Copies of the 22 kb duplication unit appear to be a portion of a larger region that is shared by multiple chromosomes, since clones lacking this unit still cross-hybridize to many of the same locations. The duplicated zones appear to have been further modified by rearrangements and/or insertions of relatively unique DNA. In some cases, these sequences were included in subsequent duplications (such as the MLCK unit at the 3[prime] end of the 3p13-contig). It is also possible that some initial duplication units were large and were later whittled down to a dense array of key functional elements [analogous to some postulated processes in intrachromosomal gene amplification (33)].
The subterminal location of many OR-containing blocks may facilitate this process, because large-scale duplications and rearrangements could occur without affecting dosage-sensitive genes in the chromosome's interior. A disproportionate number (42%) of OR-containing duplicates is located in the 10% of the genome that lies closest to the ends of chromosomes. Our observations are consistent with the notion that terminal bands could serve as a nursery for generating diversity among a subset of this gene family (7).
Large genomic duplications provide a mechanism for changing the repertoire of ORs. Functionally redundant copies could evolve through mutations to new specificities or expression patterns. Many duplicates would degenerate into pseudogenes. It is not yet clear whether the selective pressure to expand or change the OR gene repertoire is the driving force behind the large-scale duplications we observe or if OR sequences were captured inadvertently in regions with a propensity to duplicate and are innocent bystanders in the process. Sequence comparisons among the OR sequences at each location and attempts to match expressed sequences to their chromosomal location should provide some answers.
Many locations appear to contain a combination of pseudogenes and genes. We showed previously a bias for intact ORFs on chromosomes 7, 16 and 17 (6), but none of the genomic clones analyzed here maps uniquely to 16 or 7; all cross-hybridize to other locations. The duplications in common with the 11-PAC and the 3p13 contig contain apparent pseudogenes. The fact that related pseudogenes in the two regions have specific mutations in common is an indication that clusters of pseudogenes are subject to duplication. It is not known if pseudogenes serve a function; some are transcribed (34). In any case, the duplications of pseudogenes may be a negligible genomic burden when weighed against the benefits of a process that can provide new substrates for diversifying the OR repertoire.
A subset of OR-clusters is significantly diverged from other OR-containing regions
We identified a number of OR-containing clones that map to single sites in the genome. The sequences in these clones must be sufficiently diverged that they fail to cross-hybridize to other OR-containing locations. Three of the locations detected with two or more independent clones, 17p, 19p and 11p, are locations of OR genes that are already the subject of detailed genomic analysis (32,35-38). The fact that these OR clusters are among the first to be characterized in detail is probably a consequence of the fact that mapping them was uncomplicated by extensive homology with other locations in the genome. Our results also indicate that 1q21-22, 1q44, 3p14-21 and 9q32-q34 also harbor relatively unique representatives of the OR family. Our results confirm the findings of Rouquier et al. (6), who detected OR sequences at three of these locations. To our knowledge, our study is the first to point to 1q21-22 as a location of OR genes.
OR-locations missing from this collection of genomic clones
Several sites that are known to harbor OR genes are not represented in the collection of BACs that we have characterized here. Our strategy missed clones encompassing a block of OR genes that is duplicated near many telomeres including 3qter, 15qter and 19pter (7). Also notably lacking from the collection are BACs representing the 6p21 site of OR genes amidst the MHC gene cluster (39) and the 14q11.2 site near the TCR-[alpha] genes (6,40; GenBank accession no. U85195). Clones encompassing these OR genes do not cross-react by FISH with other OR-containing regions (C. Friedman and B.J. Trask, unpublished data). It is therefore interesting to note that two OR-BACs were identified that hybridize near the TCR-[alpha] locus on 14q11.2, but also to a second site at 15q11.2. This result lends further support to the notion that clusters of OR genes are a complex combination of single-copy sequence and segments with homology to two or more sites, reflecting an evolutionary history of duplications and rearrangements.
Fraction of the human genome devoted to the OR family
Our results suggest that large genomic duplications are responsible for the distribution, size and diversity of much of the OR gene family. Differences among primates suggest that the portion of the genome devoted to the OR family has changed during primate evolution, possibly as a mechanism for evolving the repertoire of OR genes. If we conservatively estimate that each of the paralogous regions encompasses 125 kb (and the results with the 3p13 contig suggest that the regions in some locations may be even larger), then more than ~3 Mb or 0.1% of the human genome is occupied by the duplications and unique segments that encompass members of the OR gene family. The similarity of many of these large repeats may be maintained by sequence exchanges among the copies (e.g. by unequal cross-over or gene conversion). The multiplicity of these large regions of homology poses a severe challenge to sequencing the human genome. However, the large-insert clones we identify here should prove useful for mapping the clusters in more detail in order to determine how the exquisite transcriptional control of these diverse receptors is achieved.
Materials and Methods
Isolation of OR-containing genomic clones
A pool of OR-specific sequences was generated by PCR from genomic DNA using degenerate primers designed in the highly conserved regions of OR genes. The OR probe pool used to screen approximately four genomic equivalents of a total genomic BAC library (CITB 978SK) was generated using primers 12464 [5[prime]-A(G/C)(A/T/C/G)TATGACCGCTATGTGGCCATCTG] and 12465 [5[prime]-CACCACAG(A/T)CAGGTGGGA(G/T)(C/G)CACAGG] recognizing a sequence [(T/S/R)YDRYVAI] in the second cytoplasmic loop, and C(A/G)SHL(T/S)VV in TM6. The probe used to screen the chromosome 3 library (LL03NC01) was amplified using primers in the conserved domains TM2 [OR5B; PMY(F/L)FL(S/A/T/G/C)NLS] and TM7 [OR3B; M(L/F/V/I)NPF(I/M)Y(S/C)L] (8; V. Brand-Arpon et al., submitted for publication). Both libraries were spotted in duplicate on high-density filters. Probes were radiolabeled by random hexamer priming, and filters were washed after hybridization to a stringency of 65°C in 0.1× SSC, or in the case of the chromosome 3 filters as described by V. Brand-Arpon et al. (submitted for publication). BAC clones were secondarily screened by PCR using primers 12464 and 12465. PCR products of 33 of the 52 BACs that passed this secondary screen and were FISH-mapped were cloned and sequenced to confirm OR homology. A comparison of the resulting sequences will be published elsewhere (K. Chan and H. Shizuya, unpublished data). Chromosome 3 clones were confirmed to contain OR sequences by PCR using OR3B/OR5B primers or by Southern blot analyses of EcoRI-restricted cosmid DNAs using the OR3B/5B PCR product amplified from total genomic DNA, as described elsewhere (V. Brand-Arpon et al., submitted for publication). Clones were selected for FISH analyses randomly, with some bias to include representatives from the cosmid collection with different restriction enzyme patterns (V. Brand-Arpon et al., submitted for publication).
Additional clones were identified for FISH analyses by cosmid walking in the chromosome 3 specific library from the ends of OR-containing clones. The resulting 250 kb contig from 3p13 is described elsewhere (V. Brand-Arpon et al., submitted for publication). Additional clones, which may derive from any of three locations on chromosome 3, were isolated by screening the chromosome 3 library with the T3 end of cos 88 and the T7 end of PAC169 by hybridization (V. Brand-Arpon et al., submitted for publication). These clones were tested for OR content by PCR and/or Southern blot analyses as described above.
The LLNL-library designations (LL03NC01) for the chromosome 3 cosmids used in this study are 3-4B17, 5-4L5, 8-6A20, 11-10A14, 13-10K1, 21-19B17, 26-21N3, 28-23I22, 32-24I3, 36-24P4, 39-25D13, 40-25H11, 45-28K17, 48-31P3, 54-37B9, 57-39H23, 65-44F15, 66-45J14, 78-51I19, 79-52K19, 80-53J7, 81-53O7, 88-24I21, 96-54B7, 97-55A16, 111-18O2, 113-26H13, 117-39K8, 202-4A19, 204-4L5, 209-23E11, 215-28I23 and 219-47A22.
FISH
Metaphase spreads were prepared from PHA-stimulated peripheral blood lymphocytes of healthy humans. A few clones were additionally mapped on human lymphoblast cell cultures GM11525 and GM10977 (Coriell). Cell lines from ATCC served as representatives of chimpanzee (Pan troglodytes, CRL-1847), gorilla (Gorilla gorilla, CRL-1854), orangutan (Pongo pygmaeus, CRL-1850), gibbon (Hylobates lar, TIB-201) and sacred baboon (Papio hamadryas, CRL-1495).
FISH was performed as described elsewhere (41). Briefly, cosmid or BAC DNA was isolated on Qiagen tip-100 columns or an Autogen system, respectively, biotinylated by nick-translation and hybridized to metaphase spreads in the presence of excess unlabeled Cot1 DNA (added to suppress hybridization of interspersed repeats throughout the genome). The same stringency conditions were used for all experiments, i.e. hybridization at 37°C in 50% formamide, 2× SSC, 10% dextran sulfate, pH 7.0; washing for 15 min at 42°C in 3 changes of 50% formamide/2× SSC followed by 15 min at 42°C in 2× SSC. Hybridization sites were detected with two layers of avidin-FITC separated by a layer of biotinylated anti-avidin. The chromosomes were QFH-banded by DAPI staining. DAPI and FITC images were digitized separately, but in registration, using a Zeiss Axiophot microscope, a Ludl filter wheel equipped with separate DAPI and FITC excitation filters, a ChromaTechnology multi-band pass emission filter, 100× NeoFluar objective, and a Princeton cooled CCD camera (Kodak 1000×1317) operated via a custom script in ScanAlytics IPLab Spectrum 3.0 software. For the image shown, the FITC image was pseudocolored red for clarity, and the DAPI banding was displayed as grey values. Hybridization signals were analyzed in at least five and more typically 10 metaphase cells per probe. Clones localizing to a single location were scored at the microscope, and only a few images were collected. For clones mapping to multiple sites, signals were scored from digitized images to avoid missing dim signals due to photobleaching during the analysis. The contrast and gain were varied interactively so that both dim and bright signals could be scored. Signal intensity was scored on a scale of 1-4, with four corresponding to a large bright signal. The average signal intensity at each location was calculated. This value measures a combination of signal intensity and labeling efficiency, because a chromosome showing no signal at a particular location was assigned the value 0.
Eleven BAC clones were blindly mapped twice using independent DNA isolates, allowing us to assess the accuracy of clone tracking and chromosomal assignments. The same result was obtained for seven duplicate pairs, including three clones with multi-site patterns. Imprecision in our assignments to 9q32-33 and 9q33-34 were identified with three duplicates. Therefore, we have combined all clones mapping to these locations into a single group. Finally, inter-sample variation in signal intensity was noted with 302B6, which produced signals at ~17 sites in one sample (as plotted in Fig.
These karyotypes shown in Figures
PCR typing of monochromosomal somatic cell hybrid panel
The Coriell monochromosomal somatic cell hybrid panel #2 was PCR typed using the following primer pairs designed from the ends of OR-containing cosmids: the T7 end of cosmid 32 (32MF, 5[prime]-ATCTCATGATCTGTTCCATCC; 32IR, 5[prime]-ATTCCAGTTAAAGGCATAACG), the T3 end of cosmid 32 (32CF, 5[prime]-TCCTAATATCACCGTGGCTC; 32CR, 5[prime]-AATATGATCACAGGGTGTACC), the T3 end of cos 5 (5AF, 5[prime]-TCTGGCCATGTGATTGCTGC; 5AR, 5[prime]-TATAAAGCGATAACTGAGGCT), the T7 end of cos 5 (5BF, 5[prime]-ACCGGGAACCTTTGGAATGC; 5BR, 5[prime]-GAGAGGTAGGCCACCGGC) and the T3 end of cos 11 (11AF, 5[prime]-TGAACCAGCTGTCCTCAGC; 11AR, 5[prime]-TGTTCCCAAGCCACGGACC). In addition, several pairs of primers were designed in the repeat- and OR-free regions of clone ends or in the 106 kb sequenced contig from 3p13 (V. Brand-Arpon et al., submitted for publication) and typed on the hybrid panel. These PCR assays are denoted A-G in Figure
ACKNOWLEDGEMENTS
We thank Lee Rowen for advice on the sequence analysis. This research was supported by grants from the NIH (GM57070-01 and HG01464), the US Department of Energy (DE-FG03-96ER62173), the Programme Génome du CNRS, the Fondation pour la Recherche Médicale, and NATO.
REFERENCES
This article has been cited by other articles:
This page is run by Oxford University Press, Great Clarendon Street, Oxford OX2 6DP, as part of the OUP Journals
Comments and feedback: www-admin{at}oup.co.uk
Last modification: 13 Nov 1998
Copyright©Oxford University Press, 1998.
![]()
CiteULike
Connotea
Del.icio.us What's this?
![]()
![]()

![]()
![]()
![]()
A. Moreno-Estrada, F. Casals, A. Ramirez-Soriano, B. Oliva, F. Calafell, J. Bertranpetit, and E. Bosch
Signatures of Selection in the Human Olfactory Receptor OR5I1 Gene
Mol. Biol. Evol.,
January 1, 2008;
25(1):
144 - 154.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
Y. Gilad, O. Man, and G. Glusman
A comparison of the human and chimpanzee olfactory receptor gene repertoires
Genome Res.,
February 1, 2005;
15(2):
224 - 230.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
B. Malnic, P. A. Godfrey, and L. B. Buck
The human olfactory receptor gene family
PNAS,
February 24, 2004;
101(8):
2584 - 2589.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
C. Amadou, R. M. Younger, S. Sims, L. H. Matthews, J. Rogers, A. Kumanovics, A. Ziegler, S. Beck, and K. Fischer Lindahl
Co-duplication of olfactory receptor and MHC class I genes in the mouse major histocompatibility complex
Hum. Mol. Genet.,
November 15, 2003;
12(22):
3025 - 3040.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
T. Newman and B. J. Trask
Complex Evolution of 7E Olfactory Receptor Genes in Segmental Duplications
Genome Res.,
May 1, 2003;
13(5):
781 - 793.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
D. P. Locke, R. Segraves, L. Carbone, N. Archidiacono, D. G. Albertson, D. Pinkel, and E. E. Eichler
Large-Scale Variation Among Human and Great Ape Genomes Determined by Array Comparative Genomic Hybridization
Genome Res.,
March 1, 2003;
13(3):
347 - 357.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
J.A. BAILEY and E.E. EICHLER
Genome-wide Detection and Analysis of Recent Segmental Duplications within Mammalian Organisms
Cold Spring Harb Symp Quant Biol,
January 1, 2003;
68(0):
115 - 124.
[Abstract]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
M. P. Gygi, M. D. Ferguson, H. C. Mefford, K. P. Lund, C. O'Day, P. Zhou, C. Friedman, G. van den Engh, M. L. Stolowitz, and B. J. Trask
Use of fluorescent sequence-specific polyamides to discriminate human chromosomes by microscopy and flow cytometry
Nucleic Acids Res.,
July 1, 2002;
30(13):
2790 - 2799.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
J. M. Young and B. J. Trask
The sense of smell: genomics of vertebrate odorant receptors
Hum. Mol. Genet.,
May 15, 2002;
11(10):
1153 - 1160.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
J. M. Young, C. Friedman, E. M. Williams, J. A. Ross, L. Tonnes-Priddy, and B. J. Trask
Different evolutionary processes shaped the mouse and human olfactory receptor gene families
Hum. Mol. Genet.,
March 1, 2002;
11(5):
535 - 546.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
R. P. Lane, J. C. Roach, I. Y. Lee, C. Boysen, A. Smit, B. J. Trask, and L. Hood
Genomic Analysis of the Olfactory Receptor Region of the Mouse and Human T-Cell Receptor alpha /delta Loci
Genome Res.,
January 1, 2002;
12(1):
81 - 87.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
E. E. Eichler, M. E. Johnson, C. Alkan, E. Tuzun, C. Sahinalp, D. Misceo, N. Archidiacono, and M. Rocchi
Divergent Origins and Concerted Expansion of Two Segmental Duplications on Chromosome 16
J. Hered.,
November 1, 2001;
92(6):
462 - 468.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
H. Kouros-Mehr, S. Pintchovski, J. Melnyk, Y.-J. Chen, C. Friedman, B. Trask, and H. Shizuya
Identification of Non-functional Human VNO Receptor Genes Provides Evidence for Vestigiality of the Human VNO
Chem Senses,
November 1, 2001;
26(9):
1167 - 1174.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
E. Linardopoulou, H. C. Mefford, O. Nguyen, C. Friedman, G. van den Engh, D. G. Farwell, M. Coltrera, and B. J. Trask
Transcriptional activity of multiple copies of a subtelomerically located olfactory receptor gene that is polymorphic in number and location
Hum. Mol. Genet.,
October 1, 2001;
10(21):
2373 - 2383.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
R. P. Lane, T. Cutforth, J. Young, M. Athanasiou, C. Friedman, L. Rowen, G. Evans, R. Axel, L. Hood, and B. J. Trask
Genomic analysis of orthologous mouse and human olfactory receptor loci
PNAS,
June 19, 2001;
98(13):
7390 - 7395.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
G. Glusman, I. Yanai, I. Rubin, and D. Lancet
The Complete Human Olfactory Subgenome
Genome Res.,
May 1, 2001;
11(5):
685 - 702.
[Abstract]
[Full Text]
![]()
![]()
![]()

![]()
![]()
![]()
T. Shiina, A. Ando, Y. Suto, F. Kasai, A. Shigenari, N. Takishima, E. Kikkawa, K. Iwata, Y. Kuwano, Y. Kitamura, et al.
Genomic Anatomy of a Premier Major Histocompatibility Complex Paralogous Region on Chromosome 1q21-q22
Genome Res.,
May 1, 2001;
11(5):
789 - 802.
[Abstract]
[Full Text]
![]()
![]()
![]()

![]()
![]()
![]()
J. E. Horvath, S. Schwartz, and E. E. Eichler
The Mosaic Structure of Human Pericentromeric DNA: A Strategy for Characterizing Complex Regions of the Human Genome
Genome Res.,
June 1, 2000;
10(6):
839 - 852.
[Abstract]
[Full Text]
![]()
![]()
![]()

![]()
![]()
![]()
Y. Ji, E. E. Eichler, S. Schwartz, and R. D. Nicholls
Structure of Chromosomal Duplicons and their Role in Mediating Human Genomic Disorders
Genome Res.,
May 1, 2000;
10(5):
597 - 610.
[Abstract]
[Full Text]
![]()
![]()
![]()

![]()
![]()
![]()
E. E. Eichler, N. Archidiacono, and M. Rocchi
CAGGG Repeats and the Pericentromeric Duplication of the Hominoid Genome
Genome Res.,
November 1, 1999;
9(11):
1048 - 1058.
[Abstract]
[Full Text]
![]()
![]()
![]()

![]()
![]()
![]()
P. Mombaerts
Seven-Transmembrane Proteins as Odorant and Chemosensory Receptors
Science,
October 22, 1999;
286(5440):
707 - 711.
[Abstract]
[Full Text]
![]()
![]()
![]()

![]()
![]()
![]()
J. R. Korenberg, X.-N. Chen, Z. Sun, Z.-Y. Shi, S. Ma, E. Vataru, D. Yimlamai, J. S. Weissenbach, H. Shizuya, M. I. Simon, et al.
Human Genome Anatomy: BACs Integrating the Genetic and Cytogenetic Maps for Bridging Genome and Biomedicine
Genome Res.,
October 1, 1999;
9(10):
994 - 1001.
[Abstract]
[Full Text]
![]()
![]()
![]()

![]()
![]()
![]()
S. Rouquier, A. Blancher, and D. Giorgi
The olfactory receptor gene repertoire in primates and mouse: Evidence for reduction of the functional fraction in primates
PNAS,
March 14, 2000;
97(6):
2870 - 2874.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
R. M. Younger, C. Amadou, G. Bethel, A. Ehlers, K. F. Lindahl, S. Forbes, R. Horton, S. Milne, A. J. Mungall, J. Trowsdale, et al.
Characterization of Clustered MHC-Linked Olfactory Receptor Genes in Human and Mouse
Genome Res.,
April 1, 2001;
11(4):
519 - 530.
[Abstract]
[Full Text]
![]()
![]()
![]()

![]()
![]()
![]()
J. A. Bailey, A. M. Yavor, H. F. Massa, B. J. Trask, and E. E. Eichler
Segmental Duplications: Organization and Impact Within the Current Human Genome Project Assembly
Genome Res.,
June 1, 2001;
11(6):
1005 - 1017.
[Abstract]
[Full Text]
[PDF]
![]()
This Article ![]()
![]()
Abstract
![]()
FREE Full Text (PDF)
![]()
Alert me when this article is cited
![]()
Alert me if a correction is posted
![]()
Services ![]()
![]()
Email this article to a friend
![]()
Similar articles in this journal
![]()
Similar articles in ISI Web of Science
![]()
Similar articles in PubMed
![]()
Alert me to new issues of the journal
![]()
Add to My Personal Archive
![]()
Download to citation manager
![]()
Search for citing articles in:
ISI Web of Science (76)
![]()
Request Permissions ![]()
Google Scholar ![]()
![]()
Articles by Trask, B. J.
![]()
Articles by Giorgi, D.
![]()
Search for Related Content
![]()
PubMed ![]()
![]()
PubMed Citation
![]()
Articles by Trask, B. J.
![]()
Articles by Giorgi, D.
![]()
Social Bookmarking ![]()
![]()
What's this?