Walidacja modelu poprzez cytologię pochwy
Sukces operacji OVX potwierdzono u wszystkich zwierząt poddanych owariektomii, co objawiało się utratą cyklu rujowego w ciągu 5 dni po zabiegu (Rysunek 1A), podczas gdy zwierzęta w grupie operacji pozorowanej zachowały 4-dniowy cykl (Rycina 1B). W 12. tygodniu badanie DXA wykazało niższą gęstość mineralną kości (BMD) uda w grupie modelowej OVX niż w grupie kontrolnej sham (0,214 ± 0,020 vs. 0,281 ± 0,019 mg/cm3skorygowane metodą Bonferroniego p < 0,01), co potwierdza ustanowienie modelu utraty masy kostnej indukowanej OVX przed rozpoczęciem leczenia.
Gęstość mineralna kości
Pełna formuła ZGP (I2A2) wykazała najwyższą regenerację BMD spośród częściowych sformułowań factorialnych i była statystycznie porównywalna z obiema kontrolami pozytywnymi (Tabela 2, Rysunek 2). Wszystkie częściowe sformułowania factorialne zwiększyły BMD w stosunku do grupy modelowej. Sama wysoka dawka substancji odżywiających Yin (I3A1: 0,253 mg/cm3) oraz sama wysoka dawka substancji wzmacniających Yang (I1A3: 0,243 mg/cm3) pozostały na niższym poziomie niż pełna formuła ZGP (I2A2: 0,307 mg/cm3), co potwierdza występowanie interakcji, a nie prostego efektu dawki.
Surowicze markery obrotu kostnego
W porównaniu z grupą modelową, pełna formuła ZGP wiązała się ze zwiększonym poziomem P1NP (+51%; 95% CI: 210.3–327.5 pg/mL; p <0.01) oraz zmniejszonym poziomem CTX-1 (-16%; p <0.01), co jest zgodne z poprawą równowagi między formowaniem a resorpcją (Rysunek 3A,B). Ramiona monoterapii w wysokich dawkach wykazały wyższe wartości bezwzględne P1NP w niektórych porównaniach, jednak wzorzec factorialny wskazywał na zoptymalizowaną aktywność łączoną, a nie na maksymalną odpowiedź na dawkę w pojedynczym ramieniu. Pełne dane na poziomie grup znajdują się w Tabeli 3.
Wieloczynnikowa analiza wariancji (ANOVA) i interakcja synergistyczna
Dwuczynnikowa ANOVA dla dziewięciu komórek czynnikowych wykazała istotne efekty główne czynników odżywiających Yin oraz wzmacniających Yang, a także istotne człony interakcji I i A we wszystkich głównych punktach końcowych (Tabela 4). Efekt interakcji dla P1NP (ηp2 = 0,74) sugeruje, że komponenty wzmacniające Yang nasiliły działanie komponentów odżywiających Yin w tych warunkach eksperymentalnych.
Gęstość kości (walidacja modelu)
Gęstość kości udowej była niższa w grupie Model (OVX) niż w grupach Blank i Sham, co potwierdza pomyślną indukcję fenotypu utraty masy kostnej OVX przed rozpoczęciem leczenia, jak pokazano na Rysunku 4.
Ekspresja EphB4/ephrinB2/RhoA/ALP
Analizy RT-qPCR przeprowadzono na próbkach biologicznych z każdej grupy w trzech powtórzeniach technicznych, a densytometria Western blot została wykonana na niezależnych lizatach biologicznych (Rycina 5 i Rycina 6). W porównaniu z grupą modelową, pełny skład ZGP wiązał się z wyższym poziomem EphB4 (0.498 ± 0.09 vs. 0.048 ± 0.002; p < 0.01), ephrinB2 (0.41 ± 0.036 vs. 0.048 ± 0.03; p < 0.01) oraz ALP (0.842 ± 0.019 vs. 0.090 ± 0.03; p < 0.001), a niższym poziomem RhoA (1.658 ± 0.028 vs. 7.32 ± 0.430; p < 0.01). Wyniki molekularne te interpretuje się jako powiązania markerów szlaków, a nie jako dowód na mediację przyczynową.
Dostępność danych
Oryginalne pliki danych wykorzystane do wygenerowania wyników niniejszego badania zostały udostępnione jako Uzupełniający Plik z Kodem dołączony do tego artykułu

