Whole-Genome Characterization of Human and Simian Immunodeficiency Virus Intrahost Diversity by Ultradeep Pyrosequencing

Rapid evolution and high intrahost sequence diversity are hallmarks of human and simian immunodeficiency virus (HIV/SIV) infection. Minor viral variants have important implications for drug resistance, receptor tropism, and immune evasion. Here, we used ultradeep pyrosequencing to sequence complete HIV/SIV genomes, detecting variants present at a frequency as low as 1%. This approach provides a more complete characterization of the viral population than is possible with conventional methods, revealing low-level drug resistance and detecting previously hidden changes in the viral population. While this work applies pyrosequencing to immunodeficiency viruses, this approach could be applied to virtually any viral pathogen.

Mapping of the healthy immunoglobulin repertoire over the course of B cell maturation by deep pyrosequencing

Healthy immunoglobulin repertoire has not been extensively evaluated reflecting in part the challenge of generating sufficiently robust data sets by conventional clonal sequencing. Deep sequencing has revolutionized the capacity to evaluate the depth and breadth of the Ig repertoire along the B cell developmental pathway, and can be used to pin point defect(s) of primary or acquired B-cell associated diseases. In this study healthy IgM and IgG repertoires were studied by 454-pyrosequencing to establish the healthy controls for diseased repertoires. (...)

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.

Clinical pharmacogenomic testing of KRAS, BRAF and EGFR mutations by high resolution melting analysis and ultra-deep pyrosequencing

BACKGROUND: Epidermal growth factor receptor (EGFR) and its downstream factors KRAS and BRAF are mutated in several types of cancer, affecting the clinical response to EGFR inhibitors. Mutations in the EGFR kinase domain predict sensitivity to the tyrosine kinase inhibitors gefitinib and erlotinib in lung adenocarcinoma, while activating point mutations in KRAS and BRAF confer resistance to the anti-EGFR monoclonal antibody cetuximab in colorectal cancer. The development of new generation methods for systematic mutation screening of these genes will allow more appropriate therapeutic choices. METHODS: We describe a high resolution melting (HRM) assay for mutation detection in EGFR exons 19-21, KRAS codon 12/13 and BRAF V600 using formalin-fixed paraffin-embedded samples. Somatic variation of KRAS exon 2 was also analysed by massively parallel pyrosequencing of amplicons with the GS Junior 454 platform. RESULTS: We tested 120 routine diagnostic specimens from patients with colorectal or lung cancer. Mutations in KRAS, BRAF and EGFR were observed in 41.9%, 13.0% and 11.1% of the overall samples, respectively, being mutually exclusive. For KRAS, six types of substitutions were detected (17 G12D, 9 G13D, 7 G12C, 2 G12A, 2 G12V, 2 G12S), while V600E accounted for all the BRAF activating mutations. Regarding EGFR, two cases showed exon 19 deletions (delE746-A750 and delE746-T751insA) and another two substitutions in exon 21 (one showed L858R with the resistance mutation T590M in exon 20, and the other had P848L mutation). Consistent with earlier reports, our results show that KRAS and BRAF mutation frequencies in colorectal cancer were 44.3% and 13.0%, respectively, while EGFR mutations were detected in 11.1% of the lung cancer specimens. Ultra-deep amplicon pyrosequencing successfully validated the HRM results and allowed detection and quantitation of KRAS somatic mutations. CONCLUSIONS: HRM is a rapid and sensitive method for moderate-throughput cost-effective screening of oncogene mutations in clinical samples. Rather than Sanger sequence validation, next-generation sequencing technology results in more accurate quantitative results in somatic variation and can be achieved at a higher throughput scale.

Ultra-Deep Pyrosequencing Detects Conserved Genomic Sites and Quantifies Linkage of Drug-Resistant Amino Acid Changes in the Hepatitis B Virus Genome

Selection of amino acid substitutions associated with resistance to nucleos(t)ide-analog (NA) therapy in the hepatitis B virus (HBV) reverse transcriptase (RT) and their combination in a single viral genome complicates treatment of chronic HBV infection and may affect the overlapping surface coding region. In this study, the variability of an overlapping polymerase-surface region, critical for NA resistance, is investigated before treatment and under antiviral therapy, with assessment of NA-resistant amino acid changes simultaneously occurring in the same genome (linkage analysis) and their influence on the surface coding region.

