Characterization of Greater Middle Eastern genetic variation for enhanced disease gene discovery

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.

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.

Virtual Gene Panels Have a Superior Diagnostic Yield for Inherited Rare Diseases Relative to Static Panels

BACKGROUND

Exome- or genome-based panels-also known as slices or virtual panels-are now a popular approach that involves comprehensive genomic sequencing while restricting analysis to subsets of genes based on patients' phenotypes. This flexible strategy enables frequent gene updates based on novel disease associations as well as reflexing to analyzing other genes up to the whole exome or genome. With recent improvements addressing limitations associated with virtual panels, the advantages of this approach, relative to static custom-based panels, remain to be systematically characterized.

METHODS

Here we perform slice testing on 1014 patients (50.5% females; average age 17 years) referred from multiple pediatric clinics within a single center in the Middle East (83% Arab population).

RESULTS

Initial analysis uncovered molecular diagnoses for 235 patients for a diagnostic yield of 23% (235/1014). "On the fly" focused analysis in most negative cases (N = 779) identified clinically significant variants correlating with patients' presentations in genes outside the originally ordered panel for another 35 patients (3.5% or 35/1024) increasing the overall diagnostic yield to 27%. The pathogenic variants underlying the additional cases (13% of all positive cases) were excluded from the original "panel" gene list, mainly as result of issues related to panel selection, novel gene-disease associations, phenotype spectrum broadening, or gene lists variability. The additional findings led to changes in clinical management in most patients (94%).

CONCLUSIONS

Our findings support slice testing as an efficient and flexible platform that facilitates updates to gene lists to achieve high clinical sensitivity and utility.

Somatic genetic alterations of oncogenes and tumor suppressor genes in cholesteatoma

BACKGROUND

Cholesteatoma is a proliferative disease that affects the tympanic cavity and temporal bone. Despite many studies and various theories, the etiopathogenesis of cholesteatoma has not been fully elucidated. Features such as invasion, migration, uncontrolled proliferation, and lack of differentiation are observed in both cholesteatoma and neoplasia.

AIMS/OBJECTIVES

The aim of this study is to investigate somatic genetic alterations in known proto-oncogenes and tumor suppressor genes in cholesteatoma.

MATERIAL AND METHODS

60 different known proto-oncogenes and tumor suppressor genes were comparatively analyzed in cholesteatoma and peripheric blood samples from 15 middle ear cholesteatoma patients using next-generation sequencing.

RESULTS

JAK3 c.2164G > A, TP53 c.284delC, and KRAS c.377A > T alterations were observed in cholesteatoma tissue but not in normal tissue. In addition, 12 different germline variants were also identified in 8 patients.

CONCLUSIONS AND SIGNIFICANCE

In this study, the presence of changes in cancer-related genes in cholesteatoma was determined and these changes were discussed in terms of possible clinical applications. We hope that the genetic alterations that emerged in this study, will be beneficial in guiding future research in this field.

Ensemble and consensus approaches to prediction of recessive inheritance for missense variants in human disease

Mode of inheritance (MOI) is necessary for clinical interpretation of pathogenic variants; however, the majority of variants lack this information. Furthermore, variant effect predictors are fundamentally insensitive to recessive-acting diseases. Here, we present MOI-Pred, a variant pathogenicity prediction tool that accounts for MOI, and ConMOI, a consensus method that integrates variant MOI predictions from three independent tools. MOI-Pred integrates evolutionary and functional annotations to produce variant-level predictions that are sensitive to both dominant-acting and recessive-acting pathogenic variants. Both MOI-Pred and ConMOI show state-of-the-art performance on standard benchmarks. Importantly, dominant and recessive predictions from both tools are enriched in individuals with pathogenic variants for dominant- and recessive-acting diseases, respectively, in a real-world electronic health record (EHR)-based validation approach of 29,981 individuals. ConMOI outperforms its component methods in benchmarking and validation, demonstrating the value of consensus among multiple prediction methods. Predictions for all possible missense variants are provided in the "Data and code availability" section.

The genetic landscape of autism spectrum disorder in an ancestrally diverse cohort

Autism spectrum disorder (ASD) comprises neurodevelopmental disorders with wide variability in genetic causes and phenotypes, making it challenging to pinpoint causal genes. We performed whole exome sequencing on a modest, ancestrally diverse cohort of 195 families, including 754 individuals (222 with ASD), and identified 38,834 novel private variants. In 68 individuals with ASD (~30%), we identified 92 potentially pathogenic variants in 73 known genes, including BCORL1, CDKL5, CHAMP1, KAT6A, MECP2, and SETD1B. Additionally, we identified 158 potentially pathogenic variants in 120 candidate genes, including DLG3, GABRQ, KALRN, KCTD16, and SLC8A3. We also found 34 copy number variants in 31 individuals overlapping known ASD loci. Our work expands the catalog of ASD genetics by identifying hundreds of variants across diverse ancestral backgrounds, highlighting convergence on nervous system development and signal transduction. These findings provide insights into the genetic underpinnings of ASD and inform molecular diagnosis and potential therapeutic targets.

Understanding the variant landscape, and genetic epidemiology of Multiple Endocrine Neoplasia in India

PURPOSE

Multiple Endocrine Neoplasia (MEN) is a group of familial cancer syndromes that encompasses several types of endocrine tumors differentiated by genetic mutations in RET, MEN1 and CDKN1B genes. Accurate diagnosis of MEN subtypes can thus be performed through genetic testing. However, MEN variants remain largely understudied in Indian populations. Additionally, few dedicated resources to understand these disorders currently exist.

METHODS

Using the gold-standard ACMG/AMP guidelines, we systematically classified variants reported across the three genes in the IndiGen dataset, and established the genetic epidemiology of MEN in the Indian population. We further classified ClinVar and Mastermind variants and compiled all into a database. Finally, we designed a multiplex primer panel for rapid variant identification.

RESULTS

We have established the genetic prevalence of MEN as the following: 1 in 1026 individuals is likely to be afflicted with MEN linked with pathogenic RET mutations. We have further created the MAPVar database containing 3280 ACMG-classified variants freely accessible at: https://clingen.igib.res.in/MAPVar/ . Finally, our NGS primer panel covers 33 exonic regions across two pools through 38 amplicons with a total amplified region of 65 kb.

CONCLUSION

Our work establishes that MEN is a prevalent disorder in India. The rare nature of Indian variants underscores the need of genomic and functional studies to establish a more comprehensive variant landscape. Additionally, our panel offers a means of cost-effective genetic testing, and the MAPVar database a ready reference to aid in a better understanding of variant pathogenicity in clinical as well as research settings.

A reference quality, fully annotated diploid genome from a Saudi individual

We have used multiple sequencing approaches to sequence the genome of a volunteer from Saudi Arabia. We use the resulting data to generate a de novo assembly of the genome, and use different computational approaches to refine the assembly. As a consequence, we provide a contiguous assembly of the complete genome of an individual from Saudi Arabia for all chromosomes except chromosome Y, and label this assembly KSA001. We transferred genome annotations from reference genomes to fully annotate KSA001, and we make all primary sequencing data, the assembly, and the genome annotations freely available in public databases using the FAIR data principles. KSA001 is the first telomere-to-telomere-assembled genome from a Saudi individual that is freely available for any purpose.

Genetic heterogeneity in familial forms of genetic generalized epilepsy: from mono- to oligogenism

