Restricted MHC class I A locus diversity in olive and hybrid olive/yellow baboons from the Southwest National Primate Research Center.

Baboons are valuable models for complex human diseases due to their genetic and physiologic similarities to humans. Deep sequencing methods to characterize full-length major histocompatibility complex (MHC) class I (MHC-I) alleles in different nonhuman primate populations were used to identify novel MHC-I alleles in baboons. We combined data from Illumina MiSeq sequencing and Roche/454 sequencing to characterize novel full-length MHC-I transcripts in a cohort of olive and hybrid olive/yellow baboons from the Southwest National Primate Research Center (SNPRC). We characterized 57 novel full-length alleles from 24 baboons and found limited genetic diversity at the MHC-I A locus, with significant sharing of two MHC-I A lineages between 22 out of the 24 animals characterized. These shared alleles provide the basis for development of tools such as MHC:peptide tetramers for studying cellular immune responses in this important animal model.

Addition of plant-growth-promoting Bacillus subtilis PTS-394 on tomato rhizosphere has no durable impact on composition of root microbiome.

BACKGROUND: Representatives of the genus Bacillus are increasingly used in agriculture to promote plant growth and to protect against plant pathogens. Unfortunately, hitherto the impact of Bacillus inoculants on the indigenous plant microbiota has been investigated exclusively for the species Bacillus amyloliquefaciens and was limited to prokaryotes, whilst eukaryotic member of this community, e.g. fungi, were not considered. RESULTS: The root-colonizing Bacillus subtilis PTS-394 supported growth of tomato plants and suppressed soil-borne diseases. Roche 454 pyrosequencing revealed that PTS-394 has only a transient impact on the microbiota community of the tomato rhizosphere. The impact on eukaryota could last up to 14 days, while that on bacterial communities lasted for 3 days only. CONCLUSIONS: Ecological adaptation and microbial community-preserving capacity are important criteria when assessing suitability of bio-inoculants for commercial development. As shown here, B. subtilis PTS-394 is acting as an environmentally compatible plant protective agent without permanent effects on rhizosphere microbial community.

Clinical Oral Condition Analysis and the Influence of Highly Active Antiretroviral Therapy on Human Salivary Microbial Community Diversity in HIV-Infected/AIDS Patients.

The purpose of this study was to assess the clinical oral status and investigate the effect of highly active antiretroviral therapy (HAART) on oral flora diversity in human immunodeficiency virus (HIV)-infected/acquired immune deficiency syndrome (AIDS) patients. We first recorded and analyzed the demographic indicators of 108 HIV-infected patients and assessed their periodontal health, dental health and oral lesion status by oral examination. Besides, we compared the changes in salivary microbial communities of healthy controls, before and after treatment of HAART-processed AIDS patients by Roche 454 sequencing and RT-qPCR. In HIV-infected/AIDS patients, age, sex, marital status, income level, smoking and oral health behaviors had an effect on periodontal clinical indicators; age and marital status were correlated with dental clinical indicators; most of them were accompanied by oral manifestations, mainly including candidiasis albicans, salivary gland disease, AIDS-associated periodontitis, and oral ulcers. Besides, a total of 487 species were detected in the saliva of AIDS patients. The microbial communities of HAART-unprocessed AIDS patients significantly differed from those processed patients, with 112 unique microbial species. More importantly, a large number of conditioned pathogens were also detected in the saliva samples of AIDS patients, which may be associated with opportunistic infections. Therefore, HAART might have a crucial role in salivary microecological balance in AIDS patients. And these patients should pay attention to the maintenance of oral health, and the early initiation of HAART may be important for the development of oral lesions.

Sequencing of E. coli strain UTI89 on multiple sequencing platforms.

OBJECTIVES: The availability of matched sequencing data for the same sample across different sequencing platforms is a necessity for validation and effective comparison of sequencing platforms. A commonly sequenced sample is the lab-adapted MG1655 strain of Escherichia coli; however, this strain is not fully representative of more complex and dynamic genomes of pathogenic E. coli strains. DATA DESCRIPTION: We present six new sequencing data sets for another E. coli strain, UTI89, which is an extraintestinal pathogenic strain isolated from a patient suffering from a urinary tract infection. We now provide matched whole genome sequencing data generated using the PacBio RSII, Oxford Nanopore MinION R9.4, Ion Torrent, ABI SOLiD, and Illumina NextSeq sequencers. Together with other publically available datasets, UTI89 has a nearly complete suite of data generated on most second- and third-generation sequencers. These data can be used as an additional validation set for new sequencing technologies and analytical methods. More than being another E. coli strain, however, UTI89 is pathogenic, with a 10% larger genome, additional pathogenicity islands, and a large plasmid, features that are common among other naturally occurring and disease-causing E. coli isolates. These data therefore provide a more medically relevant test set for development of algorithms.

