Amino acid sequences of stomach and nonstomach lysozymes of ruminants

Summary Complete amino acid sequences are presented for lysozymesc from camel and goat stomachs and compared to sequences of other lysozymesc. Tree analysis suggests that the rate of amino acid replacement went up as soon as lysozyme was recruited for the stomach function in early ruminants. The two lysozymes from goat stomach are the products of a gene duplication that probably took place before the divergence of cow, goat, and deer about 25 million years ago. Partial sequences of three lysozymes from goat tears indicated that (a) the goat tear family of lysozymes may have diverged from the stomach lysozyme family by an ancient duplication and (b) later duplications are probably responsible for the multiple forms of tear and milk lysozymes in ruminants.     

Molecular genetics and evolution of stomach and nonstomach lysozymes in the hoatzin

Multiple genes of the hoatzin encoding stomach lysozyme c and closely related members of this calcium-binding lysozyme c group were cloned from a genomic DNA library and sequenced. There are a minimum of five genes represented among these sequences that encode two distinct groups of protein sequences. One group of three genes corresponds to the stomach lysozyme amino acid sequences, and the remaining genes encode predicted proteins that are more basic in character and share several sequence identities with the pigeon egg-white lysozyme rather than with the hoatzin stomach lysozymes. Despite these structural similarities between some of the hoatzin gene products and the pigeon lysozyme, phylogenetic analyses indicate that all of the hoatzin sequences are closely related to one another. This is borne out by the relatively small genetic distances even in the intronic regions, which are not subject to the selective pressures operating on the coding regions of the stomach lysozymes. These results suggest that multiple gene duplication events have occurred during the evolution of hoatzin lysozymes.     

cDNA and amino acid sequences of rainbow trout (Oncorhynchus mykiss) lysozymes and their implications for the evolution of lysozyme and lactalbumin

Summary The complete 129-amino-acid sequences of two rainbow trout lysozymes (I and II) isolated from kidney were established using protein chemistry microtechniques. The two sequences differ only at position 86, I having aspartic acid and II having alanine. A cDNA clone coding for rainbow trout lysozyme was isolated from a cDNA library made from liver mRNA. Sequencing of the cloned cDNA insert, which was 1 kb in length, revealed a 432-bp open reading frame encoding an amino-terminal peptide of 15 amino acids and a mature enzyme of 129 amino acids identical in sequence to II. Forms I and II from kidney and liver were also analyzed using enzymatic amplification via PCR and direct sequencing; both organs contain mRNA encoding the two lysozymes. Evolutionary trees relating DNA sequences coding for lysozymesc and α-lactalbumins provide evidence that the gene duplication giving rise to conventional vertebrate lysozymesc and to lactalbumin preceded the divergence of fishes and tetrapods about 400 Myr ago. Evolutionary analysis also suggests that amino acid replacements may have accumulated more slowly on the lineage leading to fish lysozyme than on those leading to mammal and bird lysozymes.     

Detecting lineage-specific adaptive evolution of brain-expressed genes in human using rhesus macaque as outgroup

Comparative genetic analysis between human and chimpanzee may detect genetic divergences responsible for human-specific characteristics. Previous studies have identified a series of genes that potentially underwent Darwinian positive selection during human evolution. However, without a closely related species as outgroup, it is difficult to identify human-lineage-specific changes, which is critical in delineating the biological uniqueness of humans. In this study, we conducted phylogeny-based analyses of 2633 human brain-expressed genes using rhesus macaque as the outgroup. We identified 47 candidate genes showing strong evidence of positive selection in the human lineage. Genes with maximal expression in the brain showed a higher evolutionary rate in human than in chimpanzee. We observed that many immune-defense-related genes were under strong positive selection, and this trend was more prominent in chimpanzee than in human. We also demonstrated that rhesus macaque performed much better than mouse as an outgroup in identifying lineage-specific selection in humans.     

Episodic positive selection during the evolution of naphthalene dioxygenase to nitroarene dioxygenase

