Carlos Vullo

A GEP-ISFG collaborative study on the optimization of an X-STR decaplex: data on 15 Iberian and Latin American populations

In a collaborative work carried out by the Spanish and Portuguese ISFG Working Group (GEP-ISFG), a polymerase chain reaction multiplex was optimized in order to type ten X-chromosome short tandem repeats (STRs) in a single reaction, including: DXS8378, DXS9902, DXS7132, DXS9898, DXS6809, DXS6789, DXS7133, GATA172D05, GATA31E08, and DXS7423. Using this X-decaplex, each 17 of the participating laboratories typed a population sample of approximately 200 unrelated individuals (100 males and 100 females). In this work, we report the allele frequencies for the ten X-STRs in 15 samples from Argentina (Buenos Aires, Córdoba, Río Negro, Entre Ríos, and Misiones), Brazil (São Paulo, Rio de Janeiro, Paraná, and Mato Grosso do Sul), Colombia (Antioquia), Costa Rica, Portugal (Northern and Central regions), and Spain (Galicia and Cantabria). Gene diversities were calculated for the ten markers in each population and all values were above 56%. The average diversity per locus varied between 66%, for DXS7133, and 82%, for DXS6809. For this set of STRs, a high discrimination power was obtained in all populations, both in males (≥1 in 5 × 105) and females (≥1 in 3 × 109), as well as high mean exclusion chance in father/daughter duos (≥99.953%) and in father/mother/daughter trios (≥99.999%). Genetic distance analysis showed no significant differences between northern and central Portugal or between the two Spanish samples from Galicia and Cantabria. Inside Brazil, significant differences were found between Rio de Janeiro and the other three populations, as well as between São Paulo and Paraná. For the five Argentinean samples, significant distances were only observed when comparing Misiones with Entre Ríos and with Río Negro, the only two samples that do not differ significantly from Costa Rica. Antioquia differed from all other samples, except the one from Río Negro.

The Genetic Legacy of the Pre-Colonial Period in Contemporary Bolivians

Only a few genetic studies have been carried out to date in Bolivia. However, some of the most important (pre)historical enclaves of South America were located in these territories. Thus, the (sub)-Andean region of Bolivia was part of the Inca Empire, the largest state in Pre-Columbian America. We have genotyped the first hypervariable region (HVS-I) of 720 samples representing the main regions in Bolivia, and these data have been analyzed in the context of other pan-American samples (>19,000 HVS-I mtDNAs). Entire mtDNA genome sequencing was also undertaken on selected Native American lineages. Additionally, a panel of 46 Ancestry Informative Markers (AIMs) was genotyped in a sub-set of samples. The vast majority of the Bolivian mtDNAs (98.4%) were found to belong to the main Native American haplogroups (A: 14.3%, B: 52.6%, C: 21.9%, D: 9.6%), with little indication of sub-Saharan and/or European lineages; however, marked patterns of haplogroup frequencies between main regions exist (e.g. haplogroup B: Andean [71%], Sub-Andean [61%], Llanos [32%]). Analysis of entire genomes unraveled the phylogenetic characteristics of three Native haplogroups: the pan-American haplogroup B2b (originated ∼21.4 thousand years ago [kya]), A2ah (∼5.2 kya), and B2o (∼2.6 kya). The data suggest that B2b could have arisen in North California (an origin even in the north most region of the American continent cannot be disregarded), moved southward following the Pacific coastline and crossed Meso-America. Then, it most likely spread into South America following two routes: the Pacific path towards Peru and Bolivia (arriving here at about ∼15.2 kya), and the Amazonian route of Venezuela and Brazil southwards. In contrast to the mtDNA, Ancestry Informative Markers (AIMs) reveal a higher (although geographically variable) European introgression in Bolivians (25%). Bolivia shows a decreasing autosomal molecular diversity pattern along the longitudinal axis, from the Altiplano to the lowlands. Both autosomes and mtDNA revealed a low impact (1–2%) of a sub-Saharan component in Bolivians.

Frequency data for 12 mini STR loci in Argentina

Allele frequencies and forensic parameters for twelve miniSTR autosomal loci (D10S1248, D14S1434, D22S1045, D4S2364, D2S441, D1S1677, D20S480, D6S2439, D6S1056, D9S1118, D4S2639 and D17S1290) were calculated from a sample of 506 unrelated individuals from the Central-East Region of Argentina. No significant deviations from Hardy-Weinberg expectations were found. Furthermore, comparisons with other previously studied populations were made. These twelve miniSTR markers may help forensic laboratories in solving parentage testing as well as in typing degraded DNA samples.

The complete mitogenome of a 500-year-old Inca child mummy

