Skip to main content
  • Original investigation
  • Open access
  • Published:

Circadian rhythm of autonomic activity in non diabetic offsprings of type 2 diabetic patients

Abstract

The aim of the present study was to evaluate, by heart rate variability (HRV) with 24-hours ECG Holter (HRV), the circadian autonomic activity in offspring of type 2 diabetic subjects and the relation with insulin-resistance. METHODS: 50 Caucasian offsprings of type 2 diabetic subjects were divided in two groups: insulin-resistant offsprings (IR) and non insulin-resistant offsprings (NIR). Autonomic nervous activity was studied by HRV. Time domain and spectral analysis (low frequency, LF, and high frequency, HF, provide markers of sympathetic and parasympathetic modulation when assessed in normalized units) were evaluated. RESULTS. Time domain showed a reduction of total SDNN in IR (p < 0.001) and NIR (p 0.047) versus controls. Spectral analysis showed a total and night LF higher in IR and NIR than in control group (all p < 0.001). CONCLUSION. In frequency domain, the analysis of sympathetic (LF) and parasympathetic (HF) component evidenced an association between the offspring of type 2 diabetic subjects and a sympathetic overactivity. A global reduction and alteration of circadian rhythm of autonomic activity are present in offspring of type 2 diabetic patients with and without insulin resistance. The data of our study suggested that an autonomic impairment is associated with the familiarity for type 2 diabetes independently to insulin resistance and that an impairment of autonomic system activity could precede the insulin resistance.

Background

Sympathetic and parasympathetic components of neurovegetative system regulate cardiac activity. Spectral analysis of heart rate variability (HRV) is a non invasive metodic used to assess cardiac autonomic activity. The autonomic activities can be assessed in HRV by the relative distribution (evaluated in normalized units) of low frequency (LF), as an index of sympathetic modulation and of high frequency (HF), as an index of parasympathetic modulation. The analysis of HRV has been used for evaluate the autonomic cardiac activity in numerous pathophisyologic conditions [1–9]: an impaired heart rate variation is a marker of autonomic neuropathy [10, 11] as in diabetic subjects [12–16].

Several studies have showed a sympathetic overactivity also in non diabetic insulin resistant group [17–19]. However, little is known about the association between the familiarity for type 2 diabetes mellitus and the autonomic activity. An increased LF/ HF ratio (low frequency/high frequency ratio) has identified in only insulin resistant offsprings of type 2 diabetic subjects [20]. De Angelis et al [21] has showed this condition at rest but not under stimulated conditions. Laitinen et al [22] has identified a sympathetic overactivity during acute hyperinsulinemia both in insulin resistant and non insulin resistant offsprings of type 2 diabetic subjects.

The aim of the present study was to evaluate, by heart rate variability with 24-hours ECG Holter registration, the autonomic activity and circadian autonomic rhythm in offspring of type 2 diabetic subjects and to evaluate the possible impact of sympathetic activity on insulin resistance. We tested the hypothesis that in non insulin resistant offsprings of type 2 diabetic subjects an alteration of circadian rhythm of autonomic activity are present and that the impairment of autonomic system activity could precede the insulin resistance.

Subjects and methods

70 consecutive caucasian offsprings of type 2 diabetic subjects, admitted to our department, were screened. In all subjects, after an overnight fast, oral glucose tolerance tests (OGTTs) was performed: samples blood for glucose and plasma insulin were collected before and 2 h after a glucose load consisting of 75 g glucose anhydrate in 300 ml of water ingested over the course of 5 min. Also, fasting plasma insulin was measured to evaluate the insulin-resistance by the homeostasis model assessment-index (HOMA-I). Among them, 50 caucasian subjects (age: 47.71 ± 9.96 years, 32 men and 18 women), with normal OGTTs, were admitted in this study and were divided in two groups: offsprings with insulin resistance and offsprings without insulin resistance:

Subjects with hypertension [23], diabetes mellitus, impaired fasting glycemia, impaired glucose tolerance [24], obesity, dyslipidemia, cardiac arrhythmias, microalbuminuria and with drug treatment or diseases that could potentially disturbs carbohydrate metabolism (glucocorticoids, furosemide, beta-blockers, etc.) and cardiac autonomic activity (beta-blockers, anti-arrhythmics, ACE-inhibitors) were excluded.