Genetic generalized epilepsy (GGE) including childhood absence epilepsy, juvenile absence epilepsy, juvenile myoclonic epilepsy (JME), and GGE with tonic-clonic seizures (TCS) (GGE-TCS), is genetically influenced with a two- to four- fold increased risk in the first-degree relatives of patients. Since large families with GGE are very rare, international studies have focused on sporadic GGE patients using whole exome sequencing, suggesting that GGE are highly genetically heterogeneous and rather involve rare or ultra-rare variants. Moreover, a polygenic mode of inheritance is suspected in most cases. We performed SNP microarrays and whole exome sequencing in 20 families from Sudan, focusing on those with at least four affected members. Standard genetic filters and Endeavour algorithm for functional prioritization of genes selected likely susceptibility variants in FAT1, DCHS1 or ASTN2 genes. FAT1 and DCHS1 are adhesion transmembrane proteins interacting during brain development, while ASTN2 is involved in dendrite development. Our approach on familial forms of GGE is complementary to large-scale collaborative consortia studies of sporadic cases. Our study reinforces the hypothesis that GGE is genetically heterogeneous, even in a relatively limited geographic area, and mainly oligogenic, as supported by the low familial penetrance of GGE and by the Bayesian algorithm that we developed in a large pedigree with JME. Since populations with founder effect and endogamy are appropriate to study autosomal recessive pathologies, they would be also adapted to decipher genetic components of complex diseases, using the reported bayesian model.

Arab founder variants: Contributions to clinical genomics and precision medicine

BACKGROUND

Founder variants are ancestral variants shared by individuals who are not closely related. The large effect size of some of these variants in the context of Mendelian disorders offers numerous precision medicine opportunities.

METHODS

Using one of the largest datasets on Mendelian disorders in the Middle East, we identified 2,908 medically relevant founder variants derived from 18,360 exomes and genomes and investigated their contribution to the clinical annotation of the human genome.

FINDINGS

Strikingly, ∼34% of Arab founder variants are absent in gnomAD. We found a strong contribution of Arab founder variants to the identification of novel gene-disease links (n = 224) and the support/dispute (n = 81 support, n = 101 dispute) of previously reported candidate gene-disease links. The powerful segregation evidence generated by Arab founder variants allowed many ClinVar and Human Gene Mutation Database variants to be reclassified. Overall, 39.5% of diagnostic reports from our clinical lab are based on founder variants, and 19.41% of tested individuals carry at least one pathogenic founder variant. The presumptive loss-of-function mechanism that typically underlies autosomal recessive diseases means that Arab founder variants also offer unique opportunities in "druggable genome" research. Arab founder variants were also informative of migration patterns in the Middle East consistent with documented historical accounts.

CONCLUSIONS

We highlight the contribution of founder variants from an under-represented population group to precision medicine and inform future prevention programs. Our study also sheds light on the added value of these variants in supplementing other lines of research in tracing population history.

FUNDING

There is no funding for this work.

A novel missense variant in the ATPase domain of ATP8A2 and review of phenotypic variability of ATP8A2-related disorders caused by missense changes

ATPase, class 1, type 8 A, member 2 (ATP8A2) is a P4-ATPase with a critical role in phospholipid translocation across the plasma membrane. Pathogenic variants in ATP8A2 are known to cause cerebellar ataxia, impaired intellectual development, and disequilibrium syndrome 4 (CAMRQ4) which is often associated with encephalopathy, global developmental delay, and severe motor deficits. Here, we present a family with two siblings born from a consanguineous, first-cousin union from Sudan presenting with global developmental delay, intellectual disability, spasticity, ataxia, nystagmus, and thin corpus callosum. Whole exome sequencing revealed a homozygous missense variant in the nucleotide binding domain of ATP8A2 (p.Leu538Pro) that results in near complete loss of protein expression. This is in line with other missense variants in the same domain leading to protein misfolding and loss of ATPase function. In addition, by performing diffusion-weighted imaging, we identified bilateral hyperintensities in the posterior limbs of the internal capsule suggesting possible microstructural changes in axon tracts that had not been appreciated before and could contribute to the sensorimotor deficits in these individuals.

Overcoming barriers to single-cell RNA sequencing adoption in low- and middle-income countries

The advent of single-cell resolution sequencing and spatial transcriptomics has enabled the delivery of cellular and molecular atlases of tissues and organs, providing new insights into tissue health and disease. However, if the full potential of these technologies is to be equitably realised, ancestrally inclusivity is paramount. Such a goal requires greater inclusion of both researchers and donors in low- and middle-income countries (LMICs). In this perspective, we describe the current landscape of ancestral inclusivity in genomic and single-cell transcriptomic studies. We discuss the collaborative efforts needed to scale the barriers to establishing, expanding, and adopting single-cell sequencing research in LMICs and to enable globally impactful outcomes of these technologies.

Germline mutations in a G protein identify signaling cross-talk in T cells

Humans with monogenic inborn errors responsible for extreme disease phenotypes can reveal essential physiological pathways. We investigated germline mutations in GNAI2, which encodes Gαi2, a key component in heterotrimeric G protein signal transduction usually thought to regulate adenylyl cyclase-mediated cyclic adenosine monophosphate (cAMP) production. Patients with activating Gαi2 mutations had clinical presentations that included impaired immunity. Mutant Gαi2 impaired cell migration and augmented responses to T cell receptor (TCR) stimulation. We found that mutant Gαi2 influenced TCR signaling by sequestering the guanosine triphosphatase (GTPase)-activating protein RASA2, thereby promoting RAS activation and increasing downstream extracellular signal-regulated kinase (ERK)/mitogen-activated protein kinase (MAPK) and phosphatidylinositol 3-kinase (PI3K)-AKT S6 signaling to drive cellular growth and proliferation.

NOTCH3 p.Arg1231Cys is markedly enriched in South Asians and associated with stroke

The genetic factors of stroke in South Asians are largely unexplored. Exome-wide sequencing and association analysis (ExWAS) in 75 K Pakistanis identified NM_000435.3(NOTCH3):c.3691 C > T, encoding the missense amino acid substitution p.Arg1231Cys, enriched in South Asians (alternate allele frequency = 0.58% compared to 0.019% in Western Europeans), and associated with subcortical hemorrhagic stroke [odds ratio (OR) = 3.39, 95% confidence interval (CI) = [2.26, 5.10], p = 3.87 × 10-9), and all strokes (OR [CI] = 2.30 [1.77, 3.01], p = 7.79 × 10-10). NOTCH3 p.Arg231Cys was strongly associated with white matter hyperintensity on MRI in United Kingdom Biobank (UKB) participants (effect [95% CI] in SD units = 1.1 [0.61, 1.5], p = 3.0 × 10-6). The variant is attributable for approximately 2.0% of hemorrhagic strokes and 1.1% of all strokes in South Asians. These findings highlight the value of diversity in genetic studies and have major implications for genomic medicine and therapeutic development in South Asian populations.

Genome Tunisia Project: paving the way for precision medicine in North Africa

BACKGROUND

Key discoveries and innovations in the field of human genetics have led to the foundation of molecular and personalized medicine. Here, we present the Genome Tunisia Project, a two-phased initiative (2022-2035) which aims to deliver the reference sequence of the Tunisian Genome and to support the implementation of personalized medicine in Tunisia, a North African country that represents a central hub of population admixture and human migration between African, European, and Asian populations. The main goal of this initiative is to develop a healthcare system capable of incorporating omics data for use in routine medical practice, enabling medical doctors to better prevent, diagnose, and treat patients.

METHODS

