Search

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

Abstract Numerous speculations have continually emerged, trying to explore the association between COVID-19 infection and a varied range of demographic and clinical factors. Frontline healthcare workers have been the primary group exposed to this infection, and there have been limited global research that examine this cohort. However, while there are a few large studies conducted on Indian healthcare professionals to investigate their potential risk and predisposing factors to COVID-19 infection, to our knowledge there are no studies evaluating the development of long COVID in this population. This cross-sectional study systematically utilized the demographic and clinical data of 3329 healthcare workers (HCW) from a tertiary hospital in India to gain significant insights into the associations between disease prevalence, severity of SARS-Cov-2 infection and long COVID. Most of the study population was found to be vaccinated (2,615, 78.5%), while 654 (19.65%) HCWs were found to be SARS-CoV-2 positive at least once. Of the infected HCWs, 75.1% (491) did not require hospitalization, whereas the rest were hospitalized for an average duration of 9 days. A total of 206 (6.19%) individuals were found to be suffering from long COVID. Persistent weakness/tiredness was the most experienced long-COVID symptom, while females (1.79, 1.25–2.57), individuals who consumed alcohol (1.85, 1.3–2.64) or had blood group B (1.9, 1.33–2.7) were at a significantly higher risk for developing long COVID.

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.

LncRNA-based control affects cardiac pathophysiologies like myocardial infarction, coronary artery disease, hypertrophy, and myotonic muscular dystrophy. This study used a gene-break transposon (GBT) to screen zebrafish (Danio rerio) for insertional mutagenesis. We identified three insertional mutants where the GBT captured a cardiac gene. One of the adult viable GBT mutants had bradycardia (heart arrhythmia) and enlarged cardiac chambers or hypertrophy; we named it “bigheart.” Bigheart mutant insertion maps to grin2bb or N-methyl D-aspartate receptor (NMDAR2B) gene intron 2 in reverse orientation. Rapid amplification of adjacent cDNA ends analysis suggested a new insertion site transcript in the intron 2 of grin2bb. Analysis of the RNA sequencing of wild-type zebrafish heart chambers revealed a possible new transcript at the insertion site. As this putative lncRNA transcript satisfies the canonical signatures, we called this transcript grin2bb associated RNA transcript (grin2bbART). Using in situ hybridization, we confirmed localized grin2bbART expression in the heart, central nervous system, and muscles in the developing embryos and wild-type adult zebrafish atrium and bulbus arteriosus. The bigheart mutant had reduced Grin2bbART expression. We showed that bigheart gene trap insertion excision reversed cardiac-specific arrhythmia and atrial hypertrophy and restored grin2bbART expression. Morpholino-mediated antisense downregulation of grin2bbART in wild-type zebrafish embryos mimicked bigheart mutants; this suggests grin2bbART is linked to bigheart. Cardiovascular tissues use Grin2bb as a calcium-permeable ion channel. Calcium imaging experiments performed on bigheart mutants indicated calcium mishandling in the heart. The bigheart cardiac transcriptome showed differential expression of calcium homeostasis, cardiac remodeling, and contraction genes. Western blot analysis highlighted Camk2d1 and Hdac1 overexpression. We propose that altered calcium activity due to disruption of grin2bbART, a putative lncRNA in bigheart, altered the Camk2d-Hdac pathway, causing heart arrhythmia and hypertrophy in zebrafish.

BackgroundAdvances 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.MethodsWe 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.ResultsUsing 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.ConclusionWe 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.

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

Abstract Background Retinitis pigmentosa (RP) is a hereditary retinal disease which leads to visual impairment. The onset and progression of RP has physiological consequences that affects the ocular environment. Some of the key non-genetic factors which hasten the retinal degeneration in RP include oxidative stress, hypoxia and ocular inflammation. In this study, we investigated the status of the ocular immune privilege during retinal degeneration and the effect of ocular immune changes on the peripheral immune system in RP. We assessed the peripheral blood mononuclear cell stimulation by retinal antigens and their immune response status in RP patients. Subsequently, we examined alterations in ocular immune privilege machineries which may contribute to ocular inflammation and disease progression in rd1 mouse model. Results In RP patients, we observed a suppressed anti-inflammatory response to self-retinal antigens, thereby indicating a deviated response to self-antigens. The ocular milieu in rd1 mouse model indicated a significant decrease in immune suppressive ligands and cytokine TGF-B1, and higher pro-inflammatory ocular protein levels. Further, blood–retinal-barrier breakdown due to decrease in the expression of tight junction proteins was observed. The retinal breach potentiated pro-inflammatory peripheral immune activation against retinal antigens and caused infiltration of the peripheral immune cells into the ocular tissue. Conclusions Our studies with RP patients and rd1 mouse model suggest that immunological consequences in RP is a contributing factor in the progression of retinal degeneration. The ocular inflammation in the RP alters the ocular immune privilege mechanisms and peripheral immune response. These aberrations in turn create an auto-reactive immune environment and accelerate retinal degeneration.

