Polymorphic transposable elements provide new insights on high-altitude adaptation in the Tibetan Plateau.

Several studies demonstrated that populations living in the Tibetan plateau are genetically and physiologically adapted to high-altitude conditions, showing genomic signatures ascribable to the action of natural selection. However, so far most of them relied solely on inferences drawn from the analysis of coding variants and point mutations. To fill this gap, we focused on the possible role of polymorphic transposable elements in influencing the adaptation of Tibetan and Sherpa highlanders. To do so, we compared high-altitude and middle/low-lander individuals of East Asian ancestry by performing in silico analyses and differentiation tests on 118 modern and ancient samples. We detected several transposable elements associated with high altitude, which map genes involved in cardiovascular, hematological, chem-dependent and respiratory conditions, suggesting that metabolic and signaling pathways taking part in these functions are disproportionately impacted by the effect of environmental stressors in high-altitude individuals. To our knowledge, our study is the first hinting to a possible role of transposable elements in the adaptation of Tibetan and Sherpa highlanders.

Dietary, Cultural, and Pathogens-Related Selective Pressures Shaped Differential Adaptive Evolution among Native Mexican Populations.

Native American genetic ancestry has been remarkably implicated with increased risk of diverse health issues in several Mexican populations, especially in relation to the dramatic changes in environmental, dietary, and cultural settings they have recently undergone. In particular, the effects of these ecological transitions and Westernization of lifestyles have been investigated so far predominantly on Mestizo individuals. Nevertheless, indigenous groups, rather than admixed Mexicans, have plausibly retained the highest proportions of genetic components shaped by natural selection in response to the ancient milieu experienced by Mexican ancestors during their pre-Columbian evolutionary history. These formerly adaptive variants have the potential to represent the genetic determinants of some biological traits that are peculiar to Mexican people, as well as a reservoir of loci with possible biomedical relevance. To test such a hypothesis, we used genome-wide genotype data to infer the unique adaptive evolution of Native Mexican groups selected as reasonable descendants of the main pre-Columbian Mexican civilizations. A combination of haplotype-based and gene-network analyses enabled us to detect genomic signatures ascribable to polygenic adaptive traits plausibly evolved by the main genetic clusters of Mexican indigenous populations to cope with local environmental and/or cultural conditions. Some of these adaptations were found to play a role in modulating the susceptibility/resistance of these groups to certain pathological conditions, thus providing new evidence that diverse selective pressures have contributed to shape the current biological and disease-risk patterns of present-day Native and Mestizo Mexican populations.

Genomic history of the Italian population recapitulates key evolutionary dynamics of both Continental and Southern Europeans.

BACKGROUND: The cline of human genetic diversity observable across Europe is recapitulated at a micro-geographic scale by variation within the Italian population. Besides resulting from extensive gene flow, this might be ascribable also to local adaptations to diverse ecological contexts evolved by people who anciently spread along the Italian Peninsula. Dissecting the evolutionary history of the ancestors of present-day Italians may thus improve the understanding of demographic and biological processes that contributed to shape the gene pool of European populations. However, previous SNP array-based studies failed to investigate the full spectrum of Italian variation, generally neglecting low-frequency genetic variants and examining a limited set of small effect size alleles, which may represent important determinants of population structure and complex adaptive traits. To overcome these issues, we analyzed 38 high-coverage whole-genome sequences representative of population clusters at the opposite ends of the cline of Italian variation, along with a large panel of modern and ancient Euro-Mediterranean genomes. RESULTS: We provided evidence for the early divergence of Italian groups dating back to the Late Glacial and for Neolithic and distinct Bronze Age migrations having further differentiated their gene pools. We inferred adaptive evolution at insulin-related loci in people from Italian regions with a temperate climate, while possible adaptations to pathogens and ultraviolet radiation were observed in Mediterranean Italians. Some of these adaptive events may also have secondarily modulated population disease or longevity predisposition. CONCLUSIONS: We disentangled the contribution of multiple migratory and adaptive events in shaping the heterogeneous Italian genomic background, which exemplify population dynamics and gene-environment interactions that played significant roles also in the formation of the Continental and Southern European genomic landscapes.

