Low-Cost Ultra-Wide Genotyping Using Roche/454 Pyrosequencing for Surveillance of HIV Drug Resistance

Background: Great efforts have been made to increase accessibility of HIV antiretroviral therapy (ART) in low and middle-income countries. The threat of wide-scale emergence of drug resistance could severely hamper ART scale-up efforts. Population-based surveillance of transmitted HIV drug resistance ensures the use of appropriate first-line regimens to maximize efficacy of ART programs where drug options are limited. However, traditional HIV genotyping is extremely expensive, providing a cost barrier to wide-scale and frequent HIV drug resistance surveillance. Methods/Results: We have developed a low-cost laboratory-scale next-generation sequencing-based genotyping method to monitor drug resistance. We designed primers specifically to amplify protease and reverse transcriptase from Brazilian HIV subtypes and developed a multiplexing scheme using multiplex identifier tags to minimize cost while providing more robust data than traditional genotyping techniques. Using this approach, we characterized drug resistance from plasma in 81 HIV infected individuals collected in Sao Paulo, Brazil. We describe the complexities of analyzing next-generation sequencing data and present a simplified open-source workflow to analyze drug resistance data. From this data, we identified drug resistance mutations in 20% of treatment naive individuals in our cohort, which is similar to frequencies identified using traditional genotyping in Brazilian patient samples. Conclusion: The developed ultra-wide sequencing approach described here allows multiplexing of at least 48 patient samples per sequencing run, 4 times more than the current genotyping method. This method is also 4-fold more sensitive (5% minimal detection frequency vs. 20%) at a cost 3-5 x less than the traditional Sanger-based genotyping method. Lastly, by using a benchtop next-generation sequencer (Roche/454 GS Junior), this approach can be more easily implemented in low-resource settings. This data provides proof-of-concept that next-generation HIV drug resistance genotyping is a feasible and low-cost alternative to current genotyping methods and may be particularly beneficial for in-country surveillance of transmitted drug resistance.

Peptide fingerprinting of the neurotoxic fractions isolated from the secretions of sea anemones Stichodactyla helianthus and Bunodosoma granulifera. New members of the APETx-like family identified by a 454 pyrosequencing approach

Sea anemones are known to contain a wide diversity of biologically active peptides, mostly unexplored according to recent peptidomic and transcriptomic studies. In the present work, the neurotoxic fractions from the exudates of Stichodactyla helianthus and Bunodosoma granulifera were analyzed by reversed-phase chromatography and mass spectrometry. The first peptide fingerprints of these sea anemones were assessed, revealing the largest number of peptide components (156) so far found in sea anemone species, as well as the richer peptide diversity of B. granulifera in relation to S. helianthus. The transcriptomic analysis of B. granulifera, performed by massive cDNA sequencing with 454 pyrosequencing approach allowed the discovery of five new APETx-like peptides (U-AITX-Bg1a-e - including the full sequences of their precursors for four of them), which together with type 1 sea anemone sodium channel toxins constitute a very distinguishable feature of studied sea anemone species belonging to genus Bunodosoma. The molecular modeling of these new APETx-like peptides showed a distribution of positively charged and aromatic residues in putative contact surfaces as observed in other animal toxins. On the other hand, they also showed variable electrostatic potentials, thus suggesting a docking onto their targeted channels in different spatial orientations. Moreover several crab paralyzing toxins (other than U-AITX-Bg1a-e), which induce a variety of symptoms in crabs, were isolated. Some of them presumably belong to new classes of crab-paralyzing peptide toxins, especially those with molecular masses below 2 kDa, which represent the smallest peptide toxins found in sea anemones. (C) 2011 Elsevier Inc. All rights reserved.

Analysis of Hepatitis C Virus Intrahost Diversity across the Coding Region by Ultradeep Pyrosequencing

Hepatitis C virus (HCV) is the leading cause of liver disease worldwide. In this study, we analyzed four treatment-naive patients infected with subtype 1a and performed Roche/454 pyrosequencing across the coding region. We report the presence of low-level drug resistance mutations that would most likely have been missed using conventional sequencing methods. The approach described here is broadly applicable to studies of viral diversity and could help to improve the efficacy of direct-acting antiviral agents (DAA) in the treatment of HCV-infected patients.