Using different maximum-likelihood models of adaptive evolution, signatures of natural selective pressure, operating across the naphthalene family of dioxygenases, were examined. A lineage- and branch-site specific combined analysis revealed that purifying selection pressure dominated the evolutionary history of the enzyme family. Specifically, episodic positive Darwinian selection pressure, affecting only a few sites in a subset of lineages, was found to be responsible for the evolution of nitroarene dioxygenases (NArDO) from naphthalene dioxygenase (NDO). Site-specific analysis confirmed the absence of diversifying selection pressure at any particular site. Different sets of positively selected residues, obtained from branch-site specific analysis, were detected for the evolution of each NArDO. They were mainly located around the active site, the catalytic pocket and their adjacent regions, when mapped onto the 3D structure of the α-subunit of NDO. The present analysis enriches the current understanding of adaptive evolution and also broadens the scope for rational alteration of substrate specificity of enzyme by directed evolution.     

啮齿目核糖酸酶5基因(RNase5)的分子进化

RNase5是RNASE A基因超家族中的一个重要成员,是分子进化研究的理想模型之一。基于基因组水平,我们对啮齿目的3个进化枝10科17个物种开展RNase5的分子进化研究。利用TBlastN及BlastN方法鉴定每个基因组的RNase5基因,发现该基因在啮齿目的Ctenohystrica所有物种发生丢失,时间是在Ctenohystrica形成之后;邻接法和最大似然法构建的系统发育树均支持RNase5在“与小家鼠相关的进化枝”的小家鼠、褐家鼠和拉布拉多白足鼠发生三次独立基因复制事件;利用PAML软件的枝模型、位点模型及枝-位点模型计算选择压力,均检测到RNase5基因受到强烈的正选择作用。总之,我们的研究深入系统开展了RNase5在啮齿目中的分子进化,增加了该基因研究的多样性,为进一步系统认识该基因在动物的适应性进化遗传机制奠定了基础。     

Comparison of goose-type,chicken-type,and phage-type lysozymes illustrates the changes that occur in both amino acid sequence and three-dimensional structure during evolution

Summary The three-dimensional structure of goose-type lysozyme (GEWL), determined by x-ray crystallography and refined at high resolution, has similarities to the structures of hen (chicken) eggwhite lysozyme (HEWL) and bacteriophage T4 lysozyme (T4L). The nature of the structural correspondence suggests that all three classes of lysozyme diverged from a common evolutionary precursor, even though their amino acid sequences appear to be unrelated (Grütter et al. 1983).In this paper we make detailed comparisons of goose-type, chicken-type, and phage-type lysozymes. The lysozymes have undergone conformational changes at both the blobal and the local level. As in the globins, there are corresponding -helices that have rigid-body displacements relative to each other, but in some cases corresponding helices have increased or decreased in length, and in other cases there are helices in one structure that have no counterpart in another.Independent of the overall structural correspondence among the three lysozyme backbones is another, distinct correspondence between a set of three consecutive -helices in GEWL and three consecutive -helices in T4L. This structural correspondence could be due, in part, to a common energetically favorable contact between the first and the third helices.There are similarities in the active sites of the three lysozymes, but also one striking difference. Glu 73 (GEWL) spatially corresponds to Glu 35 (HEWL) and to Glu 11 (T4L). On the other hand, there are two aspartates in the GEWL active site, Asp 86 and Asp 97, neither of which corresponds exactly to Asp 52 (HEWL) or Asp 20 (T4L). (The discrepancy in the location of the carboxyl groups is about 10 Å for Asp 86 and 4 Å for Asp 97.) This lack of structural correspondence may reflect some differences in the mechanisms of action of three lysozymes. When the amino acid sequences of the three lysozyme types are aligned according to their structural correspondence, there is still no apparent relationship between the sequences except for possible weak matching in the vicinity of the active sites.     

Episodic sitewise positive selection on the signal recognition particle protein Ffh in Actinobacteria

In this work, we report a case of episodic sitewise positive selection acting on the highly conserved SRP protein Ffh in Actinobacteria. An elevated non-synonymous to synonymous mutation ratio (ω) was observed along the non-terminal branches of the species tree, which contained 11 Actinobacteria species, where positively selected residues were frequently observed within the signal sequence-binding domain. Together with the estimated lineage-specific ω ratios for several core components in the major protein translocation systems, our data suggest that this positive selection might be partially driven by the diversification of signal sequences.     

