Qatar Genome: Insights on Genomics from the Middle East

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.

Genetic background of primary and familial HLH in Qatar: registry data and population study

Background

Familial hemophagocytic lymphohistiocytosis (FHLH) is an inherited life-threatening disease. Five types are identified, with the addition of congenital immunodeficiency syndromes in which HLH is a typical manifestation. The literature on this disease is very scarce in the Middle East, with only a few scattered reports.

Methods

We report detailed demographic, clinical, and genomic data from 28 patients diagnosed with primary and familial HLH over the last decade in Qatar. An evaluation was performed of allele frequencies of deleterious variants from 12 primary and familial HLH causative genes on the Qatar Genome Programme (QGP) cohort of 14,669 Qatari individuals.

Results

The genetic diagnosis was obtained in 15 patients, and four novel mutations in Perforin 1 (PRF1), UNC13D, LYST, and RAB27A genes were found. We identified 22,945 low/high/moderate/modifier impact variants significantly enriched in the QGP in those 12 genes. The variants rs1271079313 in PRF1 and rs753966933 in RAB27A found in our patient cohort were significantly more prevalent in the QGP compared to the Genome Aggregation Database (gnomAD) database, with a high carrier frequency in the Qatari population.

Conclusions

We established the first primary and familial HLH Registry in the Gulf Region and identified novel possibly pathogenic variants present at higher frequency in the Qatari population, which could be used for screening purposes. Raising awareness about primary and familial HLH and implementing screening activities in the Qatari highly inbred population could stem into more comprehensive premarital and prenatal evaluations and faster diagnosis.

Genome-wide association study and trans-ethnic meta-analysis identify novel susceptibility loci for type 2 diabetes mellitus

BACKGROUND

The genetic basis of type 2 diabetes (T2D) is under-investigated in the Middle East, despite the rapidly growing disease prevalence. We aimed to define the genetic determinants of T2D in Qatar.

METHODS

Using whole genome sequencing of 11,436 participants (2765 T2D cases and 8671 controls) from the population-based Qatar Biobank (QBB), we conducted a genome-wide association study (GWAS) of T2D with and without body mass index (BMI) adjustment.

RESULTS

We replicated 93 known T2D-associated loci in a BMI-unadjusted model, while 96 known loci were replicated in a BMI-adjusted model. The effect sizes and allele frequencies of replicated SNPs in the Qatari population generally concurred with those from European populations. We identified a locus specific to our cohort located between the APOBEC3H and CBX7 genes in the BMI-unadjusted model. Also, we performed a transethnic meta-analysis of our cohort with a previous GWAS on T2D in multi-ancestry individuals (180,834 T2D cases and 1,159,055 controls). One locus in DYNC2H1 gene reached genome-wide significance in the meta-analysis. Assessing polygenic risk scores derived from European- and multi-ancestries in the Qatari population showed higher predictive performance of the multi-ancestry panel compared to the European panel.

CONCLUSION

Our study provides new insights into the genetic architecture of T2D in a Middle Eastern population and identifies genes that may be explored further for their involvement in T2D pathogenesis.

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.

Forging the path to precision medicine in Qatar: a public health perspective on pharmacogenomics initiatives

Pharmacogenomics (PGx) is an important component of precision medicine that promises tailored treatment approaches based on an individual's genetic information. Exploring the initiatives in research that help to integrate PGx test into clinical setting, identifying the potential barriers and challenges as well as planning the future directions, are all important for fruitful PGx implementation in any population. Qatar serves as an exemplar case study for the Middle East, having a small native population compared to a diverse immigrant population, advanced healthcare system, national genome program, and several educational initiatives on PGx and precision medicine. This paper attempts to outline the current state of PGx research and implementation in Qatar within the global context, emphasizing ongoing initiatives and educational efforts. The inclusion of PGx in university curricula and healthcare provider training, alongside precision medicine conferences, showcase Qatar's commitment to advancing this field. However, challenges persist, including the requirement for population specific implementation strategies, complex genetic data interpretation, lack of standardization, and limited awareness. The review suggests policy development for future directions in continued research investment, conducting clinical trials for the feasibility of PGx implementation, ethical considerations, technological advancements, and global collaborations to overcome these barriers.

Cross-population applications of genomics to understand the risk of multifactorial traits involving inflammation and immunity

