Strong evolutionary conservation of broadly expressed protein isoforms in the troponin I gene family and other vertebrate gene families

Functional buffering via cell-specific gene expression promotes tissue homeostasis and cancer robustness

Functional buffering that ensures biological robustness is critical for maintaining tissue homeostasis, organismal survival, and evolution of novelty. However, the mechanism underlying functional buffering, particularly in multicellular organisms, remains largely elusive. Here, we proposed that functional buffering can be mediated via expression of buffering genes in specific cells and tissues, by which we named Cell-specific Expression-BUffering (CEBU). We developed an inference index (C-score) for CEBU by computing C-scores across 684 human cell lines using genome-wide CRISPR screens and transcriptomic RNA-seq. We report that C-score-identified putative buffering gene pairs are enriched for members of the same duplicated gene family, pathway, and protein complex. Furthermore, CEBU is especially prevalent in tissues of low regenerative capacity (e.g., bone and neuronal tissues) and is weakest in highly regenerative blood cells, linking functional buffering to tissue regeneration. Clinically, the buffering capacity enabled by CEBU can help predict patient survival for multiple cancers. Our results suggest CEBU as a potential buffering mechanism contributing to tissue homeostasis and cancer robustness in humans.

Mapping Selection within Drosophila melanogaster Embryo's Anatomy

We present a survey of selection across Drosophila melanogaster embryonic anatomy. Our approach integrates genomic variation, spatial gene expression patterns, and development with the aim of mapping adaptation over the entire embryo's anatomy. Our adaptation map is based on analyzing spatial gene expression information for 5,969 genes (from text-based annotations of in situ hybridization data directly from the BDGP database, Tomancak et al. 2007) and the polymorphism and divergence in these genes (from the project DGRP, Mackay et al. 2012).The proportion of nonsynonymous substitutions that are adaptive, neutral, or slightly deleterious are estimated for the set of genes expressed in each embryonic anatomical structure using the distribution of fitness effects-alpha method (Eyre-Walker and Keightley 2009). This method is a robust derivative of the McDonald and Kreitman test (McDonald and Kreitman 1991). We also explore whether different anatomical structures differ in the phylogenetic age, codon usage, or expression bias of the genes they express and whether genes expressed in many anatomical structures show more adaptive substitutions than other genes.We found that: 1) most of the digestive system and ectoderm-derived structures are under selective constraint, 2) the germ line and some specific mesoderm-derived structures show high rates of adaptive substitution, and 3) the genes that are expressed in a small number of anatomical structures show higher expression bias, lower phylogenetic ages, and less constraint.

Evolutionary Consequences of DNA Methylation in a Basal Metazoan

Gene body methylation (gbM) is an ancestral and widespread feature in Eukarya, yet its adaptive value and evolutionary implications remain unresolved. The occurrence of gbM within protein-coding sequences is particularly puzzling, because methylation causes cytosine hypermutability and hence is likely to produce deleterious amino acid substitutions. We investigate this enigma using an evolutionarily basal group of Metazoa, the stony corals (order Scleractinia, class Anthozoa, phylum Cnidaria). We show that patterns of coral gbM are similar to other invertebrate species, predicting wide and active transcription and slower sequence evolution. We also find a strong correlation between gbM and codon bias, resulting from systematic replacement of CpG bearing codons. We conclude that gbM has strong effects on codon evolution and speculate that this may influence establishment of optimal codons.

Sexual dimorphism and the evolution of sex-biased gene expression in the brown alga ectocarpus

Males and females often have marked phenotypic differences, and the expression of these dissimilarities invariably involves sex differences in gene expression. Sex-biased gene expression has been well characterized in animal species, where a high proportion of the genome may be differentially regulated in males and females during development. Male-biased genes tend to evolve more rapidly than female-biased genes, implying differences in the strength of the selective forces acting on the two sexes. Analyses of sex-biased gene expression have focused on organisms that exhibit separate sexes during the diploid phase of the life cycle (diploid sexual systems), but the genetic nature of the sexual system is expected to influence the evolutionary trajectories of sex-biased genes. We analyze here the patterns of sex-biased gene expression in Ectocarpus, a brown alga with haploid sex determination (dioicy) and a low level of phenotypic sexual dimorphism. In Ectocarpus, female-biased genes were found to be evolving as rapidly as male-biased genes. Moreover, genes expressed at fertility showed faster rates of evolution than genes expressed in immature gametophytes. Both male- and female-biased genes had a greater proportion of sites experiencing positive selection, suggesting that their accelerated evolution is at least partly driven by adaptive evolution. Gene duplication appears to have played a significant role in the generation of sex-biased genes in Ectocarpus, expanding previous models that propose this mechanism for the resolution of sexual antagonism in diploid systems. The patterns of sex-biased gene expression in Ectocarpus are consistent both with predicted characteristics of UV (haploid) sexual systems and with the distinctive aspects of this organism's reproductive biology.

Pleiotropic mutations are subject to strong stabilizing selection

The assumption that pleiotropic mutations are more deleterious than mutations with more restricted phenotypic effects is an important premise in models of evolution. However, empirical evidence supporting this assumption is limited. Here, we estimated the strength of stabilizing selection on mutations affecting gene expression in male Drosophila serrata. We estimated the mutational variance (VM) and the standing genetic variance (VG) from two well-matched panels of inbred lines: a panel of mutation accumulation (MA) lines derived from a single inbred ancestral line and a panel of inbred lines derived from an outbred population. For 855 gene-expression traits, we estimated the strength of stabilizing selection as s = VM/VG. Selection was observed to be relatively strong, with 17% of traits having s > 0.02, a magnitude typically associated with life-history traits. Randomly assigning expression traits to five-trait sets, we used factor analytic mixed modeling in the MA data set to identify covarying traits that shared pleiotropic mutations. By assigning traits to the same trait sets in the outbred line data set, we then estimated s for the combination of traits affected by pleiotropic mutation. For these pleiotropic combinations, the median s was three times greater than s acting on the individual component traits, and 46% of the pleiotropic trait combinations had s > 0.02. Although our analytical approach was biased toward detecting mutations with relatively large effects, likely overestimating the average strength of selection, our results provide widespread support for the prediction that stronger selection can act against mutations with pleiotropic effects.

Nested insertions and accumulation of indels are negatively correlated with abundance of mutator-like transposable elements in maize and rice

Mutator-like transposable elements (MULEs) are widespread in plants and were first discovered in maize where there are a total of 12,900 MULEs. In comparison, rice, with a much smaller genome, harbors over 30,000 MULEs. Since maize and rice are close relatives, the differential amplification of MULEs raised an inquiry into the underlying mechanism. We hypothesize this is partly attributed to the differential copy number of autonomous MULEs with the potential to generate the transposase that is required for transposition. To this end, we mined the two genomes and detected 530 and 476 MULEs containing transposase sequences (candidate coding-MULEs) in maize and rice, respectively. Over 1/3 of the candidate coding-MULEs harbor nested insertions and the ratios are similar in the two genomes. Among the maize elements with nested insertions, 24% have insertions in coding regions and over half of them harbor two or more insertions. In contrast, only 12% of the rice elements have insertions in coding regions and 19% have multiple insertions, suggesting that nested insertions in maize are more disruptive. This is because most nested insertions in maize are from LTR retrotransposons, which are large in size and are prevalent in the maize genome. Our results suggest that the amplification of retrotransposons may limit the amplification of DNA transposons but not vice versa. In addition, more indels are detected among maize elements than rice elements whereas defects caused by point mutations are comparable between the two species. Taken together, more disruptive nested insertions combined with higher frequency of indels resulted in few (6%) coding-MULEs that may encode functional transposases in maize. In contrast, 35% of the coding-MULEs in rice retain putative intact transposase. This is in addition to the higher expression frequency of rice coding-MULEs, which may explain the higher occurrence of MULEs in rice than that in maize.

