Comparison of the complete protein sets of worm and yeast: orthology and divergence

Construction of protein-protein interaction network in sulfate-reducing bacteria: Unveiling of global response to Hg

Sulfate-reducing bacteria (SRB) play pivotal roles in the biotransformation of mercury (Hg). However, unrevealed global responses of SRB to Hg have restricted our understanding of details of Hg biotransformation processes. The absence of protein-protein interaction (PPI) network under Hg stimuli has been a bottleneck of proteomic analysis for molecular mechanisms of Hg transformation. This study constructed the first comprehensive PPI network of SRB in response to Hg, encompassing 67 connected nodes, 26 independent nodes, and 121 edges, covering 93% of differentially expressed proteins from both previous studies and this study. The network suggested that proteomic changes of SRB in response to Hg occurred globally, including microbial metabolism in diverse environments, carbon metabolism, nucleic acid metabolism and translation, nucleic acid repair, transport systems, nitrogen metabolism, and methyltransferase activity, partial of which could cover the known knowledge. Antibiotic resistance was the original response revealed by this network, providing insights into of Hg biotransformation mechanisms. This study firstly provided the foundational network for a comprehensive understanding of SRB's responses to Hg, convenient for exploration of potential targets for Hg biotransformation. Furthermore, the network indicated that Hg enhances the metabolic activities and modification pathways of SRB to maintain cellular activities, shedding light on the influences of Hg on the carbon, nitrogen, and sulfur cycles at the cellular level.

AtGCS promoter-driven clustered regularly interspaced short palindromic repeats/Cas9 highly efficiently generates homozygous/biallelic mutations in the transformed roots by Agrobacterium rhizogenes-mediated transformation

Agrobacterium rhizogenes-mediated (ARM) transformation is an efficient and powerful tool to generate transgenic roots to study root-related biology. For loss-of-function studies, transgenic-root-induced indel mutations by CRISPR/Cas9 only with homozygous/biallelic mutagenesis can exhibit mutant phenotype(s) (excluding recessive traits). However, a low frequency of homozygous mutants was produced by a constitutive promoter to drive Cas9 expression. Here, we identified a highly efficient Arabidopsis thaliana gamma-glutamylcysteine synthetase promoter, termed AtGCSpro, with strong activity in the region where the root meristem will initiate and in the whole roots in broad eudicots species. AtGCSpro achieved higher homozygous/biallelic mutation efficiency than the most widely used CaMV 35S promoter in driving Cas9 expression in soybean, Lotus japonicus, and tomato roots. Using the pAtGCSpro-Cas9 system, the average homozygous/biallelic mutation frequency is 1.7-fold and 8.3-fold higher than the p2 × 35Spro-Cas9 system for single and two target site(s) in the genome, respectively. Our results demonstrate the advantage of the pAtGCSpro-Cas9 system used in ARM transformation, especially its great potential in diploids with multiple-copy genes targeted mutations and polyploid plants with multiplex genome editing. AtGCSpro is conservatively active in various eudicots species, suggesting that AtGCSpro might be applied in a wide range of dicots species.

Defining characteristics and conservation of poorly annotated genes in Caenorhabditis elegans using WormCat 2.0

Omics tools provide broad datasets for biological discovery. However, the computational tools for identifying important genes or pathways in RNA-seq, proteomics, or GWAS (Genome-Wide Association Study) data depend on Gene Ontogeny annotations and are biased toward well-described pathways. This limits their utility as poorly annotated genes, which could have novel functions, are often passed over. Recently, we developed an annotation and category enrichment tool for Caenorhabditis elegans genomic data, WormCat, which provides an intuitive visualization output. Unlike Gene Ontogeny-based enrichment tools, which exclude genes with no annotation information, WormCat 2.0 retains these genes as a special UNASSIGNED category. Here, we show that the UNASSIGNED gene category enrichment exhibits tissue-specific expression patterns and can include genes with biological functions identified in published datasets. Poorly annotated genes are often considered to be potentially species-specific and thus, of reduced interest to the biomedical community. Instead, we find that around 3% of the UNASSIGNED genes have human orthologs, including some linked to human diseases. These human orthologs themselves have little annotation information. A recently developed method that incorporates lineage relationships (abSENSE) indicates that the failure of BLAST to detect homology explains the apparent lineage specificity for many UNASSIGNED genes. This suggests that a larger subset could be related to human genes. WormCat provides an annotation strategy that allows the association of UNASSIGNED genes with specific phenotypes and known pathways. Building these associations in C. elegans, with its robust genetic tools, provides a path to further functional study and insight into these understudied genes.

OUP accepted manuscript

Extensive amounts of data from next-generation sequencing and omics studies have led to the accumulation of information that provides insight into the evolutionary landscape of related proteins. Here, we present OrthoQuantum, a web server that allows for time-efficient analysis and visualization of phylogenetic profiles of any set of eukaryotic proteins. It is a simple-to-use tool capable of searching large input sets of proteins. Using data from open source databases of orthologous sequences in a wide range of taxonomic groups, it enables users to assess coupled evolutionary patterns and helps define lineage-specific innovations. The web interface allows to perform queries with gene names and UniProt identifiers in different phylogenetic clades and supplement presence with an additional BLAST search. The conservation patterns of proteins are coded as binary vectors, i.e., strings that encode the presence or absence of orthologous proteins in other genomes. These strings are used to calculate top-scoring correlation pairs needed for finding co-inherited proteins which are simultaneously present or simultaneously absent in specific lineages. Profiles are visualized in combination with phylogenetic trees in a JavaScript-based interface. The OrthoQuantum v1.0 web server is freely available at http://orthoq.bioinf.fbb.msu.ru along with documentation and tutorial.

Transcriptional Profiling of Three Pseudomonas syringae pv. actinidiae Biovars Reveals Different Responses to Apoplast-Like Conditions Related to Strain Virulence on the Host

Pseudomonas syringae pv. actinidiae is a phytopathogen that causes devastating bacterial canker in kiwifruit. Among five biovars defined by genetic, biochemical, and virulence traits, P. syringae pv. actinidiae biovar 3 (Psa3) is the most aggressive and is responsible for the most recent reported outbreaks; however, the molecular basis of its heightened virulence is unclear. Therefore, we designed the first P. syringae multistrain whole-genome microarray, encompassing biovars Psa1, Psa2, and Psa3 and the well-established model P. syringae pv. tomato, and analyzed early bacterial responses to an apoplast-like minimal medium. Transcriptomic profiling revealed i) the strong activation in Psa3 of all hypersensitive reaction and pathogenicity (hrp) and hrp conserved (hrc) cluster genes, encoding components of the type III secretion system required for bacterial pathogenicity and involved in responses to environmental signals; ii) potential repression of the hrp/hrc cluster in Psa2; and iii) activation of flagellum-dependent cell motility and chemotaxis genes in Psa1. The detailed investigation of three gene families encoding upstream regulatory proteins (histidine kinases, their cognate response regulators, and proteins with diguanylate cyclase or phosphodiesterase domains) indicated that cyclic di-GMP may be a key regulator of virulence in P. syringae pv. actinidiae biovars. The gene expression data were supported by the quantification of biofilm formation. Our findings suggest that diverse early responses to the host apoplast, even among bacteria belonging to the same pathovar, can lead to different virulence strategies and may explain the differing outcomes of infections. Based on our detailed structural analysis of hrp operons, we also propose a revision of hrp cluster organization and operon regulation in P. syringae.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Finding the generalized molecular principles of protein thermal stability

