Qualitative and quantitative dermatoglyphics of chronic kidney disease of unknown origin (CKDu) in Sri Lanka

Abstract

Background

Dermatoglyphics has been used widely in fields of medicine as a non-invasive diagnostic tool and an early assessment of risk for certain medical conditions. It reflects disturbances in fetal development during early prenatal weeks 14–22 when fingerprints develop. Dermatoglyphic asymmetry has been used to measure developmental instability during a specific period of human fetal development. Thus, the present study was planned to investigate whether digital and palmar dermatoglyphics of chronic kidney disease of unknown origin (CKDu) patients in Sri Lanka are different from healthy people.

Methods

A case control study was carried out among CKDu patients (90 males, 90 females) from a CKDu endemic area and gender-matched two control groups; one group from a CKDu endemic region (90 males, 90 females) and another group from a CKDu non-endemic region (90 males, 90 females). Dermatoglyphics were obtained using photographic methods. Both qualitative and quantitative dermatoglyphic variables were defined and analyzed according to standard criteria. Both directional (DA) and fluctuating asymmetry (FA) were assessed.

Results

Several qualitative dermatoglyphic variables had significant association with CKDu. The triradii a¹ variable was less evident in palms of CKDu cases in both genders when compared to both control groups. The FA of pattern discordance (right vs left hands) between CKDu cases and control group were significant in several digits. The FA of the ridge count was found significant in several digits, and also significant for A-B ridge count and total ridge count.

Conclusion

Based on these results, it is proposed that the mechanisms responsible for the development of CKDu might be associated with those responsible for FA observed in CKDu patients. Accordingly, a diagnostic tool based on FA could be developed for predicting risk prior to the development of CKDu.

Background

Chronic kidney disease (CKD) is a global public health problem and progressively more common in developed as well as developing nations [1]. According to a Global Burden of Disease study in 2016, the age-standardized death rate for CKD was 18.2 per 100,000 people and disability-adjusted life-years (DALYs) was 473.9 [2]. The leading underlying risk factors for CKD are diabetes and hypertension [2]. However, a devastating form of CKD has been observed recently in several regions including Sri Lanka, Central America, India, and Egypt which was not attributed to conventional risk factors [3]. It has been referred to as CKD of unknown etiology (CKDu), and recent evidence suggested causative factors could be of environmental and/or occupational origin [4]. Due to the epidemiological pattern and histopathological similarities of the disease, a new entity called chronic interstitial nephropathy of agricultural communities (CINAC) was proposed to identify this condition [4, 5]. Despite extensive investigations, its exact etiology is still a mystery.

Dermatoglyphics is the study of epidermal ridge patterns (fingerprints) on the skin of the fingers, palms, toes, and soles that commence during embryological development between the sixth and seventh week of intrauterine life, and are fully formed by the 21st week [6]. It has been used widely in fields of anthropology, genetics, and medicine and as a valuable non-invasive diagnostic tool and early assessment of risk for certain medical conditions [7]. The relationship between different dermatoglyphic traits and various medical diseases have been widely evaluated, and the main hypothesis for support of this association is “if growth of the limbs is disturbed in very early fetal life, changes in the epidermal ridge configurations are likely” [8,9,10]. It should be added, however, that both environmental and genetic factors do influence the development of dermatoglyphics [6, 10,11,12,13]. Wilms’ tumor (WT) is the most common childhood renal tumor, and Curró et al. [14] study showed a significantly lower incidence of radial loops and whorls in WT patients. Gutjahr et al [15], showed a lower occurrence of digital arch patterns in affected cases and a slightly higher frequency of whorls in WT patients. A recent systematic review revealed an association of both qualitative and quantitative dermatoglyphic traits with several kidney diseases [7]. However, it was concluded that methodological issues may limit any interpretation of significant findings [7].

Symmetry is universal in the living world, and organisms display bilateral symmetries at every turn [16]. The extent to which the average individual departs from perfect symmetry is fluctuating asymmetry (FA) [17]. More specifically, FA has been defined as random differences between the right (R) and left (L) sides of a morphological trait [18]. Values of the right minus their corresponding values on the left are in fairly normal distribution with a mean of zero or close to zero, and any variance in the distribution of R-L difference is a measure of FA [16,17,18,19]. For distinction, directional asymmetry (DA) shows a significant departure from zero in the mean of R-L differences [16, 17]. Furthermore, anti-symmetry has a zero mean (or close to it) although the distribution around the mean is either platykurtic or bimodal. In essence, FA is then a population level measure of developmental instability, robustness, and developmental noise that is influenced by genetic factors along with their interaction with environmental stressors [17, 20,21,22]. Based upon these considerations, dermatoglyphic asymmetry can be used to measure developmental instability during a specific period of human fetal development [17]. Fluctuating asymmetry has been seen in dermatoglyphics as primarily concerned with the impact of environmental stress or noise that disrupts usual genetic expression [20, 22, 23].

We hypothesized that, if CKDu is partly due to a prenatal environmental exposure among genetically predisposed individuals, it might be associated with epidermal ridge formation and possibly lead to dermatoglyphic asymmetry. Thus, the present study was planned to investigate whether digital and palmar dermatoglyphics of CKDu patients in Sri Lanka are different from healthy people.

Methods

Study design

A case control study was carried out among CKDu patients from a CKDu endemic area and gender-matched two control groups; one group was selected from CKDu endemic region and another group from CKDu non-endemic region in Sri Lanka during 2014 to 2016. Ethical clearance was obtained from the Ethics Review Committee of Faculty of Medicine and Allied Sciences, Rajarata University of Sri Lanka (ERC/2013/31). Informed written consent to participate was obtained from all participants.

Study settings

Cases

The highest number of CKDu cases are reported from the Anuradhapura district in North Central Province (NCP) of Sri Lanka and considered as endemic for CKDu [24]. CKDu patients for the current study were selected from two CKDu clinics; kidney research and treatment unit of teaching hospital Anuradhapura and base hospital Padaviya (Fig. 1a). These patients were from all divisional secretariat (DS) divisions of the Anuradhapura districts, and a few patients were from the Adjacent Sri pura DS division in the eastern province.

Fig. 1

figure1

