Y-STR analysis may in certain cases complement the autosomal STR markers in forensic, investigations, kinship testing and human identification. Hence, it would be informative to assess the probability of occurrence of the Y-STR haplotype in the Lebanese population where. This study aimed at estimating the Lebanese Y-STR allele and haplotype frequencies in 502 non-related males using the Y-filer Kit (Applied Biosystems) which includes the nine minimal Y-STR haplotype: DYS19, DYS390, DYS391, DYS392, DYS393, DYS389I/II, DYS385a/b plus 14 other Y-STR loci: DYS437, DYS438, DYS439, DYS456, DYS456, DYS448, Y-GATA-H4, DYS576, DYS570, DYS549, DYS643, DYS533 and DYS481. All 23 Y-STR loci were highly polymorphic with the marker DYS385a/b being the most polymorphic. 489 different haplotypes were defined where 476 (94.8%) carried a unique haplotype and the most common haplotypes appeared twice. This database has a discrimination capacity of 97.4% and a haplotype diversity of 0.9999%. The obtained data serves as an essential pre-requisite for using Y-chromosomal STR in routine forensic practice in the Lebanese population where the endogamy rate average is 88%.

Keywords: y-chromosome; y-str; haplotype; allele; frequency; lebanese population

Forensic investigation and human identification using DNA testing have come into popularity in the late 1980’s.^1-7 with the most commonly used polymorphic markers in forensic casework lay on the autosomal short tandem repeats (STR). In some case work, such as in rape cases or DNA mixture, autosomal STRs fail to be informative.^8,9 The amplification of Y chromosome STR markers provides an attractive alternative.^10,11 which will amplify only the male DNA in the mixture. There is a wide range of other applications to Y-STR profiling namely patrinlineal relationships, familial search.¹², disaster victim identification.¹³, studies of geographical or ethnic origins, archeogenetics.^14,15, genealogical studies, reconstruction of human history and investigation of cases involving mixtures.¹⁶ Certain limitations have constrained the use of Y-STR since it is always in a haploid state.², most of the polymorphisms lay in the non-recombining region of the Y-chromosome (NRY).¹⁷ and the NRY is inherited unchanged through paternal lineages unless a meiotic mutation occurs. Consequently, when a reference Y-STR profile matches a profile obtained from a crime scene trace, all members of the paternal lineage would probably match too, in addition to any male who shares a more distant paternal ancestry with the person to whom belongs the reference sample.¹⁸ The probability of sharing Y-STR profiles increases in isolated populations.¹⁹

According to Balding and Nichols.²⁰, it is possible that within a suspect’s sub-population, the crime scene DNA profile may be more common than in the wider population. Lebanon is a country with a population of around 4.5 million and a geographical area of 10,452 Km² situated on the eastern coast of the Mediterranean Sea. Several waves of immigration were observed in this area such as the Muslim expansion in the 7^th century, the Crusades between the 11^th and 13^th centuries and Ottoman Empire expansion in the 16^th century followed by the Armenian exode, French protectorate and the Palestinian exode in the 20^th century. All these historical events led to the diversity in terms of religious belongings to 18 different religious communities with an average rate of endogamy of 88%.²¹ Endogamy is a widely common practice and has an effect on differentiation resulting from genetic stratification by genetic relatedness. Consequently, population allele and haplotype frequencies are essential for a more accurate use of Y-STR since they provide the basis for random match probability calculations in forensics and human identification.² Inter-population variability seems to be more pronounced in the case of Y-chromosome than for the unlinked autosomal markers, which makes the definition of local databases essential for the use of Y-specific markers.²² In a crime or paterliniary study involving a suspect from an isolated sub-population group or geographically isolated region, the appropriate frequencies can be used in generating a random match probability.¹⁹ Such data are not available for the Lebanese population, and it would be highly informative to assess these frequencies where endogamy is a widely spread practice. The subject of the present work is to compile Y-chromosome allele and haplotype frequencies of natives from Lebanon using 23 Y-chromosome STR markers, calculate rate of occurrence and determine match probability in forensic casework.