Detecting Variants in the NBN Gene While Testing for Hereditary Breast Cancer: What to Do Next?

The NBN gene has been included in breast cancer (BC) multigene panels based on early studies suggesting an increased BC risk for carriers, though not confirmed by recent research. To evaluate the impact of NBN analysis, we assessed the results of NBN sequencing in 116 BRCA-negative BC patients and reviewed the literature. Three patients (2.6%) carried potentially relevant variants: two, apparently unrelated, carried the frameshift variant c.156_157delTT and another one the c.628G>T variant. The latter was subsequently found in 4/1390 (0.3%) BC cases and 8/1580 (0.5%) controls in an independent sample, which, together with in silico predictions, provided evidence against its pathogenicity. Conversely, the rare c.156_157delTT variant was absent in the case-control set; moreover, a 50% reduction of NBN expression was demonstrated in one carrier. However, in one family it failed to co-segregate with BC, while the other carrier was found to harbor also a probably pathogenic TP53 variant that may explain her phenotype. Therefore, the c.156_157delTT, although functionally deleterious, was not supported as a cancer-predisposing defect. Pathogenic/likely pathogenic NBN variants were detected by multigene panels in 31/12314 (0.25%) patients included in 15 studies. The risk of misinterpretation of such findings is substantial and supports the exclusion of NBN from multigene panels.

Dissecting the Pre-Columbian Genomic Ancestry of Native Americans along the Andes-Amazonia Divide.

Extensive European and African admixture coupled with loss of Amerindian lineages makes the reconstruction of pre-Columbian history of Native Americans based on present-day genomes extremely challenging. Still open questions remain about the dispersals that occurred throughout the continent after the initial peopling from the Beringia, especially concerning the number and dynamics of diffusions into South America. Indeed, if environmental and historical factors contributed to shape distinct gene pools in the Andes and Amazonia, the origins of this East-West genetic structure and the extension of further interactions between populations residing along this divide are still not well understood. To this end, we generated new high-resolution genome-wide data for 229 individuals representative of one Central and ten South Amerindian ethnic groups from Mexico, Peru, Bolivia, and Argentina. Low levels of European and African admixture in the sampled individuals allowed the application of fine-scale haplotype-based methods and demographic modeling approaches. These analyses revealed highly specific Native American genetic ancestries and great intragroup homogeneity, along with limited traces of gene flow mainly from the Andes into Peruvian Amazonians. Substantial amount of genetic drift differentially experienced by the considered populations underlined distinct patterns of recent inbreeding or prolonged isolation. Overall, our results support the hypothesis that all non-Andean South Americans are compatible with descending from a common lineage, while we found low support for common Mesoamerican ancestors of both Andeans and other South American groups. These findings suggest extensive back-migrations into Central America from non-Andean sources or conceal distinct peopling events into the Southern Continent.

Inferring patterns of folktale diffusion using genomic data.

Observable patterns of cultural variation are consistently intertwined with demic movements, cultural diffusion, and adaptation to different ecological contexts [Cavalli-Sforza and Feldman (1981) Cultural Transmission and Evolution: A Quantitative Approach; Boyd and Richerson (1985) Culture and the Evolutionary Process]. The quantitative study of gene-culture coevolution has focused in particular on the mechanisms responsible for change in frequency and attributes of cultural traits, the spread of cultural information through demic and cultural diffusion, and detecting relationships between genetic and cultural lineages. Here, we make use of worldwide whole-genome sequences [Pagani et al. (2016) Nature 538:238-242] to assess the impact of processes involving population movement and replacement on cultural diversity, focusing on the variability observed in folktale traditions (n = 596) [Uther (2004) The Types of International Folktales: A Classification and Bibliography. Based on the System of Antti Aarne and Stith Thompson] in Eurasia. We find that a model of cultural diffusion predicted by isolation-by-distance alone is not sufficient to explain the observed patterns, especially at small spatial scales (up to [Formula: see text]4,000 km). We also provide an empirical approach to infer presence and impact of ethnolinguistic barriers preventing the unbiased transmission of both genetic and cultural information. After correcting for the effect of ethnolinguistic boundaries, we show that, of the alternative models that we propose, the one entailing cultural diffusion biased by linguistic differences is the most plausible. Additionally, we identify 15 tales that are more likely to be predominantly transmitted through population movement and replacement and locate putative focal areas for a set of tales that are spread worldwide.

