Monthly Archives: January 2019

Complete Occlusion of Aortic Arch Branch Vessels Diagnosed by Intraoperative Point-of-Care Ultrasonography: A Case Report

DOI: 10.31038/JCRM.2019211

Abstract

The use of point-of-care ultrasonography (POCUS) is becoming increasingly widespread and clinically important. Here, we report a case in which the use of POCUS in the operating room identified a previously undiagnosed bilateral subclavian artery stenosis. A 70-year-old woman with a diagnosis of meningioma was scheduled to undergo a craniotomy. After the induction of general anesthesia, an abnormal radial artery waveform was observed bilaterally. A blood flow deficiency in both upper extremities was suspected, and POCUS was immediately performed. Using POCUS, the anesthesiologist identified stenosis and hypoplasia of the bilateral subclavian arteries, complete occlusion of the bilateral common carotid arteries, and remarkably dilated bilateral vertebral arteries. By identifying these vascular abnormalities with POCUS, the anesthesiologist was able to adjust the patient’s perioperative circulatory management and lead to a positive postoperative outcome.

Background

Point-of-care ultrasonography (POCUS) is becoming progressively more widespread and clinically valuable [1]. The use of POCUS by anesthesiologists and critical care physicians may help improve outcomes in surgical and critically ill patients [2,3]. During the perioperative period, POCUS is a very efficient and beneficial tool that can be used to diagnose previously unrecognized clinical conditions. However, to the best of our knowledge, there are very few reports in which previously undiagnosed vascular abnormalities of aortic arch branch vessels were identified by POCUS performed by anesthesiologists in the operating room. Here, we report a case in which we suspected bilateral subclavian artery stenosis based on abnormal bilateral radial artery pressure waveforms after induction of anesthesia, immediately performed POCUS, and identified a previously undiagnosed bilateral subclavian artery stenosis.

Case presentation

A 70-year-old woman diagnosed with a meningioma was scheduled to undergo a craniotomy. In addition to the meningioma, preoperative head magnetic resonance imaging (MRI) revealed that bilateral common carotid arteries were abnormally narrow and bilateral vertebral arteries were abnormally dilated; the blood vessels below the neck were not evaluated. She was alert and oriented with no neurological abnormalities, and the results of the Allen’s test were negative.

Her noninvasive blood pressure, measured from the left upper arm at admission, was 139/111 mm Hg. After induction of general anesthesia, a 20-gauge catheter was inserted in her left radial artery, and the arterial blood pressure was 79/54 mm Hg. The waveform ascended slowly, and the pulse pressure was unusually narrow. Phenylephrine (0.1 mg) was given intravenously 3 times, and her systolic pressure increased by 20 mm Hg after each dose; however, the abnormal wave form did not improve. The arterial line circuit was intact and exhibited no bending, kinking, or air bubbles. We conducted a flash test, which was normal. Because we initially suspected stenosis of the left radial artery, we inserted a second arterial catheter in her right radial artery; however, the pressure tracing revealed a similarly abnormal waveform (Fig. 1). The patient’s bilateral ulnar pulses were palpated; however, her radial pulses were stronger and more easily felt. We suspected that the patient may have had a perfusion deficiency in the bilateral upper extremities; consequently, POCUS was immediately performed.

JCRM 2019-101 - Masafumi Idei Japan_Fig. 1

Figure 1. Image shows similarly abnormal waveform of arterial blood pressure from bilateral radial arteries.

Using POCUS, we identified hypoplasia and stenosis of bilateral subclavian arteries and complete obstruction of bilateral common carotid arteries (Fig. 2). At around 3 mm in diameter, the bilateral subclavian arteries were around a quarter of the size of the accompanying subclavian veins. The diagnosis of bilateral common carotid artery obstruction was easily confirmed by color doppler imaging. Furthermore, bilateral vertebral arteries were significantly dilated with a diameter of approximately 7 mm. Because both subclavian arteries were observed to be abnormally narrow via imaging, we believed that the blood pressure measurement in the bilateral upper extremities was not reliable. Therefore, we inserted a third arterial catheter in the patient’s left dorsalis pedis artery. When the arterial blood pressure measured from her left radial artery was 81/56 (mean 66) mm Hg, it was 66/47 (mean 55) mm Hg from her right radial artery and 160/69 (mean 96) mm Hg from her left dorsalis pedis artery (Fig. 3). Based on the POCUS findings, we determined that the blood pressure measured from her dorsalis pedis artery was most accurate for perioperative circulatory management. The patient emerged from anesthesia immediately after the operation was concluded, and no neurological deficits were observed.

JCRM 2019-101 - Masafumi Idei Japan_Fig. 2

Figure 2. Ultrasonography of the neck. Hypoplasia and stenosis of bilateral subclavian arteries, complete obstruction of bilateral common carotid arteries, and dilated bilateral vertebral arteries were observed.

JCRM 2019-101 - Masafumi Idei Japan_Fig. 3

Figure 3. Image shows arterial blood pressure measured from the left radial artery and left dorsalis pedis artery.

Due to the POCUS findings, we conducted a detailed examination by enhanced computed tomography (CT) postoperatively. All 3 branch vessels of her aortic arch (brachiocephalic trunk, left common carotid artery, and left subclavian artery) were completely occluded at their regions of origin, and well-developed collateral pathways supplied blood to her arms and brain (Fig. 4). A history of aortic arch syndrome (Takayasu’s arteritis) was suspected.

JCRM 2019-101 - Masafumi Idei Japan_Fig. 4

Figure 4. Computed tomography (CT) and schema. All three branch vessels of the aortic arch were totally occluded and well-developed collateral pathways supplied blood flow to her arms and brain.

Discussion

POCUS is being utilized more often in clinical settings [1–3]. In this case, bilateral subclavian artery stenosis was suspected based on abnormal pressure waveforms of bilateral radial arteries. POCUS was immediately performed, which revealed a previously undiagnosed vascular abnormality that could affect the patient’s brain and upper extremities. Postoperative enhanced CT confirmed the findings identified by POCUS.

In POCUS findings, we could easily identify vascular abnormalities. The vessel walls of bilateral common carotid arteries were thickened circumferentially, and color doppler imaging showed complete occlusion. The bilateral subclavian arteries were only about 3 mm in diameter, and were clearly smaller than accompanying subclavian veins. Bilateral vertebral arteries, which are normally about 3 mm in diameter [4,5], were significantly dilated to approximately 7 mm in diameter.

In order to accurately interpret the POCUS findings, sufficient experience is required [1]. The common carotid and subclavian arteries are routinely identified by anesthesiologists during ultrasound-guided central venous catheter insertion, and the vertebral artery is frequently identified during ultrasound-guided brachial plexus blocks. Moreover, anesthesiologists often evaluate stenosis of the common carotid artery as part of perioperative POCUS. Because we routinely use ultrasonography and are accustomed to POCUS, we immediately noticed the abnormality of the arteries, which contributed to our ability to make this diagnosis.

Based on the POCUS findings, we relied on the lower extremity blood pressure for perioperative circulatory management. Because POCUS and preoperative head MRI results showed obstruction of the bilateral common carotid arteries and dilation of the bilateral vertebral arteries, the patient’s cerebral blood flow was thought to depend on the vertebral arteries. We judged that only the upper extremities were hypotensive, and that the blood pressure measured from the lower extremities more accurately reflected the perfusion pressure of the patient’s major organs. We investigated the regional oxygen saturation (rSO2) as an index of cerebral blood flow. However, it would have been difficult to monitor the patient’s rSO2 in this case due to the location of the surgical field.

This case highlights how circulatory management that blindly gives priority to radial artery pressure could lead to end-organ hypertension, which could increase bleeding from the surgical field. In this case, the radial artery blood pressure was significantly lower than the blood pressure of the lower extremities (about 30–40 mm Hg less in mean arterial pressure). We thought that tight blood pressure control was warranted, because excessive hypotension could have caused ischemia of the upper extremities, and hypertension could have induced surgical bleeding. Thus, we sought to maintain a mean arterial pressure in the lower extremities of approximately 80 mm Hg (130 mm Hg in systolic blood pressure), because at this pressure, the mean arterial pressure of the right radial artery stayed around 50 mm Hg. Although setting an appropriate target blood pressure was extremely difficult in the limited environment and information during surgery, we deemed this target blood pressure as appropriate because the patient had no definite history of hypertension and it would enable us to avoid excessive hypertension, which is important during craniotomies.

The patient was similarly managed in the intensive care unit postoperatively, and she experienced no bleeding, neurological problems, or other complications including major organ failure. A CT scan performed 6 days after the operation showed complete occlusion of her brachiocephalic trunk, left common carotid artery, and left subclavian artery at their regions of origin from the aortic arch. Numerous collateral pathways from arteries of the superior mediastinum and chest wall, which branched from the aorta, supplied blood to the aortic arch branch vessels. Collateral pathways also supplied blood to bilateral vertebral arteries.

Although there is no conclusive evidence that the patient’s lower extremity blood pressure reflected her cerebral blood pressure and that our target arterial pressure was exactly appropriate, our careful circulatory management that avoided excessive hypertension and hypotension during the perioperative period, might have contributed to the patient’s good postoperative neurological outcome. Because preoperative head MRI data revealed abnormalities in the bilateral common carotid arteries and bilateral vertebral arteries, further detailed evaluation of other blood vessels below the neck by enhanced CT should have been considered before operation and may be important for similar situations in the future.

After consulting with a cardiothoracic surgeon, the patient subsequently confirmed that she had a weak pulse during a medical exam when she was in her twenties. Since this issue had been identified in the patient’s youth and because she had few risk factors of arteriosclerosis such as hypertension, diabetes mellitus, and smoking, a history of aortic arch syndrome (Takayasu’s arteritis) was suspected. There were, however, few findings and symptoms suggesting active arteritis such as an elevation of inflammatory markers, and the patient denied further evaluation and treatment; therefore, it is unknown whether the complete occlusion of the aortic arch branch vessels was caused by aortic arch syndrome.

Some reports have described anesthetic management of patients diagnosed with stenosis and/or obstruction of the aortic arch branch vessels that occurs with aortic arch syndrome (Takayasu’s arteritis) [6,7]. These reports highlight the usefulness of ultrasonography to determine the best site for blood pressure monitoring [6] and the efficacy of rSO2 monitoring as an index of cerebral blood flow [7]. However, to the best of our knowledge, there are few reports in which a previously undiagnosed stenosis and/or obstruction of the aortic arch branch vessels was identified by POCUS during a perioperative period. Because POCUS has the potential to significantly improve perioperative care, anesthesiologists should strive to acquire this skill. We expect future clinical studies and cases that involve POCUS during surgery will further demonstrate how its use can improve perioperative outcomes.

Acknowledgements: We would like to thank our patient for providing to publish this case report. We would also like to thank Editage (www.editage.jp) for English language editing.

Funding: The authors declare that they have no funding.

Authors’ contributions: MI and RY wrote the draft of the manuscript. TY made the figures. IW reviewed the manuscript. The final version of the manuscript was approved by all authors.

Ethics approval and consent to participate: Not applicable.

Consent for publication: Written informed consent was obtained from the patients for publication of this case report and accompanying images.

Competing interests: The authors declare that they have no competing interests.

Author details: Department of Anesthesia, Yokohama Minami Kyosai Hospital, 1-21-1 Mutsuurahigashi, Kanazawa-ku, Yokohama 2360037, Japan

Abbreviations

POCUS: Point-of-care ultrasound

MRI: magnetic resonance imaging

ICU: intensive care unit

CT: Computed tomography

rSO2: regional oxygen saturation

References

  1. Moore CL, Copel JA. (2011) Point-of-Care Ultrasonography. N Engl J Med. 364: 749–757. [Crossref]
  2. Johnson DW, Oren-Grinberg A. (2011) Perioperative point-of-care ultrasonography: the past and the future are in anesthesiologists’ hands. Anesthesiology. 115(3): 460–2. [Crossref]
  3. Holm JH, Frederiksen CA, Juhl-Olsen P, Sloth E. (2012) Perioperative Use of Focus Assessed Transthoracic Echocardiography (FATE). Anesth Analg. 115(5): 1029–32. [Crossref]
  4. Park JH, Kim JM, Roh JK. (2007) Hypoplastic vertebral artery: frequency and associations with ischaemic stroke territory. J Neurol Neurosurg Psychiatry. 78(9): 954–8. [Crossref]
  5. Ogeng’o J, Olabu B, Sinkeet R, Ogeng’o NM, Elbusaid H. (2014) Vertebral Artery Hypoplasia in a Black Kenyan Population. Int Sch Res Notices. 2014: 934510.
  6. Narasimha PK, Chaudhuri S, Joseph TT. (2013) Utility of intra-operative ultrasound in choosing the appropriate site for blood pressure monitoring in Takayasu’s arteritis. Indian J Anaesth. 57(1): 66–8. [Crossref]
  7. Xiao W, Wang T, Fu W, Wang F, Zhao L. (2016) Regional cerebral oxygen saturation guided cerebral protection in a parturient with Takayasu’s arteritis undergoing cesarean section: a case report. J Clin Anesth. 33: 168–72. [Crossref]

Validation of an Automated Extraction Procedure for Amino Acids and Acylcarnitines for Use with Tandem Mass Spectrometry for Newborn Screening

DOI: 10.31038/EDMJ.2019314

Abstract

A certified reagent kit for newborn screening was transferred on a fully automated dried blood spot platform. The dried blood spot cards are directly eluted and the extract is online guided to tandem mass spectrometry instrument, where the amino acid and acyl carnitine panel is detected. The method takes 2 minutes per sample and requires no human interaction for up to 500 samples. The method is fully standardized through the automation and the usage of only certified consumables and reference material. The manual reagent kit was first modified to fit the automated platform, secondly validated and third, successfully transferred into a routine newborn screening laboratory.

Keywords

Dried blood spot, Newborn, Screening, Amino acid, Carnitine, Automation

Introduction

Newborn screening (NBS) is a public health program provided by most of the countries around the world aimed at screening newborns for a list of serious genetic and metabolic disorders. Early diagnosis of these conditions can help prevent their further development, which if untreated often results in brain damage, organ damage, and even death [1–4]. A routine neonatal screening procedure requires that a health professional takes a few drops of blood from the baby’s heel, applies them onto special filter paper and sends such prepared samples to a laboratory for a number of analytical tests [5].

The amino acids and acyl carnitines are detected in modern methods with tandem mass spectrometry (MS/MS). MS/MS is a fundamentally different technology than systems previously used by most newborn screening laboratories, such as bacterial inhibition assays. It is a versatile and modular system that can be easily adapted to the preferred testing approach by the user. This has led to numerous variations of newborn screening by MS/MS, and it became a challenge to compare results between laboratories. There is a recognized need to develop consensus solutions to provide more consistency between MS/MS screening programs [6–9].

This validation addresses a certified screening method using fully automated analysis equipment from sample recognition towards extraction and analysis. The method was based on the commercial MassChrom® Reagent Kit from Chromsystems to allow full standardization of the complete process. The method development was performed at the Shimadzu laboratory (Reinach, Switzerland) and the application has been transferred to a routine NBS laboratory in Switzerland (Childrens Hospital Zurich). After the method development, a validation was performed, followed by a transfer into a routine laboratory focusing on inter-, intra-day variations, correlation and robustness. To standardize the sample preparation, handling and extraction, a fully automated DBS system from CAMAG (DBS- MS 500) was integrated into the workflow. The main goal of this study was to validate this change in analytical procedure.

Application of automated DBS card handling systems, which are connected to mass spectrometry analyzers, offers a modern and fast approach where a circular area of the DBS is directly eluted from the filter paper card without any punching [9–12].

Materials and Methods

Chemicals

The following MassChrom Reagents (Chromsystems Munich, Germany) were used; Mobile Phase (No. 57001), Internal Standard (No. 57004), Extraction Buffer (No. 57008) and Mass Check Controls propionylcarnitine, decanoylcarnitine, lauroylcarnitine, myristoylcarnitine, palmitoylcarnitine and stearoylcarnitine, were purchased from Sigma-Aldrich (St. Louis, USA). MassChrom rinsing solution (No. 55007) was used for the rinsing process and also purchased from Chromsystems (Munich, Germany). Dried blood spot cards (903, TFN, MN818 and 2992) were provided by CAMAG (Muttenz, Switzerland). Fresh whole blood was obtained from the local blood donation center (Basel, Switzerland). The blood was previously tested for infectious diseases.

Analytical Methods

A DBS-MS 500 unit (CAMAG, Muttenz, Switzerland) was attached as front end to a modularHPLC system from Shimadzu (Kyoto, Japan), containing a system controller (CBM-20A), a Nexera X2pump and a degasser (DGU-20ASR). The loop outlet of the DBS-MS 500 system was connected with a1.8 m PEEK tubing (yellow 1/16» OD x .007» ID) using a KrudKatcher Ultra (KrudKatcher Ultra,Phenomenex, Torrance, CA, USA) inline filter at the mass spectrometer inlet. Analysis was performedin positive multiple reaction monitoring (MRM) mode on an electrospray ionization tandem massspectrometry system 8060 in Reinach and 8050 in Zurich (Shimadzu, Kyoto, Japan). The extract wasdirectly injected into the mass spectrometer without an analytical column.

The elution was performed isocratic with the mobile phase from the kit using a flow gradient starting at 0.2 mL/min to 0.6 mL/min in 0.3 min to 1.9 min, 1.91 min back to 0.2 mL/min, 2.0 min controller stop. The following m/z transitions were programmed for the mass spectrometry detection (Table 1);

Table 1. m/z transitions

Name

Precursor  m/z

Product  m/z

Dwell  (msec)

Q1 Pre Bias (V)

CE

Q3 Pre Bias (V)

Alanine

90.2

44.2

10

-10

-12

-19

Alanine-2H 4

94.2

48.2

10

-10

-12

-19

Arginine

175.2

70.2

10

-10

-24

-16

Arginine-2H 7

182.2

77.2

10

-10

-24

-16

Aspartic acid

134.2

134.2

10

-15

-12

-13

Aspartic acid-2H 3

137.2

75

10

-16

-16

-27

Citrulline

176.1

113.1

10

-10

-16

-25

Citrulline-2H 2

178.1

115.1

10

-10

-16

-25

Glutamic acid

148.15

84.1

10

-17

-17

-14

Glutamic acid-2H 5

153.1

88.2

10

-17

-18

-19

Glycine

76

30

10

-11

-12

-30

Glycine-13C3/ 15N1

79

32

10

-10

-16

-28

Leucine

132

86.2

10

-16

-12

-19

Leucine-2H 3

135

89.2

10

-16

-12

-19

Methionine

150.1

104.1

10

-18

-14

-22

Methionine-2H 3

153.1

107

10

-18

-13

-18

Ornithine

133.2

133.2

10

-16

-12

-27

Omithine-2H 6

139.2

76

10

-15

-19

-15

Phenylalanine

166.2

120.2

10

-18

-14

-28

Phenylalanine-2H 5

171.2

125.2

10

-18

-14

-28

Proline

116.2

70.1

10

-14

-18

-23

Proline-2H 7

123.2

77.1

10

-13

-18

-16

Tyrosine

182.1

123.1

10

-10

-18

-24

Tyrosine-2H 4

186.1

127.1

10

-10

-18

-24

Valine

118.2

72.1

10

-14

-13

-17

Valine-2H 8

126.2

80.2

10

-14

-13

-17

Carnitine

162

85

10

-21

-23

-20

Carnitine-2H 9

171

85.1

10

-10

-23

-17

Acetylcarnitine

204

84.9

10

-20

-25

-20

Acetylcarnitine-2H 3

207

85.1

10

-11

-19

-18

Propionylcarnitine

217.9

85

10

-20

-25

-20

Propionylcarnitine-2H 3

221

85.1

10

-12

-23

-18

Butyrylcarnitine

231.9

85

10

-20

-25

-20

Butyrylcarnitine-2H 3

235

85.1

10

-13

-23

-17

Valerylcarnitine

246

85.1

10

-20

-14

-21

Valerylcarnitine-2H 9

255

85.1

10

-14

-24

-17

C5DC-carnitine

276.2

85

10

-46

C5DC-carnitine-2H 6

282.6

85

10

-15

-26

-19

Hexanoylcarnitine

260

85

10

-20

-25

-20

Hexanoylcarnitine-2H 3

263

85.1

10

-13

-22

-14

Octanoylcarnitine

288

85

10

-20

-30

-20

Octanoylcarnitine-2H 3

291

85.1

10

-15

-23

-14

Decanoylcarnitine

316

85.1

10

-10

-24

-17

Decanoylcarnitine-2H 3

319

85.1

10

-10

-24

-17

Lauroylcarnitine

344

85

10

-20

-28

-20

Lauroylcarnitine-2H 3

347

85.1

10

-11

-25

-19

Myristoylcarnitine

372.2

85

10

-19

-30

-20

Myristoylcarnitine-2H 3

375

85.1

10

-11

-27

-17

Palmitovlcarnitine

399.9

85.1

10

-12

-28

-18

Palmitoylcarnitine-2H 3

402.9

85.1

10

-12

-28

-18

Stearoylcarnitine

427.9

85.1

10

-13

-29

-17

Stearoylcarnitine-2H 3

430.9

85.1

10

-13

-29

-17

ESI mode, Nebulizing gas; 1.5L/min, Heating gas Flow; 10 L/min, Interface Temperature; 300 °C,   DL Temperature 300°C, Heat Block Temperature 400°C, Drying gas Flow; 10 L/min

The data analysis uses a linear curve type using the internal standard for the area calculation. The two Chromsystems MassCheck levels are used according to the specified concentrations. All peaks are  as default integrated from 0.05 to 1.5 min with a width of 25 sec.

DBS-MS 500 Instrumentation and Settings

The MassChrom internal standard 57004 was dissolved in 25 mL extraction buffer 57008 according to the Chromsystems procedure and connected to elution bottle 1, this solution was used for extracting the DBS. The MassChrom rinsing solution was connected to rinsing bottle 1 (R1). The extraction head was cleaned in an ultra sound bath at 40 °C for 10 min prior to a large set of analyses. The extraction solvent was primed for 5 cycles and the rinsing solvents were flushed for 1 min (this process is an automated system prime method). The DBS cards were photographed with the built-in camera of the DBS-MS 500 before and after each extraction to check for the presence of a blood spot and to adjust the extraction head to the center of each spot. The Chronos for CAMAG software automatically recognized inadequate dried blood spots based on their roundness, diameter, and area. Inadequate DBS were excluded from analysis. The samples were extracted with a 4 mm diameter clamp and a volume of 60 μL and a 200 μL/min flow rate into a 20 μl loop (the 40 μl upfront volume is directed to the waste). The extraction solvent passes the sealed area on the DBS card horizontally from the bottom back to the bottom into a sample loop, which is online guided to the mass spectrometer after the elution step. The area next to the 4 mm extraction ring is not affected by the solvent and could be reused if needed. To complete the automated DBS extraction cycle, the system was rinsed for 20 s with R1 [11].

Standard and DBS Sample Preparation

DBS calibration samples for method development and implementation were prepared manually in the laboratory (described as calibrators). Later on for method validation, external quality control material from Chromsystems (Masscheck) was integrated into the analysis workflow and used for reference (described as controls).

Five calibration points were prepared for the validation process. Here, 10 mg of each analyte was dissolved in 10 ml of water for the amino acids and in 10 ml methanol for the acyl carnitines. For each level, a mix was prepared according to reference values provided by the MassCheck Controls from Chromsystems. EDTA stabilized blood was pooled and centrifuged at 1300 rcf for 5 minutes. The plasma and buffy coat layer was removed and replaced with saline (0.9% NaCl) [13]. After gently mixing, this washing procedure was repeated twice. After removing the wash solution, a spiked saline with amino acids and acyl carnitines in five different levels (A-E, Supplementary data) was added and gently mixed [14].

Aliquots of the five calibrator blood levels were spotted in 15, 30 and 50 μl (standard 30 μl) droplets by an Eppendorf pipette onto TFN, 226, 903, 818 and 2992 filter paper cards from different vendors and dried in a horizontal position for a minimum of three hours. After drying, the calibrators were placed in a plastic bag with desiccant and stored at -20 °C. Calibrators B and E were used to derive data for intra- and interday imprecision. Those levels were used to compare to the high and low control levels from the MassChrom kit which will be used for the routine afterwards.

To measure the external MassCheck control card, a DBS frame in the standardized format of 84.67 mm × 53.2 mm (w × h) was fastened to the previously cut out reference DBS spots (Figure 1). By using the X-offset of the DBS-MS 500, several extractions per reference spot can be performed.

EDMJ 2019-102 - Gaugler Stefan Switzerland_F1

Figure 1. Manually prepared QC card.

All DBS were prepared from the same blood source, prepared with the same procedure to neutralize potential hematocrit effects [11], [15], [16] , which were not investigated in this study. All samples were prepared at one site and shipped to the other laboratory under controlled environment, to also eliminate the inaccuracy of different sample preparation.

The five calibration points were used for method validation and transfer, afterwards the routine measurement will be referred to the two Chromsystems MassCheck control levels high and low.

Validation Procedure

The inter- and intra-assay precision was obtained by measuring calibrator levels B and E at 3 consecutive days six fold. Precision was evaluated within a single validation run (intra-day) as well as between three runs recorded on different days (inter-day). The precision was calculated as the percentage relative standard deviation (CV, %) within an analytical run (intra-day precision, n = 6) and over all three runs (inter-day precision, n = 18). Depending on the analyte, CVs for within-run precision can range from 15% to 25% and for inter-run precision can range from 20% to 35% for newborn screening MS/MS assays [14]. The accuracy was assessed from the overall mean of each MassCheck control concentration divided by its nominal value (bias, %). Since all signals are related to official and certified reference material from Chromsystems [17], no other accuracy assays were conducted.

The extraction recovery of the DBS-MS 500 autosampler was investigated for DBS samples at level D. The five compounds methionine, proline, valine, carnitine and lauroyl-carnitine were chosen for this experiment. Each spot was extracted five times in triplicates. Between two extractions, a drying time of approximately 15 min was programmed. Using the built-in camera of the DBS-MS 500, the extraction head automatically locked onto the same area in the center of the blood spot. The recovery was finally estimated as the percentage ratio of the analyte peak area of the first extraction to the sum of the peak areas of all subsequently conducted extractions [11, 14].

Correlation and robustness was performed on each three days in two different laboratories. The correlation is determined using the 5 calibrator levels on 903 type filter paper, measured four fold. Each day, the run included a double measurement of each MassCheck level high and low as reference.

The carryover was assessed to study the possible effects from specimens with a high concentration of an analyte on the result of the subsequent specimen or specimens. Although there are several potential sources for carryover within the MS/MS system, the most common source is the autosampler injection port and the tubing leading to the MS/MS electrospray unit. To determine the carry over in the automated DBS-MS 500 system, an extraction with a high concentration at level E is performed, followed by three extractions of blank DBS cards. For no carryover, the analyte concentrations of the blank sample should be below the LOD concentration previously determined. This should be repeated a minimum of five times [14].

Although no interferences are known or have been documented [17], the blank filter papers used for specimen collection were checked for possible interferences of common m/z for any of the analytes measured, especially since new filter paper sources and new lot numbers are used. Here, we investigated the most commonly used filter paper types from different vendors used for manufacturing the DBS cards. Paper types are TFN and 226 from Ahlstrom, 903 from EBF, MN 818 from Macherey Nagel and 2992 from Hahnemühle. Calibrator level B, D and F were measured as triplicates on each paper type to detect precision and any significant off-set.

The robustness is determined during the interference and repeatability study by varying the following parameters; technician, card type, eluent preparation (using solvent and buffer from different lots), elapsed times (with or without breaks between the measurements) and laboratory environment (location, temperature and humidity). All cards used for the interference study were spotted with level a calibrator.

Matrix effects, an extended interference study and stability tests were not performed, since this data is given by the MassChrom kit and its process remain unchanged [17].

Documentation

The documentation process is given by the DBS-MS 500 instrument which takes a picture of the DBS card before and after each extraction to assure sample traceability. All pressures on the DBS card sampler as well as on the LC-MS/MS unit are monitored and documented.

The camera documentation system checks for preset values, where quality control parameters can be integrated. If the quality criteria are not met, the system automatically checks the next spot and continues the preset program. In addition, the camera detects already extracted spots and blocks those to prevent reanalysis. The system checks the xy shift of the circle on the card and the blood spot to center the extraction spot to the middle of each DBS (to avoid inhomogeneous distribution effects). Further, the inbuilt pressure sensors are monitoring the extraction and rinsing pressure (in this study, the extraction pressure was 0.5 bar and the rinsing pressure 45 bar). The maximum extraction pressure is 1.2 bar, depending on the age of the blood sample, and the system can be rinsed with pressure up to 100 bar to prevent carry over [9].

Results

MS/MS Method

The method from the MassChrom kit was transferred to the automated DBS-MS 500 platform allowing full automation of this process. The method was successfully installed at a routine NBS laboratory. Five calibrator levels of blood spotted on DBS cards were established for the method validation and installation. For routine, the two levels high and low from the MassCheck controls were used. Since the extract is guided online to the tandem mass spectrometry system, it only takes 2 minutes per sample. All processes are overlapping, where the DBS-MS 500 runs a wash program and extracts the consecutive sample while the MS/MS is detecting and reporting the target analytes of the previous sample.

Method validation

First, validation including additionally carry-over, interference, volume effects and recovery was performed after setting up the method on the automated platform. Evaluation of accuracy was not part of this method validation, since control material from the MassChrom kit was used as reference in all experiments. The intra- and inter-day variations as well as the coefficients of determination (R2) of the calibration standards are summarized in Table 2. All factors were within the general acceptance criteria for NBS screening methods [17]. Intra-day variations are all below 15% with the exception of aspartic acid and glutamic acid which are slightly above but still acceptable since the criteria sets a value of < 25%. Also, glutamic acid gave a relative high variation for the inter-day variation (still within the criteria), compared to the other values which are all below 25%. Therefore the optimization of the MRM transition of glutamic acid was carefully monitored for the method installation at the routine site.

Table 2. Selected recovery comparison to MassCheck controls

MassCheck

Methionine

Proline

Valine

Carnitine

Laurovl-carnitine

Extraction 1

0.6159

1.421

1.597

15.499

56.966

Extraction 2

0.1284

0.460

0.527

2.416

9.854

Extraction 3

0.0557

0.154

0.184

0.653

4.461

Extraction 4

0.0276

0.063

0.083

0.293

3.981

Extraction 5

0.0154

0.031

0.035

0.152

3.368

[%]

[%]

[%]

[%]

[%]

MassCheck recovery

73.1

66.7

65.8

81.5

72.4

Sample recovery

79.9

70.3

61.5

80.0

76.4

DBS spots using 15, 30 and 50 μl blood were prepared at calibrator level A, C and E and compared as triplicates. The deviations from 50 μl to 30 μl spots were between 90.8 – 102.2 % and the deviations from 30 μl to 15 μl were between 78.2 – 88.8 %. There is a trend of smaller signals towards the 15 μl spots, however there is no obvious trend between 30 μl and 50 μl spots.

Carry-over from the high concentration to a blank DBS card passed the ICH guideline criteria and all investigated filter paper types were feasible for being used with this method. The best results were accomplished with the TFN filter paper (2.6% variation), paper types 226, 903, 229 were below 5% variation and the 818 paper moved with 6.0% on the last rank (supplementary data). The standard deviation of the complete panel from the triplicate measurement was taken into account.

The recovery of the compounds from the DBS card lies within 40 to 80%. This recovery remains constant for each analyte with the chosen extraction parameters. Since the results are always considered relative to the used MassChrom quality control card values with the same consistent recovery, this has no effect to the screening method, (see Table 2 and supplementary data).

Also the extraction behavior remains constant. This was investigated by comparison of freshly spotted blood and blood spots with 4 days of age stored at -20 °C (see supplementary data).

DBS spots using 15, 30 and 50 μl blood were prepared at calibrator level A, C and E and compared as triplicates. The deviations from 50 μl to 30 μl spots were between 90.8 – 102.2 % and the deviations from 30 μl to 15 μl were between 78.2 – 88.8 %. There is a trend of smaller signals towards the 15 μl spots, however there is no obvious trend between 30 μl and 50 μl spots.

Method Transfer into a Routine Environment

The previously developed and validated method was transferred into a routine environment. A short validation focusing on correlation and robustness was performed, where all five calibration levels were measured four-fold on 903 filter paper on three consecutive days. Each day was referred to high and low Mass Check control measured before the sample run.

Table 3. Intra-day, inter-day precision and R2 of calibrators level B (L1) and E (L2) from the spiked DBS samples (*Xle refers leucine/ Isoleucine)

Intraday day 1

Correlation

Intraday day 2

Correlation

Intraday day 3

Correlation

Inter day

LI [%]

L2 [%]

R2

LI [%]

L2 [%]

R2

LI [%]

L2 [%]

R2

LI [%]

L2 [%]

Alanine

6.8

7.6

0.994

3.7

7.1

0.981

8.6

5.4

0.998

24.5

13.0

Arginine

6.1

4.1

0.998

8.1

4.6

0.998

7.8

1.9

0.997

8.5

5.7

Aspartic acid

15.7

6.8

0.964

13.1

6.8

0.987

5.3

10.9

0.925

13.9

10.4

