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.

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.

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.

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.

Allele frequency data for 23 aSTR for different ethnic groups from Republic of Zimbabwe.

In order to determine the population allele frequencies of autosomal STR markers of forensic interest in the Zimbabwean population, we analyzed a sample of 478 individuals from 19 different ethnic groups using the PowerPlex® Fusion 6C Kit (Promega Corp, Madison, Wisconsin). The data obtained were compared among the different Zimbabwean ethnic groups as well as with several African populations to establish whether significant differences exist among them. No significant differences were found among the ethnic groups in Zimbabwe. Statistically significant differences were observed between allele frequencies in Zimbabwe and some other African populations, although FST with neighboring Bantu populations from South and Southeast regions were low (below 0.005 in most single locus comparisons).

Searching for the roots of the first free African American community.

San Basilio de Palenque is an Afro-descendant community near Cartagena, Colombia, founded in the sixteenth century. The recognition of the historical and cultural importance of Palenque has promoted several studies, namely concerning the African roots of its first inhabitants. To deepen the knowledge of the origin and diversity of the Palenque parental lineages, we analysed a sample of 81 individuals for the entire mtDNA Control Region as well as 92 individuals for 27 Y-STRs and 95 for 51 Y-SNPs. The results confirmed the strong isolation of the Palenque, with some degree of influx of Native American maternal lineages, and a European admixture exclusively mediated by men. Due to the high genetic drift observed, a pairwise FST analysis with available data on African populations proved to be inadequate for determining population affinities. In contrast, when a phylogenetic approach was used, it was possible to infer the phylogeographic origin of some lineages in Palenque. Contradicting previous studies indicating a single African origin, our results evidence parental genetic contributions from widely different African regions.

Twenty Years Later: A Comprehensive Review of the X Chromosome Use in Forensic Genetics.

The unique structure of the X chromosome shaped by evolution has led to the present gender-specific genetic differences, which are not shared by its counterpart, the Y chromosome, and neither by the autosomes. In males, recombination between the X and Y chromosomes is limited to the pseudoautosomal regions, PAR1 and PAR2; therefore, in males, the X chromosome is (almost) entirely transmitted to female offspring. On the other hand, the X chromosome is present in females with two copies that recombine along the whole chromosome during female meiosis and that is transmitted to both female and male descendants. These transmission characteristics, besides the obvious clinical impact (sex chromosome aneuploidies are extremely frequent), make the X chromosome an irreplaceable genetic tool for population genetic-based studies as well as for kinship and forensic investigations. In the early 2000s, the number of publications using X-chromosomal polymorphisms in forensic and population genetic applications increased steadily. However, nearly 20 years later, we observe a conspicuous decrease in the rate of these publications. In light of this observation, the main aim of this article is to provide a comprehensive review of the advances and applications of X-chromosomal markers in population and forensic genetics over the last two decades. The foremost relevant topics are addressed as: (i) developments concerning the number and types of markers available, with special emphasis on short tandem repeat (STR) polymorphisms (STR nomenclatures and practical concerns); (ii) overview of worldwide population (frequency) data; (iii) the use of X-chromosomal markers in (complex) kinship testing and the forensic statistical evaluation of evidence; (iv) segregation and mutation studies; and (v) current weaknesses and future prospects.

Evaluation of the Precision of Ancestry Inferences in South American Admixed Populations.

Ancestry informative markers (AIMs) are used in forensic genetics to infer biogeographical ancestry (BGA) of individuals and may also have a prominent role in future police and identification investigations. In the last few years, many studies have been published reporting new AIM sets. These sets include markers (usually around 100 or less) selected with different purposes and different population resolutions. Regardless of the ability of these sets to separate populations from different continents or regions, the uncertainty associated with the estimates provided by these panels and their capacity to accurately report the different ancestral contributions in individuals of admixed populations has rarely been investigated. This issue is addressed in this study by evaluating different AIM sets. Ancestry inference was carried out in admixed South American populations, both at population and individual levels. The results of ancestry inferences using AIM sets with different numbers of markers among admixed reference populations were compared. To evaluate the performance of the different ancestry panels at the individual level, expected and observed estimates among families and their offspring were compared, considering that (1) the apportionment of ancestry in the offspring should be closer to the average ancestry of the parents, and (2) full siblings should present similar ancestry values. The results obtained illustrate the importance of having a good balance/compromise between not only the number of markers and their ability to differentiate ancestral populations, but also a balanced differentiation among reference groups, to obtain more precise values of genetic ancestry. This work also highlights the importance of estimating errors associated with the use of a limited number of markers. We demonstrate that although these errors have a moderate effect at the population level, they may have an important impact at the individual level. Considering that many AIM-sets are being described for inferences at the individual level and not at the population level, e.g., in association studies or the determination of a suspect's BGA, the results of this work point to the need of a more careful evaluation of the uncertainty associated with the ancestry estimates in admixed populations, when small AIM-sets are used.

The first GHEP-ISFG collaborative exercise on forensic applications of massively parallel sequencing.

One of the main goals of the Spanish and Portuguese-Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) is to promote and contribute to the development and dissemination of scientific knowledge in the field of forensic genetics. The GHEP-ISFG supports several Working Commissions which develop different scientific activities. One of them, the Working Commission on "Massively Parallel Sequencing (MPS): Forensic Applications", organized its first collaborative exercise on forensic applications of MPS technology in 2019. The aim of this exercise was to assess the concordance between the MPS results and those obtained with conventional technologies (capillary electrophoresis and Sanger sequencing), as well as to compare the results obtained within the different MPS platforms and/or the different kits/panels and analysis software packages (commercial and open-access) available on the market. The seven participating laboratories analyzed some samples of the annual GHEP-ISFG proficiency test (EIADN No. 27 (2019)), using Ion Torrent™ or MiSeq FGx® platforms. Six of them sent autosomal STR sequence data, five laboratories performed MPS analysis of individual identification SNPs, four laboratories reported MPS data of Y-chromosomal STRs, and X-chromosomal STRs, three laboratories performed MPS analysis of ancestry informative SNPs and phenotype informative SNPs, two labs performed MPS analysis of the mitochondrial DNA control region, and only one lab produced MPS data of lineage informative SNPs. Autosomal STR sequencing results were highly concordant to the consensus obtained by capillary electrophoresis in the EIADN No. 27 (2019) exercise. Furthermore, in general, a high level of concordance was observed between the results of the participating laboratories, regardless of the platform used. The main discordances were due to errors during the analysis process or from sequence data obtained with low depth of coverage. In this paper we highlight some issues that still arise, such as standardization of the nomenclature for STRs analyzed by sequencing with MPS, the universal uptake of a nomenclature framework by the analysis software, and well established validation and accreditation of the new MPS platforms for use in routine forensic case-work.