Global Expression Profiling Identifies a Novel Hyaluronan ... - MDPI

Journal of

Clinical Medicine Article

Global Expression Profiling Identifies a Novel Hyaluronan Synthases 2 Gene in the Pathogenesis of Lower Extremity Varicose Veins Chia-Shan Hsieh 1,2 , Chia-Ti Tsai 3,4 , Yau-Hung Chen 5 , Sheng-Nan Chang 6 , Juey-Jen Hwang 3,6 , Eric Y. Chuang 1,2, * and I-Hui Wu 4,7, * 1 2 3 4 5 6 7

*

Department of Life Science, Genome and Systems Biology Degree Program, National Taiwan University, Taipei 10617, Taiwan; [email protected] Bioinformatics and Biostatistics Core, Center of Genomic Medicine, National Taiwan University, Taipei 10055, Taiwan Division of Cardiology, Department of Internal Medicine, National Taiwan University College of Medicine and Hospital, Taipei 10002, Taiwan; [email protected] (C.-T.T.); [email protected] (J.-J.H.) Graduate Institute of Clinical Medicine, College of Medicine, National Taiwan University, Taipei 10051, Taiwan Department of Chemistry, Tamkang University, Taipei 25137, Taiwan; [email protected] Department of Internal Medicine, National Taiwan University Hospital Yun-Lin Branch, Yun-Lin 64041, Taiwan; [email protected] Department of Surgery, National Taiwan University Hospital, Taipei, 10002, Taiwan Correspondence: [email protected] (E.Y.C.); [email protected] (I-H.W.); Tel.: +886-2-3366-3660 (E.Y.C.); +886-2-23123456 (ext. 65735) (I-H.W.)

Received: 24 October 2018; Accepted: 7 December 2018; Published: 11 December 2018







Abstract: Lower extremities varicose veins (VV) are among the most easily recognized venous abnormalities. The genetic mechanism of VV is largely unknown. In this study, we sought to explore the global expressional change of VV and identify novel genes that might play a role in VV. We used next-generation ribonucleic acid (RNA) sequence (RNA seq) technology to study the global messenger RNA expressional change in the venous samples of five diseased and five control patients. We identified several differentially expressed genes, which were further confirmed by conventional reverse transcription polymerase chain reaction (RT-PCR). Using these significant genes we performed in silico pathway analyses and found distinct transcriptional networks, such as angiogenesis, cell adhesion, vascular injury, and carbohydrate metabolisms that might be involved in the mechanism of VV. Among these significant genes, we also found hyaluronan synthases 2 gene (HAS2) played a pivotal role and governed all these pathways. We further confirmed that HAS2 expression was decreased in the venous samples of patients with VV. Finally, we used a zebrafish model with fluorescence emitting vasculature and red blood cells to see the morphological changes of the venous system and blood flow. We found that HAS2 knockdown in zebrafish resulted in dilated venous structural with static venous flow. HAS2 may modulate the transcriptional networks of angiogenesis, cell adhesion, vascular injury, and carbohydrate metabolisms in venous tissues and downregulation of HAS2 may underlie the mechanism of VV. Keywords: varicose veins; HAS2; RNA seq

1. Introduction Varicose veins (VVs) of lower extremities are among the most easily recognized vascular abnormalities with superficial venous tortuosity and enlargement [1]. Clinical manifestations can

J. Clin. Med. 2018, 7, 537; doi:10.3390/jcm7120537

www.mdpi.com/journal/jcm

J. Clin. Med. 2018, 7, 537

2 of 13

vary from leg edema to chronic, disabling venous ulceration and it affects around 25% of the adult population [2]. Age, female sex, family history, pregnancy, obesity, and prolonged standing have been shown to be predisposing factors in several epidemiological studies [3–8]. The complex pathophysiology of greater saphenous vein (GSV) valvular reflux and altered extracellular matrix proteins with increased proteolytic enzyme have been addressed as potential mechanisms of VV [9–11]. Recently, it has been shown that matrix metalloproteinases may play an important role in the remodeling of venous extracellular matrix and mechanism of VV) [8,10,12–15]. However, these processes account for only some of the tangled mechanisms, and many more remain to be identified. Genetic study is a powerful tool to identify the novel mechanism of human diseases. Several candidate genes have been identified as the possible susceptibility genes for VV, such as COL1A2, FOXC2, NOTCH3, TGF-β, thrombomodulin, and VEGF1 [16–22]. However, whether these genes are true susceptibility genes for VV remains questioned. The study recently provides novel insight into the physiology of lncRNAs and the pathogenesis of varicose veins, but still need to do more for further investigation [2]. Genome-wide association studies (GWASs), a non-candidate-driven approach, investigate the entire genome and can identify unknown or novel disease-associated single nucleotide polymorphisms (SNPs) or genes. [23–26]. Recently, GWASs have identified several novel genes that may be involved in the mechanism of venous pathogenesis [13,27]. However, it is difficult to correlate the identified significant SNPs by GWAS to the mechanism of disease. Therefore, instead of traditional forward genetic approach, in the present study we used a reverse genetic approach. To this end, we first systematically investigated and compared the global expressional changes in the venous samples of diseased and control patients using the next-generation ribonucleic acid (RNA) sequence (RNA seq) technology. We identified several differentially expressed genes that might be involved in the molecular mechanism of VV. We then used a novel zebrafish model to see the phenotypic changes when these genes were perturbed. 2. Methods 2.1. Patient Selection Five patients with VV were consecutively enrolled. All patients were diagnosed and assessed by a vascular surgeon and underwent duplex ultrasonographic scanning to assess suitability for treatment and entry into the study. Inclusion criteria were an age of 18 years or older, the presence of unilateral or bilateral primary symptomatic varicose veins at CEAP classification grade 3 or higher [28] and reflux of GSV of more than 0.5 s on duplex ultrasonography after calf compression [29]. Patients with current deep-vein thrombosis, acute superficial vein thrombosis, bleeding diathesis, pregnancy, and contraindications to the use of general or regional anesthesia were excluded. This study was approved by the institutional review board of the National Taiwan University Hospital and informed consent was obtained from participating subjects. 2.2. Surgical Removal of Venous Samples GSV ablation using the ClosureFAST (VNUS Medical Technologies Inc, San Jose, CA, USA) with concomitant stab phlebectomy were done. All the procedures were performed in a two-day inpatient setting. The GSV was cannulated at the most distal point of venous reflux, and the catheter tip was positioned 2 cm distal to the saphenofemoral junction. The GSV was ablated with the concomitant phlebectomy [30,31]. The specimen collected from the avulsion phlebectomy were used for analysis. Five control patients were consecutively enrolled. The controls patients had coronary artery disease and ischemic cardiomyopathy. Patients taking steroids or with conditions that could modify leukocyte activity (e.g., rheumatoid arthritis, vasculitis, or other collagen diseases) were excluded. Their GSV samples were harvested during elective coronary artery bypass to serve as venous grafts and the redundant part was used for study. The samples were put into and stored in liquid nitrogen immediately after they were harvested. GSVs competency and diameter were screened by the

