BACKGROUND
deciphering the evolution of eukaryotic genes is a key to understanding their function in different species.
as nearly all mammals are diploid organisms, most of their genes are expressed from both parental alleles.
nevertheless, a number of genes do not obey this rule, notably so-called imprinted genes.
they acquire specific epigenetic marks in the parental germ lines which are the cause of monoallelic expression after fertilization.
consequently, it depends on the parental origin which allele of an imprinted gene gets inactivated.
genomic imprinting has been observed in the endosperm of flowering plants and in the animal kingdom.
in the latter, imprinted genes have so far only been discovered in therian species  <cit> .
monoallelic expression may thus be evolutionary advantageous for particular functions.
based on the finding that genomic imprinting is a specific feature of species in which embryo and mother are in direct contact to each other, imprinted genes were hypothesized to regulate maternal nutrient supply during embryonic development.
the kinship theory names the different interests of the parents in the wellbeing of the embryo as a possible evolutionary origin of imprinting  <cit> .
according to this model, the mother aims at saving resources for subsequent pregnancies that are possibly fathered by different males, whereas the father is interested in the maximal exploitation of maternal resources in favor of his own offspring.
this would result in a selective pressure towards silencing of growth promoting genes in the female germ line, and of growth suppressing genes in the male germ line.
in both cases the result would be monoallelic expression of the respective genes after fertilization.
due to their monoallelic silencing, imprinted genes might react differently on natural selection than biallelically expressed genes.
in fact, their functional haploidy is reminiscent of the scenario for genes on the x chromosome, on which both positive and negative selection act more efficiently than on the autosomes .
deleterious mutations affecting the expressed allele will be subject to selective elimination whereas beneficial ones would provide immediate advantages.
on the other hand, the inactive allele may accumulate mutations that remain unexposed as long as the sex of the transmitting parent does not switch  <cit> .
imprinted genes have been supposed to play key roles in the mammalian embryo and placenta.
hence, evolution of imprinting and speciation of placental mammals might be linked.
evolution of imprinting might have subjected the affected genes to strong purifying selection, or might have triggered a strong positive selection towards species-specific functions.
fitting into the second scenario it has been reported that the evolution of imprinting of the medea gene in arabidopsis species might have been initiated by a duplication event and that the evolution of the imprinted duplicate was accompanied by strong positive selection  <cit> .
interestingly, also the evolution of mammalian imprinted genes was suggested to be influenced or even triggered by duplication events  <cit> .
studies on a limited number of imprinted genes in mouse and rat did not provide evidence for conspicuous mutation rates or positive selection in the rodent lineage  <cit> .
due to their functions as growth regulators, special attention has been paid to the imprinted genes igf <dig> and igf2r.
whereas dna sequences that encode the interface region of the igf2r and igf <dig> proteins are highly conserved, the signal sequence of igf2r that determines the protein's location in the cell is strikingly divergent between mouse and rat as well as between human and cow  <cit> .
indication of positive selection on igf <dig> in viviparous fish species that developed placenta-like structures suggests that evolutionary adaptations of growth factor encoding genes might predate their imprinting  <cit> .
to date, about  <dig> imprinted genes have been identified in human and mouse  <cit> .
in order to assess their evolutionary history in placental mammals, we analyzed sequence conservation of maternally and paternally expressed genes and their paralogs in a genome-wide comparison.
we show that the evolution of imprinted genes is characterized by an accelerated divergence at dna and protein level in rodent ancestors.
moreover, the existence of paralogous genes seems to have facilitated divergence of imprinted genes.
RESULTS
maternally expressed genes show reduced conservation
in order to assess the evolutionary conservation of imprinted genes, we selected  <dig> pairs of orthologous protein encoding human and mouse genes for which imprinting has been reported in at least one of the two species in the otago catalogue of imprinted genes  <cit>  and in the literature .
these imprinted genes were compared on a genome-wide scale to all genes annotated in the homologene database  <cit> .
this database provides information derived from alignments of orthologous refseq cdna and protein sequences.
we extracted identity values for the alignments of human cdna and amino acid sequences, respectively, to their orthologs of mouse, rat, chimpanzee, dog, cow, and chicken.
compared to genome-wide data, human and mouse ortholog pairs of maternally expressed genes are less conserved in terms of cdna identity  and the encoded proteins show a trend towards reduced identity .
additionally, a prominent feature was the increased number of gaps - corresponding to insertions or deletions of amino acids - in the sequence alignments of proteins encoded by maternally expressed genes .
this might hint at an increase of exon insertions or deletions, respectively.
a ka/ks rates are undetermined for  <dig> genes.
* p <  <dig> , ** p <  <dig>
comparing the human to dog, cow, and chimpanzee, maternally expressed genes exhibit the lowest identities on dna and protein levels, thereby supporting an increased divergence of maternally expressed genes .
however, the reduced conservation between human and non-rodent species is not statistically significant .
paternally expressed genes have a similar level of sequence conservation as the genomic background in all pairwise comparisons of the human to other mammalian species.
as imprinted genes are associated with particular, allele-specific dna methylation patterns, the observed divergence of their protein-encoding sequences might be due to an increased rate of cpg to tpg transitions.
addressing this issue we investigated the frequencies of silent cpg to tpg transitions at positions in the cdna where these mutations would not change the encoded amino acid  <cit> .
we found that both maternally expressed genes and the whole imprinted group exhibit insignificantly lower levels of silent cpg-tpg mismatches in the cdna alignments compared to genome-wide data .
hence, the increased divergence is not caused by increased deamination of methylated cpgs.
divergence of imprinted genes between rodents and other mammals
the low sequence identity of maternally expressed orthologs results apparently from an accelerated divergence in mouse or rodents.
further evidence for the latter hypothesis is provided by comparison of human and rat genes.
the obtained data are consistent with those from mouse .
here, the reduced conservation of maternally expressed genes is even more significant , which can be attributed to a higher number of substitutions in the rat  <cit> .
also between mouse and cow the conservation of all imprinted and the maternally expressed genes is lower than that of non-imprinted genes on the protein and cdna levels .
comparisons with other species, for which fewer sequences of imprinted genes are available, show a tendency towards increased divergence of murine imprinted genes as well .
interestingly, comparing mammalian genes to their orthologs in chicken, a species without imprinting effects, did not reveal any changes in the conservation of imprinted genes that might be associated with the evolution of imprinting in early mammals.
increased divergence of maternally expressed genes might have its cause in increased mutation rates, reduced purifying selection, or positive selection.
increased mutation rates are indicated by an elevated rate of synonymous substitutions   <cit> .
for human-mouse  and mouse-cow  gene pairs, ks rates are essentially similar in all groups, thereby contradicting a major influence of mutation rates.
a commonly used method for estimating selection is determining the ratio of synonymous substitutions  and non-synonymous substitutions  per site in pairwise alignments of coding dna  <cit> .
in general, ka/ks ratios of below  <dig>  indicate purifying selection  <cit>  whereas ka/ks values larger than  <dig> are indicative of positive selection.
the genome-wide median of the ka/ks ratio is  <dig>  for mouse-human gene pairs.
ka/ks tends to be elevated for the group of  <dig> maternally expressed genes  but this is not the case for the  <dig> paternally expressed ones .
the imprinted genes with the most elevated ka/ks values are the maternally expressed genes cdkn1c  and phlda <dig> , and the paternally expressed usp <dig> .
as indicated by very low ka/ks values, the strongest purifying selection appears to act on the paternally expressed genes snrpn, mest, wt <dig>  and copg <dig> .
the ka/ks ratios of the  <dig> imprinted orthologs available for mouse and cow are tentatively elevated as well .
in summary, the set of imprinted genes studied here does not contain genes with exceptionally high ka/ks ratios that indicate recent positive selection.
strong conservation of imprinted genes in modern rodents
next, we wanted to know whether the reduced conservation of imprinted genes is due to divergence in a common ancestor of mouse and rat, or whether this evolutionary process is still ongoing in modern rodents.
intriguingly, comparison of mouse and rat orthologs  shows that the conservation levels of imprinted genes are similar to those of all genes on the protein level , and even higher on dna level .
the latter fact is mostly caused by the paternally expressed genes whereas maternally expressed genes are similarly conserved as the genomic background.
the ka/ks ratio of imprinted genes  is not significantly higher than genome-wide .
murine imprinted genes contain a median of  <dig>  single nucleotide polymorphisms  per  <dig> kb of coding sequence whereas the genome-wide median is  <dig>  .
this depletion, which is not seen for human imprinted genes  strongly argues for ongoing purifying selection on the protein-coding regions of murine imprinted genes.
aka/ks rates are undetermined for  <dig> genes.
** p <  <dig> , *** p <  <dig>
contrary to the pattern observed for mouse-human gene pairs, mouse-rat imprinted orthologs have a decreased rate of synonymous substitutions  whereas the rate of nonsynonymous changes is not significantly elevated .
this finding holds for both paternally and maternally expressed genes and is in agreement with the results of smith and hurst  <cit> .
reduced ks rates may hint at selection on silent sites related to alternative splicing and rna secondary structure requirements  <cit> .
alternatively, a special chromatin structure of imprinted genes might reduce overall mutation rates in the germ lines, thereby leading to the observed lower synonymous substitution rates.
a similar connection is not seen for human-mouse or human-chimpanzee orthologs , which might point at a rodent-specific evolutionary pattern.
reconstructing patterns of ancient evolution
the elevated conservation of imprinted genes between the two rodent species is in contrast to their divergence from their orthologs in other mammals.
hence, most of the discerning dna changes must have taken place before the split of rat and mouse.
afterwards, imprinted genes were possibly subjected to strict purifying selection in mouse and rat .
since both mechanisms counteract and thus obscure the pattern in the sequences of extant species, we aimed at estimating conservation and evolutionary rates before the split of mouse and rat.
from pairwise alignment data on human-mouse , human-rat  and mouse-rat , we estimated protein and cdna identity as well as ka and ks rates between human and a common ancestor of mouse and rat  using formula  for protein or cdna identity given in percent and formula  for ka or ks rates.
formula  can be applied for sequences where mm_rn >hs_mm and mm_rn >hs_rn, which applies to the vast majority of all genes that are present for all three species in homologene.
in fact, we detected only  <dig> genes where the identity of alignments between human and rodent proteins was higher than that between mouse and rat.
this can be due to differences in the alignments caused by events such as changes in exon usage during murine evolution.
likewise, formula  is not applicable for  <dig> cases in which the ka or ks rates are higher between the two rodent species than between human and rodent sequences.
with regard to the reconstructed rodent ancestor, sequence identities become higher but more discriminative between imprinted and all genes.
comparisons of the obtained values of  <dig> imprinted genes to  <dig>  genome-wide reconstructions gave results consistent with the observations described above.
increased divergence of dna and protein sequences  and elevation of ka  of maternally expressed genes became slightly more significant.
in addition we observed trends toward higher ka values in the whole imprinted set , and higher ks and ka/ks rates of maternally expressed genes .
an alternative option to trace different evolutionary constraints in the mammalian lineages is provided by the paml package  <cit> .
we constructed branch models including the cow as fourth mammalian species.
in the imprinted group, the lineage leading to the rodent ancestor has similar ka/ks ratios as other lineages for  <dig> out of  <dig> genes.
for the remaining eleven genes, the two-ratios model assuming a different ka/ks ratio is significantly more likely than the one-ratio model assuming the same ka/ks ratio for all branches .
for four of these genes , the ka/ks ratio calculated by the two-ratios model is elevated in the lineage leading to the rodent ancestor.
all four genes were found to have ka/ks rates above the median values for human-mouse  indicating that these genes - interestingly members of prominent imprinted regions - have been subject to species-specific shifts in evolutionary patterns.
since the paml estimations of ka/ks are always lower than  <dig> for all lineages, the imprinted genes studied here apparently did not experience strong positive selection but rather relaxation of purifying constraints.
interacting evolution of imprinted genes and their paralogs
identification of general evolutionary principles on the basis of simple sequence features crucially depends on the careful construction of suitable control groups.
having originated from the same ancestral gene, paralogs are supposed to fulfill similar  functions  <cit> .
therefore, comparing imprinted genes to their paralogs seemed an attractive opportunity to reveal differences in evolution that are not related to the biochemical functions of these genes per se but rather to their particular, monoallelic expression pattern.
first, we analyzed whether there is an enrichment of paralogs in the imprinted gene group according to the ensembl release  <dig>  <cit>  annotations for  <dig>  human autosomal protein-coding genes.
to avoid a bias of genes with many paralogs, we chose for each gene the one that is listed first as its representative paralog.
this paralog is the evolutionary most recent and in most cases also the one with the highest identity.
for some imprinted genes  the approach used here gave different results than previously reported in the literature  <cit>  .
this includes duplicates that were reported to originate from retrotransposition events from the x chromosome  <cit> .
nevertheless,  <dig> % of the genes in the imprinted group possess a paralog, a slightly higher percentage than for all genes on human autosomes .
interestingly, most duplication events predate the origin of mammals, also in the imprinted group.
with the exception of ins <dig> and ins <dig>  the imprinted genes studied here were not involved in duplication events after the split of the human and rodent lineages.
several imprinted genes have been linked to paralogs on the x chromosome  <cit>  and it has been speculated that imprinting and x chromosome inactivation may have co-evolved since they require similar mechanisms of gene regulation  <cit> .
we found that the first paralogs of three imprinted genes  are located on the x chromosome.
additional three  have an x-linked paralog that is not the highest scoring one.
the resulting  <dig> % is no significant enrichment compared to the autosomal ratio of  <dig> % .
also in the mouse there is no significant enrichment of x-linked paralogs of imprinted genes.
ain cases where more than one paralog was available, protein identities were analyzed only for the gene and its evolutionary most recent paralog.
* p <  <dig> , ** p <  <dig> , *** p <  <dig>
compared to genome-wide data, human imprinted genes show less identity with their most recent paralogs on protein level .
in the mouse this relaxation in paralog conservation is more pronounced , and is probably caused by the stronger divergence of imprinted genes in the rodent ancestor as described above.
taking the average identity of all paralogs per gene as a measure yields lower identity values without affecting the significance of the differences.
to investigate whether the existence of paralogs might influence the evolution of protein-coding imprinted genes, we investigated the conservation of gene orthologs in presence or absence of paralogs.
in the entire genome and also in the case of imprinted genes, orthologs that possess paralogs are significantly more conserved between human and mouse, or mouse and rat than those without a paralog and have lower ka/ks ratios .
comparing imprinted genes with or without paralogs, respectively, to the corresponding groups of autosomal genes, reveals that imprinted genes with paralogs are subject to decreased conservation between human and rodents  and tend towards a higher ka/ks ratio  whereas there is no significant difference between genes without paralogs in both groups .
between mouse and rat, genes with paralogs behave similarly in both groups  but imprinted ones without paralogs show increased conservation on dna level and a lower ks ratio .
in all comparisons maternally and paternally expressed genes behaved similarly.
overall, genes with high numbers of paralogs did not behave differently from genes with only one or few paralogs.
athe second number refers to sequences available in the homologene database for ka/ks analyses.
* p <  <dig> , ** p <  <dig> , *** p <  <dig>
interestingly, the paralogs show a higher conservation than their imprinted counterparts and all genes in homologene and remarkably lower ka/ks ratios between human-mouse, human-rat, and mouse-rat .
no differences were found between maternally and paternally expressed genes.
purifying selection thus seems to act far more strictly on the paralogs of rodent imprinted genes than on the imprinted genes themselves.
* p <  <dig> , ** p <  <dig> , *** p <  <dig>  .
data refer to the evolutionary most recent paralog of each imprinted gene.
discussion
divergence of imprinted genes in early eutherian evolution
detailed genome-wide analyses on cdna and protein level revealed an intriguing evolutionary pattern of imprinted genes.
compared to their  chicken orthologs, imprinted genes show similar sequence identities and ka/ks distributions as all genes in the homologene database.
this suggests that in early mammals there was no specific evolutionary pressure on the protein-coding sequences of imprinted genes.
in contrast, imprinted genes of mouse and rat show an elevated divergence of coding dna sequences and proteins in comparison to their orthologs in other mammalian species, especially in the human, suggesting that there has been an increased divergence of imprinted genes in a common ancestor of mouse and rat.
although the reduced conservation was only significant for the rodent lineage, we cannot exclude that imprinted genes underwent a relaxation of selective constraints also in other species.
the observed differences might be more pronounced in rodents due to their shorter generation times and effects of their larger population sizes  <cit> .
one also has to keep in mind that the analysis of genes in species such as dog, cow, and chicken that have not been studied as long as human and mouse is hampered by a still low number and quality of available sequences.
in extant rodents, evolutionary processes seem to have shifted towards purifying selection.
at present, murine imprinted genes are apparently under strong purifying selection as suggested by the reduction of snps.
this pattern of evolution, i.e.
an initial divergence followed by fixation, is believed to be typical for the evolution of new functions, for example of duplicated gene copies  <cit>  or for the evolution of new species.
remarkably, the increased conservation in modern rodents coincides with a stricter conservation of typical dna elements in imprinted genes in the mouse than in the human, such as intronic cpg islands  <cit> .
antagonistic feedback effects of maternally expressed genes
interestingly, the increased divergence of imprinted genes affects mostly maternally expressed genes.
especially for rodents, the reduction of embryonic growth by genes expressed from the maternal alleles may have provided evolutionary advantages along with changes in placentation  <cit> .
indeed, the maternally expressed genes cdkn1c and phlda <dig> that fulfill important functions in the mouse placenta have the highest ka/ks ratios of all imprinted genes.
thus, a possible explanation for the low conservation of maternally expressed genes might be the evolution of new functions in different lineages.
however, the rather modest elevation of the ka/ks ratios argues rather for relaxed constraints than for positive selection.
an alternative scenario is the degeneration of  growth repressors.
these are likely to have different, counteracting effects on maternal and embryonic fitness: on the one hand a growth repressor might reduce the fitness of the embryo, on the other hand small embryos might increase the number of offspring.
the counteracting effects on embryonic and maternal fitness should result in a reduced selective pressure on maternally expressed growth regulators, resulting in the observed relaxed conservation among mammalian species.
since there is no physical interaction between the embryo and the father, a similar effect of paternally expressed genes on paternal fitness is hardly conceivable.
in support of this we found that the set of paternally expressed genes shows a similar or even higher level of conservation as the genome-wide background.
paralogs may facilitate divergence of rodent imprinted genes
we did not find evidence that gene duplication events originating on the x chromosomes are a key factor in the evolution of imprinted genes as it has been suggested previously  <cit> .
however, we cannot exclude that the duplication of a few genes on the x chromosome and their translocation to the autosomes might have initiated the evolution of imprinted genomic regions.
elements that initially regulated only imprinting of the duplicated genes may have strengthened and started to influence also the expression of neighboring genes.
interestingly, the six genes in our study that possess x-chromosomal paralogs are distributed over different imprinting domains, thus representing potential originators.
although duplication events might not be directly linked to evolution of imprinted gene regulation, the existence of paralogs may have enabled a greater divergence of imprinted genes by relaxation of purifying selection.
genes that possess paralogs show evolutionary patterns that differ from those of singleton genes  <cit> .
notably, there is complementary divergence between pairs of paralogs, where imprinted genes seem to be concentrated in the group of the faster evolving duplicates .
the strong conservation of their non-imprinted paralogs suggests that these genes may have compensated for the divergence by maintaining the original functions of the ancestral genes.
CONCLUSIONS
the evolutionary fate of maternally expressed genes appears to be influenced by antagonistic effects on maternal and embryonic fitness, resulting in a relaxation of selective pressure compared to biallelically expressed genes.
the observed divergence of imprinted genes might have been compensated by a restrictive conservation of their paralogs.
in most cases the duplication events predated the split of the fish and mammalian lineages and thus the evolution of imprinting.
long-lasting functional redundancy may have allowed genes that later became imprinted to specialize in functions related to placental and embryonic development.
apparent purifying selection on paternally expressed genes in modern rodents might be linked with the intensification of imprinting in species with a pronounced conflict over maternal resources.

 deciphering the evolution of eukaryotic genes is a key to understanding their function in different species. as nearly all mammals are diploid organisms, most of their genes are expressed from both parental alleles. nevertheless, a number of genes do not obey this rule, notably so-called imprinted genes. they acquire specific epigenetic marks in the parental germ lines which are the cause of monoallelic expression after fertilization. consequently, it depends on the parental origin which allele of an imprinted gene gets inactivated. genomic imprinting has been observed in the endosperm of flowering plants and in the animal kingdom. in the latter, imprinted genes have so far only been discovered in therian species. monoallelic expression may thus be evolutionary advantageous for particular functions. based on the finding that genomic imprinting is a specific feature of species in which embryo and mother are in direct contact to each other, imprinted genes were hypothesized to regulate maternal nutrient supply during embryonic development. the kinship theory names the different interests of the parents in the wellbeing of the embryo as a possible evolutionary origin of imprinting. according to this model, the mother aims at saving resources for subsequent pregnancies that are possibly fathered by different males, whereas the father is interested in the maximal exploitation of maternal resources in favor of his own offspring. this would result in a selective pressure towards silencing of growth promoting genes in the female germ line, and of growth suppressing genes in the male germ line. in both cases the result would be monoallelic expression of the respective genes after fertilization. due to their monoallelic silencing, imprinted genes might react differently on natural selection than biallelically expressed genes. in fact, their functional haploidy is reminiscent of the scenario for genes on the x chromosome, on which both positive and negative selection act more efficiently than on the autosomes. deleterious mutations affecting the expressed allele will be subject to selective elimination whereas beneficial ones would provide immediate advantages. on the other hand, the inactive allele may accumulate mutations that remain unexposed as long as the sex of the transmitting parent does not switch. imprinted genes have been supposed to play key roles in the mammalian embryo and placenta. hence, evolution of imprinting and speciation of placental mammals might be linked. evolution of imprinting might have subjected the affected genes to strong purifying selection, or might have triggered a strong positive selection towards species-specific functions. fitting into the second scenario it has been reported that the evolution of imprinting of the medea gene in arabidopsis species might have been initiated by a duplication event and that the evolution of the imprinted duplicate was accompanied by strong positive selection. interestingly, also the evolution of mammalian imprinted genes was suggested to be influenced or even triggered by duplication events. studies on a limited number of imprinted genes in mouse and rat did not provide evidence for conspicuous mutation rates or positive selection in the rodent lineage. due to their functions as growth regulators, special attention has been paid to the imprinted genes igf and igf2r. whereas dna sequences that encode the interface region of the igf2r and igf proteins are highly conserved, the signal sequence of igf2r that determines the protein's location in the cell is strikingly divergent between mouse and rat as well as between human and cow. indication of positive selection on igf in viviparous fish species that developed placenta-like structures suggests that evolutionary adaptations of growth factor encoding genes might predate their imprinting. to date, about imprinted genes have been identified in human and mouse. in order to assess their evolutionary history in placental mammals, we analyzed sequence conservation of maternally and paternally expressed genes and their paralogs in a genome-wide comparison. we show that the evolution of imprinted genes is characterized by an accelerated divergence at dna and protein level in rodent ancestors. moreover, the existence of paralogous genes seems to have facilitated divergence of imprinted genes.  maternally expressed genes show reduced conservation in order to assess the evolutionary conservation of imprinted genes, we selected pairs of orthologous protein encoding human and mouse genes for which imprinting has been reported in at least one of the two species in the otago catalogue of imprinted genes and in the literature. these imprinted genes were compared on a genome-wide scale to all genes annotated in the homologene database. this database provides information derived from alignments of orthologous refseq cdna and protein sequences. we extracted identity values for the alignments of human cdna and amino acid sequences, respectively, to their orthologs of mouse, rat, chimpanzee, dog, cow, and chicken. compared to genome-wide data, human and mouse ortholog pairs of maternally expressed genes are less conserved in terms of cdna identity and the encoded proteins show a trend towards reduced identity. additionally, a prominent feature was the increased number of gaps - corresponding to insertions or deletions of amino acids - in the sequence alignments of proteins encoded by maternally expressed genes. this might hint at an increase of exon insertions or deletions, respectively. a ka/ks rates are undetermined for genes. * p <, ** p < comparing the human to dog, cow, and chimpanzee, maternally expressed genes exhibit the lowest identities on dna and protein levels, thereby supporting an increased divergence of maternally expressed genes. however, the reduced conservation between human and non-rodent species is not statistically significant. paternally expressed genes have a similar level of sequence conservation as the genomic background in all pairwise comparisons of the human to other mammalian species. as imprinted genes are associated with particular, allele-specific dna methylation patterns, the observed divergence of their protein-encoding sequences might be due to an increased rate of cpg to tpg transitions. addressing this issue we investigated the frequencies of silent cpg to tpg transitions at positions in the cdna where these mutations would not change the encoded amino acid. we found that both maternally expressed genes and the whole imprinted group exhibit insignificantly lower levels of silent cpg-tpg mismatches in the cdna alignments compared to genome-wide data. hence, the increased divergence is not caused by increased deamination of methylated cpgs. divergence of imprinted genes between rodents and other mammals the low sequence identity of maternally expressed orthologs results apparently from an accelerated divergence in mouse or rodents. further evidence for the latter hypothesis is provided by comparison of human and rat genes. the obtained data are consistent with those from mouse. here, the reduced conservation of maternally expressed genes is even more significant, which can be attributed to a higher number of substitutions in the rat. also between mouse and cow the conservation of all imprinted and the maternally expressed genes is lower than that of non-imprinted genes on the protein and cdna levels. comparisons with other species, for which fewer sequences of imprinted genes are available, show a tendency towards increased divergence of murine imprinted genes as well. interestingly, comparing mammalian genes to their orthologs in chicken, a species without imprinting effects, did not reveal any changes in the conservation of imprinted genes that might be associated with the evolution of imprinting in early mammals. increased divergence of maternally expressed genes might have its cause in increased mutation rates, reduced purifying selection, or positive selection. increased mutation rates are indicated by an elevated rate of synonymous substitutions. for human-mouse and mouse-cow gene pairs, ks rates are essentially similar in all groups, thereby contradicting a major influence of mutation rates. a commonly used method for estimating selection is determining the ratio of synonymous substitutions and non-synonymous substitutions per site in pairwise alignments of coding dna. in general, ka/ks ratios of below indicate purifying selection whereas ka/ks values larger than are indicative of positive selection. the genome-wide median of the ka/ks ratio is for mouse-human gene pairs. ka/ks tends to be elevated for the group of maternally expressed genes but this is not the case for the paternally expressed ones. the imprinted genes with the most elevated ka/ks values are the maternally expressed genes cdkn1c and phlda, and the paternally expressed usp. as indicated by very low ka/ks values, the strongest purifying selection appears to act on the paternally expressed genes snrpn, mest, wt and copg. the ka/ks ratios of the imprinted orthologs available for mouse and cow are tentatively elevated as well. in summary, the set of imprinted genes studied here does not contain genes with exceptionally high ka/ks ratios that indicate recent positive selection. strong conservation of imprinted genes in modern rodents next, we wanted to know whether the reduced conservation of imprinted genes is due to divergence in a common ancestor of mouse and rat, or whether this evolutionary process is still ongoing in modern rodents. intriguingly, comparison of mouse and rat orthologs shows that the conservation levels of imprinted genes are similar to those of all genes on the protein level, and even higher on dna level. the latter fact is mostly caused by the paternally expressed genes whereas maternally expressed genes are similarly conserved as the genomic background. the ka/ks ratio of imprinted genes is not significantly higher than genome-wide. murine imprinted genes contain a median of single nucleotide polymorphisms per kb of coding sequence whereas the genome-wide median is. this depletion, which is not seen for human imprinted genes strongly argues for ongoing purifying selection on the protein-coding regions of murine imprinted genes. aka/ks rates are undetermined for genes. ** p <, *** p < contrary to the pattern observed for mouse-human gene pairs, mouse-rat imprinted orthologs have a decreased rate of synonymous substitutions whereas the rate of nonsynonymous changes is not significantly elevated. this finding holds for both paternally and maternally expressed genes and is in agreement with the results of smith and hurst. reduced ks rates may hint at selection on silent sites related to alternative splicing and rna secondary structure requirements. alternatively, a special chromatin structure of imprinted genes might reduce overall mutation rates in the germ lines, thereby leading to the observed lower synonymous substitution rates. a similar connection is not seen for human-mouse or human-chimpanzee orthologs, which might point at a rodent-specific evolutionary pattern. reconstructing patterns of ancient evolution the elevated conservation of imprinted genes between the two rodent species is in contrast to their divergence from their orthologs in other mammals. hence, most of the discerning dna changes must have taken place before the split of rat and mouse. afterwards, imprinted genes were possibly subjected to strict purifying selection in mouse and rat. since both mechanisms counteract and thus obscure the pattern in the sequences of extant species, we aimed at estimating conservation and evolutionary rates before the split of mouse and rat. from pairwise alignment data on human-mouse, human-rat and mouse-rat, we estimated protein and cdna identity as well as ka and ks rates between human and a common ancestor of mouse and rat using formula for protein or cdna identity given in percent and formula for ka or ks rates. formula can be applied for sequences where mm_rn >hs_mm and mm_rn >hs_rn, which applies to the vast majority of all genes that are present for all three species in homologene. in fact, we detected only genes where the identity of alignments between human and rodent proteins was higher than that between mouse and rat. this can be due to differences in the alignments caused by events such as changes in exon usage during murine evolution. likewise, formula is not applicable for cases in which the ka or ks rates are higher between the two rodent species than between human and rodent sequences. with regard to the reconstructed rodent ancestor, sequence identities become higher but more discriminative between imprinted and all genes. comparisons of the obtained values of imprinted genes to genome-wide reconstructions gave results consistent with the observations described above. increased divergence of dna and protein sequences and elevation of ka of maternally expressed genes became slightly more significant. in addition we observed trends toward higher ka values in the whole imprinted set, and higher ks and ka/ks rates of maternally expressed genes. an alternative option to trace different evolutionary constraints in the mammalian lineages is provided by the paml package. we constructed branch models including the cow as fourth mammalian species. in the imprinted group, the lineage leading to the rodent ancestor has similar ka/ks ratios as other lineages for out of genes. for the remaining eleven genes, the two-ratios model assuming a different ka/ks ratio is significantly more likely than the one-ratio model assuming the same ka/ks ratio for all branches. for four of these genes, the ka/ks ratio calculated by the two-ratios model is elevated in the lineage leading to the rodent ancestor. all four genes were found to have ka/ks rates above the median values for human-mouse indicating that these genes - interestingly members of prominent imprinted regions - have been subject to species-specific shifts in evolutionary patterns. since the paml estimations of ka/ks are always lower than for all lineages, the imprinted genes studied here apparently did not experience strong positive selection but rather relaxation of purifying constraints. interacting evolution of imprinted genes and their paralogs identification of general evolutionary principles on the basis of simple sequence features crucially depends on the careful construction of suitable control groups. having originated from the same ancestral gene, paralogs are supposed to fulfill similar functions. therefore, comparing imprinted genes to their paralogs seemed an attractive opportunity to reveal differences in evolution that are not related to the biochemical functions of these genes per se but rather to their particular, monoallelic expression pattern. first, we analyzed whether there is an enrichment of paralogs in the imprinted gene group according to the ensembl release annotations for human autosomal protein-coding genes. to avoid a bias of genes with many paralogs, we chose for each gene the one that is listed first as its representative paralog. this paralog is the evolutionary most recent and in most cases also the one with the highest identity. for some imprinted genes the approach used here gave different results than previously reported in the literature. this includes duplicates that were reported to originate from retrotransposition events from the x chromosome. nevertheless, % of the genes in the imprinted group possess a paralog, a slightly higher percentage than for all genes on human autosomes. interestingly, most duplication events predate the origin of mammals, also in the imprinted group. with the exception of ins and ins the imprinted genes studied here were not involved in duplication events after the split of the human and rodent lineages. several imprinted genes have been linked to paralogs on the x chromosome and it has been speculated that imprinting and x chromosome inactivation may have co-evolved since they require similar mechanisms of gene regulation. we found that the first paralogs of three imprinted genes are located on the x chromosome. additional three have an x-linked paralog that is not the highest scoring one. the resulting % is no significant enrichment compared to the autosomal ratio of %. also in the mouse there is no significant enrichment of x-linked paralogs of imprinted genes. ain cases where more than one paralog was available, protein identities were analyzed only for the gene and its evolutionary most recent paralog. * p <, ** p <, *** p < compared to genome-wide data, human imprinted genes show less identity with their most recent paralogs on protein level. in the mouse this relaxation in paralog conservation is more pronounced, and is probably caused by the stronger divergence of imprinted genes in the rodent ancestor as described above. taking the average identity of all paralogs per gene as a measure yields lower identity values without affecting the significance of the differences. to investigate whether the existence of paralogs might influence the evolution of protein-coding imprinted genes, we investigated the conservation of gene orthologs in presence or absence of paralogs. in the entire genome and also in the case of imprinted genes, orthologs that possess paralogs are significantly more conserved between human and mouse, or mouse and rat than those without a paralog and have lower ka/ks ratios. comparing imprinted genes with or without paralogs, respectively, to the corresponding groups of autosomal genes, reveals that imprinted genes with paralogs are subject to decreased conservation between human and rodents and tend towards a higher ka/ks ratio whereas there is no significant difference between genes without paralogs in both groups. between mouse and rat, genes with paralogs behave similarly in both groups but imprinted ones without paralogs show increased conservation on dna level and a lower ks ratio. in all comparisons maternally and paternally expressed genes behaved similarly. overall, genes with high numbers of paralogs did not behave differently from genes with only one or few paralogs. athe second number refers to sequences available in the homologene database for ka/ks analyses. * p <, ** p <, *** p < interestingly, the paralogs show a higher conservation than their imprinted counterparts and all genes in homologene and remarkably lower ka/ks ratios between human-mouse, human-rat, and mouse-rat. no differences were found between maternally and paternally expressed genes. purifying selection thus seems to act far more strictly on the paralogs of rodent imprinted genes than on the imprinted genes themselves. * p <, ** p <, *** p <. data refer to the evolutionary most recent paralog of each imprinted gene. discussion divergence of imprinted genes in early eutherian evolution detailed genome-wide analyses on cdna and protein level revealed an intriguing evolutionary pattern of imprinted genes. compared to their chicken orthologs, imprinted genes show similar sequence identities and ka/ks distributions as all genes in the homologene database. this suggests that in early mammals there was no specific evolutionary pressure on the protein-coding sequences of imprinted genes. in contrast, imprinted genes of mouse and rat show an elevated divergence of coding dna sequences and proteins in comparison to their orthologs in other mammalian species, especially in the human, suggesting that there has been an increased divergence of imprinted genes in a common ancestor of mouse and rat. although the reduced conservation was only significant for the rodent lineage, we cannot exclude that imprinted genes underwent a relaxation of selective constraints also in other species. the observed differences might be more pronounced in rodents due to their shorter generation times and effects of their larger population sizes. one also has to keep in mind that the analysis of genes in species such as dog, cow, and chicken that have not been studied as long as human and mouse is hampered by a still low number and quality of available sequences. in extant rodents, evolutionary processes seem to have shifted towards purifying selection. at present, murine imprinted genes are apparently under strong purifying selection as suggested by the reduction of snps. this pattern of evolution, i.e. an initial divergence followed by fixation, is believed to be typical for the evolution of new functions, for example of duplicated gene copies or for the evolution of new species. remarkably, the increased conservation in modern rodents coincides with a stricter conservation of typical dna elements in imprinted genes in the mouse than in the human, such as intronic cpg islands. antagonistic feedback effects of maternally expressed genes interestingly, the increased divergence of imprinted genes affects mostly maternally expressed genes. especially for rodents, the reduction of embryonic growth by genes expressed from the maternal alleles may have provided evolutionary advantages along with changes in placentation. indeed, the maternally expressed genes cdkn1c and phlda that fulfill important functions in the mouse placenta have the highest ka/ks ratios of all imprinted genes. thus, a possible explanation for the low conservation of maternally expressed genes might be the evolution of new functions in different lineages. however, the rather modest elevation of the ka/ks ratios argues rather for relaxed constraints than for positive selection. an alternative scenario is the degeneration of growth repressors. these are likely to have different, counteracting effects on maternal and embryonic fitness: on the one hand a growth repressor might reduce the fitness of the embryo, on the other hand small embryos might increase the number of offspring. the counteracting effects on embryonic and maternal fitness should result in a reduced selective pressure on maternally expressed growth regulators, resulting in the observed relaxed conservation among mammalian species. since there is no physical interaction between the embryo and the father, a similar effect of paternally expressed genes on paternal fitness is hardly conceivable. in support of this we found that the set of paternally expressed genes shows a similar or even higher level of conservation as the genome-wide background. paralogs may facilitate divergence of rodent imprinted genes we did not find evidence that gene duplication events originating on the x chromosomes are a key factor in the evolution of imprinted genes as it has been suggested previously. however, we cannot exclude that the duplication of a few genes on the x chromosome and their translocation to the autosomes might have initiated the evolution of imprinted genomic regions. elements that initially regulated only imprinting of the duplicated genes may have strengthened and started to influence also the expression of neighboring genes. interestingly, the six genes in our study that possess x-chromosomal paralogs are distributed over different imprinting domains, thus representing potential originators. although duplication events might not be directly linked to evolution of imprinted gene regulation, the existence of paralogs may have enabled a greater divergence of imprinted genes by relaxation of purifying selection. genes that possess paralogs show evolutionary patterns that differ from those of singleton genes. notably, there is complementary divergence between pairs of paralogs, where imprinted genes seem to be concentrated in the group of the faster evolving duplicates. the strong conservation of their non-imprinted paralogs suggests that these genes may have compensated for the divergence by maintaining the original functions of the ancestral genes.  the evolutionary fate of maternally expressed genes appears to be influenced by antagonistic effects on maternal and embryonic fitness, resulting in a relaxation of selective pressure compared to biallelically expressed genes. the observed divergence of imprinted genes might have been compensated by a restrictive conservation of their paralogs. in most cases the duplication events predated the split of the fish and mammalian lineages and thus the evolution of imprinting. long-lasting functional redundancy may have allowed genes that later became imprinted to specialize in functions related to placental and embryonic development. apparent purifying selection on paternally expressed genes in modern rodents might be linked with the intensification of imprinting in species with a pronounced conflict over maternal resources. 