Aneuploid acute myeloid leukemia exhibits a signature of genomic alterations in the cell cycle and protein degradation machinery.

BACKGROUND: Aneuploidy occurs in more than 20% of acute myeloid leukemia (AML) cases and correlates with an adverse prognosis. METHODS: To understand the molecular bases of aneuploid acute myeloid leukemia (A-AML), this study examined the genomic profile in 42 A-AML cases and 35 euploid acute myeloid leukemia (E-AML) cases. RESULTS: A-AML was characterized by increased genomic complexity based on exonic variants (an average of 26 somatic mutations per sample vs 15 for E-AML). The integration of exome, copy number, and gene expression data revealed alterations in genes involved in DNA repair (eg, SLX4IP, RINT1, HINT1, and ATR) and the cell cycle (eg, MCM2, MCM4, MCM5, MCM7, MCM8, MCM10, UBE2C, USP37, CK2, CK3, CK4, BUB1B, NUSAP1, and E2F) in A-AML, which was associated with a 3-gene signature defined by PLK1 and CDC20 upregulation and RAD50 downregulation and with structural or functional silencing of the p53 transcriptional program. Moreover, A-AML was enriched for alterations in the protein ubiquitination and degradation pathway (eg, increased levels of UHRF1 and UBE2C and decreased UBA3 expression), response to reactive oxygen species, energy metabolism, and biosynthetic processes, which may help in facing the unbalanced protein load. E-AML was associated with BCOR/BCORL1 mutations and HOX gene overexpression. CONCLUSIONS: These findings indicate that aneuploidy-related and leukemia-specific alterations cooperate to tolerate an abnormal chromosome number in AML, and they point to the mitotic and protein degradation machineries as potential therapeutic targets.

The genetic legacy of the Yaghnobis: A witness of an ancient Eurasian ancestry in the historically reshuffled central Asian gene pool.

OBJECTIVES: The Yaghnobis are an ethno-linguistic minority historically settled along the Yaghnob River in the Upper-Zarafshan Valley in Tajikistan. They speak a language of Old Sogdian origin, which is the only present-day witness of the Lingua Franca used along the Silk Road in Late Antiquity. The aim of this study was to reconstruct the genetic history of this community in order to shed light on its isolation and genetic ancestry within the Euro-Asiatic context. MATERIALS AND METHODS: A total of 100 DNA samples were collected in the Yaghnob and Matcha Valleys during several expeditions and their mitochondrial, Y-chromosome and autosomal genome-wide variation were compared with that from a large set of modern and ancient Euro-Asiatic samples. RESULTS: Findings from uniparental markers highlighted the long-term isolation of the Yaghnobis. Mitochondrial DNA ancestry traced an ancient link with Middle Eastern populations, whereas Y-chromosome legacy showed more tight relationships with Central Asians. Admixture, outgroup-f3, and D-statistics computed on autosomal variation corroborated Y-chromosome evidence, pointing respectively to low Anatolian Neolithic and high Steppe ancestry proportions in Yaghnobis, and to their closer affinity with Tajiks than to Iranians. DISCUSSION: Although the Yaghnobis do not show evident signs of recent admixture, they could be considered a modern proxy for the source of gene flow for many Central Asian and Middle Eastern groups. Accordingly, they seem to retain a peculiar genomic ancestry probably ascribable to an ancient gene pool originally wide spread across a vast area and subsequently reshuffled by distinct demographic events occurred in Middle East and Central Asia.

A novel founder MYO15A frameshift duplication is the major cause of genetic hearing loss in Oman.

