Skip to main content

Unusual linkage patterns of ligands and their cognate receptors indicate a novel reason for non-random gene order in the human genome



Prior to the sequencing of the human genome it was typically assumed that, tandem duplication aside, gene order is for the most part random. Numerous observers, however, highlighted instances in which a ligand was linked to one of its cognate receptors, with some authors suggesting that this may be a general and/or functionally important pattern, possibly associated with recombination modification between epistatically interacting loci. Here we ask whether ligands are more closely linked to their receptors than expected by chance.


We find no evidence that ligands are linked to their receptors more closely than expected by chance. However, in the human genome there are approximately twice as many co-occurrences of ligand and receptor on the same human chromosome as expected by chance. Although a weak effect, the latter might be consistent with a past history of block duplication. Successful duplication of some ligands, we hypothesise, is more likely if the cognate receptor is duplicated at the same time, so ensuring appropriate titres of the two products.


While there is an excess of ligands and their receptors on the same human chromosome, this cannot be accounted for by classical models of non-random gene order, as the linkage of ligands/receptors is no closer than expected by chance. Alternative hypotheses for non-random gene order are hence worth considering.


One of the most striking discoveries in the post-genomic age has been the amount of non-random gene positioning in eukaryotic genomes [1]. In the human genome, for instance, highly/broadly expressed genes cluster [25]. Likewise in yeast co-expressed genes tend to reside together [6] and such pairs tend also to be retained together over evolutionary time more than expected, given the intergene distance between them [7]. Blocks of broadly expressed mammalian genes also seem to be preserved over evolutionary time more than expected [8]. In Caenorhabditis [911], Drosophila [1214] and Arabidopsis [15], to name but three, there exists further evidence for expression clusters of some variety. These results all suggest that eukaryotic genomes are organised in a manner that permits co-expression or co-ordinate expression. Evidence also suggests linkage of functionally related genes, although on this issue the evidence is more equivocal, not least because of an ambiguity as to what "functionally related" can mean. On the one hand, in numerous eukaryotic genomes, genes from the same metabolic pathway cluster more than expected by chance [16] (for detailed case history see [17]). Likewise, linked co-expressed genes in yeast often fall within the same MIPs (Munich Information Centre For Protein Sequences) category [6] or the same Gene Ontology (GO) classification [18].

These results are so striking because they so profoundly overturn the long held assumption that genes are randomly located around eukaryotic genomes. This is not to say that possible exceptions were not considered prior to the sequencing of the complete genome. They were, however, typically dismissed as being unrepresentative or uninteresting either because they were clearly the product of tandem duplication (hox cluster, globin cluster) or were associated with weird genetics (imprinted clusters) or genes that are otherwise exceptional (e.g. clustering of rRNAs). Not all such suggestive examples could so easily be dismissed however. Here we concentrate on one class, linkage of ligands to their cognate receptors. This issue is worth systematic analysis, not least because in yeast it has recently been shown that genes whose proteins interact to form stable complexes are linked more often than expected by chance [19].

That ligands and their receptors may be linked was observed independently by several workers. Cooper [20], noting that the linkage of ligands to receptors may be common, highlights the examples of transferrin and transferrin receptor on chromosome 3q, as well as apolipoprotein E and the low density lipoprotein receptor both on chromosome 19. He also rightly cautions, however, than one can find numerous cases where ligands and receptors are not linked. Similarly, Lennard et al. [21] note the linkage of the three ligands in the interleukin 1 cluster (IL1 alpha, beta and receptor antagonist) to the two receptors [21]. The linkage of ligands to receptors has even proven to have some predictive power. Wang et al. [22] noticed that hepatocyte growth factor (HGF) and its MET receptor were both on 7q. Noting too the presence in 3p21 of both macrophage stimulating factor (MST1, a member of the same gene family as HGF), and RON (a member of the MET receptor family), they hypothesised that RON might be MST1's receptor [22]. This, in turn, they demonstrated to be the case (RON's alias is now MST1R) [22]. Popovici et al. note several of the above examples and also point to a total of 14 incidences of linked genes involved in the same pathway, not necessarily as ligand-receptor couplings [23].

This evidence prompts two questions. First, is it true that there is something odd about the linkage patterns of ligands and their cognate receptors? Second, if it is true, why might this be so? Prior authors have also suggested that linkage of ligands to receptors might be functionally important. Haig [24], observing two of the above cases (interleukin 1 and transferrin), notes that close proximity could enable linkage disequilibrium between alleles at the ligands and receptors. This linkage disequilibrium would potentially enable the spread of rare allele combinations for which there exist particular epistatic interactions. These may, Haig suggests, act as selfish maternal effect lethals, an example of which has been described in mice [25, 26]. This theory may be seen as a special case of a more general theory for linkage based on preservation of linkage disequilibrium under epistasis [2730]. One might also conjecture that ligands and receptors might at times need to be co-expressed [see e.g. [31]], so very close linkage might be beneficial for this reason as well.

If selection does act on the location of ligands and receptors (either to permit co-expression or to maintain linkage disequilibrium), then we should predict, from the above models, that when two such genes are on the same chromosome they should also be, on average, physically closer than would be expected by chance. To this end we ask two questions. First, is the mean distance between ligands and their linked receptors shorter than expected by chance? As this mode of analysis could miss an excess of cases with very tight linkage, we additionally ask whether the number of incidences of linkage within a given window size (1 Mb, 2 Mb etc.) is higher than expected by chance.

Results and discussion

No evidence for close proximity of ligands and receptors

If ligands and their cognate receptors were under selection to be in close physical proximity, we should find that the mean distance between them should be smaller than expected by chance. To test this, we examined ligand-receptor pairs from the DLRP database [32, 33]. All analyses were also performed for an augmented dataset, which additionally contains two 'cherry-picked' cases highlighted by Cooper [20] (see Methods). Contrary to expectations, the mean distance between ligand and receptor in the non-augmented data set (64.579 Mb) is higher than that found in randomized genomes (56.344 Mb; P = 0.733). The same pattern is found in the augmented data set (real = 63.466 Mb, randomized = 56.165 Mb, P = 0.709). Note that a ligand can have many receptors and that this is factored into the analysis through the randomization protocol.

