Variation in Constraint Versus Positive Selection as an Explanation for Evolutionary Rate Variation Among Anthocyanin Genes

It has been argued that downstream enzymes in metabolic pathways are expected to be subject to reduced selective constraint, while upstream enzymes, particularly those at pathway branch points, are expected to exhibit more frequent adaptive substitution than downstream enzymes. We examined whether these expectations are met for enzymes in the anthocyanin biosynthetic pathway in Ipomoea. Previous investigations have demonstrated that downstream enzymes in this pathway have substantially higher rates of nonsynonymous substitution than upstream enzymes. We demonstrate here that the difference in rates between the most upstream enzyme (CHS) and the two most downstream enzymes (ANS and UFGT) is explained almost entirely by differences in levels of selective constraint. Adaptive substitutions were not detected in any of these genes. Our results are consistent with suggestions that constraint is greater on enzymes with greater connectivity.     

Positive and negative selection on noncoding DNA close to protein-coding genes in wild house mice

During the past two decades, evidence has accumulated of adaptive evolution within protein-coding genes in a variety of species. However, with the exception of Drosophila and humans, little is known about the extent of adaptive evolution in noncoding DNA. Here, we study regions upstream and downstream of protein-coding genes in the house mouse Mus musculus castaneus, a species that has a much larger effective population size (N(e)) than humans. We analyze polymorphism data for 78 genes from 15 wild-caught M. m. castaneus individuals and divergence to a closely related species, Mus famulus. We find high levels of nucleotide diversity and moderate levels of selective constraint in upstream and downstream regions compared with nonsynonymous sites of protein-coding genes. From the polymorphism data, we estimate the distribution of fitness effects (DFE) of new mutations and infer that most new mutations in upstream and downstream regions behave as effectively neutral and that only a small fraction is strongly negatively selected. We also estimate the fraction of substitutions that have been driven to fixation by positive selection (α) and the ratio of adaptive to neutral divergence (ω(α)). We find that α for upstream and downstream regions (～ 10%) is much lower than α for nonsynonymous sites (～ 50%). However, ω(α) estimates are very similar for nonsynonymous sites (～ 10%) and upstream and downstream regions (～ 5%). We conclude that negative selection operating in upstream and downstream regions of M. m. castaneus is weak and that the low values of α for upstream and downstream regions relative to nonsynonymous sites are most likely due to the presence of a higher proportion of neutrally evolving sites and not due to lower absolute rates of adaptive substitution.     

The evolution of the major histocompatibility complex in upstream versus downstream river populations of the longnose dace

Populations in upstream versus downstream river locations can be exposed to vastly different environmental and ecological conditions and can thus harbor different genetic resources due to selection and neutral processes. An interesting question is how upstream–downstream directionality in rivers affects the evolution of immune response genes. We used next‐generation amplicon sequencing to identify eight alleles of the major histocompatibility complex (MHC) class II β exon 2 in the cyprinid longnose dace (Rhinichthys cataractae) from three rivers in Alberta, upstream and downstream of municipal and agricultural areas along contaminant gradients. We used these data to test for directional and balancing selection on the MHC. We also genotyped microsatellite loci to examine neutral population processes in this system. We found evidence for balancing selection on the MHC in the form of increased nonsynonymous variation relative to neutral expectations, and selection occurred at more amino acid residues upstream than downstream in two rivers. We found this pattern despite no population structure or isolation by distance, based on microsatellite data, at these sites. Overall, our results suggest that MHC evolution is driven by upstream–downstream directionality in fish inhabiting this system.     

The Streptomyces peucetius dpsC Gene Determines the Choice of Starter Unit in Biosynthesis of the Daunorubicin Polyketide