Geographical locations of cases (a) and controls (b)

Full size image

Control groups

Community-based control groups were selected from both CKDu endemic (Anuradhapura District) and non-endemic (Kandy District) regions (Fig. 1b). Endemic controls (EC) were recruited from 9 DS divisions from Anuradhapura district while non-endemic controls (NEC) were recruited from 7 DS divisions in Kandy districts.

Participant selection

All participants in cases and controls lived in their specified region for more than 10 years and belonged to the Sinhalese ethnic group for two generations.

Case definition

Cases included all CKDu patients that were diagnosed and followed up in two selected renal clinics. A medical doctor obtained the history and reviewed all medical records of CKDu patients and included only those patients who fulfilled the following criteria [24].

*

No past history of glomerulonephritis, pyelonephritis, renal calculi

*

Not on treatment for diabetes or (normal glycosylated hemoglobin (HbA1c); < 6.5%, normal fasting blood sugar)

*

If on treatment for hypertension, blood pressure below < 140/90 mmHg; if not on treatment for hypertension, blood pressure below < 160/100 mmHg.

We also excluded patients with a past history of snake bites when envenomation required hospital admission for Russell’s viper and cobra bites.

Control groups

A medical doctor obtained a detailed history and examination, and participants were excluded if they were found to have a past history of:

*

Kidney diseases (glomerular, tubular interstitial and cystic)

*

Hypertension or on treatment for hypertension

*

Diabetes or on treatment for diabetes

*

History of urinary infection or infectious disease affecting kidney (leptospirosis or schistosomiasis)

*

Prolong use of analgesics

*

Snake bite with envenomation, required hospital admission for Russell’s viper and cobra bites

If participants were found to a have systolic blood pressure > 120 mmHg or diastolic blood pressure > 90 mmHg, they were reexamined after a 30-min rest, if they still had high BP they were not included.

In addition, the following laboratory investigations were also carried out for control selection.

*

Estimated glomerular filtration rate (eGFR) based on CKD-EPI equation (> 60 mL/min/1.73m2) (If participant with eGFR between 60 and 90, urine albumin to creatinine ratio (ACR) and UFR values were checked for any abnormality and excluded if values were positive)

*

ACR (< 30 mg/g)

*

Normal UFR (without red cell cast)

*

HbA1C < 6.5%

The palms and fingers of both hands of CKDu patients and healthy controls were examined during initial screening to exclude any participant with a medical condition such as psoriasis or significant scaring that altered the epidermal ridge pattern.

Sample size calculation

The sample size was calculated by Kelsey formula for case control studies in OpenEpi software version 3.1 (OpenEpi, Atlanta, GA, USA) with α  =  0.05, power  =  80, ratio of cases to controls  =  1.0, hypothetical proportion of exposure among controls  =  50, and an odds ratio (OR) of 2 as a minimum difference between groups to be detected. This yielded a minimum sample size of 138 cases and 138 controls.

Obtaining fingerprints

Participants were asked to wash their hands with soap and water thoroughly to remove dirt and dry them before obtaining fingerprints. If hands were moist, each finger was wiped with rubbing alcohol and then allowed to dry to avoid light reflection of the recorded image.

First, all fingers were examined for level 1 details (patterns) by principle investigator. We used a newly developed photographic method of fingerprint recording to avoid participant discomfort and to prevent a low quality of recording associated with the conventional ink-based method. We used a Canon EOS 60D camera with a Canon EF-S 18-135 mm f/3.5–5.6 IS STM lens. Images of digit D2 (index finger) to D5 (little finger) were initially obtained while placing the hand on a stage at a right angle. Then, both palms were photographed while resting them on a stage right angles to camera. Finally, thumbs were photographed at a right angle while keeping them in an upright position. If there were extra-limital patterns, each finger or palm was photographed from their sides. All images for each participant were checked on the camera screen to avoid any recording of out of focus (blurred) images, and the procedure was repeated if necessary.

All images were stored securely in two separate hard drives. Firstly, all ten fingers of the images were screened for clarity. Then, the RAW images were opened within Canon’s Digital Photo Professional software and adjustments were made on brightness, white balance, contrast, gamma correction, and sharpness. These adjustments allowed clear visibility of ridges of the recorded prints. Then images were opened through Microsoft Paint software, and each finger pattern area was cropped. These were arranged in a single file and saved as TIFF (Tagged Image File Format) (> 300 DPI) (Fig. 2) to preserve the high-quality data.

Fig. 2

figure2

Plates with all digits and both palms of a person who has undergone analysis: D Digit. UL Ulnar loop. CPL Central pocket loop. PW plain whorl. RC ridge count. H higher count. L lower count.

Full size image

Dermatoglyphics classification

The digital fingerprint patterns were classified into eight types: ulnar loop (UL), radial loop (RL), plain arch (PA), tented arch (TA), plain whorl (PW), double loop (DL), central pocket loop (CPL), and accidental whorl (A), according to the pattern classification described by the Federal Bureau of Investigation, USA (Fig. 3) [25]. Any pattern that could not be classified was marked as undefined. Patterns that were in severely scared areas where details could not be observed or where amputation had occurred were marked as missing.

Fig. 3

figure3

The digital fingerprint patterns: a Ulnar loop. b: Radial loop. c Plain arch. d Tented arch. e Plain whorl. f Double loop. g Central pocket loop. h Accidental whorl

Full size image

Palmar dermatoglyphics variables were classified according to the Penrose Topological classification (Fig. 4) [11, 26]. Loops are specified according to the configuration area in which it occurs along with the direction of its core [26]. The configuration areas roughly correspond to fetal mounds [26]. Only true patterns, loops, and triradii are included in the description [11, 26]. Loops are designated by roman numerals according to the configuration area in which they were located and to the main direction of their cores, either distal (peripheral) or proximal (central) [11, 26]. Triradii were termed according to letters of the English alphabet [11, 26].

Fig. 4

figure4

Palmar dermatoglyphics variables: a, b Loop patterns. c Triradii

Full size image

Ridge counting and pattern intensity

