Estudo de alterações na função salivar de um modelo de ratos pré-diabéticos

A Diabetes Mellitus é uma doença metabólica com uma distribuição mundial cada vez mais acentuada, acarretando consigo, elevada morbilidade e mortalidade. Deste modo, é de suma importância a compreensão e estudo em fases precoces desta doença. Nos últimos anos a saliva tem sido usada como fonte de biomarcadores para detecção e monitorização de doenças. Para além disso, há referencia a alterações na secreção salivar e percepção dos alimentos em indivíduos diabéticos, tornando-se importante um maior conhecimento da função salivar associada a esta patologia. O objectivo geral deste trabalho é estudar a função salivar num modelo de roedores intolerantes à glicose, avaliando-se a concentração proteica total da saliva, morfologia e morfometria das glândulas salivares maiores e expressão da amilase na saliva e glândulas salivares. Não se observam diferenças entre os grupos, para a concentração em proteína total, nem para a expressão de α-amilase, na saliva. Na análise histológica das glândulas salivares maiores é possível observar dimensões significativamente maiores dos ácinos das glândulas analisadas (parótida, submandibular e sublingual), em relação aos animais normoglicémicos, sugerindo o que pode ser o início de uma patologia comum em diabéticos, a sialose (hipertrofia e hiperplasia das glândulas salivares). Através de imunomarcação para a actina observa-se um aumento da expressão de células mioepiteliais na seguinte ordem: parótida, submandibular e sublingual. Não se observam diferenças entre os grupos. A imunomarcação para a α-amilase é mais intensa nos ácinos da parótida, nos ductos granulares da submandibular e nas meias luas serosas da sublingual, em ambos os grupos. Verifica-se uma sobre-expressão nos ductos intralobulares da parótida e nos ductos granulares da sumandibular, nos animais pré-diabéticos, sugerindo alguma endocitose da enzima por parte das células dos ductos. Os resultados obtidos sugerem alterações na função salivar, numa fase prévia ao desenvolvimento da Diabetes Mellitus, as quais merecem ser exploradas.

Filodinâmica do vírus da tristeza dos citrinos: epidemiologia molecular e biologia evolutiva

O Citrus tristeza virus (CTV), responsável por várias doenças em citrinos, é um dos maiores condicionantes da citricultura a nível mundial. Existem diversos isolados de CTV com diferentes características biológicas e moleculares, sendo que os sintomas causados pelo vírus dependem essencialmente do isolado viral e da combinação variedade/porta-enxerto. A implementação de medidas de controlo da doença depende, em grande parte, do tipo de isolados presentes numa dada região. No Capítulo 2, efetuou-se uma análise comparativa entre dois métodos de tipificação de isolados de CTV e verificou-se que a caracterização por PCR assimétrico-ELISA, que considera a existência de sete grupos, é mais adequada à descrição da estrutura genética de CTV. Estes resultados foram complementados com o estudo da dinâmica de colonização de cada grupo filogenético através de um imuno-ensaio in situ (Capítulo 3). Os resultados obtidos sugerem que os isolados de CTV diferem na quantidade de células infetadas e que essa diferença parece estar relacionada com a severidade do isolado. No Capítulo 4, o estudo da variabilidade genómica da região 3’ terminal permitiu verificar que a estrutura de grupos obtida para o gene da proteína da cápside (CP) é extensível a toda a região 3’ terminal que contém os genes mais fortemente implicados na interação com o hospedeiro. A estabilidade da estrutura genética nesta região foi também inferida a partir da pesquisa de eventos de recombinação. Os resultados sugerem uma baixa frequência de recombinação entre isolados de CTV, mesmo em isolados contendo mistura de haplótipos e mantidos há mais de 12 anos no mesmo hospedeiro. Adicionalmente, foi estimada a taxa de evolução de CTV através de um método estatístico Bayesiano (Capítulo 5). Para tal, foram usadas sequências do gene da CP de isolados de diversas regiões do mundo, pertencentes a diferentes grupos filogenéticos e obtidas entre 1990 e 2010. A taxa média de evolução estimada foi de 1,58 X 10-4 substituições nucleotídicas / ano. No geral, os resultados destes dois capítulos mostram que os isolados de CTV mantêm uma elevada estabilidade genética ao longo do tempo. Finalmente, no Capítulo 6, foi estudada a situação epidemiológica de CTV em Portugal continental a partir de isolados de CTV recolhidos no campo, onde se verificou que a maioria das árvores infetadas era composta por isolados de CTV pertencentes ao grupo M, ou seja isolados considerados suaves e que não provocam sintomas severos.

