Continental Origin for Q Haplogroup Patrilineages in Argentina and Paraguay.

Haplogroup Q originated in Eurasia around 30,000 years ago. It is present in Y-chromosomes from Asia and Europe at rather low frequencies. Since America is undoubtedly one of the continents where this haplogroup is highly represented, it has been defined as one of the founding haplogroups. Its M3 clade has been early described as the most frequent, with pan-American representation. However, it was also possible to find several other haplogroup Q clades at low frequencies. Numerous mutations have been described for haplogroup Q, allowing analysis of its variability and assignment of its geographic origin. We have analyzed 442 samples of unrelated men from Argentina and Paraguay belonging to haplogroup Q; here we report specifically on 27 Q (xM3) lineages. We tested 3 single-nucleotide polymorphisms (SNPs) by amplified product-length polymorphism (APLP) analysis, 3 SNPs for restriction fragment length polymorphism (RFLP) analysis, 15 SNPs by Sanger sequencing, and 17 short tandem repeats (STRs). Our approach allowed us to identify five subhaplogroups. Q-M3 and Q-CTS2730/Z780 are undoubtedly autochthonous lineages and represent the most frequent subhaplogroups, with significant representation in self-defined aboriginal populations, and their autochthonous status has been previously described. The aim of present work was to identify the continental origin of the remaining Q lineages. Thus, we analyzed the STR haplotypes for the samples and compared them with haplotypes described by other authors for the rest of the world. Even when haplogroup Q lineages have been extensively studied in America, some of them could have their origin in post-Columbian human migration from Europe and Middle East.

GBStools: A Statistical Method for Estimating Allelic Dropout in Reduced Representation Sequencing Data.

Reduced representation sequencing methods such as genotyping-by-sequencing (GBS) enable low-cost measurement of genetic variation without the need for a reference genome assembly. These methods are widely used in genetic mapping and population genetics studies, especially with non-model organisms. Variant calling error rates, however, are higher in GBS than in standard sequencing, in particular due to restriction site polymorphisms, and few computational tools exist that specifically model and correct these errors. We developed a statistical method to remove errors caused by restriction site polymorphisms, implemented in the software package GBStools. We evaluated it in several simulated data sets, varying in number of samples, mean coverage and population mutation rate, and in two empirical human data sets (N = 8 and N = 63 samples). In our simulations, GBStools improved genotype accuracy more than commonly used filters such as Hardy-Weinberg equilibrium p-values. GBStools is most effective at removing genotype errors in data sets over 100 samples when coverage is 40X or higher, and the improvement is most pronounced in species with high genomic diversity. We also demonstrate the utility of GBS and GBStools for human population genetic inference in Argentine populations and reveal widely varying individual ancestry proportions and an excess of singletons, consistent with recent population growth.

Reconstructing Native American migrations from whole-genome and whole-exome data.

There is great scientific and popular interest in understanding the genetic history of populations in the Americas. We wish to understand when different regions of the continent were inhabited, where settlers came from, and how current inhabitants relate genetically to earlier populations. Recent studies unraveled parts of the genetic history of the continent using genotyping arrays and uniparental markers. The 1000 Genomes Project provides a unique opportunity for improving our understanding of population genetic history by providing over a hundred sequenced low coverage genomes and exomes from Colombian (CLM), Mexican-American (MXL), and Puerto Rican (PUR) populations. Here, we explore the genomic contributions of African, European, and especially Native American ancestry to these populations. Estimated Native American ancestry is 48% in MXL, 25% in CLM, and 13% in PUR. Native American ancestry in PUR is most closely related to populations surrounding the Orinoco River basin, confirming the Southern American ancestry of the Taíno people of the Caribbean. We present new methods to estimate the allele frequencies in the Native American fraction of the populations, and model their distribution using a demographic model for three ancestral Native American populations. These ancestral populations likely split in close succession: the most likely scenario, based on a peopling of the Americas 16 thousand years ago (kya), supports that the MXL Ancestors split 12.2kya, with a subsequent split of the ancestors to CLM and PUR 11.7kya. The model also features effective populations of 62,000 in Mexico, 8,700 in Colombia, and 1,900 in Puerto Rico. Modeling Identity-by-descent (IBD) and ancestry tract length, we show that post-contact populations also differ markedly in their effective sizes and migration patterns, with Puerto Rico showing the smallest effective size and the earlier migration from Europe. Finally, we compare IBD and ancestry assignments to find evidence for relatedness among European founders to the three populations.

[Association between NAT2 polymorphisms with non-syndromic cleft lip with or without cleft palate in Argentina]. / Asociación entre polimorfismos del gen NAT2 y fisura labiopalatina no sindrómica en Argentina.

