Charting the Y-chromosome ancestry of present-day Argentinean Mennonites.

Old Order Mennonite communities initially arose in Northern Europe (centered in the Netherlands) and derived from the Anabaptist movement of the 16th century. Mennonites migrated to the New World in the early 18th century, first to North America, and more recently to Mesoamerica and South America. We analyzed Y-chromosome short tandem repeats (STRs) and single nucleotide polymorphisms in males from a community of Mennonites, 'La Nueva Esperanza', which arrived to Argentina in 1985 from colonies in Bolivia and Mexico. Molecular diversity indices coupled with demographic simulations show that Mennonites have a reduced variability when compared with local Argentinean populations and 69 European population samples. Mennonite Y-STR haplotypes were mainly observed in Central Europe. In agreement, multidimensional scaling analyses based on RST genetic distances indicate that Mennonite Y-chromosomes are closely related to Central/Northern Europeans (the Netherlands, Switzerland and Denmark). In addition, statistical inferences made on the most likely geographic origin of Y-chromosome haplotypes point more specifically to the Netherlands as the populations that best represent the majority of the Mennonite Y-chromosomes. Overall, Y-chromosome variation of Mennonites shows the signatures of moderate reduction of variability when compared with source populations, which is in good agreement with their lifestyle in small endogamous demes. These genetic singularities could also help to understand disease conditions that are more prevalent among Mennonites.

Male lineages in South American native groups: evidence of M19 traveling south.

With this study, we aimed to determine the different male ancestral components of two Native American communities from Argentina, namely Toba and Colla. The analysis of 27 Y-chromosome SNPs allowed us to identify seven different haplogroups in both samples. Chromosomes carrying the M3 mutation, which typically defines the Native American haplogroup Q1a3a, were seen most frequently in the Toba community (90%). Conversely, Q1a3a was represented in 34% of the Colla Y-chromosomes, whereas haplogroup R1b1, the main representative of western European populations, exhibited the highest frequency in this population (41%). Different M3 sublineages in the Toba community could be identified by observing point mutations at both DYS385 and M19 loci. A microvariant at DYS385, named 16.1, has been characterized, which helps to further subdivide Q1a3a. It is the first time the M19 mutated allele is described in a population from Argentina. This finding supports the old age of the lineages carrying the M19 mutation, but it contradicts the previous hypothesis that the M19 mutated allele is confined to only two Equatorial-Tucano population groups from the north region of South America. The detection of M19 further south than previously thought allows questioning of the hypothesis that this lineage serves as an example of isolation after colonization. This observation also affirms the strong genetic drift to which Native Americans have been subjected. Moreover, our study illustrates a heterogeneous contribution of Europeans to these populations and supports previous studies showing that most Native American groups were subjected to European admixture that primarily involved immigrant men.

Second GHEP-ISFG exercise for DVI: "DNA-led" victims' identification in a simulated air crash.

The Spanish and Portuguese-Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) has organized a second collaborative exercise on a simulated case of Disaster Victim Identification (DVI), with the participation of eighteen laboratories. The exercise focused on the analysis of a simulated plane crash case of medium-size resulting in 66 victims with varying degrees of fragmentation of the bodies (with commingled remains). As an additional difficulty, this second exercise included 21 related victims belonging to 6 families among the 66 missings to be identified. A total number of 228 post-mortem samples were represented with aSTR and mtDNA profiles, with a proportion of partial aSTR profiles simulating charred remains. To perform the exercise, participants were provided with aSTR and mtDNA data of 51 reference pedigrees -some of which deficient-including 128 donors for identification purposes. The exercise consisted firstly in the comparison of the post-mortem genetic profiles in order to re-associate fragmented remains to the same individual and secondly in the identification of the re-associated remains by comparing aSTR and mtDNA profiles with reference pedigrees using pre-established thresholds to report a positive identification. Regarding the results of the post-mortem samples re-associations, only a small number of discrepancies among participants were detected, all of which were from just a few labs. However, in the identification process by kinship analysis with family references, there were more discrepancies in comparison to the correct results. The identification results of single victims yielded fewer problems than the identification of multiple related victims within the same family groups. Several reasons for the discrepant results were detected: a) the identity/non-identity hypotheses were sometimes wrongly expressed in the likelihood ratio calculations, b) some laboratories failed to use all family references to report the DNA match, c) In families with several related victims, some laboratories firstly identified some victims and then unnecessarily used their genetic information to identify the remaining victims within the family, d) some laboratories did not correctly use "prior odds" values for the Bayesian treatment of the episode for both post-mortem/post-mortem re-associations as well as the ante-mortem/post-mortem comparisons to evaluate the probability of identity. For some of the above reasons, certain laboratories failed to identify some victims. This simulated "DNA-led" identification exercise may help forensic genetic laboratories to gain experience and expertize for DVI or MPI in using genetic data and comparing their own results with the ones in this collaborative exercise.

