Serotonin Potentiates Transforming Growth Factor-beta3 Induced Biomechanical Remodeling in Avian Embryonic Atrioventricular Valves

Embryonic heart valve primordia (cushions) maintain unidirectional blood flow during development despite an increasingly demanding mechanical environment. Recent studies demonstrate that atrioventricular (AV) cushions stiffen over gestation, but the molecular mechanisms of this process are unknown. Transforming growth factor-beta (TGFβ) and serotonin (5-HT) signaling modulate tissue biomechanics of postnatal valves, but less is known of their role in the biomechanical remodeling of embryonic valves. In this study, we demonstrate that exogenous TGFβ3 increases AV cushion biomechanical stiffness and residual stress, but paradoxically reduces matrix compaction. We then show that TGFβ3 induces contractile gene expression (RhoA, aSMA) and extracellular matrix expression (col1α2) in cushion mesenchyme, while simultaneously stimulating a two-fold increase in proliferation. Local compaction increased due to an elevated contractile phenotype, but global compaction appeared reduced due to proliferation and ECM synthesis. Blockade of TGFβ type I receptors via SB431542 inhibited the TGFβ3 effects. We next showed that exogenous 5-HT does not influence cushion stiffness by itself, but synergistically increases cushion stiffness with TGFβ3 co-treatment. 5-HT increased TGFβ3 gene expression and also potentiated TGFβ3 induced gene expression in a dose-dependent manner. Blockade of the 5HT2b receptor, but not 5-HT2a receptor or serotonin transporter (SERT), resulted in complete cessation of TGFβ3 induced mechanical strengthening. Finally, systemic 5-HT administration in ovo induced cushion remodeling related defects, including thinned/atretic AV valves, ventricular septal defects, and outflow rotation defects. Elevated 5-HT in ovo resulted in elevated remodeling gene expression and increased TGFβ signaling activity, supporting our ex-vivo findings. Collectively, these results highlight TGFβ/5-HT signaling as a potent mechanism for control of biomechanical remodeling of AV cushions during development.


Introduction
Biomechanical remodeling is the process by which living tissues reorganize, reshape, and refit their microstructure in adaptation to changing internal and external forces. This process defines much of embryogenesis, during which initially indistinct cellular masses acquire shape and functional specificity through production and manipulation of the extracellular matrix (ECM). This is particularly important for the morphogenesis of the heart, which is critically responsible for distributing nutrients as the embryo grows. The heart transitions rapidly from a tubular structure into a multi-chambered pumping organ, simultaneously growing over 100-fold in volume [1]. The hemodynamic environment inside the heart increases dramatically in severity during this process [2][3][4], which means the biomechanical properties of the forming valves must be precisely tuned to maintain efficient unidirectional blood flow. Atrioventricular (AV) valve morphogenesis is characterized by rapid ECM accretion and turnover [5,6], which is hypothesized to be stimulated by a dynamic interaction of molecular and mechanical signaling. While numerous molecular agents important for valve morphogenesis have been identified [7][8][9][10], less is known about how these signals affect valve mechanics, which is a key readout of valve function.
The transforming growth factor-beta (TGFb) superfamily is critically important for a wide range of cellular processes [11][12][13], and is heavily involved in directing morphogenesis of AV cushions [14][15][16][17][18]. In the chick, TGFb2 and TGFb3 isoforms are necessary for the endothelial to mesenchymal transition (EMT) which initiates AV cushion development [19]. TGFb2 induces initial cellcell separation of valve endothelial cells, while TGFb3 stimulates their invasion and subsequent mesenchymal phenotype shift [15,16]. During post-EMT, these mesenchymal cells facilitate a transition in the cushion microstructure from glycosaminoglycans (GAGs) (hyaluronan, versican) toward fibrous structural proteins (collagen I, IV, V, fibronectin, periostin) [5,20,21]. This shift in ECM content translates into increased valve stiffness [22], and coincides with elevated expression of TGFb3 in the cushions and AV canal [23]. Furthermore, TGFb3 upregulates collagen I and periostin in post-EMT AV cushion explants [24], suggesting that TGFb3 is a key modulator of cushion ECM content, and consequent mechanical properties. An aim of this study is to better understand this remodeling potential of TGFb3 through a combined analysis of cushion stiffness, matrix compaction, cell proliferation, and ECM synthesis.
The capacity of TGFb3 to stimulate valvular remodeling events underscores the importance of identifying molecular signals which modulate TGFb activity. Recent studies indicate that serotonin (5-HT) interacts with TGFb signaling in adult heart valves [25,26], and can also alter valve mechanical properties [27,28]. 5-HT, which is a monoamine neurotransmitter derived from the essential amino acid tryptophan [29], increased the stiffness of porcine aortic valve cusps with the endothelial layer denuded [27], and under cyclic stretch [30]. Serotonin also increased collagen synthesis in human and sheep valve interstitial cells (VICs) [25,31]. Reports in adult VICs indicate that 5-HT can upregulate TGFb, resulting in cell differentiation and aberrant connective tissue accumulation [25,26,32]. In development, serotonin is active in key events such as cardiac progenitor patterning, left-right laterality, and migration of the neural crest [33][34][35][36][37]. Murine AV cushions express the serotonin receptors 5-HT2a and 5-HT2b, and the serotonin transporter (SERT) by the completion of EMT [38,39], which is when TGFb3 expression increases in the cushions [18,40]. Latent TGFb binding protein and serotonin binding protein are also expressed in murine post-EMT endocardial cushions [41,42], highlighting each pathway's capacity to regulate expression of their ligands. The proximity of these TGFb and 5-HT signaling components suggests that they may be interacting partners in post-EMT cushion development. Furthermore, a recent study reported TGFb1 upregulation in murine SERT KO hearts at near fetal stages, which was hypothesized to be a consequence of excess 5-HT signaling due to SERT inhibition [43]. In light of these signaling interactions in both adult and development models, we hypothesize that this mechanically relevant crosstalk of TGFb and 5-HT may play a role in modulating embryonic AV cushion biomechanics.
The objectives of this study therefore were to characterize the remodeling capacity of TGFb3 in AV cushions, and determine how TGFb3 and 5-HT may act together to regulate cushion biomechanical remodeling. Chick AV cushion biomechanics, compaction, and candidate gene expression were quantified through implementation of an ex vivo cushion culture system. We determined that TGFb3 induces AV valve stiffening through increases in cell proliferation, myofibroblastic differentiation, and collagen synthesis. 5-HT enhances the AV valve stiffening effect of TGFb3 in a dose-dependent manner. Crosstalk between TGFb3 and 5-HT signaling was investigated via molecular inhibition studies. The ex vivo results were then tested in ovo through an elevated 5-HT model. These results suggest that 5-HT may be an important potentiator of TGFb3 signaling in embryonic valve morphogenesis and biomechanical stiffening.

