Unannotated Genomic Sequences Research Articles

Impact of sequencing depth and technology on de novo RNA-Seq assembly

BackgroundRNA-Seq data is inherently nonuniform for different transcripts because of differences in gene expression. This makes it challenging to decide how much data should be generated from each sample. How much should one spend to recover the less expressed transcripts? The sequencing technology used is another consideration, as there are inevitably always biases against certain sequences. To investigate these effects, we first looked at high-depth libraries from a set of well-annotated organisms to ascertain the impact of sequencing depth on de novo assembly. We then looked at libraries sequenced from the Universal Human Reference RNA (UHRR) to compare the performance of Illumina HiSeq and MGI DNBseq™ technologies.ResultsOn the issue of sequencing depth, the amount of exomic sequence assembled plateaued using data sets of approximately 2 to 8 Gbp. However, the amount of genomic sequence assembled did not plateau for many of the analyzed organisms. Most of the unannotated genomic sequences are single-exon transcripts whose biological significance will be questionable for some users. On the issue of sequencing technology, both of the analyzed platforms recovered a similar number of full-length transcripts. The missing “gap” regions in the HiSeq assemblies were often attributed to higher GC contents, but this may be an artefact of library preparation and not of sequencing technology.ConclusionsIncreasing sequencing depth beyond modest data sets of less than 10 Gbp recovers a plethora of single-exon transcripts undocumented in genome annotations. DNBseq™ is a viable alternative to HiSeq for de novo RNA-Seq assembly.

Development of a High-Resolution Multi-Locus Microsatellite Typing Method for Colletotrichum gloeosporioides

Colletotrichum gloeosporioides is an economically important fungal pathogen causing substantial yield losses indifferent host plants. To understand the genetic diversity and molecular epidemiology of this fungus, we have developed a novel, high-resolution multi-locus microsatellite typing (MLMT) method. Bioinformatic analysis of C. gloeosporioides unannotated genome sequence yielded eight potential microsatellite loci, of which five, CG1 (GT)n, CG2 (GT1)n, CG3 (TC)n, CG4 (CT)n, and CG5 (CT1)n were selected for further study based on their universal amplification potential, reproducibility, and repeat number polymorphism. The selected microsatellites were used to analyze 31 strains of C. gloeosporioides isolated from 20 different host plants from India. All microsatellite loci were found to be polymorphic, and the approximate fragment sizes of microsatellite loci CG1, CG2, CG3, CG4, and CG5 were in ranges of 213–241, 197–227, 231–265, 209–275, and 132–188, respectively. Among the 31 isolates, 55 different genotypes were identified. The Simpson’s index of diversity (D) values for the individual locus ranged from 0.79 to 0.92, with the D value of all combined five microsatellite loci being 0.99. Microsatellite data analysis revealed that isolates from Ocimum sanctum, Capsicum annuum (chili pepper), and Mangifera indica (mango) formed distinct clusters, therefore exhibited some level of correlation between certain genotypes and host. The developed MLMT method would be a powerful tool for studying the genetic diversity and any possible genotype-host correlation in C. gloeosporioides.

GI-SVM: A sensitive method for predicting genomic islands based on unannotated sequence of a single genome.

Genomic islands (GIs) are clusters of functionally related genes acquired by lateral genetic transfer (LGT), and they are present in many bacterial genomes. GIs are extremely important for bacterial research, because they not only promote genome evolution but also contain genes that enhance adaption and enable antibiotic resistance. Many methods have been proposed to predict GI. But most of them rely on either annotations or comparisons with other closely related genomes. Hence these methods cannot be easily applied to new genomes. As the number of newly sequenced bacterial genomes rapidly increases, there is a need for methods to detect GI based solely on sequences of a single genome. In this paper, we propose a novel method, GI-SVM, to predict GIs given only the unannotated genome sequence. GI-SVM is based on one-class support vector machine (SVM), utilizing composition bias in terms of k-mer content. From our evaluations on three real genomes, GI-SVM can achieve higher recall compared with current methods, without much loss of precision. Besides, GI-SVM allows flexible parameter tuning to get optimal results for each genome. In short, GI-SVM provides a more sensitive method for researchers interested in a first-pass detection of GI in newly sequenced genomes.

BlastXtract2: Improving early exploration of (meta) genomes.

