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.

The genomic landscape of CYP2D6 variation in the Indian population.

Aim: The CYP2D6 gene is highly polymorphic, causing large interindividual variability in the metabolism of several clinically important drugs. Materials & methods: The authors investigated the diversity and distribution of CYP2D6 alleles in Indians using whole genome sequences (N = 1518). Functional consequences were assessed using pathogenicity scores and molecular dynamics simulations. Results: The analysis revealed population-specific CYP2D6 alleles (*86, *7, *111, *112, *113, *99) and remarkable differences in variant and phenotype frequencies with global populations. The authors observed that one in three Indians could benefit from a dose alteration for psychiatric drugs with accurate CYP2D6 phenotyping. Molecular dynamics simulations revealed large conformational fluctuations, confirming the predicted reduced function of *86 and *113 alleles. Conclusion: The findings emphasize the utility of comprehensive CYP2D6 profiling for aiding precision public health.

Development of pathophysiologically relevant models of sickle cell disease and ß-thalassemia for therapeutic studies.

Ex vivo cellular system that accurately replicates sickle cell disease and ß-thalassemia characteristics is a highly sought-after goal in the field of erythroid biology. In this study, we present the generation of erythroid progenitor lines with sickle cell disease and ß-thalassemia mutation using CRISPR/Cas9. The disease cellular models exhibit similar differentiation profiles, globin expression and proteome dynamics as patient-derived hematopoietic stem/progenitor cells. Additionally, these cellular models recapitulate pathological conditions associated with both the diseases. Hydroxyurea and pomalidomide treatment enhanced fetal hemoglobin levels. Notably, we introduce a therapeutic strategy for the above diseases by recapitulating the HPFH3 genotype, which reactivates fetal hemoglobin levels and rescues the disease phenotypes, thus making these lines a valuable platform for studying and developing new therapeutic strategies. Altogether, we demonstrate our disease cellular systems are physiologically relevant and could prove to be indispensable tools for disease modeling, drug screenings and cell and gene therapy-based applications.

Whole genome sequencing followed by functional analysis of genomic deletion encompassing ERCC8 and NDUFAF2 genes in a non-consanguineous Indian family reveals dysfunctional mitochondrial bioenergetics leading to infant mortality.

Genomic investigations on an infant who presented with a putative mitochondrial disorder led to identification of compound heterozygous deletion with an overlapping region of ∼142 kb encompassing two nuclear encoded genes namely ERCC8 and NDUFAF2. Investigations on fetal-derived fibroblast culture demonstrated impaired bioenergetics and mitochondrial dysfunction, which explains the phenotype and observed infant mortality in the present study. The genetic findings from this study extended the utility of whole-genome sequencing as it led to development of a MLPA-based assay for carrier screening in the extended family and the prenatal testing aiding in the birth of two healthy children.

Scalable noninvasive amplicon-based precision sequencing (SNAPseq) for genetic diagnosis and screening of ß-thalassemia and sickle cell disease using a next-generation sequencing platform.

ß-hemoglobinopathies such as ß-thalassemia (BT) and Sickle cell disease (SCD) are inherited monogenic blood disorders with significant global burden. Hence, early and affordable diagnosis can alleviate morbidity and reduce mortality given the lack of effective cure. Currently, Sanger sequencing is considered to be the gold standard genetic test for BT and SCD, but it has a very low throughput requiring multiple amplicons and more sequencing reactions to cover the entire HBB gene. To address this, we have demonstrated an extraction-free single amplicon-based approach for screening the entire ß-globin gene with clinical samples using Scalable noninvasive amplicon-based precision sequencing (SNAPseq) assay catalyzing with next-generation sequencing (NGS). We optimized the assay using noninvasive buccal swab samples and simple finger prick blood for direct amplification with crude lysates. SNAPseq demonstrates high sensitivity and specificity, having a 100% agreement with Sanger sequencing. Furthermore, to facilitate seamless reporting, we have created a much simpler automated pipeline with comprehensive resources for pathogenic mutations in BT and SCD through data integration after systematic classification of variants according to ACMG and AMP guidelines. To the best of our knowledge, this is the first report of the NGS-based high throughput SNAPseq approach for the detection of both BT and SCD in a single assay with high sensitivity in an automated pipeline.

Genetic landscape of human platelet antigen variants in the Indian population analysed from 1029 whole genomes.

