Background

Genome-wide expression data of gene microarrays can be used to infer gene networks. At a cellular level, a gene network provides a picture of the modules in which genes are densely connected, and of the hub genes, which are highly connected with other genes. A gene network is useful to identify the genes involved in the same pathway, in a protein complex or that are co-regulated. In this study, we used different methods to find gene networks in the ciliate Tetrahymena thermophila, and describe some important properties of this network, such as modules and hubs.

Methodology/Principal Findings

Using 67 single channel microarrays, we constructed the Tetrahymena gene network (TGN) using three methods: the Pearson correlation coefficient (PCC), the Spearman correlation coefficient (SCC) and the context likelihood of relatedness (CLR) algorithm. The accuracy and coverage of the three networks were evaluated using four conserved protein complexes in yeast. The CLR network with a Z-score threshold 3.49 was determined to be the most robust. The TGN was partitioned, and 55 modules were found. In addition, analysis of the arbitrarily determined 1200 hubs showed that these hubs could be sorted into six groups according to their expression profiles. We also investigated human disease orthologs in Tetrahymena that are missing in yeast and provide evidence indicating that some of these are involved in the same process in Tetrahymena as in human.

Conclusions/Significance

This study constructed a Tetrahymena gene network, provided new insights to the properties of this biological network, and presents an important resource to study Tetrahymena genes at the pathway level.

Histone H3 lysine 4 (K4) methylation is a prevalent mark associated with transcription activation and is mainly catalyzed by the MLL/SET1 family histone methyltransferases. A common feature of the mammalian MLL/SET1 complexes is the presence of three core components (RbBP5, Ash2L and WDR5) and a catalytic subunit containing a SET domain. Unlike most other histone lysine methyltransferases, all four proteins are required for efficient H3 K4 methylation. Despite extensive efforts, mechanisms for how three core components regulate MLL/SET1 methyltransferase activity remain elusive. Here we show that a heterodimer of Ash2L and RbBP5 has intrinsic histone methyltransferase activity. This activity requires the highly conserved SPRY domain of Ash2L and a short peptide of RbBP5. We demonstrate that both Ash2L and the MLL1 SET domain are capable of binding to S-adenosyl-L- [methyl-3H] methionine in the MLL1 core complex. Mutations in the MLL1 SET domain that fail to support overall H3 K4 methylation also compromise SAM binding by Ash2L. Taken together, our results show that the Ash2L/RbBP5 heterodimer plays a critical role in the overall catalysis of MLL1 mediated H3 K4 methylation. The results we describe here provide mechanistic insights for unique regulation of the MLL1 methyltransferase activity. It suggests that both Ash2L/RbBP5 and the MLL1 SET domain make direct contacts with the substrates and contribute to the formation of a joint catalytic center. Given the shared core configuration among all MLL/SET1 family HMTs, it will be interesting to test whether the mechanism we describe here can be generalized to other MLL/SET1 family members in the future.

Background

The model eukaryote, Tetrahymena thermophila, is the first ciliated protozoan whose genome has been sequenced, enabling genome-wide analysis of gene expression.

Methodology/Principal Findings

A genome-wide microarray platform containing the predicted coding sequences (putative genes) for T. thermophila is described, validated and used to study gene expression during the three major stages of the organism's life cycle: growth, starvation and conjugation.

Conclusions/Significance

Of the ∼27,000 predicted open reading frames, transcripts homologous to only ∼5900 are not detectable in any of these life cycle stages, indicating that this single-celled organism does indeed contain a large number of functional genes. Transcripts from over 5000 predicted genes are expressed at levels >5× corrected background and 95 genes are expressed at >250× corrected background in all stages. Transcripts homologous to 91 predicted genes are specifically expressed and 155 more are highly up-regulated in growing cells, while 90 are specifically expressed and 616 are up-regulated during starvation. Strikingly, transcripts homologous to 1068 predicted genes are specifically expressed and 1753 are significantly up-regulated during conjugation. The patterns of gene expression during conjugation correlate well with the developmental stages of meiosis, nuclear differentiation and DNA elimination. The relationship between gene expression and chromosome fragmentation is analyzed. Genes encoding proteins known to interact or to function in complexes show similar expression patterns, indicating that co-ordinate expression with putative genes of known function can identify genes with related functions. New candidate genes associated with the RNAi-like process of DNA elimination and with meiosis are identified and the late stages of conjugation are shown to be characterized by specific expression of an unexpectedly large and diverse number of genes not involved in nuclear functions.

Background

Tetrahymena thermophila, a widely studied model for cellular and molecular biology, is a binucleated single-celled organism with a germline micronucleus (MIC) and somatic macronucleus (MAC). The recent draft MAC genome assembly revealed low sequence repetitiveness, a result of the epigenetic removal of invasive DNA elements found only in the MIC genome. Such low repetitiveness makes complete closure of the MAC genome a feasible goal, which to achieve would require standard closure methods as well as removal of minor MIC contamination of the MAC genome assembly. Highly accurate preliminary annotation of Tetrahymena's coding potential was hindered by the lack of both comparative genomic sequence information from close relatives and significant amounts of cDNA evidence, thus limiting the value of the genomic information and also leaving unanswered certain questions, such as the frequency of alternative splicing.

Results

We addressed the problem of MIC contamination using comparative genomic hybridization with purified MIC and MAC DNA probes against a whole genome oligonucleotide microarray, allowing the identification of 763 genome scaffolds likely to contain MIC-limited DNA sequences. We also employed standard genome closure methods to essentially finish over 60% of the MAC genome. For the improvement of annotation, we have sequenced and analyzed over 60,000 verified EST reads from a variety of cellular growth and development conditions. Using this EST evidence, a combination of automated and manual reannotation efforts led to updates that affect 16% of the current protein-coding gene models. By comparing EST abundance, many genes showing apparent differential expression between these conditions were identified. Rare instances of alternative splicing and uses of the non-standard amino acid selenocysteine were also identified.

