PubMed Health⌕ Search

Biomedical subjects

Eugene V Koonin

Publications and source records attributed to Eugene V Koonin.

At least 19 recordsLinked to original sources
PubMed2026

Type IV-C CRISPR-Cas effector complexes recognize double-stranded DNA and switch on collateral cleavage of ssDNA and RNA.

Type IV-C CRISPR-Cas systems remain enigmatic compared to other class 1 systems. Here, we expand the type IV-C catalog, identifying two phylogenetically distinct clades primarily found in archaea (IV-C1) or bacteria (IV-C2), distinguishable by the Cas10IVc subunit architecture. We functionally and structurally characterize type IV-C1 systems from Thermococcus onnurineus (Ton) and Pyrococcus abyssi (Pab). Type IV-C complexes assemble with crRNAs derived from distinct CRISPR arrays and recognize a 5'-GGG-3' protospacer adjacent motif (PAM) to bind double-stranded DNA targets. Target recognition activates the HD domain of Cas10IVc, triggering metal-dependent collateral cleavage of single-stranded DNA and RNA. This behavior is explained by allosteric alignment of the HD active site, triggered by PAM-dependent R-loop formation, as revealed by cryo-EM. Together, our findings suggest that type IV-C systems provide immunity via non-specific cleavage of nucleic acids generated during mobile genetic element replication or transcription.

CP: molecular biology↗
PubMed2024

Natural history of eukaryotic DNA viruses with double jelly-roll major capsid proteins.

The phylum Preplasmiviricota (kingdom Bamfordvirae, realm Varidnaviria) is a broad assemblage of diverse viruses with comparatively short double-stranded DNA genomes (<50 kbp) that produce icosahedral capsids built from double jelly-roll major capsid proteins. Preplasmiviricots infect hosts from all cellular domains, testifying to their ancient origin, and, in particular, are associated with six of the seven supergroups of eukaryotes. Preplasmiviricots comprise four major groups of viruses, namely, polintons, polinton-like viruses (PLVs), virophages, and adenovirids. We used protein structure modeling and analysis to show that protein-primed DNA polymerases (pPolBs) of polintons, virophages, and cytoplasmic linear plasmids encompass an N-terminal domain homologous to the terminal proteins (TPs) of prokaryotic PRD1-like tectivirids and eukaryotic adenovirids that are involved in protein-primed replication initiation, followed by a viral ovarian tumor-like cysteine deubiquitinylase (vOTU) domain. The vOTU domain is likely responsible for the cleavage of the TP from the large pPolB polypeptide and is inactivated in adenovirids, in which TP is a separate protein. Many PLVs and transpovirons encode a distinct derivative of polinton-like pPolB that retains the TP, vOTU, and pPolB polymerization palm domains but lacks the exonuclease domain and instead contains a superfamily 1 helicase domain. Analysis of the presence/absence and inactivation of the vOTU domains and replacement of pPolB with other DNA polymerases in eukaryotic preplasmiviricots enabled us to outline a complete scenario for their origin and evolution.

Capsid Proteins↗
PubMed2024

Natural history of eukaryotic DNA viruses with double jelly-roll major capsid proteins.

The phylum Preplasmiviricota (kingdom Bamfordvirae, realm Varidnaviria) is a broad assemblage of diverse viruses with comparatively short double-stranded DNA genomes (<50 kbp) that produce icosahedral capsids built from double jelly-roll major capsid proteins. Preplasmiviricots infect hosts from all cellular domains, testifying to their ancient origin and, in particular, are associated with six of the seven supergroups of eukaryotes. Preplasmiviricots comprise four major groups of viruses, namely, polintons, polinton-like viruses (PLVs), virophages, and adenovirids. We employed protein structure modeling and analysis to show that protein-primed DNA polymerases (pPolBs) of polintons, virophages, and cytoplasmic linear plasmids encompass an N-terminal domain homologous to the terminal proteins (TPs) of prokaryotic PRD1-like tectivirids and eukaryotic adenovirids that are involved in protein-primed replication initiation, followed by a viral ovarian tumor-like cysteine deubiquitinylase (vOTU) domain. The vOTU domain is likely responsible for the cleavage of the TP from the large pPolB polypeptide and is inactivated in adenovirids, in which TP is a separate protein. Many PLVs and transpovirons encode a distinct derivative of polinton-like pPolB that retains the TP, vOTU and pPolB polymerization palm domains but lacks the exonuclease domain and instead contains a supefamily 1 helicase domain. Analysis of the presence/absence and inactivation of the vOTU domains, and replacement of pPolB with other DNA polymerases in eukaryotic preplasmiviricots enabled us to outline a complete scenario for their origin and evolution.

