Search

PURPOSE Tumor mutational burden (TMB) is a developing biomarker in non–small-cell lung cancer (NSCLC). Little is known regarding differences between TMB and sample location, histology, or other biomarkers. METHODS A total of 3,424 unmatched NSCLC samples, including 2,351 lung adenocarcinomas (LUADs) and 1,073 lung squamous cell carcinomas (LUSCs), underwent profiling, including next-generation sequencing of 592 cancer-related genes, programmed death ligand 1 immunohistochemistry, and TMB. The rate TMB of 10 mutations per megabase (Mb) or greater was compared between primary and metastatic LUAD and LUSC. Molecular alteration frequency was compared at a cutoff of 10 mutations/Mb. RESULTS LUAD metastases were more likely to have a TMB of 10 mutations/Mb or greater compared with primary LUADs (38% v 25%; P < .001), and this difference was most pronounced with brain metastases (61% v 35% for other metastases; P < .001). The median TMB for LUAD brain metastases was 13 mutations/Mb compared with six mutations/Mb for primary LUADs. Variability existed for other LUAD metastasis sites, with adrenal metastases most likely to meet the cutoff of 10 mutations/Mb (51%) and bone metastases least likely to meet the cutoff (19%). TMB was more commonly 10 mutations/Mb or greater for LUSC primary tumors than for LUAD primary tumors (35% v 25%, respectively; P < .001). LUSC metastases were more likely to have a TMB of 10 mutations/Mb or greater than LUSC primary tumors. Poorly differentiated disease was more likely have a TMB of 10 mutations/Mb or greater when stratified by histology and primary tumor or metastasis. Site-specific molecular differences existed at this TMB cutoff including programmed death ligand 1 positivity and STK11 and KRAS mutation rate. CONCLUSION TMB is a site-specific biomarker in NSCLC with important spatial and histologic differences. TMB is more frequently 10 mutations/Mb or greater in LUAD and LUSC metastases and highest in LUAD brain metastases. Along this TMB cutoff, clinically informative distinctions exist in other tumor profiling characteristics. Further investigation is needed to expand on these findings.

BackgroundTumor mutational burden (TMB), defined as the number of somatic mutations per megabase of interrogated genomic sequence, demonstrates predictive biomarker potential for the identification of patients with cancer most likely to respond to immune checkpoint inhibitors. TMB is optimally calculated by whole exome sequencing (WES), but next-generation sequencing targeted panels provide TMB estimates in a time-effective and cost-effective manner. However, differences in panel size and gene coverage, in addition to the underlying bioinformatics pipelines, are known drivers of variability in TMB estimates across laboratories. By directly comparing panel-based TMB estimates from participating laboratories, this study aims to characterize the theoretical variability of panel-based TMB estimates, and provides guidelines on TMB reporting, analytic validation requirements and reference standard alignment in order to maintain consistency of TMB estimation across platforms.MethodsEleven laboratories used WES data from The Cancer Genome Atlas Multi-Center Mutation calling in Multiple Cancers (MC3) samples and calculated TMB from the subset of the exome restricted to the genes covered by their targeted panel using their own bioinformatics pipeline (panel TMB). A reference TMB value was calculated from the entire exome using a uniform bioinformatics pipeline all members agreed on (WES TMB). Linear regression analyses were performed to investigate the relationship between WES and panel TMB for all 32 cancer types combined and separately. Variability in panel TMB values at various WES TMB values was also quantified using 95% prediction limits.ResultsStudy results demonstrated that variability within and between panel TMB values increases as the WES TMB values increase. For each panel, prediction limits based on linear regression analyses that modeled panel TMB as a function of WES TMB were calculated and found to approximately capture the intended 95% of observed panel TMB values. Certain cancer types, such as uterine, bladder and colon cancers exhibited greater variability in panel TMB values, compared with lung and head and neck cancers.ConclusionsIncreasing uptake of TMB as a predictive biomarker in the clinic creates an urgent need to bring stakeholders together to agree on the harmonization of key aspects of panel-based TMB estimation, such as the standardization of TMB reporting, standardization of analytical validation studies and the alignment of panel-based TMB values with a reference standard. These harmonization efforts should improve consistency and reliability of panel TMB estimates and aid in clinical decision-making.

