A mathematical framework for genetic relatedness analysis involving X chromosome aneuploidies.

The unique features of the X chromosome can be crucial to complement autosomal profiling or to disentangle complex kinship problems, providing in some cases a similar or even greater power than autosomes in paternity/maternity investigations. While theoretical and informatics approaches for pairwise X-linked kinship analyses are well established for euploid individuals, these are still lacking for individuals with an X chromosome aneuploidy. To trigger the fulfilment of this gap, this research presents a mathematical framework that enables the quantification of DNA evidence in pairwise kinship analyses, involving two non-inbred individuals, one of whom with a non-mosaic X chromosome aneuploidy: Trisomy X (47, XXX), Klinefelter (47, XXY) or Turner (45, X0) syndrome. As previously developed for a regular number of chromosomes, this approach relies on the probability of related individuals sharing identical-by-descent (IBD) alleles at one specific locus and it can be applied to any set of independently transmitted markers, with no gametic association in the population. The kinship hypotheses mostly considered in forensic casework are specifically addressed in this work, but the reasoning and procedure can be applied to virtually any pairwise kinship problem under the referred assumptions. Algebraic formulae for joint genotypic probabilities cover all the possible genotypic configurations and pedigrees. Compared with the analyses assuming individuals with a regular number of chromosomes, complicating factors rely on the different possibilities for both the parental origin of the error (either maternal or paternal), and the type of error occurred (either meiotic or post-zygotic mitotic). These imply that a non-inbred female with Triple X or a male with Klinefelter syndrome may carry two IBD alleles at the same locus. Thus, and contrarily to what occurs for the standard case, IBD partitions depend not only on the kinship hypothesis under analysis but also on the genotypic configuration of the analyzed individuals. For some cases, parameters of interest can be inferred, while for others recommended values based on the available literature are provided. This work is the starting point to analyze X-chromosomal data under the scope of kinship problems, involving individuals with aneuploidies, as it will enhance the quantification of the DNA evidence not only in forensics but also in the medical genetics field. We hope it will trigger the development of approaches including other complicating factors, as a greater number of individuals, possibility of the occurrence of mutations and/or silent alleles, as well as the analysis of linked markers.

Recommendations of the DNA Commission of the International Society for Forensic Genetics (ISFG) on short tandem repeat sequence nomenclature.

The DNA Commission of the International Society for Forensic Genetics (ISFG) has developed a set of nomenclature recommendations for short tandem repeat (STR) sequences. These recommendations follow the 2016 considerations of the DNA Commission of the ISFG, incorporating the knowledge gained through research and population studies in the intervening years. While maintaining a focus on backward compatibility with the CE data that currently populate national DNA databases, this report also looks to the future with the establishment of recommended minimum sequence reporting ranges to facilitate interlaboratory comparisons, automated solutions for sequence-based allele designations, a suite of resources to support bioinformatic development, guidance for characterizing new STR loci, and considerations for incorporating STR sequences and other new markers into investigative databases.

The IPEFA model: An initiative for online training and education as applied by the International Society for Forensic Genetics.

The IPEFA model was developed for organizing online training and education events as applied by the International Society for Forensic Genetics (ISFG). It consists of five phases: 1) Input, 2) Preparation, 3) Execution, 4) Feedback, and 5) Assessment. This document details these phases and shows IPEFA's first practical application to the 2023 edition of the virtual ISFG Summer School. Through sharing the experiences, we aim to provide transparency and engage with potential participants and teachers to (virtual) training and education events as organized by the ISFG. The model may also be useful for others organizing (online) events. We have experienced that evaluation of events with input and feedback from both the (potential) participants and teachers is essential for successful training and education. This takes time which is limited in everyone's busy agenda's and may therefore not always be performed with the care it requires. Since these aspects are crucial, however, we aim to keep following the principles as outlined in the IPEFA model.

Impact of patrilocality on contrasting patterns of paternal and maternal heritage in Central-West Africa.