Conclusion

We report here significant progress in genome closure and reannotation of Tetrahymena thermophila. Our experience to date suggests that complete closure of the MAC genome is attainable. Using the new EST evidence, automated and manual curation has resulted in substantial improvements to the over 24,000 gene models, which will be valuable to researchers studying this model organism as well as for comparative genomics purposes.

In Tetrahymena, HHT1 and HHT2 genes encode the same major histone H3; HHT3 and HHT4 encode similar minor H3 variants (H3s), H3.3 and H3.4. Green fluorescent protein (GFP)-tagged H3 is deposited onto chromatin through a DNA replication-coupled (RC) pathway. GFP-tagged H3.3 and H3.4 can be deposited both by a transcription-associated, replication-independent (RI) pathway and also weakly by an RC pathway. Although both types of H3s can be deposited by the RC pathway, DNA repair synthesis associated with meiotic recombination utilizes H3 specifically. The regions distinguishing H3 and H3.3 for their deposition pathways were identified. RC major H3 is not essential. Cells can grow without major H3 if the minor H3s are expressed at high levels. Surprisingly, cells lacking RI H3s are also viable and maintain normal nucleosome density at a highly transcribed region. The RC H3 is not detectably deposited by the RI pathway, even when there are no RI H3s available, indicating that transcription-associated RI H3 deposition is not essential for transcription. Minor H3s are also required to produce viable sexual progeny and play an unexpected role in the germ line micronuclei late in conjugation that is unrelated to transcription.

In Tetrahymena thermophila, highly phosphorylated histone H1 of growing cells becomes partially dephosphorylated when cells are starved in preparation for conjugation. To determine the effects of H1 phosphorylation on gene expression, PCR-based subtractive hybridization was used to clone cDNAs that were differentially expressed during starvation in two otherwise-isogenic strains differing only in their H1s. H1 in A5 mutant cells lacked phosphorylation, and H1 in E5 cells mimicked constitutive H1 phosphorylation. Sequences enriched in A5 cells included genes encoding proteases. Sequences enriched in E5 cells included genes encoding cdc2 kinase and a Ser/Thr kinase. These results indicate that H1 phosphorylation plays an important role in regulating the pattern of gene expression during the starvation response and that its role in transcription regulation can be either positive or negative. Treatment of starved cells with a phosphatase inhibitor caused CDC2 gene overexpression. Expression of the E5 version of H1 in starved cells containing endogenous, wild-type H1 caused the wild-type H1 to remain highly phosphorylated. These results argue that Cdc2p is the kinase that phosphorylates Tetrahymena H1, establish a positive feedback mechanism between H1 phosphorylation and CDC2 expression, and indicate that CDC2 gene expression is regulated by an H1 phosphatase.

In the ciliate Tetrahymena thermophila, approximately 15% of the germ line micronuclear DNA sequences are eliminated during formation of the somatic macronucleus. The vast majority of the internal eliminated sequences (IESs) are repeated in the micronuclear genome, and several of them resemble transposable elements. Thus, it has been suggested that DNA elimination evolved as a means for removing invading DNAs. In the present study, bacterial neo genes introduced into the germ line micronuclei were eliminated from the somatic genome. The efficiency of elimination from two different loci increased dramatically with the copy number of the neo genes in the micronuclei. The timing of neo elimination is similar to that of endogenous IESs, and they both produce bidirectional transcripts of the eliminated element, suggesting that the deletion of neo occurred by the same mechanism as elimination of endogenous IESs. These results indicate that repetition of an element in the micronucleus enhances the efficiency of its elimination from the newly formed somatic genome of Tetrahymena thermophila. The implications of these data in relation to the function and mechanism of IES elimination are discussed.

A mathematical model suggests that a HEPA/vaccine approach is viable for most buildings after a large-scale anthrax attack.

We developed a mathematical model to compare 2 indoor remediation strategies in the aftermath of an outdoor release of 1.5 kg of anthrax spores in lower Manhattan. The 2 strategies are the fumigation approach used after the 2001 postal anthrax attack and a HEPA/vaccine plan, which relies on HEPA vacuuming, HEPA air cleaners, and vaccination of reoccupants. The HEPA/vaccine approach leads to few anthrax cases among reoccupants if applied to all but the most heavily contaminated buildings, and recovery is much faster than under the decades-long fumigation plan. Only modest environmental sampling is needed. A surge capacity of 10,000 to 20,000 Hazmat workers is required to perform remediation within 6 to 12 months and to avoid permanent mass relocation. Because of the possibility of a campaign of terrorist attacks, serious consideration should be given to allowing or encouraging voluntary self-service cleaning of lightly contaminated rooms by age-appropriate, vaccinated, partially protected (through masks or hoods) reoccupants or owners.

In Tetrahymena cells, phosphorylation of linker histone H1 regulates transcription of specific genes. Phosphorylation acts by creating a localized negative charge patch and phenocopies the loss of H1 from chromatin, suggesting that it affects transcription by regulating the dissociation of H1 from chromatin. To test this hypothesis, we used FRAP of GFP-tagged H1 to analyze the effects of mutations that either eliminate or mimic phosphorylation on the binding of H1 to chromatin both in vivo and in vitro. We demonstrate that phosphorylation can increase the rate of dissociation of H1 from chromatin, providing a mechanism by which it can affect H1 function in vivo. We also demonstrate a previously undescribed ATP-dependent process that has a global effect on the dynamic binding of linker histone to chromatin.