Female Anterior Cruciate Ligaments Exhibit a Muted Mechanobiological Response to Mechanical Loading

Female Anterior Cruciate Ligaments Exhibit a Muted Mechanobiological Response to Mechanical Loading | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Female Anterior Cruciate Ligaments Exhibit a Muted Mechanobiological Response to Mechanical Loading View ORCID Profile Lauren Paschall , View ORCID Profile Maxwell Konnaris , View ORCID Profile Erdem Tabdanov , View ORCID Profile Aman Dhawan , View ORCID Profile Spencer Szczesny doi: https://doi.org/10.1101/2025.05.16.654473 Lauren Paschall 1 The Pennsylvania State University, Department of Biomedical Engineering, State College , PA 2 The National Institutes of Health, National Cancer Institute, Center for Cancer Research, Pediatric Oncology Branch , Bethesda, MD Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Lauren Paschall Maxwell Konnaris 3 The Pennsylvania State University, Department of Statistics, State College , PA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Maxwell Konnaris Erdem Tabdanov 4 The Pennsylvania State University, Department of Orthopaedics and Rehabilitation , Hershey, PA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Erdem Tabdanov Aman Dhawan 5 The Pennsylvania State University, Department of Pharmacology , Hershey, PA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Aman Dhawan Spencer Szczesny 1 The Pennsylvania State University, Department of Biomedical Engineering, State College , PA 5 The Pennsylvania State University, Department of Pharmacology , Hershey, PA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Spencer Szczesny For correspondence: ses297{at}psu.edu Abstract Full Text Info/History Metrics Preview PDF Abstract Female athletes are significantly more likely to tear their anterior cruciate ligament (ACL) compared to their male counterparts. While there are several potential reasons for this, previous data from our lab demonstrated that female ACL explants have an impaired remodeling response to loading, which may prevent the repair of fatigue damage and lead to increased ACL rupture. The objective of this study was to identify the mechanisms driving the impaired remodeling of female ACLs to cyclic loading, including the role of estrogen. ACLs were harvested from male and female New Zealand white rabbits and cyclically loaded in a tensile bioreactor followed by bulk RNA-sequencing. Additional ACL explants treated with or without estradiol were analyzed using RT-qPCR to determine the regulatory effect of estrogen on markers for tissue remodeling and inflammatory cytokines with cyclic loading. We found that female ACLs exhibited significantly fewer differentially expressed genes (DEGs) in response to loading compared to male ACLs. Additionally, multiple mechanotransduction pathways were enriched with loading only in the male ACLs. While a few estrogen-related pathways were enriched in both male and female ACLs with loading, the expression of tissue remodeling markers was not different between estrogen treatment and vehicle control. Together, our findings highlight specific mechanotransduction pathways that may be responsible for the muted biological response of female ACLs to load, which provides a potential explanation for the increased rate of ACL tears in women. Introduction While anterior cruciate ligament (ACL) tears are common in general during athletic activity 1 , female athletes are 2–8 times more likely to tear their ACLs compared to males 2 – 4 . Potential reasons for this include differences in anatomy, sex hormones, body size, and neuromuscular control 5 – 7 . However, women may also be more prone to ACL injury because of sex-specific ACL remodeling in response to fatigue loading. Recent studies suggest that ACL injuries occur from fatigue failure and accumulation of tissue damage 8 . Previous data from our lab demonstrated that female rabbit ACL explants have an impaired remodeling response to load. Specifically, we found that female ACLs fail to increase anabolic gene expression with loading and have increased catabolic gene expression compared to male ACL explants 8 , which may prevent the repair of fatigue damage explaining the increased rate of ACL tears in women. However, the cause of this sex difference in ACL mechanobiology is unknown. Differences in sex hormones (estrogen) is one obvious potential reason for the sex-dependent ACL response to loading. Estrogen receptors (ERs) are present in both male and female ACL fibroblasts 9 . Additionally, estrogen influences multiple genes responsible for soft tissue remodeling 10 – 17 ; however, the magnitude and directionality of this effect on the ACL are inconsistent 10 , 12 , 13 , 18 – 25 . Despite these discrepancies, clinical data demonstrate that estradiol concentrations correlate with increased knee laxity 26 , which alone is sufficient to explain an 8-fold increase in risk of ACL injury 27 . This effect can be attenuated through oral contraceptives that suppress spikes in estrogen levels 28 , 29 . Together, these data suggest that estrogen has a direct effect on ACL tissue mechanics and may contribute to ACL rupture. Previous data also suggest that estrogen regulates the response of cells to mechanical loading. A prior study on ACL fibroblasts found that mechanical loading completely reverses the gene expression effect of estrogen 12 . However, the mechanisms underlying this interplay between estrogen and mechanical loading are unknown. In other tissues, ER signaling has been shown to interact with multiple mechanotransduction pathways (i.e., ERK, YAP/TAZ) 30 – 38 . However, the directionality of this effect is inconsistent. Specifically, studies on cancer cells show that ER signaling both increases 36 – 38 and decreases 33 – 35 YAP activation. While previous studies demonstrate that ER signaling regulates cellular mechanotransduction, it’s not possible to predict the effect of estrogen on ACL fibroblast mechanotransduction from prior data alone. Therefore, the objective of this study was to identify potential mechanisms driving the impaired remodeling response of female ACLs. We conducted bulk RNA-sequencing on cyclically loaded male and female rabbit ACL explants to determine the innate sex differences in ACL mechanobiology in the native tissue environment. Additionally, we investigated the effect of estrogen treatment on the expression of remodeling genes in male and female loaded ACL explants. Our overall hypothesis was that female ACLs would exhibit an impaired mechanobiological response to load (increased catabolism, decreased anabolism) and this response will be influenced by estrogen. Specifically, we hypothesize that female ACLs will have altered mechanotransduction (suppression of ERK and YAP) as well as increased ER activity in response to load compared to males. We also hypothesized that estrogen treatment will increase catabolic and inflammatory gene expression in both male and female ACL explants in response to load. Materials and Methods ACL Harvest A total of 33 male and female white New Zealand rabbits (2.8 – 3.2 kg, 14 – 16 weeks old) were euthanized ( Supplementary Table 1 ) and ACLs were isolated under sterile conditions from a study approved by the Institutional Animal Care and Use Committee at Penn State University. ACLs were harvested following a previously published protocol 8 . Briefly, all surrounding tissues were removed, and the femur and tibia were cut into bone blocks to grip during mechanical loading. Cross-sectional area was determined by assuming an elliptical cross-section and major and minor diameters were determined using calipers. Mechanical Loading ACLs were placed in a custom tensile bioreactor 39 with culture media (low-glucose DMEM, 2% penicillin-streptomycin, 5% FBS, 25 mm HEPES, 4 mm GlutaMAX, 1 mm 2-Phospho-L-ascorbic acid trisodium) and kept at 37°C and 5% CO 2. A subset of the ACLs was untreated to determine the innate mechanobiology of male and female ACL explants. To investigate the effect of estrogen, ACLs were treated with 10 nM 17β-estradiol (E2) (Sigma) or ethanol (0.1% v/v) as a vehicle control. The E2 concentration was based on a previous study investigating the role of estrogen and cyclic loading on isolated ACL fibroblasts 12 . Since underloading is detrimental to tendon/ligament homeostasis 40 , 41 a 0.1 MPa static load was applied for 18 and 20 hours to acclimate the untreated and treated (vehicle/estrogen) ACLs to culture conditions respectively. The bioreactor then cyclically loaded the samples to 4 MPa at 0.5 hz for 8 h. This stress level was chosen since it represents the stress that exhibited the largest gene expression differences between male and female ACLs 8 . Control samples were maintained under a 0.1 MPa static load for the same duration. RNA Extraction Immediately following loading, samples were removed and RNA was extracted following a previously published protocol 8 . Briefly, ACLs were excised from the bone blocks and were rinsed with ice cold RNase-free water and flash frozen in liquid nitrogen. The tissue was pulverized (6775 Freezer/Mill SPEX) and total RNA was extracted with RNeasy minicolumns (RNeasy Fibrous Tissue Kit, Qiagen). RNA concentration was quantified using a Qubit 4 Fluorometer (Thermo Fisher), and RNA integrity numbers (RINs) were determined using a TapeStation (Agilent). RNA-Sequencing Samples that passed quality control (RIN > 6, total RNA > 250 ng) were sent to GENEWIZ (South Plainfield, NJ) for bulk RNA-Sequencing ( Supplementary Table 1 ). cDNA libraries were prepared using the NEBNext Ultra II RNA Library Prep for Illumina using manufacturer’s instructions (NEB, Ipswich). The sequencing libraries were validated on the Agilent TapeStation (Agilent Technologies) and quantified using a Qubit 3.0 Fluorometer (Invitrogen) as well as qPCR (KAPA biosystems). The sequencing libraries were clustered on a lane of a NovaSeq 6000 S4 flow cell. The samples were sequenced with an average of 32 million reads per sample using a 2X150 bp paired end configuration. Image analysis and base calling were conducted by the Control software. Raw sequence data (.bcl files) generated by the sequencer were converted into fastq files and de-multiplexed using Illumina’s bcl2fastq 2.17 software. One mismatch was allowed for index sequence identification. Sequenced reads were trimmed to remove possible adapter sequences and nucleotides with poor quality using Trimmomatic v.0.36. The trimmed reads were mapped to the OryCun2.0 reference genome available on ENSEMBL using the STAR aligner v.2.5.2b. Unique gene hit counts were calculated by using FeatureCount from the Subread package v.1.5.2. Differential Gene Expression Analysis Analysis and visualization of the bulk RNA-sequencing was performed utilizing the NIH Integrated Analysis Platform (NIDAP) bulk RNA-sequencing pipeline that utilizes R programs developed by a team of NCI bioinformaticians (Plantair Technologies, Denver, Colorado). Raw hit counts were uploaded and quality control measures were performed using principal component analysis and supervised hierarchical clustering to determine intergroup relationships and identify outliers. A female loaded sample (FL3) was identified as an outlier and removed ( Supplementary Fig. 1 ). After outlier removal, normalization and DEG analysis was conducted utilizing the Limma Voom R package. Five comparisons were made to compare the response of each sex to load (female load, male load) and the relationships between sex in the fresh (sex effect fresh), static (sex effect static), and loaded (sex effect load) samples. The static control and male samples were used as the baseline for the comparisons. A moderated t-statistic was used to determine p-values, and a Benjamini-Hochberg false discovery rate was used to obtain the adjusted p-values. Genes with an adjusted p-value 1 were identified as differentially expressed genes (DEGs). A list of DEGs were extracted from the model for each of the five comparisons. Pathway Analysis and Upstream Regulator Analysis After differential gene expression analysis, all DEGs for each comparison were uploaded into Ingenuity Pathway Analysis (IPA, Qiagen, Version Q4 2024). Significantly enriched pathways were determined with a -log(p-value) > 1.3. Additionally, upstream regulator analysis was performed to determine potential regulators that could explain the observed changes in gene expression. To identify significantly activated or inhibited pathways or upstream regulators, a z-score was calculated in IPA for each pathway or upstream regulator. A positive z-score suggests inhibition while a positive z-score suggests activation. The magnitude of the z-score reflects the confidence of the predicted activation or inhibition state with |z-scores| > 2 considered strong prediction and |z-scores| > 1 considered moderate prediction. Reverse Transcription quantitative PCR (RT-qPCR) The RNA extracted from each sample was diluted to 1 ng/ul, and cDNA was synthesized (High-Capacity cDNA Reverse Transcription Kit with RNase Inhibitor, Thermo Fisher). Quantitative polymerase chain reaction (qPCR) was performed using TaqMan probes ( Supplementary Table 2 ) and a StepOne Plus Real-Time PCR system to measure the expression of anabolic (C OL1A1 , COL1A2 , LOX, COL3A1, TGF β 1, ACTA, TIMP1, TIMP3) , catabolic ( MMP1, MMP2, MMP10, MMP13 ), and inflammatory ( IL-1 β , PTGS2) genes. Additionally, ER target genes ( ESR1, GPER1, PGR ) were measured along with GAPDH as the reference gene. PCR efficiency and cycle number quantification (Cq) were obtained for each individual reaction using PCR-Miner (version 4.0) 42 . Samples that exhibited undetectable fluorescence for a given were excluded from the analysis for that gene. An outlier test (ROUT) was conducted to determine if any of the reaction efficiencies for each gene were an outlier. After removal of outliers, a single efficiency value was calculated for each probe by averaging all the individual sample efficiencies. Gene expression was quantified using the delta-delta Cq method correcting for primer efficiencies 43 . The fold change for each sample was calculated as: where the probe efficiencies are listed in Supplementary Table 2 . Statistical Analysis of RT-qPCR Data The fold change for each gene was log-normalized, and Mann-Whitney tests were conducted on the normalized data to compare the gene expression between treatments and sex. To determine the effect of mechanical loading on treatment, the estrogen treated and ethanol vehicle control samples were compared to their treatment-matched statically loaded control samples (for both males and females). To determine baseline treatment-specific gene expression differences, the estrogen treated static control samples were compared to the vehicle control static control samples (for both males and females). To determine baseline sex-specific gene expression differences, the female vehicle control samples were compared to the male vehicle control samples. Additionally, Mann-Whitney tests were used on the transformed data to determine the differential expression of each gene compared to the respective control condition. Given the numerous statistical comparisons, p-values were corrected using a false discovery rate analysis within each gene category (anabolic, catabolic, inflammatory, and ER target genes) with a Q-value of 5%. Statistical significance after corrections was set as p < 0.05 with statistical trends set at p < 0.10. All statistical analysis was performed using GraphPad Prism (version 10.3.1). Results Differential Gene Expression Analysis In response to cyclic loading, male ACLs had 259 DEGs (122 upregulated, 137 downregulated) while female ACLs only had 10 DEGs (4 upregulated, 6 downregulated). When comparing the response to load between female and male ACLs (sex effect load), there are 116 DEGs (25 upregulated, 91 downregulated in female ACLs compared to males) ( Fig. 1A , Supplementary Tables 3 – 5) . When looking at the top 20 DEGs in the male samples response to load, 75% of the DEGs were downregulated. Similarly, for the sex effect load comparison, 70% of the DEGs were further downregulated in female ACLs in response to load compared to males ( Table 1 ). These differences are despite the fact that there are minimal differences between female and male ACLs prior to loading. Specifically, there were only 40 DEGs (19 upregulated, 21 downregulated) in female static ACLs compared to male samples, and only 2 DEGs (1 upregulated, 1 downregulated) in female freshly harvested ACLs compared to males ( Fig. 1B , Supplementary Table 6 - 7 ). Download figure Open in new tab Figure 1: Differentially expressed genes. (A) Response of male and female ACLs in response to cyclic loading and comparing the response between sexes (sex effect load) (n = 3 males, n = 2 females). Male and female load are relative to their sex-matched static load. Sex effect load is comparing the female ACL response to load to the male ACL response to load. (B) Comparison of baseline gene expression differences of female ACLs relative to male ACLs (n = 3 static, n = 2 fresh). All DEGs had an absolute log2FC >1 and an adjusted p-value < 0.05. View this table: View inline View popup Download powerpoint Table 1: List of top 20 DEGs in response to load in male ACLs (and the sex effect of load (female ACLs response to load relative to male ACLs response to load). Genes indicative of tissue remodeling (anabolic, catabolic, inflammatory markers) 8 were searched in the list of DEGs for all five comparisons. In response to load, male ACLs downregulated catabolic gene MMP2 and inflammatory marker IL1 β while female ACLs exhibited no differentially expressed remodeling genes ( Table 2 ). When comparing across sex, there are no differentially expressed remodeling genes for any of the comparisons. ( Table 2 ). View this table: View inline View popup Download powerpoint Table 2: Differential expression of ECM remodeling genes in all comparisons. Bold text indicates statistical significance (p < 0.05). NAs indicate that the gene was not present in RNA-sequencing data. Differentially Expressed ER Target Genes Previously reported ER target genes 44 , 45 were searched in the list of DEGs for all five comparisons. In response to load, male ACLs upregulated ESR1 and multiple genes trended towards downregulation ( IGFBP4, TPD52L1, SERPINE1 ) while female ACLs exhibited no differentially expressed ER target genes ( Table 3 ). When comparing the DEGs across sex, there were surprisingly no differentially expressed ER target genes for the freshly harvested or static samples. In response to load, ESR1 trended towards downregulation in female loaded ACLs compared to male loaded ACLs ( Table 3 ). View this table: View inline View popup Download powerpoint Table 3: Differential expression of estrogen receptor target genes in all comparisons. Bold text indicates statistical significance (p < 0.05) and italicized text indicates trends (p < 0.10). Pathway Analysis of ACLs Response to Cyclic Load In response to cyclic loading, IPA found that male ACLs had 139 signaling pathways enriched, while female ACLs only had 31 enriched pathways ( Supplementary Table 8 – 9 ). When directly comparing across sex (sex effect load), there were 82 pathways enriched ( Supplementary Table 10 ). We then investigated which estrogen-related and mechanotransduction pathways were enriched with loading. For estrogen, extra-nuclear estrogen signaling was enriched in male ACLs in response to load and was trending towards inhibition (z-score: -1), while ESR-mediated signaling was enriched in female ACLs in response to load ( Supplementary Table 8 – 9 ). For mechanotransduction, male ACLs had 15 pathways enriched in response to load ( Table 4 ) with 8 pathways predicted to be downregulated and zero pathways predicted to be upregulated ( Table 4 ). The female ACLs exhibited no enriched pathways. Additionally, there were 9 pathways enriched when directly comparing across sex (sex effect load). Prominent pathways of interest included integrin signaling, FAK signaling, MAPK signaling, PI3K/AKT signaling, Rho GTPase cycle, and G-protein coupled receptor signaling. When comparing enriched pathways between male and female ACLs in response to load, 5 pathways were enriched in both male load and sex effect of load, which were generally predicted to be inhibited ( Table 5 ). Interestingly, FAK signaling was predicted to be inhibited in the sex effect load comparison, indicating an increased inhibition of FAK in female ACLs compared to male ACLs in response to load. View this table: View inline View popup Download powerpoint Table 4: Mechanotransduction pathways enriched in response to load. Pathways with a z-score > 0 are predicted to be activated and z-score < 0 are predicted to be inhibited. View this table: View inline View popup Download powerpoint Table 5: Mechanotransduction pathways enriched in both male ACL response to load and sex effect load. Pathways with a z-score > 0 are predicted to be activated and z-score < 0 are predicted to be inhibited. Finally, we also investigated pathways associated with extracellular matrix (ECM) organization. Male ACLs had 8 ECM pathways enriched in response to load ( Table 6 ), while female ACLs had no enriched pathways. Three pathways were enriched when directly comparing across sex (sex effect load). These included collagen synthesis and degradation, activation and inhibition of matrix metalloproteinases, and elastic fiber formation. Interestingly, most of these pathways were predicted to be inhibited in the male ACLs response to load. When looking at what pathways are different between male ACLs and female ACLs in response to load, 2 pathways related to collagen synthesis were enriched in both the male load and sex effect of load comparisons ( Table 7 ). View this table: View inline View popup Download powerpoint Table 6: Extracellular matrix pathways enriched in response to load. Pathways with a z-score > 0 are predicted to be activated and z-score < 0 are predicted to be inhibited. View this table: View inline View popup Download powerpoint Table 7: Extracellular matrix pathways enriched in both male ACLs response to load and sex effect load. Pathways with a z-score > 0 are predicted to be activated and z-score < 0 are predicted to be inhibited. Upstream Regulator Analysis of ACLs Response to Cyclic Load IPA upstream regulator analysis identified several key regulators of male ACL response to load ( Table 8 ), including multiple regulators related to sex hormones (beta-estradiol, ERBB2 , dihydrotestosterone, AR, ESR1 ), growth factors (TGFβ1, VEGF, IGF1) and inflammatory cytokines (TNF, LPS, IL1β, prostaglandin E2), with the majority of these regulators predicted to be inhibited. There were no significant upstream regualtors in the female ACLs response to load ( Table 8 ). In the sex effect load, AR and TGFβ1 were predicted to be inhibited while TNF was predicted to be activated. View this table: View inline View popup Download powerpoint Table 8: List of top 20 predicted upstream regulators of response to load in female and male ACLs (relative to sex-matched static control) and the sex effect of load (female ACL response to load relative to male ACL response to load). Regulators with a z-score > 0 are predicted to be activated and z-score < 0 are predicted to be inhibited. Effects of Estrogen Treatment on the Remodeling Response of ACLs to Cyclic Load When comparing the expression of remodeling genes between male and female statically loaded ACLs cultured without estrogen (vehicle control), there were multiple anabolic ( COL1A2, LOX, TGF β 1, TIMP3 ), catabolic ( MMP1, MMP2, MMP13 ), and inflammatory ( IL1 β , PTGS2 ) markers significantly upregulated in the female ACLs compared to males and MMP1 0 was trending toward upregulation ( Supplementary Fig. 2A ). However, there was no differential expression of any of the anabolic, catabolic, or inflammatory markers with estrogen treatment in the statically loaded male or female ACLs ( Fig. 2 A-B ). This was also true when comparing the response of estrogen treatment between male and female statically loaded ACLs ( Fig. 2C ). Download figure Open in new tab Figure 2: Basal remodeling gene expression response of estrogen treatment. RT-qPCR analysis of statically loaded estrogen treated ACLs relative to their sex-matched ethanol vehicle control (represented by the red box and whisker plots). (A) Statically loaded male ACL response. (B) Statically loaded female ACL response. (C) Comparison of statically loaded male and female ACL response to estrogen treatment (n = 6 for all male samples, n =7 for female vehicle, n = 5 – 6 for female estrogen treated). Data represented as box and whisker plot with the whiskers representing the min and max data. No differential gene expression was observed with cyclical loading in the male ACLs compared to static controls with vehicle or estrogen treatment ( Fig. 3 A-B ). Direct comparisons between the vehicle control and estrogen-treated male samples found that there was no effect of estrogen on gene expression ( Supplementary Fig. 3A ). Consistently, the female cyclically loaded ACLs exhibited no differential gene expression compared to static controls with either vehicle or estrogen treatment ( Fig. 3 C-D ). Direct comparisons between vehicle and estrogen-treated cyclically loaded female samples found no effect of estrogen on gene expression ( Supplementary Fig. 3B ). Finally, direct comparison between the male and female cyclically loaded ACLs found that there were no sex differences in the vehicle or estrogen-treated groups ( Supplementary Fig. 3 C-D ). Download figure Open in new tab Figure 3: Effect of estrogen treatment on ACL remodeling in response to load. RT-qPCR analysis of cyclically loaded male and female ACLs relative to their sex and treatment-matched 0.1 MPa static load (represented by the red box and whisker plots). (A) Gene expression changes in vehicle control male ACLs. (B) Gene expression changes in estrogen treated male ACLs. (C) Gene expression changes in vehicle control female ACLs. (D) Gene expression changes in female estrogen treated ACLs (n = 6 for all male samples, n =7 for female vehicle, n = 6 for female estrogen treated). Data represented as box and whisker plot with the whiskers representing the min and max data. Effects of Estrogen Treatment on Estrogen Target Genes of ACLs to Cyclic Load In contrast to our RNA-sequencing data (sex effect static), we found GPER1 to be significantly upregulated in the statically loaded female ACLs compared to males without estrogen treatment ( Supplementary Fig. 2B ). However, there was no differential expression of any of the ER genes in the statically loaded male or female ACLs in response to estrogen treatment ( Fig. 4 A-B ). Additionally, there were no differences in gene expression when comparing the response of estrogen treatment between male and female statically loaded ACLs ( Fig. 4C ). Download figure Open in new tab Figure 4: Basal estrogen receptor target genes response to estrogen treatment. RT-qPCR analysis of statically loaded estrogen treated ACLs relative to their sex-matched ethanol vehicle control (represented by the red box and whisker plots). (A) Statically loaded male ACL response. (B) Statically loaded female ACL response. (C) Comparison of statically loaded male and female ACLs response to estrogen treatment (n = 5 for male vehicle, n = 6 for male estrogen treated, n = 6 - 7 for female vehicle, n = 6 for female estrogen treated). Data represented as box and whisker plot with the whiskers representing the min and max data. Consistent with our findings with the markers of tissue remodeling, there was no differential expression of the estrogen target genes with cyclic loading in the male ( Fig. 5 A-B ) or female samples ( Fig. 5 C-D ) with or without estrogen treatment. Similarly, there was no effect of estrogen ( Supplementary Fig. 4A-B ) or sex ( Supplementary Fig. 4C-D ) on the expression of estrogen target genes in response to loading. Download figure Open in new tab Figure 5: Effect of estrogen treatment on estrogen receptor target genes in response to load. RT-qPCR analysis of cyclically loaded male and female ACLs relative to their sex and treatment-matched 0.1 MPa static load. (A) Gene expression changes in vehicle control male ACLs. (B) Gene expression changes in estrogen treated male ACLs. (C) Gene expression changes in vehicle control female ACLs. (D) Gene expression changes in female estrogen treated ACLs (n = 6 for all male samples, n =6 - 7 for female vehicle, n = 6 for female estrogen treated). Data represented as box and whisker plot with the whiskers representing the min and max data. Discussion The objective of this work was to investigate potential mechanisms driving the impaired remodeling response of female ACLs in response to load. Consistent with our hypothesis, our RNA-sequencing data demonstrated that female ACLs had a muted response to mechanical loading compared to male ACLs. Specifically, female ACLs had significantly fewer DEGs in response to load compared to male ACLs ( Fig. 1A ). This is despite the fact that there were minimal differences in gene expression between freshly harvested female and male ACLs ( Fig. 1B ). Suggesting, that the differences in the number of DEGs in response to load between male and female ACLs are due to differences in mechanobiology rather than basal biological differences. Furthermore, multiple mechanotransduction pathways were significantly enriched in male ACLs in response to load while none were enriched in female ACLs ( Table 4 ). Surprisingly, half of these mechanotransduction pathways were predicted to be downregulated in the male samples. However, a number of prominent mechanotransduction pathways, including MAPK (i.e., ERK) and PI3K/AKT, were enriched without likely being downregulated, suggesting that these pathways may be important for the mechanobiological response of male ACLs. Additionally, when comparing the response to load between male and female ACLs, FAK (focal adhesion kinase) signaling was predicted to be further inhibited with load in female ACLs compared to male ACLs ( Table 4 , 5 ). Therefore, FAK signaling may also help explain sex differences in ACL mechanobiology given that FAK plays a prominent role in cell-matrix mechanotransduction and regulates the mechanoresponse of tendon cells 46 . Together, our RNA-sequencing data support the hypothesis that female ACLs have a muted mechanobiological response to load potentially inhibiting their ability to repair tissue damage that predisposes the ACL to rupture. Our secondary hypothesis was that the impaired mechanobiological response in female ACLs is due to estrogen signaling. Surprisingly, there were no basal differences in ER target genes between male and female ACLs ( Table 3 ). This could be due to the fact that we did not measure estrogen concentrations to confirm that estrogen levels were higher in female rabbits at the time of euthanasia. However, since rabbits are induced ovulators, female rabbits exhibit relatively stable serum estrogen concentrations over time 47 . Additionally, previous studies suggest that serum estrogen concentration is 10 times higher in female rabbits (20 – 50 pg/ml) 19 , 20 compared to male rabbits (3 pg/ml) 48 . These data suggest that the estrogen concentration should be consistently higher in the female ACL samples, suggesting that the chosen estrogen targets are not representative of estrogen signaling in rabbits. Downstream targets of estrogen activation are tissue and species specific 49 , 50 , and there are no data available on ER target genes in the ACL or rabbit, providing difficulties in validating ER activation. Additionally, we did not observe differential expression of the chosen ER targets in the female ACLs in response to load. However, the ESR-mediated signaling pathway was enriched in the female ACL response to load ( Supplementary Table 9 ). Additionally, the extra-nuclear estrogen signaling was enriched and trending toward inhibition in the male ACL response to load ( Supplementary Table 8 ). Furthermore, upstream regulator analysis revealed inhibition of ESR1 in the male ACL response to load ( Table 8 ). Therefore, while rabbit-specific markers of estrogen signaling may need to be identified, our RNA-sequencing data suggest that estrogen-related signaling may help explain sex differences in ACL mechanobiology. To further investigate this, we cyclically loaded male and female ACLs in the presence of exogenous estradiol and measured gene expression of markers of ECM remodeling. Unfortunately, we weren’t to replicate our previous work that utilized the identical loading protocol and identified a sex dependent response in these genes with loading 8 . Given that there was no effect of cyclic loading on our chosen remodeling genes, we were unable to evaluate the effect of estrogen ( Fig. 3 , Supplementary Fig. 3 A-B ). It’s possible that the lack of effect with cyclic loading is due to a poor choice of gene markers for tissue remodeling. Specifically, while our RNA-sequencing data found 259 DEGs with loading in male ACLs, only two of those genes were chosen as markers of tissue remodeling for our PCR experiments ( MMP2 , IL1 β) ( Table 2 ). However, when looking at the pathway analysis related to extracellular organization, there were only three pathways enriched when comparing the response to load between female and male ACLs ( Table 6 ). Furthermore, only a few genes enriched in these pathways are relevant to remodeling of the bulk ECM (MMP2, MMP14, COL5A3), whereas the others are more relevant to the local pericellular matrix (ADAM10, ADAM12, COL4A1, COL4A2, LAMB1) ( Supplementary Table 11 ). Together, this suggests that the mechanobiological difference between male and female ACLs may not involve bulk ECM remodeling. Still, the RNA-sequencing data suggest that ACLs do have a sex-dependent mechanobiological response, and future work is needed to identify the role of estrogen in driving this differential response. It’s also possible that the absence of an effect of estrogen on ACL explant gene expression was due to insufficient estradiol treatment. However, we used the same concentration (10 -8 M) that increased expression of COL1 and COL3 in a previous study of isolated ACL fibroblasts 12 . Additionally, our chosen concentration is three orders of magnitudes greater than serum estrogen concentrations found in rabbits (3.5 x 10 -11 M) 19 , 20 . Another possibility is that there was sufficient estradiol in the culture media to saturate estrogen signaling. While we did not use charcoal-stripped FBS in our culture media, there is only approximately 2 pg/ml (< 10 -11 M) of estradiol in 10% FBS, which was not sufficient to saturate the response of isolated ACL fibroblasts to estradiol 12 . A final possibility is that the effect of estrogen is no longer present at the end of our culture period (28 h), given that the effect of estrogen treatment on collagen gene expression decreased after 24 h 12 . A final limitation to this study was that our RNA-sequencing sample size was small for the loaded female samples due to the removal of an outlier (n=2) ( Supplementary Fig. S1 ). While it’s possible that having a smaller sample size could reduce statistical power and lead to fewer DEGs, this is likely not the reason for the low number of DEGs in the female ACLs response to load. Specifically, it’s clear from the PCA plot of the samples included in the sequencing analysis (Supplementary Fig. S1B ) that the male loaded samples separated from the male static samples to a greater amount than the female samples, which is consistent with the fewer number of DEGs in the female ACLs in response to load. In conclusion, our findings suggest that female ACLs have a muted mechanobiological response to load compared to male ACLs. This may be due to decreased activation of mechanotransduction pathways, including FAK, MAPK (i.e., ERK), and PI3K/AKT, as well as increased estrogen signaling in female ACLs in response to load. These data support the novel hypothesis that females exhibit an impaired mechanobiological response to load making female athletes more susceptible to ACL tears. Future work could compare the structural and mechanical changes that occur between male and female ACLs in response to load. This will enable a direct comparison of male and female ACLs ability to repair tissue damage. Future work is needed to better understand the role of estrogen signaling in ACL mechanobiology. Together, these data can provide insight into the ability of female ACLs to repair tissue damage and possible approaches for reducing the risk of ACL ruptures in females by modulating estrogen levels (e.g., contraceptive use). Figure Captions Download figure Open in new tab Supplementary Figure 1: Principal component analysis. Visual distribution of samples based on two principal components. Each point represents an individual sample, labeled with the sample number. (A) Principal component analysis without the outlier removal (FL3). (B) Principal component analysis with the outlier removal. Download figure Open in new tab Supplementary Figure 2: Comparison of vehicle control female and male ACL gene expression at baseline. RT-qPCR analysis of statically loaded female samples relative to male statically loaded samples. (A) Remodeling genes. (B) Estrogen receptor target genes (n = 6 for male samples and n = 7 for female samples for remodeling genes, n = 6 for male samples and n = 6-7 for female samples for estrogen receptor target genes) *p < 0.05, #p < 0.1. Data represented as box and whiskers plot with the whiskers representing the min and max data. Download figure Open in new tab Supplementary Figure 3: Comparisons of the male and female ACLs remodeling response to load. RT-qPCR analysis of cyclically loaded male and female ACLs relative to their sex and treatment-matched 0.1 MPa static load. (A) Comparison of vehicle control and estrogen treatment in male samples. (B) Comparison of vehicle control and estrogen treatment in female samples. (C) Comparison of male and female ACLs response to vehicle treatment. (D) Comparison of male and female ACLs response to estrogen treatment (n = 6 for all male samples, n =7 for female vehicle, n = 6 for female estrogen treated). Data represented as box and whisker plot with the whiskers representing the min and max data. Download figure Open in new tab Supplementary Figure 4: Comparisons of the male and female ACLs estrogen receptor target genes in response to load. RT-qPCR analysis of cyclically loaded male and female ACLs relative to their sex and treatment-matched 0.1 MPa static load. (A) Comparison of vehicle control and estrogen treatment in male samples. (B) Comparison of vehicle control and estrogen treatment in female samples. (C) Comparison of male and female ACL response to vehicle treatment. (D) Comparison of male and female ACLs response to estrogen treatment (n = 6 for all male samples, n =6 - 7 for female vehicle, n = 6 for female estrogen treated). Data represented as box and whisker plot with the whiskers representing the min and max data. Acknowledgements We would like to acknowledge the Genomics Core at Penn State for running the PCR plates. Funding for this project was provided by the Orthopaedic Research and Education Foundation (234995) and the Congressionally Directed Medical Research Program (W81XWH2110152). Funder Information Declared Orthopaedic Research and Education Foundation , 234995 Congressionally Directed Medical Research Program , W81XWH2110152 References 1. ↵ Sanders TL , Maradit Kremers H , Bryan AJ , et al. 2016 . Incidence of Anterior Cruciate Ligament Tears and Reconstruction: A 21-Year Population-Based Study . Am. J. Sports Med . 44 ( 6 ): 1502 – 1507 . OpenUrl CrossRef PubMed 2. ↵ Arendt EA , Agel J , Dick R . 1999 . Anterior cruciate ligament injury patterns among collegiate men and women . J. Athl. Train . 34 ( 2 ): 86 – 92 . OpenUrl PubMed Web of Science 3. Prodromos CC , Han Y , Rogowski J , et al. 2007 . A meta-analysis of the incidence of anterior cruciate ligament tears as a function of gender, sport, and a knee injury-reduction regimen . Arthrosc. J. Arthrosc. Relat. Surg. Off. Publ. Arthrosc. Assoc. N. Am. Int. Arthrosc. Assoc . 23 ( 12 ): 1320 – 1325 .e6. OpenUrl CrossRef 4. ↵ Joseph AM , Collins CL , Henke NM , et al. 2013 . A multisport epidemiologic comparison of anterior cruciate ligament injuries in high school athletics . J. Athl. Train . 48 ( 6 ): 810 – 817 . OpenUrl CrossRef PubMed 5. ↵ Putera KH , Kim J , Baek SY , et al. 2023 . Fatigue-driven compliance increase and collagen unravelling in mechanically tested anterior cruciate ligament. Commun . Biol . 6 ( 1 ): 564 . OpenUrl 6. Grodman LH , Beaulieu ML , Ashton-Miller JA , Wojtys EM . 2023 . Levels of ACL-straining activities increased in the six months prior to non-contact ACL injury in a retrospective survey: evidence consistent with ACL fatigue failure . Front. Physiol . 14 : 1166980 . OpenUrl CrossRef PubMed 7. ↵ Beaulieu ML , DeClercq MG , Rietberg NT , et al. 2021 . The Anterior Cruciate Ligament Can Become Hypertrophied in Response to Mechanical Loading: A Magnetic Resonance Imaging Study in Elite Athletes . Am. J. Sports Med . 49 ( 9 ): 2371 – 2378 . OpenUrl CrossRef PubMed 8. ↵ Paschall L , Carrozzi S , Tabdanov E , et al. 2024 . Cyclic loading induces anabolic gene expression in ACLs in a load-dependent and sex-specific manner . J. Orthop. Res. Off. Publ. Orthop. Res. Soc . 42 ( 2 ): 267 – 276 . OpenUrl CrossRef 9. ↵ Liu SH , al-Shaikh R , Panossian V , et al. 1996 . Primary immunolocalization of estrogen and progesterone target cells in the human anterior cruciate ligament . J. Orthop. Res. Off. Publ. Orthop. Res. Soc . 14 ( 4 ): 526 – 533 . OpenUrl CrossRef 10. ↵ Liu SH , Al-Shaikh RA , Panossian V , et al. 1997 . Estrogen affects the cellular metabolism of the anterior cruciate ligament. A potential explanation for female athletic injury . Am. J. Sports Med . 25 ( 5 ): 704 – 709 . OpenUrl CrossRef PubMed Web of Science 11. Hama H , Yamamuro T , Takeda T . 1976 . Experimental studies on connective tissue of the capsular ligament. Influences of aging and sex hormones . Acta Orthop. Scand . 47 ( 5 ): 473 – 479 . OpenUrl CrossRef PubMed Web of Science 12. ↵ Lee C-Y , Liu X , Smith CL , et al. 2004 . The combined regulation of estrogen and cyclic tension on fibroblast biosynthesis derived from anterior cruciate ligament . Matrix Biol. J. Int. Soc. Matrix Biol . 23 ( 5 ): 323 – 329 . OpenUrl CrossRef 13. ↵ Lee CA , Lee-Barthel A , Marquino L , et al. 2015 . Estrogen inhibits lysyl oxidase and decreases mechanical function in engineered ligaments . J. Appl. Physiol. Bethesda Md 1985 118 ( 10 ): 1250 – 1257 . OpenUrl 14. Wahl LM , Blandau RJ , Page RC . 1977 . Effect of hormones on collagen metabolism and collagenase activity in the pubic symphysis ligament of the guinea pig . Endocrinology 100 ( 2 ): 571 – 579 . OpenUrl CrossRef PubMed Web of Science 15. Rajabi MR , Dodge GR , Solomon S , Poole AR . 1991 . Immunochemical and immunohistochemical evidence of estrogen-mediated collagenolysis as a mechanism of cervical dilatation in the guinea pig at parturition . Endocrinology 128 ( 1 ): 371 – 378 . OpenUrl CrossRef PubMed Web of Science 16. Salamonsen LA , Woolley DE . 1996 . Matrix metalloproteinases in normal menstruation . Hum. Reprod. Oxf. Engl . 11 Suppl 2 : 124 – 133 . OpenUrl CrossRef 17. ↵ Schneikert J , Peterziel H , Defossez PA , et al. 1996 . Androgen receptor-Ets protein interaction is a novel mechanism for steroid hormone-mediated down-modulation of matrix metalloproteinase expression . J. Biol. Chem . 271 ( 39 ): 23907 – 23913 . OpenUrl Abstract / FREE Full Text 18. ↵ Slauterbeck J , Clevenger C , Lundberg W , Burchfield DM . 1999 . Estrogen level alters the failure load of the rabbit anterior cruciate ligament . J. Orthop. Res. Off. Publ. Orthop. Res. Soc . 17 ( 3 ): 405 – 408 . OpenUrl CrossRef 19. ↵ Komatsuda T , Sugita T , Sano H , et al. 2006 . Does estrogen alter the mechanical properties of the anterior cruciate ligament? An experimental study in rabbits . Acta Orthop . 77 ( 6 ): 973 – 980 . OpenUrl CrossRef PubMed Web of Science 20. ↵ Hattori K , Sano H , Komatsuda T , et al. 2010 . Effect of estrogen on tissue elasticity of the ligament proper in rabbit anterior cruciate ligament: measurements using scanning acoustic microscopy . J. Orthop. Sci. Off. J. Jpn. Orthop. Assoc . 15 ( 4 ): 584 – 588 . OpenUrl 21. Strickland SM , Belknap TW , Turner SA , et al. 2003 . Lack of hormonal influences on mechanical properties of sheep knee ligaments . Am. J. Sports Med . 31 ( 2 ): 210 – 215 . OpenUrl CrossRef PubMed Web of Science 22. Wentorf FA , Sudoh K , Moses C , et al. 2006 . The effects of estrogen on material and mechanical properties of the intra- and extra-articular knee structures . Am. J. Sports Med . 34 ( 12 ): 1948 – 1952 . OpenUrl CrossRef PubMed Web of Science 23. Warden SJ , Saxon LK , Castillo AB , Turner CH . 2006 . Knee ligament mechanical properties are not influenced by estrogen or its receptors . Am. J. Physiol. Endocrinol. Metab . 290 ( 5 ): E1034 – 1040 . OpenUrl CrossRef PubMed Web of Science 24. Yu WD , Liu SH , Hatch JD , et al. 1999 . Effect of estrogen on cellular metabolism of the human anterior cruciate ligament . Clin. Orthop . ( 366 ): 229 – 238 . 25. ↵ Yu WD , Panossian V , Hatch JD , et al. 2001 . Combined effects of estrogen and progesterone on the anterior cruciate ligament . Clin. Orthop . ( 383 ): 268 – 281 . 26. ↵ Shultz SJ , Kirk SE , Johnson ML , et al. 2004 . Relationship between sex hormones and anterior knee laxity across the menstrual cycle . Med. Sci. Sports Exerc . 36 ( 7 ): 1165 – 1174 . OpenUrl CrossRef PubMed Web of Science 27. ↵ Myer GD , Ford KR , Paterno MV , et al. 2008 . The effects of generalized joint laxity on risk of anterior cruciate ligament injury in young female athletes . Am. J. Sports Med . 36 ( 6 ): 1073 – 1080 . OpenUrl CrossRef PubMed Web of Science 28. ↵ Martineau PA , Al-Jassir F , Lenczner E , Burman ML . 2004 . Effect of the oral contraceptive pill on ligamentous laxity . Clin. J. Sport Med. Off. J. Can. Acad. Sport Med . 14 ( 5 ): 281 – 286 . OpenUrl CrossRef 29. ↵ Lee H , Petrofsky JS , Yim J . 2015 . Do Oral Contraceptives Alter Knee Ligament Damage with Heavy Exercise? Tohoku J. Exp. Med . 237 ( 1 ): 51 – 56 . OpenUrl CrossRef PubMed 30. ↵ Dos Santos EG , Dieudonne MN , Pecquery R , et al. 2002 . Rapid nongenomic E2 effects on p42/p44 MAPK, activator protein-1, and cAMP response element binding protein in rat white adipocytes . Endocrinology 143 ( 3 ): 930 – 940 . OpenUrl CrossRef PubMed Web of Science 31. Levin ER . 2005 . Integration of the extranuclear and nuclear actions of estrogen . Mol. Endocrinol. Baltim. Md 19 ( 8 ): 1951 – 1959 . OpenUrl CrossRef 32. Marino M , Galluzzo P , Ascenzi P . 2006 . Estrogen signaling multiple pathways to impact gene transcription . Curr. Genomics 7 ( 8 ): 497 – 508 . OpenUrl CrossRef PubMed Web of Science 33. ↵ Lachowski D , Cortes E , Matellan C , et al. 2020 . G Protein-Coupled Estrogen Receptor Regulates Actin Cytoskeleton Dynamics to Impair Cell Polarization . Front. Cell Dev. Biol . 8 : 592628 . OpenUrl CrossRef PubMed 34. Rice A , Cortes E , Lachowski D , et al. 2020 . GPER Activation Inhibits Cancer Cell Mechanotransduction and Basement Membrane Invasion via RhoA . Cancers 12 ( 2 ): 289 . OpenUrl CrossRef PubMed 35. ↵ Jeon Y , Yoo JE , Rhee H , et al. 2021 . YAP inactivation in estrogen receptor alpha-positive hepatocellular carcinoma with less aggressive behavior . Exp. Mol. Med . 53 ( 6 ): 1055 – 1067 . OpenUrl CrossRef PubMed 36. ↵ Wang Z , Sun L , Liang S , et al. 2019 . GPER stabilizes F-actin cytoskeleton and activates TAZ via PLCβ-PKC and Rho/ROCK-LIMK-Cofilin pathway . Biochem. Biophys. Res. Commun . 516 ( 3 ): 976 – 982 . OpenUrl CrossRef PubMed 37. Muhammad JS , Guimei M , Jayakumar MN , et al. 2021 . Estrogen-induced hypomethylation and overexpression of YAP1 facilitate breast cancer cell growth and survival . Neoplasia N. Y. N 23 ( 1 ): 68 – 79 . OpenUrl CrossRef 38. ↵ Zhou X , Wang S , Wang Z , et al. 2015 . Estrogen regulates Hippo signaling via GPER in breast cancer . J. Clin. Invest . 125 ( 5 ): 2123 – 2135 . OpenUrl CrossRef PubMed 39. ↵ Pedaprolu K , Szczesny SE . 2022 . A Novel, Open-Source, Low-Cost Bioreactor for Load-Controlled Cyclic Loading of Tendon Explants . J. Biomech. Eng . 144 ( 8 ): 084505 . OpenUrl CrossRef PubMed 40. ↵ Arnoczky SP , Tian T , Lavagnino M , Gardner K . 2004 . Ex vivo static tensile loading inhibits MMP-1 expression in rat tail tendon cells through a cytoskeletally based mechanotransduction mechanism . J. Orthop. Res. Off. Publ. Orthop. Res. Soc . 22 ( 2 ): 328 – 333 . OpenUrl CrossRef 41. ↵ Wang T , Lin Z , Ni M , et al. 2015 . Cyclic mechanical stimulation rescues achilles tendon from degeneration in a bioreactor system . J. Orthop. Res. Off. Publ. Orthop. Res. Soc . 33 ( 12 ): 1888 – 1896 . OpenUrl CrossRef 42. ↵ Zhao S , Fernald RD . 2005 . Comprehensive algorithm for quantitative real-time polymerase chain reaction . J. Comput. Biol. J. Comput. Mol. Cell Biol . 12 ( 8 ): 1047 – 1064 . OpenUrl CrossRef 43. ↵ Pfaffl M. [date unknown]. Relative Quantification . In: Real-time PCR. International University Line . p 63 – 82 . 44. ↵ Lin C-Y , Ström A , Vega VB , et al. 2004 . Discovery of estrogen receptor alpha target genes and response elements in breast tumor cells . Genome Biol . 5 ( 9 ): R66 . OpenUrl CrossRef PubMed 45. ↵ Verhaegh W , van Ooijen H , Inda MA , et al. 2014 . Selection of personalized patient therapy through the use of knowledge-based computational models that identify tumor-driving signal transduction pathways . Cancer Res . 74 ( 11 ): 2936 – 2945 . OpenUrl Abstract / FREE Full Text 46. ↵ Leahy TP , Chenna SS , Soslowsky LJ , Dyment NA . 2024 . Focal adhesion kinase regulates tendon cell mechanoresponse and physiological tendon development . FASEB J. Off. Publ. Fed. Am. Soc. Exp. Biol . 38 ( 17 ): e70050 . OpenUrl 47. ↵ Chen X , Jin R , Yang A , et al. 2023 . Behavioral and Physiological Differences in Female Rabbits at Different Stages of the Estrous Cycle . Anim. Open Access J. MDPI 13 ( 21 ): 3414 . OpenUrl 48. ↵ Celayir S . 2003 . Effects of different sex hormones on male rabbit urodynamics: an experimental study . Horm. Res . 60 ( 5 ): 215 – 220 . OpenUrl CrossRef PubMed 49. ↵ Tee MK , Rogatsky I , Tzagarakis-Foster C , et al. 2004 . Estradiol and selective estrogen receptor modulators differentially regulate target genes with estrogen receptors alpha and beta . Mol. Biol. Cell 15 ( 3 ): 1262 – 1272 . OpenUrl Abstract / FREE Full Text 50. ↵ Leitman DC , Paruthiyil S , Yuan C , et al. 2012 . Tissue-specific regulation of genes by estrogen receptors . Semin. Reprod. Med . 30 ( 1 ): 14 – 22 . OpenUrl CrossRef PubMed View the discussion thread. Back to top Previous Next Posted May 21, 2025. Download PDF Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Female Anterior Cruciate Ligaments Exhibit a Muted Mechanobiological Response to Mechanical Loading Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Female Anterior Cruciate Ligaments Exhibit a Muted Mechanobiological Response to Mechanical Loading Lauren Paschall , Maxwell Konnaris , Erdem Tabdanov , Aman Dhawan , Spencer Szczesny bioRxiv 2025.05.16.654473; doi: https://doi.org/10.1101/2025.05.16.654473 Share This Article: Copy Citation Tools Female Anterior Cruciate Ligaments Exhibit a Muted Mechanobiological Response to Mechanical Loading Lauren Paschall , Maxwell Konnaris , Erdem Tabdanov , Aman Dhawan , Spencer Szczesny bioRxiv 2025.05.16.654473; doi: https://doi.org/10.1101/2025.05.16.654473 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Bioengineering Subject Areas All Articles Animal Behavior and Cognition (7621) Biochemistry (17645) Bioengineering (13867) Bioinformatics (41872) Biophysics (21416) Cancer Biology (18549) Cell Biology (25443) Clinical Trials (138) Developmental Biology (13360) Ecology (19866) Epidemiology (2067) Evolutionary Biology (24289) Genetics (15587) Genomics (22470) Immunology (17706) Microbiology (40314) Molecular Biology (17142) Neuroscience (88456) Paleontology (666) Pathology (2826) Pharmacology and Toxicology (4815) Physiology (7634) Plant Biology (15111) Scientific Communication and Education (2042) Synthetic Biology (4285) Systems Biology (9812) Zoology (2268)