Transcriptome assembly, gene annotation and tissue gene expression atlas of the rainbow trout

Efforts to obtain a comprehensive genome sequence for rainbow trout are ongoing and will be complemented by transcriptome information that will enhance genome assembly and annotation. Previously, transcriptome reference sequences were reported using data from different sources. Although the previous work added a great wealth of sequences, a complete and well-annotated transcriptome is still needed. In addition, gene expression in different tissues was not completely addressed in the previous studies. In this study, non-normalized cDNA libraries were sequenced from 13 different tissues of a single doubled haploid rainbow trout from the same source used for the rainbow trout genome sequence. A total of ~1.167 billion paired-end reads were de novo assembled using the Trinity RNA-Seq assembler yielding 474,524 contigs > 500 base-pairs. Of them, 287,593 had homologies to the NCBI non-redundant protein database. The longest contig of each cluster was selected as a reference, yielding 44,990 representative contigs. A total of 4,146 contigs (9.2%), including 710 full-length sequences, did not match any mRNA sequences in the current rainbow trout genome reference. Mapping reads to the reference genome identified an additional 11,843 transcripts not annotated in the genome. A digital gene expression atlas revealed 7,678 housekeeping and 4,021 tissue-specific genes. Expression of about 16,000–32,000 genes (35–71% of the identified genes) accounted for basic and specialized functions of each tissue. White muscle and stomach had the least complex transcriptomes, with high percentages of their total mRNA contributed by a small number of genes. Brain, testis and intestine, in contrast, had complex transcriptomes, with a large numbers of genes involved in their expression patterns. This study provides comprehensive de novo transcriptome information that is suitable for functional and comparative genomics studies in rainbow trout, including annotation of the genome.

The C. elegans healthspan and stress-resistance assay toolkit

A wealth of knowledge on the genetic mechanisms that govern aging has emerged from the study of mutants that exhibit enhanced longevity and exceptional resilience to adverse environmental conditions. In these studies, lifespan has been an excellent proxy for establishing the rate of aging, but it is not always correlated with qualitative measures of healthy aging or 'healthspan'. Although the attributes of healthspan have been challenging to define, they share some universal features that are increasingly being incorporated into aging studies. Here we describe methods used to determine Caenorhabditis elegans healthspan. These include assessments of tissue integrity and functionality and resistance to a variety of biotic and abiotic stressors. We have chosen to include simple, rapid assays in this collection that can be easily undertaken in any C. elegans laboratory, and can be relied on to provide a preliminary but thorough insight into the healthspan of a population.

The significance of alternative transcripts for Caenorhabditis elegans transcription factor genes, based on expression pattern analysis

Background

Sequence-specific DNA-binding proteins, with their paramount importance in the regulation of expression of the genetic material, are encoded by approximately 5% of the genes in an animal’s genome. But it is unclear to what extent alternative transcripts from these genes may further increase the complexity of the transcription factor complement.

Results

Of the 938 potential C. elegans transcription factor genes, 197 were annotated in WormBase as encoding at least two distinct isoforms. Evaluation of prior evidence identified, with different levels of confidence, 50 genes with alternative transcript starts, 23 with alternative transcript ends, 35 with alternative splicing and 34 with alternative transcripts generated by a combination of mechanisms, leaving 55 that were discounted. Expression patterns were determined for transcripts for a sample of 29 transcription factor genes, concentrating on those with alternative transcript starts for which the evidence was strongest. Seamless fosmid recombineering was used to generate reporter gene fusions with minimal modification to assay expression of specific transcripts while maintaining the broad genomic DNA context and alternative transcript production. Alternative transcription factor gene transcripts were typically expressed with identical or substantially overlapping distributions rather than in distinct domains.

Conclusions

Increasingly sensitive sequencing technologies will reveal rare transcripts but many of these are clearly non-productive. The majority of the transcription factor gene alternative transcripts that are productive may represent tolerable noise rather than encoding functionally distinct isoforms.

The Caenorhabditis elegans homeobox gene ceh-19 is required for MC motorneuron function

