Humanization of yeast genes with multiple human orthologs reveals functional divergence between paralogs

Functional Divergence in Orthologous Transcription Factors: Insights from AtCBF2/3/1 and OsDREB1C

Despite traditional beliefs of orthologous genes maintaining similar functions across species, growing evidence points to their potential for functional divergence. C-repeat binding factors/dehydration-responsive element binding protein 1s (CBFs/DREB1s) are critical in cold acclimation, with their overexpression enhancing stress tolerance but often constraining plant growth. In contrast, a recent study unveiled a distinctive role of rice OsDREB1C in elevating nitrogen use efficiency (NUE), photosynthesis, and grain yield, implying functional divergence within the CBF/DREB1 orthologs across species. Here, we delve into divergent molecular mechanisms of OsDREB1C and AtCBF2/3/1 by exploring their evolutionary trajectories across rice and Arabidopsis genomes, regulatomes, and transcriptomes. Evolutionary scrutiny shows discrete clades for OsDREB1C and AtCBF2/3/1, with the Poaceae-specific DREB1C clade mediated by a transposon event. Genome-wide binding profiles highlight OsDREB1C's preference for GCCGAC compared to AtCBF2/3/1's preference for A/GCCGAC, a distinction determined by R12 in the OsDREB1C AP2/ERF domain. Cross-species multiomic analyses reveal shared gene orthogroups (OGs) and underscore numerous specific OGs uniquely bound and regulated by OsDREB1C, implicated in NUE, photosynthesis, and early flowering, or by AtCBF2/3/1, engaged in hormone and stress responses. This divergence arises from gene gains/losses (∼16.7% to 25.6%) and expression reprogramming (∼62.3% to 66.2%) of OsDREB1C- and AtCBF2/3/1-regulated OGs during the extensive evolution following the rice-Arabidopsis split. Our findings illustrate the regulatory evolution of OsDREB1C and AtCBF2/3/1 at a genomic scale, providing insights on the functional divergence of orthologous transcription factors following gene duplications across species.

From beer to breadboards: yeast as a force for biological innovation

The history of yeast Saccharomyces cerevisiae, aka brewer's or baker's yeast, is intertwined with our own. Initially domesticated 8,000 years ago to provide sustenance to our ancestors, for the past 150 years, yeast has served as a model research subject and a platform for technology. In this review, we highlight many ways in which yeast has served to catalyze the fields of functional genomics, genome editing, gene-environment interaction investigation, proteomics, and bioinformatics-emphasizing how yeast has served as a catalyst for innovation. Several possible futures for this model organism in synthetic biology, drug personalization, and multi-omics research are also presented.

Parallel Nonfunctionalization of CK1δ/ε Kinase Ohnologs Following a Whole-Genome Duplication Event

Whole-genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in Saccharomyces cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post-WGD. In S. cerevisiae, HRR25 encodes the casein kinase 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post-WGD and non-WGD species, and have nonconserved cellular localization, consistent with their ongoing loss of function. The analysis in Naumovozyma castellii shows that the noncomplementing ohnolog is expressed at a lower level and has become nonessential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

Comparative Proteomic Analysis of Two Commonly Used Laboratory Yeast Strains: W303 and BY4742

The yeast Saccharomyces cerevisiae is a powerful model system that is often used to expand our understanding of cellular processes and biological functions. Although many genetically well-characterized laboratory strains of S. cerevisiae are available, they may have different genetic backgrounds which can confound data interpretation. Here, we report a comparative whole-proteome analysis of two common laboratory yeast background strains, W303 and BY4742, in both exponential and stationary growth phases using isobaric-tag-based mass spectrometry to highlight differences in proteome complexity. We quantified over 4400 proteins, hundreds of which showed differences in abundance between strains and/or growth phases. Moreover, we used proteome-wide protein abundance to profile the mating type of the strains used in the experiment, the auxotrophic markers, and associated metabolic pathways, as well as to investigate differences in particular classes of proteins, such as the pleiotropic drug resistance (PDR) proteins. This study is a valuable resource that offers insight into mechanistic differences between two common yeast background strains and can be used as a guide to select a background that is best suited for addressing a particular biological question.