Despite their ancient past and high diversity, African populations are the least represented in human population genetic studies. In this study, uniparental markers (mtDNA and Y chromosome) were used to investigate the impact of sociocultural factors on the genetic diversity and inter-ethnolinguistic gene flow in the three major Nigerian groups: Hausa (n = 89), Yoruba (n = 135) and Igbo (n = 134). The results show a distinct history from the maternal and paternal perspectives. The three Nigerian groups present a similar substrate for mtDNA, but not for the Y chromosome. The two Niger-Congo groups, Yoruba and Igbo, are paternally genetically correlated with populations from the same ethnolinguistic affiliation. Meanwhile, the Hausa is paternally closer to other Afro-Asiatic populations and presented a high diversity of lineages from across Africa. When expanding the analyses to other African populations, it is observed that language did not act as a major barrier to female-mediated gene flow and that the differentiation of paternal lineages is better correlated with linguistic than geographic distances. The results obtained demonstrate the impact of patrilocality, a common and well-established practice in populations from Central-West Africa, in the preservation of the patrilineage gene pool and in the affirmation of identity between groups.

Microsatellites' mutation modeling through the analysis of the Y-chromosomal transmission: Results of a GHEP-ISFG collaborative study.

The Spanish and Portuguese Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) organized a collaborative study on mutations of Y-chromosomal short tandem repeats (Y-STRs). New data from 2225 father-son duos and data from 44 previously published reports, corresponding to 25,729 duos, were collected and analyzed. Marker-specific mutation rates were estimated for 33 Y-STRs. Although highly dependent on the analyzed marker, mutations compatible with the gain or loss of a single repeat were 23.2 times more likely than those involving a greater number of repeats. Longer alleles (relatively to the modal one) showed to be nearly twice more mutable than the shorter ones. Within the subset of longer alleles, the loss of repeats showed to be nearly twice more likely than the gain. Conversely, shorter alleles showed a symmetrical trend, with repeat gains being twofold more frequent than reductions. A positive correlation between the paternal age and the mutation rate was observed, strengthening previous findings. The results of a machine learning approach, via logistic regression analyses, allowed the establishment of algebraic formulas for estimating the probability of mutation depending on paternal age and allele length for DYS389I, DYS393 and DYS627. Algebraic formulas could also be established considering only the allele length as predictor for DYS19, DYS389I, DYS389II-I, DYS390, DYS391, DYS393, DYS437, DYS439, DYS449, DYS456, DYS458, DYS460, DYS481, DYS518, DYS533, DYS576, DYS626 and DYS627 loci. For the remaining Y-STRs, a lack of statistical significance was observed, probably as a consequence of the small effective size of the subsets available, a common difficulty in the modeling of rare events as is the case of mutations. The amount of data used in the different analyses varied widely, depending on how the data were reported in the publications analyzed. This shows a regrettable waste of produced data, due to inadequate communication of the results, supporting an urgent need of publication guidelines for mutation studies.

SNPtotree-Resolving the Phylogeny of SNPs on Non-Recombining DNA.

Genetic variants on non-recombining DNA and the hierarchical order in which they accumulate are commonly of interest. This variant hierarchy can be established and combined with information on the population and geographic origin of the individuals carrying the variants to find population structures and infer migration patterns. Further, individuals can be assigned to the characterized populations, which is relevant in forensic genetics, genetic genealogy, and epidemiologic studies. However, there is currently no straightforward method to obtain such a variant hierarchy. Here, we introduce the software SNPtotree v1.0, which uniquely determines the hierarchical order of variants on non-recombining DNA without error-prone manual sorting. The algorithm uses pairwise variant comparisons to infer their relationships and integrates the combined information into a phylogenetic tree. Variants that have contradictory pairwise relationships or ambiguous positions in the tree are removed by the software. When benchmarked using two human Y-chromosomal massively parallel sequencing datasets, SNPtotree outperforms traditional methods in the accuracy of phylogenetic trees for sequencing data with high amounts of missing information. The phylogenetic trees of variants created using SNPtotree can be used to establish and maintain publicly available phylogeny databases to further explore genetic epidemiology and genealogy, as well as population and forensic genetics.

Pitfalls and challenges with population assignments of individuals from admixed populations: Applying Genogeographer on Brazilian individuals.

