Combining RT-PCR-seq and RNA-seq to catalog all genic elements encoded in the human genome.

Within the ENCODE Consortium, GENCODE aimed to accurately annotate all protein-coding genes, pseudogenes, and noncoding transcribed loci in the human genome through manual curation and computational methods. Annotated transcript structures were assessed, and less well-supported loci were systematically, experimentally validated. Predicted exon-exon junctions were evaluated by RT-PCR amplification followed by highly multiplexed sequencing readout, a method we called RT-PCR-seq. Seventy-nine percent of all assessed junctions are confirmed by this evaluation procedure, demonstrating the high quality of the GENCODE gene set. RT-PCR-seq was also efficient to screen gene models predicted using the Human Body Map (HBM) RNA-seq data. We validated 73% of these predictions, thus confirming 1168 novel genes, mostly noncoding, which will further complement the GENCODE annotation. Our novel experimental validation pipeline is extremely sensitive, far more than unbiased transcriptome profiling through RNA sequencing, which is becoming the norm. For example, exon-exon junctions unique to GENCODE annotated transcripts are five times more likely to be corroborated with our targeted approach than with extensive large human transcriptome profiling. Data sets such as the HBM and ENCODE RNA-seq data fail sampling of low-expressed transcripts. Our RT-PCR-seq targeted approach also has the advantage of identifying novel exons of known genes, as we discovered unannotated exons in ~11% of assessed introns. We thus estimate that at least 18% of known loci have yet-unannotated exons. Our work demonstrates that the cataloging of all of the genic elements encoded in the human genome will necessitate a coordinated effort between unbiased and targeted approaches, like RNA-seq and RT-PCR-seq.

TLR3 and Rig-like receptor on myeloid dendritic cells and Rig-like receptor on human NK cells are both mandatory for production of IFN-gamma in response to double-stranded RNA.

Cross-talk between NK cells and dendritic cells (DCs) is critical for the potent therapeutic response to dsRNA, but the receptors involved remained controversial. We show in this paper that two dsRNAs, polyadenylic-polyuridylic acid and polyinosinic-polycytidylic acid [poly(I:C)], similarly engaged human TLR3, whereas only poly(I:C) triggered human RIG-I and MDA5. Both dsRNA enhanced NK cell activation within PBMCs but only poly(I:C) induced IFN-gamma. Although myeloid DCs (mDCs) were required for NK cell activation, induction of cytolytic potential and IFN-gamma production did not require contact with mDCs but was dependent on type I IFN and IL-12, respectively. Poly(I:C) but not polyadenylic-polyuridylic acid synergized with mDC-derived IL-12 for IFN-gamma production by acting directly on NK cells. Finally, the requirement of both TLR3 and Rig-like receptor (RLR) on mDCs and RLRs but not TLR3 on NK cells for IFN-gamma production was demonstrated using TLR3- and Cardif-deficient mice and human RIG-I-specific activator. Thus, we report the requirement of cotriggering TLR3 and RLR on mDCs and RLRs on NK cells for a pathogen product to induce potent innate cell activation.

Detection of Leishmania RNA virus in Leishmania parasites.