Parallel nonfunctionalization of CK1δ/ε kinase ohnologs following a whole-genome duplication event

Whole genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in S. cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post WGD. In S. cerevisiae, HRR25 encodes the casein kinase (CK) 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post and non-WGD species, and have non-conserved cellular localization, consistent with their ongoing loss of function. The analysis in N. castelli shows that the non-complementing ohnolog is expressed at a lower level and has become non-essential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

Species-specific protein-protein interactions govern the humanization of the 20S proteasome in yeast

Yeast and humans share thousands of genes despite a billion years of evolutionary divergence. While many human genes can functionally replace their yeast counterparts, nearly half of the tested shared genes cannot. For example, most yeast proteasome subunits are "humanizable," except subunits comprising the β-ring core, including β2c (HsPSMB7, a constitutive proteasome subunit). We developed a high-throughput pipeline to humanize yeast proteasomes by generating a large library of Hsβ2c mutants and screening them for complementation of a yeast β2 (ScPup1) knockout. Variants capable of replacing ScPup1 included (1) those impacting local protein-protein interactions (PPIs), with most affecting interactions between the β2c C-terminal tail and the adjacent β3 subunit, and (2) those affecting β2c proteolytic activity. Exchanging the full-length tail of human β2c with that of ScPup1 enabled complementation. Moreover, wild-type human β2c could replace yeast β2 if human β3 was also provided. Unexpectedly, yeast proteasomes bearing a catalytically inactive HsPSMB7-T44A variant that blocked precursor autoprocessing were viable, suggesting an intact propeptide stabilizes late assembly intermediates. In contrast, similar modifications in human β2i (HsPSMB10), an immunoproteasome subunit and the co-ortholog of yeast β2, do not enable complementation in yeast, suggesting distinct interactions are involved in human immunoproteasome core assembly. Broadly, our data reveal roles for specific PPIs governing functional replaceability across vast evolutionary distances.

Yeast Deletomics to Uncover Gadolinium Toxicity Targets and Resistance Mechanisms

Among the rare earth elements (REEs), a crucial group of metals for high-technologies. Gadolinium (Gd) is the only REE intentionally injected to human patients. The use of Gd-based contrasting agents for magnetic resonance imaging (MRI) is the primary route for Gd direct exposure and accumulation in humans. Consequently, aquatic environments are increasingly exposed to Gd due to its excretion through the urinary tract of patients following an MRI examination. The increasing number of reports mentioning Gd toxicity, notably originating from medical applications of Gd, necessitates an improved risk-benefit assessment of Gd utilizations. To go beyond toxicological studies, unravelling the mechanistic impact of Gd on humans and the ecosystem requires the use of genome-wide approaches. We used functional deletomics, a robust method relying on the screening of a knock-out mutant library of Saccharomyces cerevisiae exposed to toxic concentrations of Gd. The analysis of Gd-resistant and -sensitive mutants highlighted the cell wall, endosomes and the vacuolar compartment as cellular hotspots involved in the Gd response. Furthermore, we identified endocytosis and vesicular trafficking pathways (ESCRT) as well as sphingolipids homeostasis as playing pivotal roles mediating Gd toxicity. Finally, tens of yeast genes with human orthologs linked to renal dysfunction were identified as Gd-responsive. Therefore, the molecular and cellular pathways involved in Gd toxicity and detoxification uncovered in this study underline the pleotropic consequences of the increasing exposure to this strategic metal.

Supervised biological network alignment with graph neural networks

MOTIVATION

