Functional region prediction with a set of appropriate homologous sequences-an index for sequence selection by integrating structure and sequence information with spatial statistics

ABSTRACT: BACKGROUND: The detection of conserved residue clusters on a protein structure is one of the effective strategies for the prediction of functional protein regions. Various methods, such as Evolutionary Trace, have been developed based on this strategy. In such approaches, the conserved residues are identified through comparisons of homologous amino acid sequences. Therefore, the selection of homologous sequences is a critical step. It is empirically known that a certain degree of sequence divergence in the set of homologous sequences is required for the identification of conserved residues. However, the development of a method to select homologous sequences appropriate for the identification of conserved residues has not been sufficiently addressed. An objective and general method to select appropriate homologous sequences is desired for the efficient prediction of functional regions. RESULTS: We have developed a novel index to select the sequences appropriate for the identification of conserved residues, and implemented the index within our method to predict the functional regions of a protein. The implementation of the index improved the performance of the functional region prediction. The index represents the degree of conserved residue clustering on the tertiary structure of the protein. For this purpose, the structure and sequence information were integrated within the index by the application of spatial statistics. Spatial statistics is a field of statistics in which not only the attributes but also the geometrical coordinates of the data are considered simultaneously. Higher degrees of clustering generate larger index scores. We adopted the set of homologous sequences with the highest indexscore, under the assumption that the best prediction accuracy is obtained when the degree of clustering is the maximum. The set of sequences selected by the index led to higher functional region prediction performance than the sets of sequences selected by other sequence-based methods. CONCLUSIONS: Appropriate homologous sequences are selected automatically and objectively by the index. Such sequence selection improved the performance of functional region prediction. As far as we know, this is the first approach in which spatial statistics have been applied t o protein analyses. Such integration of structure and sequence information would be useful for other bioinformatics problems.     

Random insertion mutagenesis of the intracellular pathogen Rhodococcus equi using transposomes

The identification of virulence factors in Rhodococcus equi has been severely hampered by the lack of a method for in vivo random insertion mutagenesis. This study reports the use of transposomes to generate random insertions of a gene conferring kanamycin resistance into the genome of R. equi ATCC 33701. Southern hybridisation using the kanamycin resistance gene as probe showed that insertion of transposome is random. This was confirmed following nucleotide sequence analysis of the junction between the transposome and chromosomal DNA. The presence of a 9 bp duplication of the target sequence showed that random integration of the transposome was due to a bona fide Tn5 transposition event.     

Event-specific qualitative and quantitative PCR detection of the GMO carnation (Dianthus caryophyllus) variety Moonlite based upon the 5'-transgene integration sequence

To ensure the implementation of genetically modified organism (GMO)-labeling regulations, an event-specific detection method was developed based on the junction sequence of an exogenous integrant in the transgenic carnation variety Moonlite. The 5'-transgene integration sequence was isolated by thermal asymmetric interlaced PCR. Based upon the 5'-transgene integration sequence, the event-specific primers and TaqMan probe were designed to amplify the fragments, which spanned the exogenous DNA and carnation genomic DNA. Qualitative and quantitative PCR assays were developed employing the designed primers and probe. The detection limit of the qualitative PCR assay was 0.05% for Moonlite in 100 ng total carnation genomic DNA, corresponding to about 79 copies of the carnation haploid genome; the limit of detection and quantification of the quantitative PCR assay were estimated to be 38 and 190 copies of haploid carnation genomic DNA, respectively. Carnation samples with different contents of genetically modified components were quantified and the bias between the observed and true values of three samples were lower than the acceptance criterion (<25%) of the GMO detection method. These results indicated that these event-specific methods would be useful for the identification and quantification of the GMO carnation Moonlite.     

Characterisation of IS901 integration sites in the Mycobacterium avium genome

Data are presented on the identification and characterisation of 17 chromosomal integration loci of the insertion element IS901 in the Mycobacterium avium (cervine strain JD88/118) genome. Thirteen of these integration loci have been mapped to their corresponding positions on the M. avium strain 104 (an IS901(-) strain) genome (The Institute for Genome Research (TIGR) unfinished genome-sequencing project). Sequence data for both upstream and downstream sequence flanking regions were obtained for 12 insertion loci, while upstream sequence was obtained for five others. A consensus IS901 insertion target sequence compiled from all 17 integration sites was in broad agreement with earlier reports that were based on only two such loci. Analysis of IS901 integration site flanking sequences revealed that, like IS900 in M. avium subspecies paratuberculosis, IS901 inserts preferentially between a putative ribosome-binding sequence (RBS) and the translational start codon of an open reading frame (ORF). In BLAST X and BLAST P searches of the GenBank database, these ORFs were shown to share significant homologies with a number of other prokaryotic genes.     