High-throughput sequence analysis of turbot (Scophthalmus maximus) transcriptome using 454-pyrosequencing for the discovery of antiviral immune genes

Turbot (Scophthalmus maximus L.) is an important aquacultural resource both in Europe and Asia. However, there is little information on gene sequences available in public databases. Currently, one of the main problems affecting the culture of this flatfish is mortality due to several pathogens, especially viral diseases which are not treatable. In order to identify new genes involved in immune defense, we conducted 454-pyrosequencing of the turbot transcriptome after different immune stimulations.

Use of pyrosequencing to identify streptococci and to detect mutations causing antimicrobial resistance

Rapid identification and resistance determination of pathogens in clinical specimens is vital for accurate treatment and monitoring of infectious diseases. Antimicrobial drug resistance is increasing globally and healthcare settings are facing this cost-intensive and even life-threatening problem. The incidence of resistant pathogens in Finland has remained relatively steady and manageable at least for the time being. DNA sequencing is the gold standard method for genotyping, mutation analysis, and identification of bacteria. Due to significant cost decrease in recent years, this technique is available to many research and clinical laboratories. Pyrosequencing technique, a rapid real-time DNA sequencing method especially suitable for analyzing fairly short stretches of DNA, was used in this study. Due to its robustness and versatility, pyrosequencing was applied in this study for identification of streptococci and detection of certain mutations causing antimicrobial resistance in different bacteria. Certain streptococcal species such as S. pneumoniae and S. pyogenes are significantly important clinical pathogens. S. pneumoniae causes e.g. pneumonia and otitis media and is one of the most important community-acquired pathogens. S. pyogenes, also known as group A streptococcus, causes e.g. angina and erysipelas. In contrast, the socalled alpha-haemolytic streptococci, such as S. mitis and S. oralis, belong to the normal microbiota, which are regarded to be non-pathogenic and are nearly impossible to identify by phenotypic methods. In this thesis, a pyrosequencing method was developed for identification of streptococcal species based on the 16S rRNA sequences. Almost all streptococcal species could be differentiated from one another by the developed method, including S. pneumoniae from its close relatives S. mitis and S. oralis . New resistance genes and their variants are constantly discovered and reported. In this study, new methods for detecting certain mutations causing macrolide resistance or extended spectrum beta-lactamase (ESBL) phenotype were developed. These resistance detection approaches are not only suitable for surveillance of mechanisms causing antimicrobial resistance but also for routine analysis of clinical samples particularly in epidemic settings. In conclusion, pyrosequencing was found to be an accurate, versatile, cost-effective, and rapid DNA sequencing method that is especially suitable for mutation analysis of short DNA fragments and identification of certain bacteria.

Development of 18 Polymorphic Microsatellite Markers for Vinca minor (Apocynaceae) via 454 Pyrosequencing

• Premise of the study: Polymorphic microsatellite markers were developed in Vinca minor (Apocynaceae) to evaluate the level of clonality, population structure, and genetic diversity of the species within its native and introduced range. • Methods and Results: A total of 1371 microsatellites were found in 43,565 reads from 454 pyrosequencing of genomic V. minor DNA. Additional microsatellite loci were mined from publicly available cDNA sequences. After several rounds of screening, 18 primer pairs flanking di-, tri-, or tetranucleotide repeats were identified that revealed high levels of genetic diversity in two native Italian populations, with two to 11 alleles per locus. Clonal growth predominated in two populations from the introduced range in Germany. Five loci successfully cross-amplified in three additional Vinca species. • Conclusions: The novel polymorphic microsatellite markers are promising tools for studying clonality and population genetics of V. minor and for assessing the historical origin of Central European populations.

Transcriptome analysis for discovering candidate genes involve in embryogenesis in coconut (Cocos nucifera L.) through 454 pyrosequencing

