Potential pitfalls in the accuracy of analysis of natural sense-antisense RNA pairs by reverse transcription-PCR

Machine Learning-Assisted Direct RNA Sequencing with Epigenetic RNA Modification Detection via Quantum Tunneling

RNA sequence information holds immense potential as a drug target for diagnosing various RNA-mediated diseases and viral/bacterial infections. Massively parallel complementary DNA (c-DNA) sequencing helps but results in a loss of valuable information about RNA modifications, which are often associated with cancer evolution. Herein, we report machine learning (ML)-assisted high throughput RNA sequencing with inherent RNA modification detection using a single-molecule, long-read, and label-free quantum tunneling approach. The ML tools achieve classification accuracy as high as 100% in decoding RNA and 98% for decoding both RNA and RNA modifications simultaneously. The relationships between input features and target values have been well examined through Shapley additive explanations. Furthermore, transmission and sensitivity readouts enable the recognition of RNA and its modifications with good selectivity and sensitivity. Our approach represents a starting point for ML-assisted direct RNA sequencing that can potentially decode RNA and its epigenetic modifications at the molecular level.

Leptospira transcriptome sequencing using long-read technology reveals unannotated transcripts and potential polyadenylation of RNA molecules

IMPORTANCE

Leptospirosis, caused by the spirochete bacteria Leptospira, is a zoonotic disease of humans and animals, accounting for over 1 million annual human cases and over 60,000 deaths. We have characterized operon transcriptional units, identified novel RNA coding regions, and reported evidence of potential posttranscriptional polyadenylation in the Leptospira transcriptomes for the first time using Oxford Nanopore Technology RNA sequencing protocols. The newly identified RNA coding regions and operon transcriptional units were detected only in the pathogenic Leptospira transcriptomes, suggesting their significance in virulence-related functions. This article integrates bioinformatics, infectious diseases, microbiology, molecular biology, veterinary sciences, and public health. Given the current knowledge gap in the regulation of leptospiral pathogenicity, our findings offer valuable insights to researchers studying leptospiral pathogenicity and provide both a basis and a tool for researchers focusing on prokaryotic molecular studies for the understanding of RNA compositions and prokaryotic polyadenylation for their organisms of interest.

Three Novel Antisense Overlapping Genes in E. coli O157:H7 EDL933

The abundance of long overlapping genes in prokaryotic genomes is likely to be significantly underestimated. To date, only a few examples of such genes are fully established. Using RNA sequencing and ribosome profiling, we found expression of novel overlapping open reading frames in Escherichia coli O157:H7 EDL933 (EHEC). Indeed, the overlapping candidate genes are equipped with typical structural elements required for transcription and translation, i.e., promoters, transcription start sites, as well as terminators, all of which were experimentally verified. Translationally arrested mutants, unable to produce the overlapping encoded protein, were found to have a growth disadvantage when grown competitively against the wild type. Thus, the phenotypes found imply biological functionality of the genes at the level of proteins produced. The addition of 3 more examples of prokaryotic overlapping genes to the currently limited, yet constantly growing pool of such genes emphasizes the underestimated coding capacity of bacterial genomes. IMPORTANCE The abundance of long overlapping genes in prokaryotic genomes is likely to be significantly underestimated, since such genes are not allowed in genome annotations. However, ribosome profiling catches mRNA in the moment of being template for protein production. Using this technique and subsequent experiments, we verified 3 novel overlapping genes encoded in antisense of known genes. This adds more examples of prokaryotic overlapping genes to the currently limited, yet constantly growing pool of such genes.

Epigenetic Regulation of HIV-1 Sense and Antisense Transcription in Response to Latency-Reversing Agents

Nucleosomes positioned on the HIV-1 5' long terminal repeat (LTR) regulate sense transcription as well as the establishment and maintenance of latency. A negative-sense promoter (NSP) in the 3' LTR expresses antisense transcripts with coding and non-coding activities. Previous studies identified cis-acting elements that modulate NSP activity. Here, we used the two chronically infected T cell lines, ACH-2 and J1.1, to investigate epigenetic regulation of NSP activity. We found that histones H3 and H4 are present on the 3' LTR in both cell lines. Following treatment with histone deacetylase inhibitors (HDACi), the levels of H3K27Ac increased and histone occupancy declined. HDACi treatment also led to increased levels of RNA polymerase II (RNPII) at NSP, and antisense transcription was induced with similar kinetics and to a similar extent as 5' LTR-driven sense transcription. We also detected H3K9me2 and H3K27me3 on NSP, along with the enzymes responsible for these epigenetic marks, namely G9a and EZH2, respectively. Treatment with their respective inhibitors had little or no effect on RNPII occupancy at the two LTRs, but it induced both sense and antisense transcription. Moreover, the increased expression of antisense transcripts in response to treatment with a panel of eleven latency-reversing agents closely paralleled and was often greater than the effect on sense transcripts. Thus, HIV-1 sense and antisense RNA expression are both regulated via acetylation and methylation of lysine 9 and 27 on histone H3. Since HIV-1 antisense transcripts act as non-coding RNAs promoting epigenetic silencing of the 5' LTR, our results suggest that the limited efficacy of latency-reversing agents in the context of 'shock and kill' cure strategies may be due to concurrent induction of antisense transcripts thwarting their effect on sense transcription.

Trunk injection delivery of dsRNA for RNAi-based pest control in apple trees

BACKGROUND

RNA interference (RNAi) is a promising new approach for controlling insect pests without the use of synthetic pesticides. Trunk injection is a delivery system for woody plants that harnesses the vascular system of the tree to transport materials to the tree canopy. Full size apple trees were injected with double-stranded RNA (dsRNA), and season-long leaf samples were taken to measure the vascular mobility and temporal persistence of dsRNA, using quantitative reverse transcription polymerase chain reaction (qRT-PCR).

RESULTS

The qRT-PCR results revealed that the quantities of dsRNA in the apple leaves of treated trees were significantly greater than those in the leaves of untreated trees for both 2019 and 2020 studies. The peak dsRNA concentration in 2019 was 242 pg/30 mg of leaf tissue, and in 2020 was 16.4 pg/30 mg. The persistence of dsRNA in the apple tree canopy in 2019 was at least 84 days, and in 2020 was at least 141 days.

CONCLUSIONS

The highest mean measurement of dsRNA on a single date in 2019 was 242 pg, which is approximately equivalent to 8 ng/1 g leaf tissue. The projection using the highest replicate concentration from the same date is approximately equivalent to 27 ng/1 g leaf tissue, which may be sufficient to be considered biologically active. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Selective ablation of 3' RNA ends and processive RTs facilitate direct cDNA sequencing of full-length host cell and viral transcripts

Alternative splicing (AS) is necessary for viral proliferation in host cells and a critical regulatory component of viral gene expression. Conventional RNA-seq approaches provide incomplete coverage of AS due to their short read lengths and are susceptible to biases and artifacts introduced in prevailing library preparation methodologies. Moreover, viral splicing studies are often conducted separately from host cell transcriptome analysis, precluding an assessment of the viral manipulation of host splicing machinery. To address current limitations, we developed a quantitative full-length direct cDNA sequencing strategy to simultaneously profile viral and host cell transcripts. This nanopore-based approach couples processive reverse transcriptases with a novel one-step chemical ablation of 3' RNA ends (termed CASPR), which decreases ribosomal RNA reads and enriches polyadenylated coding sequences. We extensively validate our approach using synthetic reference transcripts and show that CASPR doubles the breadth of coverage per transcript and increases detection of long transcripts (>4 kb), while being functionally equivalent to PolyA+ selection for transcript quantification. We used our approach to interrogate host cell and HIV-1 transcript dynamics during viral reactivation and identified novel putative HIV-1 host factors containing exon skipping or novel intron retentions and delineated the HIV-1 transcriptional state associated with these differentially regulated host factors.

