Effects of gene duplication, epistasis, recombination and gene conversion on the fixation time of compensatory mutations.

Gene duplication is one of the major mechanisms of molecular evolution. Gene duplication enables copies of a gene to accumulate mutations through functional redundancy. If a gene encodes a specific protein that interacts with other proteins, RNA, or DNA, the relaxation of selective constraints caused by gene duplication might contribute to the fixation of compensatory mutations that occur at the interacting sites. In this study, we investigate the effect of gene duplication, epistasis among the duplicated copies and gene conversion on the fixation time of compensatory mutations by extending the original model of compensatory evolution proposed by Kimura in 1985. Our simulation results reveal that the time to fixation of compensatory mutations can be decreased remarkably by gene duplication if one of the duplicated loci can completely mask the deleterious effects of a mutation that occurs at the other locus. Conversely, the fixation time can be increased by gene duplication if such functional compensation is weak. We also show that the combination of the degree of functional compensation and the rate of gene conversion between duplicate loci have contrasting effects on the time to fixation of compensatory mutations.

A model of compensatory molecular evolution involving multiple sites in RNA molecules.

Consider two sites under compensatory fitness interaction, such as a Watson-Crick base pair in an RNA helix or two interacting residues in a protein. A mutation at any one of these two sites may reduce the fitness of an individual. However, fitness may be restored by the occurrence of a second mutation at the other site. Kimura modeled this process using a two-locus haploid fitness scheme with two alleles at each locus. He predicted that compensatory evolution following this model is very rare unless selection against the deleterious single mutations is weak and linkage between the interacting sites is tight. Here we investigate the question whether the rate of compensatory evolution increases if we take the context of the two directly interacting sites into account. By "context", we mean the effect of neighboring sites in an RNA helix. Interaction between the focal pair of sites under consideration and the context may lead to so-called indirect compensation. Thus, extending Kimura's classical model of compensatory evolution, we study the effects of both direct and indirect compensation on the rate of compensatory evolution. It is shown that the effects of indirect compensation are very strong. We find that recombination does not slow down the rate of compensatory evolution as predicted by the classical model. Instead, compensatory substitutions may be relatively frequent, even if linkage between the focal interacting sites is loose, selection against deleterious mutations is strong, and mutation rate is low. We compare our theoretical results with data on RNA secondary structures from vertebrate introns.

Models of compensatory molecular evolution: effects of back mutation.

Compensatory mutations are individually deleterious but appropriate combinations of mutants are harmless. For several models of compensatory molecular evolution, we consider the effects of back mutation. It is shown that the effects of back mutation on the rate of compensatory molecular evolution are weak. Further we estimate the values of selection parameter of deleterious single mutants for the models of compensatory molecular evolution both with and without back mutation using sequence data of folded RNA molecules and compare them with the previous results.

Effective population size of a population with stochastically varying size.

For a Wright-Fisher model with mutation whose population size fluctuates stochastically from generation to generation, a heterozygosity effective population size is defined by means of the equilibrium average heterozygosity of the population. It is shown that this effective population size is equal to the harmonic mean of population size if and only if the stochastic changes of population size are uncorrelated. The effective population size is larger (resp. smaller) than the harmonic mean when the stochastic changes of population size are positively (resp. negatively) autocorrelated. These results and those obtained so far for other stochastic models with fluctuating population size suggest that the property that effective population sizes are always larger than the harmonic mean under the fluctuation of population size holds only for continuous time models such as diffusion and coalescent models, whereas effective population sizes can be equal to or smaller than the harmonic mean for discrete time models.

Coordinated expression of functionally diverse fructosyltransferase genes is associated with fructan accumulation in response to low temperature in perennial ryegrass.

* Fructan is the major nonstructural carbohydrate reserve in temperate grasses. To understand regulatory mechanisms in fructan synthesis and adaptation to cold environments, the isolation, functional characterization and genetic mapping of fructosyltransferase (FT) genes in perennial ryegrass (Lolium perenne) are described. * Six cDNAs (prft1-prft6) encoding FTs were isolated from cold-treated ryegrass plants, and three were positioned on a perennial ryegrass linkage map. Recombinant proteins were produced in Pichia pastoris and enzymatic activity was characterized. Changes in carbohydrate levels and mRNA levels of FT genes during cold treatment were also analysed. * One gene encodes sucrose-sucrose 1-fructosyltransferase (1-SST), and two gene encode fructan-fructan 6G-fructosyltransferase (6G-FFT). Protein sequences for the other genes (prfts 1, 2 and 6) were similar to sucrose-fructan 6-fructosyltransferase (6-SFT). The 1-SST and prft1 genes were colocalized with an invertase gene on the ryegrass linkage map. The mRNA levels of prft1 and prft2 increased gradually during cold treatment, while those of the 1-SST and 6G-FFT genes first increased, but then decreased before increasing again during a longer period of cold treatment. * Thus at least two different patterns of gene expression have developed during the evolution of functionally diverse FT genes, which are associated in a coordinated way with fructan synthesis in a cold environment.

Compensatory evolution in diploid populations.