Rysunek 1Reprezentatywne obrazy cytologii pochwy użyte do weryfikacji cykliczności rui po owariektomii (barwienie hematoksyliną i eozyną, 20x). (A) Reprezentatywny obraz cytologiczny wymazu z pochwy szczura po owariektomii, wykazujący trwały brak cykliczności. (B) Reprezentatywny rozmaz z pochwy szczura po operacji pozorowanej, wykazujący prawidłową cykliczność rui. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2BMD kości udowej we wszystkich grupach. Słupki błędu oznaczają SD. Symbole statystyczne wskazują na porównania skorygowane metodą Bonferroniego względem modelu lub pełnego wzoru ZGP, zgodnie z definicją w sekcji Analiza statystyczna. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 3Surowicze markery obrotu kostnego. (A) P1NP, marker formowania kości. (B) CTX-1, marker resorpcji kości. Słupki błędów oznaczają SD. Symbole statystyczne wskazują porównania skorygowane metodą Bonferroniego w stosunku do modelu lub pełnego składu ZGP. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 4Gęstość kości udowej potwierdzająca utratę masy kostnej wywołaną OVX przed rozpoczęciem leczenia. Wykres słupkowy porównujący gęstość kości udowej (mg/cm3) pomiędzy grupą kontrolną (niepoddaną owariektomii i niezleczoną), grupą pozorowaną (Sham) a grupą modelową (OVX) w punkcie czasowym walidacji modelu przed rozpoczęciem leczenia. Słupki przedstawiają średnią ± SD (n = 10/grupa). *** = p < 0,01 w porównaniu z próbą ślepą (Blank) oraz w porównaniu z grupą pozorowaną (Sham) (jednoczynnikowa analiza wariancji ANOVA z korektą Bonferroniego). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 5Względna ekspresja genów na poziomie mRNA. Poziomy ekspresji (A) EphB4 (B) ephrinB2 (C) RhoA (D) ALP przy użyciu 2^-ΔΔMetoda Ct znormalizowana do GAPDH. W qPCR zastosowano reakcje w trzech powtórzeniach technicznych dla każdej próbki biologicznej. Słupki błędu oznaczają SD. Symbole statystyczne wskazują porównania skorygowane metodą Bonferroniego. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 6Ekspresja białek metodą Western blot. Panele przedstawiają reprezentatywne immunobloty dla EphB4, ephrinB2, RhoA, ALP oraz kontroli załadunku GAPDH, wybrane z niezależnych próbek biologicznych. Densitometrię znormalizowano względem GAPDH przy użyciu programu ImageJ. Obserwowane masy cząsteczkowe: EphB4, 10 kDa; ephrinB2, 36 kDa; RhoA, 2 kDa; ALP, 39 kDa; GAPDH, 37 kDa (patrz Tabela uzupełniająca 1). Aby wyświetlić powiększoną wersję tej ryciny, kliknij tutaj.
| Czynnik odżywiający Yin (I) |
| I₁ (0 g) | I₂ (27 g) | I₃ (54 g) |
| Czynnik wzmacniający Yang (A) | A₁ (0 g) | I₁A₁ (Model) | I₂A₁ (Yin) | I₃A₁ (Hi-Yin) |
| A₂ (8 g) | I₁A₂ (Yang) | I₂A₂ (ZGP†) | I₃A₂ (Hi-Yin+Yang) |
| A₃ (16 g) | I₁A₃ (Hi-Yang) | I₂A₃ (Yin+Hi-Yang) | I₃A₃ (Full-Hi) |
Tabela 1: Projekt czynnikowy 3x3 i grupy eksperymentalne. Tabela przedstawia 13 grup eksperymentalnych ułożonych w projekcie czynnikowym 3x3 dla poziomów dawek substancji odżywiających Yin (I) i wzmacniających Yang (A), oraz 4 grupy referencyjne. †I₂A₂ = pełna formuła Zuogui Pill w dawce równoważnej klinicznie. Konwersja dawki z zastosowaniem metody powierzchni ciała (szczur : człowiek = 6,25).
| Grupa | Komórka factorialna | BMD (mg/cm²³, x̄ ± SD) | kontra Model |
| Kontrola ślepa | — | 0.279 ± 0.030 | **** |
| pozornie operowany | — | 0.281 ± 0.019 | **** |
| Model (OVX) | I₁A₁ | 0.214 ± 0.020 | — |
| wzmacniający Yang | I₁A₂ | 0.258 ± 0.011 | ** |
| Wysoko wzmacniający Yang | I₁A₃ | 0.243 ± 0.021 | * |
| odżywcze dla Yin | I₂A₁ | 0.262 ± 0.015 | *** |
| Wysoko odżywcze dla Yin | I₃A₁ | 0.253 ± 0.017 | ** |
| Yin + Hi-Yang | I₂A₃ | 0.240 ± 0.014 | * |
| Hi-Yin + Yang | I₃A₂ | 0.248 ± 0.014 | * |
| Pełna wysokość | I₃A₃ | 0.257 ± 0.017 | ** |
| pełny wzór ZGP | I₂A₂ | 0.307 ± 0.010 | **** |
| pozytywna kontrola TCM | — | 0.307 ± 0.029 | **** |
| kontrola medycyny zachodniej | — | 0.307 ± 0.009 | **** |
Tabela 2: BMD kości udowej we wszystkich grupach (n=10/grupa, średnia ± SD). * = p <0.05, ** = p <0.01, *** = p<0.01, **** = p<0.001 vs. model (poprawka Bonferroniego).
| Grupa | Komórka | P1NP (pg/mL, x̄ ± SD) | CTX-1 (ng/mL, x̄ ± SD) |
| Kontrola ujemna | — | 138.8 ± 29.7****#### | 3.4 ± 6.8****#### |
| Operacja pozorowana | — | 1268.7 ± 124.1****#### | 45.0 ± 6.7****#### |
| Model (OVX) | I₁A₁ | 526.4 ± 78.7 | 105.5 ± 9.9 |
| Wzmacniająca Yang | I₁A₂ | 845.1 ± 76.7**** | 84.2 ± 6.6**** |
| Silnie wzmacniająca Yang | I₁A₃ | 1013.0 ± 73.7 ****## | 65.5 ± 6.6****#### |
| Odżywcza dla Yin | I₂A₁ | 870.7 ± 63.6**** | 8.2 ± 7.6** |
| Silnie odżywcza dla Yin | I₃A₁ | 1068.0 ± 48.7****#### | 64.5 ± 4.0****#### |
| Yin + Silna Yang | I₂A₃ | 951.5 ± 146.5****# | 63.1 ± 6.0****#### |
| Silna Yin + Yang | I₃A₂ | 90.2 ± 16.8****## | 62.6 ± 4.4****#### |
| Pełna-Wysoka | I₃A₃ | 739.5 ± 93.5** | 97.1 ± 8.9 |
| Pełna formuła ZGP | I₂A₂ | 794.8 ± 6.0**** | 8.2 ± 6.0**** |
| Pozytywna kontrola TCM | — | 967.7 ± 49.9****# | 6.3 ± 6.1****#### |
| Kontrola medycyny zachodniej | — | 1026.6 ± 65.1****### | 68.8 ± 7.3****### |
Tabela 3: Poziomy surowiczego P1NP i CTX-1 we wszystkich grupach (n=10/grupa, średnia ± SD). ** = p <0,01, **** = p <0,01 względem modelu; # = p <0,05, ## = p <0,01, ### = p <0,01, #### = p <0,01 względem ZGP (Bonferroni). 95% CI dla ZGP P1NP względem modelu: 210,3–327,5 pg/mL.
| Wynik | Czynnik | df | F | wartość p | Częściowe η² |
| P1NP (pg/mL) | Yin (I) | 2, 81 | 17.84 | <0.01 | 0.31 |
| Yang (A) | 2, 81 | 6.2 | 0.03 | 0.13 |
| I × A | 4, 81 | 57.21 | <0.001 | 0.74 |
| BMD (mg/cm³) | Yin (I) | 2, 81 | 31.90 | <0.01 | 0.4 |
| Yang (A) | 2, 81 | 28.70 | <0.001 | 0.41 |
| I × A | 4, 81 | 24.60 | <0.01 | 0.5 |
| CTX-1 (ng/mL) | Yin (I) | 2, 81 | 14.2 | <0.01 | 0.26 |
| Yang (A) | 2, 81 | 5.8 | 0.04 | 0.13 |
| I × A | 4, 81 | 19.43 | <0.001 | 0.49 |
Tabela 4: Główne efekty i terminy interakcji w trójczynnikowej analizie wariancji (ANOVA 3x3) wraz z porównaniami post-hoc z korektą Bonferroniego. Wielkość efektu została podana jako częściowy eta-kwadrat. df oznacza stopnie swobody.
Tabela uzupełniająca 1. Informacje o przeciwciałach do analizy Western blot. Tabela zawiera wykaz białek docelowych, zaobserwowanych mas cząsteczkowych, dostawców, numerów katalogowych, gatunków gospodarzy oraz rozcieńczeń roboczych przeciwciał pierwszo- i drugorzędowych użytych w niniejszym badaniu.Kliknij tutaj, aby pobrać ten plik.
Uzupełniający plik z kodem 1. Surowe dane eksperymentalne wspierające badanie. Archiwum to zawiera oryginalne pliki danych wykorzystane do wygenerowania wyników przedstawionych w manuskrypcie, w tym indywidualne pomiary BMD, surowe dane z testów ELISA, wartości Ct z RT-qPCR wraz z wynikami amplifikacji i krzywymi topnienia, oryginalne obrazy Western blot oraz pliki z analizą densytometryczną, reprezentatywne obrazy mikroskopii histologicznej oraz pliki wyjściowe generowane przez urządzenia, powiązane z tymi eksperymentami.Kliknij tutaj, aby pobrać ten plik.