Rad53 arrests leading and lagging strand DNA synthesis via distinct mechanisms in response to DNA replication stress.

DNA replication stress threatens ordinary DNA synthesis. The evolutionarily conserved DNA replication stress response pathway involves sensor kinase Mec1/ATR, adaptor protein Mrc1/Claspin, and effector kinase Rad53/Chk1, which spurs a host of changes to stabilize replication forks and maintain genome integrity. DNA replication forks consist of largely distinct sets of proteins at leading and lagging strands that function autonomously in DNA synthesis in vitro. In this article, we discuss eSPAN and BrdU-IP-ssSeq, strand-specific sequencing technologies that permit analysis of protein localization and DNA synthesis at individual strands in budding yeast. Using these approaches, we show that under replication stress Rad53 stalls DNA synthesis on both leading and lagging strands. On lagging strands, it stimulates PCNA unloading, and on leading strands, it attenuates the replication function of Mrc1-Tof1. We propose that in doing so, Rad53 couples leading and lagging strand DNA synthesis during replication stress, thereby preventing the emergence of harmful ssDNA.

Epigenetic regulation of human-specific gene expression in the prefrontal cortex.

BACKGROUND: Changes in gene expression levels during brain development are thought to have played an important role in the evolution of human cognition. With the advent of high-throughput sequencing technologies, changes in brain developmental expression patterns, as well as human-specific brain gene expression, have been characterized. However, interpreting the origin of evolutionarily advanced cognition in human brains requires a deeper understanding of the regulation of gene expression, including the epigenomic context, along the primate genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) to measure the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 acetylation (H3K27ac), both of which are associated with transcriptional activation in the prefrontal cortex of humans, chimpanzees, and rhesus macaques. RESULTS: We found a discrete functional association, in which H3K4me3HP gain was significantly associated with myelination assembly and signaling transmission, while H3K4me3HP loss played a vital role in synaptic activity. Moreover, H3K27acHP gain was enriched in interneuron and oligodendrocyte markers, and H3K27acHP loss was enriched in CA1 pyramidal neuron markers. Using strand-specific RNA sequencing (ssRNA-seq), we first demonstrated that approximately 7 and 2% of human-specific expressed genes were epigenetically marked by H3K4me3HP and H3K27acHP, respectively, providing robust support for causal involvement of histones in gene expression. We also revealed the co-activation role of epigenetic modification and transcription factors in human-specific transcriptome evolution. Mechanistically, histone-modifying enzymes at least partially contribute to an epigenetic disturbance among primates, especially for the H3K27ac epigenomic marker. In line with this, peaks enriched in the macaque lineage were found to be driven by upregulated acetyl enzymes. CONCLUSIONS: Our results comprehensively elucidated a causal species-specific gene-histone-enzyme landscape in the prefrontal cortex and highlighted the regulatory interaction that drove transcriptional activation.

RNAseq analysis of mutants in coding and non-coding transcription of phosphate genes in the yeast Saccharomyces cerevisiae.

In the yeast Saccharomyces cerevisiae phosphate starvation induces the expression of PHO genes, including PHO84, encoding an high-affinity phosphate transporter, and SPL2, encoding a regulatory protein. PHO84 is down-regulated by antisense transcription. Here, using strand-specific RNAseq the effect is studied of mutations related to sense and antisense transcription of phosphate genes. Replacement of the transcriptional terminator of PHO84 by that of CYC1 resulted, unexpectedly, in an increased antisense transcription and a strongly reduced sense transcription of PHO84 and a strongly reduced SPL2 expression. The expression of unrelated genes was altered as well. The data suggest that antisense transcription of PHO84 and not the Pho84 transporter affects the expression of SPL2. Deletion of the two putative binding sites for Ume6 in the SPL2 promoter or deletion of UME6 differently affected SPL2 expression, suggesting that Ume6 regulates SPL2 by a mechanism different from a simple binding to the putative Ume6 binding sites.

