Search

Objective: Primary antiphospholipid syndrome (PAPS) is a rare autoimmune disease characterized by the presence of antiphospholipid antibodies and the occurrence of thrombotic events and pregnancy complications. Our study aimed to identify novel genetic susceptibility loci associated with PAPS., Methods: We performed a genome-wide association study comprising 5,485 individuals (482 affected individuals) of European ancestry. Significant and suggestive independent variants from a meta-analysis of approximately 7 million variants were evaluated for functional and biological process enrichment. The genetic risk variability for PAPS in different populations was also assessed. Hierarchical clustering, Mahalanobis distance, and Dirichlet Process Mixtures with uncertainty clustering methods were used to assess genetic similarities between PAPS and other immune-mediated diseases., Results: We revealed genetic associations with PAPS in a regulatory locus within the HLA class II region near HLA-DRA and in STAT1-STAT4 with a genome-wide level of significance; 34 additional suggestive genetic susceptibility loci for PAPS were also identified. The disease risk allele near HLA-DRA is associated with overexpression of HLA-DRB6, HLA-DRB9, HLA-DQA2, and HLA-DQB2 in immune cells, vascular tissue, and nervous tissue. This association is independent of the association between PAPS and HLA-DRB1*1302. Functional analyses highlighted immune-related pathways in PAPS-associated loci. The comparison with other immune-mediated diseases revealed a close genetic relatedness to neuromyelitis optica, systemic sclerosis, and Sjögren syndrome, suggesting co-localized causal variations close to STAT1-STAT4, TNPO3, and BLK., Conclusion: This study represents a comprehensive large-scale genetic analysis for PAPS and provides new insights into the genetic basis and pathophysiology of this rare disease., (© 2024 The Author(s). Arthritis & Rheumatology published by Wiley Periodicals LLC on behalf of American College of Rheumatology.)

Genome-wide association studies (GWASs) have been successful at finding associations between genetic variants and human traits, including the immune-mediated diseases (IMDs). However, the requirement of large sample sizes for discovery poses a challenge for learning about less common diseases, where increasing volunteer numbers might not be feasible. An example of this is myositis (or idiopathic inflammatory myopathies [IIM]s), a group of rare, heterogeneous autoimmune diseases affecting skeletal muscle and other organs, severely impairing life quality. Here, we applied a feature engineering method to borrow information from larger IMD GWASs to find new genetic associations with IIM and its subgroups. Combining this approach with two clustering methods, we found 17 IMDs genetically close to IIM, including some common comorbid conditions, such as systemic sclerosis and Sjögren's syndrome, as well as hypo- and hyperthyroidism. All IIM subtypes were genetically similar within this framework. Next, we colocalized IIM signals that overlapped IMD signals, and found seven potentially novel myositis associations mapped to immune-related genes, including BLK, IRF5/TNPO3, and ITK/HAVCR2, implicating a role for both B and T cells in IIM. This work proposes a new paradigm of genetic discovery in rarer diseases by leveraging information from more common IMD, and can be expanded to other conditions and traits beyond IMD., Competing Interests: Declaration of interests Dr. Wallace receives research funding from GSK and MSD and is a part-time employee of GSK. Neither company had any influence on this work or its publication. Dr. Radstake is an employee of Abbvie and may hold stock. Abbvie had no influence on the content of this work or its publication. Dr. Chinoy has received fees as a speaker for GSK and UCB; consulting for PTC Therapeutics; advisory board member for Astra Zeneca, Pfizer, Argenx, and Galapagos; and as a data and science monitoring board chair for Horizon Therapeutics. Dr. Maurer has grants from Novartis, consulting fees from Novartis, Boehringer Ingelheim, Jannsen-Cilag, and GSK; speaker fees from Boehringer-Ingelheim, GSK, Novartis, Otsuka, and MSD; congress support from Medtalk, Pfizer, Roche, Actelion, Mepha, and MSD; and has a patent mir-29 for the treatment of SSc (US8247389, EP2331143). Dr. Wedderburn has received speaker and consultancy fees from Pfizer paid to UCL, unrelated to this work., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Phenome-wide association studies (PheWAS) facilitate the discovery of associations between a single genetic variant with multiple phenotypes. For variants which impact a specific protein, this can help identify additional therapeutic indications or on-target side effects of intervening on that protein. However, PheWAS is restricted by an inability to distinguish confounding due to linkage disequilibrium (LD) from true pleiotropy. Here we describe CoPheScan (Coloc adapted Phenome-wide Scan), a Bayesian approach that enables an intuitive and systematic exploration of causal associations while simultaneously addressing LD confounding. We demonstrate its performance through simulation, showing considerably better control of false positive rates than a conventional approach not accounting for LD. We used CoPheScan to perform PheWAS of protein-truncating variants and fine-mapped variants from disease and pQTL studies, in 2275 disease phenotypes from the UK Biobank. Our results identify the complexity of known pleiotropic genes such as APOE, and suggest a new causal role for TGM3 in skin cancer., (© 2024. The Author(s).)

