Genetic variation in a Japanese population,using the multiplex 24 STRs analysis system

Allele frequency distributions for 24 short tandem repeat (STR) loci were determined using the PowerPlex^R Fusion System (Promega) in 407 Japanese samples. The most informative locus among the 22 STR loci, excluding Amelogenin and DYS391, was Penta E (power of discrimination (PD) = 0.98), while the least informative was TPOX(PD = 0.831). The 22 loci combined matching probability (MP) was calculated to be 4.13 × 10⁻²⁶. These parameters indicated the usefulness of this 24 STR analysis in forensic personal identification and parentage testing among Japanese population.     

Population genetic data of 38 autosomal InDels in San Basilio de Palenque,the first free town in America

In order to increase the information about Indels, we report allele frequencies and statistical parameters of forensic efficiency obtained typing a sample of 114 unrelated healthy individuals living in San Basilio de Palenque – Colombia using a panel of 38 autosomal InDels. No significant deviations from Hardy–Weinberg expectations were found except in the marker rs10629077 (p = 0.0002). The present database will be useful for forensic and paternity purposes for the region studied. Moreover, these additional markers can help forensic laboratories to solve parentage testing as well as to improve the analysis of degraded DNA samples.     

Forensic oral imaging quality of hand-held dental X-ray devices: Comparison of two image receptors and two devices

Recently, different portable hand-held and battery-powered dental X-ray units have become available. Especially for forensic odontological purposes, they offer diverse advantages such as for use in disaster areas and crime-scene locations as also in autopsy rooms and mortuaries. For any application, the most important feature of these hand-held devices is the delivered image quality. The aim of this study is to evaluate the radiographic image quality acquired by two portable X-ray devices in combination with two types of image receptors and to compare the findings with the image quality of a standard intra-oral X-ray device.Eleven samples consisting of eight teeth, two dry skeletal specimens and one formalin-fixed mandible part were mounted on blocks for standardised (re)positioning. Radiological images were acquired with two hand-held (AnyRay^® 60 kVp, 0.02–4.00 mAs and NOMAD^® 60 kVp, 0.023–2.277 mAs) and one wall-mounted (MinRay^® 60/70 kVp 0.14–22.4 mAs) X-ray device combined with two image receptor systems (VistaScan^® phosphor storage plate (PSP) and SIGMA^® M CMOS Active Pixel technology sensor). The effect of X-ray source-to-object distance (SOD) was checked at 20 cm in conjunction with object to image receptor distances (OIDs) of 0.8 and 2.5 cm. For each parameter setup, the exposure times were run from low till high. An expert consent statement was achieved by agreement of four expert observers selecting the optimal images based on a developed four point quality rating system. Next, a selection of the images was assembled in a set of 198 observation screens and scored by seven observers. The observation screens were designed to compare observer scores, relations between devices, receptors and OIDs and images obtained from the different devices at equal exposure levels (mAs). All results were statistically analysed.Radiological image quality was significantly higher for phosphor plate compared with the CMOS digital receptor system (p < 0.0001). Furthermore, a significantly superior image quality was obtained for OID = 0.8 than for OID = 2.5 (p = 0.039). A significant difference in image quality between the three devices was also established (p = 0.02). The present study demonstrated the feasibility of portable X-ray systems for forensic odontological applications based on rendering optimal image quality, provided an in vitro guideline of optimal parameter settings and offered a radiological image database usable in further research.     

Identification of a forensic case using microscopy and forensically informative nucleotide sequencing (FINS): A case study of small Indian civet (Viverricula indica)

The exhibits obtained in wildlife offence cases quite often present a challenging situation for the forensic expert. The selection of proper approach for analysis is vital for a successful analysis. A generalised forensic analysis approach should proceed from the use of non-destructive techniques (morphological and microscopic examination) to partially destructive and finally destructive techniques (DNA analysis). The findings of non-destructive techniques may sometime be inconclusive but they definitely help in steering further forensic analysis in a proper direction. We describe a recent case where a very small dried skin piece (< 0.05 mg) with just one small trimmed guard hair (0.4 cm) on it was received for species identification. The single guard hair was examined microscopically to get an indication of the type of species. We also describe the extraction procedure with a lower amount of sample, using an automated extraction method (Qiagen Biorobot EZ1^®) and PCR amplification of three mitochondrial genes (16s rRNA, 12s rRNA and cytochrome b) for species identification. Microscopic examination of the single hair indicated a viverrid species but the initial DNA analysis with 16s rRNA (through NCBI BLAST) showed the highest homology (93%) with a hyaenid species (Hyaena hyaena). However, further DNA analysis based on 12s rRNA and cytochrome b gene proved that the species was indeed a viverrid i.e. Viverricula indica (small Indian civet). The highest homology shown with a Hyaenid species by the 16s rRNA sequence from the case sample was due to lack of a 16s rRNA sequence for Viverricula indica in the NCBI data base. The case highlights the importance of morphological and microscopic examinations in wildlife offence cases. With respect to DNA extraction technology we found that automatic extraction method of Biorobot EZ1^® (Qiagen) is quite useful with less amount of sample (much below recommended amount).     