Unraveling the lncRNA-miRNA-mRNA Regulatory Network Involved in Poplar Coma Development through High-Throughput Sequencing.

Poplar coma, the fluff-like appendages of seeds originating from the differentiated surface cells of the placenta and funicle, aids in the long-distance dispersal of seeds in the spring. However, it also poses hazards to human safety and causes pollution in the surrounding environment. Unraveling the regulatory mechanisms governing the initiation and development of coma is essential for addressing this issue comprehensively. In this study, strand-specific RNA-seq was conducted at three distinct stages of coma development, revealing 1888 lncRNAs and 52,810 mRNAs. The expression profiles of lncRNAs and mRNAs during coma development were analyzed. Subsequently, potential target genes of lncRNAs were predicted through co-localization and co-expression analyses. Integrating various types of sequencing data, lncRNA-miRNA-TF regulatory networks related to the initiation of coma were constructed. Utilizing identified differentially expressed genes encoding kinesin and actin, lncRNA-miRNA-mRNA regulatory networks associated with the construction and arrangement of the coma cytoskeleton were established. Additionally, relying on differentially expressed genes encoding cellulose synthase, sucrose synthase, and expansin, lncRNA-miRNA-mRNA regulatory networks related to coma cell wall synthesis and remodeling were developed. This study not only enhances the comprehension of lncRNA but also provides novel insights into the molecular mechanisms governing the initiation and development of poplar coma.

A highly discriminatory RNA strand-specific assay to facilitate analysis of the role of cis-acting elements in foot-and-mouth disease virus replication.

Foot-and-mouth-disease virus (FMDV), the aetiological agent responsible for foot-and-mouth disease (FMD), is a member of the genus Aphthovirus within the family Picornavirus. In common with all picornaviruses, replication of the single-stranded positive-sense RNA genome involves synthesis of a negative-sense complementary strand that serves as a template for the synthesis of multiple positive-sense progeny strands. We have previously employed FMDV replicons to examine viral RNA and protein elements essential to replication, but the factors affecting differential strand production remain unknown. Replicon-based systems require transfection of high levels of RNA, which can overload sensitive techniques such as quantitative PCR, preventing discrimination of specific strands. Here, we describe a method in which replicating RNA is labelled in vivo with 5-ethynyl uridine. The modified base is then linked to a biotin tag using click chemistry, facilitating purification of newly synthesised viral genomes or anti-genomes from input RNA. This selected RNA can then be amplified by strand-specific quantitative PCR, thus enabling investigation of the consequences of defined mutations on the relative synthesis of negative-sense intermediate and positive-strand progeny RNAs. We apply this new approach to investigate the consequence of mutation of viral cis-acting replication elements and provide direct evidence for their roles in negative-strand synthesis.

Analysis of long noncoding RNAs expression profiles in the human cardiac fibroblasts with cardiac fibrosis.

Cardiac fibrosis is a common pathological feature of cardiac remodelling process with disordered expression of multiple genes and eventually lead to heart failure. Emerging evidence suggests that long noncoding RNAs (lncRNAs) have emerged as critical regulators of various biological processes. However, the exact mechanisms of lncRNAs as mediators in cardiac fibrosis have not been fully elucidated. This study aimed to profile the lncRNA expression pattern in human cardiac fibroblasts (HCFs) with cardiac fibrosis. We treated HCFs with transforming growth factor-ß (TGF-ß) to induce their activation. Then, strand-specific RNA-seq was performed to profile and classify lncRNAs; and perform functional analysis in HCFs. We study the transformation of HCFs with molecular and cell biology methods. Among all identified lncRNA candidates, 176 and 526 lncRNAs were upregulated and downregulated respectively in TGF-ß-stimulated HCFs compared with controls. Functional analyses revealed that the target genes of differentially expressed lncRNAs were mainly related to focal adhesion, metabolic pathways, Hippo signaling pathway, PI3K-Akt signaling pathway, regulation of actin cytoskeleton, and hypertrophic cardiomyopathy. As a representative, novel lncRNAs NONHSAG005537 and NONHSAG017620 inhibited the proliferation, migration, invasion, and transformation of HCFs induced by TGF-ß. Collectively, our study established the expression signature of lncRNAs in cardiac fibrosis and demonstrated the cardioprotective role of NONHSAG005537 and NONHSAG017620 in cardiac fibrosis, providing a promising target for anti-fibrotic therapy.