BACKGROUND: Patients suffering from cutaneous leishmaniasis (CL) caused by New World Leishmania (Viannia) species are at high risk of developing mucosal (ML) or disseminated cutaneous leishmaniasis (DCL). After the formation of a primary skin lesion at the site of the bite by a Leishmania-infected sand fly, the infection can disseminate to form secondary lesions. This metastatic phenotype causes significant morbidity and is often associated with a hyper-inflammatory immune response leading to the destruction of nasopharyngeal tissues in ML, and appearance of nodules or numerous ulcerated skin lesions in DCL. Recently, we connected this aggressive phenotype to the presence of Leishmania RNA virus (LRV) in strains of L. guyanensis, showing that LRV is responsible for elevated parasitaemia, destructive hyper-inflammation and an overall exacerbation of the disease. Further studies of this relationship and the distribution of LRVs in other Leishmania strains and species would benefit from improved methods of viral detection and quantitation, especially ones not dependent on prior knowledge of the viral sequence as LRVs show significant evolutionary divergence. METHODOLOGY/PRINCIPAL FINDINGS: This study reports various techniques, among which, the use of an anti-dsRNA monoclonal antibody (J2) stands out for its specific and quantitative recognition of dsRNA in a sequence-independent fashion. Applications of J2 include immunofluorescence, ELISA and dot blot: techniques complementing an arsenal of other detection tools, such as nucleic acid purification and quantitative real-time-PCR. We evaluate each method as well as demonstrate a successful LRV detection by the J2 antibody in several parasite strains, a freshly isolated patient sample and lesion biopsies of infected mice. CONCLUSIONS/SIGNIFICANCE: We propose that refinements of these methods could be transferred to the field for use as a diagnostic tool in detecting the presence of LRV, and potentially assessing the LRV-related risk of complications in cutaneous leishmaniasis.

Hepatitis C virus RNA replication requires a conserved structural motif within the transmembrane domain of the NS5B RNA-dependent RNA polymerase.

Hepatitis C virus (HCV) nonstructural protein 5B (NS5B), the viral RNA-dependent RNA polymerase (RdRp), is a tail-anchored protein with a highly conserved C-terminal transmembrane domain (TMD) that is required for the assembly of a functional replication complex. Here, we report that the TMD of the HCV RdRp can be functionally replaced by a newly identified analogous membrane anchor of the GB virus B (GBV-B) NS5B RdRp. Replicons with a chimeric RdRp consisting of the HCV catalytic domain and the GBV-B membrane anchor replicated with reduced efficiency. Compensatory amino acid changes at defined positions within the TMD improved the replication efficiency of these chimeras. These observations highlight a conserved structural motif within the TMD of the HCV NS5B RdRp that is required for RNA replication.

Aminoterminal amphipathic α-helix AH1 of hepatitis C virus nonstructural protein 4B possesses a dual role in RNA replication and virus production.

Nonstructural protein 4B (NS4B) is a key organizer of hepatitis C virus (HCV) replication complex formation. In concert with other nonstructural proteins, it induces a specific membrane rearrangement, designated as membranous web, which serves as a scaffold for the HCV replicase. The N-terminal part of NS4B comprises a predicted and a structurally resolved amphipathic α-helix, designated as AH1 and AH2, respectively. Here, we report a detailed structure-function analysis of NS4B AH1. Circular dichroism and nuclear magnetic resonance structural analyses revealed that AH1 folds into an amphipathic α-helix extending from NS4B amino acid 4 to 32, with positively charged residues flanking the helix. These residues are conserved among hepaciviruses. Mutagenesis and selection of pseudorevertants revealed an important role of these residues in RNA replication by affecting the biogenesis of double-membrane vesicles making up the membranous web. Moreover, alanine substitution of conserved acidic residues on the hydrophilic side of the helix reduced infectivity without significantly affecting RNA replication, indicating that AH1 is also involved in virus production. Selective membrane permeabilization and immunofluorescence microscopy analyses of a functional replicon harboring an epitope tag between NS4B AH1 and AH2 revealed a dual membrane topology of the N-terminal part of NS4B during HCV RNA replication. Luminal translocation was unaffected by the mutations introduced into AH1, but was abrogated by mutations introduced into AH2. In conclusion, our study reports the three-dimensional structure of AH1 from HCV NS4B, and highlights the importance of positively charged amino acid residues flanking this amphipathic α-helix in membranous web formation and RNA replication. In addition, we demonstrate that AH1 possesses a dual role in RNA replication and virus production, potentially governed by different topologies of the N-terminal part of NS4B.

