A phosphamide nucleotide analog: a substrate for polymerase synthesis of DNA.

A type of modified nucleotide, deoxynucleotide γ-amidotriphosphates (dNTPγNH2s), exhibited around five times higher stability than dNTPs. These phosphamide nucleotides can be utilized by several DNA polymerases, and the amplification of a 10 kb DNA fragment through the polymerase chain reaction (PCR) can be accomplished even under conditions of high temperature, extended storage, or repeated freeze-thaw cycles. However, the control PCR with standard dNTPs was unsuccessful. These results indicate that dNTPγNH2s have the potential to substitute dNTPs in PCR.

Enzymatic Preparation of DNA with an Expanded Genetic Alphabet Using Terminal Deoxynucleotidyl Transferase and Its Applications.

Efficient preparation of DNA oligonucleotides containing unnatural nucleobases (UBs) that can pair with their cognates to form unnatural base pairs (UBPs) is an essential prerequisite for the application of UBPs in vitro and in vivo. Traditional preparation of oligonucleotides containing unnatural nucleobases largely relies on solid-phase synthesis, which needs to use unstable nucleoside phosphoramidites and a DNA synthesizer, and is environmentally unfriendly and limited in product length. To overcome these limitations of solid-phase synthesis, we developed enzymatic methods for daily laboratory preparation of DNA oligonucleotides containing unnatural nucleobase dNaM, dTPT3, or one of the functionalized dTPT3 derivatives, which can be used for orthogonal DNA labeling or the preparation of DNAs containing UBP dNaM-dTPT3, one of the most successful UBPs to date, based on the template-independent polymerase terminal deoxynucleotidyl transferase (TdT). Here, we first provide a detailed procedure for the TdT-based preparation of DNA oligonucleotides containing 3'-nucleotides of dNaM, dTPT3, or one of dTPT3 derivatives. We then present the procedures for enzyme-linked oligonucleotide assay (ELONA) and imaging of bacterial cells using DNA oligonucleotides containing 3'-nucleotides of dTPT3 derivatives with different functional groups. The procedure for enzymatic synthesis of DNAs containing an internal UBP dNaM-dTPT3 is also described. Hopefully, these methods will greatly facilitate the application of UBPs and the construction of semi-synthetic organisms with an expanded genetic alphabet.

C and G are frequently mutated into T and A in coding regions of human genes.

Nucleotide mutations in human genes have long been a hot subject for study because some of them may lead to severe human diseases. Understanding the general mutational process and evolutionary trend of human genes could help answer such questions as why certain diseases occur and what challenges we face in protecting human health. In this study, we conducted statistics on 89,895 single-nucleotide variations identified in coding regions of 18,339 human genes. The results show that C and G are frequently mutated into T and A in human genes. C/G (C or G)-to-T/A mutations lead to reduction of hydrogen bonds in double-stranded DNA because C-G and T-A base pairs are maintained by three and two hydrogen bonds respectively. C-to-T and G-to-A mutations occur predominantly in human genes because they not only reduce hydrogen bonds but also belong to transition mutation. Reduction of hydrogen bonds could reduce energy consumption not only in separating double strands of mutated DNA for transcription and replication but also in disrupting stem-loop structure of mutated mRNA for translation. It is thus considered that to reduce hydrogen bonds (and thus to reduce energy consumption in gene expression) is one of the driving forces for nucleotide mutation. Moreover, codon mutation is positively correlated to its content, suggesting that most mutations are not targeted on changing any specific codons (amino acids) but are merely for reducing hydrogen bonds. Our study provides an example of utilizing single-nucleotide variation data to infer evolutionary trend of human genes, which can be referenced to conduct similar studies in other organisms.

Resolving unknown nucleotides in the IPD-IMGT/HLA database by extended and full-length sequencing of HLA class I and II alleles.