Apolipoprotein M gene single nucleotide polymorphisms discovery in patients with chronic obstructive pulmonary disease and determined by the base-quenched probe technique.

BACKGROUND: It has been reported increased serum apoM levels seen in the patients suffered from obstructive sleep apnoea and chronic obstructive pulmonary diseases. In the present, we further examine the prevalence of apoM gene SNPs in COPD patients. And a new method base-quenched probe technique is established. METHODS: In the present study, we first used the Roche 454 GS Junior high-throughput sequencer to analyze 6 COPD samples and 6 control samples, in these samples we found 3 interesting SNPs (rs805264, rs707922 and rs707921) and then we designed primers and probes to establish a simple and quick screening method that is a base-quench probe technique and the genotype was confirmed by melting curves. With this new technique, we further determined 252 COPD samples and 248 normal subjects were applied as controls. RESULTS: A total of 19 high-confidence mutations were detected in the Roche 454, 6 mutations among them were not been reported in NCBI, but the mutation frequency was <20%. Four mutations occurred only in COPD patients, rs751064723 is located in the first exon of transcript 1 and the rest are located in either apoM gene promoter or intron region. The results of melting curve showed that the wild-type and homozygous mutants of rs805264, rs707922 and rs707921 presented melting valley at two different melting temperatures, and the results were consistent with those of DNA sequencing (K=1, P=0.000). CONCLUSIONS: The detection of apoM gene SNPs laid the foundation for the study of the relationship between COPD and apoM, and the base-quenched probe technique is simple, economic and accurate, and it is suitable for a large number of apoM genotyping studies.

Low-depth shotgun sequencing resolves complete mitochondrial genome sequence of Labeo rohita.

Labeo rohita, popularly known as rohu, is a widely cultured species in whole Indian subcontinent. In the present study, we used in-silico approach to resolve complete mitochondrial genome of rohu. Low-depth shotgun sequencing using Roche 454 GS FLX (Branford, Connecticut, USA) followed by de novo assembly in CLC Genomics Workbench version 7.0.4 (Aarhus, Denmark) revealed the complete mitogenome of L. rohita to be 16 606 bp long (accession No. KR185963). It comprised of 13 protein-coding genes, 22 tRNAs, 2 rRNAs and 1 putative control region. The gene order and organization are similar to most vertebrates. The mitogenome in the present investigation has 99% similarity with that of previously reported mitogenomes of rohu and this is also evident from the phylogenetic study using maximum-likelihood (ML) tree method. This study was done to determine the feasibility, accuracy and reliability of low-depth sequence data obtained from NGS platform as compared to the Sanger sequencing. Thus, NGS technology has proven to be competent and a rapid in-silico alternative to resolve the complete mitochondrial genome sequence, thereby reducing labors and time.

Chloroplast Genome Sequence of Pigeonpea (Cajanus cajan (L.) Millspaugh) and Cajanus scarabaeoides (L.) Thouars: Genome Organization and Comparison with Other Legumes.

Pigeonpea (Cajanus cajan (L.) Millspaugh), a diploid (2n = 22) legume crop with a genome size of 852 Mbp, serves as an important source of human dietary protein especially in South East Asian and African regions. In this study, the draft chloroplast genomes of Cajanus cajan and Cajanus scarabaeoides (L.) Thouars were generated. Cajanus scarabaeoides is an important species of the Cajanus gene pool and has also been used for developing promising CMS system by different groups. A male sterile genotype harboring the C. scarabaeoides cytoplasm was used for sequencing the plastid genome. The cp genome of C. cajan is 152,242bp long, having a quadripartite structure with LSC of 83,455 bp and SSC of 17,871 bp separated by IRs of 25,398 bp. Similarly, the cp genome of C. scarabaeoides is 152,201bp long, having a quadripartite structure in which IRs of 25,402 bp length separates 83,423 bp of LSC and 17,854 bp of SSC. The pigeonpea cp genome contains 116 unique genes, including 30 tRNA, 4 rRNA, 78 predicted protein coding genes and 5 pseudogenes. A 50 kb inversion was observed in the LSC region of pigeonpea cp genome, consistent with other legumes. Comparison of cp genome with other legumes revealed the contraction of IR boundaries due to the absence of rps19 gene in the IR region. Chloroplast SSRs were mined and a total of 280 and 292 cpSSRs were identified in C. scarabaeoides and C. cajan respectively. RNA editing was observed at 37 sites in both C. scarabaeoides and C. cajan, with maximum occurrence in the ndh genes. The pigeonpea cp genome sequence would be beneficial in providing informative molecular markers which can be utilized for genetic diversity analysis and aid in understanding the plant systematics studies among major grain legumes.

