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.

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.

Whole Genome DNA Methylation and Gene Expression Profiling of Oropharyngeal Cancer Patients in North-Eastern India: Identification of Epigenetically Altered Gene Expression Reveals Potential Biomarkers.

Oropharyngeal cancer is a subtype of head and neck squamous cell carcinoma that is associated with unique risk exposures like consumption of smokeless tobacco and areca nut and is highly prevalent in the northeastern region of India, especially Meghalaya. However, the underlying epigenetic and transcriptomic changes in this cancer type is yet to be delineated. We have undertaken a study on genome wide somatic alterations in the DNA methylation and transcriptome in oropharyngeal cancer patients from this region using genome wide techniques in paired tumors and adjacent normal tissues. By using integrative approaches, we have identified 194 epigenetically silenced and 241 epigenetically overexpressed genes in the tumor tissue of these patients. Pathways that are significantly enriched by these genes include the pathways of immune systems, such as the interleukin signaling pathways and Toll-like receptor signaling pathway. Also, osteoclast differentiation pathway was found to be epigenetically upregulated. The pathways enriched by the epigenetically downregulated genes were found to be predominantly those involved in xenobiotic metabolism and keratinization. Two major transcription factors - SPI1 and RUNX1 were identified as epigenetically dysregulated, which further modulates 129 downstream genes. Comparison of our observations with the head and neck cancer data from TCGA revealed distinct DNA methylation and gene expression landscapes which might be specific for oropharyngeal cancer. HPV DNA sequences were not detected in any of the tumor samples in RNA-Seq data. The results obtained in this study might provide improved understanding of the disease.