Genetic variants in human platelet antigens (HPAs) considered allo- or auto antigens are associated with various disorders, including neonatal alloimmune thrombocytopenia, platelet transfusion refractoriness and post-transfusion purpura. Although global differences in genotype frequencies were observed, the distributions of HPA variants in the Indian population are largely unknown. This study aims to explore the landscape of HPA variants in India to provide a basis for risk assessment and management of related complications. Population-specific frequencies of genetic variants associated with the 35 classes of HPAs (HPA-1 to HPA-35) were estimated by systematically analysing genomic variations of 1029 healthy Indian individuals as well as from global population genome datasets. Allele frequencies of the most clinically relevant HPA systems in the Indian population were found as follows, HPA-1a - 0.884, HPA-1b - 0.117, HPA-2a - 0.941, HPA-2b - 0.059, HPA-3a - 0.653, HPA-3b - 0.347, HPA-4a - 0.999, HPA-4b - 0.0010, HPA-5a - 0.923, HPA-5b - 0.077, HPA-6a - 0.998, HPA-6b - 0.002, HPA-15a - 0.582 and HPA-15b - 0.418. This study provides the first comprehensive analysis of HPA allele and genotype frequencies using large scale representative whole genome sequencing data of the Indian population.

The genome sequence of lumpy skin disease virus from an outbreak in India suggests a distinct lineage of the virus.

Although previously confined to regions within Africa, lumpy skin disease virus (LSDV) infections have caused significantly large outbreaks in several regions of the world in recent years. In 2019, an outbreak of the disease was reported in India with low rates of morbidity and no reported mortality. However, in 2022, an ongoing outbreak of LSDV spanning over seven states in India resulted in the loss of over 80,000 cattle over a period of three months. Here, we report complete genome sequences of six isolates of LSDV collected from affected cattle during an ongoing outbreak of the disease in Rajasthan, India. Analysis of these sequences showed that the genome isolates from the 2022 outbreak have a large number of genetic variations compared to the reference strain and that they form a distinct genetic lineage. This report thus highlights the importance of genome sequencing and surveillance of transboundary infectious agents to track the prevalence and emergence of variants.

Whole-genome sequencing of 1029 Indian individuals reveals unique and rare structural variants.

Structural variants contribute to genetic variability in human genomes and they can be presented in population-specific patterns. We aimed to understand the landscape of structural variants in the genomes of healthy Indian individuals and explore their potential implications in genetic disease conditions. For the identification of structural variants, a whole genome sequencing dataset of 1029 self-declared healthy Indian individuals from the IndiGen project was analysed. Further, these variants were evaluated for potential pathogenicity and their associations with genetic diseases. We also compared our identified variations with the existing global datasets. We generated a compendium of total 38,560 high-confident structural variants, comprising 28,393 deletions, 5030 duplications, 5038 insertions, and 99 inversions. Particularly, we identified around 55% of all these variants were found to be unique to the studied population. Further analysis revealed 134 deletions with predicted pathogenic/likely pathogenic effects and their affected genes were majorly enriched for neurological disease conditions, such as intellectual disability and neurodegenerative diseases. The IndiGenomes dataset helped us to understand the unique spectrum of structural variants in the Indian population. More than half of identified variants were not present in the publicly available global dataset on structural variants. Clinically important deletions identified in IndiGenomes might aid in improving the diagnosis of unsolved genetic diseases, particularly in neurological conditions. Along with basal allele frequency data and clinically important deletions, IndiGenomes data might serve as a baseline resource for future studies on genomic structural variant analysis in the Indian population.

Genetic landscape of human neutrophil antigen variants in India from population-scale genomes.

Antibodies against human neutrophil antigens (HNAs) play a significant role in various clinical conditions such as neonatal alloimmune neutropenia, transfusion-related acute lung injury, and other nonhemolytic transfusion reactions. This study aims to identify the genotype and allele frequencies of HNAs in the healthy Indian population. Ten genetic variants in four human genes encoding alleles of HNAs class I-V approved by the International Society of Blood Transfusion-Granulocyte Immunobiology Working Party were used in the analysis. Genetic variants from whole genome sequences of 1029 healthy Indian individuals corresponding to HNA alleles were analyzed. The frequencies of the variants were compared with global population datasets using an in-house computational pipeline. In HNA class I, allele frequencies of FCGR3B*01, FCGR3B*02, and FCGR3B*03 encoding HNA-1a, HNA-1b, and HNA-1c were 0.07, 0.8, and 0.13, respectively. HNA class 3 alleles namely SLC44A2*01 (encoding HNA-3a) and SLC44A2*02 (encoding HNA-3b) were found at allele frequencies of 0.78 and 0.22, respectively. The frequencies of ITGAM*01 encoding HNA-4a and ITGAM*02 encoding HNA-4a were 0.95 and 0.05, respectively. Furthermore, allele frequencies of HNA class 5 alleles were 0.32 for ITGAL*01 (encoding HNA-5a) and 0.68 for ITGAL*02 (encoding HNA-5b). Interestingly, it was also found that rs2230433 variant deciding the HNA class 5 alleles, was highly prevalent (78.2%) in the Indian population compared with other global populations. This study presents the first comprehensive report of HNA allele and genotype frequencies in the Indian population using population genome datasets of 1029 individuals. Significant difference was observed in the prevalence of HNA5a and HNA5b in India in comparison to other global populations.