The increased risk for autosomal recessive disorders is one of the most well-known medical implications of consanguinity. In the Sultanate of Oman, a country characterized by one of the highest rates of consanguineous marriages worldwide, prevalence of genetic hearing loss (GHL) is estimated to be 6/10 000. Families of GHL patients have higher consanguinity rates than the general Omani population, indicating a major role for recessive forms. Mutations in GJB2, the most commonly mutated GHL gene, have been sporadically described. We collected 97 DNA samples of GHL probands, affected/unaffected siblings and parents from 26 Omani consanguineous families. Analyzing a first family by whole-exome sequencing, we identified a novel homozygous frameshift duplication (c.1171_1177dupGCCATCT) in MYO15A, the gene linked to the deafness locus DFNB3. This duplication was then found in a total of 8/26 (28%) families, within a 849 kb founder haplotype. Reconstruction of haplotype structure at MYO15A surrounding genomic regions indicated that the founder haplotype branched out in the past two to three centuries from a haplotype present worldwide. The MYO15A duplication emerges as the major cause of GHL in Oman. These findings have major implications for the design of GHL diagnosis and prevention policies in Oman.

Inferring chronological age from DNA methylation patterns of human teeth.

OBJECTIVE: Current methods to determine chronological age from modern and ancient remains rely on both morphological and molecular approaches. However, low accuracy and the lack of standardized protocols make the development of alternative methods for the estimation of individual's age even more urgent for several research fields, such as biological anthropology, biodemography, forensics, evolutionary genetics, and ancient DNA studies. Therefore, the aim of this study is to identify genomic regions whose DNA methylation level correlates with age in modern teeth. METHODS: We used MALDI-TOF mass spectrometry to analyze DNA methylation levels of specific CpGs located in the ELOVL2, FHL2, and PENK genes. We considered methylation data from cementum, dentin and pulp of 21 modern teeth (from 17 to 77 years old) to construct a mathematical model able to exploit DNA methylation values to predict age of the individuals. RESULTS: The median difference between the real age and that estimated using DNA methylation values is 1.20 years (SD = 1.9) if DNA is recovered from both cementum and pulp of the same modern teeth, 2.25 years (SD = 2.5) if DNA is recovered from dental pulp, 2.45 years (SD = 3.3) if DNA is extracted from cementum and 7.07 years (SD = 7.0) when DNA is recovered from dentin only. DISCUSSION: We propose for the first time the evaluation of DNA methylation at ELOVL2, FHL2, and PENK genes as a powerful tool to predict age in modern teeth for anthropological applications. Future studies are needed to apply this method also to historical and relatively ancient human teeth.

Positive selection of lactase persistence among people of Southern Arabia.

OBJECTIVE: Frequency patterns of the lactase persistence (LP)-associated -13,915 G allele and archaeological records pointing to substantial role played by southern regions in the peopling and domestication processes that involved the Arabian Peninsula suggest that Southern Arabia plausibly represented the center of diffusion of such adaptive variant. Nevertheless, a well-defined scenario for evolution of Arabian LP is still to be elucidated and the burgeoning archaeological picture of complex human migrations occurred through the peninsula is not matched by an equivalent high-resolution description of genetic variation underlying this adaptive trait. To fill this gap, we investigated diversity at a wide genomic interval surrounding the LCT gene in different Southern Arabian populations. METHODS: 40 SNPs were genotyped to characterize LCT profiles of 630 Omani and Yemeni individuals to perform population structure, linkage disequilibrium, population differentiation-based and haplotype-based analyses. RESULTS: Typical Arabian LP-related variation was found in Dhofaris and Yemenis, being characterized by private haplotypes carrying the -13,915 G allele, unusual differentiation with respect to northern groups and conserved homozygous haplotype-blocks, suggesting that the adaptive allele was likely introduced in the Arabian gene pool in southern populations and was then subjected to prolonged selective pressure. CONCLUSION: By pointing to Yemen as one of the best candidate centers of diffusion of the Arabian-specific adaptive variant, obtained results indicate the spread of indigenous groups as the main process underlying dispersal of LP along the Arabian Peninsula, supporting a refugia model for Arabian demic movements occurred during the Terminal Pleistocene and Early Holocene.

