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

Near-complete Middle Eastern genomes refine autozygosity and enhance disease-causing and population-specific variant discovery

Advances in long-read sequencing have enabled routine complete assembly of human genomes, but much remains to be done to represent broader populations and show impact on disease-gene discovery. Here, we report highly accurate, near-complete and phased genomes from six Middle Eastern (ME) family trios (n = 18) with neurodevelopmental conditions, representing ancestries from Sudan, Jordan, Syria, Qatar and Afghanistan. These genomes revealed 42.2 Mb of new sequence (13.8% impacting known genes), 75 new HLA/KIR alleles and strong signals of inbreeding, with ROH covering up to one-third of chromosomes 6 and 12 in one individual. Using assembly-based variant calling, we identified 23 de novo and recessive variants as strong candidates for causing previously unresolved symptoms in the probands. The ME genomes revealed unique variation relative to existing references, showing enhanced mappability and variant calling. These results underscore the value of de novo assembly for disease variant discovery and the need for sampled ME-specific references to better characterize population-relevant variation.

Clinical Features and Predictors of Unfavorable Outcomes in Infantile Epileptic Spasms Syndrome: Results from a Multiancestry Cohort

Background: Infantile epileptic spasms syndrome is an epileptic encephalopathy with poor prognosis. Objective: To describe clinical features and predictors of unfavorable outcomes of infantile epileptic spasms syndrome in a multiancestry cohort. Methods: Retrospective cross-sectional chart review analysis at a single tertiary pediatric hospital in 56 patients with infantile epileptic spasms syndrome. Results: Fifty-nine percent were males. Most patients (57%) were from the Middle East North African region. The mean age at onset of spasms was 6.6 months (±3.5 mo). At presentation, 80% of patients had a preexisting neurologic disorder, an abnormal development (66%), and/or an abnormal neurologic examination (73%). Hypsarrhythmia on electroencephalography (EEG) was reported in 55% of patients. Etiology was identified in 75%. Abnormalities on neuroimaging were common (70%), and 54% had abnormal genetic testing. The preferred first-line treatment was vigabatrin (66%) followed by high-dose oral steroids (16%). Only a minority (∼30%) experienced complete resolution of their epileptic spasms/hypsarrhythmia by 2 weeks. At the last follow-up, 84% of patients had an abnormal development and 40% were still on antiseizure medications. Comorbidities were common (75%). The majority of patients had unfavorable outcomes, including abnormal development (84%) and/or drug-resistant epilepsy (45%). Preexisting abnormal development/neurologic examination or neurologic conditions, perinatal risk factors, neuroimaging abnormalities, and an identified etiology were associated with increased odds ratios for unfavorable developmental outcomes or drug-resistant epilepsy. Significance: Clinical features and outcomes of infantile epileptic spasms syndrome in this multiancestry cohort confirms findings of previous epidemiologic studies, recognizes some differences, and highlights risk factors that can be used for early identification of patients with predicted unfavorable outcomes.

Genomics of rare diseases in the Greater Middle East

The Greater Middle East (GME) represents a concentrated region of unparalleled genetic diversity, characterized by an abundance of distinct alleles, founder mutations and extensive autozygosity driven by high consanguinity rates. These genetic hallmarks present a unique, yet vastly untapped resource for genomic research on Mendelian diseases. Despite this immense potential, the GME continues to face substantial challenges in comprehensive data collection and analysis. This Perspective highlights the region's unique position as a natural laboratory for genetic discovery and explores the challenges that have stifled progress thus far. Importantly, we propose strategic solutions, advocating for an all-inclusive research approach. With targeted investment and focused efforts, the latent genetic wealth in the GME can be transformed into a global hub for genomic research that will redefine and advance our understanding of the human genome.

Multiethnic prevalence of the APOL1 G1 and G2 variants among the Israeli dialysis population

Background and hypothesis