BackgroundPatients, families, the healthcare system, and society as a whole are all significantly impacted by rare diseases (RDs). According to various classifications, there are currently up to 9,000 different rare diseases that have been recognized, and new diseases are discovered every month. Although very few people are affected by each uncommon disease individually, millions of people are thought to be impacted globally when all these conditions are considered. Therefore, RDs represent an important public health concern. Although crucial for clinical care, early and correct diagnosis is still difficult to achieve in many nations, especially those with low and middle incomes. Consequently, a sizeable amount of the overall burden of RD is attributable to undiagnosed RD (URD). Existing barriers and policy aspects impacting the care of patients with RD and URD remain to be investigated.MethodsTo identify unmet needs and opportunities for patients with URD, the Developing Nations Working Group of the Undiagnosed Diseases Network International (DNWG-UDNI) conducted a survey among its members, who were from 20 different nations. The survey used a mix of multiple choice and dedicated open questions covering a variety of topics. To explore reported needs and analyze them in relation to national healthcare economical aspects, publicly available data on (a) World Bank ranking; (b) Current health expenditure per capita; (c) GDP per capita; (d) Domestic general government health expenditure (% of GDP); and (e) Life expectancy at birth, total (years) were incorporated in our study.ResultsThis study provides an in-depth evaluation of the unmet needs for 20 countries: low-income (3), middle-income (10), and high-income (7). When analyzing reported unmet needs, almost all countries (N = 19) indicated that major barriers still exist when attempting to improve the care of patients with UR and/or URD; most countries report unmet needs related to the availability of specialized care and dedicated facilities. However, while the countries ranked as low income by the World Bank showed the highest prevalence of referred unmet needs across the different domains, no specific trend appeared when comparing the high, upper, and low-middle income nations. No overt trend was observed when separating countries by current health expenditure per capita, GDP per capita, domestic general government health expenditure (% of GDP) and life expectancy at birth, total (years). Conversely, both the GDP and domestic general government health expenditure for each country impacted the presence of ongoing research.ConclusionWe found that policy characteristics varied greatly with the type of health system and country. No overall pattern in terms of referral for unmet needs when separating countries by main economic or health indicators were observed. Our findings highlight the importance of identifying actionable points (e.g., implemented orphan drug acts or registries where not available) in order to improve the care and diagnosis of RDs and URDs on a global scale.

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

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

Introduction: In India, approximately 35,000 new cleft patients are born every year. Many patients receive suboptimum, improper, little, or no treatment. The cumulative burden of cleft care is up to 1 million cases. The spectrum of problems is varied, the caseload is enormous, and the logistics of treatment delivery are complicated. The Indian Council of Medical Research (ICMR) task force project was launched in 2012 to evaluate the status of cleft care in India and develop strategies to provide comprehensive cleft care through a dynamic multidisciplinary and multidimensional tool. ICMR task force project (2012–14) was conducted in Delhi, and the National Capital Region reported that at least 50% of studied cases had complex and multiple treatment needs. The needs identified were related to surgical, orthodontic, dental, ENT and Speech, rehabilitation of mutilated dentition, and various psychological disturbances among patients with Cleft anomalies. Aim: The aim of this study was to develop and test web-based application to create a system for national data of patients with cleft anomalies and digitize the patient records in a standardized preagreed format. Methods and Results: Corresponding to contemporary digital technologies and evolutionary improvements in data collection, web-based data collection instrument, including text, photographs, X-rays, and audio files, was considered the most appropriate. The experts from varied domains in consultation with ICMR and National Informatics Centre evolved a web-based data collection instrument which is named the “IndiCleft tool.” The tool has been tested and used over the years and is presently being upgraded to dynamic version for a national data and patient care registry. Conclusion: The present article describes the process of the development of a “dynamic” web-based data collection instrument. The IndiCleft tool is the national resource on cleft data in India.

The COVID-19 pandemic caused by SARS-CoV-2 has caused millions of infections and deaths worldwide. Limited treatment options and the threat from emerging variants underline the need for novel and widely accessible therapeutics. G-quadruplexes (G4s) are nucleic acid secondary structures known to affect many cellular processes including viral replication and transcription. We identified heretofore not reported G4s with remarkably low mutation frequency across >5 million SARS-CoV-2 genomes. The G4 structure was targeted using FDA-approved drugs that can bind G4s - Chlorpromazine (CPZ) and Prochlorperazine (PCZ). We found significant inhibition in lung pathology and lung viral load of SARS-CoV-2 challenged hamsters when treated with CPZ or PCZ that was comparable to the widely used antiviral drug Remdesivir. In support, in vitro G4 binding, inhibition of reverse transcription from RNA isolated from COVID-infected humans, and attenuated viral replication and infectivity in Vero cell cultures were clear in case of both CPZ and PCZ. Apart from the wide accessibility of CPZ/PCZ, targeting relatively invariant nucleic acid structures poses an attractive strategy against viruses like SARS-CoV-2, which spread fast and accumulate mutations quickly.

IntroductionRare diseases (RD) are a health priority worldwide, overall affecting hundreds of millions of people globally. Early and accurate diagnosis is essential to support clinical care but remains challenging in many countries, especially the low- and medium-income ones. Hence, undiagnosed RD (URD) account for a significant portion of the overall RD burden.MethodsIn October 2020, the Developing Nations Working Group of the Undiagnosed Diseases Network International (DNWG-UDNI) launched a survey among its members, belonging to 20 countries across all continents, to map unmet needs and opportunities for patients with URD. The survey was based on questions with open answers and included eight different domains. Conflicting interpretations were resolved in contact with the partners involved.ResultsAll members responded to the survey. The results indicated that the scientific and medical centers make substantial efforts to respond to the unmet needs of patients. In most countries, there is a high awareness of RD issues. Scarcity of resources was highlighted as a major problem, leading to reduced availability of diagnostic expertise and research. Serious equity in accessibility to services were highlighted both within and between participating countries. Regulatory problems, including securing informed consent, difficulties in sending DNA to foreign laboratories, protection of intellectual property, and conflicts of interest on the part of service providers, remain issues of concern. Finally, most respondents stressed the need to strengthen international cooperation in terms of data sharing, clinical research, and diagnostic expertise for URD patients in low and medium income countries.DiscussionThe survey highlighted that many countries experienced a discrepancy between the growing expertise and scientific value, the level of awareness and commitment on the part of relevant parties, and funding bodies. Country-tailored public health actions, including general syllabus of medical schools and of the education of other health professionals, are needed to reduce such gaps.

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

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.

