SCAR-6 elncRNA locus epigenetically regulates PROZ and modulates coagulation and vascular function.

In this study, we characterize a novel lncRNA-producing gene locus that we name Syntenic Cardiovascular Conserved Region-Associated lncRNA-6 (scar-6) and functionally validate its role in coagulation and cardiovascular function. A 12-bp deletion of the scar-6 locus in zebrafish (scar-6gib007Δ12/Δ12) results in cranial hemorrhage and vascular permeability. Overexpression, knockdown and rescue with the scar-6 lncRNA modulates hemostasis in zebrafish. Molecular investigation reveals that the scar-6 lncRNA acts as an enhancer lncRNA (elncRNA), and controls the expression of prozb, an inhibitor of factor Xa, through an enhancer element in the scar-6 locus. The scar-6 locus suppresses loop formation between prozb and scar-6 sequences, which might be facilitated by the methylation of CpG islands via the prdm14-PRC2 complex whose binding to the locus might be stabilized by the scar-6 elncRNA transcript. Binding of prdm14 to the scar-6 locus is impaired in scar-6gib007Δ12/Δ12 zebrafish. Finally, activation of the PAR2 receptor in scar-6gib007Δ12/Δ12 zebrafish triggers NF-κB-mediated endothelial cell activation, leading to vascular dysfunction and hemorrhage. We present evidence that the scar-6 locus plays a role in regulating the expression of the coagulation cascade gene prozb and maintains vascular homeostasis.

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.

Mitochondrial calcium signaling mediated transcriptional regulation of keratin filaments is a critical determinant of melanogenesis.

Mitochondria are versatile organelles that regulate several physiological functions. Many mitochondria-controlled processes are driven by mitochondrial Ca2+ signaling. However, role of mitochondrial Ca2+ signaling in melanosome biology remains unknown. Here, we show that pigmentation requires mitochondrial Ca2+ uptake. In vitro gain and loss of function studies demonstrated that Mitochondrial Ca2+ Uniporter (MCU) is crucial for melanogenesis while the MCU rheostats, MCUb and MICU1 negatively control melanogenesis. Zebrafish and mouse models showed that MCU plays a vital role in pigmentation in vivo. Mechanistically, MCU controls activation of transcription factor NFAT2 to induce expression of three keratins (keratin 5, 7 and 8), which we report as positive regulators of melanogenesis. Interestingly, keratin 5 in turn modulates mitochondrial Ca2+ uptake thereby this signaling module acts as a negative feedback loop that fine-tunes both mitochondrial Ca2+ signaling and melanogenesis. Mitoxantrone, an FDA approved drug that inhibits MCU, decreases physiological melanogenesis. Collectively, our data demonstrates a critical role for mitochondrial Ca2+ signaling in vertebrate pigmentation and reveal the therapeutic potential of targeting MCU for clinical management of pigmentary disorders. Given the centrality of mitochondrial Ca2+ signaling and keratin filaments in cellular physiology, this feedback loop may be functional in a variety of other pathophysiological conditions.

G4-binding drugs, chlorpromazine and prochlorperazine, repurposed against COVID-19 infection in hamsters.

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.

A genome-wide circular RNA transcriptome in rat.

Circular RNAs (circRNAs) are a novel class of noncoding RNAs that back-splice from 5' donor site and 3' acceptor sites to form a circular structure. A number of circRNAs have been discovered in model organisms including human, mouse, Drosophila, among other organisms. There are a few candidate-based studies on circRNAs in rat, a well-studied model organism as well. A number of pipelines have been published to identify the back splice junctions for the discovery of circRNAs but studies comparing these tools have suggested that a combination of tools would be a better approach to identify high-confidence circRNAs. The availability of a recent dataset of transcriptomes encompassing 11 tissues, 4 developmental stages, and 2 genders motivated us to explore the landscape of circRNAs in the organism in this context. In order to understand the difference among different pipelines, we employed five different combinations of tools to identify circular RNAs from the dataset. We compared the results of the different combination of tools/pipelines with respect to alignment, total number of circRNAs identified and read-coverage. In addition, we identified tissue-specific, development-stage specific and gender-specific circRNAs and further independently validated 16 circRNA junctions out of 24 selected candidates in 5 tissue samples and estimated the quantitative expression of five circRNA candidates using real-time polymerase chain reaction and our analysis suggests three candidates as tissue-enriched. This study is one of the most comprehensive studies which provides a map of circRNAs transcriptome as well as to understand the difference among different computational pipelines in rat.

Functional long non-coding and circular RNAs in zebrafish.