To manage and intelligently mine the avalanche of genomic sequences intuitive and user-friendly graphical interfaces are required. Here we present BlastXtract2 which exclusively facilitates early exploration of un-annotated genomic and metagenomic sequences. Various formats of translated searches, including the commonly used BlastX, of multiple sequences against multiple protein databases can be uploaded to a relational database server, which can be accessed via a locally installed web-server. There, an intuitive GUI allows straightforward data-mining and enables quick detection of potential frameshifts and poorly sequenced or assembled regions, thereby contributing in making BlastXtract2 a unique and valuable tool for early exploration of (meta)genomic sequences.AvailabilitySource code, documentation and an online demo version are available at https://github.com/ ClaessonLab/BlastXtract2

Identification of genes preferentially expressed in testis and spermatogonial cells of Labeo rohita by subtractive and suppressive hybridization

Genetic information linked to the testicular development of Labeo rohita (rohu) could contribute towards understanding and improving breeding efficiency. A rohu testis-enriched cDNA library was generated by the subtractive and suppressive hybridization technique to identify differentially expressed genes. Two-hundred and forty differentially screened clones were sequenced and 77 non-redundant EST sequences with an average size of 0.5 kb were screened against GenBank entries. Homology analysis using the blast tool identified 35% as known genes, 40% as uncharacterized showing homology to unannotated ESTs or genomic sequences and 25% as unknown sequences. Four novel transcripts were identified in the reverse transcription-polymerase chain reaction expression profiles of a panel of gonadal tissues including spermatogonial and somatic tissues, three of which were transcribed exclusively in testicular spermatogonial cells, and one was expressed in both gonads. The POU50-domain-containing transcript was predominantly expressed in male germ cells. The novel transcripts could be utilized as biomarkers for fish spermatogonial cells for efficient isolation, characterization and propagation of differentiated and/or undifferentiated cells for studying the process of spermatogenesis.

Ab initio construction of a eukaryotic transcriptome by massively parallel mRNA sequencing

Defining the transcriptome, the repertoire of transcribed regions encoded in the genome, is a challenging experimental task. Current approaches, relying on sequencing of ESTs or cDNA libraries, are expensive and labor-intensive. Here, we present a general approach for ab initio discovery of the complete transcriptome of the budding yeast, based only on the unannotated genome sequence and millions of short reads from a single massively parallel sequencing run. Using novel algorithms, we automatically construct a highly accurate transcript catalog. Our approach automatically and fully defines 86% of the genes expressed under the given conditions, and discovers 160 previously undescribed transcription units of 250 bp or longer. It correctly demarcates the 5' and 3' UTR boundaries of 86 and 77% of expressed genes, respectively. The method further identifies 83% of known splice junctions in expressed genes, and discovers 25 previously uncharacterized introns, including 2 cases of condition-dependent intron retention. Our framework is applicable to poorly understood organisms, and can lead to greater understanding of the transcribed elements in an explored genome.

Metal-responsive gene regulation and metal transport in Helicobacter species.

Helicobacter species are among the most successful colonizers of the mammalian gastrointestinal and hepatobiliary tract. Colonization is usually lifelong, indicating that Helicobacter species have evolved intricate mechanisms of dealing with stresses encountered during colonization of host tissues, like restriction of essential metal ions. The recent availability of genome sequences of the human gastric pathogen Helicobacter pylori, the murine enterohepatic pathogen Helicobacter hepaticus and the unannotated genome sequence of the ferret gastric pathogen Helicobacter mustelae has allowed for comparitive genome analyses. In this review we present such analyses for metal transporters, metal-storage and metal-responsive regulators in these three Helicobacter species, and discuss possible contributions of the differences in metal metabolism in adaptation to the gastric or enterohepatic niches occupied by Helicobacter species.

The Egg apparatus 1 gene from maize is a member of a large gene family found in both monocots and dicots

The maize ZmEA1 protein was recently postulated to be involved in short-range pollen tube guidance from the embryo sac. To date, EA1-like sequences had only been identified in monocot species. Using a more conserved C-terminal motif found in the monocot species, numerous ZmEA1-like sequences were retrieved in EST databases from dicot species, as well as from unannotated genomic sequences of Arabidopsis thaliana. RT-PCR analyses were produced for these unannotated genes and showed that these were indeed expressed genes. Further structural and phylogenetic analyses revealed that all members of the EA1-like (EAL) gene family shared a conserved 27-29 amino acid motif, termed the EA box near the C-terminal end, and appear to be secretory proteins. Therefore, the EA box proteins defines a new class of small secretory proteins, some of which being possibly involved in pollen tube guidance.

Analysis of secreted proteins from Aspergillus flavus