It may, however, be the case that there exist a number of ligand-receptor pairs that are much closer to each other than expected. To examine this we compared the number of ligand-receptor pairs within some critical distance of each other and compared this with the number expected by chance. In these simulations we permuted genes only within the chromosomes within which they are found, so as to ensure that the number of ligand-receptor pairs were the same in the randomized sets as in the real data set. As can be seen (Table 2) in neither data set do we find evidence for anything other than a pattern of random linkage [see Additional file 4].

Table 2 The number of occurrences of a ligand receptor pair within a given distance of each other (in Mb) for the real genome compared with 10000 randomized genomes. In this instance the first three results columns refer to the dataset excluding the two examples identified by Cooper.

It has been previously established that broadly expressed genes cluster [3, 4]. Might our simulations have produced misleading results by permitting genes of the ligands and receptors to reside in any chromosomal location? To examine this possibility we considered randomizations in which genes are swapped exclusively with ones of the same breadth of expression. None of the above results are qualitatively affected (Table 3). Using a different bin size to classify breadth of expression appears to have no effect on the results (Table 3). Likewise, permitting ligands and receptors to be located in the genome at the locations of other ligands and receptors does not affect any conclusions [see Additional file 4].

Table 3 Results of randomizations controlling for breadth of expression. The "Aug" (Augmented) data set is that containing the two ligand -receptor sets nominated by Cooper [20]. Bin size indicates the span of breadths of expression considered to be the same in the randomizations. Bin size one implies that only genes of the same breadth were switched with each other. Distances are measured in Mb. P-values are estimated by comparison of observed data with expectations obtained from randomised genomes.

The above results indicate that there is no evidence for selection for clustering of ligand and receptor. For this reason we reject a model positing epistasis between alleles of ligand and receptor as a general force acting on genomic location of these genes. Moreover the lack of tight clustering suggests that we are not witnessing clustering to enable co-regulation (by ensuring that genes are co-localised in the same chromatin block). The model suggested by Haig [24] is not, however, necessarily falsified by the above results, as he postulates selection on disequilibrium only if the genes might be involved in maternal-foetal interactions. Such a model is hard to falsify in the absence of segregation/viability data from appropriate haplotypes. However, we can note that if we further restrict our data sets to those in which either the ligand or one of the receptors is placentally expressed, the qualitative patterns described above are unaltered [see Additional file 4]. We find, therefore, no evidence for close linkage of ligand and receptor when involvement in maternal-foetal interactions might be a possibility.

An excess of ligand-receptor pairs on the same human chromosome

Above we asked whether ligands and their receptors are more closely linked than expected by chance. We can also ask if ligands and their receptors are more commonly linked (i.e. on the same chromosome) than expected by chance? In an unbiased unaugmented human data set (i.e. without the addition of the two sets highlighted by Cooper [20], see Methods) we observe 23 such pairings but expect on average 13.71 (P = 0.015). When we include the two extra sets the P value, as expected, is reduced: we observe 25 pairs but expect on average 13.8 (P = 0.005). These results support the view that in the human genome linkage of a ligand to at least one of its cognate receptors is more common than would be expected by chance. However, the majority (approx 78%) of ligands are not linked to any of their receptors, so this excess should not be considered a strong rule (although, as already noted, in special cases it has had predictive power).

No evidence for an excess of ligands-receptor pairs on the same chromosome in mouse

To ask whether the patterns observed in the human genome are also found in the mouse genome we constructed three mouse data sets and applied the three randomization protocols to each. The first two data sets are the ortholog equivalents of our two human data sets purged of duplicates by either a) Blasting or b) Blasting and removal by physical proximity (in the human genome) of ligands or receptors. That is, if two ligands were in close proximity in the human genome, even if not identified as sequence related, we would remove one before considering the location of the mouse orthologs. However, as it is possible that some ligand clusters might be unique to mouse, we additionally purged the more stringent of the above two of any groupings of ligands or receptors seen in the mouse genome.

As it happens, no matter which data set one employs or which randomization method is performed, there is not even a remote hint that ligands and cognate receptors occur more commonly on the same mouse chromosome than expected by chance [see Additional file 5]. For example, in the equivalent of the human data set purged of duplicates by Blast alone, we observe 19 ligand-receptor pairs on the same chromosome and expect 19.38 (P = 0.56) in a randomization in which the ligands and receptors can assume any genomic position currently associated with a gene. In this analysis the mean distance between ligand and receptor is 57.3 Mb but is 43.4 in the randomizations (P = 0.937). At no specified distance do we find more pairs than expected by chance [see Additional file 5]. Controlling for breadth of expression makes no difference to this conclusion [see Additional file 6]. The list of linked genes from the data set equivalent to the human set presented in Table 1 is presented in Table 4. In this set 16 ligand-receptor pairs co-occur on the same chromosome, with 14.5 expected. Only four ligand-receptor pairs are in common in the two comparable data sets.

Table 1 Incidences of occurrence, on the same human chromosome, of a ligand with one of its cognate receptors after removal of tandem duplicates by Blast and by the physical proximity method. The distance is defined as the span between the mid-positions of the ligand and the mid-position of the receptor.
Table 4 Incidences of occurrence, on the same mouse chromosome, of a ligand with one of its cognate receptors after removal of tandem duplicates by Blast and position methods. The distance is defined as the span between the mid-positions of the ligand and the mid-position of the receptor. The receptors indicated with a Y in the conserved linkage column are those that are also on the same chromosome as the same ligand in the comparable human genome set (i.e. those in Table 1).

Explaining the data: interesting biology or statistical artefact?

