Distribution of ABCB1 polymorphisms among Brazilians: impact of population admixture

Introduction: Interethnic admixture is a source of cryptic population structure that may lead to spurious genotype-phenotype associations in pharmacogenomic studies. We studied the impact of population stratification on the distribution of ABCB1 polymorphisms (1236C > T, 2677G > T/A and 3435C > T) among Brazilians, a highly admixed population with Amerindian, European and African ancestral roots. Methods: Individual DNA from 320 healthy adults was genotyped with a panel of ancestry informative markers, and the proportions of African component of ancestry (ACA) were estimated. ABCB1 genotypes were determined by the single base extension/termination method. We describe the association between ABCB1 polymorphisms and ACA by fitting a linear proportional odds logistic regression model to the data. Results: The distribution of the ABCB1 2677G > T/A and 3435C > T, but not the 1236C > T, SNPs displayed a significant trend for decreasing frequency of the T alleles and TT genotypes from White to Intermediate to Black individuals. The same trend was observed in the frequency of the T/nonG/T haplotype at the 1236, 2677 and 3435 loci. When the population sample was proportioned in quartiles, according to the individual ACA estimates, the frequency of the T allele and TT genotype at each locus declined progressively from the lowest (< 0.25 ACA) to the highest (> 0.75 ACA) quartile. Linear proportional odds logistic regression analysis confirmed that the odds of having the T allele at each locus decreases in a continuous manner with the increase of the ACA, throughout the ACA range (0.13-0.94) observed in the overall population sample. A significant association was also detected between the individual ACA estimates and the presence of the T/nonG/T haplotype in the overall population. Conclusion: Self-identification according to the racial/color categories proposed by the Brazilian Census is insufficient to properly control for population stratification in pharmacogenomic studies of ABCB1.

A single multiplex PCR and SNaPshot minisequencing reaction of 42 SNPs to classify admixture populations into mitochondrial DNA haplogroups

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A single multiplex SNaPshot reaction of 42 SNPs to classify admixture populations into mitochondrial DNA haplogroups

SNaPshot minisequencing reaction is in increasing use because of its fast detection of many polymorphisms in a single assay. In this work we described a highly sensitive single nucleotide polymorphisms (SNPs) typing method with detection of 42 mitochondrial DNA (mtDNA) SNPs in a single PCR and SNaPshot multiplex reaction in order to allow haplogroup classification in Latin American admixture population. We validated the panel typing 160 Brazilian individuals. DNA was extracted from blood spotted on filter paper using Chelex protocol. Forty SNPs were selected targeting haplogroup-specific mutations in Europeans, Africans and Asians (only precursors of Native Americans haplogroups A2, B2, C1, and D1) and two non-coding SNPs were chosen to increase the power of discrimination between individuals (SNPs positions 16,519 and 16,362). It was done using a modified version of a previously published multiplex SNaPshot minisequencing reaction established to resolve European haplogroups, adding SNPs targeting Africans (L0, L1, L2, L3, and L*) and Asians (A, B, C, and D) haplogroups based on SNPs described at PhyloTree.org build 2. PCR primers were designed using PerlPrimer software and checked with the Autodimer program. Thirty-three primer-pairs were used to amplify 42 SNPs. Using this panel, we were able to successfully classify 160 individuals into their correct haplogroups. Complete SNP profiles were obtained from 10. pg of total DNA. We conclude that it is possible to build and genotype more than 40 mtDNA SNPs in a single multiplex PCR and SNaPshot reaction, with sensitivity and reliability, resolving haplogroup classification in admixture populations. © 2011.