A multidisciplinary partnership involving Tunisian experts from different institutions has come to discern all requirements that would be of high priority to fulfill the project's goals. One of the most urgent priorities is to determine the reference sequence of the Tunisian Genome. In addition, extensive situation analysis and revision of the education programs, community awareness, appropriate infrastructure including sequencing platforms and biobanking, as well as ethical and regulatory frameworks, have been undertaken towards building sufficient capacity to integrate personalized medicine into the Tunisian healthcare system.

RESULTS

In the framework of this project, an ecosystem with all engaged stakeholders has been implemented including healthcare providers, clinicians, researchers, pharmacists, bioinformaticians, industry, policymakers, and advocacy groups. This initiative will also help to reinforce research and innovation capacities in the field of genomics and to strengthen discoverability in the health sector.

CONCLUSIONS

Genome Tunisia is the first initiative in North Africa that seeks to demonstrate the major impact that can be achieved by Human Genome Projects in low- and middle-income countries to strengthen research and to improve disease management and treatment outcomes, thereby reducing the social and economic burden on healthcare systems. Sharing this experience within the African scientific community is a chance to turn a major challenge into an opportunity for dissemination and outreach. Additional efforts are now being made to advance personalized medicine in patient care by educating consumers and providers, accelerating research and innovation, and supporting necessary changes in policy and regulation.

IFIH1 loss of function predisposes to inflammatory and SARS-CoV-2-related infectious diseases

The IFIH1 gene, encoding melanoma differentiation-associated protein 5 (MDA5), is an indispensable innate immune regulator involved in the early detection of viral infections. Previous studies described MDA5 dysregulation in weakened immunological responses, and increased susceptibility to microbial infections and autoimmune disorders. Monoallelic gain-of-function of the IFIH1 gene has been associated with multisystem disorders, namely Aicardi-Goutieres and Singleton-Merten syndromes, while biallelic loss causes immunodeficiency. In this study, nine patients suffering from recurrent infections, inflammatory diseases, severe COVID-19 or multisystem inflammatory syndrome in children (MIS-C) were identified with putative loss-of-function IFIH1 variants by whole-exome sequencing. All patients revealed signs of lymphopaenia and an increase in inflammatory markers, including CRP, amyloid A, ferritin and IL-6. One patient with a pathogenic homozygous variant c.2807+1G>A was the most severe case showing immunodeficiency and glomerulonephritis. The c.1641+1G>C variant was identified in the heterozygous state in patients suffering from periodic fever, COVID-19 or MIS-C, while the c.2016delA variant was identified in two patients with inflammatory bowel disease or MIS-C. There was a significant association between IFIH1 monoallelic loss of function and susceptibility to infections in males. Expression analysis showed that PBMCs of one patient with a c.2016delA variant had a significant decrease in ISG15, IFNA and IFNG transcript levels, compared to normal PBMCs, upon stimulation with Poly(I:C), suggesting that MDA5 receptor truncation disrupts the immune response. Our findings accentuate the implication of rare monogenic IFIH1 loss-of-function variants in altering the immune response, and severely predisposing patients to inflammatory and infectious diseases, including SARS-CoV-2-related disorders.

Pathogenic heterozygous TRPM7 variants and hypomagnesemia with developmental delay

Background

Heterozygous variants in Transient receptor potential melastatin type 7 (TRPM7), encoding an essential and ubiquitously expressed cation channel, may cause hypomagnesemia, but current evidence is insufficient to draw definite conclusions and it is unclear whether any other phenotypes can occur.

Methods

Individuals with unexplained hypomagnesemia underwent whole-exome sequencing which identified TRPM7 variants. Pathogenicity of the identified variants was assessed by combining phenotypic, functional and in silico analyses.

Results

We report three new heterozygous missense variants in TRPM7 (p.Met1000Thr, p.Gly1046Arg, p.Leu1081Arg) in individuals with hypomagnesemia. Strikingly, autism spectrum disorder and developmental delay, mainly affecting speech and motor skills, was observed in all three individuals, while two out of three also presented with seizures. The three variants are predicted to be severely damaging by in silico prediction tools and structural modeling. Furthermore, these variants result in a clear loss-of-function of TRPM7-mediated magnesium uptake in vitro, while not affecting TRPM7 expression or insertion into the plasma membrane.

Conclusions

This study provides additional evidence for the association between heterozygous TRPM7 variants and hypomagnesemia and adds developmental delay to the phenotypic spectrum of TRPM7-related disorders. Considering that the TRPM7 gene is relatively tolerant to loss-of-function variants, future research should aim to unravel by what mechanisms specific heterozygous TRPM7 variants can cause disease.

TECPR2-related hereditary sensory and autonomic neuropathy in two siblings from Palestine

Due to the majority of currently available genome data deriving from individuals of European ancestry, the clinical interpretation of genomic variants in individuals from diverse ethnic backgrounds remains a major diagnostic challenge. Here, we investigated the genetic cause of a complex neurodevelopmental phenotype in two Palestinian siblings. Whole exome sequencing identified a homozygous missense TECPR2 variant (Chr14(GRCh38):g.102425085G>A; NM_014844.5:c.745G>A, p.(Gly249Arg)) absent in gnomAD, segregating appropriately with the inheritance pattern in the family. Variant assessment with in silico pathogenicity prediction and protein modeling tools alongside population database frequencies led to classification as a variant of uncertain significance. As pathogenic TECPR2 variants are associated with hereditary sensory and autonomic neuropathy with intellectual disability, we reviewed previously published candidate TECPR2 missense variants to clarify clinical outcomes and variant classification using current approved guidelines, classifying a number of published variants as of uncertain significance. This work highlights genomic healthcare inequalities and the challenges in interpreting rare genetic variants in populations underrepresented in genomic databases. It also improves understanding of the clinical and genetic spectrum of TECPR2-related neuropathy and contributes to addressing genomic data disparity and inequalities of the genomic architecture in Palestinian populations.

Homozygous missense variants in YKT6 result in loss of function and are associated with developmental delay, with or without severe infantile liver disease and risk for hepatocellular carcinoma

PURPOSE

YKT6 plays important roles in multiple intracellular vesicle trafficking events but has not been associated with Mendelian diseases.

METHODS

We report 3 unrelated individuals with rare homozygous missense variants in YKT6 who exhibited neurological disease with or without a progressive infantile liver disease. We modeled the variants in Drosophila. We generated wild-type and variant genomic rescue constructs of the fly ortholog dYkt6 and compared their ability in rescuing the loss-of-function phenotypes in mutant flies. We also generated a dYkt6KozakGAL4 allele to assess the expression pattern of dYkt6.

RESULTS

Two individuals are homozygous for YKT6 [NM_006555.3:c.554A>G p.(Tyr185Cys)] and exhibited normal prenatal course followed by failure to thrive, developmental delay, and progressive liver disease. Haplotype analysis identified a shared homozygous region flanking the variant, suggesting a common ancestry. The third individual is homozygous for YKT6 [NM_006555.3:c.191A>G p.(Tyr64Cys)] and exhibited neurodevelopmental disorders and optic atrophy. Fly dYkt6 is essential and is expressed in the fat body (analogous to liver) and central nervous system. Wild-type genomic rescue constructs can rescue the lethality and autophagic flux defects, whereas the variants are less efficient in rescuing the phenotypes.

CONCLUSION

The YKT6 variants are partial loss-of-function alleles, and the p.(Tyr185Cys) is more severe than p.(Tyr64Cys).

Clinical Genetic Testing for Hearing Loss: Implications for Genetic Counseling and Gene-Based Therapies