Finger ridge count was defined as the number of ridges which intersect or touch a straight line drawn from the central point of a triradius to the center or core of the adjacent pattern [11]. Two ridges that result from a bifurcation of a single epidermal ridge and both cross the straight line are counted (Fig. 5). Any ridges that were close to the straight line without touching it were excluded [11]. Loop patterns have one ridge count while whorls usually have two ridge counts. Arches and other similar configurations that are not true patterns have a zero ridge count. When there was a missing finger on one hand, the ridge count on the corresponding finger of the other hand was inserted based on the considerable symmetry for this trait (this was done for three fingers in cases, two fingers in EC, one in NEC) [11].The sum of the largest ridge count on all ten fingers was defined as total ridge count (TRC) [11]. On the palm, the number of ridges that crossed a straight line connecting triradii “A” and “B” was defined as the A-B ridge count (A-B RC) [11].The pattern intensity index (PII) was calculated using the formula: PII = {(2 × % whorl + % loop) ÷ 10} [27]. PII essentially refers to the complexity of finger patterns in a specified population [6].

Fig. 5

figure5

Finger ridge counting: a Triradius. b Core

Full size image

Statistical analysis

The main demographic and clinical data were summarized using descriptive statistics.

The measure of central tendency used was the mean, and the measure of spread used was standard deviation (SD). Digital patterns were presented separately for each finger. The frequencies of the patterns were compared between the two groups using the chi-square test or the Fisher exact probability test. A P value of < 0.05 was considered statistically significant. All data were analyzed using GraphPad Prism 7.04 and SPSS® Version 20.

Directional asymmetry (DA) of Dermatoglyphics

The main aim of assessing DA is to identify whether one side was significantly larger than the other on average. We used the factorial ANOVA test to find out if significant differences existed between mean R and mean L of each homologous finger in a sample of individuals relative to the between-sides variation after accounting for measurement error. It also tested the significance of no DA and overall trait size variation among individuals. Antisymmetry artificially inflates the values of all FA indices. If a trait exhibits anti-symmetry, some portion of the between-sides variation may have a genetic basis; hence, the between-sides variance may not purely be a product of developmental noise [28]. Skewness and kurtosis statistics were used to describe any departure from normality.

Fluctuating asymmetry (FA) of Dermatoglyphics

Qualitative FA Fingerprint patterns are not morphometric traits that have a readily measurable size and shape. However, pattern types tend to be identical on homologous fingers. Hence, the degree of pattern discordance can be used as a measure of fluctuating asymmetry [29]. Therefore, the bilateral asymmetries of the finger print patterns were assessed using pattern discordance [17, 29, 30]. For example, if the R and L thumb of a person have the same pattern such as an UL on both thumbs, then the patterns are concordant. If the digital patterns differ, for instance, an UL on the right thumb and a PW on the left thumb, then the patterns are discordant.

Quantitative FA Based on standard recommendations, we used two methods to calculate fluctuating asymmetry [18]. These two methods were used by Saha et al. [31] during their study on dermatoglyphic asymmetry in psychosis.

FA calculated using correlation method

The comparisons were made in higher ridge counts between each homologous fingers of the R and L hands using Pearson product-moment correlation coefficients (r) [32]. The difference in r value of ridge count between cases and controls was calculated using Fisher’s z-transformation [31]. The square of the r of the two variables is a measure of their common variance and 1-r2 is an estimate of error variance (coefficient of indetermination) [31, 33] This is an estimate of their unshared variance and therefore of fluctuating asymmetry [29]. This method is not affected by DA [18].

Adjusted FA difference method

FA was obtained by subtracting the L hand ridge count from the R hand count of each homolog finger and dividing by the average of the sum of ridge count in both hand, and then taking the absolute value of the quotient [34]. The presence of DA artificially inflates the values of the FA difference scores. Therefore, if a significant DA was identified, FA difference scores for all measures were adjusted by subtracting the average mean R minus L difference (mean (R-L)/2) from the side with the larger mean and adding it to the smaller side of all individuals in the sample [35].

Reliability of observation

Intra-observer reliability for pattern classification was assessed by comparing observations made by the same observer during initial and later evaluations. Thumbprint patterns of 50 individuals were initially examined and recorded by the principle investigator. The same 50 individuals’ patterns were shuffled, and the thumbprints were examined and recorded a second time one month later by the principle investigator (PI). Similarly, palm print patterns were scored on two occasions. PI was blinded to his previous observation. Intra-observer reliability agreement for the pattern type on each finger was assessed with Cohen’s kappa (κ) coefficient of agreement. Interpretation of κ was based on the criteria of Landis and Koch [36].

Intra-rater reliability was also assessed by comparing two ridge counting made by the same rater during the first and second evaluations. The D1 ridge counts of 50 individuals were initially recorded by the principle investigator. These 50 individuals’ patterns were then shuffled, and their ridge counts were recorded a second time 1 month apart by the principle investigator. A similar procedure was used for palmar A-B RC. The reliability analysis was assessed by calculating the intraclass correlation coefficient (ICC) for the two recordings [37].

Results

Sample size and composition

A total of 180 (90 females and 90 males) CKDu cases were selected. The mean age of females was 60 years (SD = 10), and the mean age of males was 61 years (SD = 10). The EC consisted of 180 (90 females and 90 males) participants. The mean age of females was 44 years (SD = 10), and the mean age of males was 44 years (SD = 11). The NEC consisted of 180 (90 females and 90 males) participants; the mean age of females was 46 years (SD = 14) and mean age of males was 48 years (SD = 14).

Reliability of the dermatoglyphics analysis

There was almost perfect agreement between the two observations in both digital dermatoglyphics (κ = 0.932 [95% CI, 0.84 to 1.02], p < .0005) and palmar dermatoglyphic patterns (κ = 0.912 [95% CI, 0.82 to 1.01], p < .0001). A nearly exact agreement also was observed between the two observations on ridge counting of D1 (ICC = 0.995 [95% CI, 0.992 to 0.997], p < .0001) and A-B RC on R palm (ICC = 0.997 [95% CI, 0.994 to 0.998], p < .0001).

Qualitative dermatoglyphics

In the following tables, results are presented with three probability values according to two-way tests between CKDu cases and EC, between CKDu cases and NEC, and also between EC and NEC.