Are there any generalized molecular principles of thermal adaptation? Here, integrating the concepts of structural bioinformatics, sequence analysis, and classical knot theory, we develop a robust computational framework that seeks for mechanisms of thermal adaptation by comparing orthologous mesophilic-thermophilic and mesophilic-hyperthermophilic proteins of remarkable structural and topological similarities, and still leads us to context-independent results. A comprehensive analysis of 4741 high-resolution, non-redundant X-ray crystallographic structures collected from 11 hyperthermophilic, 32 thermophilic and 53 mesophilic prokaryotes unravels at least five "nearly universal" signatures of thermal adaptation, irrespective of the enormous sequence, structure, and functional diversity of the proteins compared. A careful investigation further extracts a set of amino acid changes that can potentially enhance protein thermal stability, and remarkably, these mutations are overrepresented in protein crystallization experiments, in disorder-to-order transitions and in engineered thermostable variants of existing mesophilic proteins. These results could be helpful to find a precise, global picture of thermal adaptation.

The genome of a subterrestrial nematode reveals adaptations to heat

The nematode Halicephalobus mephisto was originally discovered inhabiting a deep terrestrial aquifer 1.3 km underground. H. mephisto can thrive under conditions of abiotic stress including heat and minimal oxygen, where it feeds on a community of both chemolithotrophic and heterotrophic prokaryotes in an unusual ecosystem isolated from the surface biosphere. Here we report the comprehensive genome and transcriptome of this organism, identifying a signature of adaptation: an expanded repertoire of 70 kilodalton heat-shock proteins (Hsp70) and avrRpt2 induced gene 1 (AIG1) proteins. The expanded Hsp70 genes are transcriptionally induced upon growth under heat stress, and we find that positive selection is detectable in several members of this family. We further show that AIG1 may have been acquired by horizontal gene transfer (HGT) from a rhizobial fungus. Over one-third of the genes of H. mephisto are novel, highlighting the divergence of this nematode from other sequenced organisms. This work sheds light on the genomic basis of heat tolerance in a complete subterrestrial eukaryotic genome.

The nematode Halicephalobus mephisto has been found more than 1 km underground and can tolerate high heat and low oxygen. Here Weinstein et al. report the genome and transcriptome of H. mephisto and identify genomic innovations that likely are responsible for its tolerance to heat.

Data on taxonomic annotation and diversity of 18S rRNA gene amplicon libraries derived from high throughput sequencing

This Data in Brief article is a supporting information for the research article entitled “Protistan community composition in anoxic sediments from three salinity-disparate Japanese lakes” by Kataoka and Kondo (2019) [1]. Summary of 18S rRNA gene sequences originated from anoxic sediment of three lakes in two seasons using high throughput sequencing techniques (MiSeq, Illumina) was shown in this data article. Supergroup-level taxonomy was compared between the SILVA search for SILVA database and BLASTn search for the PR2 database. Alpha diversity was calculated in each sample, and beta-diversity was calculated among the six amplicon libraries. Partial sequence length between the primer set of 574*f and 1132R Hugerth et al., 2015 was compared between the forward read and the combined read.

Structure of proteins: Evolution with unsolved mysteries

Evolution of macromolecules could be considered as a milestone in the history of life. Nucleic acids are the long stretches of nucleotides that contain all the possible codes and information of life. On the other hand, proteins are their actual translated outcomes, or reflections of modifications in their structure that have occurred at a slow, but steady rate over a very long period of evolution. Over the years of research, biophysicists, biochemists, molecular and structural biologists have unfurled several layers of the structural convolutions in these chemical molecules; however evolutionists look over their structures through a different prism, which may or may not coincide with others. There remains a need to outline several well-known, but less discussed features of protein structures, like intrinsically disordered states, degron signals and different types of ubiquitin chains providing degradation signals, which help the cellular proteolytic machinery to identify and target the proteins towards degradation pathways. There are several important factors, which are critical for folding of proteins into their native three-dimensional conformations by the cytoplasmic chaperones; but in real time how the chaperones fold the newly synthesized polypeptide sequences into a particular three-dimensional shape within a fraction of second is still a mystery for biologists as well as mathematicians. Multiple similar unsolved or unaddressed questions need to be addressed in detail so that future line of research can dig deeper into the finer details of these structures of the proteins.

The Prohormone Proinsulin as a Neuroprotective Factor: Past History and Future Prospects

Proinsulin was first identified as the primary translation product of the insulin gene in Donald Steiner’s laboratory in 1967, and was the first prohormone to be isolated and sequenced. While its role as an insulin precursor has been extensively studied in the field of endocrinology, the bioactivity of the proinsulin molecule itself has received much less attention. Insulin binds to isoforms A and B of the insulin receptor (IR) with high affinity. Proinsulin, in contrast, binds with high affinity only to IR-A, which is present in the nervous system, among other tissues and elicits antiapoptotic and neuroprotective effects in the developing and postnatal nervous system. Proinsulin specifically exerts neuroprotection in the degenerating retina in mouse and rat models of retinitis pigmentosa (RP), delaying photoreceptor and vision loss after local administration in the eye or systemic (intramuscular) administration of an adeno-associated viral (AAV) vector that induces constitutive proinsulin release. AAV-mediated proinsulin expression also decreases the expression of neuroinflammation markers in the hippocampus and sustains cognitive performance in a mouse model of precocious brain senescence. We have therefore proposed that proinsulin should be considered a functionally distinct member of the insulin superfamily. Here, we briefly review the legacy of Steiner’s research, the neural expression of proinsulin, and the tissue expression patterns and functional characteristics of IR-A. We discuss the neuroprotective activity of proinsulin and its potential as a therapeutic tool in neurodegenerative conditions of the central nervous system, particularly in retinal dystrophies.

Genome-wide identification and evolution of the PIN-FORMED (PIN) gene family in Glycine max