Effect of temperature and IRF-9 gene-knockout on dynamics of vRNA, cRNA, and mRNA of viral hemorrhagic septicemia virus (VHSV).

The replication of viral hemorrhagic septicemia virus (VHSV) in appropriate host cells depends on environmental factors and the host cell's immunity. The dynamics of each VHSV RNA strand (vRNA, cRNA, and mRNA) in different conditions can provide a clue on the viral replication strategies, which can be a base for the development of efficient control measures. As VHSV is known to be sensitive to temperature and type I interferon (IFN) responses, in this study, we analyzed the effect of temperature difference (15 °C and 20 °C) and IRF-9 gene knockout on the dynamics of the three VHSV RNA strands in Epithelioma papulosum cyprini (EPC) cells using a strand-specific RT-qPCR. The tagged primers designed in this study successfully worked to quantify the three strands of VHSV. In the results of the temperature effect, the higher speed in viral mRNA transcription and the significantly higher (more than 10 times at 12-36 h) copy number of cRNA at 20 °C compared to those at 15 °C suggested the positive effect of high temperature on VHSV replication. In the results of the IRF-9 gene knockout effect, although IRF-9 gene knockout did not bring a dramatic effect on VHSV replication compared to the temperature effect, the increase of mRNA in IRF-9 KO cells was faster than normal EPC cells, which was reflected in the copy numbers of cRNA and vRNA. The IRF-9 gene knockout effect was not dramatic even in the replication of rVHSV-ΔNV-eGFP that harbors eGFP gene ORF instead of NV gene ORF. These results suggest that VHSV may be highly susceptible to pre-activated type I IFN responses but not highly susceptible to post-infection-mediated type I IFN responses or lowered type I IFN before infection. In both experiments of temperature effect and IRF-9 gene knockout effect, the copy number of cRNA never exceeded the copy number of vRNA at all assay times, suggesting that the binding efficiency of the RNP complex to the 3' end of cRNA might be lower than that to the 3' end of vRNA. Further research is needed to elucidate the regulatory mechanism that limits the amount of cRNA at an appropriate level during VHSV replication.

Revealing the novel complexity of plant long non-coding RNA by strand-specific and whole transcriptome sequencing for evolutionarily representative plant species.

BACKGROUND: Previous studies on plant long noncoding RNAs (lncRNAs) lacked consistency and suffered from many factors like heterogeneous data sources and experimental protocols, different plant tissues, inconsistent bioinformatics pipelines, etc. For example, the sequencing of RNAs with poly(A) tails excluded a large portion of lncRNAs without poly(A), and use of regular RNA-sequencing technique did not distinguish transcripts' direction for lncRNAs. The current study was designed to systematically discover and analyze lncRNAs across eight evolutionarily representative plant species, using strand-specific (directional) and whole transcriptome sequencing (RiboMinus) technique. RESULTS: A total of 39,945 lncRNAs (25,350 lincRNAs and 14,595 lncNATs) were identified, which showed molecular features of lncRNAs that are consistent across divergent plant species but different from those of mRNA. Further, transposable elements (TEs) were found to play key roles in the origination of lncRNA, as significantly large number of lncRNAs were found to contain TEs in gene body and promoter region, and transcription of many lncRNAs was driven by TE promoters. The lncRNA sequences were divergent even in closely related species, and most plant lncRNAs were genus/species-specific, amid rapid turnover in evolution. Evaluated with PhastCons scores, plant lncRNAs showed similar conservation level to that of intergenic sequences, suggesting that most lincRNAs were young and with short evolutionary age. INDUCED BY PHOSPHATE STARVATION (IPS) was found so far to be the only plant lncRNA group with conserved motifs, which may play important roles in the adaptation of terrestrial life during migration from aquatic to terrestrial. Most highly and specially expressed lncRNAs formed co-expression network with coding genes, and their functions were believed to be closely related to their co-expression genes. CONCLUSION: The study revealed novel features and complexity of lncRNAs in plants through systematic analysis, providing important insights into the origination and evolution of plant lncRNAs.