We find no evidence in mouse or man that ligands and their receptors are more closely linked on average than expected by chance. However, we do find that there are more ligand-receptor pairs on the same human chromosome than expected, a feature not found in mouse. From the above results we are faced with two possible explanations for the human data. First, that it is just a statistical blip, possibly owing to some subtle bias in the original data set (note for example, that MST1R was analysed as a potential receptor because of its linkage to MST [22]). Second, that the excess of ligand-receptor pairs on the same chromosome is the product of some deterministic force, that for some reason does not apply, or is not strong enough, in rodents. This might either be direct selection favouring the persistence of co-occurrence or a deterministic bias in the creation of co-occurrence.

That the pattern is found in humans rather than mice argues against a direct selective benefit for co-retention on the same chromosome. This is owing to the fact that the effective population size of the human population is most probably much smaller than that of mice. As such, according to the nearly-neutral theory, the efficacy of selection should be higher in mice (see also [34, 35]). Hence, if the pattern was owing to selection directly favouring co-occurrence, it is more likely to be observable in mouse rather than human, all else being equal. Moreover, it is also hard to see what direct selective benefit might accrue from co-occurrence in weak linkage. In particular, co-regulation of ligands and receptors (which is not clearly expected in the first place) would likely require much tighter linkage than observed here. We also find no evidence for a stronger similarity in the breadth of expression of the linked ligands and receptors than those on different chromosomes. We considered the difference in breadth of expression between ligand and receptor normalised by the mean of the two. Linked genes are of no more similar breadth of expression (mean difference for linked genes 0.75 +/- 0.12, for unlinked 0.686 +/- 0.036, t-test, P = 0.58).

Segmental duplication and the balance hypothesis: an alternative hypothesis for non-random gene order

A notable feature of our data set is that there are numerous cases in which a ligand-receptor pair in linkage is matched by at least one other paralogous pair also in linkage. If we define genes belonging to the same Hovergen [36] family as paralogs, then we can identify the following linked paralogous pairs from Table 1: HGF/MST1 (ligands) with MET/RON (receptors); DLL1/DLL3 with NOTCH3/NOTCH4; FGF1/FGF2 with FGFR3/FGFR4; FGF8/FGF17 with FGFR1/FGFR2 (see also [23]). Note too that FGF18 is linked to FGFR4 and to FGF1 and is sequence related to FGF8 and FGF17.

It has been argued that co-paralogy of gene pairs involved in the same pathway (of which ligand-receptor pairs are but one example) appear to be unusually common [23]. This finding is also in accord with much recent evidence suggesting that the human genome (and the vertebrate genome more generally) may be a mosaic of old large block duplications (i.e., duplications of large chunks of DNA sequence) and/or the result of whole genome duplications [37], with several of the above paralogous groups being claimed to be the result of such duplication events: FGFR 1, 2, 3 and 4 are in paralog clusters on human chromosomes 8p, 10q, 4p, and 5q respectively [38]; NOTCH 3 and 4 also appear to be in paralog clusters on chromosomes 6 and 19, with NOTCH 1 and NOTCH 2 being two further duplications of the same block [39]; HGF/MST1 and MET/RON were also previously described as belonging to co-paralogous groups [23].

Following an earlier hypothesis [19], we would like to suggest an hypothesis to explain our results that is based on the occurrence of block duplications [40]. If some ligands and receptors require an appropriate balance in their titres, then one could expect that a mutation resulting in a block duplication containing one of the pair (e.g., the ligand) might be more likely to spread through the population if it also duplicates the other (e.g., the receptor). Such co-duplication is most likely if ligand and receptor happen to be linked, while unlinked ligand-receptor pairs are less likely to be successfully duplicated. Our hypothesis may be considered as being a form of the balance hypothesis [41, 42], which supposes that proteins involved in mutual interactions need to have their titres appropriately balanced. Direct evidence for this proposition has been described in yeast, in which it is also reported that the need for balance might explain the lack of duplicability of the genes involved in complex formation [43].

This hypothesis appears tenable in the current context for several reasons. It is, for example, suggestive that the paralogous pair sets tend to be more closely linked, although the statistic is on the edge of significance (Median Test, ChiSquared = 3.59, P = 0.058, df = 1). This would be expected if constraints exist on the upper size limit of the block duplications. It may also be notable that for the many of the above genes there is evidence for dosage sensitivity, as required by the balance hypothesis (Table 5). Moreover, as re-arrangements tends to be especially common in mice [4446], any linked pairs possibly generated by block duplication are more likely to be split up, making the genome more like random.

Table 5 Evidence for or against dosage sensitivity of ligand and receptors that were possibly the source for or the consequence of block duplication

The above hypothesis also predicts that ligands and their linked receptors might be duplicated at the same time. While this can be approximately established by phylogenetic methods, these do not constitute a perfect test, as they fail to establish whether the pairs were duplicated in a block together and furthermore, genes known to be co-duplicated are very commonly not identified as such by phylogenetic methods [47]. Nonetheless, we have surveyed the available data and prior analyses and fail to find any data that contradicts the hypothesis that when a given receptor duplicated the relevant ligand did as well [see Additional file 7]. The ligands FGF1 and FGF18 are, however, probably the result of an ancient duplication that occurred independent of the receptor [48].

In some of the incidences reported here, a case for co-duplication has already been made. The linkage of the FGFs to their receptors has previously been argued, from phylogenetic data, to be owing to block or whole genome duplications [48]. This view is supported by our inspection of the phylogenetic tree of the FGFR family as presented in Hovergen [36], which suggests that at the base of the vertebrates there was one receptor which duplicated to produce the ancestors of FGFR1/2 and FGFR3/4. Duplication of both ancestral sequences then occurred very shortly after (prior to the divergence of the fish), leaving FGFR1 and FGFR2 as nearest paralogs, and FRGR3 and FGFR4 as nearest paralogs. If there was co-duplication of the receptors, we should expect to see FGF1 and FGF2 as nearest paralogs and FRF8 and FGF17 as nearest paralogs, with, in both incidences, duplication occurring near the base of the vertebrates. The nearest paralog relationships are indeed upheld [48]. Furthermore, in both instances the duplication occurred prior to the divergence of fish, as predicted.


