1. Introduction
The low-density lipoprotein (LDL) receptor has a significant function in human lipid metabolism [
1]. Numerous mutations within the LDLR gene are pathogenic, causing hereditary familial hypercholesterolemia (FH). Patients with FH have elevated blood cholesterol levels and concomitant lipid metabolism diseases, such as atherosclerosis, coronary heart disease, and Alzheimer’s disease.
The main role in the FH pathogenesis is assigned to the LDL receptors in liver cells, which are responsible for cholesterol uptake and metabolization. However, when the LDL receptors are defective or absent, hepatocytes have impaired low-density lipoprotein absorption which ultimately leads to increased cholesterol levels in the bloodstream. Elevated blood cholesterol indirectly affects blood vessels and vascular cells through the lipid accumulation inside the intima, causing oxidative stress, inflammatory cell recruitment and local cytokine production [
2,
3,
4]. Nonetheless, it is still little known about how vascular endothelial cells themselves differ between individuals with normal and pathogenic LDLR alleles.
The LDL receptors are involved in the functioning of variety of cell types, including pluripotent progenitors, vascular smooth muscle cells, renal mesangial cells, renal tubular cells, podocytes, brain endothelial cells and umbilical artery [
5]. It is believed that сellular metabolism and metabolites are crucial in regulating three-dimensional chromatin architecture and gene expression [
6]. Therefore, we hypothesize that endothelial cells with dysfunctional LDLR and associated metabolic abnormalities in cellular metabolism might themselves have dysregulated gene expression profiles independent of excess extracellular cholesterol levels. Our study focused on identifying both the LDLR protein expression and transcriptomic patterns in endothelial cells with pathogenic allelic variants in the LDLR gene and discussed how this may potentially contribute to the FH progression. Previously, we obtained two lines of induced pluripotent stem cells (iPSCs) from patients who were compound heterozygotes for pathogenic and likely pathogenic allelic variants of the LDLR gene [
7,
8,
9]. Both lines contained alterations within the LDLR gene that lead to impaired protein maturation, reducing the rate of LDLR transport to the membrane, and causing its partial retention in the endoplasmic reticulum. We obtained endothelial derivatives by direct differentiation of iPSC lines from both a conditionally healthy donor and patients with FH. We showed that typical endothelial derivatives of human iPSCs predominantly produce mature LDLR protein, whereas its level is significantly reduced in endothelial derivatives carrying pathogenic LDLR alleles. Transcriptomic analysis of FH patient-specific iPSC-derived endothelial cells (iPSC-ECs) revealed bias in gene expression in several signaling pathways. Cellular metabolic abnormalities caused by defective LDLR are manifested at the transcriptomic level as down-regulated expression of monocarboxylic acid transporters. Altered expression was also found for genes involved in exocytosis, cell-cell adhesion, cell secretion and leukocyte activation. These data suggest that endothelial cells with damaged LDLR are themselves more predisposed to dysfunction, oxidative stress, and chronic inflammation, which may be responsible for the earlier and accelerated progression of atherosclerosis and cardiovascular disease in patients with FH.
3. Discussion
The current trend in modeling of human diseases is to collect patient-specific iPSCs from individuals with hereditary pathologies and to elucidate the mechanisms of this pathology at the molecular and cellular levels using relevant cells obtained by iPSC differentiation [
23]. Patient-specific iPSC-ECs are broadly used to model various disorders, including pulmonary arterial hypertension, Moyamoya disease, fibrodysplasia ossificans progressiva, Huntington’s disease, Kawasaki’s disease, type I diabetes mellitus, atrial or ventricular septal defects, pulmonary valve stenosis, cardiomyopathy, calcified aortic valve disease, hemophilia A, diabetic endotheliopathy, Hutchinson-Gilford Progeria Syndrome, cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy, and peripheral artery disease [
27,
28,
29,
30,
31,
32].
There are known lines of iPSCs from patients with FH and examples of their use to model of this pathology, which have only been carried out on hepatocyte-like derivatives of iPSCs [
20,
21,
22,
33,
34]. Despite the knowledge that endothelial cells are the main target in atherosclerosis and the progression of cardiovascular pathology in FH, ECs or iPSC-derived ECs from patients with FH have not yet been obtained or studied. Current studies of endothelial dysfunction in FH patients have been based on blood biomarkers, blood flow and microvascular function indicators [
35,
36,
37].
In endothelial cells, LDLR is thought to be involved in the endocytosis and transcytosis of circulating LDLs [
38,
39,
40]. LDLR is found expressed in the blood-brain barrier endothelium [
35] and in human umbilical artery endothelial cells (HUAECs) [
21]. However, we found no information on the LDLR expression in endothelial cells from FH patients, including iPSC-ECs. With this work, we fill this gap and hope that studying patient-specific iPSC-ECs may potentially contribute to the understanding of the FH pathogenesis and aid in the development of new treatment approaches.
In our study, we were the first to obtain iPSC-ECs from compound heterozygous patients with FH and demonstrated dysfunction of the LDLR protein. We have shown that both iPSCs and iPSC-ECs with pathogenic LDLR alleles, have significantly reduced levels of the mature form of LDLR and an increased level of the immature form compared to controls with normal LDLR.
The mutation in one LDLR allele of the FH 1.3.1 patient is c. 530С>T, leading to the substitution p.Ser177Leu (ClinVar ID 3686), which is considered pathogenic and classified as a transport-defective mutation [
41,
42]. The mutation in the other LDLR allele of the patient is c.1054T>C, resulting in p.Cys352Arg (ClinVar ID 251618), which is considered likely pathogenic and associated with a disruption in the EGF-like domain of the LDLR protein [
38]. The second mutation can also be classified as transport-defective mutation [
39].
In patient FH 3.2.8, an extended deletion c.2141-966_2390-330del spans introns 14-16 and exons 15-16 and, according to bioinformatic predictions, is considered pathogenic or disease causing as it leads to the loss of the O-glycosylation domain, which interferes with the production of the mature LDLR protein. The second mutation c.1327T>C leads to p.Trp443Arg (ID 998052), referred as likely pathogenic, causes a recycling defect when ligands are not released from the complex with LDLR in endosomes and, as a result, LDLR does not recirculates to the cell surface. Our results on LDLR expression allow us to conclude that all LDLR mutations in the iPSCs and iPSC-EC lines lead to the reduction in the mature form of the LDLR protein, while the combination of mutations in the second patient significantly reduces both mature and immature LDLR forms.
Pharmacological blocking of endoplasmic reticulum-associated degradation (ERAD) or using pharmacological chaperones are able to provide transport of immature LDLR to the cell membrane surface, restoring receptor function despite its incomplete maturation [
36]. Thus, iPSC-derived ECs obtained from patients with FH may potentially help in the development and screening of targeted drugs to treat the consequences of specific mutations in patients with FH.
In this study, we also compared iPSC-ECs with normal and defective LDLR at the transcriptome level. We found no difference in LDLR gene expression between FH and CTRL iPSC-ECs, suggesting that LDLR function is mainly abolished at the post-translational level, as also shown by immunoblotting results.
Transcriptome profiling revealed the downregulation of the most genes (67%) in LDLR mutant iPSC-EC. We found that biological processes related to monocarboxylic acid transport, exocytosis and cell-cell adhesion were downregulated in LDLR mutant iPSC-EC lines. On the contrary, cell secretion and leukocyte activation were upregulated in these cells. Thus, the quantitative reduction of mature LDLR in iPSC-ECs abolished endothelial function, which was manifested both in the decreased ability to uptake LDL and at the transcriptomic level.
The down-regulated expression of monocarboxylic acid transporters likely reflects the expected dysregulation of cellular metabolism in iPSC-ECs with defective LDLR at the transcriptome level. Quantitative variations in metabolites associated with these transporters, including Ac-KoA and fatty acids (organic anions), may affect chromatin and gene regulation, while decreased levels of antioxidant transport could potentially render cells more susceptible to oxidative stress [
45,
46,
47].
There was also dysregulation of angiogenesis genes (XDH, SEMA5A, and SULF1) and the Rap1 signalling pathway (ADORA2B, GNAO1, SULF1, EFNA2), which promotes endothelial homeostasis and may be involved in endothelial dysfunction-associated cardiovascular pathologies [
43]. However, the iPSC-ECs with normal and defective LDLR showed no differences in the ability to form capillary-like structures in an angiogenesis test, both in terms of total vessel length and in the number of branching points. Therefore, it is likely that the identified dysregulation in angiogenesis genes may be related to other aspects of this process. It is also noteworthy that Rap1 signalling and the XDH gene are both related to oxidative stress [
49,
50].
Additionally, the reconstructed protein-protein interaction network showed an association between LDLR and some molecules such as SYT11, CNTN4, and SEMA5A. Interestingly, semaphorin and ephrin family members are involved in axon guidance and synaptic plasticity and they are also important in endothelial cell-leukocyte interactions during atherogenesis [
46]. Synaptotagmin-11 (SYT11) and contactin 4 (CNTN4) are both essential for neural development [52,53]. SYT11 is a vesicle-trafficking protein that can suppress microglial activation by inhibiting cytokine secretion and phagocytosis [
49]. CNTN4, a crucial adhesion protein, is involved in T-cell activation and oxLDL-induced cell apoptosis and inflammation in THP-1 macrophages [55,56]. Thus, our data suggest that dysregulation of neuronal guidance molecules (chemoattraction and chemorepulsion) can be considered as a possible mechanism for endothelial dysfunction in LDLR mutant iPSC-EC lines.
Overall, the transcriptomic analysis of iPSC-ECs with defective LDLR identified a number of genes and signaling pathways that could potentially be important for the progression of endothelial dysfunction in FH patients. Endothelial cells with damaged LDLR per se appear to be predisposed to dysfunction, oxidative stress and chronic inflammation, which may facilitate early and accelerated FH progression. However, further detailed studies are needed to confirm and understand our findings.
Author Contributions
Conception and design: ISZ, AIS, AAS, MSN. Preparation of experiments, analysis and interpretation of data: ISZ, AIS, AMA, AAS, MSN, AAZ, SBS, NAT, VAS, NVZ. Contribution of reagents/materials/analysis tools/administrative and financial support: YVP, VVK, MVE, SMZ. Manuscript writing: ISZ, AIS, AAS; MSN. All authors read and approved the final manuscript.