Aorta-on-a-chip reveals impaired mitochondrial dynamics as a therapeutic target for aortic aneurysm in bicuspid aortic valve disease

Background: Bicuspid aortic valve (BAV) is the most common congenital cardiovascular disease in general population and is frequently associated with the development of thoracic aortic aneurysm (TAA). There is no effective strategy to intervene with TAA progression due to an incomplete understanding of the pathogenesis. Insufficiency of NOTCH1 expression is highly related to BAV-TAA, but the underlying mechanism remains to be clarified. Methods: A comparative proteomics analysis was used to explore the biological differences between non-diseased and BAV-TAA aortic tissues. A microfluidics-based aorta-on-a-chip model was constructed to evaluate the effect of NOTCH1 deficiency on contractile phenotype and mitochondrial dynamics of human aortic smooth muscle cells (HAoSMCs). Results: Protein analyses of human aortic tissues showed the insufficient expression of NOTCH1 and impaired mitochondrial dynamics in BAV-TAA. HAoSMCs with NOTCH1-knockdown exhibited reduced contractile phenotype and were accompanied by attenuated mitochondrial fusion. Furthermore, we identified that mitochondrial fusion activators (leflunomide and teriflunomide) or mitochondrial fission inhibitor (Mdivi-1) partially rescued the disorders of mitochondrial dynamics in HAoSMCs derived from BAV-TAA patients. Conclusions: The aorta-on-a-chip model simulates the human pathophysiological parameters of aorta biomechanics and provides a platform for molecular mechanism studies of aortic disease and related drug screening. This aorta-on-a-chip model and human tissue proteomic analysis revealed that impaired mitochondrial dynamics could be a potential therapeutic target for BAV-TAA. Funding: National Key R&D Program of China, National Natural Science Foundation of China, Shanghai Municipal Science and Technology Major Project, Shanghai Science and Technology Commission, and Shanghai Municipal Education Commission.