In the past, identification of HLA alleles was limited to sequencing the region of the gene coding for the peptide binding groove, resulting in a lack of sequence information in the HLA database, challenging HLA allele assignment software programs. We investigated full-length sequences of 19 HLA class I and 7 HLA class II alleles, and we extended another 47 HLA class I alleles with sequences of 5' and 3' UTR regions that were all not yet available in the IPD-IMGT/HLA database. We resolved 8638 unknown nucleotides in the coding sequence of HLA class I and 2139 of HLA class II. Furthermore, with full-length sequencing of the 26 alleles, more than 90 kb of sequence information was added to the non-coding sequences, whereas extension of the 47 alleles resulted in the addition of 5.5 kb unknown nucleotides to the 5' UTR and > 31.7 kb to the 3' UTR region. With this information, some interesting features were observed, like possible recombination events and lineage evolutionary origins. The continuing increase in the availability of full-length sequences in the HLA database will enable the identification of the evolutionary origin and will help the community to improve the alignment and assignment accuracy of HLA alleles.

Codon Pattern and Context Analysis in Genes Triggering Alzheimer's Disease and Latent Tau Protein Aggregation Post-Anesthesia Exhibited Unique Molecular Patterns Associated with Functional Aspects.

Background: Previous reports have demonstrated post-operative dementia and Alzheimer's disease (AD), and increased amyloid-ß levels and tau hyperphosphorylation have been observed in animal models post-anesthesia. Objective: After surgical interventions, loss in memory has been observed that has been found linked with genes modulated after anesthesia. Present study aimed to study molecular pattern present in genes modulated post anesthesia and involved in characters progressing towards AD. Methods: In the present study, 17 transcript variants belonging to eight genes, which have been found to modulate post-anesthesia and contribute to AD progression, were envisaged for their compositional features, molecular patterns, and codon and codon context-associated studies. Results: The sequences' composition was G/C rich, influencing dinucleotide preference, codon preference, codon usage, and codon context. The G/C nucleotides being highly occurring nucleotides, CpGdinucleotides were also preferred; however, CpG was highly disfavored at p3-1 at the codon junction. The nucleotide composition of Cytosine exhibited a unique feature, and unlike other nucleotides, it did not correlate with codon bias. Contrarily, it correlated with the sequence lengths. The sequences were leucine-rich, and multiple leucine repeats were present, exhibiting the functional role of neuroprotection from neuroinflammation post-anesthesia. Conclusions: The analysis pave the way to elucidate unique molecular patterns in genes modulated during anesthetic treatment and might help ameliorate the ill effects of anesthetics in the future.

Regulation of micro- and small-exon retention and other splicing processes by GRP20 for flower development.

Pre-mRNA splicing is crucial for gene expression and depends on the spliceosome and splicing factors. Plant exons have an average size of ~180 nucleotides and typically contain motifs for interactions with spliceosome and splicing factors. Micro exons (<51 nucleotides) are found widely in eukaryotes and in genes for plant development and environmental responses. However, little is known about transcript-specific regulation of splicing in plants and about the regulators for micro exon splicing. Here we report that glycine-rich protein 20 (GRP20) is an RNA-binding protein and required for splicing of ~2,100 genes including those functioning in flower development and/or environmental responses. Specifically, GRP20 is required for micro-exon retention in transcripts of floral homeotic genes; these micro exons are conserved across angiosperms. GRP20 is also important for small-exon (51-100 nucleotides) splicing. In addition, GRP20 is required for flower development. Furthermore, GRP20 binds to poly-purine motifs in micro and small exons and a spliceosome component; both RNA binding and spliceosome interaction are important for flower development and micro-exon retention. Our results provide new insights into the mechanisms of micro-exon retention in flower development.

Mixed infection of ITPase-encoding potyvirid and secovirid in Mercurialis perennis: evidences for a convergent euphorbia-specific viral counterstrike.