Samples were collected from 502 non-related male individuals selected based on their geographic and religious distribution which represents the Lebanese population according to the Ministry of Interior and Municipalities in 2009 as described by El Andari et al..²³ Samples collected were either EDTA blood (n = 346) or buccal swabs (n = 156) (from right cheek, left cheek and tongue). DNA was extracted from whole blood leukocytes using the salting out method and from buccal swabs using a modified phenol-chloroform method. Samples were quantified using Nanodrop 2000 (Thermo Fisher Scientific Inc.) and diluted accordingly to approximately 1 ng/µl.

PCR amplification

DNA amplification of the 23 Y-STR loci was performed using two commercial kits: the Applied Biosystems Y-Filer^® multiplex PCR Amplification kit (Applied Biosystems, Foster City, CA) and the Promega PowerPlex^® Y23 System (Promega, Madison, USA). The 23 Y-STR systems include the 11 core loci recommended by the SWGDAM: DYS19, DYS385a/b, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS438, DYS439 and additional markers DYS437, DYS448, DYS456, DYS458, DYS635 and Y-GATA-H4. Amplifications typically contained 0.5-1.0 ng of extracted DNA. PCR reactions were carried out on GeneAmp PCR System 9700 (Applied Biosystems) using the cycling conditions as described in the manufacturers’ instructions.^24,25

DNA typing

Electrophoretic separation and detection were performed using the ABI PRISM^® 3130 Genetic Analyzer 4-capillary array system (ABI Prism 3130 Data Collection Software version 3.0) (Applied Biosystems, Foster City, CA). Size calling was performed using the GeneScan-500 Internal Lane Size Standard (LIZ-500) (Applied Biosystems) and CC5 Internal Lane Standard 500 (Promega, Madison, USA). Genotyping was performed by comparison with the provided allelic ladder and using Genemapper v4.0 (Applied Biosystems).

Statistical analysis

Y-STR data from the GeneScan^® software was transferred to an in-house software named Forensic Information Management System (FIMS) to estimate the allele frequencies for the 23 Y-STR systems. Gene diversity (GD) was calculated for each Y-STR according to the formula supplied by Nei.²⁶: GD = (n/n-1) (1-Σp_i²). Haplotype diversity (HD) was also computed with the same equation using haplotype frequencies instead of allele frequencies. Unique haplotypes (UH), random match probability (RMD) and discrimination capacity (DC).^27,2 were also calculated for the obtained data. Y-STR alleles are inherited in haplotypes, so their individual frequencies cannot be the product of the combined frequencies.²⁹ Haplotype frequency was obtained using the counting method.³⁰ Values were confirmed using Arlequin v3.5 software.³¹

Quality control

A proficiency testing quality control check was performed in conjunction with submission to the YHRD.org database.³²

Allele frequencies for 23 Y-STR markers were estimated in the Lebanese population (Table 1). 23 Y-STR systems were amplified, yet they were designated as 22 systems because of the duplicated DYS385 system which is a multi-copy Y-STR system that represents variations at two loci simultaneously, thus were analyzed as a phenotype.³³ Table 2 represents the allele frequencies for system DYS385a/b in the Lebanese population. Micro-variants (partial alleles) were observed in several occasions in system DYS458 (Table 1) whereby individuals exhibited an addition or deletion of 2 bp (alleles 16.2, 17.2, 18.2, 19.2, 20.2, 21.2, and 23.2). New and rare alleles were also documented (Table 3) according to published results in the NIST Standard Reference Database.³⁴ Null alleles and multiple peaks were reported at a number of loci (Table 3). The null alleles will be investigated and confirmed in the near future using a second set of primers. As observed previously.^35,36, DYS19 exhibits a relatively high frequency of duplications (reported to be 0.12%).

Diversity of Y-STR alleles