Simplicity has made C. elegans pharyngeal development a particularly well-studied subject. Nevertheless, here we add the previously uncharacterized homeobox gene F20D12.6/ceh-19 to the set of transcription factor genes involved. GFP reporter assays revealed that ceh-19 is expressed in three pairs of neurons, the pharyngeal pace-maker neurons MC, the amphid neurons ADF and the phasmid neurons PHA. ceh-19(tm452) mutants are viable and fertile, but grow slightly slower, produce less progeny over a prolonged period, and live longer than the wild type. These phenotypes are likely due to the moderately reduced pharyngeal pumping speed arising from the impairment of MC activity. MC neurons are still born in the ceh-19 mutants but display various morphological defects. ceh-19 expression in MC is completely lost in progeny from animals subject to RNAi for pha-4, which encodes an organ-specifying forkhead transcription factor. CEH-19 is required for the activation in MCs of the excitatory FMRFamide-like neuropeptide-encoding gene flp-2. A regulatory pathway from pha-4 through ceh-19 to flp-2 is thereby defined. The resilience of MC identity in the absence of CEH-19 may reflect the buffering qualities of transcription factor regulatory networks. genesis 51:163–178, 2013. © 2013 Wiley Periodicals, Inc.

A simplified counter-selection recombineering protocol for creating fluorescent protein reporter constructs directly from C. elegans fosmid genomic clones

BACKGROUND

Recombineering is a genetic engineering tool that enables facile modification of large episomal clones, e.g. BACs, fosmids. We have previously adapted this technology to generate, directly from fosmid-based genomic clones, fusion gene reporter constructs designed to investigate gene expression patterns in C. elegans. In our adaptation a rpsL-tet(A) positive/negative-selection cassette (RT-cassette) is first inserted and then, under negative selection, seamlessly replaced with the desired sequence. We report here on the generation and application of a resource comprising two sets of constructs designed to facilitate this particular recombineering approach.

RESULTS

Two complementary sets of constructs were generated. The first contains different fluorescent protein reporter coding sequences and derivatives while the second set of constructs, based in the copy-number inducible vector pCC1Fos, provide a resource designed to simplify RT-cassette-based recombineering. These latter constructs are used in pairs the first member of which provides a template for PCR-amplification of an RT-cassette while the second provides, as an excised restriction fragment, the desired fluorescent protein reporter sequence. As the RT-cassette is flanked by approximately 200 bp from the ends of the reporter sequence the subsequent negative selection replacement step is highly efficient. Furthermore, use of a restriction fragment minimizes artefacts negating the need for final clone sequencing. Utilizing this resource we generated single-, double- and triple-tagged fosmid-based reporters to investigate expression patterns of three C. elegans genes located on a single genomic clone.

CONCLUSIONS

We describe the generation and application of a resource designed to facilitate counter-selection recombineering of fosmid-based C. elegans genomic clones. By choosing the appropriate pair of 'insertion' and 'replacement' constructs recombineered products, devoid of artefacts, are generated at high efficiency. Gene expression patterns for three genes located on the same genomic clone were investigated via a set of fosmid-based reporter constructs generated with the modified protocol.

A survey of putative secreted and transmembrane proteins encoded in the C. elegans genome

BACKGROUND

Almost half of the Caenorhabditis elegans genome encodes proteins with either a signal peptide or a transmembrane domain. Therefore a substantial fraction of the proteins are localized to membranes, reside in the secretory pathway or are secreted. While these proteins are of interest to a variety of different researchers ranging from developmental biologists to immunologists, most of secreted proteins have not been functionally characterized so far.

RESULTS

We grouped proteins containing a signal peptide or a transmembrane domain using various criteria including evolutionary origin, common domain organization and functional categories. We found that putative secreted proteins are enriched for small proteins and nematode-specific proteins. Many secreted proteins are predominantly expressed in specific life stages or in one of the two sexes suggesting stage- or sex-specific functions. More than a third of the putative secreted proteins are upregulated upon exposure to pathogens, indicating that a substantial fraction may have a role in immune response. Slightly more than half of the transmembrane proteins can be grouped into broad functional categories based on sequence similarity to proteins with known function. By far the largest groups are channels and transporters, various classes of enzymes and putative receptors with signaling function.

CONCLUSION

Our analysis provides an overview of all putative secreted and transmembrane proteins in C. elegans. This can serve as a basis for selecting groups of proteins for large-scale functional analysis using reverse genetic approaches.

A regulatory cascade of three transcription factors in a single specific neuron, DVC, in Caenorhabditis elegans