The utility of model organisms to understand the function of a novel transcript/genes has allowed us to delineate their molecular mechanisms in maintaining cellular homeostasis. Organisms such as zebrafish have contributed a lot in the field of developmental and disease biology. Attributable to advancement and deep transcriptomics, many new transcript isoforms and non-coding RNAs such as long noncoding RNA (lncRNA) and circular RNAs (circRNAs) have been identified and cataloged in multiple databases and many more are yet to be identified. Various methods and tools have been utilized to identify lncRNAs/circRNAs in zebrafish using deep sequencing of transcriptomes as templates. Functional analysis of a few candidates such as tie1-AS, ECAL1 and CDR1as in zebrafish provides a prospective outline to approach other known or novel lncRNA/circRNA. New genetic alteration tools like TALENS and CRISPRs have helped in probing for the molecular function of lncRNA/circRNA in zebrafish. Further latest improvements in experimental and computational techniques offer the identification of lncRNA/circRNA counterparts in humans and zebrafish thereby allowing easy modeling and analysis of function at cellular level.

Rapid and accurate nucleobase detection using FnCas9 and its application in COVID-19 diagnosis.

Rapid detection of DNA/RNA pathogenic sequences or variants through point-of-care diagnostics is valuable for accelerated clinical prognosis, as witnessed during the recent COVID-19 outbreak. Traditional methods relying on qPCR or sequencing are tough to implement with limited resources, necessitating the development of accurate and robust alternative strategies. Here, we report FnCas9 Editor Linked Uniform Detection Assay (FELUDA) that utilizes a direct Cas9 based enzymatic readout for detecting nucleobase and nucleotide sequences without trans-cleavage of reporter molecules. We also demonstrate that FELUDA is 100% accurate in detecting single nucleotide variants (SNVs), including heterozygous carriers, and present a simple web-tool JATAYU to aid end-users. FELUDA is semi-quantitative, can adapt to multiple signal detection platforms, and deploy for versatile applications such as molecular diagnosis during infectious disease outbreaks like COVID-19. Employing a lateral flow readout, FELUDA shows 100% sensitivity and 97% specificity across all ranges of viral loads in clinical samples within 1hr. In combination with RT-RPA and a smartphone application True Outcome Predicted via Strip Evaluation (TOPSE), we present a prototype for FELUDA for CoV-2 detection closer to home.

Initial Insights Into the Genetic Epidemiology of SARS-CoV-2 Isolates From Kerala Suggest Local Spread From Limited Introductions.

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.

High throughput detection and genetic epidemiology of SARS-CoV-2 using COVIDSeq next-generation sequencing.

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.

Methods for Annotation and Validation of Circular RNAs from RNAseq Data.

Circular RNAs are an emerging class of transcript isoforms created by unique back splicing of exons to form a closed covalent circular structure. While initially considered as product of aberrant splicing, recent evidence suggests unique functions and conservation across evolution. While circular RNAs could be largely attributed to have little or no potential to encode for proteins, recent evidence points to at least a small subset of circular RNAs which encode for peptides. Circular RNAs are also increasingly shown to be biomarkers for a number of diseases including neurological disorders and cancer. The advent of deep sequencing has enabled large-scale identification of circular RNAs in human and other genomes. A number of computational approaches have come up in recent years to query circular RNAs on a genome-wide scale from RNA-seq data. In this chapter, we describe the application and methodology of identifying circular RNAs using three popular computational tools: FindCirc, Segemehl, and CIRI along with approaches for experimental validation of the unique splice junctions.

A genome-wide map of circular RNAs in adult zebrafish.

Circular RNAs (circRNAs) are transcript isoforms generated by back-splicing of exons and circularisation of the transcript. Recent genome-wide maps created for circular RNAs in humans and other model organisms have motivated us to explore the repertoire of circular RNAs in zebrafish, a popular model organism. We generated RNA-seq data for five major zebrafish tissues - Blood, Brain, Heart, Gills and Muscle. The repertoire RNA sequence reads left over after reference mapping to linear transcripts were used to identify unique back-spliced exons utilizing a split-mapping algorithm. Our analysis revealed 3,428 novel circRNAs in zebrafish. Further in-depth analysis suggested that majority of the circRNAs were derived from previously well-annotated protein-coding and long noncoding RNA gene loci. In addition, many of the circular RNAs showed extensive tissue specificity. We independently validated a subset of circRNAs using polymerase chain reaction (PCR) and divergent set of primers. Expression analysis using quantitative real time PCR recapitulate selected tissue specificity in the candidates studied. This study provides a comprehensive genome-wide map of circular RNAs in zebrafish tissues.

Methods to Study Long Noncoding RNA Expression and Dynamics in Zebrafish Using RNA Sequencing.