Introduction: There is inadequate information on infections with the Severe Acute Respiratory Syndrome Virus 2 (SARS-CoV-2) in children. Their clinical, as well as pathological correlation, is poorly understood. In India, children and adolescents account for 12% of all Coronavirus Disease 2019 (COVID-19) cases reported. Children accounted for roughly 11% of those impacted globally last year. However, this year, we are seeing around 20-40% of youngsters in positive instances over the world. Even babies and infants are testing positive for COVID-19, although their illness is under control and seldom becomes fatal. Children aged 5 to 12 years, on the other hand, are at a higher risk. Aim: To study the clinical, pathological and genomic characteristics among children with SARS-CoV-2 infection. Materials and Methods: This cross-sectional study was conducted among 48 paediatric positive patients for SARS-CoV-2 at Government Institute of Medical Sciences, Noida, Uttar Pradesh, and CSIR-Institute of Genomics and Integrative Biology, New Delhi, India, from 1st April 2021 to 31st May 2021. The laboratory testing was done by the real-time Reverse Transcription-Polymerase Chain Reaction (RT-PCR) method. The patients were classified as mild, moderate, severe, or asymptomatic. Their clinical and pathological findings were recorded in the case sheet. Genomic analyses were done for identifying the genetic variant in the nine selected samples. Data entry and analysis were performed using Statistical Package for Social Sciences (SPSS) version 26.0. Chi-square test was used for categorical variables and the t-test was used for continuous variables. Results: The study group has median age of 12 years. Male:Female ratio was 2:3. Most children had acquired infection from the community and 30% had the moderate illness and were admitted. Serum Glutamic-Oxalacetic Transaminase (SGOT) and Glutamic-Pyruvic Transaminase (SGPT) were raised in six patients. Alkaline Phosphatase (ALP) was raised in 21 patients and bilirubin was raised in two patients. The average duration of hospitalisation was six days (range 2-13 days). No mortality among the 48 paediatric patients studied was identified in the hospital. Delta variant (B.1.617.2) was identified in seven patients with D614, P681R, L452R mutations and B.1.617.2 was identified in two patients. Delta variant was present in the paediatric patients, but it did not prolong the hospital stay or cause mortality. Conclusion: The findings of the study suggest that children may be a potential source of infection in the SARS-CoV-2 pandemic while having an asymptomatic to mild illness.

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.

Interpopulation differences in prevalence have been well documented for a number of diseases including genetic diseases. The availability of population scale datasets and well-annotated variant resources would therefore offer a new opportunity to understand the genetic factors linked with such population differences or disparities. In the present manuscript, we evaluated the allele frequencies of genetic variants in inherited blood cancer syndromes in population genomes to understand whether they could provide an insight into these population disparities. We analysed 10 inherited disorders, encompassing 20 different genes. Pathogenic variants were systematically collected, and reclassified using the ACMG & AMP guidelines. Variant prevalence was studied using the Genome Aggregation Database (gnomAD) database, which provided a unique population template to understand prevalence across major global populations. Our analysis highlights high-confidence Pathogenic variants that are significantly enriched across specific subpopulations. To the best of our knowledge, this is the first and comprehensive analysis of genetic variants in inherited blood cancer syndromes in global populations and one of the first to suggest allele specific disparities across global populations.

ATP7A is a critical copper transporter involved in Menkes Disease, Occipital horn Syndrome and X-linked distal spinal muscular atrophy type 3 which are X linked genetic disorders. These are rare diseases and their genetic epidemiology of the diseases is unknown. A number of genetic variants in the genes have been reported in published literature as well as databases, however, understanding the pathogenicity of variants and genetic epidemiology requires the data to be compiled in a unified format. To this end, we systematically compiled genetic variants from published literature and datasets. Each of the variants were systematically evaluated for evidences with respect to their pathogenicity and classified as per the American College of Medical Genetics and the Association of Molecular Pathologists (ACMG-AMP) guidelines into Pathogenic, Likely Pathogenic, Benign, Likely Benign and Variants of Uncertain Significance. Additional integrative analysis of population genomic datasets provides insights into the genetic epidemiology of the disease through estimation of carrier frequencies in global populations. To deliver a mechanistic explanation for the pathogenicity of selected variants, we also performed molecular modeling studies. Our modeling studies concluded that the small structural distortions observed in the local structures of the protein may lead to the destabilization of the global structure. To the best of our knowledge, ATP7A Clinical Genetics Resource is one of the most comprehensive compendium of variants in the gene providing clinically relevant annotations in gene.

An epidemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causing coronavirus diseases (C0VID-19) initially reported in Wuhan, China has rapidly emerged into a global pandemic affecting millions of people worldwide. Molecular detection of SARS-CoV-2 using reverse transcription polymerase chain reaction (RT-PCR) forms the mainstay in screening, diagnosis and epidemiology of the disease. Since the virus evolves by accumulating base substitutions, mutations in the viral genome could possibly affect the accuracy of RT-PCR-based detection assays. The recent availability of genomes of SARS-CoV-2 isolates motivated us to assess the presence and potential impact of variations in target sites of the oligonucleotide primers and probes used in molecular diagnosis. We catalogued a total of 132 primer or probe sequences from literature and data available in the public domain. Our analysis revealed that a total of 5862 unique genetic variants mapped to at least one of the 132 primer or probe binding sites in the genome. A total of 29 unique variants were present in ≥ 1% of genomes from at least one of the continents (Asia, Africa, Australia, Europe, North America, and South America) that mapped to 36 unique primers or probes binding sites. Similarly, a total of 27 primer or probe binding sites had cumulative variants frequency of ≥ 1% in the global SARS-CoV-2 genomes. These included primers or probes sites which are used worldwide for molecular diagnosis as well as approved by national and international agencies. We also found 286 SARS-CoV-2 genomic regions with low variability at a continuous stretch of ≥ 20bps that could be potentially used for primer designing. This highlights the need for sequencing genomes of emerging pathogens to enable evidence-based policies for development and approval of diagnostics.

