DNA-based watermarks using the DNA-Crypt algorithm

Background  

The aim of this paper is to demonstrate the application of watermarks based on DNA sequences to identify the unauthorized use of genetically modified organisms (GMOs) protected by patents. Predicted mutations in the genome can be corrected by the DNA-Crypt program leaving the encrypted information intact. Existing DNA cryptographic and steganographic algorithms use synthetic DNA sequences to store binary information however, although these sequences can be used for authentication, they may change the target DNA sequence when introduced into living organisms.     

DNA watermarks: A proof of concept

Background  

DNA-based watermarks are helpful tools to identify the unauthorized use of genetically modified organisms (GMOs) protected by patents. In silico analyses showed that in coding regions synonymous codons can be used to insert encrypted information into the genome of living organisms by using the DNA-Crypt algorithm.     

Extracting haplotypes from diploid organisms

Each diploid organism has two alleles at every gene locus. In sexual organisms such as most plants, animals and fungi, the two alleles in an individual may be genetically very different from each other. DNA sequence data from individual alleles (called a haplotype) can provide powerful information to address a variety of biological questions and guide many practical applications. The advancement in molecular technology and computational tools in the last decade has made obtaining large-scale haplotypes feasible. This review summarizes the two basic approaches for obtaining haplotypes and discusses the associated techniques and methods. The first approach is to experimentally obtain diploid sequence information and then use computer algorithms to infer haplotypes. The second approach is to obtain haplotype sequences directly through experimentation. The advantages and disadvantages of each approach are discussed. I then discussed a specific example on how the direct approach was used to obtain haplotype information to address several fundamental biological questions of a pathogenic yeast. With increasing sophistication in both bioinformatics tools and high-throughput molecular techniques, haplotype analysis is becoming an integrated component in biomedical research.     

Comparative ab initio prediction of gene structures using pair HMMs

We present a novel comparative method for the ab initio prediction of protein coding genes in eukaryotic genomes. The method simultaneously predicts the gene structures of two un-annotated input DNA sequences which are homologous to each other and retrieves the subsequences which are conserved between the two DNA sequences. It is capable of predicting partial, complete and multiple genes and can align pairs of genes which differ by events of exon-fusion or exon-splitting. The method employs a probabilistic pair hidden Markov model. We generate annotations using our model with two different algorithms: the Viterbi algorithm in its linear memory implementation and a new heuristic algorithm, called the stepping stone, for which both memory and time requirements scale linearly with the sequence length. We have implemented the model in a computer program called DOUBLESCAN. In this article, we introduce the method and confirm the validity of the approach on a test set of 80 pairs of orthologous DNA sequences from mouse and human. More information can be found at: http://www.sanger.ac.uk/Software/analysis/doublescan/     

A CRISPR–Cas9 system for multiple genome editing and pathway assembly in Candida tropicalis

Genetic manipulation is among the most important tools for synthetic biology; however, modifying multiple genes is extremely time-consuming and can sometimes be impossible when dealing with gene families. Here, we present a clustered regularly interspaced short palindromic repeat (CRISPR) and CRISPR-associated protein 9 (Cas9) system for use in the diploid yeast Candida tropicalis that is vastly superior to traditional techniques. This system enables the rapid and reliable introduction of multiple genetic deletions or mutations, as well as a stable expression using an integrated CRISPR–Cas9 cassette or a transient CRISPR–Cas9 cassette, together with a short donor DNA. We further show that the system can be used to promote the in vivo assembly of multiple DNA fragments and their stable integration into a target locus (or loci) in C. tropicalis. Based on this system, we present a platform for the biosynthesis of β-carotene and its derivatives. These results enable the practical application of C. tropicalis and the application of the system to other organisms.     

iRSpot-DTS: Predict recombination spots by incorporating the dinucleotide-based spare-cross covariance information into Chou's pseudo components

Meiotic recombination plays an important role in the process of genetic evolution. Previous researches have shown that the recombination rates provide important information about the mechanism of recombination study. However, at present, most methods ignore the hidden correlation and spatial autocorrelation of the DNA sequence. In this study, we proposed a predictor called iRSpot-DTS to identify hot/cold spots based on the benchmark datasets. We proposed a feature extraction method called dinucleotide-based spatial autocorrelation(DSA) which can incorporate the original DNA properties and spatial information of DNA sequence. Then it used t-SNE method to remove the noise which outperformed PCA. Finally, we used SAE softmax classifier to do classification which is based on networks and can get more hidden information of DNA sequence, our iRSpot-DTS achieved remarkable performance. Jackknife cross validation tests were done on two benchmark datasets. We achieved state-of-the-art results with 96.61% overall accuracy(OA), 93.16% Matthews correlation coefficient (MCC) and over 95% in Sn and Sp which are the best in this state.     

