Leishmania sp identification by PCR associated with sequencing of target SSU rDNA in paraffin-embedded skin samples stored for more than 30 years

Paraffin-embedded samples commonly stored at educational and research institutions constitute tissues banks for follow-up or epidemiological studies; however, the paraffin inclusion process involves the use of substances that can cause DNA degradation. In this study, a PCR protocol was applied to identify Leishmania strains in 33 paraffin-embedded skin samples of patients with American cutaneous leishmaniasis. DNA was obtained by the phenol-chloroform protocol following paraffin removal and then used in PCR or nested PCR based on the nucleotide sequence of the small subunit ribosomal RNA (SSU rDNA). The amplicons obtained were cloned and sequenced to determine the single nucleotide polymorphism that distinguishes between different Leishmania species or groups. This assay allowed to distinguish organisms belonging to the subgenus Viannia and identify L. (Leishmania) amazonensis and L. (L.) chagasi of the Leishmania subgenus. Of the 33 samples, PCR and nested PCR identified 91% of samples. After sequencing the PCR product of 26 samples, 16 were identified as L. (L.) amazonensis, the other 10 contain organisms belonging to the L. (Viannia) sub-genus. These results open a huge opportunity to study stored samples and promote relevant contributions to epidemiological studies.

RNA trafficking and stabilization elements associate with multiple brain proteins

Two of the best understood somatic cell mRNA cytoplasmic trafficking elements are those governing localization of beta-actin and myelin basic protein mRNAs. These cis-acting elements bind the trans-acting factors fibroblast ZBP-1 and hnRNP A2, respectively. It is not known whether these elements fulfil other roles in mRNA metabolism. To address this question we have used Edman sequencing and western blotting to identify six rat brain proteins that bind the beta-actin element (zipcode). All are known RNA-binding proteins and differ from ZBP-1. Comparison with proteins that bind the hnRNP A2 and AU-rich response elements, A2RE/A2RE11 and AURE, showed that AURE and zipcode bind a similar set of proteins that does not overlap with those that bind A2RE11. The zipcode-binding protein, KSRP, and hnRNP A2 were selected for further study and were shown by confocal immunolluorescence microscopy to have similar distributions in the central nervous system, but they were found in largely separate locations in cell nuclei. In the cytoplasm of cultured oligodendrocytes they were segregated into separate populations of cytoplasmic granules. We conclude that not only may there be families of trans-acting factors for the same cis-acting element, which are presumably required at different stages of mRNA processing and metabolism, but independent factors may also target different and multiple RNAs in the same cell.

Heterogeneous nuclear ribonucleoprotein A3, a novel RNA trafficking response element-binding protein

The cis-acting response element, A2RE, which is sufficient for cytoplasmic mRNA trafficking in oligodendrocytes, binds a small group of rat brain proteins. Predominant among these is heterogeneous nuclear ribonucleoprotein (hnRNP) A2, a trans-acting factor for cytoplasmic trafficking of RNAs bearing A2RE-like sequences. We have now identified the other A2RE-binding proteins as hnRNP A1/A1(B), hnRNP B1, and four isoforms of hnRNP A3. The rat and human hnRNP A3 cDNAs have been sequenced, revealing the existence of alternatively spliced mRNAs. In Western blotting, 38-, 39-, 41 -, and 41.5-kDa components were all recognized by antibodies against a peptide in the glycine-rich region of hnRNP A3, but only the 41- and 41.5-kDa bands bound antibodies to a 15-residue N-terminal peptide encoded by an alternatively spliced part of exon 1. The identities of these four proteins were verified by Edman sequencing and mass spectral analysis of tryptic fragments generated from electrophoretically separated bands. Sequence-specific binding of bacterially expressed hnRNP A3 to A2RE has been demonstrated by biosensor and UV cross-linking electrophoretic mobility shift assays. Mutational analysis and confocal microscopy data support the hypothesis that the hnRNP A3 isoforms have a role in cytoplasmic trafficking of RNA.

Molecular and functional analyses of Kunjin virus infectious cDNA clones demonstrate the essential roles for NS2A in virus assembly and for a nonconservative residue in NS3 in RNA replication

