Molecular evolution of <Emphasis Type="Italic">PKD2</Emphasis> gene family in mammals

PKD2 gene encodes a critical cation channel protein that plays important roles in various developmental processes and is usually evolutionarily conserved. In the present study, we analyzed the evolutionary patterns of PKD2 and its homologous genes (PKD2L1, PKD2L2) from nine mammalian species. In this study, we demonstrated the orthologs of PKD2 gene family evolved under a dominant purifying selection force. Our results in combination with the reported evidences from functional researches suggested the entire PKD2 gene family are conserved and perform essential biological roles during mammalian evolution. In rodents, PKD2 gene family members appeared to have evolved more rapidly than other mammalian lineages, probably resulting from relaxation of purifying selection. However, positive selection imposed on synonymous sites also potentially contributed to this case. For the paralogs, our results implied that PKD2L2 genes evolved under a weaker purifying selection constraint than PKD2 and PKD2L1 genes. Interestingly, some loop regions of transmembrane domain of PKD2L2 exhibited higher P _N/P _S ratios than expected, suggesting these regions are more functional divergent in organisms and worthy of special attention. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Chun Ye, Huan Sun have contributed equally to this work.     

Significance of Population Size on the Fixation of Nonsynonymous Mutations in Genes Under Varying Levels of Selection Pressure

Previous studies observed a higher ratio of divergences at nonsynonymous and synonymous sites (ω = d_N/d_S) in species with a small population size compared to that estimated for those with a large population size. Here we examined the theoretical relationship between ω, effective population size (N_e), and selection coefficient (s). Our analysis revealed that when purifying selection is high, ω of species with small N_e is much higher than that of species with large N_e. However the difference between the two ω reduces with the decline in selection pressure (s → 0). We examined this relationship using primate and rodent genes and found that the ω estimated for highly constrained genes of primates was up to 2.9 times higher than that obtained for their orthologous rodent genes. Conversely, for genes under weak purifying selection the ω of primates was only 17% higher than that of rodents. When tissue specificity was used as a proxy for selection pressure we found that the ω of broadly expressed genes of primates was up to 2.1-fold higher than that of their rodent counterparts and this difference was only 27% for tissue specific genes. Since most of the nonsynonymous mutations in constrained or broadly expressed genes are deleterious, fixation of these mutations is influenced by N_e. This results in a higher ω of these genes in primates compared to those from rodents. Conversely, the majority of nonsynonymous mutations in less-constrained or tissue-specific genes are neutral or nearly neutral and therefore fixation of them is largely independent of N_e, which leads to the similarity of ω in primates and rodents.     

Remarkable expansions of an X-linked reproductive homeobox gene cluster in rodent evolution

Rhox is a recently identified cluster of 12 X-linked homeobox genes in mice. The expression pattern of Rhox genes during postnatal testis development corresponds to their chromosomal position, much like the colinear gene regulation of the Hox gene clusters during animal embryonic development. We here report the identification of 18 additional Rhox genes and 3 pseudogenes in mice. Comparative analyses of the mouse, rat, human, dog, cow, opossum, and chicken genomes suggest that the Rhox cluster originated in the common ancestor of primates and rodents. It subsequently underwent two remarkable expansions, first in the common ancestor of mice and rats and then in mice. Positive selection promoting amino acid substitutions was detected in some young Rhox genes, suggesting adaptive functional diversification. The recent expansions of the Rhox cluster provide an opportunity to study the mechanism and origin of colinear gene regulation, but they may also undermine the utility of mouse models for understanding the development and physiology of the human reproductive system.     

Transcriptome analysis reveals the difference between “healthy” and “common” aging and their connection with age‐related diseases

A key goal of aging research was to understand mechanisms underlying healthy aging and develop methods to promote the human healthspan. One approach is to identify gene regulations unique to healthy aging compared with aging in the general population (i.e., “common” aging). Here, we leveraged Genotype‐Tissue Expression (GTEx) project data to investigate “healthy” and “common” aging gene expression regulations at a tissue level in humans and their interconnection with diseases. Using GTEx donors' disease annotations, we defined a “healthy” aging cohort for each tissue. We then compared the age‐associated genes derived from this cohort with age‐associated genes from the “common” aging cohort which included all GTEx donors; we also compared the “healthy” and “common” aging gene expressions with various disease‐associated gene expressions to elucidate the relationships among “healthy,” “common” aging and disease. Our analyses showed that 1. GTEx “healthy” and “common” aging shared a large number of gene regulations; 2. Despite the substantial commonality, “healthy” and “common” aging genes also showed distinct function enrichment, and “common” aging genes had a higher enrichment for disease genes; 3. Disease‐associated gene regulations were overall different from aging gene regulations. However, for genes regulated by both, their regulation directions were largely consistent, implying some aging processes could increase the susceptibility to disease development; and 4. Possible protective mechanisms were associated with some “healthy” aging gene regulations. In summary, our work highlights several unique features of GTEx “healthy” aging program. This new knowledge could potentially be used to develop interventions to promote the human healthspan.     

