LncRNA VEAL2 regulates PRKCB2 to modulate endothelial permeability in diabetic retinopathy.

Long non-coding RNAs (lncRNAs) are emerging as key regulators of endothelial cell function. Here, we investigated the role of a novel vascular endothelial-associated lncRNA (VEAL2) in regulating endothelial permeability. Precise editing of veal2 loci in zebrafish (veal2gib005Δ8/+ ) induced cranial hemorrhage. In vitro and in vivo studies revealed that veal2 competes with diacylglycerol for interaction with protein kinase C beta-b (Prkcbb) and regulates its kinase activity. Using PRKCB2 as bait, we identified functional ortholog of veal2 in humans from HUVECs and named it as VEAL2. Overexpression and knockdown of VEAL2 affected tubulogenesis and permeability in HUVECs. VEAL2 was differentially expressed in choroid tissue in eye and blood from patients with diabetic retinopathy, a disease where PRKCB2 is known to be hyperactivated. Further, VEAL2 could rescue the effects of PRKCB2-mediated turnover of endothelial junctional proteins thus reducing hyperpermeability in hyperglycemic HUVEC model of diabetic retinopathy. Based on evidence from zebrafish and hyperglycemic HUVEC models and diabetic retinopathy patients, we report a hitherto unknown VEAL2 lncRNA-mediated regulation of PRKCB2, for modulating junctional dynamics and maintenance of endothelial permeability.

Whole genome sequencing of families diagnosed with cardiac channelopathies reveals structural variants missed by whole exome sequencing.

Cardiac channelopathies are a group of heritable disorders that affect the heart's electrical activity due to genetic variations present in genes coding for ion channels. With the advent of new sequencing technologies, molecular diagnosis of these disorders in patients has paved the way for early identification, therapeutic management and family screening. The objective of this retrospective study was to understand the efficacy of whole-genome sequencing in diagnosing patients with suspected cardiac channelopathies who were reported negative after whole exome sequencing and analysis. We employed a 3-tier analysis approach to identify nonsynonymous variations and loss-of-function variations missed by exome sequencing, and structural variations that are better resolved only by sequencing whole genomes. By performing whole genome sequencing and analyzing 25 exome-negative cardiac channelopathy patients, we identified 3 pathogenic variations. These include a heterozygous likely pathogenic nonsynonymous variation, CACNA1C:NM_000719:exon19:c.C2570G:p. P857R, which causes autosomal dominant long QT syndrome in the absence of Timothy syndrome, a heterozygous loss-of-function variation CASQ2:NM_001232.4:c.420+2T>C classified as pathogenic, and a 9.2 kb structural variation that spans exon 2 of the KCNQ1 gene, which is likely to cause Jervell-Lange-Nielssen syndrome. In addition, we also identified a loss-of-function variation and 16 structural variations of unknown significance (VUS). Further studies are required to elucidate the role of these identified VUS in gene regulation and decipher the underlying genetic and molecular mechanisms of these disorders. Our present study serves as a pilot for understanding the utility of WGS over clinical exomes in diagnosing cardiac channelopathy disorders.

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.

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.

1029 genomes of self-declared healthy individuals from India reveal prevalent and clinically relevant cardiac ion channelopathy variants.

