Genomic Analyses of Human European Diversity at the Southwestern Edge: Isolation, African Influence and Disease Associations in the Canary Islands

The genomic history of the indigenous people of the Canary Islands

The indigenous population of the Canary Islands, which colonized the archipelago around the 3rd century CE, provides both a window into the past of North Africa and a unique model to explore the effects of insularity. We generate genome-wide data from 40 individuals from the seven islands, dated between the 3rd-16rd centuries CE. Along with components already present in Moroccan Neolithic populations, the Canarian natives show signatures related to Bronze Age expansions in Eurasia and trans-Saharan migrations. The lack of gene flow between islands and constant or decreasing effective population sizes suggest that populations were isolated. While some island populations maintained relatively high genetic diversity, with the only detected bottleneck coinciding with the colonization time, other islands with fewer natural resources show the effects of insularity and isolation. Finally, consistent genetic differentiation between eastern and western islands points to a more complex colonization process than previously thought.

A draft conceptual model of SLC6A1 neurodevelopmental disorder

Introduction

SLC6A1 Neurodevelopmental Disorder (SLC6A1-NDD), first described in 2015, is a rare syndrome caused by a mutation in the SLC6A1 gene which encodes for the GABA Transporter 1 (GAT-1) protein. Epilepsy is one of the most common symptoms in patients and is often the primary treatment target, though the severity of epilepsy is variable. The impact of seizures and other symptoms of SLC6A1-NDD on patients and caregivers is wide-ranging and has not been described in a formal disease concept study.

Methods

A literature search was performed using the simple search term, "SLC6A1." Papers published before 2015, and those which did not describe the human neurodevelopmental disorder were removed from analysis. Open-ended interviews on lived experiences were conducted with two patient advocate key opinion leaders. An analysis of de-identified conversations between families of people with SLC6A1-NDD on social media was performed to quantify topics of concern.

Results

Published literature described symptoms in all of the following domains: neurological, visual, motor, cognitive, communication, behavior, gastrointestinal, sleep, musculo-skeletal, and emotional in addition to epilepsy. Key opinion leaders noted two unpublished features: altered hand use in infants, and developmental regression with onset of epilepsy. Analysis of social media interactions confirmed that the core symptoms of epilepsy and autistic traits were prominent concerns, but also demonstrated that other symptoms have a large impact on family life.

Discussion

For rare diseases, analysis of published literature is important, but may not be as comprehensive as that which can be gleaned from spontaneous interactions between families and through qualitative interviews. This report reflects our current understanding of the lived experience of SLC6A1-NDD. The discrepancy between the domains of disease reported in the literature and those discussed in patient conversations suggests that a formal qualitative interview-based disease concept study of SLC6A1-NDD is warranted.

Digging into the admixture strata of current-day Canary Islanders based on mitogenomes

The conquest of the Canary Islands by Europeans began at the beginning of the 15th century and culminated in 1496 with the surrender of the aborigines. The collapse of the aboriginal population during the conquest and the arrival of settlers caused a drastic change in the demographic composition of the archipelago. To shed light on this historical process, we analyzed 896 mitogenomes of current inhabitants from the seven main islands. Our findings confirm the continuity of aboriginal maternal contributions and the persistence of their genetic footprints in the current population, even at higher levels (>60% on average) than previously evidenced. Moreover, the age estimates for most autochthonous founder lineages support a first aboriginal arrival to the islands at the beginning of the first millennium. We also revealed for the first time that the main recognizable genetic influences from Europe are from Portuguese and Galicians.

Developing CIRdb as a catalog of natural genetic variation in the Canary Islanders

The current inhabitants of the Canary Islands have a unique genetic makeup in the European diversity landscape due to the existence of African footprints from recent admixture events, especially of North African components (> 20%). The underrepresentation of non-Europeans in genetic studies and the sizable North African ancestry, which is nearly absent from all existing catalogs of worldwide genetic diversity, justify the need to develop CIRdb, a population-specific reference catalog of natural genetic variation in the Canary Islanders. Based on array genotyping of the selected unrelated donors and comparisons against available datasets from European, sub-Saharan, and North African populations, we illustrate the intermediate genetic differentiation of Canary Islanders between Europeans and North Africans and the existence of within-population differences that are likely driven by genetic isolation. Here we describe the overall design and the methods that are being implemented to further develop CIRdb. This resource will help to strengthen the implementation of Precision Medicine in this population by contributing to increase the diversity in genetic studies. Among others, this will translate into improved ability to fine map disease genes and simplify the identification of causal variants and estimate the prevalence of unattended Mendelian diseases.