Landscape of Variability in Chemosensory Genes Associated With Dietary Preferences in Indian Population: Analysis of 1029 Indian Genomes.

Perception and preferences for food and beverages determine dietary behaviour and health outcomes. Inherent differences in chemosensory genes, ethnicity, geo-climatic conditions, and sociocultural practices are other determinants. We aimed to study the variation landscape of chemosensory genes involved in perception of taste, texture, odour, temperature and burning sensations through analysis of 1,029 genomes of the IndiGen project and diverse continental populations. SNPs from 80 chemosensory genes were studied in whole genomes of 1,029 IndiGen samples and 2054 from the 1000 Genomes project. Population genetics approaches were used to infer ancestry of IndiGen individuals, gene divergence and extent of differentiation among studied populations. 137,760 SNPs including common and rare variants were identified in IndiGenomes with 62,950 novel (46%) and 48% shared with the 1,000 Genomes. Genes associated with olfaction harbored most SNPs followed by those associated with differences in perception of salt and pungent tastes. Across species, receptors for bitter taste were the most diverse compared to others. Three predominant ancestry groups within IndiGen were identified based on population structure analysis. We also identified 1,184 variants that exhibit differences in frequency of derived alleles and high population differentiation (FST ≥0.3) in Indian populations compared to European, East Asian and African populations. Examples include ADCY10, TRPV1, RGS6, OR7D4, ITPR3, OPRM1, TCF7L2, and RUNX1. This is a first of its kind of study on baseline variations in genes that could govern cuisine designs, dietary preferences and health outcomes. This would be of enormous utility in dietary recommendations for precision nutrition both at population and individual level.

Comprehensive Assessment of Indian Variations in the Druggable Kinome Landscape Highlights Distinct Insights at the Sequence, Structure and Pharmacogenomic Stratum.

India confines more than 17% of the world's population and has a diverse genetic makeup with several clinically relevant rare mutations belonging to many sub-group which are undervalued in global sequencing datasets like the 1000 Genome data (1KG) containing limited samples for Indian ethnicity. Such databases are critical for the pharmaceutical and drug development industry where diversity plays a crucial role in identifying genetic disposition towards adverse drug reactions. A qualitative and comparative sequence and structural study utilizing variant information present in the recently published, largest curated Indian genome database (IndiGen) and the 1000 Genome data was performed for variants belonging to the kinase coding genes, the second most targeted group of drug targets. The sequence-level analysis identified similarities and differences among different populations based on the nsSNVs and amino acid exchange frequencies whereas a comparative structural analysis of IndiGen variants was performed with pathogenic variants reported in UniProtKB Humsavar data. The influence of these variations on structural features of the protein, such as structural stability, solvent accessibility, hydrophobicity, and the hydrogen-bond network was investigated. In-silico screening of the known drugs to these Indian variation-containing proteins reveals critical differences imparted in the strength of binding due to the variations present in the Indian population. In conclusion, this study constitutes a comprehensive investigation into the understanding of common variations present in the second largest population in the world and investigating its implications in the sequence, structural and pharmacogenomic landscape. The preliminary investigation reported in this paper, supporting the screening and detection of ADRs specific to the Indian population could aid in the development of techniques for pre-clinical and post-market screening of drug-related adverse events in the Indian population.

Pharmacogenomic landscape of Indian population using whole genomes.

