Conservation, regulation, synteny, and introns in a large-scale C. briggsae-C. elegans genomic alignment

Making worm guts: the gene regulatory network of the Caenorhabditis elegans endoderm

The nematode Caenorhabditis elegans is a triploblastic ecdysozoan, which, although it contains too few cells during embryogenesis to create discernible germ "layers," deploys similar programs for germ layer differentiation used in animals with many more cells. The endoderm arises from a single progenitor, the E cell, and is selected from among three possible fates by a three-state combinatorial regulatory system involving intersecting cell-intrinsic and intercellular signals. The core gene regulatory cascade that drives endoderm development, extending from early maternal regulators to terminal differentiation genes, is characterized by activation of successive tiers of transcription factors, including a sequential cascade of redundant GATA transcription factors. Each tier is punctuated by a cell division, raising the possibility that intercession of one cell cycle round, or DNA replication, is required for activation of the next tier. The existence of each tier in the regulatory hierarchy is justified by the assignment of a unique task and each invariably performs at least two functions: to activate the regulators in the next tier and to perform one other activity distinct from that of the next tier. While the regulatory inputs that initiate endoderm development are highly divergent, they mobilize a gene regulatory network for endoderm development that appears to be common to all triploblastic metazoans. Genome-wide functional genomic approaches, including identification of >800 transcripts that exhibit the same regulatory patterns as a number of endoderm-specific genes, are contributing to elucidation of the complete endoderm gene regulatory network in C. elegans. Dissection of the architecture of the C. elegans endoderm network may provide insights into the evolutionary plasticity and origins of this germ layer.

Intron evolution as a population-genetic process

Debate over the mechanisms responsible for the phylogenetic and genomic distribution of introns has proceeded largely without consideration of the population-genetic forces influencing the establishment and retention of novel genetic elements. However, a simple model incorporating random genetic drift and weak mutation pressure against intron-containing alleles yields predictions consistent with a diversity of observations: (i) the rarity of introns in unicellular organisms with large population sizes, and their expansion after the origin of multicellular organisms with reduced population sizes; (ii) the relationship between intron abundance and the stringency of splice-site requirements; (iii) the tendency for introns to be more numerous and longer in regions of low recombination; and (iv) the bias toward phase-0 introns. This study provides a second example of a mechanism whereby genomic complexity originates passively as a "pathological" response to small population size, and raises difficulties for the idea that ancient introns played a major role in the origin of genes by exon shuffling.

Comparative genomics approaches to study organism similarities and differences

Comparative genomics is a large-scale, holistic approach that compares two or more genomes to discover the similarities and differences between the genomes and to study the biology of the individual genomes. Comparative studies can be performed at different levels of the genomes to obtain multiple perspectives about the organisms. We discuss in detail the type of analyses that offer significant biological insights in the comparisons of (1) genome structure including overall genome statistics, repeats, genome rearrangement at both DNA and gene level, synteny, and breakpoints; (2) coding regions including gene content, protein content, orthologs, and paralogs; and (3) noncoding regions including the prediction of regulatory elements. We also briefly review the currently available computational tools in comparative genomics such as algorithms for genome-scale sequence alignment, gene identification, and nonhomology-based function prediction.

Analysis of similarity within 142 pairs of orthologous intergenic regions of Caenorhabditis elegans and Caenorhabditis briggsae

Patterns of similarity between genomes of related species reflect the distribution of selective constraint within DNA. We analyzed alignments of 142 orthologous intergenic regions of Caenorhabditis elegans and Caenorhabditis briggsae and found a mosaic pattern with regions of high similarity (phylogenetic footprints) interspersed with non-alignable sequences. Footprints cover approximately 20% of intergenic regions, often occur in clumps and are rare within 5' UTRs but common within 3' UTRs. The footprints have a higher ratio of transitions to transversions than expected at random and a higher GC content than the rest of the intergenic region. The number of footprints and the GC content of footprints within an intergenic region are higher when genes are oriented so that their 5' ends form the boundaries of the intergenic region. Overall, the patterns and characteristics identified here, along with other comparative and experimental studies, suggest that many footprints have a regulatory function, although other types of function are also possible. These conclusions may be quite general across eukaryotes, and the characteristics of conserved regulatory elements determined from genomic comparisons can be useful in prediction of regulation sites within individual DNA sequences.

Multiple subunits of the Caenorhabditis elegans anaphase-promoting complex are required for chromosome segregation during meiosis I

Two genes, originally identified in genetic screens for Caenorhabditis elegans mutants that arrest in metaphase of meiosis I, prove to encode subunits of the anaphase-promoting complex or cyclosome (APC/C). RNA interference studies reveal that these and other APC/C subunits are essential for the segregation of chromosomal homologs during meiosis I. Further, chromosome segregation during meiosis I requires APC/C functions in addition to the release of sister chromatid cohesion.

A genomic view of alternative splicing