Ethics Statement
Leghorn avian embryos from Hamburger-Hamilton stages (HH) 17-36 were utilized in this research. All procedures in this study followed the guidelines of Cornell University and NIH policy, which state that avian embryos of these stages are not considered vertebrate animals for the purposes of IACUC regulation.

Micromechanical testing
Cushion mechanical properties were measured after 24 hour treatment in the ex vivo study and at HH25 in the in ovo study using the micromechanical pipette aspiration technique [22,46,47]. A glass micropipette (,70-100 mm in diameter) was placed adjacent to the cushion surface, and a small vacuum pressure was incrementally applied. The pressure source was a 200 mL pipetter calibrated with a custom manometer. Previous strain history was mitigated by preconditioning with ,20 cycles of low pressurization (,1 Pa). The tissue was then monotonically loaded with increasing static pressure loads, at which images were captured. Aspirated length L, measured as the length from the tip of the pipette to tip of the tissue furthest inside the pipette, was converted into an experimental ''stretch ratio'', l~L zr p r p , by normalizing to the pipette radius, r p . The cushion was assumed to be an isotropic, incompressible, hyperelastic material with an exponential free where I B is the first invariant of left Cauchy Green stretch tensor. AV cushion material isotropy at HH25 was supported by a lack of preferred matrix orientation as determined by ubiquitous protein stain 5-DTAF (50 mM Invitrogen; Figure S1). The DP vs. l data was then fit to the axial stress equation of a uni-axially loaded bar of this exponential material, specifically, s axial~a C exp a l 2 z 2 l {3 h i . From previous analysis [22], the DP vs. l curve differs from the uniaxial load expression by a scale factor, c. This scale factor was numerically determined to be a function of only the material parameter a. Due to the nonlinear nature of the data, the mechanical testing data is presented as strain energy density. This was calculated as the area under the DP vs. l curve fit from l = 1-2 ( Figure 1A), which from our assumed material model is

Compaction & Opening Angle Assays
Compaction of the AV cushions was quantified as the ratio of cross-sectional area before (A 0 ) and after (A) 24 hours of culture in the different treatment conditions, denoted A/A 0 . This ratio measures the combined biomechanical remodeling effects of cell traction, proliferation, and ECM synthesis. To isolate cell traction effects, we quantified the opening angles created by micro-slit incision in AV cushions after 24 hours of treatment. The incision was made along the centerline of the spherical cushion mass, extending approximately one radius into the cushion, and immediately created a pie-wedge with defined opening angle. Opening angles are an established indicator of tissue residual stress [48], which is primarily a function of cell traction forces in our culture system. Images were taken at 1506 magnification using Zeiss Discovery v20 stereomicroscope (Spectra Services, Inc.) and QImaging Retiga 4000R Fast camera (Spectra Services, Inc). Cross-sectional area and opening angle were measured from calibrated images using NIH ImageJ image analysis software.

Immunohistochemistry (IHC)
Proliferation was assessed through bromodeoxyuridine (BrdU) incorporation into HH25 AV cushion hanging drops. BrdU reagent (Invitrogen) was added at 1:100 dilution in culture medium 6 hours prior to completion of 24 hour culture. AV cushions were then rinsed and fixed in 4% paraformeldehyde (PFA). BrdU incorporation was assessed via immunofluorescent antibody staining and confocal microscopy using anti-BrdU 488 (1:100, Invitrogen), with DRAQ5 (1:1000, Biostatus) as a DNA counterstain. Images were processed via ImageJ, and BrdU incorporation was quantified as the ratio of BrdU positive cells to total cell count. IHC was also used to label phosphorylated Smad2/3 (pSmad2/3) complex in HH25 cushions isolated from the systemic 5-HT in ovo model. Isolated cushions were fixed in 4% PFA and then stained via standard whole mount IHC protocol. The cushions were stained with primary pSmad2/3 polyclonal goat anti-human antibody (1:50, Santa Cruz) followed with 488 fluorescent secondary (1:100, Santa Cruz) and cell nuclei counter stain DRAQ5 (1:1000). pSmad2/3 was quantified as the number of cell nuclei with localized pSmad2/3 divided by the total number of cell nuclei.

PCR quantification of gene expression
At the end of 24 hour treatment, AV cushion mRNA was isolated and purified using RNEasy Isolation Kit (Qiagen). A set of 8-10 cushions were pooled per test sample. RNA integrity was determined by NanoDrop spectrometry, using A260/A280 ratio between 1.8 and 2.2 as quality control. cDNA synthesis was completed using SuperScript III first strand RT-PCR kit (Invitrogen) with oligo(dT) primers. Amplification reactions were as follows: (95uC 15 s), (54uC 15 s), (72uC 30 s). Power Syber Green (Applied Biosystems) replication indicator was read at the completion of each 72uC stage. Standard curves for all primers (listed in Table S1) were generated from HH34 brain mRNA and normalized to 18 s ribosomal RNA. Threshold cycle count, C(t), was used to calculate gene expression via the DDCt method using 18 s rRNA as a housekeeping reference gene [49].

5-HT administration in ovo
HH17 stage fertilized leghorn chicken eggs were windowed on their blunt side. Up to 1.0 mg of serotonin (Sigma) was diluted into 100 mL of PBS and dispensed directly onto the chorionic membrane at HH17, HH25, or HH31. The max 5-HT dosage was equivalent to 18 mg/kg which is comparable to other elevated 5-HT animal models (25 mg/kg and 75 mg/kg) [50,51]. After 5-HT treatment, chicks were then sealed and cultured at 55% humidity and 38uC until HH36 (Day 10). Preliminary experiments demonstrated that 5-HT treatment sometimes resulted in an ectopic heart, so additional embryos were alternatively subjected to a thoracotomy that mimicked an ectopic heart without serotonin administration as a control. Embryos were then dissected and analyzed for gross anatomical defects. Hearts with intact great artery connections were then removed, cleared, and analyzed with 3D confocal microscopy or serial section histology using Movat's pentachrome stain. Optical fluorescence tomography (OFT) of ventricular, valve, and outflow vessel anatomy was performed as previously described [52,53]. Briefly, HH36 hearts were freshly isolated and rinsed with 1% lidocaine in PBS buffer. Following rinse, hearts were perfused with fluorescein isothiocyanate-poly-L-lysine (Sigma) via micro injection and then fixed in 4% PFA. The poly-L-lysine binds to the negatively charged endothelial glycocalyx. Hearts were then cleared using Murray's (1 ng/ml, TGFb3+) and control media (TGFb32) treated cushions, n = 4. Strain energy density was calculated from the shaded regions beneath the DP vs l curves. Inset: image of aspirated HH25 AV cushion after 24 hours of culture. The pipette radius, r p , and the aspirated length, L are indicated. Scale bar = 70 mm. B) AV cushion strain energy density increased with TGFb3 treatment, but was blocked by Alk5 inhibition (SB431542, 2.6 mM). mean 6 SEM, n$7, *p,0.0001, 2-way ANOVA. doi:10.1371/journal.pone.0042527.g001 Clear, followed by deep tissue 3D imaging via fluorescence confocal microscopy. Hearts were screened for major defects, and valve morphometry were quantified from this using ImageJ. Valve measurements included leaflet length, average thickness, and minimal thickness with control n = 3 and 5-HT treatment n = 6. Average thickness (t avg ) was calculated as t avg = A L /L, where L is the annulus-tip length of the leaflet, and A L is cross-sectional area of leaflet. The location of minimum thickness was generally the same for all specimens regardless of treatment.