BACKGROUND: NAT genes are considered candidate genes for the genetic predisposition to non-syndromic Cleft lip with or without cleft palate (NSCLP), since they codify for N-acetyltransferases, enzymes responsible for the biotransformation of arylamines, hydrazine drugs, and a great number of toxins and carcinogens present in diet, cigarette smoke, and environment. AIM: To determine the association between alleles determining slow acetylator phenotype and the risk of NSCLP. MATERIAL AND METHODS: We analyzed *5 (481C>T), *6 (590G>A) and *7 (857G>A) alleles which determine the slow acetylator phenotype and *4 (wild type) allele by polymerase chain reaction/restriction fragment length polymorphism in 97 progenitor-case trios of NSCLP in Argentinian Obstetric Wards. We evaluated the transmission disequilibrium (TDT). RESULTS: TDT showed a positive association between allele *5 and NSCLP (odds ratio = 1,6; p = 0,03). CONCLUSIONS: The presence of *5 allele is significantly higher in cases with congenital NSCLP.

Human Y-chromosome SNP characterization by multiplex amplified product-length polymorphism analysis.

We designed an allele-specific amplification protocol to optimize Y-chromosome SNP typing, which is an unavoidable step for defining the phylogenetic status of paternal lineages. It allows the simultaneous highly specific definition of up to six mutations in a single reaction by amplification fragment length polymorphism (AFLP) without the need of specialized equipment, at a considerably lower cost than that based on single-base primer extension (SNaPshot™) technology or PCR-RFLP systems, requiring as little as 0.5 ng DNA and compatible with the small fragments characteristic of low-quality DNA. By designation of two primers recognizing the derived and ancestral state for each SNP, which can be differentiated by size by the addition of a noncomplementary nucleotide tail, we could define major Y clades E, F, K, R, Q, and subhaplogroups R1, R1a, R1b, R1b1b, R1b1c, J1, J2, G1, G2, I1, Q1a3, and Q1a3a1 through amplification fragments that ranged between 60 and 158bp.

Apolipoprotein E Polymorphisms in Andean Population of Jujuy, Argentina.

The pleiotropic nature of the apolipoprotein E (APOE) gene is associated with complex diseases in different populations. We analyzed APOE polymorphisms in 76 individuals from Jujuy - Argentina using NGS technology. The observed genotypes align with the expected Hardy-Weinberg equilibrium. APOE3 was the most common allele, followed by APOE4 and APOE2. The allele distribution pattern is consistent with findings in previously studied populations of Native Americans and Asians. The E4 allele's low frequency, always observed in a heterozygous state, raises questions regarding its relevance in explaining dementia and longevity associated with this marker in the Central Andes.

Human Y chromosome sequences from Q Haplogroup reveal a South American settlement pre-18,000 years ago and a profound genomic impact during the Younger Dryas.

The settlement of the Americas has been the focus of incessant debate for more than 100 years, and open questions regarding the timing and spatial patterns of colonization still remain today. Phylogenetic studies with complete human Y chromosome sequences are used as a highly informative tool to investigate the history of human populations in a given time frame. To study the phylogenetic relationships of Native American lineages and infer the settlement history of the Americas, we analyzed Y chromosome Q Haplogroup, which is a Pan-American haplogroup and represents practically all Native American lineages in Mesoamerica and South America. We built a phylogenetic tree for Q Haplogroup based on 102 whole Y chromosome sequences, of which 13 new Argentine sequences were provided by our group. Moreover, 1,072 new single nucleotide polymorphisms (SNPs) that contribute to its resolution and diversity were identified. Q-M848 is known to be the most frequent autochthonous sub-haplogroup of the Americas. The present is the first genomic study of Q Haplogroup in which current knowledge on Q-M848 sub-lineages is contrasted with the historical, archaeological and linguistic data available. The divergence times, spatial structure and the SNPs found here as novel for Q-Z780, a less frequent sub-haplogroup autochthonous of the Americas, provide genetic support for a South American settlement before 18,000 years ago. We analyzed how environmental events that occurred during the Younger Dryas period may have affected Native American lineages, and found that this event may have caused a substantial loss of lineages. This could explain the current low frequency of Q-Z780 (also perhaps of Q-F4674, a third possible sub-haplogroup autochthonous of the Americas). These environmental events could have acted as a driving force for expansion and diversification of the Q-M848 sub-lineages, which show a spatial structure that developed during the Younger Dryas period.

Software for Y-haplogroup predictions: a word of caution.