J. Clin. Med. 2018, 7, 537

3 of 13

pre-operative duplex ultrasonography. The baseline characteristics of patients with VV and controls are shown in Table 1. Table 1. Clinical characteristics of patients of the two groups Variables

VV (n = 5)

Controls (n = 5)

Age (yr) BMI Male DM HTN HLP Smoke

69.6 ± 14.7 26.1 ± 2.8 2 0 2 1 0

75 ± 11.4 22.6 ± 1.8 5 3 3 4 2

BMI, body mass index; DM, diabetes mellitus; HPL, hyperlipidemia; HTN, hypertension; VV, varicose vein; yr, years.

2.3. Sample Preparation and RNA Sequencing Total RNA was extracted from the venous samples. The tissue was homogenized with POLYTRON® PT 2100 (Dispergierund Mischtechnik, Littau, Switzerland), and Trizol solution (Gibco BRL, Grand Island, NY, USA) was added for RNA extraction. The extracted RNA was dissolved in diethyl pyrocarbonate-treated distilled water. Spectrophotometry at 260 and 280 nm was performed to measure the amount and quality of RNA. The total RNA samples were first treated with DNase I to degrade any possible deoxyribonucleic acid (DNA) contamination. Then, the messenger RNA (mRNA) was enriched by using the oligo (dT) magnetic beads (for eukaryotes). Mixed with the fragmentation buffer, the mRNA was fragmented into short fragments (about 200 base pair). Then the first strand of complimentary DNA (cDNA) was synthesized by using random hexamer-primer. Buffer, dNTPs, RNase H, and DNA polymerase I were added to synthesize the second strand. The double strand cDNA was purified with magnetic beads. End reparation and 30 -end single nucleotide A (adenine) addition were then performed. Finally, sequencing adaptors were ligated to the fragments. The fragments were enriched by polymerase chain reaction (PCR) amplification. During the quality control step, Agilent 2100 Bioanaylzer (Agilent, Santa Clara, CA, USA) and ABI StepOnePlus Real-Time PCR System (ABI, Foster City, CA 94404, USA) was used to qualify and quantify the sample library. The library products were ready for sequencing via Illumina HiSeqTM 2000 (Illumina, San Diego, CA, USA). 2.4. Reverse Transcription Polymerase Chain Reaction The extraction and quantification of mRNA by means of reverse transcription polymerase chain reaction (RT-PCR) were performed as our standard protocols [32]. The RNA was then converted to complementary DNA by reverse transcription with random hexanucleotides and avian myeloblastosis virus reverse transcriptase (Boehringer, Mannheim, Germany). The primers are as follows: CA4: forward: 50 -CACTGGTGCTACGAGGTTCA-30 , reverse: 50 -GCAACTGTTTGGCCTGGTAT-30 ; CHI3L1: forward: 50 -AAGAACAGGAACCCCAACCT-30 , reverse: 50 -TGTCTCTCCGTCCAGGGTAG-30 ; FCN1: forward: 50 -GCTCGCTGTCCTGCTAGTCT-30 , reverse: 50 -CTCCAATGACACCTGCCTCT-30 ; HAS2: forward: 50 -GCCTCATCTGTGGAGATGGT-30 , reverse: 50 -ATGCACTGAACACACCCAAA-30 . KLK5: forward: 50 -AGTCAGAAAAGGTGCGAGGA-30 , reverse: 50 -TAATCTCCCCAGGACACGAG-30 . 2.5. Bioinformatics Analysis of RNA Sequence Data High quality reads with less than 10% unknown bases and median quality scores above 20 were extracted from the raw RNA-seq data. Cleaned reads were mapped to the reference genome hg19 using SOAPaligner/SOAP2 (Beijing Genomics Institute, Beijing, China) [33] with no more than two mismatches allowed in the alignment. The expression level for each gene was calculated using the RPKM method [34]. If there was more than one transcript for a gene, the longest transcript was used to calculate expression level and coverage. Function annotation could provide gene ontology (GO) annotation information.

J. Clin. Med. 2018, 7, 537

4 of 13