Insights into molecular and cellular events involved in normal and pathological vertebral development and fusion using zebrafish as model organism

The vertebral column and its units, the vertebrae, are fundamental features, characteristic of all vertebrates. Developmental segregation of the vertebral bodies as articulated units is an intrinsic requirement to guarantee the proper function of the spine. Whenever these units become fused either during development or postsegmentation, movement is affected in a more or less severe manner, depending on the number of vertebrae affected. Nevertheless, fusion may occur as part of regular development and as a physiological requirement, like in the tetrapod sacrum or in fish posterior vertebrae forming the urostyle. In order to meet the main objective of this PhD project, which aimed to better understand the molecular and cellular events underlying vertebral fusion under physiological and pathological conditions, a detailed characterization of the vertebral fusion occurring in zebrafish caudal fin region was conducted. This showed that fusion in the caudal fin region comprised 5 vertebral bodies, from which, only fusion between [PU1++U1] and ural2 [U2+] was still traceable during development. This involved bone deposition around the notochord sheath while fusion within the remaining vertebral bodies occur at the level of the notochord sheath, as during the early establishment of the vertebral bodies. A comparison approach between the caudal fin vertebrae and the remaining vertebral column showed conserved features such as the presence of mineralization related proteins as Osteocalcin were identified throughout the vertebral column, independently on the mineralization patterns. This unexpected presence of Osteocalcin in notochord sheath, here identified as Oc1, suggested that this gene, opposing to Oc2, generally associated with bone formation and mature osteoblast activity, is potentially associated with early mineralization events including chordacentrum formation. Nevertheless, major differences between caudal fin region and anterior vertebral bodies considering arch histology and mineralization patterns, led us to use RA as an inductive factor for vertebral fusion, allowing a direct comparison of equivalent structures under normal and fusion events. This fusion phenotype was associated with notochord sheath ectopic mineralization instead of ectopic perichordal bone formation related with increased osteoblast activity, as suggested in previous reports. Additionally, alterations in ECM content, cell adhesion and blood coagulation were discussed as potentially related with the fusion phenotype. Finally, Matrix gla protein, upregulated upon RA treatment and shown to be associated with chordacentrum mineralization sites in regular development, was further described considering its potential function in vertebral formation and pathological fusion. Therefore with this work we propose zebrafish caudal fin vertebral fusion as a potential model to study both congenital and postsegmentation fusion and we present candidate factors and genes that may be further explored in order to clarify whether we can prevent vertebral fusion.

Application of molecular tools for detection of plant viruses

Grapevine leafroll disease (GLRD) is one of the most important virus diseases of grapevines worldwide, causing major economical impact. The disease has a complex aetiology and currently eleven phloem-limited viruses, termed in general Grapevine leafroll-associated virus (GLRaVs), have been identified. Two of the GLRaVs, GLRaV-1 and GLRaV-3, are included in the European certification scheme of propagation material. However, the flawed notion that GLRaV-3 is more frequent than GLRaV-1 and that all other GLRaVs are possibly not as relevant for GLRD, has until now precluded the development of specific serological and molecular detection assays and limited the scope of molecular characterization of the viruses known to be associated with the disease. Hence, few studies have addressed the phylodynamics of GLRaVs or even characterized the genetic structure of their natural populations. This generalized lack of molecular information, in turn underlie the deficient capacity to detect the viruses. The phylogenetic analyses were conducted on the basis of the heat shock protein 70 homologue (HSP70h) and the coat protein (CP) genes for GLRaV-1 and the HSP70h, the heat shock protein 90 homologue (HSP90h) and the CP genes for GLRaV-5. The data obtained for GLRaV-1 contributed 83 new CP sequences. This information was combined with previous analysis by other authors and used for the production of new polyclonal IgG, capable of detecting CP variants from all the phylogroups observed. Successful testing of this new tool included tissue print immunoblotting (TPIB) and in situ immunoassay (ISIA). The data obtained for GLRaV-5, contributed 61 new CP and 28 new HSP90h gene sequences. Eight phylogenetic groups were identified on the basis of the CP. Characterization of the genetic structure of the isolates revealed a higher diversity than previously reported and allowed the identification of dominant virus variants. For both GLRaV-1 and GLRaV-5, the effect of vegetative propagation on the virus transmission dynamics was addressed.

