Likelihood Analysis of the Chalcone Synthase Genes Suggests the Role of Positive Selection in Morning Glories (<Emphasis Type="Italic">Ipomoea</Emphasis>)

Chalcone synthase (CHS) is a key enzyme in the biosynthesis of flavonoides, which are important for the pigmentation of flowers and act as attractants to pollinators. Genes encoding CHS constitute a multigene family in which the copy number varies among plant species and functional divergence appears to have occurred repeatedly. In morning glories (Ipomoea), five functional CHS genes (A–E) have been described. Phylogenetic analysis of the Ipomoea CHS gene family revealed that CHS A, B, and C experienced accelerated rates of amino acid substitution relative to CHS D and E. To examine whether the CHS genes of the morning glories underwent adaptive evolution, maximum-likelihood models of codon substitution were used to analyze the functional sequences in the Ipomoea CHS gene family. These models used the nonsynonymous/synonymous rate ratio ( = d_N/d_S) as an indicator of selective pressure and allowed the ratio to vary among lineages or sites. Likelihood ratio test suggested significant variation in selection pressure among amino acid sites, with a small proportion of them detected to be under positive selection along the branches ancestral to CHS A, B, and C. Positive Darwinian selection appears to have promoted the divergence of subfamily ABC and subfamily DE and is at least partially responsible for a rate increase following gene duplication.     

Selection for Antimicrobial Peptide Diversity in Frogs Leads to Gene Duplication and Low Allelic Variation

Antimicrobial peptides are highly diverse pathogen-killing molecules. In many taxa, their evolution is characterized by positive selection and frequent gene duplication. It has been proposed that genes encoding antimicrobial peptides might be subject to balancing selection and/or an enhanced mutation rate, but these hypotheses have not been well evaluated because allelic variation has rarely been studied at antimicrobial peptide loci. We present an evolutionary analysis of novel antimicrobial peptide genes from leopard frogs, Rana. Our results demonstrate that a single genome contains multiple homologous copies, among which there is an excess of nonsynonymous nucleotide site divergence relative to that expected from synonymous site divergence. Thus, we confirm the trends of recurrent duplication and positive selection. Allelic variation is quite low relative to interspecies divergence, indicating a recent positive selective sweep with no evidence of balancing selection. Repeated gene duplication, rather than a balanced maintenance of divergent allelic variants at individual loci, appears to be how frogs have responded to selection for a diverse suite of antimicrobial peptides. Our data also support a pattern of enhanced synonymous site substitution in the mature peptide region of the gene, but we cannot conclude that this is due to an elevated mutation rate.     

Evolution of the proteasome components

 A phylogenetic analysis of proteasome subunits revealed two major families (α and β) which originated by an ancient gene duplication prior to the divergence of archaebacteria and eukaryotes. Numerous gene duplications have subsequently occurred in eukaryotes; at least nine of these duplications were shown to have occurred prior to the divergence of animals and fungi. In mammals, two genes encoding proteasome subunits (LMP2 and LMP7) are located in the major histocompatibility complex (MHC) region and play a specific role in generation of peptides for presentation by class I MHC molecules. Phylogenetic analysis of LMP7 and related sequences from mammals and lower vertebrates indicated that this locus arose by gene duplication prior to the divergence of jawed and jawless vertebrates; the time of this duplication was estimated to have been about 600 million years ago. The evolutionary history of the proteasome subunits provides support for a model of the evolution of new gene function postulating that, after gene duplication, the proteins encoded by daughter loci can adapt to specialized functions previously performed by the product of a single generalized ancestral locus. Received: 19 August 1996 / Revised: 24 December 1996     