IRF6 is a Risk Factor for Nonsyndromic Cleft Lip in the Brazilian Population

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a complex disorder with a worldwide incidence estimated at 1:700. Among the putative susceptibility loci, the IRF6 gene and a region at 8q24.21 have been corroborated in different populations. To test the role of IRF6 in NSCL/P predisposition in the Brazilian population, we conducted a structured association study with the SNPs rs642961 and rs590223, respectively, located at 5' and 3' of the IRF6 gene and not in strong linkage disequilibrium (LD), in patients from five different Brazilian locations. We also evaluated the effect of these SNPs in IRF6 expression in mesenchymal stem cells (MSC). We observed association between rs642961 and cleft lip only (CLO) (P = 0.009; odds ratio (OR) for AA genotype = 1.83 [95% Confidence interval (CI), 0.64-5.31]; OR for AG genotype = 1.72 [95% CI, 1.03-2.84]). This association seems to be driven by the affected patients from Barbalha, a location which presents the highest heritability estimate (H-2 = 0.85), and the A allele at rs642961 is acting through a dominant model. No association was detected for the SNP rs590223. We did not find any correlation between expression levels and genotypes of the two loci, and it is possible that these SNPs have a functional role in some specific period of embryogenesis. (C) 2012 Wiley Periodicals, Inc.

Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

Influência do polimorfismo genético de citocinas na hepatite C crônica em uma população do Rio de Janeiro

A infecção pelo vírus da hepatite C (VHC) é uma das mais comuns infecções ao redor do mundo. Aproximadamente, 20% dos pacientes infectados eliminam espontaneamente o vírus, porém a maioria dos indivíduos infectados desenvolve infecção crônica com amplo espectro de lesões hepáticas, desde inflamação leve até cirrose. A resposta imune do hospedeiro exerce grande influência sobre o desfecho da infecção pelo VHC. O objetivo deste trabalho foi analisar a influência dos polimorfismos genéticos de citocinas na susceptibilidade ou persistência da infecção por VHC e no clareamento espontâneo em uma amostra de pacientes da população do Rio de Janeiro (Brasil). Os polimorfismos genéticos das citocinas TNFA (-308), TGFB1 (codon 10 e 25), IL10 (-1082, -592), IL6 (-174) e IFNG (+874) foram analisados por PCR-SSP em 245 pacientes com hepatite C crônica (HCC), 41 pacientes que alcançaram o clareamento viral espontâneo e 189 indivíduos controle saudáveis. Além disso, os polimorfismos próximos ao gene da citocina IL28B (rs12979860, rs12980275 e rs8099917) foram analisados por PCR em tempo real em todos os grupos. O grau de fibrose e inflamação, a resposta ao tratamento e o genótipo do vírus também foram levados em consideração quanto ao desfecho da HCC Os genótipos IL28B rs12979860 CC e CT e rs12980275 AA e AG foram significativamente associados ao clareamento espontâneo e à resposta à terapia anti-viral. Da mesma forma, o alelo C (rs12979860) e o alelo A (rs12980275) foram significativamente maior no grupo Clareamento. O alelo C de IL6 (-174) foi associado com o Clareamento. Nenhuma associação entre as demais citocinas e o desfecho da HCC foi encontrada. O Genótipo TNFA (-308) GG parece estar associado com menor grau de inflamação. Além disso, a etnia auto declarada influencia a distribuição dos polimorfismos em IL6 (-174) e IL28B rs12979860 e rs8099917. Nossas observações indicam que os polimorfismos próximos ao gene da IL28B estão associados com o clareamento viral e resposta ao tratamento na população do Rio de Janeiro. Além disso, nossos resultados podem ser úteis para futuras investigações entre os polimorfismos de citocinas e a infecção por VHC numa população heterogênea como a Brasileira.

European and Polynesian admixture in the Norfolk Island population

The Norfolk Island population in the South Pacific is primarily the product of recent admixture between a small number of British male and Polynesian female founders. We identified and genotyped 128 Ancestry Informative Markers (AIMs) spread across the autosomes, X/Y chromosomes and mitochondrial DNA genome, to explore and quantify the current levels of genetic admixture in the Norfolk Islanders. On the basis of autosomal AIMs, the population shows mean European and Polynesian ancestry proportions of 88 and 12%, respectively. However, there is a substantial variation between individuals ranging from total European ancestry to near total Polynesian origin. There is a strong correlation between individual genetic estimates of Polynesian ancestry and those derived from the extensive pedigree and genealogical records of Islanders. Also in line with historical accounts, there is a substantial asymmetry in the maternal and paternal origins of the Islanders with almost all Y-chromosomes of European origin whereas at least 25% of mtDNAs appear to have a Polynesian origin. Accurate knowledge of ancestry will be important in future attempts to use the Island population in admixture mapping approaches to find the genes that underlie differences in the risk to some diseases between Europeans and Polynesians.