Compensatory mutations are individually deleterious but harmless in appropriate combinations either at more than two sites within a gene or on separate genes. Considering that dominance effects of selection and heterodimer formation of gene products may affect the rate of compensatory evolution, we investigate compensatory neutral mutation models for diploid populations. Our theoretical analysis on the average time until fixation of compensatory mutations shows that these factors play an important role in reducing the fixation time of compensatory mutations if mutation rates are not low. Compensatory evolution of heterodimers is shown to occur more easily if the deleterious effects of single mutants are recessive.

A neutral model with fluctuating population size and its effective size.

We consider a diffusion model with neutral alleles whose population size is fluctuating randomly. For this model, the effects of fluctuation of population size on the effective size are investigated. The effective size defined by the equilibrium average heterozygosity is larger than the harmonic mean of population size but smaller than the arithmetic mean of population size. To see explicitly the effects of fluctuation of population size on the effective size, we investigate a special case where population size fluctuates between two distinct states. In some cases, the effective size is very different from the harmonic mean. For this concrete model, we also obtain the stationary distribution of the average heterozygosity. Asymptotic behavior of the effective size is obtained when the population size is large and/or autocorrelation of the fluctuation is weak or strong.

Characterization of a thermostable levansucrase from Bacillus sp. TH4-2 capable of producing high molecular weight levan at high temperature.

A thermoactive and thermostable levansucrase was purified from a newly isolated thermophilic Bacillus sp. from Thailand soil. The purification was achieved by alcohol precipitation, DEAE-Cellulose and gel filtration chromatographies. The enzyme was purified to homogeneity as determined by SDS-PAGE, and had a molecular mass of 56 kDa. This levansucrase has some interesting characteristics regarding its optimum temperature and heat stability. The optimum temperature and pH were 60 degrees C and 6.0, respectively. The enzyme was completely stable after treatment at 50 degrees C for more than 1 h, and its activity increased four folds in the presence of 5 mM Fe(2+). The optimum temperature for levan production was 50 degrees C. Contrary to other levansucrases, the one presented in this study is able to produce high molecular weight levan at 50 degrees C.

A population genetic study on the transition from Jomon people to Yayoi people.

In the study on the origin of Japanese, one of main unsolved problems is the transition from the Jomon people to the Yayoi people. The main difficulty in solving this problem has been the lack of suitable skeletal materials belonging to the time between the two periods, i.e. the final Jomon and the early Yayoi Periods. Therefore, we know few details of the transition period. It is important to know who carried out a drastic change of the Yayoi culture during this transitional period, i.e. the native Jomon people or the immigrant people. By introducing population genetic models, we show that a view that the immigrant people had a significant genetic contribution to the origin of Japanese is compatible with results from anthropological and archeological studies. This result implies that the immigrant people were mainly responsible for the drastic cultural change during the transitional period.

Degradation of raw starch granules by alpha-amylase purified from culture of Aspergillus awamori KT-11.

Raw-starch-digesting alpha-amylase (Amyl III) was purified to an electrophoretically pure state from the extract of a koji culture of Aspergillus awamori KT-11 using wheat bran in the medium. The purified Amyl III digested not only soluble starch but also raw corn starch. The major products from the raw starch using Amyl III were maltotriose and maltose, although a small amount of glucose was produced. Amyl III acted on all raw starch granules that it has been tested on. However, it was considered that the action mode of the Amyl III on starch granules was different from that of glucoamylase judging from the observation of granules under a scanning electron microscope before and after enzyme reaction, and also from the reaction products. Glucoamylase (GA I) was also isolated and it was purified to an electrophoretically pure state from the extract. It was found that the electron micrographic features of the granules after treatment with the enzymes were quite different. A synergistic effect of Amyl III and GA I was observed for the digestion of raw starch granules.

Molecular cloning and determination of the nucleotide sequence of raw starch digesting alpha-amylase from Aspergillus awamori KT-11.

Complementary DNAs encoding alpha-amylases (Amyl I, Amyl III) and glucoamylase (GA I) were cloned from Aspergillus awamori KT-11 and their nucleotide sequences were determined. The sequence of Amyl III that was a raw starch digesting alpha-amylase was found to consist of a 1,902 bp open reading frame encoding 634 amino acids. The signal peptide of the enzyme was composed of 21 amino acids. On the other hand, the sequence of Amyl I, which cannot act on raw starch, consisted of a 1,500 bp ORF encoding 499 amino acids. The signal peptide of the enzyme was composed of 21 amino acids. The sequence of GA I consisted of a 1,920 bp ORF that encoded 639 amino acids. The signal peptide was composed of 24 amino acids. The amino acid sequence of Amyl III from the N-terminus to the amino acid number 499 showed 63.3% homology with Amyl I. However, the amino acid sequence from the amino acid number 501 to C-terminus, including the raw-starch-affinity site and the TS region rich in threonine and serine, showed 66.9% homology with GA I.

New action pattern of a maltose-forming alpha-amylase from Streptomyces sp. and its possible application in bakery.