Structural analysis and chromosomal localization of the mouse Psmb5 gene coding for the constitutively expressed β-type proteasome subunit

 The proteasome is a multi-subunit protease responsible for the production of peptides presented by major histocompatibility complex class I molecules. Accumulated evidence indicates that, upon stimulation with interferon-γ (IFN-γ), three β-type subunits, designated LMP2, LMP7, and PSMB10, are incorporated into the 20S proteasome by displacing the housekeeping β-type subunits designated PSMB6, PSMB5, and PSMB7, respectively. These changes in the subunit composition appear to facilitate class I-mediated antigen presentation, presumably by altering the cleavage specificities of the proteasome. In the present study, we determined the organization of the mouse gene Psmb5, coding for the PSMB5 subunit. Psmb5 is made up of three exons, spanning ∼5 kilobases. Its exon-intron organization differs radically from those of the other IFN-γ-regulated, β-type subunit genes including Lmp7 with which Psmb5 is believed to share an immediate common ancestor. The structure of the mouse Psmb5 gene is identical to that of its recently characterized human counterpart. Thus, the unique organization of the gene coding for the PSMB5 subunit appears to have been established before mammalian radiation. As well as the Psmb5 gene, the mouse genome contains a processed pseudogene designated Psmb5-ps. Interspecific backcross mapping showed that Psmb5 maps close to the Gtrgal2 locus on chromosome 14 and that Psmb5-ps is located in the vicinity of the Psme3 locus on chromosome 11. These results were confirmed by fluorescent in situ hybridization analysis that localized Psmb5 to band C2 to proximal D1 of chromosome 14 and Psmb5-ps to band D of chromosome 11. Received: 29 May 1997 / Revised: 4 June 1997     

Molecular Evolution at the Cytochrome Oxidase Subunit 2 Gene Among Divergent Populations of the Intertidal Copepod, Tigriopus californicus

The cytochrome c oxidase subunit 2 gene (COII) encodes a highly conserved protein that is directly responsible for the initial transfer of electrons from cytochrome c to cytochrome c oxidase (COX) crucial to the production of ATP during cellular respiration. Despite its integral role in electron transport, we have observed extensive intraspecific nucleotide and amino acid variation among 26 full-length COII sequences sampled from seven populations of the marine copepod, Tigriopus californicus. Although intrapopulation divergence was virtually nonexistent, interpopulation divergence at the COII locus was nearly 20% at the nucleotide level, including 38 nonsynonymous substitutions. Given the high degree of interaction between the cytochrome c oxidase subunit 2 protein (COX2) and the nuclear-encoded subunits of COX and cytochrome c (CYC), we hypothesized that some codons in the COII gene are likely to be under positive selection in order to compensate for amino acid substitutions in other subunits. Estimates of the ratio of nonsynonymous to synonymous substitution (ω), obtained using a series of maximum likelihood models of codon substitution, indicated that the majority of codons in T. californicus COII are under strong purifying selection (ω << 1), while approximately 4% of the sites in this gene appear to evolve under relaxed selective constraint (ω = 1). A branch-site maximum likelihood model identified three sites that may have experienced positive selection within the central California sequence clade in our COII phylogeny; these results are consistent with previous studies showing functional and fitness consequences among interpopulation hybrids between central and northern California populations. [Reviewing Editor: Dr. Willie Swanson]     

Protein Evolutionary Rates Correlate with Expression Independently of Synonymous Substitutions in Helicobacter pylori

In free-living microorganisms, such as Escherichia coli and Saccharomyces cerevisiae, both synonymous and nonsynonymous substitution frequencies correlate with expression levels. Here, we have tested the hypothesis that the correlation between amino acid substitution rates and expression is a by-product of selection for codon bias and translational efficiency in highly expressed genes. To this end, we have examined the correlation between protein evolutionary rates and expression in the human gastric pathogen Helicobacter pylori, where the absence of selection on synonymous sites enables the two types of substitutions to be uncoupled. The results revealed a statistically significant negative correlation between expression levels and nonsynonymous substitutions in both H. pylori and E. coli. We also found that neighboring genes located on the same, but not on opposite strands, evolve at significantly more similar rates than random gene pairs, as expected by co-expression of genes located in the same operon. However, the two species differ in that synonymous substitutions show a strand-specific pattern in E. coli, whereas the weak similarity in synonymous substitutions for neighbors in H. pylori is independent of gene orientation. These results suggest a direct influence of expression levels on nonsynonymous substitution frequencies independent of codon bias and selective constraints on synonymous sites. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Nicolas Galtier]     

siRNA下调20S蛋白酶体β5亚基抑制立氏立克次体细胞内增殖

[目的]探究宿主20S蛋白酶体β5亚基(proteasome 20S subunit beta 5,PSMB5)对革兰氏阴性专性胞内寄生菌立氏立克次体胞内生长繁殖的影响.[方法]利用小干扰RNA转染Vero和THP-1宿主细胞,设置未转染细胞为空白对照组、无义小干扰RNA转染组为阴性对照组、PSMB5特异性小干扰RNA...     