Microsatellite instability (MSI) is a key secondary effect of a defective DNA mismatch repair mechanism resulting in incorrectly replicated microsatellites in many malignant tumors. Historically, MSI detection has been performed by fragment analysis (FA) on a panel of representative genomic markers. More recently, using next‐generation sequencing (NGS) to analyze thousands of microsatellites has been shown to improve the robustness and sensitivity of MSI detection. However, NGS‐based MSI tests can be prone to population biases if NGS results are aligned to a reference genome instead of patient‐matched normal tissue. We observed an increased rate of false positives in patients of African ancestry with an NGS‐based diagnostic for MSI status utilizing 7317 microsatellite loci. We then minimized this bias by training a modified calling model that utilized 2011 microsatellite loci. With these adjustments 100% (95% CI: 89.1% to 100%) of African ancestry patients in an independent validation test were called correctly using the updated model. This poses not only a significant technical improvement but also has an important clinical impact on directing immune checkpoint inhibitor therapy., This work presents observations of a significant population bias in an NGS‐based MSI test and minimizes this bias using a computational‐based approach. The research poses not only a significant technical improvement to the NGS‐based MSI diagnostic, but has an important clinical impact on directing immune‐check point inhibitor therapy for late‐stage cancer patients of non‐European descent.

The Hawaiian archipelago provides a natural arena for understanding adaptive radiation and speciation. The Hawaiian Drosophila are one of the most diverse endemic groups in Hawaiì with up to 1,000 species. We sequenced and analyzed entire genomes of recently diverged species of Hawaiian picture-winged Drosophila, Drosophila silvestris and Drosophila heteroneura from Hawaiì Island, in comparison with Drosophila planitibia, their sister species from Maui, a neighboring island where a common ancestor of all three had likely occurred. Genome-wide single nucleotide polymorphism patterns suggest the more recent origin of D. silvestris and D. heteroneura, as well as a pervasive influence of positive selection on divergence of the three species, with the signatures of positive selection more prominent in sympatry than allopatry. Positively selected genes were significantly enriched for functional terms related to sensory detection and mating, suggesting that sexual selection played an important role in speciation of these species. In particular, sequence variation in Olfactory receptor and Gustatory receptor genes seems to play a major role in adaptive radiation in Hawaiian pictured-winged Drosophila.

A singular genome used for inference into population-based studies is a standard method in genomics. Recent studies show that spontaneous genomic variants can propagate into new generations and these changes can contribute to individual cell aging with environmental and evolutionary elements contributing to cumulative genomic variation. However, the contribution of aging to genomic changes in tissue samples remains uncharacterized. Here, we re port the impact of aging on individual human exomes and their implications. We found the human genome to be dynamic, acquiring a varying number of mutations with age (5,000 to 50,000 in 9 to 16 years). This equates to a variation rate of 9.6x10 -7 to 8.4x10 -6 bp -1 year -1 for nonsynonymous single nucleotide variants and 2.0x10 -4 to 1.0x10 -3 locus -1 year -1 for microsatellite loci in these individuals. These mutations span across 3,000 to 13,000 genes, which commonly showed association with Wnt signaling and Gonadotropin releasing hormone receptor pathways, and indicated for individuals a specific and significant enrichment for increased risk for diabetes, kidney failure, cancer, Rheumatoid arthritis, and Alzheimer's disease- conditions usually associated with aging. The results suggest that "age" is an important variable while analyzing an individual human genome to extract individual-specific clinically significant information necessary for personalized genomics.

// Jasmin H. Bavarva 1 , Hongseok Tae 1 , Lauren McIver 1 and Harold R. Garner 1 1 Virginia Bioinformatics Institute, Virginia Tech, Blacksburg, VA, USA Correspondence: Harold R. Garner, email: // Correspondence: Jasmin H. Bavarva, email: // Keywords : Nicotine; Exome sequencing; MUC4; Biomarker; Mutation targets Received : April 14, 2014 Accepted : May 27, 2014 Published : May 28, 2014 Abstract Although the connection between cancer and cigarette smoke is well established, nicotine is not characterized as a carcinogen. Here, we used exome sequencing to identify nicotine and oxidative stress-induced somatic mutations in normal human epithelial cells and its correlation with cancer. We identified over 6,400 SNVs, indels and microsatellites in each of the stress exposed cells relative to the control, of which, 2,159 were consistently observed at all nicotine doses. These included 429 nsSNVs including 158 novel and 79 cancer-associated. Over 80% of consistently nicotine induced variants overlap with variations detected in oxidative stressed cells, indicating that nicotine induced genomic alterations could be mediated through oxidative stress. Nicotine induced mutations were distributed across 1,585 genes, of which 49% were associated with cancer. MUC family genes were among the top mutated genes. Analysis of 591 lung carcinoma tumor exomes from The Cancer Genome Atlas (TCGA) revealed that 20% of non-small-cell lung cancer tumors in smokers have mutations in at least one of the MUC4, MUC6 or MUC12 genes in contrast to only 6% in non-smokers. These results indicate that nicotine induces genomic variations, promotes instability potentially mediated by oxidative stress, implicating nicotine in carcinogenesis, and establishes MUC genes as potential targets.