MS/MS techniques in proteomics make possible the identification of proteins from organisms with little or no genome sequence information available. Peptide sequences are obtained from tandem mass spectra by matching peptide mass and fragmentation information to protein sequence information from related organisms, including unannotated genome sequence data. This peptide identification data can then be grouped and reconstructed into protein data. In this study, we have used this approach to study protein secretion by Aspergillus flavus, a filamentous fungus for which very little genome sequence information is available. A. flavus is capable of degrading the flavonoid rutin (quercetin 3-O-glycoside), as the only source of carbon via an extracellular enzyme system. In this continuing study, a proteomic analysis was used to identify secreted proteins from A. flavus when grown on rutin. The growth media glucose and potato dextrose were used to identify differentially expressed secreted proteins. The secreted proteins were analyzed by 1- and 2-DE and MS/MS. A total of 51 unique A. flavus secreted proteins were identified from the three growth conditions. Ten proteins were unique to rutin-, five to glucose- and one to potato dextrose-grown A. flavus. Sixteen secreted proteins were common to all three media. Fourteen identifications were of hypothetical proteins or proteins of unknown functions. To our knowledge, this is the first extensive proteomic study conducted to identify the secreted proteins from a filamentous fungus.

BlastXtract—a new way of exploring translated searches

Searches of translated, unannotated genomic DNA sequences against protein databases is a useful early-stage method for discovering protein homologues encoded by the sequence, but generates huge amounts of output data that quickly become impregnable. BlastXtract is a web-based tool for managing and visualizing results from large translated BLAST and FastA searches. It combines the speed and storage benefits of relational database management systems with an easy-to-use graphical navigation map, and greatly facilitates the early exploration of genomic sequence. BlastXtract can be downloaded from http://bioinfo.ucc.ie/blastxtract/.

SWIFT (sequence-wide investigation with Fourier transform): a software tool for identifying proteins of a given class from the unannotated genome sequence

The ever increasing number of sequenced genomes calls for new analysis techniques, which can benefit from the methodologies developed in the field of signal processing. The present paper addresses the question of searching a pattern of amino acids (not necessarily completely specified) by means of the cross-correlation of complex sequences, obtained after suitable coding of the original amino acid sequence. Subsequently, the proposed algorithm provides a flexible strategy in setting the border between the accepted and rejected ORFs, by means of the k-means clustering of the candidate ORFs. The search for the class of proteins specified by the pattern is carried out from the most basic level, i.e. the DNA sequence, without sifting through an ensemble of previously determined ORFs. Thus, an exhaustive examination of all the occurrences of the pattern in the genome is performed. The application of the method to the search of surface proteins in Gram-positive bacteria witnesses its efficacy, in terms of both sensitivity and specificity. The comparison with the usual (and somewhat arbitrary) choice of setting a fixed value for the threshold length of the putative ORF confirms the validity of the proposed approach.

Robust Method for Proteome Analysis by MS/MS Using an Entire Translated Genome: Demonstration on the Ciliome ofTetrahymenathermophila

To improve the utility of increasingly large numbers of available unannotated and initially poorly annotated genomic sequences for proteome analysis, we demonstrate that effective protein identification can be made on a large and unannotated genome. The strategy developed is to translate the unannotated genome sequence into amino acid sequence encoding putative proteins in all six reading frames, to identify peptides by tandem mass spectrometry (MS/MS), to localize them on the genome sequence, and to preliminarily annotate the protein via a similarity search by BLAST. These tasks have been optimized and automated. Optimization to obtain multiple peptide matches in effect extends the searchable region and results in more robust protein identification. The viability of this strategy is demonstrated with the identification of 223 cilia proteins in the unicellular eukaryotic model organism Tetrahymena thermophila, whose initial genomic sequence draft was released in November 2003. To the best of our knowledge, this is the first demonstration of large-scale protein identification based on such a large, unannotated genome. Of the 223 cilia proteins, 84 have no similarity to proteins in NCBI's nonredundant (nr) database. This methodology allows identifying the locations of the genes encoding these novel proteins, which is a necessary first step to downstream functional genomic experimentation.

Genomic peculiarity of coding sequences and metabolic potential of probiotic Escherichia coli strain Nissle 1917 inferred from raw genome data