The control group consisted of 25 sex and age matched healthy non insulin resistant subjects with normal OGTTs and without familiarity for type 2 diabetes mellitus.

Height, weight and body circumferences were measured on all subjects. Body mass index (BMI, kg/m2) was calculated as weight divided by height squared. Waist-to-hip ratio (WHR) was defined as waist circumference divided by hip circumference.

Informed consent was obtained from all participants; all the investigations were performed in accordance with the participants of the Declaration of Helsinki.

Insulin resistance

The insulin-resistance was evaluated by the homeostasis model assessment index (HOMA-I) [25–28]. The HOMA-I was calculated by the formula: fasting plasma glucose (mmol/L) x fasting plasma insulin (μU/ml)/ 22,5, as described by Matthews and coworkers [29]. Insulin-resistance was defined as the third and fourth quartiles of HOMA-I.

The index subject were subdivide into two groups based on HOMA-I: 1) group of insulin-resistant offsprings (IR); 2) group of non insulin-resistant offsprings (NIR).

HRV assessment

Autonomic nervous activity was evaluated by heart rate variability (HRV) analysis during 24-hour ECG recording. All Holter recordings were performed using a three-channel recorder. Autonomic nervous activity was analysed following the recommendations of the Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology [30]. Spectral estimates of R-R interval were obtained from stationary regions free of ectopic beats and technical artefacts. The standard deviation of normal-to-normal RR intervals (SDNN) (ms) and the square root of the mean of the sum of the squares of differences between adjacent NN intervals (RMS-SD), correlated with parasympathetic system, were calculated and were divided in two periods: night (0 am – 6 am) and day (7 am – 9 pm). Fast Fourier Transform was used to obtain power spectral estimates of HRV. Total power in the frequency range (0 – 0.40 Hz) was divided into low frequency (LF: 0.04 – 0.15 Hz, modulated by sympathetic system) and high frequency (HF: 0.15 – 0.40 Hz, modulated by parasympathetic system). The power of LF and HF components was considered in normalized units (n.u.). Subjects were analysed for 24 hours, at 10 minutes interval. Artificial data and arrhythmic were excluded. The day was divided in four periods: night (0 am – 6 am), morning (7 am – 12 am), afternoon (1 pm – 6 pm), evening (7 pm – 11 pm). Data analyses were performed with software Del Mar Avionics Accuplus 363, Irvine California, USA.

Statistical analysis

All analysis were done with SPSS 12.0 (SPSS Inc., Chicago, IL, USA) for Windows XP. Data are presented as means + SD. For data with multiple time points, variables were analysed by the general linear model ANOVA and simple regression analyses were carried out by standard techniques 95% confidence intervals (CT) were calculated for regression coefficient. Means values were considered significant at p < 0.05.

Results

Clinical characteristics

The groups had not statistically significant difference in age, sex and anthropometrics parameters (i.e. BMI, waist and hip circumference, waist to hip ratio, table 1).

Table 1 Clinical characteristics for each study-group of offsprings and controls.

Heart rate variability

Tables 2 and 3 show the means of autonomic function measures for each group.

Table 2 Means (95%CI) RR intervals in the different time periods for each offspring group and controls.
Table 3 Autonomic function measures for each offspring group and controls.

Time domain

Time domain analysis of HRV showed a reduction statistically significant of total SDNN in IR (p < 0.001) and NIR (p 0.047) groups (IR: 86.10 ± 35.96 ms, NIR: 119.70 ± 28.22 ms) versus control group (131.90 ± 29.52 ms) (fig. 1). Total SDNN were reduced in IR group when compared with (NIR) NIR group (p 0.003). These results showed a total activity reduction of autonomic system in insulin resistant and non insulin resistant offsprings of type 2 diabetic patients. The autonomic activity reduction is major in IR group than NIR group.

Figure 1
figure 1