Objective: To map geographic clusters of rare disorders and congenital anomalies reported in South America., Methods: Qualitative systematic review conducted in Medline/PubMed, Lilacs, and Scielo electronic databases to identify studies meeting eligibility criteria. The strategy resulted in 1 672 unique articles, from which 164 were selected for full reading by a pair of reviewers., Results: Fifty-five articles reported at least one cluster of genetic disorders or congenital anomalies in South American territory. From these papers, 122 clusters were identified, of which half (61) were related to autosomal recessive disorders. Sixty-five (53.3%) of the clusters were located in Brazil., Conclusions: The results of the review reinforce that rare diseases and congenital anomalies can occur in a non-random way in space, which is discussed in the perspective of the complex history of formation, social organization, and genetic structure of the South American population. Mapping clusters in population medical genetics can be an important public health tool, given that such places concentrate cases of rare diseases that frequently require multiprofessional, specialized care. Therefore, these results can support important agendas in public health related to rare diseases and congenital anomalies, such as health promotion and surveillance.

Human Leukocyte Antigen (HLA)-G participates in several biological processes, including reproduction, vascular remodeling, immune tolerance, and hypoxia response. HLA-G is a potential candidate gene for high altitude adaptation since its expression is modulated in both micro and macro environment under hypoxia and constant cellular stress. Besides the promoter region, the HLA-G 3'untranslated region (UTR) influences HLA-G expression patterns through several post-transcriptional mechanisms. Currently, the 3'UTR genetic diversity in terms of altitude adaptation of Native American populations is still unexplored, particularly at high altitude ecoregions. Here, we evaluated 288 Native Americans from 9 communities located in the Andes [highland (HL); ≥2,500 m (range = 2,838-4,433 m)] and 8 populations located in lowland (LL) regions [<2,500 m (range = 80-431 m); Amazonian tropical forest, Brazilian central plateau, and Chaco] of South America. In total, nine polymorphic sites and ten haplotypes were observed. The most frequent haplotypes (UTR-1, UTR-2, and UTR-3) accounted for ∼ 77% of haplotypes found in LL, while in the HL, the same haplotypes reach ∼ 93%. Also, a remarkable high frequency of putative ancestral UTR-5 haplotype was observed in LL (21.5%), while in HL UTR-2 reaches up to 47%. Further, UTR-2 frequency positively correlates with altitude-related variables, while a negative correlation for UTR-5 was observed. From an evolutionary perspective, we observed a tendency towards balancing selection in HL and LL populations thus suggesting that haplotypes of ancient and more derived alleles may have been co-opted for relatively recent adaptations such as those experienced by modern humans in the highland and lowland of South America. We also discuss how long-term balancing selection can be a reservoir of genetic variants that can be positively selected. Finally, our study provides some pieces of evidence that HLA-G 3'UTR haplotypes may have contributed to high altitude adaptation in the Andes., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.)