In a second step, the gene diversity values of the tested Y-STR systems in Lebanese and Caucasians males were compared (Table 4).^33,37,38 Systems were ranked according to their gene diversity values. When compared to the Caucasian population similar results were obtained for the two most polymorphic systems (DYS385a/b and DYS458). However, the system with the lowest gene diversity value in each of the Lebanese and Caucasian population groups were attributed to different markers (DYS392 in Lebanese population and DYS393 in Caucasian population), indicating that while some markers could be highly informative and polymorphic in Caucasian they could be of limited value in the analysis of the Lebanese population.

Taking into consideration the high rates of endogamy in the Lebanese population, the gene diversity for the different Lebanese sub-populations was calculated to attempt and assess whether genetic differences between different religious sub-populations existed. Results showed that a marker may have variable polymorphisms among these subpopulations (Table 5). For example DYS448 could discriminate 72.5% of the individuals in the Muslim Shiite subpopulation; however, it could only discriminate 59.1% of individuals belonging to the Druze subpopulation. Another variation was recorded in system DYS392 between Armenian Orthodox (73%) which greatly differed from Christian Orthodox (30%). The latter example particularly shows the effect of endogamous marriages whereby both communities were Orthodox, yet they still differed in terms of marker gene diversity since Armenian Orthodox and Christian Orthodox do not intermingle.

This demonstrated possible heterogeneity of allele frequencies across the different Lebanese sub-populations and is in agreement with previous population studies which showed that sub-populations exhibit greater differentiation at certain loci with the possibility of identifying new and unique alleles.³⁹

Lebanese Y-STR Haplotype Frequency

The haplotype diversity (HD) in the Lebanese population was 0.9999%. A total of 489 distinct haplotypes were observed in the total data set (n = 502) with 476 haplotypes being unique and 13 haplotypes observed more than once (Table 6). Different case scenarios were seen with the 13 common haplotypes shared by two individuals each: individuals either originated from the same religious community and geographical area, the same religious community but originated from different geographical areas or vice-versa or belonged to different religious communities as well as to different geographical origins. Results showed that Y-STR haplotypes were not restricted to members originating from the same family. Individuals from distinct geographical areas and/or religious communities shared a common haplotype. These cases could be explained by common ancestry whereby non-related individuals may have had a distant common ancestor and common haplotypes were transmitted from a distant common ancestor without mutation. Haplotype diversity (HD), discrimination capacity (DC), unique haplotypes (UH) and random match probability (RMP) calculations were performed to determine how common a Y-STR haplotype occurred in the population and how frequent a random match could occur between two non-related individuals (Table 7).

The occurrence of common haplotypes raises the question of how often a haplotype is present in the population and what are the chances of a match when used in forensic cases. The frequency of the most common haplotype (n=2) was 0.009% meaning that the probability of finding an individual sharing this haplotype is one in 111 (Table 8). When comparing a haplotype occurrence within the total population and sub-population datasets, the haplotype frequency and consequently the discrimination capacity varied. The haplotype frequency decreased in the sub-populations leading to a lower discrimination capacity and higher match probability. A match frequency of only one in 111 in the total population has a match frequency of one in 30 among Muslim Sunnites and of one in seven among Druze. The match probability varies greatly when using total population v/s sub-population datasets.

Allele and haplotype frequencies for the Lebanese population were estimated. The study showed that Y-STR markers exhibit possible genetic Y-chromosomal heterogeneity within and between populations and would seem to be very useful to trace back human evolutionary processes at a historical time-scale. While many Y-STR markers were highly discriminative in the Lebanese population; others were less discriminative such as system DYS392. In the future, if specific kits were to be developed for the Lebanese population and other closely related populations with similar socio-economic characteristics, we may omit this system (DYS392) and other similar low discriminative markers. Match probability differed when comparing the national database with the respective sub-population database, thus the question that raises itself which dataset to use in case of a match profile in a forensic case. In order to properly assess this issue, more samples from each subpopulation should be tested in order to determine whether allele and haplotype frequencies along with the match probability would change.

