The genome sequence of Caenorhabditis briggsae: a platform for comparative genomics

The soil nematodes Caenorhabditis briggsae and Caenorhabditis elegans diverged from a common ancestor roughly 100 million years ago and yet are almost indistinguishable by eye. They have the same chromosome number and genome sizes, and they occupy the same ecological niche. To explore the basis for this striking conservation of structure and function, we have sequenced the C. briggsae genome to a high-quality draft stage and compared it to the finished C. elegans sequence. We predict approximately 19,500 protein-coding genes in the C. briggsae genome, roughly the same as in C. elegans. Of these, 12,200 have clear C. elegans orthologs, a further 6,500 have one or more clearly detectable C. elegans homologs, and approximately 800 C. briggsae genes have no detectable matches in C. elegans. Almost all of the noncoding RNAs (ncRNAs) known are shared between the two species. The two genomes exhibit extensive colinearity, and the rate of divergence appears to be higher in the chromosomal arms than in the centers. Operons, a distinctive feature of C. elegans, are highly conserved in C. briggsae, with the arrangement of genes being preserved in 96% of cases. The difference in size between the C. briggsae (estimated at approximately 104 Mbp) and C. elegans (100.3 Mbp) genomes is almost entirely due to repetitive sequence, which accounts for 22.4% of the C. briggsae genome in contrast to 16.5% of the C. elegans genome. Few, if any, repeat families are shared, suggesting that most were acquired after the two species diverged or are undergoing rapid evolution. Coclustering the C. elegans and C. briggsae proteins reveals 2,169 protein families of two or more members. Most of these are shared between the two species, but some appear to be expanding or contracting, and there seem to be as many as several hundred novel C. briggsae gene families. The C. briggsae draft sequence will greatly improve the annotation of the C. elegans genome. Based on similarity to C. briggsae, we found strong evidence for 1,300 new C. elegans genes. In addition, comparisons of the two genomes will help to understand the evolutionary forces that mold nematode genomes.     

Adaptive gene expression divergence inferred from population genomics

Detailed studies of individual genes have shown that gene expression divergence often results from adaptive evolution of regulatory sequence. Genome-wide analyses, however, have yet to unite patterns of gene expression with polymorphism and divergence to infer population genetic mechanisms underlying expression evolution. Here, we combined genomic expression data—analyzed in a phylogenetic context—with whole genome light-shotgun sequence data from six Drosophila simulans lines and reference sequences from D. melanogaster and D. yakuba. These data allowed us to use molecular population genetics to test for neutral versus adaptive gene expression divergence on a genomic scale. We identified recent and recurrent adaptive evolution along the D. simulans lineage by contrasting sequence polymorphism within D. simulans to divergence from D. melanogaster and D. yakuba. Genes that evolved higher levels of expression in D. simulans have experienced adaptive evolution of the associated 3′ flanking and amino acid sequence. Concomitantly, these genes are also decelerating in their rates of protein evolution, which is in agreement with the finding that highly expressed genes evolve slowly. Interestingly, adaptive evolution in 5′ cis-regulatory regions did not correspond strongly with expression evolution. Our results provide a genomic view of the intimate link between selection acting on a phenotype and associated genic evolution.     

Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans

The population genetic perspective is that the processes shaping genomic variation can be revealed only through simultaneous investigation of sequence polymorphism and divergence within and between closely related species. Here we present a population genetic analysis of Drosophila simulans based on whole-genome shotgun sequencing of multiple inbred lines and comparison of the resulting data to genome assemblies of the closely related species, D. melanogaster and D. yakuba. We discovered previously unknown, large-scale fluctuations of polymorphism and divergence along chromosome arms, and significantly less polymorphism and faster divergence on the X chromosome. We generated a comprehensive list of functional elements in the D. simulans genome influenced by adaptive evolution. Finally, we characterized genomic patterns of base composition for coding and noncoding sequence. These results suggest several new hypotheses regarding the genetic and biological mechanisms controlling polymorphism and divergence across the Drosophila genome, and provide a rich resource for the investigation of adaptive evolution and functional variation in D. simulans.     

The complete genome sequence and comparative genome analysis of the high pathogenicity Yersinia enterocolitica strain 8081

