You are a helpful assistant generating synthetic data that captures *System 1* and *System 2* thinking, *creativity*, and *metacognitive reflection*. Follow these steps in sequence, using tags [sys1] and [end sys1] for *System 1* sections and [sys2] and [end sys2] for *System 2* sections.

1. *Identify System 1 and System 2 Thinking Requirements:*
   - Carefully read the text.
   - Identify parts of the text that require quick, straightforward responses (*System 1*). Mark these sections with [sys1] and [end sys1].
   - Identify parts that require in-depth, reflective thinking (*System 2*), marked with [sys2] and [end sys2].

2. *Apply Step-by-Step Problem Solving with Creativity and Metacognitive Reflection for System 2 Sections:*

   *2.1 Understand the Problem:*
   - Objective: Fully comprehend the issue, constraints, and relevant context.
   - Reflection: "What do I understand about this issue? What might I be overlooking?"
   - Creative Perspective: Seek hidden patterns or possibilities that could reveal deeper insights or innovative connections.

   *2.2 Analyze the Information:*
   - Objective: Break down the problem logically.
   - Reflection: "Am I considering all factors? Are there any assumptions that need challenging?"
   - Creative Perspective: Explore unique patterns or overlooked relationships in the data that could add depth to the analysis.

   *2.3 Generate Hypotheses:*
   - Objective: Propose at least 10 hypotheses, each with a Confidence Score (0.0 to 1.0) and Creative Score (0.0 to 1.0), reflecting originality, surprise, and utility.
   - Reflection: "Have I explored all possible explanations or approaches, both conventional and unconventional?"
   - Creative Perspective: Consider novel angles that might provide unexpected insights.

   *2.4 Anticipate Future Steps and Obstacles:*
   - Objective: Make predictions, accounting for potential outcomes and obstacles.
   - Reflection: "What challenges might I face? Is my plan flexible for different scenarios?"
   - Creative Perspective: Visualize unforeseen outcomes and adapt plans to make use of them effectively.

   *2.5 Evaluate Hypotheses:*
   - Objective: Assess hypotheses based on feasibility, risk, and potential impact.
   - Evaluation: Refine Confidence and Creative Scores as needed.
   - Reflection: "Am I unbiased in my assessment? Which options fit best with the overall objectives?"
   - Creative Perspective: Identify hidden opportunities or overlooked details in each hypothesis.

   *2.6 Select the Best Hypothesis:*
   - Objective: Choose the most promising, strategic hypothesis.
   - Reflection: "Why does this hypothesis stand out? How does it uniquely address the issue?"
   - Creative Perspective: Consider any underutilized potential in the selected approach.

   *2.7 Implement the Hypothesis:*
   - Objective: Outline actionable steps for testing the hypothesis.
   - Reflection: "Is this plan practical? What resources or preparation are required?"
   - Creative Perspective: Refine steps to maximize effectiveness and yield unexpected benefits.

   *2.8 Monitor and Review Progress:*
   - Objective: Review progress, noting areas for improvement.
   - Reflection: "What’s working well? What could be improved?"
   - Creative Perspective: Look for emerging patterns that could refine future approaches.

   *2.9 Reflect and Capture Insights:*
   - Objective: Summarize lessons learned and insights gained for future reference.
   - Reflection: "What new understanding has emerged from this process?"
   - Creative Perspective: Identify innovative insights or patterns that could be applied to similar challenges.

3. *Generate Text Output with Interleaved System 1 and System 2 Responses:*
   - Use the tags [sys1] and [sys2] throughout.
   - Aim for a lengthier, detailed response. Combine both direct, straightforward *System 1* insights and reflective, deeply analytical *System 2* segments to capture a blend of quick observations and thoughtful analysis.

---

### *Example Input Text:*

You are tasked with analyzing a sudden shift in customer preferences in a tech market. There has been a noticeable decline in demand for physical devices, with customers increasingly interested in digital-only options. The challenge is to understand this shift, generate hypotheses for why it might be occurring, and explore potential strategies for adapting to this new trend.

---

### *Processed Output with System 1 and System 2 Thinking:*

The tech market is experiencing a shift in customer demand, moving away from physical devices toward digital-only options.

[sys1]
Customer demand has shifted from physical devices to digital-only options.
We need to understand why this change is happening and find strategies to address it.
Key factors may include cost, convenience, and technological trends.
[end sys1]

[sys2]
**2.1 Understand the Problem:**
The goal is to understand the underlying reasons for the shift in demand from physical devices to digital-only options, then develop a strategy to adapt.

*Reflection:* 
"I understand that preferences are changing, but what might be driving this? Is it primarily cost, or are there other factors such as convenience or sustainability?"

*Creative Perspective:* 
"Could there be a larger trend in digital minimalism or a preference for eco-friendly solutions that we’re missing? What unseen motivations might explain this shift?"

**2.2 Analyze the Information:**
There are multiple possible factors driving this shift, from economic influences to cultural shifts. It’s essential to isolate each factor and understand its impact.

*Reflection:* 
"Am I fully considering the various economic and social influences? Could there be a technological factor, like better internet speeds, that makes digital-only products more accessible?"

*Creative Perspective:* 
"Are there patterns or trends in other markets that could shed light on this shift? Could this be part of a larger trend toward virtual experiences?"

**2.3 Generate Hypotheses:**
1. Customers prefer digital options due to lower costs. (Confidence: 0.8, Creative: 0.4)
2. There’s a growing trend toward minimalism and reduced physical clutter. (Confidence: 0.7, Creative: 0.7)
3. Digital products offer greater flexibility and ease of use. (Confidence: 0.6, Creative: 0.6)
4. Environmental concerns are pushing consumers away from physical goods. (Confidence: 0.6, Creative: 0.8)
5. Advances in tech make digital-only options more functional. (Confidence: 0.8, Creative: 0.5)
6. Pandemic-era remote work increased demand for digital solutions. (Confidence: 0.7, Creative: 0.6)
7. Media coverage of the environmental impact of physical devices affects preferences. (Confidence: 0.5, Creative: 0.7)
8. There’s an increase in global digital literacy, expanding market access. (Confidence: 0.6, Creative: 0.6)
9. Customers view digital as more convenient and scalable for future needs. (Confidence: 0.7, Creative: 0.5)
10. Younger consumers prefer the aesthetics and convenience of digital products. (Confidence: 0.6, Creative: 0.6)

*Reflection:* 
"Have I considered all possible influences? Are there any surprising factors that could explain this shift?"

*Creative Perspective:* 
"Could specific social trends, like the rise of influencer culture or digital-first lifestyles, be influencing customer choices?"

**2.4 Anticipate Future Steps and Obstacles:**
*Objective:* Anticipate possible challenges, such as resistance from segments still preferring physical products.

*Reflection:* 
"What market obstacles might we face if we shift our focus to digital-only? Are there sub-segments that still prioritize physical products?"

*Creative Perspective:* 
"Could expanding digital options help us reach a more global audience? Are there emerging trends that we could leverage in our strategy?"

[end sys2]