BACKGROUND: In cellular organisms, inosine triphosphate pyrophosphatases (ITPases) prevent the incorporation of mutagenic deaminated purines into nucleic acids. These enzymes have also been detected in the genomes of several plant RNA viruses infecting two euphorbia species. In particular, two ipomoviruses produce replicase-associated ITPases to cope with high concentration of non-canonical nucleotides found in cassava tissues. METHOD: Using high-throughput RNA sequencing on the wild euphorbia species Mercurialis perennis, two new members of the families Potyviridae and Secoviridae were identified. Both viruses encode for a putative ITPase, and were found in mixed infection with a new partitivirid. Following biological and genomic characterization of these viruses, the origin and function of the phytoviral ITPases were investigated. RESULTS: While the potyvirid was shown to be pathogenic, the secovirid and partitivirid could not be transmitted. The secovirid was found belonging to a proposed new Comovirinae genus tentatively named "Mercomovirus", which also accommodates other viruses identified through transcriptome mining, and for which an asymptomatic pollen-associated lifestyle is suspected. Homology and phylogenetic analyses inferred that the ITPases encoded by the potyvirid and secovirid were likely acquired through independent horizontal gene transfer events, forming lineages distinct from the enzymes found in cassava ipomoviruses. Possible origins from cellular organisms are discussed for these proteins. In parallel, the endogenous ITPase of M. perennis was predicted to encode for a C-terminal nuclear localization signal, which appears to be conserved among the ITPases of euphorbias but absent in other plant families. This subcellular localization is in line with the idea that nucleic acids remain protected in the nucleus, while deaminated nucleotides accumulate in the cytoplasm where they act as antiviral molecules. CONCLUSION: Three new RNA viruses infecting M. perennis are described, two of which encoding for ITPases. These enzymes have distinct origins, and are likely required by viruses to circumvent high level of cytoplasmic non-canonical nucleotides. This putative plant defense mechanism has emerged early in the evolution of euphorbias, and seems to specifically target certain groups of RNA viruses infecting perennial hosts.

Empirical Bayes single nucleotide variant-calling for next-generation sequencing data.

One of the fundamental computational problems in cancer genomics is the identification of single nucleotide variants (SNVs) from DNA sequencing data. Many statistical models and software implementations for SNV calling have been developed in the literature, yet, they still disagree widely on real datasets. Based on an empirical Bayesian approach, we introduce a local false discovery rate (LFDR) estimator for germline SNV calling. Our approach learns model parameters without prior information, and simultaneously accounts for information across all sites in the genomic regions of interest. We also propose another LFDR-based algorithm that reliably prioritizes a given list of mutations called by any other variant-calling algorithm. We use a suite of gold-standard cell line data to compare our LFDR approach against a collection of widely used, state of the art programs. We find that our LFDR approach approximately matches or exceeds the performance of all of these programs, despite some very large differences among them. Furthermore, when prioritizing other algorithms' calls by our LFDR score, we find that by manipulating the type I-type II tradeoff we can select subsets of variant calls with minimal loss of sensitivity but dramatic increases in precision.

Sequencing by avidity enables high accuracy with low reagent consumption.

We present avidity sequencing, a sequencing chemistry that separately optimizes the processes of stepping along a DNA template and that of identifying each nucleotide within the template. Nucleotide identification uses multivalent nucleotide ligands on dye-labeled cores to form polymerase-polymer-nucleotide complexes bound to clonal copies of DNA targets. These polymer-nucleotide substrates, termed avidites, decrease the required concentration of reporting nucleotides from micromolar to nanomolar and yield negligible dissociation rates. Avidity sequencing achieves high accuracy, with 96.2% and 85.4% of base calls having an average of one error per 1,000 and 10,000 base pairs, respectively. We show that the average error rate of avidity sequencing remained stable following a long homopolymer.

DNA Encoding-Based Nucleotide Pattern and Deep Features for Instance and Class-Based Image Retrieval.