Motivation: Simple tandem repeats are highly variable genetic elements and widespread in genomes of many organisms. Next-generation sequencing technologies have enabled a robust comparison of large numbers of simple tandem repeat loci; however, analysis of their variation using traditional sequence analysis approaches still remains limiting and problematic due to variants occurring in repeat sequences confusing alignment programs into mapping sequence reads to incorrect loci when the sequence reads are significantly different from the reference sequence. Results: We have developed a program, ReviSTER, which is an automated pipeline using a ‘local mapping reference reconstruction method’ to revise mismapped or partially misaligned reads at simple tandem repeat loci. RevisSTER estimates alleles of repeat loci using a local alignment method and creates temporary local mapping reference sequences, and finally remaps reads to the local mapping references. Using this approach, ReviSTER was able to successfully revise reads misaligned to repeat loci from both simulated data and real data. Availability: ReviSTER is open-source software available at http://revister.sourceforge.net. Contact: garner@vbi.vt.edu Supplementary information: Supplementary data are available at Bioinformatics online.

The human genome is 99% complete. This study contributes to filling the 1% gap by enriching previously unknown repeat regions called microsatellites (MST). We devised a Global MST Enrichment (GME) kit to enrich and nextgen sequence 2 colorectal cell lines and 16 normal human samples to illustrate its utility in identifying contigs from reads that do not map to the genome reference. The analysis of these samples yielded 790 novel extra-referential concordant contigs that are observed in more than one sample. We searched for evidence of functional elements in the concordant contigs in two ways: (1) BLAST-ing each contig against normal RNA-Seq samples, (2) Checking for predicted functional elements using GlimmerHMM. Of the 790 concordant contigs, 37 had an exact match to at least one RNA-Seq read; 15 aligned to more than 100 RNA-Seq reads. Of the 249 concordant contigs predicted by GlimmerHMM to have functional elements, 6 had at least one exact RNA-Seq match. BLAST-ing these novel contigs against all publically available sequences confirmed that they were found in human and chimpanzee BAC and FOSMID clones sequenced as part of the original human genome project. These extra-referential contigs predominantly contained pentameric repeats, especially two motifs: AATGG and GTGGA.

Sequencing data analysis remains limiting and problematic, especially for low complexity repeat sequences and transposon elements due to inherent sequencing errors and short sequence read lengths. We have developed a program, ReviSeq, which uses a hybrid method composed of iterative remapping and local assembly upon a bacterial sequence backbone. Application of this method to six Brucella suis field isolates compared to the newly revised B. suis 1330 reference genome identified on average 13, 15, 19 and 9 more variants per sample than STAMPY/SAMtools, BWA/SAMtools, iCORN and BWA/PINDEL pipelines, and excluded on average 4, 2, 3 and 19 variants per sample, respectively. In total, using this iterative approach, we identified on average 87 variants including SNVs, short INDELs and long INDELs per strain when compared to the reference. Our program outperforms other methods especially for long INDEL calling. The program is available at http://reviseq.sourceforge.net .

Several studies have demonstrated that unmapped reads in next generation sequencing data could be used to identify infectious agents or structural variants, but there has been no intensive effort to analyze and classify all non-human sequences found in individual large data sets. To identify commonality in non-human sequences by infectious agents and putative contamination events, we analyzed non-human sequences in 150 genomic sequencing data files from the 1000 Genomes Project and observed that 0.13% of reads on average showed similarities to non-human genomes. We compared results among different sample groups divided based on ethnicities, sequencing centers and enrichment methods (whole genome sequencing vs. exome sequencing) and found that sequencing centers had specific signatures of contaminating genomes as ‘time stamps’. We also observed many unmapped reads that falsely indicated contamination because of the high similarity of human sequences to sequences in non-human genome assemblies such as mouse and Nicotiana.