Anchor-based whole genome phylogeny (ABWGP): a tool for inferring evolutionary relationship among closely related microorganisms [corrected

Phenotypic behavior of a group of organisms can be studied using a range of molecular evolutionary tools that help to determine evolutionary relationships. Traditionally a gene or a set of gene sequences was used for generating phylogenetic trees. Incomplete evolutionary information in few selected genes causes problems in phylogenetic tree construction. Whole genomes are used as remedy. Now, the task is to identify the suitable parameters to extract the hidden information from whole genome sequences that truly represent evolutionary information. In this study we explored a random anchor (a stretch of 100 nucleotides) based approach (ABWGP) for finding distance between any two genomes, and used the distance estimates to compute evolutionary trees. A number of strains and species of Mycobacteria were used for this study. Anchor-derived parameters, such as cumulative normalized score, anchor order and indels were computed in a pair-wise manner, and the scores were used to compute distance/phylogenetic trees. The strength of branching was determined by bootstrap analysis. The terminal branches are clearly discernable using the distance estimates described here. In general, different measures gave similar trees except the trees based on indels. Overall the tree topology reflected the known biology of the organisms. This was also true for different strains of Escherichia coli. A new whole genome-based approach has been described here for studying evolutionary relationships among bacterial strains and species.     

Anchor-Based Whole Genome Phylogeny (ABWGP): A Tool for Inferring Evolutionary Relationship among Closely Related Microorganims

Phenotypic behavior of a group of organisms can be studied using a range of molecular evolutionary tools that help to determine evolutionary relationships. Traditionally a gene or a set of gene sequences was used for generating phylogenetic trees. Incomplete evolutionary information in few selected genes causes problems in phylogenetic tree construction. Whole genomes are used as remedy. Now, the task is to identify the suitable parameters to extract the hidden information from whole genome sequences that truly represent evolutionary information. In this study we explored a random anchor (a stretch of 100 nucleotides) based approach (ABWGP) for finding distance between any two genomes, and used the distance estimates to compute evolutionary trees. A number of strains and species of Mycobacteria were used for this study. Anchor-derived parameters, such as cumulative normalized score, anchor order and indels were computed in a pair-wise manner, and the scores were used to compute distance/phylogenetic trees. The strength of branching was determined by bootstrap analysis. The terminal branches are clearly discernable using the distance estimates described here. In general, different measures gave similar trees except the trees based on indels. Overall the tree topology reflected the known biology of the organisms. This was also true for different strains of Escherichia coli. A new whole genome-based approach has been described here for studying evolutionary relationships among bacterial strains and species.     

Editorial: Twin peaks mission

The Y human chromosome has many ancient genes whose fidelity seems to have been preserved by tandem sequences and palindromic ‘hairpins’, compared/repaired by ‘gene conversion’. That a primary function of recombination machinery is DNA repair has been suggested, and rejected, several times; this new evidence is very persuasive. The process, better called gene conservation than gene conversion, could operate in all diploid organisms, accounting for the retention of long gene sequences without ‘informational meltdown’ ('concerted evolution'). It resembles rocket‐science computer‐redundancy error‐checking, comparison of three or four sequences, not just two. If recognition of errors in ‘converted’ sequences can be followed by either repair or rejection, the rejection option can account for the vast wastage of meiotic products. The repair option might be used in Drosophila oocytes and even zygotic nuclei, possibly other oocytes, ancient asexual lineages such as mycorrhizal fungi, perhaps the Y itself. Both evolutionary stasis (conservatism) and development and deployment of complex developmental modules can be understood in these terms so both the evolution of biodiversity and the practice of systematics may have these mechanisms as their bases. The main individual‐fitness and evolutionary advantages of diploidy were not primarily cloaking of recessive al‐leles, or allelic recombination and Mendelism, but conserving long DNA sequences.     

SWORDS: A statistical tool for analysing large DNA sequences

In this article, we present some simple yet effective statistical techniques for analysing and comparing large DNA sequences. These techniques are based on frequency distributions of DNA words in a large sequence, and have been packaged into a software called SWORDS. Using sequences available in public domain databases housed in the Internet, we demonstrate how SWORDS can be conveniently used by molecular biologists and geneticists to unmask biologically important features hidden in large sequences and assess their statistical significance.     

Immortalization of human fibroblasts transformed by origin-defective simian virus 40.

Simian virus 40 (SV40)-mediated transformation of human diploid fibroblasts has provided an effective experimental system for studies of both "senescence" in cell culture and carcinogenesis. Previous interpretations may have been complicated, however, by the semipermissive virus-cell interaction. In earlier studies, we previously demonstrated that the human diploid fibroblast line HS74 can be efficiently transformed by DNA from replication-defective mutants of SV40 containing a deletion in the viral origin for DNA synthesis (SVori-). In the current study, we found that such SVori- transformants show a significantly increased life span in culture, as compared with either HS74 or an independent transformant containing an intact viral genome, but they nonetheless undergo senescence. We have clonally isolated six immortalized derivatives of one such transformant (SV/HF-5). Growth studies indicate that the immortalized cell lines do not invariably grow better than SV/HF-5 or HS74. Genetic studies involving karyotypic analysis and Southern analysis of integrated viral sequences demonstrated both random and nonrandom alterations. All immortalized derivatives conserved one of the two copies of SV40 sequences which expressed a truncated T antigen. These cloned SV40-transformed cell lines, pre- and postimmortalization, should be useful in defining molecular changes associated with immortalization.     

Determination of radiation equivalence of ethyl methanesulphonate (EMS) for induction of gene conversion in diploid yeast.

A unit Rad-Equivalent Chemical (REC) has been suggested for purposes of quantitating the mutagenic hazards of chemicals. The usefulness of this approach is demonstrated by the establishment of a constant relationship between the forward mutation frequency and haploid genome size in various organisms for both radiation and chemical EMS. However, it is necessary to determine the radiation equivalence of chemicals in as many organisms and for as many end-points as possible. For end-points we are limited to forward mutations. Another relevant genetic end-point of interest in this regard is gene conversion which can also monitor any kind of DNA damage in a suitable diploid system. Hence, we have determined the REC value for EMS in diploid yeast with gene conversion as the end-point. This agrees well with the REC values estimated in a number of organisms with forward mutation as the end-point. This finding further underlines the generality of the REC concept.     

HICLAS: a taxonomic database system for displaying and comparing biological classification and phylogenetic trees

MOTIVATION: Numerous database management systems have been developed for processing various taxonomic data bases on biological classification or phylogenetic information. In this paper, we present an integrated system to deal with interacting classifications and phylogenies concerning particular taxonomic groups. RESULTS: An information-theoretic view (taxon view) has been applied to capture taxonomic concepts as taxonomic data entities. A data model which is suitable for supporting semantically interacting dynamic views of hierarchic classifications and a query method for interacting classifications have been developed. The concept of taxonomic view and the data model can also be expanded to carry phylogenetic information in phylogenetic trees. We have designed a prototype taxonomic database system called HICLAS (HIerarchical CLAssification System) based on the concept of taxon view, and the data models and query methods have been designed and implemented. This system can be effectively used in the taxonomic revisionary process, especially when databases are being constructed by specialists in particular groups, and the system can be used to compare classifications and phylogenetic trees. AVAILABILITY: Freely available at the WWW URL: http://aims.cps.msu.edu/hiclas/ CONTACT: pramanik@cps.msu.edu; lotus@wipm.whcnc.ac.cn     

利用CRISPR/Cas9对基因组中高度同源DNA片段编辑多样性的遗传学研究

在基因组中,编码区存在许多高度相似的基因簇或基因群(多拷贝基因),非编码区也存在大量的重复序列。这些重复序列能通过改变染色体的三维结构调控基因的转录,对于生物体的遗传与进化起到了重要的作用。其高度同源的特征使得利用CRISPR/Cas9技术进行基因组编辑时面临更加复杂的状况。如果编辑的片段是二倍体或多倍体,还会产生各条染色单体上的编辑情况不相同的现象。为此本文选择了2个位于同一染色体相距11 kb的高度同源300 bp片段(L1和L2)进行CRISPR介导的DNA片段编辑。采用一对sgRNA(分别共同靶向两片段的上、下游位点)引导Cas9对HepG2细胞两个高度相似的DNA片段进行切割。片段编辑的细胞进一步单克隆化后,对获得的22个L1/L2编辑的CRISPR单克隆细胞株进行详细的基因型鉴定。结果发现除了这两个DNA片段本身被删除外,它们之间的大片段也存在被删除的现象,三个片段的各种反转组合也很频繁。该研究结果对于采用CRISPR/Cas9系统编辑多拷贝基因或重复序列,尤其是对二倍体或多倍体生物进行基因组编辑时具有重要的借鉴和参考价值。     

‘Yeast Mail’: A Novel Saccharomyces Application (NSA) to Encrypt Messages

The universal genetic code is used by all life forms to encode biological information. It can also be used to encrypt semantic messages and convey them within organisms without anyone but the sender and recipient knowing, i.e., as a means of steganography. Several theoretical, but comparatively few experimental, approaches have been dedicated to this subject, so far. Here, we describe an experimental system to stably integrate encrypted messages within the yeast genome using a polymerase chain reaction (PCR)‐based, one‐step homologous recombination system. Thus, DNA sequences encoding alphabetical and/or numerical information will be inherited by yeast propagation and can be sent in the form of dried yeast. Moreover, due to the availability of triple shuttle vectors, Saccharomyces cerevisiae can also be used as an intermediate construction device for transfer of information to either Drosophila or mammalian cells as steganographic containers. Besides its classical use in alcoholic fermentation and its modern use for heterologous gene expression, we here show that baker's yeast can thus be employed in a novel Saccharomyces application (NSA) as a simple steganographic container to hide and convey messages.     

FBIS: A regional DNA barcode archival & analysis system for Indian fishes

DNA barcode is a new tool for taxon recognition and classification of biological organisms based on sequence of a fragment of mitochondrial gene, cytochrome c oxidase I (COI). In view of the growing importance of the fish DNA barcoding for species identification, molecular taxonomy and fish diversity conservation, we developed a Fish Barcode Information System (FBIS) for Indian fishes, which will serve as a regional DNA barcode archival and analysis system. The database presently contains 2334 sequence records of COI gene for 472 aquatic species belonging to 39 orders and 136 families, collected from available published data sources. Additionally, it contains information on phenotype, distribution and IUCN Red List status of fishes. The web version of FBIS was designed using MySQL, Perl and PHP under Linux operating platform to (a) store and manage the acquisition (b) analyze and explore DNA barcode records (c) identify species and estimate genetic divergence. FBIS has also been integrated with appropriate tools for retrieving and viewing information about the database statistics and taxonomy. It is expected that FBIS would be useful as a potent information system in fish molecular taxonomy, phylogeny and genomics. AVAILABILITY: The database is available for free at http://mail.nbfgr.res.in/fbis/     

Application of a new probabilistic model for recognizing complex patterns in glycans

MOTIVATION: The study of carbohydrate sugar chains, or glycans, has been one of slow progress mainly due to the difficulty in establishing standard methods for analyzing their structures and biosynthesis. Glycans are generally tree structures that are more complex than linear DNA or protein sequences, and evidence shows that patterns in glycans may be present that spread across siblings and into further regions that are not limited by the edges in the actual tree structure itself. Current models were not able to capture such patterns. RESULTS: We have applied a new probabilistic model, called probabilistic sibling-dependent tree Markov model (PSTMM), which is able to inherently capture such complex patterns of glycans. Not only is the ability to capture such patterns important in itself, but this also implies that PSTMM is capable of performing multiple tree structure alignments efficiently. We prove through experimentation on actual glycan data that this new model is extremely useful for gaining insight into the hidden, complex patterns of glycans, which are so crucial for the development and functioning of higher level organisms. Furthermore, we also show that this model can be additionally utilized as an innovative approach to multiple tree alignment, which has not been applied to glycan chains before. This extension on the usage of PSTMM may be a major step forward for not only the structural analysis of glycans, but it may consequently prove useful for discovering clues into their function.     

Drifting Markov models with polynomial drift and applications to DNA sequences

In this article, we introduce the drifting Markov models (DMMs) which are inhomogeneous Markov models designed for modeling the heterogeneities of sequences (in our case DNA or protein sequences) in a more flexible way than homogeneous Markov chains or even hidden Markov models (HMMs). We focus here on the polynomial drift: the transition matrix varies in a polynomial way. To show the reliability of our models on DNA, we exhibit high similarities between the probability distributions of nucleotides obtained by our models and the frequencies of these nucleotides computed by using a sliding window. In a further step, these DMMs can be used as the states of an HMM: on each of its segments, the observed process can be modeled by a drifting Markov model. Search of rare words in DNA sequences remains possible with DMMs and according to the fits provided, DMMs turn out to be a powerful tool for this purpose. The software is available on request from the author. It will soon be integrated on seq++ library (http://stat.genopole.cnrs.fr/seqpp/).     

Heterogenomic recombinants from compatible nocardiae

Recombinants obtained from matings of Nocardia erythropolis x N. canicruria were tested for their genetic stability by comparing phenotypes from direct selection with the same population after unselected growth. Contraselective loci were employed in various combinations in order that all of the mapped characters might be subjected to unselected analysis. Some recombinant class types appeared as stable haploids, whereas others behaved as heterozygous diploids, segregating out new phenotypes. All regions of the parental genomes were found to be involved in segregation, implying that the entire mapped region can become merozygotic under standard mating conditions. On the basis of segregating phenotypes, the genetic potentials of these compatible nocardiae were ascertained as follows: the formation of a diploid with subsequent segregation of parental or haploid recombinant genomes or both; persistence of the diploid through many generations; continuing reassortment of genetic information by multiple matings between parental or recombinant organisms; and, very probably, second-round recombinations within the diploid. A considerable difference in the nuclear division time between the parental organisms was postulated to have significant effects on the nature of the unselected segregants.