The interplay between genetic and environmental factors plays a significant role in interindividual variation in immune and inflammatory responses. The availability of high-throughput low-cost genotyping and next-generation sequencing has revolutionized our ability to identify human genetic variation and understand how this varies within and between populations, and the relationship with disease. In this review, we explore the potential of genomics for patient benefit, specifically in the diagnosis, prognosis and treatment of inflammatory and immune-related diseases. We summarize the knowledge arising from genetic and functional genomic approaches, and the opportunity for personalized medicine. The review covers applications in infectious diseases, rare immunodeficiencies and autoimmune diseases, illustrating advances in diagnosis and understanding risk including use of polygenic risk scores. We further explore the application for patient stratification and drug target prioritization. The review highlights a key challenge to the field arising from the lack of sufficient representation of genetically diverse populations in genomic studies. This currently limits the clinical utility of genetic-based diagnostic and risk-based applications in non-Caucasian populations. We highlight current genome projects, initiatives and biobanks from diverse populations and how this is being used to improve healthcare globally by improving our understanding of genetic susceptibility to diseases and regional pathogens such as malaria and tuberculosis. Future directions and opportunities for personalized medicine and wider application of genomics in health care are described, for the benefit of individual patients and populations worldwide.

Genomic architecture of autism spectrum disorder in Qatar: The BARAKA-Qatar Study

BACKGROUND

Autism spectrum disorder (ASD) is a neurodevelopmental condition characterized by impaired social and communication skills, restricted interests, and repetitive behaviors. The prevalence of ASD among children in Qatar was recently estimated to be 1.1%, though the genetic architecture underlying ASD both in Qatar and the greater Middle East has been largely unexplored. Here, we describe the first genomic data release from the BARAKA-Qatar Study-a nationwide program building a broadly consented biorepository of individuals with ASD and their families available for sample and data sharing and multi-omics research.

METHODS

In this first release, we present a comprehensive analysis of whole-genome sequencing (WGS) data of the first 100 families (372 individuals), investigating the genetic architecture, including single-nucleotide variants (SNVs), copy number variants (CNVs), tandem repeat expansions (TREs), as well as mitochondrial DNA variants (mtDNA) segregating with ASD in local families.

RESULTS

Overall, we identify potentially pathogenic variants in known genes or regions in 27 out of 100 families (27%), of which 11 variants (40.7%) were classified as pathogenic or likely-pathogenic based on American College of Medical Genetics (ACMG) guidelines. Dominant variants, including de novo and inherited, contributed to 15 (55.6%) of these families, consisting of SNVs/indels (66.7%), CNVs (13.3%), TREs (13.3%), and mtDNA variants (6.7%). Moreover, homozygous variants were found in 7 families (25.9%), with a sixfold increase in homozygous burden in consanguineous versus non-consanguineous families (13.6% and 1.8%, respectively). Furthermore, 28 novel ASD candidate genes were identified in 20 families, 23 of which had recurrent hits in MSSNG and SSC cohorts.

CONCLUSIONS

This study illustrates the value of ASD studies in under-represented populations and the importance of WGS as a comprehensive tool for establishing a molecular diagnosis for families with ASD. Moreover, it uncovers a significant role for recessive variation in ASD architecture in consanguineous settings and provides a unique resource of Middle Eastern genomes for future research to the global ASD community.

Genome-wide association study identifies several loci for HEV seropositivity

Hepatitis E viral (HEV) infection imposes a heavy global health burden. The variability in the prevalence of serological markers of HEV infection between different ethnic groups proposes a host genetic influence. Here, we report genetic polymorphisms associated with anti-HEV antibody positivity and level using binary- and quantitative-trait genome-wide association studies (GWAS) on a population from Qatar (n = 5829). We identified a region in 12p11.1 (lead SNP: rs559856097, allele: A, p = 2.3 × 10-10) significantly associated with anti-HEV antibodies level. This intergenic variant is located near SNORD112, a small nucleolar RNA (snoRNA). Additional gene-set and pathway enrichment analyses highlighted a strong correlation with anti-viral response-related pathways, including IFNs (alpha/beta) and interleukin-21 (IL-21) signaling. This is the first GWAS on the response to HEV infection. Further replication and functional experimentation are warranted to validate these findings.

Genetics of Inborn Errors of Immunity in highly consanguineous Middle Eastern and North African populations