An a-Amylase (EC 3.2.1.1) was purified that catalyses the production of a high level of maltose from starch without the attendant production of glucose. The enzyme was produced extracellularly by thermophilic Streptomyces sp. that was isolated from Thailand's soil. Purification was achieved by alcohol precipitation, DEAE-Cellulose, and Gel filtration chromatographies. The purified enzyme exhibited maximum activity at pH 6-7 and 60 degrees C. It had a relative molecular mass of 45 kDa, as determined by SDS-PAGE. The hydrolysis products from starch had alpha-anomeric forms, as determined by 1H-NMR. This maltose-forming alpha-Amylase completely hydrolyzed the soluble starch to produce a high level of maltose, representing up to 90%. It hydrolyzed maltotetrose and maltotriose to primarily produce maltose (82% and 62% respectively) without the attendant production of glucose. The high maltose level as a final end-product from starch and maltooligosaccharides, and the unique action pattern of this enzyme, indicate an unusual maltose-forming system. After the addition of the enzyme in the bread-baking process, the bread's volume increased and kept its softness longer than when the bread had no enzyme.

High expression level of levansucrase from Bacillus licheniformis RN-01 and synthesis of levan nanoparticles.

LsRN from Bacillus licheniformis was cloned and expressed in Escherichia coli. From a 1793 bp genomic sequence, the lsRN gene was found to be composed of a single 1446 bp ORF with a putative promoter consensus boxes and a ribosome-binding site. This ORF was predicted to encode for 482 amino acid residues. The LsRN was constitutively expressed at a relatively high level without sucrose induction. The enzyme was highly purified and an apparent size of 52 kDa with an optimum temperature and pH of 50 °C and 6.0 were determined. The wide range of M(w) of levan (1-600 kDa) was synthesized in a controlled reaction with two variable parameters: temperature and ionic strength. At high temperature (50 °C), LsRN synthesized high M(w) levan (612 kDa) as a major product while at low temperature (30 °C), low M(w) levan (11 kDa) was mainly synthesized. When 0.5M NaCl was added into the reaction, the major products at both temperatures were of the size 11 kDa. Moreover we report for the first time, an enzymatic synthesis of levan nanoparticles (NPs) by a single step reaction. The LsRN synthesized levan NPs as agglomerate with average particle size of 50 nm. The encapsulation of O-acetyl-α-tocopherol was carried out to demonstrate the applicable use of levan NPs.

Evidence for the transglycosylation of complex type oligosaccharides of glycoproteins by endo-beta-N-acetylglucosaminidase HS.

Transglycosylation activity of endo-beta-N-acetylglucosaminidase HS (Endo HS) was investigated using native human transferrin as a donor of an asparagine-linked oligosaccharide and p-nitrophenyl-beta-d-glucose (PNP-beta-d-Glc) as an acceptor of the oligosaccharide. The amount of the product increased dependent on the concentration of the acceptors. Absorption spectrum, exoglycosidase digestion and matrix assisted laser desorption and ionization-time of flight (MALDI-TOF) mass analysis of the transglycosylation product indicated that the asialobiantennary complex type oligosaccharide of human transferrin was transferred to PNP-beta-d-Glc. Endo HS also transferred the oligosaccharide of human transferrin to PNP-alpha-d-Glc, PNP-alpha-d-Gal, PNP-beta-d-Gal, PNP-beta-d-Man, PNP-beta-d-Xyl, PNP-beta-d-GlcNAc, and PNP-glycerol at a different rate. No apparent difference in the K(m) value for human transferrin as an oligosaccharide donor was observed using different acceptors, PNP-beta-d-Glc and PNP-glycerol. The amount of the transglycosylation product successively increased and became constant and then very slightly decreased during the course of enzyme reaction. Endo HS was also transferred the triantennary complex type oligosaccharide of calf fetuin and the bi-, tri-, and tetrantennary complex type oligosaccharides of human alpha(1)-acid glycoprotein to PNP-beta-d-Glc. Furthermore, Endo HS transferred an asparagine-linked oligosaccharide from a hen egg glycopeptide to PNP-beta-d-Glc. The results demonstrate that Endo HS can transfer a wide variety of asparagine-linked complex type oligosaccharides to various monosaccharides. Endo HS was distinct from other enzymes in the specificity for oligosaccharide donors and acceptors.

Coalescent process with fluctuating population size and its effective size.

We consider a Wright-Fisher model whose population size is a finite Markov chain. We introduce a sequence of two-dimensional discrete time Markov chains whose components describe the coalescent process and the fluctuation of population size. For the limiting process of the sequence of Markov chains, the relationship of the expectation of coalescence time to the harmonic and the arithmetic means of population sizes is shown, and the Laplace transform of the distribution of coalescence time is calculated. We define the coalescence effective population size (cEPS) by the expectation of coalescence time. We show that cEPS is strictly larger (resp. smaller) than the harmonic (resp. arithmetic) mean. As the population size fluctuates more quickly (resp. slowly), cEPS is closer to the harmonic (resp. arithmetic) mean. For the case of a two-valued Markov chain, we show the explicit expression of cEPS and its dependency on the sample size.