Assessing determinants of exonic evolutionary rates in mammals

From studies investigating the differences in evolutionary rates between genes, gene compactness and gene expression level have been identified as important determinants of gene-level protein evolutionary rate, as represented by nonsynonymous to synonymous substitution rate (d(N)/d(S)) ratio. However, the causes of exon-level variances in d(N)/d(S) are less understood. Here, we use principal component regression to examine to what extent 13 exon features explain the variance in d(N), d(S), and the d(N)/d(S) ratio of human-rhesus macaque or human-mouse orthologous exons. The exon features were grouped into six functional categories: expression features, mRNA splicing features, structural-functional features, compactness features, exon duplicability, and other features, including G + C content and exon length. Although expression features are important for determining d(N) and d(N)/d(S) between exons of different genes, structural-functional features and splicing features explained more of the variance for exons of the same genes. Furthermore, we show that compactness features can explain only a relatively small percentage of variance in exon-level d(N) or d(N)/d(S) in either between-gene or within-gene comparison. By contrast, d(S) yielded inconsistent results in the human-mouse comparison and the human-rhesus macaque comparison. This inconsistency may suggest rapid evolutionary changes of the mutation landscape in mammals. Our results suggest that between-gene and within-gene variation in d(N)/d(S) (and d(N)) are driven by different evolutionary forces and that the role of mRNA splicing in causing the variation in evolutionary rates of coding sequences may be underappreciated.

Evolution of vertebrate tissues driven by differential modes of gene duplication

In this study, we investigated the evolution of vertebrate tissues by examining the potential association among gene expression, duplication, and base substitution patterns. In particular, we compared whole-genome duplication (WGD) with small-scale duplication (SSD), as well as tissue restricted with ubiquitously expressed genes. All patterns were also analysed in the light of gene evolutionary rates. Among those genes characterized by rapid evolution and expressed in a restricted range of tissues, SSD was represented in a larger proportion than WGD. Conversely, genes with ubiquitous expression were associated with slower evolutionary rates and a larger proportion of WGD. The results also show that evolutionary rates were faster in genes expressed in endodermal tissues and slower in ectodermal genes. Accordingly, the proportion of the SSD and WGD genes was highest in the endoderm and ectoderm, respectively. Therefore, quickly evolving SSD genes might have contributed to the faster evolution of endodermal tissues, whereas the comparatively slowly evolving WGD genes might have functioned to maintain the basic characteristics of ectodermal tissues. Mesenchymal tissues occupied an intermediate position in this regard, whereas the patterns observed for haemocytes were unique. Rapid tissue evolution could be related to a specific gene duplication mode (SSD) and faster molecular evolution in response to exposure to the external environment. These findings reveal general patterns underlying the evolution of tissues and their corresponding genes.

The odds of duplicate gene persistence after polyploidization

Background

Gene duplication is an important biological phenomenon associated with genomic redundancy, degeneration, specialization, innovation, and speciation. After duplication, both copies continue functioning when natural selection favors duplicated protein function or expression, or when mutations make them functionally distinct before one copy is silenced.

Results

Here we quantify the degree to which genetic parameters related to gene expression, molecular evolution, and gene structure in a diploid frog - Silurana tropicalis - influence the odds of functional persistence of orthologous duplicate genes in a closely related tetraploid species - Xenopus laevis. Using public databases and 454 pyrosequencing, we obtained genetic and expression data from S. tropicalis orthologs of 3,387 X. laevis paralogs and 4,746 X. laevis singletons - the most comprehensive dataset for African clawed frogs yet analyzed. Using logistic regression, we demonstrate that the most important predictors of the odds of duplicate gene persistence in the tetraploid species are the total gene expression level and evenness of expression across tissues and development in the diploid species. Slow protein evolution and information density (fewer exons, shorter introns) in the diploid are also positively correlated with duplicate gene persistence in the tetraploid.

Conclusions

Our findings suggest that a combination of factors contribute to duplicate gene persistence following whole genome duplication, but that the total expression level and evenness of expression across tissues and through development before duplication are most important. We speculate that these parameters are useful predictors of duplicate gene longevity after whole genome duplication in other taxa.

Comparative network analysis reveals that tissue specificity and gene function are important factors influencing the mode of expression evolution in Arabidopsis and rice

Microarray experiments have yielded massive amounts of expression information measured under various conditions for the model species Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). Expression compendia grouping multiple experiments make it possible to define correlated gene expression patterns within one species and to study how expression has evolved between species. We developed a robust framework to measure expression context conservation (ECC) and found, by analyzing 4,630 pairs of orthologous Arabidopsis and rice genes, that 77% showed conserved coexpression. Examples of nonconserved ECC categories suggested a link between regulatory evolution and environmental adaptations and included genes involved in signal transduction, response to different abiotic stresses, and hormone stimuli. To identify genomic features that influence expression evolution, we analyzed the relationship between ECC, tissue specificity, and protein evolution. Tissue-specific genes showed higher expression conservation compared with broadly expressed genes but were fast evolving at the protein level. No significant correlation was found between protein and expression evolution, implying that both modes of gene evolution are not strongly coupled in plants. By integration of cis-regulatory elements, many ECC conserved genes were significantly enriched for shared DNA motifs, hinting at the conservation of ancestral regulatory interactions in both model species. Surprisingly, for several tissue-specific genes, patterns of concerted network evolution were observed, unveiling conserved coexpression in the absence of conservation of tissue specificity. These findings demonstrate that orthologs inferred through sequence similarity in many cases do not share similar biological functions and highlight the importance of incorporating expression information when comparing genes across species.

Developmental decoupling of alternative phenotypes: insights from the transcriptomes of horn-polyphenic beetles

Developmental mechanisms play an important role in determining the costs, limits, and evolutionary consequences of phenotypic plasticity. One issue central to these claims is the hypothesis of developmental decoupling, where alternate morphs result from evolutionarily independent developmental pathways. We address this assumption through a microarray study that tests whether differences in gene expression between alternate morphs are as divergent as those between sexes, a classic example of developmental decoupling. We then examine whether genes with morph-biased expression are less conserved than genes with shared expression between morphs, as predicted if developmental decoupling relaxes pleiotropic constraints on divergence. We focus on the developing horns and brains of two species of horned beetles with impressive sexual- and morph-dimorphism in the expression of horns and fighting behavior. We find that patterns of gene expression were as divergent between morphs as they were between sexes. However, overall patterns of gene expression were also highly correlated across morphs and sexes. Morph-biased genes were more evolutionarily divergent, suggesting a role of relaxed pleiotropic constraints or relaxed selection. Together these results suggest that alternate morphs are to some extent developmentally decoupled, and that this decoupling has significant evolutionary consequences. However, alternative morphs may not be as developmentally decoupled as sometimes assumed and such hypotheses of development should be revisited and refined.

Gene discovery in the horned beetle Onthophagus taurus

Background