Forensic DNA laboratory automation – Principles and guidelines

There is a lack of clear guidelines for project managers, laboratory managers and forensic scientists on strategies for the automation of forensic DNA laboratory processes and operational implementation of new technologies. This is reflected in the failure rate of projects in the forensic DNA testing environment. We present a set of guidelines and concepts important for forensic laboratory automation. Some case studies from past projects are presented. These consist of partial (or modular) automation (n = 2) and full automated robotically integrated systems (n = 2).Technology Management principles and concepts are crucial to prevent failure of projects, e.g. early adoption of untried technologies, and organizational factors. The future of laboratory automation is modular until such time as new discontinuous technologies will replace the need of the traditional manual laboratory configuration in totality.     

Estimation of legal age using calcification stages of third molars in living individuals

The increased number of adolescents and young adults with unknown or inaccurately given date of birth is a current issue in justice and legal medicine. The objective of this study was to determine the extent to which third molar calcification stages assessed on panoramic X-rays could be useful as additional criteria for forensic age estimation in living individuals, focusing on the legally important ages 17 and 18.In a retrospective multi-center study, the developmental stage of each individual's third molar was analyzed using Demirjian's scale in 2360 cases. Additionally, sex, age and ancestry were assessed.Individuals with the lowest calcification stage of all present molars in stage H were ≥ 18 years with a likelihood of ≥ 99.05% in the female (n = 388), and ≥ 99.24% in the male (n = 482) population.The lowest calcification stage of all present third molars proved to be useful as an additional reliable criterion for the determination of an age ≥ 18 years.     

Development and validation of highly selective screening and confirmatory methods for the qualitative forensic analysis of organic explosive compounds with high performance liquid chromatography coupled with (photodiode array and) LTQ ion trap/Orbitrap mass spectrometric detections (HPLC-(PDA)-LTQOrbitrap)

An LTQ-Orbitrap FTMS is a new (hybrid) mass spectrometric (MS) analyzer. It allows for the acquisition of full scan MSⁿ (n-stage fragmentations, n = 1 − n) spectra with the linear ion trap detector (LTQ) at high speed and/or with the Fourier Transform-detector (Orbitrap) with ultra high mass resolution (> 60,000 at m/z < 400 amu) and high mass accuracy (≤ 1 ppm with internal calibration). In addition it may be coupled with liquid chromatography (LC) with photo diode array (PDA) detection.Two methods for the forensic screening and confirmation of all common trace explosives in post-blast residues have been developed on this instrument using atmospheric pressure chemical ionization (APCI). In one run, the nitrogen-containing explosives are analyzed with the combination of “LC-(PDA)-APCI(−)-LTQ MS²/Orbitrap FTMS” (Method 1). In another run, peroxide explosives are analyzed with “LC-APCI(+)-LTQ MS²/Orbitrap FTMS” (Method 2).The performance of both methods has been validated according to procedures defined in the EU COMMISSION DECISION implementing Council Directive 96/23/EC concerning the performance of analytical methods and the interpretation of results (DC 2002/657/EC) and other standards (NEN 17025 and NEN 7777). The methods are highly selective due to the simultaneous utilization of the Orbitrap FTMS and LTQ MS², both of which are highly selective detectors Tested explosive compounds can be detected in the molecular ion form by the Orbitrap analyzer with minimal mass interference in different matrices when using an extremely narrow mass tolerance detection window (≤ 2 ppm). The identification of a detected compound follows an identification point system. Experimental results show that almost all explosive compounds meet the confirmation criteria (minimum 4 points) required for the positive identification by the DC 2002/657/EC.     

A new 15 autosomal STRs multiplex for domestic horse genotyping

