Production of full-length cDNA sequences by sequencing and analysis of expressed sequence tags from Schistosoma mansoni

The number of sequences generated by genome projects has increased exponentially, but gene characterization has not followed at the same rate. Sequencing and analysis of full-length cDNAs is an important step in gene characterization that has been used nowadays by several research groups. In this work, we have selected Schistosoma mansoni clones for full-length sequencing, using an algorithm that investigates the presence of the initial methionine in the parasite sequence based on the positions of alignment start between two sequences. BLAST searches to produce such alignments have been performed using parasite expressed sequence tags produced by Minas Gerais Genome Network against sequences from the database Eukaryotic Cluster of Orthologous Groups (KOG). This procedure has allowed the selection of clones representing 398 proteins which have not been deposited as S. mansoni complete CDS in any public database. Dedicated sequencing of 96 of such clones with reads from both 5' and 3' ends has been performed. These reads have been assembled using PHRAP, resulting in the production of 33 full-length sequences that represent novel S. mansoni proteins. These results shall contribute to construct a more complete view of the biology of this important parasite.

Identification of hepatitis C virus subtype 2c by sequencing analysis in patients from Córdoba, Argentina

In Argentina, most information on hepatitis C virus (HCV) genotype distribution comes from studies carried out in Buenos Aires (east province). In order to identify HCV subtypes in central Argentina, nucleotide sequencing of core region was performed in samples from 36 patients living in Córdoba, the second most populated province of Argentina. The sequence analysis identified subtype 2c as the most prevalent (50%), followed by subtype 1b (33%) and to a lesser extent by subtypes 1a (11%), 3a (3%) and 4a (3%). This is the first report of circulation of HCV subtype 2c in this region of Argentina and also such high prevalence has never been found before in the genotype distribution of South America.

Genome-wide re-sequencing of multidrug-resistant Mycobacterium leprae Airaku-3.

Genotyping and molecular characterization of drug resistance mechanisms in Mycobacterium leprae enables disease transmission and drug resistance trends to be monitored. In the present study, we performed genome-wide analysis of Airaku-3, a multidrug-resistant strain with an unknown mechanism of resistance to rifampicin. We identified 12 unique non-synonymous single-nucleotide polymorphisms (SNPs) including two in the transporter-encoding ctpC and ctpI genes. In addition, two SNPs were found that improve the resolution of SNP-based genotyping, particularly for Venezuelan and South East Asian strains of M. leprae.

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.

Comparison of HPV genotyping by type-specific PCR and sequencing

Human papillomavirus (HPV) infection is the most common sexually transmitted disease worldwide and there is a strong link between certain high-risk viral types and cervical carcinogenesis. Although there are several typing methods, it is still unclear which test is the best. This study compared the effectiveness of type-specific PCR (TS-PCR) and sequencing, with a focus on their clinical application. A total of 260 cervical samples from HPV-positive patients were tested for types 6, 11, 16, 18, 31, 33 and 35 using TS-PCR and sequencing. The genotype was identified in 36% of cases by TS-PCR and in 75% by sequencing. Sequencing was four times more likely to identify the viral type in positive samples than TS-PCR (p = 0.00). Despite being more effective for virus genotyping, sequencing was unable to identify viral types in multiple infections. Combining both techniques resulted in highly sensitive detection (87% of cases), showing that they are complementary methods. HPV genotyping is an important step in HPV management, helping to identify patients with a higher risk of developing cervical cancer and contributing to the development of type-specific vaccines.

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.

Detection of Leishmania infantum in naturally infected Lutzomyia longipalpis (Diptera: Psychodidae: Phlebotominae) and Canis familiaris in Misiones, Argentina: the first report of a PCR-RFLP and sequencing-based confirmation assay

In this study, a genotypification of Leishmaniawas performed using polimerase chain reaction-restriction fragment length polymorfism (PCR-RFLP) and sequencing techniques to identify species of Leishmaniaparasites in phlebotomine sand flies and dogs naturally infected. Between January-February of 2009, CDC light traps were used to collect insect samples from 13 capture sites in the municipality of Posadas, which is located in the province of Misiones of Argentina. Sand flies identified as Lutzomyia longipalpiswere grouped into 28 separate pools for molecular biological analysis. Canine samples were taken from lymph node aspirates of two symptomatic stray animals that had been positively diagnosed with canine visceral leishmaniasis. One vector pool of 10 sand flies (1 out of the 28 pools tested) and both of the canine samples tested positively for Leishmania infantumby PCR and RFLP analysis. PCR products were confirmed by sequencing and showed a maximum identity with L. infantum. Given that infection was detected in one out of the 28 pools and that at least one infected insect was infected, it was possible to infer an infection rate at least of 0.47% for Lu. longipalpisamong the analyzed samples. These results contribute to incriminate Lu. longipalpis as the vector of L. infantumin the municipality of Posadas, where cases of the disease in humans and dogs have been reported since 2005.