Homeobox proteins are critical regulators of developmental gene transcription and cell specification. Many insights into transcriptional regulation have been gained from studies in the nematode Caenorhabditis elegans. We investigated the expression and regulation of the C. elegans homeobox gene ceh-63, which encodes a single-homeodomain transcription factor of 152 amino acids. ceh-63 is expressed in the interneuron DVC in both sexes, from late embryogenesis through adulthood, and two pairs of uterine cells in reproductive hermaphrodites only. A reporter gene fusion, encoding GFP fused to the full-length CEH-63, also drove weak inconsistent expression in additional unidentified cells in the head and tail. A potential ceh-63 null mutant had no obvious abnormalities, except for a possible increase in subtle defects of the DVC axon projection. No behavioural responses were observed upon either laser ablation of DVC or activation of DVC through light stimulation of channelrhodopsin-2 specifically expressed in this neuron. The function of DVC therefore remains enigmatic. A transcriptional regulatory cascade operating in DVC was defined from the LIM-homeodomain protein CEH-14 through CEH-63 to the helix-turn-helix transcription factor MBR-1. Both CEH-14 and CEH-63 individually bound the mbr-1 promoter in a yeast one-hybrid assay. A model is proposed suggesting that CEH-14 activates ceh-63 and then along with CEH-63 co-ordinately activates mbr-1.

Distinct mechanisms for delimiting expression of four Caenorhabditis elegans transcription factor genes encoding activators or repressors

Regulatory transcription factors operate in networks, conferring biological robustness that makes dissection of such gene control processes difficult. The nematode Caenorhabditis elegans is a powerful molecular genetic system that allows the close scrutiny needed to understand these processes in an animal, in vivo. Strikingly lower levels of gene expression were observed when a gfp reporter was inserted into C. elegans transcription factor genes, in their broader genomic context, in comparison to when the reporter was fused to just the promoter regions. The lower level of expression is more consistent with endogenous levels of the gene products, based on independent protein and transcript assays. Through successive precise manipulations of the reporter fusion genes, elements essential for the lower level of expression were localised to the protein-coding region. With a closer focus on four transcription factor genes, the expression of both genes encoding transcriptional activators was found to be restricted by a post-transcriptional mechanism while expression of both genes encoding transcriptional repressors was delimited by transcriptional repression. An element through which the transcriptional repression acts for unc-4 was localised to a 30 base-pair region of a protein-encoding exon, with potentially wider implications for how homeobox genes operate. The hypothesis that the distinction in mechanisms delimiting expression of the two types of transcription factor genes, as observed here, may apply more widely is raised. This leads to observations concerning the implications of these different mechanisms on stochastic noise in gene expression and the consequent significance for developmental decisions in general.

The production of C. elegans transgenes via recombineering with the galK selectable marker

The creation of transgenic animals is widely utilized in C. elegans research including the use of GFP fusion proteins to study the regulation and expression pattern of genes of interest or generation of tandem affinity purification (TAP) tagged versions of specific genes to facilitate their purification. Typically transgenes are generated by placing a promoter upstream of a GFP reporter gene or cDNA of interest, and this often produces a representative expression pattern. However, critical elements of gene regulation, such as control elements in the 3' untranslated region or alternative promoters, could be missed by this approach. Further only a single splice variant can be usually studied by this means. In contrast, the use of worm genomic DNA carried by fosmid DNA clones likely includes most if not all elements involved in gene regulation in vivo which permits the greater ability to capture the genuine expression pattern and timing. To facilitate the generation of transgenes using fosmid DNA, we describe an E. coli based recombineering procedure to insert GFP, a TAP-tag, or other sequences of interest into any location in the gene. The procedure uses the galK gene as the selection marker for both the positive and negative selection steps in recombineering which results in obtaining the desired modification with high efficiency. Further, plasmids containing the galK gene flanked by homology arms to commonly used GFP and TAP fusion genes are available which reduce the cost of oligos by 50% when generating a GFP or TAP fusion protein. These plasmids use the R6K replication origin which precludes the need for extensive PCR product purification. Finally, we also demonstrate a technique to integrate the unc-119 marker on to the fosmid backbone which allows the fosmid to be directly injected or bombarded into worms to generate transgenic animals. This video demonstrates the procedures involved in generating a transgene via recombineering using this method.

Specialized chromosomes and their uses in Caenorhabditis elegans

Research on Caenorhabditis elegans involves the use of a wide range of genetic and molecular tools consisting of chromosomal material captured and modified for specific purposes. These "specialized chromosomes" come in many forms ranging from relatively simple gene deletions to complex rearrangements involving endogenous chromosomes as well as transgenic constructs. In this chapter, we describe the specialized chromosomes that are available in C. elegans, their origins, practical considerations, and methods for generation and evaluation. We will summarize their uses for biological studies, and their contribution to our knowledge about chromosome biology.

Elucidating gene regulatory mechanisms for sperm function through the integration of classical and systems approaches in C. elegans