Citrulline

4.8

5.9

0.999

6.5

5.3

0.998

6.5

4.3

0.999

7.9

8.4

Glutamic acid

15.5

7.2

0.998

1.4

4.3

0.997

3.5

3.4

0.995

32.4

11.7

Glycine

9.2

9.9

0.998

9.2

9.9

0.998

7.7

12.3

0.987

8.7

18.9

Leucine (Xle*)

6.5

4.2

0.999

2.8

5.0

0.999

5.0

3.5

0.901

12.3

9.5

Methionine

10.3

4.8

0.999

9.6

6.1

0.985

11.9

4.8

0.999

24.9

5.8

Ornithine

6.9

4.0

0.999

6.5

5.7

0.998

11.1

4.6

0.996

13.3

14.8

Phenylalanine

5.7

3.9

0.999

4.3

4.9

0.999

8.6

8.3

0.998

11.7

6.0

Proline

4.0

4.1

0.999

4.6

5.0

0.999

11.0

8.8

0.999

9.4

6.9

Tyrosine

7.8

6.8

0.999

14.5

5.4

0.991

14.1

5.3

0.997

16.2

11.1

Valine

6.0

4.3

0.999

4.7

4.2

0.999

4.8

2.9

0.999

6.5

6.3

Carnitine

6.9

8.9

0.999

1.7

9.4

0.999

6.1

8.1

0.999

6.6

15.5

Acetvlcarnitine

9.3

2.8

0.978

9.7

6.3

0.986

8.0

6.2

0.983

9.4

7.4

Propionylcarnitine

2.3

2.9

0.998

1.5

3.3

0.999

3.8

3.6

0.996

12.9

3.3

Butyrylcarnitine

8.8

2.9

0.999

6.3

4.0

0.998

10.2

8.1

0.999

11.1

7.8

Valerylcarnitine

6.3

3.5

0.999

6.3

4.4

0.999

7.6

4.3

0.996

8.0

6.2

Hexanoylcarnitine

9.2

3.1

0.999

5.0

5.6

0.998

6.5

7.7

0.998

15.0

7.0

Octanoylcarnitine

5.7

3.5

0.999

6.6

3.6

0.999

5.4

7.4

0.997

10.0

6.7

Decanoylcarnitine

5.0

2.4

0.978

3.7

12.5

0.998

7.0

4.5

0.997

15.4

9.4

Laurovlcarnitine

14.0

4.2

0.999

11.5

4.2

0.999

7.8

3.5

0.999

12.8

4.2

Myristovlcarnitine

3.2

5.0

0.997

3.4

3.3

0.998

5.9

6.6

0.999

10.8

5.6

Palmitoylcarnitine

1.6

2.7

0.996

1.3

1.8

0.997

1.3

2.5

0.996

14.8

10.6

Stearoylcarnitine

2.3

2.8

0.995

1.8

2.4

0.999

2.1

3.2

0.995

18.1

14.1

At the Children’s hospital in Zurich, the Shimadzu MS/MS 8050 was coupled to a DBS-MS 500 unit as used for the method development. The criteria were met for intra- and inter-day variation following the validation procedure (supplementary data).

Discussion

The described method is fully automated and uses exclusively certified consumables and reference material. Through this degree of standardization, the application can be directly transferred in-between newborn screening laboratories and therefore following the trend towards the development of standardized programs. The method represents a high throughput application with only 2 minutes per sample. All processes are well documented by the reporting system of the DBS- MS 500, where a picture is taken of each spot before and after the extraction and all pressures for extraction, rinsing and LC pump are monitored. The DBS card picture is analyzed with image recognition software, providing results for the spot diameter, area and roundness, which can be used as quality control criteria. Also, the system detects a spot which has been previously extracted, to avoid multiple extractions of the same area. In a standard setting, the extraction is performed automatically from the center of each spot to overcome inhomogeneous distribution effects within the DBS sample [18].

We encountered a relative high endogenous concentration of the target analytes in the available donor blood. To bring this into the range of interest, the blood was washed with saline following a protocol. This procedure allowed gently spiking the washed red blood cells with the target analytes in the according concentrations. The five calibration points were used for method development, validation and installation at the routine site. After implementing the method, all results were referred to externally prepared and certified reference cards from the MassChrom kit.

The protocol of the MassChrom kit describes an extraction of a 3.2 mm disc in 100 μl extraction buffer (including the internal standard). Here, a 4 mm area was sealed on the DBS card and extracted with 60 μl into a 20 μl loop, whereas the first 40 μl are guided into waste. This was experimentally optimized using different loop volumes, extraction volumes and flow rates. This outcome could be due to an initially high portion of certain analytes, such as salts and phospholipids, causing matrix effects and ion suppression in the ESI source. However, the basic principle of this result needs to be further investigated. The method uses 60 μl of extraction buffer per sample instead of 100 μl from the protocol, which reduces the amount of solvent and analytical standard. In addition, consumables for sample preparation according to the protocol such as Eppendorf vials and pipette tips are no longer required.

Leucine, isoleucine, hydroxy l proline and allo-isoleucine are detected as sum (Xle) and are not separated within this method. An abnormal result for this parameter will have the automatic consequence of a second tier method using an analytical column prior the mass spectrometer to properly distinguish between all isobars. This is needed for the detection of maple syrup urine disease (MSUD).

A volume of 30 μl blood was chosen for the preparation of the calibrator spots. The 50 μl standard volume from literature resulted in relatively large spots on the standard filter paper card for automation with a spot-to-spot distance of 13.7 mm. The 50 μl spots were overlapping if all four positions on the DBS card were spotted. A comparison of 30 and 50 μl droplets showed no difference in the results.

Conclusion

A validated reagent kit for the extraction and analysis of dried blood spots in the field of newborn screening was transferred and validated on an automated DBS card extraction platform. First, a validation was performed and the method was then successfully transferred into a routine environment. Each process step is well documented and all analysis steps follow Good Laboratory Practice (GLP) [19]. The method can be easily modified or extended and transferred to other routine laboratories.

Author Contributions

Stefan Gaugler and Vincente Luis Cebolla Burillo participated in the concept and design of the study. Jana Rykl contributed to the optimization of the MS methods. Stefan Gaugler and Jana Rykl did the method development work, Stefan Gaugler did the practical work at the Children’s Hospital Zurich. Stefan Gaugler was responsible for drafting the manuscript. Vincente Luis Cebolla Burillo was responsible for revising the manuscript. All authors approved the manuscript as submitted.

Conflict of Interest

Jana Rykl is an employee of Shimadzu Schweiz GmbH (Reinach, Switzerland) and Stefan Gaugler is an employee of CAMAG (Muttenz, Switzerland). None of the other authors report any conflict of interest regarding this study. All instrumentations are property of the respective companies and laboratories.

References

  1. G. Nys, M. G. M. Kok, A. C. Servais, and M. Fillet, “Beyond dried blood spot: Current microsampling techniques in the context of biomedical applications,” TrAC – Trends Anal. Chem., vol. 97, pp. 326–332, 2017.
  2. J. L. Merritt, L. L. Brody, G. Pino, and P. Rinaldo, “Newborn screening for proximal urea cycle disorders: Current evidence supporting recommendations for newborn screening,” Mol. Genet. Metab., no. February, pp. 0–1, 2018.
  3. H. Lemonde, “Newborn screening for inborn errors of metabolism,” Paediatr. Child Health (Oxford)., vol. 25, no. 3, pp. 103–107, 2014.
  4. U. Groselj, M. Zerjav, A. Smon, N. Angelkova, D. Anton, I. Baric, M. Djordjevic, L. Grimci, M. Ivanova, A. Kadam, V. Mulliqi, H. Maksic, O. Marginean, O. Margineanu, O. Milijanovic, F. Moldovanu, M. Muresan, S. Murko, M. Nanu, B. Repic, M. Samardzic, V. Sarnavka, A. Savov, M. Stojiljkovic, B. Suzic, R. Tincheva, H. Tahirovic, A. Toromanovic, N. Usurelu, and T. Battelino, “Newborn screening in southeastern Europe,” Mol. Genet. Metab., vol. 113, no. 1–2, pp. 42–45, 2014.
  5. A. W. El-Hattab, M. Almannai, and V. R. Sutton, “Newborn Screening: History, Current Status, and Future Directions,” Pediatr. Clin. North Am., vol. 65, no. 2, pp. 389–405, 2018.
  6. A. Starkweather, B. Coleman, V. B. De Mendoza, M. Fu, J. Taylor, W. Henderson, C. Kenner, D. Walker, L. Amankwaa, and C. Anderson, “American Academy of Nursing on Policy Policy brief : Improve coverage of newborn genetic screening to include the Recommended Uniform Screening Panel and newborn screening registry,” Nurs. Outlook, vol. 65, no. 4, pp. 480–484, 2017.
  7. M. Resano, M. A. Belarra, E. García-Ruiz, M. Aramendía, and L. Rello, “Dried matrix spots and clinical elemental analysis. Current status, difficulties, and opportunities,” TrAC – Trends Anal. Chem., vol. 99, pp. 75–87, 2018.
  8. R. Fingerhut, M. L. Silva Polanco, G. D. J. Silva Arevalo, and M. A. Swiderska, “First experience with a fully automated extraction system for simultaneous on-line direct tandem mass spectrometric analysis of amino acids and (acyl-)carnitines in a newborn screening setting,” Rapid Commun. Mass Spectrom., vol. 28, no. 8, pp. 965–973, 2014.
  9. S. Gaugler, J. Rykl, I. Wegner, T. von Däniken, R. Fingerhut, and G. Schlotterbeck, “Extended and Fully Automated Newborn Screening Method for Mass Spectrometry Detection,” Int. J. Neonatal Screen., vol. 4, no. 1, p. 2, 2017.
  10. “CAMAG DBS,” 2018. [Online]. Available: http://www.camag.com/en/dbs/what_is_dried_blood_spot_sampling.cfm. [Accessed: 29-May-2018].
  11. M. H. U. Duthaler, B. Berger, S. Erb, M. Battegay, E. Letang, S. Gaugler, S. Krähenbühl, “Automated high throughput analysis of antiretroviral drugs in dried blood spots.,” J Mass Spectrom., vol. 52, no. 8, pp. 534–542, 2017.
  12. S. Gaugler, J. Rykl, M. Grill, and V. L. Cebolla, “Fully automated drug screening of dried blood spots using online LC-MS / MS analysis,” J. Appl. Bioanal., vol. 4, no. 1, pp. 7–15, 2018.
  13. D. H. Chace, B. W. Adam, S. J. Smith, J. R. Alexander, S. L. Hillman, and W. H. Hannon, “Validation of accuracy-based amino acid reference materials in dried-blood spots by tandem mass spectrometry for newborn screening assays.,” Clin. Chem., vol. 45, no. 8 Pt 1, pp. 1269–77, 1999.
  14. G. L. Hoffman, “CLSI guideline NBS04-A,” 2010.
  15. P. Abu-Rabie, P. Denniff, N. Spooner, B. Z. Chowdhry, and F. S. Pullen, “Investigation of different approaches to incorporating internal standard in DBS quantitative bioanalytical workflows and their effect on nullifying hematocrit-based assay bias,” Anal. Chem., vol. 87, no. 9, pp. 4996–5003, 2015.
  16. N. Kostić, Y. Dotsikas, N. Jović, G. Stevanović, A. Malenović, and M. Medenica, “Quantitation of pregabalin in dried blood spots and dried plasma spots by validated LC-MS/MS methods,” J. Pharm. Biomed. Anal., vol. 109, pp. 79–84, 2015.
  17. Chromsystems, “MassChrom Amino acids and Acylcarnitines from Dried Blood (non derivatised), 57000, EN 06/2016 R7,” 2016.
  18. G. Nys, M. Kok, A.-C. Servais, and M. Fillet, “Beyond dried blood spot: current microsampling techniques in the context of biomedical applications,” Trends Anal. Chem., vol. In press, 2017.
  19. European Medicines Agency, “Guideline on bioanalytical method validation,” 2015.

Supplementary data

Filter paper comparison                                                                     

226

Average

3.7 %

TFN

Average

2.7 %

903

Average

3.6 %

2992

Average

2.9 %

903NBS

Average

2.6 %

818

Average

6.0 %

Calibration sample preparation

Spike A

Spike B

Spike C

Spike D

Spike E

[ul]

[ul]

[ul]

[ul]

[ul]

AA, CO/I Spike

25

50

100

167

250

AC Spike

5

10

25

50

100

H20

225

200

150

83

0

MeOH

95

90

75

50

0

NaCl 0.9%

2150

2150

2150

2150

2150

Red blood cells

2500

2500

2500

2500

2500

Total volume

5000

5000

5000

5000

5000

 

Carry over (Area blank/Area E * 100%)

Recovery

Alanine

2.1

55.6

%

Arginine

0.1

63.3

%

Aspartic acid

1.8

74.4

%

Citrulline

0.1

66.5

%

Glutamic acid

0.1

44.0

%

Glycine

1.8

68.4

%

Leucine (Xle)

0.1

65.6

%

Methionine

0.1

76.9

%

Ornithine

0.1

64.5

%

Phenylalanine

0.3

68.3

%

Proline

0.3

70.3

%

Tyrosine

0.3

55.1

%

Valine

0.3

61.5

%

Carnitine

0.1

80.0

%

Acetylcarnitine

0.1

78.7

%

Propionylcarnitine

0.1

50.1

%

Butyrvlcarnitine

0.1

81.0

%

Valerylcarnitine

0.1

81.7

%

Hexanoylcarnitine

0.3

75.6

%

Octanoylcarnitine

0.6

77.4

%

Decanoylcarnitine

0.8

72.3

%

Lauroylcarnitine

0.9

76.4

%

Mvristoylcarnitine

0.1

70.8

%

Palmitoylcarnitine

0.2

40.9

%

Stearoylcarnitine

0.2

43.4

%

Intra-day, inter-day precision and R2 of calibrators level B (LI) and E (L2) from the spiked DBS samples in the routine laboratory

Intraday day 1

Correlation

Intraday day 2

Correlation

Intraday day 3

Correlation

Interday

LI [%]

L2 [%]

R2

LI [%]

L2 [%]

R2

LI [%]

L2 [%]

R2

LI [%]

L2 [%]

Alanine

13.3

4.8

0.980

6.9

6.6

0.994

3.8

10.0

0.994

21.9

9.1

Arginine

6.5

4.0

0.999

5.0

2.0

0.998

4.7

3.2

0.999

6.9

4.5

Aspartic add

7.9

7.1

0.996

13.9

9.5

0.996

11.4

5.5

0.962

19.7

13.3

Citrulline

6.0

6.8

1.000

6.7

4.0

0.999

7.4

5.0

0.998

6.9

5.5

Glutamic add

1.0

5.4

0.994

9.8

2.9

0.993

10.8

7.0

0.999

19.6

14.2

Glycine

14.7

14.7

0.994

10.8

14.5

0.980

10.3

14.1

0.991

15.0

23.6

Leucine (Xle)

7.5

3.1

0.999

2.4

2.5

0.998

4.7

2.4

0.999

11.6

9.7

Methionine

13.2

5.0

0.997

14.3

7.3

0.997

8.4

9.5

0.999

21.3

13.7

Ornithine

5.9

2.9

0.997

15.2

9.3

0.994

9.5

6.7

0.994

12.2

13.6

Phenylalanine

9.9

7.6

0.999

4.9

3.4

0.999

6.1

3.6

1.000

7.9

14.8

Proline

9.9

5.2

0.999

4.3

2.8

0.998

5.8

5.7

0.999

7.4

6.2

Tyrosine

11.1

17.0

0.945

14.3

3.0

0.999

13.3

5.7

0.998

21.6

20.7

Valine

10.2

8.6

0.997

12.4

2.1

0.998

3.9

3.0

0.997

11.7

9.1

Carnitine

3.6

4.4

0.997

5.2

8.8

0.998

7.9

11.6

0.997

7.2

27.5

Acetylcarnitine

8.1

3.9

0.989

6.5

5.7

0.995

13.2

6.1

0.993

9.7

12.7

Propionylcarnitine

1.3

2.4

0.998

3.1

3.0

1.000

3.7

3.9

0.999

11.7

3.6

Butyrylcarnitine

6.0

5.4

0.997

8.7

4.5

0.999

11.8

5.1

0.999

9.2

10.4

Valerylcarnitine

7.3

5.5

0.999

6.7

3.8

1.000

9.6

5.3

0.999

8.5

6.5

Hexanoylcarnitine

9.7

7.1

0.999

9.6

4.8

0.999

11.2

6.2

0.999

16.3

6.9

Octanoylcarnitine

8.1

4.6

0.999

5.8

5.0

0.999

5.3

3.6

0.999

7.4

4.9

Decanoylcarnitine

15.3

6.1

0.999

5.1

3.4

0.999

6.0

2.3

0.999

10.1

4.2

Lauroylcarnitine

4.7

4.5

0.999

3.7

2.1

0.999

6.3

4.3

0.999

7.5

4.6

Myristoylcarnitine

6.9

3.8

0.998

6.2

3.2

0.998

2.6

5.9

0.997

6.1

5.1

Palmitovlcarnitine

3.9

4.7

0.995

8.2

11.4

0.896

6.7

5.5

0.992

8.8

7.8

Stearoylcarnitine

4.8

2.3

0.985

17.4

14.4

0.973

19.0

15.3

0.954

34.9

18.5

The Assessment of Quality of Life Using GOHAI among Edentulous Patients

DOI: 10.31038/JDMR.2019211

Abstract

Introduction: Edentulism may compromise the quality of life (QoL) of a patient. Geriatric patients who are satisfied with their complete dentures are usually satisfied with their daily life. This research was conducted to assess quality of life using Geriatric Oral Health Assessment Index (GOHAI) focusing on functional limitation, oro-facial pain, psychological impact and behavioural impact among edentulous patient attending USIM Polyclinic.

Materials and Methods: Thirty-two geriatric edentulous patient to be treated by third year students in polyclinic USIM were recruited. Patient completed the Malay version of GOHAI questionnaire.

Result: The data showed there were 62.5% male and 37.5% female. The mean age was 63.88 ± 5.621 year. 87.5% of the samples are wearing complete denture and 56.3% of them have been wearing denture for the past five years. The lowest mean score was 1.66 ( ± 1.6) for trouble in biting and chewing. The highest mean for GOHAI was in medication used to relieve pain with the mean score of 4.56 ( ± 0.8).

Conclusion: The QoL of geriatric patient was mostly affected by the oral function. However the least affected was seen in pain and discomfort. Despite wearing complete denture, denture replacement is essential to improve QoL.

Keyword

Edentulous and Quality of life, GOHAI, Geriatric Dentistry

Introduction

According to the Gerodontology society, the life expectancy of the elderly who need complete denture pros-thesis has increased [1]. Edentulism is considered as an outcome of poor oral healthcare and may compromise the quality of life [2]. The physiological changes in edentulous arch with continuous resorption of the residual ridges leads to the need of complete denture replacements. The use of complete denture increases the quality of life of geriatric patient in term of oral function, reduce oro-facial pain, improve oro-facial appearance, and psychosocial impact [3]. Oral health is often assessed separately from patient’s general health. Sometimes it perceived as a distinct entity from any chronic conditions of the elderly patient. Edentulism may compromise the quality of life (QoL) because patients will have difficulty eating proper food and this might lead to malnourish among the elderly. Edentulous patients need dentures to restore their oral function. Satisfaction of the denture influences psychological condition which is part of quality of life (QoL). Geriatric patient who are well satisfied with their daily lives, usually also satisfied with their complete denture [4].

There is a measure called Geriatric/General Oral Health Assessment Index (GOHAI) is an instrument consists of 12-item questionnaires and it is considered a gold standard used in measuring the oral health impact of the geriatric patients. It was developed by Atchinson and Dolan [4] and has been validated over a few languages including Malay language  [5] which is useful to measure oral-health related quality of life in Malaysia. From its name, GOHAI Index is usually used for the assessment of the quality of life in geriatric patient and not suitable to assess the overall oral health status because there was a weak correlation with the clinical measure of disease. Therefore, GOHAI cannot be used to conclude a diagnosis of dental disease because it will never assess the oral examination clinically and radiographically. However, indication for a referral or further thorough oral examination need can be obtained through GOHAI assessment as it will provide necessary information in symptoms as well as the limitation in psychosocial and functional. Therefore the aim of this paper was to assess the quality of life of the elderly population among the geriatric patients that received complete dentures from the undergraduate dental students of the Faculty of Dentistry, Islamic Science University of Malaysia (USIM) Polyclinic, Kuala Lumpur.

Methods

An approval was attained before starting the study by Board of Research Ethic of Islamic Science University of Malaysia (USIM). A total of 32 community dwelling participants were recruited among complete denture patients who will receive complete dentures from the undergraduate dental students of USIM. They must satisfy all the following inclusion criterias; Malaysians aged 60 years and above, complete denture wearers, subject with con-trolled or no underlying systemic diseases, and fluent in Malay language. The consents were obtained from par-ticipants after being informed about voluntary participation in this research. Then, self-administered questionnaire was distributed to collect information of participant’s demographic and prosthetic experience.

The participants were interviewed by a single interviewer and were asked to estimate the frequency of problems in Malay Version of GOHAI questionnaire using a six point Likert scale rating (always [5], very often [4], often [3], sometimes [2], seldom [1] or never [0]). The 12-item questionnaire was classified into four major domains which were the functional limitation, pain or discomfort, psychological and behavioural impact. The GOHAI score is determined by summing the final score of each of the 12 items ranges from 0 to 60. The score for GOHAI item number 3, 5 and 7 were maintained while the rest were reversed in order to attain a positive oral health GOHAI score. The higher GOHAI score denotes better oral health status perceived by the participants themselves. If there were 3 or more missing items in the data, the data will be dropped out. Alternately, if no more than two missing items in a data, the score of the missing will be substituted with the mean. The mean GOHAI scores in relation with demographic variables were analyzed using Independent t-test and ANOVA. 0.05 was set as the level of significance. The data were analyzed using IBM SPSS Statistic version 19.

Results

The patients’ age ranges from 60–76 years with mean age of 63.88 ± 5.6 (as shown in (Table 1) below. Almost one-third of them were female and more than half received education up to secondary education level (11 years of formal education in school). Majority of them have been wearing a complete denture for more than 5 years (56.3%). Most of the patients request for new dentures because they want to improve the chewing ability (50.0%) and aesthetic reasons (40.6%) over the other health and speech reasons. (Table 2)

Table 1. Distribution of patients according to sex.

Gender

Number

Mean

SD

Male

20

64.40

5.707

Female

12

63.00

5.608

Total

32

63.88

5.621

Table 2. Characteristics of subjects.

Participant’s characteristic

(N=32)

N

 %

Mean GOHAI score

SD

p-value

GOHAI Score

Range

40.1

7.4

(27–60)

Sex a

Male

Female

20

12

62.5

37.5

39.2

41.8

6.0

9.5

.346

Age b

<55 years

55–59 years

60–69 years

70–79 years

2

5

21

4

6.3

15.7

65.5

12.5

42.50

39.40

39.48

43.25

3.5

9.2

7.6

6.7

.789

Education Level b

No schooling

Primary school

Secondary school

University

2

6

19

5

6.3

18.8

59.4

15.7

42.0

43.0

40.4

34.8

4.2

11.7

6.4

4.0

.315

History of Denture Wearing b

< 1 year

1–3 year

3–5 year

>5 year

Never

4

3

3

18

4

12.5

9.4

9.4

56.3

12.5

39.3

44.3

41.0

40.1

37.3

10.2

14.0

3.6

7.0

4.2

.820

Reason for Wearing Denture b

Chewing

Aesthetic

Health

Speech

16

13

2

1

50.0

40.6

6.3

3.1

39.1

40.2

47.5

40.0

6.0

9.1

6.4

0.0

.539

The mean and median of the summary GOHAI scores was 40.1
( ± 7.4, range 27–60) and 39.0, respectively. Besides that, the skewness of the GOHAI was 0.75. There was only one out of 32 subjects scored the maximum in this measure. About one-third of the participants (31.3%) scored 41 and above. The mean of GOHAI scores in regard to each variable varies from 37.3 (± 4.2) to 47.5 (± 6.4). The lowest and highest mean GOHAI scores can be found in university level participants and reason of wearing denture for health purpose respectively (Table 2). Based on the data exhibited in Table 2, there is no significant association between gender, age, level of education, history of wearing denture and reason of wearing denture with the oral health related quality of life (OHRQoL). The used of medication to relieve pain in pain and discomfort domain was recorded as the highest mean for GOHAI item with 4.56 (SD 0.8) followed by 4.31 (SD 1.1) for sensitive to hot, cold or sweet foods and 4.31 (SD 1.3) for limit contact with people. The lowest mean was 1.66 (SD 1.6) in trouble biting and chewing. Meanwhile, Table 3 shows the percent of participants responding ‘sometimes,’ ‘often,’ ‘very often,’ or ‘always’ to each of the GOHAI items. The percentage of subjects responding positively to each item ranged from 16% to 80%, with ten out of 12 items being reported by 20% or more and six of 12 reported by 33.3% or more (Table 3).

Table 3. Percentage of subjects responding ‘sometimes,’ ‘often,’ ‘very of-ten,’ or ‘always’ to each GOHAI and mean each of the items.

GOHAI ITEMS

Percentage of participants responding positive response (%)*

Mean (SD)

Functional limitation

2. Trouble biting and chewing

3. Able to swallow comfortably

4. Problem to speak clearly

64

80

24

1.66 (1.6)

3.19 (1.5)

3.53 (1.6)

Pain and discomfort

5. Able to eat without discomfort

8. Used medication to relieve pain

12. Sensitive to hot, cold or sweet foods

52

16

20

2.38 (1.5)

4.56 (0.8)

4.31 (1.1)

Psychological impacts

7. Pleased with the look of teeth

9. Worried about teeth, gums or dentures

10. Self-conscious of teeth, gums or dentures

11. Uncomfortable eating in front of others

56

48

16

20

2.25 (1.4)

3.50 (1.5)

4.19 (1.4)

3.69 (1.5)

Behavioural impact

1. Limit the kinds of food

6. Limit contact with people

64

20

2.41 (1.8)

4.31 (1.3)

Discussion

TThis research was intended to evaluate the OHRQoL of the edentulous patients. There are many tools available to study the OHQoL, but the two most common measures used are GOHAI and OHIP-14. According to Ikebe et al (2012), GOHAI is more sensitive than OHIP-14, hence the choice of measures in this study. However, there was no measurements has proven that one measures is superior than the others [7]. Unfortunately, the relationship between gender, age, level of education, history of wearing denture and reason of wearing denture with the oral health related quality of life (OHRQoL) could not be established due to small sample size. For the improvement, larger sample size is needed in order to look for the significant association between those variables with the OHRQoL.

In this present study, the most prevalent impact found in the functional limitation and behavioural domains with 64% of the participants had responded positive response in both negatively worded items (trouble biting and chewing as well as limit the kinds of food). On the other hand, pain and discomfort domain (used of medication to relieve pain) found as the least prevalent impact in OHRQoL. The research conducted by Adam (2006) using OHIP-EDENT exhibited functional limitation as the most prevalent impact and psychological discomfort and physical pain came after as second and third most prevalent impact [8]. However, physical pain is the most prevalent impact according to Heydecke et al (2004) as the research used OHIP-DENT to measure OHRQoL of the samples [9].

In this study, the lowest mean for GOHAI item falls in functional limitation (trouble in biting and chewing) as it contributes a significance burden on individual as well as community. The complete denture patient experience more mastication problem compared to their dentate counterpart [8]. Meanwhile, an item in pain and discomfort domain which is “used of medication to relieve the pain” has highest mean. In other words, either the participant did not get bothered by pain while wearing denture or they were able to tolerate the pain. There was no standard classification on GOHAI score which indicates good or poor QoL. However, there was another classification used but it cannot be compare with this result because it use modified rating score. Hence, as the summary mean of GOHAI score is lower than 41, it means that almost all participants (68.7%) have relatively poor QoL and need treatment [5, 8].

Lastly, the main reason to seek complete denture treatment in this group of patients was to enhance chewing capability.

Conclusion

The QoL of geriatric patient was mostly affected by oral function and the least affected was seen in pain and discomfort. Thus, despite wearing complete denture, denture replacement is essential to improve Quality of life in geriatric patient.

Acknowledgement

This study was Funded by Short Term Grant Research Universiti Sains Islam Malaysia PPP/USG-0115/FPG/30/12715. We express our particular appreciation for the assistance of Azlan Jaafar, statistics advisor for this project.

References

  1. Mojon P (2003) The world without teeth: demographic trends. Fastest clincal dentistry insight engine Pg No: 18–21.
  2. Laurina L, Soboleva U (2006) Construction faults associated with complete denture wearers’ complains. Stomatologija, Baltic Dental and Maxillofac J 8: 61–64. [crossref]
  3. John MT, Feuerstahler L, Waller N, Baba K, Larsson P, et al. (2014) Confirmatory factor analysis of the Oral Health Impact Profile. J Oral Rehabil 41: 644–652. [crossref]
  4. Atchison KA, Dolan TA (1990) Development of the Geriatric Oral Health Assessment Index. J Dent Educ 54: 680–687. [crossref]
  5. Othman WN, et al. (2012) Validation of the Geriatric Oral Health Assessment Index (GOHAI) in the Malay lan-guage. Association between self-assessment of complete dentures and oral health-related quality of life. J Oral Rehabilitation 39: 847–857.
  6. Ikebe K, Hazeyama T, Enoki K, Murai S, Okada T, Kagawa R, et al. (2012) Comparison of GOHAI and OHIP-14 measures in relation to objective values of oral function in elderly Japanese. Community dentistry and oral epidemiology 40: 406–414. [crossref]
  7. Locker D, Matear D, Stephens M, Lawrence H, Payne B (2001) Comparison of the GOHAI and OHIP-14 as measures of the oral health-related quality of life of the elderly. Community Dent Oral Epidemi-ol 29: 373–381. [crossref]
  8. Adam RZ (2006) Do complete dentures improve the quality of life of patients? (Doctoral dissertation, University of the Western Cape).
  9. Heydecke G, Tedesco LA, Kowalski C, Inglehartet MR (2004) Complete dentures and oral health-related quality of life – do coping styles matter? Community Dent Oral Epidemiol 32: 297–306.

Vulnerability of Parkinson Patients to Iatrogenic Adverse Events: Emphasis on the Perioperative Period through a Case Report

DOI: 10.31038/JNNC.2018122

Abstract

Deep brain stimulator surgery represents an interesting option for medically resistant Parkinson’s disease patients, indisposed by “wearing-off” and motor fluctuations. Considering the neurodegenerative nature of PD, after approximately five years of natural disease evolution, one half to two-thirds of patients develops this delayed complication. The occurrence of dyskinesia and motor fluctuations interferes with activities of the daily living, consequently impacting the quality of life. A high frequency chronic stimulation targeting the subthalamic nucleus, located ventral to the thalamus and involved in the basal ganglia system, generates the same outcome than an ablation of this structure. Despite the positive impact of DBS surgery on PD symptomatology, physicians should be mindful of the vulnerability of this particular population and promote thorough monitoring, especially during the perioperative phase. A thoughtful approach considering the pathophysiological mechanisms underlying PD is required in order to select an individual’s pharmacopoeia. We report the case of a 61 year-old male, known for seventeen year of evolving PD, who developed acutely a state of generalized symmetrical rigidity, with a predilection for axial and bulbar musculature. After a meticulous review of this patient’s pharmacologic profile, sufentanil was deemed the responsible agent.

Keywords

Sufentanil, Parkinson Disease; Deep Brain Stimulator, Wearing-Off, Motor Fluctuations, Dopamine

Introduction

After several years of Parkinson’s Disease (PD) evolution, most patients develop a wearing-off effect or motor fluctuations induced by chronic exposure to exogenous levodopa. Despite the introduction of dopaminergic agonists, Catechol-O-Methyl Transferase Inhibitors (COMTI) or amantadine, the efficacy of the Deep Brain Stimulator (DBS) on this delayed complication is undeniable. On the other hand, PD patients are prone to serious adverse events justifying cautious choices in the selection of an anesthetic regimen, especially in the perioperative period.

Case Description

A 61 year-old male, known for seventeen year of evolving PD, underwent bilateral deep brain stimulator implantation. The surgery comports three systematic steps. The first stage consisted of a right temporal lead implantation. This procedure is repeated on the contralateral side on the following day. Finally, the last stage consists in the stimulator internalization.

The preoperative neurologic examination revealed a severe resting tremor implicating bilateral upper and lower limbs. The speech was significantly compromised by the facial rigidity and bradykinesia. Supplemental pertinent findings included: chin tremor, diffuse and global appendicular cogwheel rigidity, postural instability and a typical parkinsonian gait. While our patient was “en route” to the operation room (OR), he was administered his regular dopaminergic regimen. An uneventful general anesthesia involved the following agents: sufentanil, propofol, rocuronium, and lidocaine. No intra-operative anti-dopaminergic medication was judged necessary. Upon arrival to ICU (Intensive Care Unit), our gentleman developed an acute severe extra-pyramidal reaction immediately following extubation. The main concern at this point was the limitation of thoracic expansion, at imminent risk to evolve towards respiratory failure. The administration of benzodiazepines was insufficient to reach stable respiratory parameters, justifying an emergent reintubation. The acute iatrogenic reaction induced dysphagia, preventing the administration of oral medication such as levodopa-carbidopa, mandating the installation of a nasogastric tube. A direct correlation was observed between the improved chest wall compliance and the expected washout period of opioids.

Discussion

PD, described initially as the “Shaky Palsy” by James Parkinson in the 19th century, affects 0.3% of the overall population [1]. The endpoint of PD consists of progressive vanishing of the dopaminergic neurons within the substantia nigra pars compacta. This process induces an imbalance between dopaminergic inhibition and cholinergic excitation of the striatal output. Dopaminergic depletion subsequently induces excessive thalamic inhibition. The phenotypic repercussion of dopamine depletion is classical “parkinsonism”, characterized by: akinesia/bradykinesia, rigidity, tremor and postural instability.

Electrophysiological mapping mandatorily requires the patient’s cooperation. Medication preventing hyperkinetic symptoms is contraindicated in the immediate postsurgical state. The priority is to depict those involuntary movements, to assure a proper positioning of the leads and confirm the success of the procedure. Leads implantation is performed in two distinct steps. Firstly, the stereotaxic frame and leads are positioned. Local analgesia such as a scalp bloc is preferred to general anaesthesia, to maintain the patient’s cooperation. General anesthesia is performed during the second stage, consisting of the internalization of the leads and implantation of the programmable impulse generator.

Every anaesthetic agent comports pros and cons, an individualized “case-by-case” approach is preferable. Opioids should be administered cautiously in general, and even more in parkinsonian individuals considering PD’s pathophysiological basis. The striatum, representing a major source of input to the basal ganglia, contains a high concentration of opioid receptors, modulating reward, addiction and motor control [2]. Opioids neuropeptides, such as enkephalin and dynorphin. contribute to neuromodulation. A divergent form of plasticity leads to the remodeling of the nigrostriatal network, as a consequence of dopamine depletion and its treatment: exogenous L-dopa. The outcome is an upgrade of opioid transmission involving the basal ganglia, possibly preventing the development of motor dyskinesias [3]. In parallel, fentanyl and sufentanil are powerful opioids reported to induce a rare complication: “The Wooden Chest Syndrome” [4]. Several hypotheses were emitted to explain the opioid-induced rigidity resulting from narcotic exposure. PD patients are particularly vulnerable, considering their unique intrinsic features, mentioning alteration in the central mu1-opioid receptors, cerulospinal noradrenergic and glutamatergic pathways, as well as spinal motoneurons [5] (Table 1). This complication can le lethal, as increased muscular tone incapacitates the closure of the vocal cords, glottis and jaw, potentially accompanied by thoraco-abdominal rigidity impairing ventilation.

Table 1. Hypothetical Pathways and Mechanisms Predisposing Parkinson’s Disease Patients to Develop “Wooden Chest Syndrome.

Structure/Pathway

Mechanism

Striatal opioid binding sites

Interaction with dopaminergic D2-receptors.

Neurotransmitters imbalance

Reduced dopamine level results in imbalance with the cholinergic neurotransmitter system; both systems impact the control of muscle tone, leads to rigidity.

Striatal dopaminergic neurons

Dopamine depletion results in nigro-striatal dysfunction

Inhibition of tyrosine hydroxylase in the nigro-striatal system

Induced opioids, tyrosine hydroxylase is an essential enzyme acting on dopamine sysnthesis.  Overall impact is increased muscle tone/rigidity.

Decline in pallidal GABA level

Due to interconnection of the dopaminergic-gabaminergic systems, decline in the pallidal GABA concentration.  Overall impact is activation of cholinergic neurons projecting to thalamus.

Dysregulation of GABA concentration in the putamen

Impact of opioids: reduced gabaminergic transmission, improved by benzodiazepines.

Central mu1 opioid receptors

Increase in efferent motor signal, leading to increased muscle contraction and rigidity.

Locus ceruleus descending signal

Cerulospinal noradrenergic and glutamatergic impact of spinal neurons.

Pons descending signal

Activation of spinal motoneurones.

Conclusion

In conclusion, DBS surgery could be a life changing procedure for medically resistant PD candidates selected thoroughly based on their symptomatology, and realistic expectations. Anaesthetic drugs should be prescribed judiciously to prevent deleterious iatrogenic events. The importance of this case relies on three main factors: 1) the increasing popularity of DBS surgery, 2) the high prevalence of PD in our society, and 3) the high likelihood for any individual to necessitate general anesthesia in the context of a surgical procedure in a lifetime. Those elements highlight the necessity to develop optimal knowledge regarding perioperative pharmacotherapy in PD context.

