Environmental surveillance and clinical diagnostics heavily rely on the polymerase chain reaction (PCR) for target detection. A growing list of microbial threats warrants new PCR-based detection methods that are highly sensitive, specific, and multiplexable. Here, we introduce a PCR-based icosaplex (20-plex) assay for detecting 18 enteropathogen and two antimicrobial resistance genes. This multiplexed PCR assay leverages the self-avoiding molecular recognition system (SAMRS) to avoid primer dimer formation, the artificially expanded genetic information system (AEGIS) for amplification specificity, and next-generation sequencing for amplicon identification. Using parallelized multitarget TaqMan Array Cards (TAC) to benchmark performance of the 20-plex assay on wastewater, soil, and human stool samples, we found 90% agreement on positive calls and 89% agreement on negative calls. Additionally, we show how long-read and short-read sequencing information from the 20-plex can be used to further classify allelic variants of genes and distinguish subspecies. The strategy presented offers sensitive, affordable, and robust multiplex detection that can be used to support efforts in wastewater-based epidemiology, environmental monitoring, and human/animal diagnostics.

By rearranging hydrogen bond donor and acceptor groups within a standard Watson–Crick geometry, DNA can add eight independently replicable nucleotides forming four additional not found in standard Terran DNA. For many applications, the orthogonal pairing of standard and nonstandard pairs offers a key advantage. However, other applications require standard and nonstandard nucleotides to communicate with each other. This is especially true when seeking to recruit high-throughput instruments (e.g., Illumina), designed to sequence standard 4-nucleotide DNA, to sequence DNA that includes added nucleotides. For this purpose, PCR workflows are needed to replace nonstandard nucleotides in (for example) a 6-letter DNA sequence by defined mixtures of standard nucleotides built from 4 nucleotides. High-throughput sequencing can then report the sequences of those mixtures to bioinformatic alignment tools, which infer the original 6-nucleotide sequence by analysis of the mixtures. Unfortunately, the intrinsic orthogonality of standard and nonstandard nucleotides often demand polymerases that violate pairing biophysics to do this replacement, leading to inefficiencies in this “transliteration” process. Thus, laboratory in vitro evolution (LIVE) using “anthropogenic evolvable genetic information systems” (AEGIS), an important “consumer” of new sequencing tools, has been slow to be democratized; robust sequencing is needed to identify the AegisBodies and AegisZymes that AEGIS-LIVE delivers. This work introduces a new way to connect synthetic and standard molecular biology: biversal nucleotides. In an example presented here, a pyrimidine analogue (pyridine-2-one, y) pairs with Watson–Crick geometry to both a nonstandard base (2-amino-8-imidazo-[1,2a]-1,3,5-triazin-[8H]-4-one, P, the Watson–Crick partner of 6-amino-5-nitro-[1H]-pyridin-2-one, Z) and a base that completes the Watson–Crick hydrogen bond pattern (2-amino-2'-deoxyadenosine, amA). PCR amplification of GACTZP DNA with dyTP delivers products where Z:P pairs are cleanly transliterated to A:T pairs. In parallel, PCR of the same GACTZP sample at higher pH delivers products where Z:P pairs are cleanly transliterated to C:G pairs. By allowing robust sequencing of 6-letter GACTZP DNA, this workflow will help democratize AEGIS-LIVE. Further, other implementations of the biversal concept can enable communication across and between standard DNA and synthetic DNA more generally.

Adding synthetic nucleotides to DNA increases the linear information density of DNA molecules. Here we report that it also can increase the diversity of their three-dimensional folds. Specifically, an additional nucleotide (dZ, with a 5-nitro-6-aminopyridone nucleobase), placed at twelve sites in a 23-nucleotides-long DNA strand, creates a fairly stable unimolecular structure (that is, the folded Z-motif, or fZ-motif) that melts at 66.5°C at pH 8.5. Spectroscopic, gel and two-dimensional NMR analyses show that the folded Z-motif is held together by six reverse skinny dZ^-:dZ base pairs, analogous to the crystal structure of the free heterocycle. Fluorescence tagging shows that the dZ^-:dZ pairs join parallel strands in a four-stranded compact down-up-down-up fold. These have two possible structures: one with intercalated dZ^-:dZ base pairs, the second without intercalation. The intercalated structure would resemble the i-motif formed by dC:dC⁺-reversed pairing at pH ≤ 6.5. This fZ-motif may therefore help DNA form compact structures needed for binding and catalysis.

With just four building blocks, low sequence information density, few functional groups, poor control over folding, and difficulties in forming compact folds, natural DNA and RNA have been disappointing platforms from which to evolve receptors, ligands, and catalysts. Accordingly, synthetic biology has created "artificially expanded genetic information systems" (AEGIS) to add nucleotides, functionality, and information density. With the expected improvements seen in AegisBodies and AegisZymes, the task for synthetic biologists shifts to developing for expanded DNA the same analytical tools available to natural DNA. Here we report one of these, an enzyme-assisted sequencing of expanded genetic alphabet (ESEGA) method to sequence six-letter AEGIS DNA. We show how ESEGA analyses this DNA at single base resolution, and applies it to optimized conditions for six-nucleotide PCR, assessing the fidelity of various DNA polymerases, and extending this to AEGIS components with functional groups. This supports the renewed exploitation of expanded DNA alphabets in biotechnology.