Natural selection on gene expression

Changes in genetic regulation contribute to adaptations in natural populations and influence susceptibility to human diseases. Despite their potential phenotypic importance, the selective pressures acting on regulatory processes in general and gene expression levels in particular are largely unknown. Studies in model organisms suggest that the expression levels of most genes evolve under stabilizing selection, although a few are consistent with adaptive evolution. However, it has been proposed that gene expression levels in primates evolve largely in the absence of selective constraints. In this article, we discuss the microarray-based observations that led to these disparate interpretations. We conclude that in both primates and model organisms, stabilizing selection is likely to be the dominant mode of gene expression evolution. An important implication is that mutations affecting gene expression will often be deleterious and might underlie many human diseases.     

The Rates of Molecular Evolution in Rodent and Primate Mitochondrial DNA

A higher rate of molecular evolution in rodents than in primates at synonymous sites and, to a lesser extent, at amino acid replacement sites has been reported previously for most nuclear genes examined. Thus in these genes the average ratio of amino acid replacement to synonymous substitution rates in rodents is lower than in primates, an observation at odds with the neutral model of molecular evolution. Under Ohta's mildly deleterious model of molecular evolution, these observations are seen as the consequence of the combined effects of a shorter generation time (driving a higher mutation rate) and a larger effective population size (resulting in more effective selection against mildly deleterious mutations) in rodents. The present study reports the results of a maximum-likelihood analysis of the ratio of amino acid replacements to synonymous substitutions for genes encoded in mitochondrial DNA (mtDNA) in these two lineages. A similar pattern is observed: in rodents this ratio is significantly lower than in primates, again consistent only with the mildly deleterious model. Interestingly the lineage-specific difference is much more pronounced in mtDNA-encoded than in nuclear-encoded proteins, an observation which is shown to run counter to expectation under Ohta's model. Finally, accepting certain fossil divergence dates, the lineage-specific difference in amino acid replacement-to-synonymous substitution ratio in mtDNA can be partitioned and is found to be entirely the consequence of a higher mutation rate in rodents. This conclusion is consistent with a replication-dependent model of mutation in mtDNA. Received: 24 September 1999 / Accepted: 18 September 2000     

Epigenetic programming of differential gene expression in development and evolution

This review covers data on changing patterns of DNA methylation and the regulation of gene expression in mouse embryonic development. Global demethylation occurs from the eight-cell stage to the blastocyst stage in pre-implantation embryos, and global de novo methylation begins at implantation. We have used X-chromosome inactivation in female embryos as a model system to study specific CpG sites in the X-linked Pgk-1 and Gópd housekeeping genes and in the imprinted regulatory Xist gene to elucidate the role of methylation in the initiation and maintenance of differential gene activity. Methyl-ation of the X-linked housekeeping genes occurs very close in time to their inactivation, thus raising the question as to whether methylation could be causal to inactivation, as well as being involved in its maintenance. A methylation difference between sperm and eggs in the promoter region of the Xist gene, located at the X-chromosome inactivation centre, is correlated with imprinted preferential inactivation of the paternal X chromosome in extra-embryonic tissues. Based on our data, a picture of the inheritance of methylation imprints and speculation on the significance of the Xist imprint in development is presented. On a more general level, an hypothesis of evolution by “adaptive epige-netic/genetic inheritance” is considered. This proposes modification of germ line DNA in response to a change in environment and mutation at the site of modification (e.g., of methylated cytosine to thymine). Epigenetic inheritance could function to shift patterns of gene expression to buffer the evolving system against changes in environment. If the altered patterns of gene activity and inactivity persist, the modifications may become “fixed” as mutations; alternatively, previously silenced gene networks might be recruited into function, thus appearing as if they are “acquired characteristics.” An extension of this hypothesis is “foreign gene acquisition and sorting” (selection or silencing of gene function according to use). “Kidnapping” and sorting of foreign genes in this way could explain the observation that increased complexity in evolution is associated with more “junk” DNA. Adaptive epigenetic/genetic inheritance challenges the “central dogma” that information is unidirectional from the DNA to protein and the idea that Darwinian random mutation and selection are the sole mechanisms of evolution. © 1995 Wiley-Liss, Inc.     