In sum, we have described a novel pattern of co-localisation on the same chromosome of genes whose products interact, which cannot obviously be accounted for either by known models for co-ordinate regulation, nor by selection for linkage disequilibrium. The pattern may in part reflect a past history of block duplication. A version of the balance hypothesis is worth considering as underpinning to explain the results.

Tests of this hypothesis should be possible in the future. We should in principle be able, with fuller knowledge of gene order in many mammals, to reconstruct the past history of duplication and gene order re-arrangements that occurred through mammalian history. The model predicts an excess of block duplications in which both ligand and cognate receptor are found, as well as excess in which neither are found, but a dearth of those with one, but not the other. The model also predicts a general weakening of this initial signal with increasing numbers of inter-chromosomal re-arrangements, as the hypothesis proposes only an initial filter of block duplications, not ongoing direct selection to maintain linkage.


Data set assembly and curation

The table of ligand-receptor partners were extracted from the DLRP database [32, 33]. This specifies for any given single ligand the corresponding receptor or set of receptors. The genes here are referred to by gene name and Unigene id number, by reference to an old release of Unigene. These entries we updated to the current release for Homo sapiens, UniGene Build #175. For each Unigene number and gene name, the relevant Unigene page was identified [49]. If the entry remained in the new build all details were left unchanged, except in three cases where there exists a gene by the same name as that in the original dataset, at the same genomic location as the given Unigene entry, but in a separate Unigene class. In these cases the Unigene entry with matching name was employed. If the old entry had been retired then a) if only one new entry is available this was used, b) if multiple entries were found (i.e. the cluster has split), then the one with the gene name identical to that of the old entry was used, c) if no entry had the same name but all entries were at the same genomic location the entry with the most abundant sequence data was used, d) if no unambiguous match could be found the entry was eliminated. If this was the ligand then the whole entry was deleted. In a few cases separate ligand receptor blocks are collapsed to the same Unigene entry in build #175 (e.g. FGFR and FGFRB). In this instance one of the two sets was eliminated. The original file also contains a number of entries in which the ligand alone is given, with no receptor. In these cases the entry was deleted.

From this new set of Unigene identities Entrez gene [50] was searched with the current Unigene id being posted. From here we recovered a) the Entrez/LocusLink gene name b) the Entrez/LocusLink id. If the LocusLink/Entrez gene name was different from the Unigene name then the pairs were examined at LocusLink to determine that the names were synonymous. In all cases this proved to be so. From here we obtained the physical location in the NC Genbank files for each chromosome. The cDNA source annotation of each Unigene entry was employed to determine whether the gene was placentally expressed.

The data set does not include the two cases highlighted by Cooper [20]: transferrin and its receptor and apolipoprotein E and its receptor. We therefore consider a second data set in which we add these two. This is, however, problematic as we are adding only "cherry picked" data. It is thus much more likely that we should find close linkage in the expanded set compared with the original set, owing to the non-random nature of the addition to the data set. Nonetheless, should we find an absence of an effect in this expanded set, this would make for stronger evidence against the hypothesis of an over-abundance of close linkage of ligands and their receptors.

The set so defined has numerous clusters of sequence-related ligands and receptors. Such clusters are likely to have arisen from tandem gene duplications, and thus individual genes cannot be treated as independently positioned. To eliminate the effects of tandem duplication we perform an all versus all blast (with E < 0.01) of the coding sequences defined from the RefSeqs for each gene. For each pair of putative duplicates on the same chromosome one of the two was randomly selected to be removed. In a tandem cluster with more than two duplicates only one gene was considered. Using E < 0.2 resolves to the same data set.

With this approach, however, a few well described duplicate clusters are not identified. For example, there remain 6 ligand-receptor pairs associated with the 2q14 cluster of three ligands (IL1 alpha, beta and the receptor antagonist) and their two receptors (IL1R1 and Il1R2) in 2q12. However, while Blast fails to reveal either of these clusters as duplicate clusters, this contradicts the conventional wisdom, based on close analysis of gene structure, function and conserved functional parts, that they are both duplicate arrays [51, 52]. The problem in this instance is most probably that interleukins and their receptors tend often to be fast evolving, hence liable to avoid detection as duplicates unless they are relatively modern duplicates. Indeed this inability of Blast to identify orthology/paralogy of fast evolving genes has recently been well demonstrated [53].

To eliminate such problems we additionally remove one of a pair of ligands within 1 Mb of each other. Likewise we remove one of a pair of receptors should the receptors occur within 1 Mb of each other. In nearly all cases the ligands in the cluster also bind the same receptors and vice versa. One exception is Insulin-like growth factor 2 (Igf2) and Insulin, which, while sequence related and very closely linked, bind different receptors. In both cases the receptors are unlinked. Given the sequence relatedness we remove one of the two. In effect we are then asking about a tendency for a ligand cluster to be linked to a receptor cluster. The final data set specifies 108 ligand-receptor sets (106 in the non-augmented set) [see Additional file 1]. Note that most ligands have more than one receptor. When then we refer to ligand receptor "pairs," we refer to incidence in which a ligand is linked to one of its receptors. We also performed the same analyses as given below on a data set in which duplicates are defined exclusively by reference to Blast scores [see Additional file 2]. A list of linked ligands and receptors in this data are provided [see Additional file 3]. Using both data sets, in both augmented and non-augmented form (i.e. with or without the two ligand-receptor pairs highlighted by Copper (1999)) we obtain qualitatively identical results [see Additional file 4].

Mouse data