Genetic factors contribute significantly to congenital hearing loss, with non-syndromic cases being more prevalent and genetically heterogeneous. Currently, 150 genes have been associated with non-syndromic hearing loss, and their identification has improved our understanding of auditory physiology and potential therapeutic targets. Hearing loss gene panels offer comprehensive genetic testing for hereditary hearing loss, and advancements in sequencing technology have made genetic testing more accessible and affordable. Currently, genetic panel tests available at a relatively lower cost are offered to patients who face financial barriers. In this study, clinical and audiometric data were collected from six pediatric patients who underwent genetic panel testing. Known pathogenic variants in MYO15A, GJB2, and USH2A were most likely to be causal of hearing loss. Novel pathogenic variants in the MYO7A and TECTA genes were also identified. Variable hearing phenotypes and inheritance patterns were observed amongst individuals with different pathogenic variants. The identification of these variants contributes to the continually expanding knowledge base on genetic hearing loss and lays the groundwork for personalized treatment options in the future.

Measuring the Efficiency of Purging by non-random Mating in Human Populations

Human populations harbor a high concentration of deleterious genetic variants. Here, we tested the hypothesis that non-random mating practices affect the distribution of these variants, through exposure in the homozygous state, leading to their purging from the population gene pool. To do so, we produced whole-genome sequencing data for two pairs of Asian populations exhibiting different alliance rules and rates of inbreeding, but with similar effective population sizes. The results show that populations with higher rates of inbred matings do not purge deleterious variants more efficiently. Purging therefore has a low efficiency in human populations, and different mating practices lead to a similar mutational load.

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.

FLT3L governs the development of partially overlapping hematopoietic lineages in humans and mice

FMS-related tyrosine kinase 3 ligand (FLT3L), encoded by FLT3LG, is a hematopoietic factor essential for the development of natural killer (NK) cells, B cells, and dendritic cells (DCs) in mice. We describe three humans homozygous for a loss-of-function FLT3LG variant with a history of various recurrent infections, including severe cutaneous warts. The patients' bone marrow (BM) was hypoplastic, with low levels of hematopoietic progenitors, particularly myeloid and B cell precursors. Counts of B cells, monocytes, and DCs were low in the patients' blood, whereas the other blood subsets, including NK cells, were affected only moderately, if at all. The patients had normal counts of Langerhans cells (LCs) and dermal macrophages in the skin but lacked dermal DCs. Thus, FLT3L is required for B cell and DC development in mice and humans. However, unlike its murine counterpart, human FLT3L is required for the development of monocytes but not NK cells.

Illuminating the function of the orphan transporter, SLC22A10, in humans and other primates

SLC22A10 is an orphan transporter with unknown substrates and function. The goal of this study is to elucidate its substrate specificity and functional characteristics. In contrast to orthologs from great apes, human SLC22A10, tagged with green fluorescent protein, is not expressed on the plasma membrane. Cells expressing great ape SLC22A10 orthologs exhibit significant accumulation of estradiol-17β-glucuronide, unlike those expressing human SLC22A10. Sequence alignments reveal a proline at position 220 in humans, which is a leucine in great apes. Replacing proline with leucine in SLC22A10-P220L restores plasma membrane localization and uptake function. Neanderthal and Denisovan genomes show proline at position 220, akin to modern humans, indicating functional loss during hominin evolution. Human SLC22A10 is a unitary pseudogene due to a fixed missense mutation, P220, while in great apes, its orthologs transport sex steroid conjugates. Characterizing SLC22A10 across species sheds light on its biological role, influencing organism development and steroid homeostasis.

WilsonGenAI a deep learning approach to classify pathogenic variants in Wilson Disease

BACKGROUND

Advances in Next Generation Sequencing have made rapid variant discovery and detection widely accessible. To facilitate a better understanding of the nature of these variants, American College of Medical Genetics and Genomics and the Association of Molecular Pathologists (ACMG-AMP) have issued a set of guidelines for variant classification. However, given the vast number of variants associated with any disorder, it is impossible to manually apply these guidelines to all known variants. Machine learning methodologies offer a rapid way to classify large numbers of variants, as well as variants of uncertain significance as either pathogenic or benign. Here we classify ATP7B genetic variants by employing ML and AI algorithms trained on our well-annotated WilsonGen dataset.

METHODS

We have trained and validated two algorithms: TabNet and XGBoost on a high-confidence dataset of manually annotated, ACMG & AMP classified variants of the ATP7B gene associated with Wilson's Disease.

RESULTS

Using an independent validation dataset of ACMG & AMP classified variants, as well as a patient set of functionally validated variants, we showed how both algorithms perform and can be used to classify large numbers of variants in clinical as well as research settings.

CONCLUSION

We have created a ready to deploy tool, that can classify variants linked with Wilson's disease as pathogenic or benign, which can be utilized by both clinicians and researchers to better understand the disease through the nature of genetic variants associated with it.

A novel missense variant in the ATPase domain of ATP8A2 and review of phenotypic variability of ATP8A2-related disorders caused by missense changes

ATPase, class 1, type 8A, member 2 (ATP8A2) is a P4-ATPase with a critical role in phospholipid translocation across the plasma membrane. Pathogenic variants in ATP8A2 are known to cause cerebellar ataxia, mental retardation, and disequilibrium syndrome 4 (CAMRQ4) which is often associated with encephalopathy, global developmental delay, and severe motor deficits. Here, we present a family with two siblings presenting with global developmental delay, intellectual disability, spasticity, ataxia, nystagmus, and thin corpus callosum. Whole exome sequencing revealed a homozygous missense variant in the nucleotide binding domain of ATP8A2 (p.Leu538Pro) that results in near complete loss of protein expression. This is in line with other missense variants in the same domain leading to protein misfolding and loss of ATPase function. In addition, by performing diffusion-weighted imaging, we identified bilateral hyperintensities in the posterior limbs of the internal capsule suggesting possible microstructural changes in axon tracts that had not been appreciated before and could contribute to the sensorimotor deficits in these individuals.

Helper T cell immunity in humans with inherited CD4 deficiency

CD4+ T cells are vital for host defense and immune regulation. However, the fundamental role of CD4 itself remains enigmatic. We report seven patients aged 5-61 years from five families of four ancestries with autosomal recessive CD4 deficiency and a range of infections, including recalcitrant warts and Whipple's disease. All patients are homozygous for rare deleterious CD4 variants impacting expression of the canonical CD4 isoform. A shorter expressed isoform that interacts with LCK, but not HLA class II, is affected by only one variant. All patients lack CD4+ T cells and have increased numbers of TCRαβ+CD4-CD8- T cells, which phenotypically and transcriptionally resemble conventional Th cells. Finally, patient CD4-CD8- αβ T cells exhibit intact responses to HLA class II-restricted antigens and promote B cell differentiation in vitro. Thus, compensatory development of Th cells enables patients with inherited CD4 deficiency to acquire effective cellular and humoral immunity against an unexpectedly large range of pathogens. Nevertheless, CD4 is indispensable for protective immunity against at least human papillomaviruses and Trophyrema whipplei.

Prevalence of the hematopoietic rare genetic diseases in Türkiye: A retrospective study

BACKGROUND

Rare genetic diseases are an important global public health problem. At present there are defined approximately 8120 genetic diseases in 15,465 epidemiological datasets and 70% of them start in childhood. Hematopoiesis is the production of all cellular components of blood and continues throughout life.

OBJECTIVE

This study aims to present prevalence of hematopoietic rare genetic diseases recorden in Turkey.

METHODS