Coconut, Cocos nucifera L. is a major plantation crop, which ensures income for millions of people in the tropical region. Detailed molecular studies on zygotic embryo development would provide valuable clues for the identification of molecular markers to improve somatic embryogenesis. Since there is no ongoing genome project for this species, coconut expressed sequence tags (EST) would be an interesting technique to identify important coconut embryo specific genes as well as other functional genes in different biochemical pathways. The goal of this study was to analyse the ESTs by examining the transcriptome data of the different embryo tissue types together with one somatic tissue. Here, four cDNA libraries from immature embryo, mature embryo, microspore derived embryo and mature leaves were constructed. cDNA was sequenced by the Roche-454 GS-FLX system and assembled into 32621 putative unigenes and 155017 singletons. Of these unigenes, 18651 had significant sequence similarities to non-redundant protein database, from which 16153 were assigned to one or more gene ontology categories. Homologue genes, which are responsible for embryo development such as chitinase, beta-1,3-glucanase, ATP synthase CF0 subunit, thaumatin-like protein and metallothionein-like protein were identified among the embryo EST collection. Of the unigenes, 6694 were mapped into 139 KEGG pathways including carbohydrate metabolism, energy metabolism, lipid metabolism, amino acid metabolism and nucleotide metabolism. This collection of 454-derived EST data generated from different tissue types provides a significant resource for genome wide studies and gene discovery of coconut, a non-model species.

Use of Pyrosequencing and Denaturing Gradient Gel Electrophoresis to Examine the Effects of Probiotics and Essential Oil Blends on Digestive Microflora in Broilers Under Mixed Eimeria Infection

A protective digestive microflora helps prevent and reduce broiler infection and colonization by enteropathogens. In the current experiment, broilers fed diets supplemented with probiotics and essential oil (EO) blends were infected with a standard mixed Eimeria spp. to determine effects of performance enhancers on ileal and cecal microbial communities (MCs). Eight treatment groups included four controls (uninfected-unmedicated [UU], unmedicated-infected, the antibiotic BMD plus the ionophore Coban as positive control, and the ionophore as negative control), and four treatments (probiotics BC-30 and Calsporin; and EO, Crina Poultry Plus, and Crina PoultryAF). Day-old broilers were raised to 14 days in floor pens on used litter and then were moved to Petersime batteries and inoculated at 15 days with mixed Eimeria spp. Ileal and cecal samples were collected at 14 days and 7 days postinfection. Digesta DNA was subjected to pyrosequencing for sequencing of individual cecal bacteria and denaturing gradient gel electrophoresis (DGGE) for determination of changes in ileal and cecal MC according to percentage similarity coefficient (%SC). Pyrosequencing is very sensitive detecting shifts in individual bacterial sequences, whereas DGGE is able to detect gross shifts in entire MC. These combined techniques offer versatility toward identifying feed additive and mild Eimeria infection modulation of broiler MC. Pyrosequencing detected 147 bacterial species sequences. Additionally, pyrosequencing revealed the presence of relatively low levels of the potential human enteropathogens Campylobacter sp. and four Shigella spp. as well as the potential poultry pathogen Clostridiun perfringens. Pre- and postinfection changes in ileal (56%SC) and cecal (78.5%SC) DGGE profiles resulted from the coccidia infection and with increased broiler age. Probiotics and EO changed MC from those seen in UU ilea and ceca. Results potentially reflect the performance enhancement above expectations in comparison to broilers not given the probiotics or the specific EO blends as feed supplements.

Development of microsatellite markers for the Neotropical endemic Brazilian Guanabara frog, Euparkerella brasiliensis, through 454 shotgun pyrosequencing

The new-generation 454 GS-FLX Titanium pyrosequencing was used to isolate microsatellite markers for the Brazilian Guanabara frog, Euparkerella brasiliensis, an Atlantic forest endemic species. Three multiplex polymerase chain reaction sets were optimized for genotyping of 11 polymorphic (di- and tetranucleotide) microsatellite markers. Genetic diversity was assessed in 21 individuals from a population (Reserva Ecológica de Guapiaçu, REGUA) locatedin the central region of the Rio de Janeiro State, in Brazil. The mean number of alleles per locus ranged from 3 to 12. Observed and expected heterozygosities ranged from 0.095 to 0.905 and from 0.094 to 0.904, respectively. After using the Bonferroni correction for multiple tests, there was no evidence of linkage disequilibrium between pairs of loci but deviations for Hardy-Weinberg equilibrium were found in 4 loci. We found no evidence for allele dropouts or stuttering, but we detected the presence of null alleles at loci Eb10 and Eb36. These markers will be useful for analyses of fine-scale population structure and determination of relative effects of habitat loss and fragmentation on population genetic variability within species. © FUNPEC-RP.

Microsatellite loci in the oceanic whitetip shark and cross-species amplification using pyrosequencing technology

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

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.