Molecular characterization of Echinococcus granulosusfrom Peru by sequencing of the mitochondrial cytochrome C oxidase subunit 1 gene

Echinococcus granulosus, the etiologic agent of cystic echinococcosis (CE) in humans and other animal species, is distributed worldwide. Ten intra-specific variants, or genotypes (G1-G10), have been defined based on genetic diversity. To determine the genotypes present in endemic areas of Peru, samples were collected from cattle (44), sheep (41) and humans (14) from Junín, Puno Huancavelica, Cusco, Arequipa and Ayacucho. DNA was extracted from protoscolex and/or germinal layers derived from 99 E. granulosus isolates and used as templates to amplify the mitochondrial cytochrome C oxidase subunit 1 gene. The resulting polymerase chain reaction products were sequenced and further examined by sequence analysis. All isolates, independent of the host, exhibited the G1 genotype. Phylogenetic analysis showed that three isolates from Ayacucho shared the same cluster with microvariant G1(4). The G1 genotype is considered the most widespread and infectious form of E. granulosusworldwide and our results confirm that the same patterns apply to this country. Therefore, these findings should be taken into consideration in developing prevention strategies and control programs for CE in Peru.

Identification of Mycobacterium tuberculosis complex based on amplification and sequencing of the oxyR pseudogene from stored Ziehl-Neelsen-stained sputum smears in Brazil

A cross-sectional analysis of stored Ziehl-Neelsen (ZN)-stained sputum smear slides (SSS) obtained from two public tuberculosis referral laboratories located in Juiz de Fora, Minas Gerais, was carried out to distinguish Mycobacterium bovis from other members of the Mycobacterium tuberculosis complex (MTC). A two-step approach was used to distinguish M. bovis from other members of MTC: (i) oxyR pseudogene amplification to detect MTC and, subsequently, (ii) allele-specific sequencing based on the polymorphism at position 285 of this gene. The oxyR pseudogene was successfully amplified in 100 of 177 (56.5%) SSS available from 99 individuals. No molecular profile of M. bovis was found. Multivariate analysis indicated that acid-fast bacilli (AFB) results and the source laboratory were associated (p < 0.05) with oxyR pseudogene amplification. SSS that were AFB++ SSS showed more oxyR pseudogene amplification than those with AFB0, possibly due to the amount of DNA. One of the two source laboratories presented a greater chance of oxyR pseudogene amplification, suggesting that differences in sputum conservation between laboratories could have influenced the preservation of DNA. This study provides evidence that stored ZN-SSS can be used for the molecular detection of MTC.

Genotyping of human parvovirus B19 in clinical samples from Brazil and Paraguay using heteroduplex mobility assay, single-stranded conformation polymorphism and nucleotide sequencing

Heteroduplex mobility assay, single-stranded conformation polymorphism and nucleotide sequencing were utilised to genotype human parvovirus B19 samples from Brazil and Paraguay. Ninety-seven serum samples were collected from individuals presenting with abortion or erythema infectiosum, arthropathies, severe anaemia and transient aplastic crisis; two additional skin samples were collected by biopsy. After the procedure, all clinical samples were classified as genotype 1.

Intelligent PCA contribution analysis for quality estimation in batch processes. Application in a sequencing batch reactor for wastewater treatment

En aquest treball, es proposa un nou mètode per estimar en temps real la qualitat del producte final en processos per lot. Aquest mètode permet reduir el temps necessari per obtenir els resultats de qualitat de les anàlisi de laboratori. S'utiliza un model de anàlisi de componentes principals (PCA) construït amb dades històriques en condicions normals de funcionament per discernir si un lot finalizat és normal o no. Es calcula una signatura de falla pels lots anormals i es passa a través d'un model de classificació per la seva estimació. L'estudi proposa un mètode per utilitzar la informació de les gràfiques de contribució basat en les signatures de falla, on els indicadors representen el comportament de les variables al llarg del procés en les diferentes etapes. Un conjunt de dades compost per la signatura de falla dels lots anormals històrics es construeix per cercar els patrons i entrenar els models de classifcació per estimar els resultas dels lots futurs. La metodologia proposada s'ha aplicat a un reactor seqüencial per lots (SBR). Diversos algoritmes de classificació es proven per demostrar les possibilitats de la metodologia proposada.