Molecular diversity and evolution of myticin-C antimicrobial peptide variants in the Mediterranean mussel, Mytilus galloprovincialis

Mussels have diverse groups of cysteine rich, cationic antimicrobial peptides (AMPs) (defensins, mytilins, myticins, and mytimycin) that constitute an important component of their innate immune defence. Despite the identification and characterization of these AMPs in mussels, the underlying genetic mechanisms that maintain high diversity among multiple variants of the myticin-C isoform are poorly understood. Using phylogeny-based models of sequence evolution and several site-by-site frequency spectrum statistical tests for neutrality, herein we report that positive selection has been the major driving force in maintaining high diversity among the allelic-variants of the myticin-C AMP of Mytilus galloprovincialis. The statistical tests rejected the hypothesis that all polymorphism within myticin-C loci is neutral. Although a majority of the codons constrained to purifying selection (rate of amino acid replacement to the silent substitution, omega < 1), approximately 8% of the codons with omega approximately equal to 5.5 are under positive selection (omega > 1), thus indicating adaptive evolution of certain amino acids. Direct interaction of these peptides with the surrounding pathogens and/or altered/new pathogens in the changing environment is the likely cause of molecular adaptation of certain amino acid sites in myticin-C variants.     

Molecular evolution in yeast of biotechnological interest

The importance of yeast in the food and beverage industries was only realized about 1860, when the role of these organisms in food manufacture became evident. Since they grow on a wide range of substrates and can tolerate extreme physicochemical conditions, yeasts, especially the genera Saccharomyces and Kluyveromyces, have been applied to many industrial processes, Industrial strains of these genera are highly specialized organisms that have evolved to utilize a range of environments and ecological niches to their full potential. This adaptation is called "domestication". This review describes the phylogenetic relationships among Saccharomyces and Kluyveromyces species and the different mechanisms involved in the adaptive evolution of industrial yeast strains.     

Estimating the fraction of invariable codons with a capture-recapture method

Summary A codon-based approach to estimating the number of variable sites in a protein is presented. When first and second positions of codons are assumed to be replacement positions, a capture-recapture model can be used to estimate the number of variable codons from every pair of homologous and aligned sequences. The capture-recapture estimate is compared to a maximum likelihood estimate of the number of variable codons and to previous approaches that estimate the number of variable sites (not codons) in a sequence. Computer simulations are presented that show under which circumstances the capture-recapture estimate can be used to correct biases in distance matrices. Analysis of published sequences of two genes, calmodulin and serum albumin, shows that distance corrections that employ a capture-recapture estimate of the number of variable sites may be considerably different from corrections that assume that the number of variable sites is equal to the total number of positions in the sequence. Offprint requests to: A. Sidow     

Uncovering adaptive evolution in the human lineage

Background

The recent increase in human polymorphism data, together with the availability of genome sequences from several primate species, provides an unprecedented opportunity to investigate how natural selection has shaped human evolution.

Results

We compared human branch-specific substitutions with variation data in the current human population to measure the impact of adaptive evolution on human protein coding genes. The use of single nucleotide polymorphisms (SNPs) with high derived allele frequencies (DAFs) minimized the influence of segregating slightly deleterious mutations and improved the estimation of the number of adaptive sites. Using DAF ≥ 60% we showed that the proportion of adaptive substitutions is 0.2% in the complete gene set. However, the percentage rose to 40% when we focused on genes that are specifically accelerated in the human branch with respect to the chimpanzee branch, or on genes that show signatures of adaptive selection at the codon level by the maximum likelihood based branch-site test. In general, neural genes are enriched in positive selection signatures. Genes with multiple lines of evidence of positive selection include taxilin beta, which is involved in motor nerve regeneration and syntabulin, and is required for the formation of new presynaptic boutons.

Conclusions

We combined several methods to detect adaptive evolution in human coding sequences at a genome-wide level. The use of variation data, in addition to sequence divergence information, uncovered previously undetected positive selection signatures in neural genes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-599) contains supplementary material, which is available to authorized users.     

Genes that determine flower color: the role of regulatory changes in the evolution of phenotypic adaptations