Horse genotyping has a wide range of applications such as identification, pedigree verification, parentage test, forensic investigation, population genetics, analysis of diversity, legitimate registration, among others. Following the recommendations of the International Society for Forensic Genetics (ISFG) regarding the use of non-human (animal) DNA in forensic genetic investigations we have developed a multiplex PCR system of 15 autosomal tetra-nucleotide STRs loci to Equus caballus. The system includes the newly described ECAC2, ECAC4, ECAC5, ECAC9, ECAC10, ECAC12, ECAC14, ECAC15, ECAC18, ECAC21, ECAC23, ECAC26, ECAC28, ECAC29 and ECAC30 loci (on chromosomes 2, 4, 5, 9, 10, 12, 14, 15, 18, 21, 23, 26, 28, 29 and 30, respectively). The polymorphism is in average 8 alleles per marker with a maximum of eleven and a minimum of five for the population studied. All markers were in Hardy–Weinberg equilibrium, except ECAC5 (p = 0.0007). The probabilities of paternity (W), exclusion (PE) and cumulative discrimination (PD) for all loci were greater than 0.9999. This work will contribute to the implementation of standardized horse genotyping systems in the forensic community and the horse industry.     

Characterisation of forward stutter in the AmpFlSTR® SGM Plus® PCR

PCR amplification of tetrameric short tandem repeats (STRs) can lead to Taq enzyme slippage and artefact products typically one repeat unit less in size than the parent STR. These back stutter or n ? 4 amplification products are low-level relative to the amplification of the parent STR but are widely seen in the forensic community where tetrameric STRs are employed in the generation of DNA profiles. To aid the interpretation of DNA mixtures where minor contributor(s) might be present in comparable amounts to the back stutter products, the typical amounts of back stutter generated have been well characterised and guidelines for interpretation are in place. However, further artefacts thought to be Taq enzyme slippage leading to products with one repeat unit greater than the parent sequence (n + 4 or forward stutter) or two repeats less (n ? 8 or double back stutter) also occur, but these have not been well characterised despite their potential influence in mixture interpretations. Here we present findings with respect to these additional artefacts from a study of 10,000 alleles and include guidelines for interpretation.     

Age estimation by pulp/tooth area ratio in canines: Cameriere's method assessed in an Indian sample using radiovisiography

Age estimation of an individual whether living or dead is an intimidating task in forensic investigations. Since teeth are more resistant to most peri- and post-mortem changes, they are frequently used for identification and age estimation when skeletal remains are in poor condition. However, most methods are destructive and warrant extraction of teeth which is not feasible in living individuals. Cameriere's et al. put forth a radiographic method of age estimation by pulp to tooth area ratio (AR) in canines and revealed a linear regression between age and the AR. In the present study, we estimated the AR in 456 canines (upper, lower and both) in an Indian sample (114 males and 114 females) using radiovisiography technique. Linear regression equations were derived for upper canine, lower canine and both using the AR to estimate chronological age. Additionally, the efficacy of these equations was also evaluated in younger age group (<45 years). The formulas derived, i.e., age = 96.795 ? 513.561x₁ (Eq. (1)) for upper canine, age = 88.308 ? 458.137x₂ (Eq. (2)) for lower canine and age = 99.190 ? 283.537x₁ ? 306.902x₂ + 400.873x₁x₂ (Eq. (3)) for both the canines were applied to predict the chronological age. The mean value of residuals using these regression equations ranged from 4.28 to 6.39 years with upper canine equation generally giving a precise result. When these equations were applied for younger ages (<45 years), the regression equation derived from both canines gave a better result (mean residual 2.70 years). Overall these equations were better able to predict the age in younger ages, i.e., up to 45 years.     

Multiple murder and criminal careers: A latent class analysis of multiple homicide offenders

PurposeTo construct an empirically rigorous typology of multiple homicide offenders (MHOs).MethodThe current study conducted latent class analysis of the official records of 160 MHOs sampled from eight states to evaluate their criminal careers.ResultsA 3-class solution best fit the data (?2LL = ?1123.61, Bayesian Information Criterion (BIC) = 2648.15, df = 81, L² = 1179.77). Class 1 (n = 64, class assignment probability = .999) was the low-offending group marked by little criminal record and delayed arrest onset. Class 2 (n = 51, class assignment probability = .957) was the severe group that represents the most violent and habitual criminals. Class 3 (n = 45, class assignment probability = .959) was the moderate group whose offending careers were similar to Class 2.ConclusionA sustained criminal career with involvement in versatile forms of crime was observed for two of three classes of MHOs. Linkages to extant typologies and recommendations for additional research that incorporates clinical constructs are proffered.     