Consanguineous marriages in Middle Eastern and North African (MENA) countries are deeply-rooted tradition and highly prevalent resulting into increased prevalence of autosomal recessive diseases including Inborn Errors of Immunity (IEIs). Molecular genetic testing is an important diagnostic tool for IEIs since it provides a definite diagnosis, genotype-phenotype correlation, and guide therapy. In this review, we will discuss the current state and challenges of genomic and variome studies in MENA region populations, as well as the importance of funding advanced genome projects. In addition, we will review the MENA underlying molecular genetic defects of over 2457 patients published with the common IEIs, where autosomal recessive mode of inheritance accounts for 76% of cases with increased prevalence of combined immunodeficiency diseases (50%). The efforts made in the last three decades in terms of international collaboration and of in situ capacity building in MENA region countries led to the discovery of more than 150 novel genes involved in IEIs. Expanding sequencing studies within the MENA will undoubtedly be a unique asset for the IEI genetics which can advance research, and support precise genomic diagnostics and therapeutics.

A GHRHR founder mutation causes isolated growth hormone deficiency type IV in a consanguineous Pakistani family

Background

Isolated growth hormone deficiency (IGHD) is caused by a severe shortage or absence of growth hormone (GH), which results in aberrant growth and development. Patients with IGHD type IV (IGHD4) have a short stature, reduced serum GH levels, and delayed bone age.

Objectives