A highly sensitive strand-specific multiplex RT-qPCR assay for quantitation of Zika virus replication

Reverse-transcription quantitative polymerase chain reaction (RT-qPCR) is widely used to quantify viral RNA genomes for diagnostics and research, yet conventional RT-qPCR protocols are unable to accurately distinguish between the different viral RNA species that exist during infection. Here we show that false-priming and self-priming occur during reverse transcription with several published Zika virus (ZIKV) primer sets. We developed a RT-qPCR assay using tagged primers and thermostable reverse transcriptase, which greatly reduced the occurrence of nonspecific cDNA products. Furthermore, we optimized the assay for use in multiplex qPCR which allows for simultaneous quantitative detection of positive-strand and negative-strand ZIKV RNA along with an internal control from both human and mosquito cells. Importantly, this assay is sensitive enough to study early stages of virus infection in vitro. Strikingly, using this assay, we detected ZIKV negative-strand RNA as early as 3 h post-infection in mammalian cell culture, at a time point prior to the onset of positive-strand RNA synthesis. Overall, the strand-specific RT-qPCR assay developed herein is a valuable tool to quantify ZIKV RNA and to study viral replication dynamics during infection. The application of these findings has the potential to increase accuracy of RNA detection methods for a variety of viral pathogens.

HIV-1 Natural Antisense Transcription and Its Role in Viral Persistence

Natural antisense transcripts (NATs) represent a class of RNA molecules that are transcribed from the opposite strand of a protein-coding gene, and that have the ability to regulate the expression of their cognate protein-coding gene via multiple mechanisms. NATs have been described in many prokaryotic and eukaryotic systems, as well as in the viruses that infect them. The human immunodeficiency virus (HIV-1) is no exception, and produces one or more NAT from a promoter within the 3’ long terminal repeat. HIV-1 antisense transcripts have been the focus of several studies spanning over 30 years. However, a complete appreciation of the role that these transcripts play in the virus lifecycle is still lacking. In this review, we cover the current knowledge about HIV-1 NATs, discuss some of the questions that are still open and identify possible areas of future research.

Molecular Epidemiology of Citrus Leprosis Virus C: A New Viral Lineage and Phylodynamic of the Main Viral Subpopulations in the Americas

Despite the importance of viral strains/variants as agents of emerging diseases, genetic and evolutionary processes affecting their ecology are not fully understood. To get insight into this topic, we assessed the population and spatial dynamic parameters of citrus leprosis virus C (CiLV-C, genus Cilevirus, family Kitaviridae). CiLV-C is the etiological agent of citrus leprosis disease, a non-systemic infection considered the main viral disorder affecting citrus orchards in Brazil. Overall, we obtained 18 complete or near-complete viral genomes, 123 complete nucleotide sequences of the open reading frame (ORF) encoding the putative coat protein, and 204 partial nucleotide sequences of the ORF encoding the movement protein, from 430 infected Citrus spp. samples collected between 1932 and 2020. A thorough examination of the collected dataset suggested that the CiLV-C population consists of the major lineages CRD and SJP, unevenly distributed, plus a third one called ASU identified in this work, which is represented by a single isolate found in an herbarium sample collected in Asuncion, Paraguay, in 1937. Viruses from the three lineages share about 85% nucleotide sequence identity and show signs of inter-clade recombination events. Members of the lineage CRD were identified both in commercial and non-commercial citrus orchards. However, those of the lineages SJP were exclusively detected in samples collected in the citrus belt of São Paulo and Minas Gerais, the leading Brazilian citrus production region, after 2015. The most recent common ancestor of viruses of the three lineages dates back to, at least, ∼1500 years ago. Since citrus plants were introduced in the Americas by the Portuguese around the 1520s, the Bayesian phylodynamic analysis suggested that the ancestors of the main CiLV-C lineages likely originated in contact with native vegetation of South America. The intensive expansion of CRD and SJP lineages in Brazil started probably linked to the beginning of the local citrus industry. The high prevalence of CiLV-C in the citrus belt of Brazil likely ensues from the intensive connectivity between orchards, which represents a potential risk toward pathogen saturation across the region.

Co-application of biochar and nitrogen fertilizer reduced nitrogen losses from soil

Combined application of biochar and nitrogen (N) fertilizer has the potential to reduce N losses from soil. However, the effectiveness of biochar amendment on N management can vary with biochar types with different physical and chemical properties. This study aimed to assess the effect of two types of hardwood biochar with different ash contents and cation exchange capacity (CEC) on soil N mineralization and nitrous oxide (N2O) production when applied alone and in combination with N fertilizer. Soil samples collected from a temperate pasture system were amended with two types of biochar (B1 and B2), urea, and urea plus biochar, and incubated for 60 days along with soil control (without biochar or urea addition). Soil nitrate N, ammonium N, ammonia-oxidizing bacteria amoA gene transcripts, and N2O production were measured during the experiment. Compared to control, addition of B1 (higher CEC and lower ash content) alone decreased nitrate N concentration by 21% to 45% during the incubation period while the addition of B2 (lower CEC and higher ash content) alone increased the nitrate N concentration during the first 10 days. Biochar B1 also reduced the abundance of amoA transcripts by 71% after 60 days. Compared to B1 + urea, B2 + urea resulted in a significantly greater initial increase in soil ammonium and nitrate N concentrations. However, B2 + urea had a significantly lower 60-day cumulative N2O emission compared to B1 + urea. Overall, when applied with urea, the biochar with higher CEC reduced ammonification and nitrification rates, while biochar with higher ash content reduced N N2O production. Our study demonstrated that biochar has the potential to enhance N retention in soil and reduce N2O emission when it is applied with urea, but the specific effects of the added biochar depend on its physical and chemical properties.

Characterization of transcripts emanating from enhancer Eβ of the murine TCRβ locus

Enhancers are well established as critical regulators of gene expression, but the mechanisms underlying the molecular basis of their specificity and activity are only partly understood. One of the most exciting recent observations is the discovery of enhancer RNA (eRNA), a class of noncoding RNAs derived from enhancer regions. Transcription of developmentally regulated enhancers has been observed to be associated with their active state. The nature of transcripts (eRNA) and their functional attributes are diverse and context dependent. The majority of eRNA are nonpolyadenylated and present in low abundance owing to their low stability, and may represent transcriptional noise. However, some eRNAs have been reported to be reasonably long and stable, are enriched in nuclei, exhibit tissue-specific expression and may contribute to enhancer function. Transcription of enhancers has been postulated to mediate enhancer function through either the act of transcription or via the transcribed RNA per se and is a useful feature to be analysed to understand mechanisms underlying enhancer activity. Enhancer Eβ at the murine TCRβ locus has been reported to exhibit enhanced occupancy of RNA polymerase II in developing thymocytes. Here, we investigated the transcriptional potential of Eβ in developing thymocytes and detected overlapping bidirectional transcripts at Eβ ranging between 0.7 and 1.7 kb. These noncoding transcripts are capped, polyadenylated, nuclear and expressed specifically in thymocytes. Delineation of these characteristics is important to further investigate their functional roles in mediating enhancer activity.

Methodologies for Transcript Profiling Using Long-Read Technologies

