Genome and protein evolution in eukaryotes

Appropriate Nitrogen form Ratio and UV-A Supplementation Increased Quality and Production in Purple Lettuce (Lactuca sativa L.)

Purple lettuce (Lactuca sativa L. cv. Zhongshu Purple Lettuce) was chosen as the trial material, and LED intelligent light control consoles were used as the light sources. The purpose was to increase the yield and quality of purple lettuce while lowering its nitrate level. By adding various ratios of NO3--N and NH4+-N to the nutrient solution and 20 µmol m-2 s-1 UV-A based on white, red, and blue light (130, 120, 30 µmol m-2 s-1), the effects of different NO3--N/NH4+-N ratios (NO3--N, NO3--N/NH4+-N = 3/1, NH4+-N) and UV-A interaction on yield, quality, photosynthetic characteristics, anthocyanins, and nitrogen assimilation of purple lettuce were studied. In order to produce purple lettuce hydroponically under controlled environmental conditions, a theoretical foundation and technological specifications were developed, taking into account an appropriate UV-A dose and NO3--N/NH4+-N ratio. Results demonstrate that adding a 20 µmol m-2 s-1 UV-A, and a NO3--N/NH4+-N treatment of 3/1, significantly reduced the nitrate level while increasing the growth, photosynthetic rate, chlorophyll, carotenoid, and anthocyanin content of purple lettuce. The purple leaf lettuce leaves have an enhanced capacity to absorb nitrogen. Furthermore, plants have an acceleration of nitrogen metabolism, which raises the concentration of free amino acids and soluble proteins and promotes biomass synthesis. Thus, based on the NO3--N/NH4+-N (3/1) treatment, adding 20 µmol m-2 s-1 UV-A will be helpful in boosting purple lettuce production and decreasing its nitrate content.

The plant s1-like nuclease family has evolved a highly diverse range of catalytic capabilities

Plant S1-like nucleases, often referred to as nuclease I enzymes, are the main class of enzymes involved in nucleic acid degradation during plant programmed cell death. The catalytically active site of these enzymes shows a significant similarity to the well-described P1 nuclease from Penicillium citrinum. Previously published studies reported that plant S1-like nucleases possess catalytic activities similar to their fungal orthologs, i.e. they hydrolyze single-stranded DNA and RNA, and less efficiently double-stranded DNA, in the presence of zinc ions. Here we describe a comprehensive study of the nucleolytic activities of all Arabidopsis S1-like paralogs. Our results revealed that different members of this family are characterized by a surprisingly large variety of catalytic properties. We found that, in addition to Zn(2+)-dependent enzymes, this family also comprises nucleases activated by Ca(2+) and Mn(2+), which implies that the apparently well-known S1 nuclease active site in plant nucleases is able to cooperate with different activatory ions. Moreover, particular members of this class differ in their optimum pH value and substrate specificity. These results shed new light on the widely accepted classification of plant nucleases which is based on the assumption that the catalytic requirements of plant nucleases reflect their phylogenetic origin. Our results imply the need to redefine the understanding of the term 'nuclease I'. Analysis of the phylogenetic relationships between S1-like enzymes shows that plant representatives of this family evolve toward an increase in catalytic diversity. The importance of this process for the biological functions of plant S1-type enzymes is discussed.

Orthology prediction methods: a quality assessment using curated protein families