A catalog of the genetic causes of hereditary angioedema in the Canary Islands (Spain)

Hereditary angioedema (HAE) is a rare disease where known causes involve C1 inhibitor dysfunction or dysregulation of the kinin cascade. The updated HAE management guidelines recommend performing genetic tests to reach a precise diagnosis. Unfortunately, genetic tests are still uncommon in the diagnosis routine. Here, we characterized for the first time the genetic causes of HAE in affected families from the Canary Islands (Spain). Whole-exome sequencing data was obtained from 41 affected patients and unaffected relatives from 29 unrelated families identified in the archipelago. The Hereditary Angioedema Database Annotation (HADA) tool was used for pathogenicity classification and causal variant prioritization among the genes known to cause HAE. Manual reclassification of prioritized variants was used in those families lacking known causal variants. We detected a total of eight different variants causing HAE in this patient series, affecting essentially SERPING1 and F12 genes, one of them being a novel SERPING1 variant (c.686-12A>G) with a predicted splicing effect which was reclassified as likely pathogenic in one family. Altogether, the diagnostic yield by assessing previously reported causal genes and considering variant reclassifications according to the American College of Medical Genetics guidelines reached 66.7% (95% Confidence Interval [CI]: 30.1-91.0) in families with more than one affected member and 10.0% (95% CI: 1.8-33.1) among cases without family information for the disease. Despite the genetic causes of many patients remain to be identified, our results reinforce the need of genetic tests as first-tier diagnostic tool in this disease, as recommended by the international WAO/EAACI guidelines for the management of HAE.

Admixture Mapping of Sepsis in European Individuals With African Ancestries

Sepsis is a severe systemic inflammatory response to infections that is accompanied by organ dysfunction. Although the ancestral genetic background is a relevant factor for sepsis susceptibility, there is a lack of studies using the genetic singularities of a recently admixed population to identify loci involved in sepsis susceptibility. Here we aimed to discover new sepsis loci by completing the first admixture mapping study of sepsis in Canary Islanders, leveraging their distinctive genetic makeup as a mixture of Europeans and African ancestries. We used a case-control approach and inferred local ancestry blocks from genome-wide data from 113,414 polymorphisms genotyped in 343 patients with sepsis and 410 unrelated controls, all ascertained for grandparental origin in the Canary Islands (Spain). Deviations in local ancestries between cases and controls were tested using logistic regressions, followed by fine-mapping analyses based on imputed genotypes, in silico functional assessments, and gene expression analysis centered on the region of interest. The admixture mapping analysis detected that local European ancestry in a locus spanning 1.2 megabases of chromosome 8p23.1 was associated with sepsis (lowest p = 1.37 × 10⁻⁴; Odds Ratio [OR] = 0.51; 95%CI = 0.40–0.66). Fine-mapping studies prioritized the variant rs13249564 within intron 1 of MFHAS1 gene associated with sepsis (p = 9.94 × 10⁻⁴; OR = 0.65; 95%CI = 0.50–0.84). Functional and gene expression analyses focused on 8p23.1 allowed us to identify alternative genes with possible biological plausibility such as defensins, which are well-known effector molecules of innate immunity. By completing the first admixture mapping study of sepsis, our results revealed a new genetic locus (8p23.1) harboring a number of genes with plausible implications in sepsis susceptibility.

Role of the SARS-CoV-2 virus in the appearance of new onset type 1 diabetes mellitus in children in Gran Canaria, Spain

OBJECTIVES

It has been hypothesized that SARS-CoV-2 may play a role in the development of different forms of diabetes mellitus (DM). The Canary Islands have the highest incidence of type 1 DM (T1DM) reported in Spain (30-35/100,000 children under 14 years/year). In 2020-2021 we observed the highest incidence so far on the island of Gran Canaria, as a result of which we decided to evaluate the possible role of COVID-19 in the increased number of onsets.

METHODS

We examined the presence of IgG antibodies against SARS-CoV-2 in children with new onset T1DM between October 2020 and August 2021. We compared recent T1DM incidence with that of the previous 10 years.

