Higher Gene Duplicabilities for Metabolic Proteins Than for Nonmetabolic Proteins in Yeast and <Emphasis Type="Italic">E. coli</Emphasis>

Although the evolutionary significance of gene duplication has long been appreciated, it remains unclear what factors determine gene duplicability. In this study we investigated whether metabolism is an important determinant of gene duplicability because cellular metabolism is crucial for the survival and reproduction of an organism. Using genomic data and metabolic pathway data from the yeast (Saccharomyces cerevisiae) and Escherichia coli, we found that metabolic proteins indeed tend to have higher gene duplicability than nonmetabolic proteins. Moreover, a detailed analysis of metabolic pathways in these two organisms revealed that genes in the central metabolic pathways and the catabolic pathways have, on average, higher gene duplicability than do other genes and that most genes in anabolic pathways are single-copy genes.Reviewing Editor: Dr. Rüdiger Cerff     

Gene complexity and gene duplicability

Eukaryotic genes are on average more complex than prokaryotic genes in terms of expression regulation, protein length, and protein-domain structure [1-5]. Eukaryotes are also known to have a higher rate of gene duplication than prokaryotes do [6, 7]. Because gene duplication is the primary source of new genes [], the average gene complexity in a genome may have been increased by gene duplication if complex genes are preferentially duplicated. Here, we test this "gene complexity and gene duplicability" hypothesis with yeast genomic data. We show that, on average, duplicate genes from either whole-genome or individual-gene duplication have longer protein sequences, more functional domains, and more cis-regulatory motifs than singleton genes. This phenomenon is not a by-product of previously known mechanisms, such as protein function [10-13], evolutionary rate [14, 15], dosage [11], and dosage balance [16], that influence gene duplicability. Rather, it appears to have resulted from the sub-neo-functionalization process in duplicate-gene evolution [11]. Under this process, complex genes are more likely to be retained after duplication because they are prone to subfunctionalization, and gene complexity is regained via subsequent neofunctionalization. Thus, gene duplication increases both gene number and gene complexity, two important factors in the origin of genomic and organismal complexity.     

Protein complexity, gene duplicability and gene dispensability in the yeast genome

Using functional genomic and protein structural data we studied the effects of protein complexity (here defined as the number of subunit types in a protein) on gene dispensability and gene duplicability. We found that in terms of gene duplicability the major distinction in protein complexity is between hetero-complexes, each of which includes at least two different types of subunits (polypeptides), and homo-complexes, which include monomers and complexes that consist of only subunits of one polypeptide type. However, gene dispensability decreases only gradually as the number of subunit types in a protein complex increases. These observations suggest that the dosage balance hypothesis can explain well gene duplicability of complex proteins, but cannot completely explain the difference in dispensabilities between hetero-complex subunits. It is likely that knocking out a gene coding for a hetero-complex subunit would disrupt the function of the whole complex, so that the deletion effect on fitness would increase with protein complexity. We also found that multi-domain polypeptide genes are less dispensable but more duplicable than single-domain polypeptide genes. Duplicate genes derived from the whole genome duplication event in yeast are more dispensable (except for ribosomal protein genes) than other duplicate genes. Further, we found that subunits of the same protein complex tend to have similar expression levels and similar effects of gene deletion on fitness. Finally, we estimated that in yeast the contribution of duplicate genes to genetic robustness against null mutation is approximately 9%, smaller than previously estimated. In yeast, protein complexity may serve as a better indicator of gene dispensability than do duplicate genes.     

Protein Connectivity and Protein Complexity Promotes Human Gene Duplicability in a Mutually Exclusive Manner

It has previously been reported that protein complexity (i.e. number of subunits in a protein complex) is negatively correlated to gene duplicability in yeast as well as in humans. However, unlike in yeast, protein connectivity in a protein–protein interaction network has a positive correlation with gene duplicability in human genes. In the present study, we have analyzed 1732 human and 1269 yeast proteins that are present both in a protein–protein interaction network as well as in a protein complex network. In the human case, we observed that both protein connectivity and protein complexity complement each other in a mutually exclusive manner over gene duplicability in a positive direction. Analysis of human haploinsufficient proteins and large protein complexes (complex size >10) shows that when protein connectivity does not have any direct association with gene duplicability, there exists a positive correlation between gene duplicability and protein complexity. The same trend, however, is not found in case of yeast, where both protein connectivity and protein complexity independently guide gene duplicability in the negative direction. We conclude that the higher rate of duplication of human genes may be attributed to organismal complexity either by increasing connectivity in the protein–protein interaction network or by increasing protein complexity.     