2.6. Screening of Differentially Expressed Genes (DEGs) This analysis included the screening of genes that were differentially expressed among samples, GO functional enrichment analysis and KEGG pathway enrichment analysis for these DEGs. We used false discovery rate (FDR) ≤0.001 and the absolute value of log2 ratio ≥1 as the threshold to judge the significance of gene expression difference. 2.7. Pathway and Network Analysis While so many components are implicated in the mechanism of a complex trait disease, it is difficult to infer which components play a determining role. In this scenario, pathway analysis is a tool used to identify possible biological pathways involving related components. Ingenuity pathway analysis (IPA) is a powerful analysis and pathway-search tool that uncovers the significance of omics data and identifies new targets or candidate biomarkers within the context of biological systems. IPA has been specifically designed for the analysis, integration, and interpretation of data derived from omics experiments. Therefore, in the present study, to derive a possible pathway that may be implicated in the mechanism of AF related thromboembolic stroke, data were analyzed through the use of IPA (QIAGEN Inc., Germantown, MD, USA) [35]. The significance values for the canonical pathways in IPA was calculated using the right-tailed Fisher’s exact test. The p-value for a given pathway annotation is calculated by considering (1) the number of DEGs that participate in that pathway and (2) the total number of genes that are known to be associated with that pathway in the selected reference set. Higher the number of DEGs involved, the more likely it is that the association is not due to random chance, and thus more significant the p-value is. Similarly, the larger the total number of genes known to be associated with the pathway, the greater the likelihood that an association is due to random chance, and the p-value accordingly becomes less significant. The p-value identifies statistically significant over-representation of DEGs in a given pathway. The null hypothesis (H0 ) states. H0 = overlap of molecules for a particular biological category is due to chance. If p-value < 0.005, we reject the null hypothesis. This implies the association between a set of DEGs in our experiment and a given pathway is not due to random chance. 2.8. Zebrafish Model Tg(fli1:egfp)xTg(gata1:dsRed) (egfp, enhanced green fluorescent protein), zebrafish with green fluorescence emitting vasculature and red fluorescence emitting red blood cells were obtained from the Zebrafish International Resource Center (Eugene, OR, USA) [36]. Both wild-type (WT; AB strain) and Tg(fli1:egfp)xTg(gata1:dsRed) zebrafish were kept under a 14-h light and 10-h dark photoperiod at 28.5 ◦ C. The procedures for zebrafish culture, embryo collection, and fluorescence observation have been described previously [37]. We followed the standard criteria to design the zebrafish developmental stages [38]. The zebrafish experiments were approved by the Institutional Animal Care and Use Committee of the National Taiwan University College of Medicine. The anesthesia and euthanasia were performed with diluted tricane methanesulfonate (160 mg/L) and benzocaine hydrochloride (250 mg/L). 2.9. Antisense Morpholino Design and Microinjection Morpholino oligonucleotides (MO) complementary to the 50 - translational start codon flanking sequences of zebrafish targeted gene mRNAs was designed and purchased (Gene Tools LLC, Philomath, OR, USA). Zebrafish has2-morpholino (has2-MO: 50 -GCTGACCGCTTTATCACATCTCATC-30 ) was established based on zebrafish hyaluronan synthase 2 (has2) cDNA sequences (Gene Tools LLC, Philomath, OR, USA). It was prepared at stocking concentrations of 1 mM and diluted with double-distilled water to concentrations of 100, 250, 375, and 500 µM (4.5 ng per 2.3 nL). Morpholino oligonucleotide stocks (3 mmol/L) were diluted to 0.3 to 1.2 mmol/L in Danieau solution (58 mmol/L

J. Clin. Med. 2018, 7, 537

5 of 13

NaCl, 0.7 mmol/L KCl, 0.4 mmol/L MgSO4, 0.6 mmol/L Ca(NO3 )2 , 5 mmol/L HEPES, pH 7.6) before injection into the yolks of freshly fertilized fish eggs. Approximately 1 nL of morpholino oligonucleotides were injected into 40 to 180 embryos at the one- to four-cell stage, and the embryos were allowed to develop at 29 ◦ C. 2.10. Statistical Methods Continuous data were presented as means ± standard deviation (SD) and analyzed by using the Mann-Whitney U test for two-group comparisons. Categorical data were compared with the Fisher’s exact test. 3. Results 3.1. Global Expressional Profiling of VV Global expressional profiling of venous samples from patients with VV and normal venous samples using the RNA seq technique was performed. The expressional profiling of the significantly up or downregulated genes with log2 ratio >7 were listed in Table 2. Among all these differentially expressed genes, many of them encode proteins functioning as extracellular matrix enzymes, such as Kallikrein-related peptidase 5 (KLK5), chitinase-3-like protein 1 (CHI3L1), hyaluronan synthase 2 (HAS2), and carbonic anhydrase 4 (CA4). These proteins may be involved in remodeling of extracellular matrix proteins and play important roles in vascular integrity and scaffolding. The others are genes encoding cytokines such as chemokine (C-X-C motif) ligand 1 (CXCL1), chemokine (C-X-C motif) ligand 8 (CXCL8), colony stimulating factor 3 (CSF3), and interleukin 6 (IL-6), implicating the importance of paracrine factors related to cell growth and inflammation in the mechanism of VV. Table 2. Highly differentially expressed gene list between venous samples from control (CTL) patients and patients with varicose veins (VV). Gene CRABP1 SNORA77 DPEP1 PSORS1C1 CST2 HLA-DQA2 CLIC3 WFDC10B KLK5 CPNE7 CHAD WFDC3 S100A8 GSTT1 CXCL1 DAB1 CHI3L1 CSF3 CXCL8 FCGR3A/FCGR3B S100A12 CD33 FCN1 PQLC2L AQP9 SELE HAS2 CA4 IL6 THEM5 ZNF385B SRGAP2B FCGR2C IL1R2 ALAS2 HHATL

Log2 Ratio

Regulation

(VV/CTL)

(VV/CTL)

11.98 10.16 9.69 9.56 9.32 8.68 8.46 8.20 7.77 7.51 7.39 7.16 −12.33 −11.46 −10.01 −9.99 −9.98 −9.96 −9.56 −9.28 −8.95 −8.70 −8.53 −8.10 −7.93 −7.85 −7.69 −7.63 −7.62 −7.58 −7.52 −7.39 −7.18 −7.13 −7.08 −7.08

Up Up Up Up Up Up Up Up Up Up Up Up Down Down Down Down Down Down Down Down Down Down Down Down Down Down Down Down Down Down Down Down Down Down Down Down

Location

Type(s)

Cytoplasm Other Cytoplasm Other Extracellular Space Plasma Membrane Nucleus Extracellular Space Extracellular Space Cytoplasm Extracellular Space Extracellular Space Cytoplasm Cytoplasm Extracellular Space Cytoplasm Extracellular Space Extracellular Space Extracellular Space Plasma Membrane Cytoplasm Plasma Membrane Extracellular Space Other Plasma Membrane Plasma Membrane extracellular matrix Plasma Membrane Extracellular Space Cytoplasm Nucleus Other Plasma Membrane Plasma Membrane Cytoplasm Cytoplasm

transporter other peptidase other other transmembrane receptor ion channel other peptidase transporter other other other enzyme cytokine other enzyme cytokine cytokine transmembrane receptor other other other other transporter transmembrane receptor enzyme enzyme cytokine enzyme other other transmembrane receptor transmembrane receptor enzyme enzyme

J. Clin. Med. 2018, 7, 537

6 of 13

Because these genes might be implicated in the mechanism of VV, we did a pathway analysis to elucidate the possible translational meaning and significance (Table 3; Figure 1). Interestingly, the results show that decreased angiogenesis, decreased cell-to-cell binding, decreased cell invasion or motility, and enhanced vascular lesion pathways are involved in the mechanism of VV. The results also show that increased carbohydrate accumulation may be involved in the mechanism of VV. Among all these pathways, HAS2 gene downregulation plays a pivotal role to regulate and govern all these pathways (Figure 1). Table 3. Pathways involved in the transcriptional regulation of varicose veins Diseases or Functions Annotation

p-Value *

Predicted Activation State

Angiogenesis

8.48 × 10−5

Decreased

Binding of cells

3.07 × 10−4

Decreased

Invasion of cells