adenovirus↗
PubMed2024

Double-stranded RNA sequencing reveals distinct riboviruses associated with thermoacidophilic bacteria from hot springs in Japan.

Metatranscriptome sequencing expanded the known diversity of the bacterial RNA virome, suggesting that additional riboviruses infecting bacterial hosts remain to be discovered. Here we employed double-stranded RNA sequencing to recover complete genome sequences of two ribovirus groups from acidic hot springs in Japan. One group, denoted hot spring riboviruses (HsRV), consists of viruses with distinct RNA-directed RNA polymerases (RdRPs) that seem to be intermediates between typical ribovirus RdRPs and viral reverse transcriptases. This group forms a distinct phylum, Artimaviricota, or even kingdom within the realm Riboviria. We identified viruses encoding HsRV-like RdRPs in marine water, river sediments and salt marshes, indicating that this group is widespread beyond extreme ecosystems. The second group, denoted hot spring partiti-like viruses (HsPV), forms a distinct branch within the family Partitiviridae. The genome architectures of HsRV and HsPV and their identification in bacteria-dominated habitats suggest that these viruses infect thermoacidophilic bacteria.

Hot Springs↗
PubMed2006

Temporal order of evolution of DNA replication systems inferred by comparison of cellular and viral DNA polymerases.

BACKGROUND: The core enzymes of the DNA replication systems show striking diversity among cellular life forms and more so among viruses. In particular, and counter-intuitively, given the central role of DNA in all cells and the mechanistic uniformity of replication, the core enzymes of the replication systems of bacteria and archaea (as well as eukaryotes) are unrelated or extremely distantly related. Viruses and plasmids, in addition, possess at least two unique DNA replication systems, namely, the protein-primed and rolling circle modalities of replication. This unexpected diversity makes the origin and evolution of DNA replication systems a particularly challenging and intriguing problem in evolutionary biology. RESULTS: I propose a specific succession for the emergence of different DNA replication systems, drawing argument from the differences in their representation among viruses and other selfish replicating elements. In a striking pattern, the DNA replication systems of viruses infecting bacteria and eukaryotes are dominated by the archaeal-type B-family DNA polymerase (PolB) whereas the bacterial replicative DNA polymerase (PolC) is present only in a handful of bacteriophage genomes. There is no apparent mechanistic impediment to the involvement of the bacterial-type replication machinery in viral DNA replication. Therefore, I hypothesize that the observed, markedly unequal distribution of the replicative DNA polymerases among the known cellular and viral replication systems has a historical explanation. I propose that, among the two types of DNA replication machineries that are found in extant life forms, the archaeal-type, PolB-based system evolved first and had already given rise to a variety of diverse viruses and other selfish elements before the advent of the bacterial, PolC-based machinery. Conceivably, at that stage of evolution, the niches for DNA-viral reproduction have been already filled with viruses replicating with the help of the archaeal system, and viruses with the bacterial system never took off. I further suggest that the two other systems of DNA replication, the rolling circle mechanism and the protein-primed mechanism, which are represented in diverse selfish elements, also evolved prior to the emergence of the bacterial replication system. This hypothesis is compatible with the distinct structural affinities of PolB, which has the palm-domain fold shared with reverse transcriptases and RNA-dependent RNA polymerases, and PolC that has a distinct, unrelated nucleotidyltransferase fold. I propose that PolB is a descendant of polymerases that were involved in the replication of genetic elements in the RNA-protein world, prior to the emergence of DNA replication. By contrast, PolC might have evolved from an ancient non-templated polymerase, e.g., polyA polymerase. The proposed temporal succession of the evolving DNA replication systems does not depend on the specific scenario adopted for the evolution of cells and viruses, i.e., whether viruses are derived from cells or virus-like elements are thought to originate from a primordial gene pool. However, arguments are presented in favor of the latter scenario as the most parsimonious explanation of the evolution of DNA replication systems. CONCLUSION: Comparative analysis of the diversity of genomic strategies and organizations of viruses and cellular life forms has the potential to open windows into the deep past of life's evolution, especially, with the regard to the origin of genome replication systems. When complemented with information on the evolution of the relevant protein folds, this comparative approach can yield credible scenarios for very early steps of evolution that otherwise appear to be out of reach. REVIEWERS: Eric Bapteste, Patrick Forterre, and Mark Ragan.

