Widespread organisation of C. elegans genes into operons: fact or function?

Long-Read Sequencing of Human Cytomegalovirus Transcriptome Reveals RNA Isoforms Carrying Distinct Coding Potentials

The human cytomegalovirus (HCMV) is a ubiquitous, human pathogenic herpesvirus. The complete viral genome is transcriptionally active during infection; however, a large part of its transcriptome has yet to be annotated. In this work, we applied the amplified isoform sequencing technique from Pacific Biosciences to characterize the lytic transcriptome of HCMV strain Towne varS. We developed a pipeline for transcript annotation using long-read sequencing data. We identified 248 transcriptional start sites, 116 transcriptional termination sites and 80 splicing events. Using this information, we have annotated 291 previously undescribed or only partially annotated transcript isoforms, including eight novel antisense transcripts and their isoforms, as well as a novel transcript (RS2) in the short repeat region, partially antisense to RS1. Similarly to other organisms, we discovered a high transcriptional diversity in HCMV, with many transcripts only slightly differing from one another. Comparing our transcriptome profiling results to an earlier ribosome footprint analysis, we have concluded that the majority of the transcripts contain multiple translationally active ORFs, and also that most isoforms contain unique combinations of ORFs. Based on these results, we propose that one important function of this transcriptional diversity may be to provide a regulatory mechanism at the level of translation.

Trans-splicing enhances translational efficiency in C. elegans

Translational efficiency is subject to extensive regulation. However, the factors influencing such regulation are poorly understood. In Caenorhabditis elegans, 62% of genes are trans-spliced to a specific spliced leader (SL1), which replaces part of the native 5' untranslated region (5' UTR). Given the pivotal role the 5' UTR plays in the regulation of translational efficiency, we hypothesized that SL1 trans-splicing functions to regulate translational efficiency. With genome-wide analysis on Ribo-seq data, polysome profiling experiments, and CRISPR-Cas9-based genetic manipulation of trans-splicing sites, we found four lines of evidence in support of this hypothesis. First, SL1 trans-spliced genes have higher translational efficiencies than non-trans-spliced genes. Second, SL1 trans-spliced genes have higher translational efficiencies than non-trans-spliced orthologous genes in other nematode species. Third, an SL1 trans-spliced isoform has higher translational efficiency than the non-trans-spliced isoform of the same gene. Fourth, deletion of trans-splicing sites of endogenous genes leads to reduced translational efficiency. Importantly, we demonstrated that SL1 trans-splicing plays a key role in enhancing translational efficiencies of essential genes. We further discovered that SL1 trans-splicing likely enhances translational efficiency by shortening the native 5' UTRs, hence reducing the presence of upstream start codons (uAUG) and weakening mRNA secondary structures. Taken together, our study elucidates the global function of trans-splicing in enhancing translational efficiency in nematodes, paving the way for further understanding the genomic mechanisms of translational regulation.

Genome-wide analysis of trans-splicing in the nematode Pristionchus pacificus unravels conserved gene functions for germline and dauer development in divergent operons

Discovery of trans-splicing in multiple metazoan lineages led to the identification of operon-like gene organization in diverse organisms, including trypanosomes, tunicates, and nematodes, but the functional significance of such operons is not completely understood. To see whether the content or organization of operons serves similar roles across species, we experimentally defined operons in the nematode model Pristionchus pacificus. We performed affinity capture experiments on mRNA pools to specifically enrich for transcripts that are trans-spliced to either the SL1- or SL2-spliced leader, using spliced leader-specific probes. We obtained distinct trans-splicing patterns from the analysis of three mRNA pools (total mRNA, SL1 and SL2 fraction) by RNA-seq. This information was combined with a genome-wide analysis of gene orientation and spacing. We could confirm 2219 operons by RNA-seq data out of 6709 candidate operons, which were predicted by sequence information alone. Our gene order comparison of the Caenorhabditis elegans and P. pacificus genomes shows major changes in operon organization in the two species. Notably, only 128 out of 1288 operons in C. elegans are conserved in P. pacificus. However, analysis of gene-expression profiles identified conserved functions such as an enrichment of germline-expressed genes and higher expression levels of operonic genes during recovery from dauer arrest in both species. These results provide support for the model that a necessity for increased transcriptional efficiency in the context of certain developmental processes could be a selective constraint for operon evolution in metazoans. Our method is generally applicable to other metazoans to see if similar functional constraints regulate gene organization into operons.

Characterizing the transcriptional regulation of let-721, a Caenorhabditis elegans homolog of human electron flavoprotein dehydrogenase