Tracing Behçet's disease origins along the Silk Road: an anthropological evolutionary genetics perspective.

OBJECTIVES: Behçet's disease is a multifactorial vasculitis that shows its highest prevalence in geographical areas historically involved in the Silk Road, suggesting that it might have originated somewhere along these ancient trade routes. This study aims to provide a first clue towards genetic evidence for this hypothesis by testing it via an anthropological evolutionary genetics approach. METHODS: Behçet's disease variation at ancestry informative mitochondrial DNA control region and haplogroup diagnostic sites was characterised in 185 disease subjects of Italian descent and set into the Eurasian mitochondrial landscape by comparison with nearly 9,000 sequences representative of diversity observable in Italy and along the main Silk Road routes. RESULTS: Dissection of the actual genetic ancestry of disease individuals by means of population structure, spatial autocorrelation and haplogroup analyses revealed their closer relationships with some Middle Eastern and Central Asian groups settled along the Silk Road than with healthy Italians. CONCLUSIONS: These findings support the hypothesis that the Behçet's disease genetic risk has migrated to western Eurasia in parallel with ancestry components typical of Silk Road-related groups. This provided new insights that are useful to improve the understanding of disease origins and diffusion, as well as to inform future association studies aimed at properly accounting for the actual genetic ancestry of the examined Behçet's disease samples in order to minimise the detection of spurious associations and to improve the identification of genetic variants with actual clinical relevance.

Inferring the genetic history of lactase persistence along the Italian peninsula from a large genomic interval surrounding the LCT gene.

OBJECTIVE: Although genetic variants related to lactase persistence in European populations were supposed to have firstly undergone positive selection in farmers from the Balkans and Central Europe, demographic and evolutionary dynamics that subsequently shaped the distribution of this adaptive trait across the continent have still to be elucidated. To deepen the knowledge about potential routes of diffusion of lactase persistence to Western Europe we investigated variation at a large genomic region surrounding the LCT gene along the Italian peninsula, a geographical area that played a key role in population movements responsible for Neolithic diffusion across Europe. METHODS: By genotyping 40 highly selected SNPs in more than 400 Italian individuals we described gradients of nucleotide and haplotype variation potentially related to lactase persistence and compared them with those observed in several European and Mediterranean human groups. RESULTS: Multiple migratory events responsible for earlier introduction of the examined alleles in Italy than in Northern European regions could be invoked. Different demic processes occurred along the western and eastern sides of the peninsula were also inferred via linkage disequilibrium and population structure analyses. CONCLUSION: The appreciable genetic continuum observed between people from Northern or Central-Western Italy and Central European populations suggested a local arrival of lactase persistence-related variants mainly via overland routes. On the contrary, diversity of Central-Eastern and Southern Italian groups entailed also gene flow from South-Eastern Mediterranean regions, in accordance to the earlier entrance of the Neolithic in Southern Italy via maritime population movements along the Mediterranean coastlines.

The epigenetic side of human adaptation: hypotheses, evidences and theories.

CONTEXT: Epigenetics represents a still unexplored research field in the understanding of micro- and macro-evolutionary mechanisms, as epigenetic changes create phenotypic diversity within both individuals and populations. OBJECTIVE: The purpose of this review is to dissect the landscape of studies focused on DNA methylation, one of the most described epigenetic mechanisms, emphasizing the aspects that could be relevant in human adaptations. METHODS: Theories and results here considered were collected from the most recent papers published. RESULTS: The matter of DNA methylation inheritance is here described as well as the recent evolutionary theories regarding the role of DNA methylation-and epigenetics in a broader sense-in human evolution. The complex relation between (1) DNA methylation and genetic variability and (2) DNA methylation and the environmental stimuli crucial in shaping genetic and phenotypic variability through the human lineage-such as diet, climate and pathogens exposure-are described. Papers about population epigenetics are also illustrated due to their high relevance in this context. CONCLUSION: Genetic, epigenetic and phenotypic variations of the species, together with cultural ones, are considerably shaped by a vast range of environmental stimuli, thus representing the foundation of all human bio-cultural adaptations.