Soybean (Glycine max) is one of the most important crop plants. Wild and cultivated soybean varieties have significant differences worth further investigation, such as plant morphology, seed size, and seed coat development; these characters may be related to auxin biology. The PIN gene family encodes essential transport proteins in cell-to-cell auxin transport, but little research on soybean PIN genes (GmPIN genes) has been done, especially with respect to the evolution and differences between wild and cultivated soybean. In this study, we retrieved 23 GmPIN genes from the latest updated G. max genome database; six GmPIN protein sequences were changed compared with the previous database. Based on the Plant Genome Duplication Database, 18 GmPIN genes have been involved in segment duplication. Three pairs of GmPIN genes arose after the second soybean genome duplication, and six occurred after the first genome duplication. The duplicated GmPIN genes retained similar expression patterns. All the duplicated GmPIN genes experienced purifying selection (Ka/Ks < 1) to prevent accumulation of non-synonymous mutations and thus remained more similar. In addition, we also focused on the artificial selection of the soybean PIN genes. Five artificially selected GmPIN genes were identified by comparing the genome sequence of 17 wild and 14 cultivated soybean varieties. Our research provides useful and comprehensive basic information for understanding GmPIN genes.

Network-assisted target identification for haploinsufficiency and homozygous profiling screens

Chemical genomic screens have recently emerged as a systematic approach to drug discovery on a genome-wide scale. Drug target identification and elucidation of the mechanism of action (MoA) of hits from these noisy high-throughput screens remain difficult. Here, we present GIT (Genetic Interaction Network-Assisted Target Identification), a network analysis method for drug target identification in haploinsufficiency profiling (HIP) and homozygous profiling (HOP) screens. With the drug-induced phenotypic fitness defect of the deletion of a gene, GIT also incorporates the fitness defects of the gene's neighbors in the genetic interaction network. On three genome-scale yeast chemical genomic screens, GIT substantially outperforms previous scoring methods on target identification on HIP and HOP assays, respectively. Finally, we showed that by combining HIP and HOP assays, GIT further boosts target identification and reveals potential drug's mechanism of action.

Reclassification of the taxonomic status of SEMIA3007 isolated in Mexico B-11A Mex as Rhizobium leguminosarum bv. viceae by bioinformatic tools

BACKGROUND

Evidence based on genomic sequences is extremely important to confirm the phylogenetic relationships within the Rhizobium group. SEMIA3007 was analyzed within the Mesorhizobium groups to define the underlying causes of taxonomic identification. We previously used biochemical tests and phenotypic taxonomic methods to identify bacteria, which can lead to erroneous classification. An improved understanding of bacterial strains such as the Mesorhizobium genus would increase our knowledge of classification and evolution of these species.

RESULTS

In this study, we sequenced the complete genome of SEMIA3007 and compared it with five other Mesorhizobium and two Rhizobium genomes. The genomes of isolated SEMIA3007 showed several orthologs with M. huakuii, M. erdmanii and M. loti. We identified SEMIA3007 as a Mesorhizobium by comparing the 16S rRNA gene and the complete genome.

CONCLUSION

Our ortholog, 16S rRNA gene and average nucleotide identity values (ANI) analysis all demonstrate SEMIA3007 is not Rhizobium leguminosarum bv. viceae. The results of the phylogenetic analysis clearly show SEMIA3007 is part of the Mesorhizobium group and suggest a reclassification is warranted.

Fourier transform infrared microspectroscopy for the analysis of the biochemical composition of C. elegans worms

Changes in intermediary metabolism have profound effects on many aspects of C. elegans biology including growth, development and behavior. However, many traditional biochemical techniques for analyzing chemical composition require relatively large amounts of starting material precluding the analysis of mutants that cannot be grown in large amounts as homozygotes. Here we describe a technique for detecting changes in the chemical compositions of C. elegans worms by Fourier transform infrared microspectroscopy. We demonstrate that the technique can be used to detect changes in the relative levels of carbohydrates, proteins and lipids in one and the same worm. We suggest that Fourier transform infrared microspectroscopy represents a useful addition to the arsenal of techniques for metabolic studies of C. elegans worms.

Insights into the origin and evolution of the plant hormone signaling machinery

Plant hormones modulate plant growth, development, and defense. However, many aspects of the origin and evolution of plant hormone signaling pathways remain obscure. Here, we use a comparative genomic and phylogenetic approach to investigate the origin and evolution of nine major plant hormone (abscisic acid, auxin, brassinosteroid, cytokinin, ethylene, gibberellin, jasmonate, salicylic acid, and strigolactone) signaling pathways. Our multispecies genome-wide analysis reveals that: (1) auxin, cytokinin, and strigolactone signaling pathways originated in charophyte lineages; (2) abscisic acid, jasmonate, and salicylic acid signaling pathways arose in the last common ancestor of land plants; (3) gibberellin signaling evolved after the divergence of bryophytes from land plants; (4) the canonical brassinosteroid signaling originated before the emergence of angiosperms but likely after the split of gymnosperms and angiosperms; and (5) the origin of the canonical ethylene signaling pathway postdates shortly the emergence of angiosperms. Our findings might have important implications in understanding the molecular mechanisms underlying the emergence of land plants.

Punctuated emergences of genetic and phenotypic innovations in eumetazoan, bilaterian, euteleostome, and hominidae ancestors

Phenotypic traits derive from the selective recruitment of genetic materials over macroevolutionary times, and protein-coding genes constitute an essential component of these materials. We took advantage of the recent production of genomic scale data from sponges and cnidarians, sister groups from eumetazoans and bilaterians, respectively, to date the emergence of human proteins and to infer the timing of acquisition of novel traits through metazoan evolution. Comparing the proteomes of 23 eukaryotes, we find that 33% human proteins have an ortholog in nonmetazoan species. This premetazoan proteome associates with 43% of all annotated human biological processes. Subsequently, four major waves of innovations can be inferred in the last common ancestors of eumetazoans, bilaterians, euteleostomi (bony vertebrates), and hominidae, largely specific to each epoch, whereas early branching deuterostome and chordate phyla show very few innovations. Interestingly, groups of proteins that act together in their modern human functions often originated concomitantly, although the corresponding human phenotypes frequently emerged later. For example, the three cnidarians Acropora, Nematostella, and Hydra express a highly similar protein inventory, and their protein innovations can be affiliated either to traits shared by all eumetazoans (gut differentiation, neurogenesis); or to bilaterian traits present in only some cnidarians (eyes, striated muscle); or to traits not identified yet in this phylum (mesodermal layer, endocrine glands). The variable correspondence between phenotypes predicted from protein enrichments and observed phenotypes suggests that a parallel mechanism repeatedly produce similar phenotypes, thanks to novel regulatory events that independently tie preexisting conserved genetic modules.

Evolution of the F-box gene family in Euarchontoglires: gene number variation and selection patterns

F-box proteins are substrate adaptors used by the SKP1–CUL1–F-box protein (SCF) complex, a type of E3 ubiquitin ligase complex in the ubiquitin proteasome system (UPS). SCF-mediated ubiquitylation regulates proteolysis of hundreds of cellular proteins involved in key signaling and disease systems. However, our knowledge of the evolution of the F-box gene family in Euarchontoglires is limited. In the present study, 559 F-box genes and nine related pseudogenes were identified in eight genomes. Lineage-specific gene gain and loss events occurred during the evolution of Euarchontoglires, resulting in varying F-box gene numbers ranging from 66 to 81 among the eight species. Both tandem duplication and retrotransposition were found to have contributed to the increase of F-box gene number, whereas mutation in the F-box domain was the main mechanism responsible for reduction in the number of F-box genes, resulting in a balance of expansion and contraction in the F-box gene family. Thus, the Euarchontoglire F-box gene family evolved under a birth-and-death model. Signatures of positive selection were detected in substrate-recognizing domains of multiple F-box proteins, and adaptive changes played a role in evolution of the Euarchontoglire F-box gene family. In addition, single nucleotide polymorphism (SNP) distributions were found to be highly non-random among different regions of F-box genes in 1092 human individuals, with domain regions having a significantly lower number of non-synonymous SNPs.