Motivation: Polygenic scores (PGS) aim to genetically predict complex traits at an individual level. PGS are typically trained on genome-wide association summary statistics and require an independent test dataset to tune parameters. More recent methods allow parameters to be tuned on the training data, removing the need for independent test data, but approaches are computationally intensive. Based on fine-mapping principles, we present RápidoPGS, a flexible and fast method to compute PGS requiring summary-level Genome-wide association studies (GWAS) datasets only, with little computational requirements and no test data required for parameter tuning., Results: We show that RápidoPGS performs slightly less well than two out of three other widely used PGS methods (LDpred2, PRScs and SBayesR) for case-control datasets, with median r2 difference: -0.0092, -0.0042 and 0.0064, respectively, but up to 17 000-fold faster with reduced computational requirements. RápidoPGS is implemented in R and can work with user-supplied summary statistics or download them from the GWAS catalog., Availability and Implementation: Our method is available with a GPL license as an R package from CRAN and GitHub., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2021. Published by Oxford University Press.)

STAT2 plays a strategic role in defending viral infection through the signaling cascade involving the immune system initiated after type I interferon release. Many flaviviruses target the inactivation or degradation of STAT2 as a strategy to impair this host's line of defense. Primates are natural reservoirs for a range of disease-causing flaviviruses (e.g., Zika, Dengue, and Yellow Fever virus), while rodents appear less susceptible. We analyzed the STAT2 coding sequence of 28 Rodentia species and 49 Primates species. Original data from 19 Platyrrhini species were sequenced for the SH2 domain of STAT2 and included in the analysis. STAT2 has many sites whose variation can be explained by positive selection, measurement by two methods (PALM indicated 12, MEME 61). Both evolutionary tests significantly marked sites 127, 731, 739, 766, and 780. SH2 is under evolutionary constraint but presents episodic positive selection events within Rodentia: in one of them, a moderately radical change (serine > arginine) at position 638 is found in Peromyscus species, and can be implicated in the difference in susceptibility to flaviviruses within Rodentia. Some other positively selected sites are functional such as 5, 95, 203, 251, 782, and 829. Sites 251 and 287 regulate the signaling mediated by the JAK-STAT2 pathway, while 782 and 829 create a stable tertiary structure of STAT2, facilitating its connection with transcriptional co-activators. Only three positively selected sites, 5, 95, and 203, are recognized members who act on the interface between STAT2 and flaviviruses NS5 protein. We suggested that due to the higher evolutionary rate, rodents are, at this moment, taking some advantage in the battle against infections for some well-known Flaviviridae, in particular when compared to primates. Our results point to dynamics that fit with a molecular evolutionary scenario shaped by a thought-provoking virus-host arms race., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Genome-wide association studies (GWAS) have identified thousands of genetic variants that are associated with complex traits. However, a stringent significance threshold is required to identify robust genetic associations. Leveraging relevant auxiliary covariates has the potential to boost statistical power to exceed the significance threshold. Particularly, abundant pleiotropy and the non-random distribution of SNPs across various functional categories suggests that leveraging GWAS test statistics from related traits and/or functional genomic data may boost GWAS discovery. While type 1 error rate control has become standard in GWAS, control of the false discovery rate can be a more powerful approach. The conditional false discovery rate (cFDR) extends the standard FDR framework by conditioning on auxiliary data to call significant associations, but current implementations are restricted to auxiliary data satisfying specific parametric distributions, typically GWAS p-values for related traits. We relax these distributional assumptions, enabling an extension of the cFDR framework that supports auxiliary covariates from arbitrary continuous distributions ("Flexible cFDR"). Our method can be applied iteratively, thereby supporting multi-dimensional covariate data. Through simulations we show that Flexible cFDR increases sensitivity whilst controlling FDR after one or several iterations. We further demonstrate its practical potential through application to an asthma GWAS, leveraging various functional genomic data to find additional genetic associations for asthma, which we validate in the larger, independent, UK Biobank data resource., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: Anna Hutchinson is a GSK-sponsored iCASE student.