The human enteropathogen, Yersinia enterocolitica, is a significant link in the range of Yersinia pathologies extending from mild gastroenteritis to bubonic plague. Comparison at the genomic level is a key step in our understanding of the genetic basis for this pathogenicity spectrum. Here we report the genome of Y. enterocolitica strain 8081 (serotype 0:8; biotype 1B) and extensive microarray data relating to the genetic diversity of the Y. enterocolitica species. Our analysis reveals that the genome of Y. enterocolitica strain 8081 is a patchwork of horizontally acquired genetic loci, including a plasticity zone of 199 kb containing an extraordinarily high density of virulence genes. Microarray analysis has provided insights into species-specific Y. enterocolitica gene functions and the intraspecies differences between the high, low, and nonpathogenic Y. enterocolitica biotypes. Through comparative genome sequence analysis we provide new information on the evolution of the Yersinia. We identify numerous loci that represent ancestral clusters of genes potentially important in enteric survival and pathogenesis, which have been lost or are in the process of being lost, in the other sequenced Yersinia lineages. Our analysis also highlights large metabolic operons in Y. enterocolitica that are absent in the related enteropathogen, Yersinia pseudotuberculosis, indicating major differences in niche and nutrients used within the mammalian gut. These include clusters directing, the production of hydrogenases, tetrathionate respiration, cobalamin synthesis, and propanediol utilisation. Along with ancestral gene clusters, the genome of Y. enterocolitica has revealed species-specific and enteropathogen-specific loci. This has provided important insights into the pathology of this bacterium and, more broadly, into the evolution of the genus. Moreover, wider investigations looking at the patterns of gene loss and gain in the Yersinia have highlighted common themes in the genome evolution of other human enteropathogens.     

Comparison of C. elegans and C. briggsae genome sequences reveals extensive conservation of chromosome organization and synteny

To determine whether the distinctive features of Caenorhabditis elegans chromosomal organization are shared with the C. briggsae genome, we constructed a single nucleotide polymorphism–based genetic map to order and orient the whole genome shotgun assembly along the six C. briggsae chromosomes. Although these species are of the same genus, their most recent common ancestor existed 80–110 million years ago, and thus they are more evolutionarily distant than, for example, human and mouse. We found that, like C. elegans chromosomes, C. briggsae chromosomes exhibit high levels of recombination on the arms along with higher repeat density, a higher fraction of intronic sequence, and a lower fraction of exonic sequence compared with chromosome centers. Despite extensive intrachromosomal rearrangements, 1:1 orthologs tend to remain in the same region of the chromosome, and colinear blocks of orthologs tend to be longer in chromosome centers compared with arms. More strikingly, the two species show an almost complete conservation of synteny, with 1:1 orthologs present on a single chromosome in one species also found on a single chromosome in the other. The conservation of both chromosomal organization and synteny between these two distantly related species suggests roles for chromosome organization in the fitness of an organism that are only poorly understood presently.     

The parasexual cycle in Candida albicans provides an alternative pathway to meiosis for the formation of recombinant strains

Candida albicans has an elaborate, yet efficient, mating system that promotes conjugation between diploid a and α strains. The product of mating is a tetraploid a/α cell that must undergo a reductional division to return to the diploid state. Despite the presence of several “meiosis-specific” genes in the C. albicans genome, a meiotic program has not been observed. Instead, tetraploid products of mating can be induced to undergo efficient, random chromosome loss, often producing strains that are diploid, or close to diploid, in ploidy. Using SNP and comparative genome hybridization arrays we have now analyzed the genotypes of products from the C. albicans parasexual cycle. We show that the parasexual cycle generates progeny strains with shuffled combinations of the eight C. albicans chromosomes. In addition, several isolates had undergone extensive genetic recombination between homologous chromosomes, including multiple gene conversion events. Progeny strains exhibited altered colony morphologies on laboratory media, demonstrating that the parasexual cycle generates phenotypic variants of C. albicans. In several fungi, including Saccharomyces cerevisiae and Schizosaccharomyces pombe, the conserved Spo11 protein is integral to meiotic recombination, where it is required for the formation of DNA double-strand breaks. We show that deletion of SPO11 prevented genetic recombination between homologous chromosomes during the C. albicans parasexual cycle. These findings suggest that at least one meiosis-specific gene has been re-programmed to mediate genetic recombination during the alternative parasexual life cycle of C. albicans. We discuss, in light of the long association of C. albicans with warm-blooded animals, the potential advantages of a parasexual cycle over a conventional sexual cycle.     

