"The skepticism surrounding eukaryotic dual coding is unwarranted: rather than being artifacts, overlapping reading frames are often hallmarks of fascinating biology." - Wen-Yu Chung, Anton Nekrutenko and colleagues, PLoS Computational Biology, 20071
Summary
A stretch of DNA can carry more than one message. Read from a different starting letter, the same letters spell a second protein in a second "reading frame." Read on the opposite strand, they can spell a third. Inside a protein-coding gene the same letters can also carry signals for splicing, for proteins that bind DNA, and for the folding of the RNA copy. Design writers point to this as information compression that only a mind could achieve. John Sanford writes that "probably all DNA sequences in the genome encrypt multiple codes (up to 12)," and that "this type of data compression can only result from the highest level of information design."2 Evolutionary biologists answer that a second frame is usually born inside an older gene by mutation and selection, a process called overprinting, and that "proteins may arise continuously de novo."3 This article counts how much DNA carries two protein messages, prices what that costs, and asks which account fits the pattern.
The short answer: two-protein coding is real, common in the smallest genomes and rare among the annotated genes of ours. Where it occurs it couples two proteins, so most changes that alter one also alter the other. Design writers are right about that coupling. Their larger claims, that most DNA carries many codes and that beneficial changes in it must be vanishingly rare, rest on assumptions, not on measurements found here. Comparisons between related species support overprinting as a route to a second protein frame, with origins inferred rather than observed.
-
The genetic code ties two overlapping frames together. Under the standard code, about a quarter of single-letter changes in a gene leave its protein sequence unchanged: 24.6% at the base composition of human genes. Where two frames overlap on the same strand, only 1.2% of changes leave both protein sequences unchanged. About 46% change both, and 8% put a stop codon into one of them.4 About 23% change one protein and leave the other's sequence alone, and another 23% the reverse. This is arithmetic on the code, not a fitness measurement. A change that leaves both proteins alone can still matter, a change to both is not always harmful, and how many of the one-protein changes would help is not known.
-
Small viruses are full of overlaps. Large genomes are not. Counted for this article from the reference sequences, as a share of each genome's whole length, hepatitis B virus reads 50.0% of its 3,182 letters in two frames, the phage phiX174 15.1% and HIV-1 14.7%.5 Of 5,976 reference virus genomes, 53% carry an overlap of 50 letters or more.6 The bacterium E. coli reads 0.2% of its genome in two frames, most of it in overlaps of one to four letters where one gene's stop codon shares letters with the next gene's start. Human mitochondrial DNA reads 0.3% in two frames.5
-
Among annotated human genes, two proteins from one stretch are uncommon. Taking one standard ("canonical") transcript per gene from GENCODE release 50, 21,568 of 33.9 million protein-coding letters, 0.064%, are read in two frames. This count includes the 13 mitochondrial protein-coding genes, which add two overlapping pairs and a few dozen letters. Only 28 pairs of genes overlap on the same strand in different frames, and 16 pairs overlap protein-coding sequence on opposite strands.7 That is one annotation choice, not a census of every working second frame. Counting every annotated transcript raises the share to 7.2%, almost all of it from alternative transcripts that read a stretch of one gene in a second frame. Those transcript models vary in how well they are supported, and the annotation does not say how many make a stable protein. A 2007 scan found 40 human genes with a second frame of 500 letters or more conserved in mouse, 15 of them high-confidence.1 A handful of human genes have the second protein shown in cells or tissue: CDKN2A (p16 and ARF), GNAS (XLalphas and ALEX), POLG (POLGARF), PRNP (AltPrP) and ATXN1 (Alt-ATXN1).
-
RNA from both strands is common. Proteins from both strands are not. Counted for this article, 55.3% of human protein-coding genes have an exon of another gene, coding or not, overlapping one of their exons on the opposite strand. Protein-coding sequence overlaps protein-coding sequence on opposite strands at only 16 gene pairs.7 Ribosome profiling, which catches ribosomes sitting on RNA, has catalogued 7,264 translated reading frames that are not annotated proteins; 1,469 of them overlap a known gene's coding sequence in another frame.8 When 95,520 proteomics experiments were searched for their products, 24.6% of the 7,264 left a detectable peptide, almost all of them fragments displayed to the immune system, and 15 met the strict standard for a new protein.8 Translation, a detected peptide and a working protein are three different findings.
-
The constraint shows up in real genes. Across 117 overlapping RNA-virus genes, the overlapping parts had shorter evolutionary trees and fewer silent substitutions than the rest of the same genes, while amino-acid-changing substitutions did not differ significantly.99 That is the code arithmetic seen in real sequences. The same paper's ratio of the two rates cannot be relied on: its Table 1 puts the ratio lower in the overlaps, while its text and supplementary table put it higher, and the usual ratio methods are known to misread overlapping genes.1011 In hepatitis B virus, drug-resistance mutations in the polymerase gene change the surface protein read from the same letters.12 A phiX174 genome redesigned with its overlaps pulled apart grows only when a host plasmid supplies one of its genes, and in one laboratory test it released 16 phages per infected cell against 145 for the natural virus.1314 That gap mixes the effect of separating the genes with the redesign's other changes.
-
Where they can be dated, second frames look like later additions. In viruses, the newer frame of a pair, identified by comparing relatives, is usually found in one lineage only, codes for floppier protein, and uses codons unlike the rest of the genome.151516 HIV-1's antisense protein frame is open in most of the pandemic group M and not in the other groups or related ape viruses.17 Of the 7,264 human translated frames, 89.6% have no detectable protein counterpart in mammals outside the primates, a comparative finding that short proteins make hard to test.18 POLGARF traces to a mobile element that landed in the gene before placental mammals spread.19 All these ages are inferred by comparing living species. None is a birth watched as it happened, and these features are consistent with overprinting without proving it.
-
The design-side claims split. The claim that overlap constrains change is right, and the data in point 5 support it. The claim that most DNA carries many overlapping codes is not measured in any source reviewed here. Sanford's "over half the human genome is read backwards as well as forwards" carries no citation,2 and his "up to 12" codes come from a 1989 paper whose body could not be checked here. The 2013 calculation by George Montañez, Robert Marks, Jorge Fernandez and Sanford is correct arithmetic on assumed inputs: in its headline case, 12 codes per letter, 99.9% of letters already optimal, preferred letters chosen independently in each code, and every code site defined as one where change matters. Its own authors write that "the exact degree to which nucleotides participate in two or more codes is still unknown."20 Its result is a conditional probability, not a measured share of beneficial mutations. No design writer found here engages the overprinting evidence.
Verdict. Overlapping protein codes are a real feature of compact genomes and an uncommon one among the annotated genes of large genomes. Where they exist they couple two messages and hold letters in place, which is the design writers' point, and it is correct. How often a change in an overlap helps has not been measured in any source reviewed here, so the coupling does not show that helpful changes are vanishingly rare. Two protein frames on one stretch are not a general property of DNA. Where comparisons between species allow them to be dated, second frames look lineage-limited and built inside older genes, which is what overprinting predicts. That pattern supports overprinting as an explanation for some second protein frames. It does not measure the extent of other regulatory layers throughout the genome, settle where the first genes came from, or establish the age of the earth.
Part 1: Reading one stretch of DNA two ways
What "overlapping" means here
A protein-coding gene is read three letters at a time. Each three-letter word, a codon, names one amino acid or says "stop." Start one letter later and the same letters fall into different words. That is a second reading frame, often called +1, and starting two letters later gives +2. The opposite strand of the double helix offers three more frames, read in the other direction. Two genes "overlap" when the same letters are read in two of these frames and both readings make a product.
Three other uses of the word need to be kept apart, because the design-side literature mixes them.
- Two genes whose RNA copies overlap, usually on opposite strands, with no claim that both make protein. A protein gene and a non-coding RNA gene facing each other is the common case.
- One gene read in two frames in different messages. Alternative splicing can drop or add a few letters, so that downstream letters are read in another frame in one version of the message. Both versions come from the same gene.
- Other signals laid over coding letters: splice sites, places where regulatory proteins bind, RNA folding, nucleosome positioning. These are often called "codes." Most are preferences rather than codes in the strict sense: a binding site tolerates many sequences, and a nucleosome prefers some letters to others without requiring one.
Only overlap in the strict sense, two protein frames on the same letters, is counted letter by letter in this article. The others are counted where a source counts them, and named for what they are.
What counts as evidence that a second frame is used
The evidence for a second frame comes in steps, and each step says less than the one after it.
- An open frame. A stretch with no stop codon in the second frame. Long open second frames are common in ordinary human genes, as Part 5 shows, so an open frame alone proves little.
- RNA. The stretch is copied into RNA. This shows the letters are copied, not that the second frame is read.
- Translation. Ribosome profiling catches ribosomes sitting on RNA, and a three-letter rhythm in where they sit shows which frame they read. This shows that a ribosome reads the frame. It does not show that the product survives.
- A detected peptide. Mass spectrometry finds a fragment of the second protein in cell extracts. Fragments found only among the peptides cells display to the immune system show the protein was made and chopped up, not that it lasts.
- A stable protein. An antibody against the second protein finds it in normal cells, at the expected size.
- A function. Removing the second protein, while leaving the first one's amino acids unchanged, harms the cell or the organism, and putting it back repairs the harm.
- Conservation. The second frame stays open across related species. This counts as evidence of function only if common descent is assumed, and it is sometimes confounded by the first frame's own constraint.
Large catalogues sit on steps 1 to 4. The handful of proven human cases reach steps 5 and 6. A number quoted without its step is not a number about function.
What the genetic code charges for a second frame
Every single-letter change in a gene does one of three things to its protein: nothing, when the new codon names the same amino acid; swaps one amino acid; or puts in a stop. The standard code is built so that about a quarter of changes do nothing. Most of these silent changes fall on the third letter of a codon.
In a same-strand overlap, the third letter of a codon in one frame is the first or second letter of a codon in the other. A change that is silent in one frame is therefore usually an amino-acid change in the other. The table counts every single-letter change in every codon context, weighted by the base composition of human coding DNA, for each way two frames can lie over each other.
| How the two frames lie | Silent in frame A | Silent in frame B | Silent in both | A stop in either | Amino acid changed in both |
|---|---|---|---|---|---|
| One frame only, no overlap | 24.6% | 24.6% | 4.1% | ||
| Same strand, shifted 1 or 2 letters | 24.6% | 24.6% | 1.2% | 8.0% | 46.4% |
| Opposite strand, codons aligned | 24.2% | 24.6% | 2.4% | 7.1% | 49.2% |
| Opposite strand, shifted 1 letter | 24.3% | 24.8% | 0.0% | 8.3% | 44.7% |
| Opposite strand, shifted 2 letters | 24.4% | 24.9% | 17.2% | 7.1% | 62.4% |
The same-strand row is the one that matters most, because it is the commonest kind of long overlap. There, only 1.2% of changes leave both protein sequences alone, against 24.6% for a gene read one way. About 23% are silent in one frame and change the other. Wen-Yu Chung and colleagues got the same shape in 2007 by the same method, averaged over codon positions: about 1.2% silent in both and about 46% changing both.1 One opposite-strand arrangement is worse than the same-strand case: with codons shifted by one letter, no change at all is silent in both frames. Another is much milder: when the third letters of both codons fall on the same base pair, 17% of changes are silent in both.4
Three limits keep this arithmetic in its place.
- It counts changes to protein sequence, not changes to fitness. A change that is silent in both frames can still alter splicing, RNA folding or how fast the message is read. The article Synonymous Mutations and Codon Usage measures how often that matters. A change that swaps an amino acid in both proteins is not automatically harmful either; most proteins tolerate many single swaps. How often a real change helps, harms or does nothing is measured in Beneficial and Harmful Mutations, and this article does not repeat it.
- It treats letters as independent. Real genes have codon preferences and neighbour effects. The shares shift a little with base composition (24.4% silent at even composition, 24.6% at human composition), and not by enough to change the picture.
- It prices one overlap of two frames. It says nothing about how many other layers sit on the same letters. That is a separate, empirical question, and Parts 2, 3 and 6 take it up.
What the table does show is that a second protein frame removes almost all of the changes that leave both protein sequences alone, and ties most changes in one protein to a change in the other. A change that alters protein A while leaving protein B's sequence alone is still available, about 23% of changes in each direction, but how many of those would help is a question the code cannot answer. Any account of overlapping genes has to live with that coupling.
Part 2: How common overlapping genes are
Our counts, from the reference sequences
A simple measure is the number of letters read in two different protein frames. It needs the annotation of every coding sequence and its frame, and it must count each letter once, however many products use it. The table gives that count for six reference genomes, counted for this article from the NCBI RefSeq records, and for human protein-coding genes from GENCODE release 50.57 The denominators differ. For the viruses, E. coli and the mitochondrion the share is of the whole genome, most of which codes for protein. For the human row it is of protein-coding letters only, because only about 1% of the human genome, 33.9 million of 3.1 billion letters, codes for protein; as a share of the whole 3.1-billion-letter genome the same count is 0.0007%.
| Genome | Length | Letters read in two frames | Where |
|---|---|---|---|
| Hepatitis B virus, strain ayw | 3,182 | 1,591 (50.0%) | polymerase over surface 1,170; polymerase over X 250; core over polymerase 146; core over X 25 |
| Phage phiX174 | 5,386 | 814 (15.1%) | B inside A, K across A and C, E inside D, and short start-stop overlaps |
| HIV-1, reference strain HXB2 | 9,181 | 1,346 (14.7%) | includes 570 letters where the antisense protein ASP lies on the opposite strand of env |
| SARS-CoV-2, reference record | 29,903 | 20 (0.1%) | the record omits ORF9b, ORF3d and other frames found later |
| Human mitochondrion | 16,569 | 54 (0.3%) | ATP8 over ATP6 46; ND4L over ND4 7; ATP6 over COX3 1 |
| E. coli K-12 | 4,641,652 | 8,694 (0.2%) | 738 gene pairs, 470 of them overlapping by 1 to 4 letters, 60 by 30 or more, 5 by 150 or more |
| Human protein-coding DNA, nuclear and mitochondrial, one canonical transcript per gene | 33.9 million coding letters | 21,568 (0.064% of coding letters, not of the genome) | 28 same-strand pairs in different frames, two of them mitochondrial; 16 opposite-strand pairs |
These counts depend on what the annotators have accepted, and the error runs both ways. The SARS-CoV-2 record shows one direction: its 0.1% leaves out frames that ribosome profiling and patient antibodies later showed to be used, such as ORF3d.21 The other direction is that some annotated coding sequences, especially recent and computational models, may not make a working protein. So a count is a lower bound on real two-frame coding only if every annotated sequence in it is real, and it is not a census of open frames. Comparisons between rows show the broad pattern, small viruses against large genomes, but because the denominators differ the percentages are not measurements of one identical quantity. The phiX174 figure shows how definitions change a number. Bradley Wright and colleagues give the overlapped share of phiX174 as 16.8% in their text and "∼26% of nucleotides involved in encoding more than one gene" in their abstract.14 Counting letters shared by two annotated genes in any frame, including A* read in the same frame as A, gives 25.8% from the same record; counting only different frames gives our 15.1%.5
Viruses
Three surveys measured how many virus species have overlaps, each on a different sample and threshold. Nicola Chirico and colleagues, 2010: "75 per cent of the approximately 2000 known virus species have at least some gene overlap."22 Robert Belshaw and colleagues, 2007, on 701 RNA virus genomes: "56% of the viruses having some gene overlap."23 Timothy Schlub and Edward Holmes, 2020, on 5,976 reference genomes: "over half the instances of gene overlap are very small, covering <10 nt, and 84 per cent are <50 nt in length. Despite this, 53 per cent of all viruses still contained a gene overlap of 50 nt or larger."6 The share with a 50-letter overlap ran from 19% of double-stranded RNA viruses to 43% of positive-strand RNA viruses and higher in DNA viruses.6
How much of a virus genome is overlap depends on its size. The share of the genome in overlap falls as genomes get longer, in Belshaw's sample, in Chirico's and in Nadav Brandes and Michal Linial's.232224 Brandes and Linial put a ceiling on the absolute amount: "although viruses span three orders of magnitude in genome length, they almost never have over 1500 overlapping nucleotides, or over four significantly overlapping genes per virus."24 Measured as a share of genome length, Chirico found 4.6% for RNA viruses and 6.6% for DNA viruses of RNA-virus size, and corrected Belshaw's lower figure: "A lower value for overlap proportion in RNA viruses published previously by one of us ( Belshaw et al . 2007 ) was an error."22
Why small viruses overlap so much is disputed. One idea is that a capsid of fixed size limits genome length, so new genes have to fit inside old ones; Chirico found the size relation only in families with rigid icosahedral capsids, and Brandes and Linial found no such difference.2224 A second idea is that high mutation rates favour compact genomes; Chirico found "no obvious relationship between overlap proportion and mutation rate (linear regression, p = 0.50, n = 13)," on a small sample.22 A third, which Belshaw's title already names, is that overlaps are how small viruses gain new genes at all. Part 5 takes that up.
Bacteria
In bacteria, overlapping genes are common and almost all tiny. Zackary Johnson and Sallie Chisholm counted 64,989 overlapping gene pairs in the microbial genomes sequenced by 2003, about one for every six genes, and found that ">85% are <30 bp."25 Most are a stop codon of one gene sharing letters with the start of the next. Long embedded overlaps are rare. Before their own find, Lea Fellner and colleagues knew of "only three overlapping gene pairs" in the "extremely well-researched" E. coli.26 Ribosome profiling has since found many candidate starts inside genes: Sezen Meydan and colleagues found 124 internal start sites shared by two E. coli strains, 74 of them out of frame with the main gene, and tested one of the out-of-frame products in living cells.27 A 2020 review by Zachary Ardern and colleagues says "hundreds of short open reading frames antisense to annotated genes show evidence of both transcription and translation" across bacteria, while its own counts of tested cases are in single and double digits.28
The pattern
Overlap is densest where genomes are smallest. Hepatitis B virus reads half its letters twice. E. coli and the human mitochondrion read well under 1% of their genomes twice. Human coding DNA, on one canonical transcript per gene, reads 0.064% of its coding letters twice, and the whole human genome far less. Any account of overlapping codes has to explain why annotated two-frame coding is concentrated in the smallest genomes and sparse in the largest.
Part 3: Overlapping codes in the human genome
The counts, by what they count
Numbers about human overlap range from a few hundredths of a percent to "over half," because they count different things. The table puts them side by side with their denominators.
| What is counted | Result | Source |
|---|---|---|
| Coding letters read in two frames, one canonical transcript per gene, mitochondrial genes included | 21,568 of 33.9 million (0.064%) | GENCODE 50, this article7 |
| The same, MANE Select transcripts only | 2,626 (0.008%) | GENCODE 50, this article |
| The same, every annotated transcript | 2.80 million of 38.7 million (7.2%) | GENCODE 50, this article |
| Genes whose own transcripts read a stretch in two frames | 12,569 (every transcript); 8,226 (full-length coding transcripts); 1,346 (transcripts matching a CCDS entry) | GENCODE 50, this article |
| Gene pairs whose coding sequences overlap in different frames | 28 same strand; 16 opposite strands | GENCODE 50, this article |
| Protein-coding genes with an exon of another gene on the opposite strand | 11,124 of 20,107 (55.3%) | GENCODE 50, this article |
| Alternatively spliced genes with coding regions in two or three frames | 179 of 2,585 (6.9%) | Liang and Landweber 200629 |
| Genes with a second frame of 500 or more letters open in human and mouse | 40 candidates, 15 high-confidence | Chung 20071 |
| Human genes in an overlapping pair, among 13,484 human-mouse orthologs | 1,219 (9.0%), 91.9% of pairs on opposite strands, 51 pairs sharing coding sequence | Sanna 200830 |
| Translated frames found by ribosome profiling, not annotated as proteins | 7,264, of which 1,469 overlap a coding sequence in another frame | Deutsch 20268 |
Three things stand out.
- Overlapping RNA is common, overlapping protein is not. Half of protein-coding genes have another gene's exon facing them on the opposite strand. Only 16 gene pairs overlap coding sequence on opposite strands, and only 28 on the same strand in different frames. About half of those 28 involve genes added to GENCODE in the last few years, many from ribosome profiling, and two are in the mitochondrion.7
- Most of the 7.2% comes from alternative transcripts, and the count depends on which transcript models are admitted. Counting every transcript, 12,569 genes have a stretch read in two frames by their own transcripts. Limited to full-length protein-coding transcripts, the figure is 8,226 genes; limited to transcripts whose coding sequence matches the curated CCDS set, 1,346 genes, of which 261 have a two-frame stretch of 150 letters or more.7 CCDS is a set of coding sequences that several annotation groups agree on. A transcript outside it is not thereby shown to be unused, and one inside it is not thereby shown to make a stable second protein. The annotation tiers sort models by how they were built and agreed, not by tested function. Han Liang and Laura Landweber's 2006 count on a smaller curated set found two-frame coding in 6.9% of alternatively spliced genes, and gave its share of exons as under 0.5%: "The space of dual or multiple coding regions in the human genome is <0.5% of exons."29
- The strict counts and the generous counts differ by a factor of about a hundred, and every step up comes from admitting more transcript models. None of them approaches "most DNA." These are counts of protein frames. They do not measure RNA-binding sites, splicing signals or the other proposed layers of Part 6, which are counted separately where a source counts them.
Translated, detected, working: the ribosome profiling catalogue
Ribosome profiling has made the human count of translated frames much larger. In 2022 Jonathan Mudge and an international group of annotators pooled several studies into a standard set of 7,264 translated reading frames that are not annotated proteins.31 The 2026 Nature paper by Eric Deutsch and the same consortium gives the categories: 3,083 upstream frames that end before the main gene starts, 688 that start upstream and run into the main coding sequence in another frame, 720 inside the coding sequence in another frame, 61 that overlap its end, 504 downstream, 1,917 in RNAs annotated as non-coding, and 291 in other transcripts.8 So 1,469 of the 7,264 sit at least partly over a known protein in another frame.
Translation is step 3 on the ladder in Part 1. Deutsch and colleagues then searched 95,520 mass-spectrometry experiments for the products. The steps fall off fast.
- In conventional protein digests: "484 peptides passing FDR thresholds that map to 183 out of the 7,264 ncORFs (around 2.5%)."8
- In peptides displayed to the immune system: "3,116 peptides mapping to 1,785 out of 7,264 Ribo-seq ncORFs (24.6%)," almost all on the class I molecules that display fragments of proteins made inside the cell.8
- Meeting the Human Proteome Project's standard for a new protein, after manual review: "we identified 20 candidates for tier 1A status, but further scrutiny reduced this list to 15 ncORFs."8
The authors give a reason the strict count is low that has nothing to do with function: small proteins yield few peptides. Only "2 out of 36 known proteins under 50 amino acids (5.6%)" pass the same benchmark.8 So 15 is a floor for stable small proteins, not a ceiling. The immune-display count shows the frames are translated and the products made and cut up. It does not show that the products last or do anything.
On age, the same set looks young. Clara Sandmann, Jorge Ruiz-Orera and colleagues found that "most of the putative 7,264 sORF-encoded microproteins (n = 6,506; 89.6%) lack significant protein homology across non-primate mammals and could therefore be classified as evolutionarily young."18 Deutsch's group measured something different, whether the start, the stop and the openness of each frame are kept across species, and found that 2,211 frames (30.4%) scored in the top tenth of untranslated frames of the same type and length, against the 10% expected by chance.8 Frames that overlap a known gene score high on this measure, which the authors say is "presumably due to constraint to preserve the CDS."8 An internal frame can stay open because the main gene's letters are held in place, which is one reason conservation alone cannot show that a second frame is used.
Function is tested by knocking frames out. Jin Chen, Jonathan Weissman and colleagues knocked out about 2,350 non-canonical frames by gene editing in stem cells and leukemia cells and found "over 500 ORF knockout 'hits'," including 78 frames that overlap a main gene's start in another frame.32 The authors note that some upstream-frame effects may act by raising the main protein's level rather than through the small protein itself.32 John Prensner and colleagues tested 553 candidates, mostly upstream, downstream and non-coding-RNA frames rather than internal ones: "57 induced viability defects when knocked out in human cancer cell lines."33 These are the best functional counts available, and few of them are internal second frames.
The cases with proven second proteins
Five human genes have a second protein, read from letters shared with the first in another frame, detected in cells or tissue and studied in its own right.
| Gene | First product | Second product and frame | What is shown | How widely shared |
|---|---|---|---|---|
| CDKN2A | p16INK4a | ARF, from a separate first exon spliced into exon 2 and read in another frame | Both proteins stop cell division; in the mouse, "two unrelated proteins capable of inducing cell cycle arrest"34 | Mammals. Chickens have no p16INK4a, and their ARF stops at the end of its own exon35 |
| GNAS | XLalphas | ALEX, in the XL exon | Detected by antibody in neuroendocrine cells; binds XLalphas36 | Mammals, by the authors' alignments37 |
| POLG | Mitochondrial DNA polymerase subunit | POLGARF, 260 codons, from a CUG start in the leader, read in the -1 frame | Ribosome profiling, peptides, a binding partner (C1QBP), conserved by sequence across placental mammals19 | Placental mammals; arose from an inserted MIR repeat19 |
| PRNP | Prion protein | AltPrP, +3 frame | Antibody detection in human cells, brain and blood cells; removed by silencing the shared RNA38 | Expressed from human, bovine, sheep and deer cDNAs38 |
| ATXN1 | Ataxin-1 | Alt-ATXN1, +3 frame, letters 30 to 587 | Antibody detection in human cerebellum; binds ataxin-139 | Not established in the paper's abstract |
In 1995 Quelle, Sherr and colleagues called the INK4a arrangement "practically without precedent in mammalian genomes."34 Martin Klemke, Wieland Huttner and colleagues said the same of ALEX in 2001: "Our results reveal a mechanism of gene usage that is without precedent in mammalian genomes."36 As late as 2005 Anton Nekrutenko and colleagues could list only three other mammalian cases with both products characterized: INK4a, X-box binding protein 1 and a region shared by 4E-BP3 and MASK.37 The list has grown since, mostly through ribosome profiling, but slowly. The proven cases show that human genes can carry two protein messages. They do not show that many do.
Part 4: What overlap costs, measured
Part 1 priced a second frame on paper. Four kinds of measurement show whether real overlaps pay that price.
Fewer silent changes in overlapping DNA
If a second frame takes away a gene's silent changes, overlapping DNA should collect fewer of them over time than the rest of the gene. Etienne Simon-Loriere, Edward Holmes and Israel Pagán tested this on 117 overlapping genes from 55 RNA virus species, comparing the overlapping part of each gene with its non-overlapping part.
Three of their results bear directly on constraint, and they agree with each other.
- Less change overall. The total length of the evolutionary tree, a measure of how much a stretch has changed among related virus strains, "was significantly smaller in OV than in NOV regions" across all genes (P = 0.037), and significantly smaller in 65 of the 117 genes taken one by one.9 The effect was clearest where one gene lies wholly inside another (25 of 31 genes), not where genes merely overlap at their ends.9
- Fewer silent substitutions. The rate of silent substitutions was higher outside the overlaps, for all genes together (P = 0.040) and in 84 of 117 genes.9
- No significant difference in amino-acid-changing substitutions. That rate did not differ significantly between the two parts (F = 0.10, P = 0.747).9 A non-significant result is not proof of no difference, and the authors report that in terminal overlaps the rate tended to be lower outside the overlap, where conserved protein domains cluster.
That pattern is what the code arithmetic predicts: a change that would be silent in one frame changes the other protein, so fewer of them survive. The authors' abstract draws the same conclusion: "gene overlapping generally results in a decrease in the rate of RNA virus evolution through a reduction in the frequency of synonymous mutations."9
The paper's ratio of the two rates, often used as a single measure of selection, does not give a usable answer, because the paper reports it two opposite ways.10
- Its Table 1 and one sentence of the text give mean ratios of 0.23 in the overlapping parts and 0.60 in the rest.9 Counted row by row for this article, Table 1 has the overlapping part lower in 95 of 117 genes.
- Two sentences later the text says the ratio "was significantly higher in overlapping than in nonoverlapping fragments" in 81 of 117 genes.9 That cannot describe Table 1, where only 22 genes have the overlapping part higher at all.
- The supplementary table gives each gene's two rates separately. Dividing them, the overlapping part has the higher ratio in 92 of 117 genes, and the difference of the two rates is higher in the overlap in 98. Table 1's direction matches the supplement's in only 35 of 117 genes.
The individual rates point the supplement's way. With amino-acid changes about the same and silent changes fewer, the ratio in the overlap should rise, not fall. A mean of ratios across genes need not follow the means of the two rates, so this does not prove Table 1 wrong, but nothing in the paper explains its direction, and the authors' own explanation, fewer silent changes, would push the ratio the other way. This article therefore does not use the paper's ratios, in either direction. The authors read their result as "stronger purifying selection in OV than in NOV fragments" linked to constraint on silent change.9 The tree lengths and silent rates above support constraint on silent change. The ratio comparison does not settle anything further.
There is also a reason to distrust the ratio in overlaps at all. The usual methods for spotting selection assume each gene is free to change on its own. Niv Sabath, Giddy Landan and Dan Graur showed that applying them to overlapping genes "leads to the false appearance of positive selection even though the gene is in reality subject to negative selection"; for one influenza gene, PB1-F2, the ratio fell from 6.75 to 0.52 when both frames were modelled together.11 A higher ratio in overlapping parts, as the supplement shows, is what that bias would produce. Silent-substitution counts and tree lengths are more direct evidence of constraint than the ratio.
The same effect is now used to find overlaps. Andrew Firth searched alignments of about 700 RNA virus species for windows where silent changes are scarcer than expected, and his method picked out 20 of 21 known overlapping genes and 16 of 17 known overlapping non-coding elements.40 In hepatitis B virus, which reads half its letters twice, Mizokami and colleagues found the fewest silent changes in the surface gene, which lies wholly inside the polymerase gene, and explained it by the overlap: "synonymous substitutions in the S gene can cause amino acid changes in its overlapping region in a different frame."41
The two frames need not share the load evenly or in a fixed way. In the human GNAS locus, Nekrutenko and colleagues found the two frames taking turns, with amino-acid change shifting from one frame to the other across primate branches; "none of the ratios is significantly different from 1."37
One mutation, two proteins
Hepatitis B virus shows the cost in a form doctors see. The polymerase gene lies wholly over the surface-protein gene. Shuping Tong and Peter Revill: "Due to the overlap between the P and envelope genes, drug resistant mutations in the P gene can cause amino acid substitutions in or near the 'a' determinant of envelope proteins," and one resistance mutation "can convert codon 172 in the S region into a stop codon, thus truncating the envelope proteins by 55 residues."12 Some of these surface changes weaken antibody binding and help the virus escape the immune system.
This is the kind of mutation the design writers have in view: helpful in one frame and changing the other. The same review shows both outcomes. The surface changes that weaken antibody binding leave the virus better off on two counts at once, resisting the drug and dodging antibodies. The stop codon made by the rtA181T change leaves it worse off on one: "This mutant was deficient in both virion and HBsAg secretion."12 So "changes the other protein" is not the same as "harms the virus." And the overlap does not stop resistance from arising. For the older drug lamivudine, "LMV resistant mutations develop in 23% of patients following 12 months of therapy, rising to 80% by five years of treatment," though resistance to the newer tenofovir is uncommon.12
Pulling a genome's overlaps apart
The most direct test available is to rebuild a genome without its overlaps and see what happens. In 2012 Paul Jaschke, Drew Endy and colleagues designed phiX174.1, a 6,302-letter version of the 5,386-letter phage in which every coding sequence and its start signals were separated. To keep the genome at its natural length they then cut the gene for coat protein F short, making phiX174.1f. They recovered living phages "containing the øX174.1f genome from E. coli C strains that independently express full-length gene F," so the rebuilt genome works only with gene F supplied from the host.13 Their own summary: "A first fully decompressed phage genome is viable and not irreducibly complex."13
In 2020 Bradley Wright, Jaschke and colleagues measured how well the rebuilt phage grows. On the same host strain, which carries gene F on a plasmid, the natural phage released "145 ± 66 phage produced per cell while decompressed ϕX174 produced only 16 ± 2," and made plaques of 29.2 against 4.1 square millimetres.14 The rebuilt phage was also less stable when heated and attached poorly to cells at 14 degrees.14 Lysis timing was the same.
The comparison has three confounds, and they stop the gap from being assigned to overlap removal alone.
- The plates were not identical. "Wild-type ϕX174 plates contained 2 mM glucose and decompressed ϕX174 plates contained 2 mM rhamnose," the sugar that switches on the plasmid's gene F.14 The rebuilt phage depends on a gene supplied from a plasmid in the host, a complementation set-up the natural phage does not need, and part of its deficit may come from that or from the different sugar.
- Separating the genes meant changing many letters to remove the second frames while keeping the first proteins, so the rebuilt genome differs from the natural one in more than its layout. Wright and colleagues point out that the packaged DNA itself differs, which may explain the heat instability.14
- These are preprint numbers. The paper was published in ACS Synthetic Biology later in 2020, and that version was not compared here.
With those limits, the redesigned phage released about one-ninth as many phages per infected cell as the natural one. The experiment shows that phiX174 can replicate in the laboratory with its overlaps separated, provided the host supplies gene F; it does not show that the redesigned genome is viable on its own. It shows that the natural, overlapping arrangement does much better than this one redesign, but not how much of the difference comes from the overlaps themselves rather than from the many letter changes, the gene F supply or the plates.
Deleting a young overlapping gene
The bacterial case runs the other way: remove the new frame and keep the old. Fellner and colleagues stopped the young overlapping gene nog1 with a single-letter change that leaves the protein of the gene it sits in, citC, unchanged. Grown alone, the mutant grew like the parent except for slower growth with added magnesium chloride. In head-to-head competition the mutant lost: "Δnog1 shows a clear loss of fitness compared to its parental strain when grown in plain LB. When using MgCl2-supplemented broth, the decrease in fitness of Δnog1 was even more pronounced."26 In their Figure 4b the mutant falls to about a fifth of the mixed culture in plain broth and to about 4% with magnesium chloride, while a control change in citC has no effect. The authors call the benefit small: "the gene product provides only a small fitness advantage for its carrier."26 A 2020 review from the same research group reports the result backwards, as a "growth advantage over the wildtype" for the mutant.28 The original figure settles it.
What these measurements do and do not show
They show that overlapping letters collect fewer silent changes than single-frame letters, that one mutation in an overlap really does change two proteins, and that the natural overlapping phiX174 far outgrows one redesign without overlaps. That is the constraint part of the design-side argument, and it stands.
They do not measure how often a mutation in an overlap is beneficial, and so they do not support the further claim that helpful changes in overlaps are vanishingly rare. None of these studies counts beneficial mutations, and none supplies the "fraction already optimal" that the Montañez model in Part 7 assumes. The share of mutations that help, harm or do nothing is measured directly in Beneficial and Harmful Mutations, and the overlap data here do not change those measurements.
Part 5: How second frames arise: the overprinting evidence
What overprinting predicts
In 1992 Paul Keese and Adrian Gibbs gave the idea its modern form. Many virus genes, they noted, are found only in one family or subfamily, including "a range of overlapping genes that occur in viruses from about 50% of all families."3 Their proposal: "proteins may arise continuously de novo," and one way is "overprinting," in which an existing stretch of DNA is translated in a different frame.3 If that is how second frames arise, several things should be true, and each can be checked. None of them is unique to overprinting, so each is a test the idea could fail, not proof that it happened.
- The newer frame should be found in fewer lineages than the older one. The older gene is shared widely; the overprint only in the group where it arose.
- Its protein should look unlike ordinary proteins. It begins as whatever the old gene's letters happen to spell in another frame. A new product may lack a compact fold, and disorder can make a second frame easier to establish or retain. The genetic code itself also biases what amino acids that frame can spell.
- Its codon usage should not match the rest of the genome, because its codons were chosen for the older frame.
- Long open frames to start from should already be common inside existing genes.
- New overprints should matter less than old ones, since they have had less time to become built into the cell.
Lineage limits
The newer frame is usually the more narrowly shared. Corinne Rancurel, David Karlin and colleagues traced the distribution of overlapping virus proteins and "were able to confirm that 17 of these had been created de novo"; "Most proteins created de novo are orphans (i.e., restricted to one species or genus)."15 Niv Sabath, Andreas Wagner and David Karlin found 12 overprinted genes whose distribution on the family tree was clear enough to place their origin.42 Specific cases repeat the pattern.
- HIV-1 ASP. The antisense frame is open in 77% of pandemic group M sequences, with a median length of 182 codons, and not in the other HIV-1 groups or the related ape viruses; Elodie Cassan, Olivier Gascuel and colleagues conclude that it "was created recently and that its emergence in HIV-1 is concomitant with the emergence of the group M."17
- SARS-CoV-2 ORF3d. Chase Nelson and colleagues found it in SARS-CoV-2 and "also present in Guangxi pangolin-CoVs but not other closely related pangolin-CoVs or bat-CoVs."21
- Bacterial nog1. "A phylogenetic analysis suggests that the orphan gene nog1 arose by overprinting after Escherichia/Shigella separated from the other γ-proteobacteria."26
- Human dual-coding regions. Liang and Landweber checked the second frames of human dual-coding regions in other species: "one of two reading frames contains one or more stop codons in 80% of the dual coding regions in chickens, 45% in mice, and 54% in rats," and concluded that "most secondary reading frames in the dual coding regions examined evolved recently in mammals."29
- ARF. Chickens have no p16INK4a at all, and their ARF stops at the end of its own exon. Yet chicken exon 2 keeps an open frame "capable of encoding 162 amino acids that could have formed the carboxyl terminus of ARF, if splicing had occurred in the register used in mammals."35 The second frame is there in the chicken letters and unused.
- POLGARF. Gary Loughran, John Atkins and colleagues traced its start to a repeat element that inserted into the gene's leader: "POLGARF evolved ∼160 million y ago due to a mammalian-wide interspersed repeat (MIR) transposition into the 5′ leader sequence of the mammalian POLG gene, which became fixed in placental mammals."19
- The ribosome profiling catalogue. 89.6% of the 7,264 human frames have no detectable protein counterpart in mammals outside the primates.18 This is a comparative finding with a known blind spot: short proteins give homology searches little to work with, so some frames may have relatives that were not detected. It marks the frames as lineage-limited by that test, not as dated by an independent clock.
A frame found in fewer lineages is not by itself proof that it was born late. The same pattern could arise if relatives lost the frame, if it diverged beyond recognition, or if it was carried in from elsewhere. Evidence supporting a direction in the cases above includes an outgroup that keeps an open but unused frame, as the chicken does for ARF; an insertion that supplied the start, as for POLGARF; or codon usage that marks which frame came first, as in the next sections. These clues strengthen an inferred direction; they are not observations of a gene's birth.
Unusual proteins
The newer frames make odd proteins. Rancurel and colleagues found that "almost two-thirds of novel, overlapping protein regions are disordered," against "fewer than one-fourth of nonoverlapping protein regions," meaning they lack a fixed folded shape.15 Erika Kovacs and colleagues found the same in human genes that code in two frames by alternative splicing. Of 67 such genes, the older and newer frames could be told apart in 62, and the newer frame was mostly disordered "in 39 out of 62 cases." Where the second frame is shifted one way, disorder rose from 16.7% in the older frame to 56.3% in the newer; where it is shifted the other way, both frames were already high, 47.2% and 58.2%.43 Rancurel's group was careful about direction: overlapping genes may be disordered because new proteins start that way, or because disorder makes an overlap easier to keep, and "our results are compatible with both hypotheses."15 Sara Willis and Joanna Masel separated the two and estimated that "birth facilitation accounts for 32 ± 13 %" of the extra disorder in overlaps, the rest coming from the genetic code and from which overlaps survive.44
Codon usage
Angelo Pavesi, Gkikas Magiorkinis and David Karlin took 25 overlapping virus genes whose older frame was known from family trees. The older frame's codon usage was closer to the rest of the genome in 21 of the 25, with a median correlation of 0.42 against 0.19 for the newer frame.16 The newer frame reads letters that were chosen for something else, which is what overprinting predicts. Codon usage and disorder can also be shaped by selection after a frame arises, so these features are consistent with overprinting, not exclusive proof of it.
Frame bias from the genetic code
Willis and Masel also found that genes born inside other genes favour the +1 frame: "23 +1 events vs. 1 +2 event, P = 3 × 10 − 6."44 They trace this to the code itself: "+1 frameshifts yield more start codons and fewer stop codons, and hence more and longer ORFs."44 The frame most often used is the frame the code leaves most open.
Long open frames are common inside ordinary genes
If overprinting is real, long open frames should already lie hidden in ordinary genes. The test run for this article took 19,167 human MANE Select coding sequences and looked for ATG-started open frames of 100 codons or more inside each gene, in the five other frames. It compared the real genes with four kinds of shuffled copies, three replicates each: bases shuffled; codons shuffled; synonymous codons shuffled within the gene, which keeps the protein exactly; and synonymous codons shuffled only within windows of ten codons, which also keeps the local base composition.45
| Frame | Real genes | Bases shuffled | Codons shuffled | Synonymous shuffle, whole gene | Synonymous shuffle, 10-codon windows |
|---|---|---|---|---|---|
| +1 | 14.5% | 8.2% to 8.5% | 8.4% to 8.8% | 9.7% to 10.1% | 13.2% to 13.5% |
| +2 | 3.0% | 8.0% to 8.4% | 2.6% to 2.7% | 3.1% to 3.2% | 3.1% to 3.3% |
| Opposite strand, aligned | 21.5% | 9.4% to 9.5% | 20.2% to 20.6% | 21.4% to 21.8% | 21.4% to 21.5% |
| Opposite strand, shifted 1 | 25.1% | 9.3% to 9.5% | 21.3% to 22.0% | 21.3% to 21.9% | 23.9% to 24.1% |
| Opposite strand, shifted 2 | 5.1% | 9.3% to 9.5% | 3.3% to 4.0% | 4.5% to 4.6% | 4.5% to 4.9% |
About one human gene in seven hides an ATG-started frame of 100 codons or more in its +1 frame, and one in four or five in the two richer opposite-strand frames. Shuffles preserving the protein sequence and local base composition reproduce most of that frequency. Shuffling synonymous codons within short windows, which keeps both, leaves 13.2% to 13.5% in the +1 frame against 14.5% in real genes, and matches or comes close to the real share in the other frames. For frames of 150 codons or more, about the length Chung's group used, the +1 share is 4.6% in real genes and 3.9% to 4.1% after the window shuffle. By this test the +1 frame shows no avoidance of stop codons beyond what the protein and the local composition produce: it carries 4.21 stop codons per 100 codons in real genes and 4.25 to 4.38 in the shuffles that keep the protein or the codons.45 The test covers open frames and stop codons under these four nulls. It does not test other sequence features, and it does not show that any hidden frame is expressed or could become a working gene.
Two readings follow. First, the raw material for overprinting is common: thousands of human genes carry long open second frames, at rates the protein-preserving shuffles nearly reproduce without any second protein. Second, an open frame proves little. This is why Chung's group asked a harder question. Their "nearly impossible by chance" referred to frames of 500 letters or more that stay open in both human and mouse alignments, and their simulations gave "Only 0.3% of simulated alignments preserved ARFs with 500 or more nucleotides."1 An open frame in one genome is common; a long open frame kept across two species is not. The two numbers measure different things, and the design-side quotations in Part 7 use the second as though it described the first.
Old overprints matter more than new ones
Sabath and colleagues collected what was known about suppressing each of their 12 overprinted genes: "Mutations in the youngest de novo genes (taxa 1–3) have little to moderate effect, whereas suppression of older de novo genes (taxa 6 and 10–12) has severe effects."42 nog1 fits the young end: a measurable but small cost when removed. This fits a picture in which new frames start with weak or no function and some become built into the virus or cell over time.
What overprinting does not settle
- No birth has been watched. Every origin in this Part is inferred by comparing living sequences. No source reviewed here observed a second frame arising and becoming useful in real time.
- The ages are mostly relative. Sabath and colleagues write that "we are unable to estimate the absolute age of de novo genes but only their age relative to the divergence of a viral housekeeping gene."42 Cassan's group gives no date in years for ASP. Loughran's "∼160 million y" for POLGARF is a date in years, but it is read off the mainstream calibration of the placental-mammal tree, not measured by an independent clock.19 "Young" in this section means "found in fewer branches of the tree," and nothing here bears on a number of years.
- The ages rest on comparing species. Saying that a frame "arose after" a split assumes the species share an ancestor. For strains of one virus or one bacterial genus, few readers on either side would deny a shared ancestor, so the inference is close to common ground there. For human and mouse, or human and chicken, it is the very point in dispute, and the lineage pattern is evidence for overprinting only to a reader who accepts that framework or finds it the best fit to the pattern.
- The first step is not free, and it is not explained by the open frame alone. A hidden frame becomes a second gene only if it is transcribed, given a start the ribosome uses, and its product does something. Those steps may need mutations, and making an unused product costs the cell something. What overprinting removes is the need to invent two unrelated full sequences at once; it does not show which mutations occurred or how the new frame became established by selection or drift.
- Overprinting explains the second frame, not the first. Every case here starts with a working gene. Where the first gene came from is the subject of Orphan Genes, How Rare Is a Working Protein? and Gene Duplication and Divergence.
- A disordered, young, lineage-limited protein can still be important. The pattern says how the frames look, not that they are useless. ARF is a major tumour suppressor in mammals and is absent as an overlap in chickens.
Part 6: Other codes laid over coding DNA
Protein frames are only one kind of layer. Coding letters also carry signals that are not read by ribosomes. These are real, and some are well measured. They are also different in kind from a second protein, and the "number of codes" claims in Part 7 depend on counting them as though they were the same.
Trifonov's codes
Edward Trifonov argued from the 1980s that DNA carries several messages at once. His 1989 paper, by its PubMed abstract, named three: the protein code with its reading frame, a chromatin code governing how DNA wraps around nucleosomes, and a putative code for DNA loops, and said that the codes overlap "so that some nucleotides belong to several messages simultaneously."46 The body of that paper could not be obtained for this article, by three routes, so whether it lists twelve codes is unchecked. His 2011 review is more sweeping in tone and gives no count: "one has to admit that the genetic sequences carry many different codes. If we are to know what the sequences are about, we have to detect and decipher these codes. The times of surrender to 'junk' and 'selfish DNA' are over."47 Trifonov is a mainstream molecular biologist, and his argument is a scientific claim about signals in sequence, not a design argument.
Room in the code for extra signals
Shalev Itzkovitz and Uri Alon asked in 2007 whether the standard genetic code is unusually good at letting coding letters carry other messages. Compared with alternative codes, it is: the title of their paper says it is "nearly optimal for allowing additional information within protein-coding sequences."48 They offered two evolutionary explanations for how that came about: "One possibility is that the ability to include parallel codes within protein-coding sequences conferred a selection advantage during the early evolution of the genetic code. Alternatively, the genetic code might have been fixed in evolution before most parallel codes existed," and they went on to link the code's structure to its resistance to frame-shift errors.48
In 2010 Itzkovitz, Eran Hodis and Eran Segal looked for extra signals in the coding regions of about 700 species by asking where codon choice departs from what the protein and codon preferences alone predict. They found RNA folding signals, avoided start and initiation signals in bacteria, and microRNA target sites in eukaryotes, and wrote that "genomes may have evolved to encode extensive overlapping information within protein-coding regions."49 Their measure is a departure from randomised codon choice, so it finds the layers that bias synonymous codons. It does not count protein frames, and it does not say how many letters are fixed by each layer.
Protein-binding sites inside exons
In 2013 Andrew Stergachis, John Stamatoyannopoulos and colleagues mapped where regulatory proteins sit on DNA inside exons across 81 cell types and reported that "~15% of human codons are dual-use codons."50 That is the largest figure for a second layer inside human coding sequence among the sources reviewed here. It counts codons covered by a protein footprint in at least one cell type, which is not the same as codons whose letters are fixed by that protein. A 2015 reanalysis disputed the claim that these sites constrain codon choice, and the article Synonymous Mutations and Codon Usage weighs the two papers. This article does not repeat that argument.
Signals in the untranslated ends
Short frames in the leader of a messenger RNA, before the main start codon, regulate how much protein is made. Sarah Calvo, Vamsi Mootha and colleagues found that 49% of human transcripts carry at least one, and in reporter tests "uORFs typically reduce protein expression by 30–80%."51 These are the most common translated "second frames" in the ribosome profiling catalogue, 3,083 of the 7,264.8 Most of them lie outside the main coding sequence, so they add a regulatory layer to the message rather than a second protein on the same letters.
What these layers do and do not add up to
- They are real and measured. RNA structure, splicing signals, binding sites and upstream frames all constrain some letters inside genes.
- They can constrain synonymous choice. They bias which of several codons for the same amino acid is used. A binding site or a folding signal usually tolerates many sequences. So a letter in a binding site may not fit the Montañez model of Part 7, in which any change at an optimal letter of a code is harmful.
- No source reviewed here counts how many layers sit on a typical letter. The figures are per layer: 15% of codons in some footprint, 49% of transcripts with an upstream frame, 0.064% of coding letters in two accepted protein frames. None of the sources read here gives a distribution of the number of layers per letter, and none supports "up to 12" as a typical count.
Part 7: The design-side claims, tested
What the design writers say
Seven writers make the overlap argument in print. Their own words come first.
- John Sanford, Genetic Entropy, 2005 and 2014. "There is abundant evidence that most DNA sequences are poly-functional, and are therefore poly-constrained. This fact has been extensively demonstrated by Trifonov (1989). For example, most human coding sequences encode for two different RNAs that read in opposite directions (i.e., both DNA strands are transcribed; Yelin et al., 2003)." Then: "Trifonov (1989) has shown that probably all DNA sequences in the genome encrypt multiple codes (up to 12). In computer science, this type of data compression can only result from the highest level of information design." And: "Logically, all poly-constrained DNA had to be designed."2 The 2005 edition makes the same three claims.52 The 2014 edition adds: "We now know that over half the human genome is read backwards as well as forwards. A given gene produces roughly 7 overlapping transcripts," and concludes that "Mutations that are both unambiguously beneficial and also impactful enough to be selectable, must be vanishingly rare."2
- George Montañez, Robert Marks, Jorge Fernandez and John Sanford, 2013. A probability model, examined in its own section below. Sanford's book summarises it: "Any change at all will diminish total information with nearly absolute certainty. My colleagues and I have demonstrated this mathematically in a recent paper (Montañez et al., 2013)."2
- Jonathan Wells, The Myth of Junk DNA, 2011. Wells reports a 2007 study's "overlapping coding regions," a feature "that the researchers concluded 'is nearly impossible by chance.'" and Itzkovitz and Alon's finding that the code is "nearly optimal for allowing additional information within protein-coding sequences."53 He quotes William Dembski: "If these systems are designed ... we can expect the information to be densely packed and multilayered," and concludes that "dense, multi-layered embedding of information is a prediction of intelligent design."53
- Fazale Rana, Reasons to Believe, 2007. "The apparent widespread occurrence of overlapping genes doesn't make much sense from an evolutionary perspective, because of the cost they represent to the organism." "For me, the only explanation for overlapping genes is the work of a Creator. It's hard to imagine how undirected evolutionary processes could produce overlapping genes."54 Reasons to Believe is an old-earth ministry, so this is a design argument, not a young-earth one.
- Jeffrey Tomkins, Institute for Creation Research, 2014. On the 2013 binding-site paper: "a new, previously hidden genetic code embedded within the same sections of genes that code for proteins—utterly defying all naturalistic explanations for its existence." And: "Some genes overlap parts of other genes in the genome, and now it has been revealed that many genes have areas that contain dual codes within the very same sequence."55
- Ann Gauger, Ola Hössjer and Colin Reeves, 2017. "In addition, this is a strong argument for design: the sophistication required to have codes within codes would not arise by simple mutation and natural selection."56
- Casey Luskin, Science and Culture Today, 2026. "Another discovery in genomics predicted by ID is the presence of overlapping genes, where the same stretch of DNA encodes multiple genes. One strand of DNA might encode one protein while the opposite strand encodes an entirely different functional protein (see, e.g., Adelman et al., 1987 or Mehravar et al., 2021). Such an astonishing degree of information compression is only expected if a mind of genius encoded the information in DNA, not random mutation and blind natural selection."57
Claim by claim
| Claim | What it would need | What the evidence shows |
|---|---|---|
| Overlaps constrain change | Measured loss of mutational freedom in overlaps | Supported. Silent changes are scarcer in overlaps; one mutation changes two proteins; the natural phiX174 released about nine times as many phages per cell as one redesign without overlaps, a gap with several causes (Part 4) |
| "Most human coding sequences encode for two different RNAs" (Sanford, citing Yelin 2003) | A share of human coding genes with antisense RNA | Overstated as cited. Yelin's abstract reports that "> or =60% of this data set, or approximately 1,600 predicted sense-antisense transcriptional units, are transcribed from both DNA strands," a data set of 2,667 pairs, not all genes.58 By GENCODE 50, 55.3% of protein-coding genes have another gene's exon on the opposite strand, which is close to Sanford's "most" for RNA. Two RNAs are not two proteins |
| "Over half the human genome is read backwards as well as forwards" | A genome-wide count of both-strand transcription | No source given. The paragraph cites nothing for it. The nearest source in the chain, Kapranov 2007, says "most of both strands of the human genome might be transcribed."59 Transcription is a different claim from coding. The Noncoding RNA article weighs it |
| "Probably all DNA sequences ... encrypt multiple codes (up to 12)" (citing Trifonov 1989) | A per-letter count of codes across the genome | Unverified, not refuted. Trifonov's abstract names three codes and says "some nucleotides belong to several messages."46 The body could not be read. No source read for this article counts codes per letter (Part 6) |
| Two protein frames are widespread | A share of coding letters read in two frames | Not supported for large genomes, by the annotation. 0.064% of human coding letters on one canonical transcript per gene, 7.2% counting every annotated transcript, against 50% of the hepatitis B genome (Parts 2 and 3). Annotation can miss real frames and include unused ones |
| "Nearly impossible by chance" (Wells, from Chung 2007) | That a second frame could not arise by mutation | Misapplied. Chung's sentence is "although dual coding is nearly impossible by chance, a number of human transcripts contain overlapping coding regions," about frames of 500 letters kept open across human and mouse, and the same abstract says dual coding "may be advantageous."1 Open second frames are common in single genomes (Part 5) |
| "Utterly defying all naturalistic explanations" (Tomkins, on Stergachis 2013) | That the source says so, or an argument that it does | Not in the source. Stergachis's paper treats the binding sites as shaped by evolution.50 Tomkins's "dual codes" sentence cites Sanna 2008, which counts overlapping gene pairs and says nothing about dual codes.30 |
| Overlapping genes were "predicted by ID" (Luskin, Wells) | A dated prediction, made before the discovery, of where overlaps would be found | Retrospective. Overlapping genes in phiX174 were reported in 1976.60 The quotations offered are from 2004 and later, and none predicts the pattern found: dense in small viruses, sparse in large genomes, young and lineage-limited where dated |
| Adelman 1987 shows "an entirely different functional protein" on the opposite strand (Luskin) | A second protein with a function | Not in the source. Adelman's second gene "is transcribed from the other DNA strand to produce RNA of undefined function."61 Luskin's other citation, Mehravar 2021, says "the advantage and biological significance of overlapping genes is still unclear."62 |
Two of these points need saying plainly, because they cut in the design writers' favour. First, their core intuition is right: an overlap does tie two messages together, and the code arithmetic, the virus comparisons and the decompressed phage all confirm it. Second, Rana's observation that overlaps are costly is the same observation Chung's group made, and it is a fair question for any evolutionary account. The overprinting literature's answer is that new frames begin as passengers on old genes, usually short and disordered, and that mutual constraint grows as both products become needed; Part 5 lists what that answer still has to show. The measurements support the coupling claim. They do not establish either a genome-wide count of overlapping codes or the claim that a mind is the only possible explanation. Protein-frame annotations do not measure every proposed regulatory layer, and the comparative history of some second frames supports addition within an older gene.
No design writer read for this article cites the overprinting literature. Searches of Sanford's two editions, the Montañez chapter, Wells's book, Rana, Tomkins, Gauger and Luskin found no mention of Keese and Gibbs, Rancurel, Sabath or Pavesi.
The Montañez model
The model is the main quantitative design argument examined here, so its assumptions and calculations need to be tested.
What it says, in its own terms. The model takes a stretch of "functional sequence" in which each letter belongs to between zero and twelve codes. Membership is defined by consequence: "A nucleotide that does not contribute to a given code is considered neutral relative to that code." The authors also "allow for rare nucleotide sites within the functional sequence that are perfectly neutral."20 In each code a letter is already the best letter with probability p(optimal); the authors use 99.9% "for the analysis and discussion" and call it "a tunable parameter."20 At a code site that is not yet optimal, "only one of the three alternative nucleotides will be an improvement" over the existing letter.20 At an optimal site, "any change is deleterious, assuming no neutral changes." And: "There are no absolutely neutral positions in any given code, because by definition such a position is not part of that code."20 So for one code, with r = 0.001 the share not yet optimal, a random change helps with probability "p(B) = (1/3)r = 0.00033."20 The preferred letter of each code "is chosen independently of the other nucleotides that are optimal for codes that may overlap at that position."20 A change that helps one code and harms another counts as neutral, and a change is beneficial overall when it helps a majority of the codes the letter belongs to. The probability of that is their equation (1), a binomial sum. For r = 0.001 and L = 12 they print "p(B)L = 4.15 × 10−22 or I[L+] = 71 bits."20
The arithmetic follows from the inputs. Recomputed for this article, the same inputs give 3.6 x 10-22, or 71.2 bits. The printed figure is what a per-code chance of 0.00034 gives, which looks like 1/3 of 0.001 rounded up in the last digit; the text itself prints 0.00033.63 The difference is about 15% in a number of order 10-22 and is immaterial to the conclusion. The two-code values printed on page 146 reproduce exactly.
| Share of letters not yet optimal, r | 1 code | 2 codes | 3 codes | 5 codes | 12 codes |
|---|---|---|---|---|---|
| 0.001 (the paper's case) | 3.3 x 10-4 | 1.1 x 10-7 | 3.3 x 10-7 | 3.7 x 10-10 | 3.6 x 10-22 |
| 0.01 | 3.3 x 10-3 | 1.1 x 10-5 | 3.3 x 10-5 | 3.7 x 10-7 | 3.6 x 10-15 |
| 0.1 | 3.3 x 10-2 | 1.1 x 10-3 | 3.3 x 10-3 | 3.5 x 10-4 | 3.1 x 10-8 |
| 0.5 | 0.17 | 2.8 x 10-2 | 7.4 x 10-2 | 3.5 x 10-2 | 1.3 x 10-3 |
The two-code column is lower than the three-code column because, with two codes, a change must help both to count as beneficial, while with three it need help only two. The authors note the same odd-even effect.20
So the result is only as good as the inputs. The algebra is not in question. What can be checked is whether each input matches what is known about real overlapping codes.
- The number of codes, L. The paper's twelve-code sequence is "a hypothetical 100 base pair sequence," and the authors write that "The exact degree to which nucleotides participate in two or more codes is still unknown."20 For protein frames, the annotated share of human coding letters in two frames is 0.064% on one canonical transcript per gene and 7.2% across every annotated transcript; three frames are rarer still. The 17 kinds of code listed on pages 154 and 155 include isochores, nucleosome binding, three-dimensional organisation and "other undiscovered coding structures."20 Several of these are preferences that tolerate many sequences (Part 6), so it is unclear how many of their letters would meet the model's definition of a code site. No source reviewed here gives a count of codes per letter to plug in.
- The share already optimal, 99.9%. This is assumed, and the paper says so.20 No measurement of the share of letters already at their best possible letter was found for any code in the sources reviewed here. The table shows the answer changing by about nineteen orders of magnitude between r = 0.001 and r = 0.5 at twelve codes.
- What counts as a code site. Within the model, every change at an optimal site of a code is harmful by construction, because a site where change does not matter is defined as outside that code. That is a modelling choice, and the algebra built on it is consistent. The open question is how many real letters meet the definition. Real outcomes are mixed letter by letter and change by change. At a protein-coding letter, some substitutions are silent and others are not, depending on which letter replaces which. In a same-strand protein overlap, 48% of single-letter changes leave at least one protein sequence unchanged, by the code arithmetic of Part 1, a protein-sequence count rather than a fitness measurement.4 The authors themselves say that "each code must be partially 'degenerate' (imperfect) to create the 'flexibility' required to allow other overlapping codes."20 How real mutations in overlaps divide among harmful, neutral and helpful is a measurement, and the measurements are weighed in Beneficial and Harmful Mutations. No site-by-site measurement of an overlap scored in both frames was found among the sources reviewed here.
- Codes choose their preferred letters independently. This assumption is about which letter each code prefers. It is not the trivial fact that overlapping codes share letters; the model already lets them share. Whether the preferred letters of real overlapping codes are correlated has not been measured in any source reviewed here, and a correlation could push the answer either way. If two codes tend to prefer the same letter, a change that helps one tends to help the other, and beneficial changes become more likely than the model says. If they tend to prefer different letters, the reverse. The code arithmetic of Part 1 shows that effects on protein sequence are correlated between frames, a change silent in one frame being usually an amino-acid change in the other, but that is a statement about effect classes, not about which letter each protein would prefer.
- Effects are of equal size and cancel. The authors call this "not realistic" and argue that relaxing it would favour their conclusion: "we can logically expect most such ambiguous mutations to have a net deleterious effect. This is because, not only is it more likely for a random change to damage an optimized system than improve it, the nature of that damage will tend to be more pronounced than any potential improvement."20 That is their argument about effect sizes. The model does not measure effect sizes, so the claim stands or falls on the fitness measurements, not on the model.
What the model does and does not show. It prices the chance that a random change improves an already-optimised multi-code letter, under its stated inputs. Its output is a conditional probability: if a letter sits in L codes whose preferred letters are independent, each 99.9% optimised, with code sites defined as places where change matters, then a beneficial change has the stated chance. It is not a measured fraction of beneficial mutations. It does not model a new code arising, which is what overprinting is: a second frame that did not exist before, built on letters chosen for the first. Sanford's book presents it as a demonstration that "any change at all will diminish total information with nearly absolute certainty."2 The paper shows that for its assumed inputs, not for any measured genome. Sanford's further claims about selection thresholds and genetic decline are weighed in Genetic Entropy.
What survives. For two protein frames, the model's qualitative point about coupling holds: most single-letter changes that alter one protein also alter the other, and only 1.2% leave both protein sequences alone. The 1.2% is not the chance of a helpful change. For a change that improves one protein's amino-acid sequence while leaving the other's sequence alone, the candidate category is the changes that alter only the first protein, about 23% of all changes in each direction by the arithmetic of Part 1.4 What share of those would help is not known from the code, and a change in protein sequence is not the same thing as a change in the organism's fitness. Changes that alter both proteins are not all harmful either, as hepatitis B virus shows: some surface-protein changes made by drug-resistance mutations help the virus escape antibodies (Part 4). So the protein arithmetic supports coupling and trade-offs between the two frames. It does not show that beneficial changes in overlaps are vanishingly rare.
Objections and responses
-
"Annotation misses most overlaps. Ribosome profiling shows thousands more." Annotation does miss some, and this article counts the ribosome-profiling frames separately: 1,469 translated frames overlap a known coding sequence in another frame, and counting every annotated transcript gives 7.2% of coding letters in two frames. Even the generous counts leave two-frame coding a small share of human coding DNA and a tiny share of the genome. Translation is also not function: of the 7,264 frames, 15 meet the strict protein standard so far, though small proteins are hard to detect. If future work proves thousands of working second proteins, the amount side of this verdict would change, and the code arithmetic says each would couple its letters to the host gene.
-
"Silent changes in both frames still matter, so the constraint is tighter than 1.2% suggests." Agreed. RNA structure, splicing signals and codon preferences can make even both-silent changes costly, as Synonymous Mutations and Codon Usage shows. That strengthens the point about constraint where overlaps exist. It does not make overlaps more common.
-
"Overprinting is a story. A new frame would harm the old gene while it was still useless, so selection would remove it." The objection is right that the first step is not free. Making a hidden frame into a gene may need mutations, for a start the ribosome uses or for transcription, and those mutations can change the old protein; even translating an existing frame costs the cell something. What overprinting avoids is the need to invent two unrelated full sequences at once, because long open second frames already lie inside about one human gene in seven (Part 5), ribosome profiling finds translation in 1,469 frames that overlap known human genes in another frame (Part 3), and suppressing the youngest overprinted virus genes has "little to moderate effect" (Part 5).42 Whether the steps of a given origin could persist through selection or drift has to be tested case by case. No source reviewed here watched a second frame arise; every origin is inferred from comparing living sequences, and the youngest cases are in viruses.
-
"Lineage limits could mean each kind was created with its own genes, not that the second frame is young." Presence in one group and absence in its relatives does not by itself prove a late birth. The frame might have been lost in the relatives, diverged beyond recognition, or arrived by transfer, and separate creation is a further option for anyone who allows it. The evidence that fixes direction in the best cases is extra: ASP, ORF3d and nog1 sit inside older genes that their close relatives share, and across 25 virus pairs the frame judged older from family trees usually matches the rest of its genome's codon usage better than the newer frame does (Part 5); chicken exon 2 keeps an open but unused ARF frame; POLGARF's start comes from an inserted repeat. For strains of one virus or one bacterial genus, few readers on either side doubt a shared ancestor, and these features favour a frame arising within the group. For human second frames, a reader who rejects common descent can read "found only in primates" as "made for primates." A design account then has to explain why those frames look like overprints: codon usage taken from the first frame, floppy proteins, frames that stay open in relatives without being used, and in POLGARF a start supplied by a mobile repeat. Design does not forbid these features. No design writer reviewed here predicts them.
-
"Small viruses are compact because they were designed to be. Compression is a mark of design." Compactness can be designed; engineers compress data. The test is whether the overlaps look built in or added on. Where comparisons between relatives date virus overlaps, the newer frame is usually limited to one lineage and has the codon usage and disorder expected of a later addition, though those features are not exclusive proof (Part 5). Compression also tracks genome size and capsid type in ways the virus-biology literature still disputes, which is a question about what limits virus genomes rather than about who wrote them.
-
"The decompressed phage shows that overlaps are essential, so phiX174's genome is irreducibly complex." The engineers who built it say the opposite: "A first fully decompressed phage genome is viable and not irreducibly complex."13 That needs a qualifier both ways. The redesign replicates only when the host supplies gene F from a plasmid, so it is not shown to be viable on its own. And it released about one-ninth as many phages per cell as the natural virus, a gap that mixes the effect of separating the genes with the redesign's many letter changes, the supplied gene F and different plates. The experiment shows the overlaps are not strictly required for replication with help from the host, and that the natural arrangement does far better than this redesign.
-
"Chung's team found dual coding 'nearly impossible by chance.'" Their phrase refers to long second frames kept open across human and mouse, against a codon-shuffling null they themselves call a possible weakness. Their conclusion is that such frames are kept by selection and are "hallmarks of fascinating biology," not that they could not arise. Open second frames in one genome are common (Part 5).
-
"The genetic code is nearly optimal for carrying extra codes. That alone points to design." Itzkovitz and Alon's result is real. It bears on the origin of the genetic code, which is argued in Is the Genetic Code Optimized?, not on how much of a genome is overlapping. The authors give two evolutionary routes to it, and the article does not settle between them and design.
-
From the other side: "The Montañez model is a straw man." It is not. Its assumptions are written out, its algebra follows from them, apart from a rounding difference in one printed figure, and its point about coupling holds for protein overlaps. Its weakness is that its decisive inputs, the number of codes per letter, the share already optimal and the independence of each code's preferred letters, are assumed rather than measured, and its definition of a code site makes every change at an optimal site harmful by construction. It gives a conditional probability, not a measured beneficial fraction.
-
From the other side: "Overlapping genes are just a quirk of viruses. Nothing here needs answering." The design writers have a real point: overlaps couple letters, and an account of genome change has to price that. The evolutionary literature does price part of it, in the silent-substitution comparisons, the decompression experiment and the overprinting studies, though the ratio statistics often used for this are unreliable in overlaps (Part 4). The design writers' errors are about amount, history and the leap from coupling to a vanishing beneficial fraction, not about the existence of coupling.
Issues and unanswered questions
- All counts here depend on annotation, in both directions. Our GENCODE counts use release 50, count stop codons as coding and include the mitochondrial genes; the canonical count is one transcript per gene by Ensembl's choice. Later releases will add more ribosome-profiling frames, and some recent GENCODE genes behind the 28 same-strand pairs came from exactly that work. Some annotated models may also prove unused. The RefSeq virus records omit frames found later, as SARS-CoV-2 shows. None of these counts measures the non-protein layers of Part 6.
- Trifonov 1989 was read only in abstract. Springer and ScienceDirect both refused automated access, and no archive copy was found. So Sanford's "up to 12" is unverified, not refuted.
- The phiX174 numbers come from the 2020 preprint, and the comparison is confounded. The published ACS Synthetic Biology version and the full text of Jaschke 2012 were not read here, and the published numbers may differ. The gap between natural and redesigned phage mixes overlap removal with other letter changes, host-supplied gene F and different plate sugars.
- No site-by-site fitness scan of an overlap in both frames was found. A deep mutational scan of an overlapping virus region, measuring the effect of every single change on both proteins and on the virus, would test the Montañez inputs directly: how many changes are harmless in each frame, how many help, how many harm. None was found among the sources reviewed here; such a study may exist.
- Simon-Loriere 2013 reports its ratio result two opposite ways. Its Table 1 and summary means put the ratio of amino-acid to silent substitution rates lower in overlapping parts; its text count of 81 of 117 genes and its supplementary Table S1 put it higher (Part 4). Which is right could not be settled from the paper, and this article uses only the tree-length and silent-rate results, which agree with each other. A corrigendum or the authors' raw output would settle it.
- The ages of overprints are inferred. Sabath's group gives only relative ages, and Loughran's 160 million years for POLGARF comes from the mainstream calibration of the mammal tree rather than an independent clock. Every origin rests on comparing species, and for human frames the comparisons depend on homology searches that miss short or fast-changing proteins.
- Several human cases have thin follow-up. A 2012 letter commenting on the AltPrP paper was not obtained. The Alt-ATXN1 paper does not establish how widely the frame is shared. Chung's list of 15 high-confidence genes, in its Table 1, was not read.
- The human functional screens mostly test other kinds of frames. Prensner's 57 viability hits are mostly upstream, downstream and non-coding-RNA frames, and Chen's 78 start-overlap hits may act partly through the main gene. A screen designed to separate internal second frames from their host genes would sharpen the human count.
- No young-earth model of overlapping genes was found. The young-earth writers read here cite overlaps as evidence of design but do not say how overlaps came to be distributed as they are, which frames are original, or how fast viral overprints arise. A model would need to say whether second frames were created with the first, arose afterwards, or both, and to predict the codon-usage and lineage patterns of Part 5. Viral cases such as ASP and ORF3d, which young-earth writers could accept as post-creation change, are the natural place to start.
Verdict
Overlapping protein codes are real, and where they occur they couple two messages. Under the standard code, a same-strand overlap leaves only 1.2% of single-letter changes silent in both proteins, and most changes that alter one protein alter the other. Real overlaps show the effect: fewer silent substitutions in overlapping parts of virus genes, mutations that change two proteins at once, and a natural phage that far outgrows a redesign without overlaps. On the existence of that coupling, Sanford, Montañez and their colleagues are right.
That coupling is not a measured beneficial fraction. About 23% of changes in a same-strand overlap alter one protein and leave the other's sequence alone, and how many of those, or of the changes that alter both, would help an organism is not measured in any source reviewed here. The Montañez model is correct algebra on assumed inputs: the number of codes per letter, the share already optimal, independent preferred letters, and a definition of a code site under which every change at an optimal site is harmful. Its output is conditional on those inputs.
The design writers' claims about amount are not supported by the counts reviewed here. Two-frame coding fills half of a hepatitis B genome and 0.064% of human coding letters on one canonical transcript per gene, 7.2% across every annotated transcript. Other "codes" laid over coding DNA are real and can constrain synonymous choice, but no source reviewed here counts how many sit on a typical letter. "Up to 12" codes per letter is unmeasured, and "over half the human genome is read backwards as well as forwards" has no cited source.
Where comparisons between species allow a second frame to be dated, it looks added: limited to one lineage, coding a floppy protein, using codons chosen for the older frame, and in POLGARF started by a mobile repeat. That is what overprinting predicts, though these features are inferred from living sequences and are not exclusive proof, and the design writers read for this article do not engage it. The fair conclusion for a skeptical reader is that overlapping genes are a real, constraining feature of compact genomes, uncommon among the annotated genes of large ones, and not evidence that most DNA is a dense multi-layer code. Whether the first genes, the ones overprinting starts from, could arise without design is a different question, taken up in Orphan Genes and How Rare Is a Working Protein?.
Sources
- 1.
Chung, Wen-Yu et al. “A first look at ARFome: dual-coding genes in mammalian genomes.” PLoS Computational Biology 3(5):e91. 2007. Genome-wide scan of human and mouse gene alignments for alternative reading frames of 500 or more letters that stay open in both species, tested against simulated alignments and by several statistics of selection. It found 40 candidate genes, 15 of them high-confidence, and argued that the arrangement is costly but can be favoured when the two proteins work together. No experiments were done.
- 1a.
Section Abstract: “the skepticism surrounding eukaryotic dual coding is unwarranted: rather than being artifacts, overlapping reading frames are often hallmarks of fascinating biology”. ↩
- 1b.
Table S1: Proportion of substitution types in each codon position of the two frames, averaged over all nucleotide contexts. For frame F0, both-silent (SS) is 3.6%, 0% and 0% at positions 1 to 3, and both-changed (NN) 25.4%, 87.2% and 25.6%: averages of 1.2% and 46.1%. ↩
- 1c.
Section Abstract: “we identified 40 candidate genes with evolutionarily conserved overlapping coding regions”. Section Results: “Intersecting the results of the tests yielded 15 dual-coding genes as high-confidence candidates”. ↩ ↩
- 1d.
Section Results: “Only 0.3% of simulated alignments preserved ARFs with 500 or more nucleotides ( Figure S2 ).” Section Abstract: “Here we show that although dual coding is nearly impossible by chance, a number of human transcripts contain overlapping coding regions.” The authors flag their simulation's assumption that codons evolve independently as a possible weakness. ↩ ↩
- 1a.
- 2.
Sanford, John C. Genetic Entropy. FMS Publications. 4th edition. 2014. Book by a Cornell plant geneticist arguing that harmful mutations accumulate in the human genome faster than selection can remove them, so that the genome is degenerating and is young. Its argument about overlapping codes is one strand of its case that beneficial mutations are too rare to build anything. The electronic edition has no page numbers; claims are located by chapter and section.
- 2a.
Chapter 9, section item 7, Poly-constrained DNA: “Trifonov (1989) has shown that probably all DNA sequences in the genome encrypt multiple codes (up to 12). In computer science, this type of data compression can only result from the highest level of information design and results in maximal information density.” ↩ ↩
- 2b.
Chapter 11, section 3a) Overlapping codes: “We now know that over half the human genome is read backwards as well as forwards. A given gene produces roughly 7 overlapping transcripts, having different start and stop positions.” No citation is given for either figure in the paragraph. ↩ ↩
- 2c.
Chapter 9, section item 7, Poly-constrained DNA: “Any change at all will diminish total information with nearly absolute certainty. My colleagues and I have demonstrated this mathematically in a recent paper (Montañez et al., 2013). There is abundant evidence that most DNA sequences are poly-functional, and are therefore poly-constrained. This fact has been extensively demonstrated by Trifonov (1989). For example, most human coding sequences encode for two different RNAs that read in opposite directions (i.e., both DNA strands are transcribed; Yelin et al., 2003).” “Logically, all poly-constrained DNA had to be designed.” Chapter 11, section 3a) Overlapping codes: “Mutations that are both unambiguously beneficial and also impactful enough to be selectable, must be vanishingly rare.” ↩ ↩
- 2a.
- 3.
Keese, Paul K., and Adrian Gibbs. “Origins of genes: "big bang" or continuous creation?” Proceedings of the National Academy of Sciences 89(20):9489–9493. 1992. The paper that named overprinting as a general source of new genes. It compared overlapping virus genes across related viruses and used each frame's distribution among lineages to say which partner was original and which new. Mirror: PMC
- 4.
Computed for this article by
research/frame-pricing.ts(output inresearch/frame-pricing-output.txt). For each way two frames can lie over each other, the script takes every context of letters in which both codons name amino acids, makes every possible single-letter change, and counts the change as silent, amino-acid-changing or stop-making in each frame. Contexts are weighted by the base composition of 19,167 human MANE Select coding sequences from GENCODE release 50 (A 25.9%, C 25.9%, G 26.3%, T 21.8%), with letters treated as independent. At even composition the single-frame silent share is 24.4% and the same-strand both-silent share is 1.2%. In a same-strand overlap, 24.6% + 24.6% - 1.2% = 48.0% of changes are silent in at least one frame. The arithmetic counts effects on protein sequence only. ↩ ↩ ↩ ↩ - 5.
Counted for this article by
research/genbank-overlaps.ts(output inresearch/genbank-overlaps-output.txt) from NCBI RefSeq records NC_003977.2 (hepatitis B virus), NC_001422.1 (phiX174), NC_001802.1 (HIV-1 HXB2), NC_045512.2 (SARS-CoV-2), NC_012920.1 (human mitochondrion) and NC_000913.3 (E. coli K-12 MG1655). Each letter covered by an annotated coding sequence gets a label for strand and codon position; a letter with two different labels is read in two frames. Pseudogenes are skipped, and letters shared by two products read in the same frame, such as phiX174 A and A*, are not counted as a second frame. Counting phiX174 letters covered by two annotated coding sequences in any frame instead gives 1,392 of 5,386 (25.8%), byresearch/phix-any-frame.ts. ↩ ↩ ↩ ↩ - 6.
Schlub, Timothy E., and Edward C. Holmes. “Properties and abundance of overlapping genes in viruses.” Virus Evolution 6(1):veaa009. 2020. Census of overlapping coding regions of different frames in 5,976 NCBI reference virus genomes, with their lengths, frames and relation to genome size and type. The largest survey of its kind; it counts annotated genes and so misses unannotated frames. Mirror: PMC
- 6a.
Section Abstract: “over half the instances of gene overlap are very small, covering <10 nt, and 84 per cent are <50 nt in length. Despite this, 53 per cent of all viruses still contained a gene overlap of 50 nt or larger”. ↩ ↩
- 6b.
Section 3.1 Total frequency of gene overlap across all viruses: “gene overlap being least common in dsRNA viruses (19%, 95% CI 15–24%), followed by negative-sense ssRNA viruses (−ssRNA; 24%, 95% CI 20–29%) and then positive-sense ssRNA viruses (+ssRNA; 43%, 95% CI 40–45%)”. ↩
- 6a.
- 7.
Counted for this article by
research/gencode-overlaps.ts(output inresearch/gencode-overlaps-output.txt) from the GENCODE release 50 comprehensive annotation (Ensembl 116), nuclear and mitochondrial chromosomes together. Stop codons count as coding. Coding letters read in two or more frames: every annotated transcript 2,796,193 of 38,661,975 (7.2%); one Ensembl canonical transcript per gene 21,568 of 33,891,098 (0.064%, or 0.0007% of a 3.1-billion-letter genome); MANE Select 2,626 of 33,438,724 (0.008%). Gene pairs with overlapping coding sequence on canonical transcripts, readthrough genes excluded: 379 same strand and same frame, mostly the shared exons of the UGT1A cluster; 28 same strand in different frames, the longest POLG and POLGARF at 728 letters, with about half involving recently added GENCODE genes and two in the mitochondrion; 16 on opposite strands, the longest TRMT2A and RANBP1 at 246 letters, median 37. Genes with a stretch read in two frames by their own alternative transcripts: 12,569 counting every transcript; 8,226 counting only full-length protein-coding transcripts (2,756 with 150 or more letters); 1,346 counting only transcripts whose coding sequence matches a CCDS entry (261 with 150 or more letters). Protein-coding genes with an exon of another gene overlapping one of their exons on the opposite strand: 11,124 of 20,107 (55.3%); the partner is itself protein-coding in 26.0%, and the partner's exon covers this gene's coding sequence in 32.1%. ↩ ↩ ↩ ↩ ↩ ↩ - 8.
Deutsch, Eric W. et al. “Expanding the human proteome with microproteins and peptideins.” Nature. May 6, 2026. Consortium reanalysis searching 413 public proteomics datasets, 95,520 mass-spectrometry runs, for the products of the 7,264 GENCODE translated frames, with manual review of spectra and a new measure of whether frames stay open across species. It is the journal version of a 2024 preprint and replaces the preprint's counts. The authors build the annotation the article relies on, so they have an interest in the frames being real; their tier system nonetheless keeps the strict count small. Mirror: Eprints.whiterose.ac.uk
- 8a.
Figure 1b: Biotypes of the 7,264 translated frames: uORFs 3,083; upstream overlapping ORFs 688; internal ORFs 720; downstream overlapping ORFs 61; dORFs 504; lncRNA ORFs 1,917; processed-transcript ORFs 291. The three overlapping classes total 1,469. ↩ ↩ ↩ ↩
- 8b.
Section Microproteins in digest MS/MS datasets: “We found 484 peptides passing FDR thresholds that map to 183 out of the 7,264 ncORFs (around 2.5%)”. Section Microproteins as HLA-I-presented peptides: “Overall, we found 3,116 peptides mapping to 1,785 out of 7,264 Ribo-seq ncORFs (24.6%)”. “Almost all of these peptides (2,937 out of 3,116; 94.3%) were found presented by HLA class-I (HLA-I) alone”. Section Evolutionary insights to interpret ncORFs: “we identified 20 candidates for tier 1A status, but further scrutiny reduced this list to 15 ncORFs due to pseudogenic sequences, a GRCh38 assembly error and insufficient Ribo-seq evidence”. ↩ ↩ ↩ ↩
- 8c.
Section Microproteins in digest MS/MS datasets: “we found that only 2 out of 36 known proteins under 50 amino acids (5.6%) satisfy benchmarks for HUPO-HPP verification”. ↩
- 8d.
Section Evolutionary insights to interpret ncORFs: “ncORFs with biotypes that overlap annotated CDSs (upstream overlapping ORFs (uoORFs), internal ORFs (intORFs) and downstream overlapping ORFs (doORFs)) have higher ORBLv scores, presumably due to constraint to preserve the CDS”. “2,211 of the 7,264 ncORFs (30.4%), including 1,335 of the 2,915 upstream ORFs (uORFs, 45.8%), had a placental mammal ORBLq > 0.9, as compared to 10% expected for untranslated ORFs”. “only 143 of the 7,264 ncORFs (2.0%), including 74 out of 2,915 uORFs (2.5%), have a PhyloCSF24 score of >10”. ↩ ↩
- 8a.
- 9.
Simon-Loriere, Etienne et al. “The effect of gene overlapping on the rate of RNA virus evolution.” Molecular Biology and Evolution 30(8):1916–1928. 2013. Comparison of the overlapping and non-overlapping parts of 117 overlapping genes from 55 RNA virus species in 19 families, measuring tree length, silent and amino-acid-changing substitution rates and their ratio, RNA structure and protein domains in each part. Its tree lengths and silent rates show constraint in the overlaps; its ratio results conflict between Table 1, the text and supplementary Table S1, which gives each gene's two rates separately. Mirror: PMC ↩
- 9a.
Section Results, Association between Selection Pressures and Gene Evolution: “Average d N / d S values were 0.23 ± 0.04 and 0.60 ± 0.05 for OV and NOV regions, respectively.” “d N / d S was significantly higher in overlapping than in nonoverlapping fragments when all types of overlap were considered together (81/117”. Table 1: The two quotations are the first and third sentences of one paragraph. Table 1 prints the same means, with standard errors of ±0.05 inside and ±0.06 outside the overlaps, and has the overlapping part lower in 95 of its 117 gene rows, so the second sentence cannot describe Table 1. ↩ ↩
- 9b.
Section Results: “Similar analysis for d N yielded no significant differences among OV and NOV regions when all genes were considered together ( F 1,234 = 0.10; P = 0.747)”. “Finally, d S was significantly higher in NOV than in OV regions considering all genes together ( F 1,234 = 4.26; P = 0.040)”. “most of the genes showed higher d S value in NOV than in OV regions (84/117”. ↩ ↩ ↩
- 9c.
Section Abstract: “We show that gene overlapping generally results in a decrease in the rate of RNA virus evolution through a reduction in the frequency of synonymous mutations.” ↩
- 9d.
Section Results, Effect of Gene Overlapping in Virus Evolution: “A one-way analysis of variance (ANOVA) using type of region (OV and NOV) as a factor indicated that t was significantly smaller in OV than in NOV regions ( F 1,234 = 4.41; P = 0.037). In agreement with these results, in the majority of genes, t was significantly smaller in OV than in NOV regions (65/117, χ 2 = 56.18, P < 1 × 10 −5 )”. t is the total length of the evolutionary tree for each part of the gene. ↩ ↩
- 9e.
Section Results, Association between Selection Pressures and Gene Evolution: “Thus, overlapping genes are generally subjected to stronger purifying selection in OV than in NOV fragments, which seems to be associated with a greater constraint against synonymous changes regardless of the type of overlap.” ↩
- 9f.
Section Results: “the majority of the genes with internal overlapping exhibited significantly lower t values in OV than in NOV regions (25/31, χ 2 = 23.29, P < 1 × 10 −5 ), but this was not so for both types of terminal overlapping”. ↩
- 9a.
- 10.
Counted for this article by
research/simon-loriere-tables.ts(output inresearch/simon-loriere-tables-output.txt) from the 117 gene rows of Simon-Loriere and colleagues' Table 1 and the 117 rows of their supplementary Table S1 9. Table 1's ratio column has the overlapping part lower in 95 rows and higher in 22, with row means of 0.60 outside and 0.23 inside the overlaps, matching the printed means. Table S1 gives point estimates of the amino-acid-changing rate dN and the silent rate dS for each part: dS is lower in the overlap in 97 genes; dN is higher in the overlap in 62 and lower in 53; dN − dS is higher in the overlap in 98; dN/dS formed from these estimates is higher in the overlap in 92. Table 1's direction agrees with the direction from Table S1 in 35 of 117 genes. The methods section says Table 1's ratios come from the SLAC method and that per-site dN − dS values, reported in Table S1, were used to confirm them, so the values need not match exactly; their directions should mostly agree, and they do not. ↩ ↩ - 11.
Sabath, Niv et al. “A method for the simultaneous estimation of selection intensities in overlapping genes.” PLoS ONE 3(12):e3996. 2008. Statistical method for estimating selection on both frames of an overlap at once, tested on simulated genes and applied to two influenza genes earlier reported to be under positive selection. Graur is a leading critic of claims that most of the human genome is functional.
- 12.
Tong, Shuping, and Peter Revill. “Overview of hepatitis B viral replication and genetic variability.” Journal of Hepatology 64(1 Suppl):S4–S16. 2016. Review of hepatitis B virus genome organisation, replication and mutation, including the clinical consequences of the polymerase gene lying wholly over the surface-protein gene. Mirror: PMC
- 12a.
Section 3.4. Drug resistant mutants of nucleos(t)ide analogues: “Due to the overlap between the P and envelope genes, drug resistant mutations in the P gene can cause amino acid substitutions in or near the “a” determinant of envelope proteins, such as F158Y, F161L, E164D, I195M, W196S. Such mutations reduced binding of anti-S antibodies suggesting immune escape”. “The rtA181T mutation can convert codon 172 in the S region into a stop codon, thus truncating the envelope proteins by 55 residues. This mutant was deficient in both virion and HBsAg secretion”. ↩ ↩ ↩
- 12b.
Section 3.4. Drug resistant mutants of nucleos(t)ide analogues: “For example, LMV resistant mutations develop in 23% of patients following 12 months of therapy, rising to 80% by five years of treatment”. “Although mutations conferring resistance to TDF seem uncommon even following 5 years of treatment”. ↩
- 12a.
- 13.
Jaschke, Paul R. et al. “A fully decompressed synthetic bacteriophage øX174 genome assembled and archived in yeast.” Virology 434(2):278–284. 2012. Synthetic-biology paper from Drew Endy's laboratory that designed and built a phiX174 genome with every gene separated from its neighbours, assembled it in yeast, and tested whether it could produce living phage. Its title and abstract are cited here.
- 14.
Wright, Bradley W. et al. “Genome modularization reveals overlapped gene topology is necessary for efficient viral reproduction.” bioRxiv. Jun 12, 2020. Laboratory comparison of natural phiX174 with a rebuilt genome in which every overlap had been pulled apart, measuring plaque size, burst size, lysis timing, heat stability, attachment and the timing of every phage protein. Both phages were grown on a host carrying gene F on a plasmid, under different sugars. Preprint, later published in ACS Synthetic Biology 9:3079 (2020), whose figures were not compared here.
- 14a.
Section Results and Discussion: “We determined the burst-size of wild-type ϕX174 to be 145 ± 66 phage produced per cell while decompressed ϕX174 produced only 16 ± 2 phage per cell.” “the decompressed ϕX174 plaque size was significantly smaller at 4.1 ± 0.9 mm2 compared to 29.2 ± 8.0 mm2 for wild-type ϕX174”. ↩ ↩
- 14c.
Section Results and Discussion: “The decompressed ϕX174 genome encodes proteins with the same amino acid sequence as the wild-type ϕX174 genome, as such, we expect any structural changes to be either from disruption of the intimate genome-capsid interactions due to packaging a genome of different sequence or from inefficient assembly due to altered protein levels during capsid assembly.” “The results of the attachment assay at 14 °C showed that decompressed ϕX174 virions could not remain stably attached to host cells”. ↩ ↩
- 14d.
Figure 2: “Wild-type ϕX174 plates contained 2 mM glucose and decompressed ϕX174 plates contained 2 mM rhamnose.” ↩
- 14a.
- 15.
Rancurel, Corinne et al. “Overlapping genes produce proteins with unusual sequence properties and offer insight into de novo protein creation.” Journal of Virology 83(20):10719–10736. 2009. Study of 36 experimentally characterised overlapping virus genes, assigning which frame of each pair is ancestral by similarity to proteins in other families, and comparing the sequence properties and predicted disorder of old and new frames. Mirror: PMC
- 15a.
Section Discussion: “our results are compatible with both hypotheses. Indeed, almost two-thirds of novel, overlapping protein regions are disordered (Fig. 4 ), compared with fewer than one-fourth of nonoverlapping protein regions”. ↩ ↩ ↩
- 15b.
Section Abstract: “By analyzing the phylogenetic distribution of overlapping proteins, we were able to confirm that 17 of these had been created de novo and to study them individually. Most proteins created de novo are orphans (i.e., restricted to one species or genus).” ↩ ↩
- 15a.
- 16.
Pavesi, Angelo et al. “Viral proteins originated de novo by overprinting can be identified by codon usage: application to the "gene nursery" of deltaretroviruses.” PLoS Computational Biology 9(8):e1003162. 2013. Test on 25 overlapping virus genes whose older frame was known from phylogeny, showing that the older frame's codon usage matches the rest of the genome more closely than the newer frame's, and turning that into a rule for spotting new frames. Mirror: PMC
- 16a.
Section Abstract: “found that their codon usage was significantly closer to that of the rest of the viral genome than the codon usage of de novo frames”. Section Results: “r sA is higher than r sN in 21 cases (i.e. d 21 >0) and lower (i.e. d 21 <0) in only 4 cases”. “the median r sA (0.42) is significantly (P<0.01) higher than the median r sN (0.19)”. ↩ ↩
- 16a.
- 17.
Cassan, Elodie et al. “Concomitant emergence of the antisense protein gene of HIV-1 and of the pandemic.” Proceedings of the National Academy of Sciences 113(41):11537–11542. 2016. Survey of the HIV-1 antisense protein frame across thousands of HIV-1 and related ape virus sequences, finding it open in most of the pandemic group M and not in other groups or in the ape viruses. Mirror: PMC
- 18.
Sandmann, Clara-L. et al. “Evolutionary origins and interactomes of human, young microproteins and small peptides translated from short open reading frames.” Molecular Cell 83(6):994–1011. 2023. Comparative-genomics and ribosome-profiling study of the 7,264 GENCODE translated frames across primates and other mammals, with protein-interaction tests on a subset. It classified about nine in ten as young by the absence of protein homology outside primates, against controls of length-matched frames sampled from untranslated regions.
- 19.
Loughran, Gary et al. “Unusually efficient CUG initiation of an overlapping reading frame in POLG mRNA yields novel protein POLGARF.” Proceedings of the National Academy of Sciences 117(40):24936–24946. 2020. Ribosome profiling, reporter assays, mass spectrometry and comparative genomics of a second protein read from the human mitochondrial polymerase gene, from John Atkins's laboratory, which has long studied non-standard translation. Mirror: PMC
- 19a.
Section Significance: “POLG mRNA possesses a conserved CUG codon which serves as a start of translation for an exceptionally long overlapping open reading frame (260 codons in human) present in all placental mammals”. Section Abstract: “Phylogenetic analysis shows that POLGARF evolved ∼160 million y ago due to a mammalian-wide interspersed repeat (MIR) transposition into the 5′ leader sequence of the mammalian POLG gene, which became fixed in placental mammals.” ↩ ↩ ↩ ↩ ↩
- 19a.
- 20.
Montañez, George et al. “Multiple overlapping genetic codes profoundly reduce the probability of beneficial mutation.” Biological Information: New Perspectives, 139–167. World Scientific. 2013. Probability model, with a crossword-puzzle illustration, of how the chance of a beneficial point mutation falls as the number of codes sharing a letter rises. Its inputs are assumed rather than measured. Published in a conference volume edited by design advocates, with Sanford, the author of Genetic Entropy, as corresponding author.
- 20a.
Page 155: “The exact degree to which nucleotides participate in two or more codes is still unknown, but if it is at all common, it should profoundly reduce the probability of mutations which are unambiguously beneficial.” ↩ ↩
- 20b.
Page 141: “For illustration, in Figure 1 we show a hypothetical 100 base pair sequence, which participates in 12 partially overlapping codes.” Page 142: “We assume a high degree of optimization within each code, although this assumption can be relaxed, and is a tunable parameter within the model. For the analysis and discussion we assume 99.9% of the nucleotide positions defining a code are already an optimal nucleotide.” ↩ ↩
- 20c.
Page 141: “We assume this sequence is not primarily “junk DNA”, but that for the most part it encodes information, yet we allow for rare nucleotide sites within the functional sequence that are perfectly neutral.” “A nucleotide that does not contribute to a given code is considered neutral relative to that code.” Starting assumptions 1 and 2. ↩
- 20d.
Page 143: “For those nucleotides that are part of a given code but are not yet the optimal nucleotide for that code, we assume that only one of the three alternative nucleotides will be an improvement relative to the existing sub-optimal nucleotide.” Starting assumption 5. ↩
- 20e.
Page 145: “Given the stated assumptions, for any single code, a mutation at a position chosen at random that mutates has a probability of being a beneficial (B) mutation equal to p(B) = (1/3)r = 0.00033.” ↩
- 20f.
Page 144: “In the case of optimal nucleotide bases, any change is deleterious, assuming no neutral changes. Therefore, only r = 1 — p(optimal) = 0.1% = 0.001 are subject to beneficial mutation. There are no absolutely neutral positions in any given code, because by definition such a position is not part of that code.” ↩
- 20g.
Page 147: “For r = 0.001, a value of L = 12, p(B)L = 4.15 × 10−22 or I[L+] = 71 bits.” ↩
- 20h.
Page 148, Figure 3 caption: “The spikes on the surface of the plot, visible near the rear plane of the figure, result from the difference between the majority of an even number of codes and the majority of an odd number of codes.” ↩
- 20i.
Pages 154–155: “17) part of other undiscovered coding structures.” The last of 17 kinds of code listed, beginning with "part of an isochore structure" and "part of a nucleosome binding site." ↩
- 20j.
Page 155: “One interesting requirement of overlapping codes is that each code must be partially “degenerate” (imperfect) to create the “flexibility” required to allow other overlapping codes.” ↩
- 20k.
Page 143: “We assume independence in the designation of different optimal codes, such that the nucleotide deemed optimal for a given code (at a position) is chosen independently of the other nucleotides that are optimal for codes that may overlap at that position.” ↩
- 20l.
Page 156: “In our analysis we have for simplicity assumed that if a mutation has a single beneficial effect and a single deleterious effect, it is counted as neutral. However this is not realistic because we can logically expect most such ambiguous mutations to have a net deleterious effect. This is because, not only is it more likely for a random change to damage an optimized system than improve it, the nature of that damage will tend to be more pronounced than any potential improvement.” ↩
- 20a.
- 21.
Nelson, Chase W. et al. “Dynamically evolving novel overlapping gene as a factor in the SARS-CoV-2 pandemic.” eLife 9:e59633. 2020. Analysis of SARS-CoV-2 and related bat and pangolin virus genomes identifying ORF3d, a frame lying inside ORF3a, with ribosome-profiling and patient-antibody evidence that it is translated, and tracking a variant that removes it. Mirror: PMC
- 22.
Chirico, Nicola et al. “Why genes overlap in viruses.” Proceedings of the Royal Society B 277(1701):3809–3817. 2010. Comparative study of gene overlap across 36 RNA and 26 DNA virus families from the RefSeq records of 2008, testing three explanations: capsid size, mutation rate and gene novelty. It found overlap tied to genome length in icosahedral families only, and no link to mutation rate in a sample of 13 families. Mirror: PMC
- 23.
Belshaw, Robert et al. “The evolution of genome compression and genomic novelty in RNA viruses.” Genome Research 17(10):1496–1504. 2007. Survey of 701 RNA virus reference genomes measuring gene overlap against genome length across families. Its overlap-share figure was later called an error by a coauthor's group (Chirico 2010); its counts of which viruses overlap stand. Mirror: PMC
- 23a.
Section Results: “Our analysis shows 819 instances of gene overlap among the 701 reference RNA viral genomes, with 56% of the viruses having some gene overlap.” ↩
- 23b.
Section Results: “The mean overlap across all virus families represents only 1% of their genome length, but the relationship with genome size results in mean overlap being 2% in families with genomes smaller than the median length”. ↩
- 23a.
- 24.
Brandes, Nadav, and Michal Linial. “Gene overlapping and size constraints in the viral world.” Biology Direct 11:26. 2016. Survey of overlap across viral families and genera, testing whether capsid volume drives compression. It found a roughly fixed ceiling on the absolute amount of overlap rather than a fixed share, and no difference between icosahedral and other capsids. Mirror: PMC
- 25.
Johnson, Zackary I., and Sallie W. Chisholm. “Properties of overlapping genes are conserved across microbial genomes.” Genome Research 14(11):2268–2272. 2004. Census of annotated overlapping gene pairs in all microbial genomes sequenced by 2003, with their frames, directions and lengths. The number of overlaps rose with the number of genes at a slope of 0.16, and the great majority were a few letters long. Mirror: PMC
- 25a.
Section Results and Discussion: “More than 70% of overlaps are <15 bp, and >85% are <30 bp”. “The slope of the linear relationship between the number of genes and overlapping gene pairs (overlaps) is 0.16”. ↩
- 25a.
- 26.
Fellner, Lea et al. “Evidence for the recent origin of a bacterial protein-coding, overlapping orphan gene by evolutionary overprinting.” BMC Evolutionary Biology 15:283. 2015. Laboratory study of nog1, a 101-codon frame lying wholly inside the gene citC on the opposite strand in a disease-causing E. coli strain: RNA sequencing in 11 conditions, promoter and protein-fusion tests, a mutant that stops nog1 without changing citC's protein, competition against the parent strain, and a survey of related bacteria. Mirror: PMC
- 26a.
Section Background: “In the extremely well-researched Escherichia coli, as far as we know, only three overlapping gene pairs have been described: htgA/yaaW [22, 23], yghW/morA [24], and tnpA/astA [25, 26].” ↩
- 26b.
Section Results: “Interestingly, Δnog1 shows a clear loss of fitness compared to its parental strain when grown in plain LB. When using MgCl2-supplemented broth, the decrease in fitness of Δnog1 was even more pronounced”. Figure 4b: Read from the figure: the mutant ends at about 21% of cells against the parent in plain LB and about 4% with MgCl2; the citC control mutant ends near 47%. ↩
- 26c.
Section Discussion: “(v) the gene product provides only a small fitness advantage for its carrier”. ↩ ↩
- 26a.
- 27.
Meydan, Sezen et al. “Retapamulin-assisted ribosome profiling reveals the alternative bacterial proteome.” Molecular Cell 74(3):481–493. 2019. Ribosome profiling of two E. coli strains with an antibiotic that freezes ribosomes at start codons, mapping start sites inside annotated genes. Most out-of-frame internal starts were mapped, not tested. Mirror: PMC
- 27a.
Section Results: “We limited our subsequent analysis to 124 iTISs conserved between the two strains, among which, 42 directed translation in frame with the main gene, whereas start codons of 74 iTISs were out of frame relative to the main ORF”. ↩
- 27a.
- 28.
Ardern, Zachary et al. “Are antisense proteins in prokaryotes functional?” Frontiers in Molecular Biosciences 7:187. 2020. Review, by members of the laboratory that described nog1, of the evidence for proteins read from the opposite strand of bacterial genes. It reverses the sign of the nog1 competition result reported in Fellner 2015. Mirror: PMC
- 28a.
Section Abstract: “Across bacterial species hundreds of short open reading frames antisense to annotated genes show evidence of both transcription and translation”. ↩
- 28b.
Section Phenotypes of Antisense Proteins: “A strand-specific deletion mutant has a growth advantage over the wildtype in LB, and a stronger advantage in medium supplemented with magnesium chloride ( Fellner et al., 2015 ).” ↩
- 28a.
- 29.
Liang, Han, and Laura F. Landweber. “A genome-wide study of dual coding regions in human alternatively spliced genes.” Genome Research 16(2):190–196. 2006. Scan of 2,585 well-annotated human alternatively spliced genes for coding regions read in two or three frames by different transcripts, with a comparison against mouse, rat and chicken to see whether the second frame is open in those species. Mirror: PMC
- 29a.
Section Results: “Among 2585 alternatively spliced genes, we found 173 genes that contain coding sequences for two reading frames”. Section Discussion: “The space of dual or multiple coding regions in the human genome is <0.5% of exons.” The Results add six genes coding in all three frames: 179 of 2,585, or 6.9%. ↩ ↩
- 29b.
Section Results: “one of two reading frames contains one or more stop codons in 80% of the dual coding regions in chickens, 45% in mice, and 54% in rats, respectively. These results indicate that most secondary reading frames in the dual coding regions examined evolved recently in mammals”. ↩
- 29a.
- 30.
Sanna, Chaitanya R. et al. “Overlapping genes in the human and mouse genomes.” BMC Genomics 9:169. 2008. Count of genes overlapping another gene among 13,484 human-mouse orthologous pairs, by strand and by whether coding sequence is shared, against a chromosome-matched random draw.
- 31.
Mudge, Jonathan M. et al. “Standardized annotation of translated open reading frames.” Nature Biotechnology 40(7):994–999. 2022. Consensus by GENCODE annotators and ribosome-profiling groups that merged seven published studies into one catalogue of 7,264 translated frames outside annotated proteins, with rules for naming and a small pilot search of mass-spectrometry data.
- 31a.
“We find single-mapping peptide-spectrum matches (PSMs) for 13 Ribo-seq ORFs ( Supplementary Table 11 ); all but one are supported by a single PSM each”. The 2022 pilot search, before the 2026 consortium reanalysis. ↩
- 31a.
- 32.
Chen, Jin et al. “Pervasive functional translation of noncanonical human open reading frames.” Science 367(6482):1140–1146. 2020. Ribosome profiling of human stem cells and fibroblasts followed by pooled gene-editing knockout screens of 2,353 non-canonical frames in stem cells and leukemia cells, from Jonathan Weissman's laboratory, with follow-up of a few hits. Mirror: PMC
- 32a.
Section Results: “In iPSCs, our screen identified over 500 ORF knockout “hits” that result in statistically significant phenotypes. The hits include 169 genes that are variants of annotated proteins, 78 “start overlap” hits, 230 uORF hits, 91 lncRNA CDS hits, and 2 downstream CDS hits.” “A fraction of the uORF hits do not have main, canonical CDSs with fitness defects upon knockout, suggesting an independent function of the uORFs or that disruption of the uORFs result in increases in main CDS expression that results in the growth phenotype”. ↩ ↩
- 32a.
- 33.
Prensner, John R. et al. “Noncanonical open reading frames encode functional proteins essential for cancer cell survival.” Nature Biotechnology 39(6):697–704. 2021. Knockout, overexpression and tagging tests of 553 candidate non-canonical frames, mostly in non-coding RNAs and upstream or downstream of known genes, in human cancer cell lines.
- 33a.
Section Abstract: “Here, we experimentally interrogated 553 candidates selected from non-canonical ORF datasets. Of these, 57 induced viability defects when knocked out in human cancer cell lines.” ↩
- 33a.
- 34.
Quelle, D. E. et al. “Alternative reading frames of the INK4a tumor suppressor gene encode two unrelated proteins capable of inducing cell cycle arrest.” Cell 83(6):993–1000. 1995. The discovery paper for ARF: cloning of a mouse INK4a transcript read in a second frame and tests showing its product stops the cell cycle.
- 34a.
Section Abstract: “An unrelated protein (p19ARF) arises in major part from an alternative reading frame of the mouse INK4a gene, and its ectopic expression in the nucleus of rodent fibroblasts induces G1 and G2 phase arrest. Economical reutilization of coding sequences in this manner is practically without precedent in mammalian genomes”. ↩ ↩
- 34a.
- 35.
Kim, Soo-Hyun et al. “Absence of p16INK4a and truncation of ARF tumor suppressors in chickens.” Proceedings of the National Academy of Sciences 100(1):211–216. 2003. Cloning and sequencing of the chicken INK4b-ARF-INK4a region and tests of the chicken proteins in cells. It found no p16INK4a in chickens and an ARF protein encoded by one exon only. Mirror: PMC
- 35a.
Section Abstract: “chickens do not encode an equivalent of INK4a. Moreover, the reading frame for chicken ARF does not extend into exon 2, because splicing occurs in a different register to that used in mammals.” Section Results: “the −1 reading frame in this exon is capable of encoding 162 amino acids that could have formed the carboxyl terminus of ARF, if splicing had occurred in the register used in mammals”. ↩ ↩
- 35a.
- 36.
Klemke, Martin et al. “Two overlapping reading frames in a single exon encode interacting proteins—a novel way of gene usage.” EMBO Journal 20(14):3849–3860. 2001. Discovery of ALEX, a protein read in a second frame of the XL exon of the rat Gnas gene, detected with antibodies in neuroendocrine cells and shown to bind the protein made from the first frame. Mirror: PMC
- 37.
Nekrutenko, Anton et al. “Oscillating evolution of a mammalian locus with overlapping reading frames: an XLalphas/ALEX relay.” PLoS Genetics 1(2):e18. 2005. Sequence comparison of the overlapping XLalphas and ALEX frames across eight primates and other mammals, measuring amino-acid change in each frame and asking whether one frame is held still while the other changes. Mirror: PMC
- 37a.
Section Discussion: “In each case the alternative reading frame is conserved in all known mammalian orthologs of a gene.” Said of the authors' unpublished list of candidate genes with second frames. ↩
- 37b.
Section Discussion: “Only three additional cases are known where protein products of both reading frames were biochemically characterized. These include genes for the cyclin D-dependent kinase inhibitor INK4a [ 28 ], X-box protein 1 [ 29 ], and a region of overlap between 4E-BP3 and MASK [ 30 ].” ↩
- 37c.
Section Results: “none of the ratios is significantly different from 1 (p-values from Fisher's exact test are between 0.14 and 0.77), supporting our hypothesis that the two proteins constantly co-evolve and maintain XL K A/ALEX K A of approximately 1”. ↩
- 37a.
- 38.
Vanderperre, Benoît et al. “An overlapping reading frame in the PRNP gene encodes a novel polypeptide distinct from the prion protein.” FASEB Journal 25(7):2373–2386. 2011. Laboratory study of a second frame in the prion protein gene, using tagged constructs from four species, antibodies against the second protein, and RNA silencing in human cells. A published letter and reply followed in 2012.
- 38a.
Section Abstract: “AltPrP is constitutively coexpressed with human, bovine, sheep, and deer PrP.” “By transfecting cells with siRNA against PrP mRNA, we repressed expression of both PrP and AltPrP, confirming endogenous expression of AltPrP from PRNP. AltPrP was also detected in human brain homogenate, primary neurons, and peripheral blood mononuclear cells.” ↩ ↩
- 38a.
- 39.
Bergeron, Danny et al. “An out-of-frame overlapping reading frame in the ataxin-1 coding sequence encodes a novel ataxin-1 interacting protein.” Journal of Biological Chemistry 288(30):21824–21835. 2013. Laboratory study of a second frame inside the human ataxin-1 coding sequence, using tagged constructs in cultured cells and antibodies against the second protein in human brain tissue. Mirror: PMC
- 39a.
Section Abstract: “We found an alternative translation initiation ATG codon in the +3 reading frame of human ATXN1 starting 30 nucleotides downstream of the initiation codon for ATXN1 and ending at nucleotide 587.” “Polyclonal antibodies raised against Alt-ATXN1 confirmed the expression of Alt-ATXN1 in human cerebellum expressing ATXN1.” ↩
- 39a.
- 40.
Firth, Andrew E. “Mapping overlapping functional elements embedded within the protein-coding regions of RNA viruses.” Nucleic Acids Research 42(20):12425–12439. 2014. Scan of sequence alignments for about 700 RNA virus species for coding windows with fewer silent changes than expected, tested first on known overlapping genes and known non-coding elements. It found hundreds of candidate regions, many not described before. Mirror: PMC
- 40a.
Section Results: “Synplot2 detected 20/21 overlapping genes (Table 1 ) and 16/17 overlapping non-coding features (Table 2 ).” ↩
- 40a.
- 41.
Mizokami, M. et al. “Constrained evolution with respect to gene overlap of hepatitis B virus.” Journal of Molecular Evolution 44(Suppl 1):S83–S90. 1997. Comparison of silent and amino-acid-changing substitutions in each gene of 27 hepatitis B virus strains, read here from the published abstract.
- 41a.
Section Abstract: “the number of synonymous substitutions per nucleotide site for the S gene was quite low and appeared a strongly constrained evolution. Because the P gene overlaps the S gene in a different frame, the low rate of synonymous substitution for the S gene can be explained by the evolutionary constraints which are imposed on the overlapping gene region. In other words, synonymous substitutions in the S gene can cause amino acid changes in its overlapping region in a different frame.” ↩
- 41a.
- 42.
Sabath, Niv et al. “Evolution of viral proteins originated de novo by overprinting.” Molecular Biology and Evolution 29(12):3767–3780. 2012. Phylogenetic study of overlapping virus genes, identifying 12 overprinted genes whose origin could be placed on the family tree, ranking them by relative age, and comparing their rates of change and the known effects of suppressing them. Mirror: PMC
- 42a.
Section Materials and Methods: “By manual examination of these trees, we identified 12 cases of de novo origin ( tables 1 and 3 ) that met the criterion of monophyly.” Section Discussion: “One important caveat of our analysis is that we are unable to estimate the absolute age of de novo genes but only their age relative to the divergence of a viral housekeeping gene that encodes the RNA-dependent RNA polymerase.” ↩ ↩
- 42b.
Section Results: “Mutations in the youngest de novo genes (taxa 1–3) have little to moderate effect, whereas suppression of older de novo genes (taxa 6 and 10–12) has severe effects.” Table 3. ↩ ↩
- 42a.
- 43.
Kovacs, Erika et al. “Dual coding in alternative reading frames correlates with intrinsic protein disorder.” Proceedings of the National Academy of Sciences 107(12):5429–5434. 2010. Computational study of 67 human genes whose splice variants read a stretch of at least 75 letters in two frames, with protein evidence for both forms, comparing predicted disorder in the older and newer frame. Mirror: PMC
- 43a.
Section Abstract: “we identified 67 human genes with alternative splice variants comprising a dual-coding region at least 75 nucleotides in length”. “in the case of -1 frameshifts, there is a significant tendency to become more disordered upon shifting the frame (16.7% in the ancestral frame, 56.3% in the derived frame)”. “The regions encoded by the derived frame are mostly disordered (disorder percentage > 50%) in 39 out of 62 cases”. For +1 shifts the paper gives 47.2% in the ancestral frame and 58.2% in the derived frame. ↩
- 43a.
- 44.
Willis, Sara, and Joanna Masel. “Gene birth contributes to structural disorder encoded by overlapping genes.” Genetics 210(1):303–313. 2018. Study of 92 overlapping virus gene pairs from 80 species, 39 of them with the older partner known from earlier work, and 150 non-overlapping controls, separating how much of the extra disorder in overlaps comes from gene birth and how much from the genetic code and selection, and testing which frame new genes are born in. Mirror: PMC
- 44a.
Section Abstract: “Overlapping gene regions encode higher levels of intrinsic structural disorder (ISD) than nonoverlapping genes (39% vs. 25% in our viral dataset).” Section Results: “By this calculation, birth facilitation accounts for 32 ± 13 % of the estimated total difference in ISD between overlapping and nonoverlapping genes.” ↩
- 44b.
Section Results: “This unexpected result is stronger for “internal overlaps,” in which one gene is completely contained within its overlapping partner (23 +1 events vs. 1 +2 event, P = 3 × 10 − 6 )”. “One artifact of the genetic code is that +1 frameshifts yield more start codons and fewer stop codons, and hence more and longer ORFs”. ↩ ↩
- 44a.
- 46.
Trifonov, Edward N. “The multiple codes of nucleotide sequences.” Bulletin of Mathematical Biology 51(4):417–432. 1989. Theoretical review arguing that DNA carries several overlapping codes besides the protein code, and proposing a regulatory role for tandem repeats. It is the source Sanford cites for "up to 12" codes; only its abstract is cited here.
- 46a.
Section Abstract: “Several codes of nucleotide sequences are discussed including: (1) the translation framing code, responsible for correct triplet counting by the ribosome during protein synthesis; (2) the chromatin code, which provides instructions on appropriate placement of nucleosomes along the DNA molecules and their spatial arrangement; (3) a putative loop code for single-stranded RNA-protein interactions. The codes are degenerate and corresponding messages are not only interspersed but actually overlap, so that some nucleotides belong to several messages simultaneously.” ↩ ↩
- 46a.
- 47.
Trifonov, Edward N. “Thirty years of multiple sequence codes.” Genomics, Proteomics & Bioinformatics 9(1–2):1–6. 2011. Short retrospective by the originator of the multiple-codes idea, listing the codes proposed since 1980 and arguing that sequence analysis should look for them. It gives no count of codes per letter. Mirror: PMC
- 47a.
“The conclusion from the above is obvious: one has to admit that the genetic sequences carry many different codes. If we are to know what the sequences are about, we have to detect and decipher these codes. The times of surrender to “junk” and “selfish DNA” are over, and “non-coding” becomes a misnomer.” ↩
- 47a.
- 48.
Itzkovitz, Shalev, and Uri Alon. “The genetic code is nearly optimal for allowing additional information within protein-coding sequences.” Genome Research 17(4):405–412. 2007. Computational comparison of the standard genetic code with large sets of alternative codes, measuring how well each lets coding sequences carry arbitrary short sequences and stop codons in shifted frames. The standard code ranks near the top. Mirror: PMC
- 48a.
Section Robustness to translational frame-shift errors: “One possibility is that the ability to include parallel codes within protein-coding sequences conferred a selection advantage during the early evolution of the genetic code. Alternatively, the genetic code might have been fixed in evolution before most parallel codes existed.” ↩ ↩
- 48a.
- 49.
Itzkovitz, Shalev et al. “Overlapping codes within protein-coding sequences.” Genome Research 20(11):1582–1589. 2010. Comparison of short-sequence counts in the coding regions of over 700 species against randomisations that keep amino-acid sequence and codon bias, to find signals that bias synonymous codon choice. Mirror: PMC
- 49a.
Section Abstract: “Our results suggest that genomes may have evolved to encode extensive overlapping information within protein-coding regions.” ↩
- 49a.
- 50.
Stergachis, Andrew B. et al. “Exonic transcription factor binding directs codon choice and affects protein evolution.” Science 342(6164):1367–1372. 2013. DNase footprinting of protein binding inside human exons in 81 cell types, from John Stamatoyannopoulos's laboratory, reporting how many codons lie under a regulatory-protein footprint. Its conclusion about codon constraint was disputed in 2015. Mirror: PMC
- 51.
Calvo, Sarah E. et al. “Upstream open reading frames cause widespread reduction of protein expression and are polymorphic among humans.” Proceedings of the National Academy of Sciences 106(18):7507–7512. 2009. Count of upstream frames in human and mouse transcripts, reporter tests of 25 of them, and a survey of human variants that create or remove them.
- 51a.
Section Abstract: “Using reporter constructs to test 25 selected uORFs, we estimate that uORFs typically reduce protein expression by 30–80%”. Section Results: “we find that 49% of human and 44% of mouse transcripts contain at least 1 uORF”. ↩
- 51a.
- 52.
Sanford, John C. Genetic Entropy & the Mystery of the Genome. Elim Publishing. 2005. First edition of Sanford's book. Its overlapping-codes argument already appears here, before the 2013 model and the ENCODE claims added in later editions.
- 52a.
Page 132: “For example, most human coding sequences encode for two different RNAs, read in opposite directions (i.e. both DNA strands are transcribed - Yelin et al., 2003).” Page 133: “Trifonov (1989), has shown that probably all DNA sequences in the genome encrypt multiple "codes" (up to 12 codes). In computer science, this type of "data compression" can only result from the highest level of information design, and results in maximal information density.” The 2005 edition already makes the twelve-code claim; the later 2013 probability model is not its source. ↩
- 52a.
- 53.
Wells, Jonathan. The Myth of Junk DNA. Discovery Institute Press. 2011. Book by a Discovery Institute biologist arguing that research since 2000 has found function across the non-coding genome and that this fits intelligent design better than evolution. Its overlapping-codes section rests on two 2007 papers.
- 53a.
Page 105: ““overlapping coding regions,” a feature that the researchers concluded “is nearly impossible by chance.”” “The Israeli scientists concluded that the specification of amino acids by three-nucleotide “codons” in DNA is “nearly optimal for allowing additional information within protein-coding sequences.”” The first quotation's words match Chung et al. 2007, whose sentence reads "although dual coding is nearly impossible by chance, a number of human transcripts contain overlapping coding regions." The Israeli scientists are Itzkovitz and Alon 2007, named in Wells's endnote. ↩
- 53b.
Page 106: ““If these systems are designed,” he wrote in 2004, “we can expect the information to be densely packed and multilayered.” Thus “dense, multi-layered embedding of information is a prediction of intelligent design.” "He" is William Dembski. ↩
- 53a.
- 54.
Rana, Fazale. “Overlapping genes evince intelligent design.” Reasons to Believe. Jul 26, 2007. Short article by the biochemist of an old-earth design ministry, written on the 2007 dual-coding study.
- 54a.
“The apparent widespread occurrence of overlapping genes doesn't make much sense from an evolutionary perspective, because of the cost they represent to the organism. This cost is only worth it if there is a rationale for overlapping genes.” “For me, the only explanation for overlapping genes is the work of a Creator. It's hard to imagine how undirected evolutionary processes could produce overlapping genes.” ↩
- 54a.
- 55.
Tomkins, Jeffrey P. “Duons: parallel gene code defies evolution.” Institute for Creation Research, Creation Science Update. Jan 6, 2014. Short news commentary by an ICR geneticist on the 2013 binding-site paper.
- 55a.
Paragraph 1: “Researchers have just characterized a new, previously hidden genetic code embedded within the same sections of genes that code for proteins—utterly defying all naturalistic explanations for its existence.” “Some genes overlap parts of other genes in the genome, and now it has been revealed that many genes have areas that contain dual codes within the very same sequence.” The second sentence's citations are Stergachis 2013 and Sanna 2008. ↩
- 55a.
- 56.
Gauger, Ann K. et al. “Evidence for human uniqueness.” Theistic Evolution: A Scientific, Philosophical, and Theological Critique. Crossway. Edited by J. P. Moreland, Stephen C. Meyer, Christopher Shaw, Ann K. Gauger and Wayne Grudem. 2017. Chapter in a design-movement volume answering the genetic arguments for human-chimpanzee common ancestry, including the argument from shared codons.
- 56a.
Page 500: “In addition, this is a strong argument for design: the sophistication required to have codes within codes would not arise by simple mutation and natural selection.” ↩
- 56a.
- 57.
Luskin, Casey. “30 years of intelligent design predicting the unexpected genome.” Science and Culture Today. Aug 11, 2026. Article from the Discovery Institute's news site listing genome findings it counts as successful predictions of intelligent design.
- 57a.
“Another discovery in genomics predicted by ID is the presence of overlapping genes, where the same stretch of DNA encodes multiple genes. One strand of DNA might encode one protein while the opposite strand encodes an entirely different functional protein (see, e.g., Adelman et al., 1987 or Mehravar et al., 2021). Such an astonishing degree of information compression is only expected if a mind of genius encoded the information in DNA, not random mutation and blind natural selection.” ↩
- 57a.
- 58.
Yelin, Rodrigo et al. “Widespread occurrence of antisense transcription in the human genome.” Nature Biotechnology 21(4):379–386. 2003. Search of public expressed-sequence databases for RNAs from both strands of the same human locus, checked with strand-specific microarrays; cited here from its abstract.
- 58a.
Section Abstract: “The resulting data set of 2,667 sense-antisense pairs was evaluated by microarrays containing strand-specific oligonucleotide probes derived from the region of overlap.” “We conclude that > or =60% of this data set, or approximately 1,600 predicted sense-antisense transcriptional units, are transcribed from both DNA strands.” ↩
- 58a.
- 59.
Kapranov, Philipp et al. “Genome-wide transcription and the implications for genomic organization.” Nature Reviews Genetics 8(6):413–423. 2007. Review from the Affymetrix group that ran early tiling-array transcription maps, arguing that much of the genome is transcribed; cited here from its abstract.
- 59a.
Section Abstract: “Evidence indicates that most of both strands of the human genome might be transcribed, implying extensive overlap of transcriptional units and regulatory elements.” ↩
- 59a.
- 60.
Barrell, B. G. et al. “Overlapping genes in bacteriophage φX174.” Nature 264:34–41. 1976. The first report of overlapping protein-coding genes, from Frederick Sanger's laboratory during the sequencing of phiX174.
- 60a.
Section Abstract: “Bacteriophage φX174 genes D and E are translated from the same DNA sequence but in different reading frames.” ↩
- 60a.
- 61.
Adelman, J. P. et al. “Two mammalian genes transcribed from opposite strands of the same DNA locus.” Science 235(4795):1514–1517. 1987. Description of a rat locus where the gonadotropin-releasing hormone gene and a second gene are transcribed from opposite strands; cited here from its abstract.
- 61a.
Section Abstract: “A second gene, SH, is transcribed from the other DNA strand to produce RNA of undefined function.” ↩
- 61a.
- 62.
Mehravar, Majid et al. “Exon and intron sharing in opposite direction—an undocumented phenomenon in human genome—between Pou5f1 and Tcf19 genes.” BMC Genomics 22:718. 2021. Computational re-analysis of earlier expression data describing one human gene pair whose splice variants share two exons and an intron on opposite strands. Mirror: PMC
- 62a.
Section Abstract: “However, the advantage and biological significance of overlapping genes is still unclear.” ↩
- 62a.
- 63.
Recomputed for this article by
research/montanez-check.ts(output inresearch/montanez-check-output.txt) from the chapter's equation (1): the sum over k from 7 to 12 of C(12,k) pk(1-p)12-k with p = 0.001/3 gives 3.616 × 10-22, or 71.2 bits. With p = 0.00034 the same sum gives 4.154 × 10-22, the printed value. The two-code values printed on page 146, 1.11 × 10-7 and 6.66 × 10-4, reproduce exactly. ↩