HLA-DQA1 polymorphism in autochthonous Basques from Navarre (Spain): genetic position within European and Mediterranean scopes

Allele frequencies of 15 STR loci (Identifiler™ kit) in Basque-Americans

Individuals with Basque ancestry form a historically and culturally important minority of the population of the western United States. Allele frequencies for the 15 autosomal STRs in the AmpFlSTR® Identifiler® PCR Amplification Kit (Applied Biosystems) from 156 unrelated self-identified Basque individuals born in the United States are presented. Allele frequencies were used to calculate parameters commonly used in genetics and forensics including power of discrimination (PD), power of exclusion (PE), polymorphic information content (PIC), and expected heterozygosity (He). The sample population was also compared with the European Basque population and the major American ethnicities.

Ancestry markers from the human chromosome 6: Alu repeats from the MHC in autochthonous Basques

Polymorphic Alu insertions from the MHC class I region were analyzed in 215 autochthonous Basques from Guipuzcoa and Navarre provinces, with the aim of contributing new MHC Alu data in European ancestry populations. We also seek to assess both the genetic position of native Basques among worldwide samples and the efficiency of the MHC Alu elements as ancestry informative markers (AIMs). According to the MDS and AMOVA results, worldwide populations included in the comparative analyses were grouped in three major clusters defined by genetic ancestry (Africans, Asians and Europeans). The δ values (differences in weighted allele frequencies) among ancestry groups indicated that Alu elements within the alpha-block (AluHF, AluHJ and AluHG) showed an adequate resolving power to discriminate appropriately between some of the major ancestry groups. Alpha block Alu were also revealing of the exceptionality of Basques, as they allowed for the detection of genetic heterogeneity even between Basques and the other Iberian collection considered in the analysis (Valencia). Thus, analysis of the Alu loci within the alpha-block may represent a reliable, informative and cost-effective method to explore the ancestry, geographic origins and demographic history of human populations, which can be very helpful for studies into epidemiological, forensic or evolutionary perspectives.

Autosomal short tandem repeat genetic variation of the Basques in Spain

Aim

To examine population genetic structure and hypotheses of the origin of the modern Basque population in Spain using autosomal short tandem repeat (STR) data from individuals living in 27 mountain villages in the provinces of Alava, Vizcaya, Guipuzcoa, and Navarre, by comparing Basque autosomal STR variation with that of neighboring populations in Europe, as well as proposed ancestral populations in North Africa and the Caucasus.

Methods

Allele frequencies for 9 autosomal STR loci (D3S1358, D5S818, D7S820, D8S1179, D13S317, D18S51, D21S11, FGA, and vWA) and several population genetic parameters were determined for the 4 provinces in the Basque region of Spain (n = 377). Heterozygosity within the Basque population was measured using a locus-by-locus analysis of molecular variance. Relationships between the Basques and other populations were examined using a multidimensional scaling (MDS) plot of Shriver’s D_SW distance matrix.

Results

Heterozygosity levels in the Basque provinces were on the low end of the European distribution (0.805-0.812). The MDS plot of genetic distances revealed that the Basques differed from both the Caucasian and North African populations with respect to autosomal STR variation.

Conclusions

Autosomal STR analysis does not support the hypotheses of a recent common ancestor between the Basques and populations either from the Caucasus or North Africa.

Microsatellites and Alu elements from the human MHC in Valencia (Spain): analysis of genetic relationships and linkage disequilibrium

Two different sets of noncoding markers (microsatellites and Alu elements) from the human chromosome six were analysed in 106 individuals from Valencia (Spain), with the aim of exploring the effect of evolutionary forces on the genetic variability of the major histocompatibility complex (MHC) and assessing the potential usefulness of these genetic loci in phylogenetic studies. Linkage disequilibrium (LD) analyses revealed statistically significant associations among markers located in the MHC class I region, and also between the microsatellite D6S2792 and several genetic loci from MHC class I, II and III regions. Results of the Ewens-Watterson test indicated that only D6S2792 showed significant departure from selective neutrality. Despite the paucity of haplotype data in the literature, results of the phylogenetic analyses at world scale (Alu elements) showed that the genetic relationships of Valencia were mainly determined by the ethnic ancestry of the populations considered, whereas at European scale (microsatellites) population affinities were strongly influenced by geography. Our findings suggest that noncoding markers from the MHC such as Alu and microsatellite loci might have a potential value as lineage (ancestry) markers in investigations into evolutionary, medical and forensic perspectives.