Determination of amphetamine-type stimulants,ketamine and metabolites in fingernails by gas chromatography–mass spectrometry

A gas chromatography–mass spectrometry (GC–MS) method was developed and validated for the simultaneous qualification and quantification of methamphetamine (MA), amphetamine (AP), 3,4-methylenedioxy-N-methylamphetamine (MDMA), 3,4-methylenedioxy-N-amphetamine (MDA), ketamine (KET) and norketamine (NKT) in fingernails. Fingernail samples (20 mg) were washed with distilled water and methanol, digested with 1.0 M sodium hydroxide at 95 °C for 30 min, and then extracted with ethyl acetate. Extract solutions were evaporated to dryness, derivatized using heptafluorobutyric anhydride (HFBA) at 60 °C for 30 min, and analyzed by GC–MS. The linear ranges were 0.1–20.0 ng/mg for AP, MDMA and NKT, 0.2–20.0 ng/mg for MA and MDA, and 0.4–20.0 ng/mg for KET, with the coefficients of determination (r² ≥ 0.9989). The intra- and inter-day precisions were within 7.1% and 10.6%, respectively. The intra- and inter-day accuracies were ?10.9% to 0.8% and ?4.3% to 4.5%, respectively. The limits of detections (LODs) and the limits of quantifications (LOQs) for each analyte were lower than 0.094 ng/mg and 0.314 ng/mg, respectively. The recoveries were in the range of 72.3–94.9%. The average fingernail growth rates of two subjects for three years and six subjects for two months were 3.12 mm/month and 3.16 mm/month, respectively. The method proved to be suitable also for the simultaneous detection and quantification of MA, MDMA, KET and their metabolites in fingernails.     

Pressing need for national governmental recognition of forensic anthropology in South Africa as illustrated in a medico-legal case

Forensic anthropology in South Africa is well developed in the higher education sector, with advanced training and research programmes. Despite this and decades of academic involvement in casework, forensic anthropology still lacks a defined framework and mandate at a governmental level. Therefore, the involvement of forensic anthropologists’ expertise varies markedly between cases, provinces, and among various stakeholders within the country, to the detriment of dispensation of social and criminal justice. The lack of clearly defined guidelines for the rendering of the service was exemplified and demonstrated through a recent forensic case. Here, contextual information was absent, and the remains posed a challenge to analyse, ostensibly due to missing information. Numerous questions were raised during the analysis of the remains, and broader concerns about the investigative involvement of a forensic anthropologist within South African casework were brought to the fore. Through the analysis of this case, we describe the deductive processes that led to the formation of an opinion that the skeletal linear defects were the result of taphonomic changes. In addition, we highlight how these efforts where constrained and each step in the process unnecessarily hindered. Finally, we demonstrate the capacity and willingness of forensic anthropology practitioners to be involved, and how, without governmental support, it is a great potential lamentably untapped.     

SNP genotyping of forensic casework samples using the 52 SNPforID markers

The analysis of degraded DNA is one of the biggest challenges in forensic casework. SNPs, which can be amplified using small amplicons, have previously been successfully applied to the profiling of forensic evidence that could not be analyzed using conventional STRs. Here we selected the 52 SNPforID SNP markers, with amplicons that ranged in size from 59 bp to 115 bp, and used them to profile a range of casework samples from Malaysia. DNA degradation is a common problem in Malaysia due to the high temperatures and humidity. To carry out the study we modified the 52 SNPforID markers into four 13-plex SNaPshot assays to enable easier interpretation of profiles on the ABI PRISM^® 310 and 3500.Fifty-one crime samples comprising bloodstains on cloth, swabs, and a mat and 2 swabs of trace DNA from 10 crime scenes in Malaysia were profiled after DNA extraction using a phenol–chloroform method. The samples were also subjected to STR analysis using the Powerplex^® 16 system (Promega), which resulted in only 17 full profiles and 9 partial profiles; using SNPs, 36 full profiles and 5 partial profiles could be generated.     

The effect of frame rate on the ability of experienced gait analysts to identify characteristics of gait from closed circuit television footage