The two apolipoprotein L1 (APOL1) variants, G1 and G2, are common in populations of sub-Saharan African ancestry. Individuals with two of these alleles (G1 or G2) have an increased risk for a spectrum of non-diabetic chronic kidney diseases. However, these variants are typically not observed outside of populations that self-identify as current continental Africans or having clear recent African ancestry such as, most notably, African Americans, and other large population groups in the Americas and several European countries. We hypothesized that the diverse ethnic groups within the Israeli population may exhibit varying levels of recent African ancestry. Therefore, it is plausible that APOL1 risk alleles might be present even in individuals who do not self-identify as being of sub-Saharan African descent.

Methods

We non-selectively screened people with kidney failure across Israel for APOL1 risk variants using restriction fragment length polymorphism.

Results

We recruited 1744 individuals from 38 dialysis units in Israel. We identified eight patients of Moroccan Jewish, Bedouin, or Muslim Arab ancestry, who carry at least one G1 or G2 allele. None of the eight patients carried the protective APOL1 p.N264K variant. Furthermore, despite all Bedouin individuals being G2 heterozygous, the G2 minor allele frequency was significantly enriched in kidney failure cases compared to ethnically matched controls (P = .006).

Conclusions

These findings show that APOL1 G1 and G2 allelic variants are present in populations previously not appreciated to possess recent sub-Saharan ancestry and suggest that a single G2 risk variant may confer increased risk for chronic kidney disease in certain population contexts.

Patterns of population structure and genetic variation within the Saudi Arabian population

The Arabian Peninsula is considered the initial site of historic human migration out of Africa. The modern-day indigenous Arabians are believed to be the descendants who remained from the ancient split of the migrants into Eurasia. Here, we investigated how the population history and cultural practices such as endogamy have shaped the genetic variation of the Saudi Arabians. We genotyped 3,352 individuals and identified twelve genetic sub-clusters that corresponded to the geographical distribution of different tribal regions, differentiated by distinct components of ancestry based on comparisons to modern and ancient DNA references. These sub-clusters also showed variation across ranges of the genome covered in runs of homozygosity, as well as differences in population size changes over time. Using 25,488,981 variants found in whole genome sequencing data (WGS) from 302 individuals, we found that the Saudi tend to show proportionally more deleterious alleles than neutral alleles when compared to Africans/African Americans from gnomAD (e.g. a 13% increase of deleterious alleles annotated by AlphaMissense between 0.5 - 5% frequency in Saudi, compared to 7% decrease of the benign alleles; P < 0.001). Saudi sub-clusters with greater inbreeding and lower effective population sizes showed greater enrichment of deleterious alleles as well. Additionally, we found that approximately 10% of the variants discovered in our WGS data are not observed in gnomAD; these variants are also enriched with deleterious annotations. To accelerate studying the population-enriched deleterious alleles and their health consequences in this population, we made available the allele frequency estimates of 25,488,981 variants discovered in our samples. Taken together, our results suggest that Saudi's population history impacts its pattern of genetic variation with potential consequences to the population health. It further highlights the need to sequence diverse and unique populations so to provide a foundation on which to interpret medical- and pharmaco- genomic findings from these populations.

Genetic Variant Analyses Identify Novel Candidate Autism Risk Genes from a Highly Consanguineous Cohort of 104 Families from Oman

Deficits in social communication, restricted interests, and repetitive behaviours are hallmarks of autism spectrum disorder (ASD). Despite high genetic heritability, the majority of clinically diagnosed ASD cases have unknown genetic origins. We performed genome sequencing on mothers, fathers, and affected individuals from 104 families with ASD in Oman, a Middle Eastern country underrepresented in international genetic studies. This approach identified 48 novel candidate genes significantly associated with ASD in Oman. In particular, 35 of these genes have been previously implicated in neurodevelopmental disorders (NDDs) in other populations, underscoring the conserved genetic basis of ASD across ethnicities. Genetic variants within these candidate genes that would impact the encoded protein included 1 insertion, 4 frameshift, 6 splicing, 12 nonsense, and 67 missense changes. Notably, 61% of the SNVs were homozygous, suggesting a prominent recessive genetic architecture for ASD in this unique population. The scarcity of genetic studies on ASD in the Arabian Peninsula has impeded the understanding of the unique genetic landscape of ASD in this region. These findings help bridge this knowledge gap and provide valuable insights into the complex genetic basis of ASD in Oman.