Digital dermatoglyphics of males—Additional file 1: Table S1

There were no statistically significant differences in frequencies of digital patterns observed in the R and L hands between CKDu cases and control groups (EC and NEC). However, PW frequency was significantly higher on L D1 in EC compared to NEC.

Digital dermatoglyphics of females—Additional file 2: Table S2

There was a significantly higher percentage of UL pattern observed on L D5 in females between CKDu patients, while PW frequency was significantly higher on L D1 of the EC compared to CKDu cases.

When comparing CKDu cases with NEC, a significantly higher percentage of DL was observed on L D3 while a higher percentage of CPL was seen on L D2. The percentage of PW was significantly lower on L D2. In contrast, a significantly lower percentage of PA was observed on R D3 in CKDu cases.

A comparison of EC and NEC showed a significantly lower percentage of PA in EC on R D3, while a significantly higher PW frequency was found on L D4 of EC.

Palmar dermatoglyphics of males—Additional file 3: Table S3 and Additional file 4: Table S4

No significant differences between cases and controls were found regarding male frequencies of palmar loop patterns.

Regarding triradii, a significantly higher percentage of triradii a1 on the R hand was observed in both control groups and cases.

Conversely, there is a significantly higher percentage of triradii t on the L hand of cases versus NEC. However, there was also a significantly higher frequency on the L hand triradii t in EC when compared with the other control group NEC.

Palmar dermatoglyphics of females—Additional file 5: Table S5 and Additional file 6: Table S6

For the palm, a comparison of CKDu cases with NEC yielded these results; the palmar loop IIIT was found to be significantly more common on the R hand of cases. For control groups, there were three palmar patterns (II and H on L hand, and IV on R hand) in which EC consistently showed significantly higher loop frequencies than NEC.

With regard to triradii, Triradii a1 was significantly less frequent in cases against both EC and NEC on the R hand while triradii t11 was significantly more frequent in cases versus NEC for both hands. Triradius c1 showed a mixed result, CKDu cases had a significantly lower c1 frequency than EC for both hands, but EC also had a significantly higher c1 frequency than NEC for the L hand.

FA of patterns—Table 1

Fluctuating asymmetry of pattern discordance was found to be significant on D3 among males and on D2 among females when comparing CKDu cases with NEC.

Table 1 Fluctuating asymmetry of pattern discordance

Full size table

Quantitative dermatoglyphics—Table 2

General results of quantitative variables of the between cases and control groups descriptive statistics and significance tests are shown in Table 2. As can be seen, only A-B RC was found to be significantly different in the paired case control tests, in which CKDu cases for both sexes had higher values, or higher palmar ridge counts, than EC, as well as for NEC females, and marginally non-significant in males. There was also a fairly consistent tendency for CKDu cases to have higher PII values along with larger finger ridge counts than control groups. This would weakly signal more complex patterns developed in CKDu patients.

Table 2 Quantitative variables of cases and controls

Full size table

Asymmetry of pattern distributions—test for normality

Since platykurtosis was not present in the distributions of right to left ridge count differences, the occurrence of anti-symmetry is unlikely. Thus, there is some support for interpreting FA in terms of environmental effect rather than due to genetic effect.

Directional asymmetry—Additional file 7: Table S7

Significant DA findings were observed for females on D1RC in both EC and NEC groups, and also on D1 RC in males of CKDu cases.

Fluctuating asymmetry—difference method Table 3

Significant differences in FA were observed for D2 RC in cases versus both control groups in females. In females, TRC was also significant in the cases versus EC test, along with that for A-B RC for the cases versus NEC test. For males, significant differences can be seen with regard to CKDu cases for A-B RC and the EC group, and for D3 RC, TRC, and A-B RC in the tests against the NEC control group. These results signal a moderate degree for distinguishing between cases and controls. However, significant differences also were found for six variables in tests between EC and NEC control groups, namely, for D5 RC and TRC in females and D3 RC, D4 RC, TRC, and A-B RC in males. Consequently, there appears to be nearly an even mix of FA differences between CKDu and control groups, that is, for the FA difference method, CKDu shows greater FA values in half of the significant test results, while the control groups themselves show higher FA in the remaining half.

Table 3 Fluctuating asymmetry—difference method

Full size table

Fluctuating asymmetry correlation method—Table 4

The FA-correlation method only identified significant differences of FA in D3 RC and D4 RC between cases and NEC in males, and in D2 RC in females. Males alone accounted for all of the significant FA-correlation tests of EC versus NEC control groups for the three variables of D3 RC, D4 RC, and A-B RC.

Table 4 Fluctuating asymmetry—correlation method

Full size table

The summary of positive findings—Table 5

There were many significant dermatoglyphic variables observed between cases and NEC as well as between EC and NEC. Importantly, significant FA results of the qualitative variables were only observed between cases and NEC. There were significant FA (FA correlation method) results observed in D3 and D4 of the males between cases and NEC, and also between EC and NEC. Similarly, significant FA (FA-adjusted method) results were observed in D3, TRC, and ABR of males between cases and NEC, as well as between EC and NEC.

Table 5 Summary of positive findings

Full size table

Discussion

This study showed that several qualitative dermatoglyphic variables had significant association with CKDu. Also, the FA of pattern discordance (R vs L hands) between CKDu cases and control group were significant in several digits. It is of interest to note that the overall digital dermatoglyphic pattern frequencies in each of the three groups were in line with those of the Sinhalese population [27]. There were few previous studies that provided evidence to support an association between qualitative dermatoglyphics and kidney diseases [7].

In our study, the FA of the ridge count was found to be significant in several digits as well as for A-B RC and TRC. Here is our basis for interpreting these findings. The development of dermal ridge patterns is controlled by the process of initial appearance and later regression of volar pads [6]. The development of ridges occurs at the dermal-epidermal junction [38]. The formation of volar pads first appear on the fingertips on the sixth to seventh week of fetal development and are prominent over subsequent weeks [6, 38]. Volar pads start to diminish from the fifth month and disappear by the sixth month [6, 38]. According to Bonnevie [6, 39] the position and size of the volar pad is responsible for pattern configuration. Thus, a small, low pad forms arch patterns, while an elevated, large symmetrical-shaped pad forms whorls, and asymmetrical pads form loops.

