Indigenous Arabs are descendants of the earliest split from ancient Eurasian populations

Ancient genomes illuminate Eastern Arabian population history and adaptation against malaria

The harsh climate of Arabia has posed challenges in generating ancient DNA from the region, hindering the direct examination of ancient genomes for understanding the demographic processes that shaped Arabian populations. In this study, we report whole-genome sequence data obtained from four Tylos-period individuals from Bahrain. Their genetic ancestry can be modeled as a mixture of sources from ancient Anatolia, Levant, and Iran/Caucasus, with variation between individuals suggesting population heterogeneity in Bahrain before the onset of Islam. We identify the G6PD Mediterranean mutation associated with malaria resistance in three out of four ancient Bahraini samples and estimate that it rose in frequency in Eastern Arabia from 5 to 6 kya onward, around the time agriculture appeared in the region. Our study characterizes the genetic composition of ancient Arabians, shedding light on the population history of Bahrain and demonstrating the feasibility of studies of ancient DNA in the region.

North and East African mitochondrial genetic variation needs further characterization towards precision medicine

INTRODUCTION

Mitochondria are maternally inherited cell organelles with their own genome, and perform various functions in eukaryotic cells such as energy production and cellular homeostasis. Due to their inheritance and manifold biological roles in health and disease, mitochondrial genetics serves a dual purpose of tracing the history as well as disease susceptibility of human populations across the globe. This work requires a comprehensive catalogue of commonly observed genetic variations in the mitochondrial DNAs for all regions throughout the world. So far, however, certain regions, such as North and East Africa have been understudied.

OBJECTIVES

To address this shortcoming, we have created the most comprehensive quality-controlled North and East African mitochondrial data set to date and use it for characterizing mitochondrial genetic variation in this region.

METHODS

We compiled 11 published cohorts with novel data for mitochondrial genomes from 159 Sudanese individuals. We combined these 641 mitochondrial sequences with sequences from the 1000 Genomes (n = 2504) and the Human Genome Diversity Project (n = 828) and used the tool haplocheck for extensive quality control and detection of in-sample contamination, as well as Nanopore long read sequencing for haplogroup validation of 18 samples.

RESULTS

Using a subset of high-coverage mitochondrial sequences, we predict 15 potentially novel haplogroups in North and East African subjects and observe likely phylogenetic deviations from the established PhyloTree reference for haplogroups L0a1 and L2a1.

CONCLUSION

Our findings demonstrate common hitherto unexplored variants in mitochondrial genomes of North and East Africa that lead to novel phylogenetic relationships between haplogroups present in these regions. These observations call for further in-depth population genetic studies in that region to enable the prospective use of mitochondrial genetic variation for precision medicine.

Eurasian back-migration into Northeast Africa was a complex and multifaceted process

Recent studies have identified Northeast Africa as an important area for human movements during the Holocene. Eurasian populations have moved back into Northeastern Africa and contributed to the genetic composition of its people. By gathering the largest reference dataset to date of Northeast, North, and East African as well as Middle Eastern populations, we give new depth to our knowledge of Northeast African demographic history. By employing local ancestry methods, we isolated the Non-African parts of modern-day Northeast African genomes and identified the best putative source populations. Egyptians and Sudanese Copts bore most similarities to Levantine populations whilst other populations in the region generally had predominantly genetic contributions from the Arabian peninsula rather than Levantine populations for their Non-African genetic component. We also date admixture events and investigated which factors influenced the date of admixture and find that major linguistic families were associated with the date of Eurasian admixture. Taken as a whole we detect complex patterns of admixture and diverse origins of Eurasian admixture in Northeast African populations of today.

The Qatari population's genetic structure and gene flow as revealed by the Y chromosome

The Y-chromosome has been widely used in forensic genetic applications and human population genetic studies due to its uniparental origins. A large database on the Qatari population was created for comparison with other databases from the Arabian Peninsula, the Middle East, and Africa. We provide a study of 23 Y-STR loci included in PowerPlex Y23 (Promega, USA) that were genotyped to produce haplotypes in 379 unrelated males from Qatar, a country at the crossroads of migration patterns. Overall, the most polymorphic locus provided by the Promega kit was DYS458, with a genetic diversity value of 0.85 and a haplotype diversity of 0.998924. Athey's Haplogroup Predictor tool was used to predict haplogroups from Y-STR haplotypes in the Qatari population. In a median-joining network, the haplogroup J1 predominance (49%) in Qatar generated a star-like expansion cluster. The graph of population Q-matrix was developed using Y-STR data from 38 Middle Eastern and 97 African populations (11,305 individuals), and it demonstrated a stronger sub-grouping of countries within each ethnic group and showed the effect of Arabs on the indigenous Berbers of North Africa. The estimated migration rate between the Qatari and other Arabian populations was inferred using Bayesian coalescence theory in the Migrate-n program. According to the Gene Flow study, the main migration route was from Yemen to Kuwait through Qatar. Our research, using the PowerPlex Y23 database, shows the importance of gene diversity, as well as regional and social structuring, in determining the utility of demographic and forensic databases.

Investigation of germline variants in Bahraini women with breast cancer using next-generation sequencing based-multigene panel

Germline variants in BRCA1 and BRCA2 (BRCA1/2) genes are the most common cause of hereditary breast cancer. However, a significant number of cases are not linked to these two genes and additional high-, moderate- and low-penetrance genes have been identified in breast cancer. The advent of next-generation sequencing (NGS) allowed simultaneous sequencing of multiple cancer-susceptibility genes and prompted research in this field. So far, cancer-predisposition genes other than BRCA1/2 have not been studied in the population of Bahrain. We performed a targeted NGS using a multi-panel covering 180 genes associated with cancer predisposition to investigate the spectrum and frequency of germline variants in 54 women with a positive personal and/or family history of breast cancer. Sequencing analysis revealed germline variants in 29 (53.7%) patients. Five pathogenic/likely pathogenic variants in four DNA repair pathway-related genes were identified in five unrelated patients (9.3%). Two BRCA1 variants, namely the missense variant c.287A>G (p.Asp96Gly) and the truncating variant c.1066C>T (p.Gln356Ter), were detected in two patients (3.7%). Three variants in non-BRCA1/2 genes were detected in three patients (1.85% each) with a strong family history of breast cancer. These included a monoallelic missense variant c.1187G>A (p.Gly396Asp) in MUTYH gene, and two truncating variants namely c.3343C>T (p.Arg1115Ter) in MLH3 gene and c.1826G>A (p.Trp609Ter) in PMS1 gene. Other variants of uncertain significance (VUS) were also detected, and some of them were found together with the deleterious variants. In this first application of NGS-based multigene testing in Bahraini women with breast cancer, we show that multigene testing can yield additional genomic information on low-penetrance genes, although the clinical significance of these genes has not been fully appreciated yet. Our findings also provide valuable epidemiological information for future studies and highlight the importance of genetic testing, and an NGS-based multigene analysis may be applied supplementary to traditional genetic counseling.

A whole-genome reference panel of 14,393 individuals for East Asian populations accelerates discovery of rare functional variants

Underrepresentation of non-European (EUR) populations hinders growth of global precision medicine. Resources such as imputation reference panels that match the study population are necessary to find low-frequency variants with substantial effects. We created a reference panel consisting of 14,393 whole-genome sequences including more than 11,000 Asian individuals. Genome-wide association studies were conducted using the reference panel and a population-specific genotype array of 72,298 subjects for eight phenotypes. This panel yields improved imputation accuracy of rare and low-frequency variants within East Asian populations compared with the largest reference panel. Thirty-nine previously unidentified associations were found, and more than half of the variants were East Asian specific. We discovered genes with rare protein-altering variants, including LTBP1 for height and GPR75 for body mass index, as well as putative regulatory mechanisms for rare noncoding variants with cell type-specific effects. We suggest that this dataset will add to the potential value of Asian precision medicine.