RESULTS

Forty-two patients were diagnosed with T1DM (48.1/100,000 patients/year), representing a nonsignificant 25.7% increase from the expected incidence. Of the 33 patients who consented to the study, 32 presented negative IgG values, with only one patient reflecting undiagnosed past infection. Forty-four percent of patients presented with ketoacidosis at onset, which was similar to previous years.

CONCLUSIONS

We conclude that there is no direct relationship between the increased incidence of T1DM and SARS-CoV-2 in the region. The COVID-19 pandemic did not result in an increased severity of T1DM presentation.

Targeted analysis of genomic regions enriched in African ancestry reveals novel classical HLA alleles associated with asthma in Southwestern Europeans

Despite asthma has a considerable genetic component, an important proportion of genetic risks remain unknown, especially for non-European populations. Canary Islanders have the largest African genetic ancestry observed among Southwestern Europeans and the highest asthma prevalence in Spain. Here we examined broad chromosomal regions previously associated with an excess of African genetic ancestry in Canary Islanders, with the aim of identifying novel risk variants associated with asthma susceptibility. In a two-stage cases-control study, we revealed a variant within HLA-DQB1 significantly associated with asthma risk (rs1049213, meta-analysis p = 1.30 × 10^–7, OR [95% CI] = 1.74 [1.41–2.13]) previously associated with asthma and broad allergic phenotype. Subsequent fine-mapping analyses of classical HLA alleles revealed a novel allele significantly associated with asthma protection (HLA-DQA1*01:02, meta-analysis p = 3.98 × 10^–4, OR [95% CI] = 0.64 [0.50–0.82]) that had been linked to infectious and autoimmune diseases, and peanut allergy. HLA haplotype analyses revealed a novel haplotype DQA1*01:02-DQB1*06:04 conferring asthma protection (meta-analysis p = 4.71 × 10^–4, OR [95% CI] = 0.47 [0.29– 0.73]).

First Census of Patients with Hereditary Angioedema in the Canary Islands

Hereditary angioedema (HAE) is a rare genetic condition whose main symptoms are recurrent swelling in the skin, mucosa, and internal organs. Recent studies suggested that the regulation of the inflammatory response and the complement cascade are two of the pathways significantly enriched in the Canary Islands, Spain. Here, we describe the first HAE patient series in this region. Forty-one patients (33 F, 8 M) and nine healthy relatives belonging to twenty-nine families were recruited for this study, obtaining their clinical and demographic features using a data collection form, as well as blood samples for biochemical analysis. The mean age of patients was 36.8 years (ranging from 4 to 72 years). Positive family history of HAE was reported in 13 patients (32.5%), and a mean diagnosis delay of 7.9 (±12.5) years was estimated, ranging from months to 50 years. Cutaneous edema was the most common symptom (53.6%), while airway symptoms was present in 11 patients. Prophylactic treatment was indicated for 23 patients, while 14 also require on-demand rescue treatment. We estimate a minimum prevalence of 1.25:100,000 for HAE due to C1-INH deficiency or dysfunction in the Canary Islands, which is higher than the estimates for mainland Spanish populations. HAE continues to be a disease poorly recognized by health care professionals due to its confusing symptoms, leading to longer diagnosis delay. Altogether, the evidence reinforces the need for a rapid and accurate diagnosis and precision medicine-based studies to improve the patient's quality of life.

Genetic Ancestry Inference and Its Application for the Genetic Mapping of Human Diseases

Admixed populations arise when two or more ancestral populations interbreed. As a result of this admixture, the genome of admixed populations is defined by tracts of variable size inherited from these parental groups and has particular genetic features that provide valuable information about their demographic history. Diverse methods can be used to derive the ancestry apportionment of admixed individuals, and such inferences can be leveraged for the discovery of genetic loci associated with diseases and traits, therefore having important biomedical implications. In this review article, we summarize the most common methods of global and local genetic ancestry estimation and discuss the use of admixture mapping studies in human diseases.

Human Genomic Diversity Where the Mediterranean Joins the Atlantic