Long noncoding RNAs (lncRNAs) belong to a class of RNA transcripts that do not have the potential to code for proteins. LncRNAs were largely discovered in the transcriptomes of human and several model organisms, using next-generation sequencing (NGS) approaches, which have enabled a comprehensive genome scale annotation of transcripts. LncRNAs are known to have dynamic expression status and have the potential to orchestrate gene regulation at the epigenetic, transcriptional, and posttranscriptional levels. Here we describe the experimental methods involved in the discovery of lncRNAs from the transcriptome of a popular model organism zebrafish (Danio rerio). A structured and well-designed computational analysis pipeline subsequent to the RNA sequencing can be instrumental in revealing the diversity of the lncRNA transcripts. We describe one such computational pipeline used for the discovery of novel lncRNA transcripts in zebrafish. We also detail the validation of the putative novel lncRNA transcripts using qualitative and quantitative assays in zebrafish.

A G-quadruplex motif at the 3' end of sgRNAs improves CRISPR-Cas9 based genome editing efficiency.

Originating as a component of prokaryotic adaptive immunity, the type II CRISPR/Cas9 system has been repurposed for targeted genome editing in various organisms. Although Cas9 can bind and cleave DNA efficiently under in vitro conditions, its activity inside a cell can vary dramatically between targets owing to the differences between genomic loci and the availability of enough Cas9/sgRNA (single guide RNA) complex molecules for cleavage. Most methods have so far relied on Cas9 protein engineering or base modifications in the sgRNA sequence to improve CRISPR/Cas9 activity. Here we demonstrate that a structure based rational design of sgRNAs can enhance the efficiency of Cas9 cleavage in vivo. By appending a naturally forming RNA G-quadruplex motif to the 3' end of sgRNAs we can improve its stability and target cleavage efficiency in zebrafish embryos without inducing off-target activity, thereby underscoring its value in the design of better and optimized genome editing triggers.

SAGE: a comprehensive resource of genetic variants integrating South Asian whole genomes and exomes.

South Asia is home to $\sim $20% of the world population and characterized by distinct ethnic, linguistic, cultural and genetic lineages. Only limited representative samples from the region have found its place in large population-scale international genome projects. The recent availability of genome scale data from multiple populations and datasets from South Asian countries in public domain motivated us to integrate the data into a comprehensive resource. In the present study, we have integrated a total of six datasets encompassing 1213 human exomes and genomes to create a compendium of 154 814 557 genetic variants and adding a total of 69 059 255 novel variants. The variants were systematically annotated using public resources and along with the allele frequencies are available as a browsable-online resource South Asian genomes and exomes. As a proof of principle application of the data and resource for genetic epidemiology, we have analyzed the pathogenic genetic variants causing retinitis pigmentosa. Our analysis reveals the genetic landscape of the disease and suggests subset of genetic variants to be highly prevalent in South Asia.

Pharmacogenetic landscape of DPYD variants in south Asian populations by integration of genome-scale data.

AIM: Adverse drug reactions to 5-Fluorouracil(5-FU) is frequent and largely attributable to genetic variations in the DPYD gene, a rate limiting enzyme that clears 5-FU. The study aims at understanding the pharmacogenetic landscape of DPYD variants in south Asian populations. MATERIALS & METHODS: Systematic analysis of population scale genome wide datasets of over 3000 south Asians was performed. Independent evaluation was performed in a small cohort of patients. RESULTS: Our analysis revealed significant differences in the the allelic distribution of variants in different ethnicities. CONCLUSIONS: This is the first and largest genetic map the DPYD variants associated with adverse drug reaction to 5-FU in south Asian population. Our study highlights ethnic differences in allelic frequencies.

The Zebrafish GenomeWiki: a crowdsourcing approach to connect the long tail for zebrafish gene annotation.

A large repertoire of gene-centric data has been generated in the field of zebrafish biology. Although the bulk of these data are available in the public domain, most of them are not readily accessible or available in nonstandard formats. One major challenge is to unify and integrate these widely scattered data sources. We tested the hypothesis that active community participation could be a viable option to address this challenge. We present here our approach to create standards for assimilation and sharing of information and a system of open standards for database intercommunication. We have attempted to address this challenge by creating a community-centric solution for zebrafish gene annotation. The Zebrafish GenomeWiki is a 'wiki'-based resource, which aims to provide an altruistic shared environment for collective annotation of the zebrafish genes. The Zebrafish GenomeWiki has features that enable users to comment, annotate, edit and rate this gene-centric information. The credits for contributions can be tracked through a transparent microattribution system. In contrast to other wikis, the Zebrafish GenomeWiki is a 'structured wiki' or rather a 'semantic wiki'. The Zebrafish GenomeWiki implements a semantically linked data structure, which in the future would be amenable to semantic search. Database URL: http://genome.igib.res.in/twiki.