Skip Navigation

Europace 2007 9(Supplement 6):vi71-vi76; doi:10.1093/europace/eum209
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Ruchat, P.
Right arrow Articles by Kappenberger, L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ruchat, P.
Right arrow Articles by Kappenberger, L.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The European Society of Cardiology 2007. All rights reserved. For permissions, please e-mail: journals.permissions@oxfordjournals.org

A biophysical model of atrial fibrillation ablation: what can a surgeon learn from a computer model?

Patrick Ruchat1,*, Nathalie Virag2, Lam Dang3, Jürg Schlaepfer4, Etienne Pruvot4 and Lukas Kappenberger5

1 Department of Cardiovascular Surgery, University Hospital, Lausanne, Switzerland; 2 Medtronic Europe, Tolochenaz, Switzerland; 3 Signal Processing Institute, École Polytechnique Fédérale de Lausanne, Switzerland; 4 Department of Cardiology, University Hospital, Lausanne, Switzerland; 5 Cardiomet, Medical Faculty, Lausanne, Switzerland

* Corresponding author: Service de Chirurgie Cardio-Vasculaire, CHUV Rue du Bugnon 46, CH-1011 Lausanne, Switzerland. Tel: +41 21 314 2280; fax: +41 21 314 2278.E-mail address: patrick.ruchat{at}chuv.ch


    Abstract
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Conclusion
 Funding
 Acknowledgements
 References
 
Aims: Surgical ablation procedures for treating atrial fibrillation have been shown to be highly successful. However, the ideal ablation pattern still remains to be determined. This article reports on a systematic study of the effectiveness of the performance of different ablation line patterns.

Methods and results: This study of ablation line patterns was performed in a biophysical model of human atria by combining basic lines: (i) in the right atrium: isthmus line, line between vena cavae and appendage line and (ii) in the left atrium: several versions of pulmonary vein isolation, connection of pulmonary veins, isthmus line, and appendage line. Success rates and the presence of residual atrial flutter were documented. Basic patterns yielded conversion rates of only 10–25 and 10–55% in the right and the left atria, respectively. The best result for pulmonary vein isolation was obtained when a single closed line encompassed all veins (55%). Combination of lines in the right/left atrium only led to a success rate of 65/80%. Higher rates, up to 90–100%, could be obtained if right and left lines were combined. The inclusion of a left isthmus line was found to be essential for avoiding uncommon left atrial flutter.

Conclusion: Some patterns studied achieved a high conversion rate, although using a smaller number of lines than those of the Maze III procedure. The biophysical atrial model is shown to be effective in the search for promising alternative ablation strategies.

Key Words: Atrial fibrillation, Computer modelling, Ablation, Surgery


    Introduction
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Conclusion
 Funding
 Acknowledgements
 References
 
Atrial fibrillation (AF) is a growing problem that needs to be treated to prevent its impact on mortality because of severe complications such as haemodynamic dysfunctions (loss of atrial mechanical activity and irregular ventricular response that can lead to low cardiac output and heart failure) and stroke.1Go This problem is approached differently by electrophysiologists and surgeons. From the electrophysiologist's point of view, propagation of electrical activity makes the heart beat about 4 x 107 per year. From the surgeon's point of view, the heart is a pump circulating about 4 x 106 litres of blood per year. Electrophysiologists are concerned about rhythm disturbances, whereas surgeons by a decrease in the transport of blood. However, whatever the focus of interest might be, there is no reliable biological model available for a pre-evaluation of the effectiveness of therapeutic interventions. Today, most treatments for AF are based on empirical considerations, or clinical studies and many underlying mechanisms remain unclear.

