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.

Identification and characterisation of microRNAs in young adults of Angiostrongylus cantonensis via a deep-sequencing approach

Angiostrongylus cantonensis is an important causative agent of eosinophilic meningitis and eosinophilic meningoencephalitis in humans. MicroRNAs (miRNAs) are small non-coding RNAs that participate in a wide range of biological processes. This study employed a deep-sequencing approach to study miRNAs from young adults of A. cantonensis. Based on 16,880,456 high-quality reads, 252 conserved mature miRNAs including 10 antisense miRNAs that belonging to 90 families, together with 10 antisense miRNAs were identified and characterised. Among these sequences, 53 miRNAs from 25 families displayed 50 or more reads. The conserved miRNA families were divided into four groups according to their phylogenetic distribution and a total of nine families without any members showing homology to other nematodes or adult worms were identified. Stem-loop real-time polymerase chain reaction analysis of aca-miR-1-1 and aca-miR-71-1 demonstrated that their level of expression increased dramatically from infective larvae to young adults and then decreased in adult worms, with the male worms exhibiting significantly higher levels of expression than female worms. These findings provide information related to the regulation of gene expression during the growth, development and pathogenesis of young adults of A. cantonensis.

Genetics and Prognostication in Splenic Marginal Zone Lymphoma: Revelations from Deep Sequencing

Purpose: Mounting evidence supports the clinical significance of gene mutations and immunogenetic features in common mature B-cell malignancies.

Experimental Design: We undertook a detailed characterization of the genetic background of splenic marginal zone lymphoma (SMZL), using targeted resequencing and explored potential clinical implications in a multinational cohort of 175 patients with SMZL.

Results: We identified recurrent mutations in TP53 (16%), KLF2 (12%), NOTCH2 (10%), TNFAIP3 (7%), MLL2 (11%), MYD88 (7%), and ARID1A (6%), all genes known to be targeted by somatic mutation in SMZL. KLF2 mutations were early, clonal events, enriched in patients with del(7q) and IGHV1-2*04 B-cell receptor immunoglobulins, and were associated with a short median time to first treatment (0.12 vs. 1.11 years; P = 0.01). In multivariate analysis, mutations in NOTCH2 [HR, 2.12; 95% confidence interval (CI), 1.02–4.4; P = 0.044] and 100% germline IGHV gene identity (HR, 2.19; 95% CI, 1.05–4.55; P = 0.036) were independent markers of short time to first treatment, whereas TP53 mutations were an independent marker of short overall survival (HR, 2.36; 95 % CI, 1.08–5.2; P = 0.03).

Conclusions: We identify key associations between gene mutations and clinical outcome, demonstrating for the first time that NOTCH2 and TP53 gene mutations are independent markers of reduced treatment-free and overall survival, respectively.

Dinoflagellate Genomic Organization and Phylogenetic Marker Discovery Utilizing Deep Sequencing Data

Dinoflagellates possess large genomes in which most genes are present in many copies. This has made studies of their genomic organization and phylogenetics challenging. Recent advances in sequencing technology have made deep sequencing of dinoflagellate transcriptomes feasible. This dissertation investigates the genomic organization of dinoflagellates to better understand the challenges of assembling dinoflagellate transcriptomic and genomic data from short read sequencing methods, and develops new techniques that utilize deep sequencing data to identify orthologous genes across a diverse set of taxa. To better understand the genomic organization of dinoflagellates, a genomic cosmid clone of the tandemly repeated gene Alchohol Dehydrogenase (AHD) was sequenced and analyzed. The organization of this clone was found to be counter to prevailing hypotheses of genomic organization in dinoflagellates. Further, a new non-canonical splicing motif was described that could greatly improve the automated modeling and annotation of genomic data. A custom phylogenetic marker discovery pipeline, incorporating methods that leverage the statistical power of large data sets was written. A case study on Stramenopiles was undertaken to test the utility in resolving relationships between known groups as well as the phylogenetic affinity of seven unknown taxa. The pipeline generated a set of 373 genes useful as phylogenetic markers that successfully resolved relationships among the major groups of Stramenopiles, and placed all unknown taxa on the tree with strong bootstrap support. This pipeline was then used to discover 668 genes useful as phylogenetic markers in dinoflagellates. Phylogenetic analysis of 58 dinoflagellates, using this set of markers, produced a phylogeny with good support of all branches. The Suessiales were found to be sister to the Peridinales. The Prorocentrales formed a monophyletic group with the Dinophysiales that was sister to the Gonyaulacales. The Gymnodinales was found to be paraphyletic, forming three monophyletic groups. While this pipeline was used to find phylogenetic markers, it will likely also be useful for finding orthologs of interest for other purposes, for the discovery of horizontally transferred genes, and for the separation of sequences in metagenomic data sets.

Ultra-Deep Sequencing Reveals the microRNA Expression Pattern of the Human Stomach