Systems biology approach to integrative comparative genomics

As more and more systems biology approaches are used to investigate the different types of biological macromolecules, increasing numbers of whole genomic studies are now available for a large array of organisms. Whether it is genomics, transcriptomics, proteomics, interactomics or metabolomics, the full complement of genomic information on all different levels can be juxtaposed between different organisms to reveal similarities or differences, and even to provide consensus models. At the intersection of comparative genomics and systems biology lies great possibility for discovery, analysis and prediction. This paper explores this nexus and the relationship from four general levels: DNA, RNA, protein and extragenomic. For each level, we provide an overview of the methods, discuss the potential challenges and survey the current research. Finally, we suggest some organizing principles and make proposals for new areas that will be important for future research.

Phylogenetic portrait of the Saccharomyces cerevisiae functional genome

The genome of budding yeast (Saccharomyces cerevisiae) contains approximately 5800 protein-encoding genes, the majority of which are associated with some known biological function. Yet the extent of amino acid sequence conservation of these genes over all phyla has only been partially examined. Here we provide a more comprehensive overview and visualization of the conservation of yeast genes and a means for browsing and exploring the data in detail, down to the individual yeast gene, at http://yeast-phylogroups.princeton.edu. We used data from the OrthoMCL database, which has defined orthologs from approximately 150 completely sequenced genomes, including diverse representatives of the archeal, bacterial, and eukaryotic domains. By clustering genes based on similar patterns of conservation, we organized and visualized all the protein-encoding genes in yeast as a single heat map. Most genes fall into one of eight major clusters, called "phylogroups." Gene ontology analysis of the phylogroups revealed that they were associated with specific, distinct trends in gene function, generalizations likely to be of interest to a wide range of biologists.

Common and pathogen-specific virulence factors are different in function and structure

In the process of host–pathogen interactions, bacterial pathogens always employ some special genes, e.g., virulence factors (VFs) to interact with host and cause damage or diseases to host. A number of VFs have been identified in bacterial pathogens that confer upon bacterial pathogens the ability to cause various types of damage or diseases. However, it has been clarified that some of the identified VFs are also encoded in the genomes of nonpathogenic bacteria, and this finding gives rise to considerable controversy about the definition of virulence factor.

Here 1988 virulence factors of 51 sequenced pathogenic bacterial genomes from the virulence factor database (VFDB) were collected, and an orthologous comparison to a non-pathogenic bacteria protein database was conducted using the reciprocal-best-BLAST-hits approach. Six hundred and twenty pathogen-specific VFs and 1368 common VFs (present in both pathogens and nonpathogens) were identified, which account for 31.19% and 68.81% of the total VFs, respectively. The distribution of pathogen-specific VFs and common VFs in pathogenicity islands (PAIs) was systematically investigated, and pathogen-specific VFs were more likely to be located in PAIs than common VFs. The function of the two classes of VFs were also analyzed and compared in depth. Our results indicated that most but not all T3SS proteins are pathogen-specific. T3SS effector proteins tended to be distributed in pathogen-specific VFs, whereas T3SS translocation proteins, apparatus proteins, and chaperones were inclined to be distributed in common VFs. We also observed that exotoxins were located in both pathogen-specific and common VFs. In addition, the architecture of the two classes of VFs was compared, and the results indicated that common VFs had a higher domain number and lower domain coverage value, revealed that common VFs tend to be more complex and less compact proteins.

A study on the distribution of 37 well conserved families of C2H2 zinc finger genes in eukaryotes

Background

The C2H2 zinc-finger (ZNF) containing gene family is one of the largest and most complex gene families in metazoan genomes. These genes are known to exist in almost all eukaryotes, and they constitute a major subset of eukaryotic transcription factors. The genes of this family usually occur as clusters in genomes and are thought to have undergone a massive expansion in vertebrates by multiple tandem duplication events (BMC Evol Biol 8:176, 2008).

Results

In this study, we combined two popular approaches for homolog detection, Reciprocal Best Hit (RBH) (Proc Natl Acad Sci USA 95:6239–6244, 1998) and Hidden–Markov model (HMM) profiles search (Bioinformatics 14:755-763, 1998), on a diverse set of complete genomes of 124 eukaryotic species ranging from excavates to humans to identify all detectable members of 37 C2H2 ZNF gene families. We succeeded in identifying 3,890 genes as distinct members of 37 C2H2 gene families. These 37 families are distributed among the eukaryotes as progressive additions of gene blocks with increasing complexity of the organisms. The first block featuring the protists had 7 families, the second block featuring plants had 2 families, the third block featuring the fungi had 2 families (one of which was also present in plants) and the final block consisted of metazoans with 25 families. Among the metazoans, the simpler unicellular metazoans had just 15 of the 25 families while most of the bilaterians had all 25 families making up a total of 37 families. Multiple potential examples of lineage-specific gene duplications and gene losses were also observed.

Conclusions

Our hybrid approach combines features of the both RBH and HMM methods for homolog detection. This largely automated technique is much faster than manual methods and is able to detect homologs accurately and efficiently among a diverse set of organisms. Our analysis of the 37 evolutionarily conserved C2H2 ZNF gene families revealed a stepwise appearance of ZNF families, agreeing well with the phylogenetic relationship of the organisms compared and their presumed stepwise increase in complexity (Science 300:1694, 2003).

Congratulations, you have been carefully chosen to represent an important developmental regulator!

BACKGROUND

Studying a process in a new species often relies on focusing our attention to a candidate gene, encoding a protein similar to one with a known function. Not all the choices seem to be prudent.

SCOPE

This Viewpoint includes an overview of issues that are encountered during research of candidate genes. Defining a match for a gene of interest, deciding whether variation in ESTs or RNAseq data for a certain transcript, represent more than one gene. The problem of incorrect annotation of genes due to incorrect in-silico splicing, is also mentioned. The author's humble opinion on how to deal with these issues is provided.

CONCLUSIONS

The vast amount of new sequence data provides us with great possibilities for giant leaps in our understanding. Still, we cannot afford to skip over the tedious steps required to confirm that we are indeed studying the correct gene, and try to be sure that the complex expression pattern we observe is not a composite of several genes.

Localization and substrate selectivity of sea urchin multidrug (MDR) efflux transporters