Complexity, connectivity, and duplicability as barriers to lateral gene transfer

Background

Lateral gene transfer is a major force in microbial evolution and a great source of genetic innovation in prokaryotes. Protein complexity has been claimed to be a barrier for gene transfer, due to either the inability of a new gene's encoded protein to become a subunit of an existing complex (lack of positive selection), or from a harmful effect exerted by the newcomer on native protein assemblages (negative selection).

Results

We tested these scenarios using data from the model prokaryote Escherichia coli. Surprisingly, the data did not support an inverse link between membership in protein complexes and gene transfer. As the complexity hypothesis, in its strictest sense, seemed valid only to essential complexes, we broadened its scope to include connectivity in general. Transferred genes are found to be less involved in protein-protein interactions, outside stable complexes, and this is especially true for genes recently transferred to the E. coli genome. Thus, subsequent to transfer, new genes probably integrate slowly into existing protein-interaction networks. We show that a low duplicability of a gene is linked to a lower chance of being horizontally transferred. Notably, many essential genes in E. coli are conserved as singletons across multiple related genomes, have high connectivity and a highly vertical phylogenetic signal.

Conclusion

High complexity and connectivity generally do not impede gene transfer. However, essential genes that exhibit low duplicability and high connectivity do exhibit mostly vertical descent.     

Relationship between gene duplicability and diversifiability in the topology of biochemical networks

Background

Selective gene duplicability, the extensive expansion of a small number of gene families, is universal. Quantitatively, the number of genes (P_(K)) with K duplicates in a genome decreases precipitously as K increases, and often follows a power law (P_(k)∝k^-α). Functional diversification, either neo- or sub-functionalization, is a major evolution route for duplicate genes.

Results

Using three lines of genomic datasets, we studied the relationship between gene duplicability and diversifiability in the topology of biochemical networks. First, we explored scenario where two pathways in the biochemical networks antagonize each other. Synthetic knockout of respective genes for the two pathways rescues the phenotypic defects of each individual knockout. We identified duplicate gene pairs with sufficient divergences that represent this antagonism relationship in the yeast S. cerevisiae. Such pairs overwhelmingly belong to large gene families, thus tend to have high duplicability. Second, we used distances between proteins of duplicate genes in the protein interaction network as a metric of their diversification. The higher a gene’s duplicate count, the further the proteins of this gene and its duplicates drift away from one another in the networks, which is especially true for genetically antagonizing duplicate genes. Third, we computed a sequence-homology-based clustering coefficient to quantify sequence diversifiability among duplicate genes – the lower the coefficient, the more the sequences have diverged. Duplicate count (K) of a gene is negatively correlated to the clustering coefficient of its duplicates, suggesting that gene duplicability is related to the extent of sequence divergence within the duplicate gene family.

Conclusion

Thus, a positive correlation exists between gene diversifiability and duplicability in the context of biochemical networks – an improvement of our understanding of gene duplicability.     

Trehalose metabolism genes in Caenorhabditis elegans and filarial nematodes

The sugar trehalose is claimed to be important in the physiology of nematodes where it may function in sugar transport, energy storage and protection against environmental stresses. In this study we investigated the role of trehalose metabolism in nematodes, using Caenorhabditis elegans as a model, and also identified complementary DNA clones putatively encoding genes involved in trehalose pathways in filarial nematodes. In C. elegans two putative trehalose-6-phosphate synthase (tps) genes encode the enzymes that catalyse trehalose synthesis and five putative trehalase (tre) genes encode enzymes catalysing hydrolysis of the sugar. We showed by RT-PCR or Northern analysis that each of these genes is expressed as mRNA at all stages of the C. elegans life cycle. Database searches and sequencing of expressed sequence tag clones revealed that at least one tps gene and two tre genes are expressed in the filarial nematode Brugia malayi, while one tps gene and at least one tre gene were identified for Onchocerca volvulus. We used the feeding method of RNA interference in C. elegans to knock down temporarily the expression of each of the tps and tre genes. Semiquantitative RT-PCR analysis confirmed that expression of each gene was silenced by RNA interference. We did not observe an obvious phenotype for any of the genes silenced individually but gas-chromatographic analysis showed >90% decline in trehalose levels when both tps genes were targeted simultaneously. This decline in trehalose content did not affect viability or development of the nematodes.     

Gene Duplicability of Core Genes Is Highly Consistent across All Angiosperms