Positive and Negative Selection in the β-Esterase Gene Cluster of the Drosophila melanogaster Subgroup

We examine the pattern of molecular evolution of the β-esterase gene cluster, including the Est-6 and ψEst-6 genes, in eight species of the Drosophila melanogaster subgroup. Using maximum likelihood estimates of nonsynonymous/synonymous rate ratios, we show that the majority of Est-6 sites evolves under strong (48% of sites) or moderate (50% of sites) negative selection and a minority of sites (1.5%) is under significant positive selection. Est-6 sites likely to be under positive selection are associated with increased intraspecific variability. One positively selected site is responsible for the EST-6 F/S allozyme polymorphism; the same site is responsible for the EST-6 functional divergence between species of the melanogaster subgroup. For ψEst-6 83.7% sites evolve under negative selection, 16% sites evolve neutrally, and 0.3% sites are under positive selection. The positively selected sites of ψEst-6 are located at the beginning and at the end of the gene, where there is reduced divergence between D. melanogaster and D. simulans; these regions of ψEst-6 could be involved in regulation or some other function. Branch-site-specific analysis shows that the evolution of the melanogaster subgroup underwent episodic positive selection. Collating the present data with previous results for the β-esterase genes, we propose that positive and negative selection are involved in a complex relationship that may be typical of the divergence of duplicate genes as one or both duplicates evolve a new function. [Reviewing Editor: Dr. Martin Kreitman]     

The<Emphasis Type="Italic"> Trox-2</Emphasis> Hox/ParaHox gene of<Emphasis Type="Italic"> Trichoplax</Emphasis> (Placozoa) marks an epithelial boundary

Hox and ParaHox genes are implicated in axial patterning of cnidarians and bilaterians, and are thought to have originated by tandem duplication of a single ProtoHox gene followed by duplication of the resultant gene cluster. It is unclear what the ancestral role of Hox/ParaHox genes was before the divergence of Cnidaria and Bilateria, or what roles the postulated ProtoHox gene(s) played. Here we describe the full coding region, spatial expression and function of Trox-2, the single Hox/ParaHox-type gene identified in Trichoplax adhaerens (phylum Placozoa) and either a candidate ProtoHox or a ParaHox gene. Trox-2 is expressed in a ring around the periphery of Trichoplax, in small cells located between the outer margins of the upper and lower epithelial cell layers. Inhibition of Trox-2 function, either by uptake of morpholino antisense oligonucleotides or by RNA interference, causes complete cessation of growth and binary fission. We speculate that Trox-2 functions within a hitherto unrecognized population of possibly multipotential peripheral stem cells that contribute to differentiated cells at the epithelial boundary of Trichoplax.Edited by D. Tautz     

Molecular Evolution of a Small Gene Family of Wound Inducible Kunitz Trypsin Inhibitors in Populus

Maximum likelihood models of codon substitutions were used to analyze the molecular evolution of a Kunitz trypsin inhibitor (KTI) gene family in Populus and Salix. The methods support previous assertions that the KTI genes comprise a rapidly evolving gene family. Models that allow for codon specific estimates of the ratio of nonsynonymous to synonymous substitutions (ω) among sites detect positive Darwinian selection at several sites in the KTI protein. In addition, branch-specific maximum likelihood models show that there is significant heterogeneity in ω among branches of the KTI phylogeny. In particular, ω is substantially higher following duplication than speciation. There is also evidence for significant rate heterogeneity following gene duplication, suggesting different evolutionary rates in newly arisen gene duplicates. The results indicate uneven evolutionary rates both between sites in the KTI protein and among different lineages in the KTI phylogeny, which is incompatible with a neutral model of sequence evolution. [Reviewing Editor: Dr. Willie J. Swanson]     

Expression Patterns of Duplicate Genes in the Developing Root in Arabidopsis thaliana