For analysis of the patterns in mice we identified the orthologs of the human ligands and receptors by reference to the MGI curated set of mouse-human orthologs [54]. For each human gene, the locus link id was cross referenced to the mouse ortholog. Thirteen genes lacking an ortholog were removed. Mouse locus link numbers were employed to access RefSeq numbers, unigene references and the chromosomal locations. Breadth of expression was derived from Unigene cDNA source annotations. The compilation of all mouse genes and their position used in the randomization was derived from MartView [55] at Ensembl requesting those with described LocusLink ids. The positions of the ligands and receptors were found by cross-referencing their LocusLink ids to this Ensembl data set. The few that failed to be resolved by this method were ascribed a position by Blasting their RefSeq against the complete mouse genome [55]. For the randomizations only those genes with well resolved genomic locations were employed (24742 genes).

Randomization and statistics

To ask whether there are more ligand-receptor pairs on the same chromosome than expected by chance, we calculate the observed number and compare this with simulants in which we randomly permute the positions of all ligands and receptors. It is unclear on a priori grounds, however, what should be the null model for the randomization. We consider three possible models.

First we suppose that a ligand or receptor can occur in any location in the genome currently occupied by a protein coding gene. In this instance, in the human genome, the positions permitted in the randomizations correspond to the annotated positions of the 24,300 protein coding genes in the NC_0000n files for the human genome (n from 01–23).

Second, we assume that a ligand or receptor can occur in any location in the genome currently occupied by a protein coding gene with the same or comparable expression breadth. For each of the genes in the complete human and mouse sets we identified the Unigene id by following the LocusLink page pertinent to each gene and identified the breadth of expression in the same way as for the ligand-receptor set. Genes were placed in bins of 0–4, 5–9 tissues etc in which they were expressed. We consider two bin classification systems. In both, ligands and receptors in these randomizations were permitted to reside in the same location as any gene in the same bin from the complete human set.

Third, we suppose that there is something unique about ligands and receptors, such that each ligand or receptor can only be relocated to the position of another ligand or receptor. Breadth of expression is ignored as this too greatly constrains the randomizations.

For each of the above randomization protocols we determined the number of ligand-receptor pairs on the same chromosome. Significance (P) was determined from P = (r+1)/(n+1), where r is the number of simulants with the same or greater number of ligand-receptor pairs than observed in the real data and n is the number of simulants (10,000 in all instances), this being the unbiased estimator [56, 57].

To determine whether the ligand-receptor pairs that we observe on a given chromosome are more closely linked than expected we perform two analogous sets of simulations. In the first we calculate the mean distance between these pairs and compare this with the mean of the simulants described above. In the second we ask about the number of ligand-receptor pairs within a given distance of each other (e.g., within 1 Mb). We then compare this number to the mean number found after permuting all genes within the chromosomes within which they are found. By permuting on the same chromosome we control for the number of ligand-receptor pairs on the same chromosome.

All results prove to be insensitive to which of the three randomization null models is employed. Unless stated otherwise a result in the text relates to the first model in the text. All results can be found in attached files, for humans [see Additional file 4] and for mice [see Additional file 5] [see Additional file 6].