Gene duplication is an important mechanism for adding to genomic novelty. Hence, which genes undergo duplication and are preserved following duplication is an important question. It has been observed that gene duplicability, or the ability of genes to be retained following duplication, is a nonrandom process, with certain genes being more amenable to survive duplication events than others. Primarily, gene essentiality and the type of duplication (small-scale versus large-scale) have been shown in different species to influence the (long-term) survival of novel genes. However, an overarching view of “gene duplicability” is lacking, mainly due to the fact that previous studies usually focused on individual species and did not account for the influence of genomic context and the time of duplication. Here, we present a large-scale study in which we investigated duplicate retention for 9178 gene families shared between 37 flowering plant species, referred to as angiosperm core gene families. For most gene families, we observe a strikingly consistent pattern of gene duplicability across species, with gene families being either primarily single-copy or multicopy in all species. An intermediate class contains gene families that are often retained in duplicate for periods extending to tens of millions of years after whole-genome duplication, but ultimately appear to be largely restored to singleton status, suggesting that these genes may be dosage balance sensitive. The distinction between single-copy and multicopy gene families is reflected in their functional annotation, with single-copy genes being mainly involved in the maintenance of genome stability and organelle function and multicopy genes in signaling, transport, and metabolism. The intermediate class was overrepresented in regulatory genes, further suggesting that these represent putative dosage-balance-sensitive genes.     

The Enrichment of TATA Box and the Scarcity of Depleted Proximal Nucleosome in the Promoters of Duplicated Yeast Genes

Population genetic theory of gene duplication suggests that the preservation of duplicate copies requires functional divergence upon duplication. Genes that can be readily modified to produce new gene expression patterns may thus be duplicated often. In yeast, genes exhibit dichotomous expression patterns based on their promoter architectures. The expression of genes that contain TATA box or occupied proximal nucleosome (OPN) tends to be variable and respond to external signals. On the other hand, genes without TATA box or with depleted proximal nucleosome (DPN) are expressed constitutively. We find that recent duplicates in the yeast genome are heavily biased to be TATA box containing genes and not to be DPN genes. This suggests that variably expressed genes, due to the functional organization in their promoters, have higher duplicability than constitutively expressed genes.     

Evidence of selection at melanin synthesis pathway loci during silkworm domestication

The domesticated silkworm (Bombyx mori) was domesticated from wild silkworm (Bombyx mandarina) more than 5,000 years ago. During domestication, body color between B. mandarina and B. mori changed dramatically. However, the molecular mechanism of the silkworm body color transition is not known. In the present study, we examined within- and between-species nucleotide diversity for eight silkworm melanin synthesis pathway genes, which play a key role in cuticular pigmentation of insects. Our results showed that the genetic diversity of B. mori was significantly lower than that of B. mandarina and 40.7% of the genetic diversity of wild silkworm was lost in domesticated silkworm. We also examined whether position effect exists among melanin synthesis pathway genes in B. mandarina and B. mori. We found that the upstream genes have significantly lower levels of genetic diversity than the downstream genes, supporting a functional constraint hypothesis (FCH) of metabolic pathway, that is, upstream enzymes are under greater selective constraint than downstream enzymes because upstream enzymes participate in biosynthesis of a number of metabolites. We also investigated whether some of the melanin synthesis pathway genes experienced selection during domestication. Neutrality test, coalescent simulation, as well as network and phylogenetic analyses showed that tyrosine hydroxylase (TH) gene was a domestication locus. Sequence analysis further suggested that a putative expression enhancer (Abd-B-binding site) in the intron of TH gene might be disrupted during domestication. TH is the rate-limiting enzyme of melanin synthesis pathway in insects. Real-time polymerase chain reaction assay did show that the relative expression levels of TH gene in B. mori were significantly lower than that in B. mandarina at three different developmental stages, which is consistent with light body color of domesticated silkworm relative to wild silkworm. Therefore, we speculated that expression change of TH gene may contribute to the body color transition from B. mandarina to B. mori. Our results emphasize the exceptional role of gene expression regulation in morphological transition of domesticated animals.     

Differential gene expression in Arabidopsis following infection by plant-parasitic nematodes Meloidogyne incognita and Heterodera schachtii

Whole genome microarrays were used to study plant gene expression in mature Meloidogyne incognita -induced galls in Arabidopsis. We found 959 genes to be significantly differentially expressed, and two-thirds of these were down-regulated. Microarray results were confirmed by qRT-PCR. The temporal and spatial responses of four differentially expressed genes were analysed using GUS reporter plants following infection with M. incognita and the cyst nematode Heterodera schachtii . The ammonium transporter gene AtAMT1;2 was consistently and locally repressed in response to both nematodes at all developmental stages. The lateral organ boundary domain gene LBD41 showed up-regulation in the feeding sites of both nematode species, although there was variation in expression in saccate H. schachtii female feeding sites. Expression of an actin depolymerizing factor ADF3 and a lipid transfer protein was induced in feeding sites of both nematodes at the fusiform stage and this persisted in feeding sites of saccate M. incognita . These results contribute to the knowledge of how plant gene expression responds to parasitism by these nematodes as well as highlighting further differences in the mechanisms of development and maintenance of these feeding site structures.     