Prolactin (PRL) is a pleiotropic neurohormone secreted by the mammalian pituitary gland into the blood, thus reaching many tissues and organs beyond the brain. PRL binds to its receptor, PRLR, eliciting a molecular signaling cascade. This system modulates essential mammalian behaviors and promotes notable modifications in the reproductive female tissues and organs. Here, we explore how the intracellular domain of PRLR (PRLR-ICD) modulates the expression of the PRLR gene. Despite differences in the reproductive strategies between eutherian and metatherian mammals, there is no clear distinction between PRLR-ICD functional motifs. However, we found selection signatures that showed differences between groups, with many conserved functional elements strongly maintained through purifying selection across the class Mammalia. We observed a few residues under relaxed selection, the levels of which were more pronounced in Eutheria and particularly striking in primates (Simiiformes), which could represent a pre-adaptive genetic element protected from purifying selection. Alternative, new motifs, such as YLDP (318-321) and others with residues Y283 and Y290, may already be functional. These motifs would have been co-opted in primates as part of a complex genetic repertoire related to some derived adaptive phenotypes, but these changes would have no impact on the primordial functions that characterize the mammals as a whole and that are related to the PRL-PRLR system., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Takayasu arteritis is a rare inflammatory disease of large arteries. We performed a genetic study in Takayasu arteritis comprising 6,670 individuals (1,226 affected individuals) from five different populations. We discovered HLA risk factors and four non-HLA susceptibility loci in VPS8, SVEP1, CFL2, and chr13q21 and reinforced IL12B, PTK2B, and chr21q22 as robust susceptibility loci shared across ancestries. Functional analysis proposed plausible underlying disease mechanisms and pinpointed ETS2 as a potential causal gene for chr21q22 association. We also identified >60 candidate loci with suggestive association (p < 5 × 10 -5 ) and devised a genetic risk score for Takayasu arteritis. Takayasu arteritis was compared to hundreds of other traits, revealing the closest genetic relatedness to inflammatory bowel disease. Epigenetic patterns within risk loci suggest roles for monocytes and B cells in Takayasu arteritis. This work enhances understanding of the genetic basis and pathophysiology of Takayasu arteritis and provides clues for potential new therapeutic targets., (Copyright © 2020 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

We carried out an exhaustive review regarding human skin color variation and how much it may be related to vitamin D metabolism and other photosensitive molecules. We discuss evolutionary contexts that modulate this variability and hypotheses postulated to explain them; for example, a small amount of melanin in the skin facilitates vitamin D production, making it advantageous to have fair skin in an environment with little radiation incidence. In contrast, more melanin protects folate from degradation in an environment with a high incidence of radiation. Some Native American populations have a skin color at odds with what would be expected for the amount of radiation in the environment in which they live, a finding challenging the so-called "vitamin D-folate hypothesis." Since food is also a source of vitamin D, dietary habits should also be considered. Here we argue that a gene network approach provides tools to explain this phenomenon since it indicates potential alleles co-evolving in a compensatory way. We identified alleles of the vitamin D metabolism and pigmentation pathways segregated together, but in different proportions, in agriculturalists and hunter-gatherers. Finally, we highlight how an evolutionary approach can be useful to understand current topics of medical interest., (© 2020 Wiley Periodicals LLC.)

Background: Genome-wide association studies (GWAS) have identified pervasive sharing of genetic architectures across multiple immune-mediated diseases (IMD). By learning the genetic basis of IMD risk from common diseases, this sharing can be exploited to enable analysis of less frequent IMD where, due to limited sample size, traditional GWAS techniques are challenging., Methods: Exploiting ideas from Bayesian genetic fine-mapping, we developed a disease-focused shrinkage approach to allow us to distill genetic risk components from GWAS summary statistics for a set of related diseases. We applied this technique to 13 larger GWAS of common IMD, deriving a reduced dimension "basis" that summarised the multidimensional components of genetic risk. We used independent datasets including the UK Biobank to assess the performance of the basis and characterise individual axes. Finally, we projected summary GWAS data for smaller IMD studies, with less than 1000 cases, to assess whether the approach was able to provide additional insights into genetic architecture of less common IMD or IMD subtypes, where cohort collection is challenging., Results: We identified 13 IMD genetic risk components. The projection of independent UK Biobank data demonstrated the IMD specificity and accuracy of the basis even for traits with very limited case-size (e.g. vitiligo, 150 cases). Projection of additional IMD-relevant studies allowed us to add biological interpretation to specific components, e.g. related to raised eosinophil counts in blood and serum concentration of the chemokine CXCL10 (IP-10). On application to 22 rare IMD and IMD subtypes, we were able to not only highlight subtype-discriminating axes (e.g. for juvenile idiopathic arthritis) but also suggest eight novel genetic associations., Conclusions: Requiring only summary-level data, our unsupervised approach allows the genetic architectures across any range of clinically related traits to be characterised in fewer dimensions. This facilitates the analysis of studies with modest sample size by matching shared axes of both genetic and biological risk across a wider disease domain, and provides an evidence base for possible therapeutic repurposing opportunities.

Platyrrhini (New World monkeys, NWm) are a group of primates characterized by behavioral and reproductive traits that are otherwise uncommon among primates, including social monogamy, direct paternal care, and twin births. As a consequence, the study of Platyrrhine primates is an invaluable tool for the discovery of the genetic repertoire underlying these taxon-specific traits. Recently, high conservation of vasopressin (AVP) sequence, in contrast with high variability of oxytocin (OXT), has been described in NWm. AVP and OXT functions are possible due to interaction with their receptors: AVPR1a, AVPR1b, AVPR2, and OXTR; and the variability in this system is associated with the traits mentioned above. Understanding the variability in the receptors is thus fundamental to understand the function and evolution of the system as a whole. Here we describe the variability of AVPR1b coding region in 20 NWm species, which is well-known to influence behavioral traits such as aggression, anxiety, and stress control in placental mammals. Our results indicate that 4% of AVPR1b sites may be under positive selection and a significant number of sites under relaxed selective constraint. Considering the known role of AVPR1b, we suggest that some of the changes described here for the Platyrrhini may be a part of the genetic repertoire connected with the complex network of neuroendocrine mechanisms of AVP-OXT system in the modulation of the HPA axis. Thus, these changes may have promoted the emergence of social behaviors such as direct paternal care in socially monogamous species that are also characterized by small body size and twin births., (© 2019 Wiley Periodicals, Inc.)

Objectives: The major aim of this article was to estimate the demographic impact of European arrival and colonization over Native American populations from southern Brazil and Uruguay. We also compared the mitochondrial DNA (mtDNA) genetic diversity, structure, and demography of Native American lineages present in current indigenous (Natives) and nonindigenous admixed (Admixed) populations to estimate the effective population size (N e ) of contemporary and ancestral (pre-Columbian) Native American populations., Methods: We retrieved published mtDNA sequences from Native (n = 396) and Admixed (n = 309) populations from southern Brazil, Uruguay, and surrounding areas. We conducted genetic diversity, structure, and demographic analyses. Finally, we used Approximate Bayesian Computation to estimate the N e for current Native, Admixed, and pre-Columbian Native American populations., Results: We found higher Native American mtDNA genetic diversity in admixed rather than in indigenous populations (131/309 vs 27/396 different haplotypes, respectively). Only Admixed populations maintained ancient signals of the Native American population expansion approximately 14 to 17 kya, which have decayed in Natives. Our N e estimates suggest that Natives represent only 0.33% (0.18%-1.19%) of the N e for ancestral pre-Columbian indigenous populations., Conclusions: Admixed populations represent an important genetic reservoir of Native American lineages, many of which are extinct in contemporary indigenous populations. In addition, the Native American lineages present in Admixed populations retain part of the past demographic history of Native Americans. The intensity of the reduction is congruent with historical accounts of strong indigenous depopulation during the colonization process., (© 2019 Wiley Periodicals, Inc.)

Traits that undergo massive natural selection pressure, with multiple events of positive selection, are hard to find. Social behaviour, in social animals, is crucial for survival, and genetic networks involved in behaviour, such as those of serotonin (5-HT) and other neurotransmitters, must be the target of natural selection. Here, we used molecular analyses to search for signals of positive selection in the 5-HT system and found such signals in the M3-M4 intracellular domain of the 5-HT3A serotonin receptor subunit (HTR3A) in primates. We detected four amino acid sites with signs of putatively positive selection (398, 403, 432 and 416); the first three showed indications of being selected in New World monkeys (NWM, Platyrrhini), specifically in the Callitrichinae branch. Additionally, we searched for associations of these amino acid variants with social behavioural traits (i.e. sex-biased dispersal, dominance and social monogamy) using classical and Bayesian methods, and found statistically significant associations for unbiased sex dispersal (398L and 416S), unbiased sex dominance (416S) and social monogamy (416S), as well as significant positive correlation between female dispersal and 403G. Furthermore, we found putatively functional protein motifs determined by three selected sites, of which we highlight a ligand motif to GSK3 in the 416S variant, appearing only in Platyrrhini. 5-HT, 5-HT3A receptor and GSK3 are part of a network that participates in neurodevelopment and regulates behaviour, among other functions. We suggest that these genetic variations, together with those found in other neurotransmitter systems, must contribute to adaptive behaviours and consequently to fitness in NWMs., (© 2018 European Society For Evolutionary Biology. Journal of Evolutionary Biology © 2018 European Society For Evolutionary Biology.)

Objectives: To determine genetic differences between agriculturalist and hunter-gatherer southern Native American populations for selected metabolism-related markers and to test whether Neel's thrifty genotype hypothesis (TGH) could explain the genetic patterns observed in these populations., Materials and Methods: 375 Native South American individuals from 17 populations were genotyped using six markers (APOE rs429358 and rs7412; APOA2 rs5082; CD36 rs3211883; TCF7L2 rs11196205; and IGF2BP2 rs11705701). Additionally, APOE genotypes from 39 individuals were obtained from the literature. AMOVA, main effects, and gene-gene interaction tests were performed., Results: We observed differences in allele distribution patterns between agriculturalists and hunter-gatherers for some markers. For instance, between-groups component of genetic variance (F CT ) for APOE rs429358 showed strong differences in allelic distributions between hunter-gatherers and agriculturalists (p = 0.00196). Gene-gene interaction analysis indicated that the APOE E4/CD36 TT and APOE E4/IGF2BP2 A carrier combinations occur at a higher frequency in hunter-gatherers, but this combination is not replicated in archaic (Neanderthal and Denisovan) and ancient (Anzick, Saqqaq, Ust-Ishim, Mal'ta) hunter-gatherer individuals., Discussion: A complex scenario explains the observed frequencies of the tested markers in hunter-gatherers. Different factors, such as pleotropic alleles, rainforest selective pressures, and population dynamics, may be collectively shaping the observed genetic patterns. We conclude that although TGH seems a plausible hypothesis to explain part of the data, other factors may be important in our tested populations., (© 2017 Wiley Periodicals, Inc.)

