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.

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.

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.

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.

Screening for obesity in the offspring of first-cousin consanguineous couples: A Phase-I study in Saudi Arabia

Consanguineous or cousin marriages are very common in Saudi Arabia. However, owing to limited studies and insufficient knowledge about genetic diseases/disorders, many couples are unaware of the increased health risks for their offspring. Among the inherited and complex diseases from parents’ consanguinity, obesity is common; therefore, we examined the prevalence of obesity in the offspring of first-cousin consanguineous couples in Saudi Arabia. In this questionnaire-based study, 657 individuals (mean age = 18.7 ± 10.2 years; age range = 2–65 years) who were residing in Riyadh, Saudi Arabia participated. Among them, almost 90% were native Saudis. Participants mean body mass index (BMI) was 24.5 ± 9.1 kg/m2. Sex- stratified demographic details confirmed a significant association between age and BMI (p < .001). We confirmed that adolescents and adults were more prone to develop obesity. Adults and non-Saudi participants were three times more likely to develop obesity if they had first-cousin consanguineous parents than those who did not. Of the 30% of participants who were obese, 100 will be selected for Phase II, in which we plan to perform exome sequencing.

Two SNPs in the promoter region of Toll-like receptor 4 gene are not associated with smoking in Saudi Arabia

Defects in the innate immune system, particularly in Toll-like receptors (TLRs), have been reported in several cigarette smoke-promoted diseases. The aim of this study was to examine the impact of tobacco smoke on allelic frequencies of TLR4 single-nucleotide polymorphisms (SNPs) and to compare the genotypic distribution of these SNPs in a Saudi Arabian population with that in previously studied populations. DNA was extracted from 303 saliva samples collected from smokers and nonsmokers. Two transitional SNPs in the promoter region of TLR4 were selected, rs2770150 (T/C) and rs10759931 (G/A). Genotype frequencies were determined using quantitative polymerase chain reaction. Our results showed a slight effect of smoking on the distribution of rs2770150 and rs10759931. However, the differences were not significant. Thus, we conclude that the SNPs selected for this study were independent of smoking and may not be related to smoking-induced diseases.

Association between TLR-9 polymorphisms and colon cancer susceptibility in Saudi Arabian female patients

Objective

The authors aimed to explore the relationship between the expression/polymorphisms of TLR-9 and susceptibility to colon cancer development in the Saudi Arabian population.

Methods

In total, blood samples from 115 patients with colon cancer and 102 participants without colon cancer were analyzed in this study. Three single-nucleotide polymorphisms (SNPs) were selected from the TLR-9 gene, including two sites within the TLR-9 gene’s promoter region (rs352144 and rs187084) and one site in a TLR-9 intron region (rs5743839). Odds ratios (ORs) and 95% confidence intervals (CIs) were computed from logistic regression models after adjusting for age, gender, and tumor localization. To investigate the differential expression of TLR-9 in colon cancer, TLR-9 expression was evaluated using quantitative real-time reverse transcription polymerase chain reaction on 40 matched normal and colon tissues.

Results

The authors found that TLR-9 expression was decreased in colon cancer tissues as compared with that in normal tissues. Moreover, significant associations between the TLR-9 rs187084 SNP and colon cancer risk were observed in female patients only. In rs187084, the T allele had a significantly lower frequency (2.8 times) in female cancer patients than in controls (0.27 vs 0.41). The TLR-9 rs352139 and rs352144 SNPs were significantly associated with colon cancer development when the tumor was located in the rectal area.

Conclusion

The findings support the hypothesis that TLR-9 has an anticancer role in colon cancer development. Furthermore, genetic variation may influence colon cancer development, and SNPs in TLR-9 could serve as biomarkers for decision making in the treatment of females with rectal cancer.

Effect of smoking on the genetic makeup of toll-like receptors 2 and 6

Background

Cigarette smoking is a major risk factor for lung cancer, asthma, and oral cancer, and is central to the altered innate immune responsiveness to infection. Many hypotheses have provided evidence that cigarette smoking induces more genetic changes in genes involved in the development of many cigarette-related diseases. This alteration may be from single-nucleotide polymorphisms (SNPs) in innate immunity genes, especially the toll-like receptors (TLRs).

Objective

In this study, the genotype frequencies of TLR2 and TLR6 in smoking and nonsmoking population were examined.

Methods

Saliva samples were collected from 177 smokers and 126 nonsmokers. The SNPs used were rs3804100 (1350 T/C, Ser450Ser) and rs3804099 (597 T/C, Asn199Asn) for TLR2 and rs3796508 (979 G/A, Val327Met) and rs5743810 (745 T/C, Ser249Pro) for TLR6.

Results

Results showed that TLR2 rs3804100 has a significant effect in short-term smokers (OR =2.63; P=0.04), and this effect is not observed in long-term smokers (>5 years of smoking). Therefore, this early mutation may be repaired by the DNA repair system. For TLR2 rs3804099, the variation in genotype frequencies between the smokers and control patients was due to a late mutation, and its protective role appears only in long-term smokers (OR =0.40, P=0.018). In TLR6 rs5743810, the TT genotype is significantly higher in smokers than in nonsmokers (OR =6.90). The effect of this SNP is observed in long-term smokers, regardless of the smoking regime per day.

Conclusion

TLR2 (rs3804100 and rs3804099) and TLR6 (rs5743810) can be used as a potential index in the diagnosis and prevention of more diseases caused by smoking.

Expression and Polymorphism of Toll-Like Receptor 4 and Effect on NF-κB Mediated Inflammation in Colon Cancer Patients

Our aim was to evaluate the association between the expression and the polymorphism of TLR4/NF-κB pathways and colon cancer. TLR4 (rs4986790, rs10759932, rs10759931 and rs2770150) were genotyped in blood samples from Colorectal patients and healthy controls. TLR4 and cytokines inflammatory expression were evaluated by real time PCR on 40 matching normal and colon tissues and the protein level by Immunohistochemistry. The high level of TLR4 expression in colon cancer tissues is mainly due to infections by bacteria in the human colon and leads to induction of an acute secretion of inflammatory cytokines mediated by NF-κB. Also, we report here a clear evidence for an association between TLR4 rs10759931 polymorphism (OR = 0.086, CI: 0.04–0.18, P = <0.00001). This polymorphism affects the entire population without being specific to either gender or to any age group. In contrast, the rs2770150 is associated with colon cancer in women aged over 50 years and is closely linked with the decreased levels of female sex hormones during the post-menopausal period (OR = 0.188, CI: 0.074–0.48, P = <0.00084). rs10759932 and rs4986790 appear to have any association with colon cancer. Our data suggest that TLR4 SNPs could possibly serve as biomarkers for decision making in colon cancer treatment.