ManiNetCluster: a novel manifold learning approach to reveal the functional links between gene networks.

BACKGROUND: The coordination of genomic functions is a critical and complex process across biological systems such as phenotypes or states (e.g., time, disease, organism, environmental perturbation). Understanding how the complexity of genomic function relates to these states remains a challenge. To address this, we have developed a novel computational method, ManiNetCluster, which simultaneously aligns and clusters gene networks (e.g., co-expression) to systematically reveal the links of genomic function between different conditions. Specifically, ManiNetCluster employs manifold learning to uncover and match local and non-linear structures among networks, and identifies cross-network functional links. RESULTS: We demonstrated that ManiNetCluster better aligns the orthologous genes from their developmental expression profiles across model organisms than state-of-the-art methods (p-value <2.2×10-16). This indicates the potential non-linear interactions of evolutionarily conserved genes across species in development. Furthermore, we applied ManiNetCluster to time series transcriptome data measured in the green alga Chlamydomonas reinhardtii to discover the genomic functions linking various metabolic processes between the light and dark periods of a diurnally cycling culture. We identified a number of genes putatively regulating processes across each lighting regime. CONCLUSIONS: ManiNetCluster provides a novel computational tool to uncover the genes linking various functions from different networks, providing new insight on how gene functions coordinate across different conditions. ManiNetCluster is publicly available as an R package at https://github.com/daifengwanglab/ManiNetCluster.

High-Resolution Profiling of a Synchronized Diurnal Transcriptome from Chlamydomonas reinhardtii Reveals Continuous Cell and Metabolic Differentiation.

The green alga Chlamydomonas reinhardtii is a useful model organism for investigating diverse biological processes, such as photosynthesis and chloroplast biogenesis, flagella and basal body structure/function, cell growth and division, and many others. We combined a highly synchronous photobioreactor culture system with frequent temporal sampling to characterize genome-wide diurnal gene expression in Chlamydomonas. Over 80% of the measured transcriptome was expressed with strong periodicity, forming 18 major clusters. Genes associated with complex structures and processes, including cell cycle control, flagella and basal bodies, ribosome biogenesis, and energy metabolism, all had distinct signatures of coexpression with strong predictive value for assigning and temporally ordering function. Importantly, the frequent sampling regime allowed us to discern meaningful fine-scale phase differences between and within subgroups of genes and enabled the identification of a transiently expressed cluster of light stress genes. Coexpression was further used both as a data-mining tool to classify and/or validate genes from other data sets related to the cell cycle and to flagella and basal bodies and to assign isoforms of duplicated enzymes to their cognate pathways of central carbon metabolism. Our diurnal coexpression data capture functional relationships established by dozens of prior studies and are a valuable new resource for investigating a variety of biological processes in Chlamydomonas and other eukaryotes.

Nitrogen-Sparing Mechanisms in Chlamydomonas Affect the Transcriptome, the Proteome, and Photosynthetic Metabolism.

Nitrogen (N) is a key nutrient that limits global primary productivity; hence, N-use efficiency is of compelling interest in agriculture and aquaculture. We used Chlamydomonas reinhardtii as a reference organism for a multicomponent analysis of the N starvation response. In the presence of acetate, respiratory metabolism is prioritized over photosynthesis; consequently, the N-sparing response targets proteins, pigments, and RNAs involved in photosynthesis and chloroplast function over those involved in respiration. Transcripts and proteins of the Calvin-Benson cycle are reduced in N-deficient cells, resulting in the accumulation of cycle metabolic intermediates. Both cytosolic and chloroplast ribosomes are reduced, but via different mechanisms, reflected by rapid changes in abundance of RNAs encoding chloroplast ribosomal proteins but not cytosolic ones. RNAs encoding transporters and enzymes for metabolizing alternative N sources increase in abundance, as is appropriate for the soil environmental niche of C. reinhardtii. Comparison of the N-replete versus N-deplete proteome indicated that abundant proteins with a high N content are reduced in N-starved cells, while the proteins that are increased have lower than average N contents. This sparing mechanism contributes to a lower cellular N/C ratio and suggests an approach for engineering increased N-use efficiency.

Genome-wide analysis on Chlamydomonas reinhardtii reveals the impact of hydrogen peroxide on protein stress responses and overlap with other stress transcriptomes.

Reactive oxygen species (ROS) are produced by and have the potential to be damaging to all aerobic organisms. In photosynthetic organisms, they are an unavoidable byproduct of electron transfer in both the chloroplast and mitochondrion. Here, we employ the reference unicellular green alga Chlamydomonas reinhardtii to identify the effect of H2O2 on gene expression by monitoring the changes in the transcriptome in a time-course experiment. Comparison of transcriptomes from cells sampled immediately prior to the addition of H2O2 and 0.5 and 1 h subsequently revealed 1278 differentially abundant transcripts. Of those transcripts that increase in abundance, many encode proteins involved in ROS detoxification, protein degradation and stress responses, whereas among those that decrease are transcripts encoding proteins involved in photosynthesis and central carbon metabolism. In addition to these transcriptomic adjustments, we observe that addition of H2O2 is followed by an accumulation and oxidation of the total intracellular glutathione pool, and a decrease in photosynthetic O2 output. Additionally, we analyze our transcriptomes in the context of changes in transcript abundance in response to singlet O2 (O2*), and relate our H2O2 -induced transcripts to a diurnal transcriptome, where we demonstrate enrichments of H2O2 -induced transcripts early in the light phase, late in the light phase and 2 h prior to light. On this basis several genes that are highlighted in this work may be involved in previously undiscovered stress remediation pathways or acclimation responses.

Deciphering the Translation Initiation Factor 5A Modification Pathway in Halophilic Archaea.

Translation initiation factor 5A (IF5A) is essential and highly conserved in Eukarya (eIF5A) and Archaea (aIF5A). The activity of IF5A requires hypusine, a posttranslational modification synthesized in Eukarya from the polyamine precursor spermidine. Intracellular polyamine analyses revealed that agmatine and cadaverine were the main polyamines produced in Haloferax volcanii in minimal medium, raising the question of how hypusine is synthesized in this halophilic Archaea. Metabolic reconstruction led to a tentative picture of polyamine metabolism and aIF5A modification in Hfx. volcanii that was experimentally tested. Analysis of aIF5A from Hfx. volcanii by LC-MS/MS revealed it was exclusively deoxyhypusinylated. Genetic studies confirmed the role of the predicted arginine decarboxylase gene (HVO_1958) in agmatine synthesis. The agmatinase-like gene (HVO_2299) was found to be essential, consistent with a role in aIF5A modification predicted by physical clustering evidence. Recombinant deoxyhypusine synthase (DHS) from S. cerevisiae was shown to transfer 4-aminobutyl moiety from spermidine to aIF5A from Hfx. volcanii in vitro. However, at least under conditions tested, this transfer was not observed with the Hfx. volcanii DHS. Furthermore, the growth of Hfx. volcanii was not inhibited by the classical DHS inhibitor GC7. We propose a model of deoxyhypusine synthesis in Hfx. volcanii that differs from the canonical eukaryotic pathway, paving the way for further studies.

Exploiting algal NADPH oxidase for biophotovoltaic energy.

Photosynthetic microbes exhibit light-dependent electron export across the cell membrane, which can generate electricity in biological photovoltaic (BPV) devices. How electrons are exported remains to be determined; the identification of mechanisms would help selection or generation of photosynthetic microbes capable of enhanced electrical output. We show that plasma membrane NADPH oxidase activity is a significant component of light-dependent generation of electricity by the unicellular green alga Chlamydomonas reinhardtii. NADPH oxidases export electrons across the plasma membrane to form superoxide anion from oxygen. The C. reinhardtii mutant lacking the NADPH oxidase encoded by RBO1 is impaired in both extracellular superoxide anion production and current generation in a BPV device. Complementation with the wild-type gene restores both capacities, demonstrating the role of the enzyme in electron export. Monitoring light-dependent extracellular superoxide production with a colorimetric assay is shown to be an effective way of screening for electrogenic potential of candidate algal strains. The results show that algal NADPH oxidases are important for superoxide anion production and open avenues for optimizing the biological component of these devices.

Systems-level analysis of nitrogen starvation-induced modifications of carbon metabolism in a Chlamydomonas reinhardtii starchless mutant.

To understand the molecular basis underlying increased triacylglycerol (TAG) accumulation in starchless (sta) Chlamydomonas reinhardtii mutants, we undertook comparative time-course transcriptomics of strains CC-4348 (sta6 mutant), CC-4349, a cell wall-deficient (cw) strain purported to represent the parental STA6 strain, and three independent STA6 strains generated by complementation of sta6 (CC-4565/STA6-C2, CC-4566/STA6-C4, and CC-4567/STA6-C6) in the context of N deprivation. Despite N starvation-induced dramatic remodeling of the transcriptome, there were relatively few differences (5 × 10(2)) observed between sta6 and STA6, the most dramatic of which were increased abundance of transcripts encoding key regulated or rate-limiting steps in central carbon metabolism, specifically isocitrate lyase, malate synthase, transaldolase, fructose bisphosphatase and phosphoenolpyruvate carboxykinase (encoded by ICL1, MAS1, TAL1, FBP1, and PCK1 respectively), suggestive of increased carbon movement toward hexose-phosphate in sta6 by upregulation of the glyoxylate pathway and gluconeogenesis. Enzyme assays validated the increase in isocitrate lyase and malate synthase activities. Targeted metabolite analysis indicated increased succinate, malate, and Glc-6-P and decreased Fru-1,6-bisphosphate, illustrating the effect of these changes. Comparisons of independent data sets in multiple strains allowed the delineation of a sequence of events in the global N starvation response in C. reinhardtii, starting within minutes with the upregulation of alternative N assimilation routes and carbohydrate synthesis and subsequently a more gradual upregulation of genes encoding enzymes of TAG synthesis. Finally, genome resequencing analysis indicated that (1) the deletion in sta6 extends into the neighboring gene encoding respiratory burst oxidase, and (2) a commonly used STA6 strain (CC-4349) as well as the sequenced reference (CC-503) are not congenic with respect to sta6 (CC-4348), underscoring the importance of using complemented strains for more rigorous assignment of phenotype to genotype.

Activation of Autophagy by Metals in Chlamydomonas reinhardtii.

Autophagy is an intracellular self-degradation pathway by which eukaryotic cells recycle their own material in response to specific stress conditions. Exposure to high concentrations of metals causes cell damage, although the effect of metal stress on autophagy has not been explored in photosynthetic organisms. In this study, we investigated the effect of metal excess on autophagy in the model unicellular green alga Chlamydomonas reinhardtii. We show in cells treated with nickel an upregulation of ATG8 that is independent of CRR1, a global regulator of copper signaling in Chlamydomonas. A similar effect on ATG8 was observed with copper and cobalt but not with cadmium or mercury ions. Transcriptome sequencing data revealed an increase in the abundance of the protein degradation machinery, including that responsible for autophagy, and a substantial overlap of that increased abundance with the hydrogen peroxide response in cells treated with nickel ions. Thus, our results indicate that metal stress triggers autophagy in Chlamydomonas and suggest that excess nickel may cause oxidative damage, which in turn activates degradative pathways, including autophagy, to clear impaired components and recover cellular homeostasis.

The archaeal COG1901/DUF358 SPOUT-methyltransferase members, together with pseudouridine synthase Pus10, catalyze the formation of 1-methylpseudouridine at position 54 of tRNA.

The methylation of pseudouridine (Ψ) at position 54 of tRNA, producing m(1)Ψ, is a hallmark of many archaeal species, but the specific methylase involved in the formation of this modification had yet to be characterized. A comparative genomics analysis had previously identified COG1901 (DUF358), part of the SPOUT superfamily, as a candidate for this missing methylase family. To test this prediction, the COG1901 encoding gene, HVO_1989, was deleted from the Haloferax volcanii genome. Analyses of modified base contents indicated that while m(1)Ψ was present in tRNA extracted from the wild-type strain, it was absent from tRNA extracted from the mutant strain. Expression of the gene encoding COG1901 from Halobacterium sp. NRC-1, VNG1980C, complemented the m(1)Ψ minus phenotype of the ΔHVO_1989 strain. This in vivo validation was extended with in vitro tests. Using the COG1901 recombinant enzyme from Methanocaldococcus jannaschii (Mj1640), purified enzyme Pus10 from M. jannaschii and full-size tRNA transcripts or TΨ-arm (17-mer) fragments as substrates, the sequential pathway of m(1)Ψ54 formation in Archaea was reconstituted. The methylation reaction is AdoMet dependent. The efficiency of the methylase reaction depended on the identity of the residue at position 55 of the TΨ-loop. The presence of Ψ55 allowed the efficient conversion of Ψ54 to m(1)Ψ54, whereas in the presence of C55, the reaction was rather inefficient and no methylation reaction occurred if a purine was present at this position. These results led to renaming the Archaeal COG1901 members as TrmY proteins.

Biochemical characterization of Fsa16295Glu from "Fervidibacter sacchari," the first hyperthermophilic GH50 with ß-1,3-endoglucanase activity and founding member of the subfamily GH50_3.

The aerobic hyperthermophile "Fervidibacter sacchari" catabolizes diverse polysaccharides and is the only cultivated member of the class "Fervidibacteria" within the phylum Armatimonadota. It encodes 117 putative glycoside hydrolases (GHs), including two from GH family 50 (GH50). In this study, we expressed, purified, and functionally characterized one of these GH50 enzymes, Fsa16295Glu. We show that Fsa16295Glu is a ß-1,3-endoglucanase with optimal activity on carboxymethyl curdlan (CM-curdlan) and only weak agarase activity, despite most GH50 enzymes being described as ß-agarases. The purified enzyme has a wide temperature range of 4-95°C (optimal 80°C), making it the first characterized hyperthermophilic representative of GH50. The enzyme is also active at a broad pH range of at least 5.5-11 (optimal 6.5-10). Fsa16295Glu possesses a relatively high kcat/KM of 1.82 × 107 s-1 M-1 with CM-curdlan and degrades CM-curdlan nearly completely to sugar monomers, indicating preferential hydrolysis of glucans containing ß-1,3 linkages. Finally, a phylogenetic analysis of Fsa16295Glu and all other GH50 enzymes revealed that Fsa16295Glu is distant from other characterized enzymes but phylogenetically related to enzymes from thermophilic archaea that were likely acquired horizontally from "Fervidibacteria." Given its functional and phylogenetic novelty, we propose that Fsa16295Glu represents a new enzyme subfamily, GH50_3.

A hemoprotein with a zinc-mirror heme site ties heme availability to carbon metabolism in cyanobacteria.

Heme has a critical role in the chemical framework of the cell as an essential protein cofactor and signaling molecule that controls diverse processes and molecular interactions. Using a phylogenomics-based approach and complementary structural techniques, we identify a family of dimeric hemoproteins comprising a domain of unknown function DUF2470. The heme iron is axially coordinated by two zinc-bound histidine residues, forming a distinct two-fold symmetric zinc-histidine-iron-histidine-zinc site. Together with structure-guided in vitro and in vivo experiments, we further demonstrate the existence of a functional link between heme binding by Dri1 (Domain related to iron 1, formerly ssr1698) and post-translational regulation of succinate dehydrogenase in the cyanobacterium Synechocystis, suggesting an iron-dependent regulatory link between photosynthesis and respiration. Given the ubiquity of proteins containing homologous domains and connections to heme metabolism across eukaryotes and prokaryotes, we propose that DRI (Domain Related to Iron; formerly DUF2470) functions at the molecular level as a heme-dependent regulatory domain.

Quinolinate salvage and insights for targeting NAD biosynthesis in group A streptococci.

The essential coenzyme NAD plays important roles in metabolic reactions and cell regulation in all organisms. As such, NAD synthesis has been investigated as a source for novel antibacterial targets. Cross-species genomics-based reconstructions of NAD metabolism in group A streptococci (GAS), combined with focused experimental testing in Streptococcus pyogenes, led to a better understanding of NAD metabolism in the pathogen. The predicted niacin auxotrophy was experimentally verified, as well as the essential role of the nicotinamidase PncA in the utilization of nicotinamide (Nm). PncA is dispensable in the presence of nicotinate (Na), ruling it out as a viable antibacterial target. The function of the "orphan" NadC enzyme, which is uniquely present in all GAS species despite the absence of other genes of NAD de novo synthesis, was elucidated. Indeed, the quinolinate (Qa) phosphoribosyltransferase activity of NadC from S. pyogenes allows the organism to sustain growth when Qa is present as a sole pyridine precursor. Finally, the redundancy of functional upstream salvage pathways in GAS species narrows the choice of potential drug targets to the two indispensable downstream enzymes of NAD synthesis, nicotinate adenylyltransferase (NadD family) and NAD synthetase (NadE family). Biochemical characterization of NadD confirmed its functional role in S. pyogenes, and its potential as an antibacterial target was supported by inhibition studies with previously identified class I inhibitors of the NadD enzyme family. One of these inhibitors efficiently inhibited S. pyogenes NadD (sp.NadD) in vitro (50% inhibitory concentration [IC(50)], 15 µM), exhibiting a noncompetitive mechanism with a K(i) of 8 µM.

Three acyltransferases and nitrogen-responsive regulator are implicated in nitrogen starvation-induced triacylglycerol accumulation in Chlamydomonas.

Algae have recently gained attention as a potential source for biodiesel; however, much is still unknown about the biological triggers that cause the production of triacylglycerols. We used RNA-Seq as a tool for discovering genes responsible for triacylglycerol (TAG) production in Chlamydomonas and for the regulatory components that activate the pathway. Three genes encoding acyltransferases, DGAT1, DGTT1, and PDAT1, are induced by nitrogen starvation and are likely to have a role in TAG accumulation based on their patterns of expression. DGAT1 and DGTT1 also show increased mRNA abundance in other TAG-accumulating conditions (minus sulfur, minus phosphorus, minus zinc, and minus iron). Insertional mutants, pdat1-1 and pdat1-2, accumulate 25% less TAG compared with the parent strain, CC-4425, which demonstrates the relevance of the trans-acylation pathway in Chlamydomonas. The biochemical functions of DGTT1 and PDAT1 were validated by rescue of oleic acid sensitivity and restoration of TAG accumulation in a yeast strain lacking all acyltransferase activity. Time course analyses suggest than a SQUAMOSA promoter-binding protein domain transcription factor, whose mRNA increases precede that of lipid biosynthesis genes like DGAT1, is a candidate regulator of the nitrogen deficiency responses. An insertional mutant, nrr1-1, accumulates only 50% of the TAG compared with the parental strain in nitrogen-starvation conditions and is unaffected by other nutrient stresses, suggesting the specificity of this regulator for nitrogen-deprivation conditions.

Pseudouridine formation in archaeal RNAs: The case of Haloferax volcanii.

Pseudouridine (Ψ), the isomer of uridine, is commonly found at various positions of noncoding RNAs of all organisms. Ψ residues are formed by a number of single- or multisite specific Ψ synthases, which generally act as stand-alone proteins. In addition, in Eukarya and Archaea, specific ribonucleoprotein complexes, each containing a distinct box H/ACA guide RNA and four core proteins, can produce Ψ at many sites of different cellular RNAs. Cbf5 is the core Ψ synthase in these complexes. Using Haloferax volcanii as an archaeal model organism, we show that, contrary to eukaryotes, the Cbf5 homolog (HVO_2493) is not essential in this archaeon. The Cbf5-deleted strain of H. volcanii completely lacks Ψ at positions 1940, 1942, 2605, and 2591 (Escherichia coli positions 1915, 1917, 2572, and 2586) of its 23S rRNA, and contains reduced steady-state levels of some box H/ACA RNAs. Archaeal Cbf5 is known to have tRNA Ψ55 synthase activity in vitro but we could not confirm this activity in vivo in H. volcanii. Conversely, the Pus10 (previously PsuX) homolog (HVO_1979), which can produce tRNA Ψ55, as well as Ψ54 in vitro, is shown here to be essential in H. volcanii, whereas the corresponding tRNA Ψ55 synthases, Pus4 and TruB, are not essential in yeast and E. coli, respectively. Finally, we demonstrate that HVO_1852, the TruA/Pus3 homolog, is responsible for the pseudouridylation of position 39 in H. volcanii tRNAs and that the corresponding gene is not essential.

A role for tetrahydrofolates in the metabolism of iron-sulfur clusters in all domains of life.

Iron-sulfur (Fe/S) cluster enzymes are crucial to life. Their assembly requires a suite of proteins, some of which are specific for particular subsets of Fe/S enzymes. One such protein is yeast Iba57p, which aconitase and certain radical S-adenosylmethionine enzymes require for activity. Iba57p homologs occur in all domains of life; they belong to the COG0354 protein family and are structurally similar to various folate-dependent enzymes. We therefore investigated the possible relationship between folates and Fe/S cluster enzymes using the Escherichia coli Iba57p homolog, YgfZ. NMR analysis confirmed that purified YgfZ showed stereoselective folate binding. Inactivating ygfZ reduced the activities of the Fe/S tRNA modification enzyme MiaB and certain other Fe/S enzymes, although not aconitase. When successive steps in folate biosynthesis were ablated, folE (lacking pterins and folates) and folP (lacking folates) mutants mimicked the ygfZ mutant in having low MiaB activities, whereas folE thyA mutants supplemented with 5-formyltetrahydrofolate (lacking pterins and depleted in dihydrofolate) and gcvP glyA mutants (lacking one-carbon tetrahydrofolates) had intermediate MiaB activities. These data indicate that YgfZ requires a folate, most probably tetrahydrofolate. Importantly, the ygfZ mutant was hypersensitive to oxidative stress and grew poorly on minimal media. COG0354 genes of bacterial, archaeal, fungal, protistan, animal, or plant origin complemented one or both of these growth phenotypes as well as the MiaB activity phenotype. Comparative genomic analysis indicated widespread functional associations between COG0354 proteins and Fe/S cluster metabolism. Thus COG0354 proteins have an ancient, conserved, folate-dependent function in the activity of certain Fe/S cluster enzymes.

A predicted CRISPR-mediated symbiosis between uncultivated archaea.

CRISPR-Cas systems defend prokaryotic cells from invasive DNA of viruses, plasmids and other mobile genetic elements. Here, we show using metagenomics, metatranscriptomics and single-cell genomics that CRISPR systems of widespread, uncultivated archaea can also target chromosomal DNA of archaeal episymbionts of the DPANN superphylum. Using meta-omics datasets from Crystal Geyser and Horonobe Underground Research Laboratory, we find that CRISPR spacers of the hosts Candidatus Altiarchaeum crystalense and Ca. A. horonobense, respectively, match putative essential genes in their episymbionts' genomes of the genus Ca. Huberiarchaeum and that some of these spacers are expressed in situ. Metabolic interaction modelling also reveals complementation between host-episymbiont systems, on the basis of which we propose that episymbionts are either parasitic or mutualistic depending on the genotype of the host. By expanding our analysis to 7,012 archaeal genomes, we suggest that CRISPR-Cas targeting of genomes associated with symbiotic archaea evolved independently in various archaeal lineages.

Experimental evolution of a facultative thermophile from a mesophilic ancestor.

Experimental evolution via continuous culture is a powerful approach to the alteration of complex phenotypes, such as optimal/maximal growth temperatures. The benefit of this approach is that phenotypic selection is tied to growth rate, allowing the production of optimized strains. Herein, we demonstrate the use of a recently described long-term culture apparatus called the Evolugator for the generation of a thermophilic descendant from a mesophilic ancestor (Escherichia coli MG1655). In addition, we used whole-genome sequencing of sequentially isolated strains throughout the thermal adaptation process to characterize the evolutionary history of the resultant genotype, identifying 31 genetic alterations that may contribute to thermotolerance, although some of these mutations may be adaptive for off-target environmental parameters, such as rich medium. We undertook preliminary phenotypic analysis of mutations identified in the glpF and fabA genes. Deletion of glpF in a mesophilic wild-type background conferred significantly improved growth rates in the 43-to-48°C temperature range and altered optimal growth temperature from 37°C to 43°C. In addition, transforming our evolved thermotolerant strain (EVG1064) with a wild-type allele of glpF reduced fitness at high temperatures. On the other hand, the mutation in fabA predictably increased the degree of saturation in membrane lipids, which is a known adaptation to elevated temperature. However, transforming EVG1064 with a wild-type fabA allele had only modest effects on fitness at intermediate temperatures. The Evolugator is fully automated and demonstrates the potential to accelerate the selection for complex traits by experimental evolution and significantly decrease development time for new industrial strains.