Throughout the past few years, a lively debate emerged about the timing and magnitude of the human migrations between the Iberian Peninsula and the Maghreb. Several pieces of evidence, including archaeological, anthropological, historical, and genetic data, have pointed to a complex and intermingled evolutionary history in the western Mediterranean area. To study to what extent connections across the Strait of Gibraltar and surrounding areas have shaped the present-day genomic diversity of its populations, we have performed a screening of 2.5 million single-nucleotide polymorphisms in 142 samples from southern Spain, southern Portugal, and Morocco. We built comprehensive data sets of the studied area and we implemented multistep bioinformatic approaches to assess population structure, demographic histories, and admixture dynamics. Both local and global ancestry inference showed an internal substructure in the Iberian Peninsula, mainly linked to a differential African ancestry. Western Iberia, from southern Portugal to Galicia, constituted an independent cluster within Iberia characterized by an enriched African genomic input. Migration time modeling showed recent historic dates for the admixture events occurring both in Iberia and in the North of Africa. However, an integrative vision of both paleogenomic and modern DNA data allowed us to detect chronological transitions and population turnovers that could be the result of transcontinental migrations dating back from Neolithic times. The present contribution aimed to fill the gaps in the modern human genomic record of a key geographic area, where the Mediterranean and the Atlantic come together.

The incidence of skin melanoma in Gran Canaria (Canary Islands, Spain) is lower than expected in Southern Europe despite high-risk environmental conditions: an island-wide cross-sectional study

BACKGROUND

The Canary Islands are a leading European touristic destination. The ultraviolet index (UVI) in the region is the highest in Spain, and similar to indexes registered in Australia and New Zealand, which hold the highest incidence of skin melanoma worldwide. Yet according to cancer registry data, the incidence in the Canary Islands in the late 1990s was the lowest in Spain (among the lowest in Europe) and about six times lower than in New Zealand.

PURPOSE

To analyze the incidence rates of skin melanoma in Gran Canaria island between 2007 and 2018.

METHODS

The study was based in the two centres of the Canary Islands' Healthcare Service centralizing melanoma care in Gran Canaria. We analyzed crude and age-standardized (ASR) incidence rates of invasive cutaneous melanoma for the period 2007-2018 following the inclusion criteria of the International Agency for Research on Cancer (IARC). Clinical and histological characteristics of melanoma patients were assessed.

RESULTS

A total of 1058 patients were included. The incidence rates obtained matched the latest available Canary Islands' cancer registry data, confirming its reliability (ASR, Segi-Doll world standard population: 6.4 cases per 100,000 habitants for 2008-2012). The incidence was also below the latest IARC predictions for Southern Europe (GLOBOCAN 2018). Histological characteristics of patients were similar to other Southern European series.

CONCLUSIONS

The incidence of skin melanoma in Gran Canaria is unexpectedly low for a Southern European population exposed to such a high UVI. Further research in the Canary Islands could provide insight into a better understanding of melanoma pathogenesis.

Admixture mapping of asthma in southwestern Europeans with North African ancestry influences

The prevalence of asthma symptoms in Canary Islanders, a southwestern European population from Spain, is almost three times higher than the country average. Because the genetic risks identified so far explain <5% of asthma heritability, here we aimed to discover new asthma loci by completing the first admixture mapping study in Canary Islanders leveraging their distinctive genetic makeup, where significant northwest African influences coexist in the European genetic diversity landscape. A 2-stage study was conducted in 1,491 unrelated individuals self-declaring having a Canary Islands origin for the 4 grandparents. Local ancestry estimates were obtained for the shared positions with reference data from putative ancestral populations from Europe, North Africa, and sub-Saharan Africa. Case-control deviations in local ancestry were tested for each ancestry separately using logistic regressions adjusted for principal components, followed by fine-mapping analyses based on imputed genotypes and analyses of the likely deleterious exonic variants. The admixture mapping analysis of asthma detected that local North African ancestry in a locus spanning 365 kb of chromosome 16q23.3 was associated with asthma risk at study-wide significance [lowest P = 1.12 × 10^-4; odds ratio (OR) = 2.05; 95% confidence interval (CI) = 1.41-3.00]. Fine-mapping studies identified a variant associated with asthma, and results were replicated in independent samples (rs3852738, OR = 1.34; 95% CI = 1.13-1.59, P = 7.58 × 10^-4). Whole exome sequencing data from a subset of individuals revealed an enrichment of likely deleterious variants among asthma cases in 16q23.3, particularly in the phospholipase Cγ2 (PLCG2) gene (P = 3.67 × 10^-4). By completing the first mapping study of asthma in admixed populations from Europe, our results revealed a new plausible asthma locus.