Importantly, unknown environmental factors could be a common cause for exerting stress on the development of both ridge patterns and kidneys. However the differences, there were few previous studies that provided evidence to support an association between qualitative dermatoglyphics and kidney diseases [7]. In our study, the FA of the ridge count was found to be significant in several digits as well as for A-B RC and TRC. Here is our basis for interpreting these findings. The development of dermal ridge patterns is controlled by the process of initial appearance and later regression of volar pads [6]. The development of ridges occurs at the dermal-epidermal junction [40]. The formation of volar pads first appear on the fingertips on the sixth to seventh week of fetal development and are prominent over subsequent weeks [6, 40]. Volar pads start to diminish from the fifth month and disappear by the sixth month [6, 40]. According to Bonnevie [6, 38], the position and size of the volar pad are responsible for pattern configuration. Thus, a small, low pad forms arch patterns, while an elevated, large symmetrical-shaped pad forms whorls, and asymmetrical pads form loops.

Importantly, unknown environmental factors could be a common cause for exerting stress on the development of both ridge patterns and kidneys. However, the difference for TRC in CKDu cases is only observed between EC female and NEC males. TRC has been shown to be less influenced by developmental insult arising from environmental factors [12].

The triradii a1 variable was less evident in CKDu cases in both genders when compared to both control groups. Among males who lived in the endemic region, FA of A-B RC was found to be significant. Both of these variables occur in interdigital area two (ID 2). Triradii a1 is an accessory pattern that appears in the vicinity of triradii A, while A-B RC is a measure of the size of the ID 2 [11]. With regard to developmental timing, the fetal pad for the ID 2 area appears first followed by the pads on fingers, and palmar ridges form earlier and develop over a longer period than digital ridges [39, 41]. The A-B RC asymmetry is reported to be ideal for indexing developmental canalization [41, 42]. The A-B RC could be highly sensitive to environmental noises, and these results in poor canalization in CKDu cases who are exposed to unknown environmental factors during development. Environmental stressors during the second-trimester were proposed to underlie an association between familial schizophrenia and dermatoglyphics [43]. Similarly, the buffering capacity for environmental stress could be lessened in CKDu patients. Therefore, these dermatoglyphic variables could be useful in the early detection of communities who are at risk of developing CKDu.

Could certain environmental factors simultaneously affect the development of CKDu and dermatoglyphics in the North Central Province? Environment has been convincingly shown to have had an impact on the development of dermatoglyphics during the prenatal period [13]. Furthermore, several environmental factors were found to be significantly associated with CKDu patients such as pesticides, heavy metals exposure, and contaminated drinking water [4]. Quite possibly, fetal exposure to these identified agents or as yet unknown factors could alter the development of both dermatoglyphics and kidneys that make them more vulnerable to disease later in life.

A study by Kahn et al. [44] showed L and R hand RC differences between D1 and D5 in the offspring of mothers who were exposed to hunger during gestation. Further, they have showed that the same dermatoglyphic variables were significantly associated with diabetes in offspring who were exposed to environmental stress during their gestation [23]. Related to kidney development, maternal dietary imbalance was reported to be a cause of an increased risk for insufficient renal function [45]. These studies clearly indicate the potential application of dermatoglyphic markers as disease risk factors.

Previous studies have reported significant DA in thumbs [31, 46], and we assessed DA for all quantitative variables. We used adjusted FA difference that overcame the confounding effect of DA [18]. This would make FA a suitable variable for investigating possible environmental effects on dermatoglyphic development. Furthermore, we were able to recruit participants who belonged to the Sinhalese ethnic group who were identified as Sinhalese for at least two generations without miscegenation. The analyses of prints were done in sex-wise subcategories thereby removing any sex-related confounders. Measurement error could have been a limitation of this study, although it was minimized because prints were examined by a single rater under strict supervision.

Despite our careful sampling procedure and methodology considerations, there were several variables found to be statistically significant between the control groups, EC, and NEC. We did exclude the presence of CKD in those individuals in EC through both clinical and laboratory investigations. However, the overall prevalence of CKDu in Anuradhapura is 15.1%. Therefore, there could be individuals in the EC group who could have been at risk for developing CKDu in the future. Therefore, our control group from EC might not have been an ideal complement in testing dermatoglyphic variables with the other control group from NEC. In essence, since our study provides evidence that dermatoglyphics of some individuals in EC could have been altered similar to that of CKDu cases, there is a possibility that fetal developmental insult from unknown underlying causative factors of CKDu also occurred in some individuals in the EC group.

In summing up this report, some measures of digital and palmar pattern variation, and also some of the derived measures of directional and fluctuating asymmetry, were found to be significant. What then, is the application for these results? Renal damage was evident among children living in CKDu endemic region [47]; therefore, it is vital to develop methods to identify those individuals and populations who are at risk of developing CKD later in life. Furthermore, methods should be cost effective and easy to administer. Since dermatoglyphics fully develop in early fetal life and thereafter remain unchanged, these variables can be used to identify individual children or a group of children who might beat greater risk of developing CKD in their future.

When dermatoglyphic findings are combined with clinical features and results of other investigations, they can strengthen a risk of disease diagnosis so that preventive methods can commence at an early age. The photographic method for taking dermatoglyphic prints can be conducted as a point-of-care appointment, where it would be of low cost, rapid, and non-invasive, and thus appropriate for use in low socioeconomic countries that could have less advanced diagnostic facilities. This method is so convenient and practical that it could be done during the neonatal period on a routine basis. Due to rapid advancement in technology, high-quality cameras are now available in mobile phones that could be used even in busy clinics.

Future studies should be directed toward developing methods to predict future occurrence of CKDu in risk populations. Dermatoglyphic variables such as and A-B RC found in ID 2 were significant in CKDu patients. By following similar methods described by Acree [48], it is possible to identify ridge density in a predefined square in ID 2 in CKDu group, and then the probabilities of occurrence of CKDu for a given dermal ridge count could be based on Bayesian inference.

