The biochemical composition of coronavirus Nucleocapsid proteins drives the bicistronic properties of their sub-genomic mRNAs

The nucleocapsid (N) gene of coronaviruses (CoVs) encodes a structural RNA-binding protein required for viral assembly. In SARS- and MERS-CoV viruses, reading frame 2 of N contains an overlapping open reading frame (ORF) encoding an innate immune antagonist accessory protein. However, the extent of this dual coding capacity across the CoV family, and the genetic features driving overlapping ORF formation, remain unexamined. Here, we developed an in-silico tool (ORFan) to screen structural gene sub-genomic mRNA sequences from 152 species of coronavirus for the presence of overlapping accessory gene candidates. We found that frame 2 of the N gene in Alpha-, Beta- and Delta-CoVs frequently contains a potential accessory ORF (AcORF), conserved at the subgenus level, unlike the spike, envelope and membrane genes. We also show that CoV nucleotide biases in conjunction with codon biases integral to N protein function drive the relatively high abundances of start codons in frame two and stop codons in frame three, explaining the tendency towards frame 2 AcORFs. We further suggest that viral nucleotide biases and the codon composition required for RNA-binding function facilitate overlapping ORF formation in other virus families