Evaluating possible Lebanese genetic sub-structures will be essential. If sub-structures exists, they should be accounted for when producing the strengths of the DNA profile evidence in Y-STR analysis whereby theta values will be calculated and incorporated in calculations. These criteria should be used to validate and evaluate Y-STR haplotype frequencies used in match probability calculations in forensic cases, human identification and kinship studies. Results in this study showed that there is a certain level of endogamy in the Lebanese population. Hence, for better assessment of the effect of endogamy a further Y-23 study will be performed on random villages with known high endogamous rates. This article follows the guidelines for publication of population data requested by the journal.⁴⁰

The authors wish to thank and acknowledge all the Lebanese volunteers for participating in this study, and Applied Biosystems and Promega for granting and supporting this project.

The author declares there are no conflicts of interest.

1. Edwards A, Civitello A, Hammond HA, Caskey CT. DNA typing and genetic mapping with trimeric and tetrameric tandem repeats. Am J Hum Genet. 1991;49(4): 746-756.
2. Weber JL, May PE. Abundant class of human DNA polymorphisms which can be typed using the polymerase chain reaction. Am J Hum Genet. 1989;44(3): 388-396.
3. Decker AE, Kline MC, Vallone PM, et al. The impact of additional Y-STR loci on resolving common haplotypes and closely related individuals. Forensic Sci Int Genet. 2007;1(2): 215-217.
4. Litt M, Luty JA. A hypervariable microsatellite revealed by in vitro amplification of a dinucleotide repeat within the cardiac muscle actin gene. Am J Hum Genet. 1989;44(3): 397-401.
5. Chakraborty R, Kidd K. The utility of DNA typing in forensic work. Science. 1991;254 (5039): 1735-1739.
6. Caskey CT, Hammond HA. Forensic use of short tandem repeats via PCR. Advances in forensic haemogenetics. (4th edn), Springer, Germany, p. 1992;18-25.
7. Roewer L Epplen JT. Rapid and sensitive typing of forensic stains by PCR amplifications of polymorphic simple repeat sequences in case work. Forensic Sci Int. 1992;53(2): 163-171.
8. Koblinsky LF, Liottie TF, Oeser Sweat J. Forensic DNA Analysis Methods. In DNA: Forensic and Legal Applications. John Wiley and Sons, USA, pp. 2005;1-114.
9. Parson W, Niederstatter H, Kochl S, et al. When autosomal short tandem repeats fail: optimized primer and reaction design for Y-chromosome short tandem repeat analysis in forensic casework. Croat Med J. 2001;42(3): 285-287.
10. Martin P, Albarran C, Garcia O, et al.  Application of Y-STR analysis to rape cases that cannot be solved by autosomal analysis. Prog in Forensic Gen. 2000;526-528.
11. Jobling MA, Tyler Smith C. Fathers and sons: the Y chromosome and human evolution. Trends Genet. 1995;11(11): 449-456.
12. Casanova M, Leroy P, Boucekkine C,  et al.  A human Y-linked DNA polymorphism and its potential for estimating genetic and evolutionary distance. Science. 1985;230(4732): 1403-1406.
13. Calacal GC, Delfin FC, Tan MM, et al.  Identification of exhumed remains of fire tragedy victims using conventional methods and autosomal/Y-chromosomal short tandem repeat DNA profiling. Am J ForensicMed Pathol. 2005;26(3):285-291.
14. Underhill PA, Shen P, Lin AA, et al.  Y chromosome sequence variation and history of human populations. Nat genet. 2000;26(3): 358-361.
15. Daniels DL, Hall AM, Ballantyne J. SWGDAM Development Validation of a 19-Locus Y-STR System for Forensic Casework. J Forensic Sci. 2004;49(4): 668-683.
16. White PS, Tatum OL, Deaven LL, et al. New male-specific microsatellite markers from the human Y chromosome. Genomics. 1999;57(3): 433-437.
17. Jobling MA, Pandya A, Tyler Smith C. The Y chromosome in forensic analysis and paternity testing. Int J Legal Med. 1997;110(3): 118-124.
18. Prinz M. Advantages and disadvantages of Y-short tandem repeat testing in forensic casework. Forensic Sci Rev. 2003;15(2): 191-198.
19. Butler JM. Biology, technology, and genetics of STR. In Forensic DNA Typing. (2nd edn), Elsevier Academic Press, UK, pp. 2005;1-660.
20. Balding DJ, Nichols RA. DNA profile match probability calculations: how to allow for population stratification, relatedness, database and single bands. Forensic Sci Int. 1994;64(2-3): 125-140.
21. Khlat M, A Khudr. Religious endogamy and consanguinity in marriage patterns in Beirut, Lebanon. Soc Biol. 1984;33(1-2): 138-145.
22. Kayser M, Krawezak M, Excoffier L, et al.  An extensive analysis of Y-chromosomal microsatellite haplotypes in globally dispersed human populations. Am J Hum Genet. 2001;68(4): 990-1018.
23. El Andari A, Othman H, Taroni F, et al. Population genetic data for 23 STR markers from Lebanon. Forensic Sci Int Genet. 2013;7(4): e108-e113.
24.  AmpFlSTR Yfiler PCR Amplification Kit User’s Manual. Applied Biosysytems, USA, pp. 2006;1-234.
25. PowerPlex® Y23 System Kit MSDS Manual.
26. Nei M. Molecular Evolutionary Genetics. (1st edn), Columbia University Press, USA. 2006.
27. Kayser M, Caglià A, Corach D, et al. Evaluation of Y-chromosomal STRs: a multicenter study. Int J Legal Med. 1997;110(3): 125-133, 141-149.
28. Schoske R, Vallone PM, Kline MC,  et al.  High throughput Y-STR typing of U.S. populations with 27 regions of the Y chromosome using two multiplex PCR assays. Forensic Sci Int. 2004;139(2-3): 107-121.
29. Charaborty R, Jin L, Zhong Y, et al. On allele frequency computation from DNA typing data. Int J Legal Med. 1993;106(2): 103-106.
30. Holland MM, Parsons TJ. Mitochondrial DNA sequence analysis - validation and use for forensic casework. Forensic Sci Rev. 1999;11(1): 21-50.
31. Excoffier L, Laval G, Schneider S. Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinform Online 2005;1: 47-50.
32. Willuweit S, Roewer L. Y chromosome haplotype reference database (YHRD): Update. Forensic Sci Int Genet. 2007;1(2): 83-87.
33. Mulero JJ, Chang CW, Calandro LM, et al.  Development and validation of the AmpFISTR Yfiler PCR amplification kit: a male specific, single amplification 17 Y-STR multiplex system. J Forensic Sci. 2006;51(1): 64-75.
34. Butler JM, reeder DJ. Short Tandem Repeat DNA Internet DataBase. National Institute of Standards and technology, USA. 2017.
35. Butler JM, Schoske R. Duplication of DYS19 flanking regions in other parts of the Y-chromosome. Int J Legal Med. 2004;183(3): 178-183.
36. Kayser M, Sajantila A. Mutations at Y-STR loci: Implications for paternity testing and forensic analysis. Forensic Sci Int. 2001;118(2-3): 116-121.
37. John MB, Decker AE, Vallone PM, et al. Allele frequencies for 27 Y-STR loci with U.S. Caucasian, African American, and Hispanic samples. Forensic Sci Int. 2006;156(2-3): 250-260.
38. Hanson EK, Berdos PN, Ballantyne J. Testing and evaluation of 43 “Noncore” Y chromosome markers for forensic casework applications. J Forensic Sci. 2006;51(6): 1298-1314.
39. Zhivotovsky LA, Ahmed S, Wang W, et al.  The forensic DNA implications of genetic differentiation between endogmaous communities. Forensic Sci Int. 2001;119(3): 269-272.
40. Lincolin P, Carracedo A. Publication of population data of human polymorphisms. Forensic Sci Int. 2000;110(1): 3-5.