3.43 ×

10−4

Decreased

Vascular lesion Synthesis of carbohydrate Accumulation of carbohydrate

7.74 × 10−4

CHI3L1, CSF3, CXCL1, CXCL8, HAS2, IL6, S100A12, S100A8, SELE CXCL1, CXCL8, HAS2, IL1R2, IL6, SELE CHI3L1, CSF3, CXCL1, CXCL8, HAS2, IL6, S100A12, S100A8 HAS2, IL6, S100A12, SELE

1.32 × 10−3

AQP9, CXCL8, HAS2, IL1R2, IL6

1.93 × 10−3

CXCL1, CXCL8, HAS2

Molecules

* p-value was calculated using the right-tailed Fisher’s exact test. J. Clin. Med. 2018, 7, x FOR PEER REVIEW

7 of 14

Figure 1.Figure Pathway or network analysis derived thedifferentially differentially expressed identified by 1. Pathway or network analysis derivedfrom from the expressed genesgenes identified by next-generation RNA sequence technology (RNA seq) between venous samples of control patients and next-generation RNA sequence technology (RNA seq) between venous samples of control patients those of patients with varicose N = each 5 for each group. those of and patients with varicose vein. Nvein. = 5 for group.

We then validated several differentially expressed genes identified in the RNA seq by conventional RT-PCR using VV and control venous samples. Among all the up or downregulated genes in VV samples, there were several genes encoding enzymes in the extracellular matrix such as KLK5, CHI3L1, HAS2, and CA4. We have validated the expressional changes by conventional RT-PCR for these genes (Figure 2). The result showed decreased mRNA levels of CHI3L1, HAS2, and CA4 genes and increased mRNA level of KLK5 gene in VV samples, which was similar to the result of RNA seq. We also found

J. Clin. Med. 2018, 7, 537 Figure 1. Pathway or network analysis derived from the differentially expressed genes identified by next-generation RNA sequence technology (RNA seq) between venous samples of control patients geneand Ficolin encoding space protein the mRNA level of which those 1 of(FCN1) patients with varicoseextracellular vein. N = 5 for each group.

another decreased in the VV samples.

7 of 13

was also

Figure 2. Validation of the differentially expressed genes identified by next-generation RNA sequence Figure 2. Validation of thereverse differentially expressed genes identifiedchain by next-generation RNA sequence technology (RNA seq) using transcription polymerase reaction (RT-PCR) in venous technology seq)(CL) using reverse transcription chainvein reaction (RT-PCR) in venous samples of control (RNA patients and those of patientspolymerase with varicose (VV). N = 5 for each group. samples of control patients (CL) and those of patients with varicose vein (VV). N = 5 for each group. CA4, carbonic anhydrase 4; CHI3L1, chitinase-3-like protein 1; FCN1, Ficolin-1; HAS2, hyaluronan CA4, carbonic anhydrase 4; CHI3L1, chitinase-3-like protein 1; FCN1, Ficolin-1; HAS2, hyaluronan synthase 2; KLK5, kallikrein-related peptidase. synthase 2; KLK5, kallikrein-related peptidase.

3.2. Downregulation of Hyaluronan Synthases 2 in Venous Tissues from Patients with VV 3.2. Downregulation of Hyaluronan Synthases 2 in Venous Tissues from Patients with VV

Among Among all these extracellular matrix genes, we were particularly all these extracellular matrix enzyme enzyme genes, we were particularly interestedinterested in the HAS2 in the J. Clin. Med. 2018, 7, x FOR PEER REVIEW 8 of 14 HAS2 gene played a pivotal in the pathway analysis, which modulated genebecause because ititplayed a pivotal role inrole the pathway analysis, which modulated all the VV-relatedall the VV-related pathways (Figure 1). of The role of HASbiology in vascular biology has never been addressed pathways (Figure 1). The role HAS in vascular has never been addressed before. Because before. Because of its important possible important role in the of mechanism of VV,the weexperiment repeated the experiment of its possible role in the mechanism VV, we repeated again and confirmed that the mRNA level waslevel significantly lowered in thelowered venous samples patients with from again and confirmed that the mRNA was significantly in the from venous samples VV (p VV < 0.0001) patients with (p < (Figure 0.0001)3). (Figure 3).

Figure 3.Figure Expression of theofhyaluronan synthases (HAS2)ininvenous venous samples 3. Expression the hyaluronan synthases22 gene gene (HAS2) samples from from controlcontrol patients (CL) and patients varicose vein (VV).(A) (A)Expression Expression ofof HAS2 gene evaluated by reverse patients (CL) and patients with with varicose vein (VV). HAS2 gene evaluated by reverse transcription polymerase reaction (RT-PCR); Dot plotofofquantification quantification ofofRT-PCR product transcription polymerase chain chain reaction (RT-PCR); (B)(B) Dot plot RT-PCR product by by densitometry. N = 5 for each group. * p < 0.05. densitometry. N = 5 for each group. * p < 0.05. 3.3. Knockdown of HAS2 Results in Venous Dilation and Blood Flow Stasis in Zebrafish Tg(fli1:egfp)xTg(gata1:dsRed) zebrafish with green fluorescence emitting vasculature and red fluorescence emitting red blood cells were used to observe the anatomic change of vasculature (green fluorescence) and blood flow (red fluorescence) and to explore the functional role of the identified gene. We first compared the homology of amino acid sequence of HAS2 among species using Clustal W [39]. HAS2 is evolutionally conservative and expressed in both low and high vertebrates, indicating its important biological function (Figure S1). The similarity of amino acid sequence of

J. Clin. Med. 2018, 7, 537

8 of 13

3.3. Knockdown of HAS2 Results in Venous Dilation and Blood Flow Stasis in Zebrafish Tg(fli1:egfp)xTg(gata1:dsRed) zebrafish with green fluorescence emitting vasculature and red fluorescence emitting red blood cells were used to observe the anatomic change of vasculature (green fluorescence) and blood flow (red fluorescence) and to explore the functional role of the identified J.gene. Clin. Med. x FOR PEER REVIEW 9 of 14 We 2018, first 7, compared the homology of amino acid sequence of HAS2 among species using Clustal W [39]. HAS2 is evolutionally conservative and expressed in both low and high vertebrates, indicating 1). In those with severefunction phenotypes, tangled veins were found in the tail which resulted in its important biological (Figure S1). The similarity of amino aciddistal sequence of HAS2 between malformation tail (Figure 4E,F; Supplementary movie 3). zebrafish and of humans was up to 81.3% (Figure S2).