In this study, we cloned, expressed and functionally characterized Stronglycentrotus purpuratus (Sp) ATP-binding cassette (ABC) transporters. This screen identified three multidrug resistance (MDR) transporters with functional homology to the major types of MDR transporters found in humans. When overexpressed in embryos, the apical transporters Sp-ABCB1a, ABCB4a, and ABCG2a can account for as much as 87% of the observed efflux activity, providing a robust assay for their substrate selectivity. Using this assay, we found that sea urchin MDR transporters export canonical MDR susbtrates such as calcein-AM, bodipy-verapamil, bodipy-vinblastine, and mitoxantrone. In addition, we characterized the impact of nonconservative substitutions in the primary sequences of drug binding domains of sea urchin versus murine ABCB1 by mutation of Sp-ABCB1a and treatment of embryos with stereoisomeric cyclic peptide inhibitors (QZ59 compounds). The results indicated that two substitutions in transmembrane helix 6 reverse stereoselectivity of Sp-ABCB1a for QZ59 enantiomers compared with mouse ABCB1a. This suggests that subtle changes in the primary sequence of transporter drug binding domains could fine-tune substrate specificity through evolution.

large-scale screening for targeted knockouts in the Caenorhabditis elegans genome

The nematode Caenorhabditis elegans is a powerful model system to study contemporary biological problems. This system would be even more useful if we had mutations in all the genes of this multicellular metazoan. The combined efforts of the C. elegans Deletion Mutant Consortium and individuals within the worm community are moving us ever closer to this goal. At present, of the 20,377 protein-coding genes in this organism, 6764 genes with associated molecular lesions are either deletions or null mutations (WormBase WS220). Our three laboratories have contributed the majority of mutated genes, 6841 mutations in 6013 genes. The principal method we used to detect deletion mutations in the nematode utilizes polymerase chain reaction (PCR). More recently, we have used array comparative genome hybridization (aCGH) to detect deletions across the entire coding part of the genome and massively parallel short-read sequencing to identify nonsense, splicing, and missense defects in open reading frames. As deletion strains can be frozen and then thawed when needed, these strains will be an enduring community resource. Our combined molecular screening strategies have improved the overall throughput of our gene-knockout facilities and have broadened the types of mutations that we and others can identify. These multiple strategies should enable us to eventually identify a mutation in every gene in this multicellular organism. This knowledge will usher in a new age of metazoan genetics in which the contribution to any biological process can be assessed for all genes.

Beyond 'furballs' and 'dumpling soups' - towards a molecular architecture of signaling complexes and networks

The molecular architectures of intracellular signaling networks are largely unknown. Understanding their design principles and mechanisms of processing information is essential to grasp the molecular basis of virtually all biological processes. This is particularly challenging for human pathologies like cancers, as essentially each tumor is a unique disease with vastly deranged signaling networks. However, even in normal cells we know almost nothing. A few 'signalosomes', like the COP9 and the TCR signaling complexes have been described, but detailed structural information on their architectures is largely lacking. Similarly, many growth factor receptors, for example EGF receptor, insulin receptor and c-Met, signal via huge protein complexes built on large platform proteins (Gab, Irs/Dok, p130Cas[BCAR1], Frs families etc.), which are structurally not well understood. Subsequent higher order processing events remain even more enigmatic. We discuss here methods that can be employed to study signaling architectures, and the importance of too often neglected features like macromolecular crowding, intrinsic disorder in proteins and the sophisticated cellular infrastructures, which need to be carefully considered in order to develop a more mature understanding of cellular signal processing.

OrthoDisease: tracking disease gene orthologs across 100 species

Orthology is one of the most important tools available to modern biology, as it allows making inferences from easily studied model systems to much less tractable systems of interest, such as ourselves. This becomes important not least in the study of genetic diseases. We here review work on the orthology of disease-associated genes and also present an updated version of the InParanoid-based disease orthology database and web site OrthoDisease, with 14-fold increased species coverage since the previous version. Using this resource, we survey the taxonomic distribution of orthologs of human genes involved in different disease categories. The hypothesis that paralogs can mask the effect of deleterious mutations predicts that known heritable disease genes should have fewer close paralogs. We found large-scale support for this hypothesis as significantly fewer duplications were observed for disease genes in the OrthoDisease ortholog groups.

A survey of yeast genomic assays for drug and target discovery

Over the past decade, the development and application of chemical genomic assays using the model organism Saccharomyces cerevisiae has provided powerful methods to identify the mechanism of action of known drugs and novel small molecules in vivo. These assays identify drug target candidates, genes involved in buffering drug target pathways and also help to define the general cellular response to small molecules. In this review, we examine current yeast chemical genomic assays and summarize the potential applications of each approach.

The protein meta-structure: a novel concept for chemical and molecular biology

The ultimate goal of bioinformatics or computational chemical biology is the sequence-based prediction of protein functionality. However, due to the degeneracy of the primary sequence code there is no unambiguous relationship. The degeneracy can be partly lifted by going to higher levels of abstraction and, for example, incorporating 3D structural information. However, sometimes even at this conceptual level functional ambiguities often remain. Here a novel conceptual framework is described (the protein meta-structure). At this level of abstraction, the protein structure is viewed as an intricate network of interacting residues. This novel conception offers unique possibilities for chemical (molecular) biology, structural genomics and drug discovery. In this review some prototypical applications will be presented that serve to illustrate the potential of the methodology.

ANNIE: integrated de novo protein sequence annotation

Function prediction of proteins with computational sequence analysis requires the use of dozens of prediction tools with a bewildering range of input and output formats. Each of these tools focuses on a narrow aspect and researchers are having difficulty obtaining an integrated picture. ANNIE is the result of years of close interaction between computational biologists and computer scientists and automates an essential part of this sequence analytic process. It brings together over 20 function prediction algorithms that have proven sufficiently reliable and indispensable in daily sequence analytic work and are meant to give scientists a quick overview of possible functional assignments of sequence segments in the query proteins. The results are displayed in an integrated manner using an innovative AJAX-based sequence viewer. ANNIE is available online at: http://annie.bii.a-star.edu.sg. This website is free and open to all users and there is no login requirement.

A novel hybrid yeast-human network analysis reveals an essential role for FNBP1L in antibacterial autophagy

Autophagy is a conserved cellular process required for the removal of defective organelles, protein aggregates, and intracellular pathogens. We used a network analysis strategy to identify novel human autophagy components based upon the yeast interactome centered on the core yeast autophagy proteins. This revealed the potential involvement of 14 novel mammalian genes in autophagy, several of which have known or predicted roles in membrane organization or dynamics. We selected one of these membrane interactors, FNBP1L (formin binding protein 1-like), an F-BAR-containing protein (also termed Toca-1), for further study based upon a predicted interaction with ATG3. We confirmed the FNBP1L/ATG3 interaction biochemically and mapped the FNBP1L domains responsible. Using a functional RNA interference approach, we determined that FNBP1L is essential for autophagy of the intracellular pathogen Salmonella enterica serovar Typhimurium and show that the autophagy process serves to restrict the growth of intracellular bacteria. However, FNBP1L appears dispensable for other forms of autophagy induced by serum starvation or rapamycin. We present a model where FNBP1L is essential for autophagy of intracellular pathogens and identify FNBP1L as a differentially used molecule in specific autophagic contexts. By using network biology to derive functional biological information, we demonstrate the utility of integrated genomics to novel molecule discovery in autophagy.