高通量测序检测金针菇RNAi转化子插入位点及拷贝数

本研究对金针菇Flammulina velutipes的一个RNAi转化子菌株1382R3进行了高通量测序,以本实验室先前获得的野生型W23基因组数据为参考,分析了该转化子的基因插入位点以及拷贝数。转化子菌株1382R3是通过农杆菌介导将fv-hmg1-RNAi载体转化至金针菇菌株并通过PCR检测筛选标记而得到。通过BLAST将转化子测序的reads对外源载体和基因组定位,找到具有基因组序列（GS）和外源载体序列（ES）两种序列的临界reads,并据此使用PERL语言程序成功在转化子1382R3菌株中找到两个插入位点。对两个插入位置的序列分析表明：在插入位点1,T-DNA片段部分插入;在插入位点2,T-DNA全部插入到基因组。两个插入位点都对基因组内源基因的表达造成了一定的干扰。此方法拓宽了高通量测序技术的应用范围,将其运用到遗传转化插入位置和拷贝数的研究中,有利于食用菌的功能基因组及基因工程研究。     

SynBrowse: a synteny browser for comparative sequence analysis

MOTIVATION: The recent efforts of various sequence projects to sequence deeply into various phylogenies provide great resources for comparative sequence analysis. A generic and portable tool is essential for scientists to visualize and analyze sequence comparisons. RESULTS: We have developed SynBrowse, a synteny browser for visualizing and analyzing genome alignments both within and between species. It is intended to help scientists study macrosynteny, microsynteny and homologous genes between sequences. It can also aid with the identification of uncharacterized genes, putative regulatory elements and novel structural features of a species. SynBrowse is a GBrowse (the Generic Genome Browser) family software tool that runs on top of the open source BioPerl modules. It consists of two components: a web-based front end and a set of relational database back ends. Each database stores pre-computed alignments from a focus sequence to reference sequences in addition to the genome annotations of the focus sequence. The user interface lets end users select a key comparative alignment type and search for syntenic blocks between two sequences and zoom in to view the relationships among the corresponding genome annotations in detail. SynBrowse is portable with simple installation, flexible configuration, convenient data input and easy integration with other components of a model organism system. AVAILABILITY: The software is available at http://www.gmod.org CONTACT: vbrendel@iastate.edu     

ProPred1: prediction of promiscuous MHC Class-I binding sites

SUMMARY: ProPred1 is an on-line web tool for the prediction of peptide binding to MHC class-I alleles. This is a matrix-based method that allows the prediction of MHC binding sites in an antigenic sequence for 47 MHC class-I alleles. The server represents MHC binding regions within an antigenic sequence in user-friendly formats. These formats assist user in the identification of promiscuous MHC binders in an antigen sequence that can bind to large number of alleles. ProPred1 also allows the prediction of the standard proteasome and immunoproteasome cleavage sites in an antigenic sequence. This server allows identification of MHC binders, who have the cleavage site at the C terminus. The simultaneous prediction of MHC binders and proteasome cleavage sites in an antigenic sequence leads to the identification of potential T-cell epitopes. AVAILABILITY: Server is available at http://www.imtech.res.in/raghava/propred1/. Mirror site of this server is available at http://bioinformatics.uams.edu/mirror/propred1/ Supplementary information: Matrices and document on server are available at http://www.imtech.res.in/raghava/propred1/page2.html     

16S rRNA sequencing in routine bacterial identification: a 30-month experiment

Accurate identification of bacterial isolates is an essential task in clinical microbiology. Phenotypic methods are time-consuming and either fail to identify some bacteria such as Gram-positive rods entirely or at least fail to do so in some clinical situations. 16S rDNA sequencing is a recent method of identification which offers a useful alternative. In this study, we investigate the usefulness of this method for identifying a range of bacteria in a clinical laboratory under routine conditions. Over a period of 30 months, 683 isolates were obtained from clinical specimens, sequenced and analysed. For 568 of these isolates (83.1%), the sequence provided species level identification. For 108 isolates (15.8%), the identification was limited to the genus level, and for 7 isolates (1%), the sequence remained unidentifiable by 16S rDNA sequence analysis. For the isolates identified only to the genus level, the 16S rDNA approach failed to identify bacteria to the taxonomic level for 3 reasons: failure to differentiate between species in 72 isolates (66%), the lack of any closely related sequence in the database for 15 isolates (13.8%) and the presence of more than 1% of undetermined position in the sequence for 13 isolates (12%).     