Aesthetic considerations for treating the Middle Eastern patient: Thriving in Diversity international roundtable series

BACKGROUND

The Middle East has a significant influence on the global aesthetic market. Within the United States and globally, patients with Middle Eastern heritage have a wide range of ethnic and cultural backgrounds that affect their perceptions of beauty and motivations to seek cosmetic treatment.

AIMS

The aim of this roundtable was to discuss similarities and differences in anatomy and treatment preferences of Middle Eastern patients and explore how these differences may influence aesthetic practices.

PATIENTS/METHODS

In support of clinicians who wish to serve a diverse patient population, a 6-part international roundtable series focused on diversity in aesthetics was conducted from August 24, 2021, to May 16, 2022.

RESULTS

The results of the fourth roundtable in the series, the Middle Eastern Patient, are described here. A discussion of treatment preferences is included, and specific procedural information is provided for  commonly treated areas in this population (forehead, infraorbital area, and jawline).

CONCLUSIONS

Middle Eastern patients have a variety of aesthetic preferences, which are influenced by a wide range of cultural backgrounds, making it difficult to develop general statements about this demographic. There is an unmet need for research into this diverse group of patients to help physicians understand and incorporate their unique needs and desires into clinical practice.

Population Genome Programs across the Middle East and North Africa: Successes, Challenges, and Future Directions

In this review, we discuss the current state of population genome programs (PGPs) conducted in the Middle East and North African (MENA) region. This region has high prevalence of genetic diseases and significant health challenges as well as being a significantly underrepresented population in public genetic databases. The majority of ongoing PGPs represent regions in Europe, North and South America, South Asia, Australia, and Africa, with little to no descriptive information highlighted only on the MENA Region when it comes to genome programs databases, outcomes, or the challenges that MENA region countries may face establishing their own national programs. This review has identified 6 PGPs currently underway in the MENA region, namely in the Kingdom of Saudi Arabia, Qatar, Egypt, the United Arab Emirates, Bahrain, and Iran. Due to the rapidly growing involvement of the MENA region in national-scale genomic data collection, an increase in representation in public genetic databases is to be expected to occur in the near future. Whilst significant progress is being made in some MENA countries, future initiatives as well as ongoing programs will be facing several challenges related to collaboration, finance, infrastructure and institutional data access, data analysis, sustainability, health records, and biobanks. The review also reiterates the need for ensuring ethical and regulated genomic initiatives which can drive developments in personalized medicine treatments to improve patient prognosis and quality of life.

Validation of Polygenic Risk Scores for Coronary Heart Disease in a Middle Eastern Cohort Using Whole Genome Sequencing

BACKGROUND

Enthusiasm for using polygenic risk scores (PRSs) in clinical practice is tempered by concerns about their portability to diverse ancestry groups, thus motivating genome-wide association studies in non-European ancestry cohorts.

METHODS

We conducted a genome-wide association study for coronary heart disease in a Middle Eastern cohort using whole genome sequencing and assessed the performance of 6 PRSs developed with methods including LDpred (PGS000296), metaGRS (PGS000018), Pruning and Thresholding (PGS000337), and an EnsemblePRS we developed. Additionally, we evaluated the burden of rare variants in lipid genes in cases and controls. Whole genome sequencing at 30× coverage was performed in 1067 coronary heart disease cases (mean age=59 years; 70.3% males) and 6170 controls (mean age=40 years; 43.5% males).

RESULTS

The majority of PRSs performed well; odds ratio (OR) per 1 SD increase (OR_1sd) was highest for PGS000337 (OR_1sd=1.81, 95% CI [1.66-1.98], P=3.07×10^-41). EnsemblePRS performed better than individual PRSs (OR_1sd=1.8, 95% CI [1.66-1.96], P=5.89×10^-44). The OR for the 10th decile versus the remaining deciles was >3.2 for PGS000337, PGS000296, PGS000018, and reached 4.58 for EnsemblePRS. Of 400 known genome-wide significant loci, 33 replicated at P<10^-4. However, the 9p21 locus did not replicate. Six suggestive (P<10^-5) new loci/genes with plausible biological function were identified (eg, CORO7, RBM47, PDE4D). The burden of rare functional variants in LDLR, APOB, PCSK9, and ANGPTL4 was greater in cases than controls.

CONCLUSIONS

Overall, we demonstrate that PRSs derived from European ancestry genome-wide association studies performed well in a Middle Eastern cohort, suggesting these could be used in the clinical setting while ancestry-specific PRSs are developed.

Frequency of functional exonic single-nucleotide polymorphisms and haplotype distribution in the SLCO1B1 gene across genetic ancestry groups in the Qatari population

Organic anion transporting polypeptides (OATP), which are encoded by SLCO genes, participate in the hepatic elimination of drugs and xenobiotics. SLCO1B1 is an important pharmacogenomic gene (encoding OATP1B1) associated with response to the uptake of endogenous compounds, such as statin and bilirubin. Ethnicity of the patient modulates the response to these drugs; the frequency and haplotype data for SLCO1B1 genetic variants in the Arab population is lacking. Therefore, we determined the frequencies of two well-characterized SLCO1B1 single nucleotide polymorphisms (SNP) and haplotypes that affect the OATP1B1 drugs transportation activity in Qatari population. Genotyping data for two SLCO1B1 SNPs (c.388A > G, c.521 T > C) were extracted from whole exome data of 1050 Qatari individuals, who were divided into three ancestry groups, namely Bedouins, Persians/South Asians, and Africans. By way of using Fisher's exact and Chi-square tests, we evaluated the differences in minor allele frequency (MAF) of the two functional SNPs and haplotype frequencies (HF) among the three ancestry groups. The OATP1B1 phenotypes were assigned according to their function by following the guidelines from the Clinical Pharmacogenetics Implementation Consortium for SLCO1B1 and Simvastatin-Induced Myopathy.The MAF of SLCO1B1:c.388A > G was higher compared to that of SLCO1B1:c.521 T > C in the study cohort. It was significantly high in the African ancestry group compared with the other two groups, whereas SLCO1B1:c.521 T > C was significantly low in the African ancestry group compared with the other two groups. The SLCO1B1 *15 haplotype had the highest HF, followed by *1b, *1a, and *5. Only the SLCO1B1 *5 haplotype showed no significant difference in frequency across the three ancestry groups. Furthermore, we observed that the OATP1B1 normal function phenotype accounted for 58% of the Qatari individuals, the intermediate function phenotype accounted for 35% with significant differences across the ancestry groups, and the low function phenotype accounted for 6% of the total Qatari individuals with a higher trend observed in the Bedouin group.The results indicate that the phenotype frequencies of the OATP1B1 intermediate and low function in the Qatari population appear at the higher end of the frequency range seen worldwide. Thus, a pharmacogenetic screening program for SLCO1B1 variants may be necessary for the Qatari population.

Distribution of HLA-B Alleles and Haplotypes in Qatari: Recommendation for Establishing Pharmacogenomic Markers Screening for Drug Hypersensitivity

Human leukocyte antigen (HLA) proteins are present at the cellular surface of antigen-presenting cells and play a crucial role in the adaptive immune response. Class I genes, specifically certain HLA-B alleles, are associated with adverse drug reactions (ADRs) and are used as pharmacogenetic markers. Although ADRs are a common causes of hospitalization and mortality, the data on the prevalence of HLA-B pharmacogenetics markers in Arab countries are scarce. In this study, we investigated the frequencies of major HLA-B pharmacogenomics markers in the Qatari population. Next-generation sequencing data from 1,098 Qatari individuals were employed for HLA-B typing using HLA-HD version 1.4.0 and IPD-IMGT/HLA database. In addition, HLA-B pharmacogenetics markers were obtained from the HLA Adverse Drug Reaction Database. In total, 469 major HLA-B pharmacogenetic markers were identified, with HLA-B*51:01 being the most frequent pharmacogenetic marker (26.67%) in the Qatari population. Moreover, HLA-B*51:01 is associated with phenytoin- and clindamycin-induced ADRs. The second most frequent pharmacogenetic marker was the HLA-B*58:01 allele (6.56%), which is associated with allopurinol-induced ADRs. The third most frequent pharmacogenetic marker was the HLA-B*44:03 allele, which is associated with phenytoin-induced ADRs. The establishment of a pharmacogenetics screening program in Qatar for cost effective interventions aimed at preventing drug-induced hypersensitivity can be aided by the highly prevalent HLA-B pharmacogenetic markers detected here.