Single-strand specific nuclease enhances accuracy of error-corrected sequencing and improves rare mutation-detection sensitivity.

Error-corrected sequences (ECSs) that utilize double-stranded DNA sequences are useful in detecting mutagen-induced mutations. However, relatively higher frequencies of G:C > T:A (1 × 10-7 bp) and G:C > C:G (2 × 10-7 bp) errors decrease the accuracy of detection of rare G:C mutations (approximately 10-7 bp). Oxidized guanines in single-strand (SS) overhangs generated after shearing could serve as the source of these errors. To remove these errors, we first computationally discarded up to 20 read bases corresponding to the ends of the DNA fragments. Error frequencies decreased proportionately with trimming length; however, the results indicated that they were not sufficiently removed. To efficiently remove SS overhangs, we evaluated three mechanistically distinct SS-specific nucleases (S1 Nuclease, mung bean nuclease, and RecJf exonuclease) and found that they were more efficient than computational trimming. Consequently, we established Jade-Seq™, an ECS protocol with S1 Nuclease treatment, which reduced G:C > T:A and G:C > C:G errors to 0.50 × 10-7 bp and 0.12 × 10-7 bp, respectively. This was probably because S1 Nuclease removed SS regions, such as gaps and nicks, depending on its wide substrate specificity. Subsequently, we evaluated the mutation-detection sensitivity of Jade-Seq™ using DNA samples from TA100 cells exposed to 3-methylcholanthrene and 7,12-dimethylbenz[a]anthracene, which contained the rare G:C > T:A mutation (i.e., 2 × 10-7 bp). Fold changes of G:C > T:A compared to the vehicle control were 1.2- and 1.3-times higher than those of samples without S1 Nuclease treatment, respectively. These findings indicate the potential of Jade-Seq™ for detecting rare mutations and determining the mutagenicity of environmental mutagens.

In planta transcriptome analysis reveals tissue-specific expression of pathogenicity genes and microRNAs during rice-Magnaporthe interactions.

Transcriptional re-programming in host and pathogen upon leaf and neck infection is an evolving area of research for the rice blast community. Analysis of in planta rice transcriptome in leaf and neck tissues revealed tissue-specific and infection-specific expression of rice and Magnaporthe oryzae genes in host and pathogen. The glycosyl hydrolase, isocitrate lyase, cupin domain containing protein, TF2, CMPG1, CHIT17 and OsCML14 genes were uniquely expressed in leaf infection. Genes like cytochrome P450, inhibitor I family protein, GSTU6, abscisic stress ripening, and cupin domain containing protein were up-regulated during neck infection. In our microRNA sequencing study, Osa-miR166n-3p was highly expressed in upon Magnaporthe leaf infection, whereas osa-miR1661-3p, osa-miR166n-3p and osa-miR159b were overexpressed in neck infection. Here we report several transcripts being targeted by up and down regulated microRNAs during infection. The putative genes expressed upon infection in leaf and neck could be used in understanding the dual-epidemics of blast disease.

The Genome-Wide Identification of Long Non-Coding RNAs Involved in Floral Thermogenesis in Nelumbo nucifera Gaertn.