A pyrosequencing assay for the rapid discrimination of mitochondrial lineages in the Salmo trutta species complex

The European trout (Salmo trutta species complex) is genetically very diverse consisting of five distinct mitochondrial lineages that probably originated in the Pleistocene. Here, we describe a novel pyrosequencing protocol to generate two short sequence reads from the mitochondrial control region, which allow the unambiguous identification of all five lineages. The approach was found to be easily transferable between laboratories and should be a valuable tool for the assessment of genetic diversity in trout. Pyrosequencing-based assays for molecular species identification are expected to be generally useful whenever multiple positions in a short DNA sequence need to be assessed.

Characterization of CDKN2A(p16) methylation and impact in colorectal cancer: systematic analysis using pyrosequencing

Background The aim of this study is to analyse CDKN2A methylation using pyrosequencing on a large cohort of colorectal cancers and corresponding non-neoplastic tissues. In a second step, the effect of methylation on clinical outcome is addressed. Methods Primary colorectal cancers and matched non-neoplastic tissues from 432 patients underwent CDKN2A methylation analysis by pyrosequencing (PyroMarkQ96). Methylation was then related to clinical outcome, microsatellite instability (MSI), and BRAF and KRAS mutation. Different amplification conditions (35 to 50 PCR cycles) using a range of 0-100% methylated DNA were tested. Results Background methylation was at most 10% with ≥35 PCR cycles. Correlation of observed and expected values was high, even at low methylation levels (0.02%, 0.6%, 2%). Accuracy of detection was optimal with 45 PCR cycles. Methylation in normal mucosa ranged from 0 to >90% in some cases. Based on the maximum value of 10% background, positivity was defined as a ≥20% difference in methylation between tumor and normal tissue, which occurred in 87 cases. CDKN2A methylation positivity was associated with MSI (p = 0.025), BRAF mutation (p < 0.0001), higher tumor grade (p < 0.0001), mucinous histology (p = 0.0209) but not with KRAS mutation. CDKN2A methylation had an independent adverse effect (p = 0.0058) on prognosis. Conclusion The non-negligible CDKN2A methylation of normal colorectal mucosa may confound the assessment of tumor-specific hypermethylation, suggesting that corresponding non-neoplastic tissue should be used as a control. CDKN2A methylation is robustly detected by pyrosequencing, even at low levels, suggesting that this unfavorable prognostic biomarker warrants investigation in prospective studies.

Spatial Distribution of Cryptic Species Diversity in European Freshwater Amphipods (Gammarus fossarum) as Revealed by Pyrosequencing

In order to understand and protect ecosystems, local gene pools need to be evaluated with respect to their uniqueness. Cryptic species present a challenge in this context because their presence, if unrecognized, may lead to serious misjudgement of the distribution of evolutionarily distinct genetic entities. In this study, we describe the current geographical distribution of cryptic species of the ecologically important stream amphipod Gammarus fossarum (types A, B and C). We use a novel pyrosequencing assay for molecular species identification and survey 62 populations in Switzerland, plus several populations in Germany and eastern France. In addition, we compile data from previous publications (mainly Germany). A clear transition is observed from type A in the east (Danube and Po drainages) to types B and, more rarely, C in the west (Meuse, Rhone, and four smaller French river systems). Within the Rhine drainage, the cryptic species meet in a contact zone which spans the entire G. fossarum distribution range from north to south. This large-scale geographical sorting indicates that types A and B persisted in separate refugia during Pleistocene glaciations. Within the contact zone, the species rarely co-occur at the same site, suggesting that ecological processes may preclude long-term coexistence. The clear phylogeographical signal observed in this study implies that, in many parts of Europe, only one of the cryptic species is present.

Low-frequency drug-resistant HIV-1 and risk of virological failure to first-line NNRTI-based ART: a multicohort European case-control study using centralized ultrasensitive 454 pyrosequencing.