Recently, DNA encoding has shown its potential to store the vital information of the image in the form of nucleotides, namely A, C, T , and G , with the entire sequence following run-length and GC-constraint. As a result, the encoded DNA planes contain unique nucleotide strings, giving more salient image information using less storage. In this paper, the advantages of DNA encoding have been inherited to uplift the retrieval accuracy of the content-based image retrieval (CBIR) system. Initially, the most significant bit-plane-based DNA encoding scheme has been suggested to generate DNA planes from a given image. The generated DNA planes of the image efficiently capture the salient visual information in a compact form. Subsequently, the encoded DNA planes have been utilized for nucleotide patterns-based feature extraction and image retrieval. Simultaneously, the translated and amplified encoded DNA planes have also been deployed on different deep learning architectures like ResNet-50, VGG-16, VGG-19, and Inception V3 to perform classification-based image retrieval. The performance of the proposed system has been evaluated using two corals, an object, and a medical image dataset. All these datasets contain 28,200 images belonging to 134 different classes. The experimental results confirm that the proposed scheme achieves perceptible improvements compared with other state-of-the-art methods.

ALPK1 mutants causing ROSAH syndrome or Spiradenoma are activated by human nucleotide sugars.

The atypical protein kinase ALPK1 is activated by the bacterial nucleotide sugar ADP-heptose and phosphorylates TIFA to switch on a signaling pathway that combats microbial infection. In contrast, ALPK1 mutations cause two human diseases: the ALPK1[T237M] and ALPK1[Y254C] mutations underlie ROSAH syndrome (retinal dystrophy, optic nerve oedema, splenomegaly, anhidrosis, and migraine headache), while the ALPK1[V1092A] mutation accounts for 45% of spiradenoma and 30% of spiradenocarcinoma cases studied. In this study, we demonstrate that unlike wild-type (WT) ALPK1, the disease-causing ALPK1 mutants trigger the TIFA-dependent activation of an NF-κB/activator protein 1 reporter gene in the absence of ADP-heptose, which can be suppressed by either of two additional mutations in the ADP-heptose binding site that prevent the activation of WT ALPK1 by ADP-heptose. These observations are explained by our key finding that although ALPK1[T237M] and ALPK1[V1092A] are activated by bacterial ADP-heptose, they can also be activated by nucleotide sugars present in human cells (UDP-mannose, ADP-ribose, and cyclic ADP-ribose) which can be prevented by disruption of the ADP-heptose binding site. The ALPK1[V1092A] mutant was also activated by GDP-mannose, which did not activate ALPK1[T237M]. These are new examples of disease-causing mutations permitting the allosteric activation of an enzyme by endogenous molecules that the WT enzyme does not respond to. We propose that the loss of the specificity of ALPK1 for bacterial ADP-heptose underlies ROSAH syndrome and spiradenoma/spiradenocarcinoma caused by ALPK1 mutation.

Chargaff's second parity rule lies at the origin of additive genetic interactions in quantitative traits to make omnigenic selection possible.

Background: Francis Crick's central dogma provides a residue-by-residue mechanistic explanation of the flow of genetic information in living systems. However, this principle may not be sufficient for explaining how random mutations cause continuous variation of quantitative highly polygenic complex traits. Chargaff's second parity rule (CSPR), also referred to as intrastrand DNA symmetry, defined as near-exact equalities G ≈ C and A ≈ T within a single DNA strand, is a statistical property of cellular genomes. The phenomenon of intrastrand DNA symmetry was discovered more than 50 years ago; at present, it remains unclear what its biological role is, what the mechanisms are that force cellular genomes to comply strictly with CSPR, and why genomes of certain noncellular organisms have broken intrastrand DNA symmetry. The present work is aimed at studying a possible link between intrastrand DNA symmetry and the origin of genetic interactions in quantitative traits. Methods: Computational analysis of single-nucleotide polymorphisms in human and mouse populations and of nucleotide composition biases at different codon positions in bacterial and human proteomes. Results: The analysis of mutation spectra inferred from single-nucleotide polymorphisms observed in murine and human populations revealed near-exact equalities of numbers of reverse complementary mutations, indicating that random genetic variations obey CSPR. Furthermore, nucleotide compositions of coding sequences proved to be statistically interwoven via CSPR because pyrimidine bias at the 3rd codon position compensates purine bias at the 1st and 2nd positions. Conclusions: According to Fisher's infinitesimal model, we propose that accumulation of reverse complementary mutations results in a continuous phenotypic variation due to small additive effects of statistically interwoven genetic variations. Therefore, additive genetic interactions can be inferred as a statistical entanglement of nucleotide compositions of separate genetic loci. CSPR challenges the neutral theory of molecular evolution-because all random mutations participate in variation of a trait-and provides an alternative solution to Haldane's dilemma by making a gene function diffuse. We propose that CSPR is symmetry of Fisher's infinitesimal model and that genetic information can be transferred in an implicit contactless manner.