In conclusion, to some extent, dermatoglyphics of CKDu patients in this study were altered compared to non-diseased individuals. Therefore, they may be used as an easily accessible tool to assist in the early detection of groups of people who are at increased risk for developing CKDu. Of course, further studies with larger sample sizes are needed to verify and extend our findings so that eventually highly reliable dermatoglyphic biomarkers may be discovered.

Availability of data and materials

Please contact the corresponding author for data requests.

Abbreviations

A:

Accidental whorl

A-B RC:

A-B ridge count

ACR:

Albumin to creatinine ratio

CINAC:

Chronic interstitial nephritis in agricultural communities

CKD:

Chronic kidney disease

CKDu:

Chronic kidney disease of unknown etiology

CPL:

Central pocket loop

DA:

Directional asymmetry

DL:

Double loop

DS:

Divisional secretariat

EC:

Endemic control

eGFR:

Estimated glomerular filtration rate

FA:

Fluctuating asymmetry

HbA1c:

Hemoglobin A1 C

ICC:

Intraclass correlation coefficient

ID:

Inter digital

IS:

Image stabilizer

LH:

Left hand

NCP:

North central province

NEC:

Non-endemic control

OR:

Odds ratio

PA:

Plain arc

PII:

Pattern intensity index

PW:

Plain whorl

RH:

Right hand

RL:

Radial loop

SD:

Standard deviation

STM:

Stepper Motor

TA:

Tented arch

TIFF:

Tagged Image File Format

TRC:

Total ridge count

UL:

Ulnar loop

References

1. 1.

Levey AS, Atkins R, Coresh J, Cohen EP, Collins AJ, Eckardt K-U, et al. Chronic kidney disease as a global public health problem: approaches and initiatives – a position statement from Kidney Disease Improving Global Outcomes. Kidney Int. 2007;72:247–59. https://doi.org/10.1038/sj.ki.5002343.

* CAS

* Article

* PubMed

* Google Scholar

2. 2.

Xie Y, Bowe B, Mokdad AH, Xian H, Yan Y, Li T, et al. Analysis of the Global Burden of Disease study highlights the global, regional, and national trends of chronic kidney disease epidemiology from 1990 to 2016. Kidney Int. 2018;94:567–81. https://doi.org/10.1016/j.kint.2018.04.011.

* Article

* PubMed

* Google Scholar

3. 3.

Weaver VM, Fadrowski JJ, Jaar BG. Global dimensions of chronic kidney disease of unknown etiology (CKDu): a modern era environmental and/or occupational nephropathy? BMC Nephrol. 2015;16:145. https://doi.org/10.1186/s12882-015-0105-6.

* Article

* PubMed

* PubMed Central

* Google Scholar

4. 4.

Jayasumana C, Orantes C, Herrera R, Almaguer M, Lopez L, Silva LC, et al. Chronic interstitial nephritis in agricultural communities: a worldwide epidemic with social, occupational and environmental determinants. Nephrol Dial Transplant. 2016;NDT Advanc October:gfw346. https://doi.org/10.1093/ndt/gfw346.

5. 5.

Jayasumana C, Orantes C, De Broe ME. Chronic interstitial nephritis in agricultural communities: Oxford University Press; 2018. https://doi.org/10.1093/med/9780199592548.003.0366_update_001.

6. 6.

Schaumann BA, Alter M. Dermatoglyphics in medical disorders. 1st ed. Berlin, Heidelberg: Springer-Verlag; 1976. https://doi.org/10.1007/978-3-642-51620-7.

* Google Scholar

7. 7.

Wijerathne BTB, Meier RJ, Salgado SS, Agampodi SB. Dermatoglyphics in kidney diseases: a review. Springerplus. 2016;5:290. https://doi.org/10.1186/s40064-016-1783-7.

* Article

* PubMed

* PubMed Central

* Google Scholar

8. 8.

Blackwell D. A dermatoglyphic investigation of selected skin disorders: Durham University; 1994. http://etheses.dur.ac.uk/5536/

9. 9.

Schaumann B, Johnson SB. Medical applications of dermatoglyphics. Prog Clin Biol Res. 1982;84:33–44 http://www.ncbi.nlm.nih.gov/pubmed/7100197.

* CAS

* PubMed

* Google Scholar

10. 10.

Babler WJ. Embryologic development of epidermal ridges and their configurations. Birth Defects Orig Artic Ser. 1991;27:95–112. https://pdfs.semanticscholar.org/cbb8/634eb80bcc413cd86282bb04dc000a5eee5e.pdf.

11. 11.

Loesch DZ. Quantitative dermatoglyphics: classification, genetics, and pathology. London: Oxford University Press; 1983.

* Google Scholar

12. 12.

Holt SB, Penrose LS. The genetics of dermal ridges: Thomas Springfield; 1968.

13. 13.

King S, Mancini-Marïe A, Brunet A, Walker E, Meaney MJ, Laplante DP. Prenatal maternal stress from a natural disaster predicts dermatoglyphic asymmetry in humans. Dev Psychopathol. 2009;21:343–53. https://doi.org/10.1017/S0954579409000364.

* Article

* PubMed

* Google Scholar

14. 14.

Curró V, Mastroiacovo P, Castello M, Romagnoli C, Mastrangelo R, Segni G. Palmar dermatoglyphics in Wilms’ tumor. Prog Clin Biol Res. 1982;84:385–91 http://www.ncbi.nlm.nih.gov/pubmed/6285387.

* PubMed

* Google Scholar

15. 15.

Gutjahr P, Wolffram T, Emmrich P. Dermatoglyphische Untersuchungen bei Kindern mit embryonalen Tumoren. Z Kinderheilkd. 1975;120:101–10. https://doi.org/10.1007/BF00445159.

* CAS

* Article

* PubMed

* Google Scholar

16. 16.

Graham JH, Raz S, Hel-Or H, Nevo E. Fluctuating asymmetry: methods, theory, and applications. Symmetry (Basel). 2010;2:466–540. https://doi.org/10.3390/sym2020466.

* Article

* Google Scholar

17. 17.

Graham J, Özener B. Fluctuating asymmetry of human populations: a review. Symmetry (Basel). 2016;8:154. https://doi.org/10.3390/sym8120154.

* Article

* Google Scholar

18. 18.