Figure 4. Knockdown of hyaluronan synthases 2 gene (HAS2) in zebrafish embryos leads to defective Figure 4. Knockdown of hyaluronan synthases 2 gene (HAS2) in zebrafish embryos leads to defective venous structure with blood flow stasis. Tg(fli1:egfp)xTg(gata1:dsRed) zebrafish with green fluorescence venous structure with blood flow stasis. Tg(fli1:egfp)xTg(gata1:dsRed) zebrafish with green emitting vasculature and red fluorescence emitting red blood cells are used to observe the anatomic fluorescence emitting vasculature and red fluorescence emitting red blood cells are used to observe change of vasculature (green fluorescence) and blood flow (red fluorescence). (A) Morphology of the the anatomic change of vasculature (green fluorescence) and blood flow (red fluorescence). (A) tail of a representative wild-type zebrafish (48 hpf) under bright field microscope; (B) Vascular structure Morphology of the tail of a representative wild-type zebrafish (48 hpf) under bright field microscope; of the same zebrafish tail under fluorescence microscope; (C) Morphology of the tail of a representative (B) Vascular structure of the same zebrafish tail under fluorescence microscope; (C) Morphology of HAS2-morpholino (MO) injected zebrafish (48 hpf) under bright field microscope. A protruding bulb the tail of a representative HAS2-morpholino (MO) injected zebrafish (48 hpf) under bright field in the tail is noted (black arrow), mimicking human phenotype of varicose vein; (D) Vascular structure microscope. A protruding bulb in the tail is noted (black arrow), mimicking human phenotype of of the same HAS2-MO injected zebrafish tail under fluorescence microscope. A corresponding local varicose vein; (D) Vascular structure of the same HAS2-MO injected zebrafish tail under fluorescence dilation of venous structure in the caudal vein plexus (CVP) is found (white arrow); (E) Morphology microscope. A corresponding local dilation of venous structure in the caudal vein plexus (CVP) is of a representative HAS2-morpholino (MO) injected zebrafish (48 hpf) with tail malformation under found (white arrow); (E) Morphology of a representative HAS2-morpholino (MO) injected zebrafish bright field microscope; (F) Vascular structure of the same HAS2-MO injected zebrafish tail under (48 hpf) with tail malformation under bright field microscope; (F) Vascular structure of the same fluorescence microscope. Tangled veins are found in the mal-formatted tail (white stars). Magnification HAS2-MO injected zebrafish tail under fluorescence microscope. Tangled veins are found in the malis 100× in (A–D) and 200× in (E,F). Scale bar = 100 µm. Hpf indicates hours post-fertilization. formatted tail (white stars). Magnification is 100× in (A–D) and 200× in (E,F). Scale bar = 100 μm. Hpf indicates hours post-fertilization.

J. Clin. Med. 2018, 7, 537

9 of 13

Morpholino experiment was quantitatively performed to further define the role of HAS2 in the mechanism of VV. In fish with mild phenotypes, a local dilation of venous structure in the caudal vein plexus (CVP) was found, which resulted in a protruding bulb in the tail, mimicking human VV phenotype of protruding venous dilatation (Figure 4A–D). In fish with moderate phenotypes, dilation of veins, such as intersegmental vein (ISV) or sub-intestinal vein (SIV), was found (Figure 5). Sluggish blood flow in CVP was also noted (Supplementary Movie 1 and Movie 2), mimicking the human VV phenotype of venous blood flow stasis and leg edema. Furthermore, venous branches were also noted in the SIV of the MO-injected fish, which was a consequence of either a dysregulated angiogenetic process or a delay of retraction of leading venous buds during development [40]. It implies that this is a dysregulated angiogenetic process since HAS2 plays a pivotal role in the angiogenesis network, which is one of the key altered pathways in our transcriptome analysis (Figure 1). In those with severe phenotypes, tangled veins were found in the distal tail which resulted in malformation of tail (Figure 4E,F; Movie 3). J. Clin.Supplementary Med. 2018, 7, x FOR PEER REVIEW 10 of 14