DNA multigene sequencing of topotypic specimens of the fascioliasis vector Lymnaea diaphana and phylogenetic analysis of the genus Pectinidens (Gastropoda)

Freshwater lymnaeid snails are crucial in defining transmission and epidemiology of fascioliasis. In South America, human endemic areas are related to high altitudes in Andean regions. The species Lymnaea diaphana has, however, been involved in low altitude areas of Chile, Argentina and Peru where human infection also occurs. Complete nuclear ribosomal DNA 18S, internal transcribed spacer (ITS)-2 and ITS-1 and fragments of mitochondrial DNA 16S and cytochrome c oxidase (cox)1 genes of L. diaphana specimens from its type locality offered 1,848, 495, 520, 424 and 672 bp long sequences. Comparisons with New and Old World Galba/Fossaria, Palaearctic stagnicolines, Nearctic stagnicolines, Old World Radix and Pseudosuccinea allowed to conclude that (i) L. diaphana shows sequences very different from all other lymnaeids, (ii) each marker allows its differentiation, except cox1 amino acid sequence, and (iii) L. diaphana is not a fossarine lymnaeid, but rather an archaic relict form derived from the oldest North American stagnicoline ancestors. Phylogeny and large genetic distances support the genus Pectinidens as the first stagnicoline representative in the southern hemisphere, including colonization of extreme world regions, as most southern Patagonia, long time ago. The phylogenetic link of L. diaphana with the stagnicoline group may give light to the aforementioned peculiar low altitude epidemiological scenario of fascioliasis.

Molecular diagnosis of eosinophilic meningitis due to Angiostrongylus cantonensis (Nematoda: Metastrongyloidea) by polymerase chain reaction-DNA sequencing of cerebrospinal fluids of patients

Cerebrospinal fluid (CSF) samples from clinically diagnosed patients with detectable Angiostrongylus canto-nensis-specific antibodies (n = 10), patients with clinically suspected cases that tested negative for A. cantonensis-an-tibodies (n = 5) and patients with cerebral gnathostomiasis (n = 2) and neurocysticercosis (n = 2) were examined by a single-step polymerase chain reaction (PCR) method using the AC primers for the 66-kDa native protein gene. The PCR method detected A. cantonensis DNA in CSF samples from four of 10 serologically confirmed angiostrongyliasis cases. The PCR results were negative for the remaining CSF samples. The nucleotide sequences of three positive CSF-PCR samples shared 98.8-99.2% similarity with the reference sequence of A. cantonensis. These results indicate the potential application of this PCR assay with clinical CSF samples for additional support in the confirmation of eosinophilic meningitis due to A. cantonensis.

Ultra high throughput sequencing in human DNA variation detection: a comparative study on the NDUFA3-PRPF31 region.

BACKGROUND: Ultra high throughput sequencing (UHTS) technologies find an important application in targeted resequencing of candidate genes or of genomic intervals from genetic association studies. Despite the extraordinary power of these new methods, they are still rarely used in routine analysis of human genomic variants, in part because of the absence of specific standard procedures. The aim of this work is to provide human molecular geneticists with a tool to evaluate the best UHTS methodology for efficiently detecting DNA changes, from common SNPs to rare mutations. METHODOLOGY/PRINCIPAL FINDINGS: We tested the three most widespread UHTS platforms (Roche/454 GS FLX Titanium, Illumina/Solexa Genome Analyzer II and Applied Biosystems/SOLiD System 3) on a well-studied region of the human genome containing many polymorphisms and a very rare heterozygous mutation located within an intronic repetitive DNA element. We identify the qualities and the limitations of each platform and describe some peculiarities of UHTS in resequencing projects. CONCLUSIONS/SIGNIFICANCE: When appropriate filtering and mapping procedures are applied UHTS technology can be safely and efficiently used as a tool for targeted human DNA variations detection. Unless particular and platform-dependent characteristics are needed for specific projects, the most relevant parameter to consider in mainstream human genome resequencing procedures is the cost per sequenced base-pair associated to each machine.