Genome Sequencing Identifies 13 Novel Candidate Risk Genes for Autism Spectrum Disorder in a Qatari Cohort

Autism spectrum disorder (ASD) is a neurodevelopmental condition characterized by deficits in social communication, restricted interests, and repetitive behaviors. Despite considerable research efforts, the genetic complexity of ASD remains poorly understood, complicating diagnosis and treatment, especially in the Arab population, with its genetic diversity linked to migration, tribal structures, and high consanguinity. To address the scarcity of ASD genetic data in the Middle East, we conducted genome sequencing (GS) on 50 ASD subjects and their unaffected parents. Our analysis revealed 37 single-nucleotide variants from 36 candidate genes and over 200 CGG repeats in the FMR1 gene in one subject. The identified variants were classified as uncertain, likely pathogenic, or pathogenic based on in-silico algorithms and ACMG criteria. Notably, 52% of the identified variants were homozygous, indicating a recessive genetic architecture to ASD in this population. This finding underscores the significant impact of high consanguinity within the Qatari population, which could be utilized in genetic counseling/screening program in Qatar. We also discovered single nucleotide variants in 13 novel genes not previously associated with ASD: ARSF, BAHD1, CHST7, CUL2, FRMPD3, KCNC4, LFNG, RGS4, RNF133, SCRN2, SLC12A8, USP24, and ZNF746. Our investigation categorized the candidate genes into seven groups, highlighting their roles in cognitive development, including the ubiquitin pathway, transcription factors, solute carriers, kinases, glutamate receptors, chromatin remodelers, and ion channels.

The impact of consanguinity on the design of iPSC banks

The effect of consanguinity on identifying universal induced pluripotent stem cell (iPSC) donors, i.e., homozygous for the major human leukocyte antigen (HLA) loci, is unknown. The discovery sample size was calculated in a consanguineous population using a method (1qF) based on the inbreeding coefficient. The result was orders of magnitude smaller compared to the standard method.

Qatar's genetic counseling landscape: Current insights and future prospects

Genetic counseling is a dynamic and rapidly growing field. In Qatar, the significance of genetic counseling is underscored by the distinctive demographic characteristics of the population, including elevated rates of consanguinity and larger family sizes, contributing to the increased incidence of many genetic conditions. This emphasizes the crucial role of genetic counseling in addressing the specific needs of the community. Over the past decade, key health care institutions in Qatar, such as Hamad Medical Corporation and Sidra Medicine, have significantly expanded genetic counseling services encompassing premarital, reproductive, prenatal, pediatric, adult, and cancer care. This multifaceted approach reflects Qatar's health care system's commitment to addressing various aspects of genetic health and well-being across different life stages. A pivotal milestone in the field's development in Qatar was the establishment of a genetic counseling master's program at Qatar University in 2018, showcasing the country's dedication to fostering indigenous expertise in genetic counseling with the necessary competencies and cultural sensitivity to address the unique genetic counseling needs of the population. The recognition of genetic counseling as a profession and the licensure by the Ministry of Public Health in Qatar is another key achievement to ensure the high quality of service and protection of the profession. Contributing to global genetic knowledge, various academic and research entities in Qatar are conducting genetic/genomic/genetic counseling research toward advancing precision medicine in the country, and initiatives such as the Qatar Biobank and Qatar Genome Program have played a major role in catalyzing these efforts.

QPGx-CARES: Qatar pharmacogenetics clinical applications and research enhancement strategies