mtDNA variation in East Africa unravels the history of Afro-Asiatic groups.

East Africa (EA) has witnessed pivotal steps in the history of human evolution. Due to its high environmental and cultural variability, and to the long-term human presence there, the genetic structure of modern EA populations is one of the most complicated puzzles in human diversity worldwide. Similarly, the widespread Afro-Asiatic (AA) linguistic phylum reaches its highest levels of internal differentiation in EA. To disentangle this complex ethno-linguistic pattern, we studied mtDNA variability in 1,671 individuals (452 of which were newly typed) from 30 EA populations and compared our data with those from 40 populations (2970 individuals) from Central and Northern Africa and the Levant, affiliated to the AA phylum. The genetic structure of the studied populations--explored using spatial Principal Component Analysis and Model-based clustering--turned out to be composed of four clusters, each with different geographic distribution and/or linguistic affiliation, and signaling different population events in the history of the region. One cluster is widespread in Ethiopia, where it is associated with different AA-speaking populations, and shows shared ancestry with Semitic-speaking groups from Yemen and Egypt and AA-Chadic-speaking groups from Central Africa. Two clusters included populations from Southern Ethiopia, Kenya and Tanzania. Despite high and recent gene-flow (Bantu, Nilo-Saharan pastoralists), one of them is associated with a more ancient AA-Cushitic stratum. Most North-African and Levantine populations (AA-Berber, AA-Semitic) were grouped in a fourth and more differentiated cluster. We therefore conclude that EA genetic variability, although heavily influenced by migration processes, conserves traces of more ancient strata.

IL-23A, IL-23R, IL-17A and IL-17R polymorphisms in different psoriatic arthritis clinical manifestations in the northern Italian population.

To investigate the genetic variability of IL-17A, IL17-RA, IL-23A and IL-23R genes on an in-depth phenotypically characterized northern Italian Psoriatic arthritis (PsA) case-control cohort, in search for associations specific to different PsA clinical sub-phenotypes. We examined 118 patients with PsA according to CASPAR criteria (mean age 57 ± 13, female 38.4 %, mean disease duration 13.9 ± 8.6 years, peripheral disease 83.8 %, axial manifestations 34.5 %, radiological erosive disease 49 %) compared with 248 controls of the same ethnic origin matched for age and sex. The presence of axial disease was defined by the clinical axial involvement and/or the presence of radiological alteration consistent with spondyloarthropathy according to New York criteria. The presence of peripheral disease (arthritis and/or enthesitis) was defined only on clinical basis. A total of 40 SNPs, mapping within the genes mentioned above, were genotyped in both groups and used to perform association analyses by subdividing the PsA sample into subgroups according to different clinical manifestations on the basis of axial and peripheral involvements. No differences between patients and controls were found in the distribution of the IL-17A, IL17-RA, IL-23A and IL-23R genes allelic variants. Comparing patients with axial disease versus those without, we found that axial manifestations were significantly associated with the presence of IL-23R rs12401432 GG homozygosity (26.8 % vs. 5.3 %, p corr = 0.019, OR 2.63 [95 % CI 1.13-6.16]). No differences in distribution of the allelic variants were found comparing patients with versus those without peripheral disease or patients with versus without radiological peripheral erosions. In PA patients of northern Italian origin, IL-17A, IL17-RA, IL-23A and IL-23R genes allelic variants are not associated with disease susceptibility. However, a strong association with the IL-23RA rs12401432 GG genotype is associated with axial involvement of the disease.

Epigenetic aging differences between Wichí and Criollos from Argentina: Insights from genomic history and ecology.