A number of full-length cDNA clones of Kunjin virus (KUN) were previously prepared; it was shown that two of them, pAKUN and FLSDX, differed in specific infectivities of corresponding in vitro transcribed RNAs by similar to100,000-fold (A. A. Khromykh et al., J. Virol. 72:7270-7279, 1998). In this study, we analyzed a possible genetic determinant(s) of the observed differences in infectivity initially by sequencing the entire cDNAs of both clones and comparing them with the published sequence of the parental KUN strain MRM61C. We found six common amino acid residues in both cDNA clones that were different from those in the published MRM61C sequence but were similar to those in the published sequences of other flaviviruses from the same subgroup. pAKUN clone had four additional codon changes, i.e., Ile59 to Asn and Arg175 to Lys in NS2A and Tyr518 to His and Ser557 to Pro in NS3. Three of these substitutions except the previously shown marker mutation, Arg175 to Lys in NS2A, reverted to the wild-type sequence in the virus eventually recovered from pAKUN RNA-transfected BHK cells, demonstrating the functional importance of these residues in viral replication and/or viral assembly. Exchange of corresponding DNA fragments between pAKUN and FLSDX clones and site-directed mutagenesis revealed that the Tyr518-to-His mutation in NS3 was responsible for an similar to5-fold decrease in specific infectivity of transcribed RNA, while the Ile59-to-Asn mutation in NS2A completely blocked virus production. Correction of the Asn59 in pAKUN NS2A to the wild-type lie residue resulted in complete restoration of RNA infectivity. Replication of KUN replicon RNA with an Ile59-to-Asn substitution in NS2A and with a Ser557-to-Pro substitution in NS3 was not affected, while the Tyr518-to-His substitution in NS3 led to severe inhibition of RNA replication. The impaired function of the mutated NS2A in production of infectious virus was complemented in trans by the helper wild-type NS2A produced from the KUN replicon RNA. However, replicon RNA with mutated NS2A could not be packaged in trans by the KUN structural proteins. The data demonstrated essential roles for the KUN nonstructural protein NS2A in virus assembly and for NS3 in RNA replication and identified specific single-amino-acid residues involved in these functions.

High prevalence of GB Virus C/Hepatitis G Virus RNA among Brazilian blood donors

The aim of this study was to investigate the presence of the Hepatitis G Virus on a population of blood donors from São Paulo, Brazil and to evaluate its association to sociodemographic variables. Two RT-PCR systems targeting the putative 5'NCR and NS3 regions were employed and the former has shown a higher sensitivity. The observed prevalence of HGV-RNA on 545 blood donors was 9.7% (CI 95% 7.4;12.5). Statistical analysis depicted an association with race/ethnicity, black and mulatto donors being more frequently infected; and also with years of education, less educated donors presenting higher prevalences. No association was observed with other sociodemographic parameters as age, gender, place of birth and of residence. DNA sequencing of nine randomly chosen isolates demonstrated the presence of genotypes 1, 2 and 3 among our population but clustering of these Brazilian isolates was not detected upon phylogenetic analysis.

RNA-seq analysis of the Quercus suber root response to drought

Dissertação de mestrado em Plant Molecular Biology, Biotechnology and Bioentrepreneurship

An atlas of human gene expression from massively parallel signature sequencing (MPSS).

We have used massively parallel signature sequencing (MPSS) to sample the transcriptomes of 32 normal human tissues to an unprecedented depth, thus documenting the patterns of expression of almost 20,000 genes with high sensitivity and specificity. The data confirm the widely held belief that differences in gene expression between cell and tissue types are largely determined by transcripts derived from a limited number of tissue-specific genes, rather than by combinations of more promiscuously expressed genes. Expression of a little more than half of all known human genes seems to account for both the common requirements and the specific functions of the tissues sampled. A classification of tissues based on patterns of gene expression largely reproduces classifications based on anatomical and biochemical properties. The unbiased sampling of the human transcriptome achieved by MPSS supports the idea that most human genes have been mapped, if not functionally characterized. This data set should prove useful for the identification of tissue-specific genes, for the study of global changes induced by pathological conditions, and for the definition of a minimal set of genes necessary for basic cell maintenance. The data are available on the Web at http://mpss.licr.org and http://sgb.lynxgen.com.

Characteristic bimodal profiles of RNA polymerase II at thousands of active mammalian promoters.

BACKGROUND: In mammals, ChIP-seq studies of RNA polymerase II (PolII) occupancy have been performed to reveal how recruitment, initiation and pausing of PolII may control transcription rates, but the focus is rarely on obtaining finely resolved profiles that can portray the progression of PolII through sequential promoter states. RESULTS: Here, we analyze PolII binding profiles from high-coverage ChIP-seq on promoters of actively transcribed genes in mouse and humans. We show that the enrichment of PolII near transcription start sites exhibits a stereotypical bimodal structure, with one peak near active transcription start sites and a second peak 110 base pairs downstream from the first. Using an empirical model that reliably quantifies the spatial PolII signal, gene by gene, we show that the first PolII peak allows for refined positioning of transcription start sites, which is corroborated by mRNA sequencing. This bimodal signature is found both in mouse and humans. Analysis of the pausing-related factors NELF and DSIF suggests that the downstream peak reflects widespread pausing at the +1 nucleosome barrier. Several features of the bimodal pattern are correlated with sequence features such as CpG content and TATA boxes, as well as the histone mark H3K4me3. CONCLUSIONS: We thus show how high coverage DNA sequencing experiments can reveal as-yet unnoticed bimodal spatial features of PolII accumulation that are frequent at individual mammalian genes and reminiscent of transcription initiation and pausing. The initiation-pausing hypothesis is corroborated by evidence from run-on sequencing and immunoprecipitation in other cell types and species.