RNA sequencing using next-generation sequencing technologies (NGS) is currently the standard approach for gene expression profiling, particularly for large-scale high-throughput studies. NGS technologies comprise high throughput, cost efficient short-read RNA-Seq, while emerging single molecule, long-read RNA-Seq technologies have enabled new approaches to study the transcriptome and its function. The emerging single molecule, long-read technologies are currently commercially available by Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT), while new methodologies based on short-read sequencing approaches are also being developed in order to provide long range single molecule level information-for example, the ones represented by the 10x Genomics linked read methodology. The shift toward long-read sequencing technologies for transcriptome characterization is based on current increases in throughput and decreases in cost, making these attractive for de novo transcriptome assembly, isoform expression quantification, and in-depth RNA species analysis. These types of analyses were challenging with standard short sequencing approaches, due to the complex nature of the transcriptome, which consists of variable lengths of transcripts and multiple alternatively spliced isoforms for most genes, as well as the high sequence similarity of highly abundant species of RNA, such as rRNAs. Here we aim to focus on single molecule level sequencing technologies and single-cell technologies that, combined with perturbation tools, allow the analysis of complete RNA species, whether short or long, at high resolution. In parallel, these tools have opened new ways in understanding gene functions at the tissue, network, and pathway levels, as well as their detailed functional characterization. Analysis of the epi-transcriptome, including RNA methylation and modification and the effects of such modifications on biological systems is now enabled through direct RNA sequencing instead of classical indirect approaches. However, many difficulties and challenges remain, such as methodologies to generate full-length RNA or cDNA libraries from all different species of RNAs, not only poly-A containing transcripts, and the identification of allele-specific transcripts due to current error rates of single molecule technologies, while the bioinformatics analysis on long-read data for accurate identification of 5' and 3' UTRs is still in development.

Strand-specific detection of overlapping transcripts via purification involving denaturation of biotinylated cDNA

Reverse transcription-PCR (RT-PCR) is the most widely employed technique for gene expression analysis owing to its high sensitivity, easy reproducibility and fast output. It has been conceived that priming RT reactions with gene-specific primers generates cDNA only from the specific RNA. However, several reports have revealed that cDNA is synthesized even without addition of exogenous primers in RT reactions. Owing to such self-priming activity, the signals from specific strands cannot be accurately detected and can confound the expression analysis, especially in context of overlapping bidirectional transcripts. Here, we demonstrate that purification of biotin-tagged cDNA in conjunction with alkaline denaturation can obviate the problem of background priming and enable accurate strand-specific detection of overlapping transcripts.

Coronavirus endoribonuclease targets viral polyuridine sequences to evade activating host sensors

Cells carry sensors that are primed to detect invading viruses. To avoid being recognized, coronaviruses express factors that interfere with host immune sensing pathways. Previous studies revealed that a coronavirus endoribonuclease (EndoU) delays activation of the host sensor system, but the mechanism was not known. Here, we report that EndoU cleaves a viral polyuridine sequence that would otherwise activate host immune sensors. This information may be used in developing inhibitors that target EndoU activity and prevent diseases caused by coronaviruses.

Coronaviruses (CoVs) are positive-sense RNA viruses that can emerge from endemic reservoirs and infect zoonotically, causing significant morbidity and mortality. CoVs encode an endoribonuclease designated EndoU that facilitates evasion of host pattern recognition receptor MDA5, but the target of EndoU activity was not known. Here, we report that EndoU cleaves the 5′-polyuridines from negative-sense viral RNA, termed PUN RNA, which is the product of polyA-templated RNA synthesis. Using a virus containing an EndoU catalytic-inactive mutation, we detected a higher abundance of PUN RNA in the cytoplasm compared to wild-type−infected cells. Furthermore, we found that transfecting PUN RNA into cells stimulates a robust, MDA5-dependent interferon response, and that removal of the polyuridine extension on the RNA dampens the response. Overall, the results of this study reveal the PUN RNA to be a CoV MDA5-dependent pathogen-associated molecular pattern (PAMP). We also establish a mechanism for EndoU activity to cleave and limit the accumulation of this PAMP. Since EndoU activity is highly conserved in all CoVs, inhibiting this activity may serve as an approach for therapeutic interventions against existing and emerging CoV infections.

Comparative Transcriptome Profiling of Disruptive Technology, Single- Molecule Direct RNA Sequencing

Background:

The standard approach for transcriptomic profiling involves high throughput short-read sequencing technology, mainly dominated by Illumina. However, the short reads have limitations in transcriptome assembly and in obtaining full-length transcripts due to the complex nature of transcriptomes with variable length and multiple alternative spliced isoforms. Recent advances in long read sequencing by the Oxford Nanopore Technologies (ONT) offered both cDNA as well as direct RNA sequencing and has brought a paradigm change in the sequencing technology to greatly improve the assembly and expression estimates. ONT enables molecules to be sequenced without fragmentation resulting in ultra-long read length enabling the entire genes and transcripts to be fully characterized. The direct RNA sequencing method, in addition, circumvents the reverse transcription and amplification steps.

Objective:

In this study, RNA sequencing methods were assessed by comparing data from Illumina (ILM), ONT cDNA (OCD) and ONT direct RNA (ODR).

Methods:

The sensitivity & specificity of the isoform detection was determined from the data generated by Illumina, ONT cDNA and ONT direct RNA sequencing technologies using Saccharomyces cerevisiae as model. Comparative studies were conducted with two pipelines to detect the isoforms, novel genes and variable gene length.

Results:

Mapping metrics and qualitative profiles for different pipelines are presented to understand these disruptive technologies. The variability in sequencing technology and the analysis pipeline were studied.

Regulation of myosin heavy chain antisense long noncoding RNA in human vastus lateralis in response to exercise training

Alterations to muscle activity or loading state can induce changes in expression of myosin heavy chain (MHC). For example, sedentary individuals that initiate exercise training can induce a pronounced shift from IIx to IIa MHC. We sought to examine the regulatory response of MHC RNA in human subjects in response to exercise training. In particular, we examined how natural antisense RNA transcripts (NATs) are regulated throughout the MHC gene locus that includes MYH2 (IIa), MYH1 (IIx), MYH4 (IIb), and MYH8 (Neonatal) in vastus lateralis before and after a 5-wk training regime that consisted of a combination of aerobic and resistance types of exercise. The exercise program induced a IIx to IIa MHC shift that was associated with a corresponding increase in transcription on the antisense strand of the IIx MHC gene and a decrease in antisense transcription of the IIa MHC gene, suggesting an inhibitory mechanism mediated by NATs. We also report that the absence of expression of IIb MHC in human limb muscle is associated with the abundant expression of antisense transcript overlapping the IIb MHC coding gene, which is the opposite expression pattern as compared with that previously observed in rats. The NAT provides a possible regulatory mechanism for the suppressed expression of IIb MHC in humans. These data indicate that NATs may play a regulatory role with regard to the coordinated shifts in MHC gene expression that occur in human muscle in response to exercise training.

Improved measurement of absolute mRNA quantity without reverse transcription

Gene expression studies using microarrays have provided important insights into understanding the mechanisms of transcriptional regulation in a variety of biological and disease phenomena. In a previous study, we developed Photo-DEAN, a universal-microarray-based RNA quantification method that enabled reverse transcription-free multiplex measurement of the absolute amount of RNA. Photo-DEAN promotes high-throughput and bias-less transcriptome analysis without the need for common controls or additional complicated normalization steps. In this study, we empirically identified two conditions (individual specificity and uniform duplex stability) necessary for in silico design of probe sequences, allowing the Photo-DEAN method to accurately measure the absolute amount of target RNA in total RNA. We then demonstrated that using the modified probe design conditions, the Photo-DEAN method successfully measured the absolute amount of pgi mRNA spiked into E. coli total RNA. The measurement was performed at five different sites in the coding region of pgi mRNA, exhibiting no significant site dependence. Theoretical considerations suggested that probe sequences longer than the previously used 30-bases better satisfy the necessary design conditions.

Comparative transcriptome analysis reveals genes related to the yolk ratio of duck eggs