SDNN value in offsprings of type 2 diabetic patients and controls. SDNN: standard deviation of all sinus rhythm RR intervals, IR: insulin resistant offsprings of type 2 diabetic patients, NIR: non insulin resistant offsprings of type 2 diabetic patients. #p:0.003IR vs NIR; *p:0.047 controls vs NIR and † p:0.001 controls vs IR

RMS-SD in night time not increased in IR group. Therefore insulin-resistance was associated with alteration of circadian rhythm of parasympathetic component (at the night).

Frequency domain

Frequency domain analysis of HRV showed a total LF higher (all p < 0.001) in IR (70.88 ± 6.2 n.u.) and NIR (67.59 ± 4.46 n.u.) groups than in control group (55.46 ± 7.49 n.u.). No difference was noticed between IR and NIR groups (p: 0.157). HF were lower in IR (23.63 ± 6.67 n.u) and NIR (27.09 ± 3.74 n.u.) than in control group (37.65 ± 7.53) (all p < 0.001), while not statistically significant difference were between the groups of offsprings type 2 diabetic patients, p: 0.127. LF/HF total is higher in IR (3.34 ± 1.37) and NIR (3.72 ± 0.72) groups versus control group (2.44 ± 1.09), respectively p 0.016 and p 0.002.

Night (0 am – 6 am)

LF value are higher in IR (67.80 ± 9.70 n.u.) and NIR (59.60 ± 1.75 n.u.) groups than in control groups (35.21 ± 5.97) (all p < 0.001) and are higher in IR group when compared with NIR group, p 0.001). HF are lower in IR (25.72 ± 8.19 n.u.) and NIR group (36.53 ± 3.17) than in control groups (58.62 ± 6.59) (all p < 0.001 vs controls).

Morning (7 am – 12 am)

LF and HF had not statistically difference between IR (LF: 73.45 ± 9.61 n.u.; HF n.u.: 22.12 ± 7.17 n.u.), NIR (LF: 69.04 ± 7.41 n.u.; HF: 24.30 ± 8.52 n.u) and control groups (LF: 69.54 ± 8.36 n.u.; HF: 23.05 ± 7.13 n.u.).

Afternoon (1 pm – 6 pm)

LF, HF and LF/HF had not statistically difference between IR (LF: 70.59 ± 10.62 n.u.; HF: 23.21 ± 10.25 n.u.; LF/HF: 4.65 ± 1.71), NIR (LF: 70.06 ± 4.26 n.u.; HF: 21.42 ± 3.79 n.u.; LF/HF: 3.74 ± 0.46) and control groups (LF: 70.09 ± 6.62 n.u.; HF: 22.30 ± 5.09 n.u.; LF/HF: 3.67 ± 1.30).

Evening (7 pm – 11 pm)

LF were higher in IR (71.90 ± 10.70 n.u.) and NIR (71.90 ± 10.70 n.u.) groups than control group (56.75 ± 7.75 n.u.) (all p < 0.001). HF were lower in IR (23.52 ± 6.28 n.u.) and NIR (24.28 ± 9.08 n.u.) groups than control group (35.75 ± 9.14 n.u.), respectively p: 0.001 and p: 0.002. LF/HF is higher in IR (3.30 ± 1.00) and NIR (4.12 ± 2.21) groups than control group (2.00 ± 0.88) (respectively p: 0.04 and p: 0.002).

Discussion

The data of our study suggested that an autonomic impairment is associated with the familiarity for type 2 diabetes independently of insulin resistance. In frequency domain, the analysis of sympathetic (LF) and parasympathetic (HF) component and the symphatovagal balance (LF/HF) evidenced an association between the familiarity and a sympathetic overactivity, especially in nocturnal period, demonstrated by increase of LF (figure 2 and figure 3) and LF/HF ratio (figure 4). This autonomic impairment is major in insulin resistant offsprings than non insulin resistant offsprings of type 2 diabetic patients. Moreover, our study had demonstrated, in time domain analysis of HRV, a significant reduction of the total autonomic system activity in both groups, expressed by progressive decrease of SDNN value from NIR to IR groups. These results indicated that the familiarity of type 2 diabetes mellitus is related to a global reduction of autonomic nervous system and that the dysautonomia increases if offsprings are insulin resistant.