Pharmacogenetic (PGx)-informed medication prescription is a cutting-edge genomic application in contemporary medicine, offering the potential to overcome the conventional "trial-and-error" approach in drug prescription. The ability to use an individual's genetic profile to predict drug responses allows for personalized drug and dosage selection, thereby enhancing the safety and efficacy of treatments. However, despite significant scientific and clinical advancements in PGx, its integration into routine healthcare practices remains limited. To address this gap, the Qatar Genome Program (QGP) has embarked on an ambitious initiative known as QPGx-CARES (Qatar Pharmacogenetics Clinical Applications and Research Enhancement Strategies), which aims to set a roadmap for optimizing PGx research and clinical implementation on a national scale. The goal of QPGx-CARES initiative is to integrate PGx testing into clinical settings with the aim of improving patient health outcomes. In 2022, QGP initiated several implementation projects in various clinical settings. These projects aimed to evaluate the clinical utility of PGx testing, gather valuable insights into the effective dissemination of PGx data to healthcare professionals and patients, and identify the gaps and the challenges for wider adoption. QPGx-CARES strategy aimed to integrate evidence-based PGx findings into clinical practice, focusing on implementing PGx testing for cardiovascular medications, supported by robust scientific evidence. The current initiative sets a precedent for the nationwide implementation of precision medicine across diverse clinical domains.

Founder mutations and rare disease in the Arab world

Founder mutations are disease-causing variants that occur frequently in geographically or culturally isolated groups whose shared ancestor(s) carried the pathogenic variant. While some disease alleles may vanish from the genetic pool due to natural selection, variants with weaker effects may survive for a long time, thereby enhancing the prevalence of some rare diseases. These are predominantly autosomal recessive diseases but can also be autosomal dominant traits with late-onset or mild phenotypes. Cultural practices, such as endogamy and consanguinity, in these isolated groups lead to higher prevalence of such rare diseases compared to the rest of the population and worldwide. In this Perspective, we define population isolates and the underlying genetic mechanisms for accumulating founder mutations. We also discuss the current and potential scientific, clinical and public-health implications of studying founder mutations in population isolates around the world, with a particular focus on the Arab population.

Imputation accuracy across global human populations

Genotype imputation is now fundamental for genome-wide association studies but lacks fairness due to the underrepresentation of references from non-European ancestries. The state-of-the-art imputation reference panel released by the Trans-Omics for Precision Medicine (TOPMed) initiative improved the imputation of admixed African-ancestry and Hispanic/Latino samples, but imputation for populations primarily residing outside of North America may still fall short in performance due to persisting underrepresentation. To illustrate this point, we imputed the genotypes of over 43,000 individuals across 123 populations around the world and identified numerous populations where imputation accuracy paled in comparison to that of European-ancestry populations. For instance, the mean imputation r-squared (Rsq) for variants with minor allele frequencies between 1% and 5% in Saudi Arabians (n = 1,061), Vietnamese (n = 1,264), Thai (n = 2,435), and Papua New Guineans (n = 776) were 0.79, 0.78, 0.76, and 0.62, respectively, compared to 0.90-0.93 for comparable European populations matched in sample size and SNP array content. Outside of Africa and Latin America, Rsq appeared to decrease as genetic distances to European-ancestry reference increased, as predicted. Using sequencing data as ground truth, we also showed that Rsq may over-estimate imputation accuracy for non-European populations more than European populations, suggesting further disparity in accuracy between populations. Using 1,496 sequenced individuals from Taiwan Biobank as a second reference panel to TOPMed, we also assessed a strategy to improve imputation for non-European populations with meta-imputation, but this design did not improve accuracy across frequency spectra. Taken together, our analyses suggest that we must ultimately strive to increase diversity and size to promote equity within genetics research.

Burden of Mendelian disorders in a large Middle Eastern biobank

BACKGROUND

Genome sequencing of large biobanks from under-represented ancestries provides a valuable resource for the interrogation of Mendelian disease burden at world population level, complementing small-scale familial studies.

METHODS

Here, we interrogate 6045 whole genomes from Qatar-a Middle Eastern population with high consanguinity and understudied mutational burden-enrolled at the national Biobank and phenotyped for 58 clinically-relevant quantitative traits. We examine a curated set of 2648 Mendelian genes from 20 panels, annotating known and novel pathogenic variants and assessing their penetrance and impact on the measured traits.