Horned beetles, in particular in the genus Onthophagus, are important models for studies on sexual selection, biological radiations, the origin of novel traits, developmental plasticity, biocontrol, conservation, and forensic biology. Despite their growing prominence as models for studying both basic and applied questions in biology, little genomic or transcriptomic data are available for this genus. We used massively parallel pyrosequencing (Roche 454-FLX platform) to produce a comprehensive EST dataset for the horned beetle Onthophagus taurus. To maximize sequence diversity, we pooled RNA extracted from a normalized library encompassing diverse developmental stages and both sexes.

Results

We used 454 pyrosequencing to sequence ESTs from all post-embryonic stages of O. taurus. Approximately 1.36 million reads assembled into 50,080 non-redundant sequences encompassing a total of 26.5 Mbp. The non-redundant sequences match over half of the genes in Tribolium castaneum, the most closely related species with a sequenced genome. Analyses of Gene Ontology annotations and biochemical pathways indicate that the O. taurus sequences reflect a wide and representative sampling of biological functions and biochemical processes. An analysis of sequence polymorphisms revealed that SNP frequency was negatively related to overall expression level and the number of tissue types in which a given gene is expressed. The most variable genes were enriched for a limited number of GO annotations whereas the least variable genes were enriched for a wide range of GO terms directly related to fitness.

Conclusions

This study provides the first large-scale EST database for horned beetles, a much-needed resource for advancing the study of these organisms. Furthermore, we identified instances of gene duplications and alternative splicing, useful for future study of gene regulation, and a large number of SNP markers that could be used in population-genetic studies of O. taurus and possibly other horned beetles.

Pathogenic peptide deviations support a model of adaptive evolution of chordate cardiac performance by troponin mutations

In cardiac muscle, the troponin (cTn) complex is a key regulator of myofilament calcium sensitivity because it serves as a molecular switch required for translating myocyte calcium fluxes into sarcomeric contraction and relaxation. Studies of several species suggest that ectotherm chordates have myofilaments with heightened calcium responsiveness. However, genetic polymorphisms in cTn that cause increased myofilament sensitivity to activating calcium in mammals result in cardiac disease including arrhythmias, diastolic dysfunction, and increased susceptibility to sudden cardiac death. We hypothesized that specific residue modifications in the regulatory arm of troponin I (TnI) were critical in mediating the observed decrease in myofilament calcium sensitivity within the mammalian taxa. We performed large-scale phylogenetic analysis, atomic resolution molecular dynamics simulations and modeling, and computational alanine scanning. This study provides evidence that a His to Ala substitution within mammalian cardiac TnI (cTnI) reduced the thermodynamic potential at the interface between cTnI and cardiac TnC (cTnC) in the calcium-saturated state by disrupting a strong intermolecular electrostatic interaction. This key residue modification reduced myofilament calcium sensitivity by making cTnI molecularly untethered from cTnC. To meet the requirements for refined mammalian adult cardiac performance, we propose that compensatory evolutionary pressures favored mutations that enhanced the relaxation properties of cTn by decreasing its sensitivity to activating calcium.

Elevated cardiac troponin I in a dog with an intracranial meningioma and evidence of myocardial necrosis

A 10-year-old, spayed female Weimaraner was referred for an acute onset of generalized tremors, seizures, and obtundation. Neurological examination revealed severe obtundation and a right-sided menace response deficit. Neuroanatomical diagnosis was consistent with a left prosencephalic lesion. The serum cardiac troponin I level was high, indicative of acute myocardial necrosis. With magnetic resonance imaging, a mass was observed in the left olfactory bulb and tract, with extensive edema in the white matter of the left cerebrum. The hippocampus was hyperintense on T2-weighted and T2-weighted fluid-attenuated inversion recovery images. At necropsy, a meningioma of the left olfactory bulb and ischemic cell change in the neurons of the hippocampus were identified. In the heart, microscopic lesions consistent with myocardial necrosis were observed. This is the first case to document an elevated cardiac troponin I level in a dog with intracranial disease and myocardial necrosis.

The molecular structures and expression patterns of zebrafish troponin I genes

Troponin I (TnnI), a constituent of the troponin complex on the thin filament, providers a calcium-sensitive switch for striated muscle contraction. Cardiac TnnI is, therefore, a highly sensitive and specific marker of myocardial injury in acute coronary syndromes. The TnnI gene, which has been identified in birds and mammals , encodes the isoforms expressed in cardiac muscle fast skeletal muscle and slow skeletal muscle. However, very little is known about the TnnI gene in lower vertebrates. Here, we cloned and characterized the molecular structures and expression patterns of three types of zebrafish tnni genes: tnni1, tnni2, and tnn-HC (Heart and craniofacial). Based on the unrooted radial gene tree analysis of the TnnI gene among vertebrates, the zebrafish Tnni1 and TnnI2 we cloned were homologous of the slow muscle TnnI1 and fast muscle TnnI2 of other vertebrates, respectively. In addition, reverse transcription-polymerase chain reaction (RT-PCR) and whole-mount in situ hybridization demonstrated that zebrafish tnni1 and tnni2 transcripts were not detectable in the somites until 16 h post-fertilization (hpf), after which they were identified as slow-and fast muscle-specific, respectively . Interestingly, tnni-HC, a novel tnni isoform of zebrafish was expressed exclusive in heart during early cardiogenesis as 16 hpf, but then extended its expression in craniofacial muscle after 48 hpf. Thus, using zebrafish as our system model, it is suggested that the results, as noted above, may provide more insight into the molecular structure and expression pattens of the lower vertebrate TnnI gene.

Multifunctionality dominantly determines the rate of human housekeeping and tissue specific interacting protein evolution

Elucidation of the determinants of the rate of protein sequence evolution is one of the great challenges in evolutionary biology. It has been proposed that housekeeping genes are evolutionarily slower than tissue specific genes. In the present communication, we have examined different determinants that influence the evolutionary rate variation in human housekeeping and tissue specific proteins present in protein-protein interaction network. Studies on yeast proteome, revealed a predominant role of protein connectivity in determining the rate of protein evolution. However, in human, we did not observe any significant influence of protein connectivity on its evolutionary rate. Rather, a significant impact of the proportion of protein's interacting length (amount of protein interface involved in interaction with its partners), expression level and multifunctionality has been observed in determining the rate of protein evolution. We also observed that multi interface proteins are evolutionarily conserved between housekeeping and tissue specific genes and it has been found that the average number of biological processes they associated in these two sets of genes is similar. Moreover, single interface proteins in housekeeping genes evolve more slowly as compared to tissue specific genes owing to their involvement in different number of biological processes. Partial correlation analysis suggests that the relative importance of three individual factors in determining the evolutionary rate variation between housekeeping and tissue specific proteins is in the order of protein multifunctionality>protein expression level>interacting protein length.

Differential selective constraints shaping codon usage pattern of housekeeping and tissue-specific homologous genes of rice and arabidopsis

Intra-genomic variation between housekeeping and tissue-specific genes has always been a study of interest in higher eukaryotes. To-date, however, no such investigation has been done in plants. Availability of whole genome expression data for both rice and Arabidopsis has made it possible to examine the evolutionary forces in shaping codon usage pattern in both housekeeping and tissue-specific genes in plants. In the present work, we have taken 4065 rice-Arabidopsis homologous gene pairs to study evolutionary forces responsible for codon usage divergence between housekeeping and tissue-specific genes. In both rice and Arabidopsis, it is mutational bias that regulates error minimization in highly expressed genes of both housekeeping and tissue-specific genes. Our results show that, in comparison to tissue-specific genes, housekeeping genes are under strong selective constraint in plants. However, in tissue-specific genes, lowly expressed genes are under stronger selective constraint compared with highly expressed genes. We demonstrated that constraint acting on mRNA secondary structure is responsible for modulating codon usage variations in rice tissue-specific genes. Thus, different evolutionary forces must underline the evolution of synonymous codon usage of highly expressed genes of housekeeping and tissue-specific genes in rice and Arabidopsis.