The complementary strand of the human T-cell leukemia virus type 1 RNA genome encodes a bZIP transcription factor that down-regulates viral transcription.

The RNA genome of the human T-cell leukemia virus type 1 (HTLV-1) codes for proteins involved in infectivity, replication, and transformation. We report in this study the characterization of a novel viral protein encoded by the complementary strand of the HTLV-1 RNA genome. This protein, designated HBZ (for HTLV-1 bZIP factor), contains a N-terminal transcriptional activation domain and a leucine zipper motif in its C terminus. We show here that HBZ is able to interact with the bZIP transcription factor CREB-2 (also called ATF-4), known to activate the HTLV-1 transcription by recruiting the viral trans-activator Tax on the Tax-responsive elements (TxREs). However, we demonstrate that the HBZ/CREB-2 heterodimers are no more able to bind to the TxRE and cyclic AMP response element sites. Taking these findings together, the functional inactivation of CREB-2 by HBZ is suggested to contribute to regulation of the HTLV-1 transcription. Moreover, the characterization of a minus-strand gene protein encoded by HTLV-1 has never been reported until now.

Improved detection of Rhodococcus coprophilus with a new quantitative PCR assay.

Agricultural practices, such as spreading liquid manure or the utilisation of land as animal pastures, can result in faecal contamination of water resources. Rhodococcus coprophilus is used in microbial source tracking to indicate animal faecal contamination in water. Methods previously described for detecting of R. coprophilus in water were neither sensitive nor specific. Therefore, the aim of this study was to design and validate a new quantitative polymerase chain reaction (qPCR) to improve the detection of R. coprophilus in water. The new PCR assay was based on the R. coprophilus 16S rRNA gene. The validation showed that the new approach was specific and sensitive for deoxyribunucleic acid from target host species. Compared with other PCR assays tested in this study, the detection limit of the new qPCR was between 1 and 3 log lower. The method, including a filtration step, was further validated and successfully used in a field investigation in Switzerland. Our work demonstrated that the new detection method is sensitive and robust to detect R. coprophilus in surface and spring water. Compared with PCR assays that are available in the literature or to the culture-dependent method, the new molecular approach improves the detection of R. coprophilus.

Development of a duplex real-time PCR for the detection of Rickettsia spp. and typhus group rickettsia in clinical samples.

Molecular diagnosis using real-time polymerase chain reaction (PCR) may allow earlier diagnosis of rickettsiosis. We developed a duplex real-time PCR that amplifies (1) DNA of any rickettsial species and (2) DNA of both typhus group rickettsia, that is, Rickettsia prowazekii and Rickettsia typhi. Primers and probes were selected to amplify a segment of the 16S rRNA gene of Rickettsia spp. for the pan-rickettsial PCR and the citrate synthase gene (gltA) for the typhus group rickettsia PCR. Analytical sensitivity was 10 copies of control plasmid DNA per reaction. No cross-amplification was observed when testing human DNA and 22 pathogens or skin commensals. Real-time PCR was applied to 16 clinical samples. Rickettsial DNA was detected in the skin biopsies of three patients. In one patient with severe murine typhus, the typhus group PCR was positive in a skin biopsy from a petechial lesion and seroconversion was later documented. The two other patients with negative typhus group PCR suffered from Mediterranean and African spotted fever, respectively; in both cases, skin biopsy was performed on the eschar. Our duplex real-time PCR showed a good analytical sensitivity and specificity, allowing early diagnosis of rickettsiosis among three patients, and recognition of typhus in one of them.

The long non-coding RNA NEAT1 is increased in IUGR placentas, leading to potential new hypotheses of IUGR origin/development.