Although abundant biology data have been accumulated in public biology databases, such as GenBank and PIR, few easy-interface services are provided for users to ac-cess or update them. We have developed a system, named BioStore, that is composed of several programs to aid users to not only access public data but also share their own data easily. The service can be used for maintaining a local database as a repository of raw data files of several public databases and distributing the da-ta files to other users. Currently, BioStore manipulates major bio-databases and will expand to include more databases and more useful interfaces.Availability: This system is available from http://www. ccbb.re.kr/biostoreKeywords: BioStore, FTP, GenBank, repository Introduction Biology data, especially DNA and protein sequence da-ta, are very rapidly increasing. Most data have been submitted to major bioinformatics databases, such as GenBank (Benson et al., 1994), PIR (McGarvey et al., 2000), and PDB (Berman, 2000), which are distributed via FTP services to a lot of bioinformatics sites and users. As more and more bio-databases (Jeon et al., 2005; Kim et al., 2008) are constructed and a large amount of data files is kept in each database, most bi-ology researchers have difficulty in searching for in-formation and in sharing their own data with others. CCBB (Center for Computational Biology and Bioin-formatics) of the Korea Institute of Science and Technol-ogy Information has been providing an FTP server and analysis services for major bio-databases. We have de-veloped BioStore as a repository system by implement-ing a few programs to register user data, including pub-lic bio-data at CCBB. Researchers can easily register their own data with BioStore through the system to share them with others.

Many bioinformatics sites have managed local bio-databases, including major databases such as GenBank and PIR with update load. We have developed several programs to monitor the update status of these databases and to FTP them automatically. These programs can be used for maintaining local bio-databases as recent versions and providing up-to-date databases through FTP sites. Currently, the program serves major bio-databases and will extend to accommodate many more bio-databases.

Microsatellites are composed of tandemly repeated short motifs of 1–6 nucleotides. They are common in eukaryotic genomes; in humans they make up as much as 3% of the genome (Lander et al., 2001). These repeat-containing loci tend to be hypervariable with variation occurring among individuals of the same species as well as between species. The origin and evolution of microsatellites remain a major puzzle. The birth of a tetranucleotide repeat (ATGT) in the lineage leading to African apes (gorilla, bonobo and chimpanzee) and humans was documented in the n-globin pseudogene (Messier et al., 1996). To test whether the locus is under further expansion or degeneration or if there are any variations as proposed in the life-cycle of the microsatellites (Buschiazzo and Gemmell, 2006), we analyzed positions chr11:5263801:5263831 (hg19) in 82 samples from the 1000 Genomes Project (Genomes Project et al., 2010). To obtain reliable read coverage for the locus, we selected Illumina exome sequencing data and visually inspected the INDEL length in the reads aligned to the locus. We found a repeat deletion (4 bases) of the 4-mer microsatellite locus in 46 samples (56%) relative to reference allele, as supported by at least two reads in each sample. We assume that the longer (reference) microsatellite allele is ancestral, as at least two mutations are reported in this allele (rs34312249 and rs147740082), suggesting a degeneration process of the pure microsatellite locus. The short allele had 28% frequency in analyzed samples. The genus Homo is believed to have evolved 2.0 Myr (Curnoe, 2010) and modern humans 0.2 Myr (Vigilant et al., 1991) ago. Last evolutionary change in the n-globin pseudogene locus is documented between chimpanzees and humans when the microsatellite locus was expanded in humans. Unfortunately, intraspecies variation in the locus is unknown in chimpanzees and gorillas and we were unable to find the locus in the Neanderthal genome, presumably due to its incompleteness (Green et al., 2006). The locus alteration from orangutan to gorilla took ~7.0 Myr followed by species evolution with no change in the locus (gorilla to chimpanzee). Similarly, the evolution of the chimpanzee to the genus Homo took ~3.0 Myr with an accompanying n-globin microsatellite locus expansion. As time required for subsequent stage of life cycle seems to decrease by about 50% (Birth: ~7 Myr, Expansion: ~3 Myr, Degeneration: 1.8 Myr) (Figure ​(Figure1A),1A), it appears that the plausible death of the n-globin pseudogene microsatellite locus is expected within the next ~0.9 Myr (Figure ​(Figure1B).1B). However, human demography and unknown fitness effects of the microsatellite polymorphism can influence the locus evolutionary dynamics. Figure 1 Evolution and life cycle of an n-globin pseudogene microsatellite locus. (A) Evolutionary time estimate for the species emergence and corresponding change in the n-globin pseudogene microsatellite locus, illustrating different stages of life cycle and ... To test whether the scenario (life cycle of microsatellites) can be applied to other microsatellite loci, we reanalyzed our microsatellite data from multiple species (Galindo et al., 2009). We observed the birth, expansion, and degeneration of several microsatellite loci in four microsatellite motif families in human, chimpanzee, orangutan, rhesus macaque, marmoset and mouse genomes. Although, there were few loci that showed deviations from the expected life-cycle, such as AATGG repeat in DMD and HEPH gene that showed stability for an extended period of time, and in PCNXL2 where it showed one repeat diminution in chimpanzee compared to orangutan and human, majority of the repeats arguably followed the proposed life cycle. Overall, it would be of great scientific significance to constantly monitor a number of marker microsatellite loci for genetic deviations as time progresses and more samples are being sequenced.