BACKGROUND: The prevalence and genetic spectrum of cardiac channelopathies exhibit population-specific differences. We aimed to understand the spectrum of cardiac channelopathy-associated variations in India, which is characterised by a genetically diverse population and is largely understudied in the context of these disorders. RESULTS: We utilised the IndiGenomes dataset comprising 1029 whole genomes from self-declared healthy individuals as a template to filter variants in 36 genes known to cause cardiac channelopathies. Our analysis revealed 186,782 variants, of which we filtered 470 variants that were identified as possibly pathogenic (440 nonsynonymous, 30 high-confidence predicted loss of function ). About 26% (124 out of 470) of these variants were unique to the Indian population as they were not reported in the global population datasets and published literature. Classification of 470 variants by ACMG/AMP guidelines unveiled 13 pathogenic/likely pathogenic (P/LP) variants mapping to 19 out of the 1029 individuals. Further query of 53 probands in an independent cohort of cardiac channelopathy, using exome sequencing, revealed the presence of 3 out of the 13 P/LP variants. The identification of p.G179Sfs*62, p.R823W and c.420 + 2 T > C variants in KCNQ1, KCNH2 and CASQ2 genes, respectively, validate the significance of the P/LP variants in the context of clinical applicability as well as for large-scale population analysis. CONCLUSION: A compendium of ACMG/AMP classified cardiac channelopathy variants in 1029 self-declared healthy Indian population was created. A conservative genotypic prevalence was estimated to be 0.9-1.8% which poses a huge public health burden for a country with large population size like India. In the majority of cases, these disorders are manageable and the risk of sudden cardiac death can be alleviated by appropriate lifestyle modifications as well as treatment regimens/clinical interventions. Clinical utility of the obtained variants was demonstrated using a cardiac channelopathy patient cohort. Our study emphasises the need for large-scale population screening to identify at-risk individuals and take preventive measures. However, we suggest cautious clinical interpretation to be exercised by taking other cardiac channelopathy risk factors into account.

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.

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.

IndiGenomes: a comprehensive resource of genetic variants from over 1000 Indian genomes.

With the advent of next-generation sequencing, large-scale initiatives for mining whole genomes and exomes have been employed to better understand global or population-level genetic architecture. India encompasses more than 17% of the world population with extensive genetic diversity, but is under-represented in the global sequencing datasets. This gave us the impetus to perform and analyze the whole genome sequencing of 1029 healthy Indian individuals under the pilot phase of the 'IndiGen' program. We generated a compendium of 55,898,122 single allelic genetic variants from geographically distinct Indian genomes and calculated the allele frequency, allele count, allele number, along with the number of heterozygous or homozygous individuals. In the present study, these variants were systematically annotated using publicly available population databases and can be accessed through a browsable online database named as 'IndiGenomes' http://clingen.igib.res.in/indigen/. The IndiGenomes database will help clinicians and researchers in exploring the genetic component underlying medical conditions. Till date, this is the most comprehensive genetic variant resource for the Indian population and is made freely available for academic utility. The resource has also been accessed extensively by the worldwide community since it's launch.

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.

Distribution of gonadotropin-inhibitory hormone (GnIH)-like immunoreactivity in the brain and pituitary of the frog (Pelophylax esculentus) during development.

Gonadotropin-inhibitory hormone (GnIH) is a newly discovered hypothalamic RFamide peptide that influences reproduction by regulating brain and pituitary neuroendocrine functions in vertebrates. We report here for the first time, the ontogenetic description of GnIH-like immunoreactivity in the brain, olfactory system, and pituitary of the frog, Pelophylax esculentus. GnIH-like immunoreactive (GnIH-ir) elements were first observed in larvae at stage 24 in the olfactory mucosa, ventral telencephalon, and diencephalon. GnIH-ir-positive staining progressively increased in frequency and intensity during larval growth and other ir perikarya appeared in the medial septum, anterior commissure, dorsal hypothalamus, and posterior tuberculum. A decline in GnIH-ir neurons was seen along the olfactory/vomeronasal/terminal nerve complex in the stages following the pre- and prometamorphosis, while other GnIH-ir neurons showed positivity in the ventromedial surface of the olfactory bulbs and into the habenular nuclei, but the latter are no longer observed in the following stages of development. The anterior-posterior axis in several brain areas, along with the median eminence and pars intermedia of the hypophysis had the appearance of GnIH-ir fibers from early stages, with a progressive increase in the number till metamorphosis in all major subdivisions of the brain. After premetamorphosis, GnIH-ir fibers arising from labeled neurons in the suprachiasmatic nucleus could be seen contacting the ventricular lumen. The transient appearance of GnIH-ir elements in the olfactory system may hint at an olfactory placode origin in the extracranial region. The distribution of GnIH in several brain regions throughout development suggests important involvement of GnIH in multiple brain functions during development.

Protocol for next-generation sequencing of the LSD virus genome using an amplicon-based approach.

