Flow induces barrier and glycocalyx-related genes and negative surface charge in a lab-on-a-chip human blood-brain barrier model

Abstract

Microfluidic lab-on-a-chip (LOC) devices allow the study of blood-brain barrier (BBB) properties in dynamic conditions. We studied a BBB model, consisting of human endothelial cells derived from hematopoietic stem cells in co-culture with brain pericytes, in an LOC device to study fluid flow in the regulation of endothelial, BBB and glycocalyx-related genes and surface charge. The highly negatively charged endothelial surface glycocalyx functions as mechano-sensor detecting shear forces generated by blood flow on the luminal side of brain endothelial cells and contributes to the physical barrier of the BBB. Despite the importance of glycocalyx in the regulation of BBB permeability in physiological conditions and in diseases, the underlying mechanisms remained unclear. The MACE-seq gene expression profiling analysis showed differentially expressed endothelial, BBB and glycocalyx core protein genes after fluid flow, as well as enriched pathways for the extracellular matrix molecules. We observed increased barrier properties, a higher intensity glycocalyx staining and a more negative surface charge of human brain-like endothelial cells (BLECs) in dynamic conditions. Our work is the first study to provide data on BBB properties and glycocalyx of BLECs in an LOC device under dynamic conditions and confirms the importance of fluid flow for BBB culture models.

Introduction

Due to the importance of the blood-brain barrier (BBB) in many systemic as well as central nervous system (CNS) diseases ¹ improved in vitro models are crucial as research tools for the biomedical community. In the last decade, several in vitro BBB models have been developed using cell culture inserts representing static conditions and more complex lab-on-a-chip (LOC) devices. ² , ³ Fluidic and microfluidic LOCs provide fluid flow and enable the study of BBB in a dynamic condition, ⁴ , ⁵ in contrast to static cell culture inserts.

The morphology of the BBB, the structure of brain endothelial cells (EC) connected by tight intercellular junctions (TJ) and that of surrounding pericytes, astroglia endfeet, perivascular microglia and neuronal processes are well described. ⁶ Recent advances in gene sequencing and single cell analysis led to better understanding of the molecular composition of cell populations of the brain including the cell types forming the BBB. ⁷ , ⁸ More complete lists of influx and efflux transporters, receptors and transport pathways are available for the different types of BBB cells which participate in providing nutrients for the CNS and also act as protection systems. ⁶ , ⁹ The BBB shields the CNS from toxins and pathogens but also participates in the pathomechanism of CNS injuries and neuroinflammation. ¹ , ⁶

In addition to physical defense mechanisms provided by TJs, and chemical protection by efflux pumps and BBB metabolic enzymes, the fourth major line of defense is the endothelial surface glycocalyx (ESG). Blood flow induces mechanical forces acting on the surface of ECs, which via mechanosensors, signaling pathways and gene expression, modulate endothelial morphology and function. ¹⁰ , ¹¹ ESG is a sugar-protein matrix like layer covering the EC surface mainly composed of proteoglycans, glycoproteins, and glycosaminoglycans. Its unique location, composition and structure serves as a negatively charged physical barrier on the surface of ECs. ¹¹ A recent in vivo study highlights the denser structure of glycocalyx of microvessel in the brain as compared to lung and heart, and associates it with brain EC protection as a defense component of the BBB. ¹² Indeed, the surface charge of brain ECs is more negative than that of other vascular ECs, which was related to a higher level of negatively charged phosphatidylserine and phosphatidylinositol in their cell membrane, ¹³ in addition to the negative surface charge derived from the sulfate and sialic acid residues of the ESG. ¹⁰ , ¹⁴ The structural complexity and the negative surface charge of the ESG at the BBB not only provide an extra barrier on the EC, but also regulate the penetration of large molecules ¹⁵ and charged drugs or vectors to the CNS. ¹⁴ , ¹⁶ The protective role of ESG is well known in the cardiovascular system, and damage of this surface layer was demonstrated in diseases and pathologies like atherosclerosis, ischemia and inflammation. ¹² , ¹⁷ The characterization of ESG at the level of the BBB is incomplete and its importance is not fully understood.

Since LOCs provide more physiological conditions, in the present study we characterized a human BBB model composed of ECs derived from hematopoietic stem cells co-cultured with brain pericytes ¹⁸ in our microfluidic and microelectronic device. ⁵ Furthermore, our goal was to investigate the effect of fluid flow on gene expression of brain-like endothelial cells (BLECs) by in-depth massive analysis of cDNA ends sequencing (MACE-seq) with a focus on general endothelial, BBB-related and ESG-related genes. To confirm the results endothelial morphology, immunostaining of selected BBB proteins, lectin and immunostaining of ESG and surface charge measurements by laser Doppler velocimetry were performed. Our present work, together with our previous study revealing the importance of brain EC surface charge in the permeability of charged molecules, ¹⁶ draws the attention to ESG as a flow-regulated essential part of BBB models.

Material and methods

Materials

All materials used in the study were purchased from Sigma-Aldrich, Hungary Ltd. (a subsidiary of Merck, Germany), unless otherwise indicated.

Cell culture

