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.

A Complex Intrachromosomal Rearrangement Disrupting IRF6 in a Family with Popliteal Pterygium and Van der Woude Syndromes

Clefts of the lip and/or palate (CL/P) are considered the most common form of congenital anomalies occurring either in isolation or in association with other clinical features. Van der woude syndrome (VWS) is associated with about 2% of all CL/P cases and is further characterized by having lower lip pits. Popliteal pterygium syndrome (PPS) is a more severe form of VWS, normally characterized by orofacial clefts, lower lip pits, skin webbing, skeletal anomalies and syndactyly of toes and fingers. Both syndromes are inherited in an autosomal dominant manner, usually caused by heterozygous mutations in the Interferon Regulatory Factor 6 (IRF6) gene. Here we report the case of a two-generation family where the index presented with popliteal pterygium syndrome while both the father and sister had clinical features of van der woude syndrome, but without any point mutations detected by re-sequencing of known gene panels or microarray testing. Using whole genome sequencing (WGS) followed by local de novo assembly, we discover and validate a copy-neutral, 429 kb complex intra-chromosomal rearrangement in the long arm of chromosome 1, disrupting the IRF6 gene. This variant is copy-neutral, novel against publicly available databases, and segregates in the family in an autosomal dominant pattern. This finding suggests that missing heritability in rare diseases may be due to complex genomic rearrangements that can be resolved by WGS and de novo assembly, helping deliver answers to patients where no genetic etiology was identified by other means.

A recessive variant in SIM2 in a child with complex craniofacial anomalies and global developmental delay

Rare deletions and duplications on the long arm of Chromosome 21 have previously been reported in many patients with craniofacial and developmental phenotypes. However, this Down Syndrome Critical Region (DSCR) contains multiple genes, making identifying a single causative gene difficult. Here, we report a case of a boy with bicoronal craniosynostosis, facial dysmorphism, developmental delay, and intellectual impairment who was found by whole genome sequencing to have a homozygous missense mutation in the Single-Minded Homolog 2 (SIM2) gene (c.461 A  >  G, p.Tyr154Cys) within the DSCR. SIM2 encodes an essential bHLH and PAS domain transcription factor expressed during fetal brain development and acts as a master regulator of neurogenesis. This variant is globally very rare, segregates in the family, and is predicted to be highly deleterious by in silico analysis, 3D molecular modeling of protein structure, and functional analysis of zebrafish models. Zebrafish expressing the human SIM2^p.Y154C variant displayed a progressed microcephaly-like phenotype and head shape abnormalities. When combined with careful phenotyping of the patient vis-à-vis previously reported cases harboring structural variants in this critical 21q22 region, the data support a pathogenic role of SIM2 in this complex syndrome and demonstrates the utility of next-generation sequencing in prioritizing genes in contiguous deletions/duplications syndromes and diagnosing microarray-negative patients in the craniofacial clinic.