Widespread discordance of gene trees with species tree in Drosophila: evidence for incomplete lineage sorting

The phylogenetic relationship of the now fully sequenced species Drosophila erecta and D. yakuba with respect to the D. melanogaster species complex has been a subject of controversy. All three possible groupings of the species have been reported in the past, though recent multi-gene studies suggest that D. erecta and D. yakuba are sister species. Using the whole genomes of each of these species as well as the four other fully sequenced species in the subgenus Sophophora, we set out to investigate the placement of D. erecta and D. yakuba in the D. melanogaster species group and to understand the cause of the past incongruence. Though we find that the phylogeny grouping D. erecta and D. yakuba together is the best supported, we also find widespread incongruence in nucleotide and amino acid substitutions, insertions and deletions, and gene trees. The time inferred to span the two key speciation events is short enough that under the coalescent model, the incongruence could be the result of incomplete lineage sorting. Consistent with the lineage-sorting hypothesis, substitutions supporting the same tree were spatially clustered. Support for the different trees was found to be linked to recombination such that adjacent genes support the same tree most often in regions of low recombination and substitutions supporting the same tree are most enriched roughly on the same scale as linkage disequilibrium, also consistent with lineage sorting. The incongruence was found to be statistically significant and robust to model and species choice. No systematic biases were found. We conclude that phylogenetic incongruence in the D. melanogaster species complex is the result, at least in part, of incomplete lineage sorting. Incomplete lineage sorting will likely cause phylogenetic incongruence in many comparative genomics datasets. Methods to infer the correct species tree, the history of every base in the genome, and comparative methods that control for and/or utilize this information will be valuable advancements for the field of comparative genomics.     

Helicobacter pylori evolution: lineage- specific adaptations in homologs of eukaryotic Sel1-like genes

Geographic partitioning is postulated to foster divergence of Helicobacter pylori populations as an adaptive response to local differences in predominant host physiology. H. pylori's ability to establish persistent infection despite host inflammatory responses likely involves active management of host defenses using bacterial proteins that may themselves be targets for adaptive evolution. Sequenced H. pylori genomes encode a family of eight or nine secreted proteins containing repeat motifs that are characteristic of the eukaryotic Sel1 regulatory protein, whereas the related Campylobacter and Wolinella genomes each contain only one or two such “Sel1-like repeat” (SLR) genes (“slr genes”). Signatures of positive selection (ratio of nonsynonymous to synonymous mutations, d_N/d_S = ω > 1) were evident in the evolutionary history of H. pylori slr gene family expansion. Sequence analysis of six of these slr genes (hp0160, hp0211, hp0235, hp0519, hp0628, and hp1117) from representative East Asian, European, and African H. pylori strains revealed that all but hp0628 had undergone positive selection, with different amino acids often selected in different regions. Most striking was a divergence of Japanese and Korean alleles of hp0519, with Japanese alleles having undergone particularly strong positive selection (ω_J > 25), whereas alleles of other genes from these populations were intermingled. Homology-based structural modeling localized most residues under positive selection to SLR protein surfaces. Rapid evolution of certain slr genes in specific H. pylori lineages suggests a model of adaptive change driven by selection for fine-tuning of host responses, and facilitated by geographic isolation. Characterization of such local adaptations should help elucidate how H. pylori manages persistent infection, and potentially lead to interventions tailored to diverse human populations.     

Genome sequence of Rickettsia bellii illuminates the role of amoebae in gene exchanges between intracellular pathogens

The recently sequenced Rickettsia felis genome revealed an unexpected plasmid carrying several genes usually associated with DNA transfer, suggesting that ancestral rickettsiae might have been endowed with a conjugation apparatus. Here we present the genome sequence of Rickettsia bellii, the earliest diverging species of known rickettsiae. The 1,552,076 base pair–long chromosome does not exhibit the colinearity observed between other rickettsia genomes, and encodes a complete set of putative conjugal DNA transfer genes most similar to homologues found in Protochlamydia amoebophila UWE25, an obligate symbiont of amoebae. The genome exhibits many other genes highly similar to homologues in intracellular bacteria of amoebae. We sought and observed sex pili-like cell surface appendages for R. bellii. We also found that R. bellii very efficiently multiplies in the nucleus of eukaryotic cells and survives in the phagocytic amoeba, Acanthamoeba polyphaga. These results suggest that amoeba-like ancestral protozoa could have served as a genetic “melting pot” where the ancestors of rickettsiae and other bacteria promiscuously exchanged genes, eventually leading to their adaptation to the intracellular lifestyle within eukaryotic cells.     