degrees of freedom


  1. 1.

    Hurst LD, Pal C, Lercher MJ: The evolutionary dynamics of eukaryotic gene order. Nature reviews Genetics. 2004, 5: 299-310. 10.1038/nrg1319.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Caron H, van Schaik B, van der Mee M, Baas F, Riggins G, van Sluis P, Hermus MC, van Asperen R, Boon K, Voute PA, Heisterkamp S, van Kampen A, Versteeg R: The human transcriptome map: Clustering of highly expressed genes in chromosomal domains. Science. 2001, 291: 1289-1292. 10.1126/science.1056794.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Lercher MJ, Urrutia AO, Hurst LD: Clustering of housekeeping genes provides a unified model of gene order in the human genome. Nature Genetics. 2002, 31: 180-183. 10.1038/ng887.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Lercher MJ, Urrutia AO, Pavlicek A, Hurst LD: A unification of mosaic structures in the human genome. Hum Mol Genet. 2003, 12: 2411-2415. 10.1093/hmg/ddg251.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Versteeg R, van Schaik BDC, van Batenburg MF, Roos M, Monajemi R, Caron H, Bussemaker HJ, van Kampen AHC: The human transcriptome map reveals extremes in gene density, intron length, GC content, and repeat pattern for domains of highly and weakly expressed genes. Genome Res. 2003, 13: 1998-2004. 10.1101/gr.1649303.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  6. 6.

    Cohen BA, Mitra RD, Hughes JD, Church GM: A computational analysis of whole-genome expression data reveals chromosomal domains of gene expression. Nature Genet. 2000, 26: 183-186. 10.1038/79896.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Hurst LD, Williams EJ, Pal C: Natural selection promotes the conservation of linkage of co-expressed genes. Trends Genet. 2002, 18: 604-606. 10.1016/S0168-9525(02)02813-5.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Singer GAC, Lloyd AT, Huminiecki LB, Wolfe KH: Clusters of co-expressed genes in mammalian genomes are conserved by natural selection. Mol Biol Evol. 2005, 22: 767-775. 10.1093/molbev/msi062.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Lercher MJ, Blumenthal T, Hurst LD: Coexpression of neighboring genes in Caenorhabditis elegans is mostly due to operons and duplicate genes. Genome Res. 2003, 13: 238-243. 10.1101/gr.553803.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  10. 10.

    Roy PJ, Stuart JM, Lund J, Kim SK: Chromosomal clustering of muscle-expressed genes in Caenorhabditis elegans. Nature. 2002, 418: 975-979.

    CAS  PubMed  Google Scholar 

  11. 11.

    Miller MA, Cutter AD, Yamamoto I, Ward S, Greenstein D: Clustered organization of reproductive genes in the C. elegans genome. Curr Biol. 2004, 14: 1284-1290. 10.1016/j.cub.2004.07.025.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Spellman PT, Rubin GM: Evidence for large domains of similarly expressed genes in the Drosophila genome. J Biol. 2002, 1: 5-10.1186/1475-4924-1-5.

    PubMed Central  Article  PubMed  Google Scholar 

  13. 13.

    Boutanaev AM, Kalmykova AI, Shevelyou YY, Nurminsky DI: Large clusters of co-expressed genes in the Drosophila genome. Nature. 2002, 420: 666-669. 10.1038/nature01216.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Thygesen H, Zwinderman A: Modelling the correlation between the activities of adjacent genes in drosophila. BMC Bioinformatics. 2005, 6: 10-10.1186/1471-2105-6-10.

    PubMed Central  Article  PubMed  Google Scholar 

  15. 15.

    Williams EJB, Bowles DJ: Coexpression of neighboring genes in the genome of Arabidopsis thaliana. Genome Res. 2004, 14: 1060-1067. 10.1101/gr.2131104.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  16. 16.

    Lee JM, Sonnhammer EL: Genomic gene clustering analysis of pathways in eukaryotes. Genome Res. 2003, 13: 875-882. 10.1101/gr.737703.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  17. 17.

    Wong S, Wolfe KH: Birth of a metabolic gene cluster in yeast by adaptive gene relocation. Nature Genet. 2005, 37: 777-782. 10.1038/ng1584.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Fukuoka Y, Inaoka H, Kohane IS: Inter-species differences of co-expression of neighboring genes in eukaryotic genomes. BMC Genomics. 2004, 5: art. no.-4. 10.1186/1471-2164-5-4.

    Article  Google Scholar 

  19. 19.

    Teichmann SA, Veitia RA: Genes encoding subunits of stable complexes are clustered on the yeast chromosomes: An interpretation from a dosage balance perspective. Genetics. 2004, 167: 2121-2125. 10.1534/genetics.103.024505.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  20. 20.

    Cooper DN: Human Gene Evolution. 1999, Oxford, BIOS Scientific

    Google Scholar 

  21. 21.

    Lennard A, Gorman P, Carrier M, Griffiths S, Scotney H, Sheer D, Solari R: Cloning and chromosome mapping of the human interleukin-1 receptor antagonist gene. Cytokine. 1992, 4: 83-89. 10.1016/1043-4666(92)90041-O.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Wang MH, Ronsin C, Gesnel MC, Coupey L, Skeel A, Leonard EJ, Breathnach R: Identification of the Ron gene product as the receptor for the human macrophage stimulating protein. Science. 1994, 266: 117-119.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Popovici C, Leveugle M, Birnbaum D, Coulier F: Coparalogy: Physical and functional clusterings in the human genome. Biochem Biophys Res Commun. 2001, 288: 362-370. 10.1006/bbrc.2001.5794.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Haig D: Gestational drive and the green-bearded placenta. Proc Nat Acad Sci USA. 1996, 93: 6547-6551. 10.1073/pnas.93.13.6547.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  25. 25.

    Peters LL, Barker JE: Novel inheritance of the murine severe combined anemia and thrombocytopenia (Scat) phenotype. Cell. 1993, 74: 135-142. 10.1016/0092-8674(93)90301-6.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Hurst LD: scat+ is a selfish gene analogous to Medea of Tribolium castaneum. Cell. 1993, 75: 407-408. 10.1016/0092-8674(93)90375-Z.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Fisher RA: The Genetical Theory of Natural Selection. 1930, Oxford, Clarendon Press

    Google Scholar 

  28. 28.

    Bodmer WF, Parsons PA: Linkage and recombination in evolution. Adv Genet. 1962, 11: 1-100.

    Article  Google Scholar 

  29. 29.

    Nei M: Modification of linkage intensity by natural selection. Genetics. 1967, 57: 625-641.

    PubMed Central  CAS  PubMed  Google Scholar 

  30. 30.

    Nei M: Evolutionary change in linkage intensity. Nature. 1968, 218: 1160-1161.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Marquardt T, Shirasaki R, Ghosh S, Andrews SE, Carter N, Hunter T, Pfaff SL: Coexpressed EphA receptors and Ephrin-A iigands mediate opposing actions on growth cone navigation from distinct membrane domains. Cell. 2005, 121: 127-139. 10.1016/j.cell.2005.01.020.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Graeber TG, Eisenberg D: Bioinformatic identification of potential autocrine signaling loops in cancers from gene expression profiles. Nat Genet. 2001, 29: 295-300. 10.1038/ng755.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    DLRP database. []

  34. 34.

    Keightley PD, Lercher MJ, Eyre-Walker A: Evidence for widespread degradation of gene control regions in hominid genomes. PLoS Biol. 2005, 3: 282-288. 10.1371/journal.pbio.0030042.

    CAS  Article  Google Scholar 

  35. 35.

    Keightley PD, Eyre-Walker A: Deleterious mutations and the evolution of sex. Science. 2000, 290: 331-333. 10.1126/science.290.5490.331.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Duret L, Mouchiroud D, Gouy M: HOVERGEN - a database of homologous vertebrate genes. Nucl Acids Res. 1994, 22: 2360-2365.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  37. 37.

    McLysaght A, Hokamp K, Wolfe KH: Extensive genomic duplication during early chordate evolution. Nat Genet. 2002, 31: 200-204. 10.1038/ng884.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Pebusque MJ, Coulier F, Birnbaum D, Pontarotti P: Ancient large-scale genome duplications: phylogenetic and linkage analyses shed light on chordate genome evolution. Mol Biol Evol. 1998, 15: 1145-1159.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Katsanis N, Fitzgibbon J, Fisher EMC: Paralogy Mapping: Identification of a Region in the Human MHC Triplicated onto Human Chromosomes 1 and 9 Allows the Prediction and Isolation of NovelPBXandNOTCHLoci. Genomics. 1996, 35: 101-108. 10.1006/geno.1996.0328.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Stankiewicz P, Shaw CJ, Withers M, Inoue K, Lupski JR: Serial segmental duplications during primate evolution result in complex human genome architecture. Genome Res. 2004, 14: 2209-2220. 10.1101/gr.2746604.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  41. 41.

    Veitia RA: Exploring the etiology of haploinsufficiency. Bioessays. 2002, 24: 175-184. 10.1002/bies.10023.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Veitia RA: Gene dosage balance: deletions, duplications and dominance. Trends Genet. 2005, 21: 33-35. 10.1016/j.tig.2004.11.002.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Papp B, Pal C, Hurst LD: Dosage sensitivity and the evolution of gene families in yeast. Nature. 2003, 424: 194-197. 10.1038/nature01771.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Bourque G, Zdobnov EM, Bork P, Pevzner PA, Tesler G: Comparative architectures of mammalian and chicken genomes reveal highly variable rates of genomic rearrangements across different lineages. Genome Res. 2005, 15: 98-110. 10.1101/gr.3002305.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  45. 45.

    Zhao SY, Shetty J, Hou LH, Delcher A, Zhu BL, Osoegawa K, de Jong P, Nierman WC, Strausberg RL, Fraser CM: Human, mouse, and rat genome large-scale rearrangements: Stability versus speciation. Genome Res. 2004, 14: 1851-1860. 10.1101/gr.2663304.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  46. 46.

    O'Brien SJ, Menotti-Raymond M, Murphy WJ, Nash WG, Wienberg J, Stanyon R, Copeland NG, Jenkins NA, Womack JE, Graves JAM: The promise of comparative genomics in mammals. Science. 1999, 286: 458-462. 10.1126/science.286.5439.458.

    Article  PubMed  Google Scholar 

  47. 47.

    Fares MA, Byrne KP, Wolfe KH: Rate Asymmetry After Genome Duplication Causes Substantial Long Branch Attraction Artifacts in the Phylogeny of Saccharomyces Species. Mol Biol Evol. 2005

    Google Scholar 

  48. 48.

    Popovici C, Roubin R, Coulier F, Birnbaum D: An evolutionary history of the FGF superfamily. Bioessays. 2005, 27: 849-857. 10.1002/bies.20261.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Unigene. []

  50. 50.

    Entrez. []

  51. 51.

    Steinkasserer A, Spurr NK, Cox S, Jeggo P, Sim RB: The human IL-1 receptor antagonist gene (IL1RN) maps to chromosome 2q14-q21, in the region of the IL-1 alpha and IL-1 beta loci. Genomics. 1992, 13: 654-657. 10.1016/0888-7543(92)90137-H.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Dale M, Nicklin MJ: Interleukin-1 receptor cluster: gene organization of IL1R2, IL1R1, IL1RL2 (IL-1Rrp2), IL1RL1 (T1/ST2), and IL18R1 (IL-1Rrp) on human chromosome 2q. Genomics. 1999, 57: 177-179. 10.1006/geno.1999.5767.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Wolfe K: Evolutionary genomics: yeasts accelerate beyond BLAST. Curr Biol. 2004, 14: R392-394. 10.1016/j.cub.2004.05.015.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    MGD Human Mouse orthologs. []

  55. 55.

    Ensembl MartView. []

  56. 56.

    Davison AC, Hinkley DV: Bootstrap methods and their application. 1997, Cambridge, United Kingdom, Cambridge University Press

    Google Scholar 

  57. 57.

    North BV, Curtis D, Sham PC: A note on the calculation of empirical P values from Monte Carlo procedures. Am J Hum Genet. 2002, 71: 439-441. 10.1086/341527.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  58. 58.

    Jin M, Chen Y, He S, Ryan SJ, Hinton DR: Hepatocyte growth factor and its role in the pathogenesis of retinal detachment. Invest Ophthalmol Vis Sci. 2004, 45: 323-329. 10.1167/iovs.03-0355.

    Article  PubMed  Google Scholar 

  59. 59.

    Schmidt L, Duh FM, Chen F, Kishida T, Glenn G, Choyke P, Scherer SW, Zhuang Z, Lubensky I, Dean M, Allikmets R, Chidambaram A, Bergerheim UR, Feltis JT, Casadevall C, Zamarron A, Bernues M, Richard S, Lips CJ, Walther MM, Tsui LC, Geil L, Orcutt ML, Stackhouse T, Zbar B: Germline and somatic mutations in the tyrosine kinase domain of the MET proto-oncogene in papillary renal carcinomas. Nat Genet. 1997, 16: 68-73. 10.1038/ng0597-68.

    CAS  Article  PubMed  Google Scholar 

  60. 60.

    Bezerra JA, Carrick TL, Degen JL, Witte D, Degen SJ: Biological effects of targeted inactivation of hepatocyte growth factor-like protein in mice. J Clin Invest. 1998, 101: 1175-1183.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  61. 61.

    Muraoka RS, Sun WY, Colbert MC, Waltz SE, Witte DP, Degen JL, Friezner Degen SJ: The Ron/STK receptor tyrosine kinase is essential for peri-implantation development in the mouse. J Clin Invest. 1999, 103: 1277-1285.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  62. 62.

    Joutel A, Corpechot C, Ducros A, Vahedi K, Chabriat H, Mouton P, Alamowitch S, Domenga V, Cecillion M, Marechal E, Maciazek J, Vayssiere C, Cruaud C, Cabanis EA, Ruchoux MM, Weissenbach J, Bach JF, Bousser MG, Tournier-Lasserve E: Notch3 mutations in CADASIL, a hereditary adult-onset condition causing stroke and dementia. Nature. 1996, 383: 707-710. 10.1038/383707a0.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Miller DL, Ortega S, Bashayan O, Basch R, Basilico C: Compensation by fibroblast growth factor 1 (FGF1) does not account for the mild phenotypic defects observed in FGF2 null mice. Mol Cell Biol. 2000, 20: 2260-2268. 10.1128/MCB.20.6.2260-2268.2000.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  64. 64.

    Kawaguchi H, Nakamura K, Tabata Y, Ikada Y, Aoyama I, Anzai J, Nakamura T, Hiyama Y, Tamura M: Acceleration of Fracture Healing in Nonhuman Primates by Fibroblast Growth Factor-2. J Clin Endocrinol Metab. 2001, 86: 875-880. 10.1210/jc.86.2.875.

    CAS  Article  PubMed  Google Scholar 

  65. 65.

    Weinstein M, Xu X, Ohyama K, Deng CX: FGFR-3 and FGFR-4 function cooperatively to direct alveogenesis in the murine lung. Development. 1998, 125: 3615-3623.

    CAS  PubMed  Google Scholar 

  66. 66.

    Zammit C, Coope R, Gomm JJ, Shousha S, Johnston CL, Coombes RC: Fibroblast growth factor 8 is expressed at higher levels in lactating human breast and in breast cancer. Br J Cancer. 2002, 86: 1097-1103. 10.1038/sj.bjc.6600213.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  67. 67.

    Muenke M, Schell U, Hehr A, Robin NH, Losken HW, Schinzel A, Pulleyn LJ, Rutland P, Reardon W, Malcolm S: A common mutation in the fibroblast growth factor receptor 1 gene in Pfeiffer syndrome. Nat Genet. 1994, 8: 269-274. 10.1038/ng1194-269.

    CAS  Article  PubMed  Google Scholar 

  68. 68.

    Xu J, Liu Z, Ornitz DM: Temporal and spatial gradients of Fgf8 and Fgf17 regulate proliferation and differentiation of midline cerebellar structures. Development. 2000, 127: 1833-1843.

    CAS  PubMed  Google Scholar 

  69. 69.

    Hu MC, Qiu WR, Wang YP, Hill D, Ring BD, Scully S, Bolon B, DeRose M, Luethy R, Simonet WS, Arakawa T, Danilenko DM: FGF-18, a novel member of the fibroblast growth factor family, stimulates hepatic and intestinal proliferation. Mol Cell Biol. 1998, 18: 6063-6074.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  70. 70.

    OMIM. []