Statistical Analyses
All data is presented as mean 6 standard error of the mean for the number of samples reported. Statistical comparisons between groups were performed using ANOVA for data sets involving more than two groups, or two-tailed t test when only two groups were compared. Defect prevalence in the in ovo model was compared using a chi-squared statistical test. In all comparisons, differences between groups was considered statistically significant for p valves smaller than 0.05.

TGFb3 increases AV cushion stiffness
Ex vivo cultured AV cushions exhibited nonlinear mechanical behavior that was well described by the exponential constitutive model ( Figure 1A). Administration of exogenous TGFb3 (1 ng/ml) increased cushion stiffness 2.5 fold over controls (W TGFb3 = 0.96560.051 vs. W Contr = 0.37860.021, p,0.0001 Figure 1B). Inhibition of canonical TGFb signaling via the TGFb type 1 receptor Alk 5 (2.6 mM SB431542 [54]) blocked the increase in cushion stiffness (W T+TI = 0.24560.043 Figure 1B). The Alk5 inhibitor alone had no effect on cushion biomechanics. TGFb3-treated cushions compacted less than controls, with compaction quantified as the ratio of cross-sectional area before and after treatment (A/A 0 = 0.92560.028 vs. A/ A 0 = 0.50860.017, p,0.0001 Figure 2A). This was unexpected because the Cytochalasin D (CytD, 1 mM) results suggested that compaction and stiffness are directly related. CytD inhibited cytoskeletal actin polymerization which resulted in a 5.3 fold decrease in strain energy density of the AV cushions relative to control (W CytD = 0.07260.016, Figure S2A). Without actin polymerization the AV cushion cells did not compact the matrix, and the cushion did not remodel into the spherical configuration observed in all other treatments. Instead, the post-treatment cushion area was significantly larger than initial area, suggesting a relaxation of pre-treatment actin forces (A/A 0 = 1.6060.03, Figure S2B). The TGFb3 results of stiffness increase with compaction decrease did not align with this trend. Alk5 inhibition did return compaction behavior to control levels (A/ A 0 = 0.57060.035 Figure S3), indicating that the stiffness and compaction results are both dependent on activation of canonical TGFb3 signaling. To better understand the relationship between stiffness and compaction, cushion opening angles were quantified to approximate differences in cell traction forces. The opening angle of TGFb3 cushions was 1.29 fold larger than controls (74.6u62.0u vs. 57.7u61.4u, p,0.001 Figure 2B), indicating that TGFb3 treated cushions did indeed have higher cell traction forces. Together, these results demonstrate that TGFb3 induces cushion stiffening through Alk5, but with a concurrent reduction in tissue compaction that suggests other processes are also affected.

