Comparison of codon usage and tRNAs in mitochondrial genomes of Candida species

To gain insight into the nature of the mitochondrial genomes (mtDNA) of different Candida species, the synonymous codon usage bias of mitochondrial protein coding genes and the tRNAs in C. albicans, C. parapsilosis, C. stellata, C. glabrata and the closely related yeast Saccharomyces cerevisiae were analyzed. Common features of the mtDNA in Candida species are a strong A+T pressure on protein coding genes, and insufficient mitochondrial tRNA species are encoded to perform protein synthesis. The wobble site of the anticodon is always U for the NNR (NNA and NNG) codon families, which are dominated by A-ending codons, and always G for the NNY (NNC and NNU) codon families, which is dominated by U-ending codons, and always U for the NNN (NNA, NNU, NNC and NNG) codon families, which are dominated by A-ending codons and U-ending codons. Patterns of synonymous codon usage of Candida species can be classified into three groups: (1) optimal codon-anticodon usage, Glu, Lys, Leu (translated by anti-codon UAA), Gln, Arg (translated by anti-codon UCU) and Trp are containing NNR codons. NNA, whose corresponding tRNA is encoded in the mtDNA, is used preferentially. (2) Non-optimal codon-anticodon usage, Cys, Asp, Phe, His, Asn, Ser (translated by anti-codon GCU) and Tyr are containing NNY codons. The NNU codon, whose corresponding tRNA is not encoded in the mtDNA, is used preferentially. (3) Combined codon-anticodon usage, Ala, Gly, Leu (translated by anti-codon UAG), Pro, Ser (translated by anti-codon UGA), Thr and Val are containing NNN codons. NNA (tRNA encoded in the mtDNA) and NNU (tRNA not encoded in the mtDNA) are used preferentially. In conclusion, we propose that in Candida species, codons containing A or U at third position are used preferentially, regardless of whether corresponding tRNAs are encoded in the mtDNA. These results might be useful in understanding the common features of the mtDNA in Candida species and patterns of synonymous codon usage.

A genomic study of the inter-ORF distances in Saccharomyces cerevisiae

The genome of eukaryotic microbes is usually quite compacted. The yeast Saccharomyces cerevisiae is one of the best-known examples. Open reading frames (ORFs) occupy about 75% of the total DNA sequence. The existence of other, non-protein coding genes and other genetic elements leaves very little space for gene promoters and terminators. We have performed an in silico study of inter-ORF distances that shows that there is a minimum distance between two adjacent ORFs that depends on the relative orientation between them. Our analyses suggest that different kinds of promoters and terminators exist with regard to their length and ability to overlap each other. The experimental testing of some putative exceptions to the minimum length model in tandemly orientated ORF pairs suggests that, in those cases, defects in promoter or terminator functionality exist that provoke transcription of polycistronic mRNAs.

Genomic analysis of the hierarchical structure of regulatory networks

A fundamental question in biology is how the cell uses transcription factors (TFs) to coordinate the expression of thousands of genes in response to various stimuli. The relationships between TFs and their target genes can be modeled in terms of directed regulatory networks. These relationships, in turn, can be readily compared with commonplace "chain-of-command" structures in social networks, which have characteristic hierarchical layouts. Here, we develop algorithms for identifying generalized hierarchies (allowing for various loop structures) and use these approaches to illuminate extensive pyramid-shaped hierarchical structures existing in the regulatory networks of representative prokaryotes (Escherichia coli) and eukaryotes (Saccharomyces cerevisiae), with most TFs at the bottom levels and only a few master TFs on top. These masters are situated near the center of the protein-protein interaction network, a different type of network from the regulatory one, and they receive most of the input for the whole regulatory hierarchy through protein interactions. Moreover, they have maximal influence over other genes, in terms of affecting expression-level changes. Surprisingly, however, TFs at the bottom of the regulatory hierarchy are more essential to the viability of the cell. Finally, one might think master TFs achieve their wide influence through directly regulating many targets, but TFs with most direct targets are in the middle of the hierarchy. We find, in fact, that these midlevel TFs are "control bottlenecks" in the hierarchy, and this great degree of control for "middle managers" has parallels in efficient social structures in various corporate and governmental settings.

Epistasis correlates to genomic complexity

Whether systematic genetic interactions (epistasis) occur at the genomic scale remains a challenging topic in evolutionary biology. Epistasis should make a significant contribution to variation in complex traits and influence the evolution of genetic systems as sex, diploidy, dominance, or the contamination of genomes with deleterious mutations. We have collected data from widely different organisms and quantified epistasis in a common, per-generation scale. Simpler genomes, such as those of RNA viruses, display antagonistic epistasis (mutations have smaller effects together than expected); bacterial microorganisms do not apparently deviate from independent effects, whereas in multicellular eukaryotes, a transition toward synergistic epistasis occurs (mutations have larger effects together than expected). We propose that antagonistic epistasis might be a property of compact genomes with few nonpleiotropic biological functions, whereas in complex genomes, synergism might emerge from mutational robustness.

Identification and comparative analysis of sixteen fungal peptidyl-prolyl cis/trans isomerase repertoires

BACKGROUND

The peptidyl-prolyl cis/trans isomerase (PPIase) class of proteins is present in all known eukaryotes, prokaryotes, and archaea, and it is comprised of three member families that share the ability to catalyze the cis/trans isomerisation of a prolyl bond. Some fungi have been used as model systems to investigate the role of PPIases within the cell, however how representative these repertoires are of other fungi or humans has not been fully investigated.

RESULTS

PPIase numbers within these fungal repertoires appears associated with genome size and orthology between repertoires was found to be low. Phylogenetic analysis showed the single-domain FKBPs to evolve prior to the multi-domain FKBPs, whereas the multi-domain cyclophilins appear to evolve throughout cyclophilin evolution. A comparison of their known functions has identified, besides a common role within protein folding, multiple roles for the cyclophilins within pre-mRNA splicing and cellular signalling, and within transcription and cell cycle regulation for the parvulins. However, no such commonality was found with the FKBPs. Twelve of the 17 human cyclophilins and both human parvulins, but only one of the 13 human FKBPs, identified orthologues within these fungi. hPar14 orthologues were restricted to the Pezizomycotina fungi, and R. oryzae is unique in the known fungi in possessing an hCyp33 orthologue and a TPR-containing FKBP. The repertoires of Cryptococcus neoformans, Aspergillus fumigatus, and Aspergillus nidulans were found to exhibit the highest orthology to the human repertoire, and Saccharomyces cerevisiae one of the lowest.