Motivation: Inferring lengths of inherited microsatellite alleles with single base pair resolution from short sequence reads is challenging due to several sources of noise caused by the repetitive nature of microsatellites and the technologies used to generate raw sequence data. Results: We have developed a program, GenoTan, using a discretized Gaussian mixture model combined with a rules-based approach to identify inherited variation of microsatellite loci from short sequence reads without paired-end information. It effectively distinguishes length variants from noise including insertion/deletion errors in homopolymer runs by addressing the bidirectional aspect of insertion and deletion errors in sequence reads. Here we first introduce a homopolymer decomposition method which estimates error bias toward insertion or deletion in homopolymer sequence runs. Combining these approaches, GenoTan was able to genotype 94.9% of microsatellite loci accurately from simulated data with 40x sequence coverage quickly while the other programs showed Availability: GenoTan is open-source software available at http://genotan.sourceforge.net. Contact: garner@vbi.vt.edu Supplementary Information: Supplementary data are available at Bioinformatics online

Corrigendum: Draft genome sequence and genetic transformation of the oleaginous alga Nannochloropsis gaditana

Nicotine is a known risk factor for cancer development and has been shown to alter gene expression in cells and tissue upon exposure. We used Illumina® Next Generation Sequencing (NGS) technology to gain unbiased biological insight into the transcriptome of normal epithelial cells (MCF-10A) to nicotine exposure. We generated expression data from 54,699 transcripts using triplicates of control and nicotine stressed cells. As a result, we identified 138 differentially expressed transcripts, including 39 uncharacterized genes. Additionally, 173 transcripts that are primarily associated with DNA replication, recombination, and repair showed evidence for alternative splicing. We discovered the greatest nicotine stress response by HPCAL4 (up-regulated by 4.71 fold) and NPAS3 (down-regulated by -2.73 fold); both are genes that have not been previously implicated in nicotine exposure but are linked to cancer. We also discovered significant down-regulation (-2.3 fold) and alternative splicing of NEAT1 (lncRNA) that may have an important, yet undiscovered regulatory role. Gene ontology analysis revealed nicotine exposure influenced genes involved in cellular and metabolic processes. This study reveals previously unknown consequences of nicotine stress on the transcriptome of normal breast epithelial cells and provides insight into the underlying biological influence of nicotine on normal cells, marking the foundation for future studies.

Background With the advent of next-generation sequencing (NGS) technologies, full cDNA shotgun sequencing has become a major approach in the study of transcriptomes, and several different protocols in 454 sequencing have been invented. As each protocol uses its own short DNA tags or adapters attached to the ends of cDNA fragments for labeling or sequencing, different contaminants may lead to mis-assembly and inaccurate sequence products. Results We have designed and implemented a new program for raw sequence cleaning in a graphical user interface and a batch script. The cleaning process consists of several modules including barcode trimming, sequencing adapter trimming, amplification primer trimming, poly-A tail trimming, vector screening and low quality region trimming. These modules can be combined based on various sequencing applications. Conclusions ESTclean is a software package not only for cleaning cDNA sequences, but also for helping to develop sequencing protocols by providing summary tables and figures for sequencing quality control in a graphical user interface. It outperforms in cleaning read sequences from complicated sequencing protocols which use barcodes and multiple amplification primers.