The starter unit used in the biosynthesis of daunorubicin is propionyl coenzyme A (CoA) rather than acetyl-CoA, which is used in the production of most of the bacterial aromatic polyketides studied to date. In the daunorubicin biosynthesis gene cluster of Streptomyces peucetius, directly downstream of the genes encoding the beta-ketoacyl:acyl carrier protein synthase subunits, are two genes, dpsC and dpsD, encoding proteins that are believed to function as the starter unit-specifying enzymes. Recombinant strains containing plasmids carrying dpsC and dpsD, in addition to other daunorubicin polyketide synthase (PKS) genes, incorporate the correct starter unit into polyketides made by these genes, suggesting that, contrary to earlier reports, the enzymes encoded by dpsC and dpsD play a crucial role in starter unit specification. Additionally, the results of a cell-free synthesis of 21-carbon polyketides from propionyl-CoA and malonyl-CoA that used the protein extracts of recombinant strains carrying other daunorubicin PKS genes to which purified DpsC was added suggest that this enzyme has the primary role in starter unit discrimination for daunorubicin biosynthesis.     

The three tricarboxylate synthase activities of <Emphasis Type="Italic">Corynebacterium glutamicum</Emphasis> and increase of <Emphasis Type="SmallCaps">l</Emphasis>-lysine synthesis

Corynebacterium glutamicum owns a citrate synthase and two methylcitrate synthases. Characterization of the isolated enzymes showed that the two methylcitrate synthases have comparable catalytic efficiency, k _cat/K _m, as the citrate synthase with acetyl-CoA as substrate, although these enzymes are only synthesized during growth on propionate-containing media. Thus, the methylcitrate synthases have a relaxed substrate specifity, as also demonstrated by their activity with butyryl-CoA, whereas the citrate synthase does not accept acyl donors other than acetyl-CoA. A double mutant deleted of the citrate synthase gene gltA and one of the methylcitrate synthase genes, prpC1, was made unable to grow on glucose. From this mutant, a collection of suppressor mutants could be isolated which were demonstrated to have regained citrate synthase activity due to the relaxed specificity of the methylcitrate synthase PrpC2. Molecular characterization of these mutants showed that the regulator PrpR (Cg0800) located downstream of prpC1 is mutated with mutations likely to effect the secondary structure of the regulator, thus, resulting in expression of prpC2. This expression results in a citrate synthase activity, which is lower than that due to gltA in the original strain and results in increased l-lysine accumulation.     

Stronger selective constraint on downstream genes in the oxidative phosphorylation pathway of cetaceans

The oxidative phosphorylation (OXPHOS) pathway is an efficient way to produce energy via adenosine triphosphate (ATP), which is critical for sustaining an energy supply for cetaceans in a hypoxic environment. Several studies have shown that natural selection may shape the evolution of the genes involved in OXPHOS. However, how network architecture drives OXPHOS protein sequence evolution remains poorly explored. Here, we investigated the evolutionary patterns of genes in the OXPHOS pathway across six cetacean genomes within the framework of a functional network. Our results show a negative correlation between the strength of purifying selection and pathway position. This result indicates that downstream genes were subjected to stronger evolutionary constraints than upstream genes, which may be due to the dual function of ATP synthase in the OXPHOS pathway. Additionally, there was a positive correlation between codon usage bias and omega (ω = dN/dS) and a negative correlation with synonymous substitution rate (dS), indicating that the stronger selective constraint on genes (with less biased codon usage) along the OXPHOS pathway is attributable to an increase in the rate of synonymous substitution. Surprisingly, there was no significant correlation between protein–protein interactions and the evolutionary estimates, implying that highly connected enzymes may not always show greater evolutionary constraints. Compared with that observed for terrestrial mammals, we found that the signature of positive selection detected in five genes (ATP5J, LHPP, PPA1, UQCRC1 and UQCRQ) was cetacean‐specific, reflecting the importance of OXPHOS for survival in hypoxic, aquatic environments.     

Spermidine biosynthesis in Escherichia coli: promoter and termination regions of the speED operon.