RESULTS

We find that 62.5% of participants are carriers of at least 1 known pathogenic variant relating to recessive conditions, with homozygosity observed in 1 in 150 subjects (0.6%) for which Peninsular Arabs are particularly enriched versus other ancestries (5.8-fold). On average, 52.3 loss-of-function variants were found per genome, 6.5 of which affect a known Mendelian gene. Several variants annotated in ClinVar/HGMD as pathogenic appeared at intermediate frequencies in this cohort (1-3%), highlighting Arab founder effect, while others have exceedingly high frequencies (> 5%) prompting reconsideration as benign. Furthermore, cumulative gene burden analysis revealed 56 genes having gene carrier frequency > 1/50, including 5 ACMG Tier 3 panel genes which would be candidates for adding to newborn screening in the country. Additionally, leveraging 58 biobank traits, we systematically assess the impact of novel/rare variants on phenotypes and discover 39 candidate large-effect variants associating with extreme quantitative traits. Furthermore, through rare variant burden testing, we discover 13 genes with high mutational load, including 5 with impact on traits relevant to disease conditions, including metabolic disorder and type 2 diabetes, consistent with the high prevalence of these conditions in the region.

CONCLUSIONS

This study on the first phase of the growing Qatar Genome Program cohort provides a comprehensive resource from a Middle Eastern population to understand the global mutational burden in Mendelian genes and their impact on traits in seemingly healthy individuals in high consanguinity settings.

Utilizing Pharmacogenomic Data for a Safer Use of Statins among the Emirati Population

BACKGROUND

Statins are the most prescribed lipid-lowering drugs worldwide. The associated adverse events, especially muscle symptoms, have been frequently reported despite their perceived safety. Three pharmacogenes, the solute carrier organic anion transporter family member 1B1 (SLCO1B1), ATP-binding cassette subfamily G member 2 (ABCG2), and cytochrome P450 2C9 (CYP2C9) are suggested as safety biomarkers for statins. The Clinical Pharmacogenomic Implementation Consortium (CPIC) issued clinical guidelines for statin use based on these three genes.

OBJECTIVES

The present study aimed to examine variants in these pharmacogenes to predict the safety of statin use among the Emirati population.

METHODS

Analyzing 242 whole exome sequencing data at the three genes enabled the determination of the frequencies of the single nucleotide polymorphisms (SNPs), annotating the haplotypes and the predicted functions of their proteins.

RESULTS

In our cohort, 29.8% and 5.4% had SLCO1B1 decreased and poor function, respectively. The high frequency warns of the possibility of significant side effects of some statins and the importance of pharmacogenomic testing. We found a low frequency (6%) of the ABCG2:rs2231142 variant, which indicates the low probability of Emirati patients being recommended against higher rosuvastatin doses compared with other populations with higher frequencies of this variant. In contrast, we found high frequencies of the functionally impaired CYP2C9 alleles, which makes fluvastatin a less favorable choice.

CONCLUSION

Among the sparse studies available, the present one demonstrates all SLCO1B1 and CYP2C9 function-impairing alleles among Emiratis. We highlighted how population-specific pharmacogenomic data can predict safer choices of statins, especially in understudied populations.

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.

Identification of PCSK9-like human gene knockouts using metabolomics, proteomics, and whole-genome sequencing in a consanguineous population

Natural human knockouts of genes associated with desirable outcomes, such as PCSK9 with low levels of LDL-cholesterol, can lead to the discovery of new drug targets and treatments. Rare loss-of-function variants are more likely to be found in the homozygous state in consanguineous populations, and deep molecular phenotyping of blood samples from homozygous carriers can help to discriminate between silent and functional variants. Here, we combined whole-genome sequencing with proteomics and metabolomics for 2,935 individuals from the Qatar Biobank (QBB) to evaluate the power of this approach for finding genes of clinical and pharmaceutical interest. As proof-of-concept, we identified a homozygous carrier of a very rare PCSK9 variant with extremely low circulating PCSK9 levels and low LDL. Our study demonstrates that the chances of finding such variants are about 168 times higher in QBB compared with GnomAD and emphasizes the potential of consanguineous populations for drug discovery.