Patterns of positive selection in six Mammalian genomes

Genome-wide scans for positively selected genes (PSGs) in mammals have provided insight into the dynamics of genome evolution, the genetic basis of differences between species, and the functions of individual genes. However, previous scans have been limited in power and accuracy owing to small numbers of available genomes. Here we present the most comprehensive examination of mammalian PSGs to date, using the six high-coverage genome assemblies now available for eutherian mammals. The increased phylogenetic depth of this dataset results in substantially improved statistical power, and permits several new lineage- and clade-specific tests to be applied. Of approximately 16,500 human genes with high-confidence orthologs in at least two other species, 400 genes showed significant evidence of positive selection (FDR<0.05), according to a standard likelihood ratio test. An additional 144 genes showed evidence of positive selection on particular lineages or clades. As in previous studies, the identified PSGs were enriched for roles in defense/immunity, chemosensory perception, and reproduction, but enrichments were also evident for more specific functions, such as complement-mediated immunity and taste perception. Several pathways were strongly enriched for PSGs, suggesting possible co-evolution of interacting genes. A novel Bayesian analysis of the possible "selection histories" of each gene indicated that most PSGs have switched multiple times between positive selection and nonselection, suggesting that positive selection is often episodic. A detailed analysis of Affymetrix exon array data indicated that PSGs are expressed at significantly lower levels, and in a more tissue-specific manner, than non-PSGs. Genes that are specifically expressed in the spleen, testes, liver, and breast are significantly enriched for PSGs, but no evidence was found for an enrichment for PSGs among brain-specific genes. This study provides additional evidence for widespread positive selection in mammalian evolution and new genome-wide insights into the functional implications of positive selection.

Do Alu repeats drive the evolution of the primate transcriptome?

BACKGROUND

Of all repetitive elements in the human genome, Alus are unusual in being enriched near to genes that are expressed across a broad range of tissues. This has led to the proposal that Alus might be modifying the expression breadth of neighboring genes, possibly by providing CpG islands, modifying transcription factor binding, or altering chromatin structure. Here we consider whether Alus have increased expression breadth of genes in their vicinity.

RESULTS

Contrary to the modification hypothesis, we find that those genes that have always had broad expression are richest in Alus, whereas those that are more likely to have become more broadly expressed have lower enrichment. This finding is consistent with a model in which Alus accumulate near broadly expressed genes but do not affect their expression breadth. Furthermore, this model is consistent with the finding that expression breadth of mouse genes predicts Alu density near their human orthologs. However, Alus were found to be related to some alternative measures of transcription profile divergence, although evidence is contradictory as to whether Alus associate with lowly or highly diverged genes. If Alu have any effect it is not by provision of CpG islands, because they are especially rare near to transcriptional start sites. Previously reported Alu enrichment for genes serving certain cellular functions, suggested to be evidence of functional importance of Alus, appears to be partly a byproduct of the association with broadly expressed genes.

CONCLUSION

The abundance of Alu near broadly expressed genes is better explained by their preferential preservation near to housekeeping genes rather than by a modifying effect on expression of genes.

Rate of evolution in brain-expressed genes in humans and other primates

Brain-expressed genes are known to evolve slowly in mammals. Nevertheless, since brains of higher primates have evolved rapidly, one might expect acceleration in DNA sequence evolution in their brain-expressed genes. In this study, we carried out full-length cDNA sequencing on the brain transcriptome of an Old World monkey (OWM) and then conducted three-way comparisons among (i) mouse, OWM, and human, and (ii) OWM, chimpanzee, and human. Although brain-expressed genes indeed appear to evolve more rapidly in species with more advanced brains (apes > OWM > mouse), a similar lineage effect is observable for most other genes. The broad inclusion of genes in the reference set to represent the genomic average is therefore critical to this type of analysis. Calibrated against the genomic average, the rate of evolution among brain-expressed genes is probably lower (or at most equal) in humans than in chimpanzee and OWM. Interestingly, the trend of slow evolution in coding sequence is no less pronounced among brain-specific genes, vis-à-vis brain-expressed genes in general. The human brain may thus differ from those of our close relatives in two opposite directions: (i) faster evolution in gene expression, and (ii) a likely slowdown in the evolution of protein sequences. Possible explanations and hypotheses are discussed.

When calibrated against the genomic average, the rate of evolution among brain-expressed genes in humans is probably lower than or equal to that of other closely related primates.

Whether comparing morphology or cognitive ability, it is clear that the human brain has evolved rapidly relative to that of other primates. But the extent to which genes expressed in the brain also reflect this overall pattern is unclear. To address this question, it's necessary to measure any variations in the DNA sequences of these genes between human and chimpanzee. And, to do this as accurately as possible, it's also important to require an appropriate reference group to act as a benchmark against which the differences can be measured. We therefore compared publicly available genomic sequences of chimps and humans with complementary DNA sequences of several thousand genes expressed in the brain of another closely related primate—the macaque, an Old World monkey—as well as the more distantly related mouse. Our analyses of the rates of protein evolution in these species suggest that genes expressed in the human brain have in fact slowed down in their evolution since the split between human and chimpanzee, contrary to some previously published reports. We suggest that advanced brains are driven primarily by the increasing complexity in the network of gene interactions. As a result, brain-expressed genes are constrained in their sequence evolution, although their expression levels may change rapidly.

Gender-specific selection on codon usage in plant genomes

Background

Currently, there is little data available regarding the role of gender-specific gene expression on synonymous codon usage (translational selection) in most organisms, and particularly plants. Using gender-specific EST libraries (with > 4000 ESTs) from Zea mays and Triticum aestivum, we assessed whether gender-specific gene expression per se and gender-specific gene expression level are associated with selection on codon usage.

Results

We found clear evidence of a greater bias in codon usage for genes expressed in female than in male organs and gametes, based on the variation in GC content at third codon positions and the frequency of species-preferred codons. This finding holds true for both highly and for lowly expressed genes. In addition, we found that highly expressed genes have greater codon bias than lowly expressed genes for both female- and male-specific genes. Moreover, in both species, genes with female-specific expression show a greater usage of species-specific preferred codons for each of the 18 amino acids having synonymous codons. A supplemental analysis of Brassica napus suggests that bias in codon usage could also be higher in genes expressed in male gametophytic tissues than in heterogeneous (flower) tissues.

Conclusion

This study reports gender-specific bias in codon usage in plants. The findings reported here, based on the analysis of 1 497 876 codons, are not caused either by differences in the biological functions of the genes or by differences in protein lengths, nor are they likely attributable to mutational bias. The data are best explained by gender-specific translational selection. Plausible explanations for these findings and the relevance to these and other organisms are discussed.

Inference of expression-dependent negative selection based on polymorphism and divergence in the human genome

There is a mounting evidence for the correlation between the gene expression pattern and sequence divergence. However, little is known about the relationship between the gene expression pattern and polymorphism. We compiled the gene expression, polymorphism, and divergence data from the public databases of the human genome. The ratios of nonsynonymous (A) to synonymous (S) substitutions in polymorphism and divergence in the human genome were strongly influenced by the expression pattern and breadth of genes and showed strong correlations. Among the tissues we analyzed, the brain-expressed genes have the smallest and the liver-expressed genes have the largest proportion of amino acid changes both in polymorphism and divergence. The analysis implies that negative selection is the primary factor affecting expression-dependent gene evolution and the prevalent but nonuniform distribution of slightly deleterious mutations in the genome. Although the genes under relaxed negative selection evolved faster than the other genes, these genes are even more liable to slightly deleterious mutations in the population. On the other hand, nonneutral mutations in the highly conservative genes, such as brain-expressed and housekeeping genes, are largely deleterious and eliminated before they enter the population.

The two AGPase subunits evolve at different rates in angiosperms, yet they are equally sensitive to activity-altering amino acid changes when expressed in bacteria

The rate of protein evolution is generally thought to reflect, at least in part, the proportion of amino acids within the protein that are needed for proper function. In the case of ADP-glucose pyrophosphorylase (AGPase), this premise led to the hypothesis that, because the AGPase small subunit is more conserved compared with the large subunit, a higher proportion of the amino acids of the small subunit are required for enzyme activity compared with the large subunit. Evolutionary analysis indicates that the AGPase small subunit has been subject to more intense purifying selection than the large subunit in the angiosperms. However, random mutagenesis and expression of the maize (Zea mays) endosperm AGPase in bacteria show that the two AGPase subunits are equally predisposed to enzyme activity-altering amino acid changes when expressed in one environment with a single complementary subunit. As an alternative hypothesis, we suggest that the small subunit exhibits more evolutionary constraints in planta than does the large subunit because it is less tissue specific and thus must form functional enzyme complexes with different large subunits. Independent approaches provide data consistent with this alternative hypothesis.

The use of evolutionary biology concepts for genome annotation

The past decade has seen the completion of numerous whole-genome sequencing projects, began with bacterial genomes and continued with eukaryotic species from different phyla: fungi, plants and animals. Besides, more biological information are produced and are shared thanks to information exchange systems, and more biological concepts, as well as more bioinformatics tools, are available. In this article, we will describe how the evolutionary biology concepts, as well as computer science, are useful for a better understanding of biology in general and genome annotation in particular. The genome annotation process consists of taking the raw DNA produced, for example, by the genome sequencing projects, adding the layers of analysis and interpretation necessary to extract its biological significance and placing it in the context of our understanding of biological processes. Genome annotation is a multistep process falling into two broad categories: structural and functional annotation.

Relating tissue specialization to the differentiation of expression of singleton and duplicate mouse proteins

An analysis of the relationship between duplication events, the time they took place and the expression breadth of the duplicated genes supports the subfunctionalization model, in which expression divergence following gene duplication promotes the retention of a gene in multicellular species.

Background

Gene duplications have been hypothesized to be a major factor in enabling the evolution of tissue differentiation. Analyses of the expression profiles of duplicate genes in mammalian tissues have indicated that, with time, the expression patterns of duplicate genes diverge and become more tissue specific. We explored the relationship between duplication events, the time at which they took place, and both the expression breadth of the duplicated genes and the cumulative expression breadth of the gene family to which they belong.

Results

We show that only duplicates that arose through post-multicellularity duplication events show a tendency to become more specifically expressed, whereas such a tendency is not observed for duplicates that arose in a unicellular ancestor. Unlike the narrow expression profile of the duplicated genes, the overall expression of gene families tends to maintain a global expression pattern.

Conclusion

The work presented here supports the view suggested by the subfunctionalization model, namely that expression divergence in different tissues, following gene duplication, promotes the retention of a gene in the genome of multicellular species. The global expression profile of the gene families suggests division of expression between family members, whose expression becomes specialized. Because specialization of expression is coupled with an increased rate of sequence divergence, it can facilitate the evolution of new, tissue-specific functions.

The monosaccharide transporter gene family in land plants is ancient and shows differential subfamily expression and expansion across lineages

Background

In plants, tandem, segmental and whole-genome duplications are prevalent, resulting in large numbers of duplicate loci. Recent studies suggest that duplicate genes diverge predominantly through the partitioning of expression and that breadth of gene expression is related to the rate of gene duplication and protein sequence evolution.

Here, we utilize expressed sequence tag (EST) data to study gene duplication and expression patterns in the monosaccharide transporter (MST) gene family across the land plants. In Arabidopsis, there are 53 MST genes that form seven distinct subfamilies. We created profile hidden Markov models of each subfamily and searched EST databases representing diverse land plant lineages to address the following questions: 1) Are homologs of each Arabidopsis subfamily present in the earliest land plants? 2) Do expression patterns among subfamilies and individual genes within subfamilies differ across lineages? 3) Has gene duplication within each lineage resulted in lineage-specific expansion patterns? We also looked for correlations between relative EST database representation in Arabidopsis and similarity to orthologs in early lineages.

Results

Homologs of all seven MST subfamilies were present in land plants at least 400 million years ago. Subfamily expression levels vary across lineages with greater relative expression of the STP, ERD6-like, INT and PLT subfamilies in the vascular plants. In the large EST databases of the moss, gymnosperm, monocot and eudicot lineages, EST contig construction reveals that MST subfamilies have experienced lineage-specific expansions. Large subfamily expansions appear to be due to multiple gene duplications arising from single ancestral genes. In Arabidopsis, one or a few genes within most subfamilies have much higher EST database representation than others. Most highly represented (broadly expressed) genes in Arabidopsis have best match orthologs in early divergent lineages.

Conclusion

The seven subfamilies of the Arabidopsis MST gene family are ancient in land plants and show differential subfamily expression and lineage-specific subfamily expansions. Patterns of gene expression in Arabidopsis and correlation of highly represented genes with best match homologs in early lineages suggests that broadly expressed genes are often highly conserved, and that most genes have more limited expression.

Impacts of Gene Essentiality, Expression Pattern, and Gene Compactness on the Evolutionary Rate of Mammalian Proteins

Understanding the determinants of the rate of protein sequence evolution is of fundamental importance in evolutionary biology. Many recent studies have focused on the yeast because of the availability of many genome-wide expressional and functional data. Yeast studies revealed a predominant role of gene expression level and a minor role of gene essentiality in determining the rate of protein sequence evolution. Whether these rules apply to complex organisms such as mammals is unclear. Here we assemble a list of 1,138 essential and 2,341 nonessential mouse genes based on targeted gene deletion experiments and report a significant impact of gene essentiality on the rate of mammalian protein evolution. Gene expression level has virtually no effect, although tissue specificity in expression pattern has a strong influence. Unexpectedly, gene compactness, measured by average intron size and untranslated region length, has the greatest influence. Hence, the relative importance of the various factors in determining the rate of mammalian protein evolution is gene compactness > gene essentiality approximately tissue specificity > expression level. Our results suggest a considerable variation in rate determinants between unicellular organisms such as the yeast and multicellular organisms such as mammals.

Accelerated evolutionary rate may be responsible for the emergence of lineage-specific genes in ascomycota

The evolutionary origin of "orphan" genes, genes that lack sequence similarity to any known gene, remains a mystery. One suggestion has been that most orphan genes evolve rapidly so that similarity to other genes cannot be traced after a certain evolutionary distance. This can be tested by examining the divergence rates of genes with different degrees of lineage specificity. Here the lineage specificity (LS) of a gene describes the phylogenetic distribution of that gene's orthologues in related species. Highly lineage-specific genes will be distributed in fewer species in a phylogeny. In this study, we have used the complete genomes of seven ascomycotan fungi and two animals to define several levels of LS, such as Eukaryotes-core, Ascomycota-core, Euascomycetes-specific, Hemiascomycetes-specific, Aspergillus-specific, and Saccharomyces-specific. We compare the rates of gene evolution in groups of higher LS to those in groups with lower LS. Molecular evolutionary analyses indicate an increase in nonsynonymous nucleotide substitution rates in genes with higher LS. Several analyses suggest that LS is correlated with the evolutionary rate of the gene. This correlation is stronger than those of a number of other factors that have been proposed as predictors of a gene's evolutionary rate, including the expression level of genes, gene essentiality or dispensability, and the number of protein-protein interactions. The accelerated evolutionary rates of genes with higher LS may reflect the influence of selection and adaptive divergence during the emergence of orphan genes. These analyses suggest that accelerated rates of gene evolution may be responsible for the emergence of apparently orphan genes.

Low rates of expression profile divergence in highly expressed genes and tissue-specific genes during mammalian evolution

Evolutionary rates provide important information about the pattern and mechanism of evolution. Although the rate of gene sequence evolution has been well studied, the rate of gene expression evolution is poorly understood. In particular, it is unclear whether the gene expression level and tissue specificity influence the divergence of expression profiles between orthologous genes. Here we address this question using a microarray data set comprising the expression signals of 10,607 pairs of orthologous human and mouse genes from over 60 tissues per species. We show that the level of gene expression and the degree of tissue specificity are generally conserved between the human and mouse orthologs. The rate of gene expression profile change during evolution is negatively correlated with the level of gene expression, measured by either the average or the highest level among all tissues examined. This is analogous to the observation that the rate of gene (or protein) sequence evolution is negatively correlated with the gene expression level. The impacts of the degree of tissue specificity on the evolutionary rate of gene sequence and that of expression profile, however, are opposite. Highly tissue-specific genes tend to evolve rapidly at the gene sequence level but slowly at the expression profile level. Thus, different forces and selective constraints must underlie the evolution of gene sequence and that of gene expression.

Evolutionary rates and expression level in Chlamydomonas

In many biological systems, especially bacteria and unicellular eukaryotes, rates of synonymous and nonsynonymous nucleotide divergence are negatively correlated with the level of gene expression, a phenomenon that has been attributed to natural selection. Surprisingly, this relationship has not been examined in many important groups, including the unicellular model organism Chlamydomonas reinhardtii. Prior to this study, comparative data on protein-coding sequences from C. reinhardtii and its close noninterfertile relative C. incerta were very limited. We compiled and analyzed protein-coding sequences for 67 nuclear genes from these taxa; the sequences were mostly obtained from the C. reinhardtii EST database and our C. incerta EST data. Compositional and synonymous codon usage biases varied among genes within each species but were highly correlated between the orthologous genes of the two species. Relative rates of synonymous and nonsynonymous substitution across genes varied widely and showed a strong negative correlation with the level of gene expression estimated by the codon adaptation index. Our comparative analysis of substitution rates in introns of lowly and highly expressed genes suggests that natural selection has a larger contribution than mutation to the observed correlation between evolutionary rates and gene expression level in Chlamydomonas.

Gene Expression Evolves Faster in Narrowly Than in Broadly Expressed Mammalian Genes

Despite much recent interest, it remains unclear what determines the rate of evolution of gene expression. To study this issue we develop a new measure, called "Expression Conservation Index" (ECI), to quantify the degree of tissue-expression conservation between two homologous genes. Applying this measure to a large set of gene expression data from human and mouse, we show that tissue expression tends to evolve rapidly for genes that are expressed in only a limited number of tissues, whereas tissue expression can be conserved for a long time for genes expressed in a large number of tissues. Therefore, expression breadth is an important determinant for evolutionary conservation of tissue expression. In addition, we find a rapid decrease in ECI with the synonymous divergence between duplicate genes, suggesting fast divergence in tissue expression between duplicate genes.

Relationship between the tissue-specificity of mouse gene expression and the evolutionary origin and function of the proteins

A microaaray analysis of mouse gene expression combined with the proteins functional and phyletic classification suggests that phyletic age (and not function) is the dominant factor shaping the expression profle of a protein.

Background

The combination of complete genome sequence information with expression data enables us to characterize the relationship between a protein's evolutionary origin or functional category and its expression pattern. In this study, mouse proteins were assigned into functional and phyletic groups and the gene expression patterns of the different protein groupings were examined by microarray analysis in various mouse tissues.

Results

Our results suggest that the proteins that are universally distributed in all tissues are predominantly enzymes and transporters. In contrast, the tissue-specific set is dominated by regulatory proteins (signal transduction and transcription factors). An increased tendency to tissue-specificity is observed for metazoan-specific proteins. As the composition of the phyletic groups highly correlates with that of the functional groups, the data were tested in order to determine which of the two factors - function or phyletic age - is dominant in shaping the expression profile of a protein. The observed differences in expression patterns of genes between functional groups were found mainly to reflect their different phyletic origin. The connection between tissue specificity and phyletic age cannot be explained by the recent rate of evolution. Finally, although metazoan-specific proteins tend to be tissue-specific compared with phyletically conserved proteins present in all domains of life, many such 'universal' proteins are also tissue-specific.

Conclusion

The minimal cellular transcriptome of the metazoan cell differs from that of the ancestral unicellular eukaryote: new functions were added (metazoan-specific proteins), whilst other functions became specialized and no longer took place in all cells (tissue-specific pre-metazoan proteins).

Cardiac troponin I: evaluation of a biomarker for the diagnosis of heart disease in the dog

OBJECTIVES

To assess the value of measuring blood levels of the myocardial protein cardiac troponin I (cTnl) in the diagnosis of congenital and acquired heart disease in the dog and in the evaluation of the severity of heart failure.

METHODS

Serum samples obtained from healthy dogs (n = 26) and from dogs diagnosed with a variety of congenital and acquired heart conditions (n = 35) were assayed for cTnl concentration using an automated immunoassay method. Results were also analysed according to the degree of heart failure as assessed using the International Small Animal Cardiac Health Council's scheme.

RESULTS

Healthy dogs had very low or undetectable blood cTnl levels, as did dogs with congenital heart disease. However, cTnl levels were significantly elevated in dogs with acquired mitral valve disease, dilated cardiomyopathy and pericardial effusion. Blood cTnl levels also varied with severity of heart failure.

CLINICAL SIGNIFICANCE

Measurement of blood cTnl levels may be a useful aid in the diagnosis of dogs with suspected heart disease and in indicating the severity of heart failure.

Diversification and Independent Evolution of Troponin C Genes in Insects

Troponin C (TpnC), the calcium-binding subunit of the troponin regulatory complex in the muscle thin filament, is encoded by multiple genes in insects. To understand how TpnC genes have evolved, we characterized the gene number and structure in a number of insect species. The TpnC gene complement is five genes in Drosophilidae as previously reported for D. melanogaster. Gene structures are almost identical in D. pseudoobscura, D. suboboscura, and D. virilis. Developmental patterns of expression are also conserved in Drosophila subobscura and D. virilis. Similar, but not completely equivalent, TpnC gene repertoires have been identified in the Anopheles gambiae and Apis mellifera genomes. Insect TpnC sequences can be divided into three groups, allowing a systematic classification of newly identified genes. The pattern of expression of the Apis mellifera genes essentially agrees with the pattern in Drosophilidae, providing further functional support to the classification. A model for the evolution of the TpnC genes is proposed including the most likely pathway of insect TpnC diversification. Our results suggest that the rapid increase in number and sequence specialization of the adult Type III isoforms can be correlated with the evolution of the holometabolous mode of development and the acquisition of asynchronous indirect flight muscle function in insects. This evolutionarily specialization has probably been achieved independently in different insect orders.