Forensic gait analysis is increasingly being used as part of criminal investigations. A major issue is the quality of the closed circuit television (CCTV) footage used, particularly the frame rate which can vary from 25 frames per second to one frame every 4 s. To date, no study has investigated the effect of frame rate on forensic gait analysis. A single subject was fitted with an ankle foot orthosis and recorded walking at 25 frames per second. 3D motion data were also collected, providing an absolute assessment of the gait characteristics. The CCTV footage was then edited to produce a set of eight additional pieces of footage, at various frame rates. Practitioners with knowledge of forensic gait analysis were recruited and instructed to record their observations regarding the characteristics of the subject's gait from the footage. They were sequentially sent web links to the nine pieces of footage, lowest frame rate first, and a simple observation recording form, over a period of 8 months. A sample-based Pearson product-moment correlation analysis of the results demonstrated a significant positive relationship between frame rate and scores (r = 0.868, p = 0.002). The results of this study show that frame rate affects the ability of experienced practitioners to identify characteristics of gait captured on CCTV footage. Every effort should therefore be made to ensure that CCTV footage likely to be used in criminal proceedings is captured at as high a frame rate as possible.     

Cytochrome b based genetic differentiation of Indian wild pig (Sus scrofa cristatus) and domestic pig (Sus scrofa domestica) and its use in wildlife forensics

The Indian wild pig (Sus scrofa cristatus) is a protected species and listed in the Indian Wildlife (Protection) Act, 1972. The wild pig is often hunted illegally and sold in market as meat warranting punishment under law. To avoid confusion in identification of these two subspecies during wildlife forensic examinations, we describe genetic differentiation of Indian wild and domestic pigs using a molecular technique. Analysis of sequence generated from the partial fragment (421 bp) of mitochondrial DNA (mtDNA) cytochrome b (Cyt b) gene exhibited unambiguous (> 3%) genetic variation between Indian wild and domestic pigs. We observed nine forensically informative nucleotide sequence (FINS) variations between Indian wild and domestic pigs. The overall genetic variation described in this study is helpful in forensic identification of the biological samples of wild and domestic pigs. It also helped in differentiating the Indian wild pig from other wild pig races. This study indicates that domestic pigs in India are not descendent of the Indian wild pig, however; they are closer to the other wild pig races found in Asia and Europe.     

Recovery of human DNA profiles from poached deer remains part 2: Improved recovery protocol without the need for LCN analysis

Although poaching is a common wildlife crime, the high and prohibitive cost of specialised animal testing means that many cases are left un-investigated. We previously described a novel approach to wildlife crime investigation that looked at the identification of human DNA on poached animal remains (Tobe, Govan and Welch, 2011). Human DNA was successfully isolated and amplified from simulated poaching incidents, however a low template protocol was required which made this method unsuitable for use in many laboratories. We now report on an optimised recovery and amplification protocol which removes the need for low template analysis.Samples from 10 deer (40 samples total — one from each leg) analysed in the original study were re-analysed in the current study with an additional 11 deer samples. Four samples analysed using Chelex did not show any results and a new method was devised whereby the available DNA was concentrated. By combining the DNA extracts from all tapings of the same deer remains followed by concentration, the recovered quantity of human DNA was found to be 29.5 pg ± 43.2 pg, 31 × greater than the previous study. The use of the Investigator Decaplex SE (QIAGEN) STR kit provided better results in the form of more complete profiles than did the AmpF?STR® SGM Plus® kit at 30 cycles (Applied Biosystems). Re-analysis of the samples from the initial study using the new, optimised protocol resulted in an average increase of 18% of recovered alleles. Over 17 samples, 71% of the samples analysed using the optimised protocol showed sufficient amplification for comparison to a reference profile and gave match probabilities ranging from 7.7690 × 10^? 05 to 2.2706 × 10^? 14.The removal of low template analysis means this optimised method provides evidence of high probative value and is suitable for immediate use in forensic laboratories. All methods and techniques used are standard and are compatible with current SOPs. As no high cost non-human DNA analysis is required the overall process is no more expensive than the investigation of other volume crime samples. The technique is suitable for immediate use in poaching incidents.     

Evaluation of the composition of street cocaine seized in two regions of Brazil