A central goal of evolutionary genetics is to trace the causal pathway between mutations at particular genes and adaptation at the phenotypic level. The proximate objective is to identify adaptations through the analysis of molecular sequence data from specific candidate genes or their regulatory elements. In this paper, we consider the molecular evolution of floral color in the morning glory genus (Ipomoea) as a model for relating molecular and phenotypic evolution. To begin, flower color variation usually conforms to simple Mendelian transmission, thus facilitating genetic and molecular analyses. Population genetic studies of flower color polymorphisms in the common morning glory (Ipomoea purpurea) have shown that some morphs are subject to complex patterns of selection. Striking differences in floral color and morphology are also associated with speciation in the genus Ipomoea. The molecular bases for these adaptive shifts can be dissected because the biosynthetic pathways that determine floral pigmentation are well understood and many of the genes of flavonoid biosynthesis have been isolated and extensively studied. We present a comparative analysis of the level of gene expression in Ipomoea for several key genes in flavonoid biosynthesis. Specifically we ask: how frequently are adaptive shifts in flower color phenotypes associated with changes in regulation of gene expression versus mutations in structural genes? The results of this study show that most species differences in this crucial phenotype are associated with changes in the regulation of gene expression.     

Molecular population genetics and the search for adaptive evolution in plants

The first papers on plant molecular population genetics were published approximately 10 years ago. Since that time, well over 50 additional studies of plant nucleotide polymorphism have been published, and many of these studies focused on detecting the signature of balancing or positive selection at a locus. In this review, we discuss some of the theoretical and statistical issues surrounding the detection of selection, with focus on plant populations, and we also summarize the empirical plant molecular population genetics literature. At face value, the literature suggests that a history of balancing or positive selection in plant genes is rampant. In two well-studied taxa (maize and Arabidopsis) over 20% of studied genes have been interpreted as containing the signature of selection. We argue that this is probably an overstatement of the prevalence of natural selection in plant genomes, for two reasons. First, demographic effects are difficult to incorporate and have generally not been well integrated into the plant population genetics literature. Second, the genes studied to date are not a random sample, so selected genes may be overrepresented. The next generation of studies in plant molecular population genetics requires additional sampling of local populations, explicit comparisons among loci, and improved theoretical methods to control for demography. Eventually, candidate loci should be confirmed by explicit consideration of phenotypic effects.     

Evidence of selection pressures of neuraminidase gene (NA) of influenza A virus subtype H5N1 on different hosts in Guangxi Province of China

In the present study, the possible evidence of positive selection was analyzed for the neuraminidase (NA) sequences of Guangxi H5N1 strains of China. Based on an overall site-specific positive selection analysis, it was found that NA gene of H5N1 Guangxi strains underwent purifying selection and no significant positively selected sites were identified. For the branch-specific positive selection analysis, there was no positive selection evidence for the branches leading to different poultry hosts (chicken, duck and goose). Conclusively, positive selection seems not possible (if not rare) for the NA gene in influenza H5N1 subtype, at least for the samples found in Guangxi Province of China.     

基于全基因组正选择基因揭示大狐蝠和小棕蝠功能分化分子机制

小棕蝠(Myotis lucifugus)和大狐蝠(Pteropus vampyrus)在生理和行为上有诸多差异,应被诉诸于分子水平并加以系统探讨。该研究对包括小棕蝠和大狐蝠在内的7个哺乳动物全基因组同源编码序列进行了高质量比对,使用比对序列进行了分别以两个蝙蝠枝为前景枝的全基因组水平选择压力分析,并对两种蝙蝠中的正选择基因进行了富集聚类分析。结果表明,在小棕蝠枝受到正选择的基因数高于大狐蝠枝;两种蝙蝠中独立受到正选择的基因富集于不同的功能类别;且正选择基因的富集差异与小棕蝠和大狐蝠在免疫、运动协调、能量代谢和感觉器官发育等关键生物学功能方面的分化大体吻合。     

Molecular and functional evolution of human DHRS2 and DHRS4 duplicated genes