Sepsis-associated acute respiratory distress syndrome in individuals of European ancestry: a genome-wide association study

BACKGROUND

Acute respiratory distress syndrome (ARDS) is a lung inflammatory process caused mainly by sepsis. Most previous studies that identified genetic risks for ARDS focused on candidates with biological relevance. We aimed to identify novel genetic variants associated with ARDS susceptibility and to provide complementary functional evidence of their effect in gene regulation.

METHODS

We did a case-control genome-wide association study (GWAS) of 1935 European individuals, using patients with sepsis-associated ARDS as cases and patients with sepsis without ARDS as controls. The discovery stage included 672 patients admitted into a network of Spanish intensive care units between January, 2002, and January, 2017. The replication stage comprised 1345 individuals from two independent datasets from the MESSI cohort study (Sep 22, 2008-Nov 30, 2017; USA) and the VISEP (April 1, 2003-June 30, 2005) and MAXSEP (Oct 1, 2007-March 31, 2010) trials of the SepNet study (Germany). Results from discovery and replication stages were meta-analysed to identify association signals. We then used RNA sequencing data from lung biopsies, in-silico analyses, and luciferase reporter assays to assess the functionallity of associated variants.

FINDINGS

We identified a novel genome-wide significant association with sepsis-associated ARDS susceptibility (rs9508032, odds ratio [OR] 0·61, 95% CI 0·41-0·91, p=5·18 × 10^-8) located within the Fms-related tyrosine kinase 1 (FLT1) gene, which encodes vascular endothelial growth factor receptor 1 (VEGFR-1). The region containing the sentinel variant and its best proxies acted as a silencer for the FLT1 promoter, and alleles with protective effects in ARDS further reduced promoter activity (p=0·0047). A literature mining of all previously described ARDS genes validated the association of vascular endothelial growth factor A (VEGFA; OR 0·55, 95% CI 0·41-0·73; p=4·69 × 10^-5).

INTERPRETATION

A common variant within the FLT1 gene is associated with sepsis-associated ARDS. Our findings support a role for the vascular endothelial growth factor signalling pathway in ARDS pathogenesis and identify VEGFR-1 as a potential therapeutic target.

FUNDING

Instituto de Salud Carlos III, European Regional Development Funds, Instituto Tecnológico y de Energías Renovables.

Mitogenomes illuminate the origin and migration patterns of the indigenous people of the Canary Islands

The Canary Islands’ indigenous people have been the subject of substantial archaeological, anthropological, linguistic and genetic research pointing to a most probable North African Berber source. However, neither agreement about the exact point of origin nor a model for the indigenous colonization of the islands has been established. To shed light on these questions, we analyzed 48 ancient mitogenomes from 25 archaeological sites from the seven main islands. Most lineages observed in the ancient samples have a Mediterranean distribution, and belong to lineages associated with the Neolithic expansion in the Near East and Europe (T2c, J2a, X3a…). This phylogeographic analysis of Canarian ancient mitogenomes, the first of its kind, shows that some lineages are restricted to Central North Africa (H1cf, J2a2d and T2c1d3), while others have a wider distribution, including both West and Central North Africa, and, in some cases, Europe and the Near East (U6a1a1, U6a7a1, U6b, X3a, U6c1). In addition, we identify four new Canarian-specific lineages (H1e1a9, H4a1e, J2a2d1a and L3b1a12) whose coalescence dates correlate with the estimated time for the colonization of the islands (1^st millennia CE). Additionally, we observe an asymmetrical distribution of mtDNA haplogroups in the ancient population, with certain haplogroups appearing more frequently in the islands closer to the continent. This reinforces results based on modern mtDNA and Y-chromosome data, and archaeological evidence suggesting the existence of two distinct migrations. Comparisons between insular populations show that some populations had high genetic diversity, while others were probably affected by genetic drift and/or bottlenecks. In spite of observing interinsular differences in the survival of indigenous lineages, modern populations, with the sole exception of La Gomera, are homogenous across the islands, supporting the theory of extensive human mobility after the European conquest.