From worms to mammals, successful spermatogenesis depends on a gene expression profile that balances activating and repressive mechanisms. Besides developmental control of specific spermatogenic genes, male fertility requires temporal shifts in global gene expression and dramatic changes in chromatin structure and condensation. Recent studies are beginning to elucidate the molecular processes that both drive these temporal changes in gene expression and underlie fertility. In this review, we provide an overview of relevant C. elegans studies that have laid the groundwork for modern approaches. Next, we highlight recent studies that investigate how gene expression in C. elegans is modulated during spermatogenesis. These studies use large-scale genomic profiling in combination with bioinformatics, genetics, biochemistry, and in vitro methods to target specific stages or processes during sperm formation. Such studies are beginning to elucidate the multiple layers of gene regulation required during spermatogenesis, i.e., transcriptional, post-transcriptional, and epigenetic. Moreover, knowledge of how C. elegans coordinately regulates gene expression during spermatogenesis promises to provide key insights into parallel processes in mammals that are vital for fertility.

Escherichia coli MW005: lambda Red-mediated recombineering and copy-number induction of oriV-equipped constructs in a single host

Background

Escherichia coli strain EL350 contains chromosomally integrated phage lambda Red recombinase genes enabling this strain to be used for modifying the sequence of resident clones via recombineering. BAC and fosmid clones are highly suitable for modification by recombineering but, because they are present at low (1-2) copies per cell, the DNA is difficult to isolate in high yield and purity. To overcome this limitation vectors, e.g. pCC1FOS, have been constructed that contain the additional replication origin, oriV, which permits copy-number to be induced transiently when propagated in a suitable host strain, e.g. EPI300, that supplies the cognate trans-replication protein TrfA. Previously, we used EL350 and EPI300 sequentially to recombineer oriV-equipped fosmid genomic clones and, subsequently, to induce copy-number of the resulting recombinant clone. To eliminate these intervening DNA isolation and transformation steps we retrofitted EL350 with a P_BAD-driven trfA gene generating strain MW005 that supports, independently, both recombineering and copy-number induction.

Results

The P_BAD-driven copy of cre in EL350 was replaced seamlessly with a copy of trfA, PCR-amplified from EPI300 chromosomal DNA, to generate MW005. This new strain has been used to both generate, via recombineering, a number of reporter gene fusions directly from pCC1FOS-based Caenorhabditis elegans genomic clones and to transiently induce copy-number of fosmid and BAC clones prior to DNA preparation.

Conclusions

By retrofitting EL350, an established 'recombineering' E. coli strain, with a tightly regulated copy of trfA we have produced a new strain, MW005, which combines recombineering capacity with the useful ability to transiently induce copy-number of oriV-equipped clones. By coupling these two steps in a single strain, use of MW005 will enable the more rapid recombineering-mediated production of recombinant clones in the yield and quality necessary for many downstream purposes.

The Caenorhabditis elegans sirtuin gene, sir-2.1, is widely expressed and induced upon caloric restriction

INTRODUCTION

As in yeast, flies and mammals, over-expression of the Caenorhabditis elegans sirtuin gene sir-2.1 leads to extension of lifespan and deletion of the gene shortens lifespan. The sir-2.1 gene, however, is located in an operon, an organization not taken into account in previous studies of this gene's expression.

MATERIALS AND METHODS

Recombineering allowed insertion of both a mCherry and a gfp reporter gene precisely at the end of the two protein-coding regions of the 4.5kb sir-2.1 operon within a 29.3kb genomic DNA fosmid clone.

RESULTS AND DISCUSSION

In C. elegans transgenic for this recombineered fosmid, with abundant food, the sir-2.1::mCherry distribution indicated that sir-2.1 is indeed expressed in the hypodermis and many nerve cells, as previously described, but also in the intestine and in muscles. This broader expression of sir-2.1, which would fit with an expectation that SIR2.1 function in influencing lifespan might be required in most cell types, arises from transcription starting with the gene upstream of sir-2.1 in the operon. Importantly, the expression of both genes in the operon increases upon starvation, this induction also depending on the operon promoter. Furthermore, SIR-2.1::mCherry undergoes a dynamic subcellular relocalization through starvation.

Global prediction of tissue-specific gene expression and context-dependent gene networks in Caenorhabditis elegans

Tissue-specific gene expression plays a fundamental role in metazoan biology and is an important aspect of many complex diseases. Nevertheless, an organism-wide map of tissue-specific expression remains elusive due to difficulty in obtaining these data experimentally. Here, we leveraged existing whole-animal Caenorhabditis elegans microarray data representing diverse conditions and developmental stages to generate accurate predictions of tissue-specific gene expression and experimentally validated these predictions. These patterns of tissue-specific expression are more accurate than existing high-throughput experimental studies for nearly all tissues; they also complement existing experiments by addressing tissue-specific expression present at particular developmental stages and in small tissues. We used these predictions to address several experimentally challenging questions, including the identification of tissue-specific transcriptional motifs and the discovery of potential miRNA regulation specific to particular tissues. We also investigate the role of tissue context in gene function through tissue-specific functional interaction networks. To our knowledge, this is the first study producing high-accuracy predictions of tissue-specific expression and interactions for a metazoan organism based on whole-animal data.