OBJECTIVES It is still debated if pre-existing minority drug-resistant HIV-1 variants (MVs) affect the virological outcomes of first-line NNRTI-containing ART. METHODS This Europe-wide case-control study included ART-naive subjects infected with drug-susceptible HIV-1 as revealed by population sequencing, who achieved virological suppression on first-line ART including one NNRTI. Cases experienced virological failure and controls were subjects from the same cohort whose viraemia remained suppressed at a matched time since initiation of ART. Blinded, centralized 454 pyrosequencing with parallel bioinformatic analysis in two laboratories was used to identify MVs in the 1%-25% frequency range. ORs of virological failure according to MV detection were estimated by logistic regression. RESULTS Two hundred and sixty samples (76 cases and 184 controls), mostly subtype B (73.5%), were used for the analysis. Identical MVs were detected in the two laboratories. 31.6% of cases and 16.8% of controls harboured pre-existing MVs. Detection of at least one MV versus no MVs was associated with an increased risk of virological failure (OR = 2.75, 95% CI = 1.35-5.60, P = 0.005); similar associations were observed for at least one MV versus no NRTI MVs (OR = 2.27, 95% CI = 0.76-6.77, P = 0.140) and at least one MV versus no NNRTI MVs (OR = 2.41, 95% CI = 1.12-5.18, P = 0.024). A dose-effect relationship between virological failure and mutational load was found. CONCLUSIONS Pre-existing MVs more than double the risk of virological failure to first-line NNRTI-based ART.

Anthropogenic impact on diazotrophic diversity in the mangrove rhizosphere revealed by nifH pyrosequencing

Diazotrophs in the mangrove rhizosphere play a major role in providing new nitrogen to the mangrove ecosystem and their composition and activity are strongly influenced by anthropogenic activity and ecological conditions. In this study, the diversity of the diazotroph communities in the rhizosphere sediment of five tropical mangrove sites with different levels of pollution along the north and south coastline of Singapore were studied by pyrosequencing of the nifH gene. Bioinformatics analysis revealed that in all the studied locations, the diazotroph communities comprised mainly of members of the diazotrophic cluster I and cluster III. The detected cluster III diazotrophs, which were composed entirely of sulfate-reducing bacteria, were more abundant in the less polluted locations. The metabolic capacities of these diazotrophs indicate the potential for bioremediation and resiliency of the ecosystem to anthropogenic impact. In heavily polluted locations, the diazotrophic community structures were markedly different and the diversity of species was significantly reduced when compared with those in a pristine location. This, together with the increased abundance of Marinobacterium, which is a bioindicator of pollution, suggests that anthropogenic activity has a negative impact on the genetic diversity of diazotrophs in the mangrove rhizosphere.

Massively parallel sequencing and analysis of expressed sequence tags in a successful invasive plant

Background Invasive species pose a significant threat to global economies, agriculture and biodiversity. Despite progress towards understanding the ecological factors associated with plant invasions, limited genomic resources have made it difficult to elucidate the evolutionary and genetic factors responsible for invasiveness. This study presents the first expressed sequence tag (EST) collection for Senecio madagascariensis, a globally invasive plant species. Methods We used pyrosequencing of one normalized and two subtractive libraries, derived from one native and one invasive population, to generate an EST collection. ESTs were assembled into contigs, annotated by BLAST comparison with the NCBI non-redundant protein database and assigned gene ontology (GO) terms from the Plant GO Slim ontologies. Key Results Assembly of the 221 746 sequence reads resulted in 12 442 contigs. Over 50 % (6183) of 12 442 contigs showed significant homology to proteins in the NCBI database, representing approx. 4800 independent transcripts. The molecular transducer GO term was significantly over-represented in the native (South African) subtractive library compared with the invasive (Australian) library. Based on NCBI BLAST hits and literature searches, 40 % of the molecular transducer genes identified in the South African subtractive library are likely to be involved in response to biotic stimuli, such as fungal, bacterial and viral pathogens. Conclusions This EST collection is the first representation of the S. madagascariensis transcriptome and provides an important resource for the discovery of candidate genes associated with plant invasiveness. The over-representation of molecular transducer genes associated with defence responses in the native subtractive library provides preliminary support for aspects of the enemy release and evolution of increased competitive ability hypotheses in this successful invasive. This study highlights the contribution of next-generation sequencing to better understanding the molecular mechanisms underlying ecological hypotheses that are important in successful plant invasions.