LET-721 is the Caenorhabditis elegans ortholog of electron-transferring flavoprotein dehydrogenase (ETFDH). We are studying this protein in C. elegans in order to establish a tractable model system for further exploration of ETFDH structure and function. ETFDH is an inner mitochondrial membrane localized enzyme that plays a key role in the beta-oxidation of fatty acids and catabolism of amino acids and choline. ETFDH accepts electrons from at least twelve mitochondrial matrix flavoprotein dehydrogenases via an intermediate dimer protein and transfers the electrons to ubiquinone. In humans, ETFDH mutations result in the autosomal recessive metabolic disorder, multiple acyl-CoA dehydrogenase deficiency. Mutants of let-721 in C. elegans are either maternal effect lethals or semi-sterile. let-721 is transcribed in the pharynx, body wall muscle, hypoderm, intestine and somatic gonad. In addition, the subcellular localization of LET-721 agrees with predictions that it is localized to mitochondria. We identified and confirmed three cis-regulatory sequences (pha-site, rep-site, and act-site). Phylogenetic footprinting of each site indicates that they are conserved between four Caenorhabditis species. The pha-site mapped roughly 1,300 bp upstream of let-721's translational start site and is necessary for expression in pharyngeal tissues. The rep-site mapped roughly 830 bp upstream of the translational start site and represses expression of LET-721 within pharyngeal tissues. The act-site mapped roughly 800 bp upstream of the translational start site and is required for expression within spermatheca, body wall muscle, pharynx, and intestine. Taken together, we find that LET-721 is a mitochondrially expressed protein that is under complex transcriptional controls.

Evolutionary dynamics of nematode operons: easy come, slow go

Operons are widespread in prokaryotes, but are uncommon in eukaryotes, except nematode worms, where approximately 15% of genes reside in over 1100 operons in the model organism Caenorhabditis elegans. It is unclear how operons have become abundant in nematode genomes. The "one-way street" hypothesis asserts that once formed by chance, operons are very difficult to break, because the breakage would leave downstream genes in an operon without a promoter, and hence, unexpressed. To test this hypothesis, we analyzed the presence and absence of C. elegans operons in Caenorhabditis briggsae, Caenorhabditis remanei, and Caenorhabditis brenneri, using Pristionchus pacificus and Brugia malayi as outgroups, and identified numerous operon gains and losses. Coupled with experimental examination of trans-splicing patterns, our comparative genomic analysis revealed diverse molecular mechanisms of operon losses, including inversion, insertion, and relocation, but the presence of internal promoters was not found to facilitate operon losses. In several cases, the data allowed inference of mechanisms by which downstream genes are expressed after operon breakage. We found that the rate of operon gain is approximately 3.3 times that of operon loss. Thus, the evolutionary dynamics of nematode operons is better described as "easy come, slow go," rather than a "one-way street." Based on a mathematic model of operon gains and losses and additional assumptions, we projected that the number of operons in C. elegans will continue to rise by 6%-18% in future evolution before reaching equilibrium between operon gains and losses.

Characterization of gana-1, a Caenorhabditis elegans gene encoding a single ortholog of vertebrate alpha-galactosidase and alpha-N-acetylgalactosaminidase

Background

Human α-galactosidase A (α-GAL) and α-N-acetylgalactosaminidase (α-NAGA) are presumed to share a common ancestor. Deficiencies of these enzymes cause two well-characterized human lysosomal storage disorders (LSD) – Fabry (α-GAL deficiency) and Schindler (α-NAGA deficiency) diseases. Caenorhabditis elegans was previously shown to be a relevant model organism for several late endosomal/lysosomal membrane proteins associated with LSDs. The aim of this study was to identify and characterize C. elegans orthologs to both human lysosomal luminal proteins α-GAL and α-NAGA.

Results

BlastP searches for orthologs of human α-GAL and α-NAGA revealed a single C. elegans gene (R07B7.11) with homology to both human genes (α-galactosidase and α-N-acetylgalactosaminidase) – gana-1. We cloned and sequenced the complete gana-1 cDNA and elucidated the gene organization.

Phylogenetic analyses and homology modeling of GANA-1 based on the 3D structure of chicken α-NAGA, rice α-GAL and human α-GAL suggest a close evolutionary relationship of GANA-1 to both human α-GAL and α-NAGA.

Both α-GAL and α-NAGA enzymatic activities were detected in C. elegans mixed culture homogenates. However, α-GAL activity on an artificial substrate was completely inhibited by the α-NAGA inhibitor, N-acetyl-D-galactosamine.

A GANA-1::GFP fusion protein expressed from a transgene, containing the complete gana-1 coding region and 3 kb of its hypothetical promoter, was not detectable under the standard laboratory conditions. The GFP signal was observed solely in a vesicular compartment of coelomocytes of the animals treated with Concanamycin A (CON A) or NH₄Cl, agents that increase the pH of the cellular acidic compartment.

Immunofluorescence detection of the fusion protein using polyclonal anti-GFP antibody showed a broader and coarsely granular cytoplasmic expression pattern in body wall muscle cells, intestinal cells, and a vesicular compartment of coelomocytes.

Inhibition of gana-1 by RNA interference resulted in a decrease of both α-GAL and α-NAGA activities measured in mixed stage culture homogenates but did not cause any obvious phenotype.

Conclusions

GANA-1 is a single C. elegans ortholog of both human α-GAL and α-NAGA proteins. Phylogenetic, homology modeling, biochemical and GFP expression analyses support the hypothesis that GANA-1 has dual enzymatic activity and is localized in an acidic cellular compartment.