Wen-Hsiung Li

Unbiased estimation of the rates of synonymous and nonsynonymous substitution

SummaryThe current convention in estimating the number of substitutions per synonymous site (KS) and per nonsynonymous site (KA) between two protein-coding genes is to count each twofold degenerate site as one-third synonymous and two-thirds nonsynonymous because one of the three possible changes at such a site is synonymous and the other two are nonsynonymous. This counting rule can considerably overestimate theKS value because transitional mutations tend to occur more often than transversional mutations and because most transitional mutations at twofold degenerate sites are synonymous. A new method that gives unbiased estimates is proposed. An application of the new and the old method to 14 pairs of mouse and rat genes shows that the new method gives aKS value very close to the number of substitutions per fourfold degenerate site whereas the old method gives a value 30% higher. Both methods give aKA value close to the number of substitutions per nondegenerate site.

Comparative Analysis of the Receptor-Like Kinase Family in Arabidopsis and Rice

Receptor-like kinases (RLKs) belong to the large RLK/Pelle gene family, and it is known that the Arabidopsis thaliana genome contains >600 such members, which play important roles in plant growth, development, and defense responses. Surprisingly, we found that rice (Oryza sativa) has nearly twice as many RLK/Pelle members as Arabidopsis does, and it is not simply a consequence of a larger predicted gene number in rice. From the inferred phylogeny of all Arabidopsis and rice RLK/Pelle members, we estimated that the common ancestor of Arabidopsis and rice had >440 RLK/Pelles and that large-scale expansions of certain RLK/Pelle members and fusions of novel domains have occurred in both the Arabidopsis and rice lineages since their divergence. In addition, the extracellular domains have higher nonsynonymous substitution rates than the intracellular domains, consistent with the role of extracellular domains in sensing diverse signals. The lineage-specific expansions in Arabidopsis can be attributed to both tandem and large-scale duplications, whereas tandem duplication seems to be the major mechanism for recent expansions in rice. Interestingly, although the RLKs that are involved in development seem to have rarely been duplicated after the Arabidopsis–rice split, those that are involved in defense/disease resistance apparently have undergone many duplication events. These findings led us to hypothesize that most of the recent expansions of the RLK/Pelle family have involved defense/resistance-related genes.

Role of duplicate genes in genetic robustness against null mutations.

Deleting a gene in an organism often has little phenotypic effect, owing to two mechanisms of compensation. The first is the existence of duplicate genes: that is, the loss of function in one copy can be compensated by the other copy or copies. The second mechanism of compensation stems from alternative metabolic pathways, regulatory networks, and so on. The relative importance of the two mechanisms has not been investigated except for a limited study, which suggested that the role of duplicate genes in compensation is negligible. The availability of fitness data for a nearly complete set of single-gene-deletion mutants of the Saccharomyces cerevisiae genome has enabled us to carry out a genome-wide evaluation of the role of duplicate genes in genetic robustness against null mutations. Here we show that there is a significantly higher probability of functional compensation for a duplicate gene than for a singleton, a high correlation between the frequency of compensation and the sequence similarity of two duplicates, and a higher probability of a severe fitness effect when the duplicate copy that is more highly expressed is deleted. We estimate that in S. cerevisiae at least a quarter of those gene deletions that have no phenotype are compensated by duplicate genes.

Human polymorphism at microRNAs and microRNA target sites

MicroRNAs (miRNAs) function as endogenous translational repressors of protein-coding genes in animals by binding to target sites in the 3′ UTRs of mRNAs. Because a single nucleotide change in the sequence of a target site can affect miRNA regulation, naturally occurring SNPs in target sites are candidates for functional variation that may be of interest for biomedical applications and evolutionary studies. However, little is known to date about variation among humans at miRNAs and their target sites. In this study, we analyzed publicly available SNP data in context with miRNAs and their target sites throughout the human genome, and we found a relatively low level of variation in functional regions of miRNAs, but an appreciable level of variation at target sites. Approximately 400 SNPs were found at experimentally verified target sites or predicted target sites that are otherwise evolutionarily conserved across mammals. Moreover, ≈250 SNPs potentially create novel target sites for miRNAs in humans. If some variants have functional effects, they might confer phenotypic differences among humans. Although the majority of these SNPs appear to be evolving under neutrality, interestingly, some of these SNPs are found at relatively high population frequencies even in experimentally verified targets, and a few variants are associated with atypically long-range haplotypes that may have been subject to recent positive selection.

An evolutionary perspective on synonymous codon usage in unicellular organisms.

SummaryObserved patterns of synonymous codon usage are explained in terms of the joint effects of mutation, selection, and random drift. Examination of the codon usage in 165Escherichia coli genes reveals a consistent trend of increasing bias with increasing gene expression level. Selection on codon usage appears to be unidirectional, so that the pattern seen in lowly expressed genes is best explained in terms of an absence of strong selection. A measure of directional synonymous-codon usage bias, the Codon Adaptation Index, has been developed. In enterobacteria, rates of synonymous substitution are seen to vary greatly among genes, and genes with a high codon bias evolve more slowly. A theoretical study shows that the patterns of extreme codon bias observed for someE. coli (and yeast) genes can be generated by rather small selective differences. The relative plausibilities of various theoretical models for explaining nonrandom codon usage are discussed.

Uncovering Small RNA-Mediated Responses to Phosphate Deficiency in Arabidopsis by Deep Sequencing