Data on gene expression in the development of the root in Arabidopsis thaliana were used to test for expression profile differences among multi-gene families and to examine the extent to which expression differences accompanied coding sequences divergence within families. Significant differences among families were observed on two principal axes, accounting for over 80% of the variance in the expression data. The number of synonymous nucleotide substitutions per synonymous site (d_S) and the number of nonsynonymous nucleotide substitutions per nonsynonymous site (d_N) were estimated between the members of two-member families (N=428) and between phylogenetically independent sister pairs (N=190) of sequences within larger families. Ribosomal proteins and a few other proteins were exceptional in showing highly divergent expression patterns in spite of very low levels of amino acid sequence divergence, as indicated by the low d_N relative to d_S. However, the majority of gene duplicates showed relatively high levels of amino acid sequence divergence without appreciable change in expression pattern in the cell types analyzed. Reviewing Editor:Dr. Manyuan Long     

Comparative genomic analysis of the proteasome β5t subunit gene: implications for the origin and evolution of thymoproteasomes

The thymoproteasome is a recently discovered, specialized form of 20S proteasomes expressed exclusively in the thymic cortex. Although the precise molecular mechanism by which the thymoproteasome exerts its function remains to be elucidated, accumulating evidence indicates that it plays a crucial role in positive selection of T cells. In the present study, we analyzed the evolution of the β5t subunit, a β-type catalytic subunit uniquely present in thymoproteasomes. The gene coding for the β5t subunit, designated PSMB11, was identified in the cartilaginous fish, the most divergent group of jawed vertebrates compared to the other jawed vertebrates, but not in jawless vertebrates or invertebrates. Interestingly, teleost fish have two copies of apparently functional PSMB11 genes, designated PSMB11a and PSMB11b, that encode β5t subunits with distinct amino acids in the S1 pocket. BLAST searches of genome databases suggest that birds such as chickens, turkey, and zebra finch lost the PSMB11 gene, and have neither thymoproteasomes nor immunoproteasomes. In mammals, reptiles, amphibians, and teleost fishes, the PSMB11 gene (the PSMB11a gene in teleost fish) is located next to the PSMB5 gene coding for the β5 subunit of the standard 20S proteasome, indicating that the PSMB11 gene arose by tandem duplication from the evolutionarily more ancient PSMB5 gene. The general absence of introns in PSMB11 and an unusual exon–intron structure of jawed vertebrate PSMB5 suggest that PSMB5 lost introns and duplicated in tandem in a common ancestor of jawed vertebrates, with PSMB5 subsequently gaining two introns and PSMB11 remaining intronless.     

Eyespot resistance gene Pch-1 from Aegilops ventricosa is associated with a different chromosome in wheat line H-93-70 than the resistance factor in “Roazon” wheat

Summary The hexaploid wheat line H-93-70 carries a gene (Pch-1) that has been transferred from the wild grass Aegilops ventricosa and confers a high degree of resistance to eyespot diesease, caused by the fungus Pseudocercosporella herpotrichoides. Crosses of the resistant line H-93-70 with the susceptible wheat Pané 247 and with a 7D/7Ag wheat/Agropyron substitution line were carried out and F2 kernels were obtained. The kernels were cut transversally and the halves carrying the embryos were used for the resistance test, while the distal halves were used for genetic typing. Biochemical markers were used to discriminate whether the transferred Pch-1 gene was located in chromosome 7D, as is the case for a resistance factor present in Roazon wheat. In the crosses involving Pané 247, resistance was not associated with the 7D locus Pln, which determines sterol ester pattern (dominant allele in H-93-70). In the crosses with the 7D/7Ag substitution line, resistance was neither associated with protein NGE-11 (7D marker), nor alternatively inherited with respect to protein C-7 (7Ag marker). It is concluded that gene Pch-1 represents a different locus and is not an allele of the resistance factor in Roazon wheat.     

Genetic Distance in Housekeeping Genes Between <Emphasis Type="Italic">Plasmodium falciparum</Emphasis> and <Emphasis Type="Italic">Plasmodium reichenowi</Emphasis> and Within <Emphasis Type="Italic">P. falciparum</Emphasis>