The detection and identification of bio-threat agents and the study of host-pathogen interactions require a high-resolution detection platform capable of discerning closely related species. Diverse analysis methods are used to identify pathogens, specifically Brucella species or biovars. In this study, we compared four diagnostic approaches including serology-based biochemical test, PCR assay, microarray analysis using a Universal Bio-signature Detection Array (UBDA) and whole genome “deep” sequencing for Brucella organisms including a number of field isolates. We found that although there was frequent agreement among the different tests, some tests gave compound/contradictory results that were a consequence of species diversity due to mixed infections or minor contaminants as measured by UBDA and validated from whole genome sequence. By comparing these analysis techniques, we demonstrate that standard diagnostics used in the field are limited in their ability to identify genomic DNA contaminants in field isolates while UBDA and sequencing analysis are highly sensitive in tracing genomic differences among the isolates.

We conducted a genome-wide DNA methylation analysis in CD19 (+) B-cells from chronic lymphocytic leukemia (CLL) patients and normal control samples using reduced representation bisulfite sequencing (RRBS). The methylation status of 1.8-2.3 million CpGs in the CLL genome was determined; about 45% of these CpGs were located in more than 23,000 CpG islands (CGIs). While global CpG methylation was similar between CLL and normal B-cells, 1764 gene promoters were identified as being differentially methylated in at least one CLL sample when compared with normal B-cell samples. Nineteen percent of the differentially methylated genes were involved in transcriptional regulation. Aberrant hypermethylation was found in all HOX gene clusters and a significant number of WNT signaling pathway genes. Hypomethylation occurred more frequently in the gene body including introns, exons, and 3'-UTRs in CLL. The NFATc1 P2 promoter and first intron was found to be hypomethylated and correlated with upregulation of both NFATc1 RNA and protein expression levels in CLL suggesting that an epigenetic mechanism is involved in the constitutive activation of NFAT activity in CLL cells. This comprehensive DNA methylation analysis will further our understanding of the epigenetic contribution to cellular dysfunction in CLL.

Brucella suis is a causative agent of porcine brucellosis. We report the resequencing of the original sample upon which the published sequence of Brucella suis 1330 is based and describe the differences between the published assembly and our assembly at 12 loci.

The potential use of algae in biofuels applications is receiving significant attention. However, none of the current algal model species are competitive production strains. Here we present a draft genome sequence and a genetic transformation method for the marine microalga Nannochloropsis gaditana CCMP526. We show that N. gaditana has highly favourable lipid yields, and is a promising production organism. The genome assembly includes nuclear (~29 Mb) and organellar genomes, and contains 9,052 gene models. We define the genes required for glycerolipid biogenesis and detail the differential regulation of genes during nitrogen-limited lipid biosynthesis. Phylogenomic analysis identifies genetic attributes of this organism, including unique stramenopile photosynthesis genes and gene expansions that may explain the distinguishing photoautotrophic phenotypes observed. The availability of a genome sequence and transformation methods will facilitate investigations into N. gaditana lipid biosynthesis and permit genetic engineering strategies to further improve this naturally productive alga., Algae show much promise in the production of biofuels owing to their high photoautotrophic biomass and lipid production rates. In this study, the draft genome of Nannochloropsis gaditana CCMP526 and a method for the transformation of this alga are reported, facilitating the investigation of lipid synthesis and biofuel production.

MAPSI (Management and Analysis for Polyketide Synthase Type I) has been developed to offer computational analysis methods to detect type I PKS (polyketide synthase) gene clusters in genome sequences. MAPSI provides a genome analysis component, which detects PKS gene clusters by identifying domains in proteins of a genome. MAPSI also contains databases on polyketides and genome annotation data, as well as analytic components such as new PKS assembly and domain analysis. The polyketide data and analysis component are accessible through Web interfaces and are displayed with diverse information. MAPSI, which was developed to aid researchers studying type I polyketides, provides diverse components to access and analyze polyketide information and should become a very powerful computational tool for polyketide research. The system can be extended through further studies of factors related to the biological activities of polyketides.