Figure 5. Knockdown of hyaluronan synthases22 gene gene (HAS2) in zebrafish embryos leads toleads dilatedto dilated Figure 5. Knockdown of hyaluronan synthases (HAS2) in zebrafish embryos venous structure. Tg(fli1:egfp)xTg(gata1:dsRed)zebrafish zebrafish with green fluorescence emitting vasculature venous structure. Tg(fli1:egfp)xTg(gata1:dsRed) with green fluorescence emitting vasculature are used to observe the anatomic change of vasculature. (A) Inter-segmental veins (ISVs) in the tail of are used to observe the anatomic change of vasculature. (A) Inter-segmental veins (ISVs) in the tail of a a representative wild-type zebrafish (48 hpf) under fluorescence microscope (arrows); (B) Dilated representative wild-type (48 hpf) under fluorescence microscope (arrows); Dilated ISVs ISVs in the tail ofzebrafish a representative HAS2-morpholino (MO) injected zebrafish (48 hpf)(B) under in the tail fluorescence of a representative (MO) injected hpf) underwildfluorescence microscopeHAS2-morpholino (arrows); (C) The sub-intestinal vein (SIV)zebrafish in the tail of(48 a representative type(arrows); zebrafish (48 under fluorescence microscope (arrow); (D)tail The of dilated SIV (arrow) withwild-type microscope (C)hpf) The sub-intestinal vein (SIV) in the a representative protruding branches (arrowheads) in the tail of a representative HAS2-morpholino (MO) injected zebrafish (48 hpf) under fluorescence microscope (arrow); (D) The dilated SIV (arrow) with protruding zebrafish (48 hpf) under fluorescence microscope. Hpf indicates hours post-fertilization. branches Magnification (arrowheads) in the tail of a representative HAS2-morpholino (MO) injected zebrafish is 400×. Scale bar = 100 μm. (48 hpf) under fluorescence microscope. Hpf indicates hours post-fertilization. Magnification is Quantitative 400×. Scale bar = 100 data µm. after the injection of HAS2 MO in more than 1000 embryos are shown in Table

4. There was a dose-response relationship. Higher dose of MO (375 uM) was associated with a higher percentagedata of fish withthe defective vascular phenotypes lower dose MO embryos (250 uM). In theshown highestin Table 4. Quantitative after injection of HAS2 MO inthan more than 1000 are dose group (500 uM), because of the lethality of high-dose MO, most of the embryos were dead (91%). There was a dose-response relationship. Higher dose of MO (375 uM) was associated with a higher The survived embryos were probably those without or with mild MO effect and thus there was a percentage of fishlower withpercentage defectiveofvascular phenotypes lower(43%), dosecompared MO (250 In the relatively survived fish with defect than vasculature to uM). those with a highest dose group (500 uM), because of the(72%) lethality of high-dose MO, most the embryos were dead (91%). medium dose of MO injection (375 uM) (Table 4). Therefore, weof confirm that disruption of HAS2 function directly leads to vascular malformation in vertebrates.

J. Clin. Med. 2018, 7, 537

10 of 13

The survived embryos were probably those without or with mild MO effect and thus there was a relatively lower percentage of survived fish with defect vasculature (43%), compared to those with a medium dose of MO injection (72%) (375 uM) (Table 4). Therefore, we confirm that disruption of HAS2 function directly leads to vascular malformation in vertebrates. Table 4. Morphological phenotypes of zebrafish embryos (fli1/gata1) derived from fertilized eggs injected with HAS2-MO. Injection Dose (µM)

Injected Embryos

250 375 500

231 563 502

Survival Embryos 1 (Survival Rates %) 2 145 (62%) 254 (45%) 49 (9%)

GFP-Positive Embryos

Embryos without Defective Vascular Phenotypes (%) *

Embryos with Defective Vascular Phenotypes (%) *

132 240 28

84 (63%) 68 (28%) 16 (57%)

48 (37%) 172 (72%) 12 (43%)

1

Total numbers of the survival embryos at 48 hpf. 2 Survival rates indicate the numbers of survival embryos among the numbers of injected embryos at 48 hpf. * Percentages indicate the numbers of no defect or vascular defective embryos among the numbers of GFP-positive embryos at 48 hpf. hpf, hours post-fertilization. GFP, green fluorescent protein.

4. Discussion In the present study, we first used a global expressional profiling to address the molecular mechanism of VV. We further found several important pathways, such as angiogenesis, cell binding, vascular lesion and carbohydrate metabolism possibly involved in the mechanism of VV. We further used a zebrafish model to address the functional role of an identified gene HAS2 in the mechanism of VV. In the last decade, the pathophysiological theories of VV shifted from purely mechanical factors to hypotheses focusing on complex molecular and histopathologic alterations in the vessel wall and the extracellular matrix, which may be caused by the genetic inheritance [1,18,19]. In our study, we identified the global expressional change of downregulations of CHI3L1, HAS2, and CA4, and upregulation of KLK5 in the varicose venous samples with both RNA seq technique and conventional RT-PCR. All these genes encode extracellular matrix enzymes, compatible with the concept of histopathologic alterations in the extracellular matrix as a potential mechanism of VV. Extracellular matrix may provide a scaffold for the vascular structure and integrity. When the components of extracellular matrix have been altered (e.g., decomposition of extracellular matrix proteins by excessive enzyme digestion or accumulation of extracellular matrix proteins by decreased enzyme digestion), the scaffold function may deteriorate, resulting in a weakening of vascular wall support and dilatation or protrusion of the vasculature, a pathognomonic picture of VV. Furthermore, after pathway analysis using all the differentially expressed genes, we found decreased angiogenesis, decreased cellular binding and invasion, increased carbohydrate accumulation, and vascular lesion were all involved in the mechanism of VV. All these pathways are novel and have been rarely addressed as the potential mechanism of VV. Decreased angiogenesis may implicate impaired repair or proliferative activities of vessels and an injury model of VV. Furthermore, an involvement of vascular lesion network also supports the hypothesis of this injury model of VV. Decreased cell-cell binding and invasion may implicate decreased structural integrity of vessels, indicating weakening of vessel wall and possible protrusion of vessel wall, which is the prototype phenotype of VV. Our study also shows that increased carbohydrate accumulation may be involved in the pathogenesis of VV; however, further studies are required to elucidate the detailed mechanism. As mentioned before, recent studies have also highlighted the pathological changes in the extracellular matrix as the potential mechanism of VV [1,8,10,12,13,15]. Among all these extracellular matrix enzyme genes, HAS2 expression was decreased in VV samples. Hyaluronan or hyaluronic

J. Clin. Med. 2018, 7, 537

11 of 13

acid (HA) is a high molecular weight unbranched polysaccharide, and is a major constituent of the extracellular matrix and synthesized by hyaluronan synthases (HAS), which are membrane-bound enzymes and use UDP-α-N-acetyl-D-glucosamine and UDP-α-D-glucuronate as substrates to produce the glycosaminoglycan hyaluronan at the cell surface and extrude it through the membrane into the extracellular space. In our zebrafish model, we found disruption of HAS2 function directly led to venous malformation in vertebrates. Although HAS2 has been shown to play a role in vascular morphogenesis [41], its specific role in venous tissues and the pathophysiology of the most important human venous diseases—that is, VV—has never been addressed before. The mechanism by which HAS2 cause VV may be through its pivotal role to modulate and govern multiple pathways that cause VV, such as angiogenesis, cell to cell adhesion, vascular lesion, and carbohydrate metabolism. There are limitations in the present study. First, there were only five patients in each group and the differences of genetic results we found between VV patients and non-VV patients might be due to something else not related to VV. Second, the genetic changes in the VV samples might be post-transcriptional or happened due to the local environment in the legs, and the zebrafish model only showed a systemic effect. Therefore, HAS2 might not be the gene responsible for VV. Third, we did not have negative scrambled MO control group, positive control group with other MOs targeting to HAS2 and rescue studies. Therefore, we could not rule out a non-specific effect of MO. Finally, some of our control patients had diabetes that could modify leukocyte activity and might affect the genetic results. However, there is little doubt that the development and progression of VV and venous ulcer disease is multifactorial and includes a moderate-to-strong genetic component [8]. The genetic studies of VV conducted to date are few and have few functional implications. This is the first study to investigate the global gene expressional profiling in determining the candidate gene of VV. In this study, we show the possible role of HAS2 gene in the pathogenesis of VV in both global expression and functional study. Supplementary Materials: The Figure S1, Figure S2 and Supplementary Movie 1–3 are available online at http://www.mdpi.com/2077-0383/7/12/537/s1. Author Contributions: Conceptualization, C.-S.H., C.-T.T., E.Y.C., and I-H.W.; Methodology, C.-S.H., C.-T.T., and Y.-H.C.; Software, C.-S.H. and S.-N.C.; Validation, E.Y.C. and I-H.W.; Formal Analysis, C.-S.H. and I-H.W.; Investigation, C.-S.H., C.-T.T., and I-H.W.; Resources, J.-J.H., E.Y.C., and I-H.W.; Data Curation, C.-S.H. and Y.-H.C.; Writing—Original Draft Preparation, C.-S.H. and C.-T.T.; Writing—Review & Editing, C.-S.H., C.-T.T, E.Y.C. and I-H.W.; Visualization, E.Y.C. and I-H.W.; Supervision, E.Y.C. and I-H.W.; Project Administration, E.Y.C. and I-H.W.; Funding Acquisition, E.Y.C. and I-H.W. Funding: The work was supported by two grants from the Ministry of Science and Technology (MOST) in Taiwan (MOST104-2314-B-002-202 and 105-2628-B-002-050-MY2). We thank all the participating patients for their contribution to this study. We also thank the second Core Lab of the National Taiwan University Hospital for technological assistance and laboratory facility support. Conflicts of Interest: The authors declare no conflict of interest.

References 1. 2.

3. 4.

Oklu, R.; Habito, R.; Mayr, M.; Deipolyi, A.R.; Albadawi, H.; Hesketh, R.; Walker, T.G.; Linskey, K.R.; Long, C.A.; Wicky, S.; et al. Pathogenesis of varicose veins. J. Vasc. Interv. Radiol. 2012, 23, 33–39. [CrossRef] Li, X.; Jiang, X.Y.; Ge, J.; Wang, J.; Chen, G.J.; Xu, L.; Xie, D.Y.; Yuan, T.Y.; Zhang, D.S.; Zhang, H.; et al. Aberrantly expressed lncRNAs in primary varicose great saphenous veins. PLoS ONE 2014, 9, e86156. [CrossRef] Meissner, M.H.; Gloviczki, P.; Bergan, J.; Kistner, R.L.; Morrison, N.; Pannier, F.; Pappas, P.J.; Rabe, E.; Raju, S.; Villavicencio, J.L. Primary chronic venous disorders. J. Vasc. Surg. 2007, 46 (Suppl. S), 54s–67s. [CrossRef] Heit, J.A.; Rooke, T.W.; Silverstein, M.D.; Mohr, D.N.; Lohse, C.M.; Petterson, T.M.; O’Fallon, W.M.; Melton, L.J., 3rd. Trends in the incidence of venous stasis syndrome and venous ulcer: A 25-year population-based study. J. Vasc. Surg. 2001, 33, 1022–1027. [CrossRef]

J. Clin. Med. 2018, 7, 537

5. 6. 7. 8. 9.

10. 11. 12. 13. 14. 15. 16. 17.

18. 19. 20.

21.

22. 23.

24.

25. 26.

12 of 13

Brand, F.N.; Dannenberg, A.L.; Abbott, R.D.; Kannel, W.B. The epidemiology of varicose veins: The Framingham Study. Am. J. Prev. Med. 1988, 4, 96–101. [CrossRef] Bergan, J.J.; Schmid-Schonbein, G.W.; Smith, P.D.; Nicolaides, A.N.; Boisseau, M.R.; Eklof, B. Chronic venous disease. N. Engl. J. Med. 2006, 355, 488–498. [CrossRef] [PubMed] Beebe-Dimmer, J.L.; Pfeifer, J.R.; Engle, J.S.; Schottenfeld, D. The epidemiology of chronic venous insufficiency and varicose veins. Ann. Epidemiol. 2005, 15, 175–184. [CrossRef] [PubMed] Krysa, J.; Jones, G.T.; van Rij, A.M. Evidence for a genetic role in varicose veins and chronic venous insufficiency. Phlebology 2012, 27, 329–335. [CrossRef] Naik, B.; Kumar, M.; Khanna, A.K.; Suman, P.K. Clinico-histopathological study of varicose vein and role of matrix metalloproteinases-1, matrix metalloproteinases-9 and tissue inhibitor of matrix metalloproteinase-1 in varicose vein formation. Indian J. Pathol. Microbiol. 2016, 59, 25–30. [CrossRef] [PubMed] Chen, Y.; Peng, W.; Raffetto, J.D.; Khalil, R.A. Matrix Metalloproteinases in Remodeling of Lower Extremity Veins and Chronic Venous Disease. Prog. Mol. Biol. Transl. Sci. 2017, 147, 267–299. [CrossRef] Xu, Y.; Bei, Y.; Li, Y.; Chu, H. Phenotypic and functional transformation in smooth muscle cells derived from varicose veins. J. Vasc. Surg. Venous Lymphat. Disord. 2017, 5, 723–733. [CrossRef] [PubMed] Xu, H.M.; Zhao, Y.; Zhang, X.M.; Zhu, T.; Fu, W.G. Polymorphisms in MMP-9 and TIMP-2 in Chinese patients with varicose veins. J. Surg. Res. 2011, 168, e143–e148. [CrossRef] [PubMed] Jacobs, B.N.; Andraska, E.A.; Obi, A.T.; Wakefield, T.W. Pathophysiology of varicose veins. J. Vasc. Surg. Venous Lymphat. Disord. 2017, 5, 460–467. [CrossRef] [PubMed] Raffetto, J.D.; Mannello, F. Pathophysiology of chronic venous disease. Int. Angiol. J. Int. Union Angiol. 2014, 33, 212–221. Bharath, V.; Kahn, S.R.; Lazo-Langner, A. Genetic polymorphisms of vein wall remodeling in chronic venous disease: A narrative and systematic review. Blood 2014, 124, 1242–1250. [CrossRef] [PubMed] Ng, M.Y.; Andrew, T.; Spector, T.D.; Jeffery, S. Linkage to the FOXC2 region of chromosome 16 for varicose veins in otherwise healthy, unselected sibling pairs. J. Med. Genet. 2005, 42, 235–239. [CrossRef] [PubMed] Mellor, R.H.; Brice, G.; Stanton, A.W.; French, J.; Smith, A.; Jeffery, S.; Levick, J.R.; Burnand, K.G.; Mortimer, P.S. Mutations in FOXC2 are strongly associated with primary valve failure in veins of the lower limb. Circulation 2007, 115, 1912–1920. [CrossRef] Matousek, V.; Prerovsky, I. A contribution to the problem of the inheritance of primary varicose veins. Hum. Hered. 1974, 24, 225–235. Hauge, M.; Gundersen, J. Genetics of varicose veins of the lower extremities. Hum. Hered. 1969, 19, 573–580. [CrossRef] Brice, G.; Mansour, S.; Bell, R.; Collin, J.R.; Child, A.H.; Brady, A.F.; Sarfarazi, M.; Burnand, K.G.; Jeffery, S.; Mortimer, P.; et al. Analysis of the phenotypic abnormalities in lymphoedema-distichiasis syndrome in 74 patients with FOXC2 mutations or linkage to 16q24. J. Med. Genet. 2002, 39, 478–483. [CrossRef] Jin, Y.; Xu, G.; Huang, J.; Zhou, D.; Huang, X.; Shen, L. Analysis of the association between an insertion/deletion polymorphism within the 3’ untranslated region of COL1A2 and chronic venous insufficiency. Ann. Vasc. Surg. 2013, 27, 959–963. [CrossRef] [PubMed] Kowalewski, R.; Malkowski, A.; Gacko, M.; Sobolewski, K. Influence of thrombophlebitis on TGF-beta1 and its signaling pathway in the vein wall. Folia Histochem. Cytobiol. 2010, 48, 542–548. [CrossRef] [PubMed] Zhong, H.; Yang, X.; Kaplan, L.M.; Molony, C.; Schadt, E.E. Integrating pathway analysis and genetics of gene expression for genome-wide association studies. Am. J. Hum. Genet. 2010, 86, 581–591. [CrossRef] [PubMed] Schunkert, H.; Konig, I.R.; Kathiresan, S.; Reilly, M.P.; Assimes, T.L.; Holm, H.; Preuss, M.; Stewart, A.F.; Barbalic, M.; Gieger, C.; et al. Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease. Nat. Genet. 2011, 43, 333–338. [CrossRef] [PubMed] Lotta, L.A. Genome-wide association studies in atherothrombosis. Eur. J. Intern. Med. 2010, 21, 74–78. [CrossRef] Hardy, J.; Singleton, A. Genomewide association studies and human disease. N. Engl. J. Med. 2009, 360, 1759–1768. [CrossRef]

J. Clin. Med. 2018, 7, 537

27.

28.

29.

30.

31.

32.

33. 34. 35. 36. 37. 38. 39.

40.

41.

13 of 13

Shadrina, A.S.; Smetanina, M.A.; Sokolova, E.A.; Shamovskaya, D.V.; Sevost’ianova, K.S.; Shevela, A.I.; Soldatsky, E.Y.; Seliverstov, E.I.; Demekhova, M.Y.; Shonov, O.A.; et al. Allele rs2010963 C of the VEGFA gene is associated with the decreased risk of primary varicose veins in ethnic Russians. Phlebology 2018, 33, 27–35. [CrossRef] Kundu, S.; Lurie, F.; Millward, S.F.; Padberg, F., Jr.; Vedantham, S.; Elias, S.; Khilnani, N.M.; Marston, W.; Cardella, J.F.; Meissner, M.H.; et al. Recommended reporting standards for endovenous ablation for the treatment of venous insufficiency: Joint statement of The American Venous Forum and The Society of Interventional Radiology. J. Vasc. Surg. 2007, 46, 582–589. [CrossRef] Brittenden, J.; Cotton, S.C.; Elders, A.; Ramsay, C.R.; Norrie, J.; Burr, J.; Campbell, B.; Bachoo, P.; Chetter, I.; Gough, M.; et al. A Randomized Trial Comparing Treatments for Varicose Veins. N. Engl. J. Med. 2014, 371, 1218–1227. [CrossRef] Sufian, S.; Arnez, A.; Labropoulos, N.; Nguyen, K.; Satwah, V.; Marquez, J.; Chowla, A.; Lakhanpal, S. Radiofrequency ablation of the great saphenous vein, comparing one versus two treatment cycles for the proximal vein segment. Phlebology 2015, 30, 724–728. [CrossRef] Sufian, S.; Arnez, A.; Labropoulos, N.; Lakhanpal, S. Endovenous heat-induced thrombosis after ablation with 1470 nm laser: Incidence, progression, and risk factors. Phlebology 2015, 30, 325–330. [CrossRef] [PubMed] Tsai, C.T.; Hsieh, C.S.; Chang, S.N.; Chuang, E.Y.; Ueng, K.C.; Tsai, C.F.; Lin, T.H.; Wu, C.K.; Lee, J.K.; Lin, L.Y.; et al. Genome-wide screening identifies a KCNIP1 copy number variant as a genetic predictor for atrial fibrillation. Nat. Commun. 2016, 7, 10190. [CrossRef] [PubMed] Li, R.; Yu, C.; Li, Y.; Lam, T.W.; Yiu, S.M.; Kristiansen, K.; Wang, J. SOAP2: An improved ultrafast tool for short read alignment. Bioinformatics 2009, 25, 1966–1967. [CrossRef] [PubMed] Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 2008, 5, 621–628. [CrossRef] [PubMed] Kramer, A.; Green, J.; Pollard, J., Jr.; Tugendreich, S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics 2014, 30, 523–530. [CrossRef] [PubMed] Lawson, N.D.; Weinstein, B.M. In vivo imaging of embryonic vascular development using transgenic zebrafish. Dev. Biol. 2002, 248, 307–318. [CrossRef] Chen, Y.H.; Lin, Y.T.; Lee, G.H. Novel and unexpected functions of zebrafish CCAAT box binding transcription factor (NF-Y) B subunit during cartilages development. Bone 2009, 44, 777–784. [CrossRef] Kimmel, C.B.; Ballard, W.W.; Kimmel, S.R.; Ullmann, B.; Schilling, T.F. Stages of embryonic development of the zebrafish. Dev. Dyn. 1995, 203, 253–310. [CrossRef] Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [CrossRef] Hen, G.; Nicenboim, J.; Mayseless, O.; Asaf, L.; Shin, M.; Busolin, G.; Hofi, R.; Almog, G.; Tiso, N.; Lawson, N.D.; et al. Venous-derived angioblasts generate organ-specific vessels during zebrafish embryonic development. Development 2015, 142, 4266–4278. [CrossRef] Camenisch, T.D.; Spicer, A.P.; Brehm-Gibson, T.; Biesterfeldt, J.; Augustine, M.L.; Calabro, A., Jr.; Kubalak, S.; Klewer, S.E.; McDonald, J.A. Disruption of hyaluronan synthase-2 abrogates normal cardiac morphogenesis and hyaluronan-mediated transformation of epithelium to mesenchyme. J. Clin. Investig. 2000, 106, 349–360. [CrossRef] [PubMed] © 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).