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.

Characterization of the cis‐regulatory regions governing ojoplano locus

Ojoplano (opo) is a vertebrate-specific gene that was first identified in medaka fish as a recessive mutant, showing both neural crest defects and a failure of optic cup folding. In humans, this gene is associated with genetic diseases including hereditary craniofacial malformations and schizophrenia. It is localized in a 2Mb gene desert flanked by insulator sequences, between the genes SLC35B and TFAp2a. This region, syntenic between all vertebrates, represents only 2% of chromosome 6. However, it includes 23% of the all conserved cis-regulatory elements in this chromosome. Using transgenesis assays in zebrafish, we screened the enhancer activity of this locus and obtain a collection of nine enhancers. These regulatory elements were all conserved from human to teleosts and showed epigenetic marks for enhancer activity. We could associate multiple enhancers with ororfacial celfting disease and in order to explore the functionality of the enhancers, we performed a bioinformatics analysis to search for transcription factor bindings in the enhancer sequences. In terms of gene regulation we observe that H6:10137 opo enhancer has two Vsx2 binding sites and that this transcription factor regulates the expression of opo during eye development. Our findings suggest that the regulation of Vsx2 over opo is essential for optic cup folding. So far, there is no clear connection between optic cup patterning and morphogenesis. Vsx2 provides this link by controlling the expression of opo.

Structure, tissue distribution and estrogen regulation of splice variants of the sea bream estrogen receptor alpha gene

Estrogen actions are mainly mediated by specific nuclear estrogen receptors (ERs), for which different genes and a diversity of transcript variants have been identified, mainly in mammals. In this study, we investigated the presence of ER splice variants in the teleost fish gilthead sea bream (Sparus auratus), by comparison with the genomic organization of the related species Takifugu rubripes. Two exon2-deleted ERα transcript variants were isolated from liver cDNA of estradiol-treated fish. The ΔE2 variant lacks ERα exon 2, generating a premature termination codon and a putative C-terminal truncated receptor, while the ΔE2,3* variant contains an in-frame deletion of exon 2 and part of exon 3 and codes for a putative ERα protein variant lacking most of the DNA-binding domain. Both variants were expressed at very low levels in several female and male sea bream tissues, and their expression was highly inducible in liver by estradiol-17β treatment with a strong positive correlation with the typical wild-type (wt) ERα response in this tissue. These findings identify novel estrogen responsive splice variants of fish ERα, and provide the basis for future studies to investigate possible modulation of wt-ER actions by splice variants.

Effect of polyunsaturated fatty acids (PUFAs) on the proliferation and extracellular matrix mineralization and gene expression of gilthead seabream skeletal cell lines

The aquaculture industry aims at replacing significant amounts of marine fish oil by vegetable oils in fish diet. Dietary lipids have been shown to alter the fatty acid composition of bone compartments, which would impact the local production of factors controlling bone formation. Knowledge on the mechanisms underlying the nutritional regulation of bone metabolism is however scarce in fish. Two in vitro bone-derived cell systems developed from seabream (an important species for aquaculture in the Mediterranean region) vertebra, capable of in vitro mineralization and exhibiting prechondrocyte (VSa13) and pre-osteoblast (VSa16) phenotype, were used to assess the effect of certain polyunsaturated fatty acids (PUFAs; arachidonic (AA), eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids) on cell proliferation, extracellular matrix (ECM) mineralization and gene expression. While all PUFAs promoted morphological changes in both cell lines, VSa16 cell proliferation appeared to be stimulated by PUFAs in a dose dependent manner until 100M, whereas proliferation of VSa13 cells was impaired at concentrations above 10M. AA, EPA and DHA inhibited VSa13 ECM mineralization, alone and in combination, while VSa16 ECM mineralization was only inhibited by AA and EPA. DHA had the opposite effect, increasing mineralization almost by 2 fold. When EFAs were combined, DHA apparently compensated for the inhibitory effect of AA and EPA. Expression of marker genes for bone and lipid metabolisms has been investigated by qPCR and shown to be regulated in pre-osteoblasts exposed to individual PUFAs. Our results show that PUFAs are effectors of fish bone cell lines, altering cell morphology, proliferation and mineralization when added to culture medium. This work also demonstrates the suitability of our in vitro cell systems to get insights into mineralization-related effects of PUFAs in vivo and to evaluate the replacement of fish oils by vegetable oil sources in fish feeds.

The role of hypoxia in the regulation of membrane transporters

ATP binding cassette (ABC) and solute carrier (SLC) transporters are responsible for the majority of the transcellular movement of various substrates, including drugs, among epithelial cells. Despite the well characterized regulation of influx (SLC) and efflux (ABC) transporters by endogenous mediators, such as inflammatory cytokines, little is known about how changes in oxygen levels may affect expression of these transporters. In this study we showed that the expression of SLC22A4, SLC22A5, SLC22A1, SLC02B1, SLC10A2, ABCC2 and ABCC3 transporters is upregulated by hypoxia in HT29 colon carcinoma cells, but not in HepG2 hepatocarcinoma cells. Moreover, OCTN1 (SLC22A4), OCT1 (SLC22A1) and OATP-B (SLC02B1) transporter expression is also induced by inflammatory cytokines but in a smaller extent than in hypoxia. Furthermore our experiments indicate that there is no cross talk between HIF-1 and NF-κB pathways in HT-29 cells, but these two pathways act simultaneously activating common genes, such as, some SLC and ABC transporters. Our preliminary results from studies with an in vivo murine model of colitis, suggest that HIF-1is stabilized and OCTN1 is strongly induced during severe inflammation, which can be relevant for a recovery from the inflammatory process. We have also been interested in the distribution of HIF-1α variants among different ethnic groups as well as their contribution for cancer risk. Thus, we have demonstrated that there is an ethnicity-related variation in the frequency of the C1772T (P582S) single nucleotide polymorphism (SNP) in the HIF-1α gene. Furthermore, we performed a case-control study in a breast cancer population and our results suggest that there is no association between this SNP or the rare G1790A (A588T) SNP and the incidence of breast cancer. Taken together, the results obtained in this study contribute to a better knowledge of drug influx and efflux during hypoxia and inflammation as well as to the understanding of the pharmacogenetic variability of the HIF-1.

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.