Despite the advances in sequencing technology, massive proteins with known sequences remain functionally unannotated. Biological network alignment (NA), which aims to find the node correspondence between species' protein-protein interaction (PPI) networks, has been a popular strategy to uncover missing annotations by transferring functional knowledge across species. Traditional NA methods assumed that topologically similar proteins in PPIs are functionally similar. However, it was recently reported that functionally unrelated proteins can be as topologically similar as functionally related pairs, and a new data-driven or supervised NA paradigm has been proposed, which uses protein function data to discern which topological features correspond to functional relatedness.

RESULTS

Here, we propose GraNA, a deep learning framework for the supervised NA paradigm for the pairwise NA problem. Employing graph neural networks, GraNA utilizes within-network interactions and across-network anchor links for learning protein representations and predicting functional correspondence between across-species proteins. A major strength of GraNA is its flexibility to integrate multi-faceted non-functional relationship data, such as sequence similarity and ortholog relationships, as anchor links to guide the mapping of functionally related proteins across species. Evaluating GraNA on a benchmark dataset composed of several NA tasks between different pairs of species, we observed that GraNA accurately predicted the functional relatedness of proteins and robustly transferred functional annotations across species, outperforming a number of existing NA methods. When applied to a case study on a humanized yeast network, GraNA also successfully discovered functionally replaceable human-yeast protein pairs that were documented in previous studies.

AVAILABILITY AND IMPLEMENTATION

The code of GraNA is available at https://github.com/luo-group/GraNA.

Rapid, scalable, combinatorial genome engineering by marker-less enrichment and recombination of genetically engineered loci in yeast

A major challenge to rationally building multi-gene processes in yeast arises due to the combinatorics of combining all of the individual edits into the same strain. Here, we present a precise and multi-site genome editing approach that combines all edits without selection markers using CRISPR-Cas9. We demonstrate a highly efficient gene drive that selectively eliminates specific loci by integrating CRISPR-Cas9-mediated double-strand break (DSB) generation and homology-directed recombination with yeast sexual assortment. The method enables marker-less enrichment and recombination of genetically engineered loci (MERGE). We show that MERGE converts single heterologous loci to homozygous loci at ∼100% efficiency, independent of chromosomal location. Furthermore, MERGE is equally efficient at converting and combining multiple loci, thus identifying compatible genotypes. Finally, we establish MERGE proficiency by engineering a fungal carotenoid biosynthesis pathway and most of the human α-proteasome core into yeast. Therefore, MERGE lays the foundation for scalable, combinatorial genome editing in yeast.

DASH/Dam1 complex mutants stabilize ploidy in histone-humanized yeast by weakening kinetochore-microtubule attachments

Forcing budding yeast to chromatinize their DNA with human histones manifests an abrupt fitness cost. We previously proposed chromosomal aneuploidy and missense mutations as two potential modes of adaptation to histone humanization. Here, we show that aneuploidy in histone-humanized yeasts is specific to a subset of chromosomes that are defined by their centromeric evolutionary origins but that these aneuploidies are not adaptive. Instead, we find that a set of missense mutations in outer kinetochore proteins drives adaptation to human histones. Furthermore, we characterize the molecular mechanism underlying adaptation in two mutants of the outer kinetochore DASH/Dam1 complex, which reduce aneuploidy by suppression of chromosome instability. Molecular modeling and biochemical experiments show that these two mutants likely disrupt a conserved oligomerization interface thereby weakening microtubule attachments. We propose a model through which weakened microtubule attachments promote increased kinetochore-microtubule turnover and thus suppress chromosome instability. In sum, our data show how a set of point mutations evolved in histone-humanized yeasts to counterbalance human histone-induced chromosomal instability through weakening microtubule interactions, eventually promoting a return to euploidy.

A new combinatorial megaplasmid library assembly method designed to screen for minimal pathways by using SCRaMbLE