The increasing number of sequenced genomes has prompted the development of several automated orthology prediction methods. Tests to evaluate the accuracy of predictions and to explore biases caused by biological and technical factors are therefore required. We used 70 manually curated families to analyze the performance of five public methods in Metazoa. We analyzed the strengths and weaknesses of the methods and quantified the impact of biological and technical challenges. From the latter part of the analysis, genome annotation emerged as the largest single influencer, affecting up to 30% of the performance. Generally, most methods did well in assigning orthologous group but they failed to assign the exact number of genes for half of the groups. The publicly available benchmark set (http://eggnog.embl.de/orthobench/) should facilitate the improvement of current orthology assignment protocols, which is of utmost importance for many fields of biology and should be tackled by a broad scientific community.

More than 1,001 problems with protein domain databases: transmembrane regions, signal peptides and the issue of sequence homology

Large-scale genome sequencing gained general importance for life science because functional annotation of otherwise experimentally uncharacterized sequences is made possible by the theory of biomolecular sequence homology. Historically, the paradigm of similarity of protein sequences implying common structure, function and ancestry was generalized based on studies of globular domains. Having the same fold imposes strict conditions over the packing in the hydrophobic core requiring similarity of hydrophobic patterns. The implications of sequence similarity among non-globular protein segments have not been studied to the same extent; nevertheless, homology considerations are silently extended for them. This appears especially detrimental in the case of transmembrane helices (TMs) and signal peptides (SPs) where sequence similarity is necessarily a consequence of physical requirements rather than common ancestry. Thus, matching of SPs/TMs creates the illusion of matching hydrophobic cores. Therefore, inclusion of SPs/TMs into domain models can give rise to wrong annotations. More than 1001 domains among the 10,340 models of Pfam release 23 and 18 domains of SMART version 6 (out of 809) contain SP/TM regions. As expected, fragment-mode HMM searches generate promiscuous hits limited to solely the SP/TM part among clearly unrelated proteins. More worryingly, we show explicit examples that the scores of clearly false-positive hits, even in global-mode searches, can be elevated into the significance range just by matching the hydrophobic runs. In the PIR iProClass database v3.74 using conservative criteria, we find that at least between 2.1% and 13.6% of its annotated Pfam hits appear unjustified for a set of validated domain models. Thus, false-positive domain hits enforced by SP/TM regions can lead to dramatic annotation errors where the hit has nothing in common with the problematic domain model except the SP/TM region itself. We suggest a workflow of flagging problematic hits arising from SP/TM-containing models for critical reconsideration by annotation users.

Sequence homology is a fundamental principle of biology. It implies common phylogenetic ancestry of genes and, subsequently, similarity of their protein products with regard to amino acid sequence, three-dimensional structure and molecular and cellular function. Originally an esoteric concept, homology with the proxy of sequence similarity is used to justify the transfer of functional annotation from well-studied protein examples to new sequences. Yet, functional annotation via sequence similarity seems to have hit a plateau in recent years since relentless annotation transfer led to error propagation across sequence databases; thus, leading experimental follow-up work astray. It must be emphasized that the trinity of sequence, 3D structural and functional similarity has only been proven for globular segments of proteins. For non-globular regions, similarity of sequence is not necessarily a result of divergent evolution from a common ancestor but the consequence of amino acid sequence bias. In our investigation, we found that protein domain databases contain many domain models with transmembrane regions and signal peptides, non-globular segments of proteins having hydrophobic bias. Many proteins have inherited completely wrong function assignments from these domain models. We fear that future function predictions will turn out futile if this issue is not immediately addressed.

Quantifying the mechanisms of domain gain in animal proteins

Background

Protein domains are protein regions that are shared among different proteins and are frequently functionally and structurally independent from the rest of the protein. Novel domain combinations have a major role in evolutionary innovation. However, the relative contributions of the different molecular mechanisms that underlie domain gains in animals are still unknown. By using animal gene phylogenies we were able to identify a set of high confidence domain gain events and by looking at their coding DNA investigate the causative mechanisms.

Results

Here we show that the major mechanism for gains of new domains in metazoan proteins is likely to be gene fusion through joining of exons from adjacent genes, possibly mediated by non-allelic homologous recombination. Retroposition and insertion of exons into ancestral introns through intronic recombination are, in contrast to previous expectations, only minor contributors to domain gains and have accounted for less than 1% and 10% of high confidence domain gain events, respectively. Additionally, exonization of previously non-coding regions appears to be an important mechanism for addition of disordered segments to proteins. We observe that gene duplication has preceded domain gain in at least 80% of the gain events.

Conclusions

The interplay of gene duplication and domain gain demonstrates an important mechanism for fast neofunctionalization of genes.

The universal distribution of evolutionary rates of genes and distinct characteristics of eukaryotic genes of different apparent ages

The evolutionary rates of protein-coding genes in an organism span, approximately, 3 orders of magnitude and show a universal, approximately log-normal distribution in a broad variety of species from prokaryotes to mammals. This universal distribution implies a steady-state process, with identical distributions of evolutionary rates among genes that are gained and genes that are lost. A mathematical model of such process is developed under the single assumption of the constancy of the distributions of the propensities for gene loss (PGL). This model predicts that genes of different ages, that is, genes with homologs detectable at different phylogenetic depths, substantially differ in those variables that correlate with PGL. We computationally partition protein-coding genes from humans, flies, and Aspergillus fungus into age classes, and show that genes of different ages retain the universal log-normal distribution of evolutionary rates, with a shift toward higher rates in "younger" classes but also with a substantial overlap. The only exception involves human primate-specific genes that show a heavy tail of rapidly evolving genes, probably owing to gene annotation artifacts. As predicted, the gene age classes differ in characteristics correlated with PGL. Compared with "young" genes (e.g., mammal-specific human ones), "old" genes (e.g., eukaryote-specific), on average, are longer, are expressed at a higher level, possess a higher intron density, evolve slower on the short time scale, and are subject to stronger purifying selection. Thus, genome evolution fits a simple model with approximately uniform rates of gene gain and loss, without major bursts of genomic innovation.

Genome diversity in microbial eukaryotes

The genomic peculiarities among microbial eukaryotes challenge the conventional wisdom of genome evolution. Currently, many studies and textbooks explore principles of genome evolution from a limited number of eukaryotic lineages, focusing often on only a few representative species of plants, animals and fungi. Increasing emphasis on studies of genomes in microbial eukaryotes has and will continue to uncover features that are either not present in the representative species (e.g. hypervariable karyotypes or highly fragmented mitochondrial genomes) or are exaggerated in microbial groups (e.g. chromosomal processing between germline and somatic nuclei). Data for microbial eukaryotes have emerged from recent genome sequencing projects, enabling comparisons of the genomes from diverse lineages across the eukaryotic phylogenetic tree. Some of these features, including amplified rDNAs, subtelomeric rDNAs and reduced genomes, appear to have evolved multiple times within eukaryotes, whereas other features, such as absolute strand polarity, are found only within single lineages.

Correspondence of function and phylogeny of ABC proteins based on an automated analysis of 20 model protein data sets

Using our BLAST-based procedure RiPE (Retrieval-induced Phylogeny Environment), which automates the evolutionary analysis of a protein family, we assembled a set of 1138 ABC protein components [adenosine triphosphate (ATP)-binding cassette and transmembrane domain] from the protein data sets of 20 model organisms and subjected them to phylogenetic and functional analysis. For maximum speed, we based the alignment directly on a homology search with a profile of all known human ABC proteins and used neighbor-joining tree estimation. All but 11 sequences from Homo sapiens, Arabidopsis thaliana, Drosophila melanogaster, and Saccharomyces cerevisiae were placed into the correct subtree/subfamily, reproducing published classifications of the individual organisms. By following a simple "function transfer rule", our comparative phylogenetic analysis successfully predicted the known function of human ABC proteins in 19 of 22 cases. Three functional predictions did not correspond, and 10 were novel. Predictions based on BLAST alone were inferior in five cases and superior in two. Bacterial sequences were placed close to the root of most subtrees. This placement coincides with domain architecture, suggesting an early diversification of the ABC family before the kingdoms split apart. Our approach can, in principle, be used to annotate any protein family of any organism included in the study.

Characterization of an atypical antigen fromSarcoptes scabieicontaining an MADF domain

We have cloned a cDNA encoding a novel antigen from aSarcoptes scabiei(Acari) cDNA library by immunoscreening with sera fromS. scabiei-infected dogs. The antigen is encoded by a 2157 bp mRNA with a predicted open reading frame of 719 amino acids (molecular weight 79 kDa). Our sequence analysis identified the presence of a MADF domain in the N-terminus, and downstream of this domain there was a region of low sequence complexity. This latter region contained several blocks of triplets and quadruplets of polar amino acids (Asn, Gln and Ser), and these 3 amino acids represented 39·7% of all amino acids. The antigen was named Atypical Sarcoptes Antigen 1 (ASA1) since the MADF domain normally is found in proteins involved in transcriptional regulation. In addition, 15 out of 62S. scabiei-infected dogs reacted with a purified recombinant version of ASA1 in Western blot analysis. With immunohistochemistry we could show that ASA1 is expressed throughout the parasite, and that IgG specific for ASA1 binds to the inside wall of the mite's burrow. To our knowledge, this is the first description of an antigen containing an MADF domain.

Protein variety and functional diversity: Swiss-Prot annotation in its biological context

We all know that the dogma 'one gene, one protein' is obsolete. A functional protein and, likewise, a protein's ultimate function depend not only on the underlying genetic information but also on the ongoing conditions of the cellular system. Frequently the transcript, like the polypeptide, is processed in multiple ways, but only one or a few out of a multitude of possible variants are produced at a time. An overview on processes that can lead to sequence variety and structural diversity in eukaryotes is given. The UniProtKB/Swiss-Prot protein knowledgebase provides a wealth of information regarding protein variety, function and associated disorders. Examples for such annotation are shown and further ones are available at http://www.expasy.org/sprot/tutorial/examples_CRB.

Evolutionary diversity of vertebrate small heat shock proteins

All vertebrates express multiple small heat shock proteins (sHsps), which are important components of the cellular chaperoning machinery and display a spectacular diversity of functions. This ranges from remodeling the cytoskeleton and inhibiting apoptosis to serving as structural proteins in eye lens and sperm tail. Most information is available for the 10 known mammalian sHsps, formally named HspB1-B10. Only three of them (Hsp27/B1, alphaA-crystallin/B4, alphaB-crystallin/B5) have been reported from nonmammalian vertebrates, while an apparent paralog, Hsp30/B11, is found in frogs and teleost fish. To reconstruct the evolutionary diversification of the sHsps in vertebrates, we searched for additional sHsps in genome, protein, and EST databases and sequenced some avian and amphibian sHsps (HspB2, Hsp30/B11). The urochordate Ciona intestinalis was included in the search, as the outgroup of vertebrates. Orthologs of seven mammalian sHsps were now found in other vertebrate classes. Two novel sHsps, named HspB11 and HspB12, were recognized in birds, and four novel sHsps, named HspB12-B15, in teleost fish. Secondary structure predictions of orthologous sHsps from different vertebrate classes indicate conservation of the beta-sandwich structure of the functionally important C-terminal "alpha-crystallin domain," while the N-terminal domains generally have alpha-helical structures, despite their pronounced sequence variation. The constructed chordate sHsp tree is supported by shared introns, indels, and diagnostic sequences. The tree distinguishes putative orthologous and paralogous relationships, which will facilitate the functional and structural comparison of the various vertebrate sHsps. The 15 recognized paralogous vertebrate sHsps reflect the period of extensive gene duplications early in vertebrate evolution. Eleven of these sHsps are grouped in a clade that might be specific for chordates. It is inferred that at least 13 intron insertions have occurred during the evolution of chordate sHsp genes, while a single ancient intron is maintained in some lineages, in line with the general trend of massive intron gain before or during early vertebrate radiation. Interesting is the occurrence of several head-to-head located pairs of chordate sHsp genes.

Analysis of the human VPS13 gene family

The gene mutated in chorea-acanthocytosis (CHAC; approved gene symbol VPS13A) encodes chorein, a protein similar to yeast Vps13p. We detected several similar putative human proteins by BLAST analysis of chorein. We characterized the structure of three new genes encoding these CHAC-similar proteins, located on chromosomes 1p36, 8q22, and 15q21. The most similar gene in yeast to all four human genes is Vps13, and therefore the human genes were named VPS13A (CHAC, 9q21), VPS13B (8q22), VPS13C (15q21), and VPS13D (1p36). VPS13B has recently been reported as COH1, altered in Cohen syndrome. For each gene, we describe several alternative splicing variants; at least two transcripts per gene are major forms. The expression pattern of these genes is ubiquitous, with some tissue-specific differences between several transcript variants. Protein sequence comparisons suggest that intramolecular duplications have played an important role in the evolution of this gene family.

Expression analysis of ABC transporters reveals differential functions of tandemly duplicated genes in Caenorhabditis elegans

We have previously identified 60 predicted ABC transporter genes in the Caenorhabditis elegans genome and classified them into eight groups. As an initial step towards understanding how these putative ABC genes work in worms, we generated promoter-fluorescent protein fusions for the entire family to address when and where these genes are turned on in vivo. Both Aequoria green fluorescent protein (GFP) and Discosoma red fluorescent protein (RFP) were used as reporters in our transgenic assay. Observable expression is more frequently seen from fusions to genes in subfamilies B, C, D and E than those in subfamilies A and G. Sixteen worm ABC genes are found in tandem duplications, forming two four-gene clusters and four two-gene clusters. Fifteen out of the 16 duplicated gene promoters drove different or partially overlapping expression patterns, suggesting active functions for these duplicated genes. Furthermore, our results suggest that an internal promoter can cause differential expression of genes within an operon. Finally, our observations suggest that it is possible for coding sequences to function as a regulatory region for a neighbouring gene.

Evaluation of monocot and eudicot divergence using the sugarcane transcriptome

Over 40,000 sugarcane (Saccharum officinarum) consensus sequences assembled from 237,954 expressed sequence tags were compared with the protein and DNA sequences from other angiosperms, including the genomes of Arabidopsis and rice (Oryza sativa). Approximately two-thirds of the sugarcane transcriptome have similar sequences in Arabidopsis. These sequences may represent a core set of proteins or protein domains that are conserved among monocots and eudicots and probably encode for essential angiosperm functions. The remaining sequences represent putative monocot-specific genetic material, one-half of which were found only in sugarcane. These monocot-specific cDNAs represent either novelties or, in many cases, fast-evolving sequences that diverged substantially from their eudicot homologs. The wide comparative genome analysis presented here provides information on the evolutionary changes that underlie the divergence of monocots and eudicots. Our comparative analysis also led to the identification of several not yet annotated putative genes and possible gene loss events in Arabidopsis.

Lesion bypass DNA polymerases replicate across non-DNA segments

A critical feature of the robustness of the DNA replication machinery is the ability to complete its task in the presence of interfering DNA damage. A key mechanism responsible for this task is translesion replication (also termed translesion synthesis), carried out by specialized lesion bypass DNA polymerases of the Y superfamily. Here we show that in Escherichia coli, plasmids can be replicated across a segment of foreign non-DNA material, consisting of hydrocarbon chains of 3 or 12 methylene residues. This replication is carried out by DNA polymerase V and proceeds by at least two mechanisms: (i) Editing out the foreign insert, by polymerase "hopping" across it, which can be mediated by looping out of the insert, leading to its deletion, while preserving the DNA sequence. (ii) DNA synthesis through the insert, which occurs by incorporating one or two nucleotides opposite the hydrocarbon chain, yielding a net increase in the length of the DNA sequence. The remarkable ability of DNA polymerase V to insert nucleotides opposite a hydrocarbon chain shows that DNA synthesis can occur in a region of the template strand, which lacks all fundamental features of DNA, including its purine, pyrimidine, sugar, and phosphate moieties, and its hydrophilic and ionic nature. This bypass ability reflects a striking robustness of the translesion replication apparatus and is likely to contribute to its effectiveness in maintaining genome stability.

Metabolome diversity: too few genes, too many metabolites?

The multitude of metabolites found in living organisms and the calculated, unexpected small number of genes identified during genome sequencing projects discomfit biologists. Several processes on the transcription and translation level lead to the formation of isoenzymes and can therefore explain at least parts of this surprising result. However, poor enzyme specificity may also contribute to metabolome diversity. In former studies, when enzymes were isolated from natural sources, impure protein preparations were hold responsible for broad enzyme specificity. Nowadays, highly purified enzymes are available by molecular biological methods such as heterologous expression in host organisms and they can be thoroughly analyzed. During biochemical analysis of heterologously expressed enzymes poor specificity was observed for enzymes involved in fruit ripening, e.g. in flavour and color formation. Surprisingly broad specificity was shown for the reactants in the case of alcohol acyl-CoA transferase, O-methyltransferase, glucosyltransferase, P450 monooxygenases as well as polyketide synthases and for the product in the case of monoterpene synthases. Literature data confirm the assumption of limited specificity for enzymes involved in metabolism and bioformation of secondary metabolites. It is concluded that metabolome diversity is caused by low enzyme specificity but availability of suitable substrates due to compartmentation has also taken into account.