Transcriptomics approaches to identify genes and markers for production traits in giant freshwater prawn Macrobrachium rosenbergii

This study used next generation sequencing technologies to investigate growth in a cultured crustacean. The objective was to identify and characterise specific gene loci that contribute important phenotypic variation to growth as well as to develop a large set of SNP markers in candidate genes for assessing correlations between specific mutations and individual growth performance. The genomic dataset generated provides a fundamental resource for application in future crustacean stock improvement programs. Ultimately, the data can be applied to development of culture lines with improved growth performance.

Gemcitabine reactivates epigenetically silenced genes and functions as a DNA methyltransferase inhibitor

Gemcitabine is indicated in combination with cisplatin as first-line therapy for solid tumours including non-small cell lung cancer (NSCLC), bladder cancer and mesothelioma. Gemcitabine is an analogue of pyrimidine cytosine and functions as an anti-metabolite. Structurally, however, gemcitabine has similarities to 5-aza-2-deoxycytidine (decitabine/Dacogen®), a DNA methyltransferase inhibitor (DNMTi). NSCLC, mesothelioma and prostate cancer cell lines were treated with decitabine and gemcitabine. Reactivation of epigenetically silenced genes was examined by RT-PCR/qPCR. DNA methyltransferase activity in nuclear extracts and recombinant proteins was measured using a DNA methyltransferase assay, and alterations in DNA methylation status were examined using methylation-specific PCR (MS-PCR) and pyrosequencing. We observe a reactivation of several epigenetically silenced genes including GSTP1, IGFBP3 and RASSF1A. Gemcitabine functionally inhibited DNA methyltransferase activity in both nuclear extracts and recombinant proteins. Gemcitabine dramatically destabilised DNMT1 protein. However, DNA CpG methylation was for the most part unaffected by gemcitabine. In conclusion, gemcitabine both inhibits and destabilises DNA methyltransferases and reactivates epigenetically silenced genes having activity equivalent to decitabine at concentrations significantly lower than those achieved in the treatment of patients with solid tumours. This property may contribute to the anticancer activity of gemcitabine.

BDNF and TNF-α polymorphisms in memory

Here, we investigate the genetic basis of human memory in healthy individuals and the potential role of two polymorphisms, previously implicated in memory function. We have explored aspects of retrospective and prospective memory including semantic, short term, working and long-term memory in conjunction with brain derived neurotrophic factor (BDNF) and tumor necrosis factor-alpha (TNF-alpha). The memory scores for healthy individuals in the population were obtained for each memory type and the population was genotyped via restriction fragment length polymorphism for the BDNF rs6265 (Val66Met) SNP and via pyrosequencing for the TNF-alpha rs113325588 SNP. Using univariate ANOVA, a significant association of the BDNF polymorphism with visual and spatial memory retention and a significant association of the TNF-alpha polymorphism was observed with spatial memory retention. In addition, a significant interactive effect between BDNF and TNF-alpha polymorphisms was observed in spatial memory retention. In practice visual memory involves spatial information and the two memory systems work together, however our data demonstrate that individuals with the Val/Val BDNF genotype have poorer visual memory but higher spatial memory retention, indicating a level of interaction between TNF-alpha and BDNF in spatial memory retention. This is the first study to use genetic analysis to determine the interaction between BDNF and TNF-alpha in relation to memory in normal adults and provides important information regarding the effect of genetic determinants and gene interactions on human memory.

Computational epigenetic profiling of CpG islets in MTHFR

