國立台灣大學生命科學院動物學研究所 碩士論文
Graduate Institute of Zoology College of Life Science National Taiwan University
Master’s Thesis
斑馬魚水解磷酸酯受器 4 與 5 之研究
Identification of Zebrafish LPA Receptor 4 and 5
高笙詠
Milky Shen-Yung Kao
指導教授:李心予 博士 Advisor: Hsinyu Lee, PhD
中華民國九十九年八月 August, 2010
國 立 臺 灣 大 學 動 物 學 研 究 所
碩 士 論 文 斑 馬 魚 水 解 磷 酸 酯4 與5 之 研
究
高 笙詠
99 8
II
III
Acknowledgement
大學畢業三年後,終於回到學生身分,我告訴自己既然要做研究就要做好,所以凡事都 是使盡全力的去做,直到累翻了之後,睡一覺醒來再繼續。但終究實驗不是靠毅力就會順利 的,一路跌跌撞撞從零開始,我們也是把斑馬魚養了起來。其實我從來沒把魚養活過的,來 到這裡才一步一步的學著養魚,到我畢業時我的床頭也養了兩缸魚。
首先要感謝李心予老師這段時間給我的指導,放心的將斑馬魚的計畫交給我,並且全心 的推動我做任何可能的新嘗試。還有感謝張百恩老師在斑馬魚養殖和顯微注射方面的傾囊相 授,才讓我們能在一年時間之內建立斑馬魚系統。另外則是感謝廖文亮老師在飼料餵養方面 的建議,黃泉貞老師提供基因轉殖魚,以及李世傑老師和大貓學長所提供的建議,還有柏堅 學長和智涵學長的幫助。
這兩年的生活如果沒有實驗室同仁的幫忙,我一個人不可能把這個計畫執行至此,所以 我要感謝翠華,予農,和耆豪這段期間幫忙分攤養魚的工作。另外還要感謝沛翊的魚經分享,
易謙的美食資訊交流,芷欣 CoCo 所共同製造的歡樂笑點,還有仔仔、齊倫、傳恩、岳謙、
季霖這段時間的陪伴,多虧了你們,我才能夠在實驗最低潮的時刻繼續堅持,我們才能夠每 天期待回到歡樂的實驗室。
最後要感謝馥華,如果不是先幫你寫過了一遍論文,我也不會寫自己的這麼快,如果不 是妳從高中畢業死纏著我,現在我也不會有妳的支持。謝謝妳的不離不棄,謝謝妳的包容支 持,謝謝妳的樂觀開朗,謝謝妳!
最後要感謝我的父母,打小時候的放任培養,才造就了現在我的耐心和毅力。
謝謝你們相信我的能力,才讓我有機會回到學校重新開始。
IV
Abstract
Lysophosphatidic acid (LPA) is a simple lysophospholipid which regulates various important biological functions such as cell proliferation, migration, differentiation, and morphogenesis. The LPA receptor on cell membrane has been identified and, designated as LPA1-5. LPA1-3 are the major members of the endothelial differentiation gene (EDG) family. LPA4-5, structurally distinct from LPA1-3, were first discovered as purinergic receptors in P2Y family. In zebrafish (Danio rerio), LPA1 and LPA3 have been identified essential for neural tube and lymphatic vessel development, while the function of other LPA receptors are not yet identified. We searched LPA receptors using the TBLASTN algorithm in the zebrafish nucleotide database, and found zLPA4 and one sequence similar to both LPA5 and P2Y5 in mice. After aligning all the known LPA receptors, the sequence showed more similarities with LPA5 than with P2Y5. In zebrafish, zLPA4 is expressed in the gills, testis and spleen after 5 dpf, and zLPA5 is expressed in most tissues after the 18 somites stage. LPA induced Ca++ mobilization in zLPA4 or zLPA5 over-expressed B103 rat neuroblastoma was not observed. This indicates that zebrafish LPA receptors might not trigger the conventional cascade of mammalian LPA receptors. However, down-regulation of zLPA4 leads to defections in vascular development, which needs more evidence to clarify.
Keywords: Lysophosphatidic acid, Lysophosphatidic acid receptor, blood vessel, zebrafish, development
V
摘要
水解磷酸酯為一種結構簡單的酯質卻同時調控多種生理功能,包括細胞增生、移動、
分化、以及胚胎發育。目前已知之水解磷酸酯受器有五種,其中水解磷酸酯受器一至三屬於 內皮細胞分化基因群,而水解磷酸酯受器四至五則屬於另一嘌呤受器基因群。在斑馬魚系統 中,水解磷酸酯受器一已知會調控淋巴管發育,而水解磷酸酯受器三則負責神經管發育。此 研究中,我們經由 TBLASTN 演算法比對斑馬魚與小鼠之基因序列,而取得斑馬魚水解磷酸酯 受器四與五之序列。將斑馬魚水解磷酸酯受器四與五和所有已知之水解磷酸酯受器與嘌呤受 器基因群綜合比對之後,斑馬魚水解磷酸酯受器四與五分別歸屬於水解磷酸酯受器四與五之 基因群。斑馬魚水解磷酸酯受器四表現在五天大之魚鰓、睪丸以及脾臟等處。斑馬魚水解磷 酸酯受器五則在十八體節發育期之後在全身表現。斑馬魚水解磷酸酯受器四與五異位表現在 大鼠細胞株 B103 中無法接受水解磷酸酯調控細胞內鈣離子濃度。然而在斑馬魚胚胎中降低斑 馬魚水解磷酸酯受器四之表現卻明顯的抑制斑馬魚之血管發育,導致胚胎嚴重水腫。綜合以 上實驗,斑馬魚水解磷酸酯受器四應為調控斑馬魚循環發育之重要基因。
關鍵字: 水解磷酸酯、水解磷酸酯受器、斑馬魚、血管、發育
VI
Catalog
Introduction ... 1
Metabolism ... 3
EDG family LPA receptors ... 6
Discovering the “non-Edg” LPA receptor LPAR4 ... 9
LPA5 ... 10
Biological functions of LPA ... 11
Nervous system ... 12
Vascular system ... 12
Immune system ... 13
Cancer... 13
Reproduction ... 14
Gene expression and physiological function of LPA4 and LPA5 ... 14
Expression ... 15
Reproduction ... 15
Neural physiology ... 16
Platelets ... 16
VII
Interconnected mutual regulation ... 17
Zebrafish as a model for LPA receptor research ... 18
Material and Methods ... 20
Protein sequence alignment ... 20
Phylogenetic tree and percentage similarity analysis ... 20
Expression construct of zLPA4 and zLPA5 ... 20
Stage and Tissue dependent RNA expression ... 21
Whole mount in situ hybridization ... 22
LPA induced Ca++ mobilization ... 24
Bone and cartilage differential staining ... 24
Results ... 25
Hierarchical tree of all known LPA, S1P, P2RY receptors ... 25
Protein sequence alignment of human, mouse and zebrafish LPA4 and LPA5 ... 25
Percentage similarity of LPA4 and LPA5 between species ... 26
LPA induced calcium mobilization in B103 cells through zLPA4 and zLPA5 ... 26
Developmental stage dependent zLPA4 and zLPA5 RNA expression ... 27
Whole mount in situ hybridization of zLPA4 and zLPA5 ... 27
VIII
Discussion ... 28
Reference ... 33
Figures ... 40
IX
Figures
Figure.1 Heirarchical tree of LPA4 and LPA5 ... 40
Figure.2a Phylogenetic analysis of LPA4 and LPA5 in different species. ... 42
Figure.2b Phylogenetic analysis of LPA4 and LPA5 in different species. ... 43
Figure.3 LPA induced calcium mobilization in zLPA4 or zLPA5 overexpressed B103 cells.
... 44
Figure.4 Developmental stage dependant mRNA expression of zLPA4 and zLPA5 ... 45
Figure.5 Whole mount in situ hybridization of zLPA4 and zLPA5 ... 46
1
I NTRODUCTION
Lysophosphatidic acid (LPA) is a small, widespread phospholipid that acst as an extracellular signaling molecule by activating at least six known G protein-coupled receptors (GPCRs), LPA1-LPA6, and is able to exert physiological, phathophysiological and morphological effects on many cell types. Phospholipids are water-insoluble molecules that are most abundant on the cell membrane, but they can be further metabolized into several other forms, such as effective energy sources, molecules covalently attaching to proteins, and as intercellular and intracellular signaling molecules. Most bioactive lipids acts quickly, in a matter of seconds, on membrane proteins called GPCRs in an autocrine or paracrine manner (1). GPCRs act as the primary target to 50-60% of current therapeutic drugs (2). After the completion of human genome sequence, about 800 diverse genes were found to encode GPCRs, most of which are considered to bind various extracellular ligands such as lipids, amines and amino acids (3). The human genome project also revealed many GPCRs whose endogenous ligands are still unknown and are fittingly named orphan GPCRs. Many research groups are focusing on orphan GPCRs as objects of drug discovery. In fact, many GPCRs have recently been linked to their endogenous ligands (de-orphanized), yet, there are still nearly 120 GPCRs whose endogenous ligands have not been identified (3).
Among the GPCR binding lipid molecules, Lysophospholopids (LPLs) make up the majority of the GPCR-binding lipid molecules. LPLs consist of two major compounds, LPA and Sphingosine-1-phosphate (S1P), both of which are metabolites of membrane phospholipids and are stored in large quantities in platelets, which can be released into the blood upon being wounded.
LPA and S1P are structurally similar to LPLs that have distinct receptors, but the resemblance in their tissue distribution, concentration, effector pathways and broad range of biological roles make them difficult to analyze individually. Consequently, there have been many attempts to obtain the single
2
receptor gene deficient animal so as to dissect the functions of each LPL. But unfortunately, due to a large number of receptors, it is difficult to observe significant phenotypes in knockout animals, for these receptors are functionally redundant to each other (4). Thus, investigators are working on acquiring the double or triple knockout animal (4). Nonetheless, we conclude that LPA plays an important role in animal physiology for the reason that autotaxin (ATX) knockout mice uniformly die at approximately embryonic day 9.5 with pronounced neural tube and vascular defects (5-6).
LPA is present in all eukaryotic tissues at relatively low concentrations, at higher concentrations (sub-micromolar range) in blood plasma, and even higher concentrations at wounded regions. In 1996, the first high affinity, cell surface receptor for LPA was identified as LPA1 (7), this led to the discovery of two, closely related receptors, LPA2 and LPA3, and three rather divergent receptors, LPA4, LPA5 and LPA6. Each receptor is encoded by distinct genes and named LPAR1-LPAR6 in humans and Lpar1-Lpar6 in mice. Studies utilizing cloned receptors and genetic knockout mice have presented us with notable evidence about the importance of LPA in many biological systems, such as basic cellular processes as well as multiple organ systems like the nervous system. But even 14 years after the first LPA receptor was found, there still many questions left to be answered. Which is why we are putting effort on the research of zebrafish LPA receptors. Zebrafish were used in this experiment because they are easy to handl in lage number, genetically easier to manipulate than mice and have transparent embryos that develop quickly (before the 5th day of embryonic development). In addition, the large number of progeny that the fish produce allows the chance to more frequently observe embryonic lethal defects. In this research, we obtained zebrafish LPA4 and LPA5 sequences with bioinformatic and evolutional approaches, observed their expression and the impact after down-regulation of these genes with morpholino (MO) microinjection.
3
METABOLISM
The metabolism of LPA is complicated and results in numerous, chemically distinct species.
Taking LPA as a signaling molecule, the term LPA generally stands for 1-acyl-2-hydroxyl-sn-glycero-3-phosphate, rather than other forms, such as 1-alkyl- or 2-acyl-LPA (8-9), the length of acyl chain and the degree of saturation vary vastly on the precursor phospholipid.
As an example, 1-palmitoyl-posphotidylcholine is metabolized to 1-palmitoyl-LPA (16:0-LPA). The most abundant forms of LPA in human plasma are 16:0-, 18:2-, and 18:1-LPA (10). And most of the commonly available LPA in laboratory signaling studies are probably in the form of 18:1.
LPA is produced from the membrane phospholipid through two major pathways: (a) sequential activity of phospholipase D (PLD) and phospholipase A2 (PLA2) and (b) sequential activity of PLA2 and lysophospholipase D (lysoPLD) (4), which both only in the case of producing 1-acyl-LPA. For other forms of LPA production, phospholipase A1 (PLA1) and diacylglyceral kinase (DGK) are involved in the more complicated metabolism of phospholipids, diacylglycerol, and other compounds more commonly observed in serum like lysophosphatidylcholine (LPC), lydophosphaticylethanolamine (LPE) and lysophosphatidylserine (LPS). LPC, LPE and LPA are first produced by at least two mechanisms. In activated platelets LPC, LPE and LPS are produced by secreted PLA2 and phophatidylserine (PS)-specific PLA1. In plasma, LPC is produced from phosphatidylcholine (PC) in lipoprotein by lecithincholesterol acyltransferase (LCAT) and PLA1-like enzymes (8). The generated LPLs are subsequently transformed to LPA by the lysophopholipase D activity of ATX.
Autotaxin was first found to be a cancer-cell motility-stimulating factor, as a character of its assumed nucleotide phosphodiesterase activity (11). Enhanced expression of ATX is repeatedly observed in tumor tissues such as breast cancer (12), renal cancer (13), Hodgkin lymphoma (14),
4
hepatocellular carcinoma (15) and glioblastoma (16). Furthermore, ATX is identical to a cell motility-stimulating factor isolated from cell culture supernatant, and is implicated in metastatic and invasive potential in tumor tissues (17). However, the pro-migratory effect of ATX on cancer cells now known to be resulted from the autocrine signaling of LPA production. The significance of ATX was not fully appreciated until the knockout phenotype of a previous known mouse gene Enpp2 ,encoded the ATX protein that was found to have lysoPLD activity (5). The Enpp2 double null
mice died uniformly at around embryonic day 9.5 with severe neural and vascular deficient (6).
Nevertheless, Enpp2+/- heterozygotes survive to adulthood but have plasma LPA levels half of that of wild-type mice. This confirms that ATX is the major enzyme responsible for LPA production and LPA signaling is essential for development.
In the intertangled pathways of LPA production, there is an uprising bright spot that we cannot leave about, the novel synthetic pathway of LPA mediated by mPA-PLA1α/LIPH. Both mPA-PLA1α and β, also called LIPH and LIPI, respectively, belong to the lipase gene family (18-19). They have about 50% amino acid similarity with each other. Most of the time, lipase hydrolyze fatty acids of triacylglycerol (TG) at the sn-1 and sn-3 positions. There are six extracullar and three intracellular PLA1 isozymes are known in mammals. In most of the studies, the physiological substrates of these isozymes are not known. However, some of the extracellular PLA1 isozymes such as hepatic lipase and endothelial lipase are known to hydrolyze phospholipids in addition to TG. In the phylogenetic tree of the lipase of the lipase gene family, PS- PLA1, mPA-PLA1α and mPA-PLA1β cluster and form a subfamily (20), in which "m" stands for "membrane-bound" and "PA" represents "phosphatidic acid-specific" and "PS" for "phosphatidylserine-specific".
Recent research on genetics disorders of human hair loss revealed that mPA-PLA1α/LIPH could be important for hair growth through its LPA-producing activity. Two independent groups presented
5
that human hair growth deficiency is linked to a genetic defect in the phospholipase gene mPA-PLA1α/LIPH (21-22). Kazantseva et al. identified mPA-PLA1α/LIPH as a causal gene of hair loss and hair growth defect through screening about 350,000 individuals in two populations from the Volga-Ural region of Russia (22). Ali et al. also discovered mutations in the gene for mPA-PLA1α/LIPH among four generations of a Pakistani family with hereditary human hair deficiency (21). In both of the studies, affected individuals were homozygous for a deletion or a mutation in the mPA-PLA1α/LIPH gene on chromosome 3q27. LIPH is expressed in hair follicle.
Malformation of hair follicles is always associated to a hereditary disease, the results from these studies indicate that mPA-PLA1α/LIPH regulates the formation and development of hair follicles and furthermore hair growth.
Moreover, two latter studies (23-24) showed yet another causative gene for human hair loss. By mapping the causative gene, they found homozygous truncated mutations in P2Y5, which encodes an orphan GPCR. P2Y5 is the closest homologue of LPA4, with an amino acid identity of 56%.
Consequently, Pasternack et al. showed that LPA act as a ligand for P2Y5 in reporter gene and radio-labeled ligand binding experiments. Interestingly, the expression pattern of P2Y5 resemble to that of mPA-PLA1α/LIPH where mostly in inner root sheaths of hair follicles. These results of two independent results on human genetic defects indicates that LPA produced by mPA-PLA1α/LIPH activates P2Y5 in the inner root sheaths. Although the specific pathways how LPA regulates hair follicles growth are not clearly known. These studies showed that in hair follicles, LPA signaling through mPA-PLA1α/LIPH and P2Y5 , latter known as LPA6, has a critical role in hair growth.
Furthermore, mPA-PLA1α/LIPH is a actual LPA-producing enzyme in vivo.
6
EDG FAMILY LPA RECEPTORS
Since the early twentieth century, lysophospholipids have been known to have biological activity. But the effects were thought to be the consequence of nonspecific detergent-like structure of lysophospholipids that disrupts plasma membrane. The mystery stays until a receptor was found in 1996, in a search for the ventricular zone (VZ) specific expressing gene in the cerebral cortex. Which leads to the identification of ventricular zone gene 1 (vzg-1) or endothelial differentiation gene 2 (edg-2), that was proved to be a high-affinity GPCR for LPA, latter called as LPA1 (7). Thereafter, sequence similarities granted the identification of few more S1P receptors, and two more LPA receptors, LPA2 and LPA3.
LPA1
LPA1 is the first high-affinity, wide expressing receptor identified for LPA (7). The mammalian LPAR1 gene (human chromosomal locus 9q31.3) encodes an approximately 41-kDa protein
consisting of 364 amino acids with 7 transmembrane domains. Extensively expression of Lpar1 is observed in adult mice, with significant presence in at least brain, uterus, testis, lung, small intestine, heart, stomach, kidney, spleen, thymus, placenta, and skeletal muscle (25-26), which is quite similar in humans (27). LPA1 couples with and activates several types of G proteins such as Gαi/o, Gαq/11 and Gα12/13 (28-29). Most of the cellular effects of LPA1 activation are through these G proteins, such as cell proliferation and survival, cell migration, and cytoskeletal changes; altered cell-cell contact through serum-response element activation, Ca2+ mobilization, and adenylyl cyclase inhibition; and activation of mitogen-activated protein kinase, phopholipase C, Akt, and Rho pathways (25).
In the Lpar1 knockout mice, unexpected in vivo functions of this receptor were unveiled.
Lpar1-/- mice show 50% reduction of perinatal survival in a mixed (C57Bl/6J × 129) genetic background and further decreased in pure (C57Bl/6J or Balb/cByJ) genetic backgrounds (J. Chun,
7
unpublished observations). The survived mice have a reduced body size, craniofacial dysmorphism with blunted snouts, and increased apoptosis in sciatic nerve Schwann cells (30-31). Deficient in suckling, attributed to olfactory defects, likely leads to perinatal lethality.
LPA2
Lpar2 was identified from GeneBank searches of orphan GPCR genes for the 60% amino acids
identity to Lpar1. In humans, LPAR2 (chromosomal locus 19q12) encodes a protein that has as 348 residues, with calculated molecular mass of around 39kDa (32). The expression of Lpar2 is relatively restricted compared to that of Lpar1 (25, 27). In mouse, Lpar2 is clearly present in kidney, uterus, and testis and moderately expressed in lung; and with lower levels of expression in stomach, spleen, thymus, brain, and heart (25). In mouse embryo, Lpar2 is also observed in embryonic brain but decline within a week after birth (25). In human, Lpar2 is highly expressed in testis and leukocytes, with moderate expression in prostate, spleen, thymus, and pancreas (27). In the case of tumor cells, abnormal expression of Lpar2 has been reported, pairing the receptor to tumor-promoting activities.
LPA2 activates the Gαi/o, Gαq/11 and Gα12/13, which leads to downstream activation of Ras, mitogen-activated protein kinase, phosphatidylinositol 3-kinale, Rac, phospholipase C, diacylglycerol, and Rho, which is similar to LPA1 (30). LPA2 activation is generally associated with cell survival (33-34) and cell migration (35-37), which makes it a stimulating factor for cancer metastasis (38-39).
Moreover, LPA2 has been reported to have interaction with proteins other than G proteins. For example, LPA2 promotes cell migration through interactions with focal adhesion molecule TRIP6 (40-41), and several PDZ proteins and zinc finger proteins are also reported to interact directly with the carboxyl-terminal tail of LPA2 (42). Furthermore, LPA2 mediated signaling can plays as an inhibitory effects against the EGF induced migration and invasion of pancreatic cancer cells through
8
the Gα12/13/Rho pathway (43). These evidence indicates that LPA2 signaling cross-regulates between classical G protein signaling cascades and other signaling pathways to regulate the efficiency and specificity of signal transduction.
Lpar2-/- mice are viable, roughly normal, and born at normal Mendelian rations without sexual bias. But in double knockout animal of Lpar1 and Lpar2, frontal hematomas is worse than that observed in Lpar1-/-, suggesting that LPA2 is probably developmentally redundant to LPA1.
LPA3
Lpar3 was revealed with degenerate PCR-based cloning and homology searches as an orphan
GPCR gene. LPAR3 (human chromosomal locus 1p22.3-p31.1) encodes a around 40kDa GPCR which has ~50% of amino acid sequence similarity with LPA1 and LPA2. LPAR3 is expressed in human heart, testis, prostate, pancreas, lung, ovary, and brain (44-45). In mouse, expression is observed in testis, kidney, lung, small intestine, heart, stomach, spleen, brain and thymus (25).
Interestingly, Lpar3 mRNA was shown to specifically expressing in the luminal endometrial epithelium at ovum implantation (46), and the expression is regulated by progesterone and estrogen (47).
LPA3 couples with Gαi/o and Gαq but not Gα12/13, through which it stimulates LPA mediated phopholipase C activation and calcium mobilization, adenylyl cyclase inhibition and activation, and mitogen-activated protein kinase activation but not cell rounding in neuronal cells. Furthermore, LPA3 is less sensitive to saturated acyl LPA like other EDG receptors do, but has a relative higher preference for 2-acyl-LPA unsaturated fatty acids. Taking the fact that we commonly use 18:1 LPA for signaling experiments into consideration, the differences of ligand affinity between cognate LPA receptors might affect our understanding of the physiological roles of LPA receptors in vivo, unless we have further information on the natural occurrence of LPA with different isoforms.
9
Lpar3-/- mice are grossly normal and viable, but recent report showed that in female nulls, it is
important for embryo implantation, spacing and prenatal development as well (46). However, on the basis of LPA3 expression in the frontal cortex, hippocampus, and amygdala (44-45), no significant neuronal morphology were observed in the Lpar3 null mice.
DISCOVERING THE “NON-EDG”LPA RECEPTOR LPAR4
LPA4 was originally identified as an orphan GPCR from an exploration of expressed sequence tag database, and was further found to be a specific receptor for LPA through ligand screening (48).
In this research, Noguchi K et al. stably expressed orphan GPCRs in Chinese hamster ovary (CHO) cells membrane, with N-terminus hemagglutinin (HA)-tagging, and screened each of the clones with 198 lipids. By flow cytometry with anti-HA antibody, they not only confirmed significant membrane expression of each orphan GPCR but also found 1-olelyl-LPA induced both intracellular calcium mobilization and cyclic adenosine monophosphate (cAMP) formation through P2y9, also called GPR23. Competitive radioligand binding assays revealed that p2y9/GPR23 bound various LPA species but with highest affinity to 18:1-LPA. Ectopically stable expression of p2y9/GPR23 in B103 cells, a rat neuroblastoma cell line with very low endogenous LPA receptor expression, further explained LPA dependent cellular response through p2y9/GPR23 (49-50), hereafter known as LPA4.
Unlike the other three LPA receptors, LPA4 belongs to another evolutionary branch called P2Y purinergic receptors, and is encoded by a single exon. The P2Y family consists of mainly G protein coupling receptors for nucleotides and almost exists in all human tissues. However, LPA4 cannot be activated by any nucleotides or nucleosides tested, but only by LPA. The evolutionary basis of how a structurally specified purinergic receptor can binds and responds to a long chain glyceric phospholipid is still unknown, which is one of the questions we are trying to approach in this research.
10
LPAR4 gene is located on human chromosome X region q13-q21.1, with a intronless single exon
encoding 370 amino acids protein of molecular weight around 42 kDa (51-52). Intriguingly, LPA4
specifically binds to 18:1-LPA and with a structural analogs preference order of 18:1->18:0->16:0->14:0->1-akyl->1-alkenyl-LPA (48), which is similar to that reported in EDG LPA receptors. Nevertheless, in the case of human LPA4, it only shares 20-24% amino acid identity with three other EDG family LPA receptors.
In LPA4-overexpressing cells, LPA induces cellular responses such as cell rounding and stress fiber formation through the Gα12/13 and Rho/Rho-kinase pathways (49-50), as observed in LPA1, LPA2 and LPA5. LPA4 also triggers intracellular cAMP accumulation through Gαs, and calcium mobilization through Gαq/11 and Gαi (49-50), which may be a cell type dependent activity (50).
Interestingly, LPA4 is currently the only known LPA receptor that couples to Gαs.
LPA4 RNA is expressed in mouse heart, skin, thymus, ovary, developing brain and embryonic fibroblasts (49, 53). In human, similar expression pattern was reported except LPAR4 is massively expressed in the ovary (48), which may related to its genomic location of X chromosome, and refers to its role in embryonic development we will discuss later.
Recently, LPA4 double null mice have been reported, though the animal showed no apparent abnormalities. However, Lpar4-/- fibroblasts are hypersensitized to LPA induced migration response, indicating that LPA4 plays an inhibitory role in LPA induced migration (53). Which further confirmed in LPA1 expressing B103 cells, that heterologous expressing LPA4 suppresses LPA1-dependent migration and LPA-induced migration and invasion of colon cancer cells.
LPA5
GPR92/GPR93 was an orphan GPCR closely related to LPA4, with 35% amino acid identity, which also belongs to P2Y family. Human LPAR5 is located on chromosome 12q13.31 and encodes
11
372 amino acids with around 41 kDa molecular weight. Chun and colleagues showed that heterologous expression of LPA5 exerted LPA dependent activity in B103 and RH7777 cells, both with very low or no endogenous LPA response. LPA5 couples to Gαq/11 and Gα12/13, which further stimulates calcium mobilization, cAMP production and inositol phosphate elevation as well (54-55).
Moreover, LPA5 was shown to be related to LPA stimulated neurite retraction, receptor internalization and phosphatidylinositol hydrolysis (55). In murine tissue, Lpar5 is broadly expressed in embryonic brain, small intestine, skin, spleen, stomach, thymus, lung, heart, liver, and embryonic stem cells (54-55).
In a calcium mobilization assay, Fujuwara et al. screened several lipid or non-lipid compounds with site-directed mutagenesis of LPA5 residues or computational binding pocket fitting of LPA5 protein structure (56-57). They showed the rank order of LPA5 agonists is alkyl glycerol phosphate >
LPA > farnesyl phosphates >> N-arachidonoylglycine. In the research, they confirmed LPA5 to be a bona fide LPA receptor, and also found several LPA5 agonists and two non-lipid antagonists as well.
Although LPA5 has been proved to be activated by 18:1-LPA, subsequent studies also indicates that LPA5 can response to farnesyl pyrophosphate at a much higher concentrations relative to 18:1-LPA (57-58). The fact that LPAR5 has varied affinity to varied isoforms of LPA brings the question up-front: Whether distinct ligand affinity related to the different localization and expression pattern of LPA receptors, and furthermore, the interconnected regulation of one LPA receptor to each other?
BIOLOGICAL FUNCTIONS OF LPA
LPA was first known as a blood pressure regulator in 1978 (59-60). Whereas, the physiological regulative function of LPA was attributed to its role as a metabolic intermediate. Most of the physiological functions of LPA have since been discovered to be mediated by specific interaction
12
with its six cognate receptors (LPA1-6) (29, 61). Numerous gain- and loss-of-function studies have revealed a variety of biological, receptor dependent functions, such as mitogenic effects, cell differentiation, cell survival, cytoskeletal reorganization, process retraction, and cell migration (29, 61-62). These LPA-mediated activities take parts into nervous system function, vascular
development, immune system function, cancer, reproduction, fibrosis, and obesity. In this paragraph, we will have an abstract idea on those discoveries.
Nervous system
The nervous system is one of the major loci for LPA receptor expression, and LPA presents in the brain at relatively high concentrations (9). LPA1 was first found to be specifically expressed within the proliferative cortical VZ of the embryonic brain (7, 28), suggesting the importance of LPA signaling in the development of VZ cells. Moreover, LPA receptors are expressed in almost all cell types of the nervous system. And LPA signaling were known to affect several neuronal development process, such as cortical development and function (63-64), growth and folding of the cerebral cortex (65), growth cone and process retraction (66-68), survival (65), migration (69), adhesion (31), and proliferation (30, 65). Furthermore, LPA signaling may be involved in the pathophysiology of schizophrenia and neuropathic pain (70-72)
Vascular system
The regulation of vascular system involved several aspects, from the proliferation, migration, adhesion, differentiation, and assembly of vascular endothelial cells and vascular smooth muscle cells (VSMCs). Recent analysis of knockout mice showed the frequency of frontal cephalic hemorrhages in Lpar1-/- mice, increased in Lpar1-/-/ Lpar2-/- double null mice (30, 73). Moreover, ATX-null mice embryo dies with neuronal and vascular deficit at 9.5 embryonic day. These in vivo results pinpoint the value of LPA signaling in vascular development. Besides vessel structure, LPA1
13
and LPA3 were found to be the primary mediators of LPA-induced platelet activation in a pharmacological study (74). However opposing results were also reported that LPA may inhibit platelet activation in mice (75). On the regulation of VSMCs, LPA1 and LPA2 were shown to exhibit opposing effects on the migration of primary VSMCs derived from knockout mice (35). In concert with numerous bioactive mediators, including growth factors and cytokines, LPA was known to mediates wound healing (29, 61, 76-77).
Immune system
Undoubtedly, LPA is the one of the major constituent of serum that represents blood plasma upon wound. Base on the well established knowledge that LPA takes parts in wound healing, it is inevitable that people would have an interest in analyzing its relationship with the immune system.
LPA receptors are expressed in immune cells, such as lymphocytes (34, 37, 78-79), dendritic cells (DCs) (80-81), and lymphoid organs such as the spleen and thymus (29, 54, 56). In the modulation of T cells and DCs activation, LPA was shown to excite or attenuate T cell activity under different cell activation-state and varied LPA receptor combinations. Moreover, in a recent study, ATX was presented to regulate lymphocyte trafficking (82). Additional research utilizing multiple receptor deficient mice will certainly help on the detailed understanding LPA signaling in immune system
Cancer
Since LPA is known to regulate angiogenesis and lymphangiogenesis in the case of wound healing, the vessel regulating role of LPA is further tested in cancer, for elevated vessel structure is generally observed in tumors. The first clue of this correlation was reported in 1964 that LPA's metabolic precursor LPC (lysolecithin) was found to be massively increased in the serum of ovarian cancer patient (83). In a latter study, ATX was identified as a motility-stimulating factor for cancer cells (17), and which further found as a result for its lysoPLD activity in 2002 (84-85). The
14
involvement of LPA in the cancer progression and metastasis is generally due to its vessel regulating activity, which appears to be modulated at the level of receptor expression rather than ligand concentration (39, 86). Furthermore, the tumorigenic activity of LPA were shown to mostly resulted from hypoxia, which is generally observed in tumor (87).
Reproduction
Lately, numerous studies have furnished considerable evidence for the involvement of LPA signaling in the process of reproduction, such as spermatogenesis, male sexual function, ovarian function, fertilization, embryo spacing, implantation, decidualization, pregnancy maintenance, parturition, and related diseases (26). LPA is present in the follicular fluid of healthy individuals (88), and ATX activity is elevated in serum of normal pregnant women at the third trimester of pregnancy, which is further enhanced in patients at risk for preterm delivery (84, 89). In a triple gene deleted animal of LPA1-3, significant defects in germ cell survival and azoospermia were reported (90).
Moreover, LPA3-deficient female mice showed delayed implantation, embryo crowding, and reduced litter size (46), which is similar to those reported in cyclooxygenase-2 (prostaglandin producing enzyme) deficient mice. And the LPA3-dependent phenotype can be rescued by exogenous administration of prostaglandins (46), indicating that separable LPA3 signaling occur during implantation upstream of prostglandin. The specific mechanism of LPA3 regulating spacing remain to be clarified but may involve cytosolic phospholipase A2 (cPLA2α) or Wnt/β-catenin signaling for the deficient female mice show similar embryo-crowding phenotypes as that observed in Lpar3-/- mice (91-92).
GENE EXPRESSION AND PHYSIOLOGICAL FUNCTION OF LPA4 AND LPA5
Although LPA4 was originally found in 1997 (51-52), the physiological functions of LPA4 and LPA5 were still unclear, after 12 years of their discovery (4). However, as our knowledge of the EDG
15
LPA receptors mounting up, more and more groups are presenting striking results of P2Y LPA receptors, with corresponding comparison between LPA receptors.
EXPRESSION
In an expression profiling research paper, Noji et al. presented the expression patterns of LPA receptors during mouse early development in a whole-mount in situ hybridization experiment (93).
LPA4 was expressed in the prospective midbrain-hindbrain boundary (MHB), with a parallel pattern reported for ATX. At a later stage, around E12.5, the MHB expression diminished to MHB vicinity, but elevated expression was observed at hindbrain, maxilllary processes, and first bonchial arches, limb buds, developing liver and somites. Along with time, the limb bud LPA4 expression became restricted to its proximal region, and then prospective zeugopod region and the tips of the developing digits as well. On the other hand, LPA5 was shown to express in the forebrain, rostral midbrain and margins of the neural folds, just before they were to close at the midline (93). LPA5 was reported to be ubiquitously expressed at E12.5, and relatively elevated in the craniofacial region. These reports further underscored the importance of LPA4 and LPA5 in embryonic development, especially their involvement in neural folding and midbrain formation, and the digits growth.
REPRODUCTION
As Yanagida et al. presented in a review article recently (94), the notion of LPA4 and LPA5 physiological function were brought to the surface. In adult mice, LPA4 was observed in brain, heart, lung, thymus, kidney, and skeletal muscle, and intriguingly, highly expressed in ovary, uterus and placenta of female mice. The ovary centered LPA4 expression was also reported in human (48). It is therefore expected that LPA4 might plays a role in the female reproduction. Although no notable change was observed during menstrual cycle, LPA4 mRNA expression levels increased along with pregnancy at the uterus implantation site but not ovary (94). At 12.5 days post coitum (dpc), embryos
16
showed approximately 4.5-fold higher expression of LPA4 mRNA than the placenta. In combination of the previous expression profile reports, LPA4 may involve in several developmental process such as implantation and embryonic development of brain and digits. Nevertheless, the physiological functions of elevated LPA4 expression in female ovary are still unknown.
NEURAL PHYSIOLOGY
The expressed sequence tag cDNA encoding LPA4 was originally isolated from human brain (52). As discussed above, LPA was shown to stimulate neurite retraction and cell aggragation of LPA4 stably expressing B103 cells in a Gα12/13-Rho-ROCK dependent pathway (49). Moreover, LPA4 mRNA was consistently observed both in undifferentiated human embryonic stem cells and during their differentiation toward neural stem/progenitor cells (95). On the other hand, LPA5 was reported to express in the dorsal root ganglia (DRG) (55-56), which predominantly expressing transient receptor potential vanilloid type 1 (TRPV1) channel. In this research Moon MJ et al. presented that the DRG calcium mobilization stimulated farnesyl pyrophosphate (FPP) can be inhibited by transfection of siRNA of LPA5 (56). Taken together, LPA4 may be implicated in neurodevelopmental process, whereas LPA5 may play a role in the modulation of pain and nociceptive transmission.
PLATELETS
Low-density lipoprotein (LDL) with oxidative modification, which exists in the circulation or exposed after rupture of atherosclerotic plaques to blood cells, was known to activate platelets (96).
LPA is one of the main platelet-activating lipid of oxidized LDL and human atherosclerotic lesions (97), and was known to induce platelet shape change (98) or platelet aggregation at higher concentration (96). However, the platelet stimulating activity of LPA requires the synergism of adenosine diphosphate (ADP), which mediates the positive feedback of platelet activation (96). Since platelet is the main source of LPA, it is plausible that LPA may play a role as autocrine signal that
17
triggers and amplifies platelet activation response. Nevertheless, platelets act differently with varied isoforms of LPA (99), and LPA receptors expressing in platelets are mostly LPA4 and LPA5. Combined with the singular ligand affinity of LPA5, which actually parallel to those of platelets, and the fact that LPA5 agonistic activity consistently associated with the platelet shape change, LPA5
might take parts in the activation of platelets. On the contrary, platelets in approximately 20% of individuals failed to aggregate upon LPA, which is coincident with increased LPA4 mRNA (100).
These findings are potentially of clinical interest because LPA non-responders were less likely to have coronary artery disease, and the differences of biased coagulating activity may results from diverse LPA receptor expression.
INTERCONNECTED MUTUAL REGULATION
As mentioned above, LPA4 was known to negatively regulates LPA1-dependent fibroblast migration, which was the very first evidence that indicates LPA4 may plays as a brake against other LPA receptor (53). Recently, Yaworsky et al. presented that osteoblast differentiate 4 times faster with down-regulated LPA4 (101), which is also a LPA1-dependent mechanism. Moreover, LPA4
deficient mice have increased trabecular bone volume, number and thickness as well. On the other hand, more recently, Sims et al. reported LPA may also activate osteoclasts cyotosolic calcium concentration, evoke retraction and promote cell survival in a LPA1-dependent fashion (J Biol Chem. 2010 Jun 15. Epub ahead of print). However, the retraction activity of osteoclasts require the collaboration of another LPA receptor that stimulates Gα12/13, which is yet unknown, despite LPA2, LPA4 and LPA5 are known to exist in osteoclasts. Collectively, these evidences support the impression that cognate LPA receptors may regulate each other in order to exhibit delicate control of LPA activities.
18
ZEBRAFISH AS A MODEL FOR LPA RECEPTOR RESEARCH
To date, four of the known LPA receptors (LPA1-LPA4) have been separately disrupted in mice (30, 46, 53, 64, 70, 73). LPA1/LPA2 double-deficient and LPA1/LPA2/LPA3 triple-deficient strains have also been generated and analyzed (73, 90). However, none of these deficient mice exhibited a phenotype similar to ATX-deficient mice, which with embryonic lethality due to vascular defects (5-6). Therefore, different combinations of LPA receptors, especially non-Edg receptors might be necessary to recapitulate the phenotype of LPA deficient, and possibly, other subtypes of LPA receptors may remain to be involved, such as GPR87 and P2Y10 (102-103).
Basing on the knowledge of these evidences mainly observed in mammals, there still a lot yet left unknown, especially for P2Y family LPA receptors. The main obstacle in mammal LPA research is that the screening of deficient phenotypes is limited to the offspring number of null mice, despite the LPA receptor deficient may affect on the progeny survival and the difficulty of mouse in vivo assay system. Therefore, here in this study, we aim to establish a zebrafish animal model that can be genetically manipulated with large amount and analyzed in vivo, which also as a tool to reveal developmental defects that cannot be seen in the mouse uterus.
As Lee et al. presented earlier (104), zebrafish LPA1 is known to be essential for lymphatic vessel formation, and zLPA3 is critical for neural tube development (unpublished result from Lee).
Previous reports indicate that the phenotype of LPA receptor deficient in zebrafish may not necessarily parallel to those observed in mammals, which suggest that the signaling pathway might differ from fish to mice. However, the physiological function observed in zebrafish actually match those presented in ATX deficient mice, but with a much longer observation window (5 days), and outnumbered experimental repeats
19
In this study, we intended to identify and analyze the developmental and physiological function of zebrafish LPA4 and LPA5. In combination of bioinformatics and evolutionary analysis, we identified zLPA4 and zLPA5. The mRNA expression of these genes were monitored with real-time PCR and whole-mount in situ hybridization techniques. By morpholino (MO) microinjection of each gene, we observed developmental defects in zLPA4 or zLPA5 down-regulated embryos.
20
M ATERIAL AND M ETHODS
PROTEIN SEQUENCE ALIGNMENT
By using mouse LPA receptor 4 and 5 protein sequence as bait, we search for zebrafish LPA receptor 4 and 5 with TBLASTN analysis in the zebrafish nucleotide collection database. Zebrafish LPA4 is found as XP_001334713.1, and zLPA5 is found as NP_955900.1. The latter one have low similarity of about 30% with both mouse LPA5 and LPA6, but after aligning with all known the LPA, S1P receptors and P2RY receptors of different species, the sequence has highest similarity to the LPA5 family. Also, zebrafish LPA6 is identified by another group with morpholino attenuation experiment (105).
PHYLOGENETIC TREE AND PERCENTAGE SIMILARITY ANALYSIS
All the known LPA, S1P receptors and P2RY receptors are aligned with ClustalX2, and tree was drew by bootstrapping analysis with one thousand random trial. Confident reference are sited beside each branch as number of percentage. Each receptors of different species was clustered into separated group of each receptor family. The analysis is done in both protein and nucleotide sequences, and reveals similar structure that the two identified sequences belonged to LPA4 and LPA5
family. Percentage similarity of LPA4, LPA5 and LPA6 between human, mouse and zebrafish are collected and presented into chart.
EXPRESSION CONSTRUCT OF ZLPA4 AND ZLPA5
Zebrafish LPA4 is amplified from 3 dpf zebrafish WT cDNA with Phusion DNA polymerase (New England Biolab) by using forward primer: ATGGCCAGTCTTGTTCTTAA and reverse primer: TCAGAACTGAGTCTCACCAA. And zLPA5 is amplified with forward primer:
ATGACTTCAAACAACACTAC and reverse primer: TTATTGTCCAGTCCAACTCG. Both CDS
21
were ligated into pGEM T easy vector with additional adenosine to each end of the blunt end PCR product. Construct were selected and sequenced to obtain full-length CDS without mutation.
STAGE AND TISSUE DEPENDENT RNA EXPRESSION
At the night before mating morning, each pair of female and male zebrafish were transferred to a single tank with mesh bottom to allow egg passing through and keeps fishes away from eating the eggs. Females and males are separated from each other by a transparent acrylic plate which avoids fighting and helps fishes getting acquainted with each other. At the morning of mating, acrylic plate is removed as long as the lights turned on, and eggs are collected after the observation of tail chasing activity. Eggs are cultivated under 28.5°C in fish water, malformed and dead embryo are removed daily so as to keep the water clean. At the time of each developmental stage, embryo are observed under binomial microscope before picking 30 embryos for RNA extraction. Each 30 embryos are transferred into a 1.5 ml vial and left with 100 μl fish water, added 1 ml of TRIzol® (Invitrogen) and homogenized with 1 ml syringe with needles of different tip size, from 23G, 27G to 30G needles. The homogenized embryo were collected in TRIzol® and stored in -80°C till the time of RNA extraction.
Upon RNA extraction, each vial of 1ml TRIzol® is added with 100 μl chloroform and vortex for 30 seconds before standing on ice for 5 minutes, then undergo phase separation with 15 minutes of 12000 rpm centrifugation at 4°C. The upper aquarius phase is then transferred to a new vial and added with 2× volume of 100% Ethanol and stand for 15 minutes at -20°C followed by 15 minutes 12000 rpm centrifugation to precipitate RNA. Then carefully discarded the aquarius ethanol, and wash the RNA pallet with 70% ethanol, centrifugated 3000 rpm for 5 minutes. Again, discard the upper aquarius phase and air dry the RNA pallet in a laminar flow for 15 minutes at room temperature, and resuspended with 50 μl DEPC water. The resuspended RNA were heated on a hotplate at 55°C for 10 minutes before density observed with NanoDrop (Thermo Scientific). Two micro grams of embryo RNA were reverse-transcribed into cDNA with ReverTra Ace (TOYOBO), and the resulting 20 μl
22
cDNA were diluted 2× with DEPC water. Two micro liter of embryo cDNA were used as template to undergo real-time PCR by using ABsolute™ QPCR SYBR® Green Mixes (Thermo Scientific) and iCycler Thermal Cycler (BioRad). Gene specific primer for zLPA4, zLPA5, z beta Actin and zEF1α were designed by using pDraw32. And primer sequences as follows: zLPA4 intF' ATTGCTTCTTTGACCCGGTG, zLPA4 R' TCAGAACTGAGTCTCACCAA, zLPA5 F' ATGACTTCAAACAACACTAC, zLPA5 intR' CCATGTAGATCACTGGAACA, zbActin RT F'
CCAGCTGTCTTCCCATCCA, zbActin_RT_R' TCACCAACGTAGCTGTCTTTCTG,
zEF1α_RT_F' CTGGAGGCCAGCTCAAACAT, zEF1α_RT_R'
ATCAAGAAGAGTAGTACCGCTAGCATTAC. Acquired fluorescent reading were analyzed with iQ5 (BioRad), and CT values of reference genes were used to generalize deviation between samples.
ZLPA4 and zLPA5 RNA expression level were presented as relative fold to the average amount of zbActin and zEF1α.
Zebrafish RNA of different organs were collected in a similar procedure. In order to preserve as much RNA as possible, fish were sacrificed by laying on ice for 5 minutes and dissected under binomial microscope as long as fish's gill stop moving. Separated organs were homogenized in 1 ml TRIzol® with syringe with different tip opening as soon as possible. Collected organ RNA were then extracted with chloroform and precipitated by 100% ethanol, washed with 70% ethanol, air dried and resuspended in 50 μl DEPC treated water, followed by density determination. And revers transcribed into cDNA then used as template for real-time PCR, with zbActin and zEF1α as internal control.
WHOLE MOUNT IN SITU HYBRIDIZATION
The in situ hybridization experiment were carried out as in (106). In brief, probe were synthesized by using zLPA4 pGEM T and zLPA5 pGEM T vectors as template and undergo PCR
23
amplification by using VAS Taq Blue DNA polymerase (BIONOVAS) with primers containing 5' T7 promoter sequence: TAATACGACTCACTATAGGG. Each gene were amplified in two sets of primers in order to obtain sense and antisense template with 5' or 3' T7 promoter sequence. Primer sequences as follows, Sense: zLPR4 T7F' TAATACGACTCACTATAGGGgtggcattgacgactccttc,
zLPA4 R' TCAGAACTGAGTCTCACCAA, zLPA5 T7F'
TAATACGACTCACTATAGGGatgacttcaaacaacactac, zLPA5 R' TTATTGTCCAGTCCAACTCG,
Antisense: zLPA4 F' ATGGCCAGTCTTGTTCTTAA, zLPA4 T7R'
TAATACGACTCACTATAGGGtcagaactgagtctcaccaa, zLPA5 F' ATGACTTCAAACAACACTAC, zLPA5 T7R' TAATACGACTCACTATAGGGttattgtccagtccaactcg. (Capitalized letters represent the T7 promoter sequence.) The amplified template with opposite T7 promoter orientation are subjected to T7 RNA polymerase (Promega) with DIG RNA labeling mix (UTP) (Promega) so as to produce each sense and antisense strand RNA probe of each gene.
Embryo are collected at different developmental stages with 4% freshly prepared paraformaldehyde 4°C overnight, and then switch into 100% methanol and stored in -20°C till needed.
Before hybridization, the embryos are hybridized in PBS, permeabilized by proteinase K, washed in PBT and prehybridized in Hybridization Mix (HM) 70°C for 5 hours. For each 200 μl of HM, 40 ng of DIG-labeled RNA probe were added, and hybridized overnight at 70°C. After probe hybridization, embryos were washed through serial concentration of 2×SSC, 0.2×SSC, and then PBT. After 4 hours or blocking, sheep anti-digoxigenin-AP Fab fragments (Roche Diagnostics)(alkaline phosphatase conjugated) were added into blocking buffer at 1/10,000 dilution overnight at 4°C. For staining, the embryos were washed with PBT, changed into alkaline Tris buffer, and stained with NBT and BCIP in dark. The staining reaction were terminated by stop solution and mounted in 100% glycerol.
24
LPA INDUCED CA++ MOBILIZATION
The zLPA4 and zLPA5 over-expression construct were made by transferring each gene from pGEM T vecters into pIRES2 EGFP vectors through EcoRI site, the constructs were sequenced to confirm preferred orientation. Rat neuroblastoma cell line B103 were known to express very low level of intrinsic LPA receptors, which were transfected with zLPA4 pIRES, zLPA5 pIRES, or pIRES2 EGFP vector only by Lipofectamine™ 2000 (Invitrogene) for 24 hours with 500ng DNA/1.5μl Lipofectamine ratio per 90,000 cells in 96 wells plate. Fluorescence were observed 24 hours after transfection in order to confirm gene expression. About 48 hours after transfection, cells were stained with 4 uM Fluor 3-AM 37°C for 40 minutes, washed with PBS. The calcium mobilization were measured by Flexstation® 3 (Molecular Devices) by comparing time-laps fluorescent readings with different stimulations, such as 5μM LPA, 5μM S1P, vector controls or 1%
Triton X-100.
BONE AND CARTILAGE DIFFERENTIAL STAINING
Embryos were anaesthetized in 0.08% Tricaine® , and fixed in 4% paraformaldehyde/ PBS for 24 hours. Then washed in PBS and stained for cartilage in 0.1mg/ml Alcian Blue (Sigma) in ethanol/acetic acid (4:1), then rehydrated through ethanol series from 90%, 50% to 30%. The rehydrated embryos were digested overnight in 50mg/ml trypsin in 30% sodium tetraborate in water, and stained for bone with 0.4 ml of Alizarin Red S (Sigma) solution in 10 ml 0.5% KOH. At desired staining stage, embryos were destained in 1% KOH/glycerol series and stored in glycerol.
25
R ESULTS
HIERARCHICAL TREE OF ALL KNOWN LPA,S1P,P2RY RECEPTORS
In order to discover the zebrafish LPA receptor 4 and 5, we used the known mouse receptor protein sequences as bait and searched in zebrafish nucleotide sequence database. At the first glance, zebrafish LPA4 analogue XP_001334713.1 was found to has relatively high similarity to it's mammal parallel receptor. Whether the other sequence NP_955900.1 has an equally low similarity to both mouse LPA5 and LPA6, along with another zebrafish gene NP_001073524.1. Since both S1P/LPA are lysophospholipids, and LPA receptors can be divided into two families, that LPA1-3 are in the EDG family and LPA4-6 are in the P2RY family, we aligned all the known LPA, S1P and P2RY receptors protein sequences between species, and drew a hierarchical tree with random bootstrapping (Fig. 1). In the resulting figure, each receptor of different species were clustered into it's own family, with confidence reference sited at the branch. In accordance with our knowledge, S1PR1-5 are closely related to LPA1-3 for they all belong to EDG family, and LPA4/5 are similar to most other P2RY receptors. But each of the LPA receptors gathered in their own clusters, that XP_001334713.1 stays with other species' LPA4, and NP_955900.1 stays with other species' LPA5, but NP_001073524.1 were found to cluster with LPA6. The sequence NP_001073524.1 was further confirmed as zLPA6 by another group (105), and was found to be important for zebrafish embryo vessel development in a morpholino knock-down experiment. In this figure, we further verified the identity of zebrafish LPA receptor 4 and 5.
PROTEIN SEQUENCE ALIGNMENT OF HUMAN, MOUSE AND ZEBRAFISH LPA4 AND LPA5
All of the LPA receptors are GPCRs, thus as in Fig. 2a, we aligned the human, mouse and zebrafish protein sequences of LPA4 and LPA5, both alignment show conserved GPCR's seven transmembrane domains. But in zebrafish receptors, C terminus sequences are highly variated from
26
their mammalian relatives, in which cytosolic domains were known to be responsible for the binding of down stream signaling proteins such as the trimeric G protein complex (4) and PDZ scaffold proteins (42) such as NHERF2. These results suggest that the zebrafish LPA receptors 4 and 5 might have different cytosolic binding affinities compared to mammal's. In addition to cytosolic domain, the N terminus of zLPA4 is also 20 amino acids shorter than mammalian receptors, which their physiological function were not yet known. These finding reveals the conserved and variated domains of zLPA4 and zLPA5, which sheds light on the evolutionary research of LPA signalling.
PERCENTAGE SIMILARITY OF LPA4 AND LPA5 BETWEEN SPECIES
In the alignment of LPA receptors between different species, we collected the identical percentage of the same receptor, and arranged into this chart (Fig.2b). Compared to zebrafish LPA receptor 1 (104), LPA4 and LPA5 are much less conserved from fish to mammal, despite most of the LPA receptors are conserved between mouse and human. Especially in zLPA5, in contrast with the 60% similarity of zLPA4 and zLPA6, the zebrafish analogue has a low similarity to it's mammalian relatives for about 30%. In comparison with Fig. 2a, the LPA5 family is diverted from other R2RY LPA receptors such as LPA4 and LPA6, which indicates LPA5 might plays quite different roles in fishes and mammals, and it's physiological function may also differs from LPA4 due to it's broad ligand affinity (107).
LPA INDUCED CALCIUM MOBILIZATION IN B103 CELLS THROUGH ZLPA4 AND ZLPA5
All mammalian LPA receptors were known to trigger cytosolic calcium mobilization upon LPA stimulation (4). Therefore we over-expressed zLPA4 and zLPA5 in rat neuroblastoma cell line B103 along with a downstream IRES sequence and EGFP marker gene to monitor exogenous gene expression. The cells were subjected to calcium fluorophore staining, and observed for LPA stimulated calcium mobilization. But unfortunately, we are not able to observe any LPA specific
27
calcium deposition in this assay, by taking vector only transfection and vehicle treatment as negative control, which may due to the C terminus differences of zLPA4 and zLPA5 between fish and mammal.
DEVELOPMENTAL STAGE DEPENDENT ZLPA4 AND ZLPA5RNA EXPRESSION
In order to dissect the developmental roles of zLPA4 and zLPA5, first, we analyzed the expression timing and amount of these two genes. With the technique of real-time PCR, we used two reference genes zβActin and zEF1α as internal control, but we were surprised by the variation of reference gene expression along with fish development. Both genes transcription level increased with time, but zEF1α RNA peaks at 6 somites stage, however zβActin RNA level increases in a slower rate, which peak at 96 hpf stage. Taking these variations into account, the increased RNA of zLPA4 and zLPA5 before gastrulation can be seen as maternal RNA deposition. Thus, embryonic zLPA4 RNA is expressed at 5th day of development, and zLPA5 RNA is expressed since 24 somites stage. Both gene's RNA expression significantly increased after 5 dpf (data no shown).
WHOLE MOUNT IN SITU HYBRIDIZATION OF ZLPA4 AND ZLPA5
With the knowledge of the expression timing of these two genes, we further analyzed the spatial pattern of RNA expression at different developmental stages. With zHSP70 antisense probe as positive control, which zHSP70 is known to be highly expressed in 2 dpf embryo eyes, we are not able to observe any specific pattern between the comparison of antisense and sense probes of zLPA4 and zLPA5. Which suggest these two may not be specifically expressed in certain organ that can be observed under binomial microscope, however, we cannot exclude the possibility that these two LPA receptors may expressed at a smaller scale such as microvascular system or skin.