Yolk ratio is an important production index in the salted duck egg industry. Yolk constituents are deposited during development of follicles. However, the molecular mechanism responsible for different yolk ratios in laying ducks remains elusive. In this study, Shaoxing ducks laying eggs with different yolk ratios were chosen for an analysis of liver and ovary transcriptome information. Twelve libraries were constructed and generated an average of 58.5 million clean reads per library, of which 69% of clean reads from liver and 65% of clean reads from ovary were mapped to a reference genome. Between cross-phenotype groups, a total of 250 and 230 differently expressed genes (DEGs) were identified in liver and ovary respectively, of which 101 and 50 DEGs respectively were characterized. Several DEGs were detected, among which HMGCS1, HMGCR, FDFT1, (DHCR7), (STARD4), CYP46A1 and LPIN3 are involved in cholesterol metabolism-related pathways; KIAA0319, STARD4, AP1S3, SH3GL2 and CAV2 are involved in vesicular transport in the liver; and ELOVL2 and PSD2 are involved in fatty acid elongation and endocytosis in the ovary. High yolk-ratio ducks had higher activity for cholesterol synthesis and molecular trafficking. The identification of candidate genes greatly advances the understanding of the genetic basis of the formation of different yolk ratios.

DNA methylation profiling allows for characterization of atrial and ventricular cardiac tissues and hiPSC-CMs

Background

Cardiac disease modelling using human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) requires thorough insight into cardiac cell type differentiation processes. However, current methods to discriminate different cardiac cell types are mostly time-consuming, are costly and often provide imprecise phenotypic evaluation. DNA methylation plays a critical role during early heart development and cardiac cellular specification. We therefore investigated the DNA methylation pattern in different cardiac tissues to identify CpG loci for further cardiac cell type characterization.

Results

An array-based genome-wide DNA methylation analysis using Illumina Infinium HumanMethylation450 BeadChips led to the identification of 168 differentially methylated CpG loci in atrial and ventricular human heart tissue samples (n = 49) from different patients with congenital heart defects (CHD). Systematic evaluation of atrial-ventricular DNA methylation pattern in cardiac tissues in an independent sample cohort of non-failing donor hearts and cardiac patients using bisulfite pyrosequencing helped us to define a subset of 16 differentially methylated CpG loci enabling precise characterization of human atrial and ventricular cardiac tissue samples. This defined set of reproducible cardiac tissue-specific DNA methylation sites allowed us to consistently detect the cellular identity of hiPSC-CM subtypes.

Conclusion

Testing DNA methylation of only a small set of defined CpG sites thus makes it possible to distinguish atrial and ventricular cardiac tissues and cardiac atrial and ventricular subtypes of hiPSC-CMs. This method represents a rapid and reliable system for phenotypic characterization of in vitro-generated cardiomyocytes and opens new opportunities for cardiovascular research and patient-specific therapy.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0679-0) contains supplementary material, which is available to authorized users.

Detection and mitigation of spurious antisense expression with RoSA

Antisense transcription is known to have a range of impacts on sense gene expression, including (but not limited to) impeding transcription initiation, disrupting post-transcriptional processes, and enhancing, slowing, or even preventing transcription of the sense gene. Strand-specific RNA-Seq protocols preserve the strand information of the original RNA in the data, and so can be used to identify where antisense transcription may be implicated in regulating gene expression. However, our analysis of 199 strand-specific RNA-Seq experiments reveals that spurious antisense reads are often present in these datasets at levels greater than 1% of sense gene expression levels. Furthermore, these levels can vary substantially even between replicates in the same experiment, potentially disrupting any downstream analysis, if the incorrectly assigned antisense counts dominate the set of genes with high antisense transcription levels. Currently, no tools exist to detect or correct for this spurious antisense signal. Our tool, RoSA (Removal of Spurious Antisense), detects the presence of high levels of spurious antisense read alignments in strand-specific RNA-Seq datasets. It uses incorrectly spliced reads on the antisense strand and/or ERCC spikeins (if present in the data) to calculate both global and gene-specific antisense correction factors. We demonstrate the utility of our tool to filter out spurious antisense transcript counts in an Arabidopsis thaliana RNA-Seq experiment. Availability: RoSA is open source software available under the GPL licence via the Barton Group GitHub page https://github.com/bartongroup.

Direct RNA Sequencing of the Coding Complete Influenza A Virus Genome

For the first time, a coding complete genome of an RNA virus has been sequenced in its original form. Previously, RNA was sequenced by the chemical degradation of radiolabeled RNA, a difficult method that produced only short sequences. Instead, RNA has usually been sequenced indirectly by copying it into cDNA, which is often amplified to dsDNA by PCR and subsequently analyzed using a variety of DNA sequencing methods. We designed an adapter to short highly conserved termini of the influenza A virus genome to target the (-) sense RNA into a protein nanopore on the Oxford Nanopore MinION sequencing platform. Utilizing this method with total RNA extracted from the allantoic fluid of influenza rA/Puerto Rico/8/1934 (H1N1) virus infected chicken eggs (EID₅₀ 6.8 × 10⁹), we demonstrate successful sequencing of the coding complete influenza A virus genome with 100% nucleotide coverage, 99% consensus identity, and 99% of reads mapped to influenza A virus. By utilizing the same methodology one can redesign the adapter in order to expand the targets to include viral mRNA and (+) sense cRNA, which are essential to the viral life cycle, or other pathogens. This approach also has the potential to identify and quantify splice variants and base modifications, which are not practically measurable with current methods.

Capturing the 'ome': the expanding molecular toolbox for RNA and DNA library construction

All sequencing experiments and most functional genomics screens rely on the generation of libraries to comprehensively capture pools of targeted sequences. In the past decade especially, driven by the progress in the field of massively parallel sequencing, numerous studies have comprehensively assessed the impact of particular manipulations on library complexity and quality, and characterized the activities and specificities of several key enzymes used in library construction. Fortunately, careful protocol design and reagent choice can substantially mitigate many of these biases, and enable reliable representation of sequences in libraries. This review aims to guide the reader through the vast expanse of literature on the subject to promote informed library generation, independent of the application.

A multiplex RNA quantification method to determine the absolute amounts of mRNA without reverse transcription

We have developed a highly sensitive microarray-based method that determines the absolute amounts of mRNA in a total RNA sample in a multiplex manner without reverse transcription. This direct mRNA measurement promotes high-throughput testing and reduces bias in transcriptome analyses. Furthermore, quantification of the absolute amount of mRNA allows transcriptome analysis without common controls or additional, complicated normalization. The method, called Photo-DEAN, was validated using chemically synthesized RNAs of known quantities and mouse liver total RNA samples. We found that the absolute amounts of mRNA were successfully measured without the cDNA synthesis step, with a sensitivity of 15 zmol achieved in 7 h.

Variable RNA expression from recently acquired, endogenous viral elements (EVE) of white spot syndrome virus (WSSV) in shrimp

The viral accommodation hypothesis proposes that endogenous viral elements (EVE) from both RNA and DNA viruses are being continually integrated into the shrimp genome by natural host processes and that they can result in tolerance to viral infection by fortuitous production of antisense, immunospecific RNA (imRNA). Thus, we hypothesized that previously reported microarray results for the presence of white spot syndrome virus (WSSV) open reading frames (ORFs) formerly called 151, 366 and 427 in a domesticated giant tiger shrimp (Penaeus monodon) breeding stock might have represented expression from EVE, since the stock had shown uninterrupted freedom from white spot disease (WSD) for many generations. To test this hypothesis, 128 specimens from a current stock generation were confirmed for freedom from WSSV infection using two nested PCR detection methods. Subsequent nested-PCR testing revealed 33/128 specimens (26%) positive for at least one of the ORF at very high sequence identity (95-99%) to extant WSSV. Positive results for ORF 366 (now known to be a fragment of the WSSV capsid protein gene) dominated (28/33 = 84.8%), so 9 arbitrarily selected 366-positive specimens were tested by strand-specific, nested RT-PCR using DNase-treated RNA templates. This revealed variable RNA expression in individual shrimp including no RNA transcripts (n = 1), sense transcripts only (n = 1), antisense transcripts only (n = 2) or transcripts of both sense (n = 5). The latter 7 expression products indicated specimens producing putative imRNA. The variable types and numbers of the EVE and the variable RNA expression (including potential imRNA) support predictions of the viral accommodation hypothesis that EVE are randomly produced and expressed. Positive nested PCR test results for EVE of ORF 366 using DNA templates derived from shrimp sperm (germ cells), indicated that they were heritable.