Probiotic Escherichia coli strain Nissle 1917 (O6:K5:H1) is a commensal E. coli isolate that has a long tradition in medicine for the treatment of various intestinal disorders in humans. To elucidate the molecular basis of its probiotic nature, we started sequencing the genome of this organism with a whole-genome shotgun approach. A 7.8-fold coverage of the genomic sequence has been generated and is now in the finishing stage. To exploit the genome data as early as possible and to generate hypotheses for functional studies, the unfinished sequencing data were analyzed in this work using a new method [Sun, J., Zeng, A.P., 2004. IdentiCS—identification of coding sequence and in silico reconstruction of the metabolic network directly from unannotated low-coverage bacterial genome sequence. BMC Bioinformatics 5, 112] which is particularly suitable for the prediction of coding sequences (CDSs) from unannotated genome sequence. The CDSs predicted for E. coli Nissle 1917 were compared with those of all five other sequenced E. coli strains ( E. coli K-12 MG1655, E. coli K-12 W3110, E. coli CFT073, EHEC O157:H7 EDL933 and EHEC O157:H7 Sakai) published to date. Five thousand one hundred and ninety-two CDSs were predicted for E. coli Nissle 1917, of which 1065 were assigned with enzyme EC numbers. The comparison of all predicted CDSs of E. coli Nissle 1917 to the other E. coli strains revealed 108 CDSs specific for this isolate. They are organized as four big genome islands and many other smaller gene clusters. Based on CDSs with EC numbers for enzymes, the potential metabolic network of Nissle 1917 was reconstructed and compared to those of the other five E. coli strains. Overall, the comparative genomic analysis sheds light on the genomic peculiarity of the probiotic E. coli strain Nissle 1917 and is helpful for designing further functional studies long before the sequencing project is completely finished.

Multi-criterial coding sequence prediction. Combination of GeneMark with two novel, coding-character specific quantities

This work applies two recently formulated quantities, strongly correlated with the coding character of a sequence, as an additional “module” on GeneMark, in a three-criterial method. The difference in the statistical approaches implicated by the methods combined here, is expected to contribute to an efficient assignment of functionality to unannotated genomic sequences. The developed combined algorithm is used to fractionalize a collection of GeneMark-predicted exons into sub-collections of different expectation to be coding. A further modification of the algorithm allows for the assignment of an improved estimation of the probability to be coding, to GeneMark-predicted exons. This is on the basis of a suitable training set of GeneMark-predicted exons of known functionality.

Transcription of the Toxin Genes Present within the Staphylococcal Phage φSa3ms Is Intimately Linked with the Phage's Life Cycle

phiSa3ms, a lysogenic bacteriophage encoding the staphylococcal enterotoxins SEA, SEG, and SEK and the fibrinolytic enzyme staphylokinase (Sak), was identified in the unannotated genome sequence of the hypervirulent community-acquired Staphylococcus aureus strain 476. We found that mitomycin C induction of phiSa3ms led to increased transcription of all four virulence factors. The increase in sea and sak transcription was a result of read-through transcription from upstream latent phage promoters and an increase in phage copy number. The majority of the seg2 and sek2 transcripts were shown to initiate from the upstream phage cI promoter and hence were regulated by factors influencing cI transcription. The lysogeny module of phiSa3ms was shown to have some lambda-like features with divergent cI and cro genes. Band shift assays were used to identify binding sites for both CI and Cro within the region between these genes, suggesting a mechanism of control for the phiSa3ms lytic-lysogenic switch. Our findings suggest that the production of phage-encoded virulence factors in S. aureus may be regulated by processes that govern lysogeny.

Mouse Wnt9b transforming activity, tissue-specific expression, and evolution

The members of the Wnt family of secreted factors have oncogenic potential and important roles as developmental regulators. We report an analysis of mouse Wnt9b (also called Wnt15 and Wnt14b), including its cDNA sequence, chromosomal mapping, epithelial cell transforming activity, adult and embryonic tissue expression patterns, and evolution. We also deduced the full-length amino acid sequence of its close relative, Wnt9a (also called Wnt14), from unannotated genomic DNA sequences in GenBank. Full-length comparisons among Wnt amino acid sequences provide evidence that Wnt9b and Wnt9a are close paralogs of each other and are orthologs of Wnt9 genes from shark and hagfish. Mapping Wnt9b to The Jackson Laboratory BSS interspecific backcross panel places it at 63.0 cM on chromosome 11. Sequence comparisons of two pairs of linked Wnt genes (the Wnt9a– Wnt3a pair and the Wnt9b– Wnt3 pair) suggest that they arose from the relatively recent duplication of a single ancestral Wnt gene pair, confirming the close paralogous relationship of Wnt9a and Wnt9b. Wnt9b expression is primarily restricted to the kidney in the adult mouse, with lower levels detected in the preputial gland, liver, and mammary gland. Testing of staged whole mouse embryos from 9.5 to 17.5 days of gestation showed expression at all stages with a peak at day 10.5. In situ hybridization analysis showed expression in most but not all tissues of the 16.5-day embryo. No significant elevation of Wnt9b expression was detected in 29 mouse mammary tumor virus-induced tumors. Overexpression of Wnt9b in C57MG mammary epithelial cells caused small transformed foci in cell monolayers and a moderate morphological transformation in pooled colonies compared with Wnt1.