Sequence variation at two eosinophil-associated ribonuclease loci in humans

Host defense against invading pathogens is of great importance to the survival of higher organisms. We have been studying the evolution of mammalian eosinophil-associated ribonucleases (EARs), which are members of the ribonuclease A superfamily with known antipathogen activities. Earlier studies showed that positive selection promoted rapid diversification of paralogous EAR genes in both primates and rodents. Intraspecifically, however, it is unknown whether these genes also have divergent alleles. The recent discovery that the gene repertoire of the EAR family is much larger in rodents than in primates has led us to consider the possibility that primates maintain a large number of polymorphic alleles to compensate for a smaller gene repertoire. Here we present sequences of 2417 nucleotides at the two EAR loci, the eosinophil-derived neurotoxin (EDN, RNase 2) and eosinophil cationic protein (ECP, RNase 3), from >50 human individuals. Our data demonstrate that the nucleotide diversities (0.06-0.11%) at these loci are typical for human nuclear genes, thus permitting us to reject this polymorphism hypothesis. No significant departure from neutrality is noted and no signs of overdominant selection are observed. Similar patterns were observed in a preliminary study of chimpanzees. In summary, our results suggest that the antipathogen functions of the primate EARs are conserved after they are established and that these proteins are not currently undergoing rapid diversification in response to challenge from invading microorganisms.     

Differential Evolution of MAGE Genes Based on Expression Pattern and Selection Pressure

Starting from publicly-accessible datasets, we have utilized comparative and phylogenetic genome analyses to characterize the evolution of the human MAGE gene family. Our characterization of genomic structures in representative genomes of primates, rodents, carnivora, and macroscelidea indicates that both Type I and Type II MAGE genes have undergone lineage-specific evolution. The restricted expression pattern in germ cells of Type I MAGE orthologs is observed throughout evolutionary history. Unlike Type II MAGEs that have conserved promoter sequences, Type I MAGEs lack promoter conservation, suggesting that epigenetic regulation is a central mechanism for controlling their expression. Codon analysis shows that Type I but not Type II MAGE genes have been under positive selection. The combination of genomic and expression analysis suggests that Type 1 MAGE promoters and genes continue to evolve in the hominin lineage, perhaps towards functional diversification or acquiring additional specific functions, and that selection pressure at codon level is associated with expression spectrum.     

Evolution of Brain Active Gene Promoters in Human Lineage Towards the Increased Plasticity of Gene Regulation

Adaptability to a variety of environmental conditions is a prominent feature of Homo sapiens. We hypothesize that this feature can be explained by evolutionary changes in gene promoters active in the brain prefrontal cortex leading to a more flexible gene regulation network. The genotype-dependent range of gene expression can be broader in humans than in other higher primates. Thus, we searched for specific signatures of evolutionary changes in promoter architectures of multiple hominid genes, including the genes active in human cortical neurons that may indicate an increase of variability of gene expression rather than just changes in the level of expression, such as downregulation or upregulation of the genes. We performed a whole-genome search for genetic-based alterations that may impact gene regulation “flexibility” in a process of hominids evolution, such as (i) CpG dinucleotide content, (ii) predicted nucleosome-DNA dissociation constant, and (iii) predicted affinities for TATA-binding protein (TBP) in gene promoters. We tested all putative promoter regions across the human genome and especially gene promoters in active chromatin state in neurons of prefrontal cortex, the brain region critical for abstract thinking and social and behavioral adaptation. Our data imply that the origin of modern man has been associated with an increase of flexibility of promoter-driven gene regulation in brain. In contrast, after splitting from the ancestral lineages of H. sapiens, the evolution of ape species is characterized by reduced flexibility of gene promoter functioning, underlying reduced variability of the gene expression.     

Widespread positive selection in synonymous sites of mammalian genes