Identifying and genotyping transgene integration loci

The random germline integration of genetically engineered transgenes has been a powerful technique to study the role of particular genes in variety of biological processes. Although the identification of the transgene insertion site is often not essential for functional analysis of the transgene, identifying the site can have practical benefit. Enabling one to distinguish between animals that are homozygous or hemizygous for the transgene locus could facilitate breeding strategies to produce animals with a large number of genetic markers. Furthermore, founder lines generated with the same transgene construct may exhibit different phenotypes and levels of transgene expression depending on the site of integration. The goal of this report was to develop a rapid protocol for the identification and verification of transgene insertion sites. To identify host genomic sequences at the coagulation Factor X transgene integration site, DNA from a tail snip of the transgenic mouse was digested with NcoI and circularized using T4 DNA ligase. Using appropriately positioned PCR primers annealing to a transgene fragment distal to a terminal transgene restriction site (NcoI), one could amplify a fragment containing the transgene terminal region and extending into the flanking genomic sequence at the insertion site. DNA sequence determination of the amplicon permitted identification of the insertion site using a BLASTN search. FISH analysis of a metaphase spread of primary fibroblasts derived from the transgenic mouse was consistent with the identification of insertion site near the end of mouse chromosome 14. Identification of transgene insertion sites will facilitate genotyping strategies useful for the construction of mice with multiple engineered genetic markers and to distinguish among different founder lines generated by the same transgene. Furthermore, identification of the insertion site is necessary to analyze unexpected phenotypes that might be caused by insertional inactivation of an endogenous gene.     

Expanding the organismal scope of proteomics: cross-species protein identification by mass spectrometry and its implications

Due to the limited applicability of conventional protein identification methods to the proteomes of organisms with unsequenced genomes, researchers have developed approaches to identify proteins using mass spectrometry and sequence similarity database searches. Both the integration of mass spectrometry with bioinformatics and genomic sequencing drive the expanding organismal scope of proteomics.     

Selection of bead‐displayed,PNA‐encoded chemicals

The lack of efficient identification and isolation methods for specific molecular binders has fundamentally limited drug discovery. Here, we have developed a method to select peptide nucleic acid (PNA) encoded molecules with specific functional properties from combinatorially generated libraries. This method consists of three essential stages: (1) creation of a Lab‐on‐Bead? library, a one‐bead, one‐sequence library that, in turn, displays a library of candidate molecules, (2) fluorescence microscopy‐aided identification of single target‐bound beads and the extraction – wet or dry – of these beads and their attached candidate molecules by a micropipette manipulator, and (3) identification of the target‐binding candidate molecules via amplification and sequencing. This novel integration of techniques harnesses the sensitivity of DNA detection methods and the multiplexed and miniaturized nature of molecule screening to efficiently select and identify target‐binding molecules from large nucleic acid encoded chemical libraries. Beyond its potential to accelerate assays currently used for the discovery of new drug candidates, its simple bead‐based design allows for easy screening over a variety of prepared surfaces that can extend this technique's application to the discovery of diagnostic reagents and disease markers. Copyright © 2009 John Wiley & Sons, Ltd.     

Number and Accuracy of T-DNA Insertions in Transgenic Banana (Musa spp.) Plants Characterized by an Improved Anchored PCR Technique

Nineteen transgenic banana plants, produced via Agrobacterium-mediated transformation, were analyzed for the integration of T-DNA border regions using an improved anchored PCR technique. The method described is a relatively fast, three-step procedure (restriction digestion of genomic DNA, ligation of ‘vectorette’-type adaptors, and a single round of suppression PCR) for the amplification of specific T-DNA border-containing genomic fragments. Most transgenic plants carried a low number of inserts and the method was suitable for a detailed characterization of the integration events, including T-DNA border integrity as well as the insertion of non-T-DNA vector sequences, which occurred in 26% of the plants. Furthermore, the particular band pattern generated by four enzyme/primer combinations for each individual plant served as a fingerprint, allowing the identification of plants representing identical transformation events. Genomic Southern hybridization and nucleotide sequence analysis of amplification products confirmed the data obtained by anchored PCR. Sequencing of seven right or left border junction regions revealed different T-DNA processing events for each plant, indicating a relatively low frequency of precisely nicked T-DNA integration among the plants studied.     

Large-scale discovery of insertion hotspots and preferential integration sites of human transposed elements