The population of study consist of 84.680.273 people who received healthcare from the Turkish National Health Service (49.9% female, 50.1% male). TNHS collects and records electronic data which relates with illness or health information of Turkish population since 2018. All healthcare facilities utilize the Personal Electronic Health Record System (PHR), aligning with standards outlined in the Turkish National Health Data Dictionary and the Health Coding Reference Server (HCRS) established by the Ministry of Health in 2007. The data dictionary comprises essential packages such as patient application and examination records.

RESULTS

Diagnosed female population (53.04%) were higher than male (46.96%). Data shows that most of the people with rare genetic diseases were diagnosed in Marmara Region. The overall prevalence of Hematopoietic Rare Genetic Diseases higher in the years of 2021 and 2022.

CONCLUSION

The prevalence increased gradually from 2018 to 2022. The consanguinity marriage seems to be the main problem which resulted higher rate of rare genetic diseases in Türkiye.

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.

Genetic blueprint of congenital muscular dystrophies with brain malformations in Egypt: A report of 11 families

Congenital muscular dystrophies (CMDs) are a group of rare muscle disorders characterized by early onset hypotonia and motor developmental delay associated with brain malformations with or without eye anomalies in the most severe cases. In this study, we aimed to uncover the genetic basis of severe CMD in Egypt and to determine the efficacy of whole exome sequencing (WES)-based genetic diagnosis in this population. We recruited twelve individuals from eleven families with a clinical diagnosis of CMD with brain malformations that fell into two groups: seven patients with suspected dystroglycanopathy and five patients with suspected merosin-deficient CMD. WES was analyzed by variant filtering using multiple approaches including splicing and copy number variant (CNV) analysis. We identified likely pathogenic variants in FKRP in two cases and variants in POMT1, POMK, and B3GALNT2 in three individuals. All individuals with merosin-deficient CMD had truncating variants in LAMA2. Further analysis in one of the two unsolved cases showed a homozygous protein-truncating variant in Feline Leukemia Virus subgroup C Receptor 1 (FLVCR1). FLVCR1 loss of function has never been previously reported. Yet, loss of function of its paralog, FLVCR2, causes lethal hydranencephaly-hydrocephaly syndrome (Fowler Syndrome) which should be considered in the differential diagnosis for dystroglycanopathy. Overall, we reached a diagnostic rate of 86% (6/7) for dystroglycanopathies and 100% (5/5) for merosinopathy. In conclusion, our results provide further evidence that WES is an important diagnostic method in CMD in developing countries to improve the diagnostic rate, management plan, and genetic counseling for these disorders.

Unbiased phenotype and genotype matching maximizes gene discovery and diagnostic yield

PURPOSE

Widespread application of next-generation sequencing, combined with data exchange platforms, has provided molecular diagnoses for countless families. To maximize diagnostic yield, we implemented an unbiased semi-automated genematching algorithm based on genotype and phenotype matching.

METHODS

Rare homozygous variants identified in 2 or more affected individuals, but not in healthy individuals, were extracted from our local database of ∼12,000 exomes. Phenotype similarity scores (PSS), based on human phenotype ontology terms, were assigned to each pair of individuals matched at the genotype level using HPOsim.

RESULTS

33,792 genotype-matched pairs were discovered, representing variants in 7567 unique genes. There was an enrichment of PSS ≥0.1 among pathogenic/likely pathogenic variant-level pairs (94.3% in pathogenic/likely pathogenic variant-level matches vs 34.75% in all matches). We highlighted founder or region-specific variants as an internal positive control and proceeded to identify candidate disease genes. Variant-level matches were particularly helpful in cases involving inframe indels and splice region variants beyond the canonical splice sites, which may otherwise have been disregarded, allowing for detection of candidate disease genes, such as KAT2A, RPAIN, and LAMP3.

CONCLUSION

Semi-automated genotype matching combined with PSS is a powerful tool to resolve variants of uncertain significance and to identify candidate disease genes.

Hypomagnesaemia with varying degrees of extrarenal symptoms as a consequence of heterozygous CNNM2 variants

Variants in the CNNM2 gene are causative for hypomagnesaemia, seizures and intellectual disability, although the phenotypes can be variable. This study aims to understand the genotype-phenotype relationship in affected individuals with CNNM2 variants by phenotypic, functional and structural analysis of new as well as previously reported variants. This results in the identification of seven variants that significantly affect CNNM2-mediated Mg2+ transport. Pathogenicity of these variants is further supported by structural modelling, which predicts CNNM2 structure to be affected by all of them. Strikingly, seizures and intellectual disability are absent in 4 out of 7 cases, indicating these phenotypes are caused either by specific CNNM2 variant only or by additional risk factors. Moreover, in line with sporadic observations from previous reports, CNNM2 variants might be associated with disturbances in parathyroid hormone and Ca2+ homeostasis.

The genetic landscape of autism spectrum disorder in the Middle Eastern population

Introduction: Autism spectrum disorder (ASD) is characterized by aberrations in social interaction and communication associated with repetitive behaviors and interests, with strong clinical heterogeneity. Genetic factors play an important role in ASD, but about 75% of ASD cases have an undetermined genetic risk. Methods: We extensively investigated an ASD cohort made of 102 families from the Middle Eastern population of Qatar. First, we investigated the copy number variations (CNV) contribution using genome-wide SNP arrays. Next, we employed Next Generation Sequencing (NGS) to identify de novo or inherited variants contributing to the ASD etiology and its associated comorbid conditions in families with complete trios (affected child and the parents). Results: Our analysis revealed 16 CNV regions located in genomic regions implicated in ASD. The analysis of the 88 ASD cases identified 41 genes in 39 ASD subjects with de novo (n = 24) or inherited variants (n = 22). We identified three novel de novo variants in new candidate genes for ASD (DTX4, ARMC6, and B3GNT3). Also, we have identified 15 de novo variants in genes that were previously implicated in ASD or related neurodevelopmental disorders (PHF21A, WASF1, TCF20, DEAF1, MED13, CREBBP, KDM6B, SMURF1, ADNP, CACNA1G, MYT1L, KIF13B, GRIA2, CHM, and KCNK9). Additionally, we defined eight novel recessive variants (RYR2, DNAH3, TSPYL2, UPF3B KDM5C, LYST, and WNK3), four of which were X-linked. Conclusion: Despite the ASD multifactorial etiology that hinders ASD genetic risk discovery, the number of identified novel or known putative ASD genetic variants was appreciable. Nevertheless, this study represents the first comprehensive characterization of ASD genetic risk in Qatar's Middle Eastern population.

Landscape of pharmacogenetic variants associated with non-insulin antidiabetic drugs in the Indian population

INTRODUCTION

Genetic variants contribute to differential responses to non-insulin antidiabetic drugs (NIADs), and consequently to variable plasma glucose control. Optimal control of plasma glucose is paramount to minimizing type 2 diabetes-related long-term complications. India's distinct genetic architecture and its exploding burden of type 2 diabetes warrants a population-specific survey of NIAD-associated pharmacogenetic (PGx) variants. The recent availability of large-scale whole genomes from the Indian population provides a unique opportunity to generate a population-specific map of NIAD-associated PGx variants.

RESEARCH DESIGN AND METHODS

We mined 1029 Indian whole genomes for PGx variants, drug-drug interaction (DDI) and drug-drug-gene interactions (DDGI) associated with 44 NIADs. Population-wise allele frequencies were estimated and compared using Fisher's exact test.

RESULTS