A novel parallel approach to the likelihood-based estimation of admixture in population genetics

Inferring population admixture from genetic data and quantifying it is a difficult but crucial task in evolutionary and conservation biology. Unfortunately state-of-the-art probabilistic approaches are computationally demanding. Effectively exploiting the computational power of modern multiprocessor systems can thus have a positive impact to Monte Carlo-based simulation of admixture modeling. A novel parallel approach is briefly described and promising results on its message passing interface (MPI)-based C implementation are reported.

Brazilian urban population genetic structure reveals a high degree of admixture

Advances in genotyping technologies have contributed to a better understanding of human population genetic structure and improved the analysis of association studies. To analyze patterns of human genetic variation in Brazil, we used SNP data from 1129 individuals - 138 from the urban population of Sao Paulo, Brazil, and 991 from 11 populations of the HapMap Project. Principal components analysis was performed on the SNPs common to these populations, to identify the composition and the number of SNPs needed to capture the genetic variation of them. Both admixture and local ancestry inference were performed in individuals of the Brazilian sample. Individuals from the Brazilian sample fell between Europeans, Mexicans, and Africans. Brazilians are suggested to have the highest internal genetic variation of sampled populations. Our results indicate, as expected, that the Brazilian sample analyzed descend from Amerindians, African, and/or European ancestors, but intermarriage between individuals of different ethnic origin had an important role in generating the broad genetic variation observed in the present-day population. The data support the notion that the Brazilian population, due to its high degree of admixture, can provide a valuable resource for strategies aiming at using admixture as a tool for mapping complex traits in humans. European Journal of Human Genetics (2012) 20, 111-116; doi:10.1038/ejhg.2011.144; published online 24 August 2011

Conservation genetics of European wildcat (Felis silvestris silvestris): a wide and integrating analysis protocol for admixture inferences and population structure

Introgression of domestic cat genes into European wildcat (Felis silvestris silvestris) populations and reduction of wildcats’ range in Europe, leaded by habitat loss and fragmentation, are considered two of the main conservation problems for this endangered feline. This thesis addressed the questions related with the artificial hybridization and populations’ fragmentation, using a conservation genetics perspective. We combined the use of highly polymorphic loci, Bayesian statistical inferences and landscape analyses tools to investigate the origin of the geographic-genetic substructure of European wildcats (Felis silvestris silvestris) in Italy and Europe. The genetic variability of microsatellites evidenced that European wildcat populations currently distributed in Italy differentiated in, and expanded from two distinct glacial refuges during the Last Glacial Maximum. The genetic and geographic substructure detected between the eastern and western sides of the Apennine ridge, resulted by adaptation to specific ecological conditions of the Mediterranean habitats. European wildcat populations in Europe are strongly structured into 5 geographic-genetic macro clusters corresponding to: the Italian peninsular & Sicily; Balkans & north-eastern Italy; Germany eastern; central Europe; and Iberian Peninsula. Central European population might have differentiated in the extra-Mediterranean Würm ice age refuge areas (Northern Alps, Carpathians, and the Bulgarian mountain systems), while the divergence among and within the southern European populations might have resulted by the Pleistocene bio geographical framework of Europe, with three southern refugia localized in the Balkans, Italian Peninsula and Iberia Peninsula. We further combined the use of most informative autosomal SNPs with uniparental markers (mtDNA and Y-linked) for accurately detecting parental genotypes and levels of introgressive hybridization between European wild and domestic cats. A total of 11 hybrids were identified. The presence of domestic mitochondrial haplotypes shared with some wild individuals led us to hypnotize the possibility that ancient introgressive events might have occurred and that further investigation should be recommended.

Legacy of mutiny on the Bounty : founder effect and admixture on Norfolk Island

