Description

In the decades since the C. elegans genome sequence was completed (C. elegans Sequencing Consortium, 1998), ongoing curation efforts continue to further improve and annotate the genome (Harris et al., 2020). In the course of searching for alternatively-spliced genes using RNA-seq data, we happened upon an alternative splicing event in a region of the C. elegans genome with no associated gene or feature annotation. Upon closer inspection, we identified RNA-seq support for an unannotated four-exon gene, with the third exon being alternatively spliced (Fig 1A). The putative coding sequences of this gene encode proteins of either 86 or 122 amino acids, depending on the splicing choice. The longer isoform (exon three included) is the major isoform according to our RNA-seq data (Fig 1B), and this isoform encodes four predicted transmembrane domains (Fig 1C). RT-PCR of wild-type L4 stage worms further confirms that the gene is expressed at detectable levels, and that exon three is alternatively spliced (Fig 1D). We have submitted this information to WormBase, and the new gene is set to appear as WBGene00306123 (transcripts Y97E10C.2a and Y97E10C.2b) in release WS283.

While inspecting this locus, we also noted RNA-seq evidence for the gene immediately upstream (T05B11.7) that deviates from its annotated gene model. The RNA-seq data indicate that prior to reaching the annotated stop codon in exon three, the mRNA is spliced to an unannotated downstream exon (Fig 1E). Therefore, the annotated stop codon is in fact located within an intron, and the T05B11.7 gene contains four exons. This increases the coding sequence of T05B11.7 from 852 to 924 nucleotides. We also found evidence for alternative splicing in the T05B11.7 gene. Exon one harbors an alternative 5’ splice site which is used in a small fraction of the transcripts detected by RNA-seq (Fig 1F) and which is also detectable via RT-PCR (Fig 1G). These alternative isoforms (including the new addition of exon four) will appear in WS283 as transcripts T05B11.7a and T05B11.7b.

Finally, the gene T05B11.7 is currently annotated as part of an operon (CEOP5584) along with its upstream gene clic-1 (Fig 1H). Our newly-described gene Y97E10C.2, which is located immediately downstream of T05B11.7, is now included as the final gene in this operon, making CEOP5584 a 3-gene operon (Fig 1H).