Overall, we found 76 known and 52 predicted deleterious common PGx variants associated with response to type 2 diabetes therapy among Indians. We report remarkable interethnic differences in the relative cumulative counts of decreased and increased response-associated alleles across NIAD classes. Indians and South Asians showed a significant excess of decreased metformin response-associated alleles compared with other global populations. Network analysis of shared PGx genes predicts high DDI risk during coadministration of NIADs with other metabolic disease drugs. We also predict an increased CYP2C19-mediated DDGI risk for CYP3A4/3A5-metabolized NIADs, saxagliptin, linagliptin and glyburide when coadministered with proton-pump inhibitors (PPIs).

CONCLUSIONS

Indians and South Asians have a distinct PGx profile for antidiabetes drugs, marked by an excess of poor treatment response-associated alleles for various NIAD classes. This suggests the possibility of a population-specific reduced drug response in atleast some NIADs. In addition, our findings provide an actionable resource for accelerating future diabetes PGx studies in Indians and South Asians and reconsidering NIAD dosing guidelines to ensure maximum efficacy and safety in the population.

Whole-exome sequencing identifies cancer-associated variants of the endo-lysosomal ion transport channels in the Saudi population

Background

Although national efforts are underway to document the genomic variability of the Saudi population relative to other populations, such variability remains largely unexplored. Genetic variability is known to impact the fate of cells and increase or decrease the risk of a variety of complex diseases including cancer forms. Therefore, the identification of variants associated with cancer susceptibility in Saudi population may protect individuals from cancer or aid in patient-tailored therapies. The endo-lysosomal ion transport genes responsible for cationic ion homeostasis within the cell. We screened 703 single-nucleotide polymorphisms (SNPs) of the endo-lysosomal ion transporter genes in the Saudi population and identified cancer-associated variants that have been reported in other populations.

Methods

Utilizing previously derived local data of Whole-Exome Sequencing (WES), we examined SNPs of TPCN1, TPCN2, P2RX4, TRPM7, TRPV4, TRPV4, and TRPV6 genes. The SNPs were identified for those genes by our in-house database. We predicted the pathogenicity of these variants using in silico tools CADD, Polyphen-2, SIFT, PrimateAI, and FATHMM-XF. Then, we validated our findings by exploring the genetics database (VarSome, dbSNP NCB, OMIM, ClinVar, Ensembl, and GWAS Catalog) to further link cancer risk.

Results

The WES database yielded 703 SNPs found in TPCN2, P2RX4, TRPM7, TRPV4, and TRPV6 genes in 1,144 subjects. The number of variants that were found to be common in our population was 150 SNPs. We identified 13 coding-region non-synonymous variants of the endo-lysosomal genes that were most common with a minor allele frequency (MAF) of ≥ 1 %. Twelve of these variants are rs2376558, rs3750965, rs61746574, rs35264875, rs3829241, rs72928978, rs25644, rs8042919, rs17881456, rs4987682, rs4987667, and rs4987657 that were classified as cancer-associated genes.

Conclusion

Our study highlighted cancer-associated SNPs in the endo-lysosomal genes among Saudi individuals. The allelic frequencies on polymorphic variants confer susceptibility to complex diseases that are comparable to other populations. There is currently insufficient clinical data supporting the link between these SNPs and cancer risk in the Saudi population. Our data argues for initiating future cohort studies in which individuals with the identified SNPs are monitored and assessed for their likelihood of developing malignancies and therapy outcomes.

The immunopathological landscape of human pre-TCRα deficiency: From rare to common variants

We describe humans with rare biallelic loss-of-function PTCRA variants impairing pre-α T cell receptor (pre-TCRα) expression. Low circulating naive αβ T cell counts at birth persisted over time, with normal memory αβ and high γδ T cell counts. Their TCRα repertoire was biased, which suggests that noncanonical thymic differentiation pathways can rescue αβ T cell development. Only a minority of these individuals were sick, with infection, lymphoproliferation, and/or autoimmunity. We also report that 1 in 4000 individuals from the Middle East and South Asia are homozygous for a common hypomorphic PTCRA variant. They had normal circulating naive αβ T cell counts but high γδ T cell counts. Although residual pre-TCRα expression drove the differentiation of more αβ T cells, autoimmune conditions were more frequent in these patients compared with the general population.

Clinical application of next generation sequencing for Mendelian disease diagnosis in the Iranian population

Next-generation sequencing (NGS) has been proven to be one of the most powerful diagnostic tools for rare Mendelian disorders. Several studies on the clinical application of NGS in unselected cohorts of Middle Eastern patients have reported a high diagnostic yield of up to 48%, correlated with a high level of consanguinity in these populations. We evaluated the diagnostic utility of NGS-based testing across different clinical indications in 1436 patients from Iran, representing the first study of its kind in this highly consanguineous population. A total of 1075 exome sequencing and 361 targeted gene panel sequencing were performed over 8 years at a single clinical genetics laboratory, with the majority of cases tested as proband-only (91.6%). The overall diagnostic rate was 46.7%, ranging from 24% in patients with an abnormality of prenatal development to over 67% in patients with an abnormality of the skin. We identified 660 pathogenic or likely pathogenic variants, including 241 novel variants, associated with over 342 known genetic conditions. The highly consanguineous nature of this cohort led to the diagnosis of autosomal recessive disorders in the majority of patients (79.1%) and allowed us to determine the shared carrier status of couples for suspected recessive phenotypes in their deceased child(ren) when direct testing was not possible. We also highlight the observations of recessive inheritance of genes previously associated only with dominant disorders and provide an expanded genotype-phenotype spectrum for multiple less-characterized genes. We present the largest mutational spectrum of known Mendelian disease, including possible founder variants, throughout the Iranian population, which can serve as a unique resource for clinical genomic studies locally and beyond.

Recombinant IFN-γ1b Treatment in a Patient with Inherited IFN-γ Deficiency

PURPOSE

Inborn errors of IFN-γ immunity underlie Mendelian susceptibility to mycobacterial disease (MSMD). Twenty-two genes with products involved in the production of, or response to, IFN-γ and variants of which underlie MSMD have been identified. However, pathogenic variants of IFNG encoding a defective IFN-γ have been described in only two siblings, who both underwent hematopoietic stem cell transplantation (HCST).

METHODS

We characterized a new patient with MSMD by genetic, immunological, and clinical means. Therapeutic decisions were taken on the basis of these findings.

RESULTS

The patient was born to consanguineous Turkish parents and developed bacillus Calmette-Guérin (BCG) disease following vaccination at birth. Whole-exome sequencing revealed a homozygous private IFNG variant (c.224 T > C, p.F75S). Upon overexpression in recipient cells or constitutive expression in the patient's cells, the mutant IFN-γ was produced within the cells but was not correctly folded or secreted. The patient was treated for 6 months with two or three antimycobacterial drugs only and then for 30 months with subcutaneous recombinant IFN-γ1b plus two antimycobacterial drugs. Treatment with IFN-γ1b finally normalized all biological parameters. The patient presented no recurrence of mycobacterial disease or other related infectious diseases. The treatment was well tolerated, without the production of detectable autoantibodies against IFN-γ.

CONCLUSION

We describe a patient with a new form of autosomal recessive IFN-γ deficiency, with intracellular, but not extracellular IFN-γ. IFN-γ1b treatment appears to have been beneficial in this patient, with no recurrence of mycobacterial infection over a period of more than 30 months. This targeted treatment provides an alternative to HCST in patients with complete IFN-γ deficiency or at least an option to better control mycobacterial infection prior to HCST.

The impact of the Turkish population variome on the genomic architecture of rare disease traits

Purpose

The variome of the Turkish (TK) population, a population with a considerable history of admixture and consanguinity, has not been deeply investigated for insights on the genomic architecture of disease.