Highlights

  1. The endpoint of Parkinson’s disease is a progressive vanishing of the dopaminergic neurons within the “substantia nigra pars compacta”. The clinical picture resulting from this ongoing depletion is named “parkinsonism”, characterized by: akinesia/bradykinesia, rigidity, tremor and postural instability.
  2. Despite the introduction of dopaminergic agonists, catechol-O-methyl transferase inhibitors or amantadine, the efficacy of subthalamic deep brain stimulator implantation is undeniable to improve the course of dyskinesia and motor fluctuations, typically presenting after several years of disease evolution.
  3. High frequency chronic stimulation targeting the subthalamic nucleus generates the same outcome than an ablation of this structure, explaining the effectiveness of deep brain stimulator implantation in medically resistant Parkinson’s disease patients.
  4. Deep brain stimulator surgery represents an interesting option for medically resistant Parkinson Disease patients, indisposed by “wearing-off” and motor fluctuations, considering the neurodegenerative nature of Parkinson’s disease, resulting in incurable, indolent and progressive symptoms.
  5. It is essential to have in mind the pathophysiological mechanisms underlying Parkinson’s disease in order to select optimally given patient’s pharmacopoeia.
  6. Parkinson’s disease patients are particularly vulnerable during the perioperative phase. Several risk factors are linked to the modifications of the pharmacologic profile during that period. Local or generalized anaesthesia represent a supplemental risk of developing an iatrogenic adverse events or drug interaction.

Author Contribution

Catherine Maurice: Conceptualization, Drafting, Redaction, Supervision, Editing, Revision and Final Approval.

Kiran Grant: Redaction, Editing, Revision and Final Approval.

Austin M. Pereira: Redaction, Editing, Revision and Final Approval.

Yasser B. Abulhasan: Conceptualization, Drafting, Redaction, Supervision, Editing, Revision and Final Approval.

Abbreviations

DBS: Deep Brain Stimulator

PD: Parkinson’s Disease

COMT: Catechol-O-Methyltransferase

CNS: Central Nervous System

References

  1. Parkinson (1817) An essay on the Shaking Palsy. J Neuropsychiatry and Clinical Neurosciences 14: 223–236
  2. A.H. Erga, I. Dalen, A. Ushakova, J. Chung, C. Tzoulis, et al (2018) Dopaminergic and Opioid Pathways Associated with Impulse Control Disorders in Parkinson’s Disease. Front Neurol 9:109
  3. MA Cenci (2007) Dopamine Dysregulation of Movement Control in 1-DOPA-Induced Dyskinesia. Trends Neurosci 30: 236–243. [Crossref]
  4. C.P. Phua, A. Wee, A. Lim, J. Abisheganaden, A. Verma (2017)  Fentanyl-Induced Chest Wall Rigidity Syndrome in a Routine Bronchoscopy.  Respir Med Case Rep 20: 205–207.
  5. Dimitriou, I. Zogogiannis, D. Liotiri, F. Wambi, N. Tawfeeq, et al (2014) Impossible Mask Ventilation After an Unusually Dose Fentanyl-Induced Muscle Rigidity in a Patient with Essential Tremor: A Case Report and Review of the Literature. Middle East J Anaesthesiol 22: 619–622.

Development of Synthetic Cell Differentiation Agent Formulations for the Prevention and Therapy of Cancer via Targeting of Cancer Stem Cells

DOI: 10.31038/CST.2019412

Abstract

The association of methylation enzymes with telomerase constitutes a unique abnormality of cancer cells. This abnormality locks methylation enzymes in an exceptionally stable and active state so that hypomethylation of nucleic acids necessary for the cells to undergo Terminal Differentiation (TD) cannot take place. Human body produces metabolites that are able to eliminate telomerase from abnormal methylation enzymes of cancer cells to allow TD to proceed. Cell Differentiation Agent-2 (CDA-2) is a preparation of human metabolites from freshly collected urine, which has been approved for cancer therapy by the Chinese FDA. The effective components of CDA-2 are Differentiation Inducers (DIs) to target on the telomerase of abnormal methylation enzymes and Differentiation Helper Inducers (DHIs) which are the inhibitors of individual enzymes of ternary methylation enzymes. CDA-2 was very effective for the therapy of Myelodysplastic Syndrome (MDS), which is a disease attributable to Cancer Stem Cells (CSCs). We have previously carried out extensive studies on the DHIs of CDA-2.We are now focusing on the DIs of CDA-2 in order to formulate synthetic CDA for the prevention and therapy of cancer via targeting of CSCs.

DIs were purified from CDA-2 solution by procedures including differential solvent extraction, gel filtration, ion exchange chromatography, TLC, and HPLC. The mass of purified active preparation was determined by mass spectroscopy. DI activity was based on the Nitro Blue Tetrazolium (NBT) assay of HL-60 cells.

DIs of CDA-2 were found predominantly as acidic liposomal complexes extractable by dichloromethane. A good proportion of which became covalently linked to inactive carriers which were not soluble in dichloromethane, but soluble in alcohols. We have identified pregnenolone as a DHI of active liposomal complexes. After dissociation from pregnenolone, the active DIs of CDA-2 were not associated with UV absorption peaks of HPLC. We suspected that the active DIs might be acidic peptides derived from endogenous proteins, because we have previously found that acidic peptides of CDA-2 were active DIs. We, thus, randomly picked pentapeptides containing at least two acidic amino acid residues from the sequences of a- and b-hemoglobin for synthesis to test their DI activities. Indeed, acidic pentapeptides of hemoglobin were active as DIs, although the activities were not impressive. Retinoic Acid (RA) and 12-O-TetradecanoylPhorbol-13-Acetate (TPA) are well known DIs with much better activities.

In this study, we found Pyrvinium Pamoate (PP) as the best DHI, and triinosinate + tetrainosinate (I3 + I4) as an acceptable DHI.

With effective DIs and DHIs on hand, our deliberated CDA formulations were as followings: for the therapy of MDS, the CDA-MDS formulation was RA(ED25)-5P-1(ED25)-I3 + I4(RI0.5)-PP(RI0.5)-sodium pregnenolone sulfate(RI0.5); for the therapy of CSCs, the CDA-CSC formulation was RA(ED25)-TPA(ED25)-PP(RI0.5)-resveratrol(RI0.5)-curcumin(RI0.5); for the therapy of brain tumor, the CDA-BT formulation was TPA(2xED25)-PP(2xRI0.5)-sodium phenylbutyrate(RI0.5)-pyrogallol(RI0.5); and for the therapy of melanoma and pancreatic cancer, the CDA-M&P formulation was RA(ED25)-TPA(2xED25)-5P-1(ED25)-PP(2xRI0.5)-sodium tannate(RI0.5). The above CDA formulations all produced 100% NBT + on HL-60 cells.

Keywords

Cancer Prevention and Therapy, Abnormal Methylation Enzymes, Cancer Stem Cells, Synthetic CDA Formulations, Differentiation Inducers, Differentiation Helper Inducers.

Introduction

Methylation enzymes play a critical role on the regulation of cell replication and differentiation, because DNA methylation controls the expression of tissue specific genes 1] and pre-rRNA ribose methylation controls the production of ribosomes 2], which in turn dictate the commitment of cells to initiate replication 3]. If enhanced production of ribosomes is locked in place, it becomes a factor to drive carcinogenesis [4]. Biological methylation is mediated by a ternary enzyme complex consisting of Methionine AdenosylTransferase (MAT)-MethylTransferase (MT)-S-Adenosyl Homocysteine Hydrolase (SAHH) [5,6]. These enzymes must be in the ternary enzyme complex to become stable and functional. In the monomeric state, individual enzymes are quickly inactivated. SAHH is the most unstable enzyme, followed by MT, and then MAT. MTs in the monomeric state have a great tendency to be converted into nucleases to trigger apoptosis. SAHH requires a steroid factor to assume a configuration favorable for the formation of dimeric enzyme complex with MT, which is then in a position to form ternary enzyme complex with MAT. In steroid hormone target tissues such as prostate and breast, steroid hormones are the stabilizing factors of SAHH. Other tissues require similar steroid factors generated by the growth signals to stabilize SAHH [7]. In normal cells, steroid factors are the dominant factors to regulate methylation enzymes. In cancer cells and telomerase expressing primitive stem cells such as embryonic stem cells and progenitor stem cells, MAT is associated with telomerase [8], which becomes the dominant factor to regulate methylation enzymes. The association of MATL, the low Km normal isozyme of MAT, with telomerase changes the kinetic properties of MATL and the regulation of methylation enzymes. Km values of MATL and MATLT, the telomerase associated tumor isozyme, are 3 µM and 20 µM methionine, respectively, and those of SAHHL and SAHHLT are 0.3 µM and 2 µM adenosine, respectively [5, 6, 8]. The increased Km value of MATLT suggests that methylation enzymes of cancer cells have elevated levels of bound S-adenosylmethionine (AdoMet). According to Prudova et al. [9], the binding of AdoMet to a protein could protect that protein against protease digestion. It appears then that the increased pool size of AdoMet in cancer cells is very important for the stability, and therefore the activity of methylation enzymes to promote malignant growth. Chiba et al. [10] found that the pool sizes of AdoMet and S-adenosylhomocysteine (AdoHcy) shrunk greatly when cancer cells were induced to undergo TD. This finding strongly supports our arguments that the association of telomerase with methylation enzymes greatly increase the stability and the activity of methylation enzymes of cancer cells so that hypomethylation of nucleic acids required for the cell to undergo TD cannot take place [6,11]. Thus, it is very convincing that abnormal methylation enzymes play a critical role on the evolution and the progression of cancer.

Association of telomerase with methylation enzymes locks methylation enzymes in extremely stable and active state to block cell differentiation. Telomerase is actively expressed in embryonic stem cells and progenitor stem cells. Differentiation of telomerase expressing normal stem cells will be blocked like cancer cells. There is another way to achieve DNA demethylation bypassing the differentiation blockade created by abnormal methylation enzymes. Tet dioxygenases carry out oxidation of 5 mC to generate 5 hmC, 5 fC and 5 caC [12–15], and 5 caC is finally replaced with C by thymine DNA glycosylase [16,17]. 5 hmC is the stable intermediate in the oxidative demethylation of 5 mC [18]. So far three Tet dioxygenases have been identified. Tet 1 preferentially acts on the 5 mC located at the transcriptional start site, whereas Tet 2 preferentially acts on the 5 mC located in the gene body [19]. Tet 3 is expressed at very high levels in oocytes and zygotes, with rapidly declining at the two cell stage. Tet 3 is responsible for the erase of paternal 5 mC in fertilized oocytes [20]. Tet enzymes are very active in embryonic stem cells to direct extraembryonic lineage differentiation [21,22]. These enzymes are frequently mutated to become dysfunctional or silenced in cancer cells [23–26]. The expression and function of Tet enzymes marks the difference between normal primitive stem cells and cancer cells. Consequently, destabilization of abnormal methylation enzymes is the only option workable to induce TD of cancer cells.

Destabilization of abnormal methylation enzymes is a very effective strategy for cancer therapy. The therapy of acute promyelocytic leukemia (APL) with RA yielded a stunning complete remission around 90% [27]. The remission, however, was only transient. Most patients relapsed within a year, and became resistant to further treatment [28, 29]. A combination of RA and As2O3 produced a more satisfactory long lasting remission [29]. RA is a DI, and As2O3 is a DHI [7]. DI is the chemical capable of eliminating the association of telomerase from abnormal methylation enzymes, and DHI is the inhibitor of individual enzymes of ternary methylation enzymes. It appears then that a combination of DI and DHI is necessary to make a perfect drug for cancer therapy. DHI alone can be very effective for cancer therapy too. After all, human body is producing DIs, although cancer patients are unable to retain such valuable metabolites in the body. Excess DHI can salvage the loss of endogenous DIs. Imatinib mesylate is the standard care for chronic myeloid leukemia [30], which is another good example of effective cancer therapy by destabilization of abnormal methylation enzymes. Signal transduction inhibitors such as imatinib mesylate are excellent DHIs [7]. Phenylbutyrate was our creation of the first DHI which was only modestly active requiring mM concentrations to function,[31]. Nevertheless, it has demonstrated therapeutic efficacy on often untreatable brain tumors [32,33]. Brain compartment is protected by blood brain barrier. The loss of endogenous DIs is not as severe as other body compartments. Therefore, even modest DHIs can exercise good therapeutic effects. The therapeutic effect of phenylbutyrate was greatly enhanced when it was used in combination with signal transduction inhibitors [34,35], which were much better DHIs effective in µM concentrations [7]. DHIs, such as dietary polyphenols, are frequently suggested for chemoprevention of cancer [36–39].

In 1987, Liau et al. [40] brought up chemosurveillance as a natural defense mechanism against cancer. This hypothesis was based on the observation that healthy people could maintain a steady level of hydrophobic metabolites in their plasma, whereas cancer patients tended to show deficiency of such metabolites due to excessive urinary excretion [41]. Among such metabolites are chemicals capable of inducing cancer cells to undergo TD [42,43]. The evolution of cancer in the case of MDS strongly supports the validity of this hypothesis.

MDS often starts with a display of immunological disorders associated with inflammation [44], which prompts the production of inflammatory cytokines. Among such cytokines, TNF is a critical factor related to the development of MDS [45]. It causes excessive apoptosis of bone marrow stem cells, thus severely affects the ability of the patient to produce hematopoietic cells such as erythrocytes, platelets, and neutrophils. TNF is also named cachectin, because of its involvement in the symptom known as cachexia. Cachexia is a symptom commonly shared by inflammatory patients and cancer patients. A characteristic disorder of cachexia is the excessive urinary excretion of low molecular weight metabolites because of vascular hyperpermeability cause by TNF [46,47]. As a consequence, chemosurveillance normally operating in healthy people to keep a check on progenitor stem cells is disrupted under pathophysiological conditions created by TNF to allow progenitor stem cells to buildup in order to replenish unipotent stem cells wiped out by TNF. The high levels of telomerase in the peripheral and bone marrow leukocytes in MDS patients is an indication of the widespread multiplication of progenitor stem cells which express telomerase [48,49]. During the course of MDS progression, mutations on Tet2, DNMT3A, IDH1/2, ASXL1, EZH2, and RNA splicing enzymes are frequently observed [50–54], which may play a significant role on the evolution of progenitor stem cells to become CSCs [55]. As anemia symptom becomes worse, chromosomal abnormalities such as translocation and deletion characteristic of cancer cells set in to speed up replication eventually pushing MDS patients to become acute myeloid leukemia patients [56–59].

Vidaza and decitabine are the two hypomethylating agents approved for the therapy of MDS in the USA. CDA-2 is a hypomethylating agent approved by China for the therapy of cancer,[60]. Vidaza and decitabine achieve DNA hypomethylation by promoting covalent bond formation between DNA Methyltransferase (DNMT) and the azacytosine base incorporated into DNA to titrate out DNMT [61], whereas CDA-2 achieves DNA hypomethylation by converting abnormal methylation enzymes into normal enzymes [6,43]. An aborted clinical trial of CDA-2 on MDS was conducted on 117 patients in China. Based on two cycles of treatment protocols, CDA-2 yielded a slightly better therapeutic efficacy under cytological evaluation, and a marked better therapeutic efficacy under hematological improvement evaluation in comparison to vidaza and decitabine,[62, 63]. Apparently CDA-2 had a better therapeutic effect and devoid of serious adverse side effects, whereas decitabine was a proven carcinogen [64]. Since MDS is a disease attributable to CSCs [55], synthetic CDA formulations ought to do well on CSCs.

CSCs constitute only a small subpopulation within a tumor. These cells, nevertheless, are now thought to confer many of adverse characteristics that contribute to treatment failure,[65–69]. Many biological characteristics that enable cancer progression are attributable to CSCs, including angiogenesis, metastasis, recurrence, and drug resistance. Elimination of CSCs is, therefore, very critical to the success of cancer therapy. CSCs are both resistant to cytotoxic chemotherapy and radiotherapy, because these cells overexpress ATP binding cassette drug pumps, and are mostly in dormant state unresponsive to radiation [70–73]. CSCs are equivalent to progenitor stem cells of normal organs or tissues, which replicate only in response to developmental or pathological needs, e.g. growth or wound healing. Thus, CSCs are very responsive to induction of differentiation, which may be the most effective approach to target CSCs.

The objective of this study is to use CDA-2 as a model to develop synthetic CDA formulations for the prevention and therapy of cancer via targeting of CSCs. CDA-2 is a preparation of natural hydrophobic metabolites purified from freshly collected urine by reverse phase chromatograph [74]. We have carried out extensive studies on DHIs of CDA-2 [7,31,74–77]. We are now focusing on DIs of CDA-2, which are the most important active components of CDA-2. When DIs becomes available, we will be in a position to solve problems brought up by CSCs.

Methods and Materials

Chemicals and Reagents

Chemicals, chromatographic supplies, and cell culture supplies were purchased from Sigma, St. Louise, MO, unless otherwise indicated. 35×10 mm cell culture dishes were purchase from CytoOne, USA Scientific. Com. Sep-Pak C18 cartridges were purchased from Walters Associates, Milford, MA. Cosmocil C18 column was purchased from Nacalai Tesque, Kyoto, Japan. Acidic pentapeptides were purchased from GenicBio Company of Shanghai, China. 1.5 liter of CDA-2 solution, 304 mg/ml, was a gift of Mr. Zhanji Sun, the general manager of NT Pharmaceuticals, Jiangsu Co. Ltd, China.

Culture of HL-60 Cells

HL-60 cells were purchased from ATCC, Manassas, VI, which were initially maintained in ISCOVE’s modified medium, supplemented with 10% fetal bovine serum, 2 mM glutamine, 50 units/ml penicillin-50 µg/ml streptomycin for a few generations, and then transferred to RPMI 1640 medium to replace ISCOVE’s modified medium. Cells were subcultured every 3 to 4 days at an initial concentration of 5–10 × 104 cells/ml.

NBT assay

NBT assay was conducted as previously described [7]. Each 35×10 mm cell culture dish contained 2 ml of RPMI 1640 culture medium. HL-60 cells at an initial concentration of 5–10 × 104 cells were incubated with or without drugs for 3 days. Approximately 2.5×105 cells were precipitated at 600xg for 5 min. The cell pellet was suspended in 3 drops from a Pasteur pipet of NBT reagent consisting 1 mg NBT and 5 µg TPA per ml Hank balanced salt solution (HBSS), and incubated at 37o C for 30 min. The reaction was terminated by the addition of a drop from a Pasteur pipet of 4% paraformaldehyde in HBSS. NBT + cells were counted under microscope using a hemacytometer.

Determination of potency of DHIs

The potency of DHIs was assessed by the reductive index as previously described [7]. Cell culture dishes were divided into several sets of 5 dishes containing RA of different concentrations to induce between 0 to 60% NBT + . One set had RA alone as control to yield Effective Dosage50 (ED50) of RA. Other sets had different concentrations of DHIs together with RA concentrations matching the control set. After incubation at 37o C for 72 h, cell numbers from each dish were counted, and an aliquot was withdrawn for NBT assay as above described. NBT + cells in the control dishes without any drug were always below 1%. In the presence of different DHIs alone, NBT + cells in general were below 10%. The respective control value was subtracted from each experimental value to yield the actual ED value. ED50 value, defined as the dosages that induced 50% NBT + cells, were estimated from plots of NBT + values versus concentrations of RA in the absence and presence of DHIs. The Reductive Index (RI) is defined as the ED50 in the presence of DHI divided by the ED50 value of RA alone. The value is inversely related the effectiveness of the DHI agent.

Bio-Gel P2 gel filtration

DIs preparation in less than 4 ml was put onto a Bio-Gel P2 column, 2.5×95 cm for gel filtration resolution. The elution was carried out by 25 mM phosphate buffer, pH 7.8, collecting 4.2 ml/tube/5 min. An aliquot from each tube was withdrawn to dilute with 1 ml of H2O for the determination of A280 absorption, and another aliquot was withdrawn from the filtrate of 0.2 µm membrane filter for the determination of NBT + .

TLC Chromatography

DIs of different Kav fractions from Bio-Gel P2 gel filtration were recovered by C18 cartridge. The fraction was acidified to pH 2.5 and passed the solution through a C18 cartridge. The cartridge was washed twice with 5 ml H2O, and then eluted with 3 ml of 80% methanol. The methanol eluant was evaporated to dryness in a rotary evaporator. The residue was dissolved in a small amount of methanol to apply to a plate of silica gel. The chromatography was developed by ascending chromatography with BuOH-HoAc-H2O (9:2:4) for 8 h to allow the solvent to travel for 17.5 cm from the origin. The silica plate was air dried in a hood overnight. UV images were marked by a pencil. The bands were scraped off the plate with a spatula to put into centrifuge tubes for the extraction with methanol. The methanol extract was evaporated to dryness in a rotary evaporator, and the residue was dissolved in a small amount of methanol for the determination of A280 and NBT + .

Sephadex LH20 Chromatography

A methanol solution of DI subfraction in less than 4 ml was put onto a Sephadex LH20 column, 2.5×36 cm, for chromatography. The elution was carried out by methanol, collecting 2.8 ml/tube/2 min. An aliquot from each tube was withdrawn to dilute with 1 ml of methanol for the determination of A280 absorption, and another aliquot was withdrawn for the determination of NBT + .

DEAE-Sephadex Chromatography

An aqueous solution of DI subfraction with NaCl concentration less than 25 mM and pH 7.8 was put onto a DEAE-Sephadex column, 1.4×27 cm, for chromatography. The column was initially washed with H2O until A280 absorption was no longer detectable. The column was then eluted with 160 ml of a linear gradient of NaCl from 0 to 2 M, collecting 4 ml/tube/4 min. An aliquot from each tube was withdrawn to dilute with l ml of H2O for the determination of A280 absorption, and another aliquot was withdrawn from the filtrate of 0.2 µm membrane filter for the determination of NBT + .

HPLC

A 50 µl of active DI preparation was injected into a Cosmosil (5C18-RA-II) column, 4.6×250 mm, for HPLC resolution using Hewlett Packard 1050 instrument. The flow rate was set at 0.5 ml/min. The elution during the initial 10 min was carried out by a linear gradient of solution A (5% HoAc) from 100% to 0%, and solution B (80% methanol) from 0% to 100%. Thereafter, the column was eluted with solution B until no more A280 absorption was detectable. HPLC fractions were evaporated to dryness in a rotary evaporator. The residue was redissolved in a small volume of methanol for the determination of NBT + .

Mass Spectroscopy

HPLC purified DHI from CDA-2 was injected into LC-MS instrument for the determination of mass.

Results

Purification of DIs from CDA-2 solution

Obviously active DIs and DHIs of CDA-2 are making a significant contribution to protect the vast majority of healthy people from becoming cancer patients. We ought to study these active metabolites which are doing a big favor to benefit human being.

When CDA-2 solution, 304mg/ml, pH 6.7, was acidified to pH 2 with 2N HCl, 32% of A280 absorption and almost 100% of DIs were found in the sediment collected by centrifugation at 1200xg for 30 min. The sediment was 50 ml of very dark colored viscous liquid from 1.5 liter of CDA-2 solution. The sediment was first extracted 3 times with 2 volumes of dichloromethane in a flask by shaking overnight each time in a shaker. Dichloromethane extract was the orange solution of the upper phase, which was poured off the flask. Dichloromethane extracted 16.1% of the A280 absorption and 37.6% of the DI activity from the pH 2 sediment of CDA-2 solution. Next, the remaining sediment was dissolved in 5 volumes of methanol, and insoluble materials were removed by centrifugation at 1200xg for 30 min. Methanol extracted 69.7% of the A280 absorption and 60% of the DI activity from the pH 2 sediment of CDA-2 solution. Methanol insoluble residue was dissolved in dilute NaOH solution, and pH adjusted to 7.8, which constituted 14.1% of the A280 absorption and 2.4% of the DI activity of the pH 2 sediment of CDA-2 solution. DIs soluble in dichloromethane are very hydrophobic metabolites and DIs soluble in methanol are less hydrophobic.

Gel filtration profiles of DIs on Bio-Gel P2 column chromatography

Organic solvents of DI extracts were removed by rotary evaporator at temperature below 60o C. The residues were dissolved in dilute NaOH solution and pH adjusted to 7.8. 4 ml aliquot from each preparation containing approximate 100,000 A280 absorption units was put on a column of Bio-Gel P2 column for gel filtration as described in Methods and Materials. Dichloromethane extractable DIs have two major peaks at Kav = 0.43 and 0.52, and a minor peak at Kav = 1.45 as shown in Fig. 1. Kav = 0.43 and 0.52 are quite large molecular weight complexes.

CST 2019-101 - Mingliau China_F1

Figure 1. Gel filtration profile of the dichloromethane extractable DIs of CDA-2

4 ml of dichloromethane extractable DIs of CDA-2 containing approximately 100,000 A280 absorption units were put onto a column of Bio-Gel P2 column for gel filtration as described in Methods and Materials.

Gel filtration profile of methanol extractable DIs is presented in Fig. 2. Methanol extractable DIs have a major peak at Kav = 0 and two minor peaks at Kav = 0.43 and 0.52 which are probably the major components of dichloromethane extractable DIs. DIs of Kav = 0 are even bigger than those of dichloromethane extractable DIs

CST 2019-101 - Mingliau China_F2

Figure 2. Gel filtration profile of the methanol extractable DIs of CDA-2

4 ml of methanol extractable DIs of CDA-2 containing approximately 100,000 A280 absorption units were put onto a column of Bio-Gel P2 column for gel filtration as described in Methods and Materials.

The dark solution of the Kav = 0 of methanol extractable DIs produced scale-like dark precipitate and clear light color solution when the pH was acidified to 2 with 2N HCl. DI activities became dissociated from the dark inactive carrier when the precipitate was dissolved in 0.5N NaOH and incubated at 37o C as shown in Fig. 3.

CST 2019-101 - Mingliau China_F3

Figure 3. Gel filtration profiles of the DIs of Kav = 0 of Fig. 2 upon incubation with 0.5N NaOH at 37o C

DIs of Kav = 0 were obtained from Fig. 2 by adjusting the active fractions to pH 2 with 2N HCl. The active DIs became scale like dark precipitate, which was collected by centrifugation at 1200xg for 30 min. The residue was dissolved in 6 ml of 0.5N NaOH, 2 ml aliquot containing a total of 1540 A280 absorption units was withdrawn to immediately adjusted to pH 7 for gel filtration as described in Fig. l. Another 2 ml aliquot was withdrawn after incubation at 37o C or 0.5 h to be neutralized for gel filtration. And the last 2 ml was neutralized for gel filtration after incubation at 37o C for 24 h.

The major DI peak shifted from Kav = 0 to Kav = 0,09 after incubation at 37oC for 0.5 h, and then to Kav = 0.43 and 0.52 after incubation for 24 h. The Kav = 0.43 and 0.52 DIs were extractable by dichloromethane when the solution was acidified to pH 2. It appears then that DIs of Kav = 0.43 and 0.52 are the major DIs of CDA-2. A good proportion of these DIs becomes covalently linked to the inactive carrier of Kav = 0 previously identified as pigment peptides,[43].

TLC chromatography of dichloromethane extractable DIs of CDA-2

Dichloromethane extractable DIs of CDA-2 from Fig. 1 were recovered by C18 cartridges and resolved by TLC as described in Methods and Materials. Result is presented in Table 1. Significant DI activities were widespread detectable from Rf = 0.15 to 0.92. Similar widespread activities were also found on DIs of Kav = 0.52. The only difference was the relative activity of different Rf bands. DIs of Kav = 0.43 had more DIs with higher Rf values. DIs of Kav = 1.45 were much simple. Rf = 0.57 band stood out as the only major active band. TLC data serve to indicate that DIs of Kav = 0.43 and 0.52 are complexes of multiple different DIs or a few limited DIs in association with multiple different DHIs or inactive materials.

Table 1. TLC chromatography of the Kav=0.43 DIs of Fig. 1

Rf

UV images

A280/ml

% Growth

% NBT +

0.95

Bright blue fluorescent

0.05

73

5.8

0.92

Gold fluorescent

0.16

59

22.8

0.86

Dark absorption

0.23

77

7.4

0.81

Gold fluorescent

0.21

54

15.8

0.74

Gold fluorescent

0.23

68

5.3

0.62

Gold fluorescent

0.19

90

11.6

0.57

Dark absorption

0.16

81

44.0

0.44

Bright blue fluorescent

0.15

47

17.5

0.34

Multiple narrow gold fluorescent

0.10

61

17.9

0.15

Weak gold fluorescent

0.05

73

11.5

0.08

Very weak gold fluorescent

0.05

85

4.3

Methanol solution of Kav=0.43 DIs of Fig. 1 was applied to a TLC plate for chromatography as described in Methods and Materials.

Sephadex LH20 Chromatography

