ORIGINAL RESEARCH article
Improved right ventricular performance with increased tricuspid annular excursion in athlete’s heart
- 1Department of Medical Physics, University of Wisconsin-Madison, Madison, WI, USA
- 2Department of Medicine, University of Alabama at Birmingham, Birmingham, AL, USA
- 3Department of Biomedical Engineering, Wayne State University, Detroit, MI, USA
- 4Department of Electrical and Computer Engineering, Auburn University, Auburn, AL, USA
- 5Birmingham Veteran Affairs Medical Center, Birmingham, AL, USA
Background: Marathon runners (MTH) and patients with mitral regurgitation (MR) exhibit left ventricular (LV) overload, and LV geometric changes in these groups have been reported. In this study, right ventricular (RV) adaptation to chronic volume overload was evaluated in MTH and MR and normal controls together with interventricular septal remodeling and tricuspid annulus (TA) motion.
Methods: A total of 60 age-matched subjects (including 19 MTH, 17 isolated chronic compensated MR patients, and 24 normal subjects) underwent conventional cine and tagged cardiac magnetic resonance imaging. Myocardial strain and curvature were computed on the interventricular septum and RV free wall. A dual-propagation technique was applied to construct RV volume-time curves for a single cardiac cycle. Similarly, the TA was tracked throughout the cardiac cycle to create displacement over time curve.
Results: Septal curvature was significantly lower in MTH and MR compared to controls. No significant differences in RV free-wall strain or RV ejection fraction were noted among the three groups. However, longitudinal TA excursion was significantly higher in MTH compared to controls (p = 0.0061). The peak late diastolic TA velocity in MR was significantly faster than MTH (p = 0.0031) and controls (p = 0.020).
Conclusion: Increased TA kinetics allows for improved RV performance in MTH. Septal remodeling was observed in both MR and MTH, therefore a direct relationship of septal remodeling to TA kinetics in athlete’s heart could not be elucidated in this study.
The left ventricle undergoes remodeling in response to sustained changes in left ventricular (LV) pressure or volume load. In previous work, we used cardiac magnetic resonance (CMR) imaging to study LV changes in marathon runners (MTH), where the dilated left ventricle maintains an ellipsoid shape, and in mitral regurgitation (MR), where the dilated left ventricle becomes more spherical (1). Thus, the left ventricle responds differently to physiologic and pathologic chronic volume loading conditions. The stresses that lead to this remodeling are also transmitted to the right ventricle. The interventricular septum (IVS) allows for direct interaction of left and right ventricles and hence transmits systolic and diastolic forces between the ventricles. Previous studies (2–5) have shown that the deleterious effects of LV dilation in severe chronic organic MR on the structure and function of the right ventricle, leading to compression, flattening, and the consequential impairment of right ventricular (RV) systolic function. However, the influence of IVS remodeling in compensated MR on RV function has not been well studied. Furthermore, the relationship of differential LV remodeling due to physiologic versus pathologic LV volume overload to RV function has not been well described. In this study, three-dimensional (3D) geometric analysis and Lagrangian strain computation were used to define IVS remodeling and mechanics in compensated MR and MTH. We then evaluated the relationship of RV functional indices to IVS remodeling in these conditions. We hypothesized that maintenance of favorable LV and IVS geometry in physiologic LV overload noted in MTH would have favorable impact on RV function compared to pathologic LV volume overload in MR.
Materials and Methods
This study includes the same cohort of 60 subjects as previously described (1), which included 19 MTH (mean age 39 ± 10 years, 47% female), 17 patients with degenerative isolated MR (mean age 46 ± 5, 53% female), and 24 controls (mean age 45 ± 8 years, 50% female). Chronic isolated MR was defined as at least moderate severity with LV EF >60% based on echocardiographic/Doppler examination in the absence of symptoms or obstructive coronary artery diseases determined by exercise testing with nuclear perfusion. No MR patient had a history of hypertension or was taking any medication at the time of study. The study protocol was approved by the Institutional Review boards of the University of Alabama at Birmingham and Auburn University. All participants gave written informed consent.
Conventional cine MRI was acquired on a 1.5-T magnetic resonance scanner (GE Signa, Milwaukee, WI, USA) to obtain standard (2-, 3-, and 4-chamber long-axis and serial parallel short-axis) views using prospective electrocardiographically gated, breath-hold, steady state free-precession technique with following scan parameters: slice thickness = 8 mm, zero interslice gap, field of view = 40 cm × 40 cm, scan matrix = 256 × 128, flip angle = 45°, repetition/echo times = 3.8/1.6 ms. Twenty cardiac phases were reconstructed with 8–10 views per segment.
Tagged CMR was acquired with the same slice prescription as the cine acquisition. Grid tags were applied to short-axis views and stripe tagging to long-axis views using the spatial modulation of magnetization encoding method with the following scan parameters: prospective ECG gating, slice thickness = 8 mm, zero interslice gap, field of view = 40 cm × 40 cm, scan matrix = 256 × 128, flip angle = 10°, repetition/echo times = 8.0/4.2 ms, views per segment = 8–10, tag spacing = 7 mm, 20 reconstructed cardiac phases.
RV Volumetric and Tricuspid Annular Kinetics Analysis
Right ventricular endocardial contours at end-diastole (ED) (Figure 1) and end-systole (ES) were manually drawn as closed contours between the tricuspid annulus (TA) and RV apex from the short-axis views. These contours were then automatically propagated to all the other frames in the acquisition using a dual-contour propagation algorithm (6) as illustrated in Figure 2. Accurate RV segmentation in basal slices is difficult due to TA motion and partial volume effects as depicted in Figure 3. In conventional short-axis images, the RV inflow and outflow tracts are obscure. The borders among right ventricle, right atrium, and pulmonary artery were identified by viewing long-axis series and their slice projections simultaneously. To further account for this, a user-selected landmark point was specified at the TA on the RV lateral wall in both ED and ES frames of a four-chamber slice. Then, a non-rigid registration algorithm, similar to the one used in Ref. (6) for tracking mitral annular motion, was used to track this point through the cardiac cycle to form a displacement versus time curve. The displacement of this point perpendicular to the short-axis image plane was used to determine the fraction of basal short-axis slice that contributed to RV volume in that particular phase of the cardiac cycle. The validation of this propagation method on right ventricle is presented in the Data Sheet 1 in Supplementary Material.
Figure 2. The dual-propagation diagram on a mid-ventricular short-axis slice. The orange and blue arrows indicated the propagation using end diastolic (ED) and end systolic (ES) contours as the templates. These two sets of propagated contours were then combined via a weighted least-square fit to obtain the dual-propagated contour at each timeframe other than ED and ES.
Figure 3. Short-axis view of a basal RV slice. In basal slices of short-axis views, the right ventricle (red) coexists with the right atrium (RA) (magenta), aorta (blue), and pulmonary artery (PA) (green).
Once the volumes were computed at each time frame, a RV volume-time curve (VTC) was constructed and differentiated with respect to time. Early diastole and late diastole were defined as the first and second halves, respectively, of the diastolic interval. RV peak ejection rate (PER) was defined as the maximum negative time derivative during the systolic interval. RV early peak filling rate (ePFR) and atrial peak filling rate (aPFR) were defined as the maximum derivative during the early and late diastole. RV e/a ratio was computed as RVePFR over RVaPFR.
Interventricular Septal Regional Analysis
Three-dimensional (3D) IVS geometric parameters were measured from LV endocardial and epicardial contours manually traced on cine images acquired near end-diastole and end-systole as described previously (1, 7–15). The contours were traced to exclude papillary muscles. The contour data were then transformed to a coordinate system aligned along the long-axis of the left ventricle (Figure 4A) and converted to a prolate spheroidal coordinate system as described in Ref. (16). Cubic B-spline surfaces with 12 control points in the circumferential (θ) direction and 10 control points in the longitudinal direction (μ) were fit to the λ coordinates of the LV endocardial and epicardial contours (Figure 4B) for each time frame. The fit used the smoothing term described in Ref. (16), with α = 0, β = 0, and γ = 0.1.
Figure 4. Septal wall curvature computation. LV endocardial contours manually traced on magnetic resonance images were positioned in 3D (A). A surface was fit to these contours (B). Septal wall curvature was computed from the surface curvature at the septal wall (denoted by the green ball).
Three-dimensional (3D) endocardial surface curvatures were computed using standard formulas (17) at the septal wall segments defined in Ref. (18) (Figure 4B). Three-dimensional wall thickness was computed at the same segments by measuring the distance from a point on the LV epicardial surface to the closest point on the LV epicardial surface along a line perpendicular to the LV epicardial surface. The radius of curvature to wall thickness ratio (R/T) was computed by computing the reciprocal of the product of the endocardial circumferential curvature and wall thickness.
Three-dimensional LV strains were measured from tagged images at ES. Tag lines were tracked using a tag extraction technique (19) and edited by an expert user if necessary. 3D deformation and Lagrangian strain were computed by fitting a prolate spheroidal B-spline deformation model to the tag line data (20). The maximum principle strain was approximately aligned in the radial direction and is the maximum thickening strain. The maximum shortening strain was approximately aligned with circumferential direction and is the maximum contraction strain.
Interventricular septum geometric and strain parameters at the base and mid-ventricular levels were computed by averaging the LV anteroseptal and inferoseptal segments at each level. At the distal level, septal parameters were measured from the LV apical septal segment.
RV Curvature and Strain Measurements
Three-dimensional RV geometric parameters were measured using the techniques described above modified for the right ventricle. The RV lateral wall was divided into eight segments: three at base, three at mid, and two at apical level. 2D RV strains were calculated using harmonic phase (HARP) analysis (21). 2D RV mid-ventricular maximum shortening is the minimum principal strain averaged over the RV lateral wall segments at the mid-ventricular level.
One-way analysis of variance was used to compare groups for continuous and categorical variables as performed for left ventricle analysis among three subject groups. Homogeneity of variance was tested using Levene’s test. Model normality assumption was checked and if it was violated, appropriate data transformation was conducted. Tukey–Kramer procedure was performed to control the pairwise comparisons among the groups jointly in order to avoid erroneous type I error rate inflation.
Data are presented as mean ± SD. A p value <0.05 was considered statistically significant. For repeated measures, a p < 0.01 was considered statistically significant to account for correlations among parameters and locations of measurements (Table 2). All statistical analyses were performed using SAS version 9.2 (SAS Institute Inc., NC, USA).
As previously described (1), all participants had normal LVEF >55%. The RVEDV, RVESV, and stroke volume indices (normalized to body surface area) in MTH were significantly higher than MR or controls (Table 1). Mitral regurgitant volume and fraction of the MR patients was 29 ± 17 ml and 25 ± 14%. The RV mass index was significantly higher in MTH compared to MR or controls both with p < 0.0001(Table 1). This indicates that the increased demands on the right ventricle in MTH elicit greater remodeling effects than in MR. There are distinct differences in the overall RVVTCs of MTH and MR compared to controls as depicted in Figure 5. The absolute RVPER was significantly higher in MTH compared to controls (p = 0.044). However, once normalized to RVEDV, it was no longer different among the groups. Although there were differences in the RVePFR and RV e/a ratio between MTH and MR, they did not reach statistical significance.
Figure 5. RV volume-time curves for controls (solid line), marathon runners (gray line), and mitral regurgitation patients (dashed line). The solid, gray, dashed lines are the average RV volume-time curves of three groups. The error bars represent the SD at each measured time point.
The IVS radius-thickness ratio in ED was increased at distal level in MR versus controls (Table 2). The IVS circumferential curvature was significantly reduced at all three levels in MTH during both ED (p ≤ 0.0007 at all levels) and ES (p ≤ 0.0064 at all levels), whereas it was decreased at the basal (p = 0.0013) and mid (p = 0.0061) levels in ED and at the distal (p = 0.005) level in ES in MR (Table 2; Figure 6). Thus, the IVS is flatter in both MTH and MR compared to controls. The IVS wall thickness, maximum shortening, and maximum principle strain in the IVS were similar in the three groups at all levels. RV free-wall maximum shortening strains were also similar among three groups (Table 1). Despite similar thickening and strains, the mid and distal septal ES torsion were significantly less in MTH compared to controls (p = 0.0065 at mid; p = 0.0098 at distal). However, the local torsion shear angle in MTH was similar to both MR and controls.
Figure 6. End diastolic septal mid-ventricular curvature and RV tricuspid annulus displacement and strain map for a control, marathon runner, and mitral regurgitation patient. The marathon runner and the patient with MR had significantly increased septal mid curvature compared to the control. The tricuspid annulus displaced more in the marathon runner than it in the control. The RV lateral wall maximum shortening strain was similar in three groups.
Distinct differences were noted in the TA kinetics of MTH compared to controls and MR (Figure 7; Table 3). The peak TA displacement was significantly greater in MTH compared to controls (p = 0.0061). Consistent with RVVTC, the TA kinetics demonstrated significantly shorter time to peak displacement in MTH compared to controls. No difference in the TA kinetics of MR compared to controls or MTH was noted except that the peak late diastolic TA velocity was much faster than MTH (p = 0.0020) and control (p = 0.031) groups (Table 3). This along with somewhat increased atrial filling (Table 1) may indicate impairment of RV diastolic function in MR.
Figure 7. Tricuspid annulus displacement over time curve for controls (solid line), marathon runners (gray line), and mitral regurgitation patients (dashed line). The error bars represent the SD at each measured time point.
Cardiac magnetic resonance is known to yield accurate non-invasive RV volumetric measurements (22–24). In this paper, we examined structural and functional changes in the IVS and right ventricle in response to physiologic and compensated pathologic LV dilation states using standard CMR. Although the proposed modified RV short-axis series (25) can make the RV volume measurements easier and less error-prone, our technique provides a viable means to perform RV volumetric analysis in standard short-axis view and can be applied to clinical practice. RVVTC was derived with correction for TA displacement to provide more consistent systolic and diastolic parameters. Further, accurate characterization of septal and RV geometry and mechanics along with TA kinetics was performed using cine and tagged MRI.
Interventricular septum plays a critical role in the force transmission and maintenance of mechanical performance of biventricular function (4, 26–28). Here, we noted interesting adaptive changes in MTH and MR. The septal circumferential curvature was decreased in both groups compared to controls. This septal flattening is expected and reflects response to increased LV volume overload in MR and a complex interplay of pressure and volume overload on both right and left ventricles in MTH (29). There was no difference in septal principal and maximal strain measurements among the three groups. The preserved septal strains in MR likely reflect an unloading effect of compensated MR and a lack of significant pulmonary hypertension. In contrast, the normal LV geometry and lengthening maintained in MTH allows for preserved IVS strains.
Despite similar RV, EF and RV free-wall strains among the groups, the RVVTC suggests a more optimal systolic and diastolic profile in MTH compared to MR. Consistent with this, parallel changes in the TA kinetics in the three groups were also present. TA kinetics play, an important role in the overall RV systolic and diastolic function (30–32). The enhanced TA kinetics provides a mechanistic rationale for increased RV stroke volume and maintenance of RV systolic function. This finding in the MTH can all be attributed to the preload enhancing effect of increased RV volume.
In conclusion, our study described changes in the IVS mechanics and geometry due to physiologic versus pathologic left ventricle enlargement and its impact on RV function. We found that RV myocardial contractility and local curvature were similar in MTH and isolated compensated MR with preserved LVEF. There were subtle changes in septal mechanics in physiologic LV remodeling compared to controls. This combined with the findings of increased TA displacement reached in a shorter period of time in MTH suggests an overall favorable adaptive RV response in the MTH group.
WZ generated all the experimental data for validation and comparison, drew one set of contours for the inter-user variability study, and drafted the manuscript. WF developed the dual-contour propagation algorithm. CS helped in the design of the statistic testing. HG drew one set of contours served as the gold standard for validation, helped in the study design, and the revision of the manuscript. PR helped with manual contouring and study revision. SL identified manual contours to achieve consensus for validation and revised the manuscript. LD participated in study coordination. TD participated in the study design, coordination, and manuscript drafting.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by the National Heart, Lung, and Blood Institute Specialized Center for Clinically Oriented Research in Cardiac Dysfunction at the National Institutes of Health, Bethesda, MD, USA [grant number P50-HL077100 and R01-HL104018].
The Supplementary Material for this article can be found online at http://www.frontiersin.org/Journal/10.3389/fcvm.2015.00008/abstract
IVS, interventricular septum; MR, mitral regurgitation; MTH, marathon runners; RV, right ventricular; TA, tricuspid annulus; VTC, volume-time curve.
1. Schiros CG, Ahmed MI, Sanagala T, Zha W, McGiffin DC, Bamman MM, et al. Importance of three-dimensional geometric analysis in the assessment of the athlete’s heart. Am J Cardiol (2013) 111(7):1067–72. doi: 10.1016/j.amjcard.2012.12.027
2. Carlsson C, Haggstrom J, Eriksson A, Jarvinen A, Kvart C, Lord P. Size and shape of right heart chambers in mitral valve regurgitation in small-breed dogs. J Vet Intern Med (2009) 23(5):1007–13. doi:10.1111/j.1939-1676.2009.0359.x
4. Le Tourneau T, Deswarte G, Lamblin N, Foucher-Hossein C, Fayad G, Richardson M, et al. Right ventricular systolic function in organic mitral regurgitation: impact of biventricular impairment. Circulation (2013) 127(15):1597–608. doi:10.1161/CIRCULATIONAHA.112.000999
6. Feng W, Nagaraj H, Gupta H, Lloyd SG, Aban I, Perry GJ, et al. A dual propagation contours technique for semi-automated assessment of systolic and diastolic cardiac function by CMR. J Cardiovasc Magn Reson (2009) 11:30. doi:10.1186/1532-429X-11-30
7. Gupta A, Schiros CG, Gaddam KK, Aban I, Denney TS, Lloyd SG, et al. Effect of spironolactone on diastolic function in hypertensive left ventricular hypertrophy. J Hum Hypertens (2014). doi:10.1038/jhh.2014.83
8. Zheng J, Yancey DM, Ahmed MI, Wei CC, Powell PC, Shanmugam M. Increased sarcolipin expression and adrenergic drive in humans with preserved left ventricular ejection fraction and chronic isolated mitral regurgitation. Circ Heart Fail (2013) 7(1):194–202. doi:10.1161/CIRCHEARTFAILURE.113.000519
9. Schiros CG, Dell’Italia LJ, Gladden JD, Clark D III, Aban I, Gupta H, et al. Magnetic resonance imaging with 3-dimensional analysis of left ventricular remodeling in isolated mitral regurgitation: implications beyond dimensions. Circulation (2012) 125(19):2334–42. doi:10.1161/CIRCULATIONAHA.111.073239
10. Ahmed MI, Desai RV, Gaddam KK, Venkatesh BA, Agarwal S, Inusah S, et al. Relation of torsion and myocardial strains to LV ejection fraction in hypertension. JACC Cardiovasc Imaging (2012) 5(3):273–81. doi:10.1016/j.jcmg.2011.11.013
11. Ahmed MI, Aban I, Lloyd SG, Gupta H, Howard G, Inusah S, et al. A randomized controlled phase IIb trial of beta(1)-receptor blockade for chronic degenerative mitral regurgitation. J Am Coll Cardiol (2012) 60(9):833–8. doi:10.1016/j.jacc.2012.04.029
12. Gaddam K, Corros C, Pimenta E, Ahmed M, Denney T, Aban I, et al. Rapid reversal of left ventricular hypertrophy and intracardiac volume overload in patients with resistant hypertension and hyperaldosteronism: a prospective clinical study. Hypertension (2010) 55(5):1137–42. doi:10.1161/HYPERTENSIONAHA.109.141531
13. Ahmed MI, Gladden JD, Litovsky SH, Lloyd SG, Gupta H, Inusah S, et al. Increased oxidative stress and cardiomyocyte myofibrillar degeneration in patients with chronic isolated mitral regurgitation and ejection fraction >60%. J Am Coll Cardiol (2010) 55(7):671–9. doi:10.1016/j.jacc.2009.08.074
14. Ahmed MI, Sanagala T, Denney T, Inusah S, McGiffin D, Knowlan D, et al. Mitral valve prolapse with a late-systolic regurgitant murmur may be associated with significant hemodynamic consequences. Am J Med Sci (2009) 338(2):113–5. doi:10.1097/MAJ.0b013e31819d5ec6
15. Gladden JD, Ahmed MI, Litovsky SH, Schiros CG, Lloyd SG, Gupta H, et al. Oxidative stress and myocardial remodeling in chronic mitral regurgitation. Am J Med Sci (2011) 342(2):114–9. doi:10.1097/MAJ.0b013e318224ab93
18. Cerqueira MD, Weissman NJ, Dilsizian V, Jacobs AK, Kaul S, Laskey WK, et al. Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart: a statement for healthcare professionals from the cardiac imaging committee of the council on clinical cardiology of the American heart association. Circulation (2002) 105(4):539–42. doi:10.1161/hc0402.102975
19. Denney TS Jr, Gerber BL, Yan L. Unsupervised reconstruction of a three-dimensional left ventricular strain from parallel tagged cardiac images. Magn Reson Med (2003) 49(4):743–54. doi:10.1002/mrm.10434
22. Crean A, Maredia N, Ballard G, Menezes R, Wharton G, Forster J, et al. 3D Echo systematically underestimates right ventricular volumes compared to cardiovascular magnetic resonance in adult congenital heart disease patients with moderate or severe RV dilatation. J Cardiovasc Magn Reson (2011) 13(1):78. doi:10.1186/1532-429X-13-78
23. Grothues F, Moon JC, Bellenger NG, Smith GS, Klein HU, Pennell DJ. Interstudy reproducibility of right ventricular volumes, function, and mass with cardiovascular magnetic resonance. Am Heart J (2004) 147(2):218–23. doi:10.1016/j.ahj.2003.10.005
24. Sugeng L, Mor-Avi V, Weinert L, Niel J, Ebner C, Steringer-Mascherbauer R, et al. Multimodality comparison of quantitative volumetric analysis of the right ventricle. JACC Cardiovasc Imaging (2010) 3(1):10–8. doi:10.1016/j.jcmg.2009.09.017
25. Strugnell WE, Slaughter RE, Riley RA, Trotter AJ, Bartlett H. Modified RV short axis series – a new method for cardiac MRI measurement of right ventricular volumes. J Cardiovasc Magn Reson (2005) 7(5):769–74. doi:10.1080/10976640500295433
28. Banka VS, Agarwal J, Bodenheimer MM, Helfant RH. Interventricular septal motion: biventricular angiographic assessment of its relative contribution to left and right ventricular contraction. Circulation (1981) 64(5):992–6. doi:10.1161/01.CIR.64.5.992
30. Miller D, Farah MG, Liner A, Fox K, Schluchter M, Hoit BD. The relation between quantitative right ventricular ejection fraction and indices of tricuspid annular motion and myocardial performance. J Am Soc Echocardiogr (2004) 17(5):443–7. doi:10.1016/j.echo.2004.01.010
31. Haddad F, Hunt SA, Rosenthal DN, Murphy DJ. Right ventricular function in cardiovascular disease, part I: anatomy, physiology, aging, and functional assessment of the right ventricle. Circulation (2008) 117(11):1436–48. doi:10.1161/CIRCULATIONAHA.107.653584
32. Park JH, Kim JH, Lee JH, Choi SW, Jeong JO, Seong IW. Evaluation of right ventricular systolic function by the analysis of tricuspid annular motion in patients with acute pulmonary embolism. J Cardiovasc Ultrasound (2012) 20(4):181–8. doi:10.4250/jcu.2012.20.4.181
Keywords: cardiac magnetic resonance imaging, right ventricular function, tricuspid annulus displacement, interventricular septal remodeling, marathon runners, right ventricle strain, mitral regurgitation
Citation: Zha W, Schiros CG, Reddy G, Feng W, Denney TS Jr., Lloyd SG, Dell’Italia LJ and Gupta H (2015) Improved right ventricular performance with increased tricuspid annular excursion in athlete’s heart. Front. Cardiovasc. Med. 2:8. doi: 10.3389/fcvm.2015.00008
Received: 13 October 2014; Accepted: 13 February 2015;
Published online: 30 April 2015.
Edited by:Joseph B. Selvanayagam, Flinders University of South Australia, Australia
Reviewed by:Sylvia S. M. Chen, The Epworth Hospital, Australia
Christian Hamilton-Craig, University of Queensland, Australia
Copyright: © 2015 Zha, Schiros, Reddy, Feng, Denney, Lloyd, Dell’Italia and Gupta. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Himanshu Gupta, University of Alabama at Birmingham, 1808 7th Avenue South, BDB 101, CVMRI, Birmingham, AL 35294-0012, USA e-mail: firstname.lastname@example.org