Polyketides have diverse biological activities, including pharmacological functions such as antibiotic, antitumor and agrochemical properties. They are biosynthesized from short carboxylic acid precursors by polyketide synthases (PKSs). As natural polyketide products include many clinically important drugs and the volume of data on polyketides is rapidly increasing, the development of a database system to manage polyketide data is essential. MapsiDB is an integrated web database formulated to contain data on type I polyketides and their PKSs, including domain and module composition and related genome information. Data on polyketides were collected from journals and online resources and processed with analysis programs. Web interfaces were utilized to construct and to access this database, allowing polyketide researchers to add their data to this database and to use it easily.

Background Polyketides are secondary metabolites of microorganisms with diverse biological activities, including pharmacological functions such as antibiotic, antitumor and agrochemical properties. Polyketides are synthesized by serialized reactions of a set of enzymes called polyketide synthase(PKS)s, which coordinate the elongation of carbon skeletons by the stepwise condensation of short carbon precursors. Due to their importance as drugs, the volume of data on polyketides is rapidly increasing and creating a need for computational analysis methods for efficient polyketide research. Moreover, the increasing use of genetic engineering to research new kinds of polyketides requires genome wide analysis. Results We describe a system named ASMPKS (Analysis System for Modular Polyketide Synthesis) for computational analysis of PKSs against genome sequences. It also provides overall management of information on modular PKS, including polyketide database construction, new PKS assembly, and chain visualization. ASMPKS operates on a web interface to construct the database and to analyze PKSs, allowing polyketide researchers to add their data to this database and to use it easily. In addition, the ASMPKS can predict functional modules for a protein sequence submitted by users, estimate the chemical composition of a polyketide synthesized from the modules, and display the carbon chain structure on the web interface. Conclusion ASMPKS has powerful computation features to aid modular PKS research. As various factors, such as starter units and post-processing, are related to polyketide biosynthesis, ASMPKS will be improved through further development for study of the factors.

Gene structure prediction, which is to predict protein coding regions in a given nucleotide sequence, is a critical process in annotating genes and greatly affects gene analysis and genome annotation. As the gene structure of eukaryotes is much more complicated than that of prokaryotic genes, eukaryotic gene structure prediction should have more diverse and more complicated computational models. We have developed GeneChaser, a gene structure prediction program, using a duration hidden markov model. GeneChaser consists of two major processes, one of which is to train datasets to produce parameter values and the other of which is to predict protein coding regions based on the parameter values. The program predicts multiple genes rather than a single gene from a DNA sequence. To predict the gene structure for a huge chromosomal DNA sequence, it splits the sequence into overlapped fragments and performs prediction process for each fragment. A few computational models were implemented to detect signal patterns and their scanning efficiency was evaluated. Based on a few criteria, its prediction performance was compared with that of a few commonly used programs, GeneID and Morgan.

Methylation of cytosine is a post-synthesis modification that does not affect the primary DNA sequence but greatly influences gene expression level and phenotypes of an organism. As high-throughput sequencing of bisulfite-treated DNA is the most efficient method to identify methylated sites, several tools to map sequencing reads on a reference are available. But tools to visualize and to interpret the methylation level of methylation sites are currently insufficient. Herein, we present a novel tool to visualize the methylation level of CpG sites.

Microsatellites (MST), tandem repeats of 1–6 nucleotide motifs, are mutational hot-spots with a bias for insertions and deletions (INDELs) rather than single nucleotide polymorphisms (SNPs). The majority of MST instability studies are limited to a small number of loci, the Bethesda markers, which are only informative for a subset of colorectal cancers. In this paper we evaluate non-haplotype alleles present within next-gen sequencing data to evaluate somatic MST variation (SMV) within DNA repair proficient and DNA repair defective cell lines. We confirm that alleles present within next-gen data that do not contribute to the haplotype can be reliably quantified and utilized to evaluate the SMV without requiring comparisons of matched samples. We observed that SMV patterns found in DNA repair proficient cell lines without DNA repair defects, MCF10A, HEK293 and PD20 RV:D2, had consistent patterns among samples. Further, we were able to confirm that changes in SMV patterns in cell lines lacking functional BRCA2, FANCD2 and mismatch repair were consistent with the different pathways perturbed. Using this new exome sequencing analysis approach we show that DNA instability can be identified in a sample and that patterns of instability vary depending on the impaired DNA repair mechanism, and that genes harboring minor alleles are strongly associated with cancer pathways. The MST Minor Allele Caller used for this study is available at https://github.com/zalmanv/MST_minor_allele_caller. Published version

