Correlation of microsynteny conservation and disease gene distribution in mammalian genomes
© Lovell et al; licensee BioMed Central Ltd. 2009
Received: 19 December 2008
Accepted: 12 November 2009
Published: 12 November 2009
With the completion of the whole genome sequence for many organisms, investigations into genomic structure have revealed that gene distribution is variable, and that genes with similar function or expression are located within clusters. This clustering suggests that there are evolutionary constraints that determine genome architecture. However, as most of the evidence for constraints on genome evolution comes from studies on yeast, it is unclear how much of this prior work can be extrapolated to mammalian genomes. Therefore, in this work we wished to examine the constraints on regions of the mammalian genome containing conserved gene clusters.
We first identified regions of the mouse genome with microsynteny conservation by comparing gene arrangement in the mouse genome to the human, rat, and dog genomes. We then asked if any particular gene types were found preferentially in conserved regions. We found a significant correlation between conserved microsynteny and the density of mouse orthologs of human disease genes, suggesting that disease genes are clustered in genomic regions of increased microsynteny conservation.
The correlation between microsynteny conservation and disease gene locations indicates that regions of the mouse genome with microsynteny conservation may contain undiscovered human disease genes. This study not only demonstrates that gene function constrains mammalian genome organization, but also identifies regions of the mouse genome that can be experimentally examined to produce mouse models of human disease.
The availability of several mammalian genome sequences has enabled comparative genomic studies to identify regions of conserved linkage among different organisms (reviewed in [1–4]). These studies have been used to predict the genome architecture of the common mammalian ancestor [5, 6], as well as to assess recombination [7, 8] and genome evolution . Additionally, regulatory elements have been identified by the characterization of conserved regions of non-coding DNA [10–12]. Non-coding functional RNAs and microRNA targets have also been identified through comparative genomic approaches [3, 13, 14]. Gene function may be inferred from conserved proteins in other species. Likewise, comparative genomics among mammalian species is useful for predicting the functional consequences of mutations in human disease loci [15–18]. Additionally, mapping genes responsible for quantitative traits in rodents allows the prediction of locations of human quantitative traits underlying disease, based on conserved genomic structure between rodents and humans [19–22].
Although it is becoming increasingly apparent that genomes display a large degree of structural plasticity, there are nevertheless significant evolutionary constraints on genome structure. Previous studies have provided evidence for functional constraints on genome organization in prokaryotic and eukaryotic genomes . Studies in yeast demonstrate that essential genes are found in genomic clusters . The clustering of essential genes in yeast is likely driven by selection for reduced noise in gene expression levels, as essential gene clusters are localized in regions of open chromatin . Additionally, in the nematode C. elegans, essential genes are located in clusters in regions with low recombination .
Clustering of genes with similar functions has also been observed in mammalian genomes. In the human genome, genes that are in the same pathway are in closer proximity than would be expected by chance . Similarly, in the mouse genome, genes with common GO annotations are found in clusters . This is not due to tandem duplications, as most genes in the same pathway that are adjacent in the genome do not arise from duplication events . It is possible that functionally related genes are located in clusters to facilitate coordinated transcription, as many genes in clusters are co-expressed .
Many of the previous studies to detect gene clustering were based on bioinformatic analysis of genome annotation. However, there is also support for functional constraints on mammalian genome organization from experimental data. Analysis of saturation levels of mouse mutagenesis screens for lethal phenotypes directed at specific genomic regions demonstrated that mouse essential genes are disproportionately found in regions of conserved microsynteny , at least for the small number of genomic regions evaluated. To build on this prior work, we assessed microsynteny conservation on all mouse autosomes. By examining gene content in conserved and divergent genomic regions we found a significant correlation between microsynteny conservation and the density of mouse genes that are orthologous to human disease genes. As the mouse is widely used to model human disease [1, 30], the identification of this correlation will facilitate the creation of new mouse models of human disease by identifying regions of the mouse genome that contain a high density of disease gene orthologs.
Microsynteny conservation of mouse autosomes
We evaluated the level of microsynteny conservation between the mouse genome and those of human, dog and rat. First, we obtained all protein-coding genes and their genomic locations on all mouse autosomes as annotated in the Ensembl mouse genome browser [3, 31] (release 50). We also obtained protein-coding genes from the human [32, 33], rat , and dog  genomes. These were chosen because they had a sufficient level of assembly and annotation to allow comparison. The use of the dog genome as an outgroup to human, rat, and mouse improves the stringency of the study . To identify orthologs of mouse genes in the other genomes, along with their genomic locations, Ensembl BioMart homology filters were used to compile a list of orthologous genes.
Although the four mammalian genomes chosen are those with the best available annotation, the degree and quality of annotation may vary somewhat between species. In order to control for this we took additional steps to find the human, rat, and dog orthologs of mouse genes. Protein sequences of all mouse genes that did not have an annotated ortholog in another species were searched using BLAT  against the other genomes to identify orthologous sequences in the other genomes. To allow a moderately strict search with a limited number of false positives, all hits with E-value < 10-5 were retrieved. The genomic location of the best BLAT match in the other genome was used for the evaluation of microsynteny conservation. We searched a total of 6173 genes in at least one other species, finding an ortholog for 1210 of the genes in other genomes. A sensitivity analysis demonstrates that the number of genes retrieved from the BLAT searches is relatively insensitive to the choice of E-value cut-off, as changing the cut-off point from 10-3 to 10-7 results in 1266 - 1169 ortholog annotations respectively. Therefore, utilizing alternate E-value cut-off points in this range would have changed the annotation of only 0.48% of the total mouse genes analyzed in our study.
Genes were defined as having conserved microsynteny if their orthologs had the same two orthologous neighboring genes in all four species examined. Each mouse gene was queried to determine whether it met these criteria for conserved microsynteny. We then assessed the level of conservation for segments of the mouse genome, determining the percentage of conserved genes in each segment. We examined 20 Mb regions of the mouse genome, as genomic regions of this size have been analyzed experimentally through region-specific mutagenesis [37–40]. Thus, the identification of additional genomic regions with conserved microsynteny will be useful for further experimental functional genomic annotation.
We found that the conservation of microsynteny varied throughout the mouse genome. The average percentage of genes with conserved microsynteny in a 20 Mb interval was 38.29%, with a standard deviation of 11.81%. The Shapiro-Wilk test was applied to non-overlapping 20 Mb windows to determine whether the distribution of gene microsynteny was normal. A P-value of 0.16 indicated that the null hypothesis (a normal distribution) should not be rejected. However, the use of non-overlapping 20 Mb windows restricted the resolution of the study. For example, chromosome 19 has only two observations from non-overlapping windows. To improve the resolution of our study, we next examined 20 Mb intervals staggered by 5 Mb. This sliding window analysis allowed more observations on each chromosome.
Variation of microsynteny and disease gene density on all mouse autosomes.
Total for all autosomes
Comparison with sequence-based synteny blocks
Comparison of regions with Z>1 microsynteny to sequence-based synteny blocks.
Mouse Micro-synteny Interval
Dog Sequence Synteny Blocks
Rat Sequence Synteny Blocks
Human Sequence Synteny Blocks
Chr1 60 - 90 Mb
60 - 80 Mb = Chr37: 14.6 - 32.7 Mb,
80 - 90 Mb = Chr25: 40.6 - 48.1 Mb
60 - 85 Mb = Chr9: 58.5 - 84.4 Mb,
87 Mb - 90 Mb = Chr9: 84.2 - 87.0 Mb
60 - 90 Mb = Chr2: 203.6 - 234.6
Chr1 120 - 145 Mb
120 - 132 = Chr19: 31.4 - 43.3 Mb
132 - 133 Mb = Chr7: 8.7 - 9.3 Mb,
133 - 136 Mb = Chr38: 3.0 - 5.6 Mb,
136 - 140 Mb = Chr7: 3.0 - 7.1 Mb
120 - 145 Mb = Chr13: 30.9 - 57.8 Mb
120 - 132 Mb = Chr2: 125.6 - 138.4 Mb,
132 - 145 Mb = Chr1: 193.0 - 207.7 Mb
Chr1 155 - 175 Mb
155 - 168 Mb = Chr7: 19.4 - 33.9 Mb,
169 - 175 Mb = Chr38: 20.0 - 25.5 Mb
155 - 175 Mb = Chr13: 68.2 - 89.3 Mb
155 - 175 Mb = Chr1: 159.2 - 182.9 Mb
Chr2 45 - 85 Mb
45 - 53 Mb = Chr19: 52.6 - 56.7 Mb,
54 - 84 Mb = Chr36: 3.8 - 32.4 Mb,
84 - 85 Mb = Chr18: 41.4 - 41.9 Mb
45 - 85 Mb = Chr3: 25.6 - 68.4 Mb
45 - 84 Mb = Chr2: 145.2 - 188.3 Mb,
84 - 85 Mb = Chr11: 56.9 - 57.5 Mb
Chr2 115 - 140 Mb
115 - 127 = Chr30: 7.1 - 19.7 Mb,
127 - 129 Mb = Chr17: 37.5 - 40.1 Mb,
129 - 140 Mb = Chr24: 22.0 - 11.7 Mb
115 - 140 Mb = Chr3: 101.6 - 128.3 Mb
115 - 126 Mb = Chr15: 36.8 - 51.2 Mb,
126 - 127 Mb = Chr2: 95.6 - 97.0 Mb,
127 - 129 Mb = Chr2: 110.8 - 113.5 Mb,
129 - 140 Mb = Chr20: 1.4 - 13.8 Mb
Chr3 15 - 35 Mb
15 - 16 Mb = Chr24: 22.4 - 22.5 Mb,
16 - 19 Mb = Chr29: 16.2 - 18.7 Mb,
19 - 20 Mb = Chr23: 46.7 - 47.0 Mb
21 - 30 Mb = Chr34: 43.5 - 37.0 Mb
31 - 35 Mb = Chr34: 15.5 - 18.8 Mb
15 - 35 Mb = Chr2: 87.0 - 122.2 Mb
15 - 16 Mb = Chr20: 1.5 - 1.6 Mb,
16 - 19 Mb = Chr8: 64.0 - 67.0 Mb,
19 - 20 Mb = Chr3: 148.4 - 148.8 Mb,
20 - 35 Mb = Chr3: 176.7 - 182.5 Mb
Chr3 45 - 65 Mb
45 - 51 Mb = Chr19: 5.8 - 11.6 Mb,
52 - 55 Mb = Chr25: 3.5 - 8.0 Mb,
57 - 65 Mb = Chr23: 47.2 - 53.1 Mb
45 - 65 Mb = Chr2: 132.8 - 155.1 Mb
45 - 52 Mb = Chr4: 134.0 - 140.8 Mb,
52 - 55 Mb = Chr13: 36.0 - 41.1 Mb,
57 - 65 Mb = Chr3: 149.0 - 156.2 Mb
Chr3 55 - 75 Mb
55 - 57 Mb = Chr25: 7.6 - 8.0 Mb,
57 - 67 Mb = Chr23: 47.2 - 55.0 Mb,
67 - 75 Mb = Chr34: 28.0 - 35.3 Mb
55 - 75 Mb = Chr2: 144.4 - 166.1 Mb
55 - 57 Mb = Chr13: 36.0 - 36.3 Mb,
57 - 75 Mb = Chr3: 149.0 - 167.1 Mb
Chr4 30 - 65 Mb
30 - 35 Mb = Chr12: 52.6 - 49.6 Mb,
35 - 65 Mb = Chr11: 48.5 - 73.2 Mb
30 - 65 Mb = Chr5: 48.2 - 82.0 Mb
30 - 34 Mb = Chr6: 87.8 - 93.9 Mb,
34 - 45 Mb = Chr9: 27.3 - 38.4 Mb,
46 - 65 Mb = Chr9: 100.1 - 119.1 Mb
Chr4 105 - 125 Mb
105 - 107 Mb = Chr5: 57.1 - 58.9 Mb,
107 - 118 Mb = Chr15: 11.6 - 19.9 Mb,
118 - 125 Mb = Chr15: 3.2 - 7.9 Mb
105 - 125 Mb = Chr5: 127.3 - 144.8 Mb
105 - 125 Mb = Chr1: 37.2 - 55.5 Mb
Chr6 0 - 20 Mb
0 - 13 Mb = Chr14: 21.6 - 30.2 Mb
13 - 20 Mb - Chr14: 55.2 - 55.9 Mb
0 - 20 Mb = Chr4: 28.1 - 44.7 Mb
0 - 20 Mb = Chr7: 92.7 - 117.8 Mb
Chr9 40 - 70 Mb
40 - 54 Mb = Chr5: 13.7 - 27.7 Mb,
54 - 70 Mb = Chr30: 19.8 - 27.0 Mb
40 - 70 Mb = Chr8: 43.2 - 74.6 Mb
40 - 54 Mb = Chr11: 107.4 - 123.4 Mb,
54 - 70 Mb = Chr15: 59.4 - 78.3 Mb
Chr11 60 - 90 Mb
60 - 72 Mb = Chr5: 33.1 - 44.6 Mb,
72 - 90 Mb = Chr9: 33.3 - 51.0 Mb
60 - 90 Mb = Chr10: 46.4 - 78.6 Mb
60 - 76 Mb = Chr17: 1.1 - 17.7 Mb,
76 - 90 Mb = Chr17: 28.8 - 53.4 Mb
Chr11 100 - 120 Mb
100 - 120 Mb = Chr9: 3.6 - 24.5 Mb
100 - 120 Mb = Chr10: 89.0 - 109.7 Mb
100 - 104 Mb = Chr17: 39.6 - 45.2 Mb,
104 - 105 Mb = Chr7: 128.0 - 128.1 Mb,
106 - 120 Mb = Chr17: 60.5 - 79.2 Mb
Chr12 55 - 105 Mb
55 - 72 Mb = Chr8: 15.7 - 31.0 Mb,
72 - 102 Mb = Chr8: 36.6 - 65.4 Mb,
102 - 104 Mb = Chr8: 3.9 - 5.1 Mb,
104 - 105 Mb = Chr8: 65.4 - 66.4 Mb
55 - 105 Mb = Chr6: 74.3 - 127.9 Mb
55 - 72 Mb = Chr14: 34.2 - 52.3 Mb,
72 - 105 Mb = Chr14: 58.7 - 94.8 Mb
Chr15 65 - 104 Mb
65 - 76 Mb = Chr13: 31.6 - 41.2 Mb,
77 - 89 Mb = Chr10: 19.5 - 31.6 Mb,
89 - 104 Mb = Chr27: 3.8 - 18.6 Mb
65 - 104 Mb = Chr7: 103.1 - 142.6 Mb
65 - 76 Mb = Chr8: 132.9 - 146.2 Mb
77 - 89 Mb = Chr22: 36.0 - 51.2 Mb
89 - 104 Mb = Chr12: 33.5 - 54.9 Mb
Chr17 65 - 95 Mb
65 - 66 Mb = Chr3: 4.6 - 5.2 Mb,
66 - 72 Mb = Chr7: 72.2 - 78.8 Mb,
72 - 84 Mb = Chr17: 25.8 - 37.3 Mb,
84 - 90 Mb = Chr10: 48.4 - 54.8 Mb,
91 - 95 Mb = Chr7: 70.7 - 72.3 Mb
65 - 72 Mb = Chr9: 103.9 - 110.8 Mb,
72 - 77 Mb = Chr6: 17.9 - 24.1 Mb,
77 - 93 Mb = Chr6: 0.9 - 17.8 Mb,
93 - 95 Mb = Chr9: 110.7 - 112.5 Mb
65 - 66 Mb = Chr5: 109.2 - 110.0 Mb,
66 - 71 Mb = Chr18: 2.5 - 9.8 Mb,
72 - 93 Mb = Chr2: 29.0 - 51.2 Mb,
94 - 95 Mb = Chr18: 0.9 - 2.5 Mb
Chr19 30 - 65 Mb
30 - 34 Mb = Chr26: 38.4 - 41.9 Mb,
34 - 65 Mb = Chr28: 7.3 - 32.4 Mb
30 - 65 Mb = Chr1: 233.7 - 267.8 Mb
30 - 32 Mb = Chr10: 51.9 - 54.5 Mb,
32 - 65 Mb = Chr10: 89.2 - 121.2 Mb
Conserved genes are located next to other conserved genes
Although there is variation in the density of genes with conserved microsynteny across the genome, it is possible that this variation merely represents random variation within a normal distribution. To determine whether the genomic arrangement of genes with conserved microsynteny is random, we calculated the likelihood that a gene with syntenic conservation is found next to another gene with syntenic conservation. We then compared this to the frequency of conserved-synteny neighbors in a set of 10,000 randomized genomes. In each of the randomized genomes the number and position of genes is maintained, but the annotation of conservation is randomly shuffled.
Distribution of mouse orthologs of human disease-related genes
As we found that conserved genes were more likely to have conserved genes as neighbors, we investigated whether any other groups of genes were found preferentially in regions of the genome with conserved microsynteny. One group of genes that are of interest is disease genes, as they are highly relevant to human health. We therefore performed a genome-wide analysis of the mouse orthologs of human disease-related genes to assess whether they were found at a greater density in conserved regions of the mouse genome. We identified human genes with a disease-associated mutant allele from the OMIM Morbid Map database, and cross-referenced them to the mouse genome using Ensembl BioMart to identify mouse orthologs. Using the genomic locations of the mouse orthologs of human disease genes from our study on microsyntenty conservation, we determined the proportion of human disease gene orthologs in each 20 Mb sliding window of the mouse genome. The mean percentage of disease-related gene orthologs as compared to the number of total genes in a 20 Mb interval is 7.68%, with a standard deviation of 2.72%. We found variation in the distribution of disease gene orthologs in the mouse genome, with the highest percentage found on chromosome 18, and the lowest on chromosome 7 (Table 1). However, on all other autosomes, 7 - 8% of the genes are orthologs of human disease genes.
Disease gene orthlogs are located next to other disease gene orthologs
As we had found that the distribution of genes with conserved microsynteny is non-random, we examined whether that was also true for genes with disease-related orthologs. Using a similar approach, we calculated the number of mouse orthologs of human disease genes with at least one disease gene ortholog as a neighbor. As a control, we randomized which genes were annotated as disease orthologs, keeping the same total number of disease gene orthologs. From 10,000 random trails we found that the mouse orthologs of human disease genes were significantly more likely to have other disease gene orthologs as neighbors (Figure 2B, P < 0.07). This finding demonstrates that the distribution of orthologs of human disease genes in the mouse genome is not random.
Correlation between microsynteny conservation and disease gene distribution
The density of disease gene orthologs for each genomic region is shown in Figure 1 (blue lines). Z-scores are displayed to allow direct comparison between microsynteny conservation and disease orthologs. The additional calculation of Z-scores does not change the overall correlation. There is a significant correlation (P = 0.05) between microsynteny conservation and the density of disease gene orthologs for 12 of the 19 mouse autosomes. Thus, genomic regions with a high percentage of genes with conserved microsynteny also have a high percentage of disease gene orthologs. The chromosome with the best correlation between conserved microsynteny and density of disease gene orthologs is mouse chromosome 13, while the chromosome with the worst correlation is mouse chromosome 10.
To demonstrate that this correlation was not an artifact of our analysis, we randomized the annotation of disease genes. We assigned alternate genes as orthologs of human disease genes, keeping the total number of disease genes per chromosome the same as the first analysis. We then recalculated the percentage of alternate disease genes as compared to total genes in each sliding window throughout the genome, and the average and standard deviation for each sliding window. We plotted the Z-scores for each window containing these alternate disease orthologs, and compared them to the Z-scores for microsynteny conservation (Figure 1 red lines). When the chromosomal positions of orthologs of disease genes are changed to random locations, the correlation with microsynteny disappears (Pearson's R = 0.02, P < 0.58). As an additional control, we also randomized disease genes while retaining the same number of observed disease gene pairs for each chromosome. Again, we found no correlation (Pearson's R = 0.004, P < 0.93). Should the correlation between observed disease gene ortholog distribution and microsynteny conservation be an artifact of our methodology, we would also expect the randomized annotations to be correlated. This is not the case, demonstrating that the link between microsynteny correlation and density of disease gene orthologs does not arise from an artifact of the methodology.
Robustness to changes in window size
Robustness of correlation to variations in window size.
R = 0.40,
p < 4 × 10-4
R = 0.293,
p < 1 × 10-10
R = 0.256,
p < 1 × 10-10
R = 0.157,
p < 1 × 10-10
R = 0.154,
p < 1 × 10-10
R = 0.02,
p = 0.58
R = 0.018,
p = 0.59
R = 0.058,
p = 0.11
R = 0.046,
p = 0.45
R = 0.035,
p = 0.01
Actual annotation, empty windows removed
R = 0.11
p < 1 × 10-10
R = 0.11,
p < 1 × 10-10
Randomized annotation, empty windows removed
R = 0.012,
p = 0.40
R = 0.020,
p = 0.83
Number of windows in mouse genome at varying window sizes.
No annotated genes
No conserved genes
No disease genes
No disease or conserved
We have examined the relationship between microsynteny conservation and the density of orthologs of human disease genes in the mouse genome. We found a correlation between regions of conserved microsynteny and the location of mouse orthologs of human disease genes, which is consistent for variations in the window size used in our analyses. The correlation we observe suggests that regions of the mouse genome with a high density of disease gene orthologs undergo less rearrangement than regions of the genome with fewer disease gene orthologs. Genes associated with human disease are often orthologous to essential genes in other organisms . Previous studies from both mammals  and other eukaryotes [24, 26] have shown that essential genes are located in highly conserved genomic regions. Thus, disease-related genes, which perform essential functions, are more likely to be found in conserved regions of the genome.
Several studies have found that at the sequence level, human disease genes are more conserved than non-disease genes [26, 43, 44]. The sequence conservation of human disease genes with essential C. elegans orthologs is higher than those disease genes whose orthologs are not lethal when mutated . Interestingly, genes with high polymorphism among humans, but no divergence between humans and chimpanzees, are highly associated with Mendelian disease . Similarly, human disease genes with weak negative selection, where mutant alleles persist in the population, are more likely to cause diseases with Mendelian inheritance . Mendelian disease genes are more constrained evolutionarily than disease genes with non-Mendelian inheritance patterns . Together, these observations support our finding that the mouse orthologs of human disease genes are preferentially found in genomic regions with high microsynteny conservation.
Recombination may be mutagenic due to the possibility of unequal crossing-over. Thus faulty recombination events in regions with essential genes are likely to be deleterious to the survival of the organism and may thus be selected against during mammalian evolution. Studies of the human genome support this link between low recombination and essential genes. Regions of the human genome with high linkage disequilibrium, and thus low recombination, are enriched for genes associated with essential cellular functions such as response to DNA damage, cell cycle progression, or DNA and RNA metabolism . Genes that show variation in populations, such as immune response genes, are often found in regions with low linkage disequilibrium, suggesting that recombination in these regions is not deleterious to the organism . Likewise, human genes found in mutation cold spots tend to be genes involved in essential cellular processes, while those in mutation hot spots include immune response genes . These findings extend to non-coding sequences as well, as human genomic regions that are highly conserved with the pufferfish have been found to contain enhancers for developmental genes .
The correlation between disease gene density and microsynteny conservation, although significant, is not perfect. Discrepancies may come from several sources. For example, annotation of human disease genes is incomplete. Many housekeeping genes, which are likely to be essential for mammalian development, are not annotated as human disease genes, probably because mutations in these genes are lethal early in development, and thus humans with mutations are not viable . The genomic region between 55 - 75 Mb on mouse chromosome 3 shows high conservation but a low density of disease gene orthologs. However, the genes Wwtr1 and Shox2 are located in this genomic region. A mouse knock-out of Wwtr1 displays a phenotype resembling human polycystic kidney disease , and the mouse knock-out of Shox2 is lethal with cleft palate , strongly suggesting that these genes are linked to human disease, although neither is annotated as a disease gene in OMIM.
Likewise, many genes that are annotated as human disease genes may not be strictly essential for survival, and thus these genes are not expected to have conserved microsynteny. The genomic region between 85 - 105 Mb on mouse chromosome 12 has a high density of disease gene orthologs but low conservation. Mutations in the human gene SERPINA10, whose ortholog is located in this region, are associated with susceptibility to deep vein thromboses . Although SERPINA10 is annotated in OMIM as a disease gene, it is unlikely that inherited mutations in SERPINA10 present a challenge to survival of the individual, suggesting that SERPINA10 does not represent an essential gene. Finally, many diseases, especially cancers, are caused by translocation events that produce chimeric proteins. While a genomic region may have a great density of disease loci due to translocations, these regions would not show microsynteny conservation, as they are high in rearrangements.
Discrepancies between microsynteny conservation and the density of disease-related gene orthologs may also arise because other factors contribute to selective pressure on genome evolution. For example, previous studies have suggested that mammalian genes are clustered into groups based on co-expression [54, 55]. It is proposed that gene expression is therefore an evolutionary constraint on genome organization, although the effect is weak as gene clusters are found only slightly more often than by chance . There is also evidence that many over-lapping gene pairs exist in mammalian genomes, and that these gene pairs are conserved in multiple species, probably because recombination or mutation in these regions of the genome would cause deleterious mutations in both genes . Alternative mechanisms for the presence of essential genes constraining genome structure have also been proposed .
We have demonstrated the non-random distribution both of genes with conserved microsynteny and genes with disease orthologs. This observation suggests that there are constraints on genome organization in the mouse. Moreover, we have demonstrated that there is a correlation between mammalian genome architecture and gene function. It is likely that this correlation arises from gene function constraining genome organization, resulting in essential disease genes being located in regions of the mammalian genome with high conservation. The identification of a correlation between microsynteny conservation and density of disease gene orthologs suggests that additional experimental analysis of mouse genes in highly conserved genomic regions will produce new mammalian disease models by creating mutations in the orthologs of human disease genes.
All protein coding genes on mouse autosomes were retrieved using a BioMart search on the Ensembl Genome Browser (http://www.ensembl.org/index.html, release 50, NCBIM37 dataset). Orthologs of these genes in the human, rat, and dog genomes, along with their genomic positions, were retrieved using Ensembl BioMart homology filters. Any mouse genes that did not have an ortholog in another genome identified in BioMart were then subjected to a BLAT search http://genome.ucsc.edu/cgi-bin/hgBlat?command=start on the other genome using the mouse protein sequence retrieved with the BioMart sequences filter. The genomic location of the best match to the mouse protein sequence from the second genome was retrieved. Locally developed software was used to identify genes with conserved microsynteny. Each mouse gene is analyzed in turn. For each mouse gene m, its neighbors in the mouse genome (m-1 and m+1) are identified. The orientation of m-1, m and m+1 is determined from the start and stop positions. In the dog, human and rat genomes, the homologue of m is identified from both Ensembl BioMart homology filters and BLAT searches. We term these homologues d, h and r. For each of these genes their neighbors in their respective genomes (d-1, d+1, h-1, h+1, r-1 and r+1) are identified, and their orientation determined. The gene m is defined as having conserved microsynteny only if all of the following are true: (i) m, d, h, and r comprise a homologous set of genes (ii) m-1, d-1, h-1 and r-1 comprise a homologous set (iii) m+1, d+1, h+1 and r+1 comprise a homologous set (iv) the m, d, h, and r have the same orientation (v) m-1, d-1, h-1 and r-1 have the same orientation and (iv) m+1, d+1, h+1 and r+1 have the same orientation. Results are robust whether or not we limit the definition of conservation only to those genes close to each other in the genome (within 200 Kb). Results presented have conservation limited in this manner.
Where x is the proportion of genes with conserved synteny within each window, μ is the mean proportion of genes with conserved synteny in all windows and s is the standard deviation of genes with conserved synteny in all windows. Thus the Z-score represents the number of standard deviations above or below the mean.
Comparison with sequence-based synteny blocks
Each highly conserved mouse genomic region was compared to synteny maps based on sequence alignment. Maps of the mouse compared to dog, rat, and human were retrieved from the Ensembl genome browser (Comparative Genomics - Synteny, http://www.ensembl.org/Mus_musculus/Location/Synteny) database. Genomic positions of synteny blocks were retrieved from the database and tabulated.
Statistical analysis of conserved and disease gene pairs
Genes annotated as having conserved synteny or as disease orthologs were examined to determine whether the following gene on the chromosome was similarly annotated. The number of such gene pairs was determined per chromosome. By this definition a run of three genes with similar annotation would be counted as two pairs.
To determine the significance of the frequency of such gene pairs, we compared the observed value to a random distribution. The random distribution was calculated by keeping the both the total number of genes and the number of genes with a given annotation constant, but randomizing the assignment of those annotations over the total number of genes. The annotations were randomized 10,000 times, and the number of neighbors with similar annotations calculated. The significance of the observed value was determined from this simulation i.e. the likelihood of the observed value is the proportion of random values that are greater than or equal to the observed value.
Disease-related ortholog distribution
Human genes with a disease-associated mutant allele were identified from the OMIM Morbid Map database http://www.ncbi.nlm.nih.gov/Omim/getmorbid.cgi. These were cross-referenced to the mouse genome using Ensembl BioMart homology filters to identify mouse orthologs. The distribution of disease-related gene orthologs was analyzed by sliding window analysis, in the same manner as for conserved synteny. For each 20 Mb window the number of disease gene orthologs was divided by the total number of genes in a window to determine the proportion of disease gene orthologs in the window. Z-scores were then calculated for each window using the equation above.
The Pearson's correlation between the microsynteny conservation Z-scores and disease-related ortholog Z-scores was calculated, and significance calculated using the chi-squared test. As a control, we assigned new disease genes at random for each chromosome, keeping the total number of disease genes per chromosome the same. We then recalculated the density of disease genes in each sliding window throughout the genome. We plotted the Z-scores for each window containing these alternate disease orthologs, and compared them to the Z-scores for microsynteny conservation. The Pearson's correlation was then re-calculated, showing no significance for the alternate set of disease genes. Correlation analysis was also performed on window sizes of 10 Mb, staggered by 2.5 Mb, 5 Mb, staggered by 1.25 Mb, 2 Mb, staggered by 0.5 Mb, and 1 Mb, staggered by 0.25 Mb. Randomization trials were performed on 10,000 random annotations for each alternate window size as a control. To account for windows with no genes at the window sizes of 2 Mb and 1 Mb, the correlation analysis was repeated with all windows containing 0 genes removed from the actual and randomized datasets.
The authors wish to thank Stephen Oliver, Simon Hubbard, Ray Boot-Handford, David Robertson, and Stephen Brown for critical reading of the manuscript, and Sam Griffiths-Jones for advice on genomic homology searches.
- Moreno C, Lazar J, Jacob HJ, Kwitek AE: Comparative genomics for detecting human disease genes. Adv Genet. 2008, 60: 655-697.View ArticlePubMedGoogle Scholar
- Gibbs RA, Weinstock GM, Metzker ML, Muzny DM, Sodergren EJ, Scherer S, Scott G, Steffen D, Worley KC, Burch PE, et al: Genome sequence of the Brown Norway rat yields insights into mammalian evolution. Nature. 2004, 428 (6982): 493-521.View ArticlePubMedGoogle Scholar
- Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, Agarwal P, Agarwala R, Ainscough R, Alexandersson M, An P, et al: Initial sequencing and comparative analysis of the mouse genome. Nature. 2002, 420 (6915): 520-562.View ArticlePubMedGoogle Scholar
- Eppig JT, Bult CJ, Kadin JA, Richardson JE, Blake JA, Anagnostopoulos A, Baldarelli RM, Baya M, Beal JS, Bello SM: The Mouse Genome Database (MGD): from genes to mice--a community resource for mouse biology. Nucleic Acids Res. 2005, D471-475. 33 Database
- Bourque G, Pevzner PA, Tesler G: Reconstructing the genomic architecture of ancestral mammals: lessons from human, mouse, and rat genomes. Genome Res. 2004, 14 (4): 507-516.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (1): 98-110.PubMed CentralView ArticlePubMedGoogle Scholar
- Coop G, Myers SR: Live hot, die young: transmission distortion in recombination hotspots. PLoS Genet. 2007, 3 (3): e35-PubMed CentralView ArticlePubMedGoogle Scholar
- Spencer CC, Deloukas P, Hunt S, Mullikin J, Myers S, Silverman B, Donnelly P, Bentley D, McVean G: The influence of recombination on human genetic diversity. PLoS Genet. 2006, 2 (9): e148-PubMed CentralView ArticlePubMedGoogle Scholar
- Guryev V, Smits BM, Belt van de J, Verheul M, Hubner N, Cuppen E: Haplotype block structure is conserved across mammals. PLoS Genet. 2006, 2 (7): e121-PubMed CentralView ArticlePubMedGoogle Scholar
- Cooper GM, Sidow A: Genomic regulatory regions: insights from comparative sequence analysis. Curr Opin Genet Dev. 2003, 13 (6): 604-610.View ArticlePubMedGoogle Scholar
- Venkatesh B, Kirkness EF, Loh YH, Halpern AL, Lee AP, Johnson J, Dandona N, Viswanathan LD, Tay A, Venter JC, et al: Ancient noncoding elements conserved in the human genome. Science. 2006, 314 (5807): 1892-View ArticlePubMedGoogle Scholar
- Visel A, Bristow J, Pennacchio LA: Enhancer identification through comparative genomics. Semin Cell Dev Biol. 2007, 18 (1): 140-152.PubMed CentralView ArticlePubMedGoogle Scholar
- Hsu PW, Huang HD, Hsu SD, Lin LZ, Tsou AP, Tseng CP, Stadler PF, Washietl S, Hofacker IL: miRNAMap: genomic maps of microRNA genes and their target genes in mammalian genomes. Nucleic Acids Res. 2006, D135-139. 34 Database
- Torarinsson E, Yao Z, Wiklund ED, Bramsen JB, Hansen C, Kjems J, Tommerup N, Ruzzo WL, Gorodkin J: Comparative genomics beyond sequence-based alignments: RNA structures in the ENCODE regions. Genome Res. 2007, 18 (2): 242-51.View ArticlePubMedGoogle Scholar
- Springer MS, Murphy WJ: Mammalian evolution and biomedicine: new views from phylogeny. Biol Rev Camb Philos Soc. 2007, 82 (3): 375-392.View ArticlePubMedGoogle Scholar
- Coulouarn C, Factor VM, Thorgeirsson SS: Transforming growth factor-beta gene expression signature in mouse hepatocytes predicts clinical outcome in human cancer. Hepatology. 2008, 47 (6): 2059-2067.PubMed CentralView ArticlePubMedGoogle Scholar
- Borsy A, Podani J, Steger V, Balla B, Horvath A, Kosa JP, Gyurjan I, Molnar A, Szabolcsi Z, Szabo L, et al: Identifying novel genes involved in both deer physiological and human pathological osteoporosis. Mol Genet Genomics. 2009, 281 (3): 301-313.View ArticlePubMedGoogle Scholar
- Malaspina P, Picklo MJ, Jakobs C, Snead OC, Gibson KM: Comparative genomics of aldehyde dehydrogenase 5a1 (succinate semialdehyde dehydrogenase) and accumulation of gamma-hydroxybutyrate associated with its deficiency. Hum Genomics. 2009, 3 (2): 106-120.PubMed CentralView ArticlePubMedGoogle Scholar
- Stoll M, Kwitek-Black AE, Cowley AW, Harris EL, Harrap SB, Krieger JE, Printz MP, Provoost AP, Sassard J, Jacob HJ: New target regions for human hypertension via comparative genomics. Genome Res. 2000, 10 (4): 473-482.PubMed CentralView ArticlePubMedGoogle Scholar
- Twigger S, Lu J, Shimoyama M, Chen D, Pasko D, Long H, Ginster J, Chen CF, Nigam R, Kwitek A, et al: Rat Genome Database (RGD): mapping disease onto the genome. Nucleic Acids Res. 2002, 30 (1): 125-128.PubMed CentralView ArticlePubMedGoogle Scholar
- Korstanje R, DiPetrillo K: Unraveling the genetics of chronic kidney disease using animal models. Am J Physiol Renal Physiol. 2004, 287 (3): F347-352.View ArticlePubMedGoogle Scholar
- Sugiyama F, Yagami K, Paigen B: Mouse models of blood pressure regulation and hypertension. Curr Hypertens Rep. 2001, 3 (1): 41-48.View ArticlePubMedGoogle Scholar
- Hurst LD, Pal C, Lercher MJ: The evolutionary dynamics of eukaryotic gene order. Nat Rev Genet. 2004, 5 (4): 299-310.View ArticlePubMedGoogle Scholar
- Pal C, Hurst LD: Evidence for co-evolution of gene order and recombination rate. Nat Genet. 2003, 33 (3): 392-395.View ArticlePubMedGoogle Scholar
- Batada NN, Hurst LD: Evolution of chromosome organization driven by selection for reduced gene expression noise. Nat Genet. 2007, 39 (8): 945-949.View ArticlePubMedGoogle Scholar
- Kamath RS, Fraser AG, Dong Y, Poulin G, Durbin R, Gotta M, Kanapin A, Le Bot N, Moreno S, Sohrmann M, et al: Systematic functional analysis of the Caenorhabditis elegans genome using RNAi. Nature. 2003, 421 (6920): 231-237.View ArticlePubMedGoogle Scholar
- Lee JM, Sonnhammer EL: Genomic gene clustering analysis of pathways in eukaryotes. Genome Res. 2003, 13 (5): 875-882.PubMed CentralView ArticlePubMedGoogle Scholar
- Petkov PM, Graber JH, Churchill GA, DiPetrillo K, King BL, Paigen K: Evidence of a large-scale functional organization of mammalian chromosomes. PLoS Genet. 2005, 1 (3): e33-PubMed CentralView ArticlePubMedGoogle Scholar
- Hentges KE, Pollock DD, Liu B, Justice MJ: Regional variation in the density of essential genes in mice. PLoS Genet. 2007, 3 (5): e72-PubMed CentralView ArticlePubMedGoogle Scholar
- Rosenthal N, Brown S: The mouse ascending: perspectives for human-disease models. Nat Cell Biol. 2007, 9 (9): 993-999.View ArticlePubMedGoogle Scholar
- Hubbard TJ, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T: Ensembl 2007. Nucleic Acids Res. 2007, D610-617. 35 Database
- Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, et al: Initial sequencing and analysis of the human genome. Nature. 2001, 409 (6822): 860-921.View ArticlePubMedGoogle Scholar
- Venter JC, Adams MD, Myers EW, Li PW, Mural RJ, Sutton GG, Smith HO, Yandell M, Evans CA, Holt RA, et al: The sequence of the human genome. Science. 2001, 291 (5507): 1304-1351.View ArticlePubMedGoogle Scholar
- Kirkness EF, Bafna V, Halpern AL, Levy S, Remington K, Rusch DB, Delcher AL, Pop M, Wang W, Fraser CM, et al: The dog genome: survey sequencing and comparative analysis. Science. 2003, 301 (5641): 1898-1903.View ArticlePubMedGoogle Scholar
- Lunter G, Ponting CP, Hein J: Genome-wide identification of human functional DNA using a neutral indel model. PLoS Comput Biol. 2006, 2 (1): e5-PubMed CentralView ArticlePubMedGoogle Scholar
- Kent WJ: BLAT--the BLAST-like alignment tool. Genome Res. 2002, 12 (4): 656-664.PubMed CentralView ArticlePubMedGoogle Scholar
- Rinchik EM, Carpenter DA: N-ethyl-N-nitrosourea mutagenesis of a 6- to 11-cM subregion of the Fah-Hbb interval of mouse chromosome 7: Completed testing of 4557 gametes and deletion mapping and complementation analysis of 31 mutations. Genetics. 1999, 152 (1): 373-383.PubMed CentralPubMedGoogle Scholar
- Kile BT, Hentges KE, Clark AT, Nakamura H, Salinger AP, Liu B, Box N, Stockton DW, Johnson RL, Behringer RR, et al: Functional genetic analysis of mouse chromosome 11. Nature. 2003, 425 (6953): 81-86.View ArticlePubMedGoogle Scholar
- Wilson L, Ching YH, Farias M, Hartford SA, Howell G, Shao H, Bucan M, Schimenti JC: Random mutagenesis of proximal mouse chromosome 5 uncovers predominantly embryonic lethal mutations. Genome Res. 2005, 15 (8): 1095-1105.PubMed CentralView ArticlePubMedGoogle Scholar
- Hentges KE, Nakamura H, Furuta Y, Yu Y, Thompson DM, O'Brien W, Bradley A, Justice MJ: Novel lethal mouse mutants produced in balancer chromosome screens. Gene Expr Patterns. 2006, 6 (6): 653-665.View ArticlePubMedGoogle Scholar
- Goodstadt L, Ponting CP: Phylogenetic reconstruction of orthology, paralogy, and conserved synteny for dog and human. PLoS Comput Biol. 2006, 2 (9): e133-PubMed CentralView ArticlePubMedGoogle Scholar
- Furney SJ, Alba MM, Lopez-Bigas N: Differences in the evolutionary history of disease genes affected by dominant or recessive mutations. BMC Genomics. 2006, 7: 165-PubMed CentralView ArticlePubMedGoogle Scholar
- Lopez-Bigas N, Ouzounis CA: Genome-wide identification of genes likely to be involved in human genetic disease. Nucleic Acids Res. 2004, 32 (10): 3108-3114.PubMed CentralView ArticlePubMedGoogle Scholar
- Kondrashov FA, Ogurtsov AY, Kondrashov AS: Bioinformatical assay of human gene morbidity. Nucleic Acids Res. 2004, 32 (5): 1731-1737.PubMed CentralView ArticlePubMedGoogle Scholar
- Nielsen R, Hellmann I, Hubisz M, Bustamante C, Clark AG: Recent and ongoing selection in the human genome. Nat Rev Genet. 2007, 8 (11): 857-868.PubMed CentralView ArticlePubMedGoogle Scholar
- Bustamante CD, Fledel-Alon A, Williamson S, Nielsen R, Hubisz MT, Glanowski S, Tanenbaum DM, White TJ, Sninsky JJ, Hernandez RD, et al: Natural selection on protein-coding genes in the human genome. Nature. 2005, 437 (7062): 1153-1157.View ArticlePubMedGoogle Scholar
- Smith AV, Thomas DJ, Munro HM, Abecasis GR: Sequence features in regions of weak and strong linkage disequilibrium. Genome Res. 2005, 15 (11): 1519-1534.PubMed CentralView ArticlePubMedGoogle Scholar
- Chuang JH, Li H: Functional bias and spatial organization of genes in mutational hot and cold regions in the human genome. PLoS Biol. 2004, 2 (2): E29-PubMed CentralView ArticlePubMedGoogle Scholar
- Woolfe A, Goodson M, Goode DK, Snell P, McEwen GK, Vavouri T, Smith SF, North P, Callaway H, Kelly K, et al: Highly conserved non-coding sequences are associated with vertebrate development. PLoS Biol. 2005, 3 (1): e7-PubMed CentralView ArticlePubMedGoogle Scholar
- Huang H, Winter EE, Wang H, Weinstock KG, Xing H, Goodstadt L, Stenson PD, Cooper DN, Smith D, Alba MM, et al: Evolutionary conservation and selection of human disease gene orthologs in the rat and mouse genomes. Genome Biol. 2004, 5 (7): R47-PubMed CentralView ArticlePubMedGoogle Scholar
- Hossain Z, Ali SM, Ko HL, Xu J, Ng CP, Guo K, Qi Z, Ponniah S, Hong W, Hunziker W: Glomerulocystic kidney disease in mice with a targeted inactivation of Wwtr1. Proc Natl Acad Sci USA. 2007, 104 (5): 1631-1636.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu L, Gu S, Alappat S, Song Y, Yan M, Zhang X, Zhang G, Jiang Y, Zhang Z, Zhang Y, et al: Shox2-deficient mice exhibit a rare type of incomplete clefting of the secondary palate. Development. 2005, 132 (19): 4397-4406.View ArticlePubMedGoogle Scholar
- Water N, Tan T, Ashton F, O'Grady A, Day T, Browett P, Ockelford P, Harper P: Mutations within the protein Z-dependent protease inhibitor gene are associated with venous thromboembolic disease: a new form of thrombophilia. Br J Haematol. 2004, 127 (2): 190-194.View ArticlePubMedGoogle Scholar
- Singer GA, Lloyd AT, Huminiecki LB, Wolfe KH: Clusters of co-expressed genes in mammalian genomes are conserved by natural selection. Mol Biol Evol. 2005, 22 (3): 767-775.View ArticlePubMedGoogle Scholar
- Semon M, Duret L: Evolutionary origin and maintenance of coexpressed gene clusters in mammals. Mol Biol Evol. 2006, 23 (9): 1715-1723.View ArticlePubMedGoogle Scholar
- Kapranov P, Willingham AT, Gingeras TR: Genome-wide transcription and the implications for genomic organization. Nat Rev Genet. 2007, 8 (6): 413-423.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.