De novo identification of viral pathogens from cell culture hologenomes

<p>Abstract</p> <p>Background</p> <p>Fast, specific identification and surveillance of pathogens is the cornerstone of any outbreak response system, especially in the case of emerging infectious diseases and viral epidemics. This process is generally tedious and time-consuming thus making it ineffective in traditional settings. The added complexity in these situations is the non-availability of pure isolates of pathogens as they are present as mixed genomes or hologenomes. Next-generation sequencing approaches offer an attractive solution in this scenario as it provides adequate depth of sequencing at fast and affordable costs, apart from making it possible to decipher complex interactions between genomes at a scale that was not possible before. The widespread application of next-generation sequencing in this field has been limited by the non-availability of an efficient computational pipeline to systematically analyze data to delineate pathogen genomes from mixed population of genomes or hologenomes.</p> <p>Findings</p> <p>We applied next-generation sequencing on a sample containing mixed population of genomes from an epidemic with appropriate processing and enrichment. The data was analyzed using an extensive computational pipeline involving mapping to reference genome sets and <it>de-novo </it>assembly. In depth analysis of the data generated revealed the presence of sequences corresponding to <it>Japanese encephalitis </it>virus. The genome of the virus was also independently <it>de-novo </it>assembled. The presence of the virus was in addition, verified using standard molecular biology techniques.</p> <p>Conclusions</p> <p>Our approach can accurately identify causative pathogens from cell culture hologenome samples containing mixed population of genomes and in principle can be applied to patient hologenome samples without any background information. This methodology could be widely applied to identify and isolate pathogen genomes and understand their genomic variability during outbreaks.</p

Forward genetic screen using a gene-breaking trap approach identifies a novel role of grin2bb-associated RNA transcript (grin2bbART) in zebrafish heart function

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

High Resolution Methylome Map of Rat Indicates Role of Intragenic DNA Methylation in Identification of Coding Region

DNA methylation is crucial for gene regulation and maintenance of genomic stability. Rat has been a key model system in understanding mammalian systemic physiology, however detailed rat methylome remains uncharacterized till date. Here, we present the first high resolution methylome of rat liver generated using Methylated DNA immunoprecipitation and high throughput sequencing (MeDIP-Seq) approach. We observed that within the DNA/RNA repeat elements, simple repeats harbor the highest degree of methylation. Promoter hypomethylation and exon hypermethylation were common features in both RefSeq genes and expressed genes (as evaluated by proteomic approach). We also found that although CpG islands were generally hypomethylated, about 6% of them were methylated and a large proportion (37%) of methylated islands fell within the exons. Notably, we obeserved significant differences in methylation of terminal exons (UTRs); methylation being more pronounced in coding/partially coding exons compared to the non-coding exons. Further, events like alternate exon splicing (cassette exon) and intron retentions were marked by DNA methylation and these regions are retained in the final transcript. Thus, we suggest that DNA methylation could play a crucial role in marking coding regions thereby regulating alternative splicing. Apart from generating the first high resolution methylome map of rat liver tissue, the present study provides several critical insights into methylome organization and extends our understanding of interplay between epigenome, gene expression and genome stability

Silencing c-MYC expression by targeting quadruplex in P1 promoter using locked nucleic acid trap

The nuclease hypersensitive element of P1 promoter in c-MYC gene harbors a potential of unusual structure called quadruplex, which is involved in molecular recognition and function. This Hoogsteen bonded structure is in dynamic equilibrium with the usual Watson−Crick duplex structure, and these competing secondary structures undergo interconversion for execution of their respective biological roles. Herein, we investigate the sensitivity of the c-MYC quadruplex−duplex equilibrium by employing a locked nucleic acid (LNA) modified complementary strand as a pharmacological agent. Our biophysical experiments indicate that the c-MYC quadruplex under physiological conditions is stable and dominates the quadruplex−WC duplex equilibrium in both sodium and potassium buffers. This equilibrium is perturbed upon introducing the LNA modified complementary strand, which demonstrates efficient invasion of stable c-MYC quadruplex and duplex formation in contrast to the unmodified complementary strand. Our data indicate that LNA modifications confer increased thermodynamic stability to the duplex and thus favor the predominance of the duplex population over that of the quadruplex. Further, we demonstrate that this perturbation of equilibrium by a pharmacological agent results in altered gene expression. Our in vivo experiment performed using the LNA modified complementary strand suggests the influence of the quadruplex−duplex structural switch in the modulation of gene expression. We believe that this exploratory approach utilizing the selectivity and specificity of Watson−Crick base pairing of LNA bases would allow the modulation of quadruplex regulated gene expression

Human 45,X fibroblast transcriptome reveals distinct differentially expressed genes including long noncoding RNAs potentially associated with the pathophysiology of Turner syndrome.

Turner syndrome is a chromosomal abnormality characterized by the absence of whole or part of the X chromosome in females. This X aneuploidy condition is associated with a diverse set of clinical phenotypes such as gonadal dysfunction, short stature, osteoporosis and Type II diabetes mellitus, among others. These phenotypes differ in their severity and penetrance among the affected individuals. Haploinsufficiency for a few X linked genes has been associated with some of these disease phenotypes. RNA sequencing can provide valuable insights to understand molecular mechanism of disease process. In the current study, we have analysed the transcriptome profiles of human untransformed 45,X and 46,XX fibroblast cells and identified differential expression of genes in these two karyotypes. Functional analysis revealed that these differentially expressing genes are associated with bone differentiation, glucose metabolism and gonadal development pathways. We also report differential expression of lincRNAs in X monosomic cells. Our observations provide a basis for evaluation of cellular and molecular mechanism(s) in the establishment of Turner syndrome phenotypes

Validation of RNA-seq results by qRT-PCR.

<p>In total, five genes including XIST (lincRNA) were selected for real time data validation. Relative expression and FPKM values (log10) are shown for 45,X and 46,XX conditions. qRT-PCR results are found to correlate with RNA seq data, hence further confirming the expression dataset.</p

Differentially expressed long non-coding RNAs in 45,X and 46,XX fibroblast cells at P value<0.001 and FDR (q) value<0.05.

<p>Differentially expressed long non-coding RNAs in 45,X and 46,XX fibroblast cells at P value<0.001 and FDR (q) value<0.05.</p

Heat Map depicting the expression levels 116 significantly differentially expressed genes in 45,X and 46,XX fibroblast cells.

<p>Heat Map depicting the expression levels 116 significantly differentially expressed genes in 45,X and 46,XX fibroblast cells.</p

Summary of RNA-seq data generated and mapping statistics.

<p>Total number of reads generated and mapped after quality filtering for both 45,X and 46,XX human fibroblast cells.</p