low back pain is a common cause of disability with major socioeconomic consequences , affecting millions of individuals each year . although the pathophysiology of this condition is not yet fully known , intervertebral disc ( ivd ) degeneration has been shown to be a major cause of low back pain and lumbar disc herniation2 ) . since the exact pathogenesis and pathophysiology of disc degeneration is still unknown , most treatments for disc disease have focused on treating the numerous pathologic conditions arising from disc degeneration itself rather than on directly preventing , slowing , or reversing disc degeneration . recently , a number of novel therapeutic approaches have been proposed for direct management of the underlying causes of ivd degeneration . these approaches include synthetic peptide or growth factor injections1,5,10 ) , gene therapy13,14 ) , cell therapy3,11 ) , and tissue engineering4 ) . further advances in the treatment of disc degeneration will require improved understanding of the pathogenesis and pathophysiology of disc degeneration , together with testing the safety and efficacy of new approaches using suitable animal models that mimic disc degeneration in humans . previous studies have described animal models in which the ivd is stabbed to induce disc degeneration . the classic stab model by lipson and muir7 ) uses a full - thickness , ventral stab of the annulus fibrosus ( af ) to induce disc degeneration relatively quickly in rabbit lumbar discs . this model appears to produce biochemical and histological changes that are similar to those seen in disc degeneration in humans . however , the full - thickness stab model creates a macroscopic injury that induces rapid and severe disc degeneration ; thus , it does not mimic the slow , progressive disc degeneration seen in humans . degeneration after superficial annular injury is much less rapid than with the classic stab model . however , the degeneration is also less reproducible and may therefore be less useful for the assessment of treatment safety and efficacy . although both stab models result in disc degeneration , neither creates a relatively slow , progressive , and reproducible disc degeneration that allows researchers to study the safety and efficacy of new therapeutic interventions . the newly developed annulus puncture models9,17 ) , use hypodermic needles with controlled depth to produce a slower , more progressive , and more reproducible decrease in disc height than the classic stab procedure . however , this model requires an open surgical approach to the lumbar spine that is costly and bothersome to perform . the aim of the present study was to establish a simple , reproducible model of disc degeneration in rabbits through percutaneous puncturing of the annulus with needles in a minimally invasive manner . another aim was to confirm that this approach results in the slow , progressive disc degeneration that occurs in humans . to achieve these objectives , the author performed radiologic , magnetic resonance imaging ( mri ) , histological , and biological studies to examine the progress of disc degeneration . with approval from the institutional animal care and use committee , fifteen new zealand white rabbits ( 4 - 5 months old and weighing approximately 3 - 3.5 kg each ) underwent annular puncture of the l2-l3 , l3-l4 , and l4-l5 discs . five rabbits were allocated to 3 different groups according to the time at which they were sacrificed : 4 , 8 , or 20 weeks postoperatively . for the rabbits in the 20-week group , serial x - rays were performed at 0 , 2 , 4 , 8 , and 20 weeks to obtain longitudinal data on disc height changes over time . at each sacrifice time , the five rabbits were killed , and the whole lumbar spinal column was extracted ex vivo for mri analysis . after the scans , each rabbit was anesthetized with intramuscular injection of xylazine ( 5 mg / kg ) and ketamine ( 35 mg / kg ) , and the fur was shaved from the mid back and right flank . after anesthesia , a lateral plain x - ray was obtained to establish the pre - injection baseline height of the ivds . the rabbit was then placed in the lateral oblique prone position , and the injection field was sterilized with an alcohol sponge . initially , the l5-l6 disc was identified through manual palpation of the interspinous space from the mid back and pelvic rim with fluoroscopy ( vpx-200 ; toshiba corporation , tokyo , japan ) . after confirmation of the exact level , an 18-gauge angiography needle was inserted 3 - 4 cm ventral from the midline into the disc space under fluoroscopic control ( fig . after brief confirmation of the needle position in the center of the disc space by anteroposterior and lateral fluoroscopy , the needle was held in the disc space for 10 seconds . before removal of the punctured needle a constant amount of nucleus was leaked within the needle and identified by pushing the needle with an empty syringe after removal of the stylet ( fig . the l1-l2 and l5-l6 levels were designated as the non - punctured , internal controls . for each level , all procedures for identification and puncture were performed within 20 seconds . special care was taken to minimize contact with the periosteal tissues of the vertebrae because this could cause hypertrophy of the soft tissues and bony structures around the discs . the rabbits were placed in their cages after observation for recovery . at the planned time , they were sacrificed with intramuscular injection of ketamine ( 25.0 mg / kg ) followed by intravenous injection of sodium pentobarbital ( 1.2 g / kg ) . after rabbits were anesthetized , lateral plain radiographs of the lumbar spine were taken with a radiograph machine ( collimator - to film distance , 50 cm ; exposure , 5 mas ; penetration power , 44 kvp ) . during the radiographs , special care was taken to minimize axial rotation of the disc space by holding rabbits in the lateral decubitus position while ensuring the x - ray beam was straight . the lateral film was confirmed to be straight by checking for overlap of both transverse processes of the spine ( fig . the x - ray was repeated if a straight film was not obtained . to decrease the error from beam divergence in addition , each rabbit was treated with a consistent amount of anesthesia in order to provide a similar degree of muscle relaxation to minimize differences in disc height . radiographs were sent to a pacs server and analyzed with the piviewstar program ( infinitt healthcare , seoul , korea ) . the ivd height was expressed as the disc height index ( dhi ) and was analyzed quantitatively by two independent observers using the methods of l et al.8 ) with modification9 ) for evidence of changes in the vertebral body adjacent to the disc spaces of the stabbed and intact control discs . to avoid possible magnification error , changes in the dhi of injected discs were normalized to the preinjected ivd height and expressed as percent dhi ( % dhi ) . mri scans were performed on ex vivo spinal columns harvested from the rabbits . a 7-t mri machine ( biospec 70/20 usr ; bruker - biospin , billerica , ma , usa ) with a quadrature extremity coil was used . t2-weighted sections in the sagittal plane were obtained with the following settings : fast spin echo sequence with time to repetition of 3000 ms and time to echo of 70 ms ; 512256 matrix ; 10036 mm field of view ; and signal averaging of 2 . the mris were analyzed quantitatively using the modified thompson grade , in which changes in the degree and area of signal intensity were graded from 1 to 4 ( 1 , normal ; 2 , minimal decrease in signal intensity but obvious narrowing of the high signal area ; 3 , moderate decrease in signal intensity ; and 4 , severe decrease in signal intensity ) . at the three time points , the mri data for the punctured discs were pooled to assess the degree of degeneration over time . the l2-l3 and l3-l4 levels were extracted for real time quantitative polymerase chain reaction ( pcr ) analysis . the l1-l2 disc was also extracted as a non - punctured control . from each disc , the tissues were immediately placed into liquid nitrogen and frozen at -80 in preparation for pcr analysis . genes were selected from the representative forms related to the ecm component ( aggrecan ) and catabolic enzymes ( matrix metalloprotease-13 , mmp-13 ) . the frozen np samples were homogenized with a homogenizer ( mini - beadbeater ; bio spec , bartlesville , ok , usa ) in 1 ml trizol reagent ( invitrogen , carlsbad , ca , usa ) . the rt reaction was amplified in triplicate with real - time pcr ( abi prism 7900 ) in a final volume of 10 l using sybr green master mix reagent at a final concentration of 1 ( applied biosystems , foster city , ca , usa ) . glyceraldehyde phosphate dehydrogenase ( gapdh ) , a constitutively expressed housekeeping gene , was also amplified under the same conditions and used to normalize reactions . to ensure the specificity of the reaction , the size of the pcr product for each gene results were obtained with sequence detection software ( abi prism 7900 ) and evaluated with microsoft excel . the differences between the mean ct values of the gene of interest and gapdh were denoted ( ct ) , and the difference between ct of the degenerated sample and ct of the normal control sample was labeled ct . log2 ( ct ) gave the relative quantitation value of gene expression between the degenerated disc and the control disc in each animal . the primers for the rabbit - specific genes were designed in accordance with published sequences available in genbank . the sequences were as follows : aggrecan ( 5 ' gctacggagacaag gatgagttc 3 ' and 5 ' cgtaaaagacctcaccctccat 3 ' ) , mmp-13 ( 5 ' tgcccctcctcaacagtaac 3 ' and 5 ' gagcccgctgcattcttctt 3 ' ) , gapdh ( 5 ' aaggccatcaccatcttcca 3 ' and 5 ' ggatgcgttgctgacaatct 3 ' ) . the punctured disc l4 - 5 and the control disc l5 - 6 were harvested for histologic analyses . an electric saw was used to cut each disc together with the adjacent vertebral body . tissues were fixed with 10% neutral buffered formalin for 48 hours , decalcified in decalcification solution ( national diagnostics , atlanta , ga , usa ) for 3 days , and processed for paraffin sectioning . blocks embedded in paraffin were cut into mid - sagittal sections ( 4 m in thickness ) with a microtome . sections were stained with hematoxylin and eosin ( h&e ) and safranin o to determine the shape and number of cellular constituents as well as the amount of ecm , and then , analyzed under a light microscope ( nikon eclipse e800 ; nikon , melville , ny , usa ) at magnifications ranging from 40 to 200. for immunohistochemical staining , the sections were incubated with anti - collagen type ii mouse monoclonal antibody ( 1 : 100 dilution ; # ii-4c11 ; calbiochem , frankfurt , germany ) and then with hrp - labeled polymer - conjugated secondary antibody against mouse immunoglobulin g. the sections were stained with ready - to - use 3,3'-diaminobenzidine substrate - chromogen solution ( dako , produktionsvej , denmark ) and then counterstained with mayer 's hematoxylin . histologic findings included the number and nature of cells , matrix amount , and integrity of the border between the np and af . the grades were assigned 1 to 3 points for each of 4 categories , according to the severity of degeneration . the non - punctured disc at the l5-l6 level was also observed as a control in each rabbit to evaluate normal , age - related changes . longitudinal x - ray data were analyzed with repeated measures analysis of variance ( anova ) with multiple comparisons , adjusted by bonferroni method . mann - whitney u tests were used to compare degeneration for the mri , gene expression , and histology data . data were analyzed with pasw statistics 18.0 ( ibm , armonk , ny , usa ) . each rabbit was anesthetized with intramuscular injection of xylazine ( 5 mg / kg ) and ketamine ( 35 mg / kg ) , and the fur was shaved from the mid back and right flank . after anesthesia , a lateral plain x - ray was obtained to establish the pre - injection baseline height of the ivds . the rabbit was then placed in the lateral oblique prone position , and the injection field was sterilized with an alcohol sponge . initially , the l5-l6 disc was identified through manual palpation of the interspinous space from the mid back and pelvic rim with fluoroscopy ( vpx-200 ; toshiba corporation , tokyo , japan ) . after confirmation of the exact level , an 18-gauge angiography needle was inserted 3 - 4 cm ventral from the midline into the disc space under fluoroscopic control ( fig . after brief confirmation of the needle position in the center of the disc space by anteroposterior and lateral fluoroscopy , the needle was held in the disc space for 10 seconds . before removal of the punctured needle , the needle was rotated 180 degrees . a constant amount of nucleus was leaked within the needle and identified by pushing the needle with an empty syringe after removal of the stylet ( fig . the l1-l2 and l5-l6 levels were designated as the non - punctured , internal controls . for each level , all procedures for identification and puncture were performed within 20 seconds . special care was taken to minimize contact with the periosteal tissues of the vertebrae because this could cause hypertrophy of the soft tissues and bony structures around the discs . the rabbits were placed in their cages after observation for recovery . at the planned time , they were sacrificed with intramuscular injection of ketamine ( 25.0 mg / kg ) followed by intravenous injection of sodium pentobarbital ( 1.2 g / kg ) . after rabbits were anesthetized , lateral plain radiographs of the lumbar spine were taken with a radiograph machine ( collimator - to film distance , 50 cm ; exposure , 5 mas ; penetration power , 44 kvp ) . during the radiographs , special care was taken to minimize axial rotation of the disc space by holding rabbits in the lateral decubitus position while ensuring the x - ray beam was straight . the lateral film was confirmed to be straight by checking for overlap of both transverse processes of the spine ( fig . the x - ray was repeated if a straight film was not obtained . to decrease the error from beam divergence , in addition , each rabbit was treated with a consistent amount of anesthesia in order to provide a similar degree of muscle relaxation to minimize differences in disc height . radiographs were sent to a pacs server and analyzed with the piviewstar program ( infinitt healthcare , seoul , korea ) . the ivd height was expressed as the disc height index ( dhi ) and was analyzed quantitatively by two independent observers using the methods of l et al.8 ) with modification9 ) for evidence of changes in the vertebral body adjacent to the disc spaces of the stabbed and intact control discs . to avoid possible magnification error , changes in the dhi of injected discs were normalized to the preinjected ivd height and expressed as percent dhi ( % dhi ) . mri scans were performed on ex vivo spinal columns harvested from the rabbits . a 7-t mri machine ( biospec 70/20 usr ; bruker - biospin , billerica , ma , usa ) with a quadrature extremity coil was used . t2-weighted sections in the sagittal plane were obtained with the following settings : fast spin echo sequence with time to repetition of 3000 ms and time to echo of 70 ms ; 512256 matrix ; 10036 mm field of view ; and signal averaging of 2 . the mris were analyzed quantitatively using the modified thompson grade , in which changes in the degree and area of signal intensity were graded from 1 to 4 ( 1 , normal ; 2 , minimal decrease in signal intensity but obvious narrowing of the high signal area ; 3 , moderate decrease in signal intensity ; and 4 , severe decrease in signal intensity ) . at the three time points , the mri data for the punctured discs were pooled to assess the degree of degeneration over time . the l2-l3 and l3-l4 levels were extracted for real time quantitative polymerase chain reaction ( pcr ) analysis . the l1-l2 disc was also extracted as a non - punctured control . from each disc , the tissues were immediately placed into liquid nitrogen and frozen at -80 in preparation for pcr analysis . genes were selected from the representative forms related to the ecm component ( aggrecan ) and catabolic enzymes ( matrix metalloprotease-13 , mmp-13 ) . the frozen np samples were homogenized with a homogenizer ( mini - beadbeater ; bio spec , bartlesville , ok , usa ) in 1 ml trizol reagent ( invitrogen , carlsbad , ca , usa ) . the rt reaction was amplified in triplicate with real - time pcr ( abi prism 7900 ) in a final volume of 10 l using sybr green master mix reagent at a final concentration of 1 ( applied biosystems , foster city , ca , usa ) . glyceraldehyde phosphate dehydrogenase ( gapdh ) , a constitutively expressed housekeeping gene , was also amplified under the same conditions and used to normalize reactions . to ensure the specificity of the reaction , the size of the pcr product for each gene was verified by 2.0% agarose gel electrophoresis . results were obtained with sequence detection software ( abi prism 7900 ) and evaluated with microsoft excel . the differences between the mean ct values of the gene of interest and gapdh were denoted ( ct ) , and the difference between ct of the degenerated sample and ct of the normal control sample was labeled ct . log2 ( ct ) gave the relative quantitation value of gene expression between the degenerated disc and the control disc in each animal . the primers for the rabbit - specific genes were designed in accordance with published sequences available in genbank . the sequences were as follows : aggrecan ( 5 ' gctacggagacaag gatgagttc 3 ' and 5 ' cgtaaaagacctcaccctccat 3 ' ) , mmp-13 ( 5 ' tgcccctcctcaacagtaac 3 ' and 5 ' gagcccgctgcattcttctt 3 ' ) , gapdh ( 5 ' aaggccatcaccatcttcca 3 ' and 5 ' ggatgcgttgctgacaatct 3 ' ) . the punctured disc l4 - 5 and the control disc l5 - 6 were harvested for histologic analyses . an electric saw was used to cut each disc together with the adjacent vertebral body . tissues were fixed with 10% neutral buffered formalin for 48 hours , decalcified in decalcification solution ( national diagnostics , atlanta , ga , usa ) for 3 days , and processed for paraffin sectioning . blocks embedded in paraffin were cut into mid - sagittal sections ( 4 m in thickness ) with a microtome . sections were stained with hematoxylin and eosin ( h&e ) and safranin o to determine the shape and number of cellular constituents as well as the amount of ecm , and then , analyzed under a light microscope ( nikon eclipse e800 ; nikon , melville , ny , usa ) at magnifications ranging from 40 to 200. for immunohistochemical staining , the sections were incubated with anti - collagen type ii mouse monoclonal antibody ( 1 : 100 dilution ; # ii-4c11 ; calbiochem , frankfurt , germany ) and then with hrp - labeled polymer - conjugated secondary antibody against mouse immunoglobulin g. the sections were stained with ready - to - use 3,3'-diaminobenzidine substrate - chromogen solution ( dako , produktionsvej , denmark ) and then counterstained with mayer 's hematoxylin . histologic findings included the number and nature of cells , matrix amount , and integrity of the border between the np and af . the grades were assigned 1 to 3 points for each of 4 categories , according to the severity of degeneration . the non - punctured disc at the l5-l6 level was also observed as a control in each rabbit to evaluate normal , age - related changes . longitudinal x - ray data were analyzed with repeated measures analysis of variance ( anova ) with multiple comparisons , adjusted by bonferroni method . mann - whitney u tests were used to compare degeneration for the mri , gene expression , and histology data . data were analyzed with pasw statistics 18.0 ( ibm , armonk , ny , usa ) . to assess disc height change over time , % dhi was calculated in the 20-week group . when compared to baseline % dhi values obtained before the experiment , control discs showed no significant change in % dhi at any time point during the experiment . in contrast , the punctured discs showed significant narrowing of disc space at 2 weeks ( 87.2% , ) and at 4 weeks ( 80.5% ) ( p<0.001 , respectively ) ( fig . degeneration progressed gradually and was maintained over the study period ( p<0.001 , repeated measures anova ) . the mri data showed a loss of signal intensity and a decrease in area from mild to severe for each punctured disc ( l2 - 3 , l3 - 4 , and l4 - 5 ) on the t2-weighted midsagittal images ( fig . 4a ) . the two non - punctured control discs ( l1 - 2 and l5 - 6 ) maintained a consistent , normal grade over the same time period . at 4 weeks after puncture , the mri grade of the punctured discs had significantly increased compared to that of the control discs ( p<0.01 ) . the mri grade for each of the punctured discs also tended to increase over time , although this tendency was not statistically significant ( p=0.274 , the kruskal - wallis test ) ( fig . the gene expression results during the course of degeneration were calculated as relative percentages compared with the intact , non - punctured levels of mrna for each respective gene ( fig . 5 ) . representative ecm component aggrecan decreased to 45.5% at 4 weeks and to 30.7% at 20 weeks ( fig . 5a ) . the relative expression levels were not significantly different between time periods after puncture . for the catabolic enzyme mmp-13 , the levels of mrna decreased slightly at 4 weeks , increased at 8 weeks , and remained increased through week 20 ( fig . the gene expression levels were significantly increased over the time period ( p<0.05 ) . for both genes , expression levels in the punctured discs were significantly different from those of normal control discs at 20 weeks ( p<0.05 ) . representative h&e staining and safranin o staining of histologic sections from the rabbit discs showed various stages of degeneration ( fig . each healthy , non - punctured disc displayed an af with a normal pattern of fibrocartilage lamellas and a normal border between the af and np . the np was round and bloated in appearance and consisted of numerous large , notochordal cells ( fig . 6a , b depicts a disc with intermediate degeneration from a rabbit sacrificed at 4 weeks . in this disc , the annulus had a more wavy appearance , and the np was reduced in the extracellular matrix and the large vacuolated cells . the border between the annulus and nucleus was less distinct than in the normal disc . 6a , b represents a severely degenerated disc at the 20-week time point ; in this disc , most of the contents of the np had been lost , and the disc showed wavy fibrocartilage lamellas of the af with abundant associated fibroblasts . the immunohistochemical expression of type ii collagen in degenerated discs was decreased in the ecm and pericellular area of the np , compared with that of a normal control disc ( fig . the histologic score of the degenerated discs was significantly higher than that of control discs at each time point ( p<0.01 ) ( fig . the author found that the posterolateral puncture of rabbit lumbar discs with an 18-gauge angiography needle resulted in slow , progressive , and reproducible disc degeneration over 20 weeks , as confirmed by x - ray , mri , gene expression , and histologic findings . this needle puncture model appeared to induce disc degeneration more slowly than does the stab model7 ) , in which disc degeneration is achieved with a full - thickness , ventral annular lesion from a no . 11 blade . needle puncture with the selected angiography needle led to small , controlled leakage of the np and slow degeneration over time . in contrast , the classic stab model results in immediate herniation of a large bulk of the np , and studies have found almost no nuclear material in situ 1 to 2 days post stab7 ) and rapid narrowing of the disc space in 2 weeks9 ) . l et al.8 ) devised a simple radiographic technique for measurement of disc height : % dhi . for more accurate measurements , the author used masuda et al.'s9 ) modification of this technique . in addition , the author tried to properly position each rabbit , exactly align the x - ray beam , and maintain a consistent level of anesthesia in order to minimize any errors in disc height calculation . the film was repeatedly checked until both transverse processes of the spine were overlapped , suggesting that the true lateral view had been obtained ( fig . 2 ) . these methods permitted to obtain consistent results showing progressive degeneration of the discs over time . data from the rabbits in the 20-week group were used to examine serial changes in disc height over time . moreover , a repeated measures design was used to minimize any inter - subject variability in % dhi after disc puncture . a statistically significant decrease of % dhi , compared to that of control disc this result was comparable to other annular puncture model9 ) , which studied with various sized needles . the results of ex vivo mri assessment with the modified thompson grade were consistent with the radiographic findings . the results were also consistent with those of previous studies that used mri to verify disc degeneration in rabbits9,17 ) . the use of a simple numeric scale from 1 to 4 has shortcomings in terms of quantitative data analysis . this limitation can be minimized by using a quantitative , automated technique to define the region of interest with computer software17 ) . the biochemical changes associated with ivd degeneration include degradation of the ecm and an increase in catabolic enzymes . the major components of the ecm are proteoglycans and collagens , and aggrecan is the major proteoglycan in the disc . among the collagens in the disc , the catabolic involvement of mmp-13 on the ecm , especially for collagen type ii , has previously been demonstrated during degeneration of the ivd6,12,16 ) . the results of this study for changes in the aggrecan and mmp-13 are typical of those that occur during degeneration of the ivd . after early loss of vacuolated , notochordal cells within the np at 4 weeks after puncture , the af underwent infolding , the nucleus was replaced with fibrocartilage , and the ecm severely condensed . safranin o stains the proteoglycan content in ecm . in the degenerated discs , the intensity of safranin o was decreased . the expression of type ii collagen , a major ecm component in the np , decreased in degenerated discs . the histologic scale enabled the quantification of disc degeneration , and the results were consistent with those from the x - rays and mri . thus , this scale can be used in the assessment of biological treatments for disc degeneration . recently , open surgical needle puncture models have been developed , and these models appear to induce slower disc degeneration than does the classic stab model because they cause controlled leakage of the nuclear contents by annular puncture from needles of various sizes9,17 ) . the anterior surface of the disc is exposed from a posterolateral retroperitoneal approach , and the anterolateral af is punctured with a 16- or 18-gauge hypodermic needle at a controlled depth of penetration . however , the open puncture models require more technical support , time , and costs related to surgery and postsurgical care . in addition , excessive dissection or irritation of the ligaments in the perivertebral area can lead to osteophyte formation around the disc space . osteophytes make it difficult to inject therapeutic agents , such as proteins and cells . in general , when using a disc degeneration model to verify the safety and efficacy of new biological treatments , the delivery of new pharmaceutical agents should be simple and repeatable . unfortunately , in the open model , second delivery requires another open surgery , which is usually performed on the contralateral side1,10 ) . doing a repeated delivery is more complicated than an initial delivery due to the surgical stress to the experimental animals . however , in this percutaneous needle puncture model , the delivery of agents is simple and controllable through the use of fluoroscopy , and repeated injection can be performed easily at a predetermined time interval with minimal stress to the animals . to ensure that disc degeneration is reproducible and even , it is important to control the depth and location of the needle stab in the disc . previous open puncture models have controlled the depth at 5 mm with a handmade stopper9 ) or a locking forceps17 ) , but the location of the needle tip is determined by the researcher 's feel in the operative field . this model confirms the depth of the needle in real time using fluoroscopy , so obtaining a controlled depth is easy and reproducible . in addition , the exact location of the needle tip in the disc space is confirmed in the anteroposterior and lateral fluoroscopic views . preliminary experiments revealed the importance of extracting an equivalent amount of nucleus material during each puncture to create even disc degeneration ( data not shown ) . moreover , as the degree of needle rotation increases , a greater amount of disc degeneration develops . on the basis of these preliminary results , the holding time and degree of rotation of the punctured needle were maintained at 10 seconds and 180 degrees , respectively , in the present study . controlling the degree of disc degeneration is easily accomplished with this model and makes it possible to assess various hypotheses related to the grade of disc degeneration . one of the drawbacks of image - guided puncture is that it may cause repetition of the puncture . it is also technically demanding because the operator must avoid unnecessary contact with or damage to the ligamentous structures around the disc and the af itself in order to prevent uneven disc degeneration . excessive exposure to radiation is another concern , especially for researchers performing a large series of studies . thus , the use of protective devices , such as a lead gown , lead gloves , thyroid protector , and goggles , is mandatory when performing these procedures . disc degeneration from the annular puncture model has an issue that may not reflect an actual degenerative process as seen in clinics . however , controlled leakage of np using this fluoroscopic - guided needle puncture adjusted by holding time and degree of rotation may mimic various stages of disc herniation and its consequent degenerative changes . future studies using larger animals such as dogs , pigs , and non - human primates could support these findings . the author developed a rabbit model of degenerative disc disease that uses a novel percutaneous technique . the model was simple and minimally invasive and resulted in slow , progressive , and reproducible disc degeneration confirmed by radiology , biochemistry , and histology . this model can be used to study the pathogenesis of ivd degeneration and to evaluate the safety and efficacy of new biologic treatments for degenerative disc disease .
objectiveto develop a simple , reproducible model of disc degeneration in rabbits through percutaneous annular puncture and to confirm the degree of degeneration over time.methodsfifteen new zealand white rabbits ( 4 to 5 months old and weighing approximately 3 to 3.5 kg each ) underwent annular puncture of the l2-l3 , l3-l4 , and l4-l5 discs . 
 rabbits were sacrificed at 4 , 8 , or 20 weeks after puncture . for a longitudinal study to assess changes in disc height over time 
 , serial x - rays were performed at 0 , 2 , 4 , 8 , and 20 weeks for rabbits in the 20-week group . upon sacrifice 
 , the whole spinal column and discs were extracted and analyzed with magnetic resonance imaging ( mri ) , real time reverse transcriptase - polymerase chain reaction , and histological staining.resultsthe x - rays showed a slow , progressive decrease in disc height over time . 
 significant disc space narrowing compared to preoperative disc height was observed during the time period ( p<0.001 ) . 
 the mri grade , aggrecan , and matrix metalloprotease-13 mrna expression and hematoxylin and eosin / safranin o / anti - collagen ii staining were consistently indicative of degeneration , supporting the results of the x - ray data.conclusionpercutaneous annular puncture resulted in slow , reproducible disc degeneration that was confirmed by radiology , biochemistry , and histology . 
 this in vivo model can be used to study and evaluate the safety and efficacy of biologic treatments for degenerative disc disease .