Zng1 is a GTP-dependent zinc transferase needed for activation of methionine aminopeptidase.

The evolution of zinc (Zn) as a protein cofactor altered the functional landscape of biology, but dependency on Zn also created an Achilles' heel, necessitating adaptive mechanisms to ensure Zn availability to proteins. A debated strategy is whether metallochaperones exist to prioritize essential Zn-dependent proteins. Here, we present evidence for a conserved family of putative metal transferases in human and fungi, which interact with Zn-dependent methionine aminopeptidase type I (MetAP1/Map1p/Fma1). Deletion of the putative metal transferase in Saccharomyces cerevisiae (ZNG1; formerly YNR029c) leads to defective Map1p function and a Zn-deficiency growth defect. In vitro, Zng1p can transfer Zn2+ or Co2+ to apo-Map1p, but unlike characterized copper chaperones, transfer is dependent on GTP hydrolysis. Proteomics reveal mis-regulation of the Zap1p transcription factor regulon because of loss of ZNG1 and Map1p activity, suggesting that Zng1p is required to avoid a compounding effect of Map1p dysfunction on survival during Zn limitation.

Structural characterization of a soil viral auxiliary metabolic gene product - a functional chitosanase.

Metagenomics is unearthing the previously hidden world of soil viruses. Many soil viral sequences in metagenomes contain putative auxiliary metabolic genes (AMGs) that are not associated with viral replication. Here, we establish that AMGs on soil viruses actually produce functional, active proteins. We focus on AMGs that potentially encode chitosanase enzymes that metabolize chitin - a common carbon polymer. We express and functionally screen several chitosanase genes identified from environmental metagenomes. One expressed protein showing endo-chitosanase activity (V-Csn) is crystalized and structurally characterized at ultra-high resolution, thus representing the structure of a soil viral AMG product. This structure provides details about the active site, and together with structure models determined using AlphaFold, facilitates understanding of substrate specificity and enzyme mechanism. Our findings support the hypothesis that soil viruses contribute auxiliary functions to their hosts.

Synergistic use of plant-prokaryote comparative genomics for functional annotations.

BACKGROUND: Identifying functions for all gene products in all sequenced organisms is a central challenge of the post-genomic era. However, at least 30-50% of the proteins encoded by any given genome are of unknown or vaguely known function, and a large number are wrongly annotated. Many of these 'unknown' proteins are common to prokaryotes and plants. We set out to predict and experimentally test the functions of such proteins. Our approach to functional prediction integrates comparative genomics based mainly on microbial genomes with functional genomic data from model microorganisms and post-genomic data from plants. This approach bridges the gap between automated homology-based annotations and the classical gene discovery efforts of experimentalists, and is more powerful than purely computational approaches to identifying gene-function associations. RESULTS: Among Arabidopsis genes, we focused on those (2,325 in total) that (i) are unique or belong to families with no more than three members, (ii) occur in prokaryotes, and (iii) have unknown or poorly known functions. Computer-assisted selection of promising targets for deeper analysis was based on homology-independent characteristics associated in the SEED database with the prokaryotic members of each family. In-depth comparative genomic analysis was performed for 360 top candidate families. From this pool, 78 families were connected to general areas of metabolism and, of these families, specific functional predictions were made for 41. Twenty-one predicted functions have been experimentally tested or are currently under investigation by our group in at least one prokaryotic organism (nine of them have been validated, four invalidated, and eight are in progress). Ten additional predictions have been independently validated by other groups. Discovering the function of very widespread but hitherto enigmatic proteins such as the YrdC or YgfZ families illustrates the power of our approach. CONCLUSIONS: Our approach correctly predicted functions for 19 uncharacterized protein families from plants and prokaryotes; none of these functions had previously been correctly predicted by computational methods. The resulting annotations could be propagated with confidence to over six thousand homologous proteins encoded in over 900 bacterial, archaeal, and eukaryotic genomes currently available in public databases.