DIs of Kav = o.43 and 0.52 from Fig. 1 were recovered by C18 cartridges as above described. 4 ml of methanol solution containing approximately 25,000 A280 absorption units were put on a column of Sephadex LH20 column for chromatography as described in Methods and Materials. The chromatographic profile is shown in Fig. 4. Two active peaks similar to two Bio-Gel P2 peaks of Kav = 0.43 and 0.52 were found, one eluted between 150–180 ml, and the other between 180–210 ml. The separation of two peaks on Sephadex LH 20 chromatography appeared more clean cut than those on Bio-Gel P2 gel filtration. The active fractions of two peaks were pooled separately, and methanol was removed by rotary evaporator. The residue of each active fractions was suspended in 5 ml H2O, and 1N NaOH was added dropwise to bring pH up to 7.8 to dissolve the residue. NaOH soluble material gave rise to a total of 11, 500 A280 absorption units from 150–180 ml peak, and 5,500 A280 absorption units from 180–210 ml peak. These preparations were separately put on a column of DEAE-Sephadex for chromatography as described in Methods and Materials. DEAE-Sephadex chromatographic profiles are shown in Fig. 5 and 6. The 150–180 ml peak of Sephadex LH20 chromatography gave rise to two active fractions eluted between 0.55–0.75 M NaCl as a major peak and between 0.85–1.05 M NaCl as a minor shoulder peak. The 180–210 ml peak of Sephadex LH 20 chromatography produced a reverse profile: a minor shoulder peak between 0.55–0.75 M NaCl and a major peak between 0.85–1.05 M NaCl. When DEAE-Sephadex chromatographic fractions were stored at 4o C for more than 1 week, noticeable precipitate was found in active fractions from 0.55–1.05 M NaCl. The precipitate was inactive as DIs. The DI activity of the supernatant declined greatly to 22% of the original activity, which was restored to almost the original activity when enough precipitate redissolved in methanol was added back. Thus, the precipitate fits the description as DHI [31]. Upon purification by HPLC, the precipitate yielded a major UV peak at RT-19,1. The eluant of the RT-19.1 peak was evaporated to dryness by rotary evaporator, and redissolved in a small volume of methanol. Subsequent HPLC analysis showed a single RT-19.1 peak as shown in Fig. 7. The mass determination of the RT-19.l HPLC peak of Fig. 7 yielded a mass of 316.118 as shown in Fig. 8. This mass is very close to the mass of pregnenolone, which is 316.48.

CST 2019-101 - Mingliau China_F4

Figure 4. Sephadex LH20 chromatography of Kav=0.43 + 0.52 pooled material of Fig. 1

Kav=0.43 + 0.52 of Fig. 1 were recovered by C18 cartridges as described in Methods and Materials. The residue was dissolved in methanol. 4 ml of methanol solution containing approximately 25,000 A280 absorption units were put on a column of Sephadex LH20 for chromatography as described In Methods and Materials. Two fractions eluted between 150–180 ml and 180–210 ml were separately pooled, and methanol was removed by rotary evaporation under vacuum.

CST 2019-101 - Mingliau China_F5

Figure 5. DEAE-Sephadex chromatography of the 150–180 ml peak of Sephadex LH20 chromatography

The residue of the 150–180 ml peak of Fig. 4 was dissolved in dilute NaOH, and pH adjusted to 7.8. 5 ml of this solution containing a total of 11,500 A280 absorption units were put onto a column of DEAE-Sephadex for chromatography as described in Methods and Materials. Active fractions eluted between 0.55–0.75 M NaCl and 0.85–1.05 M NaCl were separately pooled.

CST 2019-101 - Mingliau China_F6

Figure 6. DEAE-Sephadex chromatography of the 180–210 ml peak of Sephadex LH20 chromatography

The residue of the 180–210 ml peak of Fig. 4 was dissolved in dilute NaOH, and pH adjusted to 7.8. 5 ml of this solution containing a total of 5,500 A280 absorption units were put onto a column of DEAE-Sephadex for chromatography as described in Methods and Materials. Active fractions eluted between 0.55–0.75 M NaCl and 0.85–1.05 M NaCl were separately pooled.

CST 2019-101 - Mingliau China_F7

Figure 7. HPLC purification of the DHI dissociated from DI of CDA-2

The HPLC RT-19.1 peak was a major UV absorption peak of the precipitate from the eluant of 0.55–0.75 M NaCl of DEAE-Sephadex chromatography shown in Fig. 5. On final analysis, it yielded a single peak at RT-19.1. HPLC was carried out as described in Methods and Materials.

CST 2019-101 - Mingliau China_F8

Figure 8. Mass spectroscopy of the RT-19.1 peak of HPLC

An aliquot of the methanol solution of the material purified to a single HPLC peak shown in Fig. 7 was injected into LC/MS instrument for the mass determination as described in Methods and Materials.

Confirmation of the HPLC RT-19.1 peak of CDA-2 as pregnenolone

Pregnenolone obtained from Sigma was dissolved in methanol to make 1.5 mg/ml, namely 4.74 mM. The A280 absorption is very low, which is 0.054/mM. The solution gave a HPLC peak at RT-19.1. When RT-19.1 peak of CDA-2 was mixed with Sigma pregnenolone, the mixture gave a single peak at RT-19.1. The potency of the HPLC RT-19.1 of CDA-2 as DHI was exactly the same as that of Sigma pregnenolone as previously reported, [7], namely 7.16 µM. It is obvious that pregnenolone is a component of active liposomal DIs of CDA-2.

Acidic pentapeptides of hemoglobin as DIs

We have tried to capture DIs in the supernatant of DEAE-Sephadex active eluants. The significant A280 absorption peaks of HPLC were all inactive as DI. We thus concluded that DIs of CDA-2 might not have A280 absorption. We have previously found that acidic peptides of the urine were active Dis [42]. We had reason to believe that plasma and urinary peptides were primarily the degradative products of endogenous proteins, because peptide profiles among different persons were similar, and those profiles were most close to the peptide profile of the spleen extract [40,41]. It is well known that erythrocytes are constantly turning over, and spleen is the organ to degrade dead erythrocytes. We assumed that the active acidic peptide of CDA-2 might be the degradative products of hemoglobin. Taking cue from the discovery of an acidic pentapeptide to modulate hemopoiesis by Laerum and Paukovits [78], we randomly picked pentapeptides containing at least two acidic amino acid residues from the sequences of α- and β-hemoglobin for synthesis to test their activities as DIs. Indeed, acidic pentapeptides of hemoglobin were active as DIs as shown in Table 2, although the activities were only modest.

Table 2. DI activity of synthetic acidic pentapeptides

Designation

Sequence

ED25, µM

5P-1

Ala-Glu-Ala-Leu-Glu

122±32

5P-2

Val-Asp-Glu-Val-Gly

154±28

5P-3

Val-Asp-Asp-Met-Pro

174±55

5P-4

Asp-Pro-Glu-Asp-Phe

207±42

PA-5P-4

PA-Asp-Pro-Glu-Asp-Phe

161±33

5P-5

Pro-Glu-Glu-Lys-Ser

415±76

Synthetic acidic pentapeptides were dissolved in 80% methanol for the test of DI activity. The dosages were selected to give ED between 20 and 30% NBT + . ED25 was obtained from the plots, and expressed as mean±S. D. from at least two experiments.

The combination of RA and acidic pentapeptides produced additive effect, whereas the combination of RA and DHIs produced synergistic effect as previously reported [7,31,76,77].

RA is a well known DI, which requires retinoic acid receptor {RAR} to activate oligoisoadenylate synthetase [79], and the product of this enzyme, oligoisoadenylate, is the responsible DI to initiate differentiation induction [6]. Oligoisoadenylate has to be synthesized inside the cell to function, because 5’ terminal triphosphate prohibits it to be absorbed from outside. Therefore, RA is only effective on cells expressing RAR. RARs are expressed on developmental stages of embryonic cells [80], in pluripotent stem cells to regulate organ and limb development [81], in CSCs [82–84], and in certain cancers such as APL [27], acute myeloid leukemia [85], neuroblastoma [86], breast cancer [87], and melanoma [88].

TPA is another well known DI, which acts as DI in some cells, but as differentiation inhibitor in other cells depending on which molecules are being affected, or as tumor promoter by increasing cell multiplication in initiated cells [89]. Biological effects of TPA are mediated through initial interaction with membrane [90]. Therefore, TPA affects almost all cancers. The tumor promotion activity of TPA can be effectively cancelled by DHIs [91]. Despite its adverse effects, TPA has been put through clinical trials for the evaluation as a cancer drug [92–94]. TPA had a very impressive DI activity on HL-60 cells as shown in Fig. 9. The maximum activity was at 0.4 nM to induce 92% NBT + . At this concentration 90% of cells were attached to the culture dish. Above this concentration, cell growth was greatly inhibited, and NBT + cells declined precipitously. ED25 of TPA was 0.17 ± 0.02 nM when tested on fast growing HL-60 cells. The growth rates measured by N3/N0 were between 13.7–15.4, almost 3 to 4 generations in 3 days. ED25 of RA was 1.17 ± 0.16 µM tested on HL-60 cells with growth rates same as those of TPA. We have previously found that HL-60 cells responded better to RA when grew slower at earlier passages. ED25 of RA was 0.78 µM when tested on HL-60 cells with N3/N0 between 6.3–7.5.

CST 2019-101 - Mingliau China_F9

Figure 9. DI activity of TPA

TPA was dissolved in methanol which was added to the culture in the amounts indicated. The volume of methanol was less than 2%. NBT assay was conducted as described in Methods and Materials.

Potential DHIs for synthetic CDA formulations

We have previously shown that RA alone could not push all HL-60 cells to complete TD [7]. 95% NBT + cells was the maximum extent achieved by RA. But in the presence of DHIs, the extent of NBT + cells could reach 100% [7]. 92% NBT + cells was the maximum extent reached by TPA in this study. Cells unable to undergo differentiation in the presence of DI alone are the cells damaged enough to prevent completion of two replication cycles needed for TD. The damaged cells if repaired become the roots for recurrence. Removal of the possible recurrent roots by DHIs is another important contribution of DHIs. The most serious damage that can be attributed to DIs are the promotion of the dissociation of methylation enzymes, and the conversion of the monomeric MTs into nucleases. MT inhibitors can effectively prevent such conversion. Uroerythrin was a good MT inhibitor to function as the most active DHI of CDA-2 [77]. This chemical is not commercially available. We need to find substitutes to prevent damages created by DIs alone. We have previously found that I3 and I4 from limited alkaline hydrolysis of poly I were potent inhibitors of nucleolar rRNA methyltransferase [2]. We, thus, prepared I3 + I4 and A3 + A4 as previously described [2] to test their DHI activities. Results shown in Table 3 indicate that I3 + I4 and A3 + A4 are active as DHIs. Included in Table 3 is PP as the most active DHI we have discovered so far. PP has been shown active as signal transduction inhibitor and inhibitor of androgen receptor,[95–98]. Both classes of inhibitors are excellent DHIs [7]. The superb activity of PP as DHI may be attributable to its dual activities to inhibit signal transduction and steroid receptor.

Table 3. Potential DHIs for synthetic CDA formulations

DHI

RI0.5

I3+I4

3.7±0.58 µM

A3+A4

4.8±1.12 µM

PP

12.2±2.2 nM

10 mg of poly I and poly A was dissolved in 1 ml of 0.3N KOH and incubated at 37o C for 0.5 h. The solution was neutralized with HCl, and put onto a column of DEAE-Sephadex for chromatography as previously described,[2]. I3 and I4 in the ratio of 1:0.9 were the major components, which were combined and recovered as previously described,[2] for the test of DHI activity. The preparations of oligonucleotides were dissolved in 80% methanol, and pp was dissolved in methanol for the determination of RI0.5 as previously described,[7]. HL-60 cells for these experiments had growth rates of N3/N0 between 7.5–9.8.

Development of synthetic CDA formulations for the prevention and therapy of cancer via targeting of CSCs

CDA-2 was effective to eliminate CSCs subpopulation [99,100], and had a very effective therapeutic efficacy on MDS [62, 63], a disease attributable to CSCs. Our immediate objective is to develop synthetic CDA formulations for the therapy of MDS, CSCs, and cancers in which CSCs play a dominant role on the pathological features and drug responses of the diseases. Cancers originated from progenitor stem cells such as MDS above described and acute myeloid leukemia resulting from MDS are enriched in CSCs. Evidence indicates that brain tumors are also originated from progenitor stem cells [101,102]. Melanoma and pancreatic cancer may not originate from progenitor stem cells, but these cancers express high levels of factors such as hypoxia inducible factors and NF-kB that can effectively drive EMT to display CSC phenotype [103 -106]. On MDS, protection of pathological cells is a necessity, because the therapy requires the differentiation of pathological cells to become functional erythrocytes. DIs and DHIs of low toxicity are preferred. On other CDA formulations, destruction of pathological cells is a preferred choice. On CSCs, toxicity is also a big concern, because CSCs naturally reject toxic chemicals. CSCs are situated in hypoxic microenvironment hard to reach by the blood, small molecules with penetration power are preferred. On brain tumor, blood brain barrier must be the major concern. Hydrophobic chemicals have better chance to cross that barrier. For melanoma and pancreatic cancer, potent DIs and DHIs are preferred. CDA formulations must consist of both DIs and DHIs to make perfect drugs. The following formulations are acceptable formulations: 3xED25 of DI + 1xRI0.5 of DHI; 2xED25 of DI + 2xRI0.5 of DHI; or 1xED25 of DI + 3xRI0.5 of DHI. 1xRI0.5 of DHI is equivalent to 1xED25 of DI [7]. The multiplicity of DI can be multiplicity of a single DI, or a combination of different DIs. Likewise, the multiplicity of DHI can be multiplicity of a single DHI, or a combination of different DHIs.

Our deliberated CDA formulations are listed in table 4. CDA-MDS, CDA-BT (brain tumor), and CDA-M&P (melanoma and pancreatic cancer) are parenteral preparations, and CDA-CSC is an oral preparation which is designed for long term application to prevent recurrence after active therapy. Thus, an oral preparation is a practical drug form. On active cancer therapy, CDA-CSC is best used in combination with drugs to inhibit growth factor receptors and signal transductions, or drugs for anti-steroid hormonal therapy. These drugs are actually good DHIs. It is also compatible with other therapeutic modalities to target on different entities, e.g. CDA-CSC to target on CSCs and cytotoxic drugs or radiation to target on non-CSCs. TPA and PP are not soluble in H2O. RA, resveratrol, and curcumin are not very soluble in H2O. Dispensing aids such as liposomal and nanoparticle technologies must be employed to increase bioavailability for clinical application. Phenylacetylglutasmine {PAG} is a major chemical constituent of CDA-2, which is also very useful to extend bioavailability of DIs and DHIs. We recommend to include 2 mM PAG in each CDA formulation for the purpose of extending bioavailability of DIs and DHIs and to cancel tumor promotion activity of TPA. We have previously shown that PAG, namely Antineoplaston A10, was effective to prevent the loss of low molecular weight metabolites including active DIs and DHIs [40]. By preventing the loss of endogenous DIs and DHIs, PAG was effective to prevent carcinogens-induced pulmonary neoplasia and hepatoma [107,108], and to achieve therapy of early stage cancer [40]. PAG is a very effective deterrent of tumor promoters. PAG is synthesized from phenylacetyl chloride and glutamine as previously described [76]. CDA formulations listed in Table 4 are plasma concentrations to achieve induction of TD of HL-60 cells above 100% on purpose to achieve maximum therapeutic effect. Multiplication of these amounts per liter of blood by a factor of 5, which is the normal blood volume of a person of 80 Kg body weight, is necessary to produce an effective dosage. Application of 3 effective dosages a day should give satisfactory therapy.

Table 4. Synthetic CDA formulations

Effect on HL-60 cells

Designation

Formulation

% Growth

%NBT+

CDA-MDS

RA(ED25)-5P-1(ED25)-I3+I4(RI0.5)-PP(RI0.5)-Na pregnenolone SO4(RI0.5)

77 ± 5.2

100

CDA-CSC

RA(ED25)-TPA(ED25)-PP(RI0.5)-resveratrol(RI0.5)-curcumin(RI0.5)

53 ± 3.6

100

CDA-BT

TPA(2xED25)-PP(2xRI0.5)-Na phenylbutyrate (RI0.5)-pyrogallol(RI0.5)

45 ± 3.9

100

CDA-M&P

TPA(2xED25)-5P-1(ED25)-PP(2xRI0.5)-Na tannate (RI0.5)

33±4.4

100

CDA-MDS is a parenteral preparation as a monotherapy of MDS. CDA-CSC is an oral preparation for use in combination with other drugs. CDA-BT is a parenteral preparation as a monotherapy of brain tumor. CDA-M&P is also a parenteral preparation as a monotherapy for melanoma and pancreatic cancer.

DIs and DHIs were methanol solution except 5P-1 which was in 80% methanol, and Na phenylbutyrate which was in H2O solution. Each component was added individually. RA(ED25) is 1.17 µM, 5P-1(ED25) is 122 µM, I3+I4(RI0.5) is 3.7 µM, PP(RI0.5) is 12.2 nM, TPA(ED25) is 0.17 nM, resveratrol (RI0.5) is 1.16 µM, curcumin(RI0.5) is 1.24 µM, Na phenylbutyrate(RI0.5) is 2 mM, and Na tannate(RI0.5) is 0.37 µM as previously reported,[7] or shown in this paper. NBT assay was carried out as described in Methods and Materials.

Discussion

Abnormal methylation enzymes are an important issue of cancer. Almost all human cancers display such an abnormality [8,109]. Abnormal methylation enzymes are the critical factor responsible for the blockade of differentiation of cancer cells [6]. The expression of telomerase fits the first hit of Knudson’s two hits theory [110]. Telomerase is also expressed in primitive stem cells such as embryonic stem cells and progenitor stem cells. Evidently blockade of differentiation is a normal biological process to build up cell mass for the development of fetus, or wound healing. The blockade of differentiation in these normal stem cells does not create problem, because Tet enzymes are functioning which can break through the blockade created by abnormal methylation enzymes to spearhead differentiation programs. The problem arises if Tet enzymes become dysfunction due to mutation or silencing such as the case of MDS [23–26]. The loss of Tet enzymes becomes the second hit of Knudson’s theory. Therefore, the evolution of cancer in the case of MDS is a perfect interpretation on the validity of the two hits theory of Knudson [110], and our chemosurveillance hypothesis [40].

Progenitor stem cells and CSCs are almost indistinguishable with only very minor differences, the functionality of Tet enzymes being a significant difference between them. Progenitor stem cells are very likely the origin of most human cancers, not just acute myeloid leukemia resulting from MDS or brain tumors. We have detected abnormal methylation enzymes to function actively in preneoplastic hyperplastic nodules induced by hepatocarcinogen [111]. Therefore, progenitor stem cells are just like time bombs, and natural DIs and DHIs are the safety devices to keep them from evolving into CSCs.

Chemicals extractable by dichloromethane are in general very hydrophobic. The detection of pregnenolone as a component of DIs of CDA-2 is a certain indication that active DIs of CDA-2 is liposomal complexes. The release of liposomal-like DIs from methanol extractable DIs by NaOH hydrolysis suggests that dichloromethane extractable DIs and methanol extractable DIs are basically the same liposomal complexes. A large proportion of these liposomal complexes become covalently linked to the inactive carrier previously identified as pigment peptides [40] either by enzymatic dehydration in the body or due to heat pyrolysis during evaporation of ethanol eluant in the process of CDA-2 preparation. The inactive pigment peptide carrier is very likely fragments of membrane. The liposomal DIs of Kav = 0.43 and 0.52 are stable in methanol, ethanol, and dichloromethane, but are not stable in butanol or high slat solution. DIs in free forms do not seem to have A280 absorption. Without A280 absorption as a probe, it is very difficult to capture DIs of CDA-2. The real DIs of CDA-2 are still up in the air.

We have identified pregnenolone as a DHI of CDA-2 in this study. Pregnenolone is the master substrate of all biologically active steroids, which is synthesized from cholesterol absorbed from diet or produced in the liver. The peak age of the production of pregnenolone is 20 years old, producing approximately 50 mg a day,[112]. The very young and the very old are the two age groups producing the least amount of pregnenolone. These are the two age groups most susceptible to cancer. Fortifying effective DHIs definitely is a good policy to amend natural insufficiency.

The therapy of APL with RA sets a classical example of the effectiveness of destabilization of abnormal methylation enzymes on cancer therapy [27]. RA alone is very effective, but imperfect. Recurrence ensues quickly [28]. This is because RA alone causes the dissociation of abnormal methylation enzymes too extensive to contribute to nuclease activity which creates damages to interrupt differentiation process. This phenomenon is analogous to the antiviral effect of interferon. Interferon is actually a DI like TPA to activate oligoisoadenylate synthetase through interaction with membrane. The product of oligoisoadenylate synthetase is a potent activator of latent nucleases to execute antiviral effect. The damage created by DI alone may affect only a small fraction of cells in the S phase when DNA is most susceptible to nuclease attack. The damages interrupt differentiation process, but when damaged cells are repaired later they become the roots for recurrence. Damages by DI alone can be prevented by DHIs. MT inhibitors can prevent protease modification to turn MT into nuclease. Steroids analogs are even better, which keep MT in dimeric complex with SAHH to protect the integrity of MT. The protection of MT from becoming nuclease is the cause that DHI is an essential partner of a perfect cancer drug.

The option for the eradication of cancer stem cells is very limited. The ideal therapeutic agents must be small molecules that are relatively non-toxic to bypass drug afflux pumps to reach adequate intracellular concentrations to trigger cancer stem cells to undergo differentiation. Synthetic CDA formulations are an attractive strategy to eradicate cancer stem cells. DIs is the most important components of CDA formulations. The choices of DIs are very limited; RA is only effective in cancers expressing RAR. TPA is universal, but is very toxic. 5Ps are universal, but are not very active. The requirement of DI is to initiate TD to reach a level above 15%, the rest can be accomplished by DHIs [7]. Therefore, even though 5Ps are not very active, they can be remedied by very active DHIs such as PP. We have developed many excellent DHIs to choose from. The mere application of effective DHIs can have excellent therapeutic effect on cancer as exemplified by imatinib mesylate. The brain compartment is a unique compartment protected by blood brain barrier, which also protect the loss of endogenous metabolites. Therefore, the deficiency of endogenous DIs and DHIs is not as severe as other body compartments. That is why DHIs have impressive therapeutic effect on brain tumor.

Destabilization of abnormal methylation enzymes is a good policy to implement at the earliest diagnosis of cancer to avoid aberrant DNA methylation no matter what cancer therapy is to follow. Aberrant DNA methylations are frequently detected in advanced cancer, which can turn ordinary responsive cancer to become very vicious cancer unresponsive to any treatment. Aberrant DNA methylations happen more readily when DNA synthesis is slowing down by cytotoxic drugs or radiation [113]. Another benefit to implement destabilization of abnormal methylation enzymes is to avoid expansion of CSC population induced by the destruction of the tumor caused by cytotoxic chemotherapy [114], or due to tumor progression [103–106]. Reconstruction of the tumor is a major mission of CSCs. Therefore, destabilization of abnormal methylation enzymes is a very good strategy to pursue for cancer therapy.

Acknowledgement

This study was supported in part by a contract awarded to CDA Therapeutics, Inc. by Xinhua Pharmaceutical Company of Zibo, Shandong China. We are indebted to the gift of CDA-2 solution from Mr. Zanji Sun of NTPharma of Tai Zhou, Jiangsu, China. We are very grateful to Christian Browder for the art work of figures.

Funding Information

The funder had no involvement in the study design; in the collection, analysis and interpretation of the data; in the writing of the report, and in the decision to submit the paper for publication. The funder had an active supervising role.