Absolute measurement of gene transcripts with Selfie-digital PCR

Absolute measurement of the number of RNA transcripts per gene is necessary to compare gene transcription among different tissues or experimental conditions and to assess transcription of genes that have a variable copy number per cell such as mitochondrial DNA. Here, we present a method called Selfie-digital PCR that measures the absolute amount of an RNA transcript produced by its own coding DNA at a particular moment. Overcoming the limitations of previous approaches, Selfie-digital PCR allows for the quantification of nuclear and mitochondrial gene transcription in a strand-specific manner that is comparable among tissues and cell types that differ in gene copy number or metabolic state. Using Selfie-digital PCR, we found that, with the exception of the liver, different organs exhibit marked variations in mitochondrial DNA copy number but similar transcription of mitochondrial DNA heavy and light chains, thus suggesting a preferential role of mitochondrial DNA abundance over its transcription in organ function. Moreover, the strand-specific analysis of mitochondrial transcription afforded by Selfie-digital PCR showed that transcription of the heavy strand was significantly higher than that of the light strand in all the tissues studied.

Multiplex, DNase-free one-step reverse transcription PCR for Plasmodium 18S rRNA and spliced gametocyte-specific mRNAs

Background

Plasmodium gametocytes are sexual stages transmitted to female Anopheles mosquitoes. While Plasmodium parasites can be differentiated microscopically on Giemsa-stained blood smears, molecular methods are increasingly used because of their increased sensitivity. Molecular detection of gametocytes requires methods that discriminate between asexual and sexual stage parasites. Commonly tested gametocyte-specific mRNAs are pfs25 and pfs230 detected by reverse transcription polymerase chain reaction (RT-PCR). However, detection of these unspliced mRNA targets requires preceding DNase treatment of nucleic acids to eliminate co-purified genomic DNA. If gametocyte-specific, spliced mRNAs could be identified, DNase treatment could be eliminated and one-step multiplexed molecular methods utilized.

Results

Expression data was used to identify highly-expressed mRNAs in mature gametocytes that were also low in antisense RNA expression in non-gametocyte stages. After testing numerous candidate mRNAs, the spliced female Pf3D7_0630000 mRNA was selected as a Plasmodium falciparum gametocyte-specific biomarker compatible with Plasmodium 18S rRNA RT-PCR. This mRNA was only detected in samples containing mature gametocytes and was absent in those containing only asexual stage parasites or uninfected human blood. PF3D7_0630000 RT-PCR detected gametocytes across a wide range of parasite densities in both spiked and clinical samples and agreed with pfs25 RT-PCR, the gold standard for RT-PCR-based gametocyte detection. PF3D7_0630000 multiplexed with Plasmodium 18S rRNA RT-PCR was more sensitive than other spliced mRNA targets for one-step RT-PCR gametocyte detection.

Conclusions

Because the spliced target does not require DNase treatment, the PF3D7_0630000 assay can be multiplexed with Plasmodium 18S rRNA for direct one-step detection of gametocytes from whole human blood.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-017-1863-3) contains supplementary material, which is available to authorized users.

Novel lincRNA Susceptibility Gene and Its Role in Etiopathogenesis of Thyrotoxic Periodic Paralysis

Thyrotoxic periodic paralysis (TPP) is a life-threatening neuromuscular complication of thyrotoxicosis characterized by muscle weakness and hypokalemia and with an unclear etiopathogenesis. However, the 17q24.3 locus had been genetically linked to TPP, in which the genetic variant rs312691 (TC genotype) in long intergenic noncoding RNA (lincRNA) CTD-2378E21.1 is located downstream of inward-rectifier potassium (Kir) channel genes [KCNJ2 and its antisense KCNJ2 (AS-KCNJ2)]. A TPP patient with a suppressed thyroid-stimulating hormone level, a high free thyroxine level of (5.8 ng/dL), and low serum potassium level of (2 mEq/L) was evaluated for Kir channel expression during and after recovery from thyrotoxicosis. We observed that circulating lincRNA and Kir expression varied in accordance with thyroid status and TC genotype. To endorse this association of a lincRNA-rs312691 variant with a genetic risk of TPP, an additional series of 37 patients with TPP and 32 patients with thyrotoxic without paralysis (TWP) were assessed. We verified that the risk of minor allele C was greater in TPP than in TWP (odds ratio, 5.289; P = 0.0062), and protective major allele T was more frequent than observed in the 1000 genome controls (odds ratio, 11.90; P < 0.0001). AS-KCNJ2 was downregulated during thyrotoxicosis in the TWP controls carrying allele T and were upregulated in those with TPP with risk allele C. Moreover, KCNJ2 (Kir2.1) expression was reduced during thyrotoxicosis and restored in euthyroid status. We further excluded any other coding variant by performing targeted exome sequencing mutational screening in 17q24.3. Our data suggest that high lincRNA AS-KCNJ2 and CDT-2378E21.1 expression, possibly driven by the triiodothyronine regulatory mechanism, reduces the Kir2.1 expression observed during thyrotoxicosis. This finding could contribute to the understanding of the reduced inward-rectifying current observed during muscle weakness in genetically susceptible TPP patients.

Expression of the Antisense-to-Latency Transcript Long Noncoding RNA in Kaposi's Sarcoma-Associated Herpesvirus

The regulation of latency is central to herpesvirus biology. Recent transcriptome-wide surveys have uncovered evidence for promiscuous transcription across the entirety of the Kaposi's sarcoma-associated herpesvirus (KSHV) genome and postulated the existence of multiple viral long noncoding RNAs (lncRNAs). Next-generation sequencing studies are highly dependent on the specific experimental approach and particular algorithms of analysis and therefore benefit from independent confirmation of the results. The antisense-to-latency transcript (ALT) lncRNA was discovered by genome-tiling microarray (Chandriani et al., J Virol 86:7934-7942, 2010, https://doi.org/10.1128/JVI.00645-10). To characterize ALT in detail, we physically isolated this lncRNA by a strand-specific hybrid capture assay and then employed transcriptome sequencing and novel reverse transcription-PCR (RT-PCR) assays to distinguish all RNA species in the KSHV latency region. These methods confirm that ALT initiates at positions 120739/121012 and encodes a single splice site, which is shared with the 3'-coterminal K14-vGPCR/ORF74 mRNA, terminating at 130873 (GenBank accession number GQ994935), resulting in an ∼10,000-nucleotide transcript. No shorter ALT isoforms were identified. This study also identified a novel intron within the LANA 5' untranslated region using a splice acceptor at 127888. In summary, ALT joins PAN/nut1/T1.1 as a bona fide lncRNA of KSHV with potentially important roles in viral gene regulation and pathogenesis.

IMPORTANCE

Increasing data support the importance of noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and lncRNAs, which have been shown to exert critical regulatory functions without coding for recognizable proteins. Defining the sequences of these ncRNAs is essential for future studies aiming to functionally characterize a specific ncRNA. Most lncRNA studies are highly dependent on high-throughput sequencing and bioinformatic analyses, few studies follow up on the initial predictions, and analyses are at times discordant. The manuscript characterizes one key viral lncRNA, ALT, by physically isolating ALT and by a sequencing-independent assay. It provides for a simple assay to monitor lncRNA expression in experimental and clinical samples. ALT is expressed antisense to the major viral latency transcripts encoding LANA as well as the viral miRNAs and thus has the potential to regulate this key part of the viral life cycle.

Assays for precise quantification of total (including short) and elongated HIV-1 transcripts