Methods

We generated and analyzed a database of variants derived from exome sequencing data of 773 TK unrelated, clinically affected individuals with various suspected Mendelian disease traits and 643 unaffected relatives.

Results

Using uniform manifold approximation and projection, we showed that the TK genomes are more similar to those of Europeans and consist of 2 main subpopulations: clusters 1 and 2 (N = 235 and 1181, respectively), which differ in admixture proportion and variome (https://turkishvariomedb.shinyapps.io/tvdb/). Furthermore, the higher inbreeding coefficient values observed in the TK affected compared with unaffected individuals correlated with a larger median span of long-sized (>2.64 Mb) runs of homozygosity (ROH) regions (P value = 2.09e-18). We show that long-sized ROHs are more likely to be formed on recently configured haplotypes enriched for rare homozygous deleterious variants in the TK affected compared with TK unaffected individuals (P value = 3.35e-11). Analysis of genotype-phenotype correlations reveals that genes with rare homozygous deleterious variants in long-sized ROHs provide the most comprehensive set of molecular diagnoses for the observed disease traits with a systematic quantitative analysis of Human Phenotype Ontology terms.

Conclusion

Our findings support the notion that novel rare variants on newly configured haplotypes arising within the recent past generations of a family or clan contribute significantly to recessive disease traits in the TK population.

Genomic disparity impacts variant classification of cancer susceptibility genes in Turkish breast cancer patients

OBJECTIVE

Turkish genome is underrepresented in large genomic databases. This study aims to evaluate the effect of allele frequency in the Turkish population in determining the clinical utility of germline findings in breast cancer, including invasive lobular carcinoma (ILC), mixed invasive ductal and lobular carcinoma (IDC-L), and ductal carcinoma (DC).

METHODS

Two clinic-based cohorts from the Umraniye Research and Training Hospital (URTH) were used in this study: a cohort consisting of 132 women with breast cancer and a non-cancer cohort consisting of 492 participants. The evaluation of the germline landscape was performed by analysis of 27 cancer genes. The frequency and type of variants in the breast cancer cohort were compared to those in the non-cancer cohort to investigate the effect of population genetics. The variant allele frequencies in Turkish Variome and gnomAD were statistically evaluated.

RESULTS

The genetic analysis identified 121 variants in the breast cancer cohort (actionable = 32, VUS = 89) and 223 variants in the non-cancer cohort (actionable = 25, VUS = 188). The occurrence of 21 variants in both suggested a possible genetic population effect. Evaluation of allele frequency of 121 variants from the breast cancer cohort showed 22% had a significantly higher value in Turkish Variome compared to gnomAD (p < 0.0001, 95% CI) with a mean difference of 60 times (ranging from 1.37-354.4). After adjusting for variant allele frequency using the ancestry-appropriate database, 6.7% (5/75) of VUS was reclassified to likely benign.

CONCLUSION

To our knowledge, this is the first study of population genetic effects in breast cancer subtypes in Turkish women. Our findings underscore the need for a large genomic database representing Turkish population-specific variants. It further highlights the significance of the ancestry-appropriate population database for accurate variant assessment in clinical settings.

Should secondary pharmacogenomic variants be actively screened and reported when diagnostic genome-wide sequencing is performed in a child?

This white paper was prepared by the Global Alliance for Genomics and Health Regulatory and Ethics Work Stream's Pediatric Task Team to review and provide perspective with respect to ethical, legal, and social issues regarding the return of secondary pharmacogenomic variants in children who have a serious disease or developmental disorder and are undergoing exome or genome sequencing to identify a genetic cause of their condition. We discuss actively searching for and reporting pharmacogenetic/genomic variants in pediatric patients, different methods of returning secondary pharmacogenomic findings to the patient/parents and/or treating clinicians, maintaining these data in the patient's health record over time, decision supports to assist using pharmacogenetic results in future treatment decisions, and sharing information in public databases to improve the clinical interpretation of pharmacogenetic variants identified in other children. We conclude by presenting a series of points to consider for clinicians and policymakers regarding whether, and under what circumstances, routine screening and return of pharmacogenomic variants unrelated to the indications for testing is appropriate in children who are undergoing genome-wide sequencing to assist in the diagnosis of a suspected genetic disease.

Ethnicity, consanguinity, and genetic architecture of hypertrophic cardiomyopathy

AIMS

Hypertrophic cardiomyopathy (HCM) is characterized by phenotypic heterogeneity that is partly explained by the diversity of genetic variants contributing to disease. Accurate interpretation of these variants constitutes a major challenge for diagnosis and implementing precision medicine, especially in understudied populations. The aim is to define the genetic architecture of HCM in North African cohorts with high consanguinity using ancestry-matched cases and controls.

METHODS AND RESULTS

Prospective Egyptian patients (n = 514) and controls (n = 400) underwent clinical phenotyping and genetic testing. Rare variants in 13 validated HCM genes were classified according to standard clinical guidelines and compared with a prospective HCM cohort of majority European ancestry (n = 684). A higher prevalence of homozygous variants was observed in Egyptian patients (4.1% vs. 0.1%, P = 2 × 10-7), with variants in the minor HCM genes MYL2, MYL3, and CSRP3 more likely to present in homozygosity than the major genes, suggesting these variants are less penetrant in heterozygosity. Biallelic variants in the recessive HCM gene TRIM63 were detected in 2.1% of patients (five-fold greater than European patients), highlighting the importance of recessive inheritance in consanguineous populations. Finally, rare variants in Egyptian HCM patients were less likely to be classified as (likely) pathogenic compared with Europeans (40.8% vs. 61.6%, P = 1.6 × 10-5) due to the underrepresentation of Middle Eastern populations in current reference resources. This proportion increased to 53.3% after incorporating methods that leverage new ancestry-matched controls presented here.

CONCLUSION

Studying consanguineous populations reveals novel insights with relevance to genetic testing and our understanding of the genetic architecture of HCM.

A Novel Truncating Mutation in PAX1 Gene Causes Otofaciocervical Syndrome Without Immunodeficiency

Otofaciocervical syndrome (OTFCS) is a rare genetic disorder of both autosomal recessive and autosomal dominant patterns of inheritance. It is caused by biallelic or monoallelic mutations in PAX1 or EYA1 genes, respectively. Here, we report an OTFCS2 female patient of 1st consanguineous healthy parents. She manifested facial dysmorphism, hearing loss, intellectual disability (ID), and delayed language development (DLD) as the main clinical phenotype. The novel homozygous variant c.1212dup (p.Gly405Argfs*51) in the PAX1 gene was identified by whole exome sequencing (WES), and family segregation confirmed the heterozygous status of the mutation in the parents using the Sanger sequencing. The study recorded a novel PAX1 variant representing the sixth report of OTFCS2 worldwide and the first Egyptian study expanding the geographic area where the disorder was confined.

Somatic Mosaicism in Amyotrophic Lateral Sclerosis and Frontotemporal Dementia Reveals Widespread Degeneration from Focal Mutations

Although mutations in dozens of genes have been implicated in familial forms of amyotrophic lateral sclerosis (fALS) and frontotemporal degeneration (fFTD), most cases of these conditions are sporadic (sALS and sFTD), with no family history, and their etiology remains obscure. We tested the hypothesis that somatic mosaic mutations, present in some but not all cells, might contribute in these cases, by performing ultra-deep, targeted sequencing of 88 genes associated with neurodegenerative diseases in postmortem brain and spinal cord samples from 404 individuals with sALS or sFTD and 144 controls. Known pathogenic germline mutations were found in 20.6% of ALS, and 26.5% of FTD cases. Predicted pathogenic somatic mutations in ALS/FTD genes were observed in 2.7% of sALS and sFTD cases that did not carry known pathogenic or novel germline mutations. Somatic mutations showed low variant allele fraction (typically <2%) and were often restricted to the region of initial discovery, preventing detection through genetic screening in peripheral tissues. Damaging somatic mutations were preferentially enriched in primary motor cortex of sALS and prefrontal cortex of sFTD, mirroring regions most severely affected in each disease. Somatic mutation analysis of bulk RNA-seq data from brain and spinal cord from an additional 143 sALS cases and 23 controls confirmed an overall enrichment of somatic mutations in sALS. Two adult sALS cases were identified bearing pathogenic somatic mutations in DYNC1H1 and LMNA, two genes associated with pediatric motor neuron degeneration. Our study suggests that somatic mutations in fALS/fFTD genes, and in genes associated with more severe diseases in the germline state, contribute to sALS and sFTD, and that mosaic mutations in a small fraction of cells in focal regions of the nervous system can ultimately result in widespread degeneration.

Identification and characterization of ATM founder mutation in BRCA-negative breast cancer patients of Arab ethnicity

Breast cancer (BC) is the most prevalent malignancy among women worldwide with germline pathogenic variants/likely pathogenic variants (PVs/LPVs) in BRCA1/2 accounting for a large portion of hereditary cases. Recently, heterozygous PVs/LPVs in the ATM serine/threonine kinase or Ataxia-telangiectasia mutated gene (ATM) has been identified as a moderate susceptibility factor for BC in diverse ethnicities. However, the prevalence of ATM PVs/LPVs in BC susceptibility in Arab populations remains largely unexplored. This study investigated the prevalence of ATM PVs/LPVs among BC patients from Saudi Arabia, employing capture-sequencing technology for ATM PVs/LPVs screening in a cohort of 715 unselected BC patients without BRCA1/2 PVs/LPVs. In addition, founder mutation analysis was conducted using the PHASE program. In our entire cohort, four unique PVs/LPVs in the ATM gene were identified in six cases (0.8%). Notably, one recurrent LPV, c.6115G > A:p.Glu2039Lys was identified in three cases, for which haplotype analysis confirmed as a novel putative founder mutation traced back to 13 generations on average. This founder mutation accounted for half of all identified mutant cases and 0.4% of total screened cases. This study further reveals a significant correlation between the presence of ATM mutation and family history of BC (p = 0.0127). These findings underscore an approximate 0.8% prevalence of ATM germline PVs/LPVs in Arab BC patients without BRCA1/2 PVs/LPVs and suggest a founder effect of specific recurrent ATM mutation. These insights can help in the design of a genetic testing strategy tailored to the local population in Saudi Arabia, thereby, enabling more accurate clinical management and risk prediction.

Genetic determinants of severe COVID-19 in young Asian and Middle Eastern patients: a case series

Studies of genetic factors associated with severe COVID-19 in young adults have been limited in non-Caucasian populations. Here, we clinically characterize a case series of patients with COVID-19, who were otherwise healthy, young adults (N = 55; mean age 34.1 ± SD 5.0 years) from 16 Asian, Middle Eastern, and North African countries. Using whole exome sequencing, we identify rare, likely deleterious variants affecting 16 immune-related genes in 17 out of 55 patients (31%), including 7 patients (41% of all carriers or 12.7% of all patients) who harbored multiple such variants mainly in interferon and toll-like receptor genes. Protein network analysis as well as transcriptomic analysis of nasopharyngeal swabs from an independent COVID-19 cohort (N = 50; 42% Asians and 22% Arabs) revealed that most of the altered genes, as identified by whole exome sequencing, and the associated molecular pathways were significantly altered in COVID-19 patients. Genetic variants tended to be associated with mortality, intensive care admission, and ventilation support. Our clinical cases series, genomic and transcriptomic findings suggest a possible role for interferon pathway genes in severe COVID-19 and highlight the importance of extending genetic studies to diverse populations to better understand the human genetics of disease.

Variants in EFCAB7 underlie nonsyndromic postaxial polydactyly

Polydactyly is the most common limb malformation that occurs in 1.6-10.6 per one thousand live births, with incidence varying with ancestry. The underlying gene has been identified for many of the ~100 syndromes that include polydactyly. While for the more common form, nonsydromic polydactyly, eleven candidate genes have been reported. We investigated the underlying genetic cause of autosomal recessive nonsyndromic postaxial polydactyly in four consanguineous Pakistani families. Some family members with postaxial polydactyly also present with syndactyly, camptodactyly, or clinodactyly. Analysis of the exome sequence data revealed two novel homozygous frameshift deletions in EFCAB7: [c.830delG;p.(Gly277Valfs*5)]; in three families and [c.1350_1351delGA;p.(Asn451Phefs*2)] in one family. Sanger sequencing confirmed that these variants segregated with postaxial polydactyly, i.e., family members with postaxial polydactyly were found to be homozygous while unaffected members were heterozygous or wild type. EFCAB7 displays expressions in the skeletal muscle and on the cellular level in cilia. IQCE-EFCAB7 and EVC-EVC2 are part of the heterotetramer EvC complex, which is a positive regulator of the Hedgehog (Hh) pathway, that plays a key role in limb formation. Depletion of either EFCAB7 or IQCE inhibits induction of Gli1, a direct Hh target gene. Variants in IQCE and GLI1 have been shown to cause nonsyndromic postaxial polydactyly, while variants in EVC and EVC2 underlie Ellis van Creveld and Weyers syndromes, which include postaxial polydactyly as a phenotype. This is the first report of the involvement of EFCAB7 in human disease etiology.

Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates

SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17β-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.

Interplay of Mendelian and polygenic risk factors in Arab breast cancer patients

BACKGROUND

Breast cancer patients from the indigenous Arab population present much earlier than patients from Western countries and have traditionally been underrepresented in cancer genomics studies. The contribution of polygenic and Mendelian risk toward the earlier onset of breast cancer in the population remains elusive.

METHODS

We performed low-pass whole genome sequencing (lpWGS) and whole-exome sequencing (WES) from 220 female breast cancer patients unselected for positive family history from the indigenous Arab population. Using publicly available resources, we imputed population-specific variants and calculated breast cancer burden-sensitive polygenic risk scores (PRS). Variant pathogenicity was also evaluated on exome variants with high coverage.

RESULTS

Variants imputed from lpWGS showed high concordance with paired exome (median dosage correlation: 0.9459, Interquartile range: 0.9410-0.9490). After adjusting the PRS to the Arab population, we found significant associations between PRS performance in risk prediction and first-degree relative breast cancer history prediction (Spearman rho=0.43, p = 0.03), where breast cancer patients in the top PRS decile are 5.53 (95% CI 1.76-17.97, p = 0.003) times more likely also to have a first-degree relative diagnosed with breast cancer compared to those in the middle deciles. In addition, we found evidence for the genetic liability threshold model of breast cancer where among patients with a family history of breast cancer, pathogenic rare variant carriers had significantly lower PRS than non-carriers (p = 0.0205, Mann-Whitney U test) while for non-carriers every standard deviation increase in PRS corresponded to 4.52 years (95% CI 8.88-0.17, p = 0.042) earlier age of presentation.

CONCLUSIONS

Overall, our study provides a framework to assess polygenic risk in an understudied population using lpWGS and identifies common variant risk as a factor independent of pathogenic variant carrier status for earlier age of onset of breast cancer among indigenous Arab breast cancer patients.