Recent studies have demonstrated the important role of plant microRNAs (miRNAs) under nutrient deficiencies. In this study, deep sequencing of Arabidopsis (Arabidopsis thaliana) small RNAs was conducted to reveal miRNAs and other small RNAs that were differentially expressed in response to phosphate (Pi) deficiency. About 3.5 million sequence reads corresponding to 0.6 to 1.2 million unique sequence tags from each Pi-sufficient or Pi-deficient root or shoot sample were mapped to the Arabidopsis genome. We showed that upon Pi deprivation, the expression of miR156, miR399, miR778, miR827, and miR2111 was induced, whereas the expression of miR169, miR395, and miR398 was repressed. We found cross talk coordinated by these miRNAs under different nutrient deficiencies. In addition to miRNAs, we identified one Pi starvation-induced DICER-LIKE1-dependent small RNA derived from the long terminal repeat of a retrotransposon and a group of 19-nucleotide small RNAs corresponding to the 5′ end of tRNA and expressed at a high level in Pi-starved roots. Importantly, we observed an increased abundance of TAS4-derived trans-acting small interfering RNAs (ta-siRNAs) in Pi-deficient shoots and uncovered an autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−) that regulates the biosynthesis of anthocyanin. This finding sheds light on the regulatory network between miRNA/ta-siRNA and its target gene. Of note, a substantial amount of miR399* accumulated under Pi deficiency. Like miR399, miR399* can move across the graft junction, implying a potential biological role for miR399*. This study represents a comprehensive expression profiling of Pi-responsive small RNAs and advances our understanding of the regulation of Pi homeostasis mediated by small RNAs.

Evolutionary analyses of the human genome.

The completion of the human genome will greatly accelerate the development of a new branch of science—evolutionary genomics. We can now directly address important questions about the evolutionary history of human genes and their regulatory sequences. Computational analyses of the human genome will reveal the number of genes and repetitive elements, the extent of gene duplication and compositional heterogeneity in the human genome, and the extent of domain shuffling and domain sharing among proteins. Here we present some first glimpses of these features.

Patterns of Nucleotide Substitution in Pseudogenes and Functional Genes

SummaryThe pattern of point mutations is inferred from nucleotide substitutions in pseudogenes. The pattern obtained suggests that transition mutations occur somewhat more frequently than transversion mutations and that mutations result more often in A or T than in G or C. Our results are discussed with respect to the predictions from Topal and Frescos model for the molecular basis of point (substitution) mutations (Nature 263:285–289, 1976). The pattern of nucleotide substitution at the first and second positions of codons in functional genes is quite similar to that in pseudogenes, but the relative frequency of the transition C→T in the sense strand is drastically reduced and those of the transversions C→G and G→C are doubled. The differences between the two patterns can be explained by the observation that in the protein evolution amino acid substitutions occur mainly between amino acids with similar biochemical properties (Grantham, Science 185:862–864, 1974). Our results for the patterns of nucleotide substitutions in pseudogenes and in functional genes lead to the prediction that both the coding and non-coding regions of protein coding genes should have high frequencies of A and T. Available data show that the non-coding regions are indeed high in A and T but the coding regions are low in T, though high in A.

An evaluation of the molecular clock hypothesis using mammalian DNA sequences.

SummaryA statistical analysis of extensive DNA sequence data from primates, rodents, and artiodacytls clearly indicates that no global molecular clock exists in mammals. Rates of nucleotide subsitution in rodents are estimated to be four to eight times higher than those in higher primates and two to four times higher than those in artiodactyls. There is strong evidence for lower substitution rates in apes and humans than in monkeys, supporting the hominoid slowdown hypothesis. There is also evidence for lower rates in humans than in apes, suggesting a further rate slowdown in the human lineage after the separation of humans from apes. By contrast, substitution rates are nearly equal in mouse and rat. These results suggest that differences in generation time or, more precisely, in the number of germline DNA replications per year are the primary cause of rate differences in mammals. Further, these differences are more in line with the neutral mutation hypothesis than if the rates are the same for short-and long-living mammals.

Rates of Nucleotide Substitution in Angiosperm Mitochondrial DNA Sequences and Dates of Divergence Between Brassica and Other Angiosperm Lineages

Abstract. We obtained 16 nucleotide sequences (∼1400 bp each) of the first intron of the mitochondrial (mt) gene for NADH subunit 4 (nad4) from 10 species of Brassicaceae. Using these new sequences and five published sequences from GenBank, we constructed a phylogenetic tree of the Brassicaceae species under study and showed that the rate of nucleotide substitution in the first intron of nad4 is very low, about 0.16–0.23 × 10−9 substitution per site per year, which is about half of the silent rate in exons of nad4. The ratios of substitution rates in this intron, ITS, and IGS are approximately 1:23:73, where ITS is the nuclear intergenic spacer between 18S and 25S rRNA genes and IGS is the intergenic spacer of 5S rRNA genes. A segment (335 bp) in the first intron of nad4 in Brassicaceae species that is absent in wheat was considered as a nonfunctional sequence and used to estimate the neutral rate (the rate of mutation) in mtDNA to be 0.5–0.7 × 10−9 substitution per site per year, which is about three times higher than the substitution rate in the rest of the first intron of nad4. We estimated that the dates of divergence are 170–235 million years (Myr) for the monocot–dicot split, 112–156 Myr for the Brassicaceae–Lettuce split, 14.5–20.4 Myr for the Brassica–Arabidopsis split, and 14.5–20.4 Myr for the Arabidopsis–Arabideae split.