Ethnic differences in pharmacogenomic (PGx) variants have been well documented in literature and could significantly impact variability in response and adverse events to therapeutics. India is a large country with diverse ethnic populations of distinct genetic architecture. India's national genome sequencing initiative (IndiGen) provides a unique opportunity to explore the landscape of PGx variants using population-scale whole genome sequences. We have analyzed the IndiGen variation dataset (N = 1029 genomes) along with global population scale databases to map the most prevalent clinically actionable and potentially deleterious PGx variants among Indians. Differential frequencies for the known and novel variants were studied and interaction of the disrupted PGx genes affecting drug responses were analyzed by performing a pathway analysis. We have highlighted significant differences in the allele frequencies of clinically actionable PGx variants in Indians when compared to the global populations. We identified 134 mostly common (allele frequency [AF] > 0.1) potentially deleterious PGx variants that could alter or inhibit the function of 102 pharmacogenes in Indians. We also estimate that on, an average, each Indian individual carried eight PGx variants (single nucleotide variants) that have a direct impact on the choice of treatment or drug dosing. We have also highlighted clinically actionable PGx variants and genes for which preemptive genotyping is most recommended for the Indian population. The study has put forward the most comprehensive PGx landscape of the Indian population from whole genomes that could enable optimized drug selection and genotype-guided prescriptions for improved therapeutic outcomes and minimizing adverse events.

An Alu insertion map of the Indian population: identification and analysis in 1021 genomes of the IndiGen project.

