Clenshaw–Curtis quadrature

A number of different Markov models of DNA sequence evolution have been proposed. These substitution models differ in terms of the parameters used to describe the rates at which one nucleotide replaces another during evolution. These models are frequently used in molecular phylogenetic analyses. In particular, they are used during the calculation of likelihood of a tree (in Bayesian and maximum likelihood approaches to tree estimation) and they are used to estimate the evolutionary distance between sequences from the observed differences between the sequences.

Introduction

These models are phenomenological descriptions of the evolution of DNA as a string of four discrete states. These Markov models do not explicitly depict the mechanism of mutation nor the action of natural selection. Rather they describe the relative rates of different changes. For example, mutational biases and purifying selection favoring conservative changes are probably both responsible for the relatively high rate of transitions compared to transversions in evolving sequences. However, the Kimura (K80) model described below merely attempts to capture the effect of both forces in a parameter that reflects the relative rate of transitions to transversions.

Evolutionary analyses of sequences are conducted on a wide variety of time scales. Thus, it is convenient to express these models in terms of the instantaneous rates of change between different states (the Q matrices below). If we are given a starting (ancestral) state at one position, the model's Q matrix and a branch length expressing the expected number of changes to have occurred since the ancestor, then we can derive the probability of the descendant sequence having each of the four states. The mathematical details of this transformation from rate-matrix to probability matrix are described in the mathematics of substitution models section of the substitution model page. By expressing models in terms of the instantaneous rates of change we can avoid estimating a large numbers of parameters for each branch on a phylogenetic tree (or each comparison if the analysis involves many pairwise sequence comparisons).

The models described on this page describe the evolution of a single site within a set of sequences. They are often used for analyzing the evolution of an entire locus by making the simplifying assumption that different sites evolve independently and are identically distributed. This assumption may be justifiable if the sites can be assumed to be evolving neutrally. If the primary effect of natural selection on the evolution of the sequences is to constrain some sites, then models of among-site rate-heterogeneity can be used. This approaches allows one to estimate only one matrix of relative rates of substitution, and another set of parameters describing the variance in the total rate of substitution across sites.

DNA evolution as a continuous-time Markov chain

Continuous-time Markov chains

Continuous-time Markov chains have the usual transition matrices which are, in addition, parameterized by time, $t$. Specifically, if $E_1,\dots ,E_4$ are the states, then the transition matrix

$P(t)=(P_{ij}(t))$ where each individual entry, $P_{ij}(t)$ refers to the probability that state $E_i$ will change to state $E_j$ in time $t$.

Example: We would like to model the substitution process in DNA sequences (i.e. Jukes–Cantor, Kimura, etc.) in a continuous-time fashion. The corresponding transition matrices will look like:

$P(t)=\left(\begin{array}{cccc}{p}_{AA}(t)& p_{GA}(t)& p_{CA}(t)& p_{TA}(t)\\ p_{AG}(t)& p_{GG}(t)& p_{CG}(t)& p_{TG}(t)\\ p_{AC}(t)& p_{GC}(t)& p_{CC}(t)& p_{TC}(t)\\ p_{AT}(t)& p_{GT}(t)& p_{CT}(t)& p_{TT}(t)\end{array}\right)$

where the top-left and bottom-right 2 × 2 blocks correspond to transition probabilities and the top-right and bottom-left 2 × 2 blocks corresponds to transversion probabilities.

Assumption: If at some time $t_0$, the Markov chain is in state $E_i$, then the probability that at time $t_0+t$, it will be in state $E_j$ depends only upon $i$, $j$ and $t$. This then allows us to write that probability as $p_{ij}(t)$.

Theorem: Continuous-time transition matrices satisfy:

$P(t+\tau )=P(t)P(\tau )$

Note: There is here a possible confusion between two meanings of the word transition. (i) In the context of Markov chains, transition is the general term that refers to the change between two states. (ii) In the context of nucleotide changes in DNA sequences, transition is a specific term that refers to the exchange between either the two purines (A ↔ G) or the two pyrimidines (C ↔ T) (for additional details, see the article about transitions in genetics). By contrast, an exchange between one purine and one pyrimidine is called a transversion.

Deriving the dynamics of substitution