The in vitro BBB model, consisting of human endothelial cells (hEC) in co-culture with bovine brain pericytes, was described by Cecchelli et al. in 2014. ¹⁸ The isolation of stem cells required the collection of human umbilical cord blood, which was approved by the Hospital ethical committee (Béthune Maternity Hospital, Béthune, France). Informed consent was obtained from the infants’ parents. The protocol was approved by the French Ministry of Higher Education and Research (CODECOH Number DC2011-1321). The study was conducted according to World Medical Association Declaration of Helsinki. The hECs were derived from CD³⁴⁺ cord blood hematopoietic stem cells by a differentiation period of 15-20 days using VEGF₁₆₅. ¹⁸ After differentiation hECs were seeded in 0.2% gelatin coated culture dishes and kept in endothelial cell culture medium (ECM, Sciencell, USA) supplemented with 5% fetal bovine serum (FBS), 1% endothelial cell growth supplement (ECGS, Sciencell), and gentamycin (50 µg/ml). After reaching confluency, hECs were gently trypsinized and 80 × 10³ cells seeded into the polyester membrane of the lab-on-a-chip (LOC) device (Figure 1(b)), coated with Matrigel (BD Biosciences, USA). The LOC was built as described in our previous publication, ⁵ a brief protocol is provided in the Supplementary Methods. Bovine brain pericytes (PC) were cultured in 0.2% gelatin coated dishes in Dulbecco’s modified Eagle’s medium (DMEM, Life Technologies, Thermo Fisher Scientific, USA) supplemented with 20% FBS, 1% Glutamax (Life Technologies) and gentamicin (50 µg/ml). When cultures reached confluency, PCs were trypsinized and 25 × 10³ cells were added to the bottom compartment of the LOC device coated with 0.2% gelatin (Figure 1(b)). During co-culture both compartments received endothelial medium. The LOC device with the cells was kept at 37 °C in a humified atmosphere and 5% CO₂. For more details, see Supplementary Methods.

Figure 1.

Characterization of the human BBB model cultured in the lab-on-a-chip (LOC). (a) Representative illustration of the device. The plastic slides carrying the gold electrodes to measure transendothelial electric resistance (TEER) are positioned at the top and bottom of the device, followed by the top and bottom channels made by PDMS and the cell culture membrane in the middle. The layers of the LOC were joined with screws. The female luer inlets were located on the top and provided easy access for both top and bottom channels. (b) Diagonal view of the LOC. Human brain endothelial cells were added to the top compartment, brain pericytes to the bottom. (c) Dynamic condition: the device was connected to a peristaltic pump and a reservoir containing cell culture medium. Fluid flow was applied at a speed of 1 ml/min for 24 hours. (d) Phase contrast images of brain endothelial cells under static and dynamic conditions. Scale bar: 100 μm. (e) TEER results were normalized to the values of the static condition which did not receive any fluid flow, values are presented as means ± SD, unpaired t-test, **p < 0.01, n = 12. (f) Apparent permeability coefficient (Papp) of the human BBB model under static and dynamic conditions, for Lucifer yellow (LY) and Evans blue labelled albumin (EBA) marker molecules. Data is shown as the % of the static condition and presented as means ± SD (unpaired t-test, ***p < 0.001, ****p < 0.0001, n = 5–8).

Dynamic culture conditions in the LOC device

To determine the importance of fluid flow on BLECs, the co-culture of BLECs with PCs lasted for 6 days under static condition (Supplementary Figure S1) then 24 hours under dynamic condition (Figure 1(c)), and compared with cells which were kept under static condition for 7 days (for details see Supplementary Methods). Cell growth was monitored by phase contrast microscopy (MotiCam 1080, MoticEurope, Barcelona, Spain), since the transparency of the gold electrodes in the LOC device enables the visualization of the entire cell monolayer. TEER was measured every day. Upon reaching the 6th day, a constant circulation of culture medium was introduced by a peristaltic pump (Masterflex, Cole-Parmer, USA) at 1 ml/min flow rate (0.4 dyne/cm²) for 24 hours for the model in dynamic conditions. A higher, physiological shear stress (1.6 dyne/cm²) was also tested by the combination of increased medium viscosity (3.5% dextran) and elevated flow rate (2.5 ml/min). Under our experimental conditions the flow in the channel was laminar based on the Reynolds’ number. Co-cultures in static conditions were kept as shown in Supplementary Figure S1. Right after this 24-hour static and dynamic condition we performed barrier integrity studies on the human BBB model, did RNA extraction of BLECs, stained for BBB proteins and glycocalyx and measured the surface charge of BLECs.

Barrier integrity measurement: Transendothelial electrical resistance and permeability

Cells were kept in co-culture for 7 days and TEER values were recorded daily. The TEER was measured as described in our previous work: ⁵ the LOC electrodes were connected to the EVOM² voltohmmeter and the TEER values were registered (Figure 1(c)). To assess the barrier integrity small molecular weight marker Lucifer yellow (LY, MW: 457 Da) ¹⁸ and Evans blue dye bound to 1% bovine serum albumin (EBA, MW: 67.5 kDa) ⁵ were used. Medium in the upper compartments of the LOC devices was replaced with 150 µl of Ringer-Hepes solution containing LY (5 µM) and EBA (165 µg/ml dye bound to 1% albumin). In the bottom compartments 350 µl of Ringer-Hepes solution was added. LOCs were incubated for 20, 40 and 60 minutes in a CO₂ incubator at 37 °C on a horizontal shaker (150 rpm/min). Samples were collected from both compartments and fluorescent intensity measured by spectrofluorometer (Horiba Fluorolog 3, Japan). For details of the permeability assay and calculation of apparent permeability coefficients (P_app), see Supplementary Methods.

Immunocytochemistry

Cells were fixed with ice cold methanol and acetone solution (1:1) for 2 minutes, washed twice with phosphate buffered saline (PBS) containing 1% FBS and blocked with 3% BSA in PBS for 1 hour. BLECs were incubated with primary antibodies for β-catenin, claudin-5, ZO-1, P-glycoprotein (PGP), glucose transporter-1 (GLUT-1), ICAM-1, VCAM-1 diluted in 3% BSA-PBS blocking buffer and incubated overnight at 4 °C. Cells were washed in PBS and incubated with secondary antibody anti-rabbit IgG-CY3 or anti-mouse IgG A488 and H33342 dye for 1 hour. The same protocol was used to stain brain PCs for α-smooth muscle actin, NG2 and PDGFR-β. Alexa fluor 488 anti-mouse IgG, anti-rabbit IgG-CY3 were used as secondary antibodies. Samples were visualized by a Leica TCS SP5 confocal laser scanning microscope. For more details, see Supplementary Methods.

Massive analysis of cDNA ends library preparation and RNA sequencing

