the kindlins represent a class of focal adhesion proteins implicated in integrin activation . they comprise three evolutionarily conserved members , kindlin-1 ( fermt,1 , c20orf42 , urp1 ; chromosome 20p12.3 ) , kindlin-2 ( fermt2 , mig2 , plekhc1 , unc112 , ; chromosome 14q22.1 ) and kindlin-3 ( fermt3 , unc-112 related protein 2 ; chromosome 11q13.1 ) , that share considerable sequence and structural similarities.1 the kindlins have a bipartite ferm ( four point one protein , ezrin , radixin , moesin ) domain interrupted by a pleckstrin homology domain and can bind directly to various classes of integrins as well as participating in inside  out integrin activation.2 the f3 sub domain of the ferm domains of all three kindlins contains a crucial phosphotyrosine binding ( ptb ) fold resembling that of talin.3,4 the relationship between kindlins and integrin signaling was explored by initial studies on c.elegans embryos harboring homozygous mutations in unc-112 , the worm ortholog of kind1 .
these embryos develop a paralyzed , arrested elongation at two - fold ( pat ) phenotype due to failure of organization of pat3/integrin in body muscle wall.5 functionally , kindlins primarily mediate protein 
protein interactions and lack catalytic domains.6 in recent years they have emerged as a key class of adaptor molecules involved in integrin activation.7 loss - of - function mutations in kindlin-1 and kindlin-3 cause kindler syndrome and leukocyte adhesion deficiency - iii syndrome , respectively . kindler syndrome was in fact the first human genetic disorder clinically associated with kindlins .
it is caused by mutation of kindlin-1 and is characterized by skin blistering , severe periodontitis and poililodermia.8,9 no human disease has yet been associated with kindlin-2 gene pathology , however kindlin-2 knockout mice die in early embryonic stage indicating the essential role it plays in development .
for instance , kindlin-1 is predominantly expressed in the epidermis and only weakly expressed in the dermis , while kindlin-3 expression is restricted exclusively to hematopoietic tissues , where it is the dominant form of kindlins expressed . on the other hand ,
kindlin-2 is ubiquitously expressed in most parts of the body.10 these differential expression patterns may in part explain the distinctive phenotypes that result from the loss of different kindlins . for instance , as noted above , the kindlin-2 homolog in c. elegans , unc-112 is essential for embryonic development.11 also , loss of kindlin-2 in mice results in pre- implantation embryonic lethality and knockdown of kindlin-2 in zebrafish reveals a strong relationship between cardiac development and the function of kindlin-2.12 consistent with the restricted expression of kindlin-3 in hematopoietic tissues , mice lacking kindlin-3 show severe osteoporosis , hemorrhage and defects in the erythrocyte membrane skeleton , and die within one week after birth.13,14 kindlin-1 shares 62% sequence similarity with kindlin-2 and 49% with kindlin-3 . until now
, the evolutionary aspects of kindlin structure and function have not been addressed . in this study
we studied the natural forces shaping the evolution of the kindlin family of proteins by comparing different evolutionary trends in different vertebrate clades .
a phylogenetic analysis of three kindlin family members illustrates the phylogenetic history of kindlin paralogs and also documents the duplication events leading to the formation of these paralogs from one single ancestral kindlin .
we show that the original kindlin arose at least as early as simple metazoans , such as hydra .
we also explored the effect of functional constraints on the evolution of these three paralogs in vertebrates and found evidence that purifying selection is a major force shaping the evolution of kindlins .
relative levels of kindlin-1 , kindlin-2 and kindlin-3 transcripts were determined by real - time rt - pcr using sybr green .
triplicate samples of each pcr mixture , each containing 4.7 l of power sybr green pcr master mixture ( applied biosystems ) , 0.3 l of a 10 pmol/l of primer mixture , 0.3 l of cdna , and water to a total volume of 10 l were transferred into a 96-well plate on an abi 7500 fast real time pcr system ( applied biosystems ) .
the samples were initially incubated at 95c for 3 min , followed by 45 cycles with 95c for 15 s , 60c for 60 s. dissociation curves were generated after each pcr run to ensure that a single , specific product was amplified .
the results were analyzed with the comparative cycle threshold ( ct ) method . for normalization , we used the expression level of -actin ( actb ) .
the pcr primers are shown in table 1 . in order to explore the evolutionary history of kindlins , sequences of the complete transcripts and the corresponding protein sequences of all three kindlins ( kindlin-1 , kindlin-2 and kindlin-3 ) from different species
were extracted from ncbi ( http://www.ncbi.nlm.nih.gov ) and ensemble ( http://www.ensembl.org ) genome browsers ( table 2 ) .
after alignment using clustalw program,15 all positions containing gaps and missing data were eliminated , because of the possible ambiguity of the alignments .
this stringent approach reduced the risk of misinterpretations the evolutionary history was inferred by using the maximum likelyhood method ( ml ) as implemented in the treefinder ( tf ) program package.16 the tf support values indicate the reliability of internal branches .
the analyses of amino acid sequences were performed using the wag2000 model17 applying eight classes of rate heterogeneity among sites ( 8 ) .
a neighbor joining analysis18 on amino acid sequences was done by the mega 4 program.19 the reliability of the nj tree was estimated by the bootstrap method , based on 1000 pseudo replicates .
a variety of methods were employed to explore the functional constraints shaping the evolution of kindlins .
pairwise comparison of the number of synonymous nucleotide substitutions per synonymous ( ds ) site and non - synonymous nucleotide substitutions per non - synonymous site ( dn ) was carried out by using nei - gojobori method.20 in addition to pairwise methods , the dn / ds ratio in different branches of the maximum - likelihood tree was estimated using the codon - based genetic algorithm implemented in the ga - branch program available at the data - monkey server ( http://www.datamonkey.org/help/gabranch.php ) .
this approach assigns each branch to an incrementally estimated class of dn / ds ratios without requiring a specification of the branches a priori and is both less parameterized and more user friendly than a fully local model.21 evolutionary distance between all possible pairs of kindlin paralogs was estimated by tajima s relative rate test.22 each pair of paralogs was compared with amphioxus protein sequences taken as an out - group .
realtime pcr analyses of the expression pattern of kindlins revealed distinct patterns of kindlin expression . as expected , kindlin-2 was expressed almost ubiquitously in all six of the tissues tested , while kindlin-1 showed a one hundred - fold lower expression in each of these tissues with the exception of human kidneys where the expression level is low but significant .
1 ) . these results are very much in agreement with existing data on the expression patterns of kindlin proteins , which also show the ubiquitous nature of kindlin-2 expression and the tissue specific expression of both kindlin-1 and kindlin-3.10 phylogenetic analysis was carried out on the amino acid sequences of kindlin proteins , with the phylogenetic tree rooted by orthologous genes from invertebrate species ( fig .
notably , an nj analysis on the same data predicted the same topology ( fig .
the resulting phylogenetic tree is very well supported for most branches except for the deep divergences of branchisotoma floridae and strongylocentrotus purpuratus and one leading to the grouping of xenopus laevis kindlin-2 with the mammals .
otherwise the vertebrate relationships for each kindlin ortholog are resolved in a manner concomitant with previous phylogenomic analyses,23 indicating a reliable evolutionary reconstruction among the kindlin families .
the phylogeny exhibits a topology of the form ( a)(bc ) ie , kindlin-1 and kindlin-3 form a cluster while kindlin-2 forms an out - group .
it appears that in invertebrates before the divergences of arthropods , a single ancestral kindlin had its ancient origin in hydra .
this gene than underwent a lineage specific duplication event in insects , giving rise to two kindlin paralogs .
although no study exists to date on the roles of kindlins in insects , it is plausible to assume that because of the diversity found in insects , the duplicated copies were maintained in response to selection pressures impacting this highly diverse group .
one of the two kindlins was lost in other higher phyla before the origin of vertebrates ie , echinoderms ( strongylocentrotus purpuratus ) urochordates ( ciona intestinalis ) , cephalochordates ( amphioxus ) and hemidchordates ( saccoglossus kowalevskii ) . however , the remaining single kindlin gene copy underwent two duplication events in vertebrates that may have occurred together with two rounds of genome duplication .
regardless of the exact mechanism , these events gave rise to three kindlin paralogs : kindlin-1 , kindlin-2 and kindlin-3 .
more data is still required to determine whether these kindlin gene duplication events in the fish were due to whole genome duplications or individual segmental gene duplication .
whatever was the cause of this duplication , it is evident from our analyses that three vertebrate kindlin paralogs originated after duplication events in the fish genome and were maintained in all subsequent vertebrate forms , probably due to the selection pressures that have promoted the diverse and complicated morphological and physiological properties of vertebrates.2427 it is important to notice that branch lengths of kindlin-2 are quite short when compared to both kindlin-1 and kindlin-3 .
very short branch lengths for kindlin-2 on the phylogenetic tree indicate that this gene has experienced a slower evolutionary rate than its paralogs .
this result is congruent with existing experimental studies on kindlin genes.9,10,13,14 for instance , kindlin-2 is a ubiquitously expressed gene playing its structural and functional roles in broad array of tissues.10 it has been indicated in various studies that ubiquitously expressed genes tend to evolve slowly compared to those with tissue specific expression.28 as noted previously , kindlin-1 and kindlin-3 are expressed in specific tissues ( epithelial tissues and the hematopoietic system , respectively ) ( fig .
1 ) , and it is therefore commensurate that the evolutionary rate of these paralogs is much higher than that of kindlin- 2 . in support of this
are kindlin knock out studies for all three kindlins which show that kindlin-2 knockout mice die during early embryogenesis.29 in contrast to the milder phenotypes of kindlin-1 and kindlin-3 , these studies support the idea that kindlin-2 is under tighter functional constraint than kindlin 1 and 3 .
comparison between non - synonymous and synonymous substitutions in orthologous transcript sequences can reveal the selective pressure that shapes the evolution of these genes .
the ratio between non - synonymous substitution per non - synonymous site to synonymous substitution per synonymous site dn / ds or omega (  ) indicates whether the evolution of genes is due to adaptive selection or due to neutral evolution .
a value of  more than 1 suggests that the gene is under positive selection . a value close to 1 suggests that a gene is under neutral selection and is experiencing neutral evolution .
however , a value of less than 1 indicates that a gene is under the influence of negative or purifying selection .
the nei - gojobori method we employed to estimate the dn / ds ratio clearly showed  values below one , ( ie , dn / ds < 1 ) for all the three kindlin paralogs in all of the vertebrate species compared for this analysis ( table 3 ) .
interestingly , the  value for kindlin-2 was much lower ( at least ten fold less ) than  values for either of the other two kindlins .
similarly the analysis within and between mammalian and non - mammalian groups clearly showed the same pattern observed in individual pairwise comparisons between different species ie , the absence of any positive selection whatsoever in all the groups .
. this difference may result from the high substitution rate often seen in fish genomes , leading to higher dn / ds ratios .
kindlin-2 in both mammals and non - mammals showed values of ds similar to kindlin-1 and kindlin-3 , but with much lower values of dn , thereby producing greatly lowered values of  ( table 4 ) . to gain further insight into the lineage specific nature of the selective pressures acting on each branch of the phylogenetic tree
this method , unlike the branch site method , does not require the manual selection of branches of interest to identify evidence for positive or negative selection . because the ga branch method does not require the user to select branches of interest , or that testing be performed one branch at a time , it experiences reduced statistical instability while also offering improved interpretability for poorly supported models .
in addition , inferences based on multiple models ( as opposed to a null - alternative pair ) are less vulnerable to model misspecification . in our study ,
ga - branch analysis selected a model with five classes of . in total , 72% of branches are assigned to a  of 0.037 , named as class d , with the remaining 28% of branches assigned to four additional classes , designated a , b , c and e , with  values of 0.102 , 0.050 , 0.037 and 0.026 , respectively .
none of the branches studied show any trend for positive selection , with the probability of positive selection being 0% for each branch of the tree .
notably , here again the very low value of  ( 0.013 ) for kindlin-2 indicates that it is evolving under the influence of much stronger negative selection than kindlin-1 and kindlin-3 ( fig . 4 and table 5 ) .
in short , whether calculated by pairwise comparison or by lineage specific analysis , dn / ds ratios consistently indicate that vertebrate kindlins have been evolving under the influence of purifying selection , with kindlin-2 under much stronger negative selection than kindlin-1 and kindlin-3 .
for instance , if adaptive selection was the main force for kindlin divergence rather than purifying selection , we would expect that the functional roles of kindlins may be diverse .
in fact , the hallmark function of kindlins is integrin activation and all the other higer order functions associated with kindlins , including cell migration , cell spreading , cell adhesion , cellular signaling and cancer promotion are associated with the ability of kindlins to activate integrin . on the other hand , these higher order processes impacted by kindlins  through integrin activation 
are so essential for organism survival and viability that very strong functional constraints exist over kindlin evolution . the more stringent functional constraint of kindlin-2 compared to either of its two counterpart paralogs is likely due to the fact that it is expressed ubiquitously in the body while kindlins 1 and 3 are tissue specific . for studying the evolutionary distances between different kindlin paralogs ,
this test involves pairwise comparison of protein sequences from kindlin paralogs of each species while using the orthologous kindlin sequences from amphioxus ( branchiostoma floridae ) as an out - group ( table 6 ) .
these tests produced intriguing results which , in correlation with phylogenetic data , show that kindlin-2 has undergone relatively limited divergence compared to both kindlin-1 and kindlin-3 , of which kindlin-3 is indicated to be most divergent .
interestingly , it seems that kindlin-2 is the representative of the original ancient amphioxus kindlin which underwent two duplications in fish giving rise to kindlin-1 and kindlin-3 .
there are two important reasons to believe that kindlin-2 is the representative of the unduplicated ancestral kindlin gene , whose duplication in vertebrates gave rise to two other paralogs .
firstly the tajima test clearly shows that kindlin-2 is closest to amphioxus kindlin , with very low levels of divergence evident compared to kindlin-1 and kindlin-3 . secondly , kindlin-2 is not only a ubiquitously expressed protein but also the only kindlin protein expressed in embryonic stem cells .
although no study has been conducted to explore the expression pattern of kindlins in amphioxus , c. elegans studies suggest that this unduplicated kindlin is expressed ubiquitously in all tissues including the embryonic stem cell.5 therefore if the amphioxus kindlin is an ortholog of c. elegans kindlin , it should also exhibit the same expression pattern , making vertebrate kindlin-2 a strong candidate as the representative of the ancestral unduplicated kindlin protein .
no species specific asymmetry was found in the divergence pattern of any of the three kindlin paralogs which is in agreement with the generally accepted principle that paralogs evolve at similar rates in different species.30 relatively high divergence of kindlin-1 and kindlin-3 relative to kindlin-2 can be explained by analyzing the expression pattern of the proteins.10 during evolution , from amphioxus to higher vertebrates , the expression pattern of the two paralogs have diverged .
our analysis of kindlin expression patterns in human tissue samples clearly shows that unlike kindlin-2 , which is expressed ubiquitously , both kindlin-1 and kindlin-3 present tissue specific expression patterns ( fig .
similarly , ussar et al have shown in their study that both kindlin-1 and kindlin-3 are expressed predominantly in epithelial tissues and the hematopoietic system respectively.10 the distinct expression patterns of all three kindlins in vertebrates thus supports the subfunctionalization model of gene duplication.31 it seems likely that that the function of the ancestral unduplicated kindlin was subfunctionalized in vertebrates in part due to the divergence of kindlin expression location .
thus , while kindlin-1 and kindlin-3 are expressed exclusively in epithelial and hematopoietic tissues , respectively , kindlin-2  being the representative of ancestral kindlin gene  is expressed in a variety of tissues , in a pattern less ubiquitous than the original unduplicated kindlin gene .
ultimately , it seems that this subfunctionalization of kindlin expression patterns may have provided a degree of selective advantage associated with the diversification of higher order functions performed by integrins in various tissues .
realtime pcr analyses of the expression pattern of kindlins revealed distinct patterns of kindlin expression . as expected , kindlin-2 was expressed almost ubiquitously in all six of the tissues tested , while kindlin-1 showed a one hundred - fold lower expression in each of these tissues with the exception of human kidneys where the expression level is low but significant .
1 ) . these results are very much in agreement with existing data on the expression patterns of kindlin proteins , which also show the ubiquitous nature of kindlin-2 expression and the tissue specific expression of both kindlin-1 and kindlin-3.10
phylogenetic analysis was carried out on the amino acid sequences of kindlin proteins , with the phylogenetic tree rooted by orthologous genes from invertebrate species ( fig .
notably , an nj analysis on the same data predicted the same topology ( fig .
3 ) . the resulting phylogenetic tree is very well supported for most branches except for the deep divergences of branchisotoma floridae and strongylocentrotus purpuratus and one leading to the grouping of xenopus laevis kindlin-2 with the mammals .
otherwise the vertebrate relationships for each kindlin ortholog are resolved in a manner concomitant with previous phylogenomic analyses,23 indicating a reliable evolutionary reconstruction among the kindlin families .
the phylogeny exhibits a topology of the form ( a)(bc ) ie , kindlin-1 and kindlin-3 form a cluster while kindlin-2 forms an out - group .
it appears that in invertebrates before the divergences of arthropods , a single ancestral kindlin had its ancient origin in hydra .
this gene than underwent a lineage specific duplication event in insects , giving rise to two kindlin paralogs .
although no study exists to date on the roles of kindlins in insects , it is plausible to assume that because of the diversity found in insects , the duplicated copies were maintained in response to selection pressures impacting this highly diverse group .
one of the two kindlins was lost in other higher phyla before the origin of vertebrates ie , echinoderms ( strongylocentrotus purpuratus ) urochordates ( ciona intestinalis ) , cephalochordates ( amphioxus ) and hemidchordates ( saccoglossus kowalevskii ) . however , the remaining single kindlin gene copy underwent two duplication events in vertebrates that may have occurred together with two rounds of genome duplication .
regardless of the exact mechanism , these events gave rise to three kindlin paralogs : kindlin-1 , kindlin-2 and kindlin-3 .
more data is still required to determine whether these kindlin gene duplication events in the fish were due to whole genome duplications or individual segmental gene duplication .
whatever was the cause of this duplication , it is evident from our analyses that three vertebrate kindlin paralogs originated after duplication events in the fish genome and were maintained in all subsequent vertebrate forms , probably due to the selection pressures that have promoted the diverse and complicated morphological and physiological properties of vertebrates.2427 it is important to notice that branch lengths of kindlin-2 are quite short when compared to both kindlin-1 and kindlin-3 .
very short branch lengths for kindlin-2 on the phylogenetic tree indicate that this gene has experienced a slower evolutionary rate than its paralogs .
this result is congruent with existing experimental studies on kindlin genes.9,10,13,14 for instance , kindlin-2 is a ubiquitously expressed gene playing its structural and functional roles in broad array of tissues.10 it has been indicated in various studies that ubiquitously expressed genes tend to evolve slowly compared to those with tissue specific expression.28 as noted previously , kindlin-1 and kindlin-3 are expressed in specific tissues ( epithelial tissues and the hematopoietic system , respectively ) ( fig .
1 ) , and it is therefore commensurate that the evolutionary rate of these paralogs is much higher than that of kindlin- 2 . in support of this
are kindlin knock out studies for all three kindlins which show that kindlin-2 knockout mice die during early embryogenesis.29 in contrast to the milder phenotypes of kindlin-1 and kindlin-3 , these studies support the idea that kindlin-2 is under tighter functional constraint than kindlin 1 and 3 .
comparison between non - synonymous and synonymous substitutions in orthologous transcript sequences can reveal the selective pressure that shapes the evolution of these genes .
the ratio between non - synonymous substitution per non - synonymous site to synonymous substitution per synonymous site dn / ds or omega (  ) indicates whether the evolution of genes is due to adaptive selection or due to neutral evolution .
a value of  more than 1 suggests that the gene is under positive selection . a value close to 1 suggests that a gene is under neutral selection and is experiencing neutral evolution .
however , a value of less than 1 indicates that a gene is under the influence of negative or purifying selection .
the nei - gojobori method we employed to estimate the dn / ds ratio clearly showed  values below one , ( ie , dn / ds < 1 ) for all the three kindlin paralogs in all of the vertebrate species compared for this analysis ( table 3 ) .
interestingly , the  value for kindlin-2 was much lower ( at least ten fold less ) than  values for either of the other two kindlins .
similarly the analysis within and between mammalian and non - mammalian groups clearly showed the same pattern observed in individual pairwise comparisons between different species ie , the absence of any positive selection whatsoever in all the groups .
this difference may result from the high substitution rate often seen in fish genomes , leading to higher dn / ds ratios .
kindlin-2 in both mammals and non - mammals showed values of ds similar to kindlin-1 and kindlin-3 , but with much lower values of dn , thereby producing greatly lowered values of  ( table 4 ) . to gain further insight into the lineage specific nature of the selective pressures acting on each branch of the phylogenetic tree
this method , unlike the branch site method , does not require the manual selection of branches of interest to identify evidence for positive or negative selection . because the ga branch method does not require the user to select branches of interest , or that testing be performed one branch at a time , it experiences reduced statistical instability while also offering improved interpretability for poorly supported models .
in addition , inferences based on multiple models ( as opposed to a null - alternative pair ) are less vulnerable to model misspecification . in our study ,
ga - branch analysis selected a model with five classes of . in total , 72% of branches are assigned to a  of 0.037 , named as class d , with the remaining 28% of branches assigned to four additional classes , designated a , b , c and e , with  values of 0.102 , 0.050 , 0.037 and 0.026 , respectively .
none of the branches studied show any trend for positive selection , with the probability of positive selection being 0% for each branch of the tree .
notably , here again the very low value of  ( 0.013 ) for kindlin-2 indicates that it is evolving under the influence of much stronger negative selection than kindlin-1 and kindlin-3 ( fig . 4 and table 5 ) .
in short , whether calculated by pairwise comparison or by lineage specific analysis , dn / ds ratios consistently indicate that vertebrate kindlins have been evolving under the influence of purifying selection , with kindlin-2 under much stronger negative selection than kindlin-1 and kindlin-3
for instance , if adaptive selection was the main force for kindlin divergence rather than purifying selection , we would expect that the functional roles of kindlins may be diverse .
in fact , the hallmark function of kindlins is integrin activation and all the other higer order functions associated with kindlins , including cell migration , cell spreading , cell adhesion , cellular signaling and cancer promotion are associated with the ability of kindlins to activate integrin . on the other hand , these higher order processes impacted by kindlins  through integrin activation 
are so essential for organism survival and viability that very strong functional constraints exist over kindlin evolution .
the more stringent functional constraint of kindlin-2 compared to either of its two counterpart paralogs is likely due to the fact that it is expressed ubiquitously in the body while kindlins 1 and 3 are tissue specific .
for studying the evolutionary distances between different kindlin paralogs , tajima s relative rate test was employed .
this test involves pairwise comparison of protein sequences from kindlin paralogs of each species while using the orthologous kindlin sequences from amphioxus ( branchiostoma floridae ) as an out - group ( table 6 ) .
these tests produced intriguing results which , in correlation with phylogenetic data , show that kindlin-2 has undergone relatively limited divergence compared to both kindlin-1 and kindlin-3 , of which kindlin-3 is indicated to be most divergent .
interestingly , it seems that kindlin-2 is the representative of the original ancient amphioxus kindlin which underwent two duplications in fish giving rise to kindlin-1 and kindlin-3 .
there are two important reasons to believe that kindlin-2 is the representative of the unduplicated ancestral kindlin gene , whose duplication in vertebrates gave rise to two other paralogs .
firstly the tajima test clearly shows that kindlin-2 is closest to amphioxus kindlin , with very low levels of divergence evident compared to kindlin-1 and kindlin-3 .
secondly , kindlin-2 is not only a ubiquitously expressed protein but also the only kindlin protein expressed in embryonic stem cells .
although no study has been conducted to explore the expression pattern of kindlins in amphioxus , c. elegans studies suggest that this unduplicated kindlin is expressed ubiquitously in all tissues including the embryonic stem cell.5 therefore if the amphioxus kindlin is an ortholog of c. elegans kindlin , it should also exhibit the same expression pattern , making vertebrate kindlin-2 a strong candidate as the representative of the ancestral unduplicated kindlin protein .
no species specific asymmetry was found in the divergence pattern of any of the three kindlin paralogs which is in agreement with the generally accepted principle that paralogs evolve at similar rates in different species.30 relatively high divergence of kindlin-1 and kindlin-3 relative to kindlin-2 can be explained by analyzing the expression pattern of the proteins.10 during evolution , from amphioxus to higher vertebrates , the expression pattern of the two paralogs have diverged .
our analysis of kindlin expression patterns in human tissue samples clearly shows that unlike kindlin-2 , which is expressed ubiquitously , both kindlin-1 and kindlin-3 present tissue specific expression patterns ( fig .
similarly , ussar et al have shown in their study that both kindlin-1 and kindlin-3 are expressed predominantly in epithelial tissues and the hematopoietic system respectively.10 the distinct expression patterns of all three kindlins in vertebrates thus supports the subfunctionalization model of gene duplication.31 it seems likely that that the function of the ancestral unduplicated kindlin was subfunctionalized in vertebrates in part due to the divergence of kindlin expression location .
thus , while kindlin-1 and kindlin-3 are expressed exclusively in epithelial and hematopoietic tissues , respectively , kindlin-2  being the representative of ancestral kindlin gene  is expressed in a variety of tissues , in a pattern less ubiquitous than the original unduplicated kindlin gene .
ultimately , it seems that this subfunctionalization of kindlin expression patterns may have provided a degree of selective advantage associated with the diversification of higher order functions performed by integrins in various tissues .
they play an important role in cell migration , cell spreading and cancer progression through a core molecular mechanism of integrin activation . in this study
we have shown that the ancestral kindlin gene was an unduplicated single gene found in organisms as primitive as hydra .
this gene then underwent a lineage specific duplication in insects , giving rise to two kindlin paralogs , while the three paralogs of kindlin found in vertebrates are the result of duplication events that occurred in fish .
of the three kindlins , kindlin-2 has undergone the least evolutionary divergence , probably due to stringent functional constraints it associated with its virtually ubiquitous expression pattern in body tissues and especially in embryonic stem cells . on the other hand
both kindlin-1 and kindlin-3 showed significantly greater divergence as a result of significantly weaker functional constraints , possibly resulting from their subfunctionalization into very specific portions of the body .
the comparison of synonymous to non - synonymous substitutions both by a pairwise method as well as a lineage specific method also indicate that all three kindlins have been evolving under strong negative selection .
 kindlin proteins represent a newly discovered family of evolutionarily conserved ferm domain - containing proteins . 
 this family includes three highly conserved proteins : kindlin-1 , kindlin-2 and kindlin-3 . 
 all three kindlin proteins are associated with focal adhesions and are involved in integrin activation . 
 the ferm domain of each kindlin is bipartite and plays a key role in integrin activation . 
 we herein explore for the first time the evolutionary history of these proteins . 
 the phylogeny of the kindlins suggests a single ancestral kindlin protein present in even the earliest metazoan ie , hydra . 
 this protein then underwent duplication events in insects and also experienced genome duplication in vertebrates , leading to the kindlin family . 
 a comparative study of the kindlin paralogs showed that kindlin-2 is the slowest evolving protein among the three family members . 
 the analysis of synonymous and non - synonymous substitutions in orthologous kindlin sequences in different species showed that all three kindlins have been evolving under the influence of purifying selection . 
 the expression pattern of kindlins along with phylogenetic studies supports the subfunctionalization model of gene duplication . 
Introduction
Material Method
Results and Discussion
Expression pattern of kindlins
Phylogeny
Estimation of the selective forces shaping the evolution of Kindlin proteins
Divergence rate between paralogs
Conclusion