An evolutionary approach to the high frequency of the Delta F508 CFTR mutation in European populations

The diffusion of the cattle pastoralism across Europe during the Neolithic period was probably accompanied by the emergence and spread of diverse contagious diseases that were unknown in the Paleolithic and that would have affected the frequency of genes directly or indirectly associated with differential susceptibility and/or resistance to infectious pathogens. We therefore propose that the high frequency of the CFTR gene, and in particular, the common Delta F508 allele mutation in current European and European-derived populations might be a consequence of the impact of selective pressures generated by the transmission of pathogenic agents from domesticated animals, mainly bovine cattle, to the man. Intestinal infectious diseases were probably a major health problem for Neolithic peoples. In such a context, a gene mutation that conferred an increased resistance to the diseases caused by pathogens transmitted by dairy cattle would have constituted a definite selective advantage, particularly in those human groups where cow's milk became an essential component of the diet. This selective advantage would be determined by an increased resistance to Cl(-)-secreting diarrheas of those individuals carrying a single copy of the Delta F508 CFTR mutation (heterozygote resistance). This hypothesis is supported by the strong association between the geography of the diffusion of cattle pastoralism (assessed indirectly by the lactase persistence distribution), the geographic distribution of a sizeable number of HLA alleles (as indicative of potential selective pressures generated by epidemic mortality) and the geographic distribution of the most common mutation causing cystic fibrosis (Delta F508). The systematic interaction of humans with infectious pathogens would have begun in northern Europe, among the carriers of the Funnel Beaker Culture, the first farmers of the North European plain, moving progressively to the south with the dissemination of the cattle pastoralism. This gradual exposure to epidemic mortality among populations located further and further south in Europe as cattle pastoralism expanded could have generated differences in CFTR gene frequencies, thereby shaping the latitudinal frequency gradients observed in present-day European populations.

Balancing selection and heterogeneity across the classical human leukocyte antigen loci: a meta-analytic review of 497 population studies

This paper presents a meta-analysis of high-resolution human leukocyte antigen (HLA) allele frequency data describing 497 population samples. Most of the datasets were compiled from studies published in eight journals from 1990 to 2007; additional datasets came from the International Histocompatibility Workshops and from the AlleleFrequencies.net database. In all, these data represent approximately 66,800 individuals from throughout the world, providing an opportunity to observe trends that may not have been evident at the time the data were originally analyzed, especially with regard to the relative importance of balancing selection among the HLA loci. Population genetic measures of allele frequency distributions were summarized across populations by locus and geographic region. A role for balancing selection maintaining much of HLA variation was confirmed. Further, the breadth of this meta-analysis allowed the ranking of the HLA loci, with DQA1 and HLA-C showing the strongest balancing selection and DPB1 being compatible with neutrality. Comparisons of the allelic spectra reported by studies since 1990 indicate that most of the HLA alleles identified since 2000 are very-low-frequency alleles. The literature-based allele-count data, as well as maps summarizing the geographic distributions for each allele, are available online.

Mitochondrial DNA haplogroup diversity in Basques: a reassessment based on HVI and HVII polymorphisms

This study provides a more complete characterization of the mitochondrial genome variability of the Basques, including data on the hypervariable segment HVII of the D-loop region, which remains relatively unknown. To that end, genomic DNA from 55 healthy men living in the Arratia Valley (Biscay province) and the Goiherri region (Guipúzcoa province) was examined by direct sequencing. Three-generation pedigree charts were compiled to ensure the collection from autochthonous individuals. The most notable findings emerging from the analysis of haplogroup composition are: (i) lack of U8a mitochondrial lineage, a rare subhaplogroup recently identified in Basques and proposed as a Paleolithic marker, (ii) low frequency of haplogroup V, which conflicts with results of earlier analyses describing high frequencies in southwestern Europe, and (iii) high frequency of haplogroup J, especially subhaplogroups J1c1 and J2a. The frequency of haplogroup J does not coincide with previous mtDNA studies in present-day Basques, but is congruent with frequencies found in prehistoric and historic Basque populations. In explaining divergence in haplogroup composition between modern Basque samples, we hypothesized spatial heterogeneity promoted by population fragmentation due to extreme limitation of dispersal opportunities during the Pleistocene glaciations. Similarities between extinct and extant Basque populations as for the high frequency of lineage J, as well as the abundance of this haplogroup in northern Spain endorse a shift in the focus of attention of mtDNA analysts. A refined dissection of haplogroup J might provide more solid evidence about the process of postglacial recolonization of Europe, and thus about the shaping of the European gene pool.