Palmer AR, Strobeck C. Fluctuating asymmetry: measurement, analysis, patterns. Annu Rev Ecol Syst. 1986;17:391–421. https://doi.org/10.1146/annurev.es.17.110186.002135.

* Article

* Google Scholar

19. 19.

Van Valen L. A study of fluctuating asymmetry. Evolution (N Y). 1962;16:125–42. https://doi.org/10.1111/j.1558-5646.1962.tb03206.x.

* Article

* Google Scholar

20. 20.

Arrieta MI, Criado B, Martinez B, Lobato MN, Gil A, Lostao CM. Fluctuating dermatoglyphic asymmetry: genetic and prenatal influences. Ann Hum Biol. 1993;20:557–63. https://doi.org/10.1080/03014469300002962

* CAS

* Article

* Google Scholar

21. 21.

Markow TA, Wandler K. Fluctuating dermatoglyphic asymmetry and the genetics of liability to schizophrenia. Psychiatry Res. 1986;19:323–8. https://doi.org/10.1016/0165-1781(86)90125-3.

* CAS

* Article

* PubMed

* Google Scholar

22. 22.

Martin NG, Jinks JL, Berry HS, Loesch DZ. A genetical analysis of diversity and asymmetry in finger ridge counts. Heredity (Edinb). 1982;48:393–405. https://doi.org/10.1038/hdy.1982.52.

* CAS

* Article

* PubMed

* Google Scholar

23. 23.

Kahn HS, Graff M, Stein AD, Lumey LH. A fingerprint marker from early gestation associated with diabetes in middle age: the Dutch Hunger Winter Families Study. Int J Epidemiol. 2009;38:101–9. https://doi.org/10.1093/ije/dyn158.

* Article

* PubMed

* Google Scholar

24. 24.

Jayatilake N, Mendis S, Maheepala P, Mehta FR. Chronic kidney disease of uncertain aetiology: prevalence and causative factors in a developing country. BMC Nephrol. 2013;14:180. https://doi.org/10.1186/1471-2369-14-180.

* Article

* PubMed

* PubMed Central

* Google Scholar

25. 25.

United States. Federal Bureau of Investigation. United States Department of Justice: The Science of Fingerprints Classification and Uses; 2006.

* Google Scholar

26. 26.

Penrose LS, Loesch D. Topological classification of palmar dermatoglyphics. J Ment Defic Res. 1970;14:111–28. https://doi.org/10.1111/j.1365-2788.1970.tb01106.x

* Article

* Google Scholar

27. 27.

Wijerathne BTB, Rathnayake GK, Adikari SC, Amarasinghe S, Abhayarathna PL, Jayasena AS. Sexual dimorphism in digital dermatoglyphic traits among Sinhalese people in Sri Lanka. J Physiol Anthropol. 2013;32:27. https://doi.org/10.1186/1880-6805-32-27.

* Article

* PubMed

* PubMed Central

* Google Scholar

28. 28.

Palmer AR, Strobeck C. Fluctuating asymmetry as a measure of developmental stability: implications of non-normal distributions and power of statistical tests. Acta Zool Fenn. 1992;191:55–70.

29. 29.

Mellor CS. Dermatoglyphic evidence of fluctuating asymmetry in schizophrenia. Br J Psychiatry. 1992;160:467–72. https://doi.org/10.1192/bjp.160.4.467.

* CAS

* Article

* PubMed

* Google Scholar

30. 30.

Ponnudurai R, Menon MS, Muthu M. Dermatoglyphic fluctuating asymmetry and symmetry in familial and non familial schizophrenia. Indian J Psychiatry. 1997;39:205–11. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2967116/.

31. 31.

Saha S, Loesch D, Chant D, Welham J, El-Saadi O, Fañanás L, et al. Directional and fluctuating asymmetry in finger and a-b ridge counts in psychosis: a case-control study. BMC Psychiatry. 2003;3:3. https://doi.org/10.1186/1471-244X-3-3

32. 32.

Holt SB. The correlations between ridge-counts on different fingers estimated from a population sample. Ann Hum Genet. 1959;23:459–60. https://doi.org/10.1111/j.1469-1809.1959.tb01487.x

* CAS

* Article

* Google Scholar

33. 33.

Micle S, Kobyliansky E. Sex differences in the intraindividual diversity of finger dermatoglyphics: pattern types and ridge counts. Hum Biol. 1988;60:123–34 http://www.jstor.org/stable/41463982.

* CAS

* PubMed

* Google Scholar

34. 34.

Sokal RR, Rohlf FJ. Biometry: The Principles and Practice of Statistics in Biological Research. 2nd ed. New York: W. H. Freeman; 1981. https://books.google.com.au/books?id=C-OTQgAACAAJ

* Google Scholar

35. 35.

Palmer RA. Fluctuating Asymmetry Analyses: A Primer. In: Markow TA, editor. Developmental Instability:Its Origins and Evolutionary Implications. London: Kluwer Academic Publishers; 1994. https://doi.org/10.1007/978-94-011-0830-0_26.

* Google Scholar

36. 36.

Landis JR, Koch GG. The measurement of observer agreement for categorical data. Biometrics. 1977;33:159. https://doi.org/10.2307/2529310.

* CAS

* Article

* PubMed

* Google Scholar

37. 37.

Rankin G, Stokes M. Reliability of assessment tools in rehabilitation: an illustration of appropriate statistical analyses. Clin Rehabil. 1998;12:187–99. https://doi.org/10.1191/026921598672178340.

* CAS

* Article

* PubMed

* Google Scholar

38. 38.

Mulvihill JJ, Smith DW. The genesis of dermatoglyphics. J Pediatr. 1969;75:579–89. https://doi.org/10.1016/S0022-3476(69)80453-1.

* CAS

* Article

* Google Scholar

39. 39.

Bonnevie K. Studies on papillary patterns of human fingers. J Genet. 1924;15:1–111. https://doi.org/10.1007/BF02983100.

* Article

* Google Scholar

40. 40.

Cummins H. The topographic history of the volar pads (walking pads; Tastballen) in the human embryo. Contrib Embryol. 1929;20:103–26.

41. 41.

