Search

In an occasional series of articles, we will be publishing autobiographical sketches from some of those working in the field of genome evolution. The series will feature both the very eminent, but also researchers closer to their start of their career, and those from underrepresented groups. The series will show the unusual paths that academics sometimes take and the obstacles they have overcome. We start this series with one of the most influential researchers in the field of molecular evolution, Wen-Hsiung Li. Wen-Hsiung has contributed enormously to the field and published on a wide diversity of topics, as described in this autobiographical sketch; he also wrote two textbooks, one of them with Dan Graur, which for many years were the bibles of the field. He was awarded the Motoo Kimura prize by the Society of Molecular Biology and Evolution in 2019 in recognition of his contributions to our subject., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

How negative selection, positive selection, and population size contribute to the large variation in nucleotide substitution rates among RNA viruses remains unclear. Here, we studied the ratios of nonsynonymous-to-synonymous substitution rates ( d N / d S ) in protein-coding genes of human RNA and DNA viruses and mammals. Among the 21 RNA viruses studied, 18 showed a genome-average d N / d S from 0.01 to 0.10, indicating that over 90% of nonsynonymous mutations are eliminated by negative selection. Only HIV-1 showed a d N / d S (0.31) higher than that (0.22) in mammalian genes. By comparing the d N / d S values among genes in the same genome and among species or strains, we found that both positive selection and population size play significant roles in the d N / d S variation among genes and species. Indeed, even in flaviviruses and picornaviruses, which showed the lowest ratios among the 21 species studied, positive selection appears to have contributed significantly to d N / d S We found the view that positive selection occurs much more frequently in influenza A subtype H3N2 than subtype H1N1 holds only for the hemagglutinin and neuraminidase genes, but not for other genes. Moreover, we found no support for the view that vector-borne RNA viruses have lower d N / d S ratios than non-vector-borne viruses. In addition, we found a correlation between d N and d S , implying a correlation between d N and the mutation rate. Interestingly, only 2 of the 8 DNA viruses studied showed a d N / d S < 0.10, while 4 showed a d N / d S > 0.22. These observations increase our understanding of the mechanisms of RNA virus evolution., Competing Interests: The authors declare no conflict of interest.

We studied the evolution of opsin genes in 59 ray-finned fish genomes. We identified the opsin genes and adjacent genes (syntenies) in each genome. Then we inferred the changes in gene copy number (N), syntenies, and tuning sites along each phylogenetic branch during evolution. The Exorh (rod opsin) gene has been retained in 56 genomes. Rh1, the intronless rod opsin gene, first emerged in ancestral Actinopterygii, and N increased to 2 by the teleost-specific whole genome duplication, but then decreased to 1 in the ancestor of Neoteleostei fishes. For cone opsin genes, the rhodopsin-like (Rh2) and long-wave-sensitive (LWS) genes showed great variation in N among species, ranging from 0 to 5 and from 0 to 4, respectively. The two short-wave-sensitive genes, SWS1 and SWS2, were lost in 23 and 6 species, respectively. The syntenies involving LWS, SWS2 and Rh2 underwent complex changes, while the evolution of the other opsin gene syntenies was much simpler. Evolutionary adaptation in tuning sites under different living environments was discussed. Our study provides a detailed view of opsin gene gains and losses, synteny changes and tuning site changes during ray-finned fish evolution.

Background: Long non-coding RNAs (lncRNAs) are important in various biological processes, but very few studies on lncRNA have been conducted in birds. To identify IncRNAs expressed during feather development, we analyzed single-stranded RNA-seq (ssRNA-seq) data from the anterior and posterior dorsal regions during zebra finch (Taeniopygia guttata) embryonic development. Using published transcriptomic data, we further analyzed the evolutionary conservation of IncRNAs in birds and amniotes., Results: A total of 1,081 lncRNAs, including 965 intergenic lncRNAs (lincRNAs), 59 intronic lncRNAs, and 57 antisense lncRNAs (lncNATs), were identified using our newly developed pipeline. These avian IncRNAs share similar characteristics with lncRNAs in mammals, such as shorter transcript length, lower exon number, lower average expression level and less sequence conservation than mRNAs. However, the proportion of lncRNAs overlapping with transposable elements in birds is much lower than that in mammals. We predicted the functions of IncRNAs based on the enriched functions of co-expressed protein-coding genes. Clusters of lncRNAs associated with natal down development were identified. The sequences and expression levels of candidate lncRNAs that shared conserved sequences among birds were validated by qPCR in both zebra finch and chicken. Finally, we identified three highly conserved lncRNAs that may be associated with natal down development., Conclusions: Our study provides the first systematical identification of avian lncRNAs using ssRNA-seq analysis and offers a resource of embryonically expressed lncRNAs in zebra finch. We also predicted the biological function of identified lncRNAs.