Background: While microRNAs (miRNAs) play important roles in tissue differentiation and in maintaining basal physiology, little is known about the miRNA expression levels in stomach tissue. Alterations in the miRNA profile can lead to cell deregulation, which can induce neoplasia. Methodology/Principal Findings: A small RNA library of stomach tissue was sequenced using high-throughput SOLiD sequencing technology. We obtained 261,274 quality reads with perfect matches to the human miRnome, and 42% of known miRNAs were identified. Digital Gene Expression profiling (DGE) was performed based on read abundance and showed that fifteen miRNAs were highly expressed in gastric tissue. Subsequently, the expression of these miRNAs was validated in 10 healthy individuals by RT-PCR showed a significant correlation of 83.97% (P<0.05). Six miRNAs showed a low variable pattern of expression (miR-29b, miR-29c, miR-19b, miR-31, miR-148a, miR-451) and could be considered part of the expression pattern of the healthy gastric tissue. Conclusions/Significance: This study aimed to validate normal miRNA profiles of human gastric tissue to establish a reference profile for healthy individuals. Determining the regulatory processes acting in the stomach will be important in the fight against gastric cancer, which is the second-leading cause of cancer mortality worldwide.

Deep sequencing of New World screw-worm transcripts to discover genes involved in insecticide resistance

Background: The New World screw-worm (NWS), Cochliomyia hominivorax, is one of the most important myiasis-causing flies, causing severe losses to the livestock industry. In its current geographical distribution, this species has been controlled by the application of insecticides, mainly organophosphate (OP) compounds, but a number of lineages have been identified that are resistant to such chemicals. Despite its economic importance, only limited genetic information is available for the NWS. Here, as a part of an effort to characterize the C. hominivorax genome and identify putative genes involved in insecticide resistance, we sampled its transcriptome by deep sequencing of polyadenylated transcripts using the 454 sequencing technology. Results: Deep sequencing on the 454 platform of three normalized libraries (larval, adult male and adult female) generated a total of 548,940 reads. Eighteen candidate genes coding for three metabolic detoxification enzyme families, cytochrome P450 monooxygenases, glutathione S transferases and carboxyl/cholinesterases were selected and gene expression levels were measured using quantitative real-time polymerase chain reaction (qRT-PCR). Of the investigated candidates, only one gene was expressed differently between control and resistant larvae with, at least, a 10-fold down-regulation in the resistant larvae. The presence of mutations in the acetylcholinesterase (target site) and carboxylesterase E3 genes was investigated and all of the resistant flies presented E3 mutations previously associated with insecticide resistance. Conclusions: Here, we provided the largest database of NWS expressed sequence tags that is an important resource, not only for further studies on the molecular basis of the OP resistance in NWS fly, but also for functional and comparative studies among Calliphoridae flies. Among our candidates, only one gene was found differentially expressed in resistant individuals, and its role on insecticide resistance should be further investigated. Furthermore, the absence of mutations in the OP target site and the high frequency of mutant carboxylesterase E3 indicate that metabolic resistance mechanisms have evolved predominantly in this species.

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.

Probing the T-cell receptor repertoire with deep sequencing.

PURPOSE OF REVIEW: To review major findings on the T-cell receptor (TCR) repertoire diversity in response to several viral infections based on conventional methods of PCR, cloning and sequencing and to discuss their limitations in light of the recent methodological advances in deep sequencing.¦RECENT FINDINGS: Direct sequencing of TCR expressed by Ag-specific T cells isolated ex vivo has revealed that the TCR repertoire is not as restricted as previously estimated. Furthermore, analyses performed independently of the T-cell clonal hierarchy have brought to light an unexpected diversity. The choice of methods is critical to characterize the complexity of the repertoire. Recent advances in deep sequencing have uncovered the diversity of the TCR repertoire and shown that the size of the repertoire in naive and Ag-experienced memory T cells is three-fold to 15-fold larger than formerly estimated. Interestingly, the TCR complementary determining region 3 sequences are not randomly selected and a certain degree of shared TCR repertoire has been observed between different individuals.¦SUMMARY: Deep sequencing is a major methodological advance allowing more accurate molecular characterization of the TCR repertoire. In the near future, such technologies will further contribute to delineate the complexity of pathogen-specific T-cell response and help defining correlates of a protective immunity.

Bacterial diversity of soil under eucalyptus assessed by 16S rDNA sequencing analysis

Studies on the impact of Eucalyptus spp. on Brazilian soils have focused on soil chemical properties and isolating interesting microbial organisms. Few studies have focused on microbial diversity and ecology in Brazil due to limited coverage of traditional cultivation and isolation methods. Molecular microbial ecology methods based on PCR amplified 16S rDNA have enriched the knowledge of soils microbial biodiversity. The objective of this work was to compare and estimate the bacterial diversity of sympatric communities within soils from two areas, a native forest (NFA) and an eucalyptus arboretum (EAA). PCR primers, whose target soil metagenomic 16S rDNA were used to amplify soil DNA, were cloned using pGEM-T and sequenced to determine bacterial diversity. From the NFA soil 134 clones were analyzed, while 116 clones were analyzed from the EAA soil samples. The sequences were compared with those online at the GenBank. Phylogenetic analyses revealed differences between the soil types and high diversity in both communities. Soil from the Eucalyptus spp. arboretum was found to have a greater bacterial diversity than the soil investigated from the native forest area.