Computational epigenetics is a new area of research focused on exploring how DNA methylation patterns affect transcription factor binding that affect gene expression patterns. The aim of this study was to produce a new protocol for the detection of DNA methylation patterns using computational analysis which can be further confirmed by bisulfite PCR with serial pyrosequencing. The upstream regulatory element and pre-initiation complex relative to CpG islets within the methylenetetrahydrofolate reductase gene were determined via computational analysis and online databases. The 1,104 bp long CpG island located near to or at the alternative promoter site of methylenetetrahydrofolate reductase gene was identified. The CpG plot indicated that CpG islets A and B, within the island, contained 62 and 75 % GC content CpG ratios of 0.70 and 0.80–0.95, respectively. Further exploration of the CpG islets A and B indicates that the transcription start sites were GGC which were absent from the TATA boxes. In addition, although six PROSITE motifs were identified in CpG B, no motifs were detected in CpG A. A number of cis-regulatory elements were found in different regions within the CpGs A and B. Transcription factors were predicted to bind to CpGs A and B with varying affinities depending on the DNA methylation status. In addition, transcription factor binding may influence the expression patterns of the methylenetetrahydrofolate reductase gene by recruiting chromatin condensation inducing factors. These results have significant implications for the understanding of the architecture of transcription factor binding at CpG islets as well as DNA methylation patterns that affect chromatin structure.

IgE-Associated IGHV Genes from Venom and Peanut Allergic Individuals Lack Mutational Evidence of Antigen Selection

Antigen selection of B cells within the germinal center reaction generally leads to the accumulation of replacement mutations in the complementarity-determining regions (CDRs) of immunoglobulin genes. Studies of mutations in IgE-associated VDJ gene sequences have cast doubt on the role of antigen selection in the evolution of the human IgE response, and it may be that selection for high affinity antibodies is a feature of some but not all allergic diseases. The severity of IgE-mediated anaphylaxis is such that it could result from higher affinity IgE antibodies. We therefore investigated IGHV mutations in IgE-associated sequences derived from ten individuals with a history of anaphylactic reactions to bee or wasp venom or peanut allergens. IgG sequences, which more certainly experience antigen selection, served as a control dataset. A total of 6025 unique IgE and 5396 unique IgG sequences were generated using high throughput 454 pyrosequencing. The proportion of replacement mutations seen in the CDRs of the IgG dataset was significantly higher than that of the IgE dataset, and the IgE sequences showed little evidence of antigen selection. To exclude the possibility that 454 errors had compromised analysis, rigorous filtering of the datasets led to datasets of 90 core IgE sequences and 411 IgG sequences. These sequences were present as both forward and reverse reads, and so were most unlikely to include sequencing errors. The filtered datasets confirmed that antigen selection plays a greater role in the evolution of IgG sequences than of IgE sequences derived from the study participants.

Identification of bacterial community composition in freshwater aquaculture system farming of Litopenaeus vannamei reveals distinct temperature-driven patterns

Change in temperature is often a major environmental factor in triggering waterborne disease outbreaks. Previous research has revealed temporal and spatial patterns of bacterial population in several aquatic ecosystems. To date, very little information is available on aquaculture environment. Here, we assessed environmental temperature effects on bacterial community composition in freshwater aquaculture system farming of Litopenaeus vannamei (FASFL). Water samples were collected over a one-year period, and aquatic bacteria were characterized by polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) and 16S rDNA pyrosequencing. Resulting DGGE fingerprints revealed a specific and dynamic bacterial population structure with considerable variation over the seasonal change, suggesting that environmental temperature was a key driver of bacterial population in the FASFL. Pyrosequencing data further demonstrated substantial difference in bacterial community composition between the water at higher (WHT) and at lower (WLT) temperatures in the FASFL. Actinobacteria, Proteobacteria and Bacteroidetes were the highest abundant phyla in the FASFL, however, a large number of unclassified bacteria contributed the most to the observed variation in phylogenetic diversity. The WHT harbored remarkably higher diversity and richness in bacterial composition at genus and species levels when compared to the WLT. Some potential pathogenenic species were identified in both WHT and WLT, providing data in support of aquatic animal health management in the aquaculture industry.