The sacred lotus (Nelumbo nucifera Gaertn.) can maintain a stable floral chamber temperature when blooming, despite ambient temperature fluctuations; however, the long non-coding RNAs (lncRNAs) involved in floral thermogenesis remain unclear. In the present study, we obtain comprehensive lncRNAs expression profiles from receptacles at five developmental stages by strand-specific RNA sequencing to reveal the lncRNAs regulatory mechanism of the floral thermogenesis of N. nucifera. A total of 22,693 transcripts were identified as lncRNAs, of which approximately 44.78% had stage-specific expression patterns. Subsequently, we identified 2579 differential expressed lncRNAs (DELs) regulating 2367 protein-coding genes mainly involved in receptacle development and reproductive process. Then, lncRNAs with floral thermogenesis identified by weighted gene co-expression network analysis (WGCNA) were mainly related to sulfur metabolism and mitochondrial electron transport chains. Meanwhile, 70 lncRNAs were predicted to act as endogenous target mimics (eTMs) for 29 miRNAs and participate in the regulation of 16 floral thermogenesis-related genes. Our dual luciferase reporter assays indicated that lncRNA LTCONS_00068702 acted as eTMs for miR164a_4 to regulate the expression of TrxL2 gene. These results deepen our understanding of the regulation mechanism of floral thermogenesis by lncRNAs and accumulate data for further research.

Identification of long noncoding natural antisense transcripts (lncNATs) correlated with drought stress response in wild rice (Oryza nivara).

BACKGROUND: Wild rice, including Oryza nivara and Oryza rufipogon, which are considered as the ancestors of Asian cultivated rice (Oryza sativa), possess high genetic diversity and serve as a crucial resource for breeding novel cultivars of cultivated rice. Although rice domestication related traits, such as seed shattering and plant architecture, have been intensively studied at the phenotypic and genomic levels, further investigation is needed to understand the molecular basis of phenotypic differences between cultivated and wild rice. Drought stress is one of the most severe abiotic stresses affecting rice growth and production. Adaptation to drought stress involves a cascade of genes and regulatory factors that form complex networks. O. nivara inhabits swampy areas with a seasonally dry climate, which is an ideal material to discover drought tolerance alleles. Long noncoding natural antisense transcripts (lncNATs), a class of long noncoding RNAs (lncRNAs), regulate the corresponding sense transcripts and play an important role in plant growth and development. However, the contribution of lncNATs to drought stress response in wild rice remains largely unknown. RESULTS: Here, we conducted strand-specific RNA sequencing (ssRNA-seq) analysis of Nipponbare (O. sativa) and two O. nivara accessions (BJ89 and BJ278) to determine the role of lncNATs in drought stress response in wild rice. A total of 1246 lncRNAs were identified, including 1091 coding-noncoding NAT pairs, of which 50 were expressed only in Nipponbare, and 77 were expressed only in BJ89 and/or BJ278. Of the 1091 coding-noncoding NAT pairs, 240 were differentially expressed between control and drought stress conditions. Among these 240 NAT pairs, 12 were detected only in Nipponbare, and 187 were detected uniquely in O. nivara. Furthermore, 10 of the 240 coding-noncoding NAT pairs were correlated with genes enriched in stress responsive GO terms; among these, nine pairs were uniquely found in O. nivara, and one pair was shared between O. nivara and Nipponbare. CONCLUSION: We identified lncNATs associated with drought stress response in cultivated rice and O. nivara. These results will improve our understanding of the function of lncNATs in drought tolerance and accelerate rice breeding.

Tight association of genome rearrangements with gene expression in conifer plastomes.

BACKGROUND: Our understanding of plastid transcriptomes is limited to a few model plants whose plastid genomes (plastomes) have a highly conserved gene order. Consequently, little is known about how gene expression changes in response to genomic rearrangements in plastids. This is particularly important in the highly rearranged conifer plastomes. RESULTS: We sequenced and reported the plastomes and plastid transcriptomes of six conifer species, representing all six extant families. Strand-specific RNAseq data show a nearly full transcription of both plastomic strands and detect C-to-U RNA-editing sites at both sense and antisense transcripts. We demonstrate that the expression of plastid coding genes is strongly functionally dependent among conifer species. However, the strength of this association declines as the number of plastomic rearrangements increases. This finding indicates that plastomic rearrangement influences gene expression. CONCLUSIONS: Our data provide the first line of evidence that plastomic rearrangements not only complicate the plastomic architecture but also drive the dynamics of plastid transcriptomes in conifers.