Figure 2
figure 2

LF value in offsprings of type 2 diabetic patients and control. IR: insulin resistant offsprings of type 2 diabetic patients, NIR: non insulin resistant offsprings of type 2 diabetic patients, LF: low frequency, hours of Holter registration: 0–6 am= night time, 7–12 am = morning. * p < 0.001 NIR and IR vs controls, #: p:0.001 IR vs NIR; & p < 0.001 controlsvs NIR and IR.

Figure 3
figure 3

Circadian variation of LF in IR, NIR and in control group. Gray circles: group of non insulin resistant offsprings of type 2 diabetic patients; black squares: group of insulin resistant offsprings of type 2 diabetic patients; white triangles: control group; IR: insulin resistant offsprings of type 2 diabetic patients, NIR: non insulin resistant offsprings of type 2 diabetic patients, LF: low frequency.

Figure 4
figure 4

LF/HF value in offsprings of type 2 diabetic patients and control. IR: insulin resistant offsprings of type 2 diabetic patients, NIR: non insulin resistant offsprings of type 2 diabetic patients, LF/HF: low frequency/high frequency, hours of Holter registration: 0–6 am= night time, 7–12 am = morning. *p.0.016 controls vs IR; †p:0.002 control vs NIR; # p:0.015 IR vs NIR; £ p:0.0001 controls vs NIR and IR; ç p:0.043 controls vs NIR.

In summary, a global reduction and alteration of circadian rhythm of autonomic activity are present in offspring of type 2 diabetic patients, without and with insulin resistance (figure 2, 3).

A long term observation can answer the question whether autonomic abnormality precede the occurence of insuline resistence and play a role in the complex pathogenesis of the insulin resistance and type 2 diabetes mellitus. Others studies occurs to explain a possible pathogenetic role of autonomic dysfunction in the development of insulin resistance and type 2 diabetes mellitus.

Limitations

A limitation in this study is the use of the HOMA index as a conventional indicator of insulin resistance. The best method for assessment of insulin resistance is the glucose clamp technique, however, the HOMA model has proved be a robust clinical and epidemiological tool in descriptions of the pathophysiology of diabetes, already quoted in > 500 publications, it has become one of standard tools in the armamentarium of the clinical physiologist [23].

Other study are necessary to determine the mechanism whereby insulin-resistance may be related to autonomic dysfunction.

Abbreviations

BMI :

body mass index

DM :

type 2 diabetes mellitus

HF :

high frequency

HOMA-I :

the homeostasis model assessment-index

HRV :

heart rate variability

LF :

low frequency

OGTTs :

oral glucose tolerance tests

RMS-SD :

the square root of the mean of the sum of the squares of differences between adjacent NN intervals

SDNN :

The standard deviation of normal-to-normal RR intervals

WHR :

waist-to-hip ratio.