Outcomes of pallidal stimulation in KMT2B dystonia have been infrequently reported prospectively. We report the six-month outcomes of bilateral GPi DBS in an Asian Indian patient with early-onset generalized dystonia associated with a novel heterozygous variant in the KMT2B gene.

Summary: 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).

Abstract Home to a culturally heterogeneous population, India is also a melting pot of genetic diversity. The population architecture characterized by multiple endogamous groups with specific marriage patterns, including the widely prevalent practice of consanguinity, not only makes the Indian population distinct from rest of the world but also provides a unique advantage and niche to understand genetic diseases. Centuries of genetic isolation of population groups have amplified the founder effects, contributing to high prevalence of recessive alleles, which translates into genetic diseases, including rare genetic diseases in India. Rare genetic diseases are becoming a public health concern in India because a large population size of close to a billion people would essentially translate to a huge disease burden for even the rarest of the rare diseases. Genomics-based approaches have been demonstrated to accelerate the diagnosis of rare genetic diseases and reduce the socio-economic burden. The Genomics for Understanding Rare Diseases: India Alliance Network (GUaRDIAN) stands for providing genomic solutions for rare diseases in India. The consortium aims to establish a unique collaborative framework in health care planning, implementation, and delivery in the specific area of rare genetic diseases. It is a nation-wide collaborative research initiative catering to rare diseases across multiple cohorts, with over 240 clinician/scientist collaborators across 70 major medical/research centers. Within the GUaRDIAN framework, clinicians refer rare disease patients, generate whole genome or exome datasets followed by computational analysis of the data for identifying the causal pathogenic variations. The outcomes of GUaRDIAN are being translated as community services through a suitable platform providing low-cost diagnostic assays in India. In addition to GUaRDIAN, several genomic investigations for diseased and healthy population are being undertaken in the country to solve the rare disease dilemma. In summary, rare diseases contribute to a significant disease burden in India. Genomics-based solutions can enable accelerated diagnosis and management of rare diseases. We discuss how a collaborative research initiative such as GUaRDIAN can provide a nation-wide framework to cater to the rare disease community of India.

Rare monogenic autoinflammatory diseases are a group of recurrent inflammatory genetic disorders caused due to genetic variants in over 37 genes. While a number of these disorders have been identified and reported in Middle Eastern populations, the carrier frequency of these genetic variants in the Middle Eastern population is not known. The availability of whole-genome and exome datasets of over 1,000 individuals from Qatar persuaded us to explore the genetic epidemiology of rare autoinflammatory genetic variants. We have systematically analyzed genetic variants in genome-scale datasets from Qatar with a compendium of variants associated with autoinflammatory diseases. The variants were systematically reclassified according to the American College of Medical Genetics and Genomics guidelines for interpretation of variant pathogenicity. Our analysis identified seven pathogenic and likely pathogenic variants with significant differences in their allele frequencies compared to the global population. The cumulative carrier frequency of these variants was found to be 2.58%. Furthermore, our analysis revealed that five genes, implicated in rare autoinflammatory diseases, were under natural selection. To the best of our knowledge, this is the first and most comprehensive study on the population-scale analysis and genetic epidemiology of genetic variants that cause rare autoinflammatory disease in Middle Eastern populations.

COVID-19, the disease caused by the novel SARS-CoV-2 coronavirus, originated as an isolated outbreak in the Hubei province of China but soon created a global pandemic and is now a major threat to healthcare systems worldwide. Following the rapid human-to-human transmission of the infection, institutes around the world have made efforts to generate genome sequence data for the virus. With thousands of genome sequences for SARS-CoV-2 now available in the public domain, it is possible to analyze the sequences and gain a deeper understanding of the disease, its origin, and its epidemiology. Phylogenetic analysis is a potentially powerful tool for tracking the transmission pattern of the virus with a view to aiding identification of potential interventions. Toward this goal, we have created a comprehensive protocol for the analysis and phylogenetic clustering of SARS-CoV-2 genomes using Nextstrain, a powerful open-source tool for the real-time interactive visualization of genome sequencing data. Approaches to focus the phylogenetic clustering analysis on a particular region of interest are detailed in this protocol.

Coronavirus disease 2019 (COVID-19) rapidly spread from a city in China to almost every country in the world, affecting millions of individuals. The rapid increase in the COVID-19 cases in the state of Kerala in India has necessitated the understanding of SARS-CoV-2 genetic epidemiology. We sequenced 200 samples from patients in Kerala using COVIDSeq protocol amplicon-based sequencing. The analysis identified 166 high-quality single-nucleotide variants encompassing four novel variants and 89 new variants in the Indian isolated SARS-CoV-2. Phylogenetic and haplotype analysis revealed that the virus was dominated by three distinct introductions followed by local spread suggesting recent outbreaks and that it belongs to the A2a clade. Further analysis of the functional variants revealed that two variants in the S gene associated with increased infectivity and five variants mapped in primer binding sites affect the efficacy of RT-PCR. To the best of our knowledge, this is the first and most comprehensive report of SARS-CoV-2 genetic epidemiology from Kerala.