Two enzymes, S-adenosylmethionine decarboxylase and spermidine synthase, are essential for the biosynthesis of spermidine in Escherichia coli. We have previously shown that the genes encoding these enzymes (speD and speE) form an operon and that the area immediately upstream from the speE gene is necessary for the expression of both the speE and speD genes. We have now studied the upstream promoter and the downstream terminator regions of this operon more completely. We have shown that the major mRNA initiation site (Ia) of the operon is located 475 base pairs (bp) upstream from the speE gene and that there is an open reading frame that encodes for a polypeptide of 115 amino acids between the Ia site and the ATG start codon for the speE gene. Downstream from the stop codon for the speD gene is a potential hairpin structure immediately followed by an mRNA termination site, t. An additional mRNA termination site, t', is present about 110 bp downstream from t and is stronger than t. By comparing our DNA fragments with those prepared from this region of the E. coli chromosome by Kohara et al., we have located the speED operon on the physical map of the E. coli chromosome. We have shown that the orientation of the speED operon is counterclockwise and that the operon is located 137.5 to 140 kbp (2.9 minutes) clockwise from the zero position of the E. coli chromosomal map.     

Patterns of evolutionary rate variation among genes of the anthocyanin biosynthetic pathway

The anthocyanin biosynthetic pathway is responsible for the production of anthocyanin pigments in plant tissues and shares a number of enzymes with other biochemical pathways. The six core structural genes of this pathway have been cloned and characterized in two taxonomically diverse plant species (maize and snapdragon). We have recently cloned these genes for a third species, the common morning glory, Ipomoea purpurea. This additional information provides an opportunity to examine patterns of evolution among genes within a single biochemical pathway. We report here that upstream genes in the anthocyanin pathway have evolved substantially more slowly than downstream genes and suggest that this difference in evolutionary rates may be explained by upstream genes being more constrained because they participate in several different biochemical pathways. In addition, regulatory genes associated with the anthocyanin pathway tend to evolve more rapidly than the structural genes they regulate, suggesting that adaptive evolution of flower color may be mediated more by regulatory than by structural genes. Finally, for individual anthocyanin genes, we found an absence of rate heterogeneity among three major angiosperm lineages. This rate constancy contrasts with an accelerated rate of evolution of three CHS-like genes in the Ipomoea lineage, indicating that these three genes have diverged without coordinated adjustment by other pathway genes.     

High-level tuber expression and sucrose inducibility of a potato Sus4 sucrose synthase gene require 5' and 3' flanking sequences and the leader intron.

The 3.6 kb of 5' flanking sequence, leader intron, and 0.7 kb of 3' sequence from the potato sucrose synthase gene Sus4-16 are sufficient to direct high-level expression in developing tubers, in basal tissues of axillary buds and shoots, and in meristems and caps of roots, and to confer sucrose inducibility in leaves. By examining a series of deletion and substitution constructs in transgenic potato plants, we found that this pattern of expression requires 5' flanking sequences both upstream and downstream of position -1500 and that sequences between positions -1500 and -267 are essential for sucrose induction. Replacement of the native 3' sequence with the nopaline synthase 3' sequence resulted in the loss of sucrose inducibility and of expression in basal tissues of axillary buds. A general decrease in expression in other tissues was also observed. Removal of the 1612-bp leader intron also had a dramatic effect on both the pattern and level of expression.     

Relative evolutionary rates of NBS-encoding genes revealed by soybean segmental duplication