The evolution of spinnable cotton fiber entailed prolonged development and a novel metabolism

A central question in evolutionary biology concerns the developmental processes by which new phenotypes arise. An exceptional example of evolutionary innovation is the single-celled seed trichome in Gossypium (“cotton fiber”). We have used fiber development in Gossypium as a system to understand how morphology can rapidly evolve. Fiber has undergone considerable morphological changes between the short, tightly adherent fibers of G. longicalyx and the derived long, spinnable fibers of its closest relative, G. herbaceum, which facilitated cotton domestication. We conducted comparative gene expression profiling across a developmental time-course of fibers from G. longicalyx and G. herbaceum using microarrays with ~22,000 genes. Expression changes between stages were temporally protracted in G. herbaceum relative to G. longicalyx, reflecting a prolongation of the ancestral developmental program. Gene expression and GO analyses showed that many genes involved with stress responses were upregulated early in G. longicalyx fiber development. Several candidate genes upregulated in G. herbaceum have been implicated in regulating redox levels and cell elongation processes. Three genes previously shown to modulate hydrogen peroxide levels were consistently expressed in domesticated and wild cotton species with long fibers, but expression was not detected by quantitative real time-PCR in wild species with short fibers. Hydrogen peroxide is important for cell elongation, but at high concentrations it becomes toxic, activating stress processes that may lead to early onset of secondary cell wall synthesis and the end of cell elongation. These observations suggest that the evolution of long spinnable fibers in cotton was accompanied by novel expression of genes assisting in the regulation of reactive oxygen species levels. Our data suggest a model for the evolutionary origin of a novel morphology through differential gene regulation causing prolongation of an ancestral developmental program.     

Chloroplast and Mitochondrial DNA Variation in HORDEUM VULGARE and HORDEUM SPONTANEUM

A survey of restriction fragment polymorphism in Hordeum vulgare and Hordeum spontaneum was made using 17 and 16 hexanucleotide restriction endonucleases on chloroplast (cp) and mitochondrial (mt) DNA, respectively. The plant accessions originated from various places throughout the Fertile Cresent and Mediterranean. The types of changes in cpDNA consisted of nucleotide substitutions and insertions and deletions on the order of 100 base pairs. In contrast, mtDNA has most likely undergone larger insertions and deletions of up to 20 kilobase pairs in addition to rearrangements. Grouping of mtDNA fragment data showed that in some cases geographical affinities existed between the two species, whereas in others there were no clear affinities. Nucleotide diversity estimates derived from the restriction fragment data were used in a number of comparisons of variability. Comparisons of overall mtDNA variability (nucleotide diversity = 9.68 x 10^-4) with cpDNA variability (nucleotide diversity = 6.38 x 10^-4 ) indicated that the former are somewhat more variable. Furthermore, there was no indication that the wild H. spontaneum (cpDNA diversity = 5.57 x 10^-4; mtDNA diversity = 6.04 x 10 ^-4) was more variable than the land races of H. vulgare (cpDNA diversity = 5.88 x 10^-4; mtDNA diversity = 9.79 x 10^-4). In fact, on the basis of mtDNA diversity, H. vulgare was the more variable species. Comparison of organelle nucleotide diversity estimates with an estimate of nuclear nucleotide diversity derived from existing isozyme data provided evidence that both organelle genomes are evolving at a slower rate than the nuclear genome.     

Asymmetrical reinforcement and Wolbachia infection in Drosophila