moreover, intensification of imprinting in modern rodents seems to have shifted the evolutionary fate of imprinted genes towards strong purifying selection. here, we aimed at assessing the evolutionary history of imprinted genes. in contrast to the majority of mammalian genes, imprinted genes are monoallelically expressed with the choice of the active allele depending on its parental origin. it is contrasted by strong conservation of paternally expressed genes in mouse and rat. in this study, we investigated the conservation of imprinted genes in vertebrate genomes and their exposition to natural selection. the apparent degeneration of maternally expressed genes may reflect a relaxation of selective pressure due to counteracting effects on maternal and embryonic fitness. in a genome-wide comparison, orthologs of imprinted genes show a stronger divergence on cdna and protein level in mammals. functional redundancy provided by the presence of highly conserved paralogs may have facilitated the divergence.
 in contrast to the majority of mammalian genes, imprinted genes are monoallelically expressed with the choice of the active allele depending on its parental origin. due to their special inheritance patterns, maternally and paternally expressed genes might be under different evolutionary pressure. here, we aimed at assessing the evolutionary history of imprinted genes.  in this study, we investigated the conservation of imprinted genes in vertebrate genomes and their exposition to natural selection. in a genome-wide comparison, orthologs of imprinted genes show a stronger divergence on cdna and protein level in mammals. this pattern is most pronounced for maternally expressed genes in rodents in comparison to their non-rodent orthologs. the divergence is not attributable to increased mutation of cpg positions. it is contrasted by strong conservation of paternally expressed genes in mouse and rat. interestingly, we found that the early divergence of imprinted genes was accompanied by an unusually strict conservation of their paralogs.  the apparent degeneration of maternally expressed genes may reflect a relaxation of selective pressure due to counteracting effects on maternal and embryonic fitness. functional redundancy provided by the presence of highly conserved paralogs may have facilitated the divergence. moreover, intensification of imprinting in modern rodents seems to have shifted the evolutionary fate of imprinted genes towards strong purifying selection. 