It is well known that nucleotide binding site (NBS)-encoding genes are duplicate-rich and fast-evolving genes. However, there is little information on the relative importance of tandem and segmental NBS duplicates and their exact evolutionary rates. The two rounds of large-scale duplication that have occurred in soybean provide a unique opportunity to investigate these issues. Comparison of NBS and non-NBS genes on segments of syntenic homoeologs shows that NBS-encoding genes evolve at least 1.5-fold faster (~1.5-fold higher synonymous and approximately 2.3-fold higher nonsynonymous substitution rates) and lose their genes approximately twofold faster than the flanking non-NBS genes. Compared with segmental duplicates, tandem NBS duplicates are more abundant in soybean, suggesting that tandem duplication is the major driving force in the expansion of NBS genes. Notably, significant sequence exchanges along with significantly positive selection were detected in most tandem-duplicated NBS gene families. The results suggest that the rapid evolution of NBS genes may be due to the combined effects of diversifying selection and frequent sequence exchanges. Interestingly, TIR-NBS-LRR genes (TNLs) have a higher nucleotide substitution rate than non-TNLs, indicating that these types of NBS genes may have a rather different evolutionary pattern. It is important to determine the exact relative evolutionary rates of TNL, non-TNL, and non-NBS genes in order to understand how fast the host plant can adjust its response to rapidly evolving pathogens in a coevolutionary context.     

Evolutionary distance estimation under heterogeneous substitution pattern among lineages

Most of the sophisticated methods to estimate evolutionary divergence between DNA sequences assume that the two sequences have evolved with the same pattern of nucleotide substitution after their divergence from their most recent common ancestor (homogeneity assumption). If this assumption is violated, the evolutionary distance estimated will be biased, which may result in biased estimates of divergence times and substitution rates, and may lead to erroneous branching patterns in the inferred phylogenies. Here we present a simple modification for existing distance estimation methods to relax the assumption of the substitution pattern homogeneity among lineages when analyzing DNA and protein sequences. Results from computer simulations and empirical data analyses for human and mouse genes are presented to demonstrate that the proposed modification reduces the estimation bias considerably and that the modified method performs much better than the LogDet methods, which do not require the homogeneity assumption in estimating the number of substitutions per site. We also discuss the relationship of the substitution and mutation rate estimates when the substitution pattern is not the same in the lineages leading to the two sequences compared.     

Biosynthesis of a natural polyketide-isoprenoid hybrid compound, furaquinocin A: identification and heterologous expression of the gene cluster

Furaquinocin (FQ) A, produced by Streptomyces sp. strain KO-3988, is a natural polyketide-isoprenoid hybrid compound that exhibits a potent antitumor activity. As a first step toward understanding the biosynthetic machinery of this unique and pharmaceutically useful compound, we have cloned an FQ A biosynthetic gene cluster by taking advantage of the fact that an isoprenoid biosynthetic gene cluster generally exists in flanking regions of the mevalonate (MV) pathway gene cluster in actinomycetes. Interestingly, Streptomyces sp. strain KO-3988 was the first example of a microorganism equipped with two distinct mevalonate pathway gene clusters. We were able to localize a 25-kb DNA region that harbored FQ A biosynthetic genes (fur genes) in both the upstream and downstream regions of one of the MV pathway gene clusters (MV2) by using heterologous expression in Streptomyces lividans TK23. This was the first example of a gene cluster responsible for the biosynthesis of a polyketide-isoprenoid hybrid compound. We have also confirmed that four genes responsible for viguiepinol [3-hydroxypimara-9(11),15-diene] biosynthesis exist in the upstream region of the other MV pathway gene cluster (MV1), which had previously been cloned from strain KO-3988. This was the first example of prokaryotic enzymes with these biosynthetic functions. By phylogenetic analysis, these two MV pathway clusters were identified as probably being independently distributed in strain KO-3988 (orthologs), rather than one cluster being generated by the duplication of the other cluster (paralogs).     

Determinants of substitution rates in mammalian genes: expression pattern affects selection intensity but not mutation rate