Reinforcement refers to the evolution of increased mating discrimination against heterospecific individuals in zones of geographic overlap and can be considered a final stage in the speciation process. One the factors that may affect reinforcement is the degree to which hybrid matings result in the permanent loss of genes from a species' gene pool. Matings between females of Drosophila subquinaria and males of D. recens result in high levels of offspring mortality, due to interspecific cytoplasmic incompatibility caused by Wolbachia infection of D. recens. Such hybrid inviability is not manifested in matings between D. recens females and D. subquinaria males. Here we ask whether the asymmetrical hybrid inviability is associated with a corresponding asymmetry in the level of reinforcement. The geographic ranges of D. recens and D. subquinaria were found to overlap across a broad belt of boreal forest in central Canada. Females of D. subquinaria from the zone of sympatry exhibit much stronger levels of discrimination against males of D. recens than do females from allopatric populations. In contrast, such reproductive character displacement is not evident in D. recens, consistent with the expected effects of unidirectional cytoplasmic incompatibility. Furthermore, there is substantial behavioral isolation within D. subquinaria, because females from populations sympatric with D. recens discriminate against allopatric conspecific males, whereas females from populations allopatric with D. recens show no discrimination against any conspecific males. Patterns of general genetic differentiation among populations are not consistent with patterns of behavioral discrimination, which suggests that the behavioral isolation within D. subquinaria results from selection against mating with Wolbachia-infected D. recens. Interspecific cytoplasmic incompatibility may contribute not only to post-mating isolation, an effect already widely recognized, but also to reinforcement, particularly in the uninfected species. The resulting reproductive character displacement not only increases behavioral isolation from the Wolbachia-infected species, but may also lead to behavioral isolation between populations of the uninfected species. Given the widespread occurrence of Wolbachia among insects, it thus appears that there are multiple ways by which these endosymbionts may directly and indirectly contribute to reproductive isolation and speciation.     

Large-scale trends in the evolution of gene structures within 11 animal genomes

We have used the annotations of six animal genomes (Homo sapiens, Mus musculus, Ciona intestinalis, Drosophila melanogaster, Anopheles gambiae, and Caenorhabditis elegans) together with the sequences of five unannotated Drosophila genomes to survey changes in protein sequence and gene structure over a variety of timescales—from the less than 5 million years since the divergence of D. simulans and D. melanogaster to the more than 500 million years that have elapsed since the Cambrian explosion. To do so, we have developed a new open-source software library called CGL (for “Comparative Genomics Library”). Our results demonstrate that change in intron–exon structure is gradual, clock-like, and largely independent of coding-sequence evolution. This means that genome annotations can be used in new ways to inform, corroborate, and test conclusions drawn from comparative genomics analyses that are based upon protein and nucleotide sequence similarities.     

Complete Mitochondrial Genome and Phylogeny of Pleistocene MammothMammuthus primigenius

Phylogenetic relationships between the extinct woolly mammoth(Mammuthus primigenius), and the Asian(Elephas maximus) and African savanna(Loxodonta africana) elephants remain unresolved. Here, we report the sequence of the complete mitochondrial genome (16,842 base pairs) of a woolly mammoth extracted from permafrost-preserved remains from the Pleistocene epoch—the oldest mitochondrial genome sequence determined to date. We demonstrate that well-preserved mitochondrial genome fragments, as long as ~1,600–1700 base pairs, can be retrieved from pre-Holocene remains of an extinct species. Phylogenetic reconstruction of the Elephantinae clade suggests thatM. primigenius andE. maximus are sister species that diverged soon after their common ancestor split from theL. africana lineage. Low nucleotide diversity found between independently determined mitochondrial genomic sequences of woolly mammoths separated geographically and in time suggests that north-eastern Siberia was occupied by a relatively homogeneous population ofM. primigenius throughout the late Pleistocene.     

Als3 is a Candida albicans invasin that binds to cadherins and induces endocytosis by host cells

Candida albicans is the most common cause of hematogenously disseminated and oropharyngeal candidiasis. Both of these diseases are characterized by fungal invasion of host cells. Previously, we have found that C. albicans hyphae invade endothelial cells and oral epithelial cells in vitro by inducing their own endocytosis. Therefore, we set out to identify the fungal surface protein and host cell receptors that mediate this process. We found that the C. albicans Als3 is required for the organism to be endocytosed by human umbilical vein endothelial cells and two different human oral epithelial lines. Affinity purification experiments with wild-type and an als3Δ/als3Δ mutant strain of C. albicans demonstrated that Als3 was required for C. albicans to bind to multiple host cell surface proteins, including N-cadherin on endothelial cells and E-cadherin on oral epithelial cells. Furthermore, latex beads coated with the recombinant N-terminal portion of Als3 were endocytosed by Chinese hamster ovary cells expressing human N-cadherin or E-cadherin, whereas control beads coated with bovine serum albumin were not. Molecular modeling of the interactions of the N-terminal region of Als3 with the ectodomains of N-cadherin and E-cadherin indicated that the binding parameters of Als3 to either cadherin are similar to those of cadherin–cadherin binding. Therefore, Als3 is a fungal invasin that mimics host cell cadherins and induces endocytosis by binding to N-cadherin on endothelial cells and E-cadherin on oral epithelial cells. These results uncover the first known fungal invasin and provide evidence that C. albicans Als3 is a molecular mimic of human cadherins.     