Despite intensive study, it is unclear which mechanisms are responsible for latent HIV infection in vivo. One potential mechanism is inhibition of HIV transcriptional elongation, which results in short abortive transcripts containing the trans-activation response (TAR) region. Because the relative levels of total (including short) and processive transcripts provide measures of HIV transcriptional initiation and elongation, there is a compelling need for techniques that accurately measure both. Nonetheless, prior assays for total transcripts have been semi-quantitative and have seen limited application to patient samples. This manuscript reports the validation of quantitative reverse transcription (RT) droplet digital PCR assays for measurement of total (TAR) and processive (R-U5/gag) HIV transcripts. Traditional RT priming strategies can efficiently detect the TAR region on long HIV transcripts but detect <4% of true short transcripts. The TAR assay presented here utilizes an initial polyadenylation step, which provides an accessible RT priming site and detects short and long transcripts with approximately equal efficiency (70%). By applying these assays to blood samples from 8 ART-treated HIV+ individuals, total HIV transcripts were detected at levels >10-fold higher than elongated transcripts, implying a substantial block to transcriptional elongation in vivo. This approach may be applied to other difficult-to-prime RNA targets.

Characterizing a novel and sensitive method to measure dsRNA in soil

Performing environmental assessments for double-stranded RNA-based agricultural products require the development of sensitive and selective methods to measure biodegradation rates of dsRNAs. We developed and characterized a novel analytical procedure that uses a molecular hybridization assay (QuantiGene(®)) to accurately measure dsRNA extracted from diverse soils. In this report, we utilize this method to demonstrate that two dsRNAs with distinct size, structure, and sequence degrade rapidly in soil with indistinguishable kinetics.

A Double-Hybridization Approach for the Transcription- and Amplification-Free Detection of Specific mRNA on a Microarray

A double-hybridization approach was developed for the enzyme-free detection of specific mRNA of a housekeeping gene. Targeted mRNA was immobilized by hybridization to complementary DNA capture probes spotted onto a microarray. A second hybridization step of Cy5-conjugated label DNA to another section of the mRNA enabled specific labeling of the target. Thus, enzymatic artifacts could be avoided by omitting transcription and amplification steps. This manuscript describes the development of capture probe molecules used in the transcription- and amplification-free analysis of RPLP0 mRNA in isolated total RNA. An increase in specific signal was found with increasing length of the target-specific section of capture probes. Unspecific signal comprising spot autofluorescence and unspecific label binding did not correlate with the capture length. An additional spacer between the specific part of the capture probe and the substrate attachment site increased the signal significantly only on a short capture probe of approximately 30 nt length.

H-NS suppresses pilE intragenic transcription and antigenic variation in Neisseria gonorrhoeae

Initially, pilE transcription in Neisseria gonorrhoeae appeared to be complicated, yet it was eventually simplified into a model where integration host factor activates a single -35/ -10 promoter. However, with the advent of high-throughput RNA sequencing, numerous small pil-specific RNAs (sense as well as antisense) have been identified at the pilE locus as well as at various pilS loci. Using a combination of in vitro transcription, site-directed mutagenesis, Northern analysis and quantitative reverse transcriptase PCR (qRT-PCR) analysis, we have identified three additional non-canonical promoter elements within the pilE gene; two are located within the midgene region (one sense and one antisense), with the third, an antisense promoter, located immediately downstream of the pilE ORF. Using strand-specific qRT-PCR analysis, an inverse correlation exists between the level of antisense expression and the amount of sense message. By their nature, promoter sequences tend to be AT-rich. In Escherichia coli, the small DNA-binding protein H-NS binds to AT-rich sequences and inhibits intragenic transcription. In N. gonorrhoeae hns mutants, pilE antisense transcription was increased twofold, with a concomitant decrease in sense transcript levels. However, most noticeably in these mutants, the absence of H-NS protein caused pilE/pilS recombination to increase dramatically when compared with WT values. Consequently, H-NS protein suppresses pilE intragenic transcription as well as antigenic variation through the pilE/pilS recombination system.

Exploring Differentially Expressed Genes and Natural Antisense Transcripts in Sheep (Ovis aries) Skin with Different Wool Fiber Diameters by Digital Gene Expression Profiling

Wool fiber diameter (WFD) is the most important economic trait of wool. However, the genes specifically controlling WFD remain elusive. In this study, the expression profiles of skin from two groups of Gansu Alpine merino sheep with different WFD (a super-fine wool group [FD = 18.0 ± 0.5 μm, n=3] and a fine wool group [FD=23.0 ± 0.5 μm, n=3]) were analyzed using next-generation sequencing-based digital gene expression profiling. A total of 40 significant differentially expressed genes (DEGs) were detected, including 9 up-regulated genes and 31 down-regulated genes. Further expression profile analysis of natural antisense transcripts (NATs) showed that more than 30% of the genes presented in sheep skin expression profiles had NATs. A total of 7 NATs with significant differential expression were detected, and all were down-regulated. Among of 40 DEGs, 3 DEGs (AQP8, Bos d2, and SPRR) had significant NATs which were all significantly down-regulated in the super-fine wool group. In total of DEGs and NATs were summarized as 3 main GO categories and 38 subcategories. Among the molecular functions, cellular components and biological processes categories, binding, cell part and metabolic process were the most dominant subcategories, respectively. However, no significant enrichment of GO terms was found (corrected P-value >0.05). The pathways that were significantly enriched with significant DEGs and NATs were mainly the lipoic acid metabolism, bile secretion, salivary secretion and ribosome and phenylalanine metabolism pathways (P < 0.05). The results indicated that expression of NATs and gene transcripts were correlated, suggesting a role in gene regulation. The discovery of these DEGs and NATs could facilitate enhanced selection for super-fine wool sheep through gene-assisted selection or targeted gene manipulation in the future.

Hypothesis: Artifacts, Including Spurious Chimeric RNAs with a Short Homologous Sequence, Caused by Consecutive Reverse Transcriptions and Endogenous Random Primers

Recent RNA-sequencing technology and associated bioinformatics have led to identification of tens of thousands of putative human chimeric RNAs, i.e. RNAs containing sequences from two different genes, most of which are derived from neighboring genes on the same chromosome. In this essay, we redefine "two neighboring genes" as those producing individual transcripts, and point out two known mechanisms for chimeric RNA formation, i.e. transcription from a fusion gene or trans-splicing of two RNAs. By our definition, most putative RNA chimeras derived from canonically-defined neighboring genes may either be technical artifacts or be cis-splicing products of 5'- or 3'-extended RNA of either partner that is redefined herein as an unannotated gene, whereas trans-splicing events are rare in human cells. Therefore, most authentic chimeric RNAs result from fusion genes, about 1,000 of which have been identified hitherto. We propose a hypothesis of "consecutive reverse transcriptions (RTs)", i.e. another RT reaction following the previous one, for how most spurious chimeric RNAs, especially those containing a short homologous sequence, may be generated during RT, especially in RNA-sequencing wherein RNAs are fragmented. We also point out that RNA samples contain numerous RNA and DNA shreds that can serve as endogenous random primers for RT and ensuing polymerase chain reactions (PCR), creating artifacts in RT-PCR.

Enhancing the speed of morpholino-DNA biosensor by electrokinetic concentration of DNA in a microfluidic chip