HLA-A, -B, -C, -DRB1 and -DQB1 allele and haplotype frequencies in Lebanese and their relatedness to neighboring and distant populations

BACKGROUND

This study examined the origin of present-day Lebanese using high-resolution HLA class I and class II allele and haplotype distributions. The study subjects comprised 152 unrelated individuals, and their HLA class I and class II alleles and two-locus and five-locus haplotypes were compared with those of neighboring and distant communities using genetic distances, neighbor-joining dendrograms, correspondence, and haplotype analyses. HLA class I (A, B, C) and class II (DRB1, DQB1) were genotyped at a high-resolution level by PCR-SSP.

RESULTS

In total, 76 alleles across the five HLA loci were detected: A*03:01 (17.1%), A*24:02 (16.5%), B*35:01 (25.7%), C*04:01 (25.3%), and C*07:01 (20.7%) were the most frequent class I alleles, while DRB1*11:01 (34.2%) and DQB1*03:01 (43.8%) were the most frequent class II alleles. All pairs of HLA loci were in significant linkage disequilibrium. The most frequent two-locus haplotypes recorded were DRB1*11:01 ~ DQB1*03:01 (30.9%), B*35:01-C*04:01 (20.7%), B*35:01 ~ DRB1*11:01 (13.8%), and A*24:02 ~ B*35:01 (10.3%). Lebanese appear to be closely related to East Mediterranean communities such as Levantines (Palestinians, Syrians, and Jordanians), Turks, Macedonians, and Albanians. However, Lebanese appear to be distinct from North African, Iberian, and Sub-Saharan communities.

CONCLUSIONS

Collectively, this indicates a limited genetic contribution of Arabic-speaking populations (from North Africa or the Arabian Peninsula) and Sub-Saharan communities to the present-day Lebanese gene pool. This confirms the notion that Lebanese population are of mixed East Mediterranean and Asian origin, with a marked European component.

Mixed model-based eQTL analysis reveals lncRNAs associated with regulation of genes involved in sex determination and spermatogenesis: The key to understanding human gender imbalance

BACKGROUND

An imbalance in the prenatal sex ratio in humans may be due to several factors affecting sperm physiology, including genetic features. In this study, we conducted a transcriptome-wide analysis of expression quantitative trait loci (eQTLs) to identify target genes associated with previously described QTLs associated with gender imbalance.

METHODS

A mixed model explaining polygenic effects by genomic covariance among individuals was used to identify the eQTLs using gene expression and genotype data from 462 European/African individuals.

RESULTS

Eight eGenes were associated with four QTLs (P < 4.00 × 10^-5), with strong associations found (P < 4.00 × 10^-8) between rs2485007 and eGenes ANKRD26P3 (P = 3.40 × 10^-9) and LINC00421 (P = 1.35 × 10^-9). ANKRD26P3 and LINC00421 are both lncRNAs associated with the control of testis-dominant genes PELP1, TAF15, NANOG, TEX14, TCF3, ZNF433, ZNF555, TEX37, FATE1, TCP11, and CYLC2 and Y-linked genes SRY and ZFY, as well as several genes with roles in spermatogenesis (ODF1, SPATC1, SPATA3, SPATA31E1, SPERT, SPATA16, MOSPD1, SPATA24, and SPO11) and sex determination (SOX family genes).

CONCLUSIONS

The above eGenes contribute directly or indirectly to gene regulation for sex determination and spermatogenesis, thereby serving as important functional clues for gender-biased selection.

Qatar Genome: Insights on Genomics from the Middle East

Despite recent biomedical breakthroughs and large genomic studies growing momentum, the Middle Eastern population, home to over 400 million people, is under-represented in the human genome variation databases. Here we describe insights from phase 1 of the Qatar Genome Program with whole genome sequenced 6,047 individuals from Qatar. We identified more than 88 million variants of which 24 million are novel and 23 million are singletons. Consistent with the high consanguinity and founder effects in the region, we found that several rare deleterious variants were more common in the Qatari population while others seem to provide protection against diseases and have shaped the genetic architecture of adaptive phenotypes. These results highlight the value of our data as a resource to advance genetic studies in the Arab and neighbouring Middle Eastern populations and will significantly boost the current efforts to improve our understanding of global patterns of human variations, human history and genetic contributions to health and diseases in diverse populations. This article is protected by copyright. All rights reserved.

The QChip1 knowledgebase and microarray for precision medicine in Qatar

Risk genes for Mendelian (single-gene) disorders (SGDs) are consistent across populations, but pathogenic risk variants that cause SGDs are typically population-private. The goal was to develop "QChip1," an inexpensive genotyping microarray to comprehensively screen newborns, couples, and patients for SGD risk variants in Qatar, a small nation on the Arabian Peninsula with a high degree of consanguinity. Over 108 variants in 8445 Qatari were identified for inclusion in a genotyping array containing 165,695 probes for 83,542 known and potentially pathogenic variants in 3438 SGDs. QChip1 had a concordance with whole-genome sequencing of 99.1%. Testing of QChip1 with 2707 Qatari genomes identified 32,674 risk variants, an average of 134 pathogenic alleles per Qatari genome. The most common pathogenic variants were those causing homocystinuria (1.12% risk allele frequency), and Stargardt disease (2.07%). The majority (85%) of Qatari SGD pathogenic variants were not present in Western populations such as European American, South Asian American, and African American in New York City and European and Afro-Caribbean in Puerto Rico; and only 50% were observed in a broad collection of data across the Greater Middle East including Kuwait, Iran, and United Arab Emirates. This study demonstrates the feasibility of developing accurate screening tools to identify SGD risk variants in understudied populations, and the need for ancestry-specific SGD screening tools.

Mutation Screening of Six Exons of ABCA4 in Iranian Stargardt Disease Patients

Purpose

Stargardt disease type 1 (STGD1) is a recessively inherited retinal disorder that can cause severe visual impairment. ABCA4 mutations are the usual cause of STGD1. ABCA4 codes a transporter protein exclusively expressed in retinal photoreceptor cells. The genecontains 50 exons. Mutations are most frequent in exons 3, 6, 12, and 13, and exons 10 and 42 each contain two common variations. We aimed to screen these exons for mutations in Iranian STGD1 patients.

Methods

Eighteen STGD1 patients were recruited for genetic analysis. Diagnosis by retina specialists was based on standard criteria, including accumulation of lipofuscin. The six ABCA4 exons were PCR amplified and sequenced by the Sanger method.

Results

One or more ABCA4-mutated alleles were identified in 5 of the 18 patients (27.8%). Five different mutations including two splice site (c.1356+1G>A and c.5836-2A>G) and three missense mutations (p.Gly1961Glu, p.Gly1961Arg, and p.Gly550Arg) were found. The p.Gly1961Glu mutation was the only mutation observed in two patients.

Conclusion

As ABCA4 mutations in exons 6, 12, 10, and 42 were identified in approximately 25% of the patients studied, these may be appropriate exons for screening projects. As in other populations, STDG1 causative ABCA4 mutations are heterogeneous among Iranian patients, and p.Gly1961Glu may be relatively frequent.

Genomewide association analysis of warfarin dose requirements in Middle Eastern and North African populations