To identify the causative mutation of IGHD in a consanguineous family comprising four affected patients with IGHD4 (MIM#618157) and explore its functional impact in silico.

Methods

Clinical and radiological studies were performed to determine the phenotypic spectrum and hormonal profile of the disease, while whole-exome sequencing (WES) and Sanger sequencing were performed to identify the disease-causing mutation. In-silico studies involved protein structural modeling and docking, and molecular dynamic simulation analyses using computational tools. Finally, data from the Qatar Genome Program (QGP) were screened for the presence of the founder variant in the Qatari population.

Results

All affected individuals presented with a short stature without gross skeletal anomalies and significantly reduced serum GH levels. Genetic mapping revealed a homozygous nonsense mutation [NM_000823:c.G214T:p.(Glu72*)] in the third exon of the growth-hormone-releasing hormone receptor gene GHRHR (MIM#139191) that was segregated in all patients. The substituted amber codon (UAG) seems to truncate the protein by deleting the C-terminus GPCR domain, thus markedly disturbing the GHRHR receptor and its interaction with the growth hormone-releasing hormone.

Conclusion

These data support that a p.Glu72* founder mutation in GHRHR perturbs growth hormone signaling and causes IGHD type IV. In-silico and biochemical analyses support the pathogenic effect of this nonsense mutation, while our comprehensive phenotype and hormonal profiling has established the genotype-phenotype correlation. Based on the current study, early detection of GHRHR may help in better therapeutic intervention.

Framework for Adoption of Next-Generation Sequencing (NGS) Globally in the Oncology Area

Radical new possibilities of improved treatment of cancer are on offer from an advanced medical technology already demonstrating its significance: next-generation sequencing (NGS). This refined testing provides unprecedentedly precise diagnoses and permits the use of focused and highly personalized treatments. However, across regions globally, many cancer patients will continue to be denied the benefits of NGS as long as some of the yawning gaps in its implementation remain unattended. The challenges at the regional and national levels are linked because putting the solutions into effect is highly dependent on cooperation between regional- and national-level cooperation, which could be hindered by shortfalls in interpretation or understanding. The aim of the paper was to define and explore the necessary conditions for NGS and make recommendations for effective implementation based on extensive exchanges with policy makers and stakeholders. As a result, the European Alliance for Personalised Medicine (EAPM) developed a maturity framework structured around demand-side and supply-side issues to enable interested stakeholders in different countries to self-evaluate according to a common matrix. A questionnaire was designed to identify the current status of NGS implementation, and it was submitted to different experts in different institutions globally. This revealed significant variability in the different aspects of NGS uptake. Within different regions globally, to ensure those conditions are right, this can be improved by linking efforts made at the national level, where patients have needs and where care is delivered, and at the global level, where major policy initiatives in the health field are underway or in preparation, many of which offer direct or indirect pathways for building those conditions. In addition, in a period when consensus is still incomplete and catching up is needed at a political level to ensure rational allocation of resources-even within individual countries-to enable the best ways to make the necessary provisions for NGS, a key recommendation is to examine where closer links between national and regional actions could complement, support, and mutually reinforce efforts to improve the situation for patients.

Genetic predisposition to porto-sinusoidal vascular disorder: A functional genomic-based, multigenerational family study

BACKGROUND AND AIMS

Porto-sinusoidal vascular disorder (PSVD) is a group of liver vascular diseases featuring lesions encompassing the portal venules and sinusoids unaccompanied by cirrhosis, irrespective of the presence/absence of portal hypertension. It can occur secondary to coagulation disorders or insult by toxic agents. However, the cause of PSVD remains unknown in most cases. Hereditary cases of PSVD are exceptionally rare, but they are of particular interest and may unveil genetic alterations and molecular mechanisms associated with the disease.

APPROACH AND RESULTS

We performed genome sequencing of four patients and two healthy individuals of a large multigenerational Lebanese family with PSVD and identified a heterozygous deleterious variant (c.547C>T, p.R183W) of FCH and double SH3 domains 1 ( FCHSD1 ), an uncharacterized gene, in patients. This variant segregated with the disease, and its pattern of inheritance was suggestive of autosomal dominant with variable expressivity. RNA structural modelling of human FCHSD1 suggests that the C-to-T substitution at position 547, corresponding to FCHSD1R183W , may increase both messenger RNA (mRNA) and protein stability and its interaction with MTOR-associated protein, LST8 homolog, a key protein of the mechanistic target of rapamycin (mTOR pathway). These predictions were substantiated by biochemical analyses, which showed that FCHSD1R183W induced high FCHSD1 mRNA stability, overexpression of FCHSD1 protein, and an increase in mTORC1 activation. This human FCHSD1 variant was introduced into mice through CRISPR/Cas9 genome editing. Nine out of the 15 mice carrying the human FCHSD1R183W variant mimicked the phenotype of human PSVD, including splenomegaly and enlarged portal vein.

CONCLUSIONS

Aberrant FCHSD1 structure and function leads to mTOR pathway overactivation and may cause PSVD.

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.

Microfluidic systems in clinical diagnosis

The use of microfluidic devices is highly attractive in the field of biomedical and clinical assessments, as their portability and fast response time have become crucial in providing opportune therapeutic treatments to patients. The applications of microfluidics in clinical diagnosis and point-of-care devices are continuously growing. The present review article discusses three main fields where miniaturized devices are successfully employed in clinical applications. The quantification of ions, sugars, and small metabolites is examined considering the analysis of bodily fluids samples and the quantification of this type of analytes employing real-time wearable devices. The discussion covers the level of maturity that the devices have reached as well as cost-effectiveness. The analysis of proteins with clinical relevance is presented and organized by the function of the proteins. The last section covers devices that can perform single-cell metabolomic and proteomic assessments. Each section discusses several strategically selected recent reports on microfluidic devices successfully employed for clinical assessments, to provide the reader with a wide overview of the plethora of novel systems and microdevices developed in the last 5 years. In each section, the novel aspects and main contributions of each reviewed report are highlighted. Finally, the conclusions and future outlook section present a summary and speculate on the future direction of the field of miniaturized devices for clinical applications.

The Genetic Spectrum of Maturity-Onset Diabetes of the Young (MODY) in Qatar, a Population-Based Study

Maturity-onset diabetes of the young (MODY) is a rare monogenic form of diabetes mellitus. In this study, we estimated the prevalence and genetic spectrum of MODY in the Middle Eastern population of Qatar using whole-genome sequencing (WGS) of 14,364 subjects from the population-based Qatar biobank (QBB) cohort. We focused our investigations on 14 previously identified genes ascribed to the cause of MODY and two potentially novel MODY-causing genes, RFX6 and NKX6-1. Genetic variations within the 16 MODY-related genes were assessed for their pathogenicity to identify disease-causing mutations. Analysis of QBB phenotype data revealed 72 subjects (0.5%) with type 1 diabetes, 2915 subjects (20.3%) with type 2 diabetes and 11,377 (79.2%) without diabetes. We identified 22 mutations in 67 subjects that were previously reported in the Human Genetic Mutation Database (HGMD) as disease-causing (DM) or likely disease causing (DM?) for MODY. We also identified 28 potentially novel MODY-causing mutations, predicted to be among the top 1% most deleterious mutations in the human genome, which showed complete (100%) disease penetrance in 34 subjects. Overall, we estimated that MODY accounts for around 2.2-3.4% of diabetes patients in Qatar. This is the first population-based study to determine the genetic spectrum and estimate the prevalence of MODY in the Middle East. Further research to characterize the newly identified mutations is warranted.

Overcoming Eurocentric bias makes for better science

To understand disease, scientists are producing comprehensive omics datasets. However, the majority of these are Eurocentric. Recently, the inclusion of patients from Asia and the Middle East in genomic analyses uncovered unique loci linked to COVID-19 severity. This demonstrates that focusing on diversity and underrepresented populations can benefit all.

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.

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.

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.

The Prevalence and Genetic Spectrum of Familial Hypercholesterolemia in Qatar Based on Whole Genome Sequencing of 14,000 Subjects

Familial hypercholesterolemia (FH) is an inherited disease characterized by reduced efficiency of low-density lipoprotein-cholesterol (LDL-C) removal from the blood and, consequently, an increased risk of life-threatening early cardiovascular complications. In Qatar, the prevalence of FH has not been determined and the disease, as in many countries, is largely underdiagnosed. In this study, we combined whole-genome sequencing data from the Qatar Genome Program with deep phenotype data from Qatar Biobank for 14,056 subjects to determine the genetic spectrum and estimate the prevalence of FH in Qatar. We used the Dutch Lipid Clinic Network (DLCN) as a diagnostic tool and scrutinized 11 FH-related genes for known pathogenic and possibly pathogenic mutations. Results revealed an estimated prevalence of 0.8% (1:125) for definite/probable cases of FH in the Qatari population. We detected 16 known pathogenic/likely pathogenic mutations in LDLR and one in PCSK9; all in a heterozygous state with high penetrance. The most common mutation was rs1064793799 (c.313+3A >C) followed by rs771019366 (p.Asp90Gly); both in LDLR. In addition, we identified 18 highly penetrant possibly pathogenic variants, of which 5 were Qatari-specific, in LDLR, APOB, PCSK9 and APOE, which are predicted to be among the top 1% most deleterious mutations in the human genome but further validations are required to confirm their pathogenicity. We did not detect any homozygous FH or autosomal recessive mutations in our study cohort. This pioneering study provides a reliable estimate of FH prevalence in Qatar based on a significantly large population-based cohort, whilst uncovering the spectrum of genetic variants associated with FH.

Patterns and distribution of de novo mutations in multiplex Middle Eastern families

While de novo mutations (DNMs) are key to genetic diversity, they are also responsible for a high number of rare disorders. To date, no study has systematically examined the rate and distribution of DNMs in multiplex families in highly consanguineous populations. Leveraging WGS profiles of 645 individuals in 146 families, we implemented a combinatorial approach using 3 complementary tools for DNM discovery in 353 unique trio combinations. We found a total of 27,168 DNMs (median: 70 single-nucleotide and 6 insertion-deletions per individual). Phasing revealed around 80% of DNMs were paternal in origin. Notably, using whole-genome methylation data of spermatogonial stem cells, these DNMs were significantly more likely to occur at highly methylated CpGs (OR: 2.03; p value = 6.62 × 10⁻¹¹). We then examined the effects of consanguinity and ethnicity on DNMs, and found that consanguinity does not seem to correlate with DNM rate, and special attention has to be considered while measuring such a correlation. Additionally, we found that Middle-Eastern families with Arab ancestry had fewer DNMs than African families, although not significant (p value = 0.16). Finally, for families with diseased probands, we examined the difference in DNM counts and putative impact across affected and unaffected siblings, but did not find significant differences between disease groups, likely owing to the enrichment for recessive disorders in this part of the world, or the small sample size per clinical condition. This study serves as a reference for DNM discovery in multiplex families from the globally under-represented populations of the Middle-East.

A de novo start-loss in EFTUD2 associated with mandibulofacial dysostosis with microcephaly: case report

Mandibulofacial dysostosis with microcephaly (MFDM) is a rare genetic disorder inherited in an autosomal dominant pattern. Major characteristics include developmental delay, craniofacial malformations such as malar and mandibular hypoplasia, and ear anomalies. Here, we report a 4.5-yr-old female patient with symptoms fitting MFDM. Using whole-genome sequencing, we identified a de novo start-codon loss (c.3G > T) in the EFTUD2. We examined EFTUD2 expression in the patient by RNA sequencing and observed a notable functional consequence of the variant on gene expression in the patient. We identified a novel variant for the development of MFDM in humans. To the best of our knowledge, this is the first report of a start-codon loss in EFTUD2 associated with MFDM.

Dispatches from Biotech beginning BeginNGS: Rapid newborn genome sequencing to end the diagnostic and therapeutic odyssey

In this Dispatch from Biotech, we briefly review the urgent need for extensive expansion of newborn screening (NBS) by genomic sequencing, and the reasons why early attempts had limited success. During the next decade transformative developments will continue in society and in the pharmaceutical, biotechnology, informatics, and medical sectors that enable prompt addition of genetic disorders to NBS by rapid whole genome sequencing (rWGS) upon introduction of new therapies that qualify them according to the Wilson and Jungner criteria (Wilson, J. M. G., & Jungner, G., World Health Organization. (1968). Principles and Practice of Screening for Disease. World Health Organization. Retrieved from https://apps.who.int/iris/handle/10665/37650). Herein we describe plans, progress, and clinical trial designs for BeginNGS (Newborn Genome Sequencing to end the diagnostic and therapeutic odyssey), a new international, pre-competitive, public-private consortium that proposes to implement a self-learning healthcare delivery system for screening all newborns for over 400 hundred genetic diseases, diagnostic confirmation, implementation of effective treatment, and acceleration of orphan drug development. We invite investigators and stakeholders worldwide to join the consortium in a prospective, multi-center, international trial of the clinical utility and cost effectiveness of BeginNGS.

Ratios of Acetaminophen Metabolites Identify New Loci of Pharmacogenetic Relevance in a Genome-Wide Association Study

Genome-wide association studies (GWAS) with non-targeted metabolomics have identified many genetic loci of biomedical interest. However, metabolites with a high degree of missingness, such as drug metabolites and xenobiotics, are often excluded from such studies due to a lack of statistical power and higher uncertainty in their quantification. Here we propose ratios between related drug metabolites as GWAS phenotypes that can drastically increase power to detect genetic associations between pairs of biochemically related molecules. As a proof-of-concept we conducted a GWAS with 520 individuals from the Qatar Biobank for who at least five of the nine available acetaminophen metabolites have been detected. We identified compelling evidence for genetic variance in acetaminophen glucuronidation and methylation by UGT2A15 and COMT, respectively. Based on the metabolite ratio association profiles of these two loci we hypothesized the chemical structure of one of their products or substrates as being 3-methoxyacetaminophen, which we then confirmed experimentally. Taken together, our study suggests a novel approach to analyze metabolites with a high degree of missingness in a GWAS setting with ratios, and it also demonstrates how pharmacological pathways can be mapped out using non-targeted metabolomics measurements in large population-based studies.

Functional Characterization of the MYO6 Variant p.E60Q in Non-Syndromic Hearing Loss Patients

Hereditary hearing loss (HHL) is a common genetic disorder accounting for at least 60% of pre-lingual deafness in children, of which 70% is inherited in an autosomal recessive pattern. The long tradition of consanguinity among the Qatari population has increased the prevalence of HHL, which negatively impacts the quality of life. Here, we functionally validated the pathogenicity of the c.178G>C, p.E60Q mutation in the MYO6 gene, which was detected previously in a Qatari HHL family, using cellular and animal models. In vitro analysis was conducted in HeLa cells transiently transfected with plasmids carrying MYO6WT or MYO6p.E60Q, and a zebrafish model was generated to characterize the in vivo phenotype. Cells transfected with MYO6WT showed higher expression of MYO6 in the plasma membrane and increased ATPase activity. Modeling the human MYO6 variants in zebrafish resulted in severe otic defects. At 72 h post-injection, MYO6p.E60Q embryos demonstrated alterations in the sizes of the saccule and utricle. Additionally, zebrafish with MYO6p.E60Q displayed super-coiled and bent hair bundles in otic hair cells when compared to control and MYO6WT embryos. In conclusion, our cellular and animal models add support to the in silico prediction that the p.E60Q missense variant is pathogenic and damaging to the protein. Since the c.178G>C MYO6 variant has a 0.5% allele frequency in the Qatari population, about 400 times higher than in other populations, it could contribute to explaining the high prevalence of hearing impairment in Qatar.