A novel single nucleotide mutation of TFL1 alters the plant architecture of Gossypium arboreum through changing the pre-mRNA splicing.

KEY MESSAGE: A single nucleotide mutation from G to A at the 201st position changed the 5' splice site and deleted 31 amino acids in the first exon of GaTFL1. Growth habit is an important agronomic trait that plays a decisive role in the plant architecture and crop yield. Cotton (Gossypium) tends to indeterminate growth, which is unsuitable for the once-over mechanical harvest system. Here, we identified a determinate growth mutant (dt1) in Gossypium arboreum by EMS mutagenesis, in which the main axis was terminated with the shoot apical meristem (SAM) converted into flowers. The map-based cloning of the dt1 locus showed a single nucleotide mutation from G to A at the 201st positions in TERMINAL FLOWER 1 (GaTFL1), which changed the alternative RNA 5' splice site and resulted in 31 amino acids deletion and loss of function of GaTFL1. Comparative transcriptomic RNA-Seq analysis identified many transporters responsible for the phytohormones, auxin, sugar, and flavonoids, which may function downstream of GaTFL1 to involve the plant architecture regulation. These findings indicate a novel alternative splicing mechanism involved in the post-transcriptional modification and TFL1 may function upstream of the auxin and sugar pathways through mediating their transport to determine the SAM fate and coordinate the vegetative and reproductive development from the SAM of the plant, which provides clues for the TFL1 mechanism in plant development regulation and provide research strategies for plant architecture improvement.

Deep learning of human polyadenylation sites at nucleotide resolution reveals molecular determinants of site usage and relevance in disease.

The genomic distribution of cleavage and polyadenylation (polyA) sites should be co-evolutionally optimized with the local gene structure. Otherwise, spurious polyadenylation can cause premature transcription termination and generate aberrant proteins. To obtain mechanistic insights into polyA site optimization across the human genome, we develop deep/machine learning models to identify genome-wide putative polyA sites at unprecedented nucleotide-level resolution and calculate their strength and usage in the genomic context. Our models quantitatively measure position-specific motif importance and their crosstalk in polyA site formation and cleavage heterogeneity. The intronic site expression is governed by the surrounding splicing landscape. The usage of alternative polyA sites in terminal exons is modulated by their relative locations and distance to downstream genes. Finally, we apply our models to reveal thousands of disease- and trait-associated genetic variants altering polyadenylation activity. Altogether, our models represent a valuable resource to dissect molecular mechanisms mediating genome-wide polyA site expression and characterize their functional roles in human diseases.

Highly Accurate Sequence- and Position-Independent Error Profiling of DNA Synthesis and Sequencing.

A comprehensive error analysis of DNA-stored data during processing, such as DNA synthesis and sequencing, is crucial for reliable DNA data storage. Both synthesis and sequencing errors depend on the sequence and the transition of bases of nucleotides; ignoring either one of the error sources leads to technical challenges in minimizing the error rate. Here, we present a methodology and toolkit that utilizes an oligonucleotide library generated from a 10-base-shifted sequence array, which is individually labeled with unique molecular identifiers, to delineate and profile DNA synthesis and sequencing errors simultaneously. This methodology enables position- and sequence-independent error profiling of both DNA synthesis and sequencing. Using this toolkit, we report base transitional errors in both synthesis and sequencing in general DNA data storage as well as degenerate-base-augmented DNA data storage. The methodology and data presented will contribute to the development of DNA sequence designs with minimal error.