Effect of Y50H and S187G substitutions on thermostability and exonuclease activity of TK1646 from Thermococcus kodakarensis.

TK1646 is a highly thermostable single strand specific 3'-5' exonuclease. Exonucleases play important role in maintaining the genome integrity at elevated temperatures. Therefore, it is important to examine the factors contributing to thermostability of these exonucleases. In this study we report on production, purification and characterization of S187G and Y50H mutants of TK1646, focusing on the factors leading to thermostability of TK1646. Characterization of the recombinant proteins indicated that these substitutions did not drastically affect the catalysis of single stranded DNA. However, both of these substitutions reduced the thermostability of the recombinant proteins. Half-lives of Y50H and S187G mutants were 95 and 155 min, respectively, at 100 °C in comparison to 180 min of the wild type. Bioinformatics analysis indicated an increase in solvent accessibility of the mutated residues and disruption of hydrogens bonds. Molecular modelling and superimposition of the 3D structures of the mutants and the wild type demonstrated that one of the active site residues, Glu145, was shifted away from the metal ion in both the mutants which may be responsible for the decrease in catalytic activity. Compact secondary structure, hydrophobicity and hydrogen bonding might be the major factors contributing to the thermostability of TK1646.

Elimination of Viroids from Tobacco Pollen Involves a Decrease in Propagation Rate and an Increase of the Degradation Processes.

Some viroids-single-stranded, non-coding, circular RNA parasites of plants-are not transmissible through pollen to seeds and to next generation. We analyzed the cause for the elimination of apple fruit crinkle viroid (AFCVd) and citrus bark cracking viroid (CBCVd) from male gametophyte cells of Nicotiana tabacum by RNA deep sequencing and molecular methods using infected and transformed tobacco pollen tissues at different developmental stages. AFCVd was not transferable from pollen to seeds in reciprocal pollinations, due to a complete viroid eradication during the last steps of pollen development and fertilization. In pollen, the viroid replication pathway proceeds with detectable replication intermediates, but is dramatically depressed in comparison to leaves. Specific and unspecific viroid degradation with some preference for (-) chains occurred in pollen, as detected by analysis of viroid-derived small RNAs, by quantification of viroid levels and by detection of viroid degradation products forming "comets" on Northern blots. The decrease of viroid levels during pollen development correlated with mRNA accumulation of several RNA-degrading factors, such as AGO5 nuclease, DICER-like and TUDOR S-like nuclease. In addition, the functional status of pollen, as a tissue with high ribosome content, could play a role during suppression of AFCVd replication involving transcription factors IIIA and ribosomal protein L5.

Strand-Specific Dual RNA Sequencing of Bronchial Epithelial Cells Infected with Influenza A/H3N2 Viruses Reveals Splicing of Gene Segment 6 and Novel Host-Virus Interactions.