Paternal and maternal mutations in X-STRs: A GHEP-ISFG collaborative study.

The GHEP-ISFG organized a collaborative study to estimate mutation rates for the markers included in the Investigator Argus X-12 QS kit Qiagen. A total of 16 laboratories gathered data from 1,612 father/mother/daughter trios, which were used to estimate both maternal and paternal mutation rates, when pooled together with other already published data. Data on fathers and mothers' age at the time of birth of the daughter were also available for ∼93 % of the cases. Population analyses were computed considering the genetic information of a subset of 1,327 unrelated daughters, corresponding to 2,654 haplotypes from residents in several regions of five countries: Argentina, Brazil, Ecuador, Portugal and Spain. Genetic differentiation analyses between the population samples from the same country did not reveal signs of significant stratification, although results from Hardy-Weinberg and linkage disequilibrium tests indicated the need of larger studies for Ecuador and Brazilian populations. The high genetic diversity of the markers resulted in a large number of haplotype combinations, showing the need of huge databases for reliable estimates of their frequencies. It should also be noted the high number of new alleles found, many of them not included in the allelic ladders provided with the kit, as very diverse populations were analyzed. The overall estimates for locus specific average mutation rates varied between 7.5E-04 (for DXS7423) and 1.1E-02 (for DXS10135), the latter being a troublesome figure for kinship analyses. Most of the found mutations (∼92 %) are compatible with the gain or loss of a single repeat. Paternal mutation rates showed to be 5.2 times higher than maternal ones. We also found that older fathers were more prone to transmit mutated alleles, having this trend not been observed in the case of the mothers.

Genetic admixture patterns in Argentinian Patagonia.

As in other Latin American populations, Argentinians are the result of the admixture amongst different continental groups, mainly from America and Europe, and to a lesser extent from Sub-Saharan Africa. However, it is known that the admixture processes did not occur homogeneously throughout the country. Therefore, considering the importance for anthropological, medical and forensic researches, this study aimed to investigate the population genetic structure of the Argentinian Patagonia, through the analysis of 46 ancestry informative markers, in 433 individuals from five different localities. Overall, in the Patagonian sample, the average individual ancestry was estimated as 35.8% Native American (95% CI: 32.2-39.4%), 62.1% European (58.5-65.7%) and 2.1% African (1.7-2.4%). Comparing the five localities studied, statistically significant differences were observed for the Native American and European contributions, but not for the African ancestry. The admixture results combined with the genealogical information revealed intra-regional variations that are consistent with the different geographic origin of the participants and their ancestors. As expected, a high European ancestry was observed for donors with four grandparents born in Europe (96.8%) or in the Central region of Argentina (85%). In contrast, the Native American ancestry increased when the four grandparents were born in the North (71%) or in the South (61.9%) regions of the country, or even in Chile (60.5%). In summary, our results showed that differences on continental ancestry contribution have different origins in each region in Patagonia, and even in each locality, highlighting the importance of knowing the origin of the participants and their ancestors for the correct interpretation and contextualization of the genetic information.

Testing for genetic structure in different urban Argentinian populations.