References

  1. Pagani M, Lombardi F, Guzzetti S, Rimoldi O, Furlan R, Pizzinelli P, Sandrone G, Malfatto G, Dell'Orto S, Piccaluga E: Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog. Circ Res. 1986, 59 (2): 178-93.

    Article  CAS  PubMed  Google Scholar 

  2. Malliani A, Pagani M, Lombardi F, Cerutti S: Cardiovascular neural regulation explored in the frequency domain. Circulation. 1991, 84 (2): 482-92.

    Article  CAS  PubMed  Google Scholar 

  3. Montano N, Ruscone TG, Porta A, Lombardi F, Pagani M, Malliani A: Power spectrum analysis of heart rate variability to assess the changes in sympathovagal balance during graded orthostatic tilt. Circulation. 1994, 90 (4): 1826-31.

    Article  CAS  PubMed  Google Scholar 

  4. Schwartz P, La Rovere M, Vanoli E: Autonomic nervous system and sudden cardiac death: experimental basis and clinical observations for post-myocardial infarction risk stratification. Circulation. 1992, 85: 177-191.

    Google Scholar 

  5. Malliani A, Schwartz P, Zanchetti A: Neural mechanism in life threatening arrhythmias. Am Heart J. 1980, 200: 705-715. 10.1016/0002-8703(80)90238-0.

    Article  Google Scholar 

  6. Tubani L, Baratta L, Giorgino F, Delfino M, Fiore G, Golluscio V, Giacovazzo M: Heart rate variability in cluster headache. Ann Ital Med Int. 2003, 18 (1): 42-46.

    PubMed  Google Scholar 

  7. Baratta L, Delfino M, Fiorentini A, Martuscelli M, Tubani L: Alterations of sympathovagal bilance by heart rate variability: a rare case of adult Still's disease. Ann Ital Med Int. 2004, 19 (3): 193-197.

    PubMed  Google Scholar 

  8. Mastrocola C, Vanacore N, Giovani A, Locuratolo N, Vella C, Alessandri A, Baratta L, Tubani L, Meco G: Twenty-four-hour heart rate variability to assess autonomic function in Parkinson's disease. Acta Neurol Scand. 1999, 99 (4): 245-247.

    Article  CAS  PubMed  Google Scholar 

  9. Laganà B, Tubani L, Maffeo N, Vella C, Makk E, Baratta L, Bonomo L: Heart rate variability and cardiac autonomic function systemic lupus erythematosus. Lupus. 1996, 5 (1): 49-55.

    PubMed  Google Scholar 

  10. Ewing DJ, Neilson JM, Shapiro CM, Stewart JA, Reid W: Twenty four hour heart rate variability: effects of posture, sleep, and time of day in healthy controls and comparison with bedside tests of autonomic function in diabetic patients. Br Heart J. 1991, 65: 239-244.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Ziegler D: Diabetic cardiovascular autonomic neuropathy: prognosis, diagnosis and treatment. Diabetes Metab Rev. 1994, 10: 339-383.

    Article  CAS  PubMed  Google Scholar 

  12. Malpas SC, Maling TJB: Heart rate variability and cardiac autonomic function in diabetes. Diabetes. 1990, 39: 1177-1181.

    Article  CAS  PubMed  Google Scholar 

  13. Pagani M, Malfatto G, Pierini S, Casati R, Masu AM, Poli M, Guzzetti S, Lombardi F, Cerutti S, Malliani A: Spectral analysis of heart rate variability in the assessment of autonomic diabetic neuropathy. J Auton Nerv Syst. 1988, 23: 143-15. 10.1016/0165-1838(88)90078-1.

    Article  CAS  PubMed  Google Scholar 

  14. Bianchi A, Bontempi B, Cerutti S, Gianoglio P, Comi G, Natali Sora MG: Spectral analysis of heart rate variability signal and respiration in diabetic subjects. Med Biol Eng Comput. 1990, 28: 205-211.

    Article  CAS  PubMed  Google Scholar 

  15. Bellavere F, Balzani I, De Masi G, Carraro M, Carenza P, Cobelli C, Thomaseth K: Power spectral analysis of heart rate variation improves assessment of diabetic cardiac autonomic neuropathy. Diabetes. 1992, 41: 633-640.

    Article  CAS  PubMed  Google Scholar 

  16. Van den Akker TJ, Koelman ASM, Hogenhuis LAH, Rompelman G: Heart rate variability and blood pressure oscillations in diabetics with autonomic neuropathy. Automedica. 1983, 4: 201-208.

    Google Scholar 

  17. Brunner EJ, Hemingway H, Walker BR, Page M, Clarke P, Juneja M, Shipley MJ, Kumari M, Andrew R, Seckl JR, Papadopoulos A, Checkley S, Rumley A, Lowe GD, Stansfeld SA, Marmot MG: Adrenocortical, autonomic, and inflammatory causes of the metabolic syndrome: nested case-control study. Circulation. 2002, 106 (21): 2659-65. 10.1161/01.CIR.0000038364.26310.BD.

    Article  CAS  PubMed  Google Scholar 

  18. Pikkujamsa SM, Huikuri HV, Airaksinen KE, Rantala AO, Kauma H, Lilja M, Savolainen MJ, Kesaniemi YA: Heart rate variability and baroreflex sensitivity in hypertensive subjects with and without metabolic features of insulin resistance syndrome. Am J Hypertens. 1998, 11 (5): 523-31. 10.1016/S0895-7061(98)00035-1.

    Article  CAS  PubMed  Google Scholar 

  19. Flanagan DE, Vaile JC, Petley GW, Moore VM, Godsland IF, Cockington RA, Robinson JS, Phillips DI: The autonomic control of heart rate and insulin resistance in young adults. J Clin Endocrinol Metab. 1999, 84 (4): 1263-7. 10.1210/jc.84.4.1263.

    CAS  PubMed  Google Scholar 

  20. Frontoni S, Bracaglia D, Baroni A, Pellegrini F, Perna M, Cicconetti E, Ciampittiello G, Menzinger G, Gambardella S: Early autonomic dysfunction in glucose-tolerant but insulin-resistant offspring of type 2 diabetic patients. Hypertension. 2003, 41 (6): 1223-7. 10.1161/01.HYP.0000073062.29546.01.

    Article  CAS  PubMed  Google Scholar 

  21. De Angelis C, Perelli P, Trezza R, Casagrande M, Biselli R, Pannitteri G, Marino B, Farrace S: Modified autonomic bilance in offsprings of diabetics detected by spectral analysis of heart rate variabilità. Metabolism. 2001, 50 (11): 1270-74. 10.1053/meta.2001.27225.

    Article  CAS  PubMed  Google Scholar 

  22. Laitinen T, Vauhkonen IK, Niskanen LK, Hartikainen JE, Lansimies EA, Uusitupa MI, Laakso M: Power spectral analysis of heart rate variability during hyperinsulinemia in nondiabetic offspring of type 2 diabetic patients: evidence for possible early autonomic dysfunction in insulin-resistant subjects. Diabetes. 48 (6): 1295-9.

  23. European Society of Hypertension-European Society of Cardiology Guidelines Committee: Guidelines for the management of arterial hypertension. J Hypertens. 2003, 21 (6): 1011-53.

    Article  Google Scholar 

  24. American Diabetes Association: Diagnosis and classification of diabetes mellitus. Diabetes Care. 2004, 27 (Suppl 1): S5-S10.

    Google Scholar 

  25. Albareda M, Rodriguez-Espinosa J, Murugo M, de Leiva A, Corcoy R: Assessment of insulin sensitivity and beta-cell function from measurements in the fasting state and during an oral glucose tolerance test. Diabetologia. 2000, 43 (12): 1507-11. 10.1007/s001250051561.

    Article  CAS  PubMed  Google Scholar 

  26. Wallace TM, Levy JC, Matthews DR: Use and abuse of HOMA modeling. Diabetes Care. 2004, 27 (6): 1487-95.

    Article  PubMed  Google Scholar 

  27. Katsuki A, Sumida Y, Gabazza EC, Murashima S, Furuta M, Araki-Sasaki R, Hori Y, Yano Y, Adachi Y: Homeostasis model assessment is a reliable indicator of insulin resistance during follow-up of patients with type 2 diabetes. Diabetes Care. 2001, 24 (2): 362-5.

    Article  CAS  PubMed  Google Scholar 

  28. Bonora E, Targher G, Alberiche M, Bonadonna RC, Saggiani F, Zenere MB, Monauni T, Muggeo M: Homeostasis model assessment closely mirrors the glucose clamp technique in the assessment of insulin sensitivity: studies in subjects with various degrees of glucose tolerance and insulin sensitivity. Diabetes Care. 2000, 23 (1): 57-63.

    Article  CAS  PubMed  Google Scholar 

  29. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC: Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985, 28 (7): 412-9. 10.1007/BF00280883.

    Article  CAS  PubMed  Google Scholar 

  30. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology: Heart rate variability: standards of measurement, physiological interpretation and clinical use. Circulation. 1996, 93 (5): 1043-65.

    Article  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to A Fiorentini.

Additional information

Competing interests

The author(s) declare that they have no competing interests.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Fiorentini, A., Perciaccante, A., Paris, A. et al. Circadian rhythm of autonomic activity in non diabetic offsprings of type 2 diabetic patients. Cardiovasc Diabetol 4, 15 (2005). https://doi.org/10.1186/1475-2840-4-15

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2840-4-15

Keywords