TGFb3 increases AV cushion proliferation and mesenchymal phenotype
Contractile phenotype markers aSMA and RhoA were significantly upregulated with TGFb3 treatment, 5.360.4 and 2.160.3 fold (6 SEM) respectively ( Figure 3A), suggesting that TGFb3 induced residual tension is partially due to an increased migratory/contractile phenotype of resident cushion mesenchyme. TGFb3 treatment also upregulated mRNA expression of col1a2 mRNA (3.860.9, p,0.05) and cyclin b2 (3.960.7 fold, p,0.05), indicative of increased collagen I synthesis and cell proliferation, respectively. BrdU incorporation confirmed that TGFb3 increased cushion cell proliferation 2.2660.36 fold over controls (p,0.0001, Figure 3B). Collectively, these results strongly suggest that while TGFb3 treated AV cushion mesenchyme are more migratory/ contractile, concomitant increases in cell proliferation and matrix synthesis work to counteract aggregate matrix compaction. This explains how the TGFb3 treated cushions are biomechanically stiffer, but appear minimally compacted. Furthermore, TGFb3 treatment increased TGFb3 transcription (2.260.6 fold, p,0.05), indicating a potential positive feedback loop for TGFb3 control of AV cushion biomechanical remodeling.

5-HT potentiates TGFb3 signaling through 5-HT2b receptor
The effect of 5-HT dose on biomechanical remodeling, independently and in combination with TGFb3, was systematically evaluated through the stiffness and compaction metrics of the AV cushion organ culture system. 5-HT administration by itself had no statistically significant effect on cushion stiffness. Combined treatment of TGFb3 with physiological 5-HT (470 nM) increased AV cushion stiffness (W T+5-HT = 1.13660.035), but high 5-HT dose (5-HT+ = 47 mM) eliminated any TGFb3 induced stiffening effect (W T+5-HT+ = 0.45760.025, Figure 4). Neither selective inhibition of the 5-HT2a (MDL100907 10 nM), 5-HT2b (SB204741 2.6 mM) receptors, nor the serotonin transporter SERT (Flouxetine 10 mM) alone affected cushion stiffness (Figure 4). Yet in combination with TGFb3, the anti-5-HT2b treatment completely blocked TGFb3 dependent stiffness and compaction behavior (Figure 4 & Figure S3). Inhibition of the 5-HT2a receptor or SERT had no measurable effect on TGFb3 induced cushion biomechanics. The compaction and stiffness changes induced by 5-HT potentiated TGFb3 followed the same trend of TGFb3 treatment alone, with compaction decreasing as stiffness increased and vice versa ( Figure 4 & Figure S3). The additional stiffening effect of 5-HT with TGFb3 was also eliminated with Alk5 inhibition, as shown through the combined treatment of TGFb3+5-HT+anti-Alk5 in Figure S4. This combined treatment generated a strain energy density similar to the TGFb3+anti-Alk5 treatment (0.20960.023 Pa vs 0.24560.16 Pa, respectively), and further supported that the effects of 5-HT signaling on AV valve remodeling is dependent on canonical TGFb signaling. Together, these findings suggest that exogenous 5-HT acts through the 5-HT2b receptor to augment or impair TGFb3 induced cushion stiffening and compaction in a dosedependent manner.  Figure 5B). Proliferation-related gene cyclin b2 was not significantly affected by 5-HT in combination with TGFb3. These results suggest that exogenous 5-HT potentiates TGFb3 more likely through interaction with upstream activation points and/or TGFb3 synthesis, rather than by interacting with TGFb3 downstream targets directly.