Throughout evolution, eukaryotic genomes have been invaded by transposable elements (TEs). Little is known about the factors leading to genomic proliferation of TEs, their preferred integration sites and the molecular mechanisms underlying their insertion. We analyzed hundreds of thousands nested TEs in the human genome, i.e. insertions of TEs into existing ones. We first discovered that most TEs insert within specific ‘hotspots’ along the targeted TE. In particular, retrotransposed Alu elements contain a non-canonical single nucleotide hotspot for insertion of other Alu sequences. We next devised a method for identification of integration sequence motifs of inserted TEs that are conserved within the targeted TEs. This method revealed novel sequences motifs characterizing insertions of various important TE families: Alu, hAT, ERV1 and MaLR. Finally, we performed a global assessment to determine the extent to which young TEs tend to nest within older transposed elements and identified a 4-fold higher tendency of TEs to insert into existing TEs than to insert within non-TE intergenic regions. Our analysis demonstrates that TEs are highly biased to insert within certain TEs, in specific orientations and within specific targeted TE positions. TE nesting events also reveal new characteristics of the molecular mechanisms underlying transposition.     

一株分离自酸梅酱的酵母菌的分离鉴定

目的:分离及鉴定来自于新疆塔城民间自制酸梅酱中的酵母菌。方法:用NL1/NL4引物对扩增酵母菌株的26S rDNAD1/D2区,测序结果进行序列分析,用Neighbour-joining(N-J)方法构建系统发育进化树,同时结合酵母的传统形态学鉴定对菌株进行鉴定。结果:26S rDNA D1/D2区序列分析表明菌株KKS与Metschnikowia aff.fructicola D3895相近,相似率为99.2%。在N-J法构建的系统发育进化树中,菌株KKS与Metschnikowia aff.fructicola聚类在同一分枝上。结论:从新疆塔城民间自制酸梅酱中分离得到一株酵母菌并将该菌株鉴定为Metschnikowia aff.fructicola。     

Evaluation of the rrn operon copy number in Bifidobacterium using real-time PCR

AIMS: A real-time PCR-based method was developed to evaluate the Bifidobacterium rRNA operon copy number. As a result of their repetitive nature, rRNA operons are very suitable targets for chromosomal integration of heterologous genes. METHODS AND RESULTS: The rrn operon multiplicity per chromosome was determined by real-time PCR quantification of the 16S rRNA amplicons obtained from genomic DNA. The values obtained in several bifidobacterial strains of human origin ranged from 1 to 5. The reliability of the method developed was confirmed by Southern hybridization technique. CONCLUSIONS: In the Bifidobacterium genus the rrn operon copies showed variability at species and strain level. The identification of Bifidobacterium strains with high rRNA multiplicity allowed the selection of potential hosts for chromosomal integration. SIGNIFICANCE AND IMPACT OF THE STUDY: The methodology here proposed represents a rapid, reliable and sensitive new tool for the quantification of rrn operon copy number in bacteria.     

HOW TO CONSTRUCT THE EXPERT SYSTEM FOR SPECIES IDENTIFICATION USING VENATION OF TORTRICINAE (LEPIDOPTERA)

Abstract  The author has reported the possibility to use venation for identification of tortricids. In the previous paper (Liu 1996), an expert system for identification using venation was also described. The present paper introduces the technique how to construct the expert system. Three original solutions in the research are published. They are: establishment of corrdinates of venation—"two-centre method", the interface to accept inputs—"model modification method" and the interface to generate outputs—"sequence method".     

Identification of metabolic network models from incomplete high-throughput datasets

MOTIVATION: High-throughput measurement techniques for metabolism and gene expression provide a wealth of information for the identification of metabolic network models. Yet, missing observations scattered over the dataset restrict the number of effectively available datapoints and make classical regression techniques inaccurate or inapplicable. Thorough exploitation of the data by identification techniques that explicitly cope with missing observations is therefore of major importance. RESULTS: We develop a maximum-likelihood approach for the estimation of unknown parameters of metabolic network models that relies on the integration of statistical priors to compensate for the missing data. In the context of the linlog metabolic modeling framework, we implement the identification method by an Expectation-Maximization (EM) algorithm and by a simpler direct numerical optimization method. We evaluate performance of our methods by comparison to existing approaches, and show that our EM method provides the best results over a variety of simulated scenarios. We then apply the EM algorithm to a real problem, the identification of a model for the Escherichia coli central carbon metabolism, based on challenging experimental data from the literature. This leads to promising results and allows us to highlight critical identification issues.