Comparative mitochondrial genome analysis of three leafhopper species of the genus Abrus Dai & Zhang (Hemiptera: Cicadellidae: Deltocephalinae) from China with phylogenetic implication.

BACKGROUND: The phylogenetic position and classification of Athysanini are poorly defined, as it includes a large group of polyphyletic genera that have historically been assigned to it mainly because they still exhibit the most typical deltocephaline genitalic and external body characters but lack the distinctive characteristics that other tribes possess. The bamboo-feeding leafhopper genus Abrus belong to the tribe Athysanini of subfamily Deltocephalinae, which currently comprises 19 valid described species, and are limited to the Oriental and Palaearctic regions in China. Although the taxonomy of Abrus are well updated, the references on comparative mitogenomic analyses of Abrus species are only known for a single species. In this study, we sequenced and analyzed the complete mitochondrial genomes (mitogenomes) of Abrus daozhenensis Chen, Yang & Li, 2012 (16,391bp) and A. yunshanensis Chen, Yang & Li, 2012 (15,768bp) (Athysanini), and compared with published mitogenome sequence of A. expansivus Xing & Li, 2014 (15,904bp). RESULTS: These Abrus species shared highly conserved mitogenomes with similar gene order to that of the putative ancestral insect with 37 typical genes and a non-coding A + T-rich region. The nucleotide composition of these genomes is highly biased toward A + T nucleotides (76.2%, 76.3%, and 74.7%), AT-skews (0.091 to 0.095, and 0.095), negative GC-skews (- 0.138, - 0.161, and - 0.138), and codon usage. All 22 tRNA genes had typical cloverleaf secondary structures, except for trnS1 (AGN) which lacks the dihydrouridine arm, and distinctively trnG in the mitogenome of A. expansivus lacks the TψC arm. Phylogenetic analyses based on 13 PCGs, 2 rRNA genes, and 22 tRNA genes consistently recovered the monophyletic Opsiini, Penthimiini, Selenocephalini, Scaphoideini, and Athysanini (except Watanabella graminea, previously sequenced species as Chlorotettix nigromaculatus) based on limited available mitogenome sequence data of 37 species. CONCLUSION: At present, Abrus belongs to the tribe Athysanini based on both morphological and molecular datasets, which is strongly supported in present phylogenetic analyses in both BI and ML methods using the six concatenated datasets: amino acid sequences and nucleotides from different combinations of protein-coding genes (PCGs), ribosomal RNA (rRNAs), and transfer RNA (tRNAs). Phylogenetic trees reconstructed herein based on the BI and ML analyses consistently recovered monophylitic Athysanini, except Watanabella graminea (Athysanini) in Opsiini with high support values.

An Analysis of Nucleotide-Amyloid Interactions Reveals Selective Binding to Codon-Sized RNA.

Interactions between RNA and proteins are the cornerstone of many important biological processes from transcription and translation to gene regulation, yet little is known about the ancient origin of said interactions. We hypothesized that peptide amyloids played a role in the origin of life and that their repetitive structure lends itself to building interfaces with other polymers through avidity. Here, we report that short RNA with a minimum length of three nucleotides binds in a sequence-dependent manner to peptide amyloids. The 3'-5' linked RNA backbone appears to be well-suited to support these interactions, with the phosphodiester backbone and nucleobases both contributing to the affinity. Sequence-specific RNA-peptide interactions of the kind identified here may provide a path to understanding one of the great mysteries rooted in the origin of life: the origin of the genetic code.

Nanopore-based direct sequencing of RNA transcripts with 10 different modified nucleotides reveals gaps in existing technology.