Electrokinetic methods that conveniently concentrate charged analytes by orders of magnitude are highly attractive for nucleic acid assays where they can bypass the complexity and costs of enzyme-based amplification. The present study demonstrates an electrokinetic concentration device incorporating charge-neutral morpholino (MO) probes: as DNA analyte is concentrated in a microfluidic channel using ion concentration polarization (ICP) it is simultaneously hybridized to spots of complementary MO probes immobilized on the channel floor. This approach is uniquely favored by the match between the optimum buffer ionic strength of approximately 10mM for both MO-DNA surface hybridization and electrokinetic concentration. The simple and easily scalable poly(dimethylsiloxane) (PDMS) microfluidic device was fabricated using soft lithography and contact printing of a conductive polymer, poly(3,4-ethylenedioxythiophene)-polystyrene sulfonate (

PEDOT

PSS) as a cation-selective membrane material. Using the microfluidic concentrator, we could increase the concentration of DNA by three orders of magnitude in less than 5 min at an electric field of 75 Vcm(-1). The 1000-fold increase in concentration of DNA led to an increase in the speed of MO-DNA hybridization by two orders of magnitude and enabled a detection sensitivity of ~1 nM within 15 min of concentration. Using the proposed microfluidic concentrator, we also demonstrated a rapid hybridization with a binary DNA mixture, containing a fully complementary and a non-complementary sequence to mimic molecular backgrounds present in real DNA samples.

Natural genetic variation impacts expression levels of coding, non-coding, and antisense transcripts in fission yeast

Our current understanding of how natural genetic variation affects gene expression beyond well-annotated coding genes is still limited. The use of deep sequencing technologies for the study of expression quantitative trait loci (eQTLs) has the potential to close this gap. Here, we generated the first recombinant strain library for fission yeast and conducted an RNA-seq-based QTL study of the coding, non-coding, and antisense transcriptomes. We show that the frequency of distal effects (trans-eQTLs) greatly exceeds the number of local effects (cis-eQTLs) and that non-coding RNAs are as likely to be affected by eQTLs as protein-coding RNAs. We identified a genetic variation of swc5 that modifies the levels of 871 RNAs, with effects on both sense and antisense transcription, and show that this effect most likely goes through a compromised deposition of the histone variant H2A.Z. The strains, methods, and datasets generated here provide a rich resource for future studies.

Identification, stability and expression of Sirt1 antisense long non-coding RNA

Natural antisense transcripts (NATs) exist ubiquitously as pivotal molecules to regulate coding gene expression. Sirtuin 1 (Sirt1) is a NAD-dependent deacetylase which is involved in myogenesis. However, whether Sirt1 transcribes NAT during C2C12 differentiation is still unknown. In this study, we identified a Sirt1 NAT which was designated as Sirt1 antisense long non-coding RNA (AS lncRNA) by sequencing and bioinformatic analysis. The level of Sirt1 AS lncRNA was greater in spleen but less in muscle tissue. The expression of both Sirt1 mRNA and Sirt1 AS lncRNA decreased during C2C12 myogenic differentiation, whereas the levels of miR-34a, which targets Sirt1, increased gradually. We further found that the half-life of Sirt1 AS lncRNA was 10h, but that of Sirt1 mRNA was 6h in C2C12 cells treated with 2 μg/ml Actinomycin D. Therefore, compared with Sirt1 mRNA, Sirt1 AS lncRNA was more stable. Overexpression of Sirt1 AS lncRNA increased the levels of Sirt1 protein, whereas overexpression of Sirt1 AS lncRNA mutant did not affect the level of Sirt1 protein in C2C12 cells. Moreover, downregulation of Sirt1 mRNA caused by miR-34a was counteracted by Sirt1 AS lncRNA in C2C12 cells. Taken together, we identified a novel NAT of Sirt1 which implicated in myogenesis through regulating Sirt1 expression.

Quantification of transgene-derived double-stranded RNA in plants using the QuantiGene nucleic acid detection platform

The expanding use of RNA interference (RNAi) in agricultural biotechnology necessitates tools for characterizing and quantifying double-stranded RNA (dsRNA)-containing transcripts that are expressed in transgenic plants. We sought to detect and quantify such transcripts in transgenic maize lines engineered to control western corn rootworm (Diabrotica virgifera virgifera LeConte) via overexpression of an inverted repeat sequence bearing a portion of the putative corn rootworm orthologue of yeast Snf7 (DvSnf7), an essential component of insect cell receptor sorting. A quantitative assay was developed to detect DvSnf7 sense strand-containing dsRNA transcripts that is based on the QuantiGene Plex 2.0 RNA assay platform from Affymetrix. The QuantiGene assay utilizes cooperative binding of multiple oligonucleotide probes with specificity for the target sequence resulting in exceptionally high assay specificity. Successful implementation of this assay required heat denaturation in the presence of the oligonucleotide probes prior to hybridization, presumably to dissociate primary transcripts carrying the duplex dsRNA structure. The dsRNA assay was validated using a strategy analogous to the rigorous enzyme-linked immunosorbent assay evaluations that are typically performed for foreign proteins expressed in transgenic plants. Validation studies indicated that the assay is sensitive (to 10 pg of dsRNA/g of fresh tissue), highly reproducible, and linear over ∼2.5 logs. The assay was validated using purified RNA from multiple maize tissue types, and studies indicate that the assay is also quantitative in crude tissue lysates. To the best of our knowledge, this is the first report of a non-polymerase chain reaction-based quantitative assay for dsRNA-containing transcripts, based on the use of the QuantiGene technology platform, and will broadly facilitate characterization of dsRNA in biological and environmental samples.

Interplay of chromatin modifications and non-coding RNAs in the heart

Precisely regulated patterns of gene expression are dependent on the binding of transcription factors and chromatin-associated determinants referred to as co-activators and co-repressors. These regulatory components function with the core transcriptional machinery to serve in critical activities to alter chromatin modification and regulate gene expression. While we are beginning to understand that cell-type specific patterns of gene expression are necessary to achieve selective cardiovascular developmental programs, we still do not know the molecular machineries that localize these determinants in the heart. With clear implications for the epigenetic control of gene expression signatures, the ENCODE (Encyclopedia of DNA Elements) Project Consortium determined that about 90% of the human genome is transcribed while only 1-2% of transcripts encode proteins. Emerging evidence suggests that non-coding RNA (ncRNA) serves as a signal for decoding chromatin modifications and provides a potential molecular basis for cell type-specific and promoter-specific patterns of gene expression. The discovery of the histone methyltransferase enzyme EZH2 in the regulation of gene expression patterns implicated in cardiac hypertrophy suggests a novel role for chromatin-associated ncRNAs and is the focus of this article.

Increase in cardiac myosin heavy-chain (MyHC) alpha protein isoform in hibernating ground squirrels, with echocardiographic visualization of ventricular wall hypertrophy and prolonged contraction

Deep hibernators such as golden-mantled ground squirrels (Callospermophilus lateralis) have multiple challenges to cardiac function during low temperature torpor and subsequent arousals. As heart rates fall from over 300 beats min(-1) to less than 10, chamber dilation and reduced cardiac output could lead to congestive myopathy. We performed echocardiography on a cohort of individuals prior to and after several months of hibernation. The left ventricular chamber exhibited eccentric and concentric hypertrophy during hibernation and thus calculated ventricular mass was ~30% greater. Ventricular ejection fraction was mildly reduced during hibernation but stroke volumes were greater due to the eccentric hypertrophy and dramatically increased diastolic filling volumes. Globally, the systolic phase in hibernation was ~9.5 times longer, and the diastolic phase was 28× longer. Left atrial ejection generally was not observed during hibernation. Atrial ejection returned weakly during early arousal. Strain echocardiography assessed the velocity and total movement distance of contraction and relaxation for regional ventricular segments in active and early arousal states. Myocardial systolic strain during early arousal was significantly greater than the active state, indicating greater total contractile movement. This mirrored the increased ventricular ejection fraction noted with early arousal. However, strain rates were slower during early arousal than during the active period, particularly systolic strain, which was 33% of active, compared with the rate of diastolic strain, which was 67% of active. As heart rate rose during the arousal period, myocardial velocities and strain rates also increased; this was matched closely by cardiac output. Curiously, though heart rates were only 26% of active heart rates during early arousal, the cardiac output was nearly 40% of the active state, suggesting an efficient pumping system. We further analyzed proportions of cardiac myosin heavy-chain (MyHC) isoforms in a separate cohort of squirrels over 5 months, including time points before hibernation, during hibernation and just prior to emergence. Hibernating individuals were maintained in both a 4°C cold room and a 20°C warm room. Measured by SDS-PAGE, relative percentages of cardiac MyHC alpha were increased during hibernation, at both hibernacula temperatures. A potential increase in contractile speed, and power, from more abundant MyHC alpha may aid force generation at low temperature and at low heart rates. Unlike many models of cardiomyopathies where the alpha isoform is replaced by the beta isoform in order to reduce oxygen consumption, ground squirrels demonstrate a potential cardioprotective mechanism to maintain cardiac output during torpor.