Here we report a genome-wide association study for non-pathological pinna morphology in over 5,000 Latin Americans. We find genome-wide significant association at seven genomic regions affecting: lobe size and attachment, folding of antihelix, helix rolling, ear protrusion and antitragus size (linear regression P values 2 × 10(-8) to 3 × 10(-14)). Four traits are associated with a functional variant in the Ectodysplasin A receptor (EDAR) gene, a key regulator of embryonic skin appendage development. We confirm expression of Edar in the developing mouse ear and that Edar-deficient mice have an abnormally shaped pinna. Two traits are associated with SNPs in a region overlapping the T-Box Protein 15 (TBX15) gene, a major determinant of mouse skeletal development. Strongest association in this region is observed for SNP rs17023457 located in an evolutionarily conserved binding site for the transcription factor Cartilage paired-class homeoprotein 1 (CART1), and we confirm that rs17023457 alters in vitro binding of CART1.

Andalusia is the most densely populated region of Spain since ancient times, and has a rich history of contacts across the Mediterranean. Earlier studies have underlined the relatively high frequency of the Sub-Saharan GM 1,17 5* haplotype in western Andalusia (Huelva province, n=252) and neighbouring Atlantic regions. Here, we provide novel data on GM/KM markers in eastern Andalusians (n=195) from Granada province, where African GM*1,17 5* frequency is relatively high (0.044). The most frequent GM haplotypes in Andalusia parallel the most common in Europe. Altogether, these data allow us to gain insight into the genetic diversity of southern Iberia. Additionally, we assess population structure by comparing our Iberian samples with 41 Mediterranean populations. GM haplotype variation across the Mediterranean reflects intense and complex interactions between North Africans and South Europeans along human history, highlighting that African influence over the Iberian Peninsula does not follow an isotropic pattern., (Copyright © 2014 Académie des sciences. Published by Elsevier SAS. All rights reserved.)