The time to the most recent common ancestor of the extant populations of Plasmodium falciparum is controversial. The controversy primarily stems from the limited availability of sequences from Plasmodium reichenowi, a chimpanzee malaria parasite closely related to P. falciparum. Since the rate of nucleotide substitution differs in different loci and DNA regions, the estimation of genetic distance between P. falciparum and P. reichenowi should be performed using orthologous sequences that are evolving neutrally. Here, we obtained full-length sequences of two housekeeping genes, sarcoplasmic and endoplasmic reticulum Ca²⁺-ATPase (serca) and lactate dehydrogenase (ldh), from 11 isolates of P. falciparum and 1 isolate of P. reichenowi and estimate the interspecific genetic distance (divergence) between the two species and intraspecific genetic distance (polymorphism) within P. falciparum. Interspecific distance and intraspecific distance at synonymous sites of interspecies-conserved regions of serca and ldh were 0.0672±0.0088 and 0.0011±0.0007, respectively, using the Nei and Gojobori method. Based on the ratio of interspecific distance to intraspecific distance, the time to the most recent common ancestor of P. falciparum was estimated to be (8.30±5.40) × 10⁴ and (11.62±7.56) × 10⁴ years ago, assuming the divergence time of the two parasite species to be 5 and 7 million years ago, respectively.This article contains an online supplementary table.Reviewing Editor: Dr. Martin Kreitman     

Organization of<Emphasis Type="Italic"> Iroquois</Emphasis> genes in fish

In mammals, a total of six iroquois (Irx) genes exist, which are organized into two clusters. Here we report on the organization of all iroquois genes present in fish, using zebrafish (Danio rerio) and pufferfish (Fugu rubripes and Tetraodon nigroviridis) as examples. A total of 10 Irx genes were found in pufferfish, and 11 in zebrafish; all but one of these genes are organized into clusters (four clusters plus one isolated gene locus). The extra fish clusters result from chromosome duplication in the fish lineage, after its divergence from tetrapod vertebrates. Two of the four fish clusters are highly conserved to the ones in mammals, with regard to similarity of genes and cluster architecture. Irx genes within the other two clusters have diverged in sequence and cluster organization, suggesting functional divergence. These results will allow us to use the zebrafish system for functional and comparative studies of iroquois genes in vertebrate development.Electronic Supplementary Material Supplementary material is available in the online version of this article at Edited by D. Tautz     

Multiple Mutations and Gene Duplications Conferring Organophosphorus Insecticide Resistance Have Been Selected at the Rop-1 Locus of the Sheep Blowfly, Lucilia cuprina

Sequences of the esterase gene E7 were compared across 41 isogenic (IV) strains of the sheep blowfly, Lucilia cuprina, and one strain of the sibling species, L. sericata. The 1.2-kb region sequenced includes sites of two insecticide resistance mutations. Gly137Asp confers resistance to organophosphorus insecticides (OPs), particularly preferring diethyl OPs such as diazinon, while Trp251Leu prefers dimethyl OPs, and particularly malathion, with the additional presence of carboxylester moieties. We found that there are just eight haplotypes among the 41 chromosomes studied: two Gly137Asp containing haplotypes, two Trp251Leu containing haplotypes, and four susceptible haplotypes, including the L. sericata sequence. While phylogenetic analysis of these haplotypes suggests that the Asp137 and Leu251 mutations each arose at least twice, evidence for recombination was detected across the region, therefore single origins for these resistance mutations cannot be ruled out. Levels of linkage disequilibrium in the data are high and significant hitchhiking is indicated by Fay and Wu s H test but not the Tajima test. A test of haplotype diversity indicates a paucity of diversity compared with neutral expectations. Both these results are consistent with a very recent selective sweep at the LcE7 locus. Interestingly, gene duplications of three different combinations of OP resistant haplotypes were identified in seven of the isogenic (IV) strains. All three types of duplication involve an Asp137 and a Trp251 haplotype. To examine whether more haplotypes existed before the hypothesised selective sweep, fragments of E7 surrounding the resistance mutations were amplified from pinned material dating back to before OPs were used. Four new sequence haplotypes, not sampled in the survey of extant haplotypes, were obtained that are all associated with susceptibility. This is suggestive of a higher historical level of susceptible allelic diversity at this locus.Reviewing Editor: Dr. Rasmus Nielsen     