Recent genome-wide analyses of alternative splicing indicate that 40-60% of human genes have alternative splice forms, suggesting that alternative splicing is one of the most significant components of the functional complexity of the human genome. Here we review these recent results from bioinformatics studies, assess their reliability and consider the impact of alternative splicing on biological functions. Although the 'big picture' of alternative splicing that is emerging from genomics is exciting, there are many challenges. High-throughput experimental verification of alternative splice forms, functional characterization, and regulation of alternative splicing are key directions for research. We recommend a community-based effort to discover and characterize alternative splice forms comprehensively throughout the human genome.

The K(A)/K(S) ratio test for assessing the protein-coding potential of genomic regions: an empirical and simulation study

Comparative genomics is a simple, powerful way to increase the accuracy of gene prediction. In this study, we show the utility of a simple test for the identification of protein-coding exons using human/mouse sequence comparisons. The test takes advantage of the fact that in the vast majority of coding regions, synonymous substitutions (K(S)) occur much more frequently than nonsynonymous ones (K(A)) and uses the K(A)/K(S) ratio as the criterion. We show the following: (1) most of the human and mouse exons are sufficiently long and have a suitable degree of sequence divergence for the test to perform reliably; (2) the test is suited for the identification of long exons and single exon genes, which are difficult to predict by current methods; (3) the test has a false-negative rate, lower than most of current gene prediction methods and a false-positive rate lower than all current methods; (4) the test has been automated and can be used in combination with other existing gene-prediction methods.

Conservation of long-range synteny and microsynteny between the genomes of two distantly related nematodes

BACKGROUND

Comparisons between the genomes of the closely related nematodes Caenorhabditis elegans and Caenorhabditis briggsae reveal high rates of rearrangement, with a bias towards within-chromosome events. To assess whether this pattern is true of nematodes in general, we have used genome sequence to compare two nematode species that last shared a common ancestor approximately 300 million years ago: the model C. elegans and the filarial parasite Brugia malayi.

RESULTS

An 83 kb region flanking the gene for Bm-mif-1 (macrophage migration inhibitory factor, a B. malayi homolog of a human cytokine) was sequenced. When compared to the complete genome of C. elegans, evidence for conservation of long-range synteny and microsynteny was found. Potential C. elegans orthologs for II of the 12 protein-coding genes predicted in the B. malayi sequence were identified. Ten of these orthologs were located on chromosome I, with eight clustered in a 2.3 Mb region. While several, relatively local, intrachromosomal rearrangements have occurred, the order, composition, and configuration of two gene clusters, each containing three genes, was conserved. Comparison of B. malayi BAC-end genome survey sequence to C. elegans also revealed a bias towards intrachromosome rearrangements.

CONCLUSIONS

We suggest that intrachromosomal rearrangement is a major force driving chromosomal organization in nematodes, but is constrained by the interdigitation of functional elements of neighboring genes.

The probability of preservation of a newly arisen gene duplicate

Newly emerging data from genome sequencing projects suggest that gene duplication, often accompanied by genetic map changes, is a common and ongoing feature of all genomes. This raises the possibility that differential expansion/contraction of various genomic sequences may be just as important a mechanism of phenotypic evolution as changes at the nucleotide level. However, the population-genetic mechanisms responsible for the success vs. failure of newly arisen gene duplicates are poorly understood. We examine the influence of various aspects of gene structure, mutation rates, degree of linkage, and population size (N) on the joint fate of a newly arisen duplicate gene and its ancestral locus. Unless there is active selection against duplicate genes, the probability of permanent establishment of such genes is usually no less than 1/(4N) (half of the neutral expectation), and it can be orders of magnitude greater if neofunctionalizing mutations are common. The probability of a map change (reassignment of a key function of an ancestral locus to a new chromosomal location) induced by a newly arisen duplicate is also generally >1/(4N) for unlinked duplicates, suggesting that recurrent gene duplication and alternative silencing may be a common mechanism for generating microchromosomal rearrangements responsible for postreproductive isolating barriers among species. Relative to subfunctionalization, neofunctionalization is expected to become a progressively more important mechanism of duplicate-gene preservation in populations with increasing size. However, even in large populations, the probability of neofunctionalization scales only with the square of the selective advantage. Tight linkage also influences the probability of duplicate-gene preservation, increasing the probability of subfunctionalization but decreasing the probability of neofunctionalization.

An extensive class of small RNAs in Caenorhabditis elegans

The lin-4 and let-7 antisense RNAs are temporal regulators that control the timing of developmental events in Caenorhabditis elegans by inhibiting translation of target mRNAs. let-7 RNA is conserved among bilaterian animals, suggesting that this class of small RNAs [microRNAs (miRNAs)] is evolutionarily ancient. Using bioinformatics and cDNA cloning, we found 15 new miRNA genes in C. elegans. Several of these genes express small transcripts that vary in abundance during C. elegans larval development, and three of them have apparent homologs in mammals and/or insects. Small noncoding RNAs of the miRNA class appear to be numerous and diverse.

Alternative splicing of the Drosophila Dscam pre-mRNA is both temporally and spatially regulated