The Arab population encompasses over 420 million people characterized by genetic admixture and a consequent rich genetic diversity. A number of genetic diseases have been reported for the first time from the population. Additionally a high prevalence of some genetic diseases including autosomal recessive disorders such as hemoglobinopathies and familial mediterranean fever have been found in the population and across the region. There is a paucity of databases cataloguing genetic variants of clinical relevance from the population. The availability of such a catalog could have implications in precise diagnosis, genetic epidemiology and prevention of disease. To fill in the gap, we have compiled DALIA, a comprehensive compendium of genetic variants reported in literature and implicated in genetic diseases reported from the Arab population. The database aims to act as an effective resource for population-scale and sub-population specific variant analyses, enabling a ready reference aiding clinical interpretation of genetic variants, genetic epidemiology, as well as facilitating rapid screening and a quick reference for evaluating evidence on genetic diseases.

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.

The rapid emergence of coronavirus disease 2019 (COVID-19) as a global pandemic affecting millions of individuals globally has necessitated sensitive and high-throughput approaches for the diagnosis, surveillance, and determining the genetic epidemiology of SARS-CoV-2. In the present study, we used the COVIDSeq protocol, which involves multiplex-PCR, barcoding, and sequencing of samples for high-throughput detection and deciphering the genetic epidemiology of SARS-CoV-2. We used the approach on 752 clinical samples in duplicates, amounting to a total of 1536 samples which could be sequenced on a single S4 sequencing flow cell on NovaSeq 6000. Our analysis suggests a high concordance between technical duplicates and a high concordance of detection of SARS-CoV-2 between the COVIDSeq as well as RT-PCR approaches. An in-depth analysis revealed a total of six samples in which COVIDSeq detected SARS-CoV-2 in high confidence which were negative in RT-PCR. Additionally, the assay could detect SARS-CoV-2 in 21 samples and 16 samples which were classified inconclusive and pan-sarbeco positive respectively suggesting that COVIDSeq could be used as a confirmatory test. The sequencing approach also enabled insights into the evolution and genetic epidemiology of the SARS-CoV-2 samples. The samples were classified into a total of 3 clades. This study reports two lineages B.1.112 and B.1.99 for the first time in India. This study also revealed 1,143 unique single nucleotide variants and added a total of 73 novel variants identified for the first time. To the best of our knowledge, this is the first report of the COVIDSeq approach for detection and genetic epidemiology of SARS-CoV-2. Our analysis suggests that COVIDSeq could be a potential high sensitivity assay for the detection of SARS-CoV-2, with an additional advantage of enabling the genetic epidemiology of SARS-CoV-2.

[This corrects the article DOI: 10.1371/journal.pone.0060066.].