The Carboxylesterase Gene Family from <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Carboxylesterases hydrolyze esters of short-chain fatty acids and have roles in animals ranging from signal transduction to xenobiotic detoxification. In plants, however, little is known of their roles. We have systematically mined the genome from the model plant Arabidopsis thaliana for carboxylesterase genes and studied their distribution in the genome and expression profile across a range of tissues. Twenty carboxylesterase genes (AtCXE) were identified. The AtCXE family shares conserved sequence motifs and secondary structure characteristics with carboxylesterases and other members of the larger / hydrolase fold superfamily of enzymes. Phylogenetic analysis of the AtCXE genes together with other plant carboxylesterases distinguishes seven distinct clades, with an Arabidopsis thaliana gene represented in six of the seven clades. The AtCXE genes are widely distributed across the genome (present in four of five chromosomes), with the exception of three clusters of tandemly duplicated genes. Of the interchromosomal duplication events, two have been mediated through newly identified partial chromosomal duplication events that also include other genes surrounding the AtCXE loci. Eighteen of the 20 AtCXE genes are expressed over a broad range of tissues, while the remaining 2 (unrelated) genes are expressed only in the flowers and siliques. Finally, hypotheses for the functional roles of the AtCXE family members are presented based on the phylogenetic relationships with other plant carboxylesterases of known function, their expression profile, and knowledge of likely esterase substrates found in plants.     

How Common Are Intragene Windows with <Emphasis Type="Italic">K</Emphasis><Subscript>A</Subscript> > <Emphasis Type="Italic">K</Emphasis><Subscript>S</Subscript> Owing to Purifying Selection on Synonymous Mutations?

One method for diagnosing the mode of sequence evolution considers the ratio of nonsynonymous substitutions per nonsynonymous site (K _A) to the corresponding figure for synonymous substitutions (K _S). A ratio (K _A/K _S) greater than unity is taken as evidence for positive selection. This, however, need not necessarily be the case. Notably, there is one instance of a high intragenic K _A/K _S peak, revealed by sliding window analysis and observed in two pairwise comparisons, better accounted for by localised purifying selection on synonymous mutations that affect splicing. Is this example exceptional? To address this we isolate intragenic domains with K _A/K _S > 1 from more than 1000 long mouse-rat orthologues. Approximately one K _A/K _S > 1 peak is found per 12–15 kb of coding sequence. Surprisingly, low synonymous substitution rates underpin more incidences than do high nonsynonymous rates. Several reasons, however, prevent us from supposing that the low synonymous rates reflect purifying selection on synonymous mutations. First, for many peaks, the null that the peak is no higher than expected given the underlying rates of evolution, cannot be rejected. Second, of 18 statistically significant incidences with unusually low K _S values, only 3 are repeatable across independent comparisons. At least two of these are within alternatively spliced exons. We conclude that repeatable statistically significant intragenic domains of low intragenic K _S are rare. As so few K _A/K _S peaks reflect increased rates of protein evolution and so few hold statistical support, we additionally conclude that sliding window analysis to infer domains of positive selection is highly error-prone.     

Mitotic gene conversion, reciprocal recombination and gene replacement at the benA, beta-tubulin, locus of Aspergillus nidulans

Summary We have developed a procedure for determining the rates of mitotic recombination of an interrupted duplication created by integration of transforming plasmid sequences at the benA, beta-tubulin, locus of Aspergillus nidulans. Transformation of a strain carrying a benomyl-resistant benA allele with plasmid AIpGM4, which carries the wild-type benA allele and the pyr4 (orotidine-5-phosphate decarboxylase) gene of Neurospora crassa, creates an interrupted duplication with plasmid sequences flanked by two benA alleles, one wild type and one benomyl resdistant. Such transformants will not grow in the presence of high levels of benomyl. Mitotic recombination causes the loss of the wild-type benA allele or conversion of the wild-type to the mutant allele resulting in nuclei carrying only the benomylresistant allele. Conidia containing such nuclei can be selected on media with high benomyl allowing easy quantitation of mitotic recombination. We found that the rate of recombination giving rise to benomyl-resistant conidia was 4.6×10^-4. Reciprocal recombination leading to benomyl-resistant conidia lacking plasmid sequences occurred at a rate of 2.0×10^-4 and gene conversion leading to benomylresistant conidia occurred at a rate of 2.6×10^-4. We selected for reciprocal recombination leading to loss of pyr4 sequences on 5-fluoro-orotic acid and used this selection for two-step gene replacement of a mutant benA allele with the wild-type allele.