Download references


We thank David N. Cooper for discussion. We thank one referee for helpful comments.

Author information



Corresponding author

Correspondence to Laurence D Hurst.

Additional information

Authors' contributions

Both authors contributed to analysis, design and write up.

Electronic supplementary material

Additional file 1: Supplement 1. The data set of ligand receptor pairs after purging of tandem duplicates by Blast score and by purging of ligands within 1 Mb of each other and receptors within 1 MB of each other. Unigene_id is the Unigene reference number. Entrez_id is the Entrez/LocusLink number. Chr is the chromosome on which the gene is found. Cytog is the cytogenetic location of the gene. Genbank_chr_file is the relevant accession number of the whole chromosome GenBank file. From, to, midpos are the ends of the gene and the middle position of the gene in base pairs, these referring to the locations in the NC files. Number of tissues is a count of the number of entries in cDNA source entry in the Unigene page. Placental: 1 means placentally expressed, 0 means no evidence for placental expression. GI: GI number. RefSeq: the ReqSeq Genbank number (XLS 35 KB)

Additional file 2: Supplement 2: The data set of ligand receptor pairs after purging of tandem duplicates by Blast scores alone. Annotation as above. (XLS 48 KB)

Additional file 3: Supplement 3: The set of ligand-receptor pairs on the same chromosome in the data set given in supplement 2 (Blast score purged alone). Annotation as for figure 1 in the table. Note: In the above three data sets the final two ligand-receptor pairs are those nominated by Cooper. (XLS 2 KB)