The assignment of individuals to a population can be of importance for the identification of mass disaster victims or criminal offenders in the field of forensic genetics. This assignment is based on biostatistical methods that process data of ancestry informative markers (AIMs), which are selected based on large allele frequency differences between the populations of interest. However, population assignments of individuals with an admixed genetic background are challenging. Admixed individuals are genetic mosaics of chromosomal segments from the parental populations, which may lead to ambiguous or no population assignment. This is problematic since admixture events are a substantial part of human history. In this study, we present challenges of interpreting the evidential weight of population assignments. We used Genogeographer for likelihood ratio (LR) calculations and Brazilians as examples of admixed individuals. Brazilians are a very heterogenous population representing a three-way admixture between Native Americans, Europeans, and Africans. Ancestry informative markers were typed in a total of 589 individuals from Brazil using the Precision ID Ancestry Panel. The Brazilians were assigned to six metapopulations (East Asia, Europe, Middle East, North Africa, South-Central Asia, Sub-Saharan Africa) defined in the Genogeographer software and LRs were calculated if the AIM profile was not an outlier in all metapopulations and simulated two-way (1:1) admixtures of the six metapopulations. Population assignments failed for 55% of the samples. These samples had significantly higher genetic contributions from East Asia, South-Central Asia and Sub-Saharan Africa, and significantly lower genetic contributions from Europe. Most of the individuals with population assignments were assigned to the metapopulations of Middle East (58%) or North Africa (36%), followed by Europe (4%), South-Central Asia (1%), and Sub-Saharan Africa (1%). For 8% of the samples, population assignments were only possible when assignments to simulated two-way (1:1) admixtures of the six metapopulations were considered. Most of these individuals were assigned to two-way admixtures of North Africa, South-Central Asia, or Sub-Saharan Africa. Relatively low median likelihood ratios (LRs<1000) were observed when comparing population likelihoods for Europe, Middle East, North Africa, South-Central Asia, or simulated 1:1 admixtures of these metapopulations. Comparisons including East Asian or Sub-Saharan African populations resulted in larger median LRs (LR>1010). The results suggested that the Precision ID Ancestry Panel provided too little information and that additional markers specifically selected for sub-continental differentiation may be required for accurate population assignment of admixed individuals. Furthermore, a Genogeographer database with additional populations including admixed populations would be advantageous for interpretation of admixed AIM profiles. It would likely increase the number of population assignments and illustrate alternatives to the most likely population, which would be valuable information for the case officer when writing the case report.

The sequence of the repetitive motif influences the frequency of multistep mutations in Short Tandem Repeats.

Microsatellites, or Short Tandem Repeats (STRs), are subject to frequent length mutations that involve the loss or gain of an integer number of repeats. This work aimed to investigate the correlation between STRs' specific repetitive motif composition and mutational dynamics, specifically the occurrence of single- or multistep mutations. Allelic transmission data, comprising 323,818 allele transfers and 1,297 mutations, were gathered for 35 Y-chromosomal STRs with simple structure. Six structure groups were established: ATT, CTT, TCTA/GATA, GAAA/CTTT, CTTTT, and AGAGAT, according to the repetitive motif present in the DNA leading strand of the markers. Results show that the occurrence of multistep mutations varies significantly among groups of markers defined by the repetitive motif. The group of markers with the highest frequency of multistep mutations was the one with repetitive motif CTTTT (25% of the detected mutations) and the lowest frequency corresponding to the group with repetitive motifs TCTA/GATA (0.93%). Statistically significant differences (α = 0.05) were found between groups with repetitive motifs with different lengths, as is the case of TCTA/GATA and ATT (p = 0.0168), CTT (p < 0.0001) and CTTTT (p < 0.0001), as well as between GAAA/CTTT and CTTTT (p = 0.0102). The same occurred between the two tetrameric groups GAAA/CTTT and TCTA/GATA (p < 0.0001) - the first showing 5.7 times more multistep mutations than the second. When considering the number of repeats of the mutated paternal alleles, statistically significant differences were found for alleles with 10 or 12 repeats, between GATA and ATT structure groups. These results, which demonstrate the heterogeneity of mutational dynamics across repeat motifs, have implications in the fields of population genetics, epidemiology, or phylogeography, and whenever STR mutation models are used in evolutionary studies in general.

The multifaceted genomic history of Ashaninka from Amazonian Peru.