To date, there has been no genomewide association study (GWAS) from the Middle East and North African (MENA) region to identify genetic variants associated with warfarin dose variability using this approach. In this study, we aimed to conduct the first GWAS of warfarin dose requirements in patients from the MENA region. A total of 132 Qatari (discovery) and 50 Egyptians (replication) were genotyped using Illumina Multi‐Ethnic Global BeadChip Array. A GWAS was performed on log‐transformed weekly warfarin dose in the studied population, adjusting for clinical characteristics and ancestry. The genomewide signals from the discovery cohort were tested in the Egyptian cohort. A GWAS meta‐analysis, including the Qatari and Egyptian cohorts, was also performed and the output from this analysis was used in a gene‐based analysis. The discovery analysis in Qatari identified five genomewide single‐nucleotide polymorphisms (SNPs) in chromosome 16. These signals were replicated in the Egyptian cohort. Combining the two data through a GWAS meta‐analysis strengthened the association in chromosome 16 with VKORC1 rs9934438 being the lead genomewide signal (β = −0.17, 6 × 10⁻¹⁵). Other SNPs were identified in chromosome 10 at a p value less than 1 × 10⁻⁵. The genetic variants within VKORC1 rs9934438 and CYP2C9 rs4086116 explained 39% and 27% of the variability in the weekly warfarin dose requirement in the Qatari and Egyptians, respectively. This is the first GWAS of warfarin dose variability in the MENA region. It confirms the importance of VKORC1 and CYP2C9 variants in warfarin dose variability among patients from the MENA region.

Thousands of Qatari genomes inform human migration history and improve imputation of Arab haplotypes

Arab populations are largely understudied, notably their genetic structure and history. Here we present an in-depth analysis of 6,218 whole genomes from Qatar, revealing extensive diversity as well as genetic ancestries representing the main founding Arab genealogical lineages of Qahtanite (Peninsular Arabs) and Adnanite (General Arabs and West Eurasian Arabs). We find that Peninsular Arabs are the closest relatives of ancient hunter-gatherers and Neolithic farmers from the Levant, and that founder Arab populations experienced multiple splitting events 12–20 kya, consistent with the aridification of Arabia and farming in the Levant, giving rise to settler and nomadic communities. In terms of recent genetic flow, we show that these ancestries contributed significantly to European, South Asian as well as South American populations, likely as a result of Islamic expansion over the past 1400 years. Notably, we characterize a large cohort of men with the ChrY J1a2b haplogroup (n = 1,491), identifying 29 unique sub-haplogroups. Finally, we leverage genotype novelty to build a reference panel of 12,432 haplotypes, demonstrating improved genotype imputation for both rare and common alleles in Arabs and the wider Middle East.

Arab populations are relatively understudied, especially their genetic architecture and historical relationship with early founders of the ancient Near East. Here, the authors examine 6,218 Qatari whole genomes, revealing insights on migration, population history and genetic structure of populations across the Middle Eastern region.

Clinical features, epidemiology, autoantibody status, HLA haplotypes and genetic mechanisms of type 1 diabetes mellitus among children in Qatar

To describe the clinical features, epidemiology, autoantibody status, HLA haplotypes and genetic mechanisms of type 1 diabetes mellitus (T1DM). Patients (0-18 years) with diabetes were recruited. Clinical data was collected, autoantibodies and c-peptide were measured. Whole Genome Sequencing was performed. Genomic data analysis was compared with the known genes linked with T1DM and HLA alleles were studied. 1096 patients had one or more antibody positivity. The incidence of T1DM in 2020 was 38.05 per 100,000 children and prevalence was 249.73. GADA was the most common autoantibody followed by IAA. Variants in GSTCD, SKAP2, SLC9B1, BANK1 were most prevalent. An association of HLA haplotypes DQA1*03:01:01G (OR = 2.46, p value = 0.011) and DQB1*03:02:01G (OR = 2.43, p value = 0.022) was identified. The incidence of T1DM in Qatar is the fourth highest in the world, IA2 autoantibody was the most specific with some patients only having ZnT8 or IA2 autoantibodies thus underlining the necessity of profiling all 4 autoantibodies. The genes associated with T1DM in the Arab population were different from those that are common in the Caucasian population. HLA-DQ was enriched in the Qatari patients suggesting that it can be considered a major risk factor at an early age.

Projecting Ancient Ancestry in Modern-Day Arabians and Iranians: A Key Role of the Past Exposed Arabo-Persian Gulf on Human Migrations

The Arabian Peninsula is strategic for investigations centered on the early structuring of modern humans in the wake of the out-of-Africa migration. Despite its poor climatic conditions for the recovery of ancient human DNA evidence, the availability of both genomic data from neighboring ancient specimens and informative statistical tools allow modeling the ancestry of local modern populations. We applied this approach to a data set of 741,000 variants screened in 291 Arabians and 78 Iranians, and obtained insightful evidence. The west-east axis was a strong forcer of population structure in the Peninsula, and, more importantly, there were clear continuums throughout time linking western Arabia with the Levant, and eastern Arabia with Iran and the Caucasus. Eastern Arabians also displayed the highest levels of the basal Eurasian lineage of all tested modern-day populations, a signal that was maintained even after correcting for a possible bias due to a recent sub-Saharan African input in their genomes. Not surprisingly, eastern Arabians were also the ones with highest similarity with Iberomaurusians, who were, so far, the best proxy for the basal Eurasians amongst the known ancient specimens. The basal Eurasian lineage is the signature of ancient non-Africans who diverged from the common European-eastern Asian pool before 50,000 years ago, prior to the later interbred with Neanderthals. Our results appear to indicate that the exposed basin of the Arabo-Persian Gulf was the possible home of basal Eurasians, a scenario to be further investigated by searching ancient Arabian human specimens.

The genomic history of the Middle East

The Middle East region is important to understand human evolution and migrations but is underrepresented in genomic studies. Here, we generated 137 high-coverage physically phased genome sequences from eight Middle Eastern populations using linked-read sequencing. We found no genetic traces of early expansions out-of-Africa in present-day populations but found Arabians have elevated Basal Eurasian ancestry that dilutes their Neanderthal ancestry. Population sizes within the region started diverging 15–20 kya, when Levantines expanded while Arabians maintained smaller populations that derived ancestry from local hunter-gatherers. Arabians suffered a population bottleneck around the aridification of Arabia 6 kya, while Levantines had a distinct bottleneck overlapping the 4.2 kya aridification event. We found an association between movement and admixture of populations in the region and the spread of Semitic languages. Finally, we identify variants that show evidence of selection, including polygenic selection. Our results provide detailed insights into the genomic and selective histories of the Middle East.

•

Middle Easterners do not have ancestry from an early out-of-Africa expansion

•

Basal Eurasian and African ancestry in Arabians deplete their Neanderthal ancestry

•

Populations experienced bottlenecks overlapping aridification events

•

Identification of recent single and polygenic signals of selection in Arabia

Analysis of HLA gene polymorphisms in East Africans reveals evidence of gene flow in two Semitic populations from Sudan

Sudan, a northeastern African country, is characterized by high levels of cultural, linguistic, and genetic diversity, which is believed to be affected by continuous migration from neighboring countries. Consistent with such demographic effect, genome-wide SNP data revealed a shared ancestral component among Sudanese Afro-Asiatic speaking groups and non-African populations, mainly from West Asia. Although this component is shared among all Afro-Asiatic speaking groups, the extent of this sharing in Semitic groups, such as Sudanese Arab, is still unknown. Using genotypes of six polymorphic human leukocyte antigen (HLA) genes (i.e., HLA-A, -C, -B, -DRB1, -DQB1, and -DPB1), we examined the genetic structure of eight East African ethnic groups with origins in Sudan, South Sudan, and Ethiopia. We identified informative HLA alleles using principal component analysis, which revealed that the two Semitic groups (Gaalien and Shokrya) constituted a distinct cluster from the other Afro-Asiatic speaking groups in this study. The HLA alleles that distinguished Semitic Arabs co-exist in the same extended HLA haplotype, and those alleles are in strong linkage disequilibrium. Interestingly, we find the four-locus haplotype "C*12:02-B*52:01-DRB1*15:02-DQB1*06:01" exclusively in non-African populations and it is widely spread across Asia. The identification of this haplotype suggests a gene flow from Asia, and likely these haplotypes were brought to Africa through back migration from the Near East. These findings will be of interest to biomedical and anthropological studies that examine the demographic history of northeast Africa.