Background and objectives: Epigenetic estimators based on DNA methylation levels have emerged as promising biomarkers of human aging. These estimators exhibit natural variations across human groups, but data about indigenous populations remain underrepresented in research. This study aims to investigate differences in epigenetic estimators between two distinct human populations, both residing in the Gran Chaco region of Argentina, the Native-American Wichí, and admixed Criollos who are descendants of intermarriages between Native Americans and the first European colonizers, using a population genetic approach. Methodology: We analyzed 24 Wichí (mean age: 39.2 ± 12.9 yo) and 24 Criollos (mean age: 41.1 ± 14.0 yo) for DNA methylation levels using the Infinium MethylationEPIC (Illumina) to calculate 16 epigenetic estimators. Additionally, we examined genome-wide genetic variation using the HumanOmniExpress BeadChip (Illumina) to gain insights into the genetic history of these populations. Results: Our results indicate that Native-American Wichí are epigenetically older compared to Criollos according to five epigenetic estimators. Analyses within the Criollos population reveal that global ancestry does not influence the differences observed, while local (chromosomal) ancestry shows positive associations between specific SNPs located in genomic regions over-represented by Native-American ancestry and measures of epigenetic age acceleration (AgeAccelHannum). Furthermore, we demonstrate that differences in population ecologies also contribute to observed epigenetic differences. Conclusions and implications: Overall, our study suggests that while the genomic history may partially account for the observed epigenetic differences, non-genetic factors, such as lifestyle and ecological factors, play a substantial role in the variability of epigenetic estimators, thereby contributing to variations in human epigenetic aging.

Imputation reliability on DNA biallelic markers for drug metabolism studies.

BACKGROUND: Imputation is a statistical process used to predict genotypes of loci not directly assayed in a sample of individuals. Our goal is to measure the performance of imputation in predicting the genotype of the best known gene polymorphisms involved in drug metabolism using a common SNP array genotyping platform generally exploited in genome wide association studies. METHODS: Thirty-nine (39) individuals were genotyped with both Affymetrix Genome Wide Human SNP 6.0 (AFFY) and Affymetrix DMET Plus (DMET) platforms. AFFY and DMET contain nearly 900000 and 1931 markers respectively. We used a 1000 Genomes Pilot + HapMap 3 reference panel. Imputation was performed using the computer program Impute, version 2. SNPs contained in DMET, but not imputed, were analysed studying markers around their chromosome regions. The efficacy of the imputation was measured evaluating the number of successfully imputed SNPs (SSNPs). RESULTS: The imputation predicted the genotypes of 654 SNPs not present in the AFFY array, but contained in the DMET array. Approximately 1000 SNPs were not annotated in the reference panel and therefore they could not be directly imputed. After testing three different imputed genotype calling threshold (IGCT), we observed that imputation performs at its best for IGCT value equal to 50%, with rate of SSNPs (MAF > 0.05) equal to 85%. CONCLUSIONS: Most of the genes involved in drug metabolism can be imputed with high efficacy using standard genome-wide genotyping platforms and imputing procedures.

Polymorphic NumtS trace human population relationships.

The human genome is constantly subjected to evolutionary forces which shape its architecture. Insertions of mitochondrial DNA sequences into nuclear genome (NumtS) have been described in several eukaryotic species, including Homo sapiens and other primates. The ongoing process of the generation of NumtS has made them valuable markers in primate phylogenetic studies, as well as potentially informative loci for reconstructing the genetic history of modern humans. Here, we report the identification of 53 human-specific NumtS by inspection of the UCSC genome browser, showing that they may be direct insertions of mitochondrial DNA into the human nuclear DNA after the human-chimpanzee split. In silico analyses allowed us to identify 14 NumtS which are polymorphic in terms of their presence/absence within the human genome in individuals of different ancestry. The allele frequencies of these polymorphic NumtS were calculated for 1000 Genomes Project sequence data from 13 populations worldwide, and principal components analysis and hierarchical clustering methods allowed the detection of strong signals of geographical structure related to the genetic diversity of these loci. All identified polymorphic human-specific NumtS together with a tandemly duplicated NumtS have also been validated by PCR amplification on a panel of 60 samples belonging to five native populations worldwide, confirming the expected NumtS variability. On the basis of these findings, we have succeeded in depicting the landscape of variation of a series of NumtS in several ethnic groups, making an advance in their identification as useful markers in the study on human population genetics.