Evolution of protein sequences is largely governed by purifying selection, with a small fraction of proteins evolving under positive selection. The evolution at synonymous positions in protein-coding genes is not nearly as well understood, with the extent and types of selection remaining, largely, unclear. A statistical test to identify purifying and positive selection at synonymous sites in protein-coding genes was developed. The method compares the rate of evolution at synonymous sites (Ks) to that in intron sequences of the same gene after sampling the aligned intron sequences to mimic the statistical properties of coding sequences. We detected purifying selection at synonymous sites in approximately 28% of the 1,562 analyzed orthologous genes from mouse and rat, and positive selection in approximately 12% of the genes. Thus, the fraction of genes with readily detectable positive selection at synonymous sites is much greater than the fraction of genes with comparable positive selection at nonsynonymous sites, i.e., at the level of the protein sequence. Unlike other genes, the genes with positive selection at synonymous sites showed no correlation between Ks and the rate of evolution in nonsynonymous sites (Ka), indicating that evolution of synonymous sites under positive selection is decoupled from protein evolution. The genes with purifying selection at synonymous sites showed significant anticorrelation between Ks and expression level and breadth, indicating that highly expressed genes evolve slowly. The genes with positive selection at synonymous sites showed the opposite trend, i.e., highly expressed genes had, on average, higher Ks. For the genes with positive selection at synonymous sites, a significantly lower mRNA stability is predicted compared to the genes with negative selection. Thus, mRNA destabilization could be an important factor driving positive selection in nonsynonymous sites, probably, through regulation of expression at the level of mRNA degradation and, possibly, also translation rate. So, unexpectedly, we found that positive selection at synonymous sites of mammalian genes is substantially more common than positive selection at the level of protein sequences. Positive selection at synonymous sites might act through mRNA destabilization affecting mRNA levels and translation.     

The genomic context of retrocopies increases their chance of functional relevancy in mammals

Described as “junk” DNA, pseudogenes are dead structures of previously active genes present in genomes. Pseudogenes are categorized into two main classes: processed pseudogenes, formed through retrotransposition, and non-processed pseudogenes, typically originated from gene decay following duplication events. The term “processed pseudogene” has changed to “retrocopy” since they are likely to evolve new functional roles and became a retrogene. Here, we surveyed 38,080 retrocopies from chimpanzee, dog, human, mouse, and rat genomes to assess their potential adaptive value. The retrocopies inserted in the same chromosome of the parental gene have higher chances of remain potentially “active” (absence of premature stop codons and frameshifts) (~26.1%), while those placed into a different chromosome have a twofold decrease chance of continuing potentially “active” (~7.52%). The genomic context of their placement seems associated with their expression. Retrocopies placed in intragenic regions and the same sense of the “host” gene have higher chances of being expressed relative to other genomic contexts. The proximity of retrocopies to their parental gene is associated with a lower decay rate, and their location likely influence their expression. Thus, despite their unclear role, retrocopies are probably involved in adaptive processes. Our results evidence natural selection acting in retrocopies.     

Molecular evolution of candidate sour taste receptor gene PKD1L3 in mammals

The PKD1L3 gene encodes an ion channel protein that can interact with the PKD2L1 protein to form a candidate sour taste receptor. In the present study, we have analyzed the evolutionary patterns of PKD1L3 genes from 10 mammalian species. The results showed that PKD1L3 genes have evolved under a dominant purifying selection force. However, for some branches and sites, PKD1L3 genes were detected to have been operated by positive selection. Moreover, some of these positive evolutionary sites are likely to participate in acid stimulus recognition. In rodents, PKD1L3 genes evolved more rapidly than other mammalian lineages. Combined with other functional research reports, our results suggest that rodents may not be the most appropriate model for functional research on the PKD1L3 gene.     

DIFFERENCES IN THE REGULATION OF GROWTH AND BIOMINERALIZATION GENES REVEALED THROUGH LONG‐TERM COMMON‐GARDEN ACCLIMATION AND EXPERIMENTAL GENOMICS IN THE PURPLE SEA URCHIN

Across heterogeneous landscapes, populations may have adaptive differences in gene regulation that adjust their physiologies to match local environments. Such differences could have origins in acclimation or in genetically fixed variation between habitats. Here we use common‐garden experiments to evaluate differences in gene expression between populations of the purple sea urchin, Strongylocentrotus purpuratus, spanning 1700 km and average temperature differences of 5°C to 8°C. Across expression profiles from 18,883 genes after 3 years of common conditions, we find highly correlated expression patterns (Pearson's r = 0.992) among most genes. However, 66 genes were differentially expressed, including many ribosomal protein and biomineralization genes, which had higher expression in urchins originally from the southern population. Gene function analyses revealed slight but pervasive expression differences in genes related to ribosomal function, metabolism, transport, “bone” development, and response to stimuli. In accord with gene expression patterns, a post‐hoc spine regrowth experiment revealed that urchins of southern origin regrew spines at a faster rate than northern urchins. These results suggest that there may be genetically controlled, potentially adaptive differences in gene regulation across habitats and that gene expression differences may be under strong enough selection to overcome high, dispersal–mediated gene flow in this marine species.