Brucella suis is the causative agent of swine brucellosis and is known to be able to infect several different hosts, including cattle, dogs, and horses, without causing disease symptoms. Here we report the complete genome sequence of Brucella suis VBI22, which was isolated from raw milk from an infected cow.

Background The reptiles, characterized by both diversity and unique evolutionary adaptations, provide a comprehensive system for comparative studies of metabolism, physiology, and development. However, molecular resources for ectothermic reptiles are severely limited, hampering our ability to study the genetic basis for many evolutionarily important traits such as metabolic plasticity, extreme longevity, limblessness, venom, and freeze tolerance. Here we use massively parallel sequencing (454 GS-FLX Titanium) to generate a transcriptome of the western terrestrial garter snake (Thamnophis elegans) with two goals in mind. First, we develop a molecular resource for an ectothermic reptile; and second, we use these sex-specific transcriptomes to identify differences in the presence of expressed transcripts and potential genes of evolutionary interest. Results Using sex-specific pools of RNA (one pool for females, one pool for males) representing 7 tissue types and 35 diverse individuals, we produced 1.24 million sequence reads, which averaged 366 bp in length after cleaning. Assembly of the cleaned reads from both sexes with NEWBLER and MIRA resulted in 96,379 contigs containing 87% of the cleaned reads. Over 34% of these contigs and 13% of the singletons were annotated based on homology to previously identified proteins. From these homology assignments, additional clustering, and ORF predictions, we estimate that this transcriptome contains ~13,000 unique genes that were previously identified in other species and over 66,000 transcripts from unidentified protein-coding genes. Furthermore, we use a graph-clustering method to identify contigs linked by NEWBLER-split reads that represent divergent alleles, gene duplications, and alternatively spliced transcripts. Beyond gene identification, we identified 95,295 SNPs and 31,651 INDELs. From these sex-specific transcriptomes, we identified 190 genes that were only present in the mRNA sequenced from one of the sexes (84 female-specific, 106 male-specific), and many highly variable genes of evolutionary interest. Conclusions This is the first large-scale, multi-organ transcriptome for an ectothermic reptile. This resource provides the most comprehensive set of EST sequences available for an individual ectothermic reptile species, increasing the number of snake ESTs 50-fold. We have identified genes that appear to be under evolutionary selection and those that are sex-specific. This resource will assist studies on gene expression and comparative genomics, and will facilitate the study of evolutionarily important traits at the molecular level.

Background Horned beetles, in particular in the genus Onthophagus, are important models for studies on sexual selection, biological radiations, the origin of novel traits, developmental plasticity, biocontrol, conservation, and forensic biology. Despite their growing prominence as models for studying both basic and applied questions in biology, little genomic or transcriptomic data are available for this genus. We used massively parallel pyrosequencing (Roche 454-FLX platform) to produce a comprehensive EST dataset for the horned beetle Onthophagus taurus. To maximize sequence diversity, we pooled RNA extracted from a normalized library encompassing diverse developmental stages and both sexes. Results We used 454 pyrosequencing to sequence ESTs from all post-embryonic stages of O. taurus. Approximately 1.36 million reads assembled into 50,080 non-redundant sequences encompassing a total of 26.5 Mbp. The non-redundant sequences match over half of the genes in Tribolium castaneum, the most closely related species with a sequenced genome. Analyses of Gene Ontology annotations and biochemical pathways indicate that the O. taurus sequences reflect a wide and representative sampling of biological functions and biochemical processes. An analysis of sequence polymorphisms revealed that SNP frequency was negatively related to overall expression level and the number of tissue types in which a given gene is expressed. The most variable genes were enriched for a limited number of GO annotations whereas the least variable genes were enriched for a wide range of GO terms directly related to fitness. Conclusions This study provides the first large-scale EST database for horned beetles, a much-needed resource for advancing the study of these organisms. Furthermore, we identified instances of gene duplications and alternative splicing, useful for future study of gene regulation, and a large number of SNP markers that could be used in population-genetic studies of O. taurus and possibly other horned beetles.