Transcriptomic analysis of the Myxococcus xanthus FruA regulon, and comparative developmental transcriptomic analysis of two fruiting body forming species, Myxococcus xanthus and Myxococcus stipitatus.

BACKGROUND: The Myxococcales are well known for their predatory and developmental social processes, and for the molecular complexity of regulation of these processes. Many species within this order have unusually large genomes compared to other bacteria, and their genomes have many genes that are unique to one specific sequenced species or strain. Here, we describe RNAseq based transcriptome analysis of the FruA regulon of Myxococcus xanthus and a comparative RNAseq analysis of two Myxococcus species, M. xanthus and Myxococcus stipitatus, as they respond to starvation and begin forming fruiting bodies. RESULTS: We show that both species have large numbers of genes that are developmentally regulated, with over half the genome showing statistically significant changes in expression during development in each species. We also included a non-fruiting mutant of M. xanthus that is missing the transcriptional regulator FruA to identify the direct and indirect FruA regulon and to identify transcriptional changes that are specific to fruiting and not just the starvation response. We then identified Interpro gene ontologies and COG annotations that are significantly up- or down-regulated during development in each species. Our analyses support previous data for M. xanthus showing developmental upregulation of signal transduction genes, and downregulation of genes related to cell-cycle, translation, metabolism, and in some cases, DNA replication. Gene expression in M. stipitatus follows similar trends. Although not all specific genes show similar regulation patterns in both species, many critical developmental genes in M. xanthus have conserved expression patterns in M. stipitatus, and some groups of otherwise unstudied orthologous genes share expression patterns. CONCLUSIONS: By identifying the FruA regulon and identifying genes that are similarly and uniquely regulated in two different species, this work provides a more complete picture of transcription during Myxococcus development. We also provide an R script to allow other scientists to mine our data for genes whose expression patterns match a user-selected gene of interest.

The time-resolved transcriptome of C. elegans.

We generated detailed RNA-seq data for the nematode Caenorhabditis elegans with high temporal resolution in the embryo as well as representative samples from post-embryonic stages across the life cycle. The data reveal that early and late embryogenesis is accompanied by large numbers of genes changing expression, whereas fewer genes are changing in mid-embryogenesis. This lull in genes changing expression correlates with a period during which histone mRNAs produce almost 40% of the RNA-seq reads. We find evidence for many more splice junctions than are annotated in WormBase, with many of these suggesting alternative splice forms, often with differential usage over the life cycle. We annotated internal promoter usage in operons using SL1 and SL2 data. We also uncovered correlated transcriptional programs that span >80 kb. These data provide detailed annotation of the C. elegans transcriptome.

Comparative analysis of the transcriptome across distant species.

The transcriptome is the readout of the genome. Identifying common features in it across distant species can reveal fundamental principles. To this end, the ENCODE and modENCODE consortia have generated large amounts of matched RNA-sequencing data for human, worm and fly. Uniform processing and comprehensive annotation of these data allow comparison across metazoan phyla, extending beyond earlier within-phylum transcriptome comparisons and revealing ancient, conserved features. Specifically, we discover co-expression modules shared across animals, many of which are enriched in developmental genes. Moreover, we use expression patterns to align the stages in worm and fly development and find a novel pairing between worm embryo and fly pupae, in addition to the embryo-to-embryo and larvae-to-larvae pairings. Furthermore, we find that the extent of non-canonical, non-coding transcription is similar in each organism, per base pair. Finally, we find in all three organisms that the gene-expression levels, both coding and non-coding, can be quantitatively predicted from chromatin features at the promoter using a 'universal model' based on a single set of organism-independent parameters.

Specific roles for the GATA transcription factors end-1 and end-3 during C. elegans E-lineage development.

end-1 and end-3 are GATA transcription factors important for specifying endoderm cell fate in Caenorhabditis elegans. Deletion of both factors together results in larval arrest, 0% survival and a fate change in the endoderm-specifying E lineage. Individual deletions of either factor, however, result in the development of viable, fertile adults, with 100% of worms developing to adults for end-1(-) and 95% for end-3(-). We sought to quantify the variable phenotypes seen in both deletions using automated cell lineaging. We quantified defects in cell lifetime, cell movement and division axis in end-3(-) embryos, while quantifying perturbations in downstream reporter gene expression in strains with homozygous deletions for either gene, showing that each deletion leads to a unique profile of downstream perturbations in gene expression and cellular phenotypes with a high correlation between early and late defects. Combining observations in both cellular and gene expression defects we found that misaligned divisions at the E2 stage resulted in ectopic expression of the Notch target ref-1 in end-3(-) embryos. Using a maximum likelihood phylogenetic approach we found end-1 and end-3 split to form two distinct clades within the Caenorhabditis lineage with distinct DNA-binding structures. These results indicate that end-1 and end-3 have each evolved into genes with unique functions during endoderm development, that end-3(-) embryos have a delay in the onset of E lineage cell fate and that end-1 has only a partially penetrant ability to activate E lineage fate.

Automated analysis of embryonic gene expression with cellular resolution in C. elegans.

We describe a system that permits the automated analysis of reporter gene expression in Caenorhabditis elegans with cellular resolution continuously during embryogenesis. We demonstrate its utility by defining the expression patterns of reporters for several embryonically expressed transcription factors. The invariant cell lineage permits the automated alignment of multiple expression profiles, allowing direct comparison of the expression of different genes' reporters. We also used this system to monitor perturbations to normal development involving changes both in cell-division timing and in cell fate. Systematic application of this system could reveal the gene activity of each cell throughout development.

Genomic mapping of single-stranded DNA in hydroxyurea-challenged yeasts identifies origins of replication.

During DNA replication one or both strands transiently become single stranded: first at the sites where initiation of DNA synthesis occurs (known as origins of replication) and subsequently on the lagging strands of replication forks as discontinuous Okazaki fragments are generated. We report a genome-wide analysis of single-stranded DNA (ssDNA) formation in the presence of hydroxyurea during DNA replication in wild-type and checkpoint-deficient rad53 Saccharomyces cerevisiae cells. In wild-type cells, ssDNA was first observed at a subset of replication origins and later 'migrated' bi-directionally, suggesting that ssDNA formation is associated with continuously moving replication forks. In rad53 cells, ssDNA was observed at virtually every known origin, but remained there over time, suggesting that replication forks stall. Telomeric regions seemed to be particularly sensitive to the loss of Rad53 checkpoint function. Replication origins in Schizosaccharomyces pombe were also mapped using our method.