In 1985, a frozen mummy was found in Cerro Aconcagua (Argentina). Archaeological studies identified the mummy as a seven-year-old Inca sacrifice victim who lived >500 years ago, at the time of the expansion of the Inca Empire towards the southern cone. The sequence of its entire mitogenome was obtained. After querying a large worldwide database of mitogenomes (>28,000) we found that the Inca haplotype belonged to a branch of haplogroup C1b (C1bi) that has not yet been identified in modern Native Americans. The expansion of C1b into the Americas, as estimated using 203 C1b mitogenomes, dates to the initial Paleoindian settlements (~18.3 thousand years ago [kya]); however, its internal variation differs between Mesoamerica and South America. By querying large databases of control region haplotypes (>150,000), we found only a few C1bi members in Peru and Bolivia (e.g. Aymaras), including one haplotype retrieved from ancient DNA of an individual belonging to the Wari Empire (Peruvian Andes). Overall, the results suggest that the profile of the mummy represents a very rare sub-clade that arose 14.3 (5–23.6) kya and could have been more frequent in the past. A Peruvian Inca origin for present-day C1bi haplotypes would satisfy both the genetic and paleo-anthropological findings.

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.

Association between Y haplogroups and autosomal AIMs reveals intra-population substructure in Bolivian populations

For the correct evaluation of the weight of genetic evidence in a forensic context, databases must reflect the structure of the population, with all possible groups being represented. Countries with a recent history of admixture between strongly differentiated populations are usually highly heterogeneous and sub-structured. Bolivia is one of these countries, with a high diversity of ethnic groups and different levels of admixture (among Native Americans, Europeans and Africans) across the territory. For a better characterization of the male lineages in Bolivia, 17 Y-STR and 42 Y-SNP loci were genotyped in samples from La Paz and Chuquisaca. Only European and Native American Y-haplogroups were detected, and no sub-Saharan African chromosomes were found. Significant differences were observed between the two samples, with a higher frequency of European lineages in Chuquisaca than in La Paz. A sample belonging to haplogroup Q1a3a1a1-M19 was detected in La Paz, in a haplotype background different from those previously found in Argentina. This result supports an old M19 North–south dispersion in South America, possibly via two routes. When comparing the ancestry of each individual assessed through his Y chromosome with the one estimated using autosomal AIMs, (a) increased European ancestry in individuals with European Y chromosomes and (b) higher Native American ancestry in the carriers of Native American Y-haplogroups were observed, revealing an association between autosomal and Y-chromosomal markers. The results of this study demonstrate that a sub-structure does exist in Bolivia at both inter- and intrapopulation levels, a fact which must be taken into account in the evaluation of forensic genetic evidence.

Ancestry informative markers: inference of ancestry in aged bone samples using an autosomal AIM-Indel multiplex.

Ancestry informative markers (AIMs) can be useful to infer ancestry proportions of the donors of forensic evidence. The probability of success typing degraded samples, such as human skeletal remains, is strongly influenced by the DNA fragment lengths that can be amplified and the presence of PCR inhibitors. Several AIM panels are available amongst the many forensic marker sets developed for genotyping degraded DNA. Using a 46 AIM Insertion Deletion (Indel) multiplex, we analyzed human skeletal remains of post mortem time ranging from 35 to 60 years from four different continents (Sub-Saharan Africa, South and Central America, East Asia and Europe) to ascertain the genetic ancestry components. Samples belonging to non-admixed individuals could be assigned to their corresponding continental group. For the remaining samples with admixed ancestry, it was possible to estimate the proportion of co-ancestry components from the four reference population groups. The 46 AIM Indel set was informative enough to efficiently estimate the proportion of ancestry even in samples yielding partial profiles, a frequent occurrence when analyzing inhibited and/or degraded DNA extracts.

Phylogeographic and genome-wide investigations of Vietnam ethnic groups reveal signatures of complex historical demographic movements

The territory of present-day Vietnam was the cradle of one of the world’s earliest civilizations, and one of the first world regions to develop agriculture. We analyzed the mitochondrial DNA (mtDNA) complete control region of six ethnic groups and the mitogenomes from Vietnamese in The 1000 Genomes Project (1000G). Genome-wide data from 1000G (~55k SNPs) were also investigated to explore different demographic scenarios. All Vietnamese carry South East Asian (SEA) haplotypes, which show a moderate geographic and ethnic stratification, with the Mong constituting the most distinctive group. Two new mtDNA clades (M7b1a1f1 and F1f1) point to historical gene flow between the Vietnamese and other neighboring countries. Bayesian-based inferences indicate a time-deep and continuous population growth of Vietnamese, although with some exceptions. The dramatic population decrease experienced by the Cham 700 years ago (ya) fits well with the Nam tiến (“southern expansion”) southwards from their original heartland in the Red River Delta. Autosomal SNPs consistently point to important historical gene flow within mainland SEA, and add support to a main admixture event occurring between Chinese and a southern Asian ancestral composite (mainly represented by the Malay). This admixture event occurred ~800 ya, again coinciding with the Nam tiến.