Despite its crucial location, the western side of Amazonia between the Andes and the source(s) of the Amazon River is still understudied from a genomic and archaeogenomic point of view, albeit possibly harboring essential information to clarify the complex genetic history of local Indigenous groups and their interactions with nearby regions,1,2,3,4,5,6,7,8 including central America and the Caribbean.9,10,11,12 Focusing on this key region, we analyzed the genome-wide profiles of 51 Ashaninka individuals from Amazonian Peru, observing an unexpected extent of genomic variation. We identified at least two Ashaninka subgroups with distinctive genomic makeups, which were differentially shaped by the degree and timing of external admixtures, especially with the Indigenous groups from the Andes and the Pacific coast. On a continental scale, Ashaninka ancestors probably derived from a south-north migration of Indigenous groups moving into the Amazonian rainforest from a southeastern area with contributions from the Southern Cone and the Atlantic coast. These ancestral populations diversified in the variegated geographic regions of interior South America, on the eastern side of the Andes, differentially interacting with surrounding coastal groups. In this complex scenario, we also revealed strict connections between the ancestors of present-day Ashaninka, who belong to the Arawakan language family,13 and those Indigenous groups that moved further north into the Caribbean, contributing to the early Ceramic (Saladoid) tradition in the islands.14,15.

Investigation of associations of European, African, Amerindian genomic ancestries and MC4R, FTO, FAIM2, BDNF loci with obesity-related traits in Rio de Janeiro, Brazil.

A complex web of causation is involved in adiposity, including environmental, social and genetic factors. We aimed to investigate associations between genetic factors such as ancestry and single nucleotide polymorphisms, and obesity-related traits in a sampled Brazilian population. A sample of 501 unrelated adults participating in 2013 at the longitudinal Pró-Saúde Study (EPS) in Rio de Janeiro, Brazil was selected. We analysed 46 AIM-InDels (insertion/deletion) as genetic ancestry markers and four single nucleotide polymorphisms located in the genes MC4R (rs17782313), FTO (rs9939609), FAIM2 (rs7138803) and BDNF (rs4074134), previously described as associated with obesity. The selected obesity-related markers were anthropometric parameters such as body mass index, waist circumference and waist-to-hip ratio, and body composition measurements namely body fat percentage, android fat mass and gynoid fat mass. The sample showed greater European ancestry (57.20%), followed by African (28.80%) and lastly Amerindian (14%). Our results suggest that the rs4074134 (BDNF) CC genotype was directly associated with gynoid fat mass, whereas body fat percentage, android fat mass and the anthropometric parameters seem not to be associated with neither ancestry nor the four polymorphisms in this population sample, most likely due to a stronger role of social, behavioural and environmental determinants.

Maternal ancestry and lineages diversity of the Santander population from Colombia.

Santander, located in the Andean region of Colombia, is one of the 32 departments of the country. Its population was shaped by intercontinental admixture between autochthonous native Americans, European settlers, and African slaves. To establish forensic databases of haplotype frequencies, the evaluation of population substructure is crucial to capture the genetic diversity in admixed populations. Total control region mitochondrial deoxyribonucleic acid haplotypes were determined for 204 individuals born in the seven provinces across the department. The maternal native heritage is highly preserved in Santander genetic background, with 90% of the haplotypes belonging to haplogroups inside A2, B4, C1, and D. Most native lineages are found broadly across the American continent, while some sub-branches are concentrated in Central America and north South America. Subtle European (6%) and African (4%) input was detected. In pairwise comparisons between provinces, relatively high FST values were found in some cases, although not statistically significant. Nonetheless, when provinces were grouped according to the principal component analysis results, significant differences were detected between groups. The database on mitochondrial deoxyribonucleic acid control region haplotype frequencies established here can be further used for populational and forensic purposes.

Tierra Del Fuego: What Is Left from the Precolonial Male Lineages?

Similar to other South American regions, Tierra del Fuego has an admixed population characterized by distinct ancestors: Native Americans who first occupied the continent, European settlers who arrived from the late 15th century onwards, and Sub-Saharan Africans who were brought to the Americas for slave labor. To disclose the paternal lineages in the current population from Tierra del Fuego, 196 unrelated males were genotyped for 23 Y-STRs and 52 Y-SNPs. Haplotype and haplogroup diversities were high, indicating the absence of strong founder or drift events. A high frequency of Eurasian haplogroups was detected (94.4%), followed by Native American (5.1%) and African (0.5%) ones. The haplogroup R was the most abundant (48.5%), with the sub-haplogroup R-S116* taking up a quarter of the total dataset. Comparative analyses with other Latin American populations showed similarities with other admixed populations from Argentina. Regarding Eurasian populations, Tierra del Fuego presented similarities with Italian and Iberian populations. In an in-depth analysis of the haplogroup R-M269 and its subtypes, Tierra del Fuego displayed a close proximity to the Iberian Peninsula. The results from this study are in line with the historical records and reflect the severe demographic change led mainly by male newcomers with paternal European origin.