Gene expression evolution occurs through changes in cis- or trans-regulatory elements or both. Interactions between transcription factors (TFs) and their binding sites (TFBSs) constitute one of the most important points where these two regulatory components intersect. In this study, we investigated the evolution of TFBSs in the promoter regions of different Saccharomyces strains and species. We divided the promoter of a gene into the proximal region and the distal region, which are defined, respectively, as the 200-bp region upstream of the transcription starting site and as the 200-bp region upstream of the proximal region. We found that the predicted TFBSs in the proximal promoter regions tend to be evolutionarily more conserved than those in the distal promoter regions. Additionally, Saccharomyces cerevisiae strains used in the fermentation of alcoholic drinks have experienced more TFBS losses than gains compared with strains from other environments (wild strains, laboratory strains, and clinical strains). We also showed that differences in TFBSs correlate with the cis component of gene expression evolution between species (comparing S. cerevisiae and its sister species Saccharomyces paradoxus) and within species (comparing two closely related S. cerevisiae strains)., (© The Author(s) 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The sequences of the untranslated regions (UTRs) of mRNAs play important roles in posttranscriptional regulation, but whether a change in UTR length can significantly affect the regulation of gene expression is not clear. In this study, we examined the connection between UTR length and Expression Correlation with cytosolic ribosomal proteins (CRP) genes (ECC), which measures the level of expression similarity of a group of genes with CRP genes under various growth conditions. We used data from the aerobic fermentation yeast Saccharomyces cerevisiae and the aerobic respiration yeast Candida albicans. To reduce statistical fluctuations, we computed the ECC for the genes in a Gene Ontology (GO) functional group. We found that in both species, ECC is strongly correlated with the 5' UTR length but not with the 3' UTR length and that the 5' UTR length is evolutionarily better conserved than the 3' UTR length. Interestingly, we found 11 GO groups that have had a substantial increase in 5' UTR length in the S. cerevisiae lineage and that the length increase was associated with a substantial decrease in ECC. Moreover, 9 of the 11 GO groups of genes are involved in mitochondrial respiration function, whose expression reprogramming has been shown to be a major factor for the evolution of aerobic fermentation. Finally, we found that an increase in 5' UTR length may decrease the +1 nucleosome occupancy. This study provides a new angle to understand the role of 5' UTR in gene expression regulation and evolution.

Many indicators of protein evolutionary rate have been proposed, but some of them are interrelated. The purpose of this study is to disentangle their correlations. We assess the strength of each indicator by controlling for the other indicators under study. We find that the number of microRNA (miRNA) types that regulate a gene is the strongest rate indicator (a negative correlation), followed by disorder content (the percentage of disordered regions in a protein, a positive correlation); the strength of disorder content as a rate indicator is substantially increased after controlling for the number of miRNA types. By dividing proteins into lowly and highly intrinsically disordered proteins (L-IDPs and H-IDPs), we find that proteins interacting with more H-IDPs tend to evolve more slowly, which largely explains the previous observation of a negative correlation between the number of protein-protein interactions and evolutionary rate. Moreover, all of the indicators examined here, except for the number of miRNA types, have different strengths in L-IDPs and in H-IDPs. Finally, the number of phosphorylation sites is weakly correlated with the number of miRNA types, and its strength as a rate indicator is substantially reduced when other indicators are considered. Our study reveals the relative strength of each rate indicator and increases our understanding of protein evolution.

It is well known that knocking out a gene in an organism often causes no phenotypic effect. One possible explanation is the existence of duplicate genes; that is, the effect of knocking out a gene is compensated by a duplicate copy. Another explanation is the existence of alternative pathways. In terms of metabolic products, the relative roles of the two mechanisms have been extensively studied in yeast but not in any multi-cellular organisms. Here, to address the functional compensation of metabolic products by duplicate genes, we quantified 35 metabolic products from 1,976 genes in knockout mutants of Arabidopsis thaliana by a high-throughput Liquid chromatography-Mass spectrometer (LC-MS) analysis. We found that knocking out either a singleton gene or a duplicate gene with distant paralogs in the genome tends to induce stronger metabolic effects than knocking out a duplicate gene with a close paralog in the genome, indicating that only duplicate genes with close paralogs play a significant role in functional compensation for metabolic products in A. thaliana. To extend the analysis, we examined metabolic products with either high or low connectivity in a metabolic network. We found that the compensatory role of duplicate genes is less important when the metabolite has a high connectivity, indicating that functional compensation by alternative pathways is common in the case of high connectivity. In conclusion, recently duplicated genes play an important role in the compensation of metabolic products only when the number of alternative pathways is small.

The genetic basis of organisms' adaptation to different environments is a central issue of molecular evolution. The budding yeast Saccharomyces cerevisiae and its relatives predominantly ferment glucose into ethanol even in the presence of oxygen. This was suggested to be an adaptation to glucose-rich habitats, but the underlying genetic basis of the evolution of aerobic fermentation remains unclear. In S. cerevisiae, the first step of glucose metabolism is transporting glucose across the plasma membrane, which is carried out by hexose transporter proteins. Although several studies have recognized that the rate of glucose uptake can affect how glucose is metabolized, the role of HXT genes in the evolution of aerobic fermentation has not been fully explored. In this study, we identified all members of the HXT gene family in 23 fully sequenced fungal genomes, reconstructed their evolutionary history to pinpoint gene gain and loss events, and evaluated their adaptive significance in the evolution of aerobic fermentation. We found that the HXT genes have been extensively amplified in the two fungal lineages that have independently evolved aerobic fermentation. In contrast, reduction of the number of HXT genes has occurred in aerobic respiratory species. Our study reveals a strong positive correlation between the copy number of HXT genes and the strength of aerobic fermentation, suggesting that HXT gene expansion has facilitated the evolution of aerobic fermentation.

Protein phosphorylation plays an important role in the regulation of protein function. Phosphorylated residues are generally assumed to be subject to functional constraint, but it has recently been suggested from a comparison of distantly related vertebrate species that most phosphorylated residues evolve at the rates consistent with the surrounding regions. To resolve the controversy, we infer the ancestral phosphoproteome of human and mouse to compare the evolutionary rates of phosphorylated and nonphosphorylated serine (S), threonine (T), and tyrosine (Y) residues. This approach enables accurate estimation of evolutionary rates as it does not assume deep conservation of phosphorylated residues. We show that phosphorylated S/T residues tend to evolve more slowly than nonphosphorylated S/T residues not only in disordered but also in ordered protein regions, indicating evolutionary conservation of phosphorylated S/T residues in mammals. Thus, phosphorylated S/T residues tend to be subject to stronger functional constraint than nonphosphorylated residues regardless of the protein regions in which they reside. In contrast, phosphorylated Y residues evolve at similar rates as nonphosphorylated ones. We also find that the human lineage has gained more phosphorylated T residues and lost fewer phosphorylated Y residues than the mouse lineage. The cause of the gain/loss imbalance remains a mystery but should be worth exploring.

The regulation of gene expression is an important determinant of organismal phenotype and evolution. However, the widespread recognition of this fact occurred long after the synthesis of evolution and genetics. Here, we give a brief sketch of thoughts regarding gene regulation in the history of evolution and genetics. We then review the development of genome-wide studies of gene regulatory variation in the context of the location and mode of action of the causative genetic changes. In particular, we review mapping of the genetic basis of expression variation through expression quantitative trait locus studies and measuring the cis/trans component of expression variation in allele-specific expression studies. We conclude by proposing a systematic integration of ideas that combines global mapping studies, cis/trans tests and modern population genetics methodologies, in order to directly estimate the forces acting on regulatory variation within and between species.

The evolution of duplicate genes has been a topic of broad interest. Here, we propose that the conservation of gene family size is a good indicator of the rate of sequence evolution and some other biological properties. By comparing the human-chimpanzee-macaque orthologous gene families with and without family size conservation, we demonstrate that genes with family size conservation evolve more slowly than those without family size conservation. Our results further demonstrate that both family expansion and contraction events may accelerate gene evolution, resulting in elevated evolutionary rates in the genes without family size conservation. In addition, we show that the duplicate genes with family size conservation evolve significantly more slowly than those without family size conservation. Interestingly, the median evolutionary rate of singletons falls in between those of the above two types of duplicate gene families. Our results thus suggest that the controversy on whether duplicate genes evolve more slowly than singletons can be resolved when family size conservation is taken into consideration. Furthermore, we also observe that duplicate genes with family size conservation have the highest level of gene expression/expression breadth, the highest proportion of essential genes, and the lowest gene compactness, followed by singletons and then by duplicate genes without family size conservation. Such a trend accords well with our observations of evolutionary rates. Our results thus point to the importance of family size conservation in the evolution of duplicate genes.

It is commonly assumed but not proven that microRNAs (miRNAs) and their targets coevolve. Under this assumption, miRNAs and targets from different species may interact adversely, resulting in reduced fitness. However, the strength of the adverse interactions may not be detectable because even outright deletions of miRNAs often manifest only subtle fitness effects. We tested and measured the strength of heterospecific interactions by carrying out transgenic experiments across Drosophila species by overexpressing the miR310s cluster of Drosophila melanogaster (Dm310s) and Drosophila pseudoobscura (Dp310s) in D. melanogaster. Flies overexpressing the heterospecific Dp310s are only one-third as viable as those overexpressing the conspecific Dm310s. The viability effect is easily detectable in comparison to the effect of the deletion of miR310s. The number of genes significantly misexpressed under the influence of Dp310s is 3-10 times greater than under Dm310s. Importantly, the numbers of predicted targets are similar between them. Expression analysis of the predicted target genes suggests that miRNAs may sometimes function to buffer fluctuations in the transcriptome output. After the buffering function has evolved, heterospecific combinations may cause adverse effects.

Both cis and trans mutations contribute to gene expression divergence within and between species. We used Saccharomyces cerevisiae as a model organism to estimate the relative contributions of cis and trans variations to the expression divergence between a laboratory (BY) and a wild (RM) strain of yeast. We examined whether genes regulated by a single transcription factor (TF; single input module, SIM genes) or genes regulated by multiple TFs (multiple input module, MIM genes) are more susceptible to trans variation. Because a SIM gene is regulated by a single immediate upstream TF, the chance for a change to occur in its trans-acting factors would, on average, be smaller than that for a MIM gene. We chose 232 genes that exhibited expression divergence between BY and RM to test this hypothesis. We examined the expression patterns of these genes in a BY-RM coculture system and in a BY-RM diploid hybrid. We found that trans variation is far more important than cis variation for expression divergence between the two strains. However, because in 75% of the genes studied, cis variation has significantly contributed to expression divergence, cis change also plays a significant role in intraspecific expression evolution. Interestingly, we found that the proportion of genes with diverged expression between BY and RM is larger for MIM genes than for SIM genes; in fact, the proportion tends to increase with the number of transcription factors that regulate the gene. Moreover, MIM genes are, on average, subject to stronger trans effects than SIM genes, though the difference between the two types of genes is not conspicuous.

There are now many known cases of orthologous or unrelated proteins in different species that have undergone parallel evolution to satisfy a similar function. However, there are no reported cases of parallel evolution for proteins that bind a common ligand but have different functions. We focused on two proteins that have different functions in steroid hormone biosynthesis and action but bind a common ligand, androgen. The first protein, androgen receptor (AR), is a nuclear hormone receptor and the second one, aromatase (cytochrome P450 19 [CYP19]), converts androgen to estrogen. We hypothesized that binding of the androgen ligand has exerted common selective pressure on both AR and CYP19, resulting in a signature of parallel evolution between these two proteins, though they perform different functions. Consistent with this hypothesis, we found that rates of amino acid change in AR and CYP19 are strongly correlated across the metazoan phylogeny, whereas no significant correlation was found in the control set of proteins. Moreover, we inferred that genomic toolkits required for steroid biosynthesis and action were present in a basal metazoan, cnidarians. The close similarities between vertebrate and sea anemone AR and CYP19 suggest a very ancient origin of their endocrine functions at the base of metazoan evolution. Finally, we found evidence supporting the hypothesis that the androgen-to-estrogen ratio determines the gonadal sex in all metazoans.

The enzymes catalyzing the first two reactions in the sulfate activation pathway, ATP-sulfurylase (S) and APS-kinase (K), are fused as 'KS' in animals but are fused as 'SK' in select bacteria and fungi. We have discovered a novel triple fusion protein of K, S, and pyrophosphatase (P) in several protozoan genomes within the Stramenopile lineage. These triple domain fusion proteins led us to hypothesize that pyrophosphatase enzymes and sulfate activation enzymes physically interact to impact the thermodynamics of the sulfate activation pathway. In support of this hypothesis, we demonstrate through biochemical assays that separately encoded KS and P proteins physically interact and that KS/P complexes activate more sulfate than KS alone. We also conclude on the basis of phylogenetic analyses that all known KS fusions originate from a single fusion event early in the eukaryotic lineage. Strikingly, our analyses support the same conclusion for all known SK fusions. These observations indicate that the patchwork of fused and nonfused S and K genes observed in modern-day eukaryotes and prokaryotes are the result of the two ancestral fusion genes evolving by an assortment of gene fissions, duplications, deletions, and horizontal transfers in different lineages. Our integrative use of genomics, phylogenetics, and biochemistry to characterize pyrophosphatase as a new member of the sulfate activation pathway should be effective at identifying new protein members and connections in other molecular pathways.

The current variance estimators for Jukes and Cantor's one-parameter model and Kimura's two-parameter model tend to underestimate the true variances when the true proportion of differences between the two sequences under study is not small. In this paper, we developed improved variance estimators, using a higher-order Taylor expansion and empirical methods. The new estimators outperform the conventional estimators and provide accurate estimates of the true variances.

The yeast Saccharomyces cerevisiae proliferates rapidly in glucose-containing media. As glucose is getting depleted, yeast cells enter the transition from fermentative to nonfermentative metabolism, known as the diauxic shift, which is associated with major changes in gene expression. To understand the expression evolution of genes involved in the diauxic shift and in nonfermentative metabolism within species, a laboratory strain (BY), a wild strain (RM), and a clinical isolate (YJM) were used in this study. Our data showed that the RM strain enters into the diauxic shift approximately 1 h earlier than the BY strain with an earlier, higher induction of many key transcription factors (TFs) involved in the diauxic shift. Our sequence data revealed sequence variations between BY and RM in both coding and promoter regions of the majority of these TFs. The key TF Cat8p, a zinc-finger cluster protein, is required for the expression of many genes in gluconeogenesis under nonfermentative growth, and its derepression is mediated by deactivation of Mig1p. Our kinetic study of CAT8 expression revealed that CAT8 induction corresponded to the timing of glucose depletion in both BY and RM and CAT8 was induced up to 50- to 90-folds in RM, whereas only 20- to 30-folds in BY. In order to decipher the relative importance of cis- and trans-variations in expression divergence in the gluconeogenic pathway during the diauxic shift, we studied the expression levels of MIG1, CAT8, and their downstream target genes in the cocultures and in the hybrid diploids of BY-RM, BY-YJM, and RM-YJM and in strains with swapped promoters. Our data showed that the differences between BY and RM in the expression of MIG1, the upstream regulator of CAT8, were affected mainly by changes in cis-elements, though also by changes in trans-acting factors, whereas those of CAT8 and its downstream target genes were predominantly affected by changes in trans-acting factors.

Despite much interest in regulatory evolution, how promoters have evolved remains poorly studied, mainly owing to paucity of data on promoter regions. Using a new set of high-quality experimentally determined core promoters of the human genome, we conducted a comparative analysis of 2,492 human and rhesus macaque promoters and their neighboring nearly neutral regions. We found that the core promoters have an average rate of nucleotide substitution substantially higher than that at 4-fold degenerate sites and only slightly lower than that for the assumed neutral controls of neighboring noncoding regions, suggesting that core promoters are subject to very weak selective constraints. Interestingly, we identified 24 core promoters (at false discovery rate = 50%) that have evolved at an accelerated rate compared with the neutral controls, suggesting that they may have undergone positive selection. The inferred positively selected genes show strong bias in molecular function. We also used population genetic approaches to examine the evolution of core promoters in human populations and found evidence of positive selection at some loci. Taken together, our results suggest that positive selection has played a substantial role in the evolution of transcriptional regulation in primates.

Background: Increasing evidence suggests that overlapping genes are much more common in eukaryotic genomes than previously thought. In this study we identified and characterized the overlapping genes in a set of 13,484 pairs of human-mouse orthologous genes., Results: About 10% of the genes under study are overlapping genes, the majority of which are different-strand overlaps. The majority of the same-strand overlaps are embedded forms, whereas most different-strand overlaps are not embedded and in the convergent transcription orientation. Most of the same-strand overlapping gene pairs show at least a tenfold difference in length, much larger than the length difference between non-overlapping neighboring gene pairs. The length difference between the two different-strand overlapping genes is less dramatic. Over 27% of the different-strand-overlap relationships are shared between human and mouse, compared to only approximately 8% conservation for same-strand-overlap relationships. More than 96% of the same-strand and different-strand overlaps that are not shared between human and mouse have both genes located on the same chromosomes in the species that does not show the overlap. We examined the causes of transition between the overlapping and non-overlapping states in the two species and found that 3' UTR change plays an important role in the transition., Conclusion: Our study contributes to the understanding of the evolutionary transition between overlapping genes and non-overlapping genes and demonstrates the high rates of evolutionary changes in the un-translated regions.

Where did vertebrate genes come from? Here we address this question by analyzing eight completely sequenced land vertebrate genomes and six completely sequenced invertebrate genomes. Approximately 70% of the vertebrate genes can be found in the six invertebrate genomes with the standard homology search criteria (denoted as V.MCL), another approximately 6% can be found with relaxed search criteria, and an additional approximately 2% can be found in sequenced fungal and bacterial genomes. Thus, a substantial proportion of vertebrate genes (approximately 22%) cannot be found in the nonvertebrate genomes studied (denoted as Vonly). Interestingly, genes in Vonly are predominantly singletons, while the majority of genes in the other three groups belong to gene families. The proteins of Vonly tend to evolve faster than those of V.MCL. Surprisingly, in many cases the family sizes in V.MCL are only as large as or even smaller than their counterparts in the invertebrates, contrary to the general perception of a larger family size in vertebrates. Interestingly, in comparison with the family size in invertebrates, vertebrate gene families involved in regulation, signal transduction, transcription, protein transport, and protein modification tend to be expanded, whereas those involved in metabolic processes tend to be contracted. Furthermore, for almost all of the functional categories with family size expansion in vertebrates, the number of gene types (i.e., the number of singletons plus the number of gene families) tends to be over-represented in Vonly, but under-represented in V.MCL. Our study suggests that gene function is a major determinant of gene family size.

There are 2 ways to infer selection pressures in the evolution of protein-coding genes, the nonsynonymous and synonymous substitution rate ratio (K(A)/K(S)) and the radical and conservative amino acid replacement rate ratio (K(R)/K(C)). Because the K(R)/K(C) ratio depends on the definition of radical and conservative changes in the classification of amino acids, we develop an amino acid classification that maximizes the correlation between K(A)/K(S) and K(R)/K(C). An analysis of 3,375 orthologous gene groups among 5 mammalian species shows that our classification gives a significantly higher correlation coefficient between the 2 ratios than those of existing classifications. However, there are many orthologous gene groups with a low K(A)/K(S) but a high K(R)/K(C) ratio. Examining the functions of these genes, we found an overrepresentation of functional categories related to development. To determine if the overrepresentation is stage specific, we examined the expression patterns of these genes at different developmental stages of the mouse. Interestingly, these genes are highly expressed in the early middle stage of development (blastocyst to amnion). It is commonly thought that developmental genes tend to be conservative in evolution, but some molecular changes in developmental stages should have contributed to morphological divergence in adult mammals. Therefore, we propose that the relaxed pressures indicated by the K(R)/K(C) ratio but not by K(A)/K(S) in the early middle stage of development may be important for the morphological divergence of mammals at the adult stage, whereas purifying selection detected by K(A)/K(S) occurs in the early middle developmental stage.

It has previously been found that, in yeast, gene essentiality is positively correlated with protein connectivity (number of interaction partners) but negatively correlated with the existence of gene duplicates and that highly connected proteins tend to have a low gene duplicability. Using data from human and mouse, we show here that, in mammals, the first of these relationships holds true, but unlike the second relationship in yeast, highly connected mammalian proteins tend to have a high gene duplicability, and there is no correlation between gene essentiality and gene duplication in mammals.

Both cis- and trans-regulatory mutations contribute to gene expression divergence within and between species. To estimate their relative contributions, we examined two yeast strains, BY (a laboratory strain) and RM (a wild strain), for their gene-expression divergence by microarray. Using these data and published ChIP-chip data, we obtained a set of single-regulator-regulated genes that showed expression divergence between BY and RM. We randomly selected 50 of these genes for further study. We developed a step-by-step approach to assess the relative contributions of cis- and trans-variations to expression divergence by using pyrosequencing to quantify the mRNA levels of the BY and RM alleles in the same culture (co-culture) and in hybrid diploids. Forty genes showed expression divergence between the two strains in co-culture, and pyrosequencing of the BY/RM hybrid diploids showed that 45% (18/40) can be attributed to differences in trans-acting factors alone, 17.5% (7/40) mainly to trans-variations, 20% (8/40) to both cis- and trans-acting factors, 7.5% (3/40) mainly to cis-variations, and 10% (4/40) to cis-acting factors alone. In addition, we replaced the BY promoter by the RM promoter in each of 10 BY genes that were found from our microarray data to have expression divergence between BY and RM, and in each case our quantitative PCR analysis revealed a cis effect of the promoter replacement on gene expression. In summary, our study suggests that trans-acting factors play the major role in expression evolution between yeast strains, but the role of cis variation is also important.

Internal duplication can enhance the function of a gene or provide raw material for the emergence of a new function in a gene. Therefore, it is interesting to see whether the frequency of internal duplication has increased during metazoan evolution. The growing number of sequenced eukaryotic genomes provides an excellent opportunity to study the change in the pattern of internal duplication in the course of metazoan evolution. We studied repeated segments in proteins in the proteomes of 11 eukaryotes. We found that the frequency of internal duplication in Caenorhabditis elegans and Drosophila melanogaster (two protostomes) is higher than that in fungi but lower than that in chordates. Moreover, the frequencies of internal duplication for the chordates studied are largely similar. We classified orthologous proteins of chordates into three antiquity groups and found that more recently derived proteins in the metazoan lineage have higher repetitiveness than older ones. Our analysis suggests that lineage-specific internal duplication in protein evolution increases with organismal complexity before the emergence of chordates but not so afterward. Proteins with repeated regions might have been preferred before the protostome-chordate split. This finding supports the suggestion that exon-shuffling occurred more frequently after the first multicellular organism appeared and might have contributed to the metazoan radiation.

The HA1 domain of HA, the major antigenic protein of influenza A viruses, contains all of the antigenic sites of HA and is under continual immune-driven selection. To resolve controversies on whether only a few or many residue sites of HA1 have undergone positive selection, whether positive selection at HA1 is continual or punctuated, and whether antigenic change is punctuated, we introduce an approach to analyze 2,248 HA1 sequences collected from 1968 to 2005. We identify 95 substitutions at 63 sites from 1968 to 2005 and show that each substitution occurred very rapidly. The rapid substitution and the fact that 57 of the 63 sites are antigenic sites indicate that hitchhiking plays a minor role and that most of these sites, many more than previously found, have undergone positive selection. Strikingly, 88 of the 95 substitutions occurred in groups, and multiple mutations at antigenic sites sped up the fixation process. Our results suggest that positive selection has been ongoing most of the time, not sporadic, and that multiple mutations at antigenic sites cumulatively enhance antigenic drift, indicating that antigenic change is less punctuated than recently proposed.

Translational selection, including gene expression, protein abundance, and codon usage bias, has been suggested as the single dominant determinant of protein evolutionary rate in yeast. Here, we show that protein structure is also an important determinant. Buried residues, which are responsible for maintaining protein structure or are located on a stable interaction surface between 2 subunits, are usually under stronger evolutionary constraints than solvent-exposed residues. Our partial correlation analysis shows that, when whole proteins are included, the variance of evolutionary rate explained by the proportion of solvent-exposed residues (P(exposed)) can reach two-thirds of that explained by translational selection, indicating that P(exposed) is the most important determinant of protein evolutionary rate next only to translational selection. Our result suggests that proteins with many residues under selective constraint (e.g., maintaining structure or intermolecular interaction) tend to evolve slowly, supporting the "fitness (functional) density" hypothesis.

It has been suggested that insertions and deletions (indels) have contributed to the sequence divergence between the human and chimpanzee genomes more than do nucleotide changes (3% vs. 1.2%). However, although there have been studies of large indels between the two genomes, no systematic analysis of small indels (i.e., indels 840,000 "small" indels, which affect >7000 UCSC-annotated human genes (>11,000 transcripts). These indels, however, amount to only approximately 0.21% sequence change in the human lineage for the regions compared, whereas in pseudogenes indels contribute to a sequence divergence of 1.40%, suggesting that most of the indels that occurred in genic regions have been eliminated. Functional analysis reveals that the genes whose coding exons have been affected by human-specific indels are enriched in transcription and translation regulatory activities but are underrepresented in catalytic and transporter activities, cellular and physiological processes, and extracellular region/matrix. This functional bias suggests that human-specific indels might have contributed to human unique traits by causing changes at the RNA and protein level.

To detect positive selection in protein-coding sequence evolution, the ratio of the nonsynonymous to synonymous substitution rate (K(A)/K(S)) is commonly used. When this ratio is higher than 1, positive selection on nonsynonymous changes is considered to have occurred. However, the question of what kinds of amino acid change are likely to be involved in positive selection has not been well studied, though intuitively it seems that radical changes frequently occur in positively selected changes. To address this question, we examined chemically radical and conservative replacements in the evolution of hepatitis C virus (HCV) protein sequences. In the envelope region, 34 positively and 440 negatively selected sites were identified by the K(A)/K(S) ratio. Radical and conservative changes were compared between the two types of selected sites using two methods. First, the numbers of radical and conservative replacements were counted at the positively and negatively selected sites according to three kinds of chemical classifications. In all three classifications, the resulting ratios of the two numbers were not statistically different for the two types of selected sites (P>0.05). Second, the distribution of chemical changes was compared between the two types of selected sites using two kinds of chemical distances. The distributions of the two chemical distances were not statistically different for the two types of selected sites (P>0.05). These results indicate that the ratio of chemically radical and conservative changes is similar for positively and negatively selected sites in the envelope protein of HCV or, in other words, there is no correlation between radical change and positive selection in the evolution of this protein.

Myostatin (GDF8) is a negative regulator of muscle growth in mammals, and loss-of-function mutations are associated with increased skeletal-muscle mass in mice, cattle, and humans. Here, we show that positive natural selection has acted on human nucleotide variation at GDF8, since the observed ratio of nonsynonymous:synonymous changes among humans is significantly greater than expected under the neutral model and is strikingly different from patterns observed across mammalian orders. Furthermore, extended haplotypes around GDF8 suggest that two amino acid variants have been subject to recent positive selection. Both mutations are rare among non-Africans yet are at frequencies of up to 31% in sub-Saharan Africans. These signatures of selection at the molecular level suggest that human variation at GDF8 is associated with functional differences.

There has been a controversy on whether alternatively spliced exons (ASEs) evolve faster than constitutively spliced exons (CSEs). Although it has been noted that ASEs are subject to weaker selective constraints than CSEs, so they evolve faster, there have also been studies that indicated slower evolution in ASEs than in CSEs. In this study, we retrieve more than 5,000 human-mouse orthologous exons and calculate the synonymous (KS) and nonsynonymous (KA) substitution rates in these exons. Our results show that ASEs have higher KA values and higher KA/KS ratios than CSEs, indicating faster amino acid-level evolution in ASEs. The faster evolution may be in part due to weaker selective constraints. It is also possible that the faster rate is in part due to faster functional evolution in ASEs. On the other hand, the majority of ASEs have lower KS values than CSEs. With reference to the substitution rate in introns, we show that the KS values in ASEs are close to the neutral substitution rate, whereas the synonymous substitution rate in CSEs has likely been accelerated. The elevated synonymous rate in CSEs is not related to CpG dinucleotides or low-complexity regions of protein but may be weakly related to codon usage bias. The overall trends of higher KA and lower KS in ASEs than in CSEs are also observed in human-rat and mouse-rat comparisons. Therefore, our observations hold for mammals of different molecular clocks.

Protein-protein interaction networks have evolved mainly through connectivity rewiring and gene duplication. However, how protein function influences these processes and how a network grows in time have not been well studied. Using protein-protein interaction data and genomic data from the budding yeast, we first examined whether there is a correlation between the age and connectivity of yeast proteins. A steady increase in connectivity with protein age is observed for yeast proteins except for those that can be traced back to Eubacteria. Second, we investigated whether protein connectivity and duplicability vary with gene function. We found a higher average duplicability for proteins interacting with external environments than for proteins localized within intracellular compartments. For example, proteins that function in the cell periphery (mainly transporters) show a high duplicability but are lowly connected. Conversely, proteins that function within the nucleus (e.g., transcription, RNA and DNA metabolisms, and ribosome biogenesis and assembly) are highly connected but have a low duplicability. Finally, we found a negative correlation between protein connectivity and duplicability.

We compared the single-nucleotide polymorphisms (SNPs) in humans in 182 housekeeping and 148 tissue-specific genes. SNPs were divided into rare and common polymorphisms based on their frequencies. We found that housekeeping genes tend to be less polymorphic than tissue-specific genes for both rare and common SNPs. Using mouse as a second species for computing sequence divergences, we found no evidence of positive selection: for both housekeeping and tissue-specific genes, the ratio of nonsynonymous to synonymous common SNPs per site showed no significant difference from that of divergence. Similarly, we observed no evidence of positive selection for the 289 and 149 genes that have orthologs available for divergence calculation between humans and chimpanzees and between humans and Old World monkeys, respectively. A comparison with previous SNP studies suggests that approximately 20% of the nonsynonymous SNPs in the human population are nearly neutral and that positive selection in the human genome might not be as frequent as previously thought.

A general picture of the role of expression divergence in the evolution of duplicate genes is emerging, thanks to the availability of completely sequenced genomes and functional genomic data, such as microarray data. It is now clear that expression divergence, regulatory-motif divergence and coding-sequence divergence all increase with the age of duplicate genes, although their exact interrelationships remain to be determined. It is also clear that gene duplication increases expression diversity and enables tissue or developmental specialization to evolve. However, the relative roles of subfunctionalization and neofunctionalization in the retention of duplicate genes remain to be clarified, especially for higher eukaryotes. In addition, the relationship between gene duplication and evolution of transcriptional regulatory networks is largely unexplored.

Transcription factors (TFs), which are central to the regulation of gene expression, are usually members of multigene families. In plants, they are involved in diverse processes such as developmental control and elicitation of defense and stress responses. To investigate if differences exist in the expansion patterns of TF gene families between plants and other eukaryotes, we first used Arabidopsis (Arabidopsis thaliana) TFs to identify TF DNA-binding domains. These DNA-binding domains were then used to identify related sequences in 25 other eukaryotic genomes. Interestingly, among 19 families that are shared between animals and plants, more than 14 are larger in plants than in animals. After examining the lineage-specific expansion of TF families in two plants, eight animals, and two fungi, we found that TF families shared among these organisms have undergone much more dramatic expansion in plants than in other eukaryotes. Moreover, this elevated expansion rate of plant TF is not simply due to higher duplication rates of plant genomes but also to a higher degree of expansion compared to other plant genes. Further, in many Arabidopsis-rice (Oryza sativa) TF orthologous groups, the degree of lineage-specific expansion in Arabidopsis is correlated with that in rice. This pattern of parallel expansion is much more pronounced than the whole-genome trend in rice and Arabidopsis. The high rate of expansion among plant TF genes and their propensity for parallel expansion suggest frequent adaptive responses to selection pressure common among higher plants.

We examined the effects of recombination on the molecular evolution of noncoding regions in pseudoautosomal regions (PARs) and recombination hotspots in hominoids. The PAR-linked regions analyzed had on average longer branch lengths than those of the recombination hotspots. Moreover, contrary to previous observations, we found no correlation between recombination rate and silent site divergence in our data set and little change in the GC content during recent hominoid evolution. This suggests that the current rate of recombination is not a good indicator of the past rates of recombination for these highly recombining regions. Furthermore, human recombination hotspots show increased AT to GC substitutions in the human lineage, while no such pattern is detected for PAR-linked regions. Together, these observations suggest that recombination hotspots in hominoids are transient in the evolutionary time-scale. Interestingly, the 16p13.3 recombination hotspot locus violates a local molecular clock, though the locus appears to be noncoding and should evolve neutrally. We hypothesize that sudden changes in recombination rate have caused the changes in substitution rate at this locus.

In eukaryotes, C5-cytosine methylation is a common mechanism associated with a variety of functions such as gene regulation or control of genomic stability. Different subfamilies of eukaryotic methyltransferases (MTases) have been identified, mainly in metazoa, plants, and fungi. In this paper, we used hidden Markov models to detect MTases in completed or almost completed eukaryotic genomes, including different species of Protozoa. A phylogenetic analysis of MTases enabled us to define six subfamilies of MTases, including two new subfamilies. The dnmt1 subfamily that includes all the known MTases with a maintenance activity seems to be absent in the Protozoa. The dnmt2 subfamily seems to be the most widespread, being present even in the nonmethylated Dictyostelium discoideum. We also found two dnmt2 members in the bacterial genus Geobacter, suggesting that horizontal transfers of MTases occurred between eukaryotes and prokaryotes. Even if the direction of transfer cannot be determined, this relationship might be useful for understanding the function of this enigmatic subfamily of MTases. Globally, our analysis reveals a great diversity of MTases in eukaryotes, suggesting the existence of different methylation systems. Our results also suggest acquisitions and losses of different MTases in every eukaryotic lineage studied and that some eukaryotes appear to be devoid of methylation.

Understanding the evolution of gene function is a primary challenge of modern evolutionary biology. Despite an expanding database from genomic and developmental studies, we are lacking quantitative methods for analyzing the evolution of some important measures of gene function, such as gene-expression patterns. Here, we introduce phylogenetic comparative methods to compare different models of gene-expression evolution in a maximum-likelihood framework. We find that expression of duplicated genes has evolved according to a nonphylogenetic model, where closely related genes are no more likely than more distantly related genes to share common expression patterns. These results are consistent with previous studies that found rapid evolution of gene expression during the history of yeast. The comparative methods presented here are general enough to test a wide range of evolutionary hypotheses using genomic-scale data from any organism.

We analyzed the completed human genome for recent segmental duplications (size > or = 1 kb and sequence similarity > or = 90%). We found that approximately 4% of the genome is covered by duplications and that the extent of segmental duplication varies from 1% to 14% among the 24 chromosomes. Intrachromosomal duplication is more frequent than interchromosomal duplication in 15 chromosomes. The duplication frequencies in pericentromeric and subtelomeric regions are greater than the genome average by approximately threefold and fourfold. We examined factors that may affect the frequency of duplication in a region. Within individual chromosomes, the duplication frequency shows little correlation with local gene density, repeat density, recombination rate, and GC content, except chromosomes 7 and Y. For the entire genome, the duplication frequency is correlated with each of the above factors. Based on known genes and Ensembl genes, the proportion of duplications containing complete genes is 3.4% and 10.7%, respectively. The proportion of duplications containing genes is higher in intrachromosomal than in interchromosomal duplications, and duplications containing genes have a higher sequence similarity and tend to be longer than duplications containing no genes. Our simulation suggests that many duplications containing genes have been selectively maintained in the genome.

Three frequently used methods for estimating the synonymous and nonsynonymous substitution rates (Ks and Ka) were evaluated and compared for their accuracies; these methods are denoted by LWL85, LPB93, and GY94, respectively. For this purpose, we used a codon-evolution model to obtain the expected Ka and Ks values for the above three methods and compared the values with those obtained by the three methods. We also proposed some modifications of LWL85 and LPB93 to increase their accuracies. Our computer simulations under the codon-evolution model showed that for sequences < or =300 codons, the performance of GY94 may not be reliable. For longer sequences, GY94 is more accurate for estimating the Ka/Ks ratio than the modified LPB93 and LWL85 in the majority of the cases studied. This is particularly so when k > or = 3, which is the transition/transversion (mutation) rate ratio. However, when k is approximately 2 and when the sequence divergence is relatively large, the modified LWL85 performed better than GY94 and the modified LPB93. The inferiority of LPB93 to LWL85 is surprising because LPB93 was intended to improve LWL85. Also, it has been thought that the codon-based method of GY94 is better than the heuristic method of LWL85, but our simulation results showed that in many cases, the opposite was true, even though our simulation was based on the codon-evolution model.

The two suborders of bats, Megachiroptera (megabats) and Microchiroptera (microbats), use different sensory modalities for perceiving their environment. Megabats are crepuscular and rely on a well-developed eyes and visual pathway, whereas microbats occupy a nocturnal niche and use acoustic orientation or echolocation more than vision as the major means of perceiving their environment. In view of the differences associated with their sensory systems, we decided to investigate the function and evolution of color vision (opsin genes) in these two suborders of bats. The middle/long wavelength (M/L) and short wavelength (S) opsin genes were sequenced from two frugivorous species of megabats, Haplonycteris fischeri and Pteropus dasymallus formosus, and one insectivorous species of microbat, Myotis velifer. Contrary to the situation in primates, where many nocturnal species have lost the functional S opsin gene, both crepuscular and strictly nocturnal species of bats that we examined have functional M/L and S opsin genes. Surprisingly, the S opsin in these bats may be sensitive to UV light, which is relatively more abundant at dawn and at dusk. The M/L opsin in these bats appears to be the L type, which is sensitive to red and may be helpful for identifying fruits among leaves or for other purposes. Most interestingly, H. fischeri has a recent duplication of the M/L opsin gene, representing to date the only known case of opsin gene duplication in non-primate mammals. Some of these observations are unexpected and may provide insights into the effect of nocturnal life on the evolution of opsin genes in mammals and the evolution of the life history traits of bats in general.

Do housekeeping genes, which are turned on most of the time in almost every tissue, evolve more slowly than genes that are turned on only at specific developmental times or tissues? Recent large-scale gene expression studies enable us to have a better definition of housekeeping genes and to address the above question in detail. In this study, we examined 1581 human-mouse orthologous gene pairs for their patterns of sequence evolution, contrasting housekeeping genes with tissue-specific genes. Our results show that, in comparison to tissue-specific genes, housekeeping genes on average evolve more slowly and are under stronger selective constraints as reflected by significantly smaller values of Ka/Ks. Besides stronger purifying selection, we explored several other factors that can possibly slow down nonsynonymous rates in housekeeping genes. Although mutational bias might slightly slow the nonsynonymous rates in housekeeping genes, it is unlikely to be the major cause of the rate difference between the two types of genes. The codon usage pattern of housekeeping genes does not seem to differ from that of tissue-specific genes. Moreover, contrary to the old textbook concept, we found that approximately 74% of the housekeeping genes in our study belong to multigene families, not significantly different from that of the tissue-specific genes ( approximately 70%). Therefore, the stronger selective constraints on housekeeping genes are not due to a lower degree of genetic redundancy.

Although the evolutionary significance of gene duplication has long been recognized, it remains unclear what determines gene duplicability. We find protein complexity to be an important determinant because the proportion of unduplicated genes (P) increases with the number of subunits in a protein. However, P is high (>or=65%) for both monomers and multimers in yeast, but

The hominoid primates (apes and humans) exhibit remarkable diversity in their social and sexual behavioral systems. This is reflected in many ways in their anatomy and physiology. For example, the testes and seminal vesicles are relatively large in species with high sperm competition like the chimpanzee and small in species with low or no sperm competition like the gorilla. Additionally, the chimpanzee is the only hominoid primate known to produce a firm copulatory plug, which presumably functions in sperm competition by blocking insemination of subsequent males. Here we examine the molecular evolution of the semenogelin genes (SEMG1 and SEMG2), which code for the predominant structural proteins in human semen. High molecular weight complexes of these proteins are responsible for the viscous gelatinous consistency of human semen; their rodent homologs are responsible for the formation of a firm copulatory plug. Chimpanzees have an expanded SEMG1 gene caused by duplications of tandem repeats, each encoding 60 amino acids, resulting in a protein nearly twice as long as that of humans. In contrast, at both SEMG1 and SEMG2 we observed several gorilla haplotypes that contain at least one premature stop codon. We suggest that these structural changes in the semenogelin proteins that have arisen since the human-chimpanzee-gorilla split may be responsible for the physiological differences between these species ejaculated semen that correlate with their sociosexual behavior.

We developed a new evolutionary method for identifying exons from genomic sequences and found 19000 potential coding exons that are absent from all existing annotations of the human genome. Of these, 13700 satisfied very stringent criteria and can with confidence be considered as novel exons. Evidently, a large number of new human genes can be identified using evolutionary approaches.

Whether nonessential genes evolve faster than essential genes has been a controversial issue. To resolve this issue, we use the data from a nearly complete set of single-gene deletions in the yeast Saccharomyces cerevisiae to assess protein dispensability. Also, instead of the nematode, which was used previously but is only distantly related to S. cerevisiae, we use another yeast, Candida albicans, as a second species to estimate the evolutionary distances between orthologous genes in two species. Our analysis reveals only a weak correlation between protein dispensability and evolutionary rate. More important, the correlation disappears when duplicate genes are removed from the analysis. And surprisingly, the average rate of nonsynonymous substitution is considerably lower than that for single-copy genes in the yeast genome. This observation suggests that structural constraints are more important in determining the rate of evolution of a protein than dispensability because duplicate genes are on average more dispensable than single-copy genes. For duplicate genes, those with only a weak effect or no effect of deletion on fitness evolve on average faster than those with a moderate or strong effect of deletion on fitness, which in turn evolve on average faster than those with a lethal effect of deletion.

Changes in technology in the past decade have had such an impact on the way that molecular evolution research is done that it is difficult now to imagine working in a world without genomics or the Internet. In 1992, GenBank was less than a hundredth of its current size and was updated every three months on a huge spool of tape. Homology searches took 30 minutes and rarely found a hit. Now it is difficult to find sequences with only a few homologs to use as examples for teaching bioinformatics. For molecular evolution researchers, the genomics revolution has showered us with raw data and the information revolution has given us the wherewithal to analyze it. In broad terms, the most significant outcome from these changes has been our newfound ability to examine the evolution of genomes as a whole, enabling us to infer genome-wide evolutionary patterns and to identify subsets of genes whose evolution has been in some way atypical.