Supplement 4:

Additional file 4: The results and analysis of randomizations, excluding randomizations controlling for breadth. The results sheet is split in two. The top half ("all genes") refers to analysis of the augmented ligand receptor data set, i.e. the one containing the two ligand-receptor pairs nominated by Cooper. The lower half excludes these two. In each, as detailed in the methods, two filtering methods were used to eliminate tandem duplicates: Blast alone (see Additional file 2) and blast plus close proximity of ligands, close proximity of receptors (see Additional file 1). Within each of these subdivisions two data sets were employed to determine where, in randomizations, a ligand and/or receptor might locate. In one case (lig/rec genes only), the positions of ligands and receptors were interchanged with ligands/receptors derived from the test set. In the second method, possible locations are anywhere in the genome where a gene, of any variety, resides ("all genes in genome"). Randomizations for each of the analysis in turn was performed either permitting a ligand or receptor to locate anywhere in the genome ("all genes") or anywhere on the same chromosome ("within chr"). Analysis was further performed on those ligand receptor sets in which at least one of the genes in the set was placentally expressed ("placental"). N is the number of genes in the test set, "sameChr", refers to the number of incidences of ligand receptor pairs on the same chromosome, "sameChrRand"", is the corresponding number for the randomizations. P_chr is the P value for the comparisons. <d> is the mean distance between ligands and receptors on the same chromosome (in kb) and <dRand> the corresponding number in the randomizations. P_d is the relevant P value. "Count[Nkb]", "countRand, P refer to the number of incidences in the real data set of ligand receptor pairs less than N kb apart, the number in the randomizations and the P value. This statistic is only relevant for the cases where randomization is done within chromosomes. The between chromosome randomization by its nature confounds the effect of the excess of pairs on the same chromosome as well as their proximity. (XLS 21 KB)

Additional file 5: Supplement 5: This supplement details the analysis (without control for breadth of expression) in the mouse data sets. (XLS 18 KB)

Additional file 6: Supplement 6: This supplement details the analysis with control for breadth of expression in the mouse. (XLS 2 KB)

Additional file 7: Supplement 7: This supplement details the evidence for co-duplication for paralogous ligand-receptor pairs (DOC 32 KB)

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Hurst, L.D., Lercher, M.J. Unusual linkage patterns of ligands and their cognate receptors indicate a novel reason for non-random gene order in the human genome. BMC Evol Biol 5, 62 (2005).

Download citation


  • Hepatocyte Growth Factor
  • Tandem Duplication
  • Cognate Receptor
  • Paralogous Pair
  • Balance Hypothesis