Consider a DNA sequence of fixed length m evolving in time by base replacement. Assume that the processes followed by the m sites are Markovian independent, identically distributed and constant in time. For a fixed site, let

$\mathbf{P}(t)=(p_A(t),p_G(t),p_C(t),p_T(t){)}^T$

be the column vector of probabilities of states $A,$ $G,$ $C,$ and $T$ at time $t$. Let

$\mathcal{ℰ}=\{A,G,C,T\}$

be the state-space. For two distinct

$x,y\in \mathcal{ℰ}$, let ${\mu }_{xy}$

be the transition rate from state $x$ to state $y$. Similarly, for any $x$, let:

${\mu }_x=\sum _{y\ne x}{\mu }_{xy}$

The changes in the probability distribution $p_A(t)$ for small increments of time $\mathrm{\Delta }t$ are given by:

$p_A(t+\mathrm{\Delta }t)=p_A(t)-p_A(t){\mu }_A\mathrm{\Delta }t+\sum _{x\ne A}{p}_x(t){\mu }_{xA}\mathrm{\Delta }t$

In other words (in frequentist language), the frequency of $A$'s at time $t+\mathrm{\Delta }t$ is equal to the frequency at time $t$ minus the frequency of the lost $A$'s plus the frequency of the newly created $A$'s.

Similarly for the probabilities $p_G(t),p_C(t),\mathrm{a}\mathrm{n}\mathrm{d}{p}_T(t)$. We can write these compactly as:

$\mathbf{P}(t+\mathrm{\Delta }t)=\mathbf{P}(t)+Q\mathbf{P}(t)\mathrm{\Delta }t$

where,

$Q=\left(\begin{array}{cccc}-{\mu }_A& {\mu }_{GA}& {\mu }_{CA}& {\mu }_{TA}\\ {\mu }_{AG}& -{\mu }_G& {\mu }_{CG}& {\mu }_{TG}\\ {\mu }_{AC}& {\mu }_{GC}& -{\mu }_C& {\mu }_{TC}\\ {\mu }_{AT}& {\mu }_{GT}& {\mu }_{CT}& -{\mu }_T\end{array}\right)$

or, alternately:

${\mathbf{P}}^{\prime }(t)=Q\mathbf{P}(t)$

where, $Q$ is the rate matrix. Note that by definition, the columns of $Q$ sum to zero.

Ergodicity

If all the transition probabilities, ${\mu }_{xy}$ are positive, i.e. if all states $x,y\in \mathcal{ℰ}$ communicate, then the Markov chain has a stationary distribution $\mathrm{\Pi }=\{{\pi }_x,x\in \mathcal{ℰ}\}$ where each ${\pi }_x$ is the proportion of time spent in state $x$ after the Markov chain has run for infinite time, and this probability does not depend upon the initial state of the process. Such a Markov chain is called, ergodic. In DNA evolution, under the assumption of a common process for each site, the stationary frequencies, ${\pi }_A,{\pi }_G,{\pi }_C,{\pi }_T$ correspond to equilibrium base compositions.

Definition A Markov process is stationary if its current distribution is the stationary distribution, i.e. $\mathbf{P}(t)=\mathrm{\Pi }$

Thus, by using the differential equation above:

$\frac{d\mathrm{\Pi }}{dt}=Q\mathrm{\Pi }=0$

Time reversibility

Definition: A stationary Markov process is time reversible if (in the steady state) the amount of change from state $x$ to $y$ is equal to the amount of change from $y$ to $x$, (although the two states may occur with different frequencies). This means that:

${\pi }_x{\mu }_{xy}={\pi }_y{\mu }_{yx}$

Not all stationary processes are reversible, however, almost all DNA evolution models assume time reversibility, which is considered to be a reasonable assumption.

Under the time reversibility assumption, let $s_{xy}={\mu }_{xy}/{\pi }_y$, then it is easy to see that:

$s_{xy}=s_{yx}$

Definition The symmetric term $s_{xy}$ is called the exchangeability between states $x$ and $y$. In other words, $s_{xy}$ is the fraction of the frequency of state $x$ that results as a result of transitions from state $y$ to state $x$.

Corollary The 12 off-diagonal entries of the rate matrix, $Q$ (note the off-diagonal entries determine the diagonal entries, since the rows of $Q$ sum to zero) can be completely determined by 9 numbers; these are: 6 exchangeability terms and 3 stationary frequencies ${\pi }_x$, (since the stationary frequencies sum to 1).