A screen for recessive speciation genes expressed in the gametes of F1 hybrid yeast

Diploid hybrids of Saccharomyces cerevisiae and its closest relative, Saccharomyces paradoxus, are viable, but the sexual gametes they produce are not. One of several possible causes of this gamete inviability is incompatibility between genes from different species—such incompatible genes are usually called “speciation genes.” In diploid F1 hybrids, which contain a complete haploid genome from each species, the presence of compatible alleles can mask the effects of (recessive) incompatible speciation genes. But in the haploid gametes produced by F1 hybrids, recessive speciation genes may be exposed, killing the gametes and thus preventing F1 hybrids from reproducing sexually. Here I present the results of an experiment to detect incompatibilities that kill hybrid gametes. I transferred nine of the 16 S. paradoxus chromosomes individually into S. cerevisiae gametes and tested the ability of each to replace its S. cerevisiae homeolog. All nine chromosomes were compatible, producing nine viable haploid strains, each with 15 S. cerevisiae chromosomes and one S. paradoxus chromosome. Thus, none of these chromosomes contain speciation genes that were capable of killing the hybrid gametes that received them. This is a surprising result that suggests that such speciation genes do not play a major role in yeast speciation.     

Odorant-binding proteins OBP57d and OBP57e affect taste perception and host-plant preference in Drosophila sechellia

Despite its morphological similarity to the other species in the Drosophila melanogaster species complex, D. sechellia has evolved distinct physiological and behavioral adaptations to its host plant Morinda citrifolia, commonly known as Tahitian Noni. The odor of the ripe fruit of M. citrifolia originates from hexanoic and octanoic acid. D. sechellia is attracted to these two fatty acids, whereas the other species in the complex are repelled. Here, using interspecies hybrids between D. melanogaster deficiency mutants and D. sechellia, we showed that the Odorant-binding protein 57e (Obp57e) gene is involved in the behavioral difference between the species. D. melanogaster knock-out flies for Obp57e and Obp57d showed altered behavioral responses to hexanoic acid and octanoic acid. Furthermore, the introduction of Obp57d and Obp57e from D. simulans and D. sechellia shifted the oviposition site preference of D. melanogaster Obp57d/e^KO flies to that of the original species, confirming the contribution of these genes to D. sechellia's specialization to M. citrifolia. Our finding of the genes involved in host-plant determination may lead to further understanding of mechanisms underlying taste perception, evolution of plant–herbivore interactions, and speciation.     

Macronuclear genome sequence of the ciliate Tetrahymena thermophila, a model eukaryote

The ciliate Tetrahymena thermophila is a model organism for molecular and cellular biology. Like other ciliates, this species has separate germline and soma functions that are embodied by distinct nuclei within a single cell. The germline-like micronucleus (MIC) has its genome held in reserve for sexual reproduction. The soma-like macronucleus (MAC), which possesses a genome processed from that of the MIC, is the center of gene expression and does not directly contribute DNA to sexual progeny. We report here the shotgun sequencing, assembly, and analysis of the MAC genome of T. thermophila, which is approximately 104 Mb in length and composed of approximately 225 chromosomes. Overall, the gene set is robust, with more than 27,000 predicted protein-coding genes, 15,000 of which have strong matches to genes in other organisms. The functional diversity encoded by these genes is substantial and reflects the complexity of processes required for a free-living, predatory, single-celled organism. This is highlighted by the abundance of lineage-specific duplications of genes with predicted roles in sensing and responding to environmental conditions (e.g., kinases), using diverse resources (e.g., proteases and transporters), and generating structural complexity (e.g., kinesins and dyneins). In contrast to the other lineages of alveolates (apicomplexans and dinoflagellates), no compelling evidence could be found for plastid-derived genes in the genome. UGA, the only T. thermophila stop codon, is used in some genes to encode selenocysteine, thus making this organism the first known with the potential to translate all 64 codons in nuclear genes into amino acids. We present genomic evidence supporting the hypothesis that the excision of DNA from the MIC to generate the MAC specifically targets foreign DNA as a form of genome self-defense. The combination of the genome sequence, the functional diversity encoded therein, and the presence of some pathways missing from other model organisms makes T. thermophila an ideal model for functional genomic studies to address biological, biomedical, and biotechnological questions of fundamental importance.     