To determine whether gene expression patterns affect mutation rates and/or selection intensity in mammalian genes, we studied the relationships between substitution rates and tissue distribution of gene expression. For this purpose, we analyzed 2,400 human/rodent and 834 mouse/rat orthologous genes, and we measured (using expressed sequence tag data) their expression patterns in 19 tissues from three development states. We show that substitution rates at nonsynonymous sites are strongly negatively correlated with tissue distribution breadth: almost threefold lower in ubiquitous than in tissue-specific genes. Nonsynonymous substitution rates also vary considerably according to the tissues: the average rate is twofold lower in brain-, muscle-, retina- and neuron-specific genes than in lymphocyte-, lung-, and liver-specific genes. Interestingly, 5' and 3' untranslated regions (UTRs) show exactly the same trend. These results demonstrate that the expression pattern is an essential factor in determining the selective pressure on functional sites in both coding and noncoding regions. Conversely, silent substitution rates do not vary with expression pattern, even in ubiquitously expressed genes. This latter result thus suggests that synonymous codon usage is not constrained by selection in mammals. Furthermore, this result also indicates that there is no reduction of mutation rates in genes expressed in the germ line, contrary to what had been hypothesized based on the fact that transcribed DNA is more efficiently repaired than nontranscribed DNA.     

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae)

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.     

Coordinate genetic control of yeast fatty acid synthase genes FAS1 and FAS2 by an upstream activation site common to genes involved in membrane lipid biosynthesis.

A systematic search for upstream controlling elements necessary for efficient expression of the yeast fatty acid synthase genes FAS1 and FAS2 revealed identical activation sites, UASFAS, in front of both FAS genes. The individual element confers, in a heterologous yeast test system, an approximately 40-fold stimulation of basal gene expression. The UASFAS motifs identified have the consensus sequence TYTTCACATGY and function in either orientation. The same sequence motif is found in the upstream regions of all so far characterized yeast genes encoding enzymes of phospholipid biosynthesis. In gel retardation assays, a protein factor, Fbf1 (FAS binding factor), was identified which interacted with UASFAS. The UASFAS motif proved to be an inositol/choline responsive element (ICRE) conferring strict repression by exogenous inositol and choline on a heterologous reporter gene. Its core sequence perfectly matches the CANNTG motif typical of basic helix-loop-helix DNA-binding proteins. In contrast to the individual UASFAS element, the intact yeast FAS promoters are not significantly influenced by inositol and choline, and thus allow nearly constitutive fatty acid synthase production. Available evidence suggests that additional cis- and trans-acting elements, other than UASFAS and Fbf1, are involved in this constitutive FAS gene expression.     

Evidence of selection at melanin synthesis pathway loci during silkworm domestication

The domesticated silkworm (Bombyx mori) was domesticated from wild silkworm (Bombyx mandarina) more than 5,000 years ago. During domestication, body color between B. mandarina and B. mori changed dramatically. However, the molecular mechanism of the silkworm body color transition is not known. In the present study, we examined within- and between-species nucleotide diversity for eight silkworm melanin synthesis pathway genes, which play a key role in cuticular pigmentation of insects. Our results showed that the genetic diversity of B. mori was significantly lower than that of B. mandarina and 40.7% of the genetic diversity of wild silkworm was lost in domesticated silkworm. We also examined whether position effect exists among melanin synthesis pathway genes in B. mandarina and B. mori. We found that the upstream genes have significantly lower levels of genetic diversity than the downstream genes, supporting a functional constraint hypothesis (FCH) of metabolic pathway, that is, upstream enzymes are under greater selective constraint than downstream enzymes because upstream enzymes participate in biosynthesis of a number of metabolites. We also investigated whether some of the melanin synthesis pathway genes experienced selection during domestication. Neutrality test, coalescent simulation, as well as network and phylogenetic analyses showed that tyrosine hydroxylase (TH) gene was a domestication locus. Sequence analysis further suggested that a putative expression enhancer (Abd-B-binding site) in the intron of TH gene might be disrupted during domestication. TH is the rate-limiting enzyme of melanin synthesis pathway in insects. Real-time polymerase chain reaction assay did show that the relative expression levels of TH gene in B. mori were significantly lower than that in B. mandarina at three different developmental stages, which is consistent with light body color of domesticated silkworm relative to wild silkworm. Therefore, we speculated that expression change of TH gene may contribute to the body color transition from B. mandarina to B. mori. Our results emphasize the exceptional role of gene expression regulation in morphological transition of domesticated animals.