References

  1. Racanelli AC, Turner FB, Kie TY, Taylor SM (2008) A mouse gene that coordinates epigenetic controls and transcriptional interference to achieve tissue-specific expression. Mol Cell Biol 281: 836–848.
  2. Liau MC, Hunt ME, Hurlbert RB (1976) Role of ribosomal RNA methylases in the regulation of ribosome production in mammalian cells. Biochemistry 15: 3158–3164. [crossref]
  3. Bernstein KA, Bleichert F, Bean JM, Cross FR, Baserga SJ (2007) Ribosome biogenesis is sensed at the Start cell cycle checkpoint. Mol Biol Cell 18: 953–964. [crossref]
  4. Justilien V, Ali SA, Jamieson L, Yin N, Cox AD, et al. (2017) Ect2-Dependent rRNA Synthesis Is Required for KRAS-TRP53-Driven Lung Adenocarcinoma. Cancer Cell 31: 256–269. [crossref]
  5. Liau MC, Chang CF, Saunders GF, Tsai YH (1981) S-Adenosylhomocysteine hydrolases as the primary target enzymes in androgen regulation of methylation complexes. Arch Biochem Biophys 208: 261–272.
  6. Liau MC (2004) Abnormal methylation enzymes: a selective molecular target for differentiation therapy of cancer. Chin Pharm J 56: 57–67.
  7. Liau MC, Kim JH, Fruehauf JP (2019) Potentiation of ATRA activity in HL-60 cells by targeting methylation enzymes. Pharmacol Pharmaceut Pharmacovig 3: 9–17
  8. Liau MC, Zhuang P, Chiou GCP (2010) Identification of the tumor factor of abnormal methylation enzymes as the catalytic subunit of telomerase. Chin Oncol Cancer Res 7: 86–96.
  9. Prudova A, Bauman Z, Braun A, Vitvitsky V, Lu SC, et al. (2006) S-adenosylmethionine stabilizes cystathionine beta-synthase and modulates redox capacity. Proc Natl Acad Sci USA 103: 6489–6494. [crossref]
  10. Chiba P, Wallner L, Kaizer E (1988) S-Adenosylmethionine metabolism in HL-60 cells: effect of cell cycle and differentiation. Biochim Biophys Acta 971: 38–45.
  11. Liau MC, Lee SS, Burzynski SR (1989) Hypomethylation of nucleic acids: a key to the induction of terminal differentiation. Intl J Exp Clin Chemother 2: 187–199.
  12. Tahiliani M, Koh KP, Shen Y, Pastor WA, Bandukwala H, et al. (2009) Conversion of 5-methylcytosine to 5-hydroxymethylcytosine in mammalian DNA by MLL partner TET1. Science 324: 930–935. [crossref]
  13. Ito S, Shen L, Dai Q, Wu SC, Collins LB, et al. (2011) Tet proteins can convert 5-methylcytosine to 5-formylcytosine and 5-carboxylcytosine. Science 333: 1300–1303. [crossref]
  14. Shen L, Song CX, He C, Zhang Y (2014) Mechanism and function of oxidative reversal of DNA and RNA methylation. Annu Rev Biochem 83: 585–614. [crossref]
  15. Wu X (2017) TET-mediated active DNA demethylation: mechanism, function and beyond. Nat Rev Genet 18: 517–534. [crossref]
  16. He YF, Li BZ, Li Z, Liu P, Wang Y, et al. (2011) Tet-mediated formation of 5-carboxylcytosine and its excision by TDG in mammalian DNA. Science 333: 1303–1307. [crossref]
  17. Maiti A, Drahat AC (2011) Thymine DNA glycosylase can rapidly excise 5-formylcytosine and 5-carboxylcytosine: potential implication for active demethylation of CpG sites. J Biol Chem 286: 35334–35338.
  18. Hu L, Lu J, Cheng J, Rao Q (2015) Structural insight into substrate preference for TET-medicated oxidation. Nature 527: 118–122.
  19. Huang Y, Charez L, Chang X, Wang X (2014) Distinct role of the methylcytosine oxidases Tet 1 and Tel 2 in mouse embryonic stem cells. Proc Natl Acad Sci USA 111: 1361–1366.
  20. Iqbal K, Jin SG, Pfeifer GP, Szabó PE (2011) Reprogramming of the paternal genome upon fertilization involves genome-wide oxidation of 5-methylcytosine. Proc Natl Acad Sci U S A 108: 3642–3647. [crossref]
  21. Ficz G, Branco MR, Seisenberger S, Santos F (2011) Dynamic regulation of 5-hydroxymethylcytosine mouse ES cells and during differentiation. Nature 473: 398–403.
  22. Monfort A, Wutz A (2013) Breathing-in epigenetic change with vitamin C. EMBO Rep 14: 337–346. [crossref]
  23. Haffiner MC, Chaux A, Meeker AK, Esopi D (2011) Global 5-hydroxymethylcytosine content is significantly reduced in tissue stem/progenitor cell compartments and in human cancers. Oncotarget 2: 627–637.
  24. Kudo Y, Tateishi K, Yamamoto K, Yamamoto S (2012) Loss of 5-hydroxymethylcytosine is accompanied with malignant cellular transformation. Cancer Sci 103: 670–676.
  25. Ficz G, Gribben JG (2014) Loss of 5-hydroxymethylcytosine in cancer: cause or consequence? Genomics 104: 352–357. [crossref]
  26. Wang KC, Kang CH, Tsai CY, Chou NH (2018) Ten-eleven translocation 1 dysfunction reduces 5-hydroxylmethylcytosine expression levels in gastric cancer cells. Oncol Lett 15: 278–284.
  27. Huang M, Ye Y, Chen S, Chai JR (1988) Use of all-trans retinoic acid in the treatment of acute promyelocytic leukemia. Blood 72: 567–572.
  28. Warrell RP, Frankel SR, Miller MA, Scheinberg DA (1991) Differentiation therapy of acute promyelocytic leukemia with tretinoin (all-trans retinoic acid). N Engl J Med 304: 1385–1393.
  29. de Thé H, Pandolfi PP, Chen Z (2017) Acute Promyelocytic Leukemia: A Paradigm for Oncoprotein-Targeted Cure. Cancer Cell 32: 552–560. [crossref]
  30. Gambacorti-Passerini C, Antolini L, Mahon FX, Guilhot F, Deininger M, et al. (2011) Multicenter independent assessment of outcomes in chronic myeloid leukemia patients treated with imatinib. J Natl Cancer Inst 103: 553–561. [crossref]
  31. Liau MC, Liau CP, Burzynski SR (1992) Potentiation of induced terminal differentiation by phenylacetic acid and related chemicals. Intl J Exptl Clin Chemother 5: 9–17.
  32. Baker MJ, Brem S, Deniels S, Sherman B (2002) Complete response of a recurrent, multicentric malignant glioma in a patient treated with phenylbutyrate. J Neurooncol 59: 239–242.
  33. Burzynski SR, Nagy-Kubove E (2011) Treatment of esthesioneuroblastoma and non-small cell-lung cancer with phenylbutyrate. J Cancer Ther 2: 518–522.
  34. Burzynski SR, Tamicki TJ, Burzynski GS, Brookman S (2014) Preliminary findings on the use of targeted therapy with pazopanib and other agents in combination with sodium phenylbutyrate in the treatment of glioblastoma multiforme. J Cancer Ther 5: 1423–1437.
  35. Burzynski SR, Burzynski GS, Brookman S (2015) A case of sustained objective response of recurrent/progressive diffuse intrinsic pontine glioma with phenylbutyrate and targeted agents. J Cancer Ther 6: 40–44.
  36. Stoner GD, Wang LS, Zikri N, Chen T, Hecht SS, et al. (2007) Cancer prevention with freeze-dried berries and berry components. Semin Cancer Biol 17: 403–410. [crossref]
  37. Amin AR, Kucuk O, Khuri FR, Shin DM (2009) Perspectives for cancer prevention with natural compounds. J Clin Oncol 27: 2712–2725. [crossref]
  38. Wang LS, Hecht SS, Carmella SG, Yu N (2009) Anthocyannin in black raspberries prevent esophageal tumors in rats. Cancer Prev Res 2: 84–93.
  39. Allegra A, Innao V, Russo S, Gorace D (2017) Anticancer activity of curcumin and its analogues: preclinical and clinical studies. Cancer Invest 35: 1–22.
  40. Liau MC, Szopa M, Burzynski B, Burzinski SR (1987) Chemosurveillance: a novel concept of the natural defense mechanism against cancer. Drug Exptl Clin Res 13: 77–82.
  41.  Liau MC, Szopa M, Burzynski B, Burzinski SR (1987) Quantitative assay of plasma and urinary peptides as an aid for the evaluation of cancer patients undergoing Antineoplaston therapy. Drug Exptl Clin Res 13: 61–70.
  42. Liau MC, Lee SS, Burzynski SR. (1988) Differentiaiton inducing components of Antineoplaston A5. Adv Exptl Clin Chemother 6/88: 9–26.
  43. Liau MC, Lee SS, Burzynski SR (1990) Modulation of cancer methylation complex isozymes as a decisive factor in the induction of terminal differentiation mediated by Antineoplaston A5. Intl J Tiss React 12: 29–36.
  44. Williamson RJ, Kruger AR, Reynolds RJ (1994) Establishing the incidence of myelodysplastic syndromes. Br J Haemato 83: 743–748.
  45. Boula A, Voulgarelis M, Giannouli S, Katrinakis G (2006) Effect of CA2 antitumor necrosis factor-a antibody therapy on hematopoiesis of patients with myelodysplastic syndromes. Clin Cancer Res 12: 3099–3108.
  46. Itkin T, Rafii S (2017) Leukemic Cells “Gas Up” Leaky Bone Marrow Blood Vessels. Cancer Cell 32: 276–278. [crossref]
  47. Passaro D, Tullio AD, Abarrategi A, Rouault-Pierre K (2017) Increased vascular permeability in the bone marrow microenvironment contributes to disease progression and drug response in acute myeloid leukemia. Cancer Cell 32: 324–341.
  48. Counter CM, Mupta J, Harley CB (1995) Telomerase activity in normal leukocytes and hematological malignancies. Blood 85: 2315–2320.
  49. Fu C, Chen Z (2002) Telomerase activity in myelodysplastic syndrome. Chin Med J (Engl) 115: 1475–1478. [crossref]
  50. Delhommeau F, Dupont S, Della Valle V, James C, Trannoy S, et al. (2009) Mutation in TET2 in myeloid cancers. N Engl J Med 360: 2289–2301. [crossref]
  51. Lasson Cam Cote G, Quintas-Cardama A (2013) The changing mutational landscape of acute myeloid leukemia and myelodysplastic syndrome. Mol Cancer Res 11: 815–827.
  52. Papaemmanuil E, Gerstuug M, Malcovati L, Tauro S (2013) Clinical and Biological implication of driver mutations in myelodysplastic syndromes. Blood 122: 3616–3627.
  53. Bejar R, Lord R, Stovenson K, Bar-Natan M (2014) TET2 mutations predict response to hypomethylating agents in myelodysplstic syndrome patient. Blood 124: 2705–2712.
  54. Kennedy JA, Ebert BL (2017) Clinical implication of genetic mutations in myelodysplastic syndrome. J Clin Oncol 35: 968–974.
  55. Woll PS, Kjällquist U, Chowdhury O, Doolittle H, Wedge DC, et al. (2014) Myelodysplastic syndromes are propagated by rare and distinct human cancer stem cells in vivo. Cancer Cell 25: 794–808. [crossref]
  56. Maeck L, Haase D, Schoch C, Hiddemann W (2000) Genetic instability in myelodysplastic syndrome: detection of microsatellite instability and loss of heterozygosity in bone marrow sample with karyotype alteration. Br J Haemato 109: 842–846.
  57. Xie D, Hefmann KW, Mori N, Miller CW (2001) Allotype analysis of the myelodysplastic syndrome. Leukemia 14: 805–810.
  58. Kumamoto K, Ban S, Oda K, Tanaka H (2002) Chromosomal instability and radiosensitivity in myelodysplastic syndrome. Leukemia 16: 2253–2258.
  59. Takeda A, Goolsby C, Yaseen NR (2006) NUP98-HOXA9 induces long-term proliferation and blocks differentiation of primary human CD34+ hematopoietic cells. Cancer Res 66: 6628–6637. [crossref]
  60. Feng F, Li Q, Ling C, Zhang Y (2005) Phase III clinical trials of the cell differentiation agent-2 (CDA-2): therapeutic efficacy on breast cancer, non-small- cell lung cancer, and primary hepatoma. Chin J Clin Oncol 2: 706–716.
  61. Santi DV, Norment A, Garrett CE (1984) Covalent bond formation between a DNA-cytosine methyltransferase and DNA containing 5-azacytosine. Proc Natl Acad Sci U S A 81: 6993–6997. [crossref]
  62. Ma J (2007) Differentiation therapy of malignant tumor and leukemia. CSCO Treaties on the Education of Clinical Oncology 480–486.
  63. Liau MC, Fruehauf JP (2015) Destabilization of abnormal methylation enzymes as a critical mechanism for CDA-2 to reverse MDS progression. The First International Forum of Myelodysplastic Syndrome in Tai Zhou, Jiangsu, China: 31–35.
  64. Prasanna P, Shack S, Wilson VL, Samid D (1995) Phenylacetate in chemoprevention: in vitro and in vivo suppression of 5-aza-2’-deoxycytidine-induced carcinogenesis. Clin Cancer Res 1: 865–871. [crossref]
  65. Jordan CT, Guzman ML, Noble M (2006) Cancer stem cells. N Engl J Med 355: 1253–1261. [crossref]
  66. Vermeulen L, Sprick MR, Kemper K, Stassi G, Medema JP (2008) Cancer stem cells–old concepts, new insights. Cell Death Differ 15: 947–958. [crossref]
  67. Zabierowski SE, Herlyn M (2008) Melanoma stem cells: the dark seed of melanoma. J Clin Oncol 26: 2890–2894. [crossref]
  68. Hemmings C (2010) The elaboration of a critical framework for understanding cancer: the cancer stem cell hypothesis. Pathology 42: 105–112. [crossref]
  69. Magee JA, Piskounova E, Morrison SJ (2012) Cancer stem cells: impact, heterogeneity, and uncertainty. Cancer Cell 21: 283–296. [crossref]
  70. Zhou S, Schuetz JD, Bunting KD, Colapietro AM, Sampath J, et al. (2001) The ABC transporter Bcrp1/ABCG2 is expressed in a wide variety of stem cells and is a molecular determinant of the side-population phenotype. Nat Med 7: 1028–1034. [crossref]
  71. Zhou M, Atkinson KL, Rosen JM (2010) Selective targeting of radiation resistant tumor initiating cells. Proc Natl Acad Sci USA 107: 3522–7.
  72. Moitra K, Lou H, Dean M (2011) Multidrug efflux pumps and cancer stem cells : in sight into multidrug resistance and therapeutic development. Clin Pharmacol Ther 89: 491–502.
  73. Frame FM, Maitland NJ (2011) Cancer stem cells, models of study and implications of therapy resistance mechanisms. Adv Exp Med Biol 720: 105–118. [crossref]
  74. Liau MC (2007) Pharmaceutical composition inducing cancer cell differentiation and the use for treatment and prevention of cancer thereof. US Patent 7232578 B2, June 19.
  75. Liau MC, Ashraf AQ, Lee SS, Burzynski SR (1990) Riboflavin as a minor active anticancer component of Antineoplaston A2 and A5. Intl J Tiss React 12: 19–26.
  76. Liau MC, Huang LJ, Lee JH, Chen SC (1998) Development of differentiation helper inducers for the differentiation therapy of cancer. Chin Pharm J 50: 289–303.
  77. Liau MC, Liau CP (2002) Methyltransferase inhibitors as excellent differentiation helper inducers for differentiation therapy of cancer. Bull Chin Cancer 166–168.
  78. Laerum OD, Paukovits WR (1984) Modulation of murine hemopoiesis in vivo by a synthetic hemoregulatory pentapeptide (HP5b). Differertiation 27: 106–112.
  79. Bourgeade MF, Besançon F (1984) Induction of 2’,5’-oligoadenylate synthetase by retinoic acid in two transformed human cell lines. Cancer Res 44: 5355–5360. [crossref]
  80. Das BC, Thapa P, Karki R, Das S, Mahapatra S, et al. (2014) Retinoic acid signaling pathways in development and diseases. Bioorg Med Chem 22: 673–683. [crossref]
  81. Cunningham TJ, Duester G (2015) Mechanisms of retinoic acid signalling and its roles in organ and limb development. Nat Rev Mol Cell Biol 16: 110–123. [crossref]
  82. Gudas LJ, Wagner JA (2011) Retinoids regulate stem cell differentiation. J Cell Physiol 226: 322–330. [crossref]
  83. Stockhausen MT, Kristoffersen K, Stobbe L, Poulsen HS (2014) Differentiation of glioblastoma multiforme stem-like cells leads to down-regulation of EGFR and EGFRviii and decreased tumorigenic and stem-like cell potential. Cancer Biol Ther 15: 216–224.
  84. Sulaiman G, Cooke A, French B (2017) MyD88 is an essential component of retinoic acid-induced differentiation in human pluripotent embryonal carcinoma cells. Cell Death Diff 24: 1975–1986.
  85. El-Houjeiri L, Saad W, Hayar B, Aouad P (2017) Antitumor effect of the atypical retinoid ST1926 in acute myeloid leukemia and nanoparticle formation prolongs lifespan and reduce tumor burden of xenograft mice. Mol Cancer Ther 16: 2047–2057.
  86. Heynen GJ, Nevedomskaya E, Palit S, Jagalur Basheer N, Lieftink C, et al. (2016) Mastermind-Like 3 Controls Proliferation and Differentiation in Neuroblastoma. Mol Cancer Res 14: 411–422. [crossref]
  87. Yang Q, Sakurai T, Kakudo K (2002) Retinoid, retinoic acid receptor beta and breast cancer. Breast Cancer Res Treat 76: 167–173. [crossref]
  88. Zhao X, Graves C, Ames SJ, Fisher DE, Spanjaard RA (2009) Mechanism of regulation and suppression of melanoma invasiveness by novel retinoic acid receptor-gamma target gene carbohydrate sulfotransferase 10. Cancer Res 69: 5218–5225. [crossref]
  89. Lotem J, Sachs L (1979) Regulation of normal differentiation in mouse and human myeloid leukemia cells by phorbol esters and the mechanism of tumor promotion. Proc Natl Acad Sci USA 76: 5158–5162.
  90. Cooper RA, Braunwald AD, Kuo AL (1982) Phorbol ester induction of leukemic cell differentiation is a membrane-mediated process. Proc Natl Acad Sci U S A 79: 2865–2869. [crossref]
  91. Weeks CE, Slaga TJ, Hemmings H, Gleason GL (1979) Inhibition of phorbol ester-induced tumor promotion in mice by vitamin A analogs and anti-inflammatory steroid. J Natl Cancer Inst 63: 401–406.
  92. Han ZT, Zhu XX, Yang RY, Sun JZ (1998) Effect if intravenous infusion of 12-O-tetradecanoyl phorbol-13-acetate (TPA) in patients with monocytic leukemia: preliminary studies on therapeutic efficacy and toxicity. Proc Natl Acad Sci USA 95: 5357–5361.
  93. Schaar DG, Goodell L, Aisner J, Cui XX (2006) Phase I clinical trial of 12-O-tetradecanoylphorbol-13-acetate for patients with relapsed/refractory malignancies. Cancer Chemother Pharmacol 57: 789–795.
  94. Bond JA, Gescher A, Verschoyle RD, Lemoine NR (2007) Cytotoxic action of phorbol esters on human pancreatic cancer cells. Int J Cancer 121: 1445–1454.
  95. Xu L, Zhang L, Hu C, Liang S (2016) WNT pathway inhibitor pyrvinium pamoate inhibits the self-renewal and metastasis of breast cancer stem cells. Intl J Oncol 48: 1175–1186.
  96. Stoddart A, Wang J, Hu C, Fernald AA, Davis EM, et al. (2017) Inhibition of WNT signaling in the bone marrow niche prevents the development of MDS in the Apcdel/+ MDS mouse model. Blood 129: 2959–2970. [crossref]
  97. Momtazi-Borojeni AA, Abdollahi E, Ghasemi F, et al. (2018) The novel role of pyrvinium in cancer therapy. J Cell Physiol 233: 2871–2881. [crossref]
  98. Barbarino M, Cesari D, Intruglio R, Indovina P (2018) Possible repurposing of pyrvinium pamoate for the treatment of mesothelioma: a preclinical assessment. J Cell Physiol 233: 7391–7401.
  99. Yao CJ, Lai GM, Chen CF (2007) Elimination of cancer stem-like “side population” cells existed in MCF-7 breast cancer cells by human urine extract CDA-2. Literature documentation on CDA-2 issued by Everlife Pharmaceutical Co, Hefei, Anhui, China, 2007; 5–11.
  100. Ren J (2009) Cancer stem cell: the new target for fighting cancer from the latest transitional research. Biomed Pharmaco Therapy 63: 626–630.
  101.  Mougeri-Sacca M, Dimarino S, De Maria R (2013) Biological and clinical implication of cancer stem cells in primary brain tumors. Front Oncol 3: 6.
  102. Lathia JD, Venere M, Rao MS, Rich JN (2011) Seeing is believing: are cancer stem cells the Loch Ness monster of tumor biology? Stem Cell Rev 7: 227–237. [crossref]
  103. Fruehauf JP, Liau MC (2011) Therapeutic targeting of HIF to reverse the cancer stem cell phenotype and epithelial-mesenchymal transition. Therapy 8: 737–740.
  104. Caramels J, Papadogeorgakis E, Hill L, Browne GJ, Richard G, et al. (2013) A switch in the expression of embryonic EMT-inducers drives the development of malignant melanoma. Cancer Cell 24: 466–480. [crossref]
  105. Smigiel JM, Parameswaran N, Jackson MW (2017) Potent EMT and CSC Phenotypes Are Induced By Oncostatin-M in Pancreatic Cancer. Mol Cancer Res 15: 478–488. [crossref]
  106. Nomura A, Gupta VK, Daner P, Sharma NS (2018) NFkB-mediated invasiveness in CD133+ pancreatic TICs is regulated by autocrine and paracrine activation of IL1 signaling. Mol Cancer Res 16: 162–72.
  107. Kampalath BN, Liau MC, Burzynski B, Burzynski SR (1987) Chemoprevention by antineoplaston A10 of benzo(a)pyrene-induced pulmonary neoplasia. Drugs Exp Clin Res 1: 51–55. [crossref]
  108. Kampalath BN, Liau MC, Burzynski B, Burzynski SR (1990) Protective effect of Antineoplaston A10 in hepatocarcinogenesis induced by aflatoxin B1. Intl J Tiss React 12: 43–50.
  109. Liau MC, Chang CF, Giovanella BC (1980) Demonstration of an altered S-adenosylmethionine synthetase in human malignant tumors xenografted into athymic nude mice. J Natl Cancer Inst 64: 1071–1075.
  110. Knudson AG Jr (1971) Mutation and cancer: statistical study of retinoblastoma. Proc Natl Acad Sci U S A 68: 820–823. [crossref]
  111. Liau MC, Chang CF, Becker FF (1979) Alteration of S-adenosylmethionine synthetase during chemical hepatocsarcinogensis and in resulting carcinoma. Cancer Res 39: 2103–2109.
  112. Morley JE (2001) Hormone, aging, and endocrines in the elderly. In: Felig P, Frohman LA, “endocrinology and Metabolism, 4th (edn), McGrow-Hill, Inc., Medical Publishing Division, New York, Pg No: 1455–1482.
  113. Liau MC, Leong Y, Liau CP, Burzynski SR (1992) Prevention of Drug-induced DNA hypermethylation by antineoplaston components. Intl J Exptl Clin Chemother 5: 19–27.
  114. Bagheri V, Memar B, Behzadi R, Aliakbarian M (2018) Isolation and identification of chemotherapy-enriched sphere forming cells from a patient with gastric cancer. J Cell Physiol 233: 7036–7046.

Mithramycin-A Induced Toxicity in HepG2 Cells is mediated by Disruption of Calcium Homeostasis

DOI: 10.31038/JPPR.2019211

Abstract

Mithramycin A is a potent inhibitor of EWS-FLI1, a transcription factor implicated in Ewings sarcoma. However, a clinical trial initiated to evaluate its efficacy revealed hepatotoxicity at doses well below those required for inhibition of EWS-FLI1 activity. In the present study, we used a chemogenomic screen against gene deletion mutants of the yeast Saccharomyces cerevisiae to generate hypothesis on its cellular mechanism of hepatotoxicity. Our findings show that it disrupts calcium homeostasis in S. cerevisiae.  Studies in HepG2 cells grown in monolayer and as spheroids confirmed that it not only induced sustained elevated cytosolic Ca2+ levels but also induced endoplasmic reticulum stress. Cells exposed to the compound were ultimately found to undergo calpain mediated apoptosis. These data suggest that mithramycin A is a direct toxin to liver cells and its toxicity is mediated, at least in part, by disruption of Ca2+ homeostasis.

Keywords

Mithramycin A; Hepatotoxicity; HepG2 cells; Calcium homeostasis

Introduction

Mithramycin A (MTA) is a microbial product first isolated from fermentation of Streptomyces plicatus [1]. In a previous work employing high throughput screening aimed at the discovery of small molecule inhibitors of the transcription factor EWS-FLI1, a therapeutic target in Ewings sarcoma, we reported the identification of MTA as the highest scoring inhibitor. EWS-FLI1, an aberrant transcription factor arising from a reciprocal chromosomal translocation mutation, is only found in Ewings sarcoma tumor cells and, as such, represents an attractive therapeutic target. As with most cancer related transcription factors, however, it has been difficult to target because of lack of surface involutions suitable for high affinity binding of small molecules. Based on results from studies in xenografts and reports of clinical activity in the early seventies [2,3], a combined phase I/II clinical trial was initiated. However, MTA was found to induce liver toxicity in patients at a plasma concentration substantially below that required to suppress EWS-FLI1 activity [4]. The discovery of MTA as a potent inhibitor of EWS-FLI1 activity has made targeting this transcription factor tractable while also opening avenues for the development of MTA analogs with improved liver toxicity profiles. We recently tested several biosynthetic analogs of MTA and identified two analogs with improved targeting of EWS-FLI1 [5]. Though its mechanism of liver toxicity is not well understood, reports of hepatotoxicity in humans along with those observed in animal models suggest that MTA may be a direct hepatotoxin [6]. A better understanding of the mechanism of liver toxicity of MTA and analogs currently under preclinical testing is needed to identify suitable analogs for further development.

Materials and Methods

Reagents, Cells and Culture Conditions

MTA (Cayman Chemical, Cat. #11434), calcium orange (ThermoFisher, Cat. #C3015), HepG2 (ATCC, Cat. #HB-8065), ionomycin (Sigma, Cat. #407953) were all commercially obtained. Caspase3/7 spheroid staining reagent was obtained from Nexcelom (Cat. #CS1-V0002-1) while luciferase reporter constructs for ATF6 (Cat. #E3661), NFAT (Cat. #E8481), and CRE (Cat. #E8471) were all obtained from Promega.

MTA Yeast Chemogenomic Screen

MTA Minimum Inhibitory Concentration (MIC) in the wild type yeast BY4743 (MATa/α his3Δ1/his3Δ1 leu2Δ0/leu2Δ0 LYS2/lys2Δ0 met15Δ0/MET15 ura3Δ0/ura3Δ0) was determined using the broth microdilution assay [7]. A range of sub-inhibitory concentrations (2.4 – 3mM) were first tested to determine a dose that causes 10–20% growth inhibition. This concentration was used for a genome-wide screening of Saccharomyces cerevisiae pooled deletion collection, composed of essential heterozygous (~1150 strains) and homozygous (~4800 strains) diploid gene deletion mutants. These were separately diluted to an OD600 of 0.06 in 900µL of YPD medium containing either 2.7mM MTA or DMSO and were grown in 48-well plates for 10 and 20 generations at 30 °C, respectively. The heterozygous and homozygous samples were then combined and 3 OD600 unit of cells for each sample were collected for genomic DNA extraction, PCR amplification of molecular barcodes, microarray hybridization and scanning as previously described [8,9]. Two independent biological replicates were performed for each sample. Intensity values for the barcodes on the TAG4 arrays are extracted by the GeneChip operating software (Affymetrix) and log2 ratio of the median normalized intensity values of each barcode in drug treated vs. DMSO control was calculated using GeneSpring software (ver. 13.0). This data has been deposited at NCBI’s Gene Expression Omnibus inder accession number GSE122458. Gene ontology (GO) enrichment analysis with BiNGO (ver. 3.0.3), a Cytoscape (ver. 3.3.0) plug-in [10], was used for functional profiling of susceptible mutants in our large-scale experiment. Significantly enriched biological processes and cellular components were evaluated using a hypergeometric test corrected for multiple hypothesis testing (P < 0.0001) using a Benjamini-Hochberg false discovery rate (FDR) correction and then visualized using the Enrichment Map (ver. 2.0.1) plug-in [11] developed for Cytoscape. Nodes in maps represent enriched biological processes or cellular components (hypergeometric test FDR < 0.0001) with edges indicating gene overlap between enriched processes/cellular components. Line width in maps is proportional to degree of overlap and edges are not shown where the overlap coefficient is less than 0.5.  GO biological processes and cellular components that were too specific (containing less than five genes) or too general (contain greater than 300 genes) were excluded from the maps.

Confirmation of MTA-sensitive Strains

Selected strains identified as MTA-sensitive in the primary screen were subsequently confirmed by growing individual isolates along with the parental strain BY4743 in 100 µL of YPD with or without 2.7mM MTA in 96-well plates for 20 hours at 30 °C and 200 rpm on a shaking incubator. The OD600 for treated and untreated samples was obtained using a BMG OmegaStar plate reader and the relative growth rate determined as the growth ratio of drug treated to untreated controls.

Generating and Imaging HepG2 Spheroids

HepG2 cells being grown in T225 flasks were harvested via trypsinization and pelleted via centrifugation. After removal of the supernatant, the cells were resuspended at 1.75×104 cells per milliliter in EMEM medium containing 10% FBS. 40 µL of this cell suspension was seeded into black-walled clear bottomed ultralow-attachment 384-well plates (Corning, #3830). Cells were incubated for 52–58 hours or until well-defined spheroids were formed. Control and MTA treated spheroids were stained with commercial Hoechst and caspase 3/7 reagents (Nexcelom, #CSK-V0003-1) following the manufacturer’s protocol. Data on spheroid size, shape, fluorescence and bright field images were all acquired on a BioTek Cytation 3 imager using a 10X objective.

JC1, ATP, and Luciferase Reporter Assays

For ATP and luciferase reporter microplate reader-based assays, HepG2 cells were seeded into 384-well plates at a density of 5000 cells per well in 27 µL of phenolred-free growth media. After incubation overnight, cells in these plates were treated with either 10 µM MTA (final concentration), 1 µM ionomycin (final concentration) or DMSO control for 24 hours prior to the addition of ATP reagent as described previously [12]. Eight replicates per condition and two independent experiments were performed. Setup of the reporter assay has been described previously [5]. But briefly, the transcription factor luciferase reporter or minimal-promoter luciferase reporter (negative control), or constitutively active (CMV) luciferase reporter (positive control) were transfected into HepG2 cells (5,000 cells/well) along with a constitutively active renilla luciferase construct in a white 96-well plate. After incubation for 24 hours, medium in each well was replaced 10 µM MTA or DMSO in either standard medium or Ca2+-free MEM spinner medium (Quality Biological, Cat. #112-021-101). After incubation for an additional 24 hours, firefly and renilla luciferase signals were read using Dual-Glo reagent (Promega, Cat. #E2920). Minimal promoter (minP) and constitutively active promoter (CMV) luciferase activity were used as controls to assess background reporter activity and general transcriptional activity, respectively. For the JC1 flow cytometry measurement of changes in mitochondrial membrane potential, HepG2 cells treated with MTA, ionomycin or DMSO for 24 hours were harvested and washed with Hank’s Balanced Salt Solution (HBSS). Pelleted cells were then stained by resuspending in 20 µM JC-1 in HBSS and incubating for 10 minutes at 37 °C. After pelleting, removal of the loading solution, and washing with HBSS, cells were finally resuspended in HBSS buffer. Samples were analyzed on an Accuri C6 flow cytometer. The data acquired were analyzed using FCS Express 4 software.  Four individual samples for each condition were tested and results are reporter as averages from two independent experiments.

Caspase-3/7 and Calpain Activation Assays

HepG2 cells were seeded at 3000 and 8000 cells per well for caspase3/7 and calpain assays, respectively, in 27 µL medium in white 384-well plates and incubated overnight. For caspase assays 3µL of 10X of MTA and controls were added and plates were incubated for 24 hours. After equilibrating to room temperature, caspase activation assay was run following the manufacturer’s protocol (Promega, #G8091) by adding 30 µL of caspase assay reagent, incubating for 30 minutes and reading luminescence using a plate reader. Six replicate wells per treatment were used. For calpain activation assay, after seeding and incubating cells overnight, cells were treated with 1 µL of BAPTA-AM (5 µM final concentration) and incubated for 30 minutes. Following this, calpain activation was measured using an assay kit (Promega, #G8501) with the following modifications in the recommended protocol. 1 µL of the synthetic calpain substrate Suc-LLVY-amino luciferin in Calpain-Glo buffer (20 µM final concentration) was first added

Western Blotting

HepG2 cells were seeded into T75 ultralow attachments Flasks (Corning, Cat. #3814) and allowed to form spheroids by placing in an incubator for 72 hours. MTA or DMSO control was added to these and flasks were allowed to incubate for an additional 48 hours. Spheroids were harvested as pellets via centrifugation. Pellets were then lysed using NP-40 lysis buffer (Thermo, # FNN0021) supplemented with PMSF (1 mM final concentration) and protease inhibitor cocktail (Thermo, #78430). Prepared Lysates were cleared by centrifugation and 20 mg of protein was loaded onto gels. Gels were blotted onto nitrocellulose membranes. These membranes were probed with GRP78 (SantaCruz, #sc-13539) antibody. Blots were then probed with secondary antibodies tagged with IRDye 680 and IRDye 800 prior to scanning with an Odyssey infra-red scanner.

Statistical Analysis

All data reported is shown as the sample mean ± the Standard Deviation (SD). Pairwise comparisons between means of controls and treatment were performed using a Student t-test (two tailed, unpaired, unpaired) where, for each couple of normally distributed populations, the null hypothesis that the means are equal were verified. Difference between control and treatment data was considered statistically significant if the Student t-test gives a significance level P (P value) less than 0.05.

Results

MTA Impacts Cation Homeostasis in S. cerevisiae

The complete pool of barcoded nnon-essential homozygous and essential heterozygous diploid deletion strains of S. cerevisiae was used to identify gene deletions that confer sensitivity to MTA. However, prior to screening this collection, a sub-inhibitory concentration of MTA that resulted in 10–20% (IC10–20) growth of the BY4743 parental strain, 2.7 mM MTA (Figure 1A), was determined under the same condition used for screening the pooled collection. Screening the pooled collection at this concentration and at a log2 ratio cut-off of 2, we found 48 homozygous gene deletion mutants whose growth rates were significantly inhibited in the presence of MTA compared to control treatment. None of the essential heterozygous deletion strains in our screen showed sensitivity to MTA suggesting lack of a specific protein target through which it exerts its inhibitory action in yeast. Hypersensitivity of non-essential homozygous deletion mutants that are involved in the cellular cation homeostasis pathway, however, pointed to the importance of cation homeostasis in resistance to the compound. An analysis of supersensitive mutants indicated a statistically significant enrichment for Gene Ontology (GO) biological process terms relevant to cation homeostasis (Table 1, Figure 1B). Hypersensitivity of deletion mutants in this major functional group to sub-inhibitory concentration of MTA was confirmed by retesting cherry-picked representative individual deletion mutants as shown in Figure 1C. Classification of MTA-sensitive mutants based on GO cellular component revealed significant enrichment for the vacuolar proton-transporting V-type ATPase complex (P-value = 6.92E-07) (Table 2, Figure 1D). Cytosolic Ca2+ homeostasis is a known constitutive function of yeast V-ATPase [13]. Taken together, the MTA chemical-genetic profile presented herein, demonstrates the compound’s impact on cation homeostasis pathway in yeast in general and Ca2+ homeostasis in particular.

Table 1. GO biological process term enrichment in Mithramycin A sensitive strains.

GO Biological Process

P-Value

Systematic Name/Standard Name

Cellular monovalent inorganic cation homeostasis

4.35E-09

YCL005W-A/VMA9, YJL129C/TRK1, YBR127C/VMA2, YEL051W/VMA8, YKL119C7VPH2, YEL027W/VMA3, YGL095C/VPS45, YLR447C/VMA6

Monovalent inorganic cation homeostasis

4.35E-09

Cellular cation homeostasis

8.52E-07

YGL167C/PMR1, YCL005W-A/VMA9, YJL129C/TRK1, YCR044C/PER1, YBR127C/VMA2, YFL051W/VMA8, YKL119C/VPH2, YEL027W/VMA3, YGL095C/VPS45, YLR447C/VMA6

Cellular Ion homeostasis

1.78E-06

Cation homeostasis

1.78E-06

Ion homeostasis

3.80E-06

Cellular chemical homeostasis

3.80E-06

Cellular homeostasis

1.56E-05

Chemical homeostasis

1.55E-05

Monovalent inorganic cation transport

1.12E-05

YCL005W-A/VMA9, YJL129C/TRK1, YBR127C/VMA2, YEL051W/VMA8, YEL027W/VMA3, YLR447C/VMA6

Intracellular pH reduction

4.35E-09

YCL005W-A/VMA9, YBR127C/VMA2, YEL051W/VMA8, YKL119C/VPH2, YEL027W/VMA3, YGL095C/VPS45, YLR447C/VMA6

pH reduction

4.35E-09

vacuolar acidification

4.35E-09

Regulation of intracellular pH

1.44E-08

Regulation of cellular pH

1.44E-08

Regulation of pH

2.00E-08

Energy coupled proton transport, against electrochemical gradient

2.60E-05

YCL005W-A/VMA9, YBR127C/VMA2, YEL027W/VMA3, YLR447C/VMA6

ATP hydrolysis coupled proton transport

2.60E-05

Proton transport

4.81E-05

YCL005W-A/VMA9, YBR127C/VMA2, YEL051W/VMA8, YEL027W/VMA3, YLR447C/VMA6

Hydrogen transport

5.02E-05

Endocytosis

3.77E-05

YNL297C/MON2, YGR167W/CLC1, YBR164C/ARL1, YEL027W/VMA3, YCR028C/FEN2, YLR240W/VPS34, YLR337C/VRP1

Vacuolar transport

5.02E-05

YLR261C/VPS63, YNL297C/MON2, YML071C/COG8, YBR164C/ARL1, YEL027W/VMA3, YGL095C/VPS45, YLR447C/VMA6, YLR372W/EL03

Table 2. GO cellular component term enrichment in MTA sensitive strains.

GO cellular component

P-Value

Systematic Name/Standard Name

vacuolar proton-transporting V-type ATPase complex

6.92E-07

YCL005W-A/VMA9, YBR127C/VMA2, YEL051W/VMA8, YEL027W/VMA3, YLR447C/VMA6

proton-transporting V-type ATPase complex

6.92E-07

proton-transporting two-sector ATPase complex

2.90E-05

vacuolar membrane

3.06E-05

YCL005W-A/VMA9, YCR044C/PER1, YBR127C/VMA2, YEL051W/VMA8, YJL154C/VPS35, YEL027W/VMA3, YCR028C/FEN2, YGL095C/VPS45, YLR240W/VPS34, YLR447C/VMA6

JPPR - 106_Girma Woldemichael_F1

Figure 1. A) An IC10–20 of 2.7mM was determined after prescreening parental (wild-type) yeast against MTA. This concentration was selected for genome-wide screening. B) Map of GO biological process enriched in the MTA chemogenomic profile. C) Confirmation of hypersensitivity of selected strains to MTA. Htb1 deletion mutant was used as negative control. Results represent Mean ± standard deviation of quadruplicates. D) MTA GO cellular component enrichment map.

MTA Induces Persistent Elevation of Cytosolic Calcium Levels in HepG2 Cells

Many proteins involved in acquiring, utilizing, storing, and regulating levels of inorganic ions are functionally conserved from yeast to man [14]. We, therefore, asked whether MTA also had an impact on Ca2+ homeostasis in mammalian liver cells. To determine effect on Ca2+ homeostasis, we used HepG2 cells grown both in monolayer and as spheroids, which are well-suited to studying direct toxicity to liver cells by small molecules [15]. We found that MTA showed a dose dependent reduction in spheroid size and cell viability in HepG2 cells grown in monolayer (Figure 2A). We also found that treatment with MTA resulted in immediate and sustained elevation of cytosolic calcium levels in HepG2 cells cultured in monolayer (Figure 2B). NFAT is known to be expressed in HepG2 cells and contributes to cell proliferation as part of the Ca2+/calcineurin/NFAT signaling pathway [16]. Therefore, a luciferase reporter for NFAT activity was also used to gauge MTA’s impact on Ca2+ levels. This was done in both Ca2+-free and Ca2+-containing medium. The results showed that while MTA treatment significantly increased NFAT reporter activity in HepG2 cells in Ca2+-containing medium, it induced an even greater reporter response in Ca2+-free medium (Figure 2C). Pretreatment with EGTA to chelate extracellular Ca2+ or treatment in Ca2+-free media only partially reversed MTA’s effect on HepG2 spheroid size and shape (Figure 2D). Because of this observation, we next examined any potential effects on the Endoplasmic Reticulum (ER), the largest intracellular Ca2+ store. Increased ATF6 expression has been shown to be a marker of ER stress in HepG2 cells [17]. An ATF6 luciferase reporter construct used to assess changes in ATF6 activity in response to MTA treatment in HepG2 cells being grown in Ca2+-containing and Ca2+-free medium showed a 4.8- and 6.9-fold increase in reporter activity, respectively, on MTA treatment relative to vehicle (Figure 2E). Changes in levels of the ER stress marker protein GRP78 were also readily apparent in lysates from MTA-treated spheroids as shown in Figure 2F [17,18]. These findings appear to implicate increase in cytosolic Ca2+ levels owing to both influx of Ca2+ from the extracellular environment and intracellular stores. Together, these findings also suggest that, similar to observations in yeast, MTA causes disruption of Ca2+ homeostasis in HepG2 cells.

JPPR - 106_Girma Woldemichael_F2

Figure 2. MTA induces rise in cytosolic Ca2+ levels in HepG2 cells. (A) MTA showed dose dependent inhibitory effect on growth of HepG2 spheroids (3D) and cells grown in monolayer (2D). Representative fluorescence imaging of control (DMSO) and MTA treated cells (B) grown as monolayer and loaded with the calcium binding dye calcium Orange AM show increased fluorescence intensity over time in MTA treated cells indicating increase in cytosolic calcium levels.  MTA treatment (10 mM) for 24 hours also resulted in increased luciferase reporter activity in HepG2 cells transiently transfected with a reporter construct for either NFAT (C) or ATF-6 (E) in both Ca2+-containing and Ca2+-free medium relative to control treatment.  (D) Representative images of spheroids of control and MTA (10 mM, 48h) treated HepG2 cells in either Ca2+-free medium or in medium pretreated with 1.8 mM EGTA to chelate extracellular Ca2+ showed partial reversal of MTA’s effect on spheroid size and shape. Western blot of lysates from spheroid treatment with MTA in both Ca2+-containing and Ca2+-free medium showed upregulation of markers of ER stress (F).

MTA Induces Cell Injury in HepG2 Cells

Dysregulation of Ca2+ homeostasis has been implicated in the induction of apoptosis via involvement in onset of the Mitochondrial Permeability Transition (MPT) [19]. As a result, we also looked at whether MTA induced MPT using JC1 to assess changes in the mitochondrial membrane potential. Treatment of HepG2 cells with MTA, however, resulted in an increase in JC1 aggregate signal (Figure 3A) suggesting hyperpolarization of the mitochondrial membrane. We sought corroboration for this finding by quantifying changes in cellular ATP content since ATP production requires the existence of mitochondrial membrane potential. Treated HepG2 cells grown either as monolayer (Figure 3B) or as spheroids (Figure 3C) were both found to have increased total ATP content consistent with findings of increased mitochondrial membrane potential.

JPPR - 106_Girma Woldemichael_F3

Figure 3. MTA increases Mitochondrial potential. A) HepG2 cells treated with 10µM MTA for 24h were stained with JC-1. 1µM of the calcium ionophore ionomycin was used as a positive control. Shown are representative plots of experiments.  Total ATP content in HepG2 cells grown as monolayer (B) and spheroids (C) was measured using ATP fueled luciferase activity 24h after treatment with 10 µM MTA.