Amphibian cardiac troponin I gene's organization, developmental expression, and regulatory properties are different from its mammalian homologue

In mammals, the expression of the troponin I-slow (TnIs) isoform is predominant in the heart during embryogenesis and, shortly after birth, is replaced by the cardiac-specific isoform, TnIc; a developmental switch thought to be mediated by thyroid hormone. Whereas, in Xenopus, TnIc is expressed at the onset of heart formation and is the only TnI isoform expressed in the heart. Herein, we demonstrate that the expression patterns of these genes appear to be common within the anuran lineage and, unlike their mammalian counterparts, are not affected by thyroid hormone. To elucidate the regulatory mechanism(s) governing the expression of the amphibian TnIc gene, we characterized the TnIc gene from Rana catesbeiana and used its 5'-flanking region to drive expression of green fluorescent protein in the Xenopus transgenic system. Our results demonstrate that a 300-bp minimal promoter containing intact GATA and CArG-box elements is sufficient to drive expression of this reporter gene in a pattern that mimics, both spatially and temporally, the expression of the endogenous Xenopus TnIc gene.

Genomic regionality in rates of evolution is not explained by clustering of genes of comparable expression profile

In mammalian genomes, linked genes show similar rates of evolution, both at fourfold degenerate synonymous sites (K4) and at nonsynonymous sites (KA). Although it has been suggested that the local similarity in the synonymous substitution rate is an artifact caused by the inclusion of disparately evolving gene pairs, we demonstrate here that this is not the case: after removal of disparately evolving genes, both (1) linked genes and (2) introns from the same gene have more similar silent substitution rates than expected by chance. What causes the local similarity in both synonymous and nonsynonymous substitution rates? One class of hypotheses argues that both may be related to the observed clustering of genes of comparable expression profile. We investigate these hypotheses using substitution rates from both human-mouse and mouse-rat comparisons, and employing three different methods to assay expression parameters. Although we confirm a negative correlation of expression breadth with both K4 and KA, we find no evidence that clustering of similarly expressed genes explains the clustering of genes of comparable substitution rates. If gene expression is not responsible, what about other causes? At least in the human-mouse comparison, the local similarity in KA can be explained by the covariation of KA and K4. As regards K4, our results appear consistent with the notion that local similarity is due to processes associated with meiotic recombination.

Effects of gene expression on molecular evolution in Arabidopsis thaliana and Arabidopsis lyrata

We analyzed the complete genome sequence of Arabidopsis thaliana and sequence data from 83 genes in the outcrossing A. lyrata, to better understand the role of gene expression on the strength of natural selection on synonymous and replacement sites in Arabidopsis. From data on tRNA gene abundance, we find a good concordance between codon preferences and the relative abundance of isoaccepting tRNAs in the complete A. thaliana genome, consistent with models of translational selection. Both EST-based and new quantitative measures of gene expression (MPSS) suggest that codon preferences derived from information on tRNA abundance are more strongly associated with gene expression than those obtained from multivariate analysis, which provides further support for the hypothesis that codon bias in Arabidopsis is under selection mediated by tRNA abundance. Consistent with previous results, analysis of protein evolution reveals a significant correlation between gene expression level and amino acid substitution rate. Analysis by MPSS estimates of gene expression suggests that this effect is primarily the result of a correlation between the number of tissues in which a gene is expressed and the rate of amino acid substitution, which indicates that the degree of tissue specialization may be an important determinant of the rate of protein evolution in Arabidopsis.

Mammalian Housekeeping Genes Evolve More Slowly than Tissue-Specific Genes

Do housekeeping genes, which are turned on most of the time in almost every tissue, evolve more slowly than genes that are turned on only at specific developmental times or tissues? Recent large-scale gene expression studies enable us to have a better definition of housekeeping genes and to address the above question in detail. In this study, we examined 1581 human-mouse orthologous gene pairs for their patterns of sequence evolution, contrasting housekeeping genes with tissue-specific genes. Our results show that, in comparison to tissue-specific genes, housekeeping genes on average evolve more slowly and are under stronger selective constraints as reflected by significantly smaller values of Ka/Ks. Besides stronger purifying selection, we explored several other factors that can possibly slow down nonsynonymous rates in housekeeping genes. Although mutational bias might slightly slow the nonsynonymous rates in housekeeping genes, it is unlikely to be the major cause of the rate difference between the two types of genes. The codon usage pattern of housekeeping genes does not seem to differ from that of tissue-specific genes. Moreover, contrary to the old textbook concept, we found that approximately 74% of the housekeeping genes in our study belong to multigene families, not significantly different from that of the tissue-specific genes ( approximately 70%). Therefore, the stronger selective constraints on housekeeping genes are not due to a lower degree of genetic redundancy.

Adaptive evolution of MRG, a neuron-specific gene family implicated in nociception

The MRG gene family (also known as SNSR) belongs to the G-protein-coupled receptor (GPCR) superfamily, is expressed specifically in nociceptive neurons, and is implicated in the modulation of nociception. Here, we show that Ka/Ks (the ratio between nonsynonymous and synonymous substitution rates) displays distinct profiles along the coding regions of MRG, with peaks (Ka/Ks>1) corresponding to extracellular domains, and valleys (Ka/Ks<1) corresponding to transmembrane and cytoplasmic domains. The extracellular domains are also characterized by a significant excess of radical amino acid changes. Statistical analysis shows that positive selection is by far the most suitable model to account for the nucleotide substitution patterns in MRG. Together, these results demonstrate that the extracellular domains of the MRG receptor family, which presumably partake in ligand binding, have experienced strong positive selection. Such selection is likely directed at altering the sensitivity and/or selectivity of nociceptive neurons to aversive stimuli. Thus, our finding suggests pain perception as an aspect of the nervous system that may have experienced a surprising level of adaptive evolution.

Translational selection and yeast proteome evolution

The primary structures of peptides may be adapted for efficient synthesis as well as proper function. Here, the Saccharomyces cerevisiae genome sequence, DNA microarray expression data, tRNA gene numbers, and functional categorizations of proteins are employed to determine whether the amino acid composition of peptides reflects natural selection to optimize the speed and accuracy of translation. Strong relationships between synonymous codon usage bias and estimates of transcript abundance suggest that DNA array data serve as adequate predictors of translation rates. Amino acid usage also shows striking relationships with expression levels. Stronger correlations between tRNA concentrations and amino acid abundances among highly expressed proteins than among less abundant proteins support adaptation of both tRNA abundances and amino acid usage to enhance the speed and accuracy of protein synthesis. Natural selection for efficient synthesis appears to also favor shorter proteins as a function of their expression levels. Comparisons restricted to proteins within functional classes are employed to control for differences in amino acid composition and protein size that reflect differences in the functional requirements of proteins expressed at different levels.

A third isoform of cytochrome c oxidase subunit VIII is present in mammals