Estimations of Mutation Rates Depend on Population Allele Frequency Distribution: The Case of Autosomal Microsatellites.

Microsatellites (or short-tandem repeats (STRs)) are widely used in anthropology and evolutionary studies. Their extensive polymorphism and rapid evolution make them the ideal genetic marker for dating events, such as the age of a gene or a population. This usage requires the estimation of mutation rates, which are usually estimated by counting the observed Mendelian incompatibilities in one-generation familial configurations (typically parent(s)-child duos or trios). Underestimations are inevitable when using this approach, due to the occurrence of mutational events that do not lead to incompatibilities with the parental genotypes ('hidden' or 'covert' mutations). It is known that the likelihood that one mutation event leads to a Mendelian incompatibility depends on the mode of genetic transmission considered, the type of familial configuration (duos or trios) considered, and the genotype(s) of the progenitor(s). In this work, we show how the magnitude of the underestimation of autosomal microsatellite mutation rates varies with the populations' allele frequency distribution spectrum. The Mendelian incompatibilities approach (MIA) was applied to simulated parent(s)/offspring duos and trios in different populational scenarios. The results showed that the magnitude and type of biases depend on the population allele frequency distribution, whatever the type of familial data considered, and are greater when duos, instead of trios, are used to obtain the estimates. The implications for molecular anthropology are discussed and a simple framework is presented to correct the naïf estimates, along with an informatics tool for the correction of incompatibility rates obtained through the MIA.

Testing the Ion AmpliSeq™ HID Y-SNP Research Panel v1 for performance and resolution in admixed South Americans of haplogroup Q.

Y haplogroups, defined by Y-SNPs, allow the reconstruction of the human Y chromosome genealogy, which is important for population, evolutionary and forensic genetics. In this study, Y-SNPs were typed and haplogroups inferred with the MPS Ion AmpliSeq™ HID Y-SNP Research Panel v1, as a high-throughput approach. Firstly, the performance of the panel was evaluated with different DNA input amounts, reagent volumes and cycle numbers. DNA-inputs from 0.5 to 1 ng generated the most balanced read depth. Combined with full reagent and 19 cycles, this offered the highest number of amplicons with a sequencing read depth of at least 20 reads. Secondly, the sub-haplogroups of 182 admixed South Americans and Greenlanders belonging to haplogroup Q were inferred and tested for potential improvement in resolution. Most samples were assigned to lineage Q-M3 with some samples assigned to lineages upstream (Q-M346, L56, L57; Q-L331, L53; Q-L54; Q-CTS11969, CTS11970) or parallel (Q-L330, L334; Q-Z780/M971) to Q-M3. Only one sample was assigned to a downstream lineage (Q-Z35615, Z35616). Most individuals of haplogroup Q with NAM ancestry could neither be distinguished from each other, nor from half of the Greenlandic samples. Typing additional, known SNPs within lineage Q-M3, Z19483 and SA05, increased the resolution of predicted haplogroups. The search for novel variants in the sequenced regions allowed the detection of 42 variants and the subdivision of lineage Q-M3 into new subclades. The variants found in six of these subclades were exclusive to certain South American countries. In light of the limited differentiation of haplogroup Q samples, the additional information on known or novel SNPs disclosed in this study when using MPS Ion AmpliSeq™ HID Y-SNP Research Panel v1 should be included in the Yleaf software, to increase the differentiation of lineage Q-M3.

The Ancestry of Eastern Paraguay: A Typical South American Profile with a Unique Pattern of Admixture.

Immigrants from diverse origins have arrived in Paraguay and produced important demographic changes in a territory initially inhabited by indigenous Guarani. Few studies have been performed to estimate the proportion of Native ancestry that is still preserved in Paraguay and the role of females and males in admixture processes. Therefore, 548 individuals from eastern Paraguay were genotyped for three marker sets: mtDNA, Y-SNPs and autosomal AIM-InDels. A genetic homogeneity was found between departments for each set of markers, supported by the demographic data collected, which showed that only 43% of the individuals have the same birthplace as their parents. The results show a sex-biased intermarriage, with higher maternal than paternal Native American ancestry. Within the native mtDNA lineages in Paraguay (87.2% of the total), most haplogroups have a broad distribution across the subcontinent, and only few are concentrated around the Paraná River basin. The frequency distribution of the European paternal lineages in Paraguay (92.2% of the total) showed a major contribution from the Iberian region. In addition to the remaining legacy of the colonial period, the joint analysis of the different types of markers included in this study revealed the impact of post-war migrations on the current genetic background of Paraguay.