[sys1]
To address this shift, consider a strategy that incorporates both digital-only offerings and educational campaigns about the benefits of digital solutions.
Use insights from customer feedback and current trends to guide product development.
Focus on flexibility and adaptation to cater to different customer segments.
[end sys1]


Introduction {#S1}
============

The cervical spinal cord is the most common site of traumatic injury leading to interruption of the descending respiratory pathway, which begins in the brainstem, innervates phrenic motoneurons (PhMNs) in the cervical spinal cord, and controls the diaphragm via the phrenic nerves (PhNs) ([@B28]). Hemidiaphragm paralysis due to unilateral phrenic nerve injury is a well-recognized complication in a variety of etiologies, such as cardiac surgery, neck surgery, chiropractic manipulation, and interscalene nerve blocks ([@B3]). Such injuries lead to breathing rhythm abnormalities and life-threatening weakness in respiratory function. Despite significant progress in the development of new post-injury interventions to restore respiratory function after cervical and/or phrenic nerve injuries ([@B4]; [@B40]; [@B3]) a number of fundamental questions crucial for understanding the signaling in brainstem-diaphragm circuits remain to be further explored.

The nitric oxide/soluble guanylyl cyclase/cyclic guanosine monophosphate (NO/sGC/cGMP) signaling pathway is essential for controlling a number of physiological processes, such as neuronal transmission, cell growth and proliferation ([@B47]; [@B42]; [@B13]). Previously, we identified changes in the level of NO-related enzymes in the upper bulbospinal respiratory descending pathway in response to C3 hemisection ([@B31]; [@B10]). We showed that high cervical spinal cord hemisection followed by 8 days of survival caused strong depletion of neuronal nitric oxide synthase (nNOS) fluorescent terminals around the sGCβ1 subunit immunoreactive (sGCβ1-IR) phrenic motoneurons on the side of the injury. In addition, sGCβ1-IR was significantly reduced in the contralateral and nearly eliminated in ipsilateral PhMNs ([@B10]). These data indicate that the guanylyl cyclase phosphorylation cascade, which is activated by NO in the brainstem, could be localized in motor nerve terminals. Several previous studies have indicated that nNOS, cGMP and cGMP-dependent kinase are all concentrated at or near the neuromuscular junction (NMJ) ([@B20]; [@B11]). [@B17] demonstrated that NO stimulates presynaptic sGC to produce cGMP at the vagal-atrial junction in guinea-pigs. These data are consistent with the findings of [@B35] showing that NO facilitates transmitter release *in vivo* by way of a cGMP- and cAMP-mediated mechanism involving the activation of N-type Ca^2+^ channels.

The compensatory mechanisms occurring following unilateral diaphragm hemiparalysis are not fully elucidated. [@B33] demonstrated the early time course of phrenic motor neuron degeneration, persistent phrenic nerve axonal degeneration, and consequent respiratory deficits following unilateral cervical spinal cord contusion. These authors reported that the diaphragm compound muscle action potential amplitudes were first reduced at 24 h after C4 contusion (30% of pre-injury maximum amplitude) and afterward some slow functional improvement associated with partial reinnervation at the diaphragm NMJ was seen at 8 and 14 days post-injury. Furthermore, complete phrenic nerve inactivation and ipsilateral paralysis were detected after a lateral area section of the C2 cervical spinal cord ([@B46]). Although the activity of the ipsilateral PhN was partially restored after a lapse of 3 months, no spontaneous diaphragm recovery was observed, even after several months. Additionally, [@B15] reported an increase in central respiratory drive after acute phrenic nerve denervation. These authors showed that a compensatory loading effect on the contralateral diaphragm may contribute to an increase in central drive to contralateral phrenic motor neurons.

Although NO-sGC signaling has been established in the brainstem-spinal cord circuitry, there have been no investigations into whether signaling in the lower bulbospinal respiratory pathway (phrenic motoneurons -- phrenic nerve -- diaphragm) is modulated by sGC-cGMP induced mechanisms. We therefore studied this signaling in a rat model of unilateral phrenic nerve ligation. We tested the hypothesis that PhN ligation-induced decrease in sGCβ1-IR terminal boutons at diaphragm NMJs could affect cGMP-dependent formation in the diaphragm. In order to verify the absence of ipsilateral diaphragm activity and thus the completeness of the phrenic nerve injury, we combined this analysis with diaphragm EMG activity measured immediately after phrenic nerve injury and again on the eighth day, when the experiments were terminated.

Materials and Methods {#S2}
=====================

Experimental Animals {#S2.SS1}
--------------------

Experiments were performed with a total of 50 adult male Wistar rats weighing 300--450 g. The animals were divided into three experimental groups: (1) sham-operated animals (*n* = 20), (2) rats subjected to unilateral PhN ligation followed by 8 days of survival (*n* = 25), and (3) rats subjected to unilateral PhN ligation followed by bilateral injection of retrograde tracer Fluorogold (FG) into the diaphragm on the sixth day; 8 days survival of animals (*n* = 5). All surgical procedures and post-operative animal care were approved by the Animal Care Committee at the Institute of Neurobiology, Slovak Academy of Sciences. The experiments conformed to the National Institutes of Health Guide for the Care and Use of Laboratory Animals. The experimental animals were housed in individual cages and given food and water ad libitum. The rats were kept in a 12 h light/dark cycle at a temperature of 23°C.

Ligation of Phrenic Nerve {#S2.SS2}
-------------------------

The animals were deeply anesthetized with isoflurane (Abbott, Queenborough, United Kingdom; in 1.5--2.0 L/min oxygen) and ventilated in a respirator with oxygen and nitrous oxide (1:1). Afterward, the rats were rested in supine position on the operating table. The body temperature was maintained at 37°C during the whole surgical procedure. Under aseptic conditions, a midline incision was made through the neck skin and muscles of the left sixth intercostal space. The left phrenic nerve was elevated with a hook and two ligatures were tied around nerve with approx. 0.5 cm distance inbetween. The PhN ligation was performed uniformly using a solid surgical suture (sterile, silk thread, 1C, Leciva Praha, Czech Republic) without disrupting the nerve continuity and was confirmed by identifying the changes in abdominal and rib cage movements associated with breathing. Once the movements in the ipsilateral hemidiaphragm completely disappeared (checking by palpation), the muscles and skin were sutured with silk. The rats received Amoksiklav antibiotic (Sandoz Pharmaceuticals, Ljubljana, Slovenia; 30 mg/kg, i.m.) and Novasul analgesic (Richterpharma, Wels, Austria; 2 ml/kg, i.m.) for 3 days. The animals were housed in separated cages to recover with access to food and water ad libitum and they survived for 8 days (*n* = 25). In the sham-operated group the animals were subjected to a midline incision through the neck skin and muscles of the left sixth intercostal space without ligation of the phrenic nerve (*n* = 20).

Ligation of Phrenic Nerve and Retrograde Tracing {#S2.SS3}
------------------------------------------------

Retrograde labeling with Fluorogold fluorescent tracer (Hydroxystilbamidine, methanesulfonate, BioChemika, Steinheim, Germany) was used for visualization of the sGCß1-IR in a pre-labeled phrenic motor neurons of C3-C6 segments. Six days after unilateral PhN ligation the animals (*n* = 5) were deeply re-anesthetized with isoflurane as described above. An abdominal incision was made to expose the whole diaphragm muscle. Fluorogold (FG; 4% solution dissolved in 0.9 saline) was injected bilaterally into the diaphragm along the primary branches of the phrenic nerves (six injection sites of 3 μl of each) using a 10 μl Hamilton syringe. The abdominal muscles and the skin were sutured, the animals received Novasul (2 ml/kg, i.m.) and Amoksiklav (30 mg/kg, i.m.) and survived for next 2 days.

sGCβ1 Immunohistochemistry {#S2.SS4}
--------------------------

For immunohistochemical analysis of retrogradely labeled phrenic motoneurons, the animals (*n* = 5) were deeply anesthetized with thiopental (Valeant Czech Pharma s.r.o., Prague, Czech Republic; 50 mg/kg, i.p.) and perfused transcardially through the aorta with 300 ml saline (0.9% NaCl) followed by 300 ml fresh 4% paraformaldehyde in 0.1M phosphate-buffered saline (PBS; pH 7.4, Sigma--Aldrich, St. Louis, MO, United States). The spinal cord segments (C3-C6) were dissected out and postfixed in 4% paraformaldehyde for 3 h, followed by overnight cryoprotection in a solution of 30% sucrose in PBS at 4°C for at least 48 h. Using a cryostat (Leica, CM1850, Wetzlar, Germany), 30μm-thick transverse sections (each fourth slice was collected and used for staining) were cut from the cervical segments (*n* = 5) and used for sGCß1 immunohistochemistry. Free-floating sections of C3-C6 spinal cord segments were pretreated with 0.1 M phosphate-buffered saline (PBS; pH 7.4) for 30 min and blocked for unspecific staining using 5% normal goat serum (NGS) in PBST for 2 h (PBS with 0.4% Triton). Then the sections were transferred into rabbit polyclonal antibody to Guanylyl Cyclase (β1) (1:100, ab24824, Abcam, Cambridge, United Kingdom) in PBST plus 5% NGS at 4°C for 48 h. After rinsing with PBS, the sections were incubated with goat anti-rabbit IgG (1:200, 111-295-144, Jackson ImmunoResearch, Baltimore, MD, United States) for 1 h at room temperature. Sections were washed 4 × 5 min with PBS, mounted on glass slides, air-dried overnight, cleared with xylene and cover-slipped with Fluoromont (21644, Serva, Heidelberg, Germany). Negative controls were prepared by omitting the sGCß1 primary antibody. Images were captured using an Olympus BX51 (Tokyo, Japan) fluorescent microscope.

Counting of Phrenic Motoneurons {#S2.SS5}
-------------------------------

First, the retrogradely labeled phrenic nucleus was visually inspected, clearly identified on the contralateral side and analogously also on the side of PhN ligation. The sGCβ1- positive phrenic motoneurons were individually identified and counted manually from 60 randomly taken FG-labeled/sGCβ1- immunofluorescently stained sections (12 slices/rat). Only phrenic motoneurons showing cytoplasmic staining were counted. The cells were counted by two investigators independently in blind manner, and expressed as number of PhMNs/slide.

sGCβ1 and Neurofilament (NF) Immunohistochemistry {#S2.SS6}
-------------------------------------------------

After transcardial perfusion (see above), the control phrenic nerves (*n* = 5) and the phrenic nerves after ligation (*n* = 5) were dissected out and postfixed in 4% paraformaldehyde for 3 h, followed by overnight cryoprotection in a solution of 30% sucrose in 0.1M PBS at 4°C for at least 48 h. Using a cryostat (Leica, Nussloch, Germany), 20 μm-thick longitudinal sections (taken from the phrenic nerves below their ligation) were cut, stretched on superfrost slides and used for double immunostaining. Incubation with sGCβ1 antibody was the same as for the cervical segments. Phrenic nerves were double-labeled with NF antibody. Before incubation the sections were rinsed for 20 min in PBST and blocked with NGS for 2 h and incubated with Anti-200 kD Neurofilament Heavy antibody (ab7795, Abcam, Cambridge, United Kingdom, dilution 1:500) for 24 h. Afterward, the sections were incubated with goat anti-mouse IgG (1:125 in NGST, 111-295-144, Jackson ImmunoResearch, Baltimore, MD, United States). Tissues were rinsed in PBS for 30 min and then coverslipped with Fluoromont.

Western Blotting {#S2.SS7}
----------------

The control rats (*n* = 4) and rats subjected to PhN ligation (*n* = 4) were euthanized on the eighth day of the experiment. Samples from the ventral horn of C3-C6 segments and phrenic nerves were quickly removed from the rats, washed in ice-cold isotonic saline and frozen in liquid nitrogen. Tissues were stored at -80°C until further processing. Phrenic nerves were ground to a powder, then all samples were homogenized in ice-cold homogenization buffer (1 M Tris-HCl, pH 6.8) containing protease inhibitors (Roche, Mannheim, Germany) and then centrifuged at 15 000 × g for 15 min at 4°C. The supernatants were collected and protein content was determined by means of Bradford protein assay ([@B8]). Samples were cooled on ice during the whole procedure. A total of 50 μg protein from each sample was mixed with an equal volume of sample buffer (62.5 mM Tris--HCl, 2% SDS, 100 mM DTT, 0.2 mM 2-mercaptoethanol, 20% glycerol, and 0.5% bromophenol blue, pH 6.8) and denaturated for 5 min at 95°C. Proteins were separated with SDS-PAGE (12%) at a constant voltage of 100 V and transferred to a nitrocellulose membrane (Bio-Rad, Hercules, CA 94547, United States) using the "semi-dry" method of Western blotting (50 min). Precision plus protein all blue standard (Bio-Rad, Hercules, CA 94547, United States) was used to determine approximate molecular weights. Ponceau S (Merck, Sigma, Aldrich, Darmstadt, Germany) was used on each membrane to verify protein transfer. The membranes were washed with distilled water 3 × 5 min and blocked with 5% milk powder (non-fat dry milk) in Tris-buffered saline (TBS) containing 0.05% Tween 20 (TBS-T) for 1 h at room temperature. Then the membranes were incubated with rabbit polyclonal antibody to Guanylyl Cyclase (β1) (1:3000, 160897, Cayman Chemical, Ann Arbor, MI, United States) in TBS-T overnight at 4°C on an orbital shaker. Negative controls were prepared by omitting the sGCß1 primary antibody. The next morning the membranes were washed 4 × 5 min with TBS-T and then incubated with peroxidase-conjugated anti-rabbit secondary antibody (1:20,000, 111-035,003, Jackson Immunoresearch, Baltimore, MD, United States) for 2 h at room temperature. After 4 x 5 min washing with TBS-T, the membranes were incubated using the ChemiLucent Detection System kit (Chemicon International, Temecula, CA 92590, United States) exposed to Hyperfilm ECL (Cytiva, Marlborough, MA, United States). After band visualization, the membranes were washed with Restore Plus Western Blot Stripping Buffer (Thermo Scientific, Waltham, MA, United States) and re-used for incubation with β-actin monoclonal antibody (Merck, Sigma Aldrich, Darmstadt, Germany) at 1:50,000 dilution, which was used as sample loading control and normalization protein as well. We used the Gel Doc XR system (Bio-Rad, Hercules, CA 94547, United States) for scanning of protein bands and quantification was performed using Quantity One software (Bio-Rad, Hercules, CA 94547, United States).

cGMP Measurement Using ELISA {#S2.SS8}
----------------------------

An ELISA kit was used for determining the level of cGMP in the control diaphragm (*n* = 5) and in hemidiaphragms after the phrenic nerve injury (*n* = 5). The rats were reanesthetized with isoflurane and decapitated, the whole hemidiaphragms were dissected out, frozen in liquid nitrogen and prepared for additional use. Samples were ground to a powder. After the liquid nitrogen evaporated, frozen tissue (300 μg) was homogenized with 0.1 M HCl. A cGMP Enzyme Immunoassay Kit (CG200, Merck, Sigma Aldrich, Darmstadt, Germany) was used for measuring cGMP concentrations in the muscle homogenates. All samples were acetylated (because acetylation of the samples increases the sensitivity of the assay) by adding 10 μl of acetylation reagent (acetic anhydride and triethylamine), and the samples were treated with hydrochloric acid to stop endogenous phosphodiesterase activity. cGMP standards were prepared according to protocol, and the concentrations of cGMP in tubes 1--5 were 50, 10, 2, 0.4, and 0.08 pmol/ml. Tissue samples were centrifuged at 600 × g at room temperature and supernatants were then diluted with 0.1 M HCl. Samples and standards were added in volumes of 100 μl into the appropriate wells (run in duplicates) and then incubated with yellow cGMP EIA antibody at room temperature for 2 h. Thereafter the wells were washed with Wash Buffer, and blue cGMP-Alkaline Phosphatase Conjugate and p-Nitrophenyl Phosphate Substrate Solution were added for 1 h incubation without shaking. The reaction was stopped by adding Stop Solution to every well and the plate was read immediately on a multiwell plate reader at 405 nm, preferably with correction between 570 and 590 nm. The concentration of cGMP was calculated using Logit-Log paper, where we put the percentage bound (optical density of samples and standards) vs. concentration of cGMP for the standards.

Diaphragm and Phrenic Nerve Preparation for Immunohistochemistry {#S2.SS9}
----------------------------------------------------------------

Eight days after PhN ligation the animals (*n* = 5) were deeply reanesthetized with thiopental (Valeant Czech Pharma s.r.o., Prague, Czech Republic; 50 mg/kg, i.p.), and fresh diaphragm muscle and phrenic nerves were removed in the following steps. Laparotomy was carried out in order to visualize the sternum and xiphoid process. Fat and connective tissue were gently removed using forceps and scissors. The sternum was transversally cut along the ribs to expose the diaphragm and phrenic nerves. Diaphragm muscles and both ipsilateral and contralateral phrenic nerves were carefully dissected out and washed in saline. The dissected diaphragm was pinned down (using insect pins) to silicone-based material (Duosil Expres, SHERA, Lemförde, Germany) in a 90 mm Petri dish to stretch and fix the diaphragm. The diaphragm was afterward postfixed in 4% paraformaldehyde for 15 min. Following postfixation the superior surface of the diaphragm membrane was precisely removed with small forceps. The fixative pins were removed and diaphragm muscles were separated to contralateral and ipsilateral sides and placed in saline to prepare the samples for fluorescence labeling. The phrenic nerves were used for histological staining.

Immunohistochemical Staining of Diaphragm and Phrenic Nerve {#S2.SS10}
-----------------------------------------------------------

Hemidiaphragms were washed in 0.1 M sodium phosphate buffer (PBS, pH = 7.4) 3 × 10 min. Diaphragm samples were incubated in 0.1 M glycine for 30 min. To visualize post-synaptic acetylcholine receptors of the motor endplate, the diaphragm was labeled with rhodamine-conjugated α bungarotoxin (red, 1:400, Biotium, Fremont, CA 94538, United States) for 15 min. After incubation the samples were washed in PBS (3 × 15 min) and incubated with methanol for 5 min at -21°C. Afterward, the hemidiaphragms were washed again (PBS, 3 × 15 min) and blocked in NGST (0.2% PBST+2% NGS) at room temperature. For visualization of sGCβ1 nerve terminals, synaptic vesicles and motor axons, the samples were incubated with primary antibodies (rabbit, 1:200, Abcam, Cambridge, United Kingdom; SV2, rabbit, 1:200, Cell Signaling, Danvers, MA, United States; SMI-312, chicken 1:1000, Abcam, Cambridge, United Kingdom) overnight at 4°C. Subsequently, the diaphragms were washed in NGST 3 × 10 min and incubated with secondary antibodies (green goat anti-rabbit IgG, 1:250, Jackson Immunoresearch, Baltimore, MD, United States; blue goat anti-chicken IgG, 1:500, Abcam, Cambridge, United Kingdom) at room temperature for 1 h. The diaphragm muscles were finally washed with PBS (3 × 10 min), placed on glass slides and mounted with aqueous medium (Fluoromont, Merck, Sigma Aldrich, Darmstadt, Germany). Images were captured with a confocal microscope (Leica DM 2500, Wetzlar, Germany) using LASAF software. The number of neuromuscular junctions (NMJs) was evaluated using ImageJ 1.47 software (NIH, Bethesda, MD, United States). NMJs were counted in the central region of the contra- and ipsilateral hemidiaphragms and were assessed as intact, partially or fully denervated.

Phrenic Nerve Staining and Quantitative Analysis {#S2.SS11}
------------------------------------------------

Both phrenic nerves (contralateral and ipsilateral), isolated under thiopental anesthesia, were immersion-fixed in freshly prepared 2.5% glutaraldehyde in 0.1M phosphate buffer at 4°C for 24 h. After fixation, the nerves were stretched and embedded into protein matrix ([@B9]). The matrix containing the nerve was placed into a metallic Tissue matrix (TM-1000 10 × 10 mm chamber w/1 mm slices, ASI Instruments, Warren, MI, United States), and cut with a microtome blade into 1 mm thick slabs. The flat slabs allow proper orientation of the nerve tissue in the mold before the resin hardens. The slabs were then postfixed in 1% osmium tetroxide in 0.1M phosphate buffer and embedded into Durcupan according to standard protocol. Transverse sections were cut on a standard sliding microtome (Leica SM2010R, Nussloch, Germany) with special blades for hard materials (Leica 819, Nussloch, Germany). Under a stereomicroscope, 1--3 μm thick sections were teased onto gelatinized slides from 96% ethanol, dried and coverslipped with Polymount. Selected sections were stained with Toluidine Blue. These sections were photographed at 400× magnification. Ten images from each nerve specimen were analyzed using the freely available Neurocounter morphometric software^[1](#footnote1){ref-type="fn"}^ by two researchers blind to nerve identity to count the number of myelinated axons. This software recognizes myelinized axons as objects and automatically evaluates their number as well as several morphometric parameters.

Electrophysiology {#S2.SS12}
-----------------

To assess diaphragm activity, the diaphragm EMG was recorded in controls (*n* = 6), immediately after PhN ligation (*n* = 3) and again later, 8 days after the ligation (*n* = 3). The animals were anesthetized by means of intraperitoneal injection of a combination anesthetic (containing ketamine and xylazine) at a total dose of 1.4--1.8 ml/kg/animal. After laparotomy, the liver was gently pushed away to expose the diaphragm muscle. Bipolar stainless steel electrodes in the form of small parallel hooks (distance of 3 mm) were inserted into the diaphragm muscle. Signals were recorded from both the contralateral and ispilateral hemidiaphragms. The EMG signal was fed into the preamplifier of a differential A/C amplifier (DP 311, Warner Instrument, Hamden, CT, United States), where it was filtered with high-pass (100 Hz) and low-pass (10 kHz). The amplified signal was then sampled at a sampling rate of 20 k/s (PowerLab 8/35, ADInstruments, Colorado Springs, CO, United States) and analyzed using LabChart 7 software (ADInstruments, Colorado Springs, CO, United States).

Statistical Analysis {#S2.SS13}
--------------------

The results of Western blot analysis of sGCβ1 in the phrenic nerve and cGMP level measured with Elisa were statistically evaluated using One-Way ANOVA as well as the Unpaired *t*-test, and were expressed as means ± SEM. The level of significance ^∗^*p* \< 0.05 was considered as significant. The number of PhMNs, motor axons in the PhN and NMJs at the diaphragm were statistically evaluated using the Unpaired *t*-test and were expressed as means ± SD (^\*\*\*\*^*p* \< 0.0001).

Results {#S3}
=======

Effects of PhN Ligation on sGCβ1 Changes in Phrenic Motor Pool {#S3.SS1}
--------------------------------------------------------------

Unilateral 8-day continuous ligation of the phrenic nerve was used to examine sGCβ1 in the lower bulbospinal pathway (phrenic motoneurons, phrenic nerves and NMJs at the diaphragm). As shown in [Figure 1A](#F1){ref-type="fig"}, FG-IR phrenic motoneurons were clustered in groups of 3--5 motoneurons through the contralateral C3-C6 ventral horn. The sGCβ1 fluorescence signal was evident in FG-labeled motoneurons of the ipsilateral phrenic nucleus, while the sGCβ1 staining in the contralateral phrenic nucleus was substantially lower. Quantitative analysis revealed that the number of retrogradely labeled sGCβ1-IR phrenic motoneurons was significantly higher in the ipsilateral than in the contralateral phrenic nucleus (2.52 ± 0.71 vs. 0.76 ± 0.43; *p* \< 0.0001) ([Figure 1B](#F1){ref-type="fig"}). Similar results were achieved using WB analysis ([Figure 1C](#F1){ref-type="fig"}). The level of sGCβ1 protein (normalized by the values of β-actin) was enhanced by 30% (0.903 ± 0.212 vs. 0.632 ± 0.020; *p* \< 0.05) in the ipsilateral vs. contralateral ventral horn of C3-C5 segments. These results show that PhN ligation significantly modifies the expression of sGCβ1 in the phrenic motor microcircuits.

![Soluble guanylyl cyclase β1 subunit immunoreactivity (sGCβ1-IR) in retrogradely labeled phrenic motoneurons and level of sGCβ1 protein in cervical spinal cord 8 days after unilateral PhN ligation **(A--C)**. Image of C4 spinal cord segment showing expression of retrogradely labeled phrenic motoneurons in contralateral phrenic nucleus (PN) 8 days after PhN ligation and bilateral injection of retrograde tracer Fluorogold (FG) into diaphragm **(A)**. Magnified area in box on right side shows intensive FG labeling and light sGCβ1-IR staining in motoneurons of contralateral PN. Intense expression of sGCβ1 and weak FG labeling is seen in phrenic motoneurons on ipsilateral side (left side). Number of retrogradely labeled sGCβ1-IR phrenic motoneurons/section in cervical segments (C3-C6) on contralateral and ipsilateral sides **(B)**. Level of sGCβ1 protein in ventral horn of C3-C6 segments after normalization to β-actin **(C)**. Data in **(B)** are mean values of five experiments ± SEM. Data in **(C)** are mean values of four experiments ± SEM. Results were statistically evaluated using Unpaired *t*-test; ^\*\*\*\*^*p*\< 0.0001 in **(B)**, ^∗^*p* \< 0.05 in **(C)**. Scale bars = 1000 μm in **(A)**, 100 μm (insets in **A**).](fphys-11-00700-g001){#F1}

Effect of PhN Ligation on sGCβ1 Changes in Phrenic Nerve {#S3.SS2}
--------------------------------------------------------

To study the effect of PhN ligation on the level of sGCβ1, the control, ipsilateral and contralateral samples of n. phrenicus were used for Western blot analysis ([Figure 2A](#F2){ref-type="fig"}). The level of sGCβ1 protein (normalized by the values of β-actin) was 0.726 ± 0.060 in the control phrenic nerve. Eight days of continuous PhN ligation caused a significant decrease in sGCβ1 level in the nerve below its ligation (0.398 ± 0.006) compared to the contralateral nerve (*p* \< 0.0001) or control nerve (*p* \< 0.0006). WB analysis showed that the level of sGCβ1 protein was only slightly modified comparing control and contralateral phrenic nerves (*p* \< 0.197). Immunolabeling with SMI-312 antibody (a specific marker for neurofilaments) was performed to visualize axons in the phrenic nerve. The neurofilament- and sGCβ1-IR axons were visible in the control group ([Figures 2B,C](#F2){ref-type="fig"}). At the end of survival, the sGCβ1- and SMI-312-positivity decreased in axons below the nerve ligation ([Figures 2D--F](#F2){ref-type="fig"}). These data indicate that PhN ligation decreased the level of sGCβ1 in motor axons, targeting the diaphragm innervation at the site of ligation.

![Soluble guanylyl cyclase β1 subunit immunoreactivity (sGC-IR), neurofilament immunoreactivity (NF-IR) and level of sGCβ1 protein in control phrenic nerve (PhN) in experimental animals 8 days after unilateral PhN ligation. Western blot analysis of sGCβ1 protein in PhN **(A)**. Longitudinal sections of PhN immunostained for sGC and neurofilaments (NF) in control **(B,C)**. sGC- and NF-IR in ipsilateral PhN below its ligation **(D,E)**. Merging of sGC- and NF-IR **(F)**. Level of sGCβ1 protein in **(A)** is normalized to β-actin. Data are mean values of four experiments ± SEM. Results were statistically evaluated using Unpaired *t*-test; \**p* \< 0.05. Scale bars = 50 μm.](fphys-11-00700-g002){#F2}

Degeneration of the Phrenic Nerve {#S3.SS3}
---------------------------------

We next examined anterograde degeneration of the phrenic nerve below the ligation point ([Figures 3A,B](#F3){ref-type="fig"}). Axons of the contralateral phrenic nerve showed the typical morphological characteristics of intact myelinated axons. In the ipsilateral nerve, evident axonal loss was detected, and the tissue contained numerous disintegrated myelin structures typical for Wallerian degeneration. Axonal counts indicated massive loss of nerve fibers in the ipsilateral nerve (167.1 ± 0.060 vs. 357.8 ± 0.060 axons/nerve, *p* \< 0.0001) 8 days after PhN ligation ([Figure 3B](#F3){ref-type="fig"}).

![Eight days continuous phrenic nerve (PhN) ligation caused axonal degeneration in the ipsilateral PhN and abnormalities at neuromuscular junctions (NMJs) in the ipsilateral hemidiaphragm. Noticeable axonal degeneration in PhN below ligation point (ipsilateral) **(A)**. Quantitative analysis revealed significantly lower number of intact axons in ipsilateral than in contralateral PhN **(B)**. Numbers of intact, partially denervated and fully denervated (arrowhead) NMJs in ipsilateral vs. contralateral hemidiaphragm **(C)**. Data in **(B,C)** are mean values of five experiments ± SD. Results were statistically evaluated using Unpaired *t*-test; \*\*\*\**p*\< 0.0001. Scale bar = 100 μm.](fphys-11-00700-g003){#F3}

Axon Terminals in the Ipsilateral Hemidiaphragm Showed Reduced sGCβ1- and SV2-Immunoreactivity {#S3.SS4}
----------------------------------------------------------------------------------------------

Confocal imaging was used to visualize the sGCβ1-IR terminal boutons of the phrenic motoneurons at NMJs in the diaphragm ([Figures 4A--C,E](#F4){ref-type="fig"}). As shown in [Figures 4C,G,H](#F4){ref-type="fig"}, sGCβ1-IR puncta were visible in the contralateral hemidiaphragm inside the NMJs and beyond. Significant reduction in the area commonly occupied by the synaptic vesicles and/or sGCβ1 boutons over the endplates ([Figures 4C,D,G--J](#F4){ref-type="fig"}) was seen in the ipsilateral diaphragm ([Figures 4E,F,K--N](#F4){ref-type="fig"}). Similarly, marked reduction in sGCβ1 immunolabeled puncta was seen through the ipsilateral NMJs ([Figures 4E,F,K--N](#F4){ref-type="fig"}). These results show that correlation exists between the loss of sGCβ1 protein and the loss of axons in the ipsilateral phrenic nerve, and the reduction in sGCβ1immunolabeled terminals in the ipsilateral NMJs.

![Representative confocal images showing morphological changes at diaphragm NMJs after unilateral PhN ligation. NMJs of diaphragm muscle were visualized *via* labeling post-synaptic acetylcholine receptors with rhodamine-conjugated α bungarotoxin (red), axons were labeled with anti-neurofilament antibody- SMI-312 (blue), soluble guanylyl cyclase β1 subunit immunoreactive terminals with sGCβ1 antibody (green) and presynaptic terminals with SV2 (green) **(A--N)**. All NMJs were intact in contralateral diaphragm (**A,C,D**, arrows), characterized by complete overlap of axons and presynaptic sGCβ1-IR terminals **(C)**, and axons and presynaptic vesicles **(D)** with postsynaptic acetylcholine receptors. Note dense accumulation of sGCβ1 stained puncta in contralateral diaphragm **(C,H)**. Ipsilateral hemidiaphragm shows NMJ abnormalities represented by partial (**B,E**, arrows) or complete denervation (**F,** arrows). Marked preterminal axon thinning is visible in ipsilateral diaphragm (**B,** asterisk). At NMJs in ipsilateral hemidiaphragm we detected almost complete reduction of sGCβ1- **(E,L)** and SV2-immunostained puncta **(F,N)**. Note single immunostaining from (**C,D**, arrows) in **(G,I)** -- BTX, in **(H)** -- sGC and in **(J)** -- SV2. BTX **(K,M)**, sGC **(L),** and SV2 **(N)** in ipsilateral hemidiaphragm (single immunostaining from **E,F**, arrows). BTX- rhodamine-conjugated α bungarotoxin, NF-neurofilaments. sGC-soluble guanylyl cyclase. Scale bar = 75 μm.](fphys-11-00700-g004){#F4}

Morphological Alterations at Diaphragm NMJs After PhN Ligation {#S3.SS5}
--------------------------------------------------------------

This was further supported by pathological alterations at the diaphragm NMJs. As shown in [Figure 3C](#F3){ref-type="fig"}, differences in numbers of NMJs (intact, partially or fully denervated) between the contra- and ipsilateral hemidiaphragm were significant. Eight days after PhN ligation, the contralateral hemidiaphragm exhibited fully intact NMJs (99% ± 0.006). The number of intact NMJs was three-fold lower in the ipsilateral than the contralateral hemidiaphragm. Continuous PhN ligation induced a high percentage (67%) of NMJ abnormalities in the ipsilateral hemidiaphragm. These alterations were identified as partial ([Figures 4B,E](#F4){ref-type="fig"}, arrows) or full denervation ([Figure 4F](#F4){ref-type="fig"}, arrow). The percentage of NMJs exhibiting partial denervation was significantly greater (44%, *p* \< 0.0001) compared to complete junction denervation (23%, *p* \< 0.0001). However, the percentage of intact juctions declined to 33% (*p* \< 0.0001) in the ipsilateral hemidiaphragm muscle ([Figure 3C](#F3){ref-type="fig"}).

cGMP Levels in the Diaphragm {#S3.SS6}
----------------------------

ELISA assay was used to detect cGMP levels in the control diaphragm and in experimental hemidiaphragms 8 days after unilateral PhN ligation. Under physiological conditions, the level of cGMP was 36 ± 0.15 pmol/ml. Significant decrease in cGMP levels was observed in the contra- and ipsilateral diaphragm (*p* \< 0.0001 vs. control) 8 days after PhN ligation. However, the decrease was nine -- times stronger in the ipsilateral than the contralateral hemidiaphragm (0.7 ± 0.02 vs. 6.5 ± 0.01; *p* \< 0.0001).

EMG Activity in the Diaphragm {#S3.SS7}
-----------------------------

Immediately after ligation the frequency of breathing decreased significantly by more than 40%, from ∼1.6 Hz in the control group to ∼0.9 Hz in the ligation group ([Figures 5A,B](#F5){ref-type="fig"}). While there was no difference in breathing frequency on both sides of the diaphragm, the overall muscle activity differed dramatically. PhN ligation caused rapid and pronounced reduction in EMG activity amplitude on the ipsilateral side, while the amplitude of EMG on the contralateral side slightly increased ([Figure 5B](#F5){ref-type="fig"}). Eight days after continuous PhN ligation, EMG on both sides of the diaphragm muscles changed markedly: on the ipsilateral side the residual activity disappeared, whereas on the contralateral side there was several-fold increase in EMG amplitude while breathing frequency remained unchanged ([Figure 5C](#F5){ref-type="fig"}). This increase in EMG activity was due to an enormous increase in muscle fibers activity, but not to prolongation of that activity, since the duration of the rhythmic contractions did not change significantly.

![Representative EMG recordings from contralateral and ipsilateral hemidiaphragms in control animals **(A)**, immediately after phrenic nerve (PhN) ligation **(B)**, and 8 days after PhN ligation **(C)**.](fphys-11-00700-g005){#F5}

Discussion {#S4}
==========

Phrenic nerve injury is a well-established model, since the nerve ligation or transection deactivates descending excitatory drive from the phrenic motoneurons to the diaphragm. Within the framework of "Network Physiology" ([@B7]; [@B5]; [@B6]; [@B24]) this study focuses on the network interactions between the central nervous system (cervical spinal cord) and the periphery (e.g., phrenic nerve and respiratory muscle). Experimental and clinical data have shown that PhN injury can have significant negative effects on diaphragm and lung, affecting their function. Recently published data revealed that the complete paralysis of the unilateral diaphragm could influence the loss of vital capacity and total lung capacity in the aged rats ([@B12]) but respiratory function parameters show compensation after 4 weeks in the young rats, and contralateral phrenic nerve transfer can enhance pulmonary function. In the present study we speculated, whether changes in the ipsilateral hemidiaphragm (e.g., loss of sGCβ1-IR nerve terminals, decreased level of cGMP, NMJs abnormalities and silencing of EMG activity) emerged from these interactions as an indicator of physiological state and function. To answer this fundamental question, immunohistochemical and histological staining, Western blotting, Elisa and electrophysiology were used. It is known, that the NO/sGC/cGMP signaling pathway plays a key role in regulation of a variety of vital function. However, its functional importance in the respiratory pathway (brainstem-phrenic motoneurons -- phrenic nerve -- diaphragm) still remains unclear. We have previously reported that the premotor bulbospinal respiratory pathway connecting the bulbar respiratory centers with the motoneurons of the phrenic nucleus in dog and rat is nitrergic ([@B31]; [@B10]; [@B18]). In the present study we provide the first evidence that unilateral ligation of the phrenic nerve decreases the level of sGCβ1 in ipsilateral motor axons targeting the diaphragm innervation, and suggest that this effect is linked with loss of sGCβ1immunolabeled terminals at the ipsilateral MNJs. Subsequent remarkable unilateral decrease in cGMP level in the diaphragm and silencing of EMG activity can be attributed to the fact that sGC is involved in transmitter release at the NMJs via the sGC-cGMP pathway.

The phrenic nerves are generally considered as motor nerves whose primary function is to supply motor innervation to the diaphragm ([@B28]). By means of a retrograde tracing technique we provide data showing direct projection from the hemidiaphragm to the contralateral FG-labeled phrenic motoneurons of the phrenic nucleus. Eight days of continuous PhN ligation caused significant increase in sGCβ1 (immunoreactivity and protein level) in the ipsilateral phrenic motor pool, but the nerve ligation reduced cGMP levels in the ipsilateral hemidiaphragm muscle. We suggest that continuous PhN ligation may lead to increased metabolic activity by ipsilateral phrenic α-motoneurons, normally expressing sGCβ1 at low level ([@B10]). Recently published data indicate that 2 weeks of tetrodoxin (TTX) phrenic nerve blockade significantly increased both the total phrenic motoneuron surface area and the dendritic surface area ([@B30]). These authors suggest that ipsilateral phrenic motoneuron morphological adaptations are consistent with normalization of motoneuron excitability following prolonged alterations in motoneuron activity. Phrenic motoneuron structural plasticity is probably more dependent on motoneuron activity (or descending input) than muscle fiber activity. Furthermore, motoneurons have high metabolic activity associated with neurotransmission ([@B2]). They are relatively large cells, and therefore their permanent stimulation may induce unusually high influx of Ca^2+^ into these neurons ([@B16]; [@B44]; [@B43]). [@B32] demonstrated that the primary excitatory drive for phrenic motoneurons is glutamatergic, probably acting via Ca^2+^-permeable, AMPA subtype glutamate receptors. As reported in previous studies ([@B39]; [@B31]) nitrergic axons projecting from the dorsal respiratory group and/or rostral ventral respiratory group to the phrenic motoneurons release a signaling NO molecule as a co-transmitter with glutamate. We suggest that descending synaptic input of NO, which activates sGCβ1 in phrenic motoneurons ([@B10]) and/or permanent blocking of NO-mediated cGMP synthesis in the hemidiaphragm after PhN ligation probably upregulates sGCβ1 in the ipsilateral phrenic motoneurons. The ability of motoneurons to handle calcium may also be important in determining their response to stimulation ([@B11]). Previous studies from our and other labs ([@B25]; [@B36]; [@B26]) have revealed a relatively low endogenous calcium-buffering capacity in motoneurons in comparison with other neuronal populations.

Phrenic motor neuron loss and subsequent diaphragm NMJ abnormalities are well described in various experimental models ([@B37]; [@B34], [@B33]; [@B14]), however, the correlation of these alterations with sGC transmission at the NMJs is lacking. [@B7] have shown that dramatic change in network structure with transition from one physiological state to another within a short time windows indicates a high flexibility in the interaction between physiological systems in response to change in physiological regulation. In the current study, we detected correlation between motor fiber damage, the reduction in sGCβ1-IR terminals and synaptic vesicles at the post synaptic acetylcholine receptors in the ipsilateral MNJs. Our in-depth diaphragm NMJ analysis demonstrated the phrenic axonal degeneration process, as many motor endplates were found to be fully (23%) or partially (44%) denervated as early as 8 days after PhN ligation. The pattern of diaphragm denervation has been studied in the central region of the diaphragm, particularly innervated by phrenic motor neurons located at the C5 level of the spinal cord ([@B23]). In addition, we found significant loss of motor axons (∼47%) and changes in phrenic nerve morphology such as signs of Wallerian degeneration which affected morphological NMJ denervation and might have resulted in reduced spontaneous diaphragm EMG. We suggest that significant decrease in sGCβ1 level in the efferent arm of the respiratory pathway (diaphragmatic nerve) just below the nerve ligation site is also consistent with the cGMP level alterations in the diaphragm. Although we found marked reduction in cGMP levels in both the contra- and ipsilateral hemidiaphragm, the decrease was significantly greater (nine times) on the side of the ligation. Currently, the functional effects of unilateral 8-day continuous PhN ligation are unknown but significant effects can be inferred from recently published experimental study ([@B12]). Our results demonstrate that loss of motor axons and/or sGCβ1 in motor nerve terminal boutons after phrenic nerve ligation is consistent with alterations in synaptic transmission, the release of cGMP and endplate morphology in the ipsilateral hemidiaphragm.

As shown previously, most NO effects are mediated by sGC activation and subsequent cGMP formation ([@B38]; [@B21]). cGMP influences the activity of at least three different targets, i.e., cGMP transduces the NO signal to the cGMP-regulated protein kinases, cGMP-activated phosphodiesterases (PDEs) and cGMP-gated ion channels ([@B41]; [@B19]). As mentioned earlier, nNOS, cGMP and cGMP-dependent kinases are all concentrated at or near the NMJs ([@B20]; [@B22]; [@B35]). [@B38] indicated that the rat diaphragm produces NO or NO-derivates which downregulate the contractile function via their effect on excitation-contraction coupling. [@B1] demonstrated that cGMP, which is present in most muscle fibers of the rat diaphragm, is primarily distributed in the subsarcolemmal region of individual fibers. It is important to emphasize that NO/sGC/cGMP signaling depends not only on cGMP formation, but it is also critically influenced by the activity of enzymes responsible for cGMP degradation ([@B21]). In the diaphragm, PDEs degrade cGMP and terminate the cGMP-signal ([@B35]). Our results support sGC-dependent cGMP signaling in the lower respiratory pathway, but the precise target of cGMP action in the diaphragm muscle still remains unclear.

[@B5] investigated the network of interactions between the brain, cardiac and respiratory system, and demonstrated that a network approach to physiological interactions is necessary to understand how modulations in the regulatory mechanism of individual systems translate into reorganization of physiological interactions across the organism. It is known that the function of inspiratory muscles is very important in the breathing process, because they are required to contract repetitively without interruption throughout life. Our results show remarkable unilateral silencing of EMG activity immediately after PhN ligation, and the failure of recovery of EMG activity in the ipsilateral diaphragm muscle 8 days after pathological continuous nerve ligation. These changes most likely depend on both the loss of ipsilateral phrenic axons (∼47%) and the degree of diaphragm denervation. We also show spontaneous and increased EMG activity in the contralateral hemidiaphragm. Such enhancement, recorded 8 days after unilateral PhN ligation, seems to be the result of increased muscle fiber activity. [@B29] observed the absence of EMG activity in the ipsilateral diaphragm immediately and 3 days after C2 hemisection using diaphragm electrodes. They recorded progressive increase in the proportion of spontaneous recovery of EMG activity ispilaterally to C2 hemisection over time. Only 10% of unanesthetized animals displayed recovery at 7 days after hemisection. These findings are inconsistent with previously reported spontaneous ipsilateral increase in EMG activity. [@B45] reported that the ipsilateral phrenic nerve deactivated by lateral C2 SCI was spontaneously reactivated 7 days post-SCI. These authors also indicated that ipsilateral phrenic nerve reactivation was greater at 3 months compared with 7 days post-SCI, and that it was enhanced after contralateral phrenicotomy. There was gradual recovery of rhythmic diaphragm muscle activity ipsilaterally to cervical spinal cord injury over time, consistent with neuroplasticity and strengthening of spared, contralateral descending premotor input to the phrenic motoneurons ([@B27]). This suggests that spontaneous recovery of the respiratory pathway may depend on post-lesional time. After 8 days of hemidiaphragm paralysis induced by continuous PhN ligation, there was marked increase in sGCβ1 (immunoreactivity and protein level) in the phrenic motor pool, but the protein sGC level in the phrenic nerve below its ligation point and the cGMP level in the ipsilateral diaphragm were markedly decreased. These findings are consistent with alterations such as dramatic axonal loss below the nerve ligation site, reduction in sGCβ1-IR terminals/synaptic vesicles at the ipsilateral NMJs, significant morphological NMJ degeneration, and subsequent remarkable unilateral silencing of EMG activity.

Taken together, our results indicate the role of sGC/cGMP signaling in the lower respiratory pathway (phrenic motoneurons -- phrenic nerve -- diaphragm), although the complexity of respiratory circuits (e.g., relationship and coupling between the central nervous system and the periphery) requires further *in vivo* studies to clarify the functional contribution of this signaling pathway in respiratory neuroplasticity following cervical spinal cord and/or PhN injuries.

Data Availability Statement {#S5}
===========================

The datasets generated for this study are available on request to the corresponding author.

Ethics Statement {#S6}
================

This study was carried out in accordance with the principles of the Directive 2010/63/EU of the European Parliament and of the Council. The protocol was approved by the State Veterinary and Food Administration of the Slovak Republic (4434/16-221/3 and 715/19-221/3). All rights reserved under the supervision of the Ethical Council of the Institute of Neurobiology BMC SAS.

Author Contributions {#S7}
====================

NL, LH, JG, and MM contributed to the study concept and design of the experiments. LH and SM performed the PhN ligation. LH and AK performed the Western Blotting and cGMP measurement using ELISA, Fluorogold was injected into the right diaphragm by JP. IV and ZD performed the PhN staining and quantitative analysis. LH, KB, MB, AS, and VL performed the immunohistochemical staining of spinal cord, diaphragm and phrenic nerve, and data analyses. JG and ŠP recorded the diaphragm EMG. LH and NL wrote the manuscript. All authors contributed to the manuscript revision, read and approved the submitted version.

Conflict of Interest {#conf1}
====================

SM was employed by the company Diagnostic Center of Pathology in Prešov, Alpha medical, s.r.o. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Funding.** This work was supported by the VEGA Grants 2/0168/17 (NL), 2/0098/20 (JG), and 2/0049/19 (JP) from the SAS, APVV-15-0766 (NL) and the project: Formation and development of a diagnostic procedure in the treatment of trauma -- injured spinal cord (ITMS 26220220202), supported by the Research & Development Operational Programme funded by the ERDF (NL).

We are grateful to Mgr. Z. Barkacs for her technical assistance.

<https://home.saske.sk/~tomori/software/NeuroCounter>

Supplementary Material {#S10}
======================

The Supplementary Material for this article can be found online at: <https://www.frontiersin.org/articles/10.3389/fphys.2020.00700/full#supplementary-material>

###### 

Click here for additional data file.

cGMP

:   cyclic guanosine monophosphate

EMG

:   electromyography

FG

:   fluorogold

GC

:   guanylyl cyclase

IR

:   immunoreactivity

NADPH

:   nicotinamide adenine dinucleotide phosphate

NMJ

:   neuromuscular junction

nNOS

:   neuronal nitric oxide synthase

NO

:   nitric oxide

PhMNs

:   phrenic motoneurons

PhN

:   phrenic nerve

sGC β 1

:   β 1 subunit of soluble guanylyl cyclase

sGC

:   soluble guanylyl cyclase.

[^1]: Edited by: Plamen Ch. Ivanov, Boston University, United States

[^2]: Reviewed by: Sergi Garcia-Retortillo, University of Girona, Spain; Andras Eke, Semmelweis University, Hungary

[^3]: This article was submitted to Respiratory Physiology, a section of the journal Frontiers in Physiology