fog-2 and the evolution of self-fertile hermaphroditism in Caenorhabditis

Somatic and germline sex determination pathways have diverged significantly in animals, making comparisons between taxa difficult. To overcome this difficulty, we compared the genes in the germline sex determination pathways of Caenorhabditis elegans and C. briggsae, two Caenorhabditis species with similar reproductive systems and sequenced genomes. We demonstrate that C. briggsae has orthologs of all known C. elegans sex determination genes with one exception: fog-2. Hermaphroditic nematodes are essentially females that produce sperm early in life, which they use for self fertilization. In C. elegans, this brief period of spermatogenesis requires FOG-2 and the RNA-binding protein GLD-1, which together repress translation of the tra-2 mRNA. FOG-2 is part of a large C. elegans FOG-2-related protein family defined by the presence of an F-box and Duf38/FOG-2 homogy domain. A fog-2-related gene family is also present in C. briggsae, however, the branch containing fog-2 appears to have arisen relatively recently in C. elegans, post-speciation. The C-terminus of FOG-2 is rapidly evolving, is required for GLD-1 interaction, and is likely critical for the role of FOG-2 in sex determination. In addition, C. briggsae gld-1 appears to play the opposite role in sex determination (promoting the female fate) while maintaining conserved roles in meiotic progression during oogenesis. Our data indicate that the regulation of the hermaphrodite germline sex determination pathway at the level of FOG-2/GLD-1/tra-2 mRNA is fundamentally different between C. elegans and C. briggsae, providing functional evidence in support of the independent evolution of self-fertile hermaphroditism. We speculate on the convergent evolution of hermaphroditism in Caenorhabditis based on the plasticity of the C. elegans germline sex determination cascade, in which multiple mutant paths yield self fertility.     

Phylogenomics of the reproductive parasite Wolbachia pipientis wMel: a streamlined genome overrun by mobile genetic elements

The complete sequence of the 1,267,782 bp genome of Wolbachia pipientis wMel, an obligate intracellular bacteria of Drosophila melanogaster, has been determined. Wolbachia, which are found in a variety of invertebrate species, are of great interest due to their diverse interactions with different hosts, which range from many forms of reproductive parasitism to mutualistic symbioses. Analysis of the wMel genome, in particular phylogenomic comparisons with other intracellular bacteria, has revealed many insights into the biology and evolution of wMel and Wolbachia in general. For example, the wMel genome is unique among sequenced obligate intracellular species in both being highly streamlined and containing very high levels of repetitive DNA and mobile DNA elements. This observation, coupled with multiple evolutionary reconstructions, suggests that natural selection is somewhat inefficient in wMel, most likely owing to the occurrence of repeated population bottlenecks. Genome analysis predicts many metabolic differences with the closely related Rickettsia species, including the presence of intact glycolysis and purine synthesis, which may compensate for an inability to obtain ATP directly from its host, as Rickettsia can. Other discoveries include the apparent inability of wMel to synthesize lipopolysaccharide and the presence of the most genes encoding proteins with ankyrin repeat domains of any prokaryotic genome yet sequenced. Despite the ability of wMel to infect the germline of its host, we find no evidence for either recent lateral gene transfer between wMel and D. melanogaster or older transfers between Wolbachia and any host. Evolutionary analysis further supports the hypothesis that mitochondria share a common ancestor with the α-Proteobacteria, but shows little support for the grouping of mitochondria with species in the order Rickettsiales. With the availability of the complete genomes of both species and excellent genetic tools for the host, the wMel–D. melanogaster symbiosis is now an ideal system for studying the biology and evolution of Wolbachia infections.