Scaling of branch lengths

By comparing extant sequences, one can determine the amount of sequence divergence. This raw measurement of divergence provides information about the number of changes that have occurred along the path separating the sequences. The simple count of differences (the Hamming distance) between sequences will often underestimate the number of substitution because of multiple hits (see homoplasy). Trying to estimate the exact number of changes that have occurred is difficult, and usually not necessary. Instead, branch lengths (and path lengths) in phylogenetic analyses are usually expressed in the expected number of changes per site. The path length is the product of the duration of the path in time and the mean rate of substitutions. While their product can be estimated, the rate and time are not identifiable from sequence divergence.

The descriptions of rate matrices on this page accurately reflect the relative magnitude of different substitutions, but these rate matrices are not scaled such that a branch length of 1 yields one expected change. This scaling can be accomplished by multiplying every element of the matrix by the same factor, or simply by scaling the branch lengths. If we use the β to denote the scaling factor, and ν to denote the branch length measured in the expected number of substitutions per site then βν is used the transition probability formulae below in place of μt. Note that ν is a parameter to be estimated from data, and is referred to as the branch length, while β is simply a number that can be calculated from the rate matrix (it is not a separate free parameter).

The value of β can be found by forcing the expected rate of flux of states to 1. The diagonal entries of the rate-matrix (the Q matrix) represent -1 times the rate of leaving each state. For time-reversible models, we know the equilibrium state frequencies (these are simply the π_i parameter value for state i). Thus we can find the expected rate of change by calculating the sum of flux out of each state weighted by the proportion of sites that are expected to be in that class. Setting β to be the reciprocal of this sum will guarantee that scaled process has an expected flux of 1:

$\beta =1/(-\sum _i{\pi }_i{\mu }_{ii})$

For example, in the Jukes-Cantor, the scaling factor would be 4/(3μ)' because the rate of leaving each state is 3μ/4.

Most common models of DNA evolution

JC69 model (Jukes and Cantor, 1969)

JC69 is the simplest substitution model. There are several assumptions. It assumes equal base frequencies $({\pi }_A={\pi }_G={\pi }_C={\pi }_T=\frac{1}{4})$ and equal mutation rates. The only parameter of this model is therefore $\mu$, the overall substitution rate. As previously mentioned, this variable becomes a constant when we normalize to the mean-rate to 1.

$Q=\left(\begin{array}{cccc}*& \frac{\mu }{4}& \frac{\mu }{4}& \frac{\mu }{4}\\ \frac{\mu }{4}& *& \frac{\mu }{4}& \frac{\mu }{4}\\ \frac{\mu }{4}& \frac{\mu }{4}& *& \frac{\mu }{4}\\ \frac{\mu }{4}& \frac{\mu }{4}& \frac{\mu }{4}& *\end{array}\right)$

$P=\left(\begin{array}{cccc}\frac{1}{4}+\frac{3}{4}{e}^{-t\mu }& \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }& \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }& \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }\\ \\ \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }& \frac{1}{4}+\frac{3}{4}{e}^{-t\mu }& \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }& \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }\\ \\ \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }& \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }& \frac{1}{4}+\frac{3}{4}{e}^{-t\mu }& \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }\\ \\ \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }& \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }& \frac{1}{4}-\frac{1}{4}{e}^{-t\mu }& \frac{1}{4}+\frac{3}{4}{e}^{-t\mu }\end{array}\right)$

When branch length, $\nu$, is measured in the expected number of changes per site then:

$P_{ij}(\nu )=\{\begin{array}{cc}\frac{1}{4}+\frac{3}{4}{e}^{-4\nu /3}& \text{ if }i=j\\ \frac{1}{4}-\frac{1}{4}{e}^{-4\nu /3}& \text{ if }i\ne j\end{array}$

It is worth notice that $\nu =\frac{3}{4}t\mu =(\frac{\mu }{4}+\frac{\mu }{4}+\frac{\mu }{4})t$ what stands for sum of any column (or row) of matrix $Q$ multiplied by time and thus means expected number of substitutions in time $t$ (branch duration) for each particular site (per site) when the rate of substitution equals $\mu$.

Given the proportion $p$ of sites that differ between the two sequences the Jukes-Cantor estimate of the evolutionary distance (in terms of the expected number of changes) between two sequences is given by

$\hat{d}=-\frac{3}{4}\ln (1-\frac{4}{3}p)=\hat{\nu }$

The $p$ in this formula is frequently referred to as the $p$-distance. It is a sufficient statistic for calculating the Jukes-Cantor distance correction, but is not sufficient for the calculation of the evolutionary distance under the more complex models that follow (also note that $p$ used in subsequent formulae is not identical to the "$p$-distance").

K80 model (Kimura, 1980)

The K80 model distinguishes between transitions (A <-> G, i.e. from purine to purine, or C <-> T, i.e. from pyrimidine to pyrimidine) and transversions (from purine to pyrimidine or vice versa). In Kimura's original description of the model the α and β were used to denote the rates of these types of substitutions, but it is now more common to set the rate of transversions to 1 and use κ to denote the transition/transversion rate ratio (as is done below). The K80 model assumes that all of the bases are equally frequent (π_T=π_C=π_A=π_G=0.25).

Rate matrix $Q=\left(\begin{array}{cccc}*& \kappa & 1& 1\\ \kappa & *& 1& 1\\ 1& 1& *& \kappa \\ 1& 1& \kappa & *\end{array}\right)$

The Kimura two-parameter distance is given by:

$\hat{d}=-\frac{1}{2}\ln (1-2p-q)-\frac{1}{4}\ln (1-2q)$

where p is the proportion of sites that show transitional differences and q is the proportion of sites that show transversional differences.

F81 model (Felsenstein 1981)

Felsenstein's 1981 model is an extension of the JC69 model in which base frequencies are allowed to vary from 0.25 (${\pi }_T\ne {\pi }_C\ne {\pi }_A\ne {\pi }_G$)

Rate matrix:

$Q=\left(\begin{array}{cccc}*& {\pi }_C& {\pi }_A& {\pi }_G\\ {\pi }_T& *& {\pi }_A& {\pi }_G\\ {\pi }_T& {\pi }_C& *& {\pi }_G\\ {\pi }_T& {\pi }_C& {\pi }_A& *\end{array}\right)$

When branch length, ν, is measured in the expected number of changes per site then:

$\beta =1/(1-{\pi }_A^2-{\pi }_C^2-{\pi }_G^2-{\pi }_T^2)$

$P_{ij}(\nu )=\{\begin{array}{cc}{e}^{-\beta \nu }+{\pi }_j(1-e^{-\beta \nu })& \text{ if }i=j\\ {\pi }_j(1-e^{-\beta \nu })& \text{ if }i\ne j\end{array}$

HKY85 model (Hasegawa, Kishino and Yano 1985)

The HKY85 model can be thought of as combining the extensions made in the Kimura80 and Felsenstein81 models. Namely, it distinguishes between the rate of transitions and transversions (using the κ parameter), and it allows unequal base frequencies (${\pi }_T\ne {\pi }_C\ne {\pi }_A\ne {\pi }_G$). [ Felsenstein described a similar (but not equivalent) model in 1984 using a different parameterization; that latter model is referred to as the F84 model. ]

Rate matrix $Q=\left(\begin{array}{cccc}*& \kappa {\pi }_C& {\pi }_A& {\pi }_G\\ \kappa {\pi }_T& *& {\pi }_A& {\pi }_G\\ {\pi }_T& {\pi }_C& *& \kappa {\pi }_G\\ {\pi }_T& {\pi }_C& \kappa {\pi }_A& *\end{array}\right)$

If we express the branch length, ν in terms of the expected number of changes per site then:

$\beta =\frac{1}{2({\pi }_A+{\pi }_G)({\pi }_C+{\pi }_T)+2\kappa [({\pi }_A{\pi }_G)+({\pi }_C{\pi }_T)]}$

$P_{AA}(\nu ,\kappa ,\pi )=[{\pi }_A({\pi }_A+{\pi }_G+({\pi }_C+{\pi }_T)e^{-\beta \nu })+{\pi }_G{e}^{-(1+({\pi }_A+{\pi }_G)(\kappa -1.0))\beta \nu }]/({\pi }_A+{\pi }_G)$

$P_{AC}(\nu ,\kappa ,\pi )={\pi }_C(1.0-e^{-\beta \nu })$

$P_{AG}(\nu ,\kappa ,\pi )=[{\pi }_G({\pi }_A+{\pi }_G+({\pi }_C+{\pi }_T)e^{-\beta \nu })-{\pi }_G{e}^{-(1+({\pi }_A+{\pi }_G)(\kappa -1.0))\beta \nu }]/({\pi }_A+{\pi }_G)$

$P_{AT}(\nu ,\kappa ,\pi )={\pi }_T(1.0-e^{-\beta \nu })$

and formula for the other combinations of states can be obtained by substituting in the appropriate base frequencies.

T92 model (Tamura 1992)

T92 is a simple mathematical method developed to estimate the number of nucleotide substitutions per site between two DNA sequences, by extending Kimura’s (1980) two-parameter method to the case where a G+C-content bias exists. This method will be useful when there are strong transition-transversion and G+C-content biases, as in the case of Drosophila mitochondrial DNA. (Tamura 1992)

One frequency only ${\pi }_{GC}$

${\pi }_G={\pi }_C=\frac{{\pi }_{GC}}{2}$

${\pi }_A={\pi }_T=\frac{(1-{\pi }_{GC})}{2}$

Rate matrix $Q=\left(\begin{array}{cccc}*& \kappa (1-{\pi }_{GC})/2& (1-{\pi }_{GC})/2& (1-{\pi }_{GC})/2\\ \kappa {\pi }_{GC}/2& *& {\pi }_{GC}/2& {\pi }_{GC}/2\\ (1-{\pi }_{GC})/2& (1-{\pi }_{GC})/2& *& \kappa (1-{\pi }_{GC})/2\\ {\pi }_{GC}/2& {\pi }_{GC}/2& \kappa {\pi }_{GC}/2& *\end{array}\right)$

The evolutionary distance between two noncoding sequences according to this model is given by

$d=-h\ln (1-\frac{p}{h}-q)-\frac{1}{2}(1-h)\ln (1-2q)$

where $h=2\theta (1-\theta )$ where $\theta \in (0,1)$ is the GC content.

TN93 model (Tamura and Nei 1993)

The TN93 model distinguishes between the two different types of transition - i.e. (A <-> G) is allowed to have a different rate to (C<->T). Transversions are all assumed to occur at the same rate, but that rate is allowed to be different from both of the rates for transitions.

TN93 also allows unequal base frequencies (${\pi }_T\ne {\pi }_C\ne {\pi }_A\ne {\pi }_G$).

Rate matrix $Q=\left(\begin{array}{cccc}*& {\kappa }_1{\pi }_C& {\pi }_A& {\pi }_G\\ {\kappa }_1{\pi }_T& *& {\pi }_A& {\pi }_G\\ {\pi }_T& {\pi }_C& *& {\kappa }_2{\pi }_G\\ {\pi }_T& {\pi }_C& {\kappa }_2{\pi }_A& *\end{array}\right)$

GTR: Generalised time-reversible (Tavaré 1986)

GTR is the most general neutral, independent, finite-sites, time-reversible model possible. It was first described in a general form by Simon Tavaré in 1986.

The GTR parameters consist of an equilibrium base frequency vector, $\mathrm{\Pi }=({\pi }_1,{\pi }_2,{\pi }_3,{\pi }_4)$, giving the frequency at which each base occurs at each site, and the rate matrix

$Q=\left(\begin{array}{cccc}-(x_1+x_2+x_3)& \frac{{\pi }_1{x}_1}{{\pi }_2}& \frac{{\pi }_1{x}_2}{{\pi }_3}& \frac{{\pi }_1{x}_3}{{\pi }_4}\\ x_1& -(\frac{{\pi }_1{x}_1}{{\pi }_2}+x_4+x_5)& \frac{{\pi }_2{x}_4}{{\pi }_3}& \frac{{\pi }_2{x}_5}{{\pi }_4}\\ x_2& x_4& -(\frac{{\pi }_1{x}_2}{{\pi }_3}+\frac{{\pi }_2{x}_4}{{\pi }_3}+x_6)& \frac{{\pi }_3{x}_6}{{\pi }_4}\\ x_3& x_5& x_6& -(\frac{{\pi }_1{x}_3}{{\pi }_4}+\frac{{\pi }_2{x}_5}{{\pi }_4}+\frac{{\pi }_3{x}_6}{{\pi }_4})\end{array}\right)$

Therefore, GTR (for four characters, as is often the case in phylogenetics) requires 6 substitution rate parameters, as well as 4 equilibrium base frequency parameters. However, this is usually eliminated down to 9 parameters plus $\mu$, the overall number of substitutions per unit time. When measuring time in substitutions ($\mu$=1) only 8 free parameters remain.

In general, to compute the number of parameters, one must count the number of entries above the diagonal in the matrix, i.e. for n trait values per site $\frac{n^2-n}{2}$, and then add n for the equilibrium base frequencies, and subtract 1 because $\mu$ is fixed. One gets

$\frac{n^2-n}{2}+n-1=\frac{1}{2}{n}^2+\frac{1}{2}n-1.$

For example, for an amino acid sequence (there are 20 "standard" amino acids that make up proteins), one would find there are 209 parameters. However, when studying coding regions of the genome, it is more common to work with a codon substitution model (a codon is three bases and codes for one amino acid in a protein). There are $4^3=64$ codons, but the rates for transitions between codons which differ by more than one base is assumed to be zero. Hence, there are $\frac{20\times 19\times 3}{2}+64-1=633$ parameters.

See also

• Molecular evolution
• Molecular clock
• UPGMA

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.

Further reading

• One of the biggest reasons investing in a Singapore new launch is an effective things is as a result of it is doable to be lent massive quantities of money at very low interest rates that you should utilize to purchase it. Then, if property values continue to go up, then you'll get a really high return on funding (ROI). Simply make sure you purchase one of the higher properties, reminiscent of the ones at Fernvale the Riverbank or any Singapore landed property Get Earnings by means of Renting
  
  In its statement, the singapore property listing - website link, government claimed that the majority citizens buying their first residence won't be hurt by the new measures. Some concessions can even be prolonged to chose teams of consumers, similar to married couples with a minimum of one Singaporean partner who are purchasing their second property so long as they intend to promote their first residential property. Lower the LTV limit on housing loans granted by monetary establishments regulated by MAS from 70% to 60% for property purchasers who are individuals with a number of outstanding housing loans on the time of the brand new housing purchase. Singapore Property Measures - 30 August 2010 The most popular seek for the number of bedrooms in Singapore is 4, followed by 2 and three. Lush Acres EC @ Sengkang
  
  Discover out more about real estate funding in the area, together with info on international funding incentives and property possession. Many Singaporeans have been investing in property across the causeway in recent years, attracted by comparatively low prices. However, those who need to exit their investments quickly are likely to face significant challenges when trying to sell their property – and could finally be stuck with a property they can't sell. Career improvement programmes, in-house valuation, auctions and administrative help, venture advertising and marketing, skilled talks and traisning are continuously planned for the sales associates to help them obtain better outcomes for his or her shoppers while at Knight Frank Singapore. No change Present Rules
  
  Extending the tax exemption would help. The exemption, which may be as a lot as $2 million per family, covers individuals who negotiate a principal reduction on their existing mortgage, sell their house short (i.e., for lower than the excellent loans), or take part in a foreclosure course of. An extension of theexemption would seem like a common-sense means to assist stabilize the housing market, but the political turmoil around the fiscal-cliff negotiations means widespread sense could not win out. Home Minority Chief Nancy Pelosi (D-Calif.) believes that the mortgage relief provision will be on the table during the grand-cut price talks, in response to communications director Nadeam Elshami. Buying or promoting of blue mild bulbs is unlawful.
  
  A vendor's stamp duty has been launched on industrial property for the primary time, at rates ranging from 5 per cent to 15 per cent. The Authorities might be trying to reassure the market that they aren't in opposition to foreigners and PRs investing in Singapore's property market. They imposed these measures because of extenuating components available in the market." The sale of new dual-key EC models will even be restricted to multi-generational households only. The models have two separate entrances, permitting grandparents, for example, to dwell separately. The vendor's stamp obligation takes effect right this moment and applies to industrial property and plots which might be offered inside three years of the date of buy. JLL named Best Performing Property Brand for second year running
  
  The data offered is for normal info purposes only and isn't supposed to be personalised investment or monetary advice. Motley Fool Singapore contributor Stanley Lim would not personal shares in any corporations talked about. Singapore private home costs increased by 1.eight% within the fourth quarter of 2012, up from 0.6% within the earlier quarter. Resale prices of government-built HDB residences which are usually bought by Singaporeans, elevated by 2.5%, quarter on quarter, the quickest acquire in five quarters. And industrial property, prices are actually double the levels of three years ago. No withholding tax in the event you sell your property. All your local information regarding vital HDB policies, condominium launches, land growth, commercial property and more
  
  There are various methods to go about discovering the precise property. Some local newspapers (together with the Straits Instances ) have categorised property sections and many local property brokers have websites. Now there are some specifics to consider when buying a 'new launch' rental. Intended use of the unit Every sale begins with 10 p.c low cost for finish of season sale; changes to 20 % discount storewide; follows by additional reduction of fiftyand ends with last discount of 70 % or extra. Typically there is even a warehouse sale or transferring out sale with huge mark-down of costs for stock clearance. Deborah Regulation from Expat Realtor shares her property market update, plus prime rental residences and houses at the moment available to lease Esparina EC @ Sengkang
• One of the biggest reasons investing in a Singapore new launch is an effective things is as a result of it is doable to be lent massive quantities of money at very low interest rates that you should utilize to purchase it. Then, if property values continue to go up, then you'll get a really high return on funding (ROI). Simply make sure you purchase one of the higher properties, reminiscent of the ones at Fernvale the Riverbank or any Singapore landed property Get Earnings by means of Renting
  
  In its statement, the singapore property listing - website link, government claimed that the majority citizens buying their first residence won't be hurt by the new measures. Some concessions can even be prolonged to chose teams of consumers, similar to married couples with a minimum of one Singaporean partner who are purchasing their second property so long as they intend to promote their first residential property. Lower the LTV limit on housing loans granted by monetary establishments regulated by MAS from 70% to 60% for property purchasers who are individuals with a number of outstanding housing loans on the time of the brand new housing purchase. Singapore Property Measures - 30 August 2010 The most popular seek for the number of bedrooms in Singapore is 4, followed by 2 and three. Lush Acres EC @ Sengkang
  
  Discover out more about real estate funding in the area, together with info on international funding incentives and property possession. Many Singaporeans have been investing in property across the causeway in recent years, attracted by comparatively low prices. However, those who need to exit their investments quickly are likely to face significant challenges when trying to sell their property – and could finally be stuck with a property they can't sell. Career improvement programmes, in-house valuation, auctions and administrative help, venture advertising and marketing, skilled talks and traisning are continuously planned for the sales associates to help them obtain better outcomes for his or her shoppers while at Knight Frank Singapore. No change Present Rules
  
  Extending the tax exemption would help. The exemption, which may be as a lot as $2 million per family, covers individuals who negotiate a principal reduction on their existing mortgage, sell their house short (i.e., for lower than the excellent loans), or take part in a foreclosure course of. An extension of theexemption would seem like a common-sense means to assist stabilize the housing market, but the political turmoil around the fiscal-cliff negotiations means widespread sense could not win out. Home Minority Chief Nancy Pelosi (D-Calif.) believes that the mortgage relief provision will be on the table during the grand-cut price talks, in response to communications director Nadeam Elshami. Buying or promoting of blue mild bulbs is unlawful.
  
  A vendor's stamp duty has been launched on industrial property for the primary time, at rates ranging from 5 per cent to 15 per cent. The Authorities might be trying to reassure the market that they aren't in opposition to foreigners and PRs investing in Singapore's property market. They imposed these measures because of extenuating components available in the market." The sale of new dual-key EC models will even be restricted to multi-generational households only. The models have two separate entrances, permitting grandparents, for example, to dwell separately. The vendor's stamp obligation takes effect right this moment and applies to industrial property and plots which might be offered inside three years of the date of buy. JLL named Best Performing Property Brand for second year running
  
  The data offered is for normal info purposes only and isn't supposed to be personalised investment or monetary advice. Motley Fool Singapore contributor Stanley Lim would not personal shares in any corporations talked about. Singapore private home costs increased by 1.eight% within the fourth quarter of 2012, up from 0.6% within the earlier quarter. Resale prices of government-built HDB residences which are usually bought by Singaporeans, elevated by 2.5%, quarter on quarter, the quickest acquire in five quarters. And industrial property, prices are actually double the levels of three years ago. No withholding tax in the event you sell your property. All your local information regarding vital HDB policies, condominium launches, land growth, commercial property and more
  
  There are various methods to go about discovering the precise property. Some local newspapers (together with the Straits Instances ) have categorised property sections and many local property brokers have websites. Now there are some specifics to consider when buying a 'new launch' rental. Intended use of the unit Every sale begins with 10 p.c low cost for finish of season sale; changes to 20 % discount storewide; follows by additional reduction of fiftyand ends with last discount of 70 % or extra. Typically there is even a warehouse sale or transferring out sale with huge mark-down of costs for stock clearance. Deborah Regulation from Expat Realtor shares her property market update, plus prime rental residences and houses at the moment available to lease Esparina EC @ Sengkang

External links

• DAWG: DNA Assembly With Gaps — free software for simulating sequence evolution

Template:MolecularEvolution I do not commonly speak about myself. I get pleasure from sharing enjoyable written content with others. I also choose to investigate the world wide web since I think you can find much to determine to choose from. When I am on the world-wide-web.
|
My name: Gladis Bickford
Age: 38 years old
Country: United States
Town: Aurora
ZIP: 80011
Address: 4182 Mckinley Avenue
|
I am 36 years old and my name is Gladis Bickford. I life in Aurora (United States).
|
I’m Gladis from Aurora doing my final year engineering in History. I did my schooling, secured 89% and hope to find someone with same interests in Fencing.
|
Im addicted to my hobby Fencing.
I try to learn Norwegian in my free time.
|
My name is Gladis (29 years old) and my hobbies are Shooting sport and Conlanging.
|
Hello! I am Gladis. I am satisfied that I can unite to the entire world. I live in United States, in the CO region. I dream to visit the different nations, to get acquainted with appealing people.
|
I'm Gladis (28) from Aurora, United States.
I'm learning Norwegian literature at a local high school and I'm just about to graduate.
I have a part time job in a backery.
|
Hi!
My name is Gladis and I'm a 21 years old girl from United States.
|
Hello! My name is Gladis.
It is a little about myself: I live in United States, my city of Aurora.
It's called often Eastern or cultural capital of CO. I've married 2 years ago.
I have two children - a son (Melva) and the daughter (Armand). We all like Fencing.
|
Hello!
I'm Norwegian male :).
I really like Sons of Anarchy!
|
Hello from United States. I'm glad to came here. My first name is Gladis.
I live in a city called Aurora in east United States.
I was also born in Aurora 36 years ago. Married in December year 2011. I'm working at the the office.
|
My name's Gladis Bickford but everybody calls me Gladis. I'm from United States. I'm studying at the university (2nd year) and I play the Trumpet for 6 years. Usually I choose songs from the famous films ;).
I have two brothers. I love Gaming, watching TV (Sons of Anarchy) and Origami.
|
I'm a 47 years old and study at the high school (History).
In my free time I'm trying to learn Norwegian. I have been twicethere and look forward to go there sometime in the future. I love to read, preferably on my ebook reader. I really love to watch American Dad and Grey's Anatomy as well as documentaries about nature. I love Fencing.
|
I'm Gladis and I live in Aurora.
I'm interested in History, Fencing and Norwegian art. I like travelling and reading fantasy.
|
I'm Gladis and I live with my husband and our 3 children in Aurora, in the CO south part. My hobbies are Lapidary, Vehicle restoration and Amateur astronomy.
|
Hello, I'm Gladis, a 22 year old from Aurora, United States.
My hobbies include (but are not limited to) Bowling, Bboying and watching Sons of Anarchy.
|
Hi there! :) My name is Gladis, I'm a student studying History from Aurora, United States.
|
I'm Gladis and I live in a seaside city in northern United States, Aurora. I'm 26 and I'm will soon finish my study at History.
|
I am Gladis from Aurora. I am learning to play the Trumpet. Other hobbies are Fencing.
|
I am Gladis and was born on 6 December 1990. My hobbies are Meteorology and Sewing.
|
My name is Gladis and I am studying Neuroscience and Business at Aurora / United States.
}

Feel free to surf to my blog: computer store elk grove village