Human genetic variation studies and HLA class II loci

Human leucocyte antigen (HLA) loci widely known for their role in generation of immune responses by encoding cell-surface heterodimers are often considered to be effective for the purpose of reconstructing human phylogenies due to high degree of polymorphism and less recombination. In the present study, we have made an attempt to study HLA class II loci (DRB1, DQA1 and DQB1) in inferring phylogenetic relationship based on both phylogenetic and haplotype approach. In the phylogenetic approach, the compiled database of 19 populations got segregated and finely resolved in three basal clusters with very high bootstrap values corresponding to four geo-ethnic groups of Africans, Orientals, Americans and Caucasians. Maximum- likelihood phylogram has placed North Indian Hindus alongside other Caucasian populations. Haplotype analysis revealed high range of haplotype diversity with nearly 144 observed haplotypes. The haplotype distribution suggested that numbers of Caucasian-specific haplotypes are frequently found among north Indian Hindus. Our results indicate that if the property of less recombination is explored to assign extended haplotypes, followed by strong interpretation based on more logistic statistical model, then there is a high possibility that HLA class II loci can infer exact and accurate phylogenetic assessments as revealed by mtDNA and Y-chromosome markers.

No association between ADRA2A polymorphisms and schizophrenia

There is evidence to suggest that the alpha(2A)-adrenergic receptor may be involved in schizophrenia. With attention directed at the upstream regulatory region of the gene which codes for this receptor (ADRA2A), we proposed that single nucleotide polymorphisms (SNPs) within this region influences susceptibility to schizophrenia by altering the expression of this receptor. We opted to test for an influence on susceptibility by association study using 112 schizophrenic/schizoaffective disorder patients and 159 controls. The region of interest was screened for SNPs using a combination of bioinformatic searches and sequencing. A total of nine SNPs were discovered, of which four (-5972-G/A, -2211-A/T, -1291-C/G and -261-G/A) were genotyped in the entire clinical sample. No associations were evident, suggesting no influence for these SNPs in susceptibility to schizophrenia.

Polymorphic Alu insertions and the genetic structure of Iberian Basques

Eight Alu sequences (ACE, TPA25, PV92, APO, FXIIIB, D1, A25 and B65) were analyzed in two samples from Navarre and Guipúzcoa provinces (Basque Country, Spain). Alu data for other European, Caucasus and North African populations were compiled from the literature for comparison purposes to assess the genetic relationships of the Basques in a broader geographic context. Results of both MDS plot and AMOVA revealed spatial heterogeneity among these three population clusters clearly defined by geography. On the contrary, no substantial genetic heterogeneity was found between the Basque samples, or between Basques and other Europeans (excluding Caucasus populations). Moreover, the genetic information obtained from Alu data conflicts with hypotheses linking the origin of Basques with populations from North Africa (Berbers) or from the Caucasus region (Georgia). In order to explain the reduced genetic heterogeneity detected by Alu insertions among Basque subpopulations, values of the Wright's F(ST )statistic were estimated for both Alu markers and a set of short tandem repeats (STRs) in terms of two geographical scales: (1) the Basque Country, (2) Europe (including Basques). In the Basque area, estimates of Wahlund's effect for both genetic markers showed no statistical difference between Basque subpopulations. However, when this analysis was performed on a European scale, F(ST) values were significantly higher for Alu insertions than for STR alleles. From these results, we suggest that the spatial heterogeneity of the Basque gene pool identified in previous polymorphism studies is relatively recent and probably caused by a differential process of genetic admixture with non-Basque neighboring populations modulated by the effect of a linguistic barrier to random mating.

Gene flow in the Iberian Peninsula determined from Y-chromosome STR loci