5-HT modulates TGFb3 regulation of AV cushion mesenchyme phenotype
We also analyzed the mRNA expression of intracellular 5-HT (i5-HT) related genes transglutaminase 2 (TGM2) and SERT. i5-HT transamidates small GTPases and matrix proteins, in a process called ''serotonylation'' [44]. TGM2 is an i5-HT binding partner which assists transamidation of RhoA [55] and fibronectin [56], altering tissue mechanics through GTPase activation and matrix protein cross-linking, respectively. SERT mRNA expres-sion was significantly increased with 5-HT treatment (1.560.2 fold, p,0.05), but was downregulated with the 5-HT+ dose (0.4660.12 fold, Figure S5A). TGFb3 treatment stimulated a 4.061.0 fold increase in TGM2, but SERT transcription remained near control levels (0.7060.11, Figure S5B). Addition of 5-HT with TGFb3 significantly decreased SERT and TGM2 mRNA, regardless of 5-HT dose. Although TGFb3 treatment did upregulate TGM2, the downregulation of SERT by 5-HT treatment and the lack of mechanical changes seen with the SERT inhibitor suggest that serotonylation is not a primary mechanism of stiffness increase in the ex vivo culture remodeling results.

Elevated 5-HT induces atrioventricular valvuloseptal defects in ovo
As the effects of TGFb signaling on valve formation are well studied [17,18,57], we here test whether exogenous 5-HT administration in ovo alters valve morphogenesis. 5-HT administration in ovo at HH17 induced a spectrum of cardiac defects by HH36 (Day 10) as summarized in Table 1. Temporal and dosage dependant viability curves ( Figure S6A) showed that a 0.7 mg dose  was over 50% lethal at HH36, but administration of the same dose of 5-HT at HH25 or HH31 did not result in further lethality or defect formation (data not shown). The only gross malformations observed were localized to the heart and chest wall. Approximately 42% (24/57) of affected embryos exhibited an ectopic heart which protruded through an incomplete chest wall closure ( Figure  S6B). To confirm that interior defects resulted specifically from 5-HT exposure and not secondarily from the ectopia, an experimental thoracotomy was performed to model the ectopic condition. We found no statistically significant occurrence of any cardiac defects with experimental ectopia, supporting that 5-HT was responsible for the cardiac defects observed. A ventricular septal defect (VSD or SVSD) occurred in 42% (24/57) of the defective embryos. Approximately 18% (10/57) of the embryos exhibited double outlet right ventricle (DORV) defects. 5-HT administration also resulted in significantly enlarged atria with thinned walls in 35% (20/57) of the defective embryos (Table 1, Figure 6A). All of the embryos with DORV also exhibited highly stenotic or atretic atrioventricular (AV) valves ( Figure 6B), with the normally muscular flap valve in the right AV canal appearing thin and fibrous like the left AV valve. Regardless of gross cardiac defect identified, the average (0.14460.009 mm, mean 6 SEM) and minimal (0.08060.007 mm) thickness of the left AV septal leaflet was thinner in 5-HT treated embryos than controls (0.19160.009 and 0.16560.023 mm respectively, Figure 6C). No differences were found in mural leaflet thickness, or in the length of either leaflet. The reduction in AV valve thickness with 5-HT treatment indicated an increase in tissue compaction, and may possibly be a recapitulation of the migratory/contractile phenotype observed ex vivo.  Figure 7B). Interestingly, the TGFb3 mRNA expression was comparable to that observed in the ex ovo organ culture treatment of TGFb3 alone (2.260.6), 5-HT alone (1.960.1), and TGFb3+5-HT (2.060.3). aSMA and col1a2 mRNA were also upregulated in ovo with 5-HT, but less than with direct TGFb3 administration in ex vivo culture (aSMA -1.5 vs 5.7, RhoA -1.7 vs 2.1). The similar mRNA profiles of the candidate genes in both models suggested that 5-HT also potentiates TGFb signaling in AV cushions in ovo. To confirm that the 5-HT treatment was indeed modulating TGFb signaling activity in ovo, we quantified nuclear pSmad2/3 expression in HH25 cushions with and without 5-HT treatment (Figure 8). 5-HT treatment increased the number of cell nuclei with localized pSmad2/3 expression 2.660.8 fold over control embryos (0.2860.04 vs. 0.1160.03, p,0.01). Together these results demonstrate that 5-HT potentiates TGFb signaling in AV cushions to control contractile differentiation, proliferation, and biomechanical remodeling.