Assessing the genetic burden of familial hypercholesterolemia in a large middle eastern biobank

BACKGROUND

The genetic architecture underlying Familial Hypercholesterolemia (FH) in Middle Eastern Arabs is yet to be fully described, and approaches to assess this from population-wide biobanks are important for public health planning and personalized medicine.

METHODS

We evaluate the pilot phase cohort (n = 6,140 adults) of the Qatar Biobank (QBB) for FH using the Dutch Lipid Clinic Network (DLCN) criteria, followed by an in-depth characterization of all genetic alleles in known dominant (LDLR, APOB, and PCSK9) and recessive (LDLRAP1, ABCG5, ABCG8, and LIPA) FH-causing genes derived from whole-genome sequencing (WGS). We also investigate the utility of a globally established 12-SNP polygenic risk score to predict FH individuals in this cohort with Arab ancestry.

RESULTS

Using DLCN criteria, we identify eight (0.1%) 'definite', 41 (0.7%) 'probable' and 334 (5.4%) 'possible' FH individuals, estimating a prevalence of 'definite or probable' FH in the Qatari cohort of ~ 1:125. We identify ten previously known pathogenic single-nucleotide variants (SNVs) and 14 putatively novel SNVs, as well as one novel copy number variant in PCSK9. Further, despite the modest sample size, we identify one homozygote for a known pathogenic variant (ABCG8, p. Gly574Arg, global MAF = 4.49E-05) associated with Sitosterolemia 2. Finally, calculation of polygenic risk scores found that individuals with 'definite or probable' FH have a significantly higher LDL-C SNP score than 'unlikely' individuals (p = 0.0003), demonstrating its utility in Arab populations.

CONCLUSION

We design and implement a standardized approach to phenotyping a population biobank for FH risk followed by systematically identifying known variants and assessing putative novel variants contributing to FH burden in Qatar. Our results motivate similar studies in population-level biobanks - especially those with globally under-represented ancestries - and highlight the importance of genetic screening programs for early detection and management of individuals with high FH risk in health systems.

Ancestry-related distribution of Runs of homozygosity and functional variants in Qatari population

Background

Describing how genetic history shapes the pattern of medically relevant variants could improve the understanding of how specific loci interact with each other and affect diseases and traits prevalence. The Qatari population is characterized by a complex history of admixture and substructure, and the study of its population genomic features would provide valuable insights into the genetic landscape of functional variants. Here, we analyzed the genomic variation of 186 newly-genotyped healthy individuals from the Qatari peninsula.

Results

We discovered an intricate genetic structure using ancestry related analyses. In particular, the presence of three different clusters, Cluster 1, Cluster 2 and Cluster 3 (with Near Eastern, South Asian and African ancestry, respectively), was detected with an additional fourth one (Cluster 4) with East Asian ancestry. These subpopulations show differences in the distribution of runs of homozygosity (ROH) and admixture events in the past, ranging from 40 to 5 generations ago. This complex genetic history led to a peculiar pattern of functional markers under positive selection, differentiated in shared signals and private signals. Interestingly we found several signatures of shared selection on SNPs in the FADS2 gene, hinting at a possible common evolutionary link to dietary intake. Among the private signals, we found enrichment for markers associated with HDL and LDL for Cluster 1(Near Eastern ancestry) and Cluster 3 (South Asian ancestry) and height and blood traits for Cluster 2 (African ancestry).

The differences in genetic history among these populations also resulted in the different frequency distribution of putative loss of function variants. For example, homozygous carriers for rs2884737, a variant linked to an anticoagulant drug (warfarin) response, are mainly represented by individuals with predominant Bedouin ancestry (risk allele frequency G at 0.48).

Conclusions

We provided a detailed catalogue of the different ancestral pattern in the Qatari population highlighting differences and similarities in the distribution of selected variants and putative loss of functions. Finally, these results would provide useful guidance for assessing genetic risk factors linked to consanguinity and genetic ancestry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01087-1.

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.