Expressional profiling study revealed unique expressional patterns and dramatic expressional divergence of maize alpha-zein super gene family

The alpha-zein super gene family encodes the most predominant storage protein in maize (Zea mays) endosperm. In maize inbred line B73, it consists of four gene families with 41 member genes. In this study, we combined quantitative real-time PCR and random clone sequencing to successfully profile the expression of alpha-zein super gene family during endosperm development. We found that only 18 of the 41 member genes were expressed, and their expression levels diverge greatly. At the gene family level, all families had characteristic "up-and-down" oscillating expressional patterns that diverged into two major groups. At the individual gene level, member genes showed dramatic divergence of expression patterns, indicating fast differentiation of their expression regulation. A comparison study among different inbred lines revealed significantly different expressed gene sets, indicating the existence of highly diverged haplotypes. Large gene families containing long gene clusters, e.g. z1A or z1C, mainly contributed the highly divergent haplotypes. In addition, allelic genes also showed significant divergence in their expressional levels. These results indicated a highly dynamic and fast evolving nature to the maize alpha-zein super gene family, which might be a common feature for other large gene families.

The origins of multicellularity and the early history of the genetic toolkit for animal development

Multicellularity appeared early and repeatedly in life's history; its instantiations presumably required the confluence of environmental, ecological, and genetic factors. Comparisons of several independently evolved pairs of multicellular and unicellular relatives indicate that transitions to multicellularity are typically associated with increases in the numbers of genes involved in cell differentiation, cell-cell communication, and adhesion. Further examination of the DNA record suggests that these increases in gene complexity are the product of evolutionary innovation, tinkering, and expansion of genetic material. Arguably, the most decisive multicellular transition was the emergence of animals. Decades of developmental work have demarcated the genetic toolkit for animal multicellularity, a select set of a few hundred genes from a few dozen gene families involved in adhesion, communication, and differentiation. Examination of the DNA records of the earliest-branching animal phyla and their closest protist relatives has begun to shed light on the origins and assembly of this toolkit. Emerging data favor a model of gradual assembly, with components originating and diversifying at different time points prior to or shortly after the origin of animals.

Description of International Caenorhabditis elegans Experiment first flight (ICE-FIRST)

Traveling, living and working in space is now a reality. The number of people and length of time in space is increasing. With new horizons for exploration it becomes more important to fully understand and provide countermeasures to the effects of the space environment on the human body. In addition, space provides a unique laboratory to study how life and physiologic functions adapt from the cellular level to that of the entire organism. Caenorhabditis elegans is a genetic model organism used to study physiology on Earth. Here we provide a description of the rationale, design, methods, and space culture validation of the ICE-FIRST payload, which engaged C. elegans researchers from four nations. Here we also show C. elegans growth and development proceeds essentially normally in a chemically defined liquid medium on board the International Space Station (10.9 day round trip). By setting flight constraints first and bringing together established C. elegans researchers second, we were able to use minimal stowage space to successfully return a total of 53 independent samples, each containing more than a hundred individual animals, to investigators within one year of experiment concept. We believe that in the future, bringing together individuals with knowledge of flight experiment operations, flight hardware, space biology, and genetic model organisms should yield similarly successful payloads.

Clustering orthologous proteins across phylogenetically distant species

The quality of orthologous protein clusters (OPCs) is largely dependent on the results of the reciprocal BLAST (basic local alignment search tool) hits among genomes. The BLAST algorithm is very efficient and fast, but it is very difficult to get optimal solution among phylogenetically distant species because the genomes with large evolutionary distance typically have low similarity in their protein sequences. To reduce the false positives in the OPCs, thresholding is often employed on the BLAST scores. However, the thresholding also eliminates large numbers of true positives as the orthologs from distant species likely have low BLAST scores. To rectify this problem, we introduce a new hybrid method combining the Recursive and the Markov CLuster (MCL) algorithms without using the BLAST thresholding. In the first step, we use InParanoid to produce n(n-1)/2 ortholog tables from n genomes. After combining all the tables into one, our clustering algorithm clusters ortholog pairs recursively in the table. Then, our method employs MCL algorithm to compute the clusters and refines the clusters by adjusting the inflation factor. We tested our method using six different genomes and evaluated the results by comparing against Kegg Orthology (KO) OPCs, which are generated from manually curated pathways. To quantify the accuracy of the results, we introduced a new intuitive similarity measure based on our Least-move algorithm that computes the consistency between two OPCs. We compared the resulting OPCs with the KO OPCs using this measure. We also evaluated the performance of our method using InParanoid as the baseline approach. The experimental results show that, at the inflation factor 1.3, we produced 54% more orthologs than InParanoid sacrificing a little less accuracy (1.7% less) than InParanoid, and at the factor 1.4, produced not only 15% more orthologs than InParanoid but also a higher accuracy (1.4% more) than InParanoid.

On the extent and origins of genic novelty in the phylum Nematoda

BACKGROUND

The phylum Nematoda is biologically diverse, including parasites of plants and animals as well as free-living taxa. Underpinning this diversity will be commensurate diversity in expressed genes, including gene sets associated specifically with evolution of parasitism.

METHODS AND FINDINGS

Here we have analyzed the extensive expressed sequence tag data (available for 37 nematode species, most of which are parasites) and define over 120,000 distinct putative genes from which we have derived robust protein translations. Combined with the complete proteomes of Caenorhabditis elegans and Caenorhabditis briggsae, these proteins have been grouped into 65,000 protein families that in turn contain 40,000 distinct protein domains. We have mapped the occurrence of domains and families across the Nematoda and compared the nematode data to that available for other phyla. Gene loss is common, and in particular we identify nearly 5,000 genes that may have been lost from the lineage leading to the model nematode C. elegans. We find a preponderance of novelty, including 56,000 nematode-restricted protein families and 26,000 nematode-restricted domains. Mapping of the latest time-of-origin of these new families and domains across the nematode phylogeny revealed ongoing evolution of novelty. A number of genes from parasitic species had signatures of horizontal transfer from their host organisms, and parasitic species had a greater proportion of novel, secreted proteins than did free-living ones.

CONCLUSIONS

These classes of genes may underpin parasitic phenotypes, and thus may be targets for development of effective control measures.

Orthology and functional conservation in eukaryotes

In recent years, it has become clear that all of the organisms on the Earth are related to each other in ways that can be documented by molecular sequence comparison. In this review, we focus on the evolutionary relationships among the proteins of the eukaryotes, especially those that allow inference of function from one species to another. Data and illustrations are derived from specific comparison of eight species: Homo sapiens, Mus musculus, Arabidopsis thaliana, Caenorhabditis elegans, Danio rerio, Saccharomyces cerevisiae, and Plasmodium falciparum.

Linking the genomes of nonmodel teleosts through comparative genomics

Recently the genomes of two more teleost species have been released: the medaka (Oryzias latipes), and the three-spined stickleback (Gasterosteus aculateus). The rapid developments in genomics of fish species paved the way to new and valuable research in comparative genetics and genomics. With the accumulation of information in model species, the genetic and genomic characterization of nonmodel, but economically important species, is now feasible. Furthermore, comparison of low coverage gene maps of aquacultured fish species against fully sequenced fish species will enhance the efficiency of candidate genes identification projected for quantitative trait loci (QTL) scans for traits of commercial interest. This study shows the syntenic relationship between the genomes of six different teleost species, including three fully sequenced model species: Tetraodon nigroviridis, Oryzias latipes, Gasterosteus aculateus, and three marine species of commercial and evolutionary interest: Sparus aurata, Dicentrarchus labrax, Oreochromis spp. All three commercial fish species belong to the order Perciformes, which is the richest in number of species (approximately 10,000) but poor in terms of available genomic information and tools. Syntenic relationships were established by using 800 EST and microsatellites sequences successfully mapped on the RH map of seabream. Comparison to the stickleback genome produced most positive BLAT hits (58%) followed by medaka (32%) and Tetraodon (30%). Thus, stickleback was used as the major stepping stone to compare seabass and tilapia to seabream. In addition to the significance for the aquaculture industry, this approach can encompass important ecological and evolutionary implications.

Functional protein divergence in the evolution of Homo sapiens

BACKGROUND

Protein-coding regions in a genome evolve by sequence divergence and gene gain and loss, altering the gene content of the organism. However, it is not well understood how this has given rise to the enormous diversity of metazoa present today.

RESULTS

To obtain a global view of human genomic evolution, we quantify the divergence of proteins by functional category at different evolutionary distances from human.

CONCLUSION

This analysis highlights some general systems-level characteristics of human evolution: regulatory processes, such as signal transducers, transcription factors and receptors, have a high degree of plasticity, while core processes, such as metabolism, transport and protein synthesis, are largely conserved. Additionally, this study reveals a dynamic picture of selective forces at short, medium and long evolutionary timescales. Certain functional categories, such as 'development' and 'organogenesis', exhibit temporal patterns of sequence divergence in eukaryotes relative to human. This framework for a grammar of human evolution supports previously postulated theories of robustness and evolvability.

Evolution of protein domain promiscuity in eukaryotes

Numerous eukaryotic proteins contain multiple domains. Certain domains show a tendency to occur in diverse domain architectures and can be considered "promiscuous." These promiscuous domains are, typically, involved in protein-protein interactions and play crucial roles in interaction networks, particularly those that contribute to signal transduction. A systematic comparative-genomic analysis of promiscuous domains in eukaryotes is described. Two quantitative measures of domain promiscuity are introduced and applied to the analysis of 28 genomes of diverse eukaryotes. Altogether, 215 domains are identified as strongly promiscuous. The fraction of promiscuous domains in animals is shown to be significantly greater than that in fungi or plants. Evolutionary reconstructions indicate that domain promiscuity is a volatile, relatively fast-changing feature of eukaryotic proteins, with few domains remaining promiscuous throughout the evolution of eukaryotes. Some domains appear to have attained promiscuity independently in different lineages, for example, animals and plants. It is proposed that promiscuous domains persist within a relatively small pool of evolutionarily stable domain combinations from which numerous rare architectures emerge during evolution. Domain promiscuity positively correlates with the number of experimentally detected domain interactions and with the strength of purifying selection affecting a domain. Thus, evolution of promiscuous domains seems to be constrained by the diversity of their interaction partners. The set of promiscuous domains is enriched for domains mediating protein-protein interactions that are involved in various forms of signal transduction, especially in the ubiquitin system and in chromatin. Thus, a limited repertoire of promiscuous domains makes a major contribution to the diversity and evolvability of eukaryotic proteomes and signaling networks.

Family expansion and gene rearrangements contributed to the functional specialization of PRDM genes in vertebrates

Background

Progressive diversification of paralogs after gene expansion is essential to increase their functional specialization. However, mode and tempo of this divergence remain mostly unclear. Here we report the comparative analysis of PRDM genes, a family of putative transcriptional regulators involved in human tumorigenesis.

Results

Our analysis assessed that the PRDM genes originated in metazoans, expanded in vertebrates and further duplicated in primates. We experimentally showed that fast-evolving paralogs are poorly expressed, and that the most recent duplicates, such as primate-specific PRDM7, acquire tissue-specificity. PRDM7 underwent major structural rearrangements that decreased the number of encoded Zn-Fingers and modified gene splicing. Through internal duplication and activation of a non-canonical splice site (GC-AG), PRDM7 can acquire a novel intron. We also detected an alternative isoform that can retain the intron in the mature transcript and that is predominantly expressed in human melanocytes.

Conclusion

Our findings show that (a) molecular evolution of paralogs correlates with their expression pattern; (b) gene diversification is obtained through massive genomic rearrangements; and (c) splicing modification contributes to the functional specialization of novel genes.

Systematic identification of SH3 domain-mediated human protein–protein interactions by peptide array target screening

Systematic identification of direct protein-protein interactions is often hampered by difficulties in expressing and purifying the corresponding full-length proteins. By taking advantage of the modular nature of many regulatory proteins, we attempted to simplify protein-protein interactions to the corresponding domain-ligand recognition and employed peptide arrays to identify such binding events. A group of 12 Src homology (SH) 3 domains from eight human proteins (Swiss-Prot ID: SRC, PLCG1, P85A, NCK1, GRB2, FYN, CRK) were used to screen a peptide target array composed of 1536 potential ligands, which led to the identification of 921 binary interactions between these proteins and 284 targets. To assess the efficiency of the peptide array target screening (PATS) method in identifying authentic protein-protein interactions, we examined a set of interactions mediated by the PLCgamma1 SH3 domain by coimmunoprecipitation and/or affinity pull-downs using full-length proteins and achieved a 75% success rate. Furthermore, we characterized a novel interaction between PLCgamma1 and hematopoietic progenitor kinase 1 (HPK1) identified by PATS and demonstrated that the PLCgamma1 SH3 domain negatively regulated HPK1 kinase activity. Compared to protein interactions listed in the online predicted human interaction protein database (OPHID), the majority of interactions identified by PATS are novel, suggesting that, when extended to the large number of peptide interaction domains encoded by the human genome, PATS should aid in the mapping of the human interactome.

Modeling tumor predisposing FH mutations in yeast: effects on fumarase activity, growth phenotype and gene expression profile