CONCLUSION

Given this data, we would hypothesize that: (i) the evolution of the fungal PPIases is driven, at least in part, by the size of the proteome, (ii) evolutionary pressures differ both between the different PPIase families and the different fungi, and (iii) whilst the cyclophilins and parvulins have evolved to perform conserved functions, the FKBPs have evolved to perform more variable roles. Also, the repertoire of Cryptococcus neoformans may represent a better model fungal system within which to study the functions of the PPIases as its genome size and genetic tractability are equal to those of Saccharomyces cerevisiae, whilst its repertoires exhibits greater orthology to that of humans. However, further experimental investigations are required to confirm this.

A model for the evolution of paralog families in genomes

We introduce and analyse a simple probabilistic model of genome evolution. It is based on three fundamental evolutionary events: gene loss, duplication and accumulated change. This is motivated by previous works which consisted in fitting the available genomic data into, what is called paralog distributions. This formalism is described by a system of infinite number of linear equations. We show that this system generates a semigroup of linear operators on the space l (1). We prove that size distribution of paralogous gene families in a genome converges to the equilibrium as time goes to infinity. Moreover we show that when probabilities of gene removal and duplication are close to each other, then the resulting distribution is close to logarithmic distribution. Some empirical results for yeast genomes are presented.

Isw1 functions in parallel with the NuA4 and Swr1 complexes in stress-induced gene repression

The packaging of DNA into chromatin allows eukaryotic cells to organize and compact their genomes but also creates an environment that is generally repressive to nuclear processes that depend upon DNA accessibility. There are several classes of enzymes that modulate the primary structure of chromatin to regulate various DNA-dependent processes. The biochemical activities of the yeast Isw1 ATP-dependent chromatin-remodeling enzyme have been well characterized in vitro, but little is known about how these activities are utilized in vivo. In this work, we sought to discern genetic backgrounds that require Isw1 activity for normal growth. We identified a three-way genetic interaction among Isw1, the NuA4 histone acetyltransferase complex, and the Swr1 histone replacement complex. Transcription microarray analysis revealed parallel functions for these three chromatin-modifying factors in the regulation of TATA-containing genes, including the repression of a large number of stress-induced genes under normal growth conditions. In contrast to a recruitment-based model, we find that the NuA4 and Swr1 complexes act throughout the genome while only a specific subset of the genome shows alterations in transcription.

Identification and functional characterization of hCLS1, a human cardiolipin synthase localized in mitochondria

In eukaryotic cells, CLS (cardiolipin synthase) is involved in the final step of cardiolipin synthesis by catalysing the transfer of a phosphatidyl residue from CDP-DAG (diacylglycerol) to PG (phosphatidylglycerol). Despite an important role of cardiolipin in regulating mitochondrial function, a gene encoding the mammalian CLS has not been identified so far. We report in the present study the identification and characterization of a human cDNA encoding the first mammalian CLS [hCLS1 (human CLS1)]. The predicted hCLS1 peptide sequence shares significant homology with the yeast and plant CLS proteins. The recombinant hCLS1 enzyme expressed in COS-7 cells catalysed efficiently the synthesis of cardiolipin in vitro using CDP-DAG and PG as substrates. Furthermore, overexpression of hCLS1 cDNA in COS-7 cells resulted in a significant increase in cardiolipin synthesis in intact COS-7 cells without any significant effects on the activity of the endogenous phosphatidylglycerophosphate synthase of the transfected COS-7 cells. Immunohistochemical analysis demonstrated that the recombinant hCLS1 protein was localized to the mitochondria when transiently expressed in COS-7 cells, which was further corroborated by results from subcellular fractionation analyses of the recombinant hCLS1 protein. Northern-blot analysis showed that the hCLS1 gene was predominantly expressed in tissues that require high levels of mitochondrial activities for energy metabolism, with the highest expression in skeletal and cardiac muscles. High levels of hCLS1 expression were also detected in liver, pancreas, kidney and small intestine, implying a functional role of hCLS1 in these tissues.

Saccharomyces cerevisiae cells have three Omega class glutathione S-transferases acting as 1-Cys thiol transferases

The Saccharomyces cerevisiae genome encodes three proteins that display similarities with human GSTOs (Omega class glutathione S-transferases) hGSTO1-1 and hGSTO2-2. The three yeast proteins have been named Gto1, Gto2 and Gto3, and their purified recombinant forms are active as thiol transferases (glutaredoxins) against HED (beta-hydroxyethyl disulphide), as dehydroascorbate reductases and as dimethylarsinic acid reductases, while they are not active against the standard GST substrate CDNB (1-chloro-2,4-dinitrobenzene). Their glutaredoxin activity is also detectable in yeast cell extracts. The enzyme activity characteristics of the Gto proteins contrast with those of another yeast GST, Gtt1. The latter is active against CDNB and also displays glutathione peroxidase activity against organic hydroperoxides such as cumene hydroperoxide, but is not active as a thiol transferase. Analysis of point mutants derived from wild-type Gto2 indicates that, among the three cysteine residues of the molecule, only the residue at position 46 is required for the glutaredoxin activity. This indicates that the thiol transferase acts through a monothiol mechanism. Replacing the active site of the yeast monothiol glutaredoxin Grx5 with the proposed Gto2 active site containing Cys46 allows Grx5 to retain some activity against HED. Therefore the residues adjacent to the respective active cysteine residues in Gto2 and Grx5 are important determinants for the thiol transferase activity against small disulphide-containing molecules.

TOS9 regulates white-opaque switching in Candida albicans

In Candida albicans, the a1-alpha2 complex represses white-opaque switching, as well as mating. Based upon the assumption that the a1-alpha2 corepressor complex binds to the gene that regulates white-opaque switching, a chromatinimmunoprecipitation-microarray analysis strategy was used to identify 52 genes that bound to the complex. One of these genes, TOS9, exhibited an expression pattern consistent with a "master switch gene." TOS9 was only expressed in opaque cells, and its gene product, Tos9p, localized to the nucleus. Deletion of the gene blocked cells in the white phase, misexpression in the white phase caused stable mass conversion of cells to the opaque state, and misexpression blocked temperature-induced mass conversion from the opaque state to the white state. A model was developed for the regulation of spontaneous switching between the opaque state and the white state that includes stochastic changes of Tos9p levels above and below a threshold that induce changes in the chromatin state of an as-yet-unidentified switching locus. TOS9 has also been referred to as EAP2 and WOR1.

Inferring direct regulatory targets from expression and genome location analyses: a comparison of transcription factor deletion and overexpression

BACKGROUND

Effects on gene expression due to environmental or genetic changes can be easily measured using microarrays. However, indirect effects on expression can be substantial. The indirect effects of a perturbation need to be distinguished from the direct effects if we are to understand the structure and behavior of regulatory networks.

RESULTS

The most direct way to perturb a transcriptional network is to alter transcription factor activity. Here, for the first time, we compare expression changes and genomic binding in a simple regulon under conditions of both low and high transcription factor activity. Specifically, we assessed the effects on expression and binding due to deletion of the yeast LEU3 transcription factor gene and effects due to elevation of Leu3 activity. Leu3 activity was elevated through overexpression and the introduction of a mutation that renders the protein constitutively active. Genes that are bound and/or regulated by Leu3 under one or both conditions were characterized in terms of their functional annotations and their predicted potential to be bound by Leu3. We also assessed the evolutionary conservation of the predicted binding potential using a novel alignment-independent method. Both perturbations yield genes that are likely to be direct targets of Leu3, including most of the classically defined targets. Additional direct targets are identified by each of the methods. However, experimental and computational criteria suggest that most genes whose expression is affected by the Leu3 genotype are unlikely to be regulated by binding of the protein.

CONCLUSION

Most genes that are differentially expressed by Leu3 are not direct targets despite the exceptional simplicity of the regulon, and the unusually direct nature of the perturbations investigated. These conclusions are reached through computational analyses that support and extend chromatin immunoprecipitation data on the identities of direct targets. These results have implications for the interpretation of expression experiments, especially in cases for which chromatin immunoprecipitation data are unavailable, incomplete, or ambiguous.

Genome-wide patterns of histone modifications in yeast

Post-translational histone modifications and histone variants generate complexity in chromatin to enable the many functions of the chromosome. Recent studies have mapped histone modifications across the Saccharomyces cerevisiae genome. These experiments describe how combinations of modified and unmodified states relate to each other and particularly to chromosomal landmarks that include heterochromatin, subtelomeric chromatin, centromeres, origins of replication, promoters and coding regions. Such patterns might be important for the regulation of heterochromatin-mediated silencing, chromosome segregation, DNA replication and gene expression.

Asymmetric positioning of nucleosomes and directional establishment of transcriptionally silent chromatin by Saccharomyces cerevisiae silencers

In Saccharomyces cerevisiae, silencers flanking the HML and HMR loci consist of various combinations of binding sites for Abf1p, Rap1p, and the origin recognition complex (ORC) that serve to recruit the Sir silencing complex, thereby initiating the establishment of transcriptionally silent chromatin. There have been seemingly conflicting reports concerning whether silencers function in an orientation-dependent or -independent manner, and what determines the directionality of a silencer has not been explored. We demonstrate that chromatin plays a key role in determining the potency and directionality of silencers. We show that nucleosomes are asymmetrically distributed around the HML-I or HMR-E silencer so that a nucleosome is positioned close to the Abf1p side but not the ORC side of the silencer. This coincides with preferential association of Sir proteins and transcriptional silencing on the Abf1p side of the silencer. Elimination of the asymmetry in nucleosome positioning at a silencer leads to comparable silencing on both sides. Asymmetric nucleosome positioning in the immediate vicinity of a silencer is independent of its orientation and genomic context, indicating that it is the inherent property of the silencer. Moreover, it is also independent of the Sir complex and thus precedes the formation of silent chromatin. Finally, we demonstrate that asymmetric positioning of nucleosomes and directional silencing by a silencer depend on ORC and Abf1p. We conclude that the HML-I and HMR-E silencers promote asymmetric positioning of nucleosomes, leading to unequal potentials of transcriptional silencing on their sides and, hence, directional silencing.

Discovery of aspoquinolones A-D, prenylated quinoline-2-one alkaloids from Aspergillus nidulans, motivated by genome mining

Motivated by the observation that the Aspergillus nidulans genome bears multiple anthranilic acid synthase gene copies, the fungal metabolome was reinvestigated under various fermentation conditions, resulting in the discovery of novel prenylated quinolin-2-one alkaloids, two of which bear unprecedented terpenoid side chains.

Generation, annotation and analysis of ESTs from Trichoderma harzianum CECT 2413

BACKGROUND

The filamentous fungus Trichoderma harzianum is used as biological control agent of several plant-pathogenic fungi. In order to study the genome of this fungus, a functional genomics project called "TrichoEST" was developed to give insights into genes involved in biological control activities using an approach based on the generation of expressed sequence tags (ESTs).

RESULTS

Eight different cDNA libraries from T. harzianum strain CECT 2413 were constructed. Different growth conditions involving mainly different nutrient conditions and/or stresses were used. We here present the analysis of the 8,710 ESTs generated. A total of 3,478 unique sequences were identified of which 81.4% had sequence similarity with GenBank entries, using the BLASTX algorithm. Using the Gene Ontology hierarchy, we performed the annotation of 51.1% of the unique sequences and compared its distribution among the gene libraries. Additionally, the InterProScan algorithm was used in order to further characterize the sequences. The identification of the putatively secreted proteins was also carried out. Later, based on the EST abundance, we examined the highly expressed genes and a hydrophobin was identified as the gene expressed at the highest level. We compared our collection of ESTs with the previous collections obtained from Trichoderma species and we also compared our sequence set with different complete eukaryotic genomes from several animals, plants and fungi. Accordingly, the presence of similar sequences in different kingdoms was also studied.

CONCLUSION

This EST collection and its annotation provide a significant resource for basic and applied research on T. harzianum, a fungus with a high biotechnological interest.

A role for TFIIIC transcription factor complex in genome organization

Eukaryotic genome complexity necessitates boundary and insulator elements to partition genomic content into distinct domains. We show that inverted repeat (IR) boundary elements flanking the fission yeast mating-type heterochromatin domain contain B-box sequences, which prevent heterochromatin from spreading into neighboring euchromatic regions by recruiting transcription factor TFIIIC complex without RNA polymerase III (Pol III). Genome-wide analysis reveals TFIIIC with Pol III at all tRNA genes, many of which cluster at pericentromeric heterochromatin domain boundaries. However, a single tRNA(phe) gene with modest TFIIIC enrichment is insufficient to serve as boundary and requires RNAi-associated element to restrain heterochromatin spreading. Remarkably, we found TFIIIC localization without Pol III at many sites located between divergent promoters. These sites appear to act as chromosome-organizing clamps by tethering distant loci to the nuclear periphery, at which TFIIIC is concentrated into several distinct bodies. Our analyses uncover a general genome organization mechanism involving conserved TFIIIC complex.

TFIIIC boxes in the genome

In this issue of Cell, Noma et al. (2006) show that B-boxes and TFIIIC limit the spread of heterochromatin at the silent mat region in the fission yeast genome. Global analysis of TFIIIC distribution revealed dispersed sites of association that coalesce at the nuclear periphery, suggesting that TFIIIC may act as a barrier throughout the genome.

A genomic code for nucleosome positioning

Eukaryotic genomes are packaged into nucleosome particles that occlude the DNA from interacting with most DNA binding proteins. Nucleosomes have higher affinity for particular DNA sequences, reflecting the ability of the sequence to bend sharply, as required by the nucleosome structure. However, it is not known whether these sequence preferences have a significant influence on nucleosome position in vivo, and thus regulate the access of other proteins to DNA. Here we isolated nucleosome-bound sequences at high resolution from yeast and used these sequences in a new computational approach to construct and validate experimentally a nucleosome-DNA interaction model, and to predict the genome-wide organization of nucleosomes. Our results demonstrate that genomes encode an intrinsic nucleosome organization and that this intrinsic organization can explain approximately 50% of the in vivo nucleosome positions. This nucleosome positioning code may facilitate specific chromosome functions including transcription factor binding, transcription initiation, and even remodelling of the nucleosomes themselves.

Is optimal gene order impossible?

Recent evidence suggests that yeast genes encoding proteins that are present in the same protein complex tend to be linked and to be co-expressed. More generally, we found that genes that are close to each other in the protein interaction network tend to be linked more often than expected and are often co-expressed. Unexpectedly, we found that linked genes in network proximity have unusually high recombination rates. Because high recombination rates are associated with high rates of genome re-organization, our findings might explain why the clustering of genes in proximity in the network is such a weak effect: there could be a co-evolutionary cycle of physical linkage for co-expression, upwards modification of the recombination rate and concomitant break-up of a cluster. Under such a model an "optimal" gene order is never stable.

Tools to study molecular mechanisms of Aspergillus pathogenicity

The unique nature of Aspergillus fungi represents a challenge for scrutinizing the attributes that render these saprophytic microorganisms pathogenic or allergenic under certain environmental circumstances. Recent publication of the genomic sequence from an isolate of the major pathogen Aspergillus fumigatus denotes enormous progress in aiming at cellular features and gene products that contribute to its pathogenicity. Latest developments to study virulence-related characteristics comprise profiling techniques, conditional gene inactivation and precise manipulation of the genome by means of gene targeting. Advances in assessing the virulence potential of particular mutant strains in alternative test systems complement these approaches.

A phylogenomic analysis of the Ascomycota

An automated procedure was developed to extract orthologous sequences from fungal genomes and incorporate them into phylogenomic analyses in a timely and efficient manner. This approach involves parsing an all versus all BLASTP search of 17 proteomes and creating a similarity matrix from e-values, which is then used to cluster proteins into related groups by means of a Markov Clustering algorithm. After performing this analysis at different stringency levels, 854 single copy protein clusters, which were ubiquitously distributed in all 17 proteomes, were identified. These clusters were culled to include only those clusters where all proteins had best hits to and received hits from a protein within the same cluster. The final data set included gapless alignments for 781 clusters of orthologous sequences that were concatenated into one super alignment containing 195,664 amino acid characters. Neighbor-joining distance and maximum likelihood analyses resulted in identical topologies and all except one node received 100% bootstrap support. The node supporting Stagonospora nodorum's position received 83% support or higher; it was also the only taxon differentially resolved in the maximum parsimony analyses. All analyses resolved the two derived subphyla Pezizomycotina and Saccharomycotina, and Schizosaccharomyces pombe as an early diverging lineage of the Ascomycota. Importantly, these analyses resolved the Leotiomycetes as the sister group to the Sordariomycetes, a region of the Ascomycota phylogeny that has remained problematic in molecular phylogenetic studies of more limited character sampling. Additional phylogenetic analyses which included orthologous sequences from an unannotated ascomycotan genome (e.g., Coccidioides immitis) and subsets of orthologs with different characteristics supported this topology. Phylogenetic analyses of the 595 orthologs which included C. immitis resulted in an identical topology to the previous 781 ortholog analysis and correctly placed C. immitis in the Eurotiomycetes. This demonstrated the correct identification of orthologs and the ability to incorporate unannotated genomic data into a common phylogenetic analysis.

Ser/Thr-rich domains are associated with genetic variation and morphogenesis inSaccharomyces cerevisiae

Proteins containing regions of amino acid bias are often found in eukaryotes and are associated with particular functional groups. We have carried out a genomic analysis of yeast proteins containing regions with a significant bias of Ser and Thr residues. Our findings reveal that a high number are cell surface proteins or regulatory proteins involved in many aspects of cell differentiation. Furthermore, in Saccharomyces-related species, a highly significant correlation exists between the frequency of Ser-rich regions and DNA repeats, indicating that their generation may rely on similar factors. Cluster analysis shows that Ser/Thr-rich regions, located within the tandem repeats of cell surface proteins, are encoded to an increased frequency by UCU (Ser) and ACU/ACC (Thr), implying that mutational events that generate iterations could involve these codons. Replication slippage is proposed to be a contributing factor, as mounting evidence suggests that repeat generation in cell surface proteins can occur independently of meiosis. To reinforce this argument, we have discovered a premeiotic association between Mre11p, a nuclease involved in DNA repair, and ORFs encoding Ser/Thr-rich regions. Several macromolecules involved in the glycosylation and phosphorylation of proteins require Ser and Thr residues as binding sites. Ser/Thr-rich regions, through polymorphisms, are associated with the evolution of functional sites, particularly in providing motifs for glycosylation and phosphorylation. These results point to a Ser/Thr-biased somatic mutation mechanism that contributes to rapid evolution in yeast.

Genetic evidence for 18S rRNA binding and an Rps19p assembly function of yeast nucleolar protein Nep1p

The nucleolar protein Nep1 and its human homologue were previously shown to be involved in the maturation of 18S rRNA and to interfere directly or indirectly with a methylation reaction. Here, we report that the loss-of-function mutation Deltasnr57 and multicopy expression of the ribosomal 40S subunit protein 19 (Rps19p) can partially suppress the Saccharomyces cerevisiae Deltanep1 growth defect. SnR57 mediates 2'-O-ribose-methylation of G(1570) in the 18S rRNA. By performing a three-hybrid screen, we isolated several short RNA sequences with strong binding affinity to Nep1p. All isolated RNAs shared a six-nucleotide consensus motif C/UUCAAC. Furthermore, one of the isolated RNAs exactly corresponded to nucleotides 1553-1577 of the 18S rRNA, which includes G(1570), the site of snR57-dependent 18S rRNA methylation. From protein-protein crosslink data and the cryo-EM map of the S. cerevisiae small ribosomal subunit, we suggest that Rps19p is localized in close vicinity to the Nep1p 18S rRNA binding site. Our results suggest that Nep1p binds adjacent to helix 47 of the 18S rRNA and possibly supports the association of Rps19p to pre-ribosomal particles.

Characterization and complementation of a Fus3/Kss1 type MAPK from Tuber borchii, TBMK

Mitogen-activated protein kinases (MAPK) are used by organisms to transduce extra cellular signals from the environment in cellular events such as proliferation and differentiation. In the present study, we have characterized the first MAPK from the ectomycorrhizal fungus Tuber borchii (TBMK) which belongs to the YERK1 (yeast extra cellular regulated kinase) subfamily. TBMK is present as a single copy in the genome and the codified protein was phosphorylated during the interaction with the host plant, Tilia americana. Complementation studies showed that TBMK restores pheromone signaling in Saccharomyces cerevisiae and partially restores invasive growth of Fusarium oxysporum that lack the fmk1 gene. This suggests a protein kinase activity and its involvement in the infection processes. Hence, TBMK could play an important role during the pre-symbiotic phase of T. borchii with its host plant in the modulation of genes necessary for the establishment of symbiosis leading to the synthesis of functional ectomycorrhizae.

Meiotic proteins bqt1 and bqt2 tether telomeres to form the bouquet arrangement of chromosomes

In many organisms, meiotic chromosomes are bundled at their telomeres to form a "bouquet" arrangement. The bouquet formation plays an important role in homologous chromosome pairing and therefore progression of meiosis. As meiotic telomere clustering occurs in response to mating pheromone signaling in fission yeast, we looked for factors essential for bouquet formation among genes induced under mating pheromone signaling. This genome-wide search identified two proteins, Bqt1 and Bqt2, that connect telomeres to the spindle-pole body (SPB; the centrosome equivalent in fungi). Neither Bqt1 nor Bqt2 alone functions as a connector, but together the two proteins form a bridge between Rap1 (a telomere protein) and Sad1 (an SPB protein). Significantly, when both Bqt1 and Bqt2 are ectopically expressed in mitotic cells, they also form a bridge between Rap1 and Sad1. Thus, a complex including Bqt1 and Bqt2 is essential for connecting telomeres to the SPB.

Comparative genomics and genome evolution in yeasts

Yeasts provide a powerful model system for comparative genomics research. The availability of multiple complete genome sequences from different fungal groups--currently 18 hemiascomycetes, 8 euascomycetes and 4 basidiomycetes--enables us to gain a broad perspective on genome evolution. The sequenced genomes span a continuum of divergence levels ranging from multiple individuals within a species to species pairs with low levels of protein sequence identity and no conservation of gene order. One of the most interesting emerging areas is the growing number of events such as gene losses, gene displacements and gene relocations that can be attributed to the action of natural selection.

Genome-scale analysis reveals Sst2 as the principal regulator of mating pheromone signaling in the yeast Saccharomyces cerevisiae

A common property of G protein-coupled receptors is that they become less responsive with prolonged stimulation. Regulators of G protein signaling (RGS proteins) are well known to accelerate G protein GTPase activity and do so by stabilizing the transition state conformation of the G protein alpha subunit. In the yeast Saccharomyces cerevisiae there are four RGS-homologous proteins (Sst2, Rgs2, Rax1, and Mdm1) and two Galpha proteins (Gpa1 and Gpa2). We show that Sst2 is the only RGS protein that binds selectively to the transition state conformation of Gpa1. The other RGS proteins also bind Gpa1 and modulate pheromone signaling, but to a lesser extent and in a manner clearly distinct from Sst2. To identify other candidate pathway regulators, we compared pheromone responses in 4,349 gene deletion mutants representing nearly all nonessential genes in yeast. A number of mutants produced an increase (sst2, bar1, asc1, and ygl024w) or decrease (cla4) in pheromone sensitivity or resulted in pheromone-independent signaling (sst2, pbs2, gas1, and ygl024w). These findings suggest that Sst2 is the principal regulator of Gpa1-mediated signaling in vivo but that other proteins also contribute in distinct ways to pathway regulation.

Annotation of unknown yeast ORFs by correlation analysis of microarray data and extensive literature searches

Changes in the expression of genes were used to elucidate the metabolic pathways and regulatory mechanisms that respond to environmental or genetic modifications. Results from previously published chemostat datasets were merged with novel data generated in the present study. ORFs displaying significant changes in expression that correlated with those of other ORFs were analysed using GO mapping tools and supplemented by literature information. The strategy developed was used to propose annotations for ORFs of unknown function. The following ORFs were assigned functions as a result of this study: YMR090w, YGL157w, YGR243w, YLR327c, YER121w, YFR017c, YGR067c, YKL187c, YGR236c (SPG1), YMR107w (SPG4), YMR206w, YER067w, YJL103c, YNL175C (NOP13) YJL200C, YDL070C (FMP16) and YGR173W.

Support vector machine for classification of meiotic recombination hotspots and coldspots in Saccharomyces cerevisiae based on codon composition

Background

Meiotic double-strand breaks occur at relatively high frequencies in some genomic regions (hotspots) and relatively low frequencies in others (coldspots). Hotspots and coldspots are receiving increasing attention in research into the mechanism of meiotic recombination. However, predicting hotspots and coldspots from DNA sequence information is still a challenging task.

Results

We present a novel method for classification of hot and cold ORFs located in hotspots and coldspots respectively in Saccharomyces cerevisiae, using support vector machine (SVM), which relies on codon composition differences. This method has achieved a high classification accuracy of 85.0%. Since codon composition is a fusion of codon usage bias and amino acid composition signals, the ability of these two kinds of sequence attributes to discriminate hot ORFs from cold ORFs was also investigated separately. Our results indicate that neither codon usage bias nor amino acid composition taken separately performed as well as codon composition. Moreover, our SVM based method was applied to the full genome: We predicted the hot/cold ORFs from the yeast genome by using cutoffs of recombination rate. We found that the performance of our method for predicting cold ORFs is not as good as that for predicting hot ORFs. Besides, we also observed a considerable correlation between meiotic recombination rate and amino acid composition of certain residues, which probably reflects the structural and functional dissimilarity between the hot and cold groups.

Conclusion

We have introduced a SVM-based novel method to discriminate hot ORFs from cold ones. Applying codon composition as sequence attributes, we have achieved a high classification accuracy, which suggests that codon composition has strong potential to be used as sequence attributes in the prediction of hot and cold ORFs.

Synthesis of Novel Porous Magnetic Silica Microspheres as Adsorbents for Isolation of Genomic DNA

An improved procedure is described for preparation of novel mesoporous microspheres consisting of magnetic nanoparticles homogeneously dispersed in a silica matrix. The method is based on a three-step process, involving (i) formation of hematite/silica composite microspheres by urea-formaldehyde polymerization, (ii) calcination of the composite particles to remove the organic constituents, and (iii) in situ transformation of the iron oxide in the composites by hydrogen reductive reaction. The as-synthesized magnetite/silica composite microspheres were nearly monodisperse, mesoporous, and magnetizable, with as typical values an average diameter of 3.5 microm, a surface area of 250 m(2)/g, a pore size of 6.03 nm, and a saturation magnetization of 9.82 emu/g. These magnetic particles were tested as adsorbents for isolation of genomic DNA from Saccharomyces cerevisiae cells and maize kernels. The results are quite encouraging as the magnetic particle based protocols lead to the extraction of genomic DNA with satisfactory integrity, yield, and purity. Being hydrophilic in nature, the porous magnetic silica microspheres are considered a good alternative to polystyrene-based magnetic particles for use in biomedical applications where nonspecific adsorption of biomolecules is to be minimized.

The Nse5-Nse6 dimer mediates DNA repair roles of the Smc5-Smc6 complex

Stabilization and processing of stalled replication forks is critical for cell survival and genomic integrity. We characterize a novel DNA repair heterodimer of Nse5 and Nse6, which are nonessential nuclear proteins critical for chromosome segregation in fission yeast. The Nse5/6 dimer facilitates DNA repair as part of the Smc5-Smc6 holocomplex (Smc5/6), the basic architecture of which we define. Nse5-Nse6 [corrected] (Nse5 and Nse6) [corrected] mutants display a high level of spontaneous DNA damage and mitotic catastrophe in the absence of the master checkpoint regulator Rad3 (hATR). Nse5/6 mutants are required for the response to genotoxic agents that block the progression of replication forks, acting in a pathway that allows the tolerance of irreparable UV lesions. Interestingly, the UV sensitivity of Nse5/6 [corrected] is suppressed by concomitant deletion of the homologous recombination repair factor, Rhp51 (Rad51). Further, the viability of Nse5/6 mutants depends on Mus81 and Rqh1, factors that resolve or prevent the formation of Holliday junctions. Consistently, the UV sensitivity of cells lacking Nse5/6 can be partially suppressed by overexpressing the bacterial resolvase RusA. We propose a role for Nse5/6 mutants in suppressing recombination that results in Holliday junction formation or in Holliday junction resolution.

Functional partitioning of yeast co-expression networks after genome duplication

Several species of yeast, including the baker's yeast Saccharomyces cerevisiae, underwent a genome duplication roughly 100 million years ago. We analyze genetic networks whose members were involved in this duplication. Many networks show detectable redundancy and strong asymmetry in their interactions. For networks of co-expressed genes, we find evidence for network partitioning whereby the paralogs appear to have formed two relatively independent subnetworks from the ancestral network. We simulate the degeneration of networks after duplication and find that a model wherein the rate of interaction loss depends on the “neighborliness” of the interacting genes produces networks with parameters similar to those seen in the real partitioned networks. We propose that the rationalization of network structure through the loss of pair-wise gene interactions after genome duplication provides a mechanism for the creation of semi-independent daughter networks through the division of ancestral functions between these daughter networks.

An analysis of how duplicated networks of genes (as a result of whole genome duplication in yeast) evolved shows that network partitioning occurred through loss of interactions, resulting in independent subnetworks.

Agrobacterium-mediated transformation of Sclerotinia sclerotiorum

Ascospores from the phytopathogenic fungus Sclerotinia sclerotiorum were transformed to hygromycin B resistance by co-cultivation with Agrobacterium tumefaciens. Transformed spores germinated and grew on PDA supplemented with 100 ug/ml hygromycin B. The presence of mitotically stable hph gene integration at random sites in the genome was confirmed by PCR and Southern blot analysis. A transformation frequency of 8 x 10(-5) was achieved in five separate experiments. This study is the first report of success co-cultivating A. tumefaciens with S. sclerotiorum. This report of a reproducible Agrobacterium-mediated transformation method should allow the development of T-DNA tagging as a system for insertional mutagenesis in S. sclerotiorum and provide a simple and reliable method for genetic manipulation.

Overlapping contributions of Msh1p and putative recombination proteins Cce1p, Din7p, and Mhr1p in large-scale recombination and genome sorting events in the mitochondrial genome of Saccharomyces cerevisiae

The mechanisms that govern mutation avoidance in the mitochondrial genome, though believed to be numerous, are poorly understood. The identification of individual genes has implicated mismatch repair and several recombination pathways in maintaining the fidelity and structural stability of mitochondrial DNA. However, the majority of genes in these pathways have not been identified and the interactions between different pathways have not been extensively studied. Additionally, the multicopy presence of the mitochondrial genome affects the occurrence and persistence of mutant phenotypes, making mitochondrial DNA transmission and sorting important factors affecting mutation accumulation. We present new evidence that the putative recombination genes CCE1, DIN7, and MHR1 have overlapping function with the mismatch repair homolog MSH1 in point mutation avoidance and suppression of aberrant recombination events. In addition, we demonstrate a novel role for Msh1p in mtDNA transmission, a role not predicted by studies of its nuclear homologs.

Revealing hidden relationships among yeast genes involved in chromosome segregation using systematic synthetic lethal and synthetic dosage lethal screens