In this work, seven multiallelic short-tandem repeat (STR) loci from the nonrecombining region of the human Y-chromosome (DYS19, DYS389 I, DYS389 II, DYS390, DYS391, DYS392, and DYS393) were typed in a sample of residents in the Basque Country (RBAS). In all, 40 different Y-STR haplotypic combinations were identified, resulting in a value of haplotypic diversity of 0.979. Y-STR data compiled from previous works were used for studying Y-chromosome diversity in the Iberian Peninsula and for assessing the effects of migratory movements on the genetic background of the population living currently in territories traditionally occupied by native (autochthonous) Basques. An analysis of the spatial distribution of allelic frequencies of the Y-STRs revealed a geographic pattern characterized by variation gradients (frequency clines) oriented for the most part in the direction southwest-northeast. Accordingly, a neighbor-joining analysis showed a relative polarization between populations located in the northeast and center of the Iberian Peninsula, and the rest of the samples considered. The study sample (RBAS) occupied an intermediate position in the population tree between the autochthonous Basques (BASQ) and the remaining samples. Interestingly, the RBAS collection only showed genetic heterogeneity with that of native Basques (PhiST = 0.013, P < 0.05). Estimates of admixture proportions in the gene pool of RBAS indicated a high level of hybridization with Basque (56%) and non-Basque (44%) genes, which could explain the genetic differentiation observed between BASQ and RBAS.

Sequence-based typing of the complete coding sequence of DQA1 and phenotype frequencies in the Dutch Caucasian population

Typing of DQA1 by sequencing has been a challenge because of a 3-nucleotide deletion in exon 2 in half of the alleles. Furthermore, 19 of the 28 alleles cannot be identified on basis of exon 2 alone, but need additional exon information. With the sequencing strategy presented here the complete exons 1-4 are sequenced heterozygously, enabling identification of all DQA1 alleles by sequence-based typing (SBT). Exons 1-4 were amplified and sequenced separately, the combined sequences were used for automated allele assignment. The method was validated by typing 21 individuals with all possible different allele group combinations. In addition 26 quality control samples were correctly typed by this method. To determine the phenotype frequencies 155 unrelated Dutch Caucasian individuals were DQA1 typed. In total 15 known and two new DQA1 alleles were identified. DQA1*0103 and *0505 were the most frequent alleles with phenotype frequencies of 30% and 29%, respectively. The SBT method presented here is an improvement compared to already existing protocols in that the complete exon sequence is obtained for all coding exons, using identical polymerase chain reaction conditions. Furthermore, all exons are sequenced heterozygously, facilitating allele assignment and reducing the number of amplification reactions.

Genetic position of Valencia (Spain) in the Mediterranean basin according toAlu insertions

In this work, eight human-specific Alu insertion polymorphisms (ACE, TPA25, PV92, APO, FXIIIB, D1, A25, and B65) were typed in 106 unrelated healthy individuals born in the province of Valencia (Spain), with the aim of analyzing the genetic relationships between this region of the Iberian Peninsula and other Mediterranean populations. To that end, Alu data on Eastern European, Western European, and North African populations were compiled from previous studies. The genetic information was stressed by means of genetic distances (R matrix method), nonmetric multidimensional scaling (MDS) and analyses of molecular variance (AMOVA). In Valencia, the most common Alu insertion was APO (0.940), and the least frequent was A25 (0.104). The average gene diversity (GD) computed for the sample examined was comparatively high (0.382). The insertion frequencies estimated for the eight Alu markers were very similar to the mean frequencies calculated for the whole set of populations included in the study, suggesting the hybrid nature of the Valencia's gene pool. MDS and AMOVA results generated from Alu data reveal that the Mediterranean has acted as a strong genetic boundary between the north (Europe) and the south (Northern Africa), resulting in significant gene diversity between the populations of the two regions. Restricted exclusively to the European scope, we suggest the possibility that the Mediterranean could have also acted as a migratory passageway, propitiating the dissemination of cultures and genes between the east and west of Europe and giving rise to some homogenization of gene frequencies among coastal dwelling populations.

Allelic frequencies of 13 STR loci in autochthonous Basques from the province of Vizcaya (Spain)

Allelic frequencies of 13 STR loci (D3S1358, VWA, FGA, D8S1179, D21S11, D18S51, D5S818, D13S317, D16S539, TH01, TPOX, CSF1PO, and D7S820) were estimated from a sample of 73 unrelated healthy donors natives of the Spanish Basque province of Vizcaya. These STR loci constitute the core of polymerase chain reaction (PCR)-based DNA genetic markers in the US Combined DNA Index System (CODIS). All STR loci analysed met Hardy-Weinberg expectations. Based upon the allelic frequencies, forensically important parameters including gene diversity (GD), polymorphism information content (PIC) and power of discrimination (PD) were calculated.

Microsatellite data support subpopulation structuring among Basques