INTRODUCTION: Intrauterine Growth Restriction (IUGR) is a multifactorial disease defined by an inability of the fetus to reach its growth potential. IUGR not only increases the risk of neonatal mortality/morbidity, but also the risk of metabolic syndrome during adulthood. Certain placental proteins have been shown to be implicated in IUGR development, such as proteins from the GH/IGF axis and angiogenesis/apoptosis processes. METHODS: Twelve patients with term IUGR pregnancy (birth weight < 10th percentile) and 12 CTRLs were included. mRNA was extracted from the fetal part of the placenta and submitted to a subtraction method (Clontech PCR-Select cDNA Subtraction). RESULTS: One candidate gene identified was the long non-coding RNA NEAT1 (nuclear paraspeckle assembly transcript 1). NEAT1 is the core component of a subnuclear structure called paraspeckle. This structure is responsible for the retention of hyperedited mRNAs in the nucleus. Overall, NEAT1 mRNA expression was 4.14 (±1.16)-fold increased in IUGR vs. CTRL placentas (P = 0.009). NEAT1 was exclusively localized in the nuclei of the villous trophoblasts and was expressed in more nuclei and with greater intensity in IUGR placentas than in CTRLs. PSPC1, one of the three main proteins of the paraspeckle, co-localized with NEAT1 in the villous trophoblasts. The expression of NEAT1_2 mRNA, the long isoform of NEAT1, was only modestly increased in IUGR vs. CTRL placentas. DISCUSSION/CONCLUSION: The increase in NEAT1 and its co-localization with PSPC1 suggests an increase in paraspeckles in IUGR villous trophoblasts. This could lead to an increased retention of important mRNAs in villous trophoblasts nuclei. Given that the villous trophoblasts are crucial for the barrier function of the placenta, this could in part explain placental dysfunction in idiopathic IUGR fetuses.

Genome-wide RNA polymerase II profiles and RNA accumulation reveal kinetics of transcription and associated epigenetic changes during diurnal cycles.

Interactions of cell-autonomous circadian oscillators with diurnal cycles govern the temporal compartmentalization of cell physiology in mammals. To understand the transcriptional and epigenetic basis of diurnal rhythms in mouse liver genome-wide, we generated temporal DNA occupancy profiles by RNA polymerase II (Pol II) as well as profiles of the histone modifications H3K4me3 and H3K36me3. We used these data to quantify the relationships of phases and amplitudes between different marks. We found that rhythmic Pol II recruitment at promoters rather than rhythmic transition from paused to productive elongation underlies diurnal gene transcription, a conclusion further supported by modeling. Moreover, Pol II occupancy preceded mRNA accumulation by 3 hours, consistent with mRNA half-lives. Both methylation marks showed that the epigenetic landscape is highly dynamic and globally remodeled during the 24-hour cycle. While promoters of transcribed genes had tri-methylated H3K4 even at their trough activity times, tri-methylation levels reached their peak, on average, 1 hour after Pol II. Meanwhile, rhythms in tri-methylation of H3K36 lagged transcription by 3 hours. Finally, modeling profiles of Pol II occupancy and mRNA accumulation identified three classes of genes: one showing rhythmicity both in transcriptional and mRNA accumulation, a second class with rhythmic transcription but flat mRNA levels, and a third with constant transcription but rhythmic mRNAs. The latter class emphasizes widespread temporally gated posttranscriptional regulation in the mouse liver.

T-Coffee: a web server for the multiple sequence alignment of protein and RNA sequences using structural information and homology extension.

This article introduces a new interface for T-Coffee, a consistency-based multiple sequence alignment program. This interface provides an easy and intuitive access to the most popular functionality of the package. These include the default T-Coffee mode for protein and nucleic acid sequences, the M-Coffee mode that allows combining the output of any other aligners, and template-based modes of T-Coffee that deliver high accuracy alignments while using structural or homology derived templates. These three available template modes are Expresso for the alignment of protein with a known 3D-Structure, R-Coffee to align RNA sequences with conserved secondary structures and PSI-Coffee to accurately align distantly related sequences using homology extension. The new server benefits from recent improvements of the T-Coffee algorithm and can align up to 150 sequences as long as 10,000 residues and is available from both http://www.tcoffee.org and its main mirror http://tcoffee.crg.cat.