No impact of soluble epoxide hydrolase rs4149243, rs2234914 and rs751142 genetic variants on the development of type II diabetes and its hypertensive complication among Jordanian patients

BACKGROUND

Human soluble epoxide hydrolase plays a major role in cardiovascular homoeostasis. Genetic variants in the EPHX2 gene among different ethnic groups are associated with cardiovascular complications, such as hypertension. However, no reports regarding the association of EPHX2 genotype with hypertension among type II diabetic (T2D) patients of Middle Eastern Jordanian origin exist.

OBJECTIVE

The current study aimed to elucidate the association of the EPHX2 allele, genotype and haplotype with T2D, hypertension and parameters of lipid profile parameters among Jordanian T2D patients.

METHODS

Ninety-three genomic DNA samples of non-diabetic controls and 97 samples from T2D patients were genotyped for EPHX2 rs4149243, rs2234914 and rs751142 genetic variants. The DNA samples were amplified using polymerase chain reaction (PCR) and then sequenced using Applied Biosystems Model (ABI3730x1). The functionality of intronic EPHX2 variants was predicted using the in silico Berkely Drosophila Genome Project software.

RESULTS

We found no significant (P >.05) association between the EPHX2 rs4149243, rs2234914 and rs751142 allele, genotype and haplotype and the incidence of T2D and hypertension. Additionally, no association (P >.05) between these EPHX2 genetic variants with the baseline total cholesterol, low- and high-density lipoproteins and triglycerides among both non-diabetic and diabetic volunteers was found. However, we found an inter-ethnic variation (χ² -test, P value ˂ .05) in the allele frequency of the EPHX2 rs4149243 and rs2234914 variants between Jordanians and other ethnic populations. Also, the in silico Berkely Drosophila Genome Project software predicted that the intronic EPHX2 rs4149243 could alter the splicing of intron 7.

CONCLUSIONS

It can be concluded from this study that EPHX2 rs4149243, rs2234914 and rs751142 genetic variants do not play a role in the development of T2D and hypertension among Jordanian T2D patients. Further genetic studies with larger sample sizes are needed to find out the association of other functional EPHX2 variants with cardiovascular diseases among T2D patients in Jordan.

Mitochondrial haplogroup J associated with higher risk of obesity in the Qatari population

Obesity, a major risk factor for metabolic disorders, is highly prevalent in Qatari population. Maternal transmission of obesity traits can be significant; for example, X haplogroup is known to be associated with lower BMI and body fat mass in Northern Europeans and T haplogroup which is a sister haplogroup of J is known to be associated with obesity in Caucasian subjects from Austria and Southern Italy. We aimed to delineate the mitochondrial haplogroups and variants associated with obesity in Qatari population. Mitochondrial genomes of 864 Qatari individuals were extracted from whole exome sequencing data with an average coverage of 77X. We distributed the participants into 2 sub-cohorts: obese (BMI ≥ 30) and non-obese (BMI < 30); the mean value of BMI from these two groups were 36.5 ± 5.7 and 26.5 ± 2.6, respectively. Mitochondrial haplogroup profiling followed by uni- and multivariant association tests adjusted for covariates were performed. Qatari individuals with mitochondrial haplogroup J had an increased (twofold) risk of obesity (odds ratio [OR] 1.925; 95% CI 1.234–3.002; P = 0.0038; the Bonferroni adjusted P value threshold is 0.0041), whereas the individuals with haplogroup X were at low risk of obesity (OR 0.387; 95% CI 0.175–0.857; P = 0.019). Further, a set of 38 mitochondrial variants were found to be associated (at P ≤ 0.05) with obesity in models adjusted for age, sex and haplogroup.

Population structure of indigenous inhabitants of Arabia

Modern day Saudi Arabia occupies the majority of historical Arabia, which may have contributed to ancient waves of migration out of Africa. This ancient history has left a lasting imprint in the genetics of the region, including the diverse set of tribes that call Saudi Arabia their home. How these tribes relate to each other and to the world's major populations remains an unanswered question. In an attempt to improve our understanding of the population structure of Saudi Arabia, we conducted genomic profiling of 957 unrelated individuals who self-identify with 28 large tribes in Saudi Arabia. Consistent with the tradition of intra-tribal unions, the subjects showed strong clustering along tribal lines with the distance between clusters correlating with their geographical proximities in Arabia. However, these individuals form a unique cluster when compared to the world's major populations. The ancient origin of these tribal affiliations is supported by analyses that revealed little evidence of ancestral origin from within the 28 tribes. Our results disclose a granular map of population structure and have important implications for future genetic studies into Mendelian and common diseases in the region.

An integrated personal and population-based Egyptian genome reference

A small number of de novo assembled human genomes have been reported to date, and few have been complemented with population-based genetic variation, which is particularly important for North Africa, a region underrepresented in current genome-wide references. Here, we combine long- and short-read whole-genome sequencing data with recent assembly approaches into a de novo assembly of an Egyptian genome. The assembly demonstrates well-balanced quality metrics and is complemented with variant phasing via linked reads into haploblocks, which we associate with gene expression changes in blood. To construct an Egyptian genome reference, we identify genome-wide genetic variation within a cohort of 110 Egyptian individuals. We show that differences in allele frequencies and linkage disequilibrium between Egyptians and Europeans may compromise the transferability of European ancestry-based genetic disease risk and polygenic scores, substantiating the need for multi-ethnic genome references. Thus, the Egyptian genome reference will be a valuable resource for precision medicine.

Whole Genome Sequencing of Four Representatives From the Admixed Population of the United Arab Emirates

Whole genome sequences (WGS) of four nationals of the United Arab Emirates (UAE) at an average coverage of 33X have been completed and described. The selection of suitable subpopulation representatives was informed by a preceding comprehensive population structure analysis. Representatives were chosen based on their central location within the subpopulation on a principal component analysis (PCA) and the degree to which they were admixed. Novel genomic variations among the different subgroups of the UAE population are reported here. Specifically, the WGS analysis identified 4,161,067-4,798,806 variants in the four individual samples, where approximately 80% were single nucleotide polymorphisms (SNPs) and 20% were insertions or deletions (indels). An average of 2.75% was found to be novel variants according to dbSNP (build 151). This is the first report of structural variants (SV) from WGS data from UAE nationals. There were 15,677-20,339 called SVs, of which around 13.5% were novel. The four samples shared 1,399,178 variants, each with distinct variants as follows: 1,085,524 (for the individual denoted as UAE S011), 1,228,559 (UAE S012), 791,072 (UAE S013), and 906,818 (UAE S014). These results show a previously unappreciated population diversity in the region. The synergy of WGS and genotype array data was demonstrated through variant annotation of the former using 2.3 million allele frequencies for the local population derived from the latter technology platform. This novel approach of combining breadth and depth of array and WGS technologies has guided the choice of population genetic representatives and provides complementary, regionalized allele frequency annotation to new genomes comprising millions of loci.

Discovering Novel Biochemical and Genetic Markers for Coronary Heart Disease in Qatari Individuals: The Initiative Qatar Cardiovascular Biorepository

Background:

We aimed to describe the creation and challenge of a DNA and plasma biorepository (Qatar Cardiovascular Biorepository) with linkage to the electronic health record of cardiovascular risk factors to facilitate discovery of novel genetic and proteomic biomarkers for coronary heart disease in Qatari individuals.