The vast accumulation of knowledge from genome sequencing projects and studies with model organisms has presented a remarkable challenge to biologists: to understand the functions of thousands of highly conserved genes and how they work together to regulate fundamental cellular processes. This challenge is compounded by the inescapable reality that most genes are 'buffered' by other genes that contribute to the same biological processes, limiting the impact of phenotypic studies with single mutants. In budding yeast, functional genomic methods have been developed for the systematic application of established genetic techniques. In particular, the Synthetic Genetic Array (SGA) method allows genome-wide synthetic lethal (SL) and synthetic dosage lethal (SDL) screens thus enabling an unbiased survey of genetic interactions. We have used genes encoding components of the yeast kinetochore as a biological testbed for assaying the utility of SGA-based SL and SDL screens for revealing new pathways and genes involved in chromosome segregation. We identified 211 nonessential deletion mutants that were unable to tolerate either overexpression or loss of function of kinetochore genes. Our study uncovered a wealth of relationships between gene products that functionally interact with the kinetochore, and also highlighted the value of performing genome-wide screens with both hypomorphic and hypermorphic alleles of query genes. Here, we will highlight our recent kinetochore SGA genomic screens, in the broader context of applying complementary genetic screening approaches in the systematic exploration of biological pathways or functional complexes.

Disentangling information flow in the Ras-cAMP signaling network

The perturbation of signal-transduction molecules elicits genomic-expression effects that are typically neither restricted to a small set of genes nor uniform. Instead there are broad, varied, and complex changes in expression across the genome. These observations suggest that signal transduction is not mediated by isolated pathways of information flow to distinct groups of genes in the genome. Rather, multiple entangled paths of information flow influence overlapping sets of genes. Using the Ras-cAMP pathway in Saccharomyces cerevisiae as a model system, we perturbed key pathway elements and collected genomic-expression data. Singular value decomposition was applied to separate the genome-wide transcriptional response into weighted expression components exhibited by overlapping groups of genes. Molecular interaction data were integrated to connect gene groups to perturbed signaling elements. The resulting series of linked subnetworks maps multiple putative pathways of information flow through a dense signaling network, and provides a set of testable hypotheses for complex gene-expression effects across the genome.

GONOME: measuring correlations between GO terms and genomic positions

Background:

Current methods to find significantly under- and over-represented gene ontology (GO) terms in a set of genes consider the genes as equally probable "balls in a bag", as may be appropriate for transcripts in micro-array data. However, due to the varying length of genes and intergenic regions, that approach is inappropriate for deciding if any GO terms are correlated with a set of genomic positions.

Results:

We present an algorithm – GONOME – that can determine which GO terms are significantly associated with a set of genomic positions given a genome annotated with (at least) the starts and ends of genes. We show that certain GO terms may appear to be significantly associated with a set of randomly chosen positions in the human genome if gene lengths are not considered, and that these same terms have been reported as significantly over-represented in a number of recent papers. This apparent over-representation disappears when gene lengths are considered, as GONOME does. For example, we show that, when gene length is taken into account, the term "development" is not significantly enriched in genes associated with human CpG islands, in contradiction to a previous report. We further demonstrate the efficacy of GONOME by showing that occurrences of the proteosome-associated control element (PACE) upstream activating sequence in the S. cerevisiae genome associate significantly to appropriate GO terms. An extension of this approach yields a whole-genome motif discovery algorithm that allows identification of many other promoter sequences linked to different types of genes, including a large group of previously unknown motifs significantly associated with the terms 'translation' and 'translational elongation'.

Conclusion:

GONOME is an algorithm that correctly extracts over-represented GO terms from a set of genomic positions. By explicitly considering gene size, GONOME avoids a systematic bias toward GO terms linked to large genes. Inappropriate use of existing algorithms that do not take gene size into account has led to erroneous or suspect conclusions. Reciprocally GONOME may be used to identify new features in genomes that are significantly associated with particular categories of genes.

Accumulation of unstable promoter-associated transcripts upon loss of the nuclear exosome subunit Rrp6p in Saccharomyces cerevisiae

Mutations in RRP6 result in the accumulation of aberrant polyadenylated transcripts from small nucleolar RNA genes. We exploited this observation to search for novel noncoding RNA genes in the yeast genome. When RNA from rrp6Delta yeast is compared with wild-type on whole-genome microarrays, numerous intergenic loci exhibit an increased mutant/wild type signal ratio. Among these loci, we found one encoding a new C/D box small nucleolar RNA, as well as a surprising number that gave rise to heterogeneous Trf4p-polyadenylated RNAs with lengths of approximately 250-500 nt. This class of RNAs is not easily detected in wild-type cells and appears associated with promoters. Fine mapping of several such transcripts shows they originate near known promoter elements but do not usually extend far enough to act as mRNAs, and may regulate the transcription of downstream mRNAs. Rather than being uninformative transcriptional "noise," we hypothesize that these transcripts reflect important features of RNA polymerase activity at the promoter. This activity is normally undetectable in wild-type cells because the transcripts are somehow distinguished from true mRNAs and are degraded in an Rrp6p-dependent fashion in the nucleus.

Agrobacterium T-DNA-mediated integration and gene replacement in the brown rot pathogen Monilinia fructicola

A transformation system utilizing Agrobacterium tumefaciens was developed for targeted gene disruption in Monilinia fructicola, a fungal pathogen that causes brown rot disease of stone fruits. Transformation with a vector containing the neomycin phosphotransferase II (nptII) cassette flanked with 4 kb cutinase gene (Mfcut1) flanking sequences resulted in an average of 13 transformants per 10(5) spores. When assayed by PCR and DNA blot analyses, more than 50% of the transformants recovered had integrated in the targeted Mfcut1 locus. Both target-gene-specific and non-specific integrations carried direct (head-to-tail) repeat T-DNA integrations. Sequence analysis of these T-DNA integrations revealed that 26 bp of the T-DNA right border were missing at the junctions between direct repeats in all cases. The recombination event during non-specific T-DNA integration in this fungus was unlike that reported in Agrobacterium-mediated transformation in plants.

[Analysis of the secreted proteins encoded by genes in genome of filamental fungus (Neurospora crassa)]

The completed genome sequence of filamental fungus Neurospora crassa has completed and 10,082 open reading frames (ORFs) were predicted in whole genome. Among them, 437 proteins were secreted ones based on the combination of prediction algorithms SignalP v3.0 and PSORT, transmembrane domains prediction algorithms TMHMM v2.0 and THUMBUP, potential GPI-anchor sites prediction algorithm big-PI predictor and the subcellular protein location prediction algorithm TargetP v1.01. The minimum ORF length of the encoding sequences for 437 secreted proteins was 252 bp, the maximum ORF length was 6 604 bp and average was 1 433 bp. The length of signal peptide of secreted proteins ranged from 15 to 59 amino acids. Among the 437 secreted proteins, 205 were provided with function description. The function of secretory proteins were involved in kinds of enzyme, cell energy, transmission, cell recovered itself and defense mechanism et al. These results implied that processes catalyzed by these proteins probably occurred outside of cells, at least partially.

The complete DNA sequence of the mitochondrial genome of the dermatophyte fungus Epidermophyton floccosum

We report here the complete nucleotide sequence of the 30.9-kb mitochondrial genome of the dermatophyte fungus Epidermophyton floccosum. All genes are encoded on the same DNA strand and include seven subunits of the reduced nicotinamide adenine dinucleotide ubiquinone oxireductase (nad1, nad2, nad3, nad4, nad4L, nad5, and nad6), three subunits of cytochrome oxidase (cox1, cox2, and cox3), apocytochrome b (cob), three subunits of ATP synthase (atp6, atp8, and atp9), the small and large ribosomal RNAs (rns and rnl), and 25 tRNAs. A ribosomal protein gene (rps5) is present as an intronic ORF in the large ribosomal subunit. The genes coding for cob and cox1 carry one intron and nad5 carries two introns with ORFs. The mtDNA of E. floccosum has the same gene order as Trichophyton rubrum mtDNA, with the exception of some tRNA genes. Maximum likelihood phylogenetic analysis confirms T. rubrum as a close relative of E. floccosum. This is the first complete mitochondrial sequence of a species of the order Onygenales. This sequence is available under GenBank accession number AY916130.

Shape-specific recognition in the structure of the Vts1p SAM domain with RNA

Although the abundant sterile alpha motif (SAM) domain was originally classified as a protein-protein interaction domain, it has recently been shown that certain SAM domains have the ability to bind RNA, defining a new type of post-transcriptional gene regulator. To further understand the function of SAM-RNA recognition, we determined the solution structures of the SAM domain of the Saccharomyces cerevisiae Vts1p (Vts1p-SAM) and the Smaug response element (SRE) stem-loop RNA as a complex and in isolation. The structures show that Vts1p-SAM recognizes predominantly the shape of the SRE rather than its sequence, with the exception of a G located at the tip of the pentaloop. Using microarray gene profiling, we identified several genes in S. cerevisiae that seem to be regulated by Vts1p and contain one or more copies of the SRE.

The dawn of fungal pathogen genomics

Recent advances in sequencing technologies have led to a remarkable increase in the number of sequenced fungal genomes. Several important plant pathogenic fungi are among those that have been sequenced or are being sequenced. Additional fungal pathogens are likely to be sequenced in the near future. Analysis of the available genomes has provided useful information about genes that may be important for plant infection and colonization. Genome features, such as repetitive sequences, telomeres, conserved syntenic blocks, and expansion of pathogenicity-related genes, are discussed in detail with Magnaporthe oryzae (M. grisea) and Fusarium graminearum as examples. Functional and comparative genomic studies in plant pathogenic fungi, although still in the early stages and limited to a few pathogens, have enormous potential to improve our understanding of the molecular mechanisms involved in host-pathogen interactions. Development of advanced genomics tools and infrastructure is critical for efficient utilization of the vast wealth of available genome sequence information and will form a solid foundation for systems biology studies of plant pathogenic fungi.

Evolution of thehemiascomycete yeasts: on life styles and the importance of inbreeding

The term 'breeding system' is used to describe the morphological and behavioural aspects of the sexual life cycle of a species. The yeast breeding system provides three alternatives that enable hapoids to return to the diploid state that is necessary for meiosis: mating of unrelated haploids (amphimixis), mating between spores from the same tetrad (intratetrad mating, automixis) and mother daughter mating upon mating type switching (haplo-selfing). The frequency of specific mating events affects the level of heterozygosity present in individuals and the genetic diversity of populations. This review discusses the reproductive strategies of yeasts, in particular S. cerevisiae (Bakers' or budding yeast). Emphasis is put on intratetrad mating, its implication for diversity, and how the particular genome structure could have evolved to ensure the preservation of a high degree of heterozygosity in conjunction with frequent intratetrad matings. I also discuss how the ability of yeast to control the number of spores that are formed accounts for high intratetrad mating rates and for enhanced transmission of genomic variation. I extend the discussion to natural genetic variation and propose that a high level of plasticity is inherent in the yeast breeding system, which may allow variation of the breeding behaviour in accordance with the needs imposed by the environment.

Balance between transcription and RNA degradation is vital for Saccharomyces cerevisiae mitochondria: reduced transcription rescues the phenotype of deficient RNA degradation

The Saccharomyces cerevisiae SUV3 gene encodes the helicase component of the mitochondrial degradosome (mtEXO), the principal 3'-to-5' exoribonuclease of yeast mitochondria responsible for RNA turnover and surveillance. Inactivation of SUV3 (suv3Delta) causes multiple defects related to overaccumulation of aberrant transcripts and precursors, leading to a disruption of mitochondrial gene expression and loss of respiratory function. We isolated spontaneous suppressors that partially restore mitochondrial function in suv3Delta strains devoid of mitochondrial introns and found that they correspond to partial loss-of-function mutations in genes encoding the two subunits of the mitochondrial RNA polymerase (Rpo41p and Mtf1p) that severely reduce the transcription rate in mitochondria. These results show that reducing the transcription rate rescues defects in RNA turnover and demonstrates directly the vital importance of maintaining the balance between RNA synthesis and degradation.

Mining yeast in silico unearths a golden nugget for mitochondrial biology

NADH:ubiquinone oxidoreductase (complex I) of the electron transport chain is a multimeric mitochondrial enzyme of approximately 1000 kDa consisting of 46 different proteins encoded by both the mitochondrial and nuclear genomes. Little is known about the cellular mechanisms and protein chaperones that guide its assembly. In this issue of the JCI, Ogilvie et al. use genomic sequence data to compare the proteins produced by yeasts with and without complex I in order to generate a list of proteins whose human orthologs might serve as complex I assembly proteins. The gene encoding one of these candidate proteins, B17.2L, was found to harbor a nonsense mutation in one of 28 patients with a deficiency of complex I. B17.2L associated with subcomplexes that are seen when complex I assembly is incomplete. The research described here combines clever model organism genomics and bioinformatics with sophisticated human molecular and biochemical genetics to identify the first mammalian protein required for the normal assembly of complex I.