Lumpy skin disease (LSD) is a viral disease predominantly affecting cattle caused by a poxvirus belonging to the capripoxvirus genus. Here, we present a protocol for next-generation sequencing of the LSD virus genome using an amplicon-based approach. We describe steps for DNA extraction, viral DNA enrichment, amplicon pooling and purification, and library preparation and pooling. We then detail procedures for sequencing and computational analysis. This protocol can be adapted to any Illumina sequencing platform as an accelerated and scalable system. For complete details on the use and execution of this protocol, please refer to Bhatt et al.1,2.

Spectrum of rare and common mitochondrial DNA variations from 1029 whole genomes of self-declared healthy individuals from India.

Mitochondrial disorders are a class of heterogeneous disorders caused by genetic variations in the mitochondrial genome (mtDNA) as well as the nuclear genome. The spectrum of mtDNA variants remains unexplored in the Indian population. In the present study, we have cataloged 2689 high confidence single nucleotide variants, small insertions and deletions in mtDNA in 1029 healthy Indian individuals. We found a major proportion (76.5 %) of the variants being rare (AF<=0.005) in the studied population. Intriguingly, we found two 'confirmed' pathogenic variants (m.1555 A>G and m.14484 T>C) with a frequency of ∼1 in 250 individuals in our dataset. The high carrier frequency underscores the need for screening of the mtDNA pathogenic mutations in newborns in India. Interestingly, our analysis also revealed 202 variants in our dataset which have been 'reported' in disease cases as per the MITOMAP database. Additionally, we found the frequency of haplogroup M (52.2 %) to be the highest among all the 18 top-level haplogroups found in our dataset. In comparison to the global population datasets, 20 unique mtDNA variants are found in the Indian population. We hope the whole genome sequencing based compendium of mtDNA variants along with their allele frequencies and heteroplasmy levels in the Indian population will drive additional genome scale studies for mtDNA. Furthermore, the identification of clinically relevant variants in our dataset will aid in better clinical interpretation of the variants in mitochondrial disorders.

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.

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 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.

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.

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.

A Rapid and Scalable Multiplex PCR-Based Next-Generation Amplicon Sequencing Method for Familial Hypercholesterolemia Genetic Screening.

BACKGROUND: Familial hypercholesterolemia (FH) is a frequently underdiagnosed genetic disorder characterized by elevated low-density lipoprotein (LDL) levels. Genetic testing of LDLR, APOB, and PCSK9 genes can identify variants in up to 80% of clinically diagnosed patients. However, limitations in time, scalability, and cost have hindered effective next-generation sequencing of these genes. Additionally, pharmacogenomic variants are associated with statin-induced adverse effects in FH patients. To address these challenges, we developed a multiplex primer-based amplicon sequencing approach for FH genetic testing. METHODS: Multiplex primers were designed for the exons of the LDLR, APOB, and PCSK9 genes, as well as for pharmacogenomic variants rs4149056 (SLCO1B1:c.521T > A), rs2306283 (SLCO1B1:c.388A > G), and rs2231142 (ABCG2:c.421C > A). Analytical validation using samples with known pathogenic variants and clinical validation with 12 FH-suspected probands were conducted. Library preparation was based on a bead-based tagmentation method, and sequencing was conducted on the NovaSeq 6000 platform. RESULTS: Our approach ensured no amplicon dropouts, with over 100× coverage on each amplicon. Known variants in 2 samples were successfully detected. Further, we identified one heterozygous LDLR (p.Glu228Ter) variant and 2 homozygous cases of LDLR (p.Lys294Ter) and LDLR (p.Ser177Leu) variants in patients. Pharmacogenomic analysis revealed that overall 3 patients may require reduced statin doses. Our approach offered reduced library preparation time (approximately 3 h), greater scalability, and lower costs (under $50) for FH genetic testing. CONCLUSIONS: Our method effectively sequences LDLR, APOB, and PCSK9 genes including pharmacogenomic variants that will guide appropriate screening and statin dosing, thus increasing both efficiency and affordability.