Actively retrotransposing primate-specific Alu repeats display insertion-deletion (InDel) polymorphism through their insertion at new loci. In the global datasets, Indian populations remain under-represented and so do their Alu InDels. Here, we report the genomic landscape of Alu InDels from the recently released 1021 Indian Genomes (IndiGen) (available at https://clingen.igib.res.in/indigen). We identified 9239 polymorphic Alu insertions that include private (3831), rare (3974) and common (1434) insertions with an average of 770 insertions per individual. We achieved an 89% PCR validation of the predicted genotypes in 94 samples tested. About 60% of identified InDels are unique to IndiGen when compared to other global datasets; 23% of sites were shared with both SGDP and HGSVC; among these, 58% (1289 sites) were common polymorphisms in IndiGen. The insertions not only show a bias for genic regions, with a preference for introns but also for the associated genes showing enrichment for processes like cell morphogenesis and neurogenesis (P-value < 0.05). Approximately, 60% of InDels mapped to genes present in the OMIM database. Finally, we show that 558 InDels can serve as ancestry informative markers to segregate global populations. This study provides a valuable resource for baseline Alu InDels that would be useful in population genomics.

Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes.

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

Clinical, immunological and genomic characteristics of children with X-linked agammaglobulinemia from Kerala, South India.

X-linked agammaglobulinemia (XLA) is an X-linked recessive primary immunodeficiency disorder caused due to a pathogenic variant in the Bruton tyrosine (BTK) gene with an incidence of 1:379,000 live births and 1:190,000 male births. Patients affected with XLA present with recurrent infections of the gastrointestinal and respiratory tracts. Here we report the first case series of 17 XLA patients of 10 South Indian families with a wide spectrum of clinical and genetic features. In our cohort, patients presented mainly with recurrent pneumonia, gastrointestinal infection, otitis media, pyoderma, abscesses, empyema, arthritis, and osteomyelitis. Using next-generation and Sanger sequencing we have identified 10 unique pathogenic and likely pathogenic variants in 17 patients. This encompasses three nonsynonymous, two stop-gain, two frameshifts, two structural, and one splicing variant, out of which two of them are novel. Based on the type of variant, patients had variable clinical features and treatment responses. We have also evaluated Btk protein expression for six patients in comparison to the healthy individuals and determined mosaic Btk expression patterns in four mothers. We have also performed family screening in 6 families using Sanger sequencing and identified 19 carriers for the variant. The diagnosis for the patients led to the proper treatment i.e. 15 patients were on intravenous immunoglobulin (IVIG) and the other two had successful hematopoietic stem cell transplantation (HSCT). Unfortunately, two of our patients died due to sepsis, while on IVIG. We envision the present study could help in better understanding of patients with XLA and help in family screening and prenatal diagnosis. To the best of our knowledge, this is the largest case series of patients affected with XLA from South India.

Variants of concern responsible for SARS-CoV-2 vaccine breakthrough infections from India.

Emerging reports of SARS-CoV-2 breakthrough infections entail methodical genomic surveillance for determining the efficacy of vaccines. This study elaborates genomic analysis of isolates from breakthrough infections following vaccination with AZD1222/Covishield and BBV152/Covaxin. Variants of concern B.1.617.2 and B.1.1.7 responsible for cases surge in April-May 2021 in Delhi, were the predominant lineages among breakthrough infections.

Genetic epidemiology of autoinflammatory disease variants in Indian population from 1029 whole genomes.

BACKGROUND: Autoinflammatory disorders are the group of inherited inflammatory disorders caused due to the genetic defect in the genes that regulates innate immune systems. These have been clinically characterized based on the duration and occurrence of unprovoked fever, skin rash, and patient's ancestry. There are several autoinflammatory disorders that are found to be prevalent in a specific population and whose disease genetic epidemiology within the population has been well understood. However, India has a limited number of genetic studies reported for autoinflammatory disorders till date. The whole genome sequencing and analysis of 1029 Indian individuals performed under the IndiGen project persuaded us to perform the genetic epidemiology of the autoinflammatory disorders in India. RESULTS: We have systematically annotated the genetic variants of 56 genes implicated in autoinflammatory disorder. These genetic variants were reclassified into five categories (i.e., pathogenic, likely pathogenic, benign, likely benign, and variant of uncertain significance (VUS)) according to the American College of Medical Genetics and Association of Molecular pathology (ACMG-AMP) guidelines. Our analysis revealed 20 pathogenic and likely pathogenic variants with significant differences in the allele frequency compared with the global population. We also found six causal founder variants in the IndiGen dataset belonging to different ancestry. We have performed haplotype prediction analysis for founder mutations haplotype that reveals the admixture of the South Asian population with other populations. The cumulative carrier frequency of the autoinflammatory disorder in India was found to be 3.5% which is much higher than reported. CONCLUSION: With such frequency in the Indian population, there is a great need for awareness among clinicians as well as the general public regarding the autoinflammatory disorder. To the best of our knowledge, this is the first and most comprehensive population scale genetic epidemiological study being reported from India.

Genomic characterization and epidemiology of an emerging SARS-CoV-2 variant in Delhi, India.

Delhi, the national capital of India, experienced multiple severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreaks in 2020 and reached population seropositivity of >50% by 2021. During April 2021, the city became overwhelmed by COVID-19 cases and fatalities, as a new variant, B.1.617.2 (Delta), replaced B.1.1.7 (Alpha). A Bayesian model explains the growth advantage of Delta through a combination of increased transmissibility and reduced sensitivity to immune responses generated against earlier variants (median estimates: 1.5-fold greater transmissibility and 20% reduction in sensitivity). Seropositivity of an employee and family cohort increased from 42% to 87.5% between March and July 2021, with 27% reinfections, as judged by increased antibody concentration after a previous decline. The likely high transmissibility and partial evasion of immunity by the Delta variant contributed to an overwhelming surge in Delhi.

An optimized, amplicon-based approach for sequencing of SARS-CoV-2 from patient samples using COVIDSeq assay on Illumina MiSeq sequencing platforms.

Sequencing of SARS-CoV-2 genomes is crucial for understanding the genetic epidemiology of the COVID-19 pandemic. It is also critical for understanding the evolution of the virus and also for the rapid development of diagnostic tools. The present protocol is a modification of the Illumina COVIDSeq test. We describe an amplicon-based next-generation sequencing approach with short turnaround time, adapted for bench-top sequencers like MiSeq, iSeq, and MiniSeq. For complete details on the use and execution of this protocol, please refer to Bhoyar et al. (2021).

Whole genome sequencing identifies novel structural variant in a large Indian family affected with X-linked agammaglobulinemia.

X-linked agammaglobulinemia (XLA, OMIM #300755) is a primary immunodeficiency disorder caused by pathogenic variations in the BTK gene, characterized by failure of development and maturation of B lymphocytes. The estimated prevalence worldwide is 1 in 190,000 male births. Recently, genome sequencing has been widely used in difficult to diagnose and familial cases. We report a large Indian family suffering from XLA with five affected individuals. We performed complete blood count, immunoglobulin assay, and lymphocyte subset analysis for all patients and analyzed Btk expression for one patient and his mother. Whole exome sequencing (WES) for four patients, and whole genome sequencing (WGS) for two patients have been performed. Carrier screening was done for 17 family members using Multiplex Ligation-dependent Probe Amplification (MLPA) and haplotype ancestry mapping using fineSTRUCTURE was performed. All patients had hypogammaglobulinemia and low CD19+ B cells. One patient who underwent Btk estimation had low expression and his mother showed a mosaic pattern. We could not identify any single nucleotide variants or small insertion/ deletions from the WES dataset that correlates with the clinical feature of the patient. Structural variant analysis through WGS data identifies a novel large deletion of 5,296 bp at loci chrX:100,624,323-100,629,619 encompassing exons 3-5 of the BTK gene. Family screening revealed seven carriers for the deletion. Two patients had a successful HSCT. Haplotype mapping revealed a South Asian ancestry. WGS led to identification of the accurate genetic mutation which could help in early diagnosis leading to improved outcomes, prevention of permanent organ damage and improved quality of life, as well as enabling genetic counselling and prenatal diagnosis in the family.