Journal Article↗
PubMed2006

Compensatory relationship between splice sites and exonic splicing signals depending on the length of vertebrate introns.

BACKGROUND: The signals that determine the specificity and efficiency of splicing are multiple and complex, and are not fully understood. Among other factors, the relative contributions of different mechanisms appear to depend on intron size inasmuch as long introns might hinder the activity of the spliceosome through interference with the proper positioning of the intron-exon junctions. Indeed, it has been shown that the information content of splice sites positively correlates with intron length in the nematode, Drosophila, and fungi. We explored the connections between the length of vertebrate introns, the strength of splice sites, exonic splicing signals, and evolution of flanking exons. RESULTS: A compensatory relationship is shown to exist between different types of signals, namely, the splice sites and the exonic splicing enhancers (ESEs). In the range of relatively short introns (approximately, < 1.5 kilobases in length), the enhancement of the splicing signals for longer introns was manifest in the increased concentration of ESEs. In contrast, for longer introns, this effect was not detectable, and instead, an increase in the strength of the donor and acceptor splice sites was observed. Conceivably, accumulation of A-rich ESE motifs beyond a certain limit is incompatible with functional constraints operating at the level of protein sequence evolution, which leads to compensation in the form of evolution of the splice sites themselves toward greater strength. In addition, however, a correlation between sequence conservation in the exon ends and intron length, particularly, in synonymous positions, was observed throughout the entire length range of introns. Thus, splicing signals other than the currently defined ESEs, i.e., potential new classes of ESEs, might exist in exon sequences, particularly, those that flank long introns. CONCLUSION: Several weak but statistically significant correlations were observed between vertebrate intron length, splice site strength, and potential exonic splicing signals. Taken together, these findings attest to a compensatory relationship between splice sites and exonic splicing signals, depending on intron length.

Animals↗
PubMed2006

Long intervals of stasis punctuated by bursts of positive selection in the seasonal evolution of influenza A virus.

BACKGROUND: The interpandemic evolution of the influenza A virus hemagglutinin (HA) protein is commonly considered a paragon of rapid evolutionary change under positive selection in which amino acid replacements are fixed by virtue of their effect on antigenicity, enabling the virus to evade immune surveillance. RESULTS: We performed phylogenetic analyses of the recently obtained large and relatively unbiased samples of the HA sequences from 1995-2005 isolates of the H3N2 and H1N1 subtypes of influenza A virus. Unexpectedly, it was found that the evolution of H3N2 HA includes long intervals of generally neutral sequence evolution without apparent substantial antigenic change ("stasis" periods) that are characterized by an excess of synonymous over nonsynonymous substitutions per site, lack of association of amino acid replacements with epitope regions, and slow extinction of coexisting virus lineages. These long periods of stasis are punctuated by shorter intervals of rapid evolution under positive selection during which new dominant lineages quickly displace previously coexisting ones. The preponderance of positive selection during intervals of rapid evolution is supported by the dramatic excess of amino acid replacements in the epitope regions of HA compared to replacements in the rest of the HA molecule. In contrast, the stasis intervals showed a much more uniform distribution of replacements over the HA molecule, with a statistically significant difference in the rate of synonymous over nonsynonymous substitution in the epitope regions between the two modes of evolution. A number of parallel amino acid replacements - the same amino acid substitution occurring independently in different lineages - were also detected in H3N2 HA. These parallel mutations were, largely, associated with periods of rapid fitness change, indicating that there are major limitations on evolutionary pathways during antigenic change. The finding that stasis is the prevailing modality of H3N2 evolution suggests that antigenic changes that lead to an increase in fitness typically result from epistatic interactions between several amino acid substitutions in the HA and, perhaps, other viral proteins. The strains that become dominant due to increased fitness emerge from low frequency strains thanks to the last amino acid replacement that completes the set of replacements required to produce a significant antigenic change; no subset of substitutions results in a biologically significant antigenic change and corresponding fitness increase. In contrast to H3N2, no clear intervals of evolution under positive selection were detected for the H1N1 HA during the same time span. Thus, the ascendancy of H1N1 in some seasons is, most likely, caused by the drop in the relative fitness of the previously prevailing H3N2 lineages as the fraction of susceptible hosts decreases during the stasis intervals. CONCLUSION: We show that the common view of the evolution of influenza virus as a rapid, positive selection-driven process is, at best, incomplete. Rather, the interpandemic evolution of influenza appears to consist of extended intervals of stasis, which are characterized by neutral sequence evolution, punctuated by shorter intervals of rapid fitness increase when evolutionary change is driven by positive selection. These observations have implications for influenza surveillance and vaccine formulation; in particular, the possibility exists that parallel amino acid replacements could serve as a predictor of new dominant strains. REVIEWERS: Ron Fouchier (nominated by Andrey Rzhetsky), David Krakauer, Christopher Lee.

Journal Article↗
PubMed2006

Evolution of glyoxylate cycle enzymes in Metazoa: evidence of multiple horizontal transfer events and pseudogene formation.

BACKGROUND: The glyoxylate cycle is thought to be present in bacteria, protists, plants, fungi, and nematodes, but not in other Metazoa. However, activity of the glyoxylate cycle enzymes, malate synthase (MS) and isocitrate lyase (ICL), in animal tissues has been reported. In order to clarify the status of the MS and ICL genes in animals and get an insight into their evolution, we undertook a comparative-genomic study. RESULTS: Using sequence similarity searches, we identified MS genes in arthropods, echinoderms, and vertebrates, including platypus and opossum, but not in the numerous sequenced genomes of placental mammals. The regions of the placental mammals' genomes expected to code for malate synthase, as determined by comparison of the gene orders in vertebrate genomes, show clear similarity to the opossum MS sequence but contain stop codons, indicating that the MS gene became a pseudogene in placental mammals. By contrast, the ICL gene is undetectable in animals other than the nematodes that possess a bifunctional, fused ICL-MS gene. Examination of phylogenetic trees of MS and ICL suggests multiple horizontal gene transfer events that probably went in both directions between several bacterial and eukaryotic lineages. The strongest evidence was obtained for the acquisition of the bifunctional ICL-MS gene from an as yet unknown bacterial source with the corresponding operonic organization by the common ancestor of the nematodes. CONCLUSION: The distribution of the MS and ICL genes in animals suggests that either they encode alternative enzymes of the glyoxylate cycle that are not orthologous to the known MS and ICL or the animal MS acquired a new function that remains to be characterized. Regardless of the ultimate solution to this conundrum, the genes for the glyoxylate cycle enzymes present a remarkable variety of evolutionary events including unusual horizontal gene transfer from bacteria to animals.

Journal Article↗
PubMed2006

Origins and impact of constraints in evolution of gene families.

Recent investigations of high-throughput genomic and phenomic data have uncovered a variety of significant but relatively weak correlations between a gene's functional and evolutionary characteristics. In particular, essential genes and genes with paralogs have a slight propensity to evolve more slowly than nonessential genes and singletons, respectively. However, given the weakness and multiplicity of these associations, their biological relevance remains uncertain. Here, we show that existence of an essential paralog can be used as a specific and strong gauge of selection. We partition gene families in several genomes into two classes: those that include at least one essential gene (E-families) and those without essential genes (N-families). We find that weaker purifying selection causes N-families to evolve in a more dynamic regime with higher rates both of duplicate fixation and pseudogenization. Because genes in E-families are subject to significantly stronger purifying selection than those in N-families, they survive longer and exhibit greater sequence divergence. Longer average survival time also allows for divergence of upstream regulatory regions, resulting in change of transcriptional context among paralogs in E-families. These findings are compatible with differential division of ancestral functions (subfunctionalization) or emergence of novel functions (neofunctionalization) being the prevalent modes of evolution of paralogs in E-families as opposed to pseudogenization (nonfunctionalization), which is the typical fate of paralogs in N-families. Unlike other characteristics of genes, such as essentiality, existence of paralogs, or expression level, membership in an E-family or an N-family strongly correlates with the level of selection and appears to be a major determinant of a gene's evolutionary fate.

Animals↗
PubMed2006

Kinase activity of overexpressed HipA is required for growth arrest and multidrug tolerance in Escherichia coli.

Overexpression of the HipA protein of the HipBA toxin/antitoxin module leads to multidrug tolerance in Escherichia coli. HipA is a "toxin" that causes reversible dormancy, whereas HipB is an antitoxin that binds HipA and acts as a transcriptional repressor of the hipBA operon. Comparative sequence analysis shows that HipA is a member of the phosphatidylinositol 3/4-kinase superfamily. The kinase activity of HipA was examined. HipA was autophosphorylated in the presence of ATP in vitro, and the purified protein appeared to carry a single phosphate group on serine 150. Thus, HipA is a serine kinase that is at least partially phosphorylated in vivo. Overexpression of HipA caused inhibition of cell growth and increase in persister formation. Replacing conserved aspartate 309 in the conserved kinase active site or aspartate 332 in the Mg2+-binding site with glutamine produced mutant proteins that lost the ability to stop cellular growth upon overexpression. Replacing serine 150 with alanine yielded a similarly inactive protein. The mutant proteins were then examined for their ability to increase antibiotic tolerance. Cells overexpressing wild-type HipA were highly tolerant to cefotaxime, a cell wall synthesis inhibitor, to ofloxacin, a fluoroquinolone inhibitor of DNA gyrase, and to topoisomerase IV and were almost completely resistant to killing by mitomycin C, which forms DNA adducts. The mutant proteins did not protect cells from cefotaxime or ofloxacin and had an impaired ability to protect from mitomycin C. Taken together, these results suggest that the protein kinase activity of HipA is essential for persister formation.

Amino Acid Sequence↗
PubMed2006

Mathematical modeling of tumor therapy with oncolytic viruses: effects of parametric heterogeneity on cell dynamics.

BACKGROUND: One of the mechanisms that ensure cancer robustness is tumor heterogeneity, and its effects on tumor cells dynamics have to be taken into account when studying cancer progression. There is no unifying theoretical framework in mathematical modeling of carcinogenesis that would account for parametric heterogeneity. RESULTS: Here we formulate a modeling approach that naturally takes stock of inherent cancer cell heterogeneity and illustrate it with a model of interaction between a tumor and an oncolytic virus. We show that several phenomena that are absent in homogeneous models, such as cancer recurrence, tumor dormancy, and others, appear in heterogeneous setting. We also demonstrate that, within the applied modeling framework, to overcome the adverse effect of tumor cell heterogeneity on the outcome of cancer treatment, a heterogeneous population of an oncolytic virus must be used. Heterogeneity in parameters of the model, such as tumor cell susceptibility to virus infection and the ability of an oncolytic virus to infect tumor cells, can lead to complex, irregular evolution of the tumor. Thus, quasi-chaotic behavior of the tumor-virus system can be caused not only by random perturbations but also by the heterogeneity of the tumor and the virus. CONCLUSION: The modeling approach described here reveals the importance of tumor cell and virus heterogeneity for the outcome of cancer therapy. It should be straightforward to apply these techniques to mathematical modeling of other types of anticancer therapy. REVIEWERS: Leonid Hanin (nominated by Arcady Mushegian), Natalia Komarova (nominated by Orly Alter), and David Krakauer.

Journal Article↗
PubMed2006

Genome-wide analysis of substrate specificities of the Escherichia coli haloacid dehalogenase-like phosphatase family.

Haloacid dehalogenase (HAD)-like hydrolases are a vast superfamily of largely uncharacterized enzymes, with a few members shown to possess phosphatase, beta-phosphoglucomutase, phosphonatase, and dehalogenase activities. Using a representative set of 80 phosphorylated substrates, we characterized the substrate specificities of 23 soluble HADs encoded in the Escherichia coli genome. We identified small molecule phosphatase activity in 21 HADs and beta-phosphoglucomutase activity in one protein. The E. coli HAD phosphatases show high catalytic efficiency and affinity to a wide range of phosphorylated metabolites that are intermediates of various metabolic reactions. Rather than following the classical "one enzyme-one substrate" model, most of the E. coli HADs show remarkably broad and overlapping substrate spectra. At least 12 reactions catalyzed by HADs currently have no EC numbers assigned in Enzyme Nomenclature. Surprisingly, most HADs hydrolyzed small phosphodonors (acetyl phosphate, carbamoyl phosphate, and phosphoramidate), which also serve as substrates for autophosphorylation of the receiver domains of the two-component signal transduction systems. The physiological relevance of the phosphatase activity with the preferred substrate was validated in vivo for one of the HADs, YniC. Many of the secondary activities of HADs might have no immediate physiological function but could comprise a reservoir for evolution of novel phosphatases.

Catalysis↗
PubMed2006

The ancient Virus World and evolution of cells.

BACKGROUND: Recent advances in genomics of viruses and cellular life forms have greatly stimulated interest in the origins and evolution of viruses and, for the first time, offer an opportunity for a data-driven exploration of the deepest roots of viruses. Here we briefly review the current views of virus evolution and propose a new, coherent scenario that appears to be best compatible with comparative-genomic data and is naturally linked to models of cellular evolution that, from independent considerations, seem to be the most parsimonious among the existing ones. RESULTS: Several genes coding for key proteins involved in viral replication and morphogenesis as well as the major capsid protein of icosahedral virions are shared by many groups of RNA and DNA viruses but are missing in cellular life forms. On the basis of this key observation and the data on extensive genetic exchange between diverse viruses, we propose the concept of the ancient virus world. The virus world is construed as a distinct contingent of viral genes that continuously retained its identity throughout the entire history of life. Under this concept, the principal lineages of viruses and related selfish agents emerged from the primordial pool of primitive genetic elements, the ancestors of both cellular and viral genes. Thus, notwithstanding the numerous gene exchanges and acquisitions attributed to later stages of evolution, most, if not all, modern viruses and other selfish agents are inferred to descend from elements that belonged to the primordial genetic pool. In this pool, RNA viruses would evolve first, followed by retroid elements, and DNA viruses. The Virus World concept is predicated on a model of early evolution whereby emergence of substantial genetic diversity antedates the advent of full-fledged cells, allowing for extensive gene mixing at this early stage of evolution. We outline a scenario of the origin of the main classes of viruses in conjunction with a specific model of precellular evolution under which the primordial gene pool dwelled in a network of inorganic compartments. Somewhat paradoxically, under this scenario, we surmise that selfish genetic elements ancestral to viruses evolved prior to typical cells, to become intracellular parasites once bacteria and archaea arrived at the scene. Selection against excessively aggressive parasites that would kill off the host ensembles of genetic elements would lead to early evolution of temperate virus-like agents and primitive defense mechanisms, possibly, based on the RNA interference principle. The emergence of the eukaryotic cell is construed as the second melting pot of virus evolution from which the major groups of eukaryotic viruses originated as a result of extensive recombination of genes from various bacteriophages, archaeal viruses, plasmids, and the evolving eukaryotic genomes. Again, this vision is predicated on a specific model of the emergence of eukaryotic cell under which archaeo-bacterial symbiosis was the starting point of eukaryogenesis, a scenario that appears to be best compatible with the data. CONCLUSION: The existence of several genes that are central to virus replication and structure, are shared by a broad variety of viruses but are missing from cellular genomes (virus hallmark genes) suggests the model of an ancient virus world, a flow of virus-specific genes that went uninterrupted from the precellular stage of life's evolution to this day. This concept is tightly linked to two key conjectures on evolution of cells: existence of a complex, precellular, compartmentalized but extensively mixing and recombining pool of genes, and origin of the eukaryotic cell by archaeo-bacterial fusion. The virus world concept and these models of major transitions in the evolution of cells provide complementary pieces of an emerging coherent picture of life's history. REVIEWERS: W. Ford Doolittle, J. Peter Gogarten, and Arcady Mushegian.

Journal Article↗
PubMed2006

Global similarity and local divergence in human and mouse gene co-expression networks.

BACKGROUND: A genome-wide comparative analysis of human and mouse gene expression patterns was performed in order to evaluate the evolutionary divergence of mammalian gene expression. Tissue-specific expression profiles were analyzed for 9,105 human-mouse orthologous gene pairs across 28 tissues. Expression profiles were resolved into species-specific coexpression networks, and the topological properties of the networks were compared between species. RESULTS: At the global level, the topological properties of the human and mouse gene coexpression networks are, essentially, identical. For instance, both networks have topologies with small-world and scale-free properties as well as closely similar average node degrees, clustering coefficients, and path lengths. However, the human and mouse coexpression networks are highly divergent at the local level: only a small fraction (<10%) of coexpressed gene pair relationships are conserved between the two species. A series of controls for experimental and biological variance show that most of this divergence does not result from experimental noise. We further show that, while the expression divergence between species is genuinely rapid, expression does not evolve free from selective (functional) constraint. Indeed, the coexpression networks analyzed here are demonstrably functionally coherent as indicated by the functional similarity of coexpressed gene pairs, and this pattern is most pronounced in the conserved human-mouse intersection network. Numerous dense network clusters show evidence of dedicated functions, such as spermatogenesis and immune response, that are clearly consistent with the coherence of the expression patterns of their constituent gene members. CONCLUSION: The dissonance between global versus local network divergence suggests that the interspecies similarity of the global network properties is of limited biological significance, at best, and that the biologically relevant aspects of the architectures of gene coexpression are specific and particular, rather than universal. Nevertheless, there is substantial evolutionary conservation of the local network structure which is compatible with the notion that gene coexpression networks are subject to purifying selection.

Animals↗
PubMed2006

Evolutionary systems biology: links between gene evolution and function.

The recent accumulation of genome-wide data on various facets of gene expression, function and evolution stimulated the emergence of a new field, evolutionary systems biology. Many significant correlations were detected between variables that characterize the functioning of a gene, such as expression level, knockout effect, connectivity of genetic and protein-protein interaction networks, and variables that describe gene evolution, such as sequence evolution rate and propensity for gene loss. The first attempts on multidimensional analysis of genomic data yielded composite variables that describe the 'status' of a gene in the genomic community. However, it remains uncertain whether different functional variables affect gene evolution synergistically or there is a single, dominant factor. The number of translation events, linked to selection for translational robustness, was proposed as a candidate for such a major determinant of protein evolution. These developments show that, although the methodological basis of evolutionary systems biology is not yet fully solidified, this area of research is already starting to yield fundamental biological insights.

Evolution, Molecular↗