Construction of bovine adenovirus type 3 E1 and E3, substitution plasmids and the sequencing analysis of DNA pol and pTP /

The relative ease to concentrate and purify adenoviruses, their well characterized mid-sized genome, and the ability to delete non-essential regions from their genome to accommodate foreign gene, made adenoviruses a suitable candidate for the construction of vectors. The use of adenoviral vectors in gene therapy, vaccination, and as a general vector system for expressing foreign genes have been documented for some time. In this study, the objective was to rescue a BAV3 E1 or E3 recombinant vector carrying the kanamycin resistant gene, a dominant selectable marker with useful applications in studying vectored gene expression in mammalian cells. To accomplish the objective of this study, more information about BAV3 DNA sequences was required in order to make the manipulation of the virus genome accessible. Therefore, sequencing of the BAV3 genome from 1 1 .7% to 30.8% was carried out. Analysis of the determined sequences revealed the primary structure of important viral gene products coded by E2 including BAV3 DNA pol and precursor to terminal protein. Comparative analysis of these proteins with their counterparts from human and non human adenoviruses revealed important insights as to the evolutionary lineage of BAV3. In order to insert the kanamycin resistance gene in either E1 or E3, it was necessary to delete BAV3 sequences to accommodate the foreign gene so as not to exceed the limit of the packaging capacity of the virus. To construct a recombinant BAV3 in which a foreign gene was inserted in the deleted E1 region, an E1 shuttle vector was constructed. This involved the deletion from the viral sequences a region between 1.3% to 9% and inserting the kanamycin resistance gene to replace the deletion. The E1 shuttle vector contained the left (0%- 53.9%) segment of the genome and was expected to generate BAV3 recombinants that can be grown and propagated in cells that can complement the missing E1 functions. To construct a similar shuttle vector for E3 deletion, DNA sequences extending from 78.9% to 82.5% (1281 bp) were deleted from within the E3 region that had been cloned into a plasmid vector. The deleted region corresponds to those that have been shown to be non-essential for viral replication in cell culture. The resulting plasmid was used to construct another recombinant plasmid with BAV3 DNA sequences extending from 37.1% to 100% and with a deletion of E3 sequences that were replaced by kanamycin resistance gene. This shuttle plasmid was used in cotransfections with digested viral DNA in an attempt to rescue a recombinant BAV3 carrying the kanamycin resistance gene to replace the deleted E3. In spite of repeated attempts of transfection, El or E3 recombinant BAV3 were not isolated. It seems that other approaches should be applied to make a final conclusion on BAV3 infectivity.

Suitability of DNA extracted from archival specimens of fruit-eating bats of the genus Artibeus (Chiroptera, Phyllostomidae) for polymerase chain reaction and sequencing analysis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Bacterial diversity of soil under eucalyptus assessed by 16S rDNA sequencing analysis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Comparative clinical sample preparation of DNA and RNA viral nucleic acids for a commercial deep sequencing system (Illumina MiSeq (R))

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Whole genome sequencing analysis of Plasmodium vivax using whole genome capture

Background: Malaria caused by Plasmodium vivax is an experimentally neglected severe disease with a substantial burden on human health. Because of technical limitations, little is known about the biology of this important human pathogen. Whole genome analysis methods on patient-derived material are thus likely to have a substantial impact on our understanding of P. vivax pathogenesis and epidemiology. For example, it will allow study of the evolution and population biology of the parasite, allow parasite transmission patterns to be characterized, and may facilitate the identification of new drug resistance genes. Because parasitemias are typically low and the parasite cannot be readily cultured, on-site leukocyte depletion of blood samples is typically needed to remove human DNA that may be 1000X more abundant than parasite DNA. These features have precluded the analysis of archived blood samples and require the presence of laboratories in close proximity to the collection of field samples for optimal pre-cryopreservation sample preparation. Results: Here we show that in-solution hybridization capture can be used to extract P. vivax DNA from human contaminating DNA in the laboratory without the need for on-site leukocyte filtration. Using a whole genome capture method, we were able to enrich P. vivax DNA from bulk genomic DNA from less than 0.5% to a median of 55% (range 20%-80%). This level of enrichment allows for efficient analysis of the samples by whole genome sequencing and does not introduce any gross biases into the data. With this method, we obtained greater than 5X coverage across 93% of the P. vivax genome for four P. vivax strains from Iquitos, Peru, which is similar to our results using leukocyte filtration (greater than 5X coverage across 96% of the genome). Conclusion: The whole genome capture technique will enable more efficient whole genome analysis of P. vivax from a larger geographic region and from valuable archived sample collections.