The Drosophila melanogaster Down syndrome cell adhesion molecule (Dscam) gene encodes an axon guidance receptor that can express 38,016 different mRNAs by virtue of alternative splicing. The Dscam gene contains 95 alternative exons that are organized into four clusters of 12, 48, 33, and 2 exons each. Although numerous Dscam mRNA isoforms can be synthesized, it remains to be determined whether different Dscam isoforms are synthesized at different times in development or in different tissues. We have investigated the alternative splicing of the Dscam exon 4 cluster, which contains 12 mutually exclusive alternative exons, and found that Dscam exon 4 alternative splicing is developmentally regulated. The most highly regulated exon, 4.2, is infrequently used in early embryos but is the predominant exon 4 variant used in adults. Moreover, the developmental regulation of exon 4.2 alternative splicing is conserved in D. yakuba. In addition, different adult tissues express distinct collections of Dscam mRNA isoforms. Given the role of Dscam in neural development, these results suggest that the regulation of alternative splicing plays an important role in determining the specificity of neuronal wiring. In addition, this work provides a framework to determine the mechanisms by which complex alternative splicing events are regulated.

Surveying Saccharomyces genomes to identify functional elements by comparative DNA sequence analysis

Comparative sequence analysis has facilitated the discovery of protein coding genes and important functional sequences within proteins, but has been less useful for identifying functional sequence elements in nonprotein-coding DNA because the relatively rapid rate of change of nonprotein-coding sequences and the relative simplicity of non-coding regulatory sequence elements necessitates the comparison of sequences of relatively closely related species. We tested the use of comparative DNA sequence analysis to aid identification of promoter regulatory elements, nonprotein-coding RNA genes, and small protein-coding genes by surveying random DNA sequences of several Saccharomyces yeast species, with the goal of learning which species are best suited for comparisons with S. cerevisiae. We also determined the DNA sequence of a few specific promoters and RNA genes of several Saccharomyces species to determine the degree of conservation of known functional elements within the genome. Our results lead us to conclude that comparative DNA sequence analysis will enable identification of functionally conserved elements within the yeast genome, and suggest a path for obtaining this information.

Introns of Entamoeba histolytica and Entamoeba dispar

The genome of Entamoeba histolytica is considered to possess very few intervening sequences (introns), as only 5 intron-containing genes from this protozoan parasite have been reported so far. However, while sequencing a number of genomic contigs as well as three independent genes coding for ribosomal protein L27a, we have identified 9 additional intron-containing genes of E. histolytica and the closely related species Entamoeba dispar, indicating that introns are more common in these organisms than previously suggested. The various amoeba introns are relatively short comprising between 46 and 115 nucleotides only and have a higher AT-content compared to the corresponding exon sequences. In contrast to higher eukaryotes, amoeba introns do not contain a well-conserved branch point consensus, and have extended donor and acceptor splice sites of the sequences G

Comparative analysis of the gene-dense ACHE/TFR2 region on human chromosome 7q22 with the orthologous region on mouse chromosome 5

Chromosome 7q22 has been the focus of many cytogenetic and molecular studies aimed at delineating regions commonly deleted in myeloid leukemias and myelodysplastic syndromes. We have compared a gene-dense, GC-rich sub-region of 7q22 with the orthologous region on mouse chromosome 5. A physical map of 640 kb of genomic DNA from mouse chromosome 5 was derived from a series of overlapping bacterial artificial chromosomes. A 296 kb segment from the physical map, spanning ACHE: to Tfr2, was compared with 267 kb of human sequence. We identified a conserved linkage of 12 genes including an open reading frame flanked by ACHE: and Asr2, a novel cation-chloride cotransporter interacting protein Cip1, Ephb4, Zan and Perq1. While some of these genes have been previously described, in each case we present new data derived from our comparative sequence analysis. Adjacent unfinished sequence data from the mouse contains an orthologous block of 10 additional genes including three novel cDNA sequences that we subsequently mapped to human 7q22. Methods for displaying comparative genomic information, including unfinished sequence data, are becoming increasingly important. We supplement our printed comparative analysis with a new, Web-based program called Laj (local alignments with java). Laj provides interactive access to archived pairwise sequence alignments via the WWW. It displays synchronized views of a dot-plot, a percent identity plot, a nucleotide-level local alignment and a variety of relevant annotations. Our mouse-human comparison can be viewed at http://web.uvic.ca/~bioweb/laj.html. Laj is available at http://bio.cse.psu.edu/, along with online documentation and additional examples of annotated genomic regions.

Alternative splicing: increasing diversity in the proteomic world

How can the genome of Drosophila melanogaster contain fewer genes than the undoubtedly simpler organism Caenorhabditis elegans? The answer must lie within their proteomes. It is becoming clear that alternative splicing has an extremely important role in expanding protein diversity and might therefore partially underlie the apparent discrepancy between gene number and organismal complexity. Alternative splicing can generate more transcripts from a single gene than the number of genes in an entire genome. However, for the vast majority of alternative splicing events, the functional significance is unknown. Developing a full catalog of alternatively spliced transcripts and determining each of their functions will be a major challenge of the upcoming proteomic era.