Establishment of HEV RNA reverse transcriptase PCR assays for the diagnosis of acute and chronic hepatitis E

Background and aim: H epatitis E v irus (HEV) infection has emerged as a c ause o f travel-related a nd autochthonous a cute hepatitis as well as chronic hepatitis in immunosuppressed patients. While t ravel-related cases a re c aused primarily b y infections w ith HEV of g enotype 1 ( HEV-1), autochthonous c ases a nd chronic cases a re d ue t o genotype 3 (HEV-3), which is s hared between humans and diverse animal species. The aim of this study was to establish HEV RNA detection assays f or q uantitative v iral load testing and genotyping. Methods: V iral RNA was p urified from plasma or s erum a nd converted to cDNA prior to (1) multiplex real-time PCR for HEV RNA quantification and (2) multiplex PCR coupled to DNA sequencing for HEV genotype determination. Real-time PCR was d esigned to match a ll known HEV genotypes available i n Genbank while PCR was designed using conserved primers flanking a variable region of the HEV RNA. Results: In a validation panel, the newly developed assays allowed for the reliable detection and genotyping of HEV-1 or HEV-3. Cases of t ravel-related and a utochthonous a cute h epatitis E a s well a s chronic hepatitis E i n immunosuppressed patients have b een identified using t hese a ssays a nd will be p resented in detail. Anti- HEV antibodies were n egative i n three well-characterized patients with chronic hepatitis E after organ transplantation. Conclusions: We developed and validated a quantitative HEV RNA detection assay that c an now be o ffered on a r outine basis (www.chuv.ch/imul/imu-collaborations-viral_hepatitis). Genotyping can also be offered on selected cases. HEV RNA detection is key in diagnosing chronic hepatitis E i n immunosuppressed patients with unexplained transaminase elevations, as serology can be negative in these patients.

Analysis of HIV-1 expression level and sense of transcription by high-throughput sequencing of the infected cell.

Next-generation sequencing offers an unprecedented opportunity to jointly analyze cellular and viral transcriptional activity without prerequisite knowledge of the nature of the transcripts. SupT1 cells were infected with a vesicular stomatitis virus G envelope protein (VSV-G)-pseudotyped HIV vector. At 24 h postinfection, both cellular and viral transcriptomes were analyzed by serial analysis of gene expression followed by high-throughput sequencing (SAGE-Seq). Read mapping resulted in 33 to 44 million tags aligning with the human transcriptome and 0.23 to 0.25 million tags aligning with the genome of the HIV-1 vector. Thus, at peak infection, 1 transcript in 143 is of viral origin (0.7%), including a small component of antisense viral transcription. Of the detected cellular transcripts, 826 (2.3%) were differentially expressed between mock- and HIV-infected samples. The approach also assessed whether HIV-1 infection modulates the expression of repetitive elements or endogenous retroviruses. We observed very active transcription of these elements, with 1 transcript in 237 being of such origin, corresponding on average to 123,123 reads in mock-infected samples (0.40%) and 129,149 reads in HIV-1-infected samples (0.45%) mapping to the genomic Repbase repository. This analysis highlights key details in the generation and interpretation of high-throughput data in the setting of HIV-1 cellular infection.

Added Value of Deep Sequencing Relative to Population Sequencing in Heavily Pre-Treated HIV-1-Infected Subjects.