Human DHRS2 and DHRS4 genes code for similar NADP-dependent short-chain carbonyl-reductase enzymes having different substrate specificity. Human DHRS2 and DHRS4 enzymes share several common sequence motives including residues responsible for coenzyme binding as well as for the intimate catalytic oxido-reductase mechanism, while their substrate-binding sequences have very low similarity. We found that DHRS2 and DHRS4 genes are syntenic outparalogues originated from a duplication of the DHRS4 gene that took place before the formation of the mammalian clade. DHRS2 gene evolved more rapidly and underwent positive selection on more sites than the DHRS4 gene. DHRS2 sites under positive selection were mainly located on the enzyme active site thus showing that substrate specificity drove the divergence from the DHRS4 enzyme. Rapid divergent evolution brought the human DHRS2 enzyme to have subcellular localization, synthesis regulation and specialized cellular functions very different from those of the human DHRS4 enzyme.     

Molecular evolution of the antiretroviral <Emphasis Type="Italic">TRIM5</Emphasis> gene

In 2004, the first report of TRIM5α as a cellular antiretroviral factor triggered intense interest among virologists, particularly because some primate orthologs of TRIM5α have activity against HIV. Since that time, a complex and eventful evolutionary history of the TRIM5 locus has emerged. A review of the TRIM5 literature constitutes a veritable compendium of evolutionary phenomena, including elevated rates of nonsynonymous substitution, divergence in subdomains due to short insertions and deletions, expansions and contractions in gene copy number, pseudogenization, balanced polymorphism, trans-species polymorphism, convergent evolution, and the acquisition of new domains by exon capture. Unlike most genes, whose history is dominated by long periods of purifying selection interspersed with rare instances of genetic innovation, analysis of restriction factor loci is likely to be complicated by the unpredictable and more-or-less constant influence of positive selection. In this regard, the molecular evolution and population genetics of restriction factor loci most closely resemble patterns that have been documented for immunity genes, such as class I and II MHC genes, whose products interact directly with microbial targets. While the antiretroviral activity encoded by TRIM5 provides plausible mechanistic hypotheses for these unusual evolutionary observations, evolutionary analyses have reciprocated by providing significant insights into the structure and function of the TRIM5α protein. Many of the lessons learned from TRIM5 should be applicable to the study of other restriction factor loci, and molecular evolutionary analysis could facilitate the discovery of new antiviral factors, particularly among the many TRIM genes whose functions remain as yet unidentified.     

Understanding the evolution and stability of the G-matrix

The G-matrix summarizes the inheritance of multiple, phenotypic traits. The stability and evolution of this matrix are important issues because they affect our ability to predict how the phenotypic traits evolve by selection and drift. Despite the centrality of these issues, comparative, experimental, and analytical approaches to understanding the stability and evolution of the G-matrix have met with limited success. Nevertheless, empirical studies often find that certain structural features of the matrix are remarkably constant, suggesting that persistent selection regimes or other factors promote stability. On the theoretical side, no one has been able to derive equations that would relate stability of the G-matrix to selection regimes, population size, migration, or to the details of genetic architecture. Recent simulation studies of evolving G-matrices offer solutions to some of these problems, as well as a deeper, synthetic understanding of both the G-matrix and adaptive radiations.     

Heterogeneous Nature and Distribution of Interruptions in Dinucleotides May Indicate the Existence of Biased Substitutions Underlying Microsatellite Evolution

Some aspects of microsatellite evolution, such as the role of base substitutions, are far from being fully understood. To examine the significance of base substitutions underlying the evolution of microsatellites we explored the nature and the distribution of interruptions in dinucleotide repeats from the human genome. The frequencies that we inferred in the repetitive sequences were statistically different from the frequencies observed in other noncoding sequences. Additionally, we detected that the interruptions tended to be towards the ends of the microsatellites and 5'-3' asymmetry. In all the estimates nucleotides forming the same repetitive motif seem to be affected by different base substitution rates in AC and AG. This tendency itself could generate patterning and similarity in flanking sequences and reconcile these phenomena with the high mutation rate found in flanking sequences without invoking convergent evolution. Nevertheless, our data suggest that there is a regional bias in the substitution pattern of microsatellites. The accumulation of random substitutions alone cannot explain the heterogeneity and the asymmetry of interruptions found in this study or the relative frequency of different compound microsatellites in the human genome. Therefore, we cannot rule out the possibility of a mutational bias leading to convergent or parallel evolution in flanking sequences.