RNA undergoes complex posttranscriptional processing including chemical modifications of the nucleotides. The resultant-modified nucleotides are an integral part of RNA sequences that must be considered in studying the biology of RNA and in the design of RNA therapeutics. However, the current "RNA-sequencing" methods primarily sequence complementary DNA rather than RNA itself, which means that the modifications present in RNA are not captured in the sequencing results. Emerging direct RNA-sequencing technologies, such as those offered by Oxford Nanopore, aim to address this limitation. In this study, we synthesized and used Nanopore technology to sequence RNA transcripts consisting of canonical nucleotides and 10 different modifications in various concentrations. The results show that direct RNA sequencing still has a baseline error rate of >10%, and although some modifications can be detected, many remain unidentified. Thus, there is a need to develop sequencing technologies and analysis methods that can comprehensively capture the total complexity of RNA. The RNA sequences obtained through this project are made available for benchmarking analysis methods.

New algorithm for the analysis of nucleotide and amino acid evolutionary relationships based on Klein four-group.

Phylogenetics is the study of ancestral relationships among biological species. Such sequence analyses are often represented as phylogenetic trees. The branching pattern of each tree and its topology reflect the evolutionary relatedness between analyzed sequences. We present a Klein four-group algorithm (K4A) for the evolutionary analysis of nucleotide and amino acid sequences. Klein four-group set of operators consists of: identity e (U), and three elements-a = transition (C), b = transversion (G) and c = transition-transversion or complementarity (A). We generated Klein four-group based distance matrices of: 1. Cayley table (CK4), 2. Table rows (K4R), 3. Table columns (K4C), and 4. Euclidean 2D distance (K4E). The performance of the matrices was tested on a dataset of RecA proteins in bacteria, eukaryotes (Rad51 homolog) and archaea (RadA homolog). RecA and its functional homologs are found in all species, and are essential for the repair and maintenance of DNA. Consequently, they represent a good model for the study of evolutionary relationship of protein and nucleotide sequences. The ancestral relationship between the sequences was correctly classified by all K4A matrices concerning general topology. All distance matrices exhibited small variations among species, and overall results of tree classification were in agreement with the general patterns obtained by standard BLOSUM and PAM substitution matrices. During the evolution of a code there is a phase of optimization of system rules, the ambiguity of a code is eliminated, and the system starts producing specific components. Klein four-group algorithm is consistent with the concept of ambiguity reduction. It also enables the use of different genetic code table variants optimized for particular transitions in evolution based on biological specificity.

Development of nucleotide signatures for common poisonous organisms provides a new strategy for food poisoning diagnosis.

DNA barcoding is widely used in toxic species authentication, but due to serious DNA degradation of forensic materials, the application of full-length barcode sequences in food poisoning diagnosis is greatly limited. Nucleotide signature, a shorter specific molecular marker, derived from traditional DNA barcoding has been proposed as an emerging tool of toxic species detection in deeply processed materials. In this study, to resolve the frequent food poisoning accidents with unknown origin, we envisioned developing a nucleotide signature data set of common poisonous organisms and combining high-throughput sequencing (HTS) to reveal the poisoning cause. Ninety-three individuals and 1093 DNA barcode sequences of twelve common poisonous plants, fish, mushrooms and their related species were collected. Through sequence alignment and screening, the nucleotide signatures were respectively developed and validated as their specific molecular markers. The sequence length varied from 19 bp to 38 bp. These fragments were conserved within the same species or genera, and the specificity between related species has been also demonstrated. To further evaluate the application potential of nucleotide signature in forensic diagnosis, simulated forensic specimens (SFS) containing different poisonous ingredients were sequenced by HTS with PCR-free libraries. As a result, the nucleotide signature was successfully captured from original HTS data without assembly and annotation, accompanied by a high detection sensitivity of 0.1 ng/µl in mixture system. Therefore, this method was suitable for the assay of forensic materials with serious DNA degradation. The present study undoubtedly provides a new perspective and strong support for the detection of toxic ingredients and the diagnosis of food poisoning.