Human proteins expressed in yeast are common to enhance protein production while the expression of functional human pathways remain challenging. Here, we propose a simple and economical high-throughput gene assembly method to create a yeast megaplasmid library from human cDNA to screen for minimal human functional pathways. We introduced artificial promoters followed by symmetric loxP sites into the megaplasmids using Golden Gate assembly coupled with streptavidin-bead-based purification. The isolated high molecular weight, randomly assembled cDNA megaplasmid library may be useful for high-throughput directed evolution experiments and may be adapted for use in other model organisms.

The First Yarrowia lipolytica Yeast Models Expressing Hepatitis B Virus X Protein: Changes in Mitochondrial Morphology and Functions

Chronic hepatitis B virus infection is the dominant cause of hepatocellular carcinoma, the main cause of cancer death. HBx protein, a multifunctional protein, is essential for pathogenesis development; however, the underlying mechanisms are not fully understood. The complexity of the system itself, and the intricate interplay of many factors make it difficult to advance in understanding the mechanisms underlying these processes. The most obvious solution is to use simpler systems by reducing the number of interacting factors. Yeast cells are particularly suitable for studying the relationships between oxidative stress, mitochondrial dynamics (mitochondrial fusion and fragmentation), and mitochondrial dysfunction involved in HBx-mediated pathogenesis. For the first time, genetically modified yeast, Y. lipolytica, was created, expressing the hepatitis B virus core protein HBx, as well as a variant fused with eGFP at the C-end. It was found that cells expressing HBx experienced stronger oxidative stress than the control cells. Oxidative stress was alleviated by preincubation with the mitochondria-targeted antioxidant SkQThy. Consistent with these data, in contrast to the control cells (pZ-0) containing numerous mitochondrial forming a mitochondrial reticulum, in cells expressing HBx protein, mitochondria were fragmented, and preincubation with SkQThy partially restored the mitochondrial reticulum. Expression of HBx had a significant influence on the bioenergetic function of mitochondria, making them loosely coupled with decreased respiratory rate and reduced ATP formation. In sum, the first highly promising yeast model for studying the impact of HBx on bioenergy, redox-state, and dynamics of mitochondria in the cell and cross-talk between these parameters was offered. This fairly simple model can be used as a platform for rapid screening of potential therapeutic agents, mitigating the harmful effects of HBx.

Evidence that conserved essential genes are enriched for pro-longevity factors

At the cellular level, many aspects of aging are conserved across species. This has been demonstrated by numerous studies in simple model organisms like Saccharomyces cerevisiae, Caenorhabdits elegans, and Drosophila melanogaster. Because most genetic screens examine loss of function mutations or decreased expression of genes through reverse genetics, essential genes have often been overlooked as potential modulators of the aging process. By taking the approach of increasing the expression level of a subset of conserved essential genes, we found that 21% of these genes resulted in increased replicative lifespan in S. cerevisiae. This is greater than the ~ 3.5% of genes found to affect lifespan upon deletion, suggesting that activation of essential genes may have a relatively disproportionate effect on increasing lifespan. The results of our experiments demonstrate that essential gene overexpression is a rich, relatively unexplored means of increasing eukaryotic lifespan.

Full humanization of the glycolytic pathway in Saccharomyces cerevisiae

Although transplantation of single genes in yeast plays a key role in elucidating gene functionality in metazoans, technical challenges hamper humanization of full pathways and processes. Empowered by advances in synthetic biology, this study demonstrates the feasibility and implementation of full humanization of glycolysis in yeast. Single gene and full pathway transplantation revealed the remarkable conservation of glycolytic and moonlighting functions and, combined with evolutionary strategies, brought to light context-dependent responses. Human hexokinase 1 and 2, but not 4, required mutations in their catalytic or allosteric sites for functionality in yeast, whereas hexokinase 3 was unable to complement its yeast ortholog. Comparison with human tissues cultures showed preservation of turnover numbers of human glycolytic enzymes in yeast and human cell cultures. This demonstration of transplantation of an entire essential pathway paves the way for establishment of species-, tissue-, and disease-specific metazoan models.

Humanized yeast to model human biology, disease and evolution