Akt signalling through GSK-3beta, mTOR and Foxo1 is involved in human skeletal muscle hypertrophy and atrophy

Skeletal muscle size is tightly regulated by the synergy between anabolic and catabolic signalling pathways which, in humans, have not been well characterized. Akt has been suggested to play a pivotal role in the regulation of skeletal muscle hypertrophy and atrophy in rodents and cells. Here we measured the amount of phospho-Akt and several of its downstream anabolic targets (glycogen synthase kinase-3beta (GSK-3beta), mTOR, p70(s6k) and 4E-BP1) and catabolic targets (Foxo1, Foxo3, atrogin-1 and MuRF1). All measurements were performed in human quadriceps muscle biopsies taken after 8 weeks of both hypertrophy-stimulating resistance training and atrophy-stimulating de-training. Following resistance training a muscle hypertrophy ( approximately 10%) and an increase in phospho-Akt, phospho-GSK-3beta and phospho-mTOR protein content were observed. This was paralleled by a decrease in Foxo1 nuclear protein content. Following the de-training period a muscle atrophy (5%), relative to the post-training muscle size, a decrease in phospho-Akt and GSK-3beta and an increase in Foxo1 were observed. Atrogin-1 and MuRF1 increased after the hypertrophy and decreased after the atrophy phases. We demonstrate, for the first time in human skeletal muscle, that the regulation of Akt and its downstream signalling pathways GSK-3beta, mTOR and Foxo1 are associated with both the skeletal muscle hypertrophy and atrophy processes

Toxicity of triclosan, penconazole and metalaxyl on Caulobacter crescentus and a freshwater microbial community as assessed by flow cytometry.

Biocides are widely used for domestic hygiene, agricultural and industrial applications. Their widespread use has resulted in their introduction into the environment and raised concerns about potential deleterious effects on aquatic ecosystems. In this study, the toxicity of the biocides triclosan, penconazole and metalaxyl were evaluated with the freshwater bacterium Caulobacter crescentus and with a freshwater microbial community using a combination of single- and double-stain flow cytometric assays. Growth of C.  crescentus and the freshwater community were repressed by triclosan but not by penconazole or metalaxyl at concentrations up to 250 μM. The repressive effect of triclosan was dependent on culture conditions. Caulobacter crescentus was more sensitive to triclosan when grown with high glucose at high cell density than when grown directly in sterilized lake water at low cell density. This suggests that the use of conventional growth conditions may overestimate biocide toxicity. Additional experiments showed that the freshwater community was more sensitive to triclosan than C.  crescentus, with 10 nM of triclosan being sufficient to repress growth and change the phylogenetic composition of the community. These results demonstrate that isolate-based assays may underestimate biocide toxicity and highlight the importance of assessing toxicity directly on natural microbial communities. Because 10 nM of triclosan is within the range of concentrations observed in freshwater systems, these results also raise concerns about the risk of introducing triclosan into the environment.

RDM1, a novel RNA recognition motif (RRM)-containing protein involved in the cell response to cisplatin in vertebrates.

A variety of cellular proteins has the ability to recognize DNA lesions induced by the anti-cancer drug cisplatin, with diverse consequences on their repair and on the therapeutic effectiveness of this drug. We report a novel gene involved in the cell response to cisplatin in vertebrates. The RDM1 gene (for RAD52 Motif 1) was identified while searching databases for sequences showing similarities to RAD52, a protein involved in homologous recombination and DNA double-strand break repair. Ablation of RDM1 in the chicken B cell line DT40 led to a more than 3-fold increase in sensitivity to cisplatin. However, RDM1-/- cells were not hypersensitive to DNA damages caused by ionizing radiation, UV irradiation, or the alkylating agent methylmethane sulfonate. The RDM1 protein displays a nucleic acid binding domain of the RNA recognition motif (RRM) type. By using gel-shift assays and electron microscopy, we show that purified, recombinant chicken RDM1 protein interacts with single-stranded DNA as well as double-stranded DNA, on which it assembles filament-like structures. Notably, RDM1 recognizes DNA distortions induced by cisplatin-DNA adducts in vitro. Finally, human RDM1 transcripts are abundant in the testis, suggesting a possible role during spermatogenesis.