Study of transcripts and reactive oxygen species involved in the defence responses of Quercus suber against Phytophthora cinnamomi

The present work has the merit of exploring an insight into the activation of defence genes of Quercus suber during response to infection by Phytophthora cinnamomi. Thus, cDNA-AFLP methodology was used to identify gene fragments differentially present in the mRNA profiles of host cells of micropropagated Q. suber plantlets roots infected with zoospores of P. cinnamomi at different post challenge time points. Six candidate genes were selected based on their interesting cDNA-AFLP expression patterns and homology to genes known to play a role in defence. These six genes encode a cinnamyl alcohol dehydrogenase 2 (QsCAD2), a protein disulphide isomerase (QsPDI), a CC-NBS-LRR resistance protein (QsRPc), thaumatin-like protein (QsTLP), chitinase (QsCHI) and a 1,3-beta glucanase (QsGLU). The current work has been successful in evaluation of the expression of these genes by qRT-PCR. Data analysis revealed that transcript levels of QsRPc, QsCHI, QsCAD2 and QsPDI increased during the early hours of inoculation, while transcript profiles of thaumatin-like protein showed decreasing. No expression was detected for 1,3-beta-glucanase (QsGLU). Furthermore, the choice of suitable reference genes in any new experimental system is absolutely crucial in qRT-PCR; for this reason in this study and for the first time a set of potential reference genes were analyzed and validated for qRT-PCR normalization in the patho-system Phytophthora-Q. suber. Four candidate reference genes polimerase II (QsRPII), eukaryotic translation initiation factor 5A(QsEIF-5A), b-tubulin (QsTUB) and a medium subunit family protein of Clathrin adaptor complexes (QsCACs) were evaluated to determine the most stable internal references in Q. suber. Analysis of stability of genes was carried out using Genex software. Results indicated all these four potential reference genes assumed stable expression. Data analysis revealed that QsRPII and QsCACs were the two most stable genes, while genes QsTUB and QsEIF-5A were the third and the fourth most stable gene, respectively. In this study, a plasmid-based quantitative PCR method was developed to measure P. cinnamomi colonization during infection process of Q. suber. Plasmid-based detection of P. cinnamomi showed a gradual accumulation of the pathogen DNA in cork oak root tips up to 24 h post infection. The higher increase in P. cinnamomi/plasmid DNA ratio occurred between 18 and 24 h. One of the primary objectives of this research was to study the effect of cinnamomins (elicitins secreted by P. cinnamomin) on inducing defence mechanism against the pathogen, as recent histological and ultra-structural studies showed that P. cinnamomi was restricted to the outer cortex root fragments pre-treated with capsicien and cryptogein, suggesting that elicitins can stimulate plant defence reactions against P. cinnamomi. To complement these studies and to have a clear view of the nature of the interaction, the role of cinnamomins in the production of the oxidative burst [ROS and ROS scavenging enzymes such as superoxide dismutase (SOD), catalase (CAT) and peroxidase (POD)] and in the defence responses was evaluated. Cork oak seedlings were pretreated with alpha-cinnamomin and then inoculated with P. cinnamomi mycelia. Results showed a significant higher production of reactive oxygen species (ROS) (H2O2 and O2•-) in elicitin and non-elicitin treated roots in interaction with P. cinnamomi in comparison to the corresponding control. The plant group inoculated with the pathogen after cinnamomin treatment showed an earlier increase in H2O2 production but this was lower as compared with that group inoculated with P. cinnamomi alone. Also, in elicitin pre-treated group generally, a lower level of O2•− production during infection was observed as compared with inoculated roots with P. cinnamomi alone without elicitin treatment. Furthermore, in this study, we evaluated activities of antioxidant enzymes upon challenge with P. cinnamomi, with and without pretreatment with alpha cinnamomin. Results indicated that the activities of defense enzymes POD, SOD and CAT increased after P. cinnamomi inoculation when compared with those in the control group. Also, in the group treated with alpha-cinnamomin followed by P. cinnamomi inoculation, a higher level of enzymatic activities was detected as compared with elicitin non-treated group, which suggest the protective effect of alpha-cinnamomin against the pathogen due to higher elevated levels of defense enzymes POD, SOD and CAT during the infection period. Furthermore, a sensitive qPCR method was applied to measure the pathogen biomass in elicited and non-elicited Q. suber roots challenged with P. cinnamomi to elucidate the effect of cinnamomins on the colonization of P. cinnamomi. Plasmid-based quantification of P. cinnamomi showed a significant decrease in accumulation of the pathogen DNA in cork oak roots after treatment with alpha and beta-cinnamomins which attest the role of cinnamomins in promoting defense responses in cork oak against P. cinnamomi invasion.