Fifteen autosomal short tandem repeat (STR) markers (D3S1358, HUMTH01, D21S11, D18S51, PENTA E, D5S818, D13S317, D7S820, D16S539, CSF1PO, PENTA D, HUMvWA, D8S1179, HUMTPOX, FGA) were analyzed in 1734 individuals living in urban areas of cities from six different Argentinian provinces (Buenos Aires, Neuquén, Tucumán, La Pampa, San Luis, Santa Cruz) in order to determine if a common urban database could be used in Argentina for forensic purposes. Frequencies estimates, Hardy-Weinberg equilibrium (HWE), and other parameters of forensic interest were computed. Comparisons between the six populations, and with published data from one Native American population from Argentina and other urban populations from Argentina and Europe were also performed. Our results reveal evidences for population structure, both when testing for genetic differentiation and when comparing frequencies distributions between different pairs of populations. Therefore, caution should be taken when using a common pooled database with general forensic purposes in Argentina.

Species identification in forensic samples using the SPInDel approach: A GHEP-ISFG inter-laboratory collaborative exercise.

DNA is a powerful tool available for forensic investigations requiring identification of species. However, it is necessary to develop and validate methods able to produce results in degraded and or low quality DNA samples with the high standards obligatory in forensic research. Here, we describe a voluntary collaborative exercise to test the recently developed Species Identification by Insertions/Deletions (SPInDel) method. The SPInDel kit allows the identification of species by the generation of numeric profiles combining the lengths of six mitochondrial ribosomal RNA (rRNA) gene regions amplified in a single reaction followed by capillary electrophoresis. The exercise was organized during 2014 by a Working Commission of the Spanish and Portuguese-Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG), created in 2013. The 24 participating laboratories from 10 countries were asked to identify the species in 11 DNA samples from previous GHEP-ISFG proficiency tests using a SPInDel primer mix and control samples of the 10 target species. A computer software was also provided to the participants to assist the analyses of the results. All samples were correctly identified by 22 of the 24 laboratories, including samples with low amounts of DNA (hair shafts) and mixtures of saliva and blood. Correct species identifications were obtained in 238 of the 241 (98.8%) reported SPInDel profiles. Two laboratories were responsible for the three cases of misclassifications. The SPInDel was efficient in the identification of species in mixtures considering that only a single laboratory failed to detect a mixture in one sample. This result suggests that SPInDel is a valid method for mixture analyses without the need for DNA sequencing, with the advantage of identifying more than one species in a single reaction. The low frequency of wrong (5.0%) and missing (2.1%) alleles did not interfere with the correct species identification, which demonstrated the advantage of using a method based on the analysis of multiple loci. Overall, the SPInDel method was easily implemented by laboratories using different genotyping platforms, the interpretation of results was straightforward and the SPInDel software was used without any problems. The results of this collaborative exercise indicate that the SPInDel method can be applied successfully in forensic casework investigations.

A comprehensive Y-STR portrait of Argentinean populations.

A study of 23 Y-STRs was conducted in 257 individuals living in urban areas from eight Argentinean provinces. The data were meta-analyzed together with 364 profiles obtained from the literature that represent other five provinces. A total of 255 different haplotypes were observed (253 singletons). Genetic structure estimated from analysis of molecular variance (AMOVA) and exploring different grouping scenarios, yielded high within population variance. Not surprisingly, analyses of genetic distances with respect to main ancestral continental populations indicated Argentinean haplotypes to be closely related to European ones. Overall, these results provide a quite complete picture of the patterns of Y chromosome variation in Argentina, notably contributing to increase the previous national database, and consequently allowing a better estimation of parameters of interest in forensic casework and parentage testing.

Genomic continuity of Argentinean Mennonites.