Host-influenza virus interplay at the transcript level has been extensively characterized in epithelial cells. Yet, there are no studies that simultaneously characterize human host and influenza A virus (IAV) genomes. We infected human bronchial epithelial BEAS-2B cells with two seasonal IAV/H3N2 strains, Brisbane/10/07 and Perth/16/09 (reference strains for past vaccine seasons) and the well-characterized laboratory strain Udorn/307/72. Strand-specific RNA sequencing (RNA-seq) of the infected BEAS-2B cells allowed for simultaneous analysis of host and viral transcriptomes, in addition to pathogen genomes, to reveal changes in mRNA expression and alternative splicing (AS). In general, patterns of global and immune gene expression induced by the three IAVs were mostly shared. However, AS of host transcripts and small nuclear RNAs differed between the seasonal and laboratory strains. Analysis of viral transcriptomes showed deletions of the polymerase components (defective interfering-like RNAs) within the genome. Surprisingly, we found that the neuraminidase gene undergoes AS and that the splicing event differs between seasonal and laboratory strains. Our findings reveal novel elements of the host-virus interaction and highlight the importance of RNA-seq in identifying molecular changes at the genome level that may contribute to shaping RNA-based innate immunity.IMPORTANCE The use of massively parallel RNA sequencing (RNA-seq) has revealed insights into human and pathogen genomes and their evolution. Dual RNA-seq allows simultaneous dissection of host and pathogen genomes and strand-specific RNA-seq provides information about the polarity of the RNA. This is important in the case of negative-strand RNA viruses like influenza virus, which generate positive (complementary and mRNA) and negative-strand RNAs (genome) that differ in their potential to trigger innate immunity. Here, we characterize interactions between human bronchial epithelial cells and three influenza A/H3N2 strains using strand-specific dual RNA-seq. We focused on this subtype because of its epidemiological importance in causing significant morbidity and mortality during influenza epidemics. We report novel elements that differ between seasonal and laboratory strains highlighting the complexity of the host-virus interplay at the RNA level.

Development of strand-specific real-time RT-PCR for the analysis of SCRV transcription and replication dynamics.

To distinguish between three types of Siniperca chuatsi rhabdovirus (SCRV) viral RNA (vRNA, cRNA, and mRNA) and investigate SCRV transcription and replication dynamics in Chinese perch brain CPB cells, a novel, strand-specific, reverse transcriptase quantitative real-time PCR (RT-qPCR) assay was established. The method is based on strand-specific reverse transcription, using tagged primers to add a 'tag' sequence at the 5' end. We used the 'tag' sequence as the forward primer and a strand-specific reverse primer to quantify the three types of RNA. Three types of synthetic viral RNA were used as reference standards for validation and quantification. These assays were optimized to produce a standard curve from 102 to 107 copies/µL, with an efficiency of 91-101% and an R2 value of 0.9949-0.9999. The coefficients of variation for repeatability and reproducibility were less than 2.85% and 5.52%, respectively. Using this method, specific target RNA was detected at a 3500-70,000 fold higher level than other types of RNA. This method was also used to evaluate the dynamics of vRNA, cRNA and mRNA synthesis in CPB cells infected with SCRV. The results indicate that the intracellular dynamics of vRNA, cRNA and mRNA are different. In the earliest phase of SCRV infection, all three types of viral RNA increased very slowly. The copy number of vRNA and mRNA increased exponentially from 4 h post infection, while cRNA increased from 6 h post infection. The amount of cRNA was lower than vRNA and mRNA throughout the infection. The novel, strand-specific RT-qPCR method developed in this study provides critical data to aid the understanding of transcription and replication during SCRV infection.

Strand-specific RNA-seq analysis of the Acidithiobacillus ferrooxidans transcriptome in response to magnesium stress.

Bioleaching is a promising process for 350 million tons Jinchuan low-grade pentlandite. But, Jinchuan pentlandite has lots of magnesium and high concentration of Mg2+ is harmful to bioleaching microorganisms. Thus, finding a way to improve the adaption of microorganisms to Mg2+ is a key for bioleaching. In the study, we found that oxidizing activity, bioleaching ability and biofilm formation of A.f were inhibited by Mg2+ stress. In addition, we analyzed mRNA and small RNA (sRNA) of Acidithiobacillus ferrooxidans (A.f) under Mg2+ stress by strand-specific RNA-sequencing (ssRNA-seq). After the bioinformatics process, 2475 coding genes were obtained, and there were 33 differential expression genes (DEGs) in 0.1 M-VS-Con, including 28 down-regulated and 5 up-regulated, whereas 52 DEGs were obtained in 0.5 M-VS-Con, including 28 down-regulated and 24 up-regulated. Gene ontology analysis showed most of DEGs were involved in catalytic activity, metabolic process and single-organism process. Furthermore, we identified 636 sRNA and some differential expression sRNA that may respond to Mg2+ stress. Further analysis of DEGs suggested that Mg2+ stress reduced biofilm formation perhaps through inhibiting Type IV Pili-related gene expression and inhibited bacterial activity perhaps through affecting carbon fixation. The study provided the foundation to understand the mechanisms of Mg2+ resistance in A.f and may be helpful to improve bioleaching ability for pentlandit.