OBJECTIVE: To explore the potential of deep HIV-1 sequencing for adding clinically relevant information relative to viral population sequencing in heavily pre-treated HIV-1-infected subjects. METHODS: In a proof-of-concept study, deep sequencing was compared to population sequencing in HIV-1-infected individuals with previous triple-class virological failure who also developed virologic failure to deep salvage therapy including, at least, darunavir, tipranavir, etravirine or raltegravir. Viral susceptibility was inferred before salvage therapy initiation and at virological failure using deep and population sequencing genotypes interpreted with the HIVdb, Rega and ANRS algorithms. The threshold level for mutant detection with deep sequencing was 1%. RESULTS: 7 subjects with previous exposure to a median of 15 antiretrovirals during a median of 13 years were included. Deep salvage therapy included darunavir, tipranavir, etravirine or raltegravir in 4, 2, 2 and 5 subjects, respectively. Self-reported treatment adherence was adequate in 4 and partial in 2; one individual underwent treatment interruption during follow-up. Deep sequencing detected all mutations found by population sequencing and identified additional resistance mutations in all but one individual, predominantly after virological failure to deep salvage therapy. Additional genotypic information led to consistent decreases in predicted susceptibility to etravirine, efavirenz, nucleoside reverse transcriptase inhibitors and indinavir in 2, 1, 2 and 1 subject, respectively. Deep sequencing data did not consistently modify the susceptibility predictions achieved with population sequencing for darunavir, tipranavir or raltegravir. CONCLUSIONS: In this subset of heavily pre-treated individuals, deep sequencing improved the assessment of genotypic resistance to etravirine, but did not consistently provide additional information on darunavir, tipranavir or raltegravir susceptibility. These data may inform the design of future studies addressing the clinical value of minority drug-resistant variants in treatment-experienced subjects.

Sequencing of E2 and NS5A regions of HCV genotype 3a in Brazilian patients with chronic hepatitis

Hepatitis C virus (HCV) is a major cause of liver disease throughout the world. The NS5A and E2 proteins of HCV genotype 1 were reported to inhibit the double-stranded (ds) RNA-dependent protein kinase (PKR), which is involved in the cellular antiviral response induced by interferon (IFN). The response to IFN therapy is quite different between genotypes, with response rates among patients infected with types 2 and 3 that are two-three-fold higher than in patients infected with type 1. Interestingly, a significant percentage of HCV genotype 3-infected patients do not respond to treatment at all. The aim of this paper was to analyse the sequences of fragments of the E2 and NS5A regions from 33 outpatients infected with genotype 3a, including patients that have responded (SVR) or not responded (NR) to treatment. HCV RNA was extracted and amplified with specific primers for the NS5A and E2 regions and the PCR products were then sequenced. The sequences obtained covered amino acids (aa) 636-708 in E2 and in NS5A [including the IFN sensitivity determining region (ISDR), PKR-binding domain and extended V3 region)]. In the E2 and NS5A regions, we did observe aa changes among patients, but these changes were not statistically significant between the SVR and NR groups. In conclusion, our results suggest that the ISDR domain is not predictive of treatment success in patients infected with HCV genotype 3a.

Direct detection of a BRAF mutation in total RNA from melanoma cells using cantilever arrays.

Malignant melanoma, the deadliest form of skin cancer, is characterized by a predominant mutation in the BRAF gene. Drugs that target tumours carrying this mutation have recently entered the clinic. Accordingly, patients are routinely screened for mutations in this gene to determine whether they can benefit from this type of treatment. The current gold standard for mutation screening uses real-time polymerase chain reaction and sequencing methods. Here we show that an assay based on microcantilever arrays can detect the mutation nanomechanically without amplification in total RNA samples isolated from melanoma cells. The assay is based on a BRAF-specific oligonucleotide probe. We detected mutant BRAF at a concentration of 500 pM in a 50-fold excess of the wild-type sequence. The method was able to distinguish melanoma cells carrying the mutation from wild-type cells using as little as 20 ng µl(-1) of RNA material, without prior PCR amplification and use of labels.

A comparison of methods for differential expression analysis of RNA-seq data.

BACKGROUND: Finding genes that are differentially expressed between conditions is an integral part of understanding the molecular basis of phenotypic variation. In the past decades, DNA microarrays have been used extensively to quantify the abundance of mRNA corresponding to different genes, and more recently high-throughput sequencing of cDNA (RNA-seq) has emerged as a powerful competitor. As the cost of sequencing decreases, it is conceivable that the use of RNA-seq for differential expression analysis will increase rapidly. To exploit the possibilities and address the challenges posed by this relatively new type of data, a number of software packages have been developed especially for differential expression analysis of RNA-seq data. RESULTS: We conducted an extensive comparison of eleven methods for differential expression analysis of RNA-seq data. All methods are freely available within the R framework and take as input a matrix of counts, i.e. the number of reads mapping to each genomic feature of interest in each of a number of samples. We evaluate the methods based on both simulated data and real RNA-seq data. CONCLUSIONS: Very small sample sizes, which are still common in RNA-seq experiments, impose problems for all evaluated methods and any results obtained under such conditions should be interpreted with caution. For larger sample sizes, the methods combining a variance-stabilizing transformation with the 'limma' method for differential expression analysis perform well under many different conditions, as does the nonparametric SAMseq method.

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver.

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes.