Large-scale comparison of fungal sequence information: mechanisms of innovation in Neurospora crassa and gene loss in Saccharomyces cerevisiae

The Neurospora crassa molybdate transporter: Characterizing a novel transporter homologous to the plant MOT1 family

Molybdenum (Mo) is an essential element for animals, plants, and fungi. To achieve biological activity in eukaryotes, Mo must be complexed into the molybdenum cofactor (Moco). Cells are known to take up Mo in the form of the oxyanion molybdate. However, molybdate transporters are scarcely characterized in the fungal kingdom. In plants and algae, molybdate is imported into the cell via two families of molybdate transporters (MOT), MOT1 and MOT2. For the filamentous fungus Neurospora crassa, a sequence homologous to the MOT1 family was previously annotated. Here we report a characterization of this molybdate-related transporter, encoded by the ncmot-1 gene. We found that the deletion of ncmot-1 leads to an accumulation of total Mo within the mycelium and a roughly 51 % higher tolerance against high molybdate levels when grown on ammonium medium. The localization of a GFP tagged NcMOT-1 was identified among the vacuolar membrane. Thereby, we propose NcMOT-1 as an exporter, transporting molybdate out of the vacuole into the cytoplasm. Lastly, the heterologous expression of NcMOT-1 in Saccharomyces cerevisiae verifies the functionality of this protein as a MOT. Our results open the way towards understanding molybdate transport as part of Mo homeostasis and Moco-biosynthesis in fungi.

Structural basis for the interaction between the first SURP domain of the SF3A1 subunit in U2 snRNP and the human splicing factor SF1

SURP domains are exclusively found in splicing-related proteins in all eukaryotes. SF3A1, a component of the U2 snRNP, has two tandem SURP domains, SURP1, and SURP2. SURP2 is permanently associated with a specific short region of SF3A3 within the SF3A protein complex whereas, SURP1 binds to the splicing factor SF1 for recruitment of U2 snRNP to the early spliceosomal complex, from which SF1 is dissociated during complex conversion. Here, we determined the solution structure of the complex of SURP1 and the human SF1 fragment using nuclear magnetic resonance (NMR) methods. SURP1 adopts the canonical topology of α1-α2-3₁₀ -α3, in which α1 and α2 are connected by a single glycine residue in a particular backbone conformation, allowing the two α-helices to be fixed at an acute angle. A hydrophobic patch, which is part of the characteristic surface formed by α1 and α2, specifically contacts a hydrophobic cluster on a 16-residue α-helix of the SF1 fragment. Furthermore, whereas only hydrophobic interactions occurred between SURP2 and the SF3A3 fragment, several salt bridges and hydrogen bonds were found between the residues of SURP1 and the SF1 fragment. This finding was confirmed through mutational studies using bio-layer interferometry. The study also revealed that the dissociation constant between SURP1 and the SF1 fragment peptide was approximately 20 μM, indicating a weak or transient interaction. Collectively, these results indicate that the interplay between U2 snRNP and SF1 involves a transient interaction of SURP1, and this transient interaction appears to be common to most SURP domains, except for SURP2.

Phylogenomic Analyses of 2,786 Genes in 158 Lineages Support a Root of the Eukaryotic Tree of Life between Opisthokonts and All Other Lineages

Advances in phylogenomics and high-throughput sequencing have allowed the reconstruction of deep phylogenetic relationships in the evolution of eukaryotes. Yet, the root of the eukaryotic tree of life remains elusive. The most popular hypothesis in textbooks and reviews is a root between Unikonta (Opisthokonta + Amoebozoa) and Bikonta (all other eukaryotes), which emerged from analyses of a single-gene fusion. Subsequent, highly cited studies based on concatenation of genes supported this hypothesis with some variations or proposed a root within Excavata. However, concatenation of genes does not consider phylogenetically-informative events like gene duplications and losses. A recent study using gene tree parsimony (GTP) suggested the root lies between Opisthokonta and all other eukaryotes, but only including 59 taxa and 20 genes. Here we use GTP with a duplication-loss model in a gene-rich and taxon-rich dataset (i.e., 2,786 gene families from two sets of 155 and 158 diverse eukaryotic lineages) to assess the root, and we iterate each analysis 100 times to quantify tree space uncertainty. We also contrasted our results and discarded alternative hypotheses from the literature using GTP and the likelihood-based method SpeciesRax. Our estimates suggest a root between Fungi or Opisthokonta and all other eukaryotes; but based on further analysis of genome size, we propose that the root between Opisthokonta and all other eukaryotes is the most likely.

Multiple independent origins for a subtelomeric locus associated with growth rate in Fusarium circinatum

Fusarium is a diverse assemblage that includes a large number of species of considerable medical and agricultural importance. Not surprisingly, whole genome sequences for many Fusarium species have been published or are in the process of being determined, the availability of which is invaluable for deciphering the genetic basis of key phenotypic traits. Here we investigated the distribution, genic composition, and evolutionary history of a locus potentially determining growth rate in the pitch canker pathogen F. circinatum. We found that the genomic region underlying this locus is highly conserved amongst F. circinatum and its close relatives, except for the presence of a 12 000 base pair insertion in all of the examined isolates of F. circinatum. This insertion encodes for five genes and our phylogenetic analyses revealed that each was most likely acquired through horizontal gene transfer from polyphyletic origins. Our data further showed that this region is located in a region low in G+C content and enriched for repetitive sequences and transposable elements, which is situated near the telomere of Chromosome 3 of F. circinatum. As have been shown for other fungi, these findings thus suggest that the emergence of the unique 12 000 bp region in F. circinatum is linked to the dynamic evolutionary processes associated with subtelomeres that, in turn, have been implicated in the ecological adaptation of fungal pathogens.

The calcium-binding protein EpANN from the lichenized fungus Endocarpon pusillum enhances stress tolerance in yeast and plants

Annexins are calcium-phospholipid binding proteins that play a significant role in the Ca²⁺signaling pathway. These proteins are essential for plants to effectively respond to abiotic stresses. However, their functions and mechanisms remain largely unknown in fungi. In this study, an annexin gene, Epann, was cloned from the lichenized fungus Endocarpon pusillum, a drought resistant organism. Our results showed that Epann was induced by several abiotic stresses in E. pusillum. Heterologous expression of the Epann gene enhanced the stress tolerance of Saccharomyces cerevisiae. Under heat-shock conditions, the EpANN proteins were significantly aggregated and the aggregation sites were located on peroxisomes. In heat-shocked cells, Epann reduced the reactive oxygen species level mainly through its intracellular peroxidase activity and regulation of stress-related genes. Transgenic Arabidopsis plants overexpressing Epann exhibited a higher germination rate under oxidative stress and stronger drought tolerance. Our results provide a mechanistic understanding of the role of annexins in abiotic stress responses and suggest that this lichenized fungal gene could be a promising resource to generate stress-tolerant transgenic organisms.

Genomic and transcriptomic analyses of the Chinese Maotai-flavored liquor yeast MT1 revealed its unique multi-carbon co-utilization

BACKGROUND

Revealing genetic mechanisms behind specific physiological characteristics of Saccharomyces cerevisiae from specific environments is important for industrial applications and requires precise understanding.

RESULTS

Maotai strain MT1 was isolated from the complicated Chinese Maotai-flavored liquor-making environment with extremely high temperatures, and acidic and ethanol stresses. Compared with the type strain S288c, MT1 can tolerate high acidity (pH 2.0), high ethanol levels (16 %) and high temperatures (44 °C). In addition, MT1 can simultaneously utilize various sugars, including glucose, sucrose, galactose, maltose, melibiose, trehalose, raffinose and turanose. Genomic comparisons identified a distinct MT1 genome, 0.5 Mb smaller than that of S288c. There are 145 MT1-specific genes that are not in S288c, including MEL1, MAL63, KHR1, BIO1 and BIO6. A transcriptional comparison indicated that HXT5 and HXT13, which are theoretically repressed by glucose, were no longer inhibited in MT1 and were highly expressed even in a medium containing 70 g/L glucose. Thus, other sugars may be co-utilized with glucose by MT1 without diauxic growth.

CONCLUSIONS

Based on a functional genomics analysis, we revealed the genetic basis and evolutionary mechanisms underlying the traits of the Chinese Maotai-flavored yeast MT1. This work provides new insights for the genetic breeding of yeast and also enriches the genetic resources of yeast.

Fungal annexins: a mini review

The large family of annexins is composed of more than a thousand members which are typically phospholipid-binding proteins. Annexins act in a number of signalling networks and membrane trafficking events which are fundamental to cell physiology. Annexins exert their functions mainly through their calcium-dependent membrane binding abilities; however, some calcium-independent interactions have been documented in the literature. Although mammalian and plant annexins have been well characterized, little is known about this family in fungi. This mini review summarizes the available data on fungal annexins.

Differential gene retention as an evolutionary mechanism to generate biodiversity and adaptation in yeasts

The evolutionary history of the characters underlying the adaptation of microorganisms to food and biotechnological uses is poorly understood. We undertook comparative genomics to investigate evolutionary relationships of the dairy yeast Geotrichum candidum within Saccharomycotina. Surprisingly, a remarkable proportion of genes showed discordant phylogenies, clustering with the filamentous fungus subphylum (Pezizomycotina), rather than the yeast subphylum (Saccharomycotina), of the Ascomycota. These genes appear not to be the result of Horizontal Gene Transfer (HGT), but to have been specifically retained by G. candidum after the filamentous fungi-yeasts split concomitant with the yeasts' genome contraction. We refer to these genes as SRAGs (Specifically Retained Ancestral Genes), having been lost by all or nearly all other yeasts, and thus contributing to the phenotypic specificity of lineages. SRAG functions include lipases consistent with a role in cheese making and novel endoglucanases associated with degradation of plant material. Similar gene retention was observed in three other distantly related yeasts representative of this ecologically diverse subphylum. The phenomenon thus appears to be widespread in the Saccharomycotina and argues that, alongside neo-functionalization following gene duplication and HGT, specific gene retention must be recognized as an important mechanism for generation of biodiversity and adaptation in yeasts.

The Neurospora crassa chr-1 gene is up-regulated by chromate and its encoded CHR-1 protein causes chromate sensitivity and chromium accumulation

The ChrA membrane protein belongs to the CHR superfamily of chromate ion transporters, which includes homologues from bacteria, archaea and eukaryotes. Bacterial ChrA homologues confer chromate resistance by exporting chromate ions from the cell's cytoplasm. The Neurospora crassa strain 74-A chr-1 gene encodes a putative CHR-1 protein of 507 amino acid residues, which belongs to the CHR superfamily. RT-PCR assays showed that expression of the chr-1 gene was up-regulated by chromate exposure of N. crassa cultures. Introduction in N. crassa of sense and antisense fragments of the chr-1 gene, as part of a silencing module within the pSilent-1 vector, produced transformants with a phenotype of resistance to chromate and diminished accumulation of chromium, as compared with the control strain containing only the vector. A chromate-resistance phenotype was also observed in N crassa strains deleted in the genomic chr-1 gene, thus confirming that the absence of CHR-1 protein confers chromate resistance to the fungus. The cDNA from N. crassa chr-1 gene (Ncchr-1) was cloned into the pYES2 vector under the control of a GAL promoter and the resulting recombinant plasmid was transferred to the yeast Saccharomyces cerevisiae. Galactose-induced S. cerevisiae transformants expressing Ncchr-1 were more sensitive to chromate and accumulated 2.5 times more chromium than the induced strain containing only the vector. Excess sulfate, a chromate analog, was unable to protect S. cerevisiae chr-1 transformants from chromate toxicity. These data indicate that the N. crassa CHR-1 protein functions as a transporter that takes up chromate; it also appears that this transport occurs in a sulfate-independent fashion. This is the first report assigning a role as a chromate transporter to a nonbacterial CHR protein.

Sequence Analysis of SSR-Flanking Regions Identifies Genome Affinities between Pasture Grass Fungal Endophyte Taxa

Fungal species of the Neotyphodium and Epichloë genera are endophytes of pasture grasses showing complex differences of life-cycle and genetic architecture. Simple sequence repeat (SSR) markers have been developed from endophyte-derived expressed sequence tag (EST) collections. Although SSR array size polymorphisms are appropriate for phenetic analysis to distinguish between taxa, the capacity to resolve phylogenetic relationships is limited by both homoplasy and heteroploidy effects. In contrast, nonrepetitive sequence regions that flank SSRs have been effectively implemented in this study to demonstrate a common evolutionary origin of grass fungal endophytes. Consistent patterns of relationships between specific taxa were apparent across multiple target loci, confirming previous studies of genome evolution based on variation of individual genes. Evidence was obtained for the definition of endophyte taxa not only through genomic affinities but also by relative gene content. Results were compatible with the current view that some asexual Neotyphodium species arose following interspecific hybridisation between sexual Epichloë ancestors. Phylogenetic analysis of SSR-flanking regions, in combination with the results of previous studies with other EST-derived SSR markers, further permitted characterisation of Neotyphodium isolates that could not be assigned to known taxa on the basis of morphological characteristics.

Cloning and characterization of Rhodotorula glutinis thymine hydroxylase

Thymine hydroxylase (TH) is a member of the alpha-ketoglutarate-dependent nonheme iron dioxygenase family that includes a series of DNA repair proteins including alkB. Substantial interest in this family of enzymes derives from their capacity to modify DNA bases and precursors by oxidation. Previously, a sequence has been published for cloned Rhodotorula glutinis TH. However, the minimal reported activity of this enzyme, coupled with inconsistencies with previously published mass spectrometry data, compelled us to reexamine TH. The sequence reported here differs from the previously reported sequence at two amino acid positions and is consistent with previously reported mass spectrometry data. The cloned enzyme characterized in this report displayed substantial activity, indicating that the sequence differences are critical for activity. The substrate selectivity of TH against a series of pyrimidine analogues is consistent with that reported for the wild-type enzyme and, in part, explains the mode of selection of uracil analogues. A preliminary model of the active site has been constructed for the purposes of comparing TH with other members of this family. TH and alkB share in common the capacity to oxidize N-methyl groups. However, TH has the added capacity to oxidize the 5-methyl group of thymine, a property that is potentially important for enzymes that could act on DNA and modify DNA-protein interactions.

Identification of potential marker genes for Trichoderma harzianum strains with high antagonistic potential against Rhizoctonia solani by a rapid subtraction hybridization approach

A rapid subtraction hybridization approach was used to isolate genes differentially expressed during mycelial contact between Trichoderma harzianum (Hypocrea lixii) and Rhizoctonia solani, and could serve as marker genes for selection of superior biocontrol strains. Putatively positive clones were evaluated by transcription analysis during mycelial contact with R. solani versus growth on glucose, and for their differential transcription between two strains with either strong or poor biocontrol capability before, at, and after contact with R. solani. Besides four clones, which had similarity to putative but as yet uncharacterized proteins, they comprised ribosomal proteins, proteins involved in transcriptional switch and regulation, amino acid and energy catabolism, multidrug resistance, and degradation of proteins and glucans. Transcription of three clones was evaluated in five T. harzianum strains under confrontation conditions with R. solani. Two clones-acetyl-xylane esterase AXE1 and endoglucanase Cel61b-showed significant upregulation during in vivo confrontation of a T. harzianum strain that successively demonstrated a very high antagonistic capability towards R. solani, while expression was progressively lower in a series of T. harzianum strains with intermediate to poor antagonistic activity. These clones are promising candidates for use as markers in the screening of improved T. harzianum biocontrol strains.

Evolution of SET-domain protein families in the unicellular and multicellular Ascomycota fungi

BACKGROUND

The evolution of multicellularity is accompanied by the occurrence of differentiated tissues, of organismal developmental programs, and of mechanisms keeping the balance between proliferation and differentiation. Initially, the SET-domain proteins were associated exclusively with regulation of developmental genes in metazoa. However, finding of SET-domain genes in the unicellular yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe suggested that SET-domain proteins regulate a much broader variety of biological programs. Intuitively, it is expected that the numbers, types, and biochemical specificity of SET-domain proteins of multicellular versus unicellular forms would reflect the differences in their biology. However, comparisons across the unicellular and multicellular domains of life are complicated by the lack of knowledge of the ancestral SET-domain genes. Even within the crown group, different biological systems might use the epigenetic 'code' differently, adapting it to organism-specific needs. Simplifying the model, we undertook a systematic phylogenetic analysis of one monophyletic fungal group (Ascomycetes) containing unicellular yeasts, Saccharomycotina (hemiascomycetes), and a filamentous fungal group, Pezizomycotina (euascomycetes).

RESULTS

Systematic analysis of the SET-domain genes across an entire eukaryotic phylum has outlined clear distinctions in the SET-domain gene collections in the unicellular and in the multicellular (filamentous) relatives; diversification of SET-domain gene families has increased further with the expansion and elaboration of multicellularity in animal and plant systems. We found several ascomycota-specific SET-domain gene groups; each was unique to either Saccharomycotina or Pezizomycotina fungi. Our analysis revealed that the numbers and types of SET-domain genes in the Saccharomycotina did not reflect the habitats, pathogenicity, mechanisms of sexuality, or the ability to undergo morphogenic transformations. However, novel genes have appeared for functions associated with the transition to multicellularity. Descendents of most of the SET-domain gene families found in the filamentous fungi could be traced in the genomes of extant animals and plants, albeit as more complex structural forms.

CONCLUSION

SET-domain genes found in the filamentous species but absent from the unicellular sister group reflect two alternative evolutionary events: deletion from the yeast genomes or appearance of novel structures in filamentous fungal groups. There were no Ascomycota-specific SET-domain gene families (i.e., absent from animal and plant genomes); however, plants and animals share SET-domain gene subfamilies that do not exist in the fungi. Phylogenetic and gene-structure analyses defined several animal and plant SET-domain genes as sister groups while those of fungal origin were basal to them. Plants and animals also share SET-domain subfamilies that do not exist in fungi.

GENOMIC INSIGHTS INTO EVOLUTIONARY RELATIONSHIPS AMONG HETEROKONT LINEAGES EMPHASIZING THE PHAEOPHYCEAE(1)

Heterokonts comprise a large and diverse group of organisms unified by the heterokont biflagellate condition. Monophyly of many of these lineages is well established, but evolutionary relationships among the various lineages remain elusive. Among these lineages, the brown algae (Phaeophyceae) are a monophyletic, taxonomically diverse, and ecologically critical group common to marine environments. Despite their biological and scientific importance, consensus regarding brown algal phylogeny and taxonomic relationships is missing. Our long-term research goal is to produce a well-resolved taxon-rich phylogeny of the class to assess evolutionary patterns and taxonomic relationships among brown algal lineages and their relationship to other closely related heterokont groups. To accomplish this goal and augment existing loci for phaeophycean-wide systematic studies, we generated expressed sequence tags (ESTs) from several major brown algal lineages and from the heterokont lineage representing the closest sister group to brown algae. To date, we have successfully constructed cDNA libraries for two lineages (Choristocarpus tenellus Zanardini and Schizocladia ischiensis E. C. Henry, Okuda et H. Kawai) and in the library test phase obtained up to 1,600 ESTs per organism. Annotation results showed a gene discovery rate of 45%-50% for each library revealing 500-700 unique genes from each organism. We have identified several potential genes for phylogenetic inference and used these loci for preliminary molecular clock analyses. Our molecular clock analysis suggests that the basal divergence in brown algae occurred around the time of the pennate-centric diatom divergence. Here we report this analysis and other uses of ESTs in brown algal phylogenomics and the utility of these data for resolving the phylogeny of this group.

Dynamic behavior of an intrinsically unstructured linker domain is conserved in the face of negligible amino acid sequence conservation

Proteins or regions of proteins that do not form compact globular structures are classified as intrinsically unstructured proteins (IUPs). IUPs are common in nature and have essential molecular functions, but even a limited understanding of the evolution of their dynamic behavior is lacking. The primary objective of this work was to test the evolutionary conservation of dynamic behavior for a particular class of IUPs that form intrinsically unstructured linker domains (IULD) that tether flanking folded domains. This objective was accomplished by measuring the backbone flexibility of several IULD homologues using nuclear magnetic resonance (NMR) spectroscopy. The backbone flexibility of five IULDs, representing three kingdoms, was measured and analyzed. Two IULDs from animals, one IULD from fungi, and two IULDs from plants showed similar levels of backbone flexibility that were consistent with the absence of a compact globular structure. In contrast, the amino acid sequences of the IULDs from these three taxa showed no significant similarity. To investigate how the dynamic behavior of the IULDs could be conserved in the absence of detectable sequence conservation, evolutionary rate studies were performed on a set of nine mammalian IULDs. The results of this analysis showed that many sites in the IULD are evolving neutrally, suggesting that dynamic behavior can be maintained in the absence of natural selection. This work represents the first experimental test of the evolutionary conservation of dynamic behavior and demonstrates that amino acid sequence conservation is not required for the conservation of dynamic behavior and presumably molecular function.

Dynamic behavior of an intrinsically unstructured linker domain is conserved in the face of negligible amino acid sequence conservation

Proteins or regions of proteins that do not form compact globular structures are classified as intrinsically unstructured proteins (IUPs). IUPs are common in nature and have essential molecular functions, but even a limited understanding of the evolution of their dynamic behavior is lacking. The primary objective of this work was to test the evolutionary conservation of dynamic behavior for a particular class of IUPs that form intrinsically unstructured linker domains (IULD) that tether flanking folded domains. This objective was accomplished by measuring the backbone flexibility of several IULD homologues using nuclear magnetic resonance (NMR) spectroscopy. The backbone flexibility of five IULDs, representing three kingdoms, was measured and analyzed. Two IULDs from animals, one IULD from fungi, and two IULDs from plants showed similar levels of backbone flexibility that were consistent with the absence of a compact globular structure. In contrast, the amino acid sequences of the IULDs from these three taxa showed no significant similarity. To investigate how the dynamic behavior of the IULDs could be conserved in the absence of detectable sequence conservation, evolutionary rate studies were performed on a set of nine mammalian IULDs. The results of this analysis showed that many sites in the IULD are evolving neutrally, suggesting that dynamic behavior can be maintained in the absence of natural selection. This work represents the first experimental test of the evolutionary conservation of dynamic behavior and demonstrates that amino acid sequence conservation is not required for the conservation of dynamic behavior and presumably molecular function.

The evolution of seminal ribonuclease: pseudogene reactivation or multiple gene inactivation events?

Two approaches, one novel, are applied to analyze the divergent evolution of ruminant seminal ribonucleases (RNases), paralogs of the well-known pancreatic RNases of mammals. Here, the goal was to identify periods of divergence of seminal RNase under functional constraints, periods of divergence as a pseudogene, and periods of divergence driven by positive selection pressures. The classical approach involves the analysis of nonsynonymous to synonymous replacements ratios (omega) for the branches of the seminal RNase evolutionary tree. The novel approach coupled these analyses with the mapping of substitutions on the folded structure of the protein. These analyses suggest that seminal RNase diverged during much of its history after divergence from pancreatic RNase as a functioning protein, followed by homoplastic inactivations to create pseudogenes in multiple ruminant lineages. Further, they are consistent with adaptive evolution only in the most recent episode leading to the gene in modern oxen. These conclusions contrast sharply with the view, cited widely in the literature, that seminal RNase decayed after its formation by gene duplication into an inactive pseudogene, whose lesions were repaired in a reactivation event. Further, the 2 approaches, omega estimation and mapping of replacements on the protein structure, were compared by examining their utility for establishing the functional status of the seminal RNase genes in 2 deer species. Hog and roe deer share common lesions, which strongly suggests that the gene was inactive in their last common ancestor. In this specific example, the crystallographic approach made the correct implication more strongly than the omega approach. Studies of this type should contribute to an integrated framework of tools to assign functional and nonfunctional episodes to recently created gene duplicates and to understand more broadly how gene duplication leads to the emergence of proteins with novel functions.

Comparative EST analysis provides insights into the basal aquatic fungus Blastocladiella emersonii

Background

Blastocladiella emersonii is an aquatic fungus of the Chytridiomycete class, which is at the base of the fungal phylogenetic tree. In this sense, some ancestral characteristics of fungi and animals or fungi and plants could have been retained in this aquatic fungus and lost in members of late-diverging fungal species. To identify in B. emersonii sequences associated with these ancestral characteristics two approaches were followed: (1) a large-scale comparative analysis between putative unigene sequences (uniseqs) from B. emersonii and three databases constructed ad hoc with fungal proteins, animal proteins and plant unigenes deposited in Genbank, and (2) a pairwise comparison between B. emersonii full-length cDNA sequences and their putative orthologues in the ascomycete Neurospora crassa and the basidiomycete Ustilago maydis.

Results

Comparative analyses of B. emersonii uniseqs with fungi, animal and plant databases through the two approaches mentioned above produced 166 B. emersonii sequences, which were identified as putatively absent from other fungi or not previously described. Through these approaches we found: (1) possible orthologues of genes previously identified as specific to animals and/or plants, and (2) genes conserved in fungi, but with a large difference in divergence rate in B. emersonii. Among these sequences, we observed cDNAs encoding enzymes from coenzyme B₁₂-dependent propionyl-CoA pathway, a metabolic route not previously described in fungi, and validated their expression in Northern blots.

Conclusion

Using two different approaches involving comparative sequence analyses, we could identify sequences from the early-diverging fungus B. emersonii previously considered specific to animals or plants, and highly divergent sequences from the same fungus relative to other fungi.

Plant pathogen forensics: capabilities, needs, and recommendations

A biological attack on U.S. crops, rangelands, or forests could reduce yield and quality, erode consumer confidence, affect economic health and the environment, and possibly impact human nutrition and international relations. Preparedness for a crop bioterror event requires a strong national security plan that includes steps for microbial forensics and criminal attribution. However, U.S. crop producers, consultants, and agricultural scientists have traditionally focused primarily on strategies for prevention and management of diseases introduced naturally or unintentionally rather than on responding appropriately to an intentional pathogen introduction. We assess currently available information, technologies, and resources that were developed originally to ensure plant health but also could be utilized for postintroduction plant pathogen forensics. Recommendations for prioritization of efforts and resource expenditures needed to enhance our plant pathogen forensics capabilities are presented.

Comparative genomic analysis of phytopathogenic fungi using expressed sequence tag (EST) collections

SUMMARY We describe the analysis of 57 727 unique expressed sequence tags (ESTs) from 15 species of phytopathogenic and three species of saprophytic fungi. This resource is held within the COGEME phytopathogen EST database (http://cogeme.ex.ac.uk/). Comparative analysis was performed to investigate the differences between pathogenic and free-living fungi based on a substantial collection of expressed gene sequences and available, completed fungal genome sequences. We report that the expressed gene inventories of pathogenic fungi were not significantly more similar to each other than to those of free-living filamentous fungi. As expected, however, filamentous fungi as a group share more sequences in common than with the free-living yeast species Saccharomyces cerevisiae. Interestingly, ESTs of the obligate biotrophic fungus Blumeria graminis f. sp. hordei were more dissimilar to those of all other fungal species assessed, having a lower number of sequences in common with filamentous ascomycetes studied to date and also possessing a larger proportion of unisequences of unknown function. Our analysis of ESTs in the COGEME database enabled identification of a set of functional groups of genes that are more highly represented in the genomes of pathogenic fungi than non-pathogenic species.

Contribution of horizontal gene transfer to the evolution of Saccharomyces cerevisiae

The genomes of the hemiascomycetes Saccharomyces cerevisiae and Ashbya gossypii have been completely sequenced, allowing a comparative analysis of these two genomes, which reveals that a small number of genes appear to have entered these genomes as a result of horizontal gene transfer from bacterial sources. One potential case of horizontal gene transfer in A. gossypii and 10 potential cases in S. cerevisiae were identified, of which two were investigated further. One gene, encoding the enzyme dihydroorotate dehydrogenase (DHOD), is potentially a case of horizontal gene transfer, as shown by sequencing of this gene from additional bacterial and fungal species to generate sufficient data to construct a well-supported phylogeny. The DHOD-encoding gene found in S. cerevisiae, URA1 (YKL216W), appears to have entered the Saccharomycetaceae after the divergence of the S. cerevisiae lineage from the Candida albicans lineage and possibly since the divergence from the A. gossypii lineage. This gene appears to have come from the Lactobacillales, and following its acquisition the endogenous eukaryotic DHOD gene was lost. It was also shown that the bacterially derived horizontally transferred DHOD is required for anaerobic synthesis of uracil in S. cerevisiae. The other gene discussed in detail is BDS1, an aryl- and alkyl-sulfatase gene of bacterial origin that we have shown allows utilization of sulfate from several organic sources. Among the eukaryotes, this gene is found in S. cerevisiae and Saccharomyces bayanus and appears to derive from the alpha-proteobacteria.

A phylogenomic approach to reconstructing the diversification of serine proteases in fungi

Using a phylogenomic approach with 10 fungi of very different virulence and habitat, we determined that there was substantial diversification of subtilase-type proteases early in ascomycete history (with subsequent loss in many lineages) but with no comparable diversification of trypsins. Patterns of intron loss and the degree of divergence between paralogues demonstrated that the proliferation of proteinase K subtilases and subtilisin type subtilases seen in pathogenic ascomycetes (Metarhizium anisopliae, Magnaporthe grisea, Fusarium graminearum) occurred after the basidiomycete/ascomycete split but predated radiation of ascomycete lineages. This suggests that the early ascomycetes had a lifestyle that selected for multiple proteases, whereas the current disparity in gene numbers between ascomycete lineages results from retention of genes in at least some pathogens that have been lost in other lineages (yeasts, Aspergillus nidulans, Neurospora crassa). A similar prevailing trend towards lineage specific gene loss of trypsins in saprophytes and some pathogens suggests that their phylogenetic breadth will have been much wider in early fungi than currently.

Evaluation of monocot and eudicot divergence using the sugarcane transcriptome

Over 40,000 sugarcane (Saccharum officinarum) consensus sequences assembled from 237,954 expressed sequence tags were compared with the protein and DNA sequences from other angiosperms, including the genomes of Arabidopsis and rice (Oryza sativa). Approximately two-thirds of the sugarcane transcriptome have similar sequences in Arabidopsis. These sequences may represent a core set of proteins or protein domains that are conserved among monocots and eudicots and probably encode for essential angiosperm functions. The remaining sequences represent putative monocot-specific genetic material, one-half of which were found only in sugarcane. These monocot-specific cDNAs represent either novelties or, in many cases, fast-evolving sequences that diverged substantially from their eudicot homologs. The wide comparative genome analysis presented here provides information on the evolutionary changes that underlie the divergence of monocots and eudicots. Our comparative analysis also led to the identification of several not yet annotated putative genes and possible gene loss events in Arabidopsis.

EST Analysis of the Cnidarian Acropora millepora Reveals Extensive Gene Loss and Rapid Sequence Divergence in the Model Invertebrates

A significant proportion of mammalian genes are not represented in the genomes of Drosophila, Caenorhabditis or Saccharomyces, and many of these are assumed to have been vertebrate innovations. To test this assumption, we conducted a preliminary EST project on the anthozoan cnidarian, Acropora millepora, a basal metazoan. More than 10% of the Acropora ESTs with strong metazoan matches to the databases had clear human homologs but were not represented in the Drosophila or Caenorhabditis genomes; this category includes a surprising diversity of transcription factors and metabolic proteins that were previously assumed to be restricted to vertebrates. Consistent with higher rates of divergence in the model invertebrates, three-way comparisons show that most Acropora ESTs match human sequences much more strongly than they do any Drosophila or Caenorhabditis sequence. Gene loss has thus been much more extensive in the model invertebrate lineages than previously assumed and, as a consequence, some genes formerly thought to be vertebrate inventions must have been present in the common metazoan ancestor. The complexity of the Acropora genome is paradoxical, given that this organism contains apparently few tissue types and the simplest extant nervous system consisting of a morphologically homogeneous nerve net.

Correlations between gene expression and gene conservation in fission yeast

Genes can be expressed at a wide range of levels, and they show different degrees of cross-species conservation. We compared gene expression levels to gene conservation by integrating microarray data from fission yeast (Schizosaccharomyces pombe) with lists of "core" genes (present in worm and budding and fission yeasts), "yeast-specific" genes (present in budding and fission yeasts, but not in worm), and "pombe-specific" genes (present in fission yeast only). Whereas a disproportionate number of core genes are highly expressed in vegetatively growing cells, many pombe-specific genes are expressed at lower levels. This bias is less pronounced in cells undergoing sexual development, when many pombe-specific genes become highly expressed. This implies that organism-specific proteins are more likely to function during specialized processes such as cellular differentiation. Accordingly, pombe-specific genes were overrepresented among genes induced during sexual development; they were particularly enriched in a group of genes induced during meiotic prophase, when homologous chromosomes pair and recombine. This raises the possibility that organism-specific genes with functions in meiotic prophase favor speciation by preventing fruitful meiosis between closely related organisms. Finally, the set of genes induced late during sexual differentiation, at the time of spore formation, was enriched in yeast-specific genes, indicating that these genes play specialized roles in ascospore development.

What's in the genome of a filamentous fungus? Analysis of the Neurospora genome sequence

The German Neurospora Genome Project has assembled sequences from ordered cosmid and BAC clones of linkage groups II and V of the genome of Neurospora crassa in 13 and 12 contigs, respectively. Including additional sequences located on other linkage groups a total of 12 Mb were subjected to a manual gene extraction and annotation process. The genome comprises a small number of repetitive elements, a low degree of segmental duplications and very few paralogous genes. The analysis of the 3218 identified open reading frames provides a first overview of the protein equipment of a filamentous fungus. Significantly, N.crassa possesses a large variety of metabolic enzymes including a substantial number of enzymes involved in the degradation of complex substrates as well as secondary metabolism. While several of these enzymes are specific for filamentous fungi many are shared exclusively with prokaryotes.

Transcriptome characterization of the dimorphic and pathogenic fungus Paracoccidioides brasiliensis by EST analysis

Paracoccidioides brasiliensis is a pathogenic fungus that undergoes a temperature-dependent cell morphology change from mycelium (22 degrees C) to yeast (36 degrees C). It is assumed that this morphological transition correlates with the infection of the human host. Our goal was to identify genes expressed in the mycelium (M) and yeast (Y) forms by EST sequencing in order to generate a partial map of the fungus transcriptome. Individual EST sequences were clustered by the CAP3 program and annotated using Blastx similarity analysis and InterPro Scan. Three different databases, GenBank nr, COG (clusters of orthologous groups) and GO (gene ontology) were used for annotation. A total of 3,938 (Y = 1,654 and M = 2,274) ESTs were sequenced and clustered into 597 contigs and 1,563 singlets, making up a total of 2,160 genes, which possibly represent one-quarter of the complete gene repertoire in P. brasiliensis. From this total, 1,040 were successfully annotated and 894 could be classified in 18 functional COG categories as follows: cellular metabolism (44%); information storage and processing (25%); cellular processes-cell division, posttranslational modifications, among others (19%); and genes of unknown functions (12%). Computer analysis enabled us to identify some genes potentially involved in the dimorphic transition and drug resistance. Furthermore, computer subtraction analysis revealed several genes possibly expressed in stage-specific forms of P. brasiliensis. Further analysis of these genes may provide new insights into the pathology and differentiation of P. brasiliensis.

A genomics approach to crop pest and disease research

Genome-wide analyses of gene function and gene expression are beginning to yield valuable information in many areas of biological research, and these genomic tools are now being applied to crop pest and disease research. DNA sequencing of cDNA libraries to generate sets of expressed sequence tags (ESTs) are allowing gene compendiums for crop diseases to be compiled. Annotation of such data collections is also providing a wealth of functional information about gene products through similarities to proteins with known function. The next phase of the functional genomics era will be to employ large-scale techniques to knock out or silence genes in order to synthesize gene-specific mutants for phenotypic analysis and to use micro-array methodology to analyze global gene expression, protein turnover and protein processing during the processes of parasitism and colonization. Application of these technologies promises to accelerate the pace that biological information relevant to crop protection accrues. The ability of researchers to assimilate this information into complex models and workable hypotheses is, thus, set to revolutionize the way we study pests and diseases of crop plants.

Deducing the origin of soluble adenylyl cyclase, a gene lost in multiple lineages

The family of eukaryotic adenylyl cyclases consists of a very large group of 12 transmembrane adenylyl cyclases and a very small group of soluble adenylyl cyclase (sAC). Orthologs of human sAC are present in rat Dictyostelium and bacteria but absent from the completely sequenced genomes of Drosophila melanogaster, Caenorhabditis elegans, Arabidopsis thaliana, and Saccharomyces cereviciae. sAC consists of two cyclase domains and a long approximately 1,000 amino acid C-terminal (sCKH) region. This sCKH region and one cyclase domain have been found in only four bacterial genes; the sCKH region was also detected in bacterial Lux-transcription factors and in complex bacterial and fungal kinases. The phylogenies of the kinase and cyclase domains are identical to the phylogeny of the corresponding sCKH domain, suggesting that the sCKH region fused with the other domains early during evolution in bacteria. The amino acid sequences of sAC proteins yield divergence times from the human lineage for rat and Dictyostelium that are close to the reported divergence times of many other proteins in these species. The combined results suggest that the sCKH region was fused with one cyclase domain in bacteria, and a second cyclase domain was added in bacteria or early eukaryotes. The sAC was retained in a few bacteria and throughout the entire evolution of the human lineage but lost independently from many bacteria and from the lineages of plants, yeast, worms, and flies. We conclude that within the family of adenylyl cyclases, soluble AC was poorly fixed during evolution, whereas membrane-bound AC has expanded to form the subgroups of prevailing adenylyl and guanylyl cyclases.

The role of lineage-specific gene family expansion in the evolution of eukaryotes

A computational procedure was developed for systematic detection of lineage-specific expansions (LSEs) of protein families in sequenced genomes and applied to obtain a census of LSEs in five eukaryotic species, the yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe, the nematode Caenorhabditis elegans, the fruit fly Drosophila melanogaster, and the green plant Arabidopsis thaliana. A significant fraction of the proteins encoded in each of these genomes, up to 80% in A. thaliana, belong to LSEs. Many paralogous gene families in each of the analyzed species are almost entirely comprised of LSEs, indicating that their diversification occurred after the divergence of the major lineages of the eukaryotic crown group. The LSEs show readily discernible patterns of protein functions. The functional categories most prone to LSE are structural proteins, enzymes involved in an organism's response to pathogens and environmental stress, and various components of signaling pathways responsible for specificity, including ubiquitin ligase E3 subunits and transcription factors. The functions of several previously uncharacterized, vastly expanded protein families were predicted through in-depth protein sequence analysis, for example, small-molecule kinases and methylases that are expanded independently in the fly and in the nematode. The functions of several other major LSEs remain mysterious; these protein families are attractive targets for experimental discovery of novel, lineage-specific functions in eukaryotes. LSEs seem to be one of the principal means of adaptation and one of the most important sources of organizational and regulatory diversity in crown-group eukaryotes.

Microevolutionary genomics of bacteria

The availability of multiple complete genome sequences from the same species can facilitate attempts to systematically address basic questions in genome evolution. We refer to such efforts as "microevolutionary genomics". We report the results of comparative analyses of complete intraspecific genome (and proteome) sequences from four bacterial species--Chlamydophila pneumoniae, Escherichia coli, Helicobacter pylori and Neisseria meningitidis. Comparisons of average synonymous (K(s)) and nonsynonymous (K(a)) substitution rates were used to assess the influence of various biological factors on the rate of protein evolution. For example, E. coli experiences the most intense purifying selection of the species analyzed, and this may be due to the relatively larger population size of this species. In addition, essential genes were shown to be more evolutionarily conserved than nonessential genes in E. coli and duplicated genes have higher rates of evolution than unique genes for all species studied except C. pneumoniae. Different functional categories of genes were shown to evolve at significantly different rates emphasizing the role of category-specific functional constraints in determining evolutionary rates. Finally, functionally characterized genes tend to be conserved between strains, while uncharacterized genes are over-represented among the unique, strain-specific genes. This suggests the possibility that nonessential genes are responsible for driving the evolutionary diversification between strains.

Predicting the distribution, conservation, and functions of SNAREs and related proteins in fungi

Hyphal tip growth, the hallmark of the fungi, requires highly polarized and localized exocytosis, but how this requirement is met is unknown. Members of conserved protein families called SNAREs and Rabs mediate vesicle trafficking and fusion at virtually every step of the intracellular pathway in all examined eukaryotes. We have searched the available nearly complete fungal genomes, established the presence or absence of members of the SNARE and Rab families in these genomes, and predicted their evolutionary relationships to one another. Comparisons with the extensively studied Saccharomyces cerevisiae indicate that, in general, most of the members of these families (including those involved in mediating exocytosis) are conserved. The presence of exceptional SNAREs and Rabs in some fungi that are not conserved in S. cerevisiae may be indicative of specialized steps that occur in these fungi. The implications of these findings for current tip growth models are discussed.

In silico sequence analysis of arylamine N-acetyltransferases: evidence for an absence of lateral gene transfer from bacteria to vertebrates and first description of paralogs in bacteria

The arylamine N-acetyltransferases (NATs) are xenobiotic-metabolizing enzymes responsible for the biotransformation of various arylamine and heterocyclic amines, including drugs and carcinogenic compounds. NAT and NAT-like genes have been identified in several vertebrate and eubacterial species. Little is known about their evolutionary history, but the horizontal transfer of NAT genes from bacteria to vertebrates was recently suggested [S. Salzberg, O. White, J. Peterson, J. Eisen, Science 292 (2001) 1903]. We used various bioinformatics-based approaches to screen eukaryotic and prokaryotic genomes. We identified Mesorhizobium loti NAT genes as the first examples of NAT paralogs in prokaryotes. As shown for vertebrate species, the existence of NAT paralogs in this bacterium may be accounted for by enzymatic specialization after gene duplication. Phylogenetic analysis following the identification of a NAT ortholog in the nonvertebrate species Ciona intestinalis indicated that NAT genes are unlikely to be examples of direct horizontal gene transfer (HGT). Our study suggests that NAT genes have evolved from a common ancestor, with a succession of nonvertebrate intermediates. The absence of NAT genes in yeast, nematode worms, fruit flies, and mustard weed may result from gene loss in these nonvertebrate lineages. These results provide new insight into the taxonomic distribution and evolutionary history of this class of drug-metabolizing enzymes.

Timeline: Neurospora: a model of model microbes

In the 1940s, studies with Neurospora pioneered the use of microorganisms in genetic analysis and provided the foundations for biochemical genetics and molecular biology. What has happened since this orange mould was used to show that genes control metabolic reactions? How did it come to be the fungal counterpart of Drosophila? We describe its continued use during the heyday of research with Escherichia coli and yeast, and its emergence as a biological model for higher fungi.

Comparative genomics and evolution of proteins involved in RNA metabolism

RNA metabolism, broadly defined as the compendium of all processes that involve RNA, including transcription, processing and modification of transcripts, translation, RNA degradation and its regulation, is the central and most evolutionarily conserved part of cell physiology. A comprehensive, genome-wide census of all enzymatic and non-enzymatic protein domains involved in RNA metabolism was conducted by using sequence profile analysis and structural comparisons. Proteins related to RNA metabolism comprise from 3 to 11% of the complete protein repertoire in bacteria, archaea and eukaryotes, with the greatest fraction seen in parasitic bacteria with small genomes. Approximately one-half of protein domains involved in RNA metabolism are present in most, if not all, species from all three primary kingdoms and are traceable to the last universal common ancestor (LUCA). The principal features of LUCA's RNA metabolism system were reconstructed by parsimony-based evolutionary analysis of all relevant groups of orthologous proteins. This reconstruction shows that LUCA possessed not only the basal translation system, but also the principal forms of RNA modification, such as methylation, pseudouridylation and thiouridylation, as well as simple mechanisms for polyadenylation and RNA degradation. Some of these ancient domains form paralogous groups whose evolution can be traced back in time beyond LUCA, towards low-specificity proteins, which probably functioned as cofactors for ribozymes within the RNA world framework. The main lineage-specific innovations of RNA metabolism systems were identified. The most notable phase of innovation in RNA metabolism coincides with the advent of eukaryotes and was brought about by the merge of the archaeal and bacterial systems via mitochondrial endosymbiosis, but also involved emergence of several new, eukaryote-specific RNA-binding domains. Subsequent, vast expansions of these domains mark the origin of alternative splicing in animals and probably in plants. In addition to the reconstruction of the evolutionary history of RNA metabolism, this analysis produced numerous functional predictions, e.g. of previously undetected enzymes of RNA modification.

Large scale analysis of sequences from Neurospora crassa

After 50 years of analysing Neurospora crassa genes one by one large scale sequence analysis has increased the number of accessible genes tremendously in the last few years. Being the only filamentous fungus for which a comprehensive genomic sequence database is publicly accessible N. crassa serves as the model for this important group of microorganisms. The MIPS N. crassa database currently holds more than 16 Mb of non-redundant data of the chromosomes II and V analysed by the German Neurospora Genome Project. This represents more than one-third of the genome. Open reading frames (ORFs) have been extracted from the sequence and the deduced proteins have been annotated extensively. They are classified according to matches in sequence databases and attributed to functional categories according to their relatives. While 41% of analysed proteins are related to known proteins, 30% are hypothetical proteins with no match to a database entry. The entire genome is expected to comprise some 13000 protein coding genes, more than twice as many as found in yeasts, and reflects the high potential of filamentous fungi to cope with various environmental conditions.

Genome and protein evolution in eukaryotes

The past year has seen the completion of the genome sequence of the flowering plant Arabidopsis thaliana and the initial sequence reports of the human genome. The availability of completely sequenced eukaryotic genomes from disparate phylogenetic lineages has opened the door to comparative analyses and a better understanding of the evolutionary processes shaping genomes. Complex many-to-many relationships between genes from different species appear to be the norm, suggesting that transfer of detailed functional annotation will not be straightforward. In addition to expansion and contraction of gene families, new genes evolve from recombination of pre-existing domains, although some domain families do appear to have evolved recently and to be specific to restricted phylogenetic lineages. The overall picture is of a huge diversity of gene content within eukaryotic genomes, reflecting different functional demands in different species.

On the species of origin: diagnosing the source of symbiotic transcripts

BACKGROUND

Most organisms have developed ways to recognize and interact with other species. Symbiotic interactions range from pathogenic to mutualistic. Some molecular mechanisms of interspecific interaction are well understood, but many remain to be discovered. Expressed sequence tags (ESTs) from cultures of interacting symbionts can help identify transcripts that regulate symbiosis, but present a unique challenge for functional analysis. Given a sequence expressed in an interaction between two symbionts, the challenge is to determine from which organism the transcript originated. For high-throughput sequencing from interaction cultures, a reliable computational approach is needed. Previous investigations into GC nucleotide content and comparative similarity searching provide provisional solutions, but a comparative lexical analysis, which uses a likelihood-ratio test of hexamer counts, is more powerful.

RESULTS

Validation with genes whose origin and function are known yielded 94% accuracy. Microbial (non-plant) transcripts comprised 75% of a Phytophthora sojae-infected soybean (Glycine max cv Harasoy) library, contrasted with 15% or less in root tissue libraries of Medicago truncatula from axenic, Phytophthora medicaginis-infected, mycorrhizal, and rhizobacterial treatments. Mycorrhizal libraries contained about 23% microbial transcripts; an axenic plant library contained a similar proportion of putative microbial transcripts.

CONCLUSIONS

Comparative lexical analysis offers numerous advantages over alternative approaches. Many of the transcripts isolated from mixed cultures were of unknown function, suggesting specificity to symbiotic metabolism and therefore candidates likely to be interesting for further functional investigation. Future investigations will determine whether the abundance of non-plant transcripts in a pure plant library indicates procedural artifacts, horizontally transferred genes, or other phenomena.

Fungal genomics and pathogenicity

The filamentous fungal genetics community has enthusiastically embraced the utilization of genomics technologies to resolve long-standing issues in fungal biology. For example, such technologies have been proposed to study the mechanics of tip growth, photoreception, gene silencing, the molecular basis of conidiation, the pathway leading to sexual reproduction, and mechanisms of pathogenesis. These studies have provided a refreshing change of pace in research on filamentous fungi, which has lagged behind that on other eukaryotes in the exploitation of genome-wide methodologies. Despite the late start, several fungal genome sequencing projects are underway. The resulting databases will allow the comprehensive analysis of developmental processes that are characteristic of fungi, including the molecular nature of pathogenicity. DNA databases underpin analyses of the fungal transcriptome, proteome, and metabolome. This combined information will contribute to our basic understanding of not only the mechanics of infection but also the evolution of pathogenicity.

Fungal Zuotin proteins evolved from MIDA1-like factors by lineage-specific loss of MYB domains

Proteins are often characterized by the presence of multiple domains, which make specific contributions to their cellular function. While the gain of domains in proteins by duplication and shuffling is well established, domain loss is poorly documented. Here, we provide evidence that domain loss has played an important role in the evolution of protein architecture and function by demonstrating that fungal Zuotin proteins evolved from MIDA1-like proteins, present in animals and plants, by complete loss of the carboxyl-terminal MYB domains. Phylogenetic analyses of the DnaJ motif (the J domain) present in both Zuotin and MIDA1 proteins were complicated by the limited length and profound differences in evolutionary rates exhibited by this domain. To rigorously examine J domain phylogeny, we combined the nonparametric bootstrap with Monte Carlo simulation. This method, which we have designated the resampled parametric bootstrap, allowed us to assess type I and type II error associated with these analyses. These results revealed significant support for domain loss rather than domain gain or gene loss involving paralogs. The absence of sequences related to the MIDA1 MYB domains in Saccharomyces cerevisiae further indicates that the domains have been completely lost, consistent with known functional differences between Zuotin and MIDA1 proteins. These analyses suggest that the description of additional examples of complete domain loss may provide a method to identify orthologous proteins exhibiting functional differences using genomic sequence data.

Genetic and molecular analysis of circadian rhythms in Neurospora

Over the course of the past 40 years Neurospora has become a well-known and uniquely tractable model system for the analysis of the molecular basis of eukaryotic circadian oscillatory systems. Molecular bases for the period length and sustainability of the rhythm, light, and temperature resetting of the circadian system and for gating of light input and light effects are becoming understood, and Neurospora promises to be a suitable system for examining the role of coupled feedback loops in the clock. Many of these insights have shown or foreshadow direct parallels in mammalian systems, including the mechanism of light entrainment, the involvement of PAS:PAS heterodimers as transcriptional activators in essential clock-associated feedback loops, and dual role of FRQ in the loop as an activator and a repressor; similarities extend to the primary sequence level in at least one case, that of WC-1 and BMAL1. Work on circadian output in Neurospora has identified more than a dozen regulated genes and has been at the forefront of studies aimed at understanding clock control of gene expression.

Microbial genes in the human genome: lateral transfer or gene loss?

The human genome was analyzed for evidence that genes had been laterally transferred into the genome from prokaryotic organisms. Protein sequence comparisons of the proteomes of human, fruit fly, nematode worm, yeast, mustard weed, eukaryotic parasites, and all completed prokaryote genomes were performed, and all genes shared between human and each of the other groups of organisms were collected. About 40 genes were found to be exclusively shared by humans and bacteria and are candidate examples of horizontal transfer from bacteria to vertebrates. Gene loss combined with sample size effects and evolutionary rate variation provide an alternative, more biologically plausible explanation.

Analysis of the pdx-1 (snz-1/sno-1) region of the Neurospora crassa genome: correlation of pyridoxine-requiring phenotypes with mutations in two structural genes

We report the analysis of a 36-kbp region of the Neurospora crassa genome, which contains homologs of two closely linked stationary phase genes, SNZ1 and SNO1, from Saccharomyces cerevisiae. Homologs of SNZ1 encode extremely highly conserved proteins that have been implicated in pyridoxine (vitamin B6) metabolism in the filamentous fungi Cercospora nicotianae and in Aspergillus nidulans. In N. crassa, SNZ and SNO homologs map to the region occupied by pdx-1 (pyridoxine requiring), a gene that has been known for several decades, but which was not sequenced previously. In this study, pyridoxine-requiring mutants of N. crassa were found to possess mutations that disrupt conserved regions in either the SNZ or SNO homolog. Previously, nearly all of these mutants were classified as pdx-1. However, one mutant with a disrupted SNO homolog was at one time designated pdx-2. It now appears appropriate to reserve the pdx-1 designation for the N. crassa SNZ homolog and pdx-2 for the SNO homolog. We further report annotation of the entire 36,030-bp region, which contains at least 12 protein coding genes, supporting a previous conclusion of high gene densities (12,000-13,000 total genes) for N. crassa. Among genes in this region other than SNZ and SNO homologs, there was no evidence of shared function. Four of the genes in this region appear to have been lost from the S. cerevisiae lineage.

Genome comparisons highlight similarity and diversity within the eukaryotic kingdoms

In 2000, the number of completely sequenced eukaryotic genomes increased to four. The addition of Drosophila and Arabidopsis into this cohort permits additional insights into the processes that have shaped evolution. Analysis and comparisons of both completed genomes and partially sequenced genomes have already shed light on mechanisms such as gene duplication and gene loss that have long been hypothesized to be major forces in speciation. Indeed, duplicate gene pairs in Saccharomyces, Arabidopsis, Caenorhabditis and Drosophila are high: 30%, 60%, 48% and 40%, respectively. Evidence of horizontal gene-transfer, thought to be a major evolutionary force in bacteria, has been found in Arabidopsis. The release of the 'first draft' of the human genome sequence in 2000 heralds a new stage of biological study. Understanding the as-yet-unannotated human genome will be largely based on conclusions, techniques and tools developed during the analysis and comparison of the genome of these four model organisms.