Study of influence of regulatory polymorphisms of expression in development of breast cancer

The human genome has millions of genetics variants that can affect gene expression. These variants are known as cis-regulatory variants and are responsible for intra-species phenotypic differences and individual susceptibility to disease. One of the diseases affected by cis-regulatory variants is breast cancer. Breast cancer is one of the most common cancers, with approximately 4500 new cases each year in Portugal. Breast cancer has many genes mutated and TP53 has been shown to be relevant for this disease. TP53 is one of the most commonly mutated genes in human cancer and it is involved in cell cycle regulation and apoptosis. Previous work by Maia et al has shown that TP53 has differential allelic expression (DAE), which suggests that this gene may be under the influence of cis-regulatory variants. Also, its DAE pattern is totally altered in breast tumours with normal copy number. We hypothesized that cis-regulatory variants affecting TP53 may have a role in breast cancer development and treatment. The present work aims to identify the cis-regulatory variants playing a role in TP53 expression, using in silico, in vitro and in vivo approaches. By bioinformatic tools we have identified candidate cis-regulatory variants and predicted the possible transcription factor binding sites that they affect. By EMSA we studied DNA-protein interactions in this region of TP53. The in silico analysis allowed us to identified three candidate cis-regulatory SNPs which may affect the binding of seven transcription factors. However, the EMSA experiments have not been conclusive and we have not yet confirmed whether any of the identified SNPs are associated with gene expression control of TP53. We will carry out further experiments to validate our findings.

The activity of mouse Cerberus like 2 during cardiogenesis - genetic and morphogenetic studies

Cardiogenesis is a delicate and complex process that requires the coordination of an intricate network of pathways and the different cell types. Therefore, understanding heart development at the morphogenetic level is an essential requirement to uncover the causes of congenital heart disease and to provide insight for disease therapies. Mouse Cerberus like 2 (Cerl2) has been defined as a Nodal antagonist in the node with an important role in the Left-Right (L/R) axis establishment, at the early embryonic development. As expected, Cerl2 knockout mice (Cerl2-/-) showed multiple laterality defects with associated cardiac failure. In order to identify the endogenous role of Cerl2 during heart formation independent of its described functions in the node, we accurately analyzed animals where laterality defects were not present. We thereby unravel the consequences of Cerl2 lossof- function in the heart, namely increased left ventricular thickness due to hyperplasia of cardiomyocytes and de-regulated expression of cardiac genes. Furthermore, the Cerl2 mutant neonates present impaired cardiac function. Once that the cardiac expression of Cerl2 is mostly observed in the left ventricle until around midgestration, this result suggest a specific regulatory role of Cerl2 during the formation of the left ventricular myoarchitecture. Here, we present two possible molecular mechanisms underlying the cardiac Cerl2 function, the regulation of Cerl2 antagonist in activation of the TGFßs/Nodal/Activin/Smad2 signaling identified by increased Smad2 phosphorilation in Cerl2-/- hearts and the negative feedback between Cerl2 and Wnt/ß-catenin signaling in heart formation. In this work and since embryonic stem cells derived from 129 mice strain is extensively used to produce targeted mutants, we also present echocardiographic reference values to progressive use of juveniles and young adult 129/Sv strain in cardiac studies. In addition, we investigate the cardiac physiology of the surviving Cerl2 mutants in 129/Sv background over time through a follow-up study using echocardiographic analysis. Our results revealed that Cerl2-/- mice are able to improve and maintain the diastolic and most of systolic cardiac physiologic parameters as analyzed until young adult age. Since Cerl2 is no longer expressed in the postnatal heart, we suggest that an intrinsic and compensatory mechanism of adaptation may be active for recovering the decreased cardiac function found in Cerl2 mutant neonates. Altogether, these data highlight the role of Cerl2 during embryonic heart development in mice. Furthermore, we also suggest that Cerl2-/- may be an interesting model to uncover the molecular, cellular and physiological mechanisms behind the improvement of the cardiac function, contributing to the development of therapeutic approaches to treat heart failures.

Evolution, gene regulation and functional analysis of BMP2 in fish

Bone morphogenetic proteins (BMPs) are multifunctional growth factors belonging to the transforming growth factor β (TGFβ) superfamily with a central role in bone formation and mineralization. BMP2, a founding member of this family, has demonstrated remarkable osteogenic properties and is clinically used to promote bone repair and fracture healing. Lack of basic data on factors regulating BMP2 expression and activity have hampered a better understanding of its role in bone formation and bone-related diseases. The objective of this work was to collect new functional data and determine spatiotemporal expression patterns in a fish system aiming towards a better understanding of BMP2 function and regulation. Transcriptional and post-transcriptional regulation of gilthead seabream BMP2 gene was inferred from luciferase reporter systems. Several bone- and cartilage-related transcription factors (e.g. RUNX3, MEF2c, SOX9 and ETS1) were found to regulate BMP2 transcription, while microRNA 20a was shown to affect stability of the BMP2 transcript and thus the mineralogenic capacity of fish bone-derived host cells. The regulation of BMP2 activity through an interaction with the matrix Gla protein (MGP) was investigated in vitro using BMP responsive elements (BRE) coupled to luciferase reporter gene. Although we demonstrated the functionality of the experimental system in a fish cell line and the activation of BMP signaling pathway by seabream BMP2, no conclusive evidence could be collected on a possible interaction beween MGP and BMP2. The evolutionary relationship among the members of BMP2/4/16 subfamily was inferred from taxonomic and phylogenetic analyses. BMP16 diverged prior to BMP2 and BMP4 and should be the result of an ancient genome duplication that occurred early in vertebrate evolution. Structural and functional data suggested that all three proteins are effectors of the BMP signaling pathway, but expression data revealed different spatiotemporal patterns in teleost fish suggesting distinct mechanisms of regulation. In this work, through the collection of novel data, we provide additional insight into the regulation, the structure and the phylogenetic relationship of BMP2 and its closely related family members.

Molecular characterization of Gla-rich protein (GRP) gene expression and function

Gla-rich protein (GRP) is a vitamin K-dependent protein related to bone and cartilage recently described. This protein is characterized by a large number of Gla (γ-carboxyglutamic acid) residues being the protein with the highest Gla content of any known protein. It was found in a widely variety of tissues but highest levels was found in skeletal and cartilaginous tissues. This small secreted protein was also expressed and accumulated in soft tissues and it was clearly associated with calcification pathologies in the same tissues. Although the biological importance of GRP remains to be elucidated, it was suggested a physiological role in cartilage development and calcification process during vertebrate skeleton formation. Using zebrafish, an accepted model to study skeletal development, we have described two grp paralog genes, grp1 and grp2, which exhibited distinct patterns of expression, suggesting different regulatory pathways for each gene. Gene synteny analysis showed that grp2 gene is more closely related to tetrapod grp, although grp1 gene was proposed to be the vertebrate ortholog by sequence comparison. In addition, we identified a functional promoter of grp2 gene and using a functional approach we confirmed the involvement of transcription factors from Sox family (Sox9b and Sox10) in the regulation of grp2 expression. In an effort to provide more information about the function of grp isoforms, we generated two zebrafish transgenic lines capable to overexpress conditionally grp genes and possible roles in the skeleton development were studied. To better understand GRP function a mammalian system was used and the analysis of knockout mice showed that GRP is involved in chondrocyte maturation and the absence of GRP is associated to proteoglycans loss in calcified articular cartilage. In addition, we detected differences in chondrogenesis markers in articular chondrocyte primary culture. Overall, our data suggest a main role for GRP on chondrocyte differentiation.

Caraterização dos supressores p20 e p23 do Citrus tristeza virus e avaliação da possibilidade da sua inibição