Hyper-IgD syndrome (HIDS, OMIM #260920) is a rare autosomal recessive autoinflammatory disorder caused by pathogenic variants in the mevalonate kinase (MVK) gene. HIDS has an incidence of 1:50,000 to 1:5,000, and is thought to be prevalent mainly in northern Europe. Here, we report a case series of HIDS from India, which includes ten patients from six families who presented with a wide spectrum of clinical features such as recurrent fever, oral ulcers, rash, arthritis, recurrent diarrhea, hepatosplenomegaly, and high immunoglobulin levels. Using whole exome sequencing (WES) and/or Sanger capillary sequencing, we identified five distinct genetic variants in the MVK gene from nine patients belonging to six families. The variants were classified as pathogenic or likely pathogenic as per the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG-AMP) guidelines for annotation of sequence variants. Over 70% of patients in the present study had two recurrent mutations in MVK gene i.e. a nonsynonymous variant p.V377I, popularly known as the 'Dutch mutation', along with a splicing variant c.226+2delT in a compound heterozygous form. Identity by descent analysis in two patients with the recurrent variants identified a 6.7 MB long haplotype suggesting a founder effect in the South Indian population. Our analysis suggests that a limited number of variants account for the majority of the patients with HIDS in South India. This has implications in clinical diagnosis, as well as in the development of cost-effective approaches for genetic diagnosis and screening. To our best knowledge, this is the first and most comprehensive case series of clinically and genetically characterized patients with HIDS from India.

Recent advances in the field of genomics have seen the successful implementation of whole exome sequencing as a rapid and efficient diagnostic strategy in several genodermatoses. The aim of this study was to explore the potential of molecular studies in dystrophic epidermolysis bullosa (DEB) in India. Whole exome sequencing was performed using genomic DNA from each case of epidermolysis bullosa, followed by massively parallel sequencing. Resulting reads were mapped to the human reference genome hg19. Sanger sequencing subsequently confirmed the potentially pathogenic mutations. Whole exome sequencing of 18 patients with DEB from 17 unrelated Indian families revealed 20 distinct sequence variants in the COL7A1 gene including 2 widely prevalent mutations. Dominant inheritance was seen in 7 patients, while 11 patients showed a highly variable recessive DEB. This preliminary study using exome sequencing is clearly encouraging and will serve as the basis for future large-scale molecular studies to actively identify and understand DEB in the Indian population.

Abstract Background RNA is known to play diverse roles in gene regulation. The clues for this regulatory function of RNA are embedded in its ability to fold into intricate secondary and tertiary structure. Results We report the transcriptome-wide RNA secondary structure in zebrafish at single nucleotide resolution using Parallel Analysis of RNA Structure (PARS). This study provides the secondary structure map of zebrafish coding and non-coding RNAs. The single nucleotide pairing probabilities of 54,083 distinct transcripts in the zebrafish genome were documented. We identified RNA secondary structural features embedded in functional units of zebrafish mRNAs. Translation start and stop sites were demarcated by weak structural signals. The coding regions were characterized by the three-nucleotide periodicity of secondary structure and display a codon base specific structural constrain. The splice sites of transcripts were also delineated by distinct signature signals. Relatively higher structural signals were observed at 3’ Untranslated Regions (UTRs) compared to Coding DNA Sequence (CDS) and 5’ UTRs. The 3′ ends of transcripts were also marked by unique structure signals. Secondary structural signals in long non-coding RNAs were also explored to better understand their molecular function. Conclusions Our study presents the first PARS-enabled transcriptome-wide secondary structure map of zebrafish, which documents pairing probability of RNA at single nucleotide precision. Our findings open avenues for exploring structural features in zebrafish RNAs and their influence on gene expression.

Abstract Background Only a handful of long noncoding RNAs have been functionally characterized. They are known to modulate regulation through interacting with other biomolecules in the cell: DNA, RNA and protein. Though there have been detailed investigations on lncRNA-miRNA and lncRNA-protein interactions, the interaction of lncRNAs with DNA have not been studied extensively. In the present study, we explore whether lncRNAs could modulate genomic regulation by interacting with DNA through the formation of highly stable DNA:DNA:RNA triplexes. Methods We computationally screened 23,898 lncRNA transcripts as annotated by GENCODE, across the human genome for potential triplex forming sequence stretches (PTS). The PTS frequencies were compared across 5′UTR, CDS, 3′UTR, introns, promoter and 1000 bases downstream of the transcription termination sites. These regions were annotated by mapping to experimental regulatory regions, classes of repeat regions and transcription factors. We validated few putative triplex mediated interactions where lncRNA-gene pair interaction is via pyrimidine triplex motif using biophysical methods. Results We identified 20,04,034 PTS sites to be enriched in promoter and intronic regions across human genome. Additional analysis of the association of PTS with core promoter elements revealed a systematic paucity of PTS in all regulatory regions, except TF binding sites. A total of 25 transcription factors were found to be associated with PTS. Using an interaction network, we showed that a subset of the triplex forming lncRNAs, have a positive association with gene promoters. We also demonstrated an in vitro interaction of one lncRNA candidate with its predicted gene target promoter regions. Conclusions Our analysis shows that PTS are enriched in gene promoter and largely associated with simple repeats. The current study suggests a major role of a subset of lncRNAs in mediating chromatin organization modulation through CTCF and NSRF proteins.

Abstract Liver plays a key role in maintaining glucose homeostasis and impaired hepatic glucose metabolism is associated with type 2 diabetes. In the present study, we used RNA sequencing to profile the transcriptome of the livers of diabetic db/db mice as compared to the normal db/+ mice and identified 218 differentially expressed genes. Amongst these, there were 3 lncRNAs that were significantly downregulated and H19 was the most altered lncRNA in the livers of db/db mice. H19 expression significantly correlated with the expression of genes of the glycolysis and gluconeogenesis pathways, which suggest that altered hepatic H19 levels can directly or indirectly modulate their expression. Inhibition of H19 using specific siRNA in HepG2 cells and primary mouse hepatocytes significantly increased the levels of gluconeogenic genes. This was subsequently accompanied by increased hepatic glucose output. Further,H19 depletion in HepG2 cells impaired insulin signaling and increased nuclear localization of FoxO1, an important transcriptional regulator of gluconeogenic gene expression. Our results reveal a novel link between decreased H19 levels and impaired gluconeogenesis via regulation of FoxO1 nuclear levels. These put forth interesting observations on the regulatory role of H19 in altering hepatic physiology during diabetes.

Adolescence is the most critical phase of human growth that radically alters physiology of the body and wherein any inconsistency can lead to serious health consequences in adulthood. The timing and pace at which various developmental events occur during adolescence is highly diverse and thus results in a drastic change in blood biochemistry. Monitoring the physiological levels of various biochemical measures in ample number of individuals during adolescence can call up for an early intervention in managing metabolic diseases in adulthood. Today, only a couple of studies in different populations have investigated blood biochemistry in a small number of adolescents however, there is no comprehensive biochemical data available worldwide. In view, we performed a cross-sectional study in a sizeable group of 7,618 Indian adolescents (3,333 boys and 4,285 girls) aged between 11-17 years to inspect the distribution of values of common biochemical parameters that generally prevails during adolescence and we observed that various parameters considerably follow the reported values from different populations being 3.56-6.49mmol/L (fasting glucose), 10.60-199.48pmol/L (insulin), 0.21-3.22nmol/L (C-peptide), 3.85-6.25% (HbA1c), 2.49-5.54mmol/L (total cholesterol), 1.16-3.69mmol/L (LDL), 0.78-1.85mmol/L (HDL), 0.33-2.24mmol/L (triglycerides), 3.56-11.45mmol/L (urea), 130.01-440.15μmol/L (uric acid) and 22.99-74.28μmol/L (creatinine). Barring LDL and triglycerides, all parameters differed significantly between boys and girls (p< 0.001). Highest difference was seen for uric acid (p = 1.3 x10-187) followed by C-peptide (p = 6.6 x10-89). Across all ages during adolescence, glycemic and nitrogen metabolites parameters varied markedly with gender. Amongst lipid parameters, only HDL levels were found to be significantly associated with gender following puberty (p< 0.001). All parameters except urea, differed considerably in obese and lean adolescents (p< 0.0001). The present study asserts that age, sex and BMI are the essential contributors to variability in blood biochemistry during adolescence. Our composite data on common blood biochemical measures will benefit future endeavors to define reference intervals in adolescents especially when the global availability is scarce.

Background: Long non-coding RNAs (lncRNAs), are being reported to be extensively involved in diverse regulatory roles and have exhibited numerous disease associations. LncRNAs modulate their function through interaction with other biomolecules in the cell including DNA, RNA, and proteins. The availability of genome-scale experimental datasets of RNA binding proteins (RBP) motivated us to understand the role of lncRNAs in terms of its interactions with these proteins. In the current report, we demonstrate a comprehensive study of interactions between RBP and lncRNAs at a transcriptome scale through extensive analysis of the crosslinking and immunoprecipitation (CLIP) experimental datasets available for 70 RNA binding proteins.Results: Our analysis suggests that density of interaction sites for these proteins was significantly higher for specific sub-classes of lncRNAs when compared to protein-coding transcripts. We also observe a positional preference of these RBPs across lncRNA and protein coding transcripts in addition to a significant co-occurrence of RBPs having similar functions, suggesting a modular organization of these elements across lncRNAs.Conclusion: The significant enrichment of RBP sites across some lncRNA classes is suggestive that these interactions might be important in understanding the functional role of lncRNA. We observed a significant enrichment of RBPs which are involved in functional roles such as silencing, splicing, mRNA processing, and transport, indicating the potential participation of lncRNAs in such processes.

Earthworms show a wide spectrum of regenerative potential with certain species like Eisenia fetida capable of regenerating more than two-thirds of their body while other closely related species, such as Paranais litoralis seem to have lost this ability. Earthworms belong to the phylum Annelida, in which the genomes of the marine oligochaete Capitella telata and the freshwater leech Helobdella robusta have been sequenced and studied. Herein, we report the transcriptomic changes in Eisenia fetida (Indian isolate) during regeneration. Following injury, E. fetida regenerates the posterior segments in a time spanning several weeks. We analyzed gene expression changes both in the newly regenerating cells and in the adjacent tissue, at early (15days post amputation), intermediate (20days post amputation) and late (30 days post amputation) by RNAseq based de novo assembly and comparison of transcriptomes. We also generated a draft genome sequence of this terrestrial red worm using short reads and mate-pair reads. An in-depth analysis of the miRNome of the worm showed that many miRNA gene families have undergone extensive duplications. Sox4, a master regulator of TGF-beta mediated epithelial-mesenchymal transition was induced in the newly regenerated tissue. Genes for several proteins such as sialidases and neurotrophins were identified amongst the differentially expressed transcripts. The regeneration of the ventral nerve cord was also accompanied by the induction of nerve growth factor and neurofilament genes. We identified 315 novel differentially expressed transcripts in the transcriptome, that have no homolog in any other species. Surprisingly, 82% of these novel differentially expressed transcripts showed poor potential for coding proteins, suggesting that novel ncRNAs may play a critical role in regeneration of earthworm.

Severe combined immunodeficiency is a large clinically heterogeneous group of disorders caused by a defect in the development of humoral or cellular immune responses. At least 13 genes are known to be involved in the pathophysiology of the disease and the mutation spectrum in SCID has been well documented. Mutations of the recombination-activating genes RAG 1 and RAG 2 are associated with a range of clinical presentations including, severe combined immunodeficiency and autoimmunity. Recently, our understanding of the molecular basis of immune dysfunction in RAG deficiency has improved tremendously with newer insights into the ultrastructure of the RAG complex. In this report, we describe the application of whole exome sequencing for arriving at a molecular diagnosis in a child suffering from B- T- NK+ severe combined immunodeficiency. Apart from making the accurate molecular diagnosis, we also add a genetic variation c.2308G>A p.E770K to the compendium of variations associated with the disease.

Syndromes of mineralocorticoid excess (SME) are closely related clinical manifestations occurring within a specific set of diseases. Overlapping clinical manifestations of such syndromes often create a dilemma in accurate diagnosis, which is crucial for disease surveillance and management especially in rare genetic disorders. Here we demonstrate the use of whole exome sequencing (WES) for accurate diagnosis of rare SME and report that p.R337C variation in the HSD11B2 gene causes progressive apparent mineralocorticoid excess (AME) syndrome in a South Indian family of Mappila origin.

X-linked agammaglobulinemia (XLA) is an extremely rare inherited primary immunodeficiency characterized by recurrent bacterial infections, decrease in number of mature B cells and low serum immunoglobulins. XLA is caused by mutations in the gene encoding Bruton's tyrosine kinase. We report a case of a young Indian boy suspected to have XLA. Immunophenotyping was performed for the affected child using CD20, CD19 and CD3 antibodies. Whole exome sequencing was performed using trio-based approach. The variants were further analyzed using capillary sequencing in the trio as well as maternal grandmother. Initial immunophenotyping in the affected child showed decreased count of CD19+ B cells. To strengthen the clinical findings and confirm the diagnosis of XLA, we performed whole exome sequencing. Our analysis identified a novel frameshift insertion (c.1325dupT) in the BTK gene, which was further validated by Sanger sequencing. Our approach shows the potential in using whole exome sequencing to pinpoint the molecular lesion, enabling timely diagnosis and genetic counseling, and potentially offering prenatal genetic testing for the family.

Pseudomonas aeruginosa (Pa) is the leading cause of chronic lung infection in Cystic Fibrosis (CF) patients. It is well recognized that CF epithelial cells fail to develop an appropriate response to infection, allowing bacterial colonization and a chronic inflammatory response. Since long non-coding RNAs (lncRNAs), are known to play a key role in regulating mammalian innate immune response, we hypothesized that CF cells exposed to Pa could express a specific lncRNA signature responsible of the maladaptative CF response. We analyzed transcriptomic datasets to compare the expression profiles of lncRNAs in primary CF and non-CF epithelial cells infected with Pa at 0, 2, 4, and 6 h of infection. Our analysis identified temporal expression signatures of 25, 73, 15, and 26 lncRNA transcripts differentially expressed at 0, 2, 4, and 6 h post-infection respectively, between CF and non-CF cells. In addition, we identified profiles specific to CF and non-CF cells. The differential expression of two candidate lncRNAs were independently validated using real-time PCR. We identified a specific CF signature of lncRNA expression in a context of Pa infection that could potentially play a role in the maladaptive immune response of CF patients.

[This corrects the article DOI: 10.1371/journal.pgen.1005954.].

Dystrophic epidermolysis bullosa simplex (DEB) is a phenotypically diverse inherited skin fragility disorder. It is majorly manifested by appearance of epidermal bullae upon friction caused either by physical or environmental trauma. The phenotypic manifestations also include appearance of milia, scarring all over the body and nail dystrophy. DEB can be inherited in a recessive or dominant form and the recessive form of DEB (RDEB) is more severe. In the present study, we identify a novel p.G2254fs mutation in COL7A1 gene causing a sporadic case of RDEB by whole exome sequencing (WES). Apart from adding a novel frameshift Collagen VII mutation to the repertoire of known mutations reported in the disease, to the best of our knowledge, this is the first report of a genetically characterized case of DEB from India.

Dystrophic epidermolysis bullosa simplex (DEB) is a phenotypically diverse inherited skin fragility disorder. It is majorly manifested by appearance of epidermal bullae upon friction caused either by physical or environmental trauma. The phenotypic manifestations also include appearance of milia, scarring all over the body and nail dystrophy. DEB can be inherited in a recessive or dominant form and the recessive form of DEB (RDEB) is more severe. In the present study, we identify a novel p.G2254fs mutation in COL7A1 gene causing a sporadic case of RDEB by whole exome sequencing (WES). Apart from adding a novel frameshift Collagen VII mutation to the repertoire of known mutations reported in the disease, to the best of our knowledge, this is the first report of a genetically characterized case of DEB from India.

We report here the ~670 Mb genome assembly of the Asian seabass (Lates calcarifer), a tropical marine teleost. We used long-read sequencing augmented by transcriptomics, optical and genetic mapping along with shared synteny from closely related fish species to derive a chromosome-level assembly with a contig N50 size over 1 Mb and scaffold N50 size over 25 Mb that span ~90% of the genome. The population structure of L. calcarifer species complex was analyzed by re-sequencing 61 individuals representing various regions across the species' native range. SNP analyses identified high levels of genetic diversity and confirmed earlier indications of a population stratification comprising three clades with signs of admixture apparent in the South-East Asian population. The quality of the Asian seabass genome assembly far exceeds that of any other fish species, and will serve as a new standard for fish genomics.

The organization of structure and function of cardiac chambers in vertebrates is defined by chamber-specific distinct gene expression. This peculiarity and uniqueness of the genetic signatures demonstrates functional resolution attributed to the different chambers of the heart. Altered expression of the cardiac chamber genes can lead to individual chamber related dysfunctions and disease patho-physiologies. Information on transcriptional repertoire of cardiac compartments is important to understand the spectrum of chamber specific anomalies. We have carried out a genome wide transcriptome profiling study of the three cardiac chambers in the zebrafish heart using RNA sequencing. We have captured the gene expression patterns of 13,396 protein coding genes in the three cardiac chambers-atrium, ventricle and bulbus arteriosus. Of these, 7,260 known protein coding genes are highly expressed (≥10 FPKM) in the zebrafish heart. Thus, this study represents nearly an all-inclusive information on the zebrafish cardiac transcriptome. In this study, a total of 96 differentially expressed genes across the three cardiac chambers in zebrafish were identified. The atrium, ventricle and bulbus arteriosus displayed 20, 32 and 44 uniquely expressing genes respectively. We validated the expression of predicted chamber-restricted genes using independent semi-quantitative and qualitative experimental techniques. In addition, we identified 23 putative novel protein coding genes that are specifically restricted to the ventricle and not in the atrium or bulbus arteriosus. In our knowledge, these 23 novel genes have either not been investigated in detail or are sparsely studied. The transcriptome identified in this study includes 68 differentially expressing zebrafish cardiac chamber genes that have a human ortholog. We also carried out spatiotemporal gene expression profiling of the 96 differentially expressed genes throughout the three cardiac chambers in 11 developmental stages and 6 tissue types of zebrafish. We hypothesize that clustering the differentially expressed genes with both known and unknown functions will deliver detailed insights on fundamental gene networks that are important for the development and specification of the cardiac chambers. It is also postulated that this transcriptome atlas will help utilize zebrafish in a better way as a model for studying cardiac development and to explore functional role of gene networks in cardiac disease pathogenesis.

An increase in the expression of estrogen receptors (ER) and the expanded population of ER-positive cells are two common phenotypes of breast cancer. Detection of the aberrantly expressed ERα in breast cancer is carried out using ERα-antibodies and radiolabelled ligands to make decisions about cancer treatment and targeted therapy. Capitalizing on the beneficial advantages of aptamer over the conventional antibody or radiolabelled ligand, we have identified a DNA aptamer that selectively binds and facilitates the detection of ERα in human breast cancer tissue sections. The aptamer is identified using the high throughput sequencing assisted SELEX screening. Biophysical characterization confirms the binding and formation of a thermodynamically stable complex between the identified DNA aptamer (ERaptD4) and ERα (Ka = 1.55±0.298×108 M(-1); ΔH = 4.32×104±801.1 cal/mol; ΔS = -108 cal/mol/deg). Interestingly, the specificity measurements suggest that the ERaptD4 internalizes into ERα-positive breast cancer cells in a target-selective manner and localizes specifically in the nuclear region. To harness these characteristics of ERaptD4 for detection of ERα expression in breast cancer samples, we performed the aptamer-assisted histochemical analysis of ERα in tissue samples from breast cancer patients. The results were validated by performing the immunohistochemistry on same samples with an ERα-antibody. We found that the two methods agree strongly in assay output (kappa value = 0.930, p-value