This work evaluates cocaine purity and the concentration ranges of adulterants and inorganic constituents for 31 street cocaine samples seized in two different regions of Brazil from July 2008 to May 2010. Cocaine and adulterants, such as caffeine, lidocaine and benzocaine, were quantified by Gas chromatography–mass spectrometry (GC–MS), and the inorganic constituents were determined by Inductively Coupled Plasma-Optical Emission Spectrometry (ICP-OES) and ion chromatography (IC). The cocaine concentrations in the samples seized in the Amazonas state (AM samples) ranged from 154 to 978 mg g^? 1, and these samples did not contain any of the adulterants studied. The cocaine concentrations in the samples seized in the Minas Gerais state (MG samples) ranged from 63.9 to 753 mg g^? 1. Caffeine was the main adulterant found in 76% of the MG samples, ranging in concentration from 5.5 to 645.3 mg g^? 1. Lidocaine was found in 66.7% of the MG samples, with concentrations ranging from 16.3 to 576.7 mg g^? 1. Benzocaine was found in only one MG sample, at a concentration of 84.8 mg g^? 1. Fourteen elements were identified by ICP-OES, and a wide variation was observed in the concentrations of Ca, Mg, Na, P, Al, Fe, Mn and Zn. Pearson Product–moment Correlations between the analytes allowed the constituents to be associated with the chemicals used in the manufacturing of cocaine and with some common diluents. The study of the purity of cocaine and the presence and concentration of adulterants and inorganic constituents is important because the latter can have deleterious effects on health.     

Freezing skeletal muscle tissue does not affect its decomposition in soil: Evidence from temporal changes in tissue mass,microbial activity and soil chemistry based on excised samples

The study of decaying organisms and death assemblages is referred to as forensic taphonomy, or more simply the study of graves. This field is dominated by the fields of entomology, anthropology and archaeology. Forensic taphonomy also includes the study of the ecology and chemistry of the burial environment. Studies in forensic taphonomy often require the use of analogues for human cadavers or their component parts. These might include animal cadavers or skeletal muscle tissue. However, sufficient supplies of cadavers or analogues may require periodic freezing of test material prior to experimental inhumation in the soil. This study was carried out to ascertain the effect of freezing on skeletal muscle tissue prior to inhumation and decomposition in a soil environment under controlled laboratory conditions. Changes in soil chemistry were also measured. In order to test the impact of freezing, skeletal muscle tissue (Sus scrofa) was frozen (?20 °C) or refrigerated (4 °C). Portions of skeletal muscle tissue (～1.5 g) were interred in microcosms (72 mm diameter × 120 mm height) containing sieved (2 mm) soil (sand) adjusted to 50% water holding capacity. The experiment had three treatments: control with no skeletal muscle tissue, microcosms containing frozen skeletal muscle tissue and those containing refrigerated tissue. The microcosms were destructively harvested at sequential periods of 2, 4, 6, 8, 12, 16, 23, 30 and 37 days after interment of skeletal muscle tissue. These harvests were replicated 6 times for each treatment. Microbial activity (carbon dioxide respiration) was monitored throughout the experiment. At harvest the skeletal muscle tissue was removed and the detritosphere soil was sampled for chemical analysis. Freezing was found to have no significant impact on decomposition or soil chemistry compared to unfrozen samples in the current study using skeletal muscle tissue. However, the interment of skeletal muscle tissue had a significant impact on the microbial activity (carbon dioxide respiration) and chemistry of the surrounding soil including: pH, electroconductivity, ammonium, nitrate, phosphate and potassium. This is the first laboratory controlled study to measure changes in inorganic chemistry in soil associated with the decomposition of skeletal muscle tissue in combination with microbial activity.     

Comparison of the three age estimation methods: Which is more reliable for Turkish children?

BackgroundThree atlases—the GÖK, the Greulich–Pyle (GP), and the Tanner–Whitehouse (TW3)—are used frequently for age determination in Turkey. The purpose of this study was to evaluate the applicability of these three methods related to the skeletal age assessment for Turkish adolescents.Materials and methodsThe conventional roentgenograms of the left hands and wrists, elbows, shoulders, and pelvises of 333 healthy Caucasian children (164 females, 169 males) who fit the study and the criteria of each atlas were obtained. The mean differences (± standard deviation [S.D.] in years) between the chronologic age (CA) and the skeletal age (BA), which were obtained by using each age estimation method, were calculated and tested using t-test.ResultsFor girls, the most accurate method was the TW3 (mean differences (d): ?0.21 (p < 0.05)), following by the GP (d: 0.66 (p < 0.001), and the GÖK (d: 2.99 (p < 0.001)). For boys, the most accurate method was the GP (d: ?0.02 (p > 0.05)), followed by the TW3 (d: ?0.18 (p < 0.05)) and GÖK (d: 1.05 (p < 0.001)).Discussion and conclusionsResults show that the TW3 (for girls) and the GP (for boys) methods are more appropriate than the GÖK atlas for estimating the BA. GÖK could be used for boys aged 11–14 years but it should not be used for girls.