Background: The APOE gene has received much attention due to the remarkable spatial variation patterns of some of its genotypes and alleles in human populations and to its relevance in biomedicine., Aim: This work was addressed to investigate the extent of APOE polymorphism between autochthonous Andalusians originating from Huelva and Granada provinces. No data on this marker in these southern Spanish coastal populations are available up to date., Subjects and Methods: This study used genomic DNA from healthy, unrelated Andalusians of both sexes (n = 322). All samples were genotyped for two SNPs, rs429358 and rs7412, which determine the three APOE alleles: ε2, ε3 and ε4. For analyses, a TaqMan-based technique was applied using a RT-PCR. Comparisons with other Mediterranean populations were performed based on multivariate analysis., Results: A relatively high frequency of ε4 in Granada (eastern Andalusia), as well as a low ε2 frequency in Huelva (western Andalusia) were observed. The finding that ε4 allele in Southern Spain and Portugal is higher than expected given its geographical location poses an interesting question for this study, given the well-established APOE-ε4 gradient in Europe., Conclusion: This population study may represent useful information for further prospective anthropological and molecular genetic studies focused on unravelling the relationship between population genetic composition and specific human diseases.

Background: The archeology and history of the ancient Mediterranean have shown that this sea has been a permeable obstacle to human migration. Multiple cultural exchanges around the Mediterranean have taken place with presumably population admixtures. A gravitational territory of those migrations has been the Iberian Peninsula. Here we present a comprehensive analysis of the maternal gene pool, by means of control region sequencing and PCR-RFLP typing, of autochthonous Andalusians originating from the coastal provinces of Huelva and Granada, located respectively in the west and the east of the region., Results: The mtDNA haplogroup composition of these two southern Spanish populations has revealed a wide spectrum of haplogroups from different geographical origins. The registered frequencies of Eurasian markers, together with the high incidence and diversification of African maternal lineages (15% of the total mitochondrial variability) among Huelva Andalusians when compared to its eastwards relatives of Granada and other Iberian populations, constitute relevant findings unknown up-to-date on the characteristics of mtDNA within Andalusia that testifies a female population substructure. Therefore, Andalusia must not be considered a single, unique population., Conclusions: The maternal legacy among Andalusians reflects distinctive local histories, pointing out the role of the westernmost territory of Peninsular Spain as a noticeable recipient of multiple and diverse human migrations. The obtained results underline the necessity of further research on genetic relationships in both sides of the western Mediterranean, using carefully collected samples from autochthonous individuals. Many studies have focused on recent North African gene flow towards Iberia, yet scientific attention should be now directed to thoroughly study the introduction of European genes in northwest Africa across the sea, in order to determine its magnitude, timescale and methods, and to compare them to those terrestrial movements from eastern Africa and southwestern Asia.