In animals, a crucial facet of any gene's function is the tissue or cell type in which that gene is expressed and the proteins that it interacts with in that cell. However, genome-wide identification of expression across the multitude of tissues of varying size and complexity is difficult to achieve experimentally. In this paper, we show that microararray data collected from whole animals can be analyzed to yield high-quality predictions of tissue-specific expression. These predictions are of better or comparable accuracy to tissue-specific expression determined from high-throughput experiments. Our results provide a global view of tissue-specific expression in Caenorhabditis elegans, allowing us to address the question of how expression patterns are regulated and to analyze how the functions of genes that are expressed in several tissues are influenced by the cellular context.

A simplified, robust, and streamlined procedure for the production of C. elegans transgenes via recombineering

Background

The nematode Caenorhabditis elegans has emerged as a powerful system to study biologic questions ranging from development to aging. The generation of transgenic animals is an important experimental tool and allows use of GFP fusion proteins to study the expression of genes of interest or generation of epitope tagged versions of specific genes. Transgenes are often generated by placing a promoter upstream of a reporter gene or cDNA. This often produces a representative expression pattern, but important exceptions have been observed. To better capture the genuine expression pattern and timing, several investigators have modified large pieces of DNA carried by BACs or fosmids for use in the construction of transgenic animals via recombineering. However, these techniques are not in widespread use despite the advantages when compared to traditional approaches. Additionally, some groups have encountered problems with employing these techniques. Hence, we sought identify ways to improve the simplicity and reliability of the procedure.

Results

We describe here several important modifications we have made to existing protocols to make the procedure simpler and more robust. Among these are the use of galK gene as a selection marker for both the positive and negative selection steps in recombineering, the use of R6K based plasmids which eliminate the need for extensive PCR product purification, a means to integrate the unc-119 marker on to the fosmid backbone, and placement of homology arms to commonly used GFP and TAP fusion genes flanking the galK cassette which reduces the cost of oligos by 50%.

Conclusion

We have made several significant changes that allow the production of C. elegans transgenes from a commercially available fosmid library in a robust and streamlined manner. These changes make the technique more attractive especially to small academic labs unfamiliar with recombineering.

ELT-2 is the predominant transcription factor controlling differentiation and function of the C. elegans intestine, from embryo to adult

Starting with SAGE-libraries prepared from C. elegans FAC-sorted embryonic intestine cells (8E-16E cell stage), from total embryos and from purified oocytes, and taking advantage of the NextDB in situ hybridization data base, we define sets of genes highly expressed from the zygotic genome, and expressed either exclusively or preferentially in the embryonic intestine or in the intestine of newly hatched larvae; we had previously defined a similarly expressed set of genes from the adult intestine. We show that an extended TGATAA-like sequence is essentially the only candidate for a cis-acting regulatory motif common to intestine genes expressed at all stages. This sequence is a strong ELT-2 binding site and matches the sequence of GATA-like sites found to be important for the expression of every intestinal gene so far analyzed experimentally. We show that the majority of these three sets of highly expressed intestinal-specific/intestinal-enriched genes respond strongly to ectopic expression of ELT-2 within the embryo. By flow-sorting elt-2(null) larvae from elt-2(+) larvae and then preparing Solexa/Illumina-SAGE libraries, we show that the majority of these genes also respond strongly to loss-of-function of ELT-2. To test the consequences of loss of other transcription factors identified in the embryonic intestine, we develop a strain of worms that is RNAi-sensitive only in the intestine; however, we are unable (with one possible exception) to identify any other transcription factor whose intestinal loss-of-function causes a phenotype of comparable severity to the phenotype caused by loss of ELT-2. Overall, our results support a model in which ELT-2 is the predominant transcription factor in the post-specification C. elegans intestine and participates directly in the transcriptional regulation of the majority (>80%) of intestinal genes. We present evidence that ELT-2 plays a central role in most aspects of C. elegans intestinal physiology: establishing the structure of the enterocyte, regulating enzymes and transporters involved in digestion and nutrition, responding to environmental toxins and pathogenic infections, and regulating the downstream intestinal components of the daf-2/daf-16 pathway influencing aging and longevity.