Molecular and clinical insights into complex genomic rearrangements related to MECP2 duplication syndrome.

MECP2 duplication syndrome (MDS) is caused by copy number variation (CNV) spanning the MECP2 gene at Xq28 and is a major cause of intellectual disability (ID) in males. Herein, we describe two unrelated males harboring non-recurrent complex Xq28 rearrangements associated with MDS. Copy number gains were initially detected by quantitative real-time polymerase chain reaction and further delineated by high-resolution array comparative genomic hybridization, familial segregation, expression analysis and X-chromosome inactivation (XCI) evaluation in a carrier mother. SNVs within the rearrangements and/or fluorescent in situ hybridization (FISH) were used to assess the parental origin of the rearrangements. Patient 1 exhibited an intrachromosomal rearrangement, whose structure is consistent with a triplicated segment presumably embedded in an inverted orientation between two duplicated sequences (DUP-TRP/INV-DUP). The rearrangement was inherited from the carrier mother, who exhibits extreme XCI skewing and subtle psychiatric symptoms. Patient 2 presented a de novo (X;Y) unbalanced translocation resulting in duplication of Xq28 and deletion of Yp, originated in the paternal gametogenesis. Neurodevelopmental trajectory and non-neurological symptoms were consistent with previous reports, with the exception of cerebellar vermis hypoplasia in patient 2. Although both patients share the core MDS phenotype, patient 1 showed MECP2 transcript levels in blood similar to controls. Understanding the molecular mechanisms related to MDS is essential for designing targeted therapeutic strategies.

Investigating genetic diversity in admixed populations from Ecuador.

OBJECTIVES: According to demographic history, Ecuador has experienced shifts in its Native American populations caused by European colonization and the African slave trade. The continuous admixture events among Europeans, Native Americans, and Africans occurred differently in each region of the country, producing a stratified population. Thus, the aim of this study was to investigate the level of genetic substructure in the Ecuadorian Mestizo population. MATERIALS AND METHODS: A total of 377 male and 209 female samples were genotyped for two sets of X-chromosomal markers (32 X-Indels and 12 X-STRs). Population analyses performed included Hardy-Weinberg equilibrium tests, LD analysis, PCA, pairwise FST s, and AMOVA. RESULTS: Significant levels of LD were observed between markers separated by distances of less than 1 cM, as well as between markers separated by distances varying from 10.891 to 163.53 cM. Among Ecuadorian regions, Amazonia showed the highest average R2 value. DISCUSSION: When X-chromosomal and autosomal differentiation values were compared, a sex-biased admixture between European men and Native American and African women was revealed, as well as between African men and Native American women. Moreover, a distinct Native American ancestry was discernible in the Amazonian population, in addition to sex-biased gene flow between Amazonia and the Andes and Pacific coast regions. Overall, these results underline the importance of integrating X chromosome information to achieve a more comprehensive view of the genetic and demographic histories of South American admixed populations.

Patterns of genetic diversity in Colombia for 38 indels used in human identification.

The current population of Colombia has a genetic heterogeneity resulting from different migrations from other continents and within the country. In addition, there are small groups in their territory that have remained isolated and therefore have a different genetic pool in relation to that of the neighbouring urban populations. This population stratification must be considered in forensic analysis, being more complex for markers with marked intercontinental differentiation. In this study, population differentiation in Colombian admixed, native, and Afro-descendant populations was evaluated for a group of 38 indels described for forensic use. Allelic frequencies and parameters of forensic relevance were determined in each of the groups defined based on population differentiation analyses. In addition to the differences found between population groups, the results show that the set of 38 indels analysed could be useful in studies of individual identification in Colombia. The exclusion power presented by this set of markers suggests the need for joint use with other markers, being able to complement the STRs in paternity cases. High levels of both power of discrimination and exclusion were found when complementing the 38 HID-indels with a second multiplex, for a total of 83 indels.