Results of An International Study of Quantitative BCR/ABL Assays Using Defined RNA Samples

Molecular monitoring of BCR/ABL transcripts by real time quantitative reverse transcription PCR (qRT-PCR) is an essential technique for clinical management of patients with BCR/ABL-positive CML and ALL. Though quantitative BCR/ABL assays are performed in hundreds of laboratories worldwide, results among these laboratories cannot be reliably compared due to heterogeneity in test methods, data analysis, reporting, and lack of quantitative standards. Recent efforts towards standardization have been limited in scope. Aliquots of RNA were sent to clinical test centers worldwide in order to evaluate methods and reporting for e1a2, b2a2, and b3a2 transcript levels using their own qRT-PCR assays. Total RNA was isolated from tissue culture cells that expressed each of the different BCR/ABL transcripts. Serial log dilutions were prepared, ranging from 100 to 10-5, in RNA isolated from HL60 cells. Laboratories performed 5 independent qRT-PCR reactions for each sample type at each dilution. In addition, 15 qRT-PCR reactions of the 10-3 b3a2 RNA dilution were run to assess reproducibility within and between laboratories. Participants were asked to run the samples following their standard protocols and to report cycle threshold (Ct), quantitative values for BCR/ABL and housekeeping genes, and ratios of BCR/ABL to housekeeping genes for each sample RNA. Thirty-seven (n=37) participants have submitted qRT-PCR results for analysis (36, 37, and 34 labs generated data for b2a2, b3a2, and e1a2, respectively). The limit of detection for this study was defined as the lowest dilution that a Ct value could be detected for all 5 replicates. For b2a2, 15, 16, 4, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. For b3a2, 20, 13, and 4 labs showed a limit of detection at the 10-5, 10-4, and 10-3 dilutions, respectively. For e1a2, 10, 21, 2, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. Log %BCR/ABL ratio values provided a method for comparing results between the different laboratories for each BCR/ABL dilution series. Linear regression analysis revealed concordance among the majority of participant data over the 10-1 to 10-4 dilutions. The overall slope values showed comparable results among the majority of b2a2 (mean=0.939; median=0.9627; range (0.399 - 1.1872)), b3a2 (mean=0.925; median=0.922; range (0.625 - 1.140)), and e1a2 (mean=0.897; median=0.909; range (0.5174 - 1.138)) laboratory results (Fig. 1-3)). Thirty-four (n=34) out of the 37 laboratories reported Ct values for all 15 replicates and only those with a complete data set were included in the inter-lab calculations. Eleven laboratories either did not report their copy number data or used other reporting units such as nanograms or cell numbers; therefore, only 26 laboratories were included in the overall analysis of copy numbers. The median copy number was 348.4, with a range from 15.6 to 547,000 copies (approximately a 4.5 log difference); the median intra-lab %CV was 19.2% with a range from 4.2% to 82.6%. While our international performance evaluation using serially diluted RNA samples has reinforced the fact that heterogeneity exists among clinical laboratories, it has also demonstrated that performance within a laboratory is overall very consistent. Accordingly, the availability of defined BCR/ABL RNAs may facilitate the validation of all phases of quantitative BCR/ABL analysis and may be extremely useful as a tool for monitoring assay performance. Ongoing analyses of these materials, along with the development of additional control materials, may solidify consensus around their application in routine laboratory testing and possible integration in worldwide efforts to standardize quantitative BCR/ABL testing.