Genomic diversity based on 13 short tandem repeat (STR) loci (D3S1358, vWA, FGA, D8S1179, D21S11, D18S51, D5S818, D13S317, D7S820, D16S539, TH01, TPOX, and CSF1PO) is reported for the first time in Basques from the provinces of Guipúzcoa and Navarre (Spain). STR data from previous studies on Basques from Alava and Vizcaya provinces were also examined using hierarchal analysis of molecular variance (AMOVA) and genetic admixture estimations to ascertain whether the Basques are genetically heterogeneous. To assess the genetic position of Basques in a broader geographic context, we conducted phylogenetic analyses based on F(ST) genetic distances [neighbor-joining trees and multidimensional scaling (MDS)] using data compiled in previous publications. The genetic profile of the Basque groups revealed distinctive regional partitioning of short tandem repeat (STR) diversity. Consistent with the above, native Basques clearly segregated from other populations from Europe (including Spain), North Africa, and the Middle East. The main line of genetic discontinuity inferred from the spatial variability of the microsatellite diversity in Basques significantly overlapped the geographic distribution of the Basque language. The genetic heterogeneity among native Basque groups correlates with the peculiar geography of peopling and marital structure in rural Basque zones and with language boundaries resulting from the uneven impact of Romance languages in the different Basque territories.

Association between dopamine D3 receptor gene polymorphisms and schizophrenia in an isolate population

There are several lines of evidence implicating the dopamine D3 receptor in the pathophysiology of schizophrenia. The Ser9Gly polymorphism of the dopamine D3 receptor gene (DRD3) has been the most extensively investigated DRD3 variant in connection with the disease but results have been inconclusive. Recent reports indicate that the Ser9Gly polymorphism is in linkage disequilibrium with other markers, but association studies between DRD3 haplotypes and schizophrenia have had mixed results. Genetic heterogeneity may be one of the causes of contradicting results. In order to clarify the role of DRD3 alterations in the aetiology of disease, we have investigated three D3 genetic variants (Ser9Gly, -205-G/A, -7685-G/C) in a sample of patients with schizophrenia or schizoaffective disorder (N=118) and controls (N=162) recruited from a human isolate from Navarra (Northern Spain) of Basque origin. Although no association was found between the Ser9Gly or the -205-A/G polymorphisms and disease, an excess of allele -7685-C was observed in patients (p=0.002 after correction for multiple analyses). Haplotype analysis shows the three markers to be in strong linkage disequilibrium (p<0.0001) and strongly associated with disease (p<1x 10(-5)). These results may suggest that these polymorphisms exert a combined or synergistic effect on susceptibility to schizophrenia, or are in linkage with an unknown causative factor. However, further replication in independent samples is required.

Genetic polymorphism and linkage disequilibrium of the HLA-DP region in Basques from Navarre (Spain)

In this work, a sample of 116 individuals from an autochthonous Basque population (northern Navarre, Spain) was typed at the DNA level for HLA-DPA1 and HLA-DPB1 loci, with the aim of analysing the genetic polymorphism and the linkage disequilibrium (LD) of the HLA-DP region. In this Basque subpopulation, the most frequent alleles were *0103 (0.767) and *0201 (0.185) for DPA1 locus, whereas for DPB1 locus the *0401 allele was predominant (0.307). Accordingly, the most frequent haplotype was DPA1*0103-DPB1*0401 (0.300), which showed a significant LD. However, the haplotypes that most differentiated the sample of Navarre from other worldwide populations already analysed were DPA1*0105-DPB1*1901 (0.011) and DPA1*0201-DPB1*7601 (0.021), both of which showed a strong LD. Analysis of the relationships between populations based on data provided by HLA-DPA1 and HLA-DPB1 loci revealed a high genetic affinity between the Basque samples (North Navarre and Guipúzcoa), which, in turn, tend to plot separately from the remaining European populations. Gene frequency clines for DPB1*01, DPB1*04 and DPB1*11 alleles among European populations are reported for the first time. These alleles showed maximum values of F(ST) (0.033, 0.034 and 0.025, respectively). Various evolutionary forces were considered in discussing the origin of the spatial structuring of the gene frequencies: (i) gene flow, argued from the hypotheses of Post-glacial recolonization from southern Europe or the demic diffusion of farmers from the Near East into Europe, and (ii) the existence of selective pressures that could have generated genetic microdifferentiation.