Surgical ablation of AF aims at creating lines of block to interrupt electrical conduction and, as a result, prevent the AF re-entrant process. The surgical Maze III procedure by Cox et al.2Go is the gold standard for patients undergoing open heart surgery for concomitant surgery such as valve replacement or coronary bypass. Although Cox et al.3Go,4Go reported excellent long-term results, the ideal AF ablation pattern still remains to be determined. This pattern should be able to terminate AF and prevent its re-occurrence with a limited number of ablation lines having minimal length, although allowing a correct mechanical activity/blood pumping during normal sinus rhythm. Furthermore, it would be desirable to apply this pattern with a method that is less invasive than surgery such as radiofrequency,5Go microwave, or cryoablation. No significant difference in the post-operative sinus conversion rate could be seen between the classical ‘cut and sew’ and the alternative sources of energy that were used to treat AF.6Go

The evaluation of different lesion patterns is usually performed in clinical studies or animal experiments. Computer modelling is a rapidly growing field, which is an additional tool in testing AF ablation patterns. The main advantages are the access to any value of interest and the repeatability of experiments in controlled conditions. Several computer models of atria have been developed over the last decades. Moe et al.7Go developed the first model of AF in the 1960s using cellular automata. Later, more sophisticated models were published, taking into account several aspects of the complex atrial anatomy.8Go–11Go Existing models differ in the tradeoff made between accuracy of the anatomical representation and the computational load. Modelling of ablation procedures requires simulations of long periods of AF on a computer modelling scale (several minutes instead of several seconds generally). The heavy computational demand of the advanced available atrial models precludes their use in this type of analysis. Therefore, a simpler model requiring far less computation effort was designed such that sufficient realism for the particular application was retained.12Go Good agreement was previously obtained between ablation success rates obtained from this computer model and those reported from clinical procedures.12Go,13Go

We have used this biophysical model of AF ablation to systematically study the performance of each ablation line generally performed by surgeons separately, as well as several of their possible combinations.


    Methods
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Conclusion
 Funding
 Acknowledgements
 References
 
Biophysical model of atrial fibrillation
The biophysical model of AF used in this study was specifically designed to take into account the main aspects of human atrial geometry, based on magnetic resonance images, while restricting computational load. As shown in Figure 1, the openings in the geometry of the connections to the valves and major vessels are represented. Next to the overall geometry, these openings are essential elements of the geometry, as atrial activation has to propagate around these orifices. The geometry was presented by a monolayer structure with homogeneous properties. The atrial tissue was discretized using a triangular mesh, the edges of which had an average length of 0.6 mm. At each node of the mesh, the electric activity was modelled by membrane kinetics models adjusted to match human atrial properties.14Go,15Go The resulting computational load (burden) is one hour of PC processor work to model one second of atrial activity.