For decades, budding yeast, a single-cellular eukaryote, has provided remarkable insights into human biology. Yeast and humans share several thousand genes despite morphological and cellular differences and over a billion years of separate evolution. These genes encode critical cellular processes, the failure of which in humans results in disease. Although recent developments in genome engineering of mammalian cells permit genetic assays in human cell lines, there is still a need to develop biological reagents to study human disease variants in a high-throughput manner. Many protein-coding human genes can successfully substitute for their yeast equivalents and sustain yeast growth, thus opening up doors for developing direct assays of human gene function in a tractable system referred to as 'humanized yeast'. Humanized yeast permits the discovery of new human biology by measuring human protein activity in a simplified organismal context. This Review summarizes recent developments showing how humanized yeast can directly assay human gene function and explore variant effects at scale. Thus, by extending the 'awesome power of yeast genetics' to study human biology, humanizing yeast reinforces the high relevance of evolutionarily distant model organisms to explore human gene evolution, function and disease.

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

Coexpression reveals conserved gene programs that co-vary with cell type across kingdoms

What makes a mouse a mouse, and not a hamster? Differences in gene regulation between the two organisms play a critical role. Comparative analysis of gene coexpression networks provides a general framework for investigating the evolution of gene regulation across species. Here, we compare coexpression networks from 37 species and quantify the conservation of gene activity 1) as a function of evolutionary time, 2) across orthology prediction algorithms, and 3) with reference to cell- and tissue-specificity. We find that ancient genes are expressed in multiple cell types and have well conserved coexpression patterns, however they are expressed at different levels across cell types. Thus, differential regulation of ancient gene programs contributes to transcriptional cell identity. We propose that this differential regulation may play a role in cell diversification in both the animal and plant kingdoms.

Benzylic Dehydroxylation of Echinocandin Antifungal Drugs Restores Efficacy against Resistance Conferred by Mutated Glucan Synthase

Each year, infections caused by fungal pathogens claim the lives of about 1.6 million people and affect the health of over a billion people worldwide. Among the most recently developed antifungal drugs are the echinocandins, which noncompetitively inhibit β-glucan synthase, a membrane-bound protein complex that catalyzes the formation of the main polysaccharide component of the fungal cell wall. Resistance to echinocandins is conferred by mutations in FKS genes, which encode the catalytic subunit of the β-glucan synthase complex. Here, we report that selective removal of the benzylic alcohol of the nonproteinogenic amino acid 3S,4S-dihydroxy-l-homotyrosine of the echinocandins anidulafungin and rezafungin, restored their efficacy against a large panel of echinocandin-resistant Candida strains. The dehydroxylated compounds did not significantly affect the viability of human-derived cell culture lines. An analysis of the efficacy of the dehydroxylated echinocandins against resistant Candida strains, which contain mutations in the FKS1 and/or FKS2 genes of the parental strains, identified amino acids of the Fks proteins that are likely to reside in proximity to the l-homotyrosine residue of the bound drug. This study describes the first example of a chemical modification strategy to restore the efficacy of echinocandin drugs, which have a critical place in the arsenal of antifungal drugs, against resistant fungal pathogens.

Fluconazole-COX Inhibitor Hybrids: A Dual-Acting Class of Antifungal Azoles

When used in combination with azole antifungal drugs, cyclooxygenase (COX) inhibitors such as ibuprofen improve antifungal efficacy. We report the conjugation of a chiral antifungal azole pharmacophore to COX inhibitors and the evaluation of activity of 24 hybrids. Hybrids derived from ibuprofen and flurbiprofen were considerably more potent than fluconazole and comparable to voriconazole against a panel of Candida species. The potencies of hybrids composed of an S-configured azole pharmacophore were higher than those with an R-configured pharmacophore. Tolerance, defined as the ability of a subpopulation of cells to grow in the presence of the drug, to the hybrids was lower than to fluconazole and voriconazole. The hybrids were active against a mutant lacking CYP51, the target of azole drugs, indicating that these agents act via a dual mode of action. This study established that azole-COX inhibitor hybrids are a novel class of potent antifungals with clinical potential.