40
Bicuspid aortic valve (BAV) disease is the most common congenital cardiovascular abnormality and is found in nearly 1.4% of the general population (Garg et al., 2005;Michelena et al., 2011; The relationship between NOTCH1 insufficiency and mitochondrial dysfunction in human Figure 1. The relationship between NOTCH1 insufficiency and mitochondrial dysfunction in BAV-

6
The online version of this article includes the following source data for figure 1: 155 Source data 1. The original raw data for western blotting.
Source data 2. List of total differential protein expression between non-diseased and BAV-TAA aortic tissues.

219
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021. ; https://doi.org/10.1101/2021.07.07.21260132 doi: medRxiv preprint  experiments. An approximately 60% reduction in NOTCH1 mRNA expression was found in the 245 NOTCH1-KD group compared with the WT and NC groups (Figure 4-figure supplement 1a).

246
Under static conditions, the mRNA expression of SM22 and CNN1 was upregulated and that of 247 OPN was downregulated in the NOTCH1-KD group compared with the WT and NC groups 248 (Figure 4-figure supplement 1a). The same tendency was found for the protein expression 249 levels by western blotting analysis (Figure 4-figure supplement 1b-c).

250
NOTCH1-KD and WT HAoSMCs were cultured on aorta-on-chip models under rhythmic 251 strain and static conditions to characterize the expression of SM22 and CNN1. As shown in 252 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted July 9, 2021. ; https://doi.org/10.1101/2021.07.07.21260132 doi: medRxiv preprint Figure 4b-c, immunofluorescent staining images showed that SM22 and CNN1 were upregulated in the NOTCH-KD group compared with the WT group under static conditions. However, the 254 opposite results were observed under rhythmic strain conditions: the expression of SM22 and 255 CNN1 in NOTCH-KD HAoSMCs was lower than that in the WT group under rhythmic strain 256 conditions. Western blotting analyses revealed similar alterations in the expression of SM22 and 257 CNN1 (Figure 4d-e). In WT HAoSMCs, rhythmic strain induced the upregulation of SM22 and 258 CNN1 expression compared with the levels observed under static conditions. However, a 259 downregulation of SM22 and CNN1 expression was detected in NOTCH-KD HAoSMCs exposed 260 to rhythmic strain. These results suggested that rhythmic strain induced different effects between 261 WT and NOTCH-KD HAoSMCs. In NOTCH-KD HAoSMCs, SM22 and CNN1 expression were

270
Representative images of the western blotting analyses of SM22 and CNN1 in the WT and NOTCH1-KD 271 groups exposed to rhythmic low or high strains for 24 hrs. (e) Quantification of the total band densities for 272 individual proteins normalized to the corresponding band density of β-actin (n = 4, **P < 0.01, ***P < 0.001, 273 ****P < 0.0001, two-way ANOVA followed by Tukey's post hoc test). All the data are expressed as the 274 means ± SDs.

275
The online version of this article includes the following figure supplement(s) and source data for figure 4:

280
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021. ; https://doi.org/10.1101/2021.07.07.21260132 doi: medRxiv preprint Furthermore, we assessed the alterations in mitochondrial function and dynamics in 281 NOTCH1-KD HAoSMCs. Specifically, western blotting was performed to detect the protein expression of MFN1 and MFN2 (mitochondrial fusion-related proteins) and DRP1 and MFF (mitochondrial fission-related proteins). No significant difference in MFN1, MFN2, DRP1 and MFF 284 protein expression was found between the WT and NC groups under static conditions, which 285 suggested that the expression of these mitochondrial dynamics-related proteins was not affected 286 by the lentivirus vector (Figure 4-figure supplement 1d-e). As shown in Figure 5a

326
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021.  independent biological replicates each with two to three technical replicates were plotted. *P < 0.05, **P < 341 0.01, two-way ANOVA followed by Tukey's post hoc test). (i) The ATP concentrations were measured using 342 an ATP Determination Kit (n=3, data from three independent biological replicates each with two technical 343 replicates were plotted. **P < 0.01, two-way ANOVA followed by Tukey's post hoc test). Quantitative 344 measurements were calculated using ImageJ software. All the data are expressed as the means ± SDs.

345
The online version of this article includes the following source data for figure 5:      . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Drugs rescued the impairment of mitochondrial dynamics in NOTCH1-insufficient
To confirm whether the decreased contractile phenotype of NOTCH1-knockdown HAoSMCs can be rescued by inhibition of mitochondrial fission or activation of mitochondrial fusion, we evaluated the phenotypic alterations and mitochondrial dynamics of NOTCH1-insufficient and p-HAoSMCs purchased from ATCC (PCS-100-012) were transfected with NOTCH1-shRNA 363 to generate NOTCH1-knockdown HAoSMCs.

364
The effect of drug testing on aorta-on-a-chip model was desirable, even though the results 365 did not exhibit all the positive responses. In the CRL1999 HAoSMCs (Figure 6b) Table V. Significantly lower NOTCH1 protein expression was 383 found in p-HAoSMCs from BAV-TAA aortic tissues compared with p-HAoSMCs from non-384 diseased aortic tissue (Figure 6-figure supplement 1). A significant reduction in MFN1 and 385 MFN2 expression was found in p-HAoSMCs from BAV-TAA aortic tissues under static and 386 rhythmic high-strain conditions. We also found that SM22 and CNN1 expression was significantly Figure 6. Screening of drugs that can rescue the cell phenotype and mitochondrial dynamics using loose of vascular tone, which is one of the major pathological process in TAA. In this study, we found that NOTCH1 insufficiency in HAoSMCs downregulated the contractile phenotype proteins 480 SM22 and CNN1 under rhythmic strain. Also, the expressions of SM22 and CNN1 were decreased in BAV-TAA aortic tissues compared to non-diseased aortic tissues. Noseda et al.

482
reported that NOTCH1 activation was needed for expression of the contractile phenotype protein 483 αSMA in VSMCs via the NOTCH/CSL axis (Noseda et al., 2006). High et al. reported that 484 NOTCH1 promotes the differentiation of the cardiac neural crest into differentiated VSMCs and 485 the expression of αSMA (High et al., 2007). However, the specific mechanism by which NOTCH1 486 leads to a reduced contractile phenotype has not been fully investigated. Based on results

487
obtained from the aorta-on-a-chip model, it comes to the assumption that NOTCH1-insufficiency 488 reduced the contractile phenotype through decreasing mitochondrial fusion.

489
Mitochondrial dynamics itself has been recognized as one of the critical factors that regulate

524
The experimental data obtained in this study encourage further studies on the application of 525 leflunomide and teriflunomide in BAV-TAA.

526
In conclusion, we constructed a aorta-on-a-chip model, which could serve as a 527 complementary tool to the current cell culture system and animal models. Using the aorta-on-a-528 chip model, we found that NOTCH1 insufficiency in HAoSMCs induced phenotypic switching from 529 a contractile to a synthetic phenotype accompanied by an impairment of mitochondrial fusion, 530 implying its potential role as a therapeutic target for BAV-TAA. At the current stage, in vitro 531 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021. ; https://doi.org/10.1101/2021.07.07.21260132 doi: medRxiv preprint microphysiological models and animal models are complementary to each other in the sense that 532 they are able to provide more comprehensive basis for preclinical assays with greater predictive 533 power.
The structure of the three-layer microfluidic aorta-on-a-chip model was designed using computer-aneurysm (BAV-TAA; mean age: 59.3 years; range: 43-72 years; 4 males). The patients' basic information is available in Supplementary File- Table I. glucose Dulbecco's modified Eagle's medium (DMEM, Gibco) with 20% fetal bovine serum (FBS, 590 Gibco) and 1% penicillin and streptomycin (p/s, Gibco) for 2-3 weeks at 37°C and 5% CO2 in a 591 humidified incubator. After approximately 10-12 days, the p-HASMCs started to migrate out of the 592 tissue pieces. When the cells reached approximately 80% confluency, first-passaged cells were 593 rinsed with PBS, digested using 0.25% trypsin (Gibco), and replated in smooth muscle cell culture conditions. The nuclei were counterstained with 4′,6-diamidino-2-phenyllindole (DAPI) (Thermo Fisher Scientific) for 10 min. The aorta-on-a-chip models were then disassembled, and images were acquired with a fluorescence microscope (Leica DMi8) and analyzed using ImageJ software. incubated for 10 min at room temperature. After incubation, the ATP concentration was measured using a luminometer. An ATP concentration standard curve was then established and used to calculate the ATP concentration of each sample.

686
The aortic samples were minced and lysed with RIPA on ice for 30 min. The extracts were 687 centrifuged at 14,000 rpm and 4°C for 25 min, and the supernatant was collected after 688 centrifugation. The protein concentrations were quantified using a BCA Protein Assay kit (Thermo 689 Fisher Scientific). Filter-aided sample preparation (FASP) was performed for protein digestion.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021. ; https://doi.org/10.1101/2021.07.07.21260132 doi: medRxiv preprint The experimental results are expressed as the means ± standard deviations (SDs). A minimum of three individual replications of each group were used for the relative analyses. The statistical significance was indicated by *P < 0.05, **P < 0.01, ***P < 0.001 and ****P < 0.0001.

The Supplementary File include the following:
Online Tables I-VI  measurements  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021. ; https://doi.org/10.1101/2021.07.07.21260132 doi: medRxiv preprint  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021. ; https://doi.org/10.1101/2021.07.07.21260132 doi: medRxiv preprint (c) Length-to-width ratio of HAoSMCs after exposure to low or high rhythmic strain for 24 hrs. (n = 3, cells were measured in three fields per sample, ****P < 0.0001, one-way ANOVA followed by Tukey's post hoc test). (d) Alignments of HAoSMCs exposed to low or high rhythmic strain for 24 hrs. (e) Representative images of the western blotting analyses of protein markers of the contractile phenotype (SM22 and CNN1) and synthetic phenotype (OPN) of HAoSMCs after exposure to low or high rhythmic strain for 24 hrs. (f) Quantification of the total band densities for individual proteins normalized to the corresponding band of β-actin (n = 4, *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, one-way ANOVA followed by Tukey's post hoc test). All the data are expressed as the means ± SDs.  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021.  Representative images of immunofluorescence staining of SM22 and CNN1 after exposure to rhythmic low or high strain for 24 hrs. The scale bar represents 200 μm. (c) Intensity of immunofluorescence staining of SM22 and CNN1 (*P < 0.05, **P < 0.01, ****P < 0.0001, two-way ANOVA followed by Tukey's post hoc test). (d) Representative images of the western blotting analyses of SM22 and CNN1 in the WT and NOTCH1-KD groups exposed to rhythmic low or high strains for 24 hrs. (e) Quantification of the total band densities for individual proteins normalized to the corresponding band density of β-actin (n = 4, **P < 0.01, ***P < 0.001, ****P < 0.0001, two-way ANOVA followed by Tukey's post hoc test). All the data are expressed as the means ± SDs.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021. ; Quantification of the total band densities for four individual proteins normalized to the corresponding band density of β-actin (n = 4, *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, two-way ANOVA followed by Tukey's post hoc test). (c) MitoTracker staining images of the mitochondrial morphologies in the WT and NOTCH1-KD groups under rhythmic low or high strain or static conditions. The scale bar represents 50 μm. (d) Quantification of the mitochondria length (n = 3, **P < 0.01, ****P < 0.0001, two-way ANOVA followed by Tukey's post hoc test). (e) TMRM staining of the mitochondrial membrane potentials in the WT and NOTCH1-KD groups under rhythmic low or high strain or static conditions. The scale bar represents 100 μm. (f) Quantification of the relative TMRM fluorescence intensity (**P < 0.01, two-way ANOVA followed by Tukey's post hoc test). (g) MitoSOX staining of mitochondrial superoxide generation in the WT and NOTCH1-KD groups. The scale bar represents 100 μm. (h) Quantification of the relative MitoSOX fluorescence intensity (*P < 0.05, **P < 0.01, two-way ANOVA followed by Tukey's post hoc test). (i) The ATP concentrations were measured using an ATP Determination Kit (**P < 0.01, two-way ANOVA followed by Tukey's post hoc test). Quantitative measurements were calculated using ImageJ software. All the data are expressed as the means ± SDs.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021.  SM22 and CNN1 were normalized to the corresponding band densities of β-actin. All the data are expressed as the means ± SDs. n = 3, *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, one-way ANOVA followed by Tukey's post hoc test.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 9, 2021. ; https://doi.org/10.1101/2021.07.07.21260132 doi: medRxiv preprint . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted July 9, 2021.  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted July 9, 2021. ; https://doi.org/10.1101/2021.07.07.21260132 doi: medRxiv preprint