Differentially expressed genes from RNA-Seq and functional enrichment results are affected by the choice of single-end versus paired-end reads and stranded versus non-stranded protocols.

BACKGROUND: RNA-Seq is now widely used as a research tool. Choices must be made whether to use paired-end (PE) or single-end (SE) sequencing, and whether to use strand-specific or non-specific (NS) library preparation kits. To date there has been no analysis of the effect of these choices on identifying differentially expressed genes (DEGs) between controls and treated samples and on downstream functional analysis. RESULTS: We undertook four mammalian transcriptomics experiments to compare the effect of SE and PE protocols on read mapping, feature counting, identification of DEGs and functional analysis. For three of these experiments we also compared a non-stranded (NS) and a strand-specific approach to mapping the paired-end data. SE mapping resulted in a reduced number of reads mapped to features, in all four experiments, and lower read count per gene. Up to 4.3% of genes in the SE data and up to 12.3% of genes in the NS data had read counts which were significantly different compared to the PE data. Comparison of DEGs showed the presence of false positives (average 5%, using voom) and false negatives (average 5%, using voom) using the SE reads. These increased further, by one or two percentage points, with the NS data. Gene ontology functional enrichment (GO) of the DEGs arising from SE or NS approaches, revealed striking differences in the top 20 GO terms, with as little as 40% concordance with PE results. Caution is therefore advised in the interpretation of such results. By comparison, there was overall consistency in gene set enrichment analysis results. CONCLUSIONS: A strand-specific protocol should be used in library preparation to generate the most reliable and accurate profile of expression. Ideally PE reads are also recommended particularly for transcriptome assembly. Whilst SE reads produce a DEG list with around 5% of false positives and false negatives, this method can substantially reduce sequencing cost and this saving could be used to increase the number of biological replicates thereby increasing the power of the experiment. As SE reads, when used in association with gene set enrichment, can generate accurate biological results, this may be a desirable trade-off.

Genome-wide analysis of rice cis-natural antisense transcription under cadmium exposure using strand-specific RNA-Seq.

BACKGROUND: The elucidation of novel transcripts and their expression in response to various stress conditions is necessary to understand the transcriptional network of plants as an adaptation to biotic and abiotic stresses. We performed strand-specific RNA-Seq (ssRNA-Seq) on rice exposed to cadmium (Cd) for 24 h and investigated the expression of cis-natural antisense transcripts (cis-NATs), a class of endogenous coding or non-protein-coding RNAs with sequence complementarity to the opposite strands of RAP transcripts. RESULTS: Many RAP transcripts possessed cis-NATs and these cis-NATs were responsive to some extent. Cis-NATs were upregulated from 26, 266 and 409 RAP gene loci, while 2054, 2501 and 2825 RAP transcripts were upregulated from 38,123 RAP loci under high Cd exposure in roots at 1, 12 and 24 h, respectively. In addition, most of the upregulated cis-NATs showed little upregulation under ABA or cold treatment. A number of cis-NATs were upregulated from less than 35 RAP gene loci in different tissue and time-point combinations under low Cd exposure, suggesting that cis-NATs respond to environmental stress. Furthermore, 409 RAP transcripts with upregulated cis-NATs were classified into three groups based on the expression of the RAP transcripts from the opposite DNA strand, including 138 upregulated, 128 invariable, and 143 downregulated transcripts, although the responses of cis-NATs and RAP transcripts were not always correlated. CONCLUSIONS: We have shown that the cis-NATs identified by ssRNA-Seq analysis are novel genes and that some of them are stress-specific and show different responses depending on the degree of stress and tissue. These results improve our understanding of the complete molecular mechanism of plant adaptation to Cd exposure.