RNA isolation and RNA integrity analysis were performed as described in Supplementary Methods. Genome-wide gene expression profiling was performed using massive analysis of cDNA ends (MACE-seq) with RNA extracted from BLECs co-cultured with PCs in static and dynamic conditions (3 and 5 biological replicates, respectively). A total of 8 libraries were constructed using the Rapid MACE-Seq kit (GenXPro GmbH, Frankfurt, Germany), according to the manufacturer’s protocol. MACE-seq performs gene expression profiling by sequencing part of the 3′ end of mRNA transcipts. While synthetizing one cDNA molecule from each mRNA transcript, MACE-seq can accurately quantify transcribed polyadenylated transcripts. For more details, see Supplementary Methods and Figure S2.

Bioinformatic analysis of MACE-seq data

Approximately 62 million MACE-seq reads were obtained across all 8 libraries, and subsequently processed and analyzed bioinformatically. PCR-duplicates were identified using the TrueQuant technology and subsequently removed from raw data. The remaining reads were further poly(A)-trimmed and low-quality reads were discarded. In the following step, the clean reads were aligned to the human reference genome (hg38, http://genome.ucsc.edu/cgi-bin/hgTables) using bowtie2 mapping tool resulting in a dataset of 28,834 different genes. The gene count data was normalized to account for differences in library size and RNA composition bias by calculating the median of gene expression ratios using DESeq2 R/Bioconductor package. ¹⁹ Testing for differential gene expression was also performed using the DESeq2 R/Bioconductor package. As a result, p-value and log2(fold change) (log2FC) were obtained for each gene in the dataset. False discovery rate (FDR) analysis was estimated to account for multiple testing. The gene expression level was calculated for each gene using the normalized value for the dynamic condition divided by the normalized value for the static condition multiplied by 100. Genes with a p-value < 0.05 and |log2FC| > 1 were considered to be differentially expressed. To perform functional profiling analysis g:Profiler was used to identify over-represented biochemical pathways from 3 databases (KEGG, Reactome and Gene Ontology) and to calculate the statistical significance of each pathway. GOplot was used to calculate the zscore from each over-represented pathway. ²⁰ The transcriptomic datasets generated and analyzed in this study are available in the Gene Expression Omnibus (GEO) repository. ²¹ GEO accession number is GSE155671 and will also be available at the BBBHub (http://bbbhub.unibe.ch).

Zeta potential measurements

The zeta potential reflecting the surface charge of brain endothelial cells was measured by laser Doppler velocimetry using a Zetasizer Nano instrument (Malvern Instruments, UK) equipped with a He-Ne laser (λ = 632.8 nm) as described previously. ¹⁶ For more details, see Supplementary Methods.

Staining of cell surface glycocalyx

After 24-hour dynamic and static condition BLECs in co-culture with PCs were washed and fixed by 1% paraformaldehyde in PBS. Cells were not permeabilized to ensure only surface labeling. Sialic acid and N-acetyl-D-glucosamine residues within the glycocalyx were labeled by A488 conjugated wheat germ agglutinin (WGA) lectin (W11261, Invitrogen; 5 µg/ml). ²² BLECs were also immunostained for chondroitin sulfate (anti-chondroitin sulfate mouse monoclonal primary antibody, Merck, Germany; 1:100; AB_476879). Pictures were taken with Leica TCS SP5 confocal laser scanning microscope at random positions (14 images/membrane; three independent experiments). Fluorescent images were analyzed for staining intensity using the FIJI (ImageJ) software. For details see Supplementary Methods.

Statistics

GraphPad Prism 5.0 Software (GraphPad Software, USA) was used for statistical analysis. Gaussian distribution of the data was tested with the KS normality test. Data showing Gaussian distribution were analyzed with unpaired t-test or one-way ANOVA followed by Bonferroni post-test. Results were considered statistically significant at p-values < 0.05. Data are represented as means ± SD. Experiments were repeated at least twice and a minimum number of 3 biological parallels were used.

Results

Human BBB model in the LOC device: Effect of flow on barrier properties

The human BBB model has been characterized and comparisons of BLECs in monoculture vs. co-culture with brain pericytes were described in static culture conditions. ¹⁸ , ²³ Our first goal was to optimize and characterize this human BBB model in the LOC device under flow. The LOC has two fluid compartments, similarly to culture inserts: a porous PET membrane separates the top and bottom channels (Figure 1(a) and (b)). One of the biggest advantages of the LOC compared to the cell culture insert setup is the possibility of applying fluid flow to induce shear stress to mimic more physiological conditions for brain EC cultures. ⁵ As in the case of cell culture inserts, ¹⁸ , ²³ BLECs were added to the top compartment, and brain PCs were seeded to the bottom compartment (Figure 1(b)). During the first 6 days feeding was applied every 8 hours (Supplementary Figure 1), followed by a 24-hour dynamic condition (flow rate: 1 ml/min; 0.4 dyne/cm²; Figure1(c)). For the static model 7 days of automatic feeding was used and no fluid flow was applied. The tightness of the barrier, which reflects the permissiveness of the monolayer for the passage of ions and molecules was characterized by measurement of TEER and permeability for LY and EBA. ⁵ As shown in Figure 1(e) 24-hour fluid flow elevated the TEER significantly by 18% (361.8 ± 166.3 Ω×cm² to 425.5 ± 188.8 Ω×cm²). The paracellular permeability for LY was significantly decreased after fluid flow by 78% (16.3 ± 3.6 to 3.6 ± 1.4 × 10⁻⁶cm/s) similarly to permeability for EBA, which was significantly decreased by 93% (1.77 ± 0.02 to 0.14 ± 0.01 ×10⁻⁶cm/s) after dynamic condition (Figure 1(f)). The morphology of BLECs also changed after dynamic as compared to static condition. After flow a more elongated cell shape was visible (Figure 1(d)) than in cells under static condition. No difference was observed between the permeability or TEER values of the low and high shear groups in dynamic conditions (Figure S3(b) and (c)).

Transcriptomic profile of the human BBB model: Effect of flow on the global gene expression profile

To further characterize the effect of fluid flow on the human BBB co-culture model, we analyzed the transcriptomic gene expression profile of BLECs using the MACE-seq. Gene expression profiling analysis was able to detect expression from 28,807 different genes (Supplementary Figure S4(a)). A total of 396 genes had a p-value lower than 0.05 and an |log2FC| >1, which were considered differentially expressed (Figure 2(b)). From all the expressed genes, 174 (0.6%) and 222 (0.8%) were up and down regulated, respectively (Figure 2(a)). Using more stringent thresholds of |log2FC| >2 or <–2, 67 (0.2%) upregulated and 73 (0.3%) downregulated genes were identified, respectively. Hierarchical clustering and principal component analysis (PCA) were used to assess relatedness between samples (Supplementary Figure 4(a) and (b)). PCA and clustering analysis grouped the expression profiles from dynamic and static condition in two separate clusters. Many of the top 50 differentially expressed, up- or downregulated genes (Supplementary Table S1-S3) are related to cell adhesion, basal membrane composition, surface glycocalyx, enzymatic or transport processes at the BBB, which we discuss in the relevant sections below.

Figure 2.

(a) Number of differentially expressed transcripts of human endothelial cells co-cultured with brain pericytes under dynamic and static conditions. The total number of differentially expressed transcripts are shown on top of the bars. Testing for differential gene expression was performed using the DESeq2 R/Bioconductor package. ¹⁹ (b) The volcano plot identifies the total changes in the data set. Black: no change, Blue: downregulated genes, Red: upregulated genes.

Transcriptomic gene expression profile of the human BBB model: Effect of flow on endothelial cell markers

During the analysis of MACE data, we first wanted to verify the influence of fluid flow on the expression of general endothelial marker and other endothelial cell-related genes. Several basic endothelial marker genes were expressed in BLECs, including von Willebrand factor (VWF), vascular endothelial growth factor genes (VEGFA, B, C) and their receptors (VEGFR1-3), angiopoietins (ANGPT1, -2), endothelial cell specific molecule 1 (ESM1) and endothelial nitric oxide synthase (NOS3), as expected. However, many of these genes were not changed by flow (Figure 3(a)). Among the endothelial genes 24-hour flow significantly upregulated endothelin 1 (EDN1), nitric oxide synthase interacting protein (NOSIP) and vascular endothelial growth factor receptor 1 (VEGFR1) (Figure 3(a)) and significantly (13-fold) downregulated the lymphatic vessel endothelial marker gene LYVE1 to a very low expression level. In dynamic condition among the adhesion molecules (Figure 3(b)), the expression of intercellular adhesion molecule-1 gene (ICAM1) was significantly elevated. Other surface molecules were not changed by flow, such as adhesion regulating molecule-1 (ADRM1), activated leukocyte cell adhesion molecule (ALCAM), ICAM2, E- and P-selectins (SELE, -P) and vascular cell adhesion molecule-1 (VCAM1). Flow had no effect on endothelial cytoskeleton genes in the human BBB model, including actinins (ACTN), filamin (FLN), talin (TLN), tensin (TNS), vinculin (VCL) and vimentin (VIM), while significantly increased the expression of α-actinin-1 (ACTN1) and the actin cross-linker transgelin (TAGLN) genes. Importantly, genes for several members of the integrin family, especially α-subunits (ITGA5, -11, -E, -V) were overexpressed after flow (Figure 3(b)), along with a number of brain endothelial basal membrane genes like collagen type IV (COL4A1), as well as fibulin 5 (FBLN5), FNDC3B and members of the matrix metallopeptidase family (MMP1, -2, -10, -14) (Figure 3(c)). From the members of the ADAM (a disintegrin and metalloproteinase) family expressed in the human BBB model only one gene, ADAM-12 was upregulated by flow. Several basal membrane genes were unchanged by flow, including dystrophin (DMD), elastin (ELN), perlecan or basement membrane-specific heparan sulfate proteoglycan core protein (HSPG2), laminins (LAM) and vitronectin (VTN). The significant upregulation of genes for cytoskeleton proteins, integrins and basal membrane molecules indicates pathways leading to better attachment of the cell layer at the abluminal side and demonstrates the importance of dynamic conditions for BBB models to better mimic physiological conditions in vitro.

Figure 3.

Transcriptomic gene expression profile of (a) Endothelial cell specific genes. (b) Endothelial cell surface molecules, cytoskeleton molecule and integrin genes; (c) Basal membrane molecule and matrix metalloproteinase (MMP) genes. Expression is shown as the relative expression (%) of the genes present in human endothelial cells co-cultured with brain pericytes in dynamic condition as compared to static condition. Testing for differential gene expression was performed using the DESeq2 R/Bioconductor package. ¹⁹ Genes with a p-value <0.05 and less than 50% or more than 200% gene expression levels were considered to be differentially expressed. Red color labels upregulation and statistically significant expression changes, blue color shows downregulation and statistically significant expression change, cream color indicates no change in the gene expression (*p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001).

Transcriptomic gene expression profile of the human BBB model: Effect of flow on BBB-related genes

Three main junction groups participate in the formation of BBB: tight junctions (TJ), adherens junctions (AJ) and gap junctions (GJ).^24,25 These intercellular connections enable brain endothelial cells to form a dynamic interface between the blood and CNS. ⁶ Several TJ transmembrane or linker proteins were present at an unchanged level after the dynamic condition, including endothelial cell-selective adhesion molecule (ESAM), occludin (OCLN), claudin-12 (CLDN12), junctional adhesion molecule 3 (JAM3), MARVELD1 and -D2, zonula occludens protein-1 and -2 (TJP1 and TJP2) (Figure 4(a)). The gene expression level of JAM1, as well as that of CLDN1, -3, -5 and -7 were decreased after flow condition (Figure 4(a)). Among AJ proteins the gene expression level of vascular endothelial or VE-cadherin (CDH5), neuronal or N-cadherin (CDH2), nectin-2 and -3, α- and β-catenins (CTNNA1, -B1) remained unchanged, while the expression levels of protocadherins-1 and -9 (PCDH1, -9) genes increased. Genes for gap junction proteins, known to be involved in intercellular communication and the regulation of TJs and AJs, connexin 43 (GJA1) and connexin 40 (GJA5) were significantly upregulated by dynamic condition (Figure 4(a)).

Figure 4.

Transcriptomic gene expression profile of (a) tight, adherens and gap junction protein genes. (b) drug efflux transporter and drug metabolizing enzyme genes. Expression is shown as the relative expression (%) of the genes present in human endothelial cells co-cultured with brain pericytes in dynamic condition as compared to static condition. Testing for differential gene expression was performed using the DESeq2 R/Bioconductor package. ¹⁹ Genes with a p-value < 0.05 and less than 50% or more than 200% gene expression levels were considered to be differentially expressed. Red color labels upregulation and statistically significant expression changes, blue color shows downregulation and statistically significant expression change, cream color indicates no change in the gene expression (*p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001).

Active efflux pumps and drug metabolizing enzymes represent another line of defense at the level of BBB by providing protection against toxic chemical compounds. As shown in Figure 4(b) the gene expression level for many of the important ATP cassette binding (ABC) and other transporters participating in drug and metabolite efflux was unchanged after flow condition: ABCB1 (PGP), ABCG2 (breast cancer resistance protein, BCRP), members of the ABCC family (ABCC1-6 or MRP1-6) and the excitatory amino acid transporter-3 belonging to the solute carrier (SLC) family (SLC1A1/EAAT3). Other transporters, like the lipid transporters ABCA3 and ABCA8, expressed in both peripheral and brain tissues, were significantly downregulated after the dynamic condition along with the organic anion-transporting polypeptide 1 (SLCO4A1/OATP1) and the organic cation transporter-2 (SLC22A5/OCTN2). The only drug efflux transporter which gene was upregulated by flow was SLC22A5 (OCTN1). The cytochrome (CYP) P450 enzymes, which participate in steroid and drug metabolism (so called drug metabolism phase I enzymes) and thus contribute to subsequent efflux of conjugated drugs at the level of BBB, are expressed in different barrier culture models. ²⁶ Many members of the CYP family were expressed in the human BBB model and CYP1A1 and CYP1B1 genes were significantly overexpressed after flow (Figure 4(b)). Additionally, drug metabolism phase II enzymes of the BBB, epoxide hydrolase-1 (EPHX1) and glutathione S-transferase π (GSTP1) were also upregulated by flow (Figure 4(b)).

SLCs including glucose, monocarboxylate, amino acid and peptide transporters were expressed either with an unchanged level or were significantly overexpressed after flow in BLECs, except for SLC7A5/LAT1, which was downregulated (Supplementary Figure S5(a)). Among ion transporters and channels expressed at the BBB, ¹ two ion transporter genes were upregulated, while potassium channel genes were downregulated by flow condition (Supplementary Figure S5(b)). The gene expression of many important BBB receptor genes was unchanged, except LRP3 and LRP5 which were increased and LDLRAD3, very low-density lipoprotein receptor (VLDLR) and caveolin-1 (CAV1) which were decreased (Supplementary Figure S5(b)).

We characterized the BLEC/PC model and verified the flow induced changes for BBB proteins by immunostaining (Figure 5). Brain PCs stained for α-smooth muscle actin and NG2 and showed a uniform morphology (Figure 5(a)). In the MACE study CLDN5 expression levels were lower after fluid flow. In the fluorescence intensity evaluation no significant change, only a trend of decrease was observed for claudin-5 immunostaining (Figure 5(b)). Flow did not change the expression of β-catenin and ZO1 genes which was supported by the analysis of immunostaining for these junctional proteins (Figure 5(c) and (d)). No change was seen for PGP and GLUT1 gene expression and immunostaining (Figure 5(e) and (f)). ICAM1 gene expression was elevated after fluid flow, which was corroborated with the image analysis of the staining (Figure 5(g)). Flow did not alter VCAM-1 expression at gene or protein level (Figure 5(h)). The immunostaining for claudin-5 and β-catenin was also compared between the two dynamic conditions, and in the higher shear group a decreased immunostaining intensity was found for both junctional proteins (Figure S2(d) to (f)).

Figure 5.

Immunocytochemistry of brain pericytes and human brain-like endothelial cells. (a) Brain pericytes were characterized by immunostaining for platelet-derived growth factor (PDGF)-β receptor, α-smooth muscle actin (α-SMA) and neuron-glial antigen 2 (NG2). Brain-like endothelial cells in static and dynamic conditions were characterized by immunostaining followed by fluorescence intensity evaluation (n = 12–14) for: (a–c) junctional proteins claudin-5, zonula occludens – 1 and β-catenin; (e-f) efflux pump P-glycoprotein (PGP) and solute carrier glucose transporter-1 (GLUT1); (g-h) cell adhesion molecules intercellular adhesion molecule-1 (ICAM-1) and vascular cell adhesion molecule-1 (VCAM-1). Scale bar: 20 μm. Unpaired t-test, **p < 0.01, n.s. non significant.

Effect of flow on the human BBB model: Endothelial surface charge and glycocalyx

The endothelial surface glycocalyx (ESG) plays an important role in the mechano-sensing and transduction functions of endothelial cells, which are essential to maintain vascular integrity and homeostasis. ¹⁰ ESG is a rich layer of carbohydrates connected through proteoglycans and glycoproteins to the surface of ECs (Figure 6(a) and (b)). ²⁷ Most of the in vitro BBB models based on cell culture inserts lack fluid flow, therefore represent a static environment. While the importance of the ESG and extracellular matrix are in general acknowledged as elements of the defense system of the BBB, this is an under researched area. This is the reason why we aimed to investigate the transcriptomic gene expression changes of the surface glycocalyx and the extracellular matrix components after static and dynamic conditions of BLECs co-cultured with PCs in our LOC. As shown in Figure 6(c), glycocalyx core proteins, like decorin (DCN), glypican-1 (GPC1), syndecan-2 (SDC2) and versican (VCAN) genes were significantly upregulated in the dynamic condition. Other ESG core protein genes like biglycan (BGN), CD44, and other syndecans (SDC1, -3, -4) were also present, but their expression levels were not changed by flow. Four types of galectins that bind specifically to β-galactoside sugars within the glycocalyx were expressed in BLECs (GAL1, -3, -8, -9), of which the GAL3 gene was decreased by flow (Figure 6(c)). The human BBB model expressed a large number of enzymes which participate in the synthesis and remodeling of the ESG (Figure 6(c)). The genes of four enzymes, the carbohydrate sulfotransferase-1 (CHST1), heparanase (HPSE), hyaluronidase-2 (HYAL2) and the heparan sulfate 6-O-endosulfatase SULF2, which selectively removes 6-O-sulfate groups from heparan sulfate, were found to be upregulated by flow, while the expression of only one enzyme, β-1,4-galactosyltransferase-5 (B4GALT5) decreased. The expression changes in the ESG core proteins and galectins were caused by flow, since the effect of co-culture with PCs in static conditions did not induce such changes (Supplementary Figure S6(a)).

Figure 6.

(a) Illustration showing the effect of fluid flow on the endothelial surface glycocalyx (ESG). (b) Drawing representing ESG core proteins, glycoproteins and glucosaminoglycans. HA: hyaluronic acid, HS: heparan sulfate, CS: chondroitin sulfate. (c) Transcriptomic gene expression profile of ESG-related genes. Expression is shown as the relative expression (%) of the genes present in human endothelial cells co-cultured with brain pericytes in dynamic condition as compared to static condition. Testing for differential gene expression was performed using the DESeq2 R/Bioconductor package. ¹⁹ Genes with a p-value < 0.05 and less than 50% or more than 200% gene expression levels were considered to be differentially expressed. Red color labels upregulation and statistically significant expression changes, blue color shows downregulation and statistically significant expression change, cream color indicates no change in the gene expression (**p < 0.01, ***p < 0.001, ****p < 0.0001). (d) Zeta potential measured by laser Doppler velocimetry (means ± SD, n = 10; unpaired t-test, * p < 0.05 compared to static condition). (e) and (f) Fluorescent intensity analysis of the images showing wheat germ agglutinin (WGA) lectin staining and chondroitin sulfate immunocytochemistry. Image analysis values are presented as means ± SD, n = 14 images/groups; unpaired t-test, **p < 0.01 compared to static condition. (g) and (h) Staining of ESG on brain endothelial cells with fluorescently labeled WGA lectin that binds to the sialic acid residues and with chondroitin sulfate immunocytochemistry. Scale bar: 20 µm.

Brain EC glycocalyx is one of the thickest within the vasculature ¹² and also cultured brain ECs have highly negative surface charge. ¹⁶ The main sources of this negative charge are the lipid head groups of the plasma membrane and the ESG, composed of highly negatively charged polysaccharide chains. Since ESG related genes were upregulated by flow, zeta potential measurement was performed to determine if the changes at gene expression influence the absolute value of the surface charge of BLECs. We measured a significant decrease in surface charge of BLECs (Figure 6(d)) in dynamic condition (–12.4 ± 1.4 mV) compared to static condition (–11.0 ± 1.0 mV). We also compared the effect of co-culture with PCs in static conditions using cell culture inserts, and found that the negative surface charge of BLECs in the co-culture model decreased (–11.7 ± 1.2 mV) compared to the monoculture (–10.6 ± 1.4 mV) (Supplementary Figure 6(b)), indicating the importance of both flow and co-culture in the negative surface charge of brain ECs.

The more negative surface was related to a denser ESG as confirmed by the fluorescently labelled WGA lectin staining: a 50% increase in the fluorescence intensity was observed in the confocal images after 24-hour flow (Figure 6(e) and (g)). A 30% increase in WGA lectin staining was also observed in BLECs in co-culture with PCs as compared to monocultures in static conditions (Supplementary Figure S6(c) and (d)). The BLECs were immunolabeled for the surface glycosaminoglycan chondroitin sulfate (CS), but no change was found in the staining intensity (Figure 6(h) and (f)). We saw no difference between the lower and the higher shear stress BLEC groups for WGA lectin staining (Figure S2(h) and (i)) or for zeta potential (Figure S2(g)).

The effect of flow on the ESG was corroborated by the gene enrichment profiles (Figure 7). We found that the most significantly upregulated pathways after the introduction of flow conditions were the extracellular matrix and structure pathways and ESG related pathways.

Figure 7.

Functional profiling analysis of extracellular matrix-related pathways in static versus dynamic condition. The x-axis represents the statistical significance calculated using g:Profiler while the zscore represents the tendency of the regulation of these pathways calculated using GOplot.

Discussion

RNA sequencing analysis provided new data on mouse brain EC transcriptomics ⁷ and on zonation of ECs in mouse brain vasculature. ⁸ Culture models of the BBB in static condition have also been characterized at the transcriptomic level ²⁶ , ²⁸ including the BLEC model we use in the present study. ²³ We found one study on the effect of fluid flow on the transcriptome of a human cell-based BBB model, ⁴ but no study analyzed transcriptomic changes of a human BBB model in an LOC-based system in dynamic condition yet.

The BBB is composed of brain ECs that have strong interaction with the neighboring astroglia endfeet and pericytes. Brain pericytes share a common basal membrane with the ECs and they have an important role in the development, maintenance, and regulation of BBB functions. ¹ , ²⁹ Co-culture of brain ECs with PCs is known to elevate the tightness of the paracellular barrier and increase the BBB properties. ¹⁸ , ³⁰ , ³¹ The human BBB model used in our study, in which hECs are derived from hematopoietic stem cells and co-cultured with bovine brain PCs, is a well characterized in vitro model. ¹⁸ , ²³ , ³² There is only one recent study which describes the adaptation of this BBB model to the use of short-term (4 and 30 minutes) fluid flow with the goal to study immune cell transmigration. ³² However, no longer term (24-hour) fluid flow and its effects on this BBB model have been investigated until now.

We evaluated the effects of two dynamic conditions on the BBB model in the LOC device, a 0.4 dynes/cm² shear stress and a four-times higher, 1.6 dynes/cm² condition, typical for brain postcapillary venules. To sum up our observations, we can conclude that the higher shear condition was well-tolerated by the BLECs, and was similarly beneficial to the barrier integrity elevation and glycocalyx expression as the lower shear condition. The 0.4 dynes/cm² shear is lower than that in brain microvessels in vivo, and is a limitation of the study. However, this lower shear stress provided sufficient dynamic condition for the BLEC model to observe improved barrier integrity and the expression of important BBB genes and characteristics. The characterization of a human co-culture BBB model under static and dynamic conditions in an LOC device and the global transcriptomic gene expression analysis may provide valuable data to further improve BBB LOC models.

First, we measured the barrier properties of the BLEC and bovine PC co-culture model in the LOC device. In static conditions in the device similar TEER values were obtained for this human BBB model as for several types of in vitro co-culture models, including primary cell based ones. ²⁶ After introducing fluid flow, TEER values significantly increased and permeability of the BLEC monolayers significantly decreased for both LY and EBA markers. Similar changes were found by previous studies using human BBB culture models in LOC devices. ⁴ , ⁵ The permeability values obtained for the marker molecule LY in dynamic conditions in the LOC device was lower than those in previous studies using culture inserts representing static conditions. ¹⁸ , ³² In contrast to the dynamic in vitro system using hollow fiber cartridges ⁴ our LOC device ⁵ allowed the morphological observation of BLECs, which showed similar immunostaining of junctional proteins as in the culture insert models. ¹⁸ , ³² Fluid flow changed cell shape to be more elongated and also realigned the cells in the direction of flow, in accordance with a previous study on bovine brain endothelial cells. ³³ These observations prove, that the LOC device is functional and show the importance of fluid flow in barrier integrity regulation of in vitro human BBB models. ⁴ , ⁵

Blood flow in capillaries and the resulting shear force are important physiological regulators of EC functions in the periphery ¹⁰ and in the brain. ¹ Many in vitro BBB models are based on culture inserts, therefore lack fluid flow. ²⁶ Hollow fiber cartridges, ⁴ microfluidic LOC devices ³ , ³⁴ as well as bioengineered cerebral vessels ³⁵ allow the study of fluid flow on BBB models. A linked organ-on-chip model of the human neurovascular unit revealed the metabolic coupling of endothelial and neuronal cells, ³⁴ while a bioengineered model of human cerebral arteries helped to study the molecular mechanism of cerebrovascular amyloid angiopathy. ³⁵ Flow-mediated regulation of endothelial genes has been studied on vascular ECs for a long time. Upregulation of genes TGF-β, EDN1/ET-1, CCL2/MCP-1 and ICAM-1 has been described previously ³⁶ and we confirmed these data in the present study. The observed changes in our BBB model are complex: while expression level of the vasoconstrictor EDN-1 gene was increased, the level of one of its receptor gene, EDNRB, was decreased, and the vasodilating factor NOSIP, which interacts with NOS3, was found to be significantly upregulated. VEGFR1, upregulated by flow in our model, mediates endothelial cytoprotection via serine/threonine-specific protein kinase AKT/protein kinase B. ³⁷ LYVE-1 is not expressed in mouse brain capillary ECs. ⁸ In static conditions BLECs expressed LYVE1 gene at low level, which was reduced by 13-fold to a very low expression level in dynamic conditions indicating that BBB-like properties were induced in the LOC device.

Among EC surface adhesion molecules flow elevated the expression of ICAM1 gene, as described in the case of vascular ECs ³⁶ and a human BBB model ⁴ which can be important for immune cell transmigration studies. ³² The gene expression level of ICAM1 after upregulation in the dynamic setup was still low and at a similar level observed in mouse brain capillary ECs. ⁸ We found no change in the gene expression for adhesion molecules ICAM2, E and P-selectins.

Endothelial cells respond to fluid flow by changing their morphology to a more elongated shape and align with the flow direction. ¹¹ Our morphological findings on the human BBB model in the LOC device are in accordance with these observations, and further supported by the upregulation of cytoskeletal genes ACTN1, known to interact with NOS3, and TAGLN, crosslinking actin filaments and participating in cytoskeletal reorganization. Flow is also known to increase the adhesion of ECs to the basal membrane via integrins, which interact with collagen, laminin, and fibronectin. ²⁷ We found that several integrin-α subunit genes were significantly increased in dynamic conditions, while integrin-β subunits were unchanged or downregulated. In a human EC and astrocyte co-culture model an increase in the gene expression for both α and β integrin subunits was observed after flow, ⁴ which may indicate a special role for astrocytes in the induction of these genes. In our study the differentially expressed basal membrane genes COL4A1, FBLN5 and the fibronectin related FNDC3B together with integrins may indicate a stronger attachment of the cells in dynamic condition. MMPs are important in different physiological and pathological processes at the BBB. In cerebral ischemia MMPs participate in both the vascular injury and the repair phase during angiogenesis and reestablishment of blood flow. ³⁸ In our study several MMPs were differentially expressed, and this upregulation may be related with basal membrane remodeling induced by flow.

The paracellular tightness of the BBB is controlled by transmembrane TJ proteins. We found the downregulation of epithelial claudins CLDN1, -3, -7 by flow. These claudins were also expressed at a low level in BBB culture models ²⁶ and in isolated brain ECs. ⁷ , ⁸ Flow did not change the expression level of important TJ genes ESAM, OCLDN, JAMs, MARVELDs and linkers TJPs. The expression level of claudin-5 gene, considered as the most important claudin in brain ECs was decreased by flow in our co-culture model with PCs, which was not observed in the astrocyte co-culture model. ⁴ In contrast, the genes of adherens junction proteins were either unchanged or upregulated together with gap junction genes. These, together with the unchanged level of several TJ genes and increased level of basal membrane protein and integrin gene expression might explain that we observed an increased barrier integrity of the BBB model by functional measurements, namely TEER and permeability for marker molecules.

The chemical protection of the CNS is maintained by active efflux pumps, mainly ABC transporters, the gene expression level of which was mostly unchanged in dynamic condition. Enzymes participating in drug metabolism were more sensitive to the effect of flow, and their level was unchanged or upregulated like in the case of phase I enzymes CYP1A1 and CYP1B1, and phase II enzymes EPHX1 and GSTP1. The upregulation of P450 enzymes in a human BBB model by fluid flow was also described. ⁴ The increased gene expression level of phase II enzymes EPHX1 and GSTP1 in brain ECs was also observed by co-culture with PCs. ²³

SLC transporters are key for the proper transport of nutrients to the CNS. ⁹ The expression level of SLC transporters was either unchanged or in many cases increased in dynamic conditions. In addition to nutrients, the BBB regulates the transport of ions by SLCs. On the human BBB model flow increased the gene expression of anion exchanger AE2 and the K⁺- Cl^– cotransporter (KCC1). While the level of voltage dependent Ca²⁺ (CACN family) and anion (VDAC) channels was unchanged by flow, several K⁺ channels were found downregulated, which were also expressed at a very low level in isolated brain capillary ECs. ⁸ In the human BBB co-culture model with astrocytes an upregulation of Ca²⁺ and K⁺ channels was observed in dynamic conditions ⁴ indicating an important role for the co-culture conditions.

Many important proteins and peptides, like insulin and holotransferrin, cross the BBB by receptor-mediated pathway. ¹ Flow did not change the expression level of many of these receptors. One of the exception is the gene of LRP5, a canonical WNT pathway signaling co-receptor, which was increased. Importantly, LRP5 participates in barrier genesis and BBB maturation. ¹ Flow also downregulated in the human BBB model the expression level of caveolin-1 gene. Since isolated brain ECs express much less CAV1 than peripheral (lung) ECs, ⁸ these changes may point to barrier maturation in the present model. We confirmed the presence of several BBB related proteins in the BLEC model by immunostaining. In concordance with the gene expression data, no changes were found between the static and dynamic conditions for junctional proteins β-catenin and ZO-1, the efflux pump PGP, the glucose transporter GLUT1, and the adhesion molecule VCAM-1. We found lower gene expression levels but no significant change for immunostaining intensity of claudin-5 after fluid flow. The gene expression of ICAM-1 adhesion molecule was elevated by flow, which was corroborated with immunostaining.

The study of Cucullo et al. described transcriptomic changes in BBB-related genes induced by flow on a human BBB model, ⁴ but no data are available on the gene expression of mechano-sensing and luminal glycocalyx components. Continuous blood flow regulates the composition of ESG, the dynamic equilibrium of glycoproteins, proteoglycans, glycosaminoglycans and the associated plasma proteins. ²⁷ ESG is not only important as an element of the physical barrier ¹⁵ with its highly negative charge but also control the stability of endothelial cells. We observed the differential upregulation of major glycocalyx core protein genes DCN, SDC2 and VCAN, as well as an increase in GPC1 expression, which act as mechano-sensors and also link the ESG to the cytoskeleton, therefore actively participate in flow induced morphological changes in ECs. ¹⁰ The increased gene expression level of glycocalyx-related enzymes heparanase and hyaluronidase-2, as well as carbohydrate sulfotransferase-1 and the heparan sulfate 6-O-endosulfatase SULF2 adding and removing sulfate groups, which are important for the overall negative charge of the glycocalyx, may indicate active remodeling of the ESG in dynamic condition. Direct measurement of the surface charge of the ECs by a laser Doppler velocimetry confirmed that cells became more negatively charged. These data together with increased WGA lectin staining support that the ESG of the human BBB model became more robust. These observations are further corroborated that the three most enriched pathways in the human BBB model after flow were related to extracellular matrix organization, extracellular matrix structure organization and collagen-containing extracellular matrix. In addition, several other pathways related to cell surface, extracellular matrix proteoglycan and dermatan sulfate biosynthesis and metabolism were also increased.

In conclusion, flow increased barrier properties, induced several key general endothelial and BBB-related genes on BLECs in the LOC device. In addition, flow not only upregulated extracellular matrix and glycocalyx-related genes and pathways, but made the brain endothelial cell surface more negatively charged and more rich in lectin binding sites. These results strongly argue for the inclusion of flow in BBB models and draw the attention to the importance of the endothelial surface glycocalyx as an element of the BBB. This human model can be used as a tool to study the role of cell surface glycocalyx in BBB physiology and pathology.

Supplemental Material

sj-pdf-1-jcb-10.1177_0271678X21992638 - Supplemental material for Flow induces barrier and glycocalyx-related genes and negative surface charge in a lab-on-a-chip human blood-brain barrier model

Supplemental material, sj-pdf-1-jcb-10.1177_0271678X21992638 for Flow induces barrier and glycocalyx-related genes and negative surface charge in a lab-on-a-chip human blood-brain barrier model by Ana R Santa-Maria, Fruzsina R Walter, Ricardo Figueiredo, András Kincses, Judit P Vigh, Marjolein Heymans, Maxime Culot, Peter Winter, Fabien Gosselet, András Dér and Mária A Deli in Journal of Cerebral Blood Flow & Metabolism