Investigation of Genetic Causes in Patients with Congenital Heart Disease in Qatar: Findings from the Sidra Cardiac Registry

Congenital heart disease (CHD) is one of the most common forms of birth defects worldwide, with a prevalence of 1–2% in newborns. CHD is a multifactorial disease partially caused by genetic defects, including chromosomal abnormalities and single gene mutations. Here, we describe the Sidra Cardiac Registry, which includes 52 families and a total of 178 individuals, and investigate the genetic etiology of CHD in Qatar. We reviewed the results of genetic tests conducted in patients as part of their clinical evaluation, including chromosomal testing. We also performed whole exome sequencing (WES) to identify potential causative variants. Sixteen patients with CHD had chromosomal abnormalities that explained their complex CHD phenotype, including six patients with trisomy 21. Moreover, using exome analysis, we identified potential CHD variants in 24 patients, revealing 65 potential variants in 56 genes. Four variants were classified as pathogenic/likely pathogenic based on the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) classification; these variants were detected in four patients. This study sheds light on several potential genetic variants contributing to the development of CHD. Additional functional studies are needed to better understand the role of the identified variants in the pathogenesis of CHD.

Middle Eastern Genetic Variation Improves Clinical Annotation of the Human Genome

Genetic variation in populations of Middle Eastern origin remains highly underrepresented in most comprehensive genomic databases. This underrepresentation hampers the functional annotation of the human genome and challenges accurate clinical variant interpretation. To highlight the importance of capturing genetic variation in the Middle East, we aggregated whole exome and genome sequencing data from 2116 individuals in the Middle East and established the Middle East Variation (MEV) database. Of the high-impact coding (missense and loss of function) variants in this database, 53% were absent from the most comprehensive Genome Aggregation Database (gnomAD), thus representing a unique Middle Eastern variation dataset which might directly impact clinical variant interpretation. We highlight 39 variants with minor allele frequency >1% in the MEV database that were previously reported as rare disease variants in ClinVar and the Human Gene Mutation Database (HGMD). Furthermore, the MEV database consisted of 281 putative homozygous loss of function (LoF) variants, or complete knockouts, of which 31.7% (89/281) were absent from gnomAD. This set represents either complete knockouts of 83 unique genes in reportedly healthy individuals, with implications regarding disease penetrance and expressivity, or might affect dispensable exons, thus refining the clinical annotation of those regions. Intriguingly, 24 of those genes have several clinically significant variants reported in ClinVar and/or HGMD. Our study shows that genetic variation in the Middle East improves functional annotation and clinical interpretation of the genome and emphasizes the need for expanding sequencing studies in the Middle East and other underrepresented populations.

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.

Poking COVID-19: Insights on Genomic Constraints among Immune-Related Genes between Qatari and Italian Populations

Host genomic information, specifically genomic variations, may characterize susceptibility to disease and identify people with a higher risk of harm, leading to better targeting of care and vaccination. Italy was the epicentre for the spread of COVID-19 in Europe, the first country to go into a national lockdown and has one of the highest COVID-19 associated mortality rates. Qatar, on the other hand has a very low mortality rate. In this study, we compared whole-genome sequencing data of 14398 adults and Qatari-national to 925 Italian individuals. We also included in the comparison whole-exome sequence data from 189 Italian laboratory-confirmed COVID-19 cases. We focused our study on a curated list of 3619 candidate genes involved in innate immunity and host-pathogen interaction. Two population-gene metric scores, the Delta Singleton-Cohort variant score (DSC) and Sum Singleton-Cohort variant score (SSC), were applied to estimate the presence of selective constraints in the Qatari population and in the Italian cohorts. Results based on DSC and SSC metrics demonstrated a different selective pressure on three genes (MUC5AC, ABCA7, FLNA) between Qatari and Italian populations. This study highlighted the genetic differences between Qatari and Italian populations and identified a subset of genes involved in innate immunity and host-pathogen interaction.