Figure 1
View larger version (44K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 1 Biophysical model of human atria derived from magnetic resonance images. Tricuspid valve, TV; mitral valve, MV; inferior vena cava, IVC; superior vena cava, SVC; pulmonary veins, PVs; left/right atrium appendages, LAA/RAA; sino-atrial node, SAN. The atrial surface is meshed using triangular elements with an average length of 0.6 mm. At each node of the mesh, an atrial action potential, like the one shown, is computed and colour-coded (red for potentials greater than 20 mV and blue for resting potentials at –80 mV). Atrial fibrillation is initiated by rapid burst pacing at the SA node in a substrate of electrically remodelled atrial cells. The type of chronic AF simulated can be seen as being expressed by multiple re-entrant wavelets.

 
The biophysical model allows the simulation of different types of arrhythmias. Various conditions for initiation and perpetuation of AF in this model have been studied and were found to be in agreement with the limited available experimental data, either from detailed mapping on small segments of atrial tissue or from spot check recordings more widely spaced, as is obtained for electro-physiological (EP) studies.16Go

The study is based on a model of chronic AF in the form of sustained multiple re-entrant wavelets. The model tissue substrate of AF was set up on electrophysiological data recorded in patients suffering from permanent AF.15Go Atrial fibrillation was initiated by rapid burst pacing in the sino-atrial node region, leading to sustained AF, which continued unless terminated by an intervention such as the application of ablation lines.

Testing of ablation patterns
In the computer model, after AF had been initiated, ablation was simulated by instantaneously assigning an infinite resistivity between the cells located on the lines of the pattern studied. In this monolayer model, this completely blocked propagation across these lines, thus simulating a perfect (transmural) ablation.12Go

For testing each of the ablation line patterns, the following procedure was carried out (Figure 2). As a first step, as described earlier, sustained chronic AF was initiated by means of burst pacing. A database comprising 40 different simulated AF states was created by randomly selecting and storing subsequent states of AF. These 40 states differ in the activity present in the tissue, such as the number of wavelets. The ablation patterns were applied instantaneously and were considered successful if AF termination was observed within an arbitrary period of 30 s. A detailed analysis of the activity following the application of lines indeed showed that termination occurring after this period was generally because of spontaneous events. When applied to the AF database, the percentage of the 40 simulations in which AF terminated was taken as the success rate of the ablation procedure. In addition to the percentage of return to sinus rhythm, the conversion to atrial flutter was also documented.


Figure 2
View larger version (29K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 2 Simulation of ablation patterns. In the first step, AF is initiated via burst pacing. A database of 40 elements is constructed by randomly selecting different time instants in the simulated electrogram. Each ablation pattern is applied abruptly to the AF states from this database and the result after 30 s is evaluated and classified as follows: return to sinus rhythm, conversion to atrial flutter, or perpetuating AF.

 
The AF ablation patterns tested included those commonly performed by surgeons (Figure 3). Different versions of the individual lines present in a full Maze III procedure, both in the right and left atria, were tested individually as well as several combinations of these lines in increasing numbers. Basic ablation lines in the right atrium were an isthmus line (RAI), a line between vena cavae (SIVC), and an appendage line (RAA). Basic ablation lines in the left atrium were several versions of pulmonary vein isolation (one by one IPV1, two by two IPV2, and four together IPV4), a connection of pulmonary veins (CPV), an isthmus line (LAI), and an appendage line (LAA).


Figure 3
View larger version (39K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 3 Systematic testing of ablation patterns with a comparison to Maze III as a reference. Patterns in the right atrium: right isthmus line between tricuspid valve and inferior vena cava (RAI), line between superior and inferior vena cava (SIVC), and right atrial appendage line (RAA). Patterns in the left atrium: isolation of pulmonary veins one by one at 5 mm distance to ostia (IPV1), isolation of pulmonary veins two by two with a 5 mm distance between ablated areas (IPV2), isolation of four pulmonary veins (IPV4), IPV4 and left atrial isthmus line between inferior pulmonary veins and mitral valve (LAI), connection of three pulmonary veins (CPV) and LAI, left atrial appendage line (LAA).

 
Each ablation pattern required a simulation time of a total of up to 20 min (40 different AF states times 30 s, less if AF terminated sooner). Generally, this took several hours of simulated AF and several days of computation time on a set of several PCs.


    Results
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Conclusion
 Funding
 Acknowledgements
 References
 
The best result and the shortest time to AF termination were obtained in the computer model with the standard Maze III procedure (100%), which was used as a reference in this study. The percentages of successful conversion to sinus rhythm are summarized in Table 1. This table contains three parts ordered by increasing number of lines in each part: (i) right lines only, (ii) left lines only, and (iii) both right and left lines.


View this table:
[in this window]
[in a new window]

 
Table 1 Success rates in converting AF to sinus rhythm for each simulated ablation pattern

 
Each basic ablation pattern presented in Figure 3 was evaluated individually. For the right atrium, the basic patterns yielded conversion rates of 10–25%. For the left atrium, conversion rate was higher by 10–55%, but with a possibility of conversion into uncommon left atrial flutter. If we compare the different possibilities for isolation of pulmonary veins, the best result was obtained when a single closed ablation line encompassed all four veins (55%), or when they were connected together with the addition of the left isthmus line (60%).

The best result in the right atrium, with 65% conversion, was obtained by combining all basic lines, whereas a combination of all left lines lead to a success rate of 80%. Higher rates could be obtained only if right and left lines were combined. Some patterns achieved the maximum conversion rate using a smaller number of ablation lines than the Maze III procedure. This was the case for the patterns combining an optimal isolation of pulmonary veins, the left isthmus line, and the line between vena cavae in the right atrium. We could also observe that the left isthmus line is required in the left atrium in order to avoid the occurrence of uncommon left atrial flutter at some stage after the ablation.


    Discussion
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Conclusion
 Funding
 Acknowledgements
 References
 
Computer modelling has some limitations, which render the comparison with clinical data sometimes difficult. The timescales of computer modelling and clinical data are very different. The longest computer-simulated AF lasts only a few minutes, whereas AF in reality can last for days, months, or years. To overcome this problem, this model of chronic AF was developed by introducing electrical remodelling as observed in atrial cells after 1 or 2 weeks of pacing-induced AF with arrhythmia becoming self-sustained and atrial effective refractory period markedly shortened.17Go Therefore, it is only the transition phase, following the application of the ablation lines, that is modelled. Long-term effects of ablation cannot be studied clearly by means of this approach.

In this type of study, the geometry and electrophysiological properties are the same in all experiments, and the comparison of the effectiveness of the different ablation line patterns is devoid of confounding factors such as differences between patients and their clinical condition. The results shown here relate to chronic AF involving multiple re-entrant wavelets. Other types of AF, such as the ones involving rapid foci originating from the pulmonary veins as observed during paroxysmal AF,18Go which may create differences in success rate for the lines located in the left atrium, could also be studied by means of similar models with different properties.

Despite some of the limitations, the biophysical model is an interesting tool to assess surgical ablation patterns for patients suffering from permanent AF. It offers the possibility to test lines in a ‘reversible’ way in a human model and to observe the AF termination process in detail. A previous comparison between simulations performed with this model and patients who underwent radiofrequency ablation showed a positive correlation for conversion rates to sinus rhythm and residual atrial flutter.13Go Therefore, conversion rates can be studied realistically for different patterns in the case of chronic substrates of AF, which is the usually the case when surgeons perform ablation procedures during open heart surgery.

The results obtained for the simulated Maze III procedure are comparable with clinical data reported for this pattern, showing a long-term success rate ranging from 80 to 99%.4Go The success rate obtained for the simulations of this pattern is very high (100%) because all the lines created are perfect and non-conducting, which may not be the case in reality. Indeed, previous simulation studies have shown that performance can decrease with Maze III from 100 to 80% for a gap of 3 mm located in the line connecting the pulmonary veins and mitral annulus.12Go This decrease in performance is linked to the apparition of a residual atypical atrial flutter around the mitral valve.

The systematic testing of the ablation lines composing the Maze III patterns confirmed the need of having lines in both the right and left atria. This study is a generalization of the study presented by Ruchat et al.,19Go in which a minimal ablation pattern was defined by successfully adding ablation lines until a full Maze III procedure, suggesting that ablation in the right atrium could be simplified to a single line joining both vena cavae. The results shown in Table 1 indicate that combination of lines only in the right/left atrium may lead to a success rate in the range 10–65/10–85%. The maximum conversion rates obtained are comparable with the published results of clinical studies: long-term success rate of atrial compartmentalization of 67%20Go and long-term success rate of pulmonary vein region ablation of 82%.21Go

As could be expected, individual lines generally have little impact on AF (10–25%). A noteworthy exception formed the isolation of pulmonary veins with a success rate in the range 20–55%. Isolation of pulmonary veins can be performed using different approaches, mostly developed to treat paroxysmal AF.22Go The situation is different in the case of permanent AF, and our study shows that better conversion rates are obtained if all pulmonary veins are encompassed by one single closed ablation line rather than ablated individually. In our case, this is due to the fact that one area will be less likely to generate wavebreaks sustaining AF.

Our simulation studies indicated that all basic lines in the left atrium can lead to atypical atrial flutter, whereas this is not the case for right lines. These atrial flutters can be stopped by the introduction of a left atrial isthmus line. This line should be transmural to be effective.12Go,13Go In the case of isolation of pulmonary veins one by one or two by two, the rate of residual atrial flutter is also dependent on the size of the area between ablated zones, as shown in a previous study.23Go

Combination of lines in both the right and left atria increased conversion rate in the range 90–100% (assuming perfect transmural lines). This biophysical model gives some insight into the ablation process for an optimal minimal pattern. One or two lines are required in the right atrium, but the line choice does not seem to be critical. For the left atrium, the critical areas to be ablated are the pulmonary vein region and the left isthmus line.


    Conclusion
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Conclusion
 Funding
 Acknowledgements
 References
 
This model-based study aimed at identifying ablation patterns with high conversion rates, similar to the best ones reported in the literature, although causing less damage to the tissue. A potential pattern that came up so far is the isolation of the pulmonary veins combined with a left isthmus line to prevent atrial flutter and a line between vena cavae or a right atrial isthmus line.


    Funding
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Conclusion
 Funding
 Acknowledgements
 References
 
This study was made possible by grants from the Theo-Rossi-Di-Montelera Foundation, Medtronic Europe, the Swiss Governmental Commission of Innovative Technologies (CTI), and the Swiss National Science Foundation (SNSF).


    Acknowledgements
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Conclusion
 Funding
 Acknowledgements
 References
 
The authors wish to thank Josée Morisette, Arthur Stillman and Ryan Lahm, who kindly furnished the atrial geometry surface model.

Conflict of interest: none declared.


    References
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Conclusion
 Funding
 Acknowledgements
 References
 
[1] Braunwald E. Shattuck lecture—cardiovascular medicine at the turn of the millennium: triumphs, concerns and opportunities. N Engl J Med (1997) 337:1360–9.[Free Full Text]

[2] Cox JL, Jaquiss RD, Schuessler RB, Boineau JP. Modification of the maze procedure for atrial flutter and atrial fibrillation. II. Surgical technique of the maze III procedure. J Thorac Cardiovasc Surg (1995) 110:485–95.[Abstract/Free Full Text]

[3] Cox JL, Ad N, Palazzo T. Impact of the maze procedure on the stroke rate in patients with atrial fibrillation. J Thorac Cardiovasc Surg (1997) 63:1070–5.

[4] Cox JL. Atrial transport function after the maze procedure for atrial fibrillation: a 10 year clinical experience. Am Heart J (1998) 136:934–6.[CrossRef][Web of Science][Medline]

[5] Ruchat P, Schlaepfer J, Delabays A, Hurni M, Milne J, von Segesser LK. Left atrial radiofrequency compartmentalization for chronic atrial fibrillation during heart surgery. Thorac Cardiovasc Surg (2002) 50:155–9.[CrossRef][Web of Science][Medline]

[6] Khargi K, Hutten BA, Lemke B, Deneke T. Surgical treatment of atrial fibrillation: a systematic review. Eur J Cardiothorac Surg (2005) 27:258–65.[Abstract/Free Full Text]

[7] Moe G, Rheinbold W, Abildskov J. A computer model of atrial fibrillation. Am Heart J (1964) 67:200–20.[CrossRef][Web of Science][Medline]

[8] Harrild D, Henriquez C. A computer model of normal conduction in the human atria. Circ Res (2000) 87:e25–36.[Abstract/Free Full Text]

[9] Vigmond EJ, Ruckdeschel R, Trayanova N. Reentry in a morphologically realistic atrial model. J Cardiovasc Electrophysiol (2001) 12:1046–54.[CrossRef][Web of Science][Medline]

[10] Kneller J, Kalifa J, Zou R, Zaitsev AV, Warren M, Berenfeld O, et al. Mechanisms of atrial fibrillation termination by pure sodium channel blockade in an ionically-realistic mathematical model. Circ Res (2005) 96:35–47.[CrossRef]

[11] Seeman G, Hoper C, Sachse FB, Dossel O, Holden AV, Zhang H. Heterogeneous three-dimensional anatomical and electrophysiological model of human atria. Philos Trans R Soc A (2006) 364:1465–81.[CrossRef]

[12] Dang L, Virag N, Ihara Z, Jacquemet V, Vesin J-M, Schlaepfer J, et al. Evaluation of ablation patterns using a biophysical model of atrial fibrillation. Ann Biomed Eng (2005) 33:465–74.[CrossRef][Web of Science][Medline]

[13] Ruchat P, Dang L, Schlaepfer J, Virag N, von Segesser LK, Kappenberger L. Use of a biophysical model of atrial fibrillation in the interpretation of the outcome of surgical ablation procedures. Eur J Cardiothorac Surg (2007) 32:90–5.[Abstract/Free Full Text]

[14] Luo C-H, Rudy Y. A model of the ventricular cardiac action potential. Circ Res (1991) 68:1501–26.[Abstract/Free Full Text]

[15] Li D, Zhang L, Knellet J, Nattel S. Potential ionic mechanism for repolarization differences between canine right and left atrium. Circ Res (2001) 88:1168–75.[Abstract/Free Full Text]

[16] Jacquemet V, Virag N, Kappenberger L. Wavelength and vulnerability to atrial fibrillation: insights from a computer model of human atria. Europace (2005) 7:S83–92.[Abstract/Free Full Text]

[17] Thijssen VLJL, Ausma J, Borgers M. Structural remodelling during chronic atrial fibrillation: act of programmed cell survival. Cardiovasc Res (2001) 52:14–24.[Abstract/Free Full Text]

[18] Haissaguerre M, Jais P, Shah DC, Takashi A, Hocini M, Quiniou G, et al. Spontaneous initiation of atrial fibrillation by ectopic beats originating in the pulmonary veins. N Engl J Med (1998) 339:659–66.[Abstract/Free Full Text]

[19] Ruchat P, Dang L, Virag N, Schlaepfer J, von Segesser LK, Kappenberger L. A biophysical model of atrial fibrillation to define the appropriate ablation pattern in modified maze. Eur J Cardiothorac Surg (2007) 31:65–9.[Abstract/Free Full Text]

[20] Garg A, Finneran W, Mollerus M, Birgersdotter-Green O, Fujimura O, Time L, et al. Right atrial compartmentalization using radiofrequency catheter ablation for management of patients with refractory atrial fibrillation. J Thorac Cardiovasc Surg (1999) 10:763–71.

[21] Kottkamp H, Hindricks G, Hammel D, Austschbach R, Mergenthaler J, Borggrefe M, et al. Intraoperative radiofrequency ablation of chronic atrial fibrillation: a left curative approach by elimination of anatomic ‘anchor’ reentrant circuit. J Cardiovasc Electrophysiol (1999) 10:772–80.[Web of Science][Medline]

[22] Morady F. Catheter ablation of supraventricular arrhythmias: state of the art. J Cardiovasc Electrophysiol (2004) 43:2044–53.

[23] Rotter M, Dang L, Jacquemet V, Virag N, Kappenbeger L, Haissaguerre M. Impact of varying ablation patterns in a simulation model of persistent atrial fibrillation. Pacing Clin Electrophysiol (2007) 30:314–21.[CrossRef][Medline]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Ruchat, P.
Right arrow Articles by Kappenberger, L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ruchat, P.
Right arrow Articles by Kappenberger, L.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?