Whole Genome Comparison of Thermus sp. NMX2.A1 Reveals Principle Carbon Metabolism Differences with Closest Relation Thermus scotoductus SA-01.

Genome sequencing of the yellow-pigmented, thermophilic bacterium Thermus sp. NMX2.A1 resulted in a 2.29 Mb draft genome that encodes for 2312 proteins. The genetic relationship between various strains from the genus Thermus was assessed based on phylogenomic analyses using a concatenated set of conserved proteins. The resulting phylogenetic tree illustrated that Thermus sp. NMX2 A.1 clusters together with Thermus scotoductus SA-01, despite being isolated from vastly different geographical locations. The close evolutionary relationship and metabolic parallels between the two strains has previously been recognized; however, neither strain's genome data were available at that point in time. Genomic comparison of the Thermus sp. NMX2.A1 and T. scotoductus SA-01, as well as other closely related Thermus strains, revealed a high degree of synteny at both the genomic and proteomic level, with processes such as denitrification and natural cell competence appearing to be conserved. However, despite this high level of similarity, analysis revealed a complete, putative Calvin-Benson-Bassham (CBB) cycle in NMX2.A1 that is absent in SA-01. Analysis of horizontally transferred gene islands provide evidence that NMX2 selected these genes due to pressure from its HCO3 (-) rich environment, which is in stark contrast to that of the deep subsurface isolated SA-01.

Transcriptomic SNP discovery for custom genotyping arrays: impacts of sequence data, SNP calling method and genotyping technology on the probability of validation success.

BACKGROUND: Single nucleotide polymorphism (SNP) discovery is an important goal of many studies. However, the number of 'putative' SNPs discovered from a sequence resource may not provide a reliable indication of the number that will successfully validate with a given genotyping technology. For this it may be necessary to account for factors such as the method used for SNP discovery and the type of sequence data from which it originates, suitability of the SNP flanking sequences for probe design, and genomic context. To explore the relative importance of these and other factors, we used Illumina sequencing to augment an existing Roche 454 transcriptome assembly for the Antarctic fur seal (Arctocephalus gazella). We then mapped the raw Illumina reads to the new hybrid transcriptome using BWA and BOWTIE2 before calling SNPs with GATK. The resulting markers were pooled with two existing sets of SNPs called from the original 454 assembly using NEWBLER and SWAP454. Finally, we explored the extent to which SNPs discovered using these four methods overlapped and predicted the corresponding validation outcomes for both Illumina Infinium iSelect HD and Affymetrix Axiom arrays. RESULTS: Collating markers across all discovery methods resulted in a global list of 34,718 SNPs. However, concordance between the methods was surprisingly poor, with only 51.0 % of SNPs being discovered by more than one method and 13.5 % being called from both the 454 and Illumina datasets. Using a predictive modeling approach, we could also show that SNPs called from the Illumina data were on average more likely to successfully validate, as were SNPs called by more than one method. Above and beyond this pattern, predicted validation outcomes were also consistently better for Affymetrix Axiom arrays. CONCLUSIONS: Our results suggest that focusing on SNPs called by more than one method could potentially improve validation outcomes. They also highlight possible differences between alternative genotyping technologies that could be explored in future studies of non-model organisms.

Precise and efficient antibody epitope determination through library design, yeast display and next-generation sequencing.