No presente trabalho pretendeu-se caracterizar a capacidade supressora das proteínas p20 e p23 de diferentes grupos filogenéticos do CTV e o possível silenciamento da p23 de forma a ser incluída numa estratégia de proteção. A atividade supressora local da proteína p23 de todos os grupos filogenéticos foi caracterizada em Nicotiana benthamiana da linha 16C. Todas as proteínas testadas foram capaz de suprimir o silenciamento local, mas não o silenciamento a curta distância. A supressão local mais eficiente verificou-se para a p23 do Gp 5 e a menos eficiente para os isolados dos Gps 2 e M. Surpreendentemente, a p23 do Gp 5 aboliu completamente o silenciamento sistémico, sugerindo que existe uma relação entre a intensidade do silenciamento local e sistémico. A capacidade supressora local conjunta das proteínas p20 e p23 foi avaliada. A coexpressão de ambas as proteínas revelou atividade supressora mais forte comparada com a capacidade de cada proteína individual, mesmo quando inoculada com metade da densidade ótica, sugerindo a existência de sinergismo entre as proteínas p20 e p23. Para analisar as propriedades supressoras a longo prazo, as proteínas p20 e p23 foram inseridas no vetor viral TRV que assegurou a sua disseminação pela planta e expressão por um período mais alargado. Foram observados sintomas em N. benthamiana para todas as modalidades testadas, tais como, nanismo da planta, lesões necróticas severas nas folhas inoculadas e nas folhas novas ligeiros sintomas de mosaico e enrolamento. Contudo, sistemicamente não foram registadas diferenças na capacidade supressora das proteínas p20 e p23. A possibilidade para silenciar sistemicamente a proteína p23 quando incluída num genoma viral foi avaliada através do uso de plantas e enxertos transgénicos para a p23. A estratégia que envolve o uso de enxertos transgénicos parece indicar resultados promissores que conduzem ao silenciamento da p23, contudo, são resultados que devem ser encarados como preliminares.

Molecular structure and functional analysis of runx3 in zebrafish

Tese de doutoramento, Ciências Biomédicas, Universidade do Algarve, Departamento de Ciências Biomédicas e Medicina, 2014

Cloning of the Bone Gla Protein gene from the teleost fish Sparus aurata (gilthead seabream). Molecular organization, developmental appearance and evolutionary implications

Tese de Doutoramento, Biologia Molecular, Faculdade de Ciências do Mar e do Ambiente, Universidade do Algarve, 2001

Fracionamento por processos de membranas do soro de queijo com vista ao uso na produção de bioetanol

Dissertação de mestrado, Engenharia Biológica, Faculdade de Ciências e da Tecnologia, Universidade do Algarve, 2015

Effects of mixture of nanoparticles under different salinity in the clams Ruditapes philippinarum and Ruditapes decussatus

Dissertação de mestrado, Biologia Marinha, Faculdade de Ciências e Tecnologia, Universidadde do Algarve, 2015

Identification of novel molecular determinants of tissue mineralization in fish

The identification of genes involved in signaling and regulatory pathways, and matrix formation is paramount to the better understanding of the complex mechanisms of bone formation and mineralization, and critical to the successful development of therapies for human skeletal disorders. To achieve this objective, in vitro cell systems derived from skeletal tissues and able to mineralize their extracellular matrix have been used to identify genes differentially expressed during mineralization and possibly new markers of bone and cartilage homeostasis. Using cell systems of fish origin and techniques such as suppression subtractive hybridization and microarray hybridization, three genes never associated with mechanisms of calcification were identified: the calcium binding protein S100-like, the short-chain dehydrogenase/reductase sdr-like and the betaine homocysteine S-methyltransferase bhmt3. Analysis of the spatial-temporal expression of these 3 genes by qPCR and in situ hybridization revealed: (1) the up-regulation of sdr-like transcript during in vitro mineralization of gilthead seabream cell lines and its specificity for calcified tissues and differentiating osteoblasts; (2) the up-regulation of S100-like and the down-regulation of bhmt3 during in vitro mineralization and the central role of both genes in cartilaginous tissues undergoing endo/perichondral mineralization in juvenile fish. While expression of S100-like and bhmt3 was restricted to calcified tissues, sdr-like transcript was also detected in soft tissues, in particular in tissues of the gastrointestinal tract. Functional analysis of gene promoters revealed the transcriptional regulation of the 3 genes by known regulators of osteoblast and chondrocyte differentiation/mineralization: RUNX2 and RAR (sdr-like), ETS1 (s100-like; bhmt3), SP1 and MEF2c (bhmt3). The evolutionary relationship of the different orthologs and paralogs identified within the scope of this work was also inferred from taxonomic and phylogenetic analyses and revealed novel protein subfamilies (S100-like and Sdr-like) and the explosive diversity of Bhmt family in particular fish groups (Neoteleostei). Altogether our results contribute with new data on SDR, S100 and BHMT proteins, evidencing for the first time the role for these three proteins in mechanisms of mineralization in fish and emphasized their potential as markers of mineralizing cartilage and bone in developing fish.