The population of Norfolk Island, located off the eastern coast of Australia, possesses an unusual and fascinating history. Most present-day islanders are related to a small number of the 'Bounty' mutineer founders. These founders consisted of Caucasian males and Polynesian females and led to an admixed present-day population. By examining a single large pedigree of 5742 individuals, spanning >200 years, we analyzed the influence of admixture and founder effect on various cardiovascular disease (CVD)-related traits. On account of the relative isolation of the population, on average one-third of the genomes of present-day islanders (single large pedigree individuals) is derived from 17 initial founders. The proportion of Polynesian ancestry in the present-day individuals was found to significantly influence total triglycerides, body mass index, systolic blood pressure and diastolic blood pressure. For various cholesterol traits, the influence of ancestry was less marked but overall the direction of effect for all CVD-related traits was consistent with Polynesian ancestry conferring greater CVD risk. Marker-derived homozygosity was computed and agreed with measures of inbreeding derived from pedigree information. Founder effect (inbreeding and marker-derived homozygosity) significantly influenced height. In conclusion, both founder effect and extreme admixture have substantially influenced the genetic architecture of a variety of CVD-related traits in this population.

Population structure, effective population size and adverse effects of stocking in the endangered Australian eastern freshwater cod Maccullochella ikei

Microsatellite markers were used to examine spatio-temporal genetic variation in the endangered eastern freshwater cod Maccullochella ikei in the Clarence River system, eastern Australia. High levels of population structure were detected. A model-based clustering analysis of multilocus genotypes identified four populations that were highly differentiated by F-statistics (FST = 0· 09 − 0· 49; P < 0· 05), suggesting fragmentation and restricted dispersal particularly among upstream sites. Hatchery breeding programmes were used to re-establish locally extirpated populations and to supplement remnant populations. Bayesian and frequency-based analyses of hatchery fingerling samples provided evidence for population admixture in the hatchery, with the majority of parental stock sourced from distinct upstream sites. Comparison between historical and contemporary wild-caught samples showed a significant loss of heterozygosity (21%) and allelic richness (24%) in the Mann and Nymboida Rivers since the commencement of stocking. Fragmentation may have been a causative factor; however, temporal shifts in allele frequencies suggest swamping with hatchery-produced M. ikei has contributed to the genetic decline in the largest wild population. This study demonstrates the importance of using information on genetic variation and population structure in the management of breeding and stocking programmes, particularly for threatened species.

Living on the edge : Population genetics of Finno-Ugric-speaking humans in North Eurasia

The major aim of this thesis was to examine the origins and distribution of uniparental and autosomal genetic variation among the Finno-Ugric-speaking human populations living in Boreal and Arctic regions of North Eurasia. In more detail, I aimed to disentangle the underlying molecular and population genetic factors which have produced the patterns of uniparental and autosomal genetic diversity in these populations. Among Finno-Ugrics the genetic amalgamation and clinal distribution of West and East Eurasian gene pools were observed within uniparental markers. This admixture indicates that North Eurasia was colonized through Central Asia/ South Siberia by human groups already carrying both West and East Eurasian lineages. The complex combination of founder effects, gene flow and genetic drift underlying the genetic diversity of the Finno-Ugric- speaking populations were emphasized by low haplotype diversity within and among uniparental and biparental markers. A high prevalence of lactase persistence allele among the North Eurasian Finno- Ugric agriculturalist populations was also shown indicating a local adaptation to subsistence change with lactose rich diet. Moreover, the haplotype background of lactase persistence allele among the Finno- Ugric-speakers strongly suggested that the lactase persistence T-13910 mutation was introduced independently more than once to the North Eurasian gene pool. A significant difference in genetic diversity, haplotype structure and LD distribution within the cytochrome P450 CYP2C and CYP2D regions revealed the unique gene pool of the Finno-Ugric Saami created mainly by population genetic processes compared to other Europeans and sub-Saharan Mandenka population. From all studied populations the Saami showed also significantly the highest allele frequency of a CYP2C19 gene mutation causing variable drug reactions. The diversity patterns observed within CYP2C and CYP2D regions emphasize the strong effect of demographic history shaping genetic diversity and LD especially among such small and constant size populations as the Finno-Ugric-speaking Saami. Moreover, the increased LD in Saami due to genetic drift and/or admixture was shown to offer an advantage for further attempts to identify alleles associated to common complex pharmacogenetic traits.