Heterozygous mutations in the fumarase (FH) gene cause the tumor predisposition syndrome hereditary leiomyomatosis and renal cell cancer (MIM 605839). While most families segregate a benign phenotype of multiple leiomyomas, others display a phenotype with early-onset renal cancer and leiomyosarcoma. Modifier genes may play a role in this, but an alternative explanation is simple genotype-phenotype association. FH mutations predisposing to cancer appear to be truncating or in fully conserved amino acids, suggesting that mutations severely affecting FH activity might predispose to malignancy. In the present study, we analyzed 2 conserved fumarase mutations in yeast. H153R has been described in 3 cancer predisposition families; whereas all 3 reported K187R families have displayed the benign phenotype. Examining H153R and K187R should clarify whether cancer-related FH mutations differ from their benign phenotype-associated counterparts. Yeast strains containing the 2 mutations, and knockout and wild type (WT) references, were created and the growth phenotypes studied on selected carbon sources to assess mitochondrial function. Additionally, Fum1 protein production and activity were measured, and the strains were subjected to transcriptional profiling. On nonfermentable lactate medium, the fumarase knockout strains did not grow, whereas the mutants showed no differences, as compared to WT yeast. Although both mutant strains produced fumarase, a considerable decrease in enzyme activity was seen in mutants with respect to WT. Transcription of the majority of Krebs cycle enzymes was downregulated in response to mutations in fumarase. In conclusion, both mutants displayed some, albeit greatly reduced, fumarase activity. This activity was sufficient to support normal growth on nonfermentable carbon source, unlike the deletion phenotype, demonstrating the significance of the residual activity. The findings support the hypothesis that modifier gene(s), rather than phenotype-genotype effects, display a major role in determining tumor phenotypes in families segregating FH mutations.

Regulation of glutamate receptor B pre-mRNA splicing by RNA editing

RNA-editing enzymes of the ADAR family convert adenosines to inosines in double-stranded RNA substrates. Frequently, editing sites are defined by base-pairing of the editing site with a complementary intronic region. The glutamate receptor subunit B (GluR-B) pre-mRNA harbors two such exonic editing sites termed Q/R and R/G. Data from ADAR knockout mice and in vitro editing assays suggest an intimate connection between editing and splicing of GluR-B pre-mRNA.

By comparing the events at the Q/R and R/G sites, we can show that editing can both stimulate and repress splicing efficiency. The edited nucleotide, but not ADAR binding itself, is sufficient to exert this effect. The presence of an edited nucleotide at the R/G site reduces splicing efficiency of the adjacent intron facilitating alternative splicing events occurring downstream of the R/G site.

Lack of editing inhibits splicing at the Q/R site. Editing of both the Q/R nucleotide and an intronic editing hotspot are required to allow efficient splicing. Inefficient intron removal may ensure that only properly edited mRNAs become spliced and exported to the cytoplasm.

Protein subcomplexes--molecular machines with highly specialized functions

Complex cellular processes are accomplished by the concerted action of hierarchically organized functional modules. Protein complexes are major components which act as highly specialized molecular machines. Here we present a statistical procedure to find insightful substructures in protein complexes based on large-scale protein complex data: we identify statistically significant common protein subcomplexes (SCs) contained in different protein complexes. We analyze recently published data of the two model organisms Saccharomyces cerevisiae (four different data sets) and Escherichia coli, as well as human protein complex data. Our method identifies well-characterized protein assemblies with known functions which act as own functional entities in the cell. In addition, we also identified hitherto unknown functional entities that should be studied experimentally in future. We discuss two typical properties of protein subcomplexes: 1) subcomplexes are enriched with essential proteins (which implies that the whole SCs may be strongly conserved) and 2) SCs are functionally and spatially more homogeneous than the experimentally found protein assemblies. The latter property is exploited to propose functions for so far unknown proteins of S. cerevisiae.

Growth-rate regulated genes have profound impact on interpretation of transcriptome profiling in Saccharomyces cerevisiae

BACKGROUND

Growth rate is central to the development of cells in all organisms. However, little is known about the impact of changing growth rates. We used continuous cultures to control growth rate and studied the transcriptional program of the model eukaryote Saccharomyces cerevisiae, with generation times varying between 2 and 35 hours.

RESULTS

A total of 5930 transcripts were identified at the different growth rates studied. Consensus clustering of these revealed that half of all yeast genes are affected by the specific growth rate, and that the changes are similar to those found when cells are exposed to different types of stress (>80% overlap). Genes with decreased transcript levels in response to faster growth are largely of unknown function (>50%) whereas genes with increased transcript levels are involved in macromolecular biosynthesis such as those that encode ribosomal proteins. This group also covers most targets of the transcriptional activator RAP1, which is also known to be involved in replication. A positive correlation between the location of replication origins and the location of growth-regulated genes suggests a role for replication in growth rate regulation.

CONCLUSION

Our data show that the cellular growth rate has great influence on transcriptional regulation. This, in turn, implies that one should be cautious when comparing mutants with different growth rates. Our findings also indicate that much of the regulation is coordinated via the chromosomal location of the affected genes, which may be valuable information for the control of heterologous gene expression in metabolic engineering.

Long-term trends in evolution of indels in protein sequences

Background

In this paper we describe an analysis of the size evolution of both protein domains and their indels, as inferred by changing sizes of whole domains or individual unaligned regions or "spacers". We studied relatively early evolutionary events and focused on protein domains which are conserved among various taxonomy groups.

Results

We found that more than one third of all domains have a statistically significant tendency to increase/decrease in size in evolution as judged from the overall domain size distribution as well as from the size distribution of individual spacers. Moreover, the fraction of domains and individual spacers increasing in size is almost twofold larger than the fraction decreasing in size.

Conclusion

We showed that the tolerance to insertion and deletion events depends on the domain's taxonomy span. Eukaryotic domains are depleted in insertions compared to the overall test set, namely, the number of spacers increasing in size is about the same as the number of spacers decreasing in size. On the other hand, ancient domain families show some bias towards insertions or spacers which grow in size in evolution. Domains from several Gene Ontology categories also demonstrate certain tendencies for insertion or deletion events as inferred from the analysis of spacer sizes.

Similarity of the domain structure of proteins as a basis for the conservation of meiosis

Meiosis is conserved in all eucaryotic kingdoms, and homologous rows of variability are revealed for the cytological traits of meiosis. To find the nature of these phenomenons, we reviewed the most-studied meiosis-specific proteins and studied them with the methods of bioinformatics. We found that synaptonemal complex proteins have no homology of amino-acid sequence, but are similar in the domain organization and three-dimensional (3D) structure of functionally important domains in budding yeast, nematode, Drosophila, Arabidopsis, and human. Recombination proteins of Rad51/Dmc1 family are conserved to the extent which permits them to make filamentous single-strand deoxyribonucleic acid (ssDNA)-protein intermediates of meiotic recombination. The same structural principles are valid for conservation of the ultrastructure of kinetochores, cell gap contacts, and nuclear pore complexes, such as in the cases when ultrastructure 3D parameters are important for the function. We suggest that self-assembly of protein molecules plays a significant role in building-up of all biological structures mentioned.