Mennonites are Anabaptist communities that originated in Central Europe about 500 years ago. They initially migrated to different European countries, and in the early 18th century they established their first communities in North America, from where they moved to other American regions. We aimed to analyze an Argentinean Mennonite congregation from a genome-wide perspective by way of investigating >580.000 autosomal SNPs. Several analyses show that Argentinean Mennonites have European ancestry without signatures of admixture with other non-European American populations. Among the worldwide datasets used for population comparison, the CEU, which is the best-subrogated Central European population existing in The 1000 Genome Project, is the dataset showing the closest genome affinity to the Mennonites. When compared to other European population samples, the Mennonites show higher inbreeding coefficient values. Argentinean Mennonites show signatures of genetic continuity with no evidence of admixture with Americans of Native American or sub-Saharan African ancestry. Their genome indicates the existence of an increased endogamy compared to other Europeans most likely mirroring their lifestyle that involve small communities and historical consanguineous marriages.

STR data for PowerPlex 16 System from Neuquen population, SW Argentina.

Allele frequencies for the 15 autosomic STR loci included in the PowerPlex 16 System kit (Promega) were estimated from a sample of 111 unrelated individuals living in Neuquen province, southwest of Argentina. Population showed to be in HWE.

STR data for PowerPlex 16 System from Buenos Aires population, Argentina.

Allele frequencies for the 15 autosomic STR loci included in the PowerPlex 16 System kit (Promega Corp.) were estimated from a sample of 143 unrelated individuals living in Capital Federal and in Buenos Aires Metropolitan Area, Argentina. Population showed to be in HWE.

Evaluating methods to correct for population stratification when estimating paternity indexes.

The statistical interpretation of the forensic genetic evidence requires the use of allelic frequency estimates in the reference population for the studied markers. Differences in the genetic make up of the populations can be reflected in statistically different allelic frequency distributions. One can easily figure out that collecting such information for any given population is not always possible. Therefore, alternative approaches are needed in these cases in order to compensate for the lack of information. A number of statistics have been proposed to control for population stratification in paternity testing and forensic casework, Fst correction being the only one recommended by the forensic community. In this study we aimed to evaluate the performance of Fst to correct for population stratification in forensics. By way of simulations, we first tested the dependence of Fst on the relative sizes of the sub-populations, and second, we measured the effect of the Fst corrections on the Paternity Index (PI) values compared to the ones obtained when using the local reference database. The results provide clear-cut evidence that (i) Fst values are strongly dependent on the sampling scheme, and therefore, for most situations it would be almost impossible to estimate real values of Fst; and (ii) Fst corrections might unfairly correct PI values for stratification, suggesting the use of local databases whenever possible to estimate the frequencies of genetic profiles and PI values.

Y chromosome microsatellite genetic variation in two Native American populations from Argentina: population stratification and mutation data.

Two Native American populations from North and northwest regions of Argentina (Toba and Colla) were analyzed for 17 Y chromosome short tandem repeat loci (Y-STRs), namely, DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS385, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635 and GATA H4.1. Over 357 allele transfers, two one-step mutations could be detected at DYS456 and GATA H4.1 loci. A new 16.1 'micro-variant' allele was observed for DYS385, characterized by an insertion at the fifth GAAA repeat. We also observed two alleles at the DYS448 locus in three samples (two from Toba and one from Colla). A total of 34 and 16 different haplotypes were detected for Toba and Colla, respectively, the former with a haplotype diversity value of 0.9769+/-0.01, whereas 0.9497+/-0.02 for the latter. Significant population differences were observed between Colla and Toba, at least in part, due to a more prevalent European input in the Colla. In agreement with this observation is the fact that the genetic distances between Colla and Iberian populations are lower than those observed between Iberian and any other Native American population. The results of multiscaling dimensional analysis and genetic distances (Rst) among Native American population samples also reflect this fact. The data show the existence of clear population stratification in the Argentina, a fact that should be taken into account in forensic casework.

Variabilidad genética en 17 microsatélites del Cromosoma Y en dos poblaciones nativas de Argentina (Tobas, Collas) / Genetic Variability in 17 Y-STR loci of a native population from Argentine (Tobas, Collas)

Se analizan 17 marcadores microsatélites del cromosoma Y del tipo Short Tandem Repeats (STRs) (DYS19,DYS3891, DYS38911, DYS390, DYS391, DYS392, DYS393, DYS385, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, y GATA H4)en dos poblaciones nativas de las regiones Norte y Noroeste de la Argentina (Tobas y Collas)