The terminal enzyme of the mitochondrial respiratory chain, cytochrome c oxidase (COX), contains three mitochondrial and ten nuclear encoded subunits in mammals. Three of the nuclear subunits (VIa, VIIa, and VIII) have muscle and non-muscle-specific isoforms, subunit IV contains a lung-specific isoform, and subunit VIb contains a testes-specific isoform. For subunit VIII, the smallest nuclear encoded COX polypeptide, we have now found a third gene (COX 8-3), which has been identified in human, lemur, rat, and mouse, suggesting that it is present in a broad range of Eutherian mammals. Sequence similarity and gene structure support the homology of COX8-3 to the other subunit VIII isoforms, indicating that all three are the product of gene duplications. COX VIII-3 protein is mitochondrially-targeted, as shown by a fluorescent COX VIII3/DsRed fusion protein. Both the mitochondrial targeting and its sequence conservation suggest that COXVIII-3 functions as part of the COX holoenzyme and could have a tissue-specific role, as is the case for the other two isoforms. Questions remain about where COX8-3 is predominantly expressed. However, detection of full-length cDNAs, lower levels of sequence divergence at the first and second codon positions compared to the third, and a conserved gene structure indicate that COX VIII-3 is an expressed gene whose origin dates to at least 91 million years ago.

Troponin I binds polycystin-L and inhibits its calcium-induced channel activation

Polycystin-L (PCL) is an isoform of polycystin-2, the product of the second gene associated with autosomal dominant polycystic kidney disease, and functions as a Ca(2+)-regulated nonselective cation channel. We recently demonstrated that polycystin-2 interacts with troponin I, an important regulatory component of the actin microfilament complex in striated muscle cells and an angiogenesis inhibitor. In this study, using the two-microelectrode voltage-clamp technique and Xenopus oocyte expression system, we showed that the calcium-induced PCL channel activation is substantially inhibited by the skeletal and cardiac troponin I (60% and 31% reduction, respectively). Reciprocal co-immunoprecipitation experiments demonstrated that PCL physically associates with the skeletal and cardiac troponin I isoforms in overexpressed Xenopus oocytes and mouse fibroblast NIH 3T3 cells. Furthermore, both native PCL and cardiac troponin I were present in human heart tissues where they indeed associate with each other. GST pull-down and microtiter binding assays showed that the C-terminus of PCL interacts with the troponin I proteins. The yeast two-hybrid assay further verified this interaction and defined the corresponding interacting domains of the PCL C-terminus and troponin I. Taken together, this study suggests that troponin I acts as a regulatory subunit of the PCL channel complex and provides the first direct evidence that PCL is associated with the actin cytoskeleton through troponin I.

Expression of enhancers is altered in Drosophila melanogaster hybrids

The molecular foundations of evolution are difficult to trace because most protein sequences are virtually identical in closely related species. The largest fraction of sequence within the genome, however, is composed of noncoding sequences where regulatory elements locate to various sites. It has been suggested that changes in the activity of these elements may trigger evolutionary change. In Drosophila, the enhancer trap procedure identifies regulatory sequences in the genome after the insertion of a P-element-based construct. We generated new insertions and characterized their expression domains in the adult eye and larval imaginal disks using the white and LacZ reporter genes. Lines with robust expression patterns in D. melanogaster were analyzed in hybrids to test the conservation of regulatory mechanisms between species. Most of the enhancers used in this study modified their expression in hybrids with the mating species D. mauritiana and D. simulans. Expression changes resulted either in gain or loss of expression and were cell-type or hybrid-genome specific. Further characterization of a limited number of enhancers in D. melanogaster showed that expression domains could adapt to changes in cell number during development but not after the completion of cell proliferation. Also, expression of some enhancers appeared to be sensitive to heterochromatin from the Y but not the X chromosome. Taken together, these results demonstrate the high sensitivity of regulatory mechanisms of gene expression as a prime source of evolutionary change and suggest quantitative changes in available transcription factors as one of the mechanisms involved.

Polycystin-2 interacts with troponin I, an angiogenesis inhibitor

Polycystin-2 (PC2), encoded by the PKD2 gene, is mutated in 10-15% of autosomal dominant polycystic kidney disease (ADPKD) patients. PC2 is a Ca(2+)-permeable nonselective cation channel and is present in kidney and many other organs. Likewise, PKD2-mutated patients and mice exhibit extrarenal abnormalities. In comparison with cysts in the kidney, liver, and pancreas, abnormalities in the heart, brain, and vascular vessels are less understood. In particular, roles of PC2 in muscle and endothelia remain largely unknown. In the present study, using a yeast two-hybrid screening, we discovered that the PC2 carboxyl terminal domain (D682-V968) interacts with the cardiac troponin I, an important regulatory component of the actin microfilament in cardiac muscle cells. This interaction was demonstrated by GST pull-down and microtiter binding assays. Dose-dependent binding between PC2 and troponin I followed a Michaelis-Menten relationship, indicating a 1:1 binding stoichiometry. The interacting domains were located to the R872-H927 segment of PC2 and the M1-V107 and K106-L158 segments of troponin I. Co-immunoprecipitation experiments demonstrated that the cardiac and two skeletal isoforms of troponin I were all associated with PC2, when coexpressed in mouse fibroblast NIH 3T3 cells and Xenopus oocytes. Furthermore, reciprocal co-immunoprecipitation verified the interaction between the native polycystin-2 and troponin I in human adult heart tissues. This study thus provides new evidence for a direct attachment of PC2 to the actin microfilament network, in addition to the recently identified association between PC2 and trypomyosin-1. Troponin I functions as an inhibitory subunit of the troponin complex for calcium-dependent regulation of muscle contraction and as an inhibitor of angiogenesis seen in ADPKD. It is possible that altered interaction due to pathogenic polycystin-1 or -2 mutations can account for angiogenesis in ADPKD and may be corrected to some extent by exogenous troponin I.

The slow isoform of Xenopus troponin I is expressed in developing skeletal muscle but not in the heart

In birds and mammals three isoforms of troponin I (TnI) exist; a slow (TnIs), a fast (TnIf) and a cardiac (TnIc). Although each of these isoforms is expressed in the adult forms of these organisms in a muscle fiber-type-specific manner, the gene encoding TnIs is also expressed within the developing heart of these vertebrates. Herein, our results demonstrate that the developing heart of Xenopus laevis, unlike its counterpart in birds and mammals, does not express the gene encoding the TnIs isoform and that the expression of this gene, as well as the one encoding the Xenopus TnIf isoform, is restricted to skeletal muscle.

Identification of a matrix-binding domain in MAGP1 and MAGP2 and intracellular localization of alternative splice forms

MAGP1 is a small molecular mass protein associated with microfibrils in the extracellular matrix (ECM). To identify the molecular basis of its interaction with other microfibrillar proteins, deletion constructs of MAGP1 were expressed in a mammalian cell system that served as a model for microfibril assembly. This study identified a 54-amino acid sequence in the carboxyl-terminal region of the protein that defines a matrix-binding domain that is sufficient to target MAGP1 to the ECM. Site-directed mutagenesis demonstrated that binding activity is dependent on the presence of 7 cysteine residues in this sequence. MAGP2 contains a sequence similar to the matrix-binding domain of MAGP1, but could not associate with the ECM because of a single amino acid change. Two naturally occurring MAGP1 splice variants, MAGP1B (human-specific) and MAGP1D (found in mice), localized intracellularly when expressed as chimeric proteins with green fluorescent protein in rat lung fibroblasts. This suggests a second action site for MAGP1.