Methods:

A prospective case–control study was conducted between October 2013 and February 2018. CHD was defined as a history of an acute coronary syndrome (myocardial infarction [MI]/unstable angina) or coronary revascularization. Controls were identified from blood donors who had no history of coronary heart diesase. After informed consent, blood samples were obtained for DNA and plasma. Demographic, laboratory, and clinical variables were derived from the electronic medical record, and information regarding history of cardiovascular diseases and risk factors was collected from surveys. Challenges in establishing the biorepository were noted, and processes to promote use of the biorepository by Qatari investigators were put in place.

Results:

During the study period, 2671 individuals were approached; of them, 2087 participants were recruited (1029 patients and 1058 controls). Relevant risk factors were ascertained from the electronic health record and surveys. The mean age was 49 ± 16 years, with 61% males. Challenges included setting up the infrastructure for qatar cardiovascular biorepository, developing an informed consent document in Arabic/English, and meeting target recruitment goals. The prevalence of diabetes mellitus, hypertension, dyslipidemia, and smoking was 41%, 44.5%, 40%, and 19%, respectively. History of myocardial infarction, percutaneous coronary intervention, and coronary artery bypass surgery was 55%, 68%, and 17%, respectively, among patients.

Conclusions:

This study addresses the challenges in setting up qatar cardiovascular biorepository, the first cardiovascular genomics biorepository in the Arab Middle Eastern region. QCBio is a unique resource for identifying genetic susceptibility variants and novel circulating markers for coronary heart disease in Qatari adults and enables individualized assessment of risk for coronary heart disease.

Molecular Blood Group Screening in Donors from Arabian Countries and Iran Using High-Throughput MALDI-TOF Mass Spectrometry and PCR-SSP

Background and Aims

Only little is known about blood groups other than ABO blood groups and Rhesus factors in Arabian countries and Iran. During the last years, increased migration to Central Europe has put a focus on the question how to guarantee blood supply for patients from these countries, particularly because hemoglobinopathies with the need of regular blood support are more frequent in patients from that region. Therefore, blood group allele frequencies should be determined in individuals from Arabian countries and Iran by molecular typing and compared to a German rare donor panel.

Methods

1,111 samples including 800 individuals from Syria, 147 from Iran, 123 from the Arabian Peninsula, and 41 from Northern African countries were included in a MALDI-TOF MS assay to detect polymorphisms coding for Kk, Fy(a/b), Fy_null, C_w, Jk(a/b), Jo(a+/a-), Lu(a/b), Lu(8/14), Ss, Do(a/b), Co(a/b), In(a/b), Js(a/b), Kp(a/b), and variant alleles RHCE*c.697C>G and RHCE *c.733C>G. Yt(a/b), S-s-U-, Vel_null, Co_null, and RHCE *c.667G>T were tested by PCR-SSP.

Results

Of the Arabian donors, 2% were homozygous for the FY *02.01N allele (Fy_null), and 15.7% carried the heterozygous mutation. However, 0.8% of the German donors also carried 1 copy of the allele. 3.6% of all and 29.3% of Northern African donors were heterozygous for the RHCE *c.733C>G substitution, 0.4% of the Syrian probands were heterozygous for DO *01/DO *01.-05, a genotype that was lacking in German donors. Whereas the KEL *02.06 allele coding for the Js(a) phenotype was missing in Germans; 0.8% of the Syrian donors carried 1 copy of this allele. 1.8% of the Syrian but only 0.3% of the German donors were negative for YT *01. One donor from Northern Africa homo-zygously carried the GYPB *270+5g>t mutation, inducing the S-s-U+^w phenotype, and in 2 German donors a GYPB *c.161G>A exchange, which induces the Mit+ phenotype, caused a GYPB *03 allele dropout in the MALDI assay. The overall failure rate of the Arabian panel was 0.4%.

Conclusions

Some blood group alleles that are largely lacking in Europeans but had been described in African individuals are present in Arabian populations at a somewhat lower frequency. In single cases, it could be challenging to provide immunized Arabian patients with compatible blood.

Qatari Genotype May Contribute to Complications in Type 2 Diabetes

OBJECTIVE

There is increasing evidence of a strong genetic component in type 2 diabetes (T2DM) that may contribute to diabetes complications. Given the high prevalence of diabetes with its associated complications in the Middle East, we sought to determine if the genotype within a Middle East population may be contributory. Therefore, three genotype-based Qatari ancestral groups, Q1 Arab Bedouin, Q2 Asian/Persian, and Q3 sub-Saharan African, with a fourth admixed group were correlated with T2DM prevalence and its complications to determine if they differed between the 4 Qatari ancestries, particularly for the SLMAP allele-associated diabetic retinopathy.

METHODS

In this cross-sectional study, 398 Qatari subjects, 220 with and 178 without T2DM, were genotyped by Affymetrix 500k SNP arrays. Ancestry was correlated with diabetes complications.

RESULTS

398 subjects were included, the mean age was 49.8 years, and 56.8% were male. The genotype-based ancestry and T2DM prevalence were as follows: 164 (41.2%) with ancestry Q1, 60.4% with T2DM; 149 (37.4%) with ancestry Q2, 49.7% with T2DM; 31 (7.8%) with ancestry Q3, 61.3% with T2DM; and 54 (13.6%) with "admixed" ancestry, 51.9% with T2DM. For patients with diabetes, hypertension (p < 0.035) and retinopathy (p < 0.016) were greater in the Q3 ancestry.

CONCLUSION

These data suggest that the genotype may contribute to complication risk, as exemplified by the increase in hypertension and retinopathy in the Q3 ancestry, though the SLMAP allele was not implicated; however, diabetes prevalence did not differ between the four Qatari ancestries.

Distinct genetic variation and heterogeneity of the Iranian population

Iran, despite its size, geographic location and past cultural influence, has largely been a blind spot for human population genetic studies. With only sparse genetic information on the Iranian population available, we pursued its genome-wide and geographic characterization based on 1021 samples from eleven ethnic groups. We show that Iranians, while close to neighboring populations, present distinct genetic variation consistent with long-standing genetic continuity, harbor high heterogeneity and different levels of consanguinity, fall apart into a cluster of similar groups and several admixed ones and have experienced numerous language adoption events in the past. Our findings render Iran an important source for human genetic variation in Western and Central Asia, will guide adequate study sampling and assist the interpretation of putative disease-implicated genetic variation. Given Iran's internal genetic heterogeneity, future studies will have to consider ethnic affiliations and possible admixture.

Genome-Wide Homozygosity Mapping Reveals Genes Associated With Cognitive Ability in Children From Saudi Arabia

Recent studies of the genetic foundations of cognitive ability rely on large samples (in extreme, hundreds of thousands) of individuals from relatively outbred populations of mostly European ancestry. Hypothesizing that the genetic foundation of cognitive ability depends on the broader population-specific genetic context, we performed a genome-wide association study and homozygosity mapping of cognitive ability estimates obtained through latent variable modeling in a sample of 354 children from a consanguineous population of Saudi Arabia. Approximately half of the sample demonstrated significantly elevated homozygosity levels indicative of inbreeding, and among those with elevated levels, homozygosity was negatively associated with cognitive ability. Further homozygosity mapping identified a specific run, inclusive of the GRIA4 gene, that survived corrections for multiple testing for association with cognitive ability. The results suggest that in a consanguineous population, a notable proportion of the variance in cognitive ability in the normal range in children might be regulated by population-specific mechanisms such as patterns of elevated homozygosity. This observation has implications for the field's understanding of the etiological bases of intelligence and its variability around the world.

The clinical and genetic characteristics of permanent neonatal diabetes (PNDM) in the state of Qatar

Background