The ability of antibodies to bind an antigen with a high degree of affinity and specificity has led them to become the largest and fastest growing class of therapeutic proteins. Clearly identifying the epitope at which they bind their cognate antigen provides insight into their mechanism of action and helps differentiate antibodies that bind the same antigen. Here, we describe a method to precisely and efficiently map the epitopes of a panel of antibodies in parallel over the course of several weeks. This method relies on the combination of rational library design, quantitative yeast surface display and next-generation DNA sequencing and was demonstrated by mapping the epitopes of several antibodies that neutralize alpha toxin from Staphylococcus aureus. The accuracy of this method was confirmed by comparing the results to the co-crystal structure of one antibody and alpha toxin and was further refined by the inclusion of a lower-affinity variant of the antibody. In addition, this method produced quantitative insight into the epitope residues most critical for the antibody-antigen interaction and enabled the relative affinities of each antibody toward alpha toxin variants to be estimated. This affinity estimate serves as a predictor of neutralizing antibody potency and was used to anticipate the ability of each antibody to effectively bind and neutralize naturally occurring alpha toxin variants secreted by strains of S. aureus, including clinically relevant strains. Ultimately this type information can be used to help select the best clinical candidate among a set of antibodies against a given antigen.

The antibody genetics of multiple sclerosis: comparing next-generation sequencing to sanger sequencing.

We previously identified a distinct mutation pattern in the antibody genes of B cells isolated from cerebrospinal fluid (CSF) that can identify patients who have relapsing-remitting multiple sclerosis (RRMS) and patients with clinically isolated syndromes who will convert to RRMS. This antibody gene signature (AGS) was developed using Sanger sequencing of single B cells. While potentially helpful to patients, Sanger sequencing is not an assay that can be practically deployed in clinical settings. In order to provide AGS evaluations to patients as part of their diagnostic workup, we developed protocols to generate AGS scores using next-generation DNA sequencing (NGS) on CSF-derived cell pellets without the need to isolate single cells. This approach has the potential to increase the coverage of the B-cell population being analyzed, reduce the time needed to generate AGS scores, and may improve the overall performance of the AGS approach as a diagnostic test in the future. However, no investigations have focused on whether NGS-based repertoires will properly reflect antibody gene frequencies and somatic hypermutation patterns defined by Sanger sequencing. To address this issue, we isolated paired CSF samples from eight patients who either had MS or were at risk to develop MS. Here, we present data that antibody gene frequencies and somatic hypermutation patterns are similar in Sanger and NGS-based antibody repertoires from these paired CSF samples. In addition, AGS scores derived from the NGS database correctly identified the patients who initially had or subsequently converted to RRMS, with precision similar to that of the Sanger sequencing approach. Further investigation of the utility of the AGS in predicting conversion to MS using NGS-derived antibody repertoires in a larger cohort of patients is warranted.

A fast and robust method for full genome sequencing of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Type 1 and Type 2.

PRRSV is a positive-sense RNA virus with a high degree of genetic variability among isolates. For diagnostic sensitivity and vaccine design it is essential to monitor PRRSV genetic diversity. However, to date only a few full genome sequences of PRRSV isolates have been made publicly available. In the present study, fast and robust methods for long range RT-PCR amplification and subsequent next generation sequencing (NGS) were developed and validated on nine Type 1 and nine Type 2 PRRSV viruses. The methods generated robust and reliable sequences both on primary material and cell culture adapted viruses and the protocols performed well on all three NGS platforms tested (Roche 454 FLX, Illumina HiSeq2000, and Ion Torrent PGM™ Sequencer). These methods will greatly facilitate the generation of more full genome PRRSV sequences globally.

Explanatory chapter: next generation sequencing.

Technological breakthroughs in sequencing technologies have driven the advancement of molecular biology and molecular genetics research. The advent of high-throughput Sanger sequencing (for information on the method, see Sanger Dideoxy Sequencing of DNA) in the mid- to late-1990s made possible the accelerated completion of the human genome project, which has since revolutionized the pace of discovery in biomedical research. Similarly, the advent of next generation sequencing is poised to revolutionize biomedical research and usher a new era of individualized, rational medicine. The term next generation sequencing refers to technologies that have enabled the massively parallel analysis of DNA sequence facilitated through the convergence of advancements in molecular biology, nucleic acid chemistry and biochemistry, computational biology, and electrical and mechanical engineering. The current next generation sequencing technologies are capable of sequencing tens to hundreds of millions of DNA templates simultaneously and generate >4 gigabases of sequence in a single day. These technologies have largely started to replace high-throughput Sanger sequencing for large-scale genomic projects, and have created significant enthusiasm for the advent of a new era of individualized medicine.