Increased mitochondrial membrane potential has been shown to be an early event in the induction of apoptosis in response to treatment by compounds in HepG2 cells [20]. Spheroid imaging and size measurements in many of the experiments conducted also revealed a decrease in size and change in shape upon MTA treatment suggesting impact on cell proliferation and survival. As a result, we first looked at whether MTA induced apoptosis. Imaging of MTA treated HepG2 spheroids stained with caspase 3/7 reagent measuring caspase 3/7 activity revealed that treated cells were positively stained (Figure  4A). Caspase 3/7 activity in MTA treated HepG2 cells grown as a monolayer were also found to be induced 3-fold 48 hours after treatment (Figure 4B). No increase in the activity of caspases 8 or 9 was observed in MTA treated cells. Since induction of apoptosis via activation of caspase 7 by calpain has been implicated in response to increase in cytosolic Ca2+ concentrations [21], calpain activation in HepG2 cells was assessed. It was found that MTA induced a 3-fold increase in calpain activity. Calpain activation by MTA was reversed by pretreatment with the cell permeable calcium chelator BAPTA-AM (Figure 4C).  Pretreatment with BAPTA-AM was also found to inhibit caspase 3/7 (Figure 4D) activation suggesting that this activation is dependent on elevation of free cytosolic Ca2+ levels.  Together, the data on disruption of mitochondrial membrane potential and induction of apoptosis by MTA in HepG2 cells point to MTA being a direct toxin.

JPPR - 106_Girma Woldemichael_F4

Figure 4. MTA induces apoptosis in HepG2 cells.  A) Control, 10 µM MTA and 20 nM staurosporine (positive control) treated HepG2 spheroids were stained with both Hoechst 33342 DNA stain (purple) and caspase 3/7 stain (green) and imaged. Shown are representative images acquired 48h after treatment.   B) Relative caspase 3/7/8/9 activity 48h after treatment with 10 µM MTA.   C) Calpain activation was detected in response to MTA treatment in HepG2 cells using a luminescence detection kit. MTA induced activation was reversed by pretreatment of cells with the intracellular calcium chelator BAPTA-AM.   D) Pretreatment of HepG2 cells with 5 µM BAPTA-AM significantly reduced caspase-3/7 activation by in HepG2 cells treated with 10 µM MTA for 48h.

Discussion

In a previous effort to find novel agents targeting EWS-FLI1 in Ewings sarcoma, we identified MTA as the highest scoring lead compound. However, liver toxicity during clinical trials at subinhibitory concentrations halted development of MTA. A better understanding of the mechanism of its liver toxicity is needed for further development of more recent analogs with improved activity profile. To the best of our knowledge, this is the first study to describe a potential cellular mechanism for its hepatotoxicity.

In the present study, we performed a chemical genomics screen of non-essential homozygous and essential heterozygous gene deletion mutants of S. cerevisiae to gain insight into the MTA’s mode of action. Screening of this “disruptome” against bioactive compounds has been shown to provide valuable insight into target genes and gene involved in resistance to these compounds [22]. Through gene ontology analysis of sensitive strains, we showed that MTA’s chemogenomic profile revealed its involvement in dysregulation of calcium homeostasis through impact on V-ATPases, which are critical for generation of a pH gradient that drives secondary transporters to maintain cellular ion homeostasis. Comparison of MTA’s chemogenomic profile with that of Amiodarone’s shows close similarity [23,24]. Amiodarone’s antifungal activity is mediated by perturbation of calcium homeostasis with hypersensitivity of vma mutants in its chemogenomic profile ascribed to defects in ion homeostasis. Hypersensitivity of multiple vma deletion mutants encoding subunits of the vacuolar membrane H+-ATPase in the MTA chemogenomic profile (i.e., vma9Δ, vma2Δ, vma8Δ, vma3Δ and vma6Δ) not only indicates disruption of the cation homeostasis pathway by MTA but also underscores the critical role of V-ATPases in resistance to MTA’s inhibitory action [25].

There are many reports that link drug-induced liver toxicity to perturbation of calcium homeostasis. Elevation in cytosolic Ca2+ concentration is implicated in toxic liver injury associated with a number of compounds including diclofenac [26], senecionine and trans-4-OH-2-hexenal [27], halothane [28,29], and antivirals efavirenz and ritonavir [30]. In HepG2 cells, we found that MTA induced a sustained increase in cytosolic calcium levels through mobilization from the extra cellular medium and intracellular stores resulting in inhibition of cell proliferation. However, unlike several reports in different cell lines where such elevation causes mitochondrial depolarization and mitochondria mediated apoptosis, we found no alterations in the levels of or activation of mitochondria and apoptosis related proteins including Bcl-2 and BAX ruling out immediate involvement in induction of apoptosis by mitochondria. On the other hand, studies have shown that calpains are activated by sustained elevation of cytosolic calcium in HepG2 cells [31]. Our findings of calpain activation and partial reversal of the apoptosis in MTA treated cells with a calpain inhibitor are consistent with this observation and further underscore the critical role of impaired calcium homeostasis pathway in drug-induced direct hepatotoxicity.

Although its exact cellular mechanism of action has not been elucidated in mammalian cells, studies have shown that MTA forms complexes with GC-rich regions of DNA [32]. Formation of this complex is, however, dependent on the presence of divalent cations [33]. This DNA binding ability is thought to be responsible for its anticancer effects where carcinogenesis and disease progression is driven by transcription factors that preferentially bind GC-rich regions. Comparison of MTA’s haploinsufficiency and homozygous deletion profiles in yeast deletion mutants with those of well-characterized DNA binding compounds such as doxorubicin and actinomycin D [34], however, did not reveal enrichments in genes involved in DNA synthesis and repair. Western blotting analysis of lysates from MTA treated monolayer or spheroid HepG2 cultures also did not show activation of proteins such as ATM and ATR involved in DNA repair. These suggest that, its effect in HepG2 cells is primarily mediated through its impact on calcium homeostasis.

In summary, we propose that MTA acts as a direct hepatotoxin and that this activity ensues subsequent to dysregulation of calcium homeostasis and induction of ER stress resulting in apoptosis. This discovery will allow for a better characterization of hepatotoxicity much earlier in the discovery phase of synthetic and biosynthetic analogs currently being generated for targeting EWS-FLI1 in Ewings sarcoma.

Acknowledgement

This work has been funded in part with Federal funds from the Frederick National Laboratory for Cancer Research, National Institutes of Health, under contract HHSN261200800001E and in part by the Intramural Research Program of NIH, Frederick National Lab, and Center for Cancer Research. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does the mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government.

References

  1. Haskell T, Ryder A, P. Frohardt R, A. Fusari S, L. Jakubowski Z, et al. (1958) The Isolation and Characterization of Three Crystalline Antibiotics from Streptomyces plicatus. J Am Chem Soc 80: 743–747.
  2. Grohar P, Woldemichael G, Griffin L, Mendoza A, Chen Q-R, et al. (2011) Identification of an Inhibitor of the EWS-FLI1 Oncogenic Transcription Factor by High-Throughput Screening. Journal of the National Cancer Institute 103: 962–978.
  3. Kofman S, Perlia CP, Economou SG (1973) Mithramycin in the treatment of metastatic Ewing’s sarcoma. Cancer 31: 889–893. [crossref]
  4. Grohar PJ, Glod J, Peer CJ, Sissung TM, et al. (2017) A phase I/II trial and pharmacokinetic study of mithramycin in children and adults with refractory Ewing sarcoma and EWS-FLI1 fusion transcript. Cancer Chemother Pharmacol 80: 645–652. [crossref]
  5. Osgood CL, Maloney N, Kidd CG, Kitchen-Goosen SM, Segars LE, et al. (2016) Identification of mithramycin analogs with improved targeting of the EWS-FLI1 transcription factor. Clinical Cancer Research 22: 4105–4118. [crossref]
  6. King PD, Perry MC (2001) Hepatotoxicity of chemotherapy. Oncologist 6: 162–176. [crossref]
  7. CLSI (2008) Reference method for broth dilution antifungal susceptibility testing of yeasts. approved standard M27- A3, 3rd edn, Clinical and Laboratory Standards Institute, Wayne, PA ,USA.
  8. Ericson E, Hoon S, Onge RPS, Giaever G, Nislow C (2010) Exploring gene function and drug action using chemogenomic dosage assays. Methods enzymol 470: 233–255 [crossref]
  9. Pierce SE, Davis RW, Nislow C, Giaever G (2007) Genome-wide analysis of barcoded Saccharomyces cerevisiae gene-deletion mutants in pooled cultures. Nat protoc 2: 2958–2974. [crossref]
  10. Maere S, Heymans K, Kuiper M (2005) BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21: 3448–3449. [crossref]
  11. Merico D, Isserlin R, Stueker O, Emili A, Bader GD (2010) Enrichment map: a network-based method for gene-set enrichment visualization and interpretation. PLoS One 5: e13984. [crossref]
  12. Caropreso V, Darvishi E1, Turbyville TJ2, Ratnayake R1, Grohar PJ3, et al. (2016) Englerin A Inhibits EWS-FLI1 DNA Binding in Ewing Sarcoma Cells. J Biol Chem 291: 10058–10066. [crossref]
  13. Förster C, Kane PM (2000) Cytosolic Ca2+ Homeostasis Is a Constitutive Function of the V-ATPase in Saccharomyces cerevisiae. Journal of Biological Chemistry 275: 38245–38253. [crossref]
  14. Cyert MS, Philpott CC (2013) Regulation of cation balance in Saccharomyces cerevisiae. Genetics 193: 677–713. [crossref]
  15. Ramaiahgari SC, den Braver MW, Herpers B, Terpstra V, Commandeur JNM, et al. (2014) A 3D in vitro model of differentiated HepG2 cell spheroids with improved liver-like properties for repeated dose high-throughput toxicity studies. Archives of Toxicology 88: 1083–1095. [crossref]
  16. Wang S, Kang X, Cao S, Cheng H, Wang D, et al. (2012) Calcineurin/NFATc1 Pathway Contributes to Cell Proliferation in Hepatocellular Carcinoma. Digestive Diseases and Sciences 57: 3184–3188. [crossref]
  17. Bosquet A, Guaita-Esteruelas S, Saavedra P, Rodríguez-Calvo R, Heras M, Girona J, et al. (2016) Exogenous FABP4 induces endoplasmic reticulum stress in HepG2 liver cells. Atherosclerosis 249: 191–199. [crossref]
  18. Woldemichael GM, Turbyville TJ, Vasselli JR, Linehan WM, McMahon JB (2012) Lack of a Functional VHL Gene Product Sensitizes Renal Cell Carcinoma Cells to the Apoptotic Effects of the Protein Synthesis Inhibitor Verrucarin A. Neoplasia 14: 771-IN728. [crossref]
  19. Schanne FA, Kane AB, Young EE, Farber JL (1979) Calcium dependence of toxic cell death: a final common pathway. Science 206: 700–702. [crossref]
  20. Cao J, Liu Y, Jia L, Zhou H-M, Kong Y, Yang G, et al. (2007) Curcumin induces apoptosis through mitochondrial hyperpolarization and mtDNA damage in human hepatoma G2 cells. Free Radical Biology and Medicine 43: 968–975. [crossref]
  21. Gafni J, Cong X, Chen SF, Gibson BW, Ellerby LM (2009) Calpain-1 cleaves and activates caspase-7. J Biol Chem 284: 25441–25449. [crossref]
  22. Giaever G1, Nislow C (2014) The yeast deletion collection: a decade of functional genomics. Genetics 197: 451–465. [crossref]
  23. Gupta SS, Ton VK, Beaudry V, Rulli S, Cunningham K, et al. (2003) Antifungal activity of amiodarone is mediated by disruption of calcium homeostasis. The Journal of biological chemistry 278: 28831–28839. [crossref]
  24. Muend S, Rao R (2008) Fungicidal activity of amiodarone is tightly coupled to calcium influx. FEMS Yeast Res 8: 425–431. [crossref]
  25. Forster C, Kane PM (2000) Cytosolic Ca2+ homeostasis is a constitutive function of the V-ATPase in Saccharomyces cerevisiae. J Biol Chem 275: 38245–38253. [crossref]
  26. Lim MS, Lim PL, Gupta R, Boelsterli UA (2006) Critical role of free cytosolic calcium, but not uncoupling, in mitochondrial permeability transition and cell death induced by diclofenac oxidative metabolites in immortalized human hepatocytes. Toxicology and applied pharmacology 217: 322–331. [crossref]
  27. Griffin DS, Segall HJ (1987) Role of cellular calcium homeostasis in toxic liver injury induced by the pyrrolizidine alkaloid senecionine and the alkenal trans-4-OH-2-hexenal. Journal of biochemical toxicology 2: 155–167. [crossref]
  28. Farrell GC, Mahoney J, Bilous M, Frost L (1988) Altered hepatic calcium homeostasis in guinea pigs with halothane-induced hepatotoxicity. J Pharmacol Exp Ther 247: 751–756. [crossref]
  29. Yu WF, Yang LQ, Zhou MT, Liu ZQ, Li Q (2005) Ca2+ cytochemical changes of hepatotoxicity caused by halothane and sevoflurane in enzyme-induced hypoxic rats. World J Gastroenterol 11: 5025–5028. [crossref]
  30. Kao E, Shinohara M, Feng M, Lau MY, Ji C (2012) Human immunodeficiency virus protease inhibitors modulate Ca2+ homeostasis and potentiate alcoholic stress and injury in mice and primary mouse and human hepatocytes. Hepatology (Baltimore, Md) 56: 594–604. [crossref]
  31. Caro AA, Cederbaum AI (2002) Role of calcium and calcium-activated proteases in CYP2E1-dependent toxicity in HEPG2 cells. J Biol Chem 277: 104–113. [crossref]
  32. Albertini V, Jain A, Vignati S, Napoli S, Rinaldi A, et al. (2006) Novel GC-rich DNA-binding compound produced by a genetically engineered mutant of the mithramycin producer Streptomyces argillaceus exhibits improved transcriptional repressor activity: implications for cancer therapy. Nucleic acids research 34: 1721–1734. [crossref]
  33. Cons BM1, Fox KR (1989) Interaction of mithramycin with metal ions and DNA. Biochem Biophys Res Commun 160: 517–524. [crossref]
  34. Hoepfner D, Helliwell SB, Sadlish H, Schuierer S, Filipuzzi I, et al. (2014) High-resolution chemical dissection of a model eukaryote reveals targets, pathways and gene functions. Microbiological Research 169: 107–120. [crossref]

Might Actinomycin be Used to Cure All Cancer?

DOI: 10.31038/CST.2019411

In addition to curing Wilm’s tumor and choriocarcinoma, I propose actinomycin to be more widely used to cure ALL cancer – i.e., curing lung, colon, breast, ovarian and other cancers – for the reasons described in this communication.

The discovery that actinomycin is a powerful anticancer agent against Wilm’s tumor and choriocarcinoma was made by Sidney Farber in the early 1950’s. This miraculous discovery is known to dramatically cure patients suffering from these embryonic tumors — and still remains the treatment of choice to this day!

Although it is not known whether actinomycin can be used to treat ALL cancers – the purpose of this letter is to call attention to the possibility that trace amounts of actinomycin given over an extended period of time could prove to be a powerful anticancer chemotherapeutic regimen. It is important that the medical community realize this, and to understand how actinomycin acts to defeat cancer.

Background Information

Actinomycin D is a naturally occurring cyclic-polypeptide containing antibiotic known to bind to DNA and inhibit RNA synthesis (see Figure 1) 1–4]. It does this by interfering with the elongation of growing RNA chains by the RNA polymerase enzyme 5]. Nucleolar 5S ribosomal RNA-synthesis is known to be particularly sensitive to the presence of actinomycin, and this accounts for its pharmacological activity as well as its extreme toxicity to mammalian cells [6, 7].

CST-2019-101_F1

Figure 1. Chemical structure of actinomycin D. (L meval: L methyl valine; sar: sarcosine; L pro: L proline; D val: D valine; L thr: L threonine)

Stereochemistry of Actinomycin-DNA Binding

Several years ago, we determined the three-dimensional structure of an actinomycin-deoxyguanosine complex by x-ray crystallography [8–11]. The stereochemical information obtained from this study suggested a model to understand the general features of how actinomycin binds to DNA. According to this model, the phenoxazone ring system on actinomycin intercalates between adjacent guanine-cytosine base-pairs, while pentapeptide chains lie in the narrow groove of the B- helix to form hydrogen bonds with guanine residues on opposite chains. Implicit in this model was the assumption that actinomycin binds to B-DNA, or a distorted B-DNA form. The possibility that actinomycin might bind to some other discretely different DNA conformational-state was not envisioned at that time.

A modification to this actinomycin-DNA binding model was subsequently proposed, which allows one to understand its mechanism of action (see Figure 2). This model is similar to the previous one; however, it predicts that actinomycin to bind to (what we have called) beta-DNA (i.e., not to B-DNA), this being a metastable and hyperflexible premelted form inferred from our wider crystallographic studies of planar drug molecules intercalated into a series of DNA-like and RNA-like self-complementary dinucleoside-monophosphates [12–14].

CST-2019-101_F2

Figure 2. Actinomycin: beta-DNA binding model.

Figure 3 shows this same (extended) beta-DNA structure “pinned” by ethidium. The complex is an organized right-handed double helical structure in which the beta-structural element plus the intercalator form the asymmetric unit of the helix. This maximally elongated and unwound DNA duplex-structure, pinned by ethidium at saturating concentrations, readily explains the well-known observation of neighbor-exclusion intercalative drug-binding [15–17].

CST-2019-101_F3

Figure 3. Ethidium: beta-DNA binding model

Mechanism of Action of Actinomycin D

I have next proposed beta-DNA to be an obligatory intermediate (i.e., a transition-state intermediate) in DNA-melting. This concept readily leads to understanding the mechanism of action of actinomycin D.

Figure 4 (a) shows an electron-micrograph of nucleolar 5-S ribosomal-RNA genes undergoing very active transcription in malignant Hela cells [18] and my interpretation of this process (b) — which indicates the mechanism of action of actinomycin D [19, 20].

CST-2019-101_F4

Figure 4. (a) Electron-micrograph of nucleolar 5-S ribosomal-RNA genes undergoing very active transcription within dividing amphibian oocytes [18]. (b) Interpretation of this photomicrograph showing how actinomycin might act to inhibit this process. Actinomycin binds to beta-DNA, a conformational intermediate that exists within the boundaries connecting double-stranded B-DNA with single-stranded DNA in the transcription complex. This immobilizes (i.e., “pins”) the complex, interfering with the elongation of growing RNA chains.

Actinomycin intercalates into beta-DNA found within the boundaries connecting double-stranded B-DNA with single-stranded DNA in the transcription-complex. This immobilizes (i.e., “pins”) the complex, interfering with the elongation of growing RNA-chains. In extremely active genes such as these, RNA polymerases lie in a close-packed arrangement along DNA. Interference with the movement of one polymerase by actinomycin is expected to inhibit the movement of other polymerases. This predicts nucleolar 5-S ribosomal RNS synthesis to be extremely sensitive to the presence of actinomycin [21–25].

Might actinomycin be expected to preferentially kill malignant cells?

Nucleoli within each nucleus in malignant cells (such as Hela-cells) are expected to contain large numbers of tandem repeats of 5-S ribosomal genes undergoing transcription. Since rapidly dividing malignant cells require increased numbers of ribosomes to carry out protein synthesis, they need many more tandem repeats of 5-S ribosomal genes per nucleoli than normal cells (or alternatively, it is possible that each nucleus within a malignant cell contains many more nucleoli than normal cells, the number of tandem repeats of 5-S ribosomal genes in each nucleolus remaining the same).

The presence of an effect such as this could allow actinomycin to preferentially kill malignant cells. For this reason, trace amounts of actinomycin given over extended periods of time can be expected to be a powerful anticancer chemotherapeutic regimen. Of course, it is first necessary to carry out the appropriate experiments in mice or in other related mammals before attempting clinical use.