GAP-43 promoter elements in transgenic zebrafish reveal a difference in signals for axon growth during CNS development and regeneration

A pivotal event in neural development is the point at which differentiating neurons become competent to extend long axons. Initiation of axon growth is equally critical for regeneration. Yet we have a limited understanding of the signaling pathways that regulate the capacity for axon growth during either development or regeneration. Expression of a number of genes encoding growth associated proteins (GAPs) accompanies both developmental and regenerative axon growth and has led to the suggestion that the same signaling pathways regulate both modes of axon growth. We have tested this possibility by asking whether a promoter fragment from a well characterized GAP gene, GAP-43, is sufficient to activate expression in both developing and regenerating neurons. We generated stable lines of transgenic zebrafish that express green fluorescent protein (GFP) under regulation of a 1 kb fragment of the rat GAP-43 gene, a fragment that contains a number of evolutionarily conserved elements. Analysis of GFP expression in these lines confirms that the rat 1 kb region can direct growth-associated expression of the transgene in differentiating neurons that extend long axons. Furthermore, this region supports developmental down-regulation of transgene expression which, like the endogenous gene, coincides with neuronal maturation. Strikingly, these same sequences are insufficient for directing expression in regenerating neurons. This finding suggests that signaling pathways regulating axon growth during development and regeneration are not the same. While these results do not exclude the possibility that pathways involved in developmental axon growth are also active in regenerative growth, they do indicate that signaling pathway(s) controlling activation of the GAP-43 gene after CNS injury differ in at least one key component from the signals controlling essential features of developmental axon growth.     

Oesophagostomum dentatum: potential as a model for genomic studies of strongylid nematodes, with biotechnological prospects

There are substantial gaps in the knowledge of the molecular processes of development and reproduction in parasitic nematodes, despite the fact that understanding such processes could lead to novel ways of treating and controlling parasitic diseases, through blocking or disrupting key biological pathways. Biotechnological advances through large-scale sequencing projects, approaches for the analysis of differential gene and protein expression and functional genomics (e.g., double-stranded RNA interference) now provide opportunities to investigate the molecular basis of developmental processes in some parasitic nematodes. The porcine nodule worm, Oesophagostomum dentatum (order Strongylida), may provide a platform for testing the function of genes from this and related nematodes, given that this species can be grown and maintained in culture in vitro for periods longer than other nematodes of the same order. In this article, we review relevant biological, biochemical and molecular biological and genomic information about O. dentatum and propose that the O. dentatum - pig system provides an attractive model for exploring molecular developmental and reproductive processes in strongylid nematodes, leading toward new intervention methods and biotechnological outcomes.     

Evolutionary significance of gene expression divergence

Recent large-scale studies of evolutionary changes in gene expression among mammalian species have led to the proposal that gene expression divergence may be neutral with respect to organismic fitness. Here, we employ a comparative analysis of mammalian gene sequence divergence and gene expression divergence to test the hypothesis that the evolution of gene expression is predominantly neutral. Two models of neutral gene expression evolution are considered: 1-purely neutral evolution (i.e., no selective constraint) of gene expression levels and patterns and 2-neutral evolution accompanied by selective constraint. With respect to purely neutral evolution, levels of change in gene expression between human-mouse orthologs are correlated with levels of gene sequence divergence that are determined largely by purifying selection. In contrast, evolutionary changes of tissue-specific gene expression profiles do not show such a correlation with sequence divergence. However, divergence of both gene expression levels and profiles are significantly lower for orthologous human-mouse gene pairs than for pairs of randomly chosen human and mouse genes. These data clearly point to the action of selective constraint on gene expression divergence and are inconsistent with the purely neutral model; however, there is likely to be a neutral component in evolution of gene expression, particularly, in tissues where the expression of a given gene is low and functionally irrelevant. The model of neutral evolution with selective constraint predicts a regular, clock-like accumulation of gene expression divergence. However, relative rate tests of the divergence among human-mouse-rat orthologous gene sets reveal clock-like evolution for gene sequence divergence, and to a lesser extent for gene expression level divergence, but not for the divergence of tissue-specific gene expression profiles. Taken together, these results indicate that gene expression divergence is subject to the effects of purifying selective constraint and suggest that it might also be substantially influenced by positive Darwinian selection.