Development of a strand-specific real-time qRT-PCR for the accurate detection and quantitation of West Nile virus RNA

Studying the tropism and replication kinetics of West Nile virus (WNV) in different cell types in vitro and in tissues in animal models is important for understanding its pathogenesis. As detection of the negative strand viral RNA is a more reliable indicator of active replication for single-stranded positive-sense RNA viruses, the specificity of qRT-PCR assays currently used for the detection of WNV positive and negative strand RNA was reassessed. It was shown that self- and falsely-primed cDNA was generated during the reverse transcription step in an assay employing unmodified primers and several reverse transcriptases. As a result, a qRT-PCR assay using the thermostable rTth in combination with tagged primers was developed, which greatly improved strand specificity by circumventing the events of self- and false-priming. The reliability of the assay was then addressed in vitro using BV-2 microglia cells as well as in C57/BL6 mice. It was possible to follow the kinetics of positive and negative-strand RNA synthesis both in vitro and in vivo; however, the sensitivity of the assay will need to be optimized in order to detect and quantify negative-strand RNA synthesis in the very early stages of infection. Overall, the strand-specific qRT-PCR assay developed in this study is an effective tool to quantify WNV RNA, reassess viral replication, and study tropism of WNV in the context of WNV pathogenesis.

Directionality of noncoding human RNAs: how to avoid artifacts

Inactivation of tumor suppressor and metastasis suppressor genes via epigenetic silencing is a frequent event in human cancers. Recent work has shown new mechanisms of epigenetic silencing, based on the occurrence of long noncoding promoter-spanning antisense and/or sense RNAs (lncRNAs), which constitute part of chromatin silencing complexes. Using reverse transcription polymerase chain reaction (RT-PCR), we have started to scan "triple negative" and Her2-overexpressing breast cancer cell lines for directional/bidirectional transcription through promoters of tumor suppressor and metastasis suppressor genes known to be epigenetically silenced in vivo. Surprisingly, we found that RT-PCR-amplified products were obtained at high frequency in the absence of exogenous primers. These amplified products resulted from RT priming via transcripts originating from promoter or upstream spanning regions. Consequently, this priming overruled directionality determination and led to false detection-identification of such lncRNAs. We show that this prevalent "no primer" artifact can be eliminated by treating the RNA preparations with periodate, performing RT reactions at highly elevated temperatures, or a combination of both. These experimental improvements enabled determination of the presence and directionality of individual promoter-spanning long noncoding RNAs with certainty. Examples for the BRMS1 metastasis suppressor gene, as well as RAR-β2 and CST6 human tumor suppressor genes, in breast carcinoma cell lines are presented.

Development of a strand specific SYBRGreen RT-PCR for a GIII.2 bovine norovirus

A strand specific SYBRGreen RT-PCR was developed for a bovine norovirus (GIII.2). HEK293 cells were transfected with a plasmid containing the complete virus genome and copy DNA was produced with viral RNA strand-specific primers that introduced nucleotide changes. Amplicons from the negative and positive viral RNA strands, and from potential transcripts made by sequence independent transcription, were separated by melting curve analysis. The RT-PCR showed high strand specificity and could be a useful tool to study virus replication in replicon and reverse genetic systems and in screening for low levels of virus replication in norovirus permissive cell lines.

A genome-wide investigation of expression characteristics of natural antisense transcripts in liver and muscle samples of pigs

Natural antisense transcripts are endogenous transcripts that are complementary to the sense-strand of DNA. These transcripts have been identified in various eukaryotic species and are involved in a broad range of regulatory events and biological processes. However, their general biological functions, expression characteristics and regulatory mechanisms are still unclear. In this study, 497 liver and 586 muscle samples were harvested from a White Duroc×Erhualian F₂ resource population. The expression profiles of sense and antisense transcripts were determined by tag-based RNA sequencing. We identified 33.7% and 20.4% of transcripts having both sense and antisense expression, and 12.5% and 6.1% of transcripts only expressing antisense transcripts in liver and muscle, respectively. More than 32.2% of imprinting or predicted imprinting genes in the geneimprint database were detected with both sense and antisense expression. The correlations between sense and antisense expression in sense-antisense pairs were diverse in both liver and muscle, showing positive, negative or absent correlation. Antisense expression increases gene expression variability. More interestingly, compared to eQTL mapping of sense transcripts in which more than one eQTL was mapped for a transcript, only one eQTL was identified for each antisense transcript, and the percentage of cis-eQTL in antisense eQTL was higher than that in sense eQTL. This suggests that the expressions of antisense transcripts tend to be cis-regulated by a single genomic locus. To our knowledge, this study is the first systematical investigation of antisense transcription in pigs. The findings improve our understanding of the complexity of porcine transcriptome.

The search for new chlorophyll-binding proteins in the cyanobacterium Synechocystis sp. PCC 6803

Light harvesting provides a major challenge in the production of biofuels from microorganisms; while sunlight provides the energy necessary for biomass/biofuel production, at the same time it damages the cells. The genome of Synechocystis sp. PCC 6803 was searched for open reading frames that might code for yet unidentified chlorophyll-binding proteins with low molecular mass that could be involved in stress-adaptation. Amongst 9167 hypothetical ORFs corresponding to potential polypeptides of 100 amino acids or less, two were identified that had the potential to be pigment-binding, because they (i) encoded a potential transmembrane region, (ii) showed sequence similarity with known chlorophyll-binding domains, (iii) were conserved in other cyanobacterial species, and (iv) their codon adaptation index indicated significant translation probability. The two ORFs were located complementary (antisense) and internal to the ferrochelatase (hemH) and the pyruvate dehydrogenase (pdh) genes and therefore were named a-fch and a-pdh, respectively. Transcription of both genes was confirmed; however, no translated proteins could be detected immunologically. Whereas mutations within a-pdh or a-fch did not lead to any obvious phenotype, it is clear that transcripts and proteins over and above the currently known set may play a role in defining the physiology of cyanobacteria and other organisms.

BeeDoctor, a versatile MLPA-based diagnostic tool for screening bee viruses

The long-term decline of managed honeybee hives in the world has drawn significant attention to the scientific community and bee-keeping industry. A high pathogen load is believed to play a crucial role in this phenomenon, with the bee viruses being key players. Most of the currently characterized honeybee viruses (around twenty) are positive stranded RNA viruses. Techniques based on RNA signatures are widely used to determine the viral load in honeybee colonies. High throughput screening for viral loads necessitates the development of a multiplex polymerase chain reaction approach in which different viruses can be targeted simultaneously. A new multiparameter assay, called "BeeDoctor", was developed based on multiplex-ligation probe dependent amplification (MLPA) technology. This assay detects 10 honeybee viruses in one reaction. "BeeDoctor" is also able to screen selectively for either the positive strand of the targeted RNA bee viruses or the negative strand, which is indicative for active viral replication. Due to its sensitivity and specificity, the MLPA assay is a useful tool for rapid diagnosis, pathogen characterization, and epidemiology of viruses in honeybee populations. "BeeDoctor" was used for screening 363 samples from apiaries located throughout Flanders; the northern half of Belgium. Using the "BeeDoctor", virus infections were detected in almost eighty percent of the colonies, with deformed wing virus by far the most frequently detected virus and multiple virus infections were found in 26 percent of the colonies.