References

  1. KIRK JM (1960) The mode of action of actinomycin D. Biochim Biophys Acta 42: 167–169.
  2. GOLDBERG IH, RABINOWITZ M (1962) Actionmycin D inhibition of deoxyribonucleic acid-dependent synthesis of ribonucleic acid. Science 136: 315–316.
  3. GOLDBERG IH, RABINOWITZ M, REICH E (1962) Basis of actinomycin action. I. DNA binding and inhibition of RNA-polymerase synthetic reactions by actinomycin. Proc Natl Acad Sci USA 48: 2094–2101.
  4. Reich E, Goldberg IH (1964) Actinomycin and nucleic acid function. Prog Nucleic Acid Res Mol Biol 3: 183–234.
  5. Sentenac A, Simon EJ, Fromageot P (1968) Initiation of chains by RNA polymerase and the effects of inhibitors studied by a direct filtration technique. Biochim Biophys Acta 161: 299–308.
  6. Perry RP (1963) Selective effects of actinomycin D on the intracellular distribution of RNA synthesis in tissue culture cells.  Exp Cell Res 29: 400–406
  7. Goldberg, IH (1975) Handb Exp Pharmacol. 38: 582–592
  8. Sobell HM, Jain SC, Sakore TD, Nordman CE (1971) Stereochemistry of actinomycin–DNA binding. Nat New Biol 231: 200–205.
  9. Jain, SC, Sobell, HM (1972) Stereochemistry of actinomycin binding to DNA. I. Refinement and further structural details of the actinomycin-deoxyguanosine crystalline complex. J Mol Biol 68: 1–20
  10. Sobell HM, Jain SC (1972) Stereochemistry of actinomycin binding to DNA. II. Detailed molecular model of actinomycin-DNA complex and its implications. J Mol Biol 68: 21–34.
  11. Sobell HM (1974) How actinomycin binds to DNA. Sci Am 231: 82–91.
  12. Sobell HM (1985) Actinomycin and DNA transcription. Proc Natl Acad Sci USA 82: 5328–5331.
  13. Sobell, HM (2009) Premeltons in DNA, Explanatory Publications, Lake Luzerne, NY ISBN 978-0-615-33828-6
  14. Sobell, HM (2013) Organization of DNA in Chromatin, Explanatory Publications, Lake Luzerne, NY ISBN 978-0-692-01974-0
  15. Crothers, DM (1968) Calculation of binding isotherms for heterogeous polymers. Biopolymers 6: 575–583
  16. Wells RD, Larson JE (1970) Studies on the binding of actinomycin D to DNA and DNA model polymers. J Mol Biol 49: 319–342.
  17. Bond PJ, Langridge R, Jennette KW, Lippard SJ (1975) X-ray fiber diffraction evidence for neighbor exclusion binding of a platinum metallointercalation reagent to DNA. Proc Natl Acad Sci USA 72: 4825–4829.
  18. Miller OL Jr, Beatty BR (1969) Visualization of nucleolar genes. Science 164: 955–957.
  19. Sobell HM1 (2016) Premeltons in DNA. J Struct Funct Genomics 17: 17–31.
  20. Sobell, HM (2016) How actinomycin binds to DNA and exerts its mechanism of action. J Struct Funct Genomics
  21. Sobell, HM (2017) World Journal of Pharmaceutical Research 6: 78–84
  22. Sobell, HM (2018) International Journal of Biopharmaceutical Sciences Boffin Access UK  Premeltons in DNA (in press)
  23. Sobell, HM (2018) International Journal of Biopharmaceutical Sciences Boffin Access UK ` Organization of DNA in Chromatin (in press)
  24. Leroy Liu and James Wang have provided a key insight into the nature of DNA supercoiling accompanying transcription that has shed additional light on this question. They have theorized that – in the presence of significant resistance to the rotational motion of the RNA polymerase and its nascent RNA chain around DNA during transcription – the advancing polymerase generates positive superhelicity in the DNA template ahead of it, and negative superhelicity behind it. In nucleolar genes, where there may as many as 200 RNA polymerases moving down the DNA template while synthesizing growing ribosomal RNA-chains, positive and negative superhelical DNA regions between them annihilate one-another, causing adjacent chains to bond-together to form “trains” of transcription complexes, these now moving synchronously along DNA. If this were the case, then the binding by one actinomycin molecule is sufficient to stop the entire “transcription-train” from moving along DNA.
  25. Liu LF1, Wang JC (1987) Supercoiling of the DNA template during transcription. Proc Natl Acad Sci USA 84: 7024–7027.

Osteomyelitis in Children, What to do?

DOI: 10.31038/IJOT.2019211

Introduction

Osteomyelitis is defined as inflammation of the bone with subsequent bone destruction [1]. Osteomyelitis in children has been a diagnostic challenge for decades. Hematogenous osteomyelitis presents clinically versatile, depending on age and causative organism. Historically osteomyelitis has been classified by pathogenesis and duration. Osteomyelitis in children is primary of haematogenous origin [2–5] and was thought to originate from the metaphysis of long bones due to low blood flow in end capillaries [6]. Metaphyseal vessels are open ended towards the physis (growth plate) in long bones. Therefore Stephen RF. et al. [7] suggested the junctions between epiphysis and metaphysis to be the origin of infection in acute osteomyelitis. Time between onset of symptoms and confirmed diagnosis defines acute (within 14 days), subacute (within 3 months) and chronic osteomyelitis (more than 3 months) [8, 9].

Osteomyelitis in children most often originates from long bones, the lower extremities being the most common site, distributed as femur (23–29%), tibia (19–26%) and fibula (4–10%) [5, 9].

Fractures and bony malignancies as Ewing’s sarcoma are important differential diagnoses to rule out during initial assesment. Complications of untreated osteomyelitis in children are septic arthritis, growth retardation due to damage to the physis, bone deformation, angular deformity, septicemia, organ failure and death [5, 10, 11].

Diagnostics

The diagnostic process is multimodal. Both biochemistry and image modalities are required to support the diagnosis sufficiently. Especially imaging modalities has evolved over recent years enabling detailed information on infection status. [3]

A systematic review on acute heamatogenous osteomyelitis [8] showed elevated C-reactive protein (CRP) in 80.5% of children on admission on presentation, 91% had raised Erythrocyte Sedimentation Rate (ESR) and 35.9% had leukocytosis. Pääkönen M. et al evaluated inflammatory markers in 265 children with septic arthritis, osteomyelitis or both. They found sensitivity of ESR and CRP to be 94% and 95% respectively. Raised CRP and ESR simultaneously have a sensitivity of 98% for osteoarticular infection [12]. CRP has short half-life and is inexpensive which makes it useful for monitoring treatment [8, 9].

Subacute osteomyelitis presents with mild symptoms and often no positive laboratory findings [13]. Several studies [13, 14] have presented case series of subacute osteomyelitis in children, with low sensitivities of CRP, ESR and leukocytes.

Standard radiographs have differential diagnostic value, whereas early diagnosis of infection is dependent on more dynamic image modalities. Visible infection on plain radiographs is seen two to three weeks after onset of symptoms [9].

Magnetic resonance imaging (MRI) is the preferred modality for detecting primary focus of infection and it can be helpful in planning surgery. Several studies have documented sensitivity of 82–100% and specificity of 75–99% [10, 15]

Ultra Sonography can reveal subperiosteal changes, evaluate soft tissue and joint effusion and hence have a role in supporting the diagnosis [5].

Bone scans can be useful in multifocal disease and when no clear origin of infection has been found. In neonates specificity is lower due to higher rate of false negative scans [5, 8]. In a recent study, white blood cell scintigraphy was found convincing in detecting post traumatic osteomyelitis [16] Another study reported sensitivity of 79% and a specificity of 97% in detecting fracture related bony infection in peripheral bones [17].

Microbiology and Antibiotics

Staphylococcus Aureus, streptococcus and gram-negative organisms are primary causative organisms in young children [8]. In recent years the facultative anaerobic Gram-negative bacillus Kingella Kingae has been recognized as among the most common pathogen in children between six months and four years of age with haematogenous osteomyelitis [4, 6]. Several published case series includes patients, in which no bacteriological causative organism is identified by microbiological tests on blood, pus or bone biopsy [11, 18].

Antibiotic treatment regimes of acute osteomyelitis in children have historically been 4–6 weeks in total [3, 5]. Two systematic reviews find short intravenous course (3–4 days) followed by 3 weeks of oral therapy as effective as longer intravenous treatment regimes, in uncomplicated cases of acute osteomyelitis.[8, 19] A recent open review concludes that definite guidelines for treatment length and route of administration are yet to be established [4].

Initial empirical treatment, if prevalence of Methicillin Sensitive Staphylococcus Aureus (MSSA) >90%, is in several studies recommended as short course intravenous antistaphylococcal penicillin. Benzylpenicillin/cephalosporin is added if patient is not immunized against Haemophilus Influenzae. When clinical improvement and lowered inflammatory markers treatment is finished by oral regime 3–4 weeks [3–5, 8]. In Methicillin Resistant Staphylococcus Aureus (MRSA) endemic areas, prevalence >10% in community Clindamycin or Vancomycin are recommended [3, 9].

Case Presentation

A nine year old boy was referred to our outpatient clinic from the general practitioner, with plain radiographic findings suspicious of bony malignancy in left proximal fibula (Figure 1). The patient was otherwise healthy and had followed standard Danish children vaccine program.

IJOT - 106_Rikke Thorninger_F1

Figure 1. Primary plain x-­ray when patient was admitted to our department showing isolated process in left proximal fibula.

The boy had 4 weeks lasting constant pain in proximal fibula of the left leg. There had been a minor contusion to the left knee prior to onset of pain. There had been no fewer, swelling or redness of the left leg at any time during period of pain.

Physical examination revealed slim inconspicuously looking legs with no difference between sides. Inspection showed no difference between healthy and affected leg. Palpation of the proximal fibula was painful.

Left leg movement was free and painless, and the neurovascular status was normal. X-ray and subsequently magnetic resonance scanning and computed tomography scan showed large sequester in the left proximal fibula (Figure 3), arising suspicion of sub-acute osteomyelitis. Blood sample inflammatory makers was all within normal range and there were no systemic signs of infection.

The sequester in left fibula was surgically removed with additional decortication above affected bone and clearing of the bone marrow canal. A Gentamicin implant was placed in the bone defect. The pathologic bone was microscopically evaluated and cultured. Cultures were all negative, microscopy showed no signs of acute or chronic inflammation. Morphology was found suspicious of osteofibrous dysplasia. Subsequent multidisciplinary team conference, at Aarhus University Hospital concluded that radiographic material and history suggested subacute osteomyelitis. Osteofibrous dysplasia and malignancy were excluded.

The patient underwent 14 days intravenous treatment postoperatively, with benzylpenicillin and dicloxacillin. After discharge treatment was completed with 4 weeks of per oral dicloxacillin. Amoxycillin/ Clauvulanic acid was primary per oral choice, but was stopped due to allergic skin reaction presenting within five days. Throughout the treatment period all blood tests were within normal range.

At three month follow up, the patient had recovered and had complete remission of pain. Radiographs showed healing in the defect area (Figure 2).

IJOT - 106_Rikke Thorninger_F2

Figure 2. Plain radiographs three months after surgery showed healing around bony defect in proximal fibula.

IJOT - 106_Rikke Thorninger_F3

Figure 3. Magnetic resonance image prior to sugery, showing sequestra in left proximal fibula.

Discussion

In the presented case, a multidisciplinary team discussed paraclinical findings and sustained subacute osteomyelitis as tentative diagnosis, although blood samples had been normal and there were no sign of inflammation in sequestral bone biopsies. In a retrospective study with 121 children diagnosed with acute haematogenous osteomyelitis, sensitivity of blood culture was 32.4% and sensitivity of biopsy culture was 46.6%. In this study 16.5% of patients included underwent surgery [11] Several studies suggest that surgery is reserved for those not responding to standard treatment, but may provide early microbiological diagnose and accelerate patient recovery [3–5]

The patient of this case was offered surgery and subsequently antibiotic therapy. The clinical follow up showed full recovery and remission of symptoms. Inflammatory markers, CRP, leukocytes and ESR was unchanged and within normal range. As pointed out in several studies [3, 4, 8, 20] there are no definite guidelines or clear consensus for assessment and treatment strategy for acute or especially subacute osteomyelitis in children. In this case treatment for suspected osteomyelitis was effectuated without final radiographic or microbiological diagnosis. Total remission of symptoms was seen within the first week after surgery, and there were no clinical sign of infection local or systemic.

The results of this treatment strategy were early recovery and minimal delay due to further diagnostic processes.

Bone abscess representing subacute osteomyelitis can clinically and radiographically mimic bony malignant tumors, being and important differential diagnosis [21]. Dhanoa A. et al. [20] presented six cases, including one child and three adolescents, with initial suspicion of bony malignancy. All cases where diagnosed with subacute osteomyelitis, confirmed by histopathological exam of needle biopsies and microbiological test. The diagnostic process represents a potential delay in relevant treatment of subacute osteomyelitis, and has been described as a clinical and diagnostic challenge [13, 20].

Our patient underwent surgical removal and in total 6 weeks of antibiotic treatment, although no microbiological diagnosis was obtained and biopsy showed no clear sign of osteomyelitis. For subacute osteomyelitis we recommend short antibiotic therapy, initial intravenously and shift to oral therapy guided primarily on clinical remission rather than laboratory testing. In cases with radiographically well defined abscesses, surgically debridement and cleansing of affected bone might support faster recovery and shorter antibiotic treatment regimes.

Further studies are needed to define and test diagnostic algorithms to support clinical decision making and minimize diagnostic delay.

References

  1. de Graaf  H, Sukhtankar P, Arch B, et al (2017) Duration of intravenous antibiotic therapy for children with acute osteomyelitis or septic arthritis: a feasibility study. Health Technol Assess 21: 1–164. [crossref]
  2. Schmitt SK (2017) Osteomyelitis. Infect Dis Clin North Am 31: 325–338. [crossref]
  3. Harik NS, MS Smeltzer (2010) Management of acute hematogenous osteomyelitis in children. Expert Rev Anti Infect Ther 8: 175- 181.
  4. Iliadis AD, Ramachandran M (2017) Paediatric bone and joint infection. EFORT Open Rev 2: 7–12. [crossref]
  5. Yeo A, Ramachandran M (2014) Acute haematogenous osteomyelitis in children. BMJ 348: 66. [crossref]
  6. Jaramillo D, et al (2017) Hematogenous Osteomyelitis in Infants and Children: Imaging of a Changing Disease. Radiology  283: 629–643.
  7. Stephen RF, MK Benson, S Nade (2012) Misconceptions about childhood acute osteomyelitis. J Child Orthop 6: 353–356.
  8. Dartnell J, M Ramachandran, M Katchburian (2012) Haematogenous acute and subacute paediatric osteomyelitis: a systematic review of the literature. J Bone Joint Surg Br 94: 584–595.
  9. Peltola H, Pääkkönen M (2014) Acute osteomyelitis in children. N Engl J Med 370: 352–360. [crossref]
  10. van Schuppen J, MM van Doorn, RR van Rijn (2012) Childhood osteomyelitis: imaging characteristics. Insights Imaging 3: 519–533.
  11. Chiappini E, et al (2017) Epidemiology and Management of Acute Haematogenous Osteomyelitis in a Tertiary Paediatric Center. Int J Environ Res Public Health 14(5).
  12. Paakkonen M, et al (2010) Sensitivity of erythrocyte sedimentation rate and C– reactive protein in childhood bone and joint infections. Clin Orthop Relat Res 468: 861–866.
  13. Spyropoulou V, et al (2016) Primary subacute hematogenous osteomyelitis in children: a clearer bacteriological etiology. J Child Orthop 10: 241–246.
  14. Foster CE, et al (2018) Brodie’s Abscess in Children: A Ten–Year Single Institution Retrospective Review. Pediatr Infect Dis J 2018.
  15. Thevenin–Lemoine C, et al (2016) MRI of acute osteomyelitis in long bones of children: Pathophysiology study. Orthop Traumatol Surg Res 102:  831–837.
  16. Govaert GA, IJpma FF, McNally M, McNally E, et al (2017) Accuracy of diagnostic imaging modalities for peripheral post-traumatic osteomyelitis – a systematic review of the recent literature. Eur J Nucl Med Mol Imaging 44: 1393–1407. [crossref]
  17. Govaert, GAM, et al (2018) High diagnostic accuracy of white blood cell scintigraphy for fracture related infections: Results of a large retrospective single–center study. Injury 49: 1085–1090.
  18. Floyed, RL, RW Steele (2003) Culture negative osteomyelitis. Pediatr Infect Dis J  22: 731–736.
  19. Howard–Jones, AR, D Isaacs (2013) Systematic review of duration and choice of systemic antibiotic therapy for acute haematogenous bacterial osteomyelitis in children. J Paediatr Child Health 49: 760–768.
  20. Dhanoa A, VA Singh (2010) Subacute osteomyelitis masquerading as primary bone sarcoma: report of six cases. Surg Infect (Larchmt) 11: 475–478.
  21. McCarville MB (2009) The child with bone pain: malignancies and mimickers. Cancer Imaging 9: 115–121. [crossref]

Predicting Voice Mutation by Larynx and Voice Modifications

DOI: 10.31038/SRR.2019212

Abstract

Objective: To determine whether the vocal folds length, along with the acoustic voice parameters measurements, can predict the moment of upcoming voice mutation and assess the process of a child’s maturation.

Study design: A cohort study started with examination of children at a premutation age, and a follow up 2.5 and 5 years later.

Setting: Referral center (Claros Otorhinolaryngology Clinic)

Subjects and methods: Children at a premutation age were examined, with a follow up at a mutation and postmutation age. During each visit a CT examination was performed to determine vocal folds length, followed by an examination of the acoustic voice parameters and a videolaryngoscopy and videostroboscopy. Obtained values were analyzed statistically to find the correlations between them and the reported age of mutation.

Results: 50 children (25 males aged 11.5, and 25 females aged 9.5, with a follow up 2.5 and 5 years later) were examined. A study started with 73 children, but 23 of them failed to attend the first or second follow up. Statistical significance was reported for a correlation between the age of mutation and loudness in boys aged 14 (r = 0.48, b = 0.31), vocal folds length in boys aged 14 (b = – 2.18), and loudness in boys aged 11.5 (b = -0.15); and for girls for a correlation between the age of mutation and decrease of fundamental frequency between ages 9.5 and 12 (r = 0.5, b = 0.01).

Conclusion: The parameters mentioned above have a correlation with the moment of mutation and might in future become an additional way of evaluating a child’s development.

Keywords

Voice mutation, Voice break, Vocal fold length, Voice acoustic parameters, Child’s development

Introduction

Proper development in a child can determine his/her educational, professional and emotional future. We watch carefully during child’s growth if this process is not disturbed. However, it is not easy to evaluate; especially during puberty, which is unique to every child and is determined by such an unpredictable and complex factor as the game of hormones [1]. Many authors agree that, while assessing the moment of puberty in girls is easy because of the presence of menarche and breast growth, it is more difficult in boys because of a lack of these concrete breaking moments and its extended character [2, 3]. It is well known that, in contrast, the situation is opposite with respect to mutation; i.e. it happens in a much more subtle way for girls, whilst for boys mutation happens suddenly, more dramatically, and more noticeably [2, 4, 5]. The interesting phenomenon of voice break during puberty has tempted many authors to evaluate the development of a child by assessing the age of voice break as a clean sign of puberty.1, 6, 7 Although exploring the subject of mutation by evaluating the acoustic parameters of the voice has received a reasonable amount of attention in the literature [6, 8–11], the other ways of predicting the time of mutation, especially by examining the length of vocal folds based on CT scans, are to the best of our knowledge underexplored. We conducted our research to address this gap.

Our Clinic is a widely known consultancy for professional opera singers of the Gran Teatro del Liceo in Barcelona, specializing in the issue of voice since 1970, and also providing medical support for a large number of Spanish children. The medical data presented in this article is the result of our work over the past five years. The objective of our research was to establish correlations between the change of vocal folds length and acoustic parameters and signs of voice break described by children and their parents, and therefore to determine which combination of parameters would be the best to evaluate a child’s development.

Changes in the vocal box over the growth of a child are the consequence of complex coordination between the respiratory, digestive and nervous systems [12], as well as anatomical, histological and neurological modifications [13]. In comparison to the adult larynx, the pediatric larynx has disadvantages in voice production on an anatomical level [4], insofar as the ratio of the membranous vocal fold length to the total vocal fold length is lower, the cartilaginous framework is less rigid, and the incidence and degree of posterior glottic chink is increased [14]. Histologically, the pediatric larynx also varies a lot compared to the adult one, which manifests mainly in increased cellularity and decreased cellular differentiation and organization [15], as well as in the lamina propria which begins as a monolayer [16, 17], and changes into a bilayer around the age of 10 and into a trilayer after puberty [18]. Some authors argue that the triple structure of lamina propria occurs already at the age of 7 [4, 19], however, it is widely accepted that the distribution and composition of the collagen and elastic fibers do not mimick those of adults until puberty [15, 16, 18, 20]. Other differences in the pediatric larynx include elevated overall subglottic pressure and recruitment of a greater percentage of pulmonary capacity [4, 21], which changes during mutation.

The reason for these changes lies mainly in the histological structure of the vocal fold – female and male vocal folds alike express androgen receptors in the cytoplasm of the laryngeal gland, progesterone receptors in the nuclei of the same cells, and estrogen receptors in the epithelial cells of the larynx [4, 22, 23], leading to muscle thickening, final development of the trilayered lamina propria, changes in elastin and collagen deposition between the layers, variable lubrication and vocal fold elongation [24]. The expression of the receptors is similar for boys and girls, however differences in the level of hormones between genders causes differences in vocal fold development. In contrast, the other histological features vary: there is more elastin in the cover than in the ligament of male vocal folds, while the elastin in the female lamina propria is more compact [20, 25]. These differences lead to enormous distinctions in the male and female mutations. Apart from the difference in its dynamics, for girls, voice break happens earlier [26–28], starting from the age of 10 and finishing about the age of 14 years, whilst for boys it happens around the age of 12–16 years [29], with some period of voice instability [30, 31]. For both genders it results in the enlargement of arytenoids, expansion of the laryngeal muscles and ligaments [32, 33], completing of the glottal closure, lengthening of the framework of the larynx, and lengthening and rounding of the vocal folds, 4 which is highly related to changes of the acoustic parameters. It is important to note that nowadays mutation occurs much earlier than in the past [2, 34]. This was widely described by Daw, who recorded the age of voice break in members of J. S. Bach’s choirs in Leipzig in 17271749 as being 18 years old [35].

Materials and Methods

The study protocol was acknowledged, reviewed and approved by the internal ethics committee of our medical center, Claros Otorhinolaryngology Clinic Institutional Review Board. All of the parents and children were informed about the examination technique and provided written informed consent. We examined children of a premutation age, with a follow up 2.5 and 5 years later (at the mutation and postmutation age). Exclusion criteria were: vocal fold pathologies, history of neck trauma, previous intubations or laryngeal, head and neck or torso surgeries that have caused changes in vocal folds structure. For the power of a test equal to 0.9 (90%), the smallest sample size was calculated for each checked independent variable, and for statistically significant variables it varied from 10 to 41.

As advised by the pediatric voice assessment guidelines and European Laryngological Society (ELS), subjective and instrumental acoustic evaluations of the voice and aerodynamic performance, as well as visual evaluation of the larynx, were performed [36–38]. During each of the three appointments that the child attended, the voice parameters were measured by a speech- language pathologist. We chose these specific parameters based on advice from the literature: fundamental frequency as the basic, classical objective parameter of the voice [13], vocal range as quite a broad parameter, and because of that a strong sign of a voice disorder if pathological, shimmer and jitter as described as non-invasive, relatively easily applicable and objective [13, 39–41], and, furthermore, highly related to voice problems and dysphonia, [36, 42–44] loudness because it is believed to be a necessary parameter to objectify the result of checked jitter and shimmer [39, 45, 46], and maximum phonation time because it is believed to be the simplest, most easily measured aerodynamic parameter of phonation [47]. To perform the examinations we used sustained vowels taking examples from the approved authors [13, 39, 48, 49]. The vocal recording was performed in accordance with the Union of European Phoniatricians recommendations, with the child in a standing position, in a silent room, with noise level no higher than 40 dB, and with a microphone placed in front of the mouth at a 30 cm distance [50]. We used a microphone from Bruel & Kjaer Rhino-larynx Stroboscope—Type 4914 (Bruel & Kjaer Sound & Vibration, Denmark). All children were examined and recorded in the same conditions. Based on approved literature we defined the norms of all checked parameters [47, 51–59], and we compared them to obtained values. Finally, an ENT consultant examined the vocal folds during every visit to exclude any pathologies, performing a videolaryngoscopy with a rigid endoscope followed by a videostroboscopy (Hopkins II telescope 70 degrees, Karl Storz, Germany).

On every single visit, every 2.5 years, after parents and children provided written informed consent again, CT scanning was performed the way confirmed to be accurate before in our different study [60], using Philips Brilliance ICT 256 (Medical Systems, Netherlands), in the supine position, from the level of the frontal to the level of the aortic arch. Acquisition parameters consisted of a tube current—250 mA, 120 kV, 128×0.625 detector collimation, 0.75-second rotation time, pitch 0.993, scan field of view of250, standard resolution, raw slice thickness – 1 mm. For laryngeal evaluation we added a set of axial reconstruction 2×2 angled through C4 C6 disc spaces. The reconstruction interval was 0.5 mm and the slice thickness was 1 mm. Using standard CT software, a radiologist measured the precise length of the vocal folds in the axial view of the glottis, the longitudinal size of the glottis was estimated in a midsagittal plane (from anterior to a posterior boundary), and in the axial plane, and the length of vocal folds was measured between the anterior commissure and the most posterior part of vocal folds.

Finally, after the third examination (five years after the initial one), the children and their parents answered a survey. The first questions included gender, current age and presumed age of mutation. The following parts of the survey included questions about signs of mutation and voice problems during voice break. The next set of questions related to the age of menarche and breast growth for girls and the age of the first signs of puberty for boys. Lastly, they were asked to complete with the speech-language pathologist the GRBAS scale, which gives scores from 0 to 3 for hoarseness, roughness, breathiness, asthenia, and vocal strain [61].

Data was then implemented into Statistica 13.1 (StatSoft Poland, Cracow) software. Statistical significance was reported at the alpha level of 0.05. P value below 0.05 was considered significant. While analyzing the data we performed the Pearson correlation coefficient test, as well as an analysis of multiple regression and simple linear regression. Correlation coefficients were interpreted to determine whether the effect size was low (correlation coefficient-O.lO), medium (correlation coefficient~0.30) or high (correlation coefficient~0.50). Hypothesis tests were designed as two-tailed. A hypothesis null was formulated as HO: there is no correlation between the change of vocal folds length or acoustic parameters and the moment of voice break (r = 0, b = 0), against the alternative hypothesis H1: there is a correlation (r≠0, b≠0). We created graphs and classification trees to present our findings. The power of the test was determined, and the confidence intervals (Cl) were established for the obtained values.

Results

50 children of a premutation age were our final study group (25 males and 25 females) with a follow up 2.5 and 5 years later (in the mutation and postmutation age). Exactly half of them were males examined at age 11.5, age 14 and age 16.5, and the other half were females examined at age 9.5, 12 and 14.5.

While analyzing the correlations between all the obtained variables and the age of mutation with the Pearson correlation coefficient test, we reported statistical significance in the correlation between the age of mutation and loudness in boys aged 14 (positive correlation coefficient r = 0.48, 0.48, 95%, CI:0, 10–0.73, P = .015, power of the test = .7). Multiple regression analysis showed statistical significance in the correlation between the age of mutation and vocal fold length in boys aged 14 (negative coefficient b = -2.18, P = .044), as well as loudness in boys aged 11.5 (negative coefficient b = -0.15, P = .047), loudness in boys aged 14 (positive coefficient b = 0.31 , P = .022), and loudness in boys aged 16.5 (negative coefficient b = -0.31, P = .040\ however, loudness in boys aged 16.5 cannot be treated as a predictor of mutation, which presumably had occurred earlier). The effect size for these coefficients was R2 = 0.88 (0.88, 95%, CI: 0.72–0.93) and the Cohen’s coefficient was f2 = 7.33. Deeper analysis of these calculations is shown in the classification trees (Figures 1 and 2).

OTO-181055.pdf

Figure 1. Classification tree for loudness in boys aged 11, 5, loudness in boys aged 14, loudness in boys aged 16.5, and vocal folds length in boys aged 14.

Loudness in boys aged 14. and subsequently vocal folds length in boys aged 14, differentiate cases the best. Combined, they are a good prediction of the age of mutation.

SRR Pedro Carlos - 2018-102_F2

Figure 2. Classification tree for all the values obtained in boys aged 11.5 and 14. In these age groups vocal folds length in boys aged 11.5 was the best parameter differentiating cases

In the same calculations for girls, we reported statistical significance in the correlation between the age of mutation and the age of first menstruation (positive coefficient r = 0.84, 0, 84, 95%, CI:0, 66–0, 92, P<.001. power of the test = 1) and increase of fundamental frequency between the age of 9.5 and 12 years old (negative coefficient r = -0.5, 0, 5, 95%, CI:0.13–0.74, P = , 011, power of the test = .75). Since the correlation coefficient between the age of the mutation and the age of the first menstruation was high, we could perform an analysis of simple linear regression, results of which are shown in Figure 3. Multiple regression analysis also showed a statistically significant correlation between the age of mutation and the age of the first menstruation (positive coefficient b = 0.7, P<.001). and increase of fundamental frequency (negative coefficient b = -0.01, P = , 039), which is in agreement with previous calculations. The effect size for these coefficients was R2 = 0.76 (0.88, 95%, CI:0.54–0.86) and the Cohen’s coefficient was f2 = 3.1. Further analysis is shown in the classification trees (Figures 4 and 5).

SRR Pedro Carlos - 2018-102_F3

Figure 3. Scatterplot showing analysis of simple regression and a prediction zone of 95% for the mutation age. We can see that, e.g., for the menstrual age of 12 (axis x), with 95% of probability the mutation will occur between the age of 11.8 and 13.2 (axis y).

SRR Pedro Carlos - 2018-102_F4

Figure 4. Classification tree for the values obtained in girls aged 9.5, 12 and 14.5. Age of the first menstruation was the best parameter differentiating cases in respect of the mutation age.

SRR Pedro Carlos - 2018-102_F5

Figure 5. Classification tree for the increases of the values obtained in girls between ages 9.5 and 12. Age of the first menstruation was the best parameter differentiating cases.

Discussion

The main aim of our study was to find the parameter which correlates the best with the age of mutation, and therefore could possibly serve to predict the age of mutation and evaluate development of the child. We took under further considerations only the values with confirmed statistical significance.

With respect to boys, our calculations showed a statistically significant correlation between the age of mutation and loudness in boys aged 14. vocal fold length in boys aged 14, and loudness in hoys aged 11.5. For loudness in hoys aged 14 the correlation coefficient was positive, which tells us that the louder a child sings at the age of 14, the later he has the mutation. For loudness in boys aged 11.5 and vocal fold length in boys aged 14 the correlation coefficient was negative, which tells us that the louder boy sings at the age of 11.5, the lower the age of mutation, and – most importantly – the longer vocal folds are at the age of 14, the sooner the mutation will start. This is especially interesting in the case of boys, because despite of the obvious vocal folds lengthening with age, male mutation is well accepted to be a sudden and steep change, [2, 4, 5] with periods of higher voice interrupted by periods of lower voice, [30,31] and is dependent on numerous systematic changes described above, therefore it is not simply related to the change of vocal folds length.

It is worth emphasizing, that the correlation between the age of mutation and loudness in hoys aged 14 was confirmed to be statistically significant by all the statistical tests that we performed, and also was the best value differentiating the cases in respect of the age of mutation in the classification tree (Figure 1), therefore it is a variable worth special attention.

Our study also revealed interesting findings in girls. The calculations showed a statistically significant negative correlation between the age of mutation and increase of fundamental frequency between the age of 9.5 and 12. Which means that the more fundamental frequency drops between the age of 9.5 and 12 years old, the later mutation occurs. We have also confirmed a statistically significant positive correlation between the age of mutation and age of the first menstruation. Which means that the sooner the first menstruation appears, the sooner the voice break starts. This is not a surprising result; however, it gave us the opportunity to deepen our statistical analysis. Figure 3 illustrates an analysis of simple regression and the prediction zone of 95% for the mutation age, which means it allows us to predict with 95% of probability the age of voice break knowing the age of menstruation. This way of illustrating the correlation has the potential to be extremely useful in everyday medical and choral practice. The effect size of our results measured by the correlation coefficients was high.

It is important to point out that we have to consider the possible lack of precision in radiological measurements, however it is worth noting that our CT examinations had especially high resolution parameters, and that the CT scans were analyzed multiple times, in different views and planes. A further limitation of our study might be uncertainty about the proper understanding of our instructions during the acoustic examinations (which uncertainty accompanies scientists in every study involving children [62]), however, the age of children involved in our study was not so low as to make this a major concern.

Although several studies have investigated the subject of mutation and ways of predicting it, there is still room to explore it further. Decoster et al. investigated changes in acoustic parameters in girls, however boys were not the subject of the study [2]. Hacki and Heitmiiller, as well as Boltezar et al., did address the subject of mutation, yet in relation to acoustic voice parameters, not vocal folds length [6, 31]. Similarly, numerous studies examined acoustic voice parameters in pediatric population, though other examinations, such as vocal folds length measurements, were beyond the scope of the research [8–11]. Rogers et al. evaluated vocal fold growth as a function of age in a large group of patients [7], however using measuring sticks in total anaesthesia, and emphasized that it might have lengthened the vocal folds [63]. There are several other studies investigating vocal folds length in relation to age, however measurements were performed post mortem, and therefore did not reflect the actual conditions of the living human being’s body [18, 64–67]. Hollien, similarly as in our study, has used radiological imaging; however he has used X-ray images, which are less precise than the CT scans used in our research [68]. Thus, we are tempted to claim that our research is original, and to the best of our knowledge explores aspects not addressed before, adding an important contribution to still not exhausted research about a child’s development and the subject of mutation.

Conclusion

In the academic pursuit of knowledge, evaluating the proper maturation of a child has a special place of a particular concern. This is unsurprising, given that childhood can determine the future of a young human being. However, as we are all different, it is also difficult to determine whether development is proper, and, at the same time, so easy to miss the red flags. Undoubtedly, the period of puberty is the most challenging, both for the human body, which goes through multiple changes, and for scientists, who try to establish reference points to make the evaluation of maturing easier. We believe this hunt is never finished. In the future, one way to asses this might be a routine examination of vocal folds length and acoustic voice parameters. Our study attempted to make our contribution in bringing this future closer.

Acknowledgement

The authors report no conflicts of interest. The authors report no financial and material support for the research and the work reported in the manuscript.

References

  1. Fenichel P (2012) Delayed puberty. Endocr Dev 22: 138–159. [crossref]
  2. Decoster W, Ghesquiere S, Van Steenberge S. Great talent (2008) excellent voices-no problem for pubertal girls? Logoped Phoniatr Vocol 33: 104–112.
  3. Mumm R, Hermanussen M, Scheffler C (2016) New reference centiles for boys’ height, weight and body mass index used voice break as the marker of biological age. Acta Paediatr 105: e459–463. [crossref]
  4. Kerschner JE, Merati AL (2008) Science of Voice Production from Infancy Through Adolescence. In: Hartnick CJ, Boseley ME. Pediatric Voice Disorders. 1st edn. San Diego, California: Plural Publishing Inc; Pg No: 23–29.
  5. McAllister A, Sjölander P (2013) Children’s voice and voice disorders. Semin Speech Lang 34: 71–79. [crossref]
  6. Hacki T, Heitmiiller S (1999) Development of the child’s voice: premutation, mutation. Int J Pediatr Otorhinolaryngol 49: 141–144.
  7. Rogers DJ, Setlur J, Raol N, Maurer R, Hartnick CJ (2014) Evaluation of true vocal fold growth as a function of age. Otolaryngol Head Neck Surg 151: 681–686. [crossref]
  8. Böhme G, Stuchlik G (1995) Voice profiles and standard voice profile of untrained children. J Voice 9: 304–307. [crossref]
  9. McAllister A, Sederholm E, Sundberg J, Gramming P (1994) Relations between voice range profiles and physiological and perceptual voice characteristics in ten-year-old children. J Voice 8: 230–239.
  10. Sorenson DN (1989) A fundamental frequency investigation of children ages 6–10 years old. J Commun Disord 22: 115–123. [crossref]
  11. Pedersen MF, Moller S, Krabbe S, Bennett P, Svenstrup B (1990) Fundamental voice frequency in female puberty measured with electroglottography during continuous speech as secondary sex characteristic. Int J Pediatr Otorhinolaryngol 20: 17–24.
  12. Baken RJ, Orlikoff RF (1997) Voice measurement: is more better. Logop Phoniatr Vocol 22: 147–151.
  13. Nicollas R, Garrel R, Ouaknine M, Giovanni A, Nazarian B, et al. (2008) Normal voice in children between 6 and 12 years of age: database and nonlinear analysis. J Voice 22: 671–675. [crossref]
  14. Pabon JPH, McAllister A, Sederholm E, Sundberg J (2000) Dynamics and voice quality information in the computer phonetograms of children’s voices. In: Wide P, (ed). Child voice Stockholm: KTH Voice Research Centre Pg No: 85–100.
  15. Hartnick CJ, Rehbar R, Prasad V (2005) Development and maturation of the pediatric human vocal fold lamina propria. Laryngoscope 115: 4–15. [crossref]
  16. Hirano M, Sato K (1993) Histological color atlas of the human larynx. San Diego, California: Singular Publishing Group, USA.
  17. Hartnick CJ, Rehbar R, Prasad V (2005) Development and maturation of the pediatric human vocal fold lamina propria. Laryngoscope 115: 4–15. [crossref]
  18. Hirano M, Kurita S, Nakashima T (1983) Growth, development, and aging of human vocal folds. In: Bless DM, Abbs JH, ed. Vocal fold physiology. San Diego, California: College-Hill.
  19. Boseley ME, Hartnick CJ (2006) Development of the human true vocal fold: depth of cell layers and quantifying cell types within the lamina propria. Ann Otol Rhinol Laryngol 115: 784–788. [crossref]
  20. Chan RW, Fu M, Young L, Tirunagari N (2007) Relative contributions of collagen and elastin to elasticity of the vocal fold under tension. Ann Biomed Eng 35: 1471–1483. [crossref]
  21. Stathopoulos ET, Sapienza CM (1997) Developmental changes in laryngeal and respiratory function with variations in sound pressure level. J Speech Lang Hear Res 40: 595–614. [crossref]
  22. Newman SR, Butler J, Hammond EH, Gray SD (2000) Preliminary report on hormone receptors in the human vocal fold. J Voice 14: 72–81. [crossref]
  23. Brunings JW, Schepens JJ, Peutz-Kootstra CJ, Kross KW (2013) The expression of estrogen and progesterone receptors in the human larynx. J Voice 27: 376–380. [crossref]
  24. Kuhn MA (2014) Histological changes in vocal fold growth and aging. Curr Opin Otolaryngol Head Neck Surg 22: 460–465. [crossref]
  25. Kahane JC (1988) Histologic structure and properties of the human vocal folds. Ear Nose Throat J 67: 322, 324–325, 329–30. [crossref]
  26. Marshall WA, Tanner JM. Puberty. In: Falkner F, Tanner JM (1986) Human Growth-A Comprehensive Treatise 2: Postnatal Growth Neurobiology. 2nd ed. Plenum, New York: Springer, Pg No: 171–210.
  27. Tanner JM (1981) Growth and maturation during adolescence. Nutr Rev 39: 43–55. [crossref]
  28. Bogin B (1999) Patterns of human growth. (2nd edn) Cambridge University Press, London.
  29. Balasubramaniam RK, NN (2017) Voice Mutation During Adolescence in Mangalore, India: Implications for the Assessment and Management of Mutational Voice Disorders. J Voice 31: 511.e29–511.e33.
  30. Fuchs M, Fröehlich M, Hentschel B, Stuermer IW, Kruse E, et al. (2007) Predicting mutational change in the speaking voice of boys. J Voice 21: 169–178. [crossref]
  31. Boltezar IH, Burger ZR, Zargi M (1997) Instability of voice in adolescence: pathologic condition or normal developmental variation? J Pediatr 130: 185–190. [crossref]
  32. Fitch WT, Giedd J (1999) Morphology and development of the human vocal tract: a study using magnetic resonance imaging. J Acoust Soc Am 106: 1511–1522. [crossref]
  33. Ferlito A (2000) In: Ferlito A., ed. Diseases of the Larynx. London: Oxford; New York: Arnold. Oxford Univeristy Press.
  34. Gackle L (2000) Understanding voice transformation in female adolescents. In: Welch G, ed. Bodymind and Voice. Minnesota, USA: The Voice Care Network; Pg No: 739–744.
  35. Daw SF (1970) Age of boys’ puberty in Leipzig, 1727–49, as indicated by voice breaking in J.S. Bach’s choir members. Hum Biol 42: 87–89. [crossref]
  36. Cohen W, Wynne DM, KubbaH, McCartney E (2012) Development of a minimum protocol for assessment in the paediatric voice clinic. Part 1: Evaluating vocal function. Logoped Phoniatr Vocol 37: 33–38.
  37. Kelchner LN, Brehm SB, de Alarcon A, Weinrich B (2012) Update on pediatric voice and airway disorders: assessment and care. Curr Opin Otolaryngol Head Neck Surg 20: 160–164. [crossref]
  38. Hartnick C, Boseley M (2008) Pediatric Voice Disorders: Diagnosis and Treatment. San Diego, CA: Plural Publishing.
  39. Brockmann-Bauser M, Beyer D, Bohlender JE (2014) Clinical relevance of speaking voice intensity effects on acoustic jitter and shimmer in children between 5;0 and 9; 11 years. Int J Pediatr. Otorhinolaryngol 78: 2121–2126.
  40. Oates J, Kirkby R (1980) An acoustical investigation of voice quality disorders in children with vocal nodules. Australian J HumanCommDis 8: 28–39.
  41. Meal lister A, Sundberg J, Hibi SR (1998) Acoustic measurements and perceptual evaluation of hoarseness in children’s voices. Logoped Phoniatr Vocol 23: 27–38.
  42. Lopes LW, Barbosa Lima IL, Alves Almeida LN, Cavalcante DP, de Almeida AA (2012) Severity of voice disorders in children: correlations between perceptual and acoustic data. J Voice 26: 819.e7–819.el2.
  43. Shah RK, Engel SH, Choi SS (2008) Relationship between voice quality and vocal nodule size. Otolaryngol Head Neck Surg 139: 723–726. [crossref]
  44. Valadez V, Ysunza A, Ocharan-Hernandez E, Garrido-Bustamante N, Sanchez-Valerio A, et al, (2012) Voice parameters and videonasolaryngoscopy in children with vocal nodules: a longitudinal study, before and after voice therapy. Int J Pediatr Otorhinolaryngol 76: 1361–1365.
  45. Brockmann-Bauser M (2012) Improving Jitter and Shimmer Measurements in Normal Voices. Schulz-Kirchner Verlag GmbH Idstein 205.
  46. Brockmann-Bauser M, Drinnan MJ (2011) Routine acoustic voice analysis: time to think again? Curr Opin Otolaryngol Head Neck Surg 19: 165–170. [crossref]
  47. Pribuisiene R, Uloza V, Kardisiene V (2011) Voice characteristics of children aged between 6 and 13 years: Impact of age, gender, and vocal training. Logopedics Phoniatrics Vocology 36: 150–155.
  48. Glaze L, Bless D, Milenkovic P, Susser R (1988) Acoustic characteristics of children’s voice. J Voice 2: 312–319.
  49. Cohen SR, Perelman N, Nimni ME, Mahnovski V, Cheung DT (1993) Whole organ evaluation of collagen in the developing human larynx structures and adjoining anatomic structures (hyoid and trachea). Ann Otol Rhinol Laryngol 102: 655–659.
  50. Shutte HK, Seidner W (1983) Recommendations by the Union of European Phoniatricians (UEP): Standardizing voice area measurement/phonetography. Folia Phoniatr 35: 286–288.
  51. Bennett S (1983) A 3-year longitudinal study of school-aged children’s fundamental frequencies. J Speech Hear Res 26: 137–141. [crossref]
  52. Stathopoulos ET, Sapienza C (1993) Respiratory and laryngeal measures of children during vocal intensity variation. J Acoust Soc Am 94: 2531–2543. [crossref]
  53. Finnegan DE (1985) Maximum phonation time for children with normal voices. Folia Phoniatr (Basel) 37: 209–215. [crossref]
  54. Horii Y (1983) Automatic analysis of voice fundamental frequency and intensity using a VisiPitch. J Speech Hear Res 25: 467–471.
  55. Banik A, Arya S, Kant A (2015) Vocal Parameters in Children between 4 To 12 Years of Age: An Attempt to Establish a Prototype Database. International Journal of Scientific and Research Publications 11: 446–453.
  56. Ting HN, Chia SY, Manap HH, Ho AH, Tiu KY, Abdul Hamid B (2012) Fundamental Frequency and Perturbation Measures of Sustained Vowels in Malaysian Malay Children Between 7 and 12 Years Old. J Voice. 26: 425–30.
  57. Bonet M, Casan P (1994) Evaluation of dysphonia in a children’s choir. Folia Phoniatr Logop 46: 27–34. [crossref]
  58. Williamson G (2006) Human communication: A linguistic introduction. 2nd ed. Billingham, UK: Speech Language Services.
  59. Maturo S, Hill C, Bunting G, Ballif C, Maurer R, et al. (2012) Establishment of a normative pediatric acoustic database. Arch Otolaryngol Head Neck Surg 138: 956–961. [crossref]
  60. Claros P, Sobolewska AZ, Domenech-Claros A, Claros-Pujol A, Pujol C, et al. (2018) CT-based Morphometric Analysis of Professional Opera Singers’ Vocal Folds. J Voice 21: 30614–30618.
  61. Sundberg J (2000) Emotive transforms. Phonetica 57: 95–112. [crossref]
  62. Brockmann-Bauser M, Beyer D, Bohlender JE (2014) Clinical relevance of speaking voice intensity effects on acoustic jitter and shimmer in children between 5;0 and 9; 11 years. Int J Pediatr Otorhinolaryngol 78: 2121–2126.
  63. Su MC, Yeh TH, Tan CT, Lin CD, Linne OC, et al. (2002) Measurement of adult vocal fold length. J Laryngol Otol 116: 447–449. [[crossref]
  64. Eckel HE, Sprinzl GM, Sittel C, Koebke J, Damm M, et al. (2000) Anatomy of the glottis and subglottis in the pediatric larynx]. HNO 48: 501–507. [crossref]
  65. Eckel HE, Sittel C (1995) Morphometry of the larynx in horizontal sections. Am J Otolaryngol 16: 40–48. [[crossref]
  66. Friedrich G, Kainz J (1988) Morphometry of the larynx in horizontal sections: normal data for the quantitative evaluation of current imaging technics in German], Laryngol Rhinol Otol (Stuttg) 67: 269–274.
  67. Kahane JC (1978) A morphological study of the human prepubertal and pubertal larynx. Am J Anat 151: 11–19. [crossref]
  68. HOLLIEN H (1960) Some laryngeal correlates of vocal pitch. J Speech Hear Res 3: 52–58. [crossref]

Prevalence, Symptoms and Treatment of Vocal Fatigue in Professional Opera Singers: Clinician’s Experiences

DOI: 10.31038/SRR.2019211

Summary

Purpose and study design: The research aims to estimate the prevalence and characteristic symptoms of vocal fatigue in professional opera singers. Also, the paper is summary of many years of clinical experience with management of singers with vocal fatigue in Clarós Clinic. The research was designed as a retrospective observational study.

Material and method: In the group of 250 professional opera singers who were examined in Clarós Clinic in 10 years period, fifty-five cases of vocal fatigue were reported and evaluated. Among subjects were 21 men and 34 women. Mean age of participants was 46,78y (range 19–72 years old, standard deviation: 12.18 year). Representation of classical voice types was: 14 tenors, 5 baritones, 2 basses, 22 sopranos, 10 mezzos, 2 contralti.

Results: Prevalence of vocal fatigue in 10 years period in the study group was 22%. The three most frequent symptoms observed in the study group were: muscle pain (87.27%), muscle fatigue (76.36%) and diffuse sharp pain in the neck (70.37%). Statistical analysis showed significance only for the relationship between female opera singers and incidence of symptoms such as tremulous voice and muscle pain. All other symptoms were not statistically related to gender. Contradictory to expectations, application of anti-inflammatory drugs was statistically associated with the longer duration of the symptoms.

Conclusions: Vocal fatigue may be underestimated but is one of the most frequent problems encountered in the ENT practice that provides care for professional voice users. In the presented research, administration of anti-inflammatory drugs has not been associated with faster recovery from vocal fatigue.

Keywords

vocal fatigue, opera singers, treatment of vocal fatigue, retrospective observational study, the prevalence of vocal fatigue.

Introduction

Professional opera singers are considered to be the highest class artists among singers. In fact, they should be perceived as Olympian athletes when comes to quality and ability of their vocal folds. Exceptional singing is achieved by years of training and vocal regime. Unfortunately, the requirements set for the opera singers’ voice entail the unique susceptibility to vocal fatigue. The research reported fatigue to be a common complaint among professional voice users [1].

Vocal fatigue presents a challenge to research and clinical practice. Despite the amount data gathered on that matter, precise definition and guidelines are still not well established.

Unlike in case of muscle fatigue, the creation of the conditions that allow the investigation be far more difficult because of complicated mechanism of voice production. Many studies which attempted to induce vocal fatigue yielded a different and inconsistent result. In case of singing, the task is especially challenging because repertoire can vary significantly. Furthermore, aspect like frequency of the performance can play a role in the development of vocal fatigue. Some types of voice are particularly vulnerable, like a soprano. However, the data provided on the subject is anecdotal, obtained from clinicians’ experience instead of large case-control studies.

Mechanisms which underlay vocal fatigue are multifaceted. Titze listed some potential physiological and biomechanical factors that may contribute: fatigue of respiratory and laryngeal muscles, fatigue of non-muscular tissues of the larynx, and changes in vocal folds’ viscosity [2]. As one of the others, the neuromuscular fatigue express inability of muscles to sustain the tension under repeated stimulation [3].

The exceptional ability of the larynx to produce sound is possible because of vocal folds which are covered with non-muscular, pliable tissue that generates multiple and rapid vibrations [4]. Furthermore, viscous properties of vocal folds epithelium allow lubrication and shock absorption [5]. Research showed that prolonged, high-pitched phonation could increase: frictional energy loss, heat dissipation and tissue viscosity. All these factors can lead to tissue fatigue [4]. The influence of tissue biomechanics on phonation makes the study about vocal fatigue more complicated than the research of fatigue involving other skeletal muscles.

The expression vocal fatigue has been applied and arbitrarily understood. For presented research, vocal fatigue definition was employed from Solomon as the self- report of an increased sense of effort with prolonged phonation, whether or not there are observable or measurable decrements in phonation function [5]. The aspect of self-reporting is particularly crucial amid opera singers who are the harshest judges of their vocal abilities.

Moreover, vocal fatigue has been described clinically by several authors. Common symptoms enlisted in literature were: husky vocal quality, breathy vocal quality, loss of voice, pitch breaks, inability to maintain regular pitch, reduce pitch range, lack of vocal carrying power, reduce loudness range, need to use greater vocal effort [1].

Opera singers are interesting subjects to study vocal fatigue because in comparison to other professional types of vocalists their vocal abilities are challenged at the highest level. Furthermore, years of vocal regime provide them with high self-awareness of the vocal capabilities which helps to notice any vocal problems at first onset.

Clarós Clinic provides medical support for a large number of professional opera singers of Liceu Opera theatre in Barcelona, Music Conservatory students and other singers since 1970’. Authors decided to analyse reports of vocal fatigue from last ten years among opera singers and describe common symptoms, frequency and interventions that were applied.

The research aims to determine the prevalence of vocal fatigue, symptoms and factors which may underlay the occurrence of that problem.

Paper is summary of our experience with treatment and management of vocal fatigue among professional opera singers.

Overall, vocal fatigue is still an intriguing and persistent problem when presented in clinical practice and is the challenge in the face of consequences which may cause for the professional singer with a tight schedule. Furthermore, the struggle with finding a substitute for a famous singer who suffers from vocal fatigue provides additional pressure for the clinician who has to offer effective treatment.

Materials and Methods

Study design

The research was designed as a retrospective observational study. Evaluation of vocal fatigue cases presented in Clarós Clinic during ten years period was conducted to obtain prevalence, symptoms and to analyse the treatment.

The research protocol was approved by ethics committee of Clarós Clinic medical centre.

Participants

In the group of 250 professional opera singers who were treated in Clarós Clinic in 10 years, fifty-five cases of vocal fatigue were reported and evaluated for this research. Among subjects were 21 men and 34 women. Mean age of participants was 46,78y (range 19–72 years old, standard deviation: 12.18 year). Representation of voice types was: 14 tenors, 5 baritones, 2 basses, 22 sopranos, 10 mezzos, 2 contralti.

Medical evidence

Eleven most common symptoms and six risk factors of vocal fatigue identified in the study group were gathered in Table 1.

Table 1. Common symptoms of vocal fatigue and risk factors with definitions.

Symptoms

Description

Hoarseness

Patient complains of hoarse voice

Breathy voice

Breathy vocal quality, running out of breath while talking

Whispery voice

Patient is only able to whisper

Tremulous voice

Unsteady voice, voice affected by trembling or tremors

Neck muscle pain

Patient complains of muscular pain in the neck

Sharp pain localised on the neck: diffused or localised

Patient complains of sharp pain in the neck during speaking or singing

Neck muscle fatigue

Muscular fatigue present on the neck

Tissue fatigue

Increased viscosity of vocal folds’ mucosa and tissue stiffness

Neck strain

Increased tension in the neck

Stiffness of the vocal folds

Increased tension of vocal folds

Changes in vibrato

Inability to maintain proper vibrato

Risk factors

Muscle overstrained

Excessive effort during performance

Overuse of voice

Inadequate vocal rest regime

Incorrect technique

Insufficient training

Singing warm-up

Inappropriate singing warm-up

Shouting

High-pitched, forced phonation in short period

Inadequate repertoire

Prolonged singing beyond the appropriate tessitura (most acceptable and comfortable vocal range for the given singer [28])

Medical records were evaluated to search for possible causes and factor which may contribute to the development of vocal fatigue. Patients reported: overuse of voice, stress, incorrect technique, vocal warm-up, improper repertoire and shouting. Interventions which were used: voice rest regime, medications (anti-inflammatory drugs, steroids: hydrocortisone). Also, duration of the symptoms was measured. In some cases, the repetitive character of vocal fatigue was noticed.

 All cases of professional opera singers were presented to the same, senior, most experienced ENT consultant. He was responsible for the patient examination and evaluation. All medical records of patients with vocal fatigue were created by the senior consultant who followed the similar protocol in every case. Table 1 gathers the common symptoms and risk factors assessed in the study with the description of researcher interpretation.

In every case, the standard medical interview was gathered. Moreover, patients underwent ENT examination consisted of endoscopic evaluation and neck palpation.

A senior most experienced ENT consultant performed video laryngoscopy with conventional equipment to examine the vocal folds (Karl Storz® 70 degrees rigid endoscope and HD camera).

The table below shows typical symptoms and risk factors which were assessed in medical records. Description represents how researcher stated or understood given symptom or risk factor. Definitions were based on literature and authors experienced [3,6].

Statistical Analysis

Data was collected in Excel sheet and implemented to Statistica 13.1 (Statsoft) software for statistical analysis. Statistical significance was accepted at the alpha level of 0.05. A p-value below 0.05 was considered significant.

Contingency tables were used to analyse quality variable obtained from patients’ medical history. Percentage values for every symptom and risk factor were calculated for men, women and the complete group. The Chi2 test was used to assess statistical dependence between gender and incidence of given symptom and risk factor.

Contingency tables were also implemented to analyse the possible influence of anti-inflammatory medication on the duration of symptoms.

Due to the qualitative nature of the data obtained from medical history more complicated analyses were not possible.

It needs to be to highlight that some types of voice are very infrequent which influence the precision of statistical analysis. In this study, the contralti, basses and countertenors occurred in the minority.

Results

The results section is divided into two parts: descriptive data (tables and chart) and the statistical analysis. The precise characteristic of descriptive data obtained from medical history is presented in table 2, 3 and 4. Percentage values of symptoms incidence are illustrated in Chart 1.

SRR-2019-101-Pedro Carlos_F1

Chart 1. Bar chart. Percentage value of symptoms’ incidence.

Table 2.
Common symptoms of vocal fatigue in the study group.

Whispery voice

Tremulous voice

Hoarseness

Muscle pain

Sharp pain

Muscle fatigue

Tissue fatigue

Vocal folds stiffness

Vibrato changes

Pitch changes

Diffuse

Localised

Men

18.18%

5.45%

21.82%

38.18%

23.63%

7.41%

32.73%

12.73%

23.64%

9.09%

9.09%

Women

45.45%

25.45%

40.00%

49.09%

40.74%

22.22%

43.64%

20.00%

27.27%

43.64%

12.73%

General

63.64%

30.91%

61.82%

87.27%

70.37%

29.63%

76.36%

32.73%

50.91%

52.73%

21.82%

Chi2*

p = 0.052

p = 0.037

p = 0.574

p = 0.026

p = 0.23

p = 0.19

p = 0.93

p = 0.19

p = 0.779

p = 0.000

*The Chi2 for Independence test. Significant values (p≤0.05). NOTE: VF- vocal fatigue.

The table shows the percentage of VF occurrence and percentage of cases in which medications were applied in the groups. The results of the Chi2 test for independence between enlisted factor and gender are also given. Significance was reported at p-level below 0.05. None of given factors had gender predilection in the presented group.

Table 3. Risk factors of vocal fatigue analysed in the study group.

Muscle overstrain

Overuse of voice

Incorrect technique

Singing warm-up

Shouting

Inadequate repertoire

Reoccurrence VF

Medication

Men

32.73%

30.91%

3.64%

18.18%

20.00%

10.91%

7.27%

10.91%

Women

49.09%

49.09%

10.91%

21.82%

29.09%

18.18%

23.65%

14.55%

General

81.82%

80.00%

14.55%

40%

49.09%

29.09%

30.91%

25.24%

Chi2*

p = 0.55

p = 0.889

p = 0.406

p = 0.34

p = 0.7

p = 0.946

p = 0.13

p = 0.676

*The Chi2 for Independence test. Significant values (p ≤ 0.05). NOTE: VF- vocal fatigue.

The table shows the percentage value of risk factors incidence in the groups and results of the Chi2 test for independence between enlisted factor and gender. Significance was reported at p-level below 0.05. None of given risk factors had gender predilection in the presented group.

Table 4. Percentage value of vocal fatigue reoccurrence.  Percentage of cases in which anti-inflammatory medication was applied.

Reoccurrence VF

Medication

Men

7.27%

10.91%

Women

23.65%

14.55%

General

30.91%

25.24%

Chi2*

p=0.13

p=0.676

*The Chi2 for Independence test. Significant values (p≤0.05). NOTE: VF- vocal fatigue.

The Chi2 for Independence test was used to assess the relationship between gender and vocal fatigue symptoms. Results of the Chi2 test are also given in tables 2, 3.

Table 5 presents the percentage values of vocal fatigue recurrence among different types of classical voices.

Table 5. Recurrence of vocal fatigue among different voice types (tessitura- defined as most acceptable and comfortable vocal range for the given singer [28]).

Recurrence of vocal fatigue in different voice types

Yes

No

General

Soprano

12.73%

27.27%

40.00%

Mezzosoprano

10.91%

7.27%

18.18%

Contralto

0,00%

3.64%

3.64%

Tenor

5.45%

20.00%

25.45%

Baritone

1.82%

7.27%

9.09%

Bass

0.00%

3.64%

3.64%

General

30.91%

69.09%

100.00%

The table shows the percentage value of risk factors incidence in the groups and results of the Chi2 test for independence between enlisted factor and gender. Significance was reported at p-level below 0.05. None of given risk factors had gender predilection in the presented group.

The table presents the percentage of vocal fatigue reoccurrence among different voice types. In presented data sopranos and tenors were most frequently affected by the reoccurrence of VF. Examined contralti and basses had not experienced vocal fatigue more than one time at the moment of evaluation.

The chart shows the percentage value of symptoms which occurred in the whole group (general), men and women groups. Bars help to illustrate which symptoms were most common and compared them between groups. Therefore, three most frequent symptoms were muscle fatigue, diffuse sharp pain and muscle pain. Most common symptom among women as well as in men was muscle pain.

Results – Summary

The Clarós Clinic provided medical support for 250 professional opera singers during the time that the data was gathered. Symptoms of vocal fatigue were reported in 55 operatic vocalists, and these cases were enlisted to the research.

Prevalence of vocal fatigue in 10 years period among opera singers examined in the Clinic was 22%.

The three most frequent symptoms observed in the study group were: muscle pain (87.27%), muscle fatigue (76.36%) and diffuse sharp pain (70.37%). Most frequent complaints were pain-related, especially in female singers group.

Most common vocal complaints were: whispery voice (63.64%) and hoarseness (61.82%). More than a half of opera singers (52.73%) had difficulties with maintaining the vocal pitch during singing and quarter (21.82%) noticed changes in vibrato

The Chi2 test showed statistical significance only for the relationship between female opera singers and incidence of symptoms such as tremulous voice and muscle pain (test Chi2: p = 0.037, p = 0.026). All other symptoms were not statistically related to gender.

Most common risk factors were: muscle overstain and overuse of voice which were present in over 80% of cases. Much less frequently, singers reported incorrect technique and inadequate repertoire as contributors to vocal fatigue (table 3).

Recurrence of vocal fatigue was noted in one-third of the singers and was more distinctive for female singers, especially sopranos and mezzo-sopranos (table 4).

Mean duration of the symptoms was 3.8day (standard deviation: +/- 1.9day). The median value was 3 day.

Contingency tables were also used to estimate the relationship between duration of the symptoms and application of anti-inflammatory medications. Results of the Chi2 test showed that shorter length of the symptoms was related to lack of administration medication (test Chi2: p = 0.009).

Discussion

Vocal fatigue is an interesting and often debilitating condition, affecting many professional voice users. It gained much attention in the field of research, yet mechanisms which underlie the onset of this condition and its pathophysiology are still not fully understood. The amounts of vocal effort and specific elements which can trigger vocal fatigue are part of the ongoing debate. More research has to be done to develop reliable guidelines for management and treatment of vocal fatigue.

Every group of professional voice users have its characteristic which helps to study individual exposure factors. Opera singers have the individual susceptibility to fatigue which may interfere with social and occupational functioning.

The purpose of this study on vocal fatigue in opera singers was to state prevalence, characterise symptoms of this condition and review management.

Prevalence of vocal fatigue in studied group was 22%, which suggests that it might concern every fifth singer. As previous research showed among other types of professional singers, VF can cause voice impairment even more often. In a study conducted on the large group of various kinds of singers (opera singers- 49.8%), Phyland reported that participants experienced vocal fatigue in the previous year in 69% of cases [7]. In the research, fatigue was the second most frequent problem reported by singers after hoarseness [7].

Among most common symptoms reported by opera singers in the study group: two were pain related. Over 87% of singers pointed out the muscle pain as a single most common symptom of vocal fatigue. Female opera singers tended to suffer more from muscle pain than male singers, which was also confirmed statistically significant. Muscle pain was usually localised in throat, jaw and neck, but also in chest and back. These findings were consistent with previous reports which stated that most common pain present in singers were a sore throat (66%), pain during speaking (41%) and neck pain (35%). However, the study mentioned above pointed out the tendency for a sore throat among male vocalists, other types of body pains had no difference according to gender [8]. The presence of pain can severely compromise singer’s performance and negatively influence the quality of life. An important factor which helps to prevent the muscle pain is the concern for proper technique and vocal rest regimen.

Neuromuscular fatigue has been widely investigated in the literature. It can be presumed that muscle of the respiratory and phonatory system can fatigue and contribute to the deterioration of phonation or the perception of increased vocal effort, especially during prolonged high-pitched phonation. Undoubtedly, the research showed that respiratory muscles are highly unlikely to experienced fatigue. More recent findings presented evidence from whole body exercise suggesting that respiratory muscle fatigue occurs only following constant high-intensity training [9]. This situation cannot occur during regular physical activity, even as challenging as prolonged singing.

The distinction between fatigue of skeletal muscle and phonation muscles is relevant, because of the different histological structure. The capability of a muscle to maintain contraction over an extended period is related to a distribution of different motor units within the muscle body. In case of the larynx, the vast majority of intrinsic laryngeal muscle have fatigue-resistant muscle fibres (type I and IIa) rather than fatigable (type IIb) [10]. More recent studies provide interesting data showing the even more complicated histological structure of human intrinsic laryngeal muscles than presented in animal models [10]. These facts help to explain why singers usually complain about fatigue of muscle and experience discomfort in areas primarily localised in throat, jaw, and neck.

One of the unique aspects is non-muscular tissue fatigue which represents mechanical exhaustion. Mechanical deterioration represents the amount of strain that material can tolerate before breaking down. The fatigue represents progressive structural damage that results from mechanical stress (force per unit) imposed by strain on the material. Titze reported in one of his research that non-muscular tissue fatigue could cause damage to the laryngeal mucosa, but the quantity and duration of the physical stress were uncertain [11]. The author also described tensile stress which is required for high pitch phonation as a most significant mechanical stress in vocal folds vibration [11]

 Tissue viscosity plays a role in response to mechanical stress because that feature refers to individual properties of vocal folds’ mucosa responsible for lubrication and shock absorption. Research demonstrated that viscosity highly depends on the systemic and superficial hydration of mucosa. Singers in a situation of reduced systemic hydration may be particularly prone to experience the vocal fatigue [12]. Factor as prolonged, high-pitched phonation without proper hydration can lead to stress and strain which placed on the tissue can provoke fatigue. Furthermore, the viscosity of vocal folds mucosa can be affected by the decreased humidity of environment as in case of oral breathing [13].

The more recent study confirmed positive influence of systemic hydration on perceptual parameters of voice quality in singers. The improvement was seen in higher fundamental frequency, less cycle to cycle variation in pitch, or longer phonation time, depending on the individual. Hydrated vocal folds allow for optimal vibration, increase ease of phonation and prevent structural damage to vocal to vocal folds mucosa [14,15,16].

Non-muscular tissue biomechanical properties which include mucosal viscosity plays a significant role in the development of fatigue. In presented study 1/3 of singers have suffered from problems related to tissue fatigue.

Further, changes in vibrato were observed in over 52% singers. These changes represent acoustic differences related to vocal fatigue and can seriously interfere with performance. In previous research, Titze noted that muscle fatigue results in decreased ability to maintain stable tension in vocal folds [2]. In another interesting study, Boucher attempted to isolate acoustic signs of fatigue in laryngeal muscles. The study showed that 12 conventional acoustic parameters that were measured neither demonstrated consistent linear relationships with the fatigue estimates. Though, average values demonstrated the consistent peaks in vocal tremor and appeared near the points of critical shifts in muscle fatigue. In sum, rises in tremor corresponding to shifts in muscle fatigue appear robust in the face of fluctuations in modal pitch [17].

A more recent study helped to differentiate the acoustic changes related to vocal fatigue. The research was the first to demonstrate a link between tremor and observed muscle fatigue that is specifically attributable to voice effort and not only to fatigue linked with waking hours. Also, the results do not support applications of F0 or other conventional acoustic parameters as manifestations of fatigue in laryngeal structures. Even though many research showed significant rises in F0 as a result of vocal effort with reference to group averages, the recent reports showed that individual or cross-subject changes in F0, as in other conventional acoustic parameters do not consistently indicate fatigue in laryngeal structures [18].

Despite years of research, no consensus has emerged that could support the elaboration of guidelines for vocal fatigue. In the most basic approach, investigations on vocal fatigue have concentrated on identifying changes in voice in tests where fatigue was provoked by tasks of various “vocal load”. The assignments varied across different research; participants were asked to read or sing at varying pitch or intensities for a variable period extending from few minutes to several hours [19]. Those arguments make cross-study comparison useless with results on suggested vocal symptoms inconsistent and sometimes contradictory.

In our study as the first line of treatment, the vocal rest regime was applied in every case. Importance of vocal rest was underlined in many research and is usually required as first line intervention when vocal fatigue is experienced by singer [20]. Stress and anxiety management is also crucial for maintaining the good psychological condition of a singer. The aim of physical and mental approach to prevention of vocal fatigue is optimisation of the performance efficacy [21]. This translate to minimising muscular activation, achievable by improving posture and relaxing muscles [22].

Professional opera singers often follow the vocal routine which usually begins under the influence of their singing teacher and speech-language pathologist. Vocal hygiene practices contain moderation in amount and type of voice use, reduction of stress, avoidance of phonotraumatic behaviours like shouting, talking over crowds, aspects like systemic hydration and humidification to improve performance and ensure voice longevity [23]. Effects of systemic hydration and vocal rest were proved to have a significant influence on the decrease of vocal fatigue and maintaining good vocal quality in general [16,24].

In literature, researchers postulated to set safety limits of vibration dose (phonatory time) for professional voice users to protect people in several occupations (singers, teachers). For instance, proper recovery time has been worked out for professional athletes who abuse their body in different ways. Importantly, Titze divided recovery for short- and long-term. The first one takes place immediately when phonation is stopped [25]. The primary benefit from short recovery is for the muscles whose chemicals get reset before next contractions. On the contrary, traumatised epithelial cells need more prolonged healing. Some of them after being heavily bombarded during vocal folds contraction can degenerate and be shed off. New cells will grow underneath, but that requires time. Furthermore, some destruction of the structural matrix of the lamina propria may be present after prolonged phonation. Fibroblasts activity is necessary for the repairing to continue constantly. The recovery process may range from several hours to 72h to complete [25]. As was presented in the study, the actual phonation for opera singers in 2–3h opera was of the series 20–30min for leading role, and their schedule performance was on the order of 3 per week [25]. Overall, these facts put opera singers in the favourable position for proper vocal recovery. Nevertheless, the type and loudness of phonation were not considered in research calculation, but they can play a tremendous role. Given this points, more research is needed to create appropriate guidelines to prevent singers from vocal fatigue.

Singers appear to be at particular risk of developing voice problems. Formal assessments of singers experiencing a voice problem at any given moment in a time range from about 20 to 50% [26]. Moreover, the impact of voice problems on quality of life was widely investigated in the literature. Many studies proved that decrease in voice quality and other voice impairments affects severely quality of life. Of course, in case of professional singers, this problem grows to the crucial role because directly concerns the source of income.

Attempts at analysing vocal fatigue in patients who already experienced this condition are difficult because of subject heterogeneity and burden with problems of data interpretation. Even individuals selected for having only symptoms of vocal fatigue usually present variable baseline and outcome data.

Professional opera singers have high stakes in sustaining excellent vocal condition but also experienced unique vocal demands, making them important population to study.

Useful tools to adapt to everyday practice are scales which help to evaluate singers’ perceptions of physical aspects of singing status. A good example is EASE scale (Evaluation of the Ability to Sing Easily) which is clinical outcome test for symptomatic aspects of compromised vocal health but was not designed as the primarily disease-specific instrument [27]. This test serves as the measure of potential changes in the singing voice which may indicate effects of vocal effort and may help to detect singers with the increased risk of possible development of voice disorders. The EASE was designed to help singers in assessing vocal load threshold, recovery time to assist performance scheduling, help to predict the development of vocal problems, evaluate therapeutic outcome in management for specific needs of the singer’s voice, lastly to provide supportive data for determining performance fitness [27].

Conclusion

Vocal fatigue may be underestimated but is one of the most frequent problems encountered in the ENT practice which provides care for professional voice users. The most frequent symptoms were muscle fatigue, diffuse sharp pain and muscle pain. In the study group, administration of anti-inflammatory drugs has not been associated with faster recovery from vocal fatigue.

References

  1. Kitch JA, Oates J (1994) The perceptual features of vocal fatigue as self-reported by a group of actors and singers. J Voice 8: 207–214.
  2. Titze I.R. (1984) Vocal fatigue: some biomechanical considerations. In: Lawrence VL (ed). (Transcripts of the Twelfth Symposium: Care of the professional voice Part One:Scientific Papers. New York: The Voice. Foundation;) Pg No: 97–104.
  3. Welham NV, Maclagan MA (2003) Vocal fatigue: current knowledge and future directions. J Voice 17: 21–30. [crossref]
  4. Cooper DS, Titze IR (1985) Generation and dissipation of heat in vocal fold tissue. J Speech Hear Res 28: 207–215. [crossref]
  5. Solomon NP (2008) Vocal fatigue and its relation to vocal hyperfunction . Int J Speech Lang Pathol 10: 254–266. [crossref]
  6. Sataloff RT () Professional Voice: The Science and Art of Clinical Care, 3-Volume Set, Fourth Edition. https://books.google.es/books?id=ImVGDwAAQBAJ&pg=PA547&dq =the+performance+voice+benninger&hl=pl&sa=X&ved= 0ahUKEwiUhqbJlLPZAhWDWxQKHRS7ARIQ6AEIcTAJ#v = onepage&q&f = false. Accessed February 20, 2018.
  7. Phyland DJ, Oates J, Greenwood KM (1999) Self reported voice problems among three groups of professional singers. 13: 602–611.
  8. Rocha C1, Moraes M, Behlau M (2012) Pain in popular singers. J Soc Bras Fonoaudiol 24: 374–380. [crossref]
  9. Faghy MA, Brown PI (2014) Thoracic load carriage-induced respiratory muscle fatigue. Eur J Appl Physiol 114: 1085–1093. [crossref]
  10. Wu YZ, Crumley RL, Armstrong WB, Caiozzo VJ (2000) New perspectives about human laryngeal muscle: single-fiber analyses and interspecies comparisons. Arch Otolaryngol Head Neck Surg 126: 857.
  11. Titze IR (1994) Mechanical stress in phonation. J Voice. 8: 99–105.
  12. Pearl Solomon N, Stemmle DiMattia M (2000) Effects of a vocally fatiguing task and systemic hydration on phonation threshold pressure. J Voice 14: 341–362
  13. Sivasankar M, Fisher KV (2002) Oral breathing increases Pth and vocal effort by superficial drying of vocal fold mucosa. J Voice 16: 172–181. [crossref]
  14. van Wyk L, Cloete M, Hattingh D, van der Linde J, Geertsema S (2017) The Effect of Hydration on the Voice Quality of Future Professional Vocal Performers. J Voice 31: 111. [crossref]
  15. Verdolini K, Min Y, Titze IR, Lemke J, Brown K, et al. (2002) Biological mechanisms underlying voice changes due to dehydration. J Speech Lang Hear Res 45: 268–281. [crossref]
  16. Leydon C, Wroblewski M, Eichorn N, Sivasankar M (2010) A meta-analysis of outcomes of hydration intervention on phonation threshold pressure. J Voice 24: 637–643.
  17. Boucher VJ (2008) Acoustic Correlates of Fatigue in Laryngeal Muscles: Findings for a Criterion-Based Prevention of Acquired Voice Pathologies. J Speech Lang Hear Res 51: 1161.
  18. Boucher VJ, Ayad T (2010) Physiological Attributes of Vocal Fatigue and Their Acoustic Effects: A Synthesis of Findings for a Criterion-Based Prevention of Acquired Voice Disorders. J Voice 24: 324–336.
  19. Caraty MJ, Montacié C (2014) Vocal fatigue induced by prolonged oral reading: Analysis and detection. Comput Speech Lang 28: 453–466.
  20. Behlau M, Oliveira G (2009) Vocal hygiene for the voice professional. Curr Opin Otolaryngol Head Neck Surg 17: 149–154. [crossref]
  21. Boucher VJ, Ahmarani C, Ayad T (2006) Physiologic features of vocal fatigue: Electromyographic spectral-compression in laryngeal muscles. Laryngoscope. 116: 959–965.
  22. Angsuwarangsee T, Morrison M (2002) Extrinsic laryngeal muscular tension in patients with voice disorders. J Voice 16: 333–343. [crossref]
  23. Timmermans B1, Vanderwegen J, De Bodt MS (2005) Outcome of vocal hygiene in singers. Curr Opin Otolaryngol Head Neck Surg 13: 138–142. [crossref]
  24. Yiu EM, Chan RM (2003) Effect of hydration and vocal rest on the vocal fatigue in amateur karaoke singers. J Voice 17: 216–227. [crossref]
  25. Titze IR (1999) Toward occupational safety criteria for vocalization. Logop Phoniatr Vocology 24: 49–54.
  26. Verdolini K, Ramig LO (2001) Review: occupational risks for voice problems. Logoped Phoniatr Vocol 26: 37–46. [crossref]
  27. Phyland DJ, Pallant JF, Benninger MS, et al. (2013) Development and preliminary validation of the EASE: A tool to measure perceived singing voice function. J Voice 27: 454–462.
  28. Smith B (2013) Choral Pedagogy. In: Brenda J, Sataloff RT (eds). (3rdedn). Plural Publishing, North America: 344.