Neonatal diabetes mellitus (NDM) is a rare condition that occurs within the first six months of life. Permanent NDM (PNDM) is caused by mutations in specific genes that are known for their expression at early and/or late stages of pancreatic beta‐ cell development, and are either involved in beta‐cell survival, insulin processing, regulation, and release. The native population in Qatar continues to practice consanguineous marriages that lead to a high level of homozygosity. To our knowledge, there is no previous report on the genomics of NDM among the Qatari population. The aims of the current study are to identify patients with NDM diagnosed between 2001 and 2016, and examine their clinical and genetic characteristics.

Methods

To calculate the incidence of PNDM, all patients with PNDM diagnosed between 2001 and 2016 were compared to the total number of live births over the 16‐year‐period. Whole Genome Sequencing (WGS) was used to investigate the genetic etiology in the PNDM cohort.

Results

PNDM was diagnosed in nine (n = 9) patients with an estimated incidence rate of 1:22,938 live births among the indigenous Qatari. Seven different mutations in six genes (PTF1A, GCK, SLC2A2, EIF2AK3, INS, and HNF1B) were identified. In the majority of cases, the genetic etiology was part of a previously identified autosomal recessive disorder. Two novel de novo mutations were identified in INS and HNF1B.

Conclusion

Qatar has the second highest reported incidence of PNDM worldwide. A majority of PNDM cases present as rare familial autosomal recessive disorders. Pancreas associated transcription factor 1a (PTF1A) enhancer deletions are the most common cause of PNDM in Qatar, with only a few previous cases reported in the literature.

Genome-Wide Characterization of Arabian Peninsula Populations: Shedding Light on the History of a Fundamental Bridge between Continents

The Arabian Peninsula (AP) was an important crossroad between Africa, Asia, and Europe, being the cradle of the structure defining these main human population groups, and a continuing path for their admixture. The screening of 741,000 variants in 420 Arabians and 80 Iranians allowed us to quantify the dominant sub-Saharan African admixture in the west of the peninsula, whereas South Asian and Levantine/European influence was stronger in the east, leading to a rift between western and eastern sides of the Peninsula. Dating of the admixture events indicated that Indian Ocean slave trade and Islamization periods were important moments in the genetic makeup of the region. The western-eastern axis was also observable in terms of positive selection of diversity conferring lactose tolerance, with the West AP developing local adaptation and the East AP acquiring the derived allele selected in European populations and existing in South Asia. African selected malaria resistance through the DARC gene was enriched in all Arabian genomes, especially in the western part. Clear European influences associated with skin and eye color were equally frequent across the Peninsula.

Analyses of Neanderthal introgression suggest that Levantine and southern Arabian populations have a shared population history

OBJECTIVES

Modern humans are thought to have interbred with Neanderthals in the Near East soon after modern humans dispersed out of Africa. This introgression event likely took place in either the Levant or southern Arabia depending on the dispersal route out of Africa that was followed. In this study, we compare Neanderthal introgression in contemporary Levantine and southern Arabian populations to investigate Neanderthal introgression and to study Near Eastern population history.

MATERIALS AND METHODS

We analyzed genotyping data on >400,000 autosomal SNPs from seven Levantine and five southern Arabian populations and compared these data to those from populations from around the world including Neanderthal and Denisovan genomes. We used f₄ and D statistics to estimate and compare levels of Neanderthal introgression between Levantine, southern Arabian, and comparative global populations. We also identified 1,581 putative Neanderthal-introgressed SNPs within our dataset and analyzed their allele frequencies as a means to compare introgression patterns in Levantine and southern Arabian genomes.

RESULTS

We find that Levantine and southern Arabian populations have similar levels of Neanderthal introgression to each other but lower levels than other non-Africans. Furthermore, we find that introgressed SNPs have very similar allele frequencies in the Levant and southern Arabia, which indicates that Neanderthal introgression is similarly distributed in Levantine and southern Arabian genomes.

DISCUSSION

We infer that the ancestors of contemporary Levantine and southern Arabian populations received Neanderthal introgression prior to separating from each other and that there has been extensive gene flow between these populations.

Population-scale genomics-Enabling precision public health

The current excitement for affordable genomics technologies and national precision medicine initiatives marks a turning point in worldwide healthcare practices. The last decade of global population sequencing efforts has defined the enormous extent of genetic variation in the human population resulting in insights into differential disease burden and response to therapy within and between populations. Population-scale pharmacogenomics helps to provide insights into the choice of optimal therapies and an opportunity to estimate, predict and minimize adverse events. Such an approach can potentially empower countries to formulate national selection and dosing policies for therapeutic agents thereby promoting public health with precision. We review the breadth and depth of worldwide population-scale sequencing efforts and its implications for the implementation of clinical pharmacogenetics toward making precision medicine a reality.

Variation and Functional Impact of Neanderthal Ancestry in Western Asia

Neanderthals contributed genetic material to modern humans via multiple admixture events. Initial admixture events presumably occurred in Western Asia shortly after humans migrated out of Africa. Despite being a focal point of admixture, earlier studies indicate lower Neanderthal introgression rates in some Western Asian populations as compared with other Eurasian populations. To better understand the genome-wide and phenotypic impact of Neanderthal introgression in the region, we sequenced whole genomes of nine present-day Europeans, Africans, and the Western Asian Druze at high depth, and analyzed available whole genome data from various other populations, including 16 genomes from present-day Turkey. Our results confirmed previous observations that contemporary Western Asian populations, on an average, have lower levels of Neanderthal-introgressed DNA relative to other Eurasian populations. Modern Western Asians also show comparatively high variability in Neanderthal ancestry, which may be attributed to the complex demographic history of the region. We further replicated the previously described depletion of putatively functional sequences among Neanderthal-introgressed haplotypes. Still, we find dozens of common Neanderthal-introgressed haplotypes in the Turkish sample associated with human phenotypes, including anthropometric and metabolic traits, as well as the immune response. One of these haplotypes is unusually long and harbors variants that affect the expression of members of the CCR gene family and are associated with celiac disease. Overall, our results paint a complex first picture of the genomic impact of Neanderthal introgression in the Western Asian populations.

Assessment of coding region variants in Kuwaiti population: implications for medical genetics and population genomics

Consanguineous populations of the Arabian Peninsula have been underrepresented in global efforts that catalogue human exome variability. We sequenced 291 whole exomes of unrelated, healthy native Arab individuals from Kuwait to a median coverage of 45X and characterised 170,508 single-nucleotide variants (SNVs), of which 21.7% were ‘personal’. Up to 12% of the SNVs were novel and 36% were population-specific. Half of the SNVs were rare and 54% were missense variants. The study complemented the Greater Middle East Variome by way of reporting many additional Arabian exome variants. The study corroborated Kuwaiti population genetic substructures previously derived using genome-wide genotype data and illustrated the genetic relatedness among Kuwaiti population subgroups, Middle Eastern, European and Ashkenazi Jewish populations. The study mapped 112 rare and frequent functional variants relating to pharmacogenomics and disorders (recessive and common) to the phenotypic characteristics of Arab population. Comparative allele frequency data and carrier distributions of known Arab mutations for 23 disorders seen among Arabs, of putative OMIM-listed causal mutations for 12 disorders observed among Arabs but not yet characterized for genetic basis in Arabs, and of 17 additional putative mutations for disorders characterized for genetic basis in Arab populations are presented for testing in future Arab studies.

Exome sequencing-based identification of novel type 2 diabetes risk allele loci in the Qatari population

BACKGROUND

Type 2 diabetes (T2D) susceptibility is influenced by genetic and lifestyle factors. To date, the majority of genetic studies of T2D have been in populations of European and Asian descent. The focus of this study is on genetic variations underlying T2D in Qataris, a population with one of the highest incidences of T2D worldwide.

RESULTS