Discussion
In this study we implemented a quantitative organ culture assay that simultaneously interrogated the contributions of cellular and  molecular signaling to drive cushion tissue-level remodeling and biomechanical strengthening. TGFb3 stimulated a 2.5 fold increase in biomechanical stiffness (Figure 1), generated in part by an increase in cell traction. This contractile phenotype is a common outcome of TGFb signaling in post-natal valve tissue. For instance, porcine aortic valves express contractile marker aSMA when stimulated by TGFb1 in situ [58]. Porcine aortic valve interstitial cells (VICs) embedded in collagen gels expressed aSMA in response to TGFb1, and demonstrated significant gel compaction over untreated gels [59]. Similarly, TGFb3 treated embryonic AV progenitors compacted collagen gels to 10% of initial area [23]. Yet in contrast to these reports, TGFb3 induced contractility did not result in hyper-compacted AV cushions (Figure 2), but instead compacted less than controls. A key distinction between these two assays is that in vitro collagen gel cultures have much lower cell densities than our ex vivo system. The effect of proliferation on volume change is virtually undetectable in these gels, and cell traction dominates the compaction behavior. In native tissues, especially in the embryo, changes in cell proliferation and/or apoptosis have a significant impact on resulting tissue volume and apparent compaction. Hanging-drop culture of AV cushions enables precise control of the biochemical environment while maintaining the natural structural and cellular composition of the cushion. The lack of compaction with TGFb3 treatment is therefore most likely due to a counterbalancing from increases in cell proliferation ( Figure 3B) and ECM synthesis. This supports a mechanism of simultaneous tissue growth, matrix reorganization, and biomechanical stiffening during embryonic valve formation that is driven by a complex coordination of cell tractions, matrix synthesis, and cell proliferation. These findings underscore that embryonic valve mechanics, which is critical for proper valve function, cannot be inferred strictly from isolated compaction, proliferation, or matrix synthesis data, but is best measured directly.
The interplay of TGFb3 and 5-HT signaling was most notably seen through the potentiation of TGFb3 gene expression by 5-HT  Thoracotomy control of ectopic heart condition was created by mechanically debriding the chest dermis and pericardium at HH17. dose ( Figure 5). The physiological 5-HT concentration upregulated TGFb3 expression, while the high concentration downregulated expression. Upregulation of TGFb expression by 5-HT has been observed in several cardiac cells and tissues, though the molecular mechanism is still unclear. Adult aortic valve interstitial cells treated with 5-HT have increased TGFb1 activity, predominantly through the 5-HT2a receptor [25,60]. Neo-natal rat cardiac fibroblasts treated with 5-HT and 5-HT2a agonists upregulated aSMA protein expression, which is a marker for fibroblast differentiation and a gene induced by TGFb signaling [61]. Similarly, TGFb1 and aSMA expression were elevated in SERT cre-lox KO mice hearts through heightened 5-HT2a signaling in late embryonic stage mice, purportedly due to excess 5-HT from SERT inhibition [43]. Other reports point to 5-HT2b as the key mechanism. 5-HT administration in adult rats increased 5-HT2b mRNA expression in both aortic and mitral valves, demonstrating a positive response to 5-HT treatment [50]. SERT mRNA was downregulated in these valves denoting a negative response to elevated 5-HT, which our results also demonstrate ( Figure S5B). The 5-HT2b receptor, TGFb receptor type I and II, and the TGFb latent binding protein were all more expressed in canine myxomatous mitral valves than normal valves, suggesting a coupling of these two pathways through 5-HT2b [26]. Long-term 5-HT treatment of rats generated valve-related echocardiographic and histology defects [62], but these defects did not occur in rats simultaneously treated with a 5-HT2b inhibitor [51]. This suggests that the 5-HT2b receptor may be a key pathway for cardiac and valve tissue remodeling. Cardiac fibroblast studies indicate that 5-HT upregulates TGFb1 through a mutual transactivation of the epidermal growth factor (EGF) pathway and the 5-HT2b receptor [45,63]. Our results support a 5-HT2b dependant mechanism, as seen by 5-HT2b inhibition effectively blocking TGFb3 stiffening. The TGFb stiffening effect was independent of 5-HT2a and SERT. Although TGFb3 upregulated TGM2 expression, 5-HT treatment mitigated this expression which suggests TGM2 activity does not contribute to the enhanced stiffening of TGFb/5-HT signaling. High 5-HT also mitigated TGFb3 stiffening, which may be due to desensitization of the 5-HT2b receptor by sustained high 5-HT exposure. 5-HT increased pSmad2/3 phosphorylation in cushion mesenchyme in ovo. This suggests that 5-HT signaling through 5HT2b may interact with Smad2/3 signaling, but further studies are warranted to clarify potential roles of other intermediate or downstream targets.
In our in ovo model, systemic 5-HT elevation induced severe heart defects, including failure of the ventricular septum to close, ballooned atria, DORV, and hyper-contracted AV valves. Variations of these defects have been observed in other TGFb and 5-HT related studies. VSDs are the most prevalent congenital heart defects observed, occurring in approximately 50% of all clinical cardiac malformations [64,65]. Selective serotonin uptake inhibitors (SSRI) taken during the first trimester of pregnancy were associated with a statistical increase in VSD prevalence in newborns [66]. Our data supports elevated extracellular 5-HT as a possible cause of this correlation. Removal of TGFb secondary messenger Smad4 causes VSDs and other lethal congenital defects, which are presumed to be the consequence of decreased TGFb signaling [67]. Yet removal of TGFb inhibitory messenger Smad7 also generates VSDs [68], indicating that exacerbated TGFb signaling can also generate significant cardiac defects. The dilated atria observed in our model are not explicitly reported in other 5-HT studies, suggesting the defect may result from secondary effects, such as altered hemodynamics from valve incompetence. For instance, enlarged atria have been induced in zebrafish embryos through mechanical obstruction of the AV  canal [69]. Our avian model exhibited a small (18%), but statistically significant, penetrance of DORV, which is a predominant congenital defect in TGFb2 KO mice (87% penetrance) [17]. Collectively these defects highlight the morphogenetic potential of 5-HT in early cardiac development, and the similar spectrum of defects generated across 5-HT and TGFb related animal models.
An interaction of TGFb and 5-HT signaling was observed in ovo through the upregulation of TGFb3 and contractile genes in the AV cushions ( Figure 7B), the increase in pSmad2/3 expression (Figure 8), and the resulting thinned valve morphology ( Figure 6C). While the pSmad2/3 and mRNA expression confirms that aspects of the ex vivo results occur in ovo, it is unclear whether elevated TGFb signaling at HH25 is solely responsible for the thinned valve morphology observed at HH36. Hyperplastic and thickened AV valves occur in TGFb2 KO (31% penetrance) [17,70], and TGFb latent binding protein KO (81% penetrance) [57] animals, which supports this hypothesis. However, systemic 5-HT administration in adult rats generates thickened valves, with treatment duration dependent remodeling. Subcutaneous 5-HT injections for 7 days in adult rats produced thickened AV valves rich in GAGs [50], while 3-month treatment increased valve thickness, but consisted primarily of collagen [62]. Thickened, collagen-rich valves are also reported in adult SERT KO mice [71], and at late embryonic stage SERT KO pups [43]. Together these results indicate that elevated 5-HT signaling can instigate valvular remodeling in vivo, but changes in valve microstructure and morphology are clearly dependent on other factors such as treatment duration, specimen age, or secondary effects from accompanying congenital malformations. Altered hemodynamic loading can also generate defects, as evidenced through the serious malformations stimulated by mechanical perturbation [72,73]. Yet hemodynamic loading is simultaneously a consequence and stimulant of molecular signaling, interacting in a cyclical rather than a linear cause-effect manner. This again emphasizes the importance of direct assessment of mechanical stiffness, because it can distinguish the influence of these microstructure and microenvironment variations on valve performance.
Embryonic valve formation and maturation utilizes multiple TGFb isoforms in spatially and temporally restricted ways that are also somewhat different between species [16,18]. We chose to focus on TGFb3 over either TGFb1 or TGFb2 because of its principal role in cell invasion during chick cushion EMT [16], and confirmed increase in expression during post-EMT [23]. Our results establish a molecular mechanism for short-term (24 hours) TGFb3 stimulation on AV cushion biomechanical remodeling, but the effects of prolonged signaling on biomechanical and morphological changes remain unclear. This could be addressed with a combined in vivo/in vitro experimentation over more time points using a system like the approach presented here. The ex vivo culture system contains both endocardial and mesenchymal cells, but the lack of chick reactive antibodies prohibited the determination of cell specific responses. Our in ovo exogenous 5-HT administration model data complements existing data on genetic mutant animal models of TGFb and 5-HT related signaling in cardiac development [17,74]. Future studies will need to investigate whether the serotonin effects of TGFb3 change with TGFb3 dose.
In conclusion, tissue mechanics, cell phenotype, and molecular signaling all simultaneously direct and control tissue morphogenesis. Our results suggest that TGFb is a potent stimulator of cushion stiffening, and that 5-HT is a key regulator of this stimulating effect. Connecting signaling networks with cell and tissue level responses will become increasingly important for understanding post-EMT valve remodeling and potentially other embryonic remodeling events. The quantitative experimental systems presented herein are an attractive approach for elucidating these multi-scale mechanisms and their downstream consequences. Table S1 RT-PCR Primer Sequences. (TIF) Figure S1 Minimal ECM organization in HH25 cushion supports use of an isotropic mechanical testing technique. A) Confocal image of a HH25 cushion with ECM labeled via 5-DTAF protein stain at 106 magnification. B) 406 magnification. Note the lack of matrix fiber density or preferential fiber orientation at this stage of development. (TIF) Figure S2 Compaction-related stiffness control. A) Molecular inhibition of actin polymerization (Cytochalasin D, 1 mM) caused an 80-85% reduction in effective modulus. mean 6 SEM, n$6 *p,0.0001, t-test B) Cushion area increased with actin inhibition, resulting in a 3 fold decrease in measured compaction compared to control. Insets: Representative images of AV cushions before and after treatment, scale bar = 100 mm. mean 6 SEM, n$12, *p,0.0001, t-test. (TIF) Figure S3 TGFb3-induced decrease in compaction was blocked through inhibition of Alk5 (SB431542, 2.6 mM) or 5-HTR2b (SB204741 35 mM, anti-5-HT2b). Neither 5-HTR2a inhibitor (MDL100907 10 nM, anti-5-HT2a) nor serotonin transporter inhibitor (Fluoxetine 10 mM, anti-SERT) affected TGFb3 compaction behavior. mean 6 SEM, n$7, *p,0.05, ttest with respect to untreated controls. (TIF) Figure S4 TGFb3 and 5-HT stiffness generation is dependent on Alk5 signaling pathway. Strain energy density (Pa) of cushions treated with TGFb3 (1 ng/ml) only, TGFb3+Alk5 inhibitor (SB431542, 2.6 mM anti-Alk5), TGFb3+5-HT (470 nM), and TGFb3+5-HT+anti-Alk5. mean 6 SEM, n$8. Different letter pairings denotes statistically significant p,0.05, 2-way ANOVA. (TIF) Figure S5 Intracellular 5-HT uptake is modulated by 5-HT dose. A) 5-HT transporter (SERT) gene expression was downregulated via high 5-HT (47 mM, 5-HT+) dose, while transglutaminase 2 (TGM2) was not affected. The physiological dose of 5-HT (470 nM, 5-HT) had no effect on either SERT or TGM2 gene expression. B) TGFb3 (1 ng/ml) stimulated 4-fold increase in TGM2, which was mitigated by either doses of 5-HT. TGFb3 had no effect on SERT expression. mean 6 SEM, n = 3-4, *p,0.05, t-test. (TIF) Figure S6 Characterization of in ovo 5-HT administration model. A) Plot of avian embryo viability as a function of time and 5-HT dose. 5-HT administration to the surface of HH17 chick embryos resulted in greater than 70% lethality at dosages above 0.75 mg. The majority of deaths occurred within 48 hours of incubation. Doses of 0.5 mg and below were over 80% viable with virtually no morphological defects. Doses administered at later incubation times (Day 5, Day 7) did not result in lethality or defects by HH36 (data not shown). 5-HT administration at the predicted 50% lethality dose (0.7 mg/100 mL) resulted in 55% lethality by Day 10. B) Representative image of ectopic heart (arrow) and unclosed chest (dashed line) observed with both 5-HT treatment and thoracotomy sham controls. (TIF)