Retention of duplicated genes in evolution

Gene duplication is a prevalent phenomenon across the tree of life. The processes that lead to the retention of duplicated genes are not well understood. Functional genomics approaches in model organisms, such as yeast, provide useful tools to test the mechanisms underlying retention with functional redundancy and divergence of duplicated genes, including fates associated with neofunctionalization, subfunctionalization, back-up compensation, and dosage amplification. Duplicated genes may also be retained as a consequence of structural and functional entanglement. Advances in human gene editing have enabled the interrogation of duplicated genes in the human genome, providing new tools to evaluate the relative contributions of each of these factors to duplicate gene retention and the evolution of genome structure.

Sequence, structural and functional conservation among the human and fission yeast ELL and EAF transcription elongation factors

BACKGROUND

Transcription elongation is a dynamic and tightly regulated step of gene expression in eukaryotic cells. Eleven nineteen Lysine rich Leukemia (ELL) and ELL Associated Factors (EAF) family of conserved proteins are required for efficient RNA polymerase II-mediated transcription elongation. Orthologs of these proteins have been identified in different organisms, including fission yeast and humans.

METHODS AND RESULTS

In the present study, we have examined the sequence, structural and functional conservation between the fission yeast and human ELL and EAF orthologs. Our computational analysis revealed that these proteins share some sequence characteristics, and were predominantly disordered in both organisms. Our functional complementation assays revealed that both human ELL and EAF proteins could complement the lack of ell1⁺ or eaf1⁺ in Schizosaccharomyces pombe respectively. Furthermore, our domain mapping experiments demonstrated that both the amino and carboxyl terminal domains of human EAF proteins could functionally complement the S. pombe eaf1 deletion phenotypes. However, only the carboxyl-terminus domain of human ELL was able to partially rescue the phenotypes associated with lack of ell1⁺ in S. pombe.

CONCLUSIONS

Collectively, our work adds ELL-EAF to the increasing list of human-yeast complementation gene pairs, wherein the simpler fission yeast can be used to further enhance our understanding of the role of these proteins in transcription elongation and human disease.

Ancestral sequences of a large promiscuous enzyme family correspond to bridges in sequence space in a network representation

Evolutionary relationships of protein families can be characterized either by networks or by trees. Whereas trees allow for hierarchical grouping and reconstruction of the most likely ancestral sequences, networks lack a time axis but allow for thresholds of pairwise sequence identity to be chosen and, therefore, the clustering of family members with presumably more similar functions. Here, we use the large family of arylsulfatases and phosphonate monoester hydrolases to investigate similarities, strengths and weaknesses in tree and network representations. For varying thresholds of pairwise sequence identity, values of betweenness centrality and clustering coefficients were derived for nodes of the reconstructed ancestors to measure the propensity to act as a bridge in a network. Based on these properties, ancestral protein sequences emerge as bridges in protein sequence networks. Interestingly, many ancestral protein sequences appear close to extant sequences. Therefore, reconstructed ancestor sequences might also be interpreted as yet-to-be-identified homologues. The concept of ancestor reconstruction is compared to consensus sequences, too. It was found that hub sequences in a network, e.g. reconstructed ancestral sequences that are connected to many neighbouring sequences, share closer similarity with derived consensus sequences. Therefore, some reconstructed ancestor sequences can also be interpreted as consensus sequences.

Reconstruction of human genome evolution in yeast: an educational primer for use with "systematic humanization of the yeast cytoskeleton discerns functionally replaceable from divergent human genes"