The development of online software designed for genetic studies has been exponentially growing, providing numerous benefits to the scientific community. However, they should be used with care, since some require adjustments. The efficiency of two programs for haplogroup prediction was tested with 119 samples of known haplotypes and haplogroups from Argentine populations. Quantitative estimates of the predictive quality of both software systems were computed with the uncertainty coefficient; and sensitivity, specificity, positive, and negative likelihood ratios were also calculated to assert the reliability of both programs, showing high probabilities of assigning an incorrect haplogroup.

Fine-scale genomic analyses of admixed individuals reveal unrecognized genetic ancestry components in Argentina.

Similarly to other populations across the Americas, Argentinean populations trace back their genetic ancestry into African, European and Native American ancestors, reflecting a complex demographic history with multiple migration and admixture events in pre- and post-colonial times. However, little is known about the sub-continental origins of these three main ancestries. We present new high-throughput genotyping data for 87 admixed individuals across Argentina. This data was combined to previously published data for admixed individuals in the region and then compared to different reference panels specifically built to perform population structure analyses at a sub-continental level. Concerning the Native American ancestry, we could identify four Native American components segregating in modern Argentinean populations. Three of them are also found in modern South American populations and are specifically represented in Central Andes, Central Chile/Patagonia, and Subtropical and Tropical Forests geographic areas. The fourth component might be specific to the Central Western region of Argentina, and it is not well represented in any genomic data from the literature. As for the European and African ancestries, we confirmed previous results about origins from Southern Europe, Western and Central Western Africa, and we provide evidences for the presence of Northern European and Eastern African ancestries.

Brief communication: Restricted geographic distribution for Y-Q* paragroup in South America.

We analyzed 21 paragroup Q* Y chromosomes from South American aboriginal and urban populations. Our aims were to evaluate the phylogenetic status, geographic distribution, and genetic diversity in these groups of chromosomes and compare the degree of genetic variation in relation to Q1a3a haplotypes. All Q* chromosomes from our series and five samples from North American Q* presented the derivate state for M346, that is present upstream to M3, and determined Q1a3* paragroup. We found a restrictive geographic distribution and low frequency of Q1a3* in South America. We assumed that this low frequency could be reflecting extreme drift effects. However, several estimates of gene diversity do not support the existence of a severe bottleneck. The mean haplotype diversity expected was similar to that for South American Q1a3* and Q1a3a (0.478 and 0.501, respectively). The analysis of previous reports from other research groups and this study shows the highest frequencies of Q* for the West Corner and the Grand Chaco regions of South America. At present, there is no information on whether the phylogenetic status of Q* paragoup described in previous reports is similar to that of Q1a3* paragroup though our results support this possibility.

Population structure in Argentina.

We analyzed 391 samples from 12 Argentinian populations from the Center-West, East and North-West regions with the Illumina Human Exome Beadchip v1.0 (HumanExome-12v1-A). We did Principal Components analysis to infer patterns of populational divergence and migrations. We identified proportions and patterns of European, African and Native American ancestry and found a correlation between distance to Buenos Aires and proportion of Native American ancestry, where the highest proportion corresponds to the Northernmost populations, which is also the furthest from the Argentinian capital. Most of the European sources are from a South European origin, matching historical records, and we see two different Native American components, one that spreads all over Argentina and another specifically Andean. The highest percentages of African ancestry were in the Center West of Argentina, where the old trade routes took the slaves from Buenos Aires to Chile and Peru. Subcontinentaly, sources of this African component are represented by both West Africa and groups influenced by the Bantu expansion, the second slightly higher than the first, unlike North America and the Caribbean, where the main source is West Africa. This is reasonable, considering that a large proportion of the ships arriving at the Southern Hemisphere came from Mozambique, Loango and Angola.

Asociación entre polimorfismos del gen NAT2 y fisura labiopalatina no sindrómica en Argentina / Association between NAT2 polymorphisms with non-syndromic cleft lip with or without cleft palate in Argentina

Background: NAT genes are considered candidate genes for the genetic predisposition to non-syndromic Cleft lip with or without cleft palate (NSCLP), since they codify for N-acetyltransferases, enzymes responsible for the biotransformation of arylamines, hydrazine drugs, and a great number of toxins and carcinogens present in diet, cigarette smoke, and environment. Aim: To determine the association between alleles determining slow acetylator phenotype and the risk of NSCLP. Material and Methods: We analyzed *5 (481C>T), *6 (590G>A) and *7 (857G>A) alleles which determine the slow acetylator phenotype and *4 (wild type) allele by polymerase chain reaction/restriction fragment length polymorphism in 97 progenitor-case trios of NSCLP in Argentinian Obstetric Wards. We evaluated the transmission disequilibrium (TDT). Results: TDT showed a positive association between allele *5 and NSCLP (odds ratio = 1,6; p = 0,03). Conclusions: The presence of *5 allele is significantly higher in cases with congenital NSCLP.