Illumina HiSeq exome sequencing was performed on 864 Qatari subjects (574 T2D cases, 290 controls). Sequence kernel association test (SKAT) gene-based analysis identified an association for low frequency potentially deleterious variants in 6 genes. However, these findings were not replicated by SKAT analysis in an independent cohort of 12,699 exomes, primarly due to the absence of low frequency potentially deleterious variants in 5 of the 6 genes. Interestingly one of the genes identified, catenin beta 1 (CTNNB1, β-catenin), is the key effector of the Wnt pathway and interacts with the nuclear receptor transcription factor 7-like 2 (TCF7L2), variants which are the most strongly associated with risk of developing T2D worldwide. Single variant analysis did not identify any associated variants, suggesting the SKAT association signal was not driven by individual variants. None of the 6 associated genes were among 634 previously described T2D genes.

CONCLUSIONS

The observation that genes not previously linked to T2D in prior studies of European and Asian populations are associated with T2D in Qatar provides new insights into the complexity of T2D pathogenesis and emphasizes the importance of understudied populations when assessing genetic variation in the pathogenesis of common disorders.

Whole-Genome-Sequence-Based Haplotypes Reveal Single Origin of the Sickle Allele during the Holocene Wet Phase

Five classical designations of sickle haplotypes are made on the basis of the presence or absence of restriction sites and are named after the ethno-linguistic groups or geographic regions from which the individuals with sickle cell anemia originated. Each haplotype is thought to represent an independent occurrence of the sickle mutation rs334 (c.20A>T [p.Glu7Val] in HBB). We investigated the origins of the sickle mutation by using whole-genome-sequence data. We identified 156 carriers from the 1000 Genomes Project, the African Genome Variation Project, and Qatar. We classified haplotypes by using 27 polymorphisms in linkage disequilibrium with rs334. Network analysis revealed a common haplotype that differed from the ancestral haplotype only by the derived sickle mutation at rs334 and correlated collectively with the Central African Republic (CAR), Cameroon, and Arabian/Indian haplotypes. Other haplotypes were derived from this haplotype and fell into two clusters, one composed of Senegal haplotypes and the other composed of Benin and Senegal haplotypes. The near-exclusive presence of the original sickle haplotype in the CAR, Kenya, Uganda, and South Africa is consistent with this haplotype predating the Bantu expansions. Modeling of balancing selection indicated that the heterozygote advantage was 15.2%, an equilibrium frequency of 12.0% was reached after 87 generations, and the selective environment predated the mutation. The posterior distribution of the ancestral recombination graph yielded a sickle mutation age of 259 generations, corresponding to 7,300 years ago during the Holocene Wet Phase. These results clarify the origin of the sickle allele and improve and simplify the classification of sickle haplotypes.

Effect of genetic polymorphisms in SREBF-SCAP pathway on therapeutic response to rosuvastatin in Saudi metabolic syndrome patients

AIM

Genetic variants contribute to statins' therapeutic variability. SREBF-SCAP pathway is a key player in lipid homeostasis. Hence, effect of SREBF-SCAP polymorphisms on therapeutic response was studied.

PATIENTS & METHODS

Metabolic syndrome patients of either sex were prescribed rosuvastatin 10 mg for 24 weeks. Clinical, anthropometric and lipid measurements were done before and after treatment. Genotyping was done by pyrosequencing.

RESULTS & CONCLUSION

No associations of SCAP and SREBF-1a genotypes with baseline lipids but significant associations with lipid reductions were observed. Significant effect of SCAP (GG; B = -8.16, p = 0.001); SREBF-1a (GG; B = -7.47, p = 0.001) and SREBF-1a (-delG; B = -7.42, p = 0.001) was observed on total cholesterol reduction. Additive trend was found between SCAP genotypes and lipid reductions. A total of 88% responders have SCAP 'G' allele (p = 0.001). Patients carrying SCAP (GG) and SREBF-1a (GG and -delG) have 9.5-, 8.6- and 14.6-times more likelihood of being responders (p < 0.05). 'G' allele in SCAP and SREBF-1a is significant predictor of rosuvastatin response.

Between Lake Baikal and the Baltic Sea: genomic history of the gateway to Europe

BACKGROUND

The history of human populations occupying the plains and mountain ridges separating Europe from Asia has been eventful, as these natural obstacles were crossed westward by multiple waves of Turkic and Uralic-speaking migrants as well as eastward by Europeans. Unfortunately, the material records of history of this region are not dense enough to reconstruct details of population history. These considerations stimulate growing interest to obtain a genetic picture of the demographic history of migrations and admixture in Northern Eurasia.

RESULTS

We genotyped and analyzed 1076 individuals from 30 populations with geographical coverage spanning from Baltic Sea to Baikal Lake. Our dense sampling allowed us to describe in detail the population structure, provide insight into genomic history of numerous European and Asian populations, and significantly increase quantity of genetic data available for modern populations in region of North Eurasia. Our study doubles the amount of genome-wide profiles available for this region. We detected unusually high amount of shared identical-by-descent (IBD) genomic segments between several Siberian populations, such as Khanty and Ket, providing evidence of genetic relatedness across vast geographic distances and between speakers of different language families. Additionally, we observed excessive IBD sharing between Khanty and Bashkir, a group of Turkic speakers from Southern Urals region. While adding some weight to the "Finno-Ugric" origin of Bashkir, our studies highlighted that the Bashkir genepool lacks the main "core", being a multi-layered amalgamation of Turkic, Ugric, Finnish and Indo-European contributions, which points at intricacy of genetic interface between Turkic and Uralic populations. Comparison of the genetic structure of Siberian ethnicities and the geography of the region they inhabit point at existence of the "Great Siberian Vortex" directing genetic exchanges in populations across the Siberian part of Asia. Slavic speakers of Eastern Europe are, in general, very similar in their genetic composition. Ukrainians, Belarusians and Russians have almost identical proportions of Caucasus and Northern European components and have virtually no Asian influence. We capitalized on wide geographic span of our sampling to address intriguing question about the place of origin of Russian Starovers, an enigmatic Eastern Orthodox Old Believers religious group relocated to Siberia in seventeenth century. A comparative reAdmix analysis, complemented by IBD sharing, placed their roots in the region of the Northern European Plain, occupied by North Russians and Finno-Ugric Komi and Karelian people. Russians from Novosibirsk and Russian Starover exhibit ancestral proportions close to that of European Eastern Slavs, however, they also include between five to 10 % of Central Siberian ancestry, not present at this level in their European counterparts.

CONCLUSIONS

Our project has patched the hole in the genetic map of Eurasia: we demonstrated complexity of genetic structure of Northern Eurasians, existence of East-West and North-South genetic gradients, and assessed different inputs of ancient populations into modern populations.

Homozygous nonsense mutation in SCHIP1/IQCJ-SCHIP1 causes a neurodevelopmental brain malformation syndrome

We report a consanguineous Arab family with 3 affected siblings who display a disorder of global developmental delay, learning difficulties, facial dysmorphism, hearing impairments, and cataract. The clinical phenotype was associated with characteristic brain magnetic resonance imaging (MRI) features of axonal guidance defects involving anterior commissure agenesis as well as scattered areas of polymicrogyria-cobblestone complex. Whole genome sequencing revealed a novel nonsense mutation (159609921C>T) that segregated in the family consistent in an autosomal recessive pattern. This mutation located in the C-terminal region shared by the Schwanomin-Interacting Protein1 (SCHIP1) isoforms including the IQCJ-SCHIP1. The in vitro expression of SCHIP1 and IQCJ-SCHIP1 truncated mutant isoforms (NM_001197109.1; p.R209* and NM_001197114.1; p.R501*, respectively) were markedly reduced as compared to their full-length versions suggesting protein stability/folding impairment. The pathogenic nature of this mutation is supported by a previously reported mouse knockout of Schip1 isoforms, which phenocopied the human axon guidance abnormality. This is the first report of a SCHIP1/IQCJ-SCHIP1 point mutation in humans associated with a neurological-developmental phenotype.