The evolution of eukaryotic organisms starting with the last eukaryotic common ancestor was accompanied by lineage-specific expansion of gene families. A paper by Garge et al. provides an excellent opportunity to have students explore how expansion of gene families via gene duplication results in protein specialization, in this case in the context of eukaryotic cytoskeletal organization . The authors tested hypotheses about conserved protein function by systematic "humanization" of the yeast cytoskeletal components while employing a wide variety of methodological approaches. We outline several exercises to promote students' ability to explore the genomic databases, perform bioinformatic analyses, design experiments for functional analysis of human genes in yeast and critically interpret results to address both specific and general questions.

Measuring Pharmacogene Variant Function at Scale Using Multiplexed Assays

As costs of next-generation sequencing decrease, identification of genetic variants has far outpaced our ability to understand their functional consequences. This lack of understanding is a central challenge to a key promise of pharmacogenomics: using genetic information to guide drug selection and dosing. Recently developed multiplexed assays of variant effect enable experimental measurement of the function of thousands of variants simultaneously. Here, we describe multiplexed assays that have been performed on nearly 25,000 variants in eight key pharmacogenes (ADRB2, CYP2C9, CYP2C19, NUDT15, SLCO1B1, TMPT, VKORC1, and the LDLR promoter), discuss advances in experimental design, and explore key challenges that must be overcome to maximize the utility of multiplexed functional data. Expected final online publication date for the Annual Review of Pharmacology and Toxicology, Volume 62 is January 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Discovery of new vascular disrupting agents based on evolutionarily conserved drug action, pesticide resistance mutations, and humanized yeast

Thiabendazole (TBZ) is an FDA-approved benzimidazole widely used for its antifungal and antihelminthic properties. We showed previously that TBZ is also a potent vascular disrupting agent and inhibits angiogenesis at the tissue level by dissociating vascular endothelial cells in newly formed blood vessels. Here, we uncover TBZ's molecular target and mechanism of action. Using human cell culture, molecular modeling, and humanized yeast, we find that TBZ selectively targets only 1 of 9 human β-tubulin isotypes (TUBB8) to specifically disrupt endothelial cell microtubules. By leveraging epidemiological pesticide resistance data and mining chemical features of commercially used benzimidazoles, we discover that a broader class of benzimidazole compounds, in extensive use for 50 years, also potently disrupt immature blood vessels and inhibit angiogenesis. Thus, besides identifying the molecular mechanism of benzimidazole-mediated vascular disruption, this study presents evidence relevant to the widespread use of these compounds while offering potential new clinical applications.

Special Care Is Needed in Applying Phylogenetic Comparative Methods to Gene Trees with Speciation and Duplication Nodes

How gene function evolves is a central question of evolutionary biology. It can be investigated by comparing functional genomics results between species and between genes. Most comparative studies of functional genomics have used pairwise comparisons. Yet it has been shown that this can provide biased results, as genes, like species, are phylogenetically related. Phylogenetic comparative methods should be used to correct for this, but they depend on strong assumptions, including unbiased tree estimates relative to the hypothesis being tested. Such methods have recently been used to test the “ortholog conjecture,” the hypothesis that functional evolution is faster in paralogs than in orthologs. Although pairwise comparisons of tissue specificity (τ) provided support for the ortholog conjecture, phylogenetic independent contrasts did not. Our reanalysis on the same gene trees identified problems with the time calibration of duplication nodes. We find that the gene trees used suffer from important biases, due to the inclusion of trees with no duplication nodes, to the relative age of speciations and duplications, to systematic differences in branch lengths, and to non-Brownian motion of tissue specificity on many trees. We find that incorrect implementation of phylogenetic method in empirical gene trees with duplications can be problematic. Controlling for biases allows successful use of phylogenetic methods to study the evolution of gene function and provides some support for the ortholog conjecture using three different phylogenetic approaches.

Genome-wide toxicogenomic study of the lanthanides sheds light on the selective toxicity mechanisms associated with critical materials

Lanthanides are a series of critical elements widely used in multiple industries, such as optoelectronics and healthcare. Although initially considered to be of low toxicity, concerns have emerged during the last few decades over their impact on human health. The toxicological profile of these metals, however, has been incompletely characterized, with most studies to date solely focusing on one or two elements within the group. In the current study, we assessed potential toxicity mechanisms in the lanthanide series using a functional toxicogenomics approach in baker's yeast, which shares many cellular pathways and functions with humans. We screened the homozygous deletion pool of 4,291 Saccharomyces cerevisiae strains with the lanthanides and identified both common and unique functional effects of these metals. Three very different trends were observed within the lanthanide series, where deletions of certain proteins on membranes and organelles had no effect on the cellular response to early lanthanides while inducing yeast sensitivity and resistance to middle and late lanthanides, respectively. Vesicle-mediated transport (primarily endocytosis) was highlighted by both gene ontology and pathway enrichment analyses as one of the main functions disturbed by the majority of the metals. Protein-protein network analysis indicated that yeast response to lanthanides relied on proteins that participate in regulatory paths used for calcium (and other biologically relevant cations), and lanthanide toxicity included disruption of biosynthetic pathways by enzyme inhibition. Last, multiple genes and proteins identified in the network analysis have human orthologs, suggesting that those may also be targeted by lanthanides in humans.

Propagation of Mitochondria-Derived Reactive Oxygen Species within the Dipodascus magnusii Cells

Mitochondria are considered to be the main source of reactive oxygen species (ROS) in the cell. It was shown that in cardiac myocytes exposed to excessive oxidative stress, ROS-induced ROS release is triggered. However, cardiac myocytes have a network of densely packed organelles that do not move, which is not typical for the majority of eukaryotic cells. The purpose of this study was to trace the spatiotemporal development (propagation) of prooxidant-induced oxidative stress and its interplay with mitochondrial dynamics. We used Dipodascus magnusii yeast cells as a model, as they have advantages over other models, including a uniquely large size, mitochondria that are easy to visualize and freely moving, an ability to vigorously grow on well-defined low-cost substrates, and high responsibility. It was shown that prooxidant-induced oxidative stress was initiated in mitochondria, far preceding the appearance of generalized oxidative stress in the whole cell. For yeasts, these findings were obtained for the first time. Preincubation of yeast cells with SkQ1, a mitochondria-addressed antioxidant, substantially diminished production of mitochondrial ROS, while only slightly alleviating the generalized oxidative stress. This was expected, but had not yet been shown. Importantly, mitochondrial fragmentation was found to be primarily induced by mitochondrial ROS preceding the generalized oxidative stress development.

The ortholog conjecture revisited: the value of orthologs and paralogs in function prediction

MOTIVATION

The computational prediction of gene function is a key step in making full use of newly sequenced genomes. Function is generally predicted by transferring annotations from homologous genes or proteins for which experimental evidence exists. The 'ortholog conjecture' proposes that orthologous genes should be preferred when making such predictions, as they evolve functions more slowly than paralogous genes. Previous research has provided little support for the ortholog conjecture, though the incomplete nature of the data cast doubt on the conclusions.

RESULTS

We use experimental annotations from over 40 000 proteins, drawn from over 80 000 publications, to revisit the ortholog conjecture in two pairs of species: (i) Homo sapiens and Mus musculus and (ii) Saccharomyces cerevisiae and Schizosaccharomyces pombe. By making a distinction between questions about the evolution of function versus questions about the prediction of function, we find strong evidence against the ortholog conjecture in the context of function prediction, though questions about the evolution of function remain difficult to address. In both pairs of species, we quantify the amount of information that would be ignored if paralogs are discarded, as well as the resulting loss in prediction accuracy. Taken as a whole, our results support the view that the types of homologs used for function transfer are largely irrelevant to the task of function prediction. Maximizing the amount of data used for this task, regardless of whether it comes from orthologs or paralogs, is most likely to lead to higher prediction accuracy.

AVAILABILITY AND IMPLEMENTATION

https://github.com/predragradivojac/oc.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.