Bogle AC, Reed T, Rose RJ. Replication of asymmetry of a-b ridge count and behavioral discordance in monozygotic twins. Behav Genet. 1994;24:65–72. https://doi.org/10.1007/BF01067929.

* CAS

* Article

* Google Scholar

42. 42.

Jantz RL, Webb RS. Dermatoglyphic asymmetry as a measure of canalization. Ann Hum Biol. 1980;7:489–93. https://doi.org/10.1080/03014468000004611.

* CAS

* Article

* Google Scholar

43. 43.

Bracha HS, Torrey EF, Gottesman II, Bigelow LB, Cunniff C. Second-trimester markers of fetal size in schizophrenia: a study of monozygotic twins. Am J Psychiatry. 1992;149:1355–61. https://doi.org/10.1176/ajp.149.10.1355.

* CAS

* Article

* PubMed

* Google Scholar

44. 44.

Kahn HS, Graff M, Stein AD, Zybert PA, McKeague IW, Lumey LH. A fingerprint characteristic associated with the early prenatal environment. Am J Hum Biol. 20:59–65. https://doi.org/10.1002/ajhb.20672.

* Article

* Google Scholar

45. 45.

Gluckman PD, Hanson MA, Cooper C, Thornburg KL. Effect of In Utero and Early-Life Conditions on Adult Health and Disease. N Engl J Med. 2008;359:61–73. https://doi.org/10.1056/NEJMra0708473.

* CAS

* Article

* PubMed

* PubMed Central

* Google Scholar

46. 46.

Jantz RL. On the levels of dermatoglyphic variation. Birth Defects Orig Artic Ser. 1979;15:53–61 http://www.ncbi.nlm.nih.gov/pubmed/552883.

* CAS

* PubMed

* Google Scholar

47. 47.

Agampodi SB, Amarasinghe GS, Naotunna PGCR, Jayasumana CS, Siribaddana SH. Early renal damage among children living in the region of highest burden of chronic kidney disease of unknown etiology (CKDu) in Sri Lanka. BMC Nephrol. 2018;19:115. https://doi.org/10.1186/s12882-018-0911-8.

* CAS

* Article

* PubMed

* PubMed Central

* Google Scholar

48. 48.

Acree MA. Is there a gender difference in fingerprint ridge density? Forensic Sci Int. 1999;102:35–44. https://doi.org/10.1016/S0379-0738(99)00037-7.

* CAS

* Article

* PubMed

* Google Scholar

Download references

Acknowledgements

We acknowledge all the support given by Professor Sisira Siribaddana, Dr. Thilini Agampodi, Dr. Rathis Deveraj, Professor Sarath Edirisinghe, Professor Malini Uduphille, Dr. Geetha Rathnayake, Mr. Indika Senavirathna, Dr. Kosala Weerakoon, Dr. Channa Jayasumana, and Dr. Anjana Silva. We specially thank Dr. Pradeep Rajasinghe, Dr. Anuradha Rathnayake, Dr. Anushitha Bandara, Mr. Sarath Wijerathne, Mr. Amila Ruwan, and Mr. Nalin. We also thank Dr. Saman Munasinghe, Dr. Rajeewa Dassanayake, Dr. Janyanatha Kahandawa, Dr. Chanaka Weerarathne, Dr. Pubudu Ranaweera, Dr. Geeth, Dr. NCD Ariyarathne, and Dr Shamitha Dassanayake. We sincerely thank all volunteer participants who agreed to participate in this study.

Funding

This work was partially funded by RJT/RP&HDC/2013/Med.&Alli.Sci./R/02.

Author information

Affiliations

1. Department of Community Medicine, Faculty of Medicine and Allied Sciences, Rajarata University of Sri Lanka, Saliyapura, 50008, Sri Lanka

* Buddhika Thilanga Bandara Wijerathne

* & Suneth Buddhika Agampodi

2. Department of Anthropology, Indiana University, Bloomington, IN, USA

* Robert John Meier

3. Department of Anatomy, Faculty of Medicine, University of Kelaniya, Ragama, Sri Lanka

* Sujatha Senadeera Salgado

Authors

1. Search for Buddhika Thilanga Bandara Wijerathne in:

* PubMed •

* Google Scholar

2. Search for Robert John Meier in:

* PubMed •

* Google Scholar

3. Search for Sujatha Senadeera Salgado in:

* PubMed •

* Google Scholar

4. Search for Suneth Buddhika Agampodi in:

* PubMed •

* Google Scholar

Contributions

BTBW conceived the research question, reviewed the data, performed the data analysis, and drafted the manuscript. BTBW, RJM, SS, and SBA participated in designing the study. BTBW collected the data. BTBW did the data analysis, and RJM, SBA, and SS the supervision. All authors have reviewed and approved the final manuscript.

Corresponding author

Correspondence to Buddhika Thilanga Bandara Wijerathne.

Ethics declarations

Ethics approval and consent to participate

This study was reviewed and approved by the Ethics Review Committee of Faculty of Medicine and Allied Sciences, Rajarata University of Sri Lanka (ERC/2013/31). Informed written consent to participate was obtained from all participants.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1. Digital dermatoglyphics of males.

Additional file 2: Table S2. Digital dermatoglyphics of females.

Additional file 3: Table S3. Palmar dermatoglyphics (loops) of males.

Additional file 4: Table S4. Palmar dermatoglyphics (triradii) of males.

Additional file 5: Table S5. Palmar dermatoglyphics (loops) of females.

Additional file 6: Table S6. Palmar dermatoglyphics (triradii) of females.

Additional file 7: Table S7. Directional asymmetry.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wijerathne, B.T.B., Meier, R.J., Salgado, S.S. et al. Qualitative and quantitative dermatoglyphics of chronic kidney disease of unknown origin (CKDu) in Sri Lanka. J Physiol Anthropol 39, 1 (2020) doi:10.1186/s40101-019-0207-0

Download citation

*

Received: 06 June 2019

*

Accepted: 11 September 2019

*

Published: 17 January 2020

*

DOI: https://doi.org/10.1186/s40101-019-0207-0

Keywords

* Chronic kidney disease

* Dermatoglyphics

* Fluctuating asymmetry

* Prenatal stress

* Sri Lanka