Dissimilatory Fe(III) and Mn(IV) reduction by Shewanella putrefaciens requires ferE, a homolog of the pulE (gspE) type II protein secretion gene.

Shewanella putrefaciens strain 200 respires anaerobically on a wide range of compounds as the sole terminal electron acceptor, including ferric iron [Fe(III)] and manganese oxide [Mn(IV)]. Previous studies demonstrated that a 23.3-kb S. putrefaciens wild-type DNA fragment conferred metal reduction capability to a set of respiratory mutants with impaired Fe(III) and Mn(IV) reduction activities (T. DiChristina and E. DeLong, J. Bacteriol. 176:1468-1474, 1994). In the present study, the smallest complementing fragment was found to contain one open reading frame (ORF) (ferE) whose translated product displayed 87% sequence similarity to Aeromonas hydrophila ExeE, a member of the PulE (GspE) family of proteins found in type II protein secretion systems. Insertional mutants E726 and E912, constructed by targeted replacement of wild-type ferE with an insertionally inactivated ferE construct, were unable to respire anaerobically on Fe(III) or Mn(IV) yet retained the ability to grow on all other terminal electron acceptors. Nucleotide sequence analysis of regions flanking ferE revealed the presence of one partial and two complete ORFs whose translated products displayed 55 to 70% sequence similarity to the PulD, -F, and -G homologs of type II secretion systems. A contiguous cluster of 12 type II secretion genes (pulC to -N homologs) was found in the unannotated genome sequence of Shewanella oneidensis (formerly S. putrefaciens) MR-1. A 91-kDa heme-containing protein involved in Fe(III) reduction was present in the peripheral proteins loosely attached to the outside face of the outer membrane of the wild-type and complemented (Fer+) B31 transconjugates yet was missing from this location in Fer mutants E912 and B31 and in uncomplemented (Fer-) B31 transconjugates. Membrane fractionation studies with the wild-type strain supported this finding: the 91-kDa heme-containing protein was detected with the outer membrane fraction and not with the inner membrane or soluble fraction. These findings provide the first genetic evidence linking dissimilatory metal reduction to type II protein secretion and provide additional biochemical evidence supporting outer membrane localization of S. putrefaciens proteins involved in anaerobic respiration on Fe(III) and Mn(IV).

A Novel Family of Chitin-binding Proteins from Insect Type 2 Peritrophic Matrix: cDNA SEQUENCES, CHITIN BINDING ACTIVITY, AND CELLULAR LOCALIZATION

The peritrophic matrix is a prominent feature of the digestive tract of most insects, but its function, formation, and even its composition remain contentious. This matrix is a molecular sieve whose toughness and elasticity are generated by glycoproteins, proteoglycans, and chitin fibrils. We now describe a small, highly conserved protein, peritrophin-15, which is an abundant component of the larval peritrophic matrices of the Old World screwworm fly, Chrysomya bezziana, and sheep blowfly, Lucilia cuprina. Their deduced amino acid sequences code for a 8-kDa secreted protein characterized by a highly conserved and novel register of six cysteines. Two Drosophila homologues have also been identified from unannotated genomic sequences. Recombinant peritrophin-15 binds strongly and specifically to chitin; however, the stoichiometry of binding is low (1:10,000 N-acetyl glucosamine). We propose that peritrophin-15 caps the ends of the chitin polymer. Immunogold studies localized peritrophin-15 to the peritrophic matrix and specific vesicles in cells of the cardia, the small organ of the foregut responsible for peritrophic matrix synthesis. The vesicular contents are disgorged at the base of microvilli underlying the newly formed peritrophic matrix. This is the first time that the process of synthesis and integration of a peritrophic matrix protein into the nascent peritrophic matrix has been observed.

Combining structural analysis of DNA with search routines for the detection of transcription regulatory elements.

Analysis of unannotated genomic sequences for regulatory regions depends on a reliable recognition of individual cis-acting elements. Since some of them have a very low conserved sequence pattern, additional criteria are required. Using molecular modelling techniques, we have created a complete database for the conversion of base sequences into profiles of structural parameters of DNA. On this basis, search routines can be developed that scan for profile matches. They may be used instead of or, probably most appropriate in most cases, in combination with conventional sequence pattern searches.