Gene duplication within the Green Lineage: the case of TEL genes

Recent years have witnessed a breathtaking increase in the availability of genome sequence data, providing evidence of the highly duplicate nature of eukaryotic genomes. Plants are exceptional among eukaryotic organisms in that duplicate loci compose a large fraction of their genomes, partly because of the frequent occurrence of polyploidy (or whole-genome duplication) events. Tandem gene duplication and transposition have also contributed to the large number of duplicated genes in plant genomes. Evolutionary analyses allowed the dynamics of duplicate gene evolution to be studied and several models were proposed. It seems that, over time, many duplicated genes were lost and some of those that were retained gained new functions and/or expression patterns (neofunctionalization) or subdivided their functions and/or expression patterns between them (subfunctionalization). Recent studies have provided examples of genes that originated by duplication with successive diversification within plants. In this review, we focused on the TEL (TERMINAL EAR1-like) genes to illustrate such mechanisms. Emerged from the mei2 gene family, these TEL genes are likely to be land plant-specific. Phylogenetic analyses revealed one or two TEL copies per diploid genome. TEL gene degeneration and loss in several Angiosperm species such as in poplar and maize seem to have occurred. In Arabidopsis thaliana, whose genome experienced at least three polyploidy events followed by massive gene loss and genomic reorganization, two TEL genes were retained and two new shorter TEL-like (MCT) genes emerged. Molecular and expression analyses suggest for these genes sub- and neofunctionalization events, but confirmation will come from their functional characterization.     

Genome comparisons highlight similarity and diversity within the eukaryotic kingdoms

In 2000, the number of completely sequenced eukaryotic genomes increased to four. The addition of Drosophila and Arabidopsis into this cohort permits additional insights into the processes that have shaped evolution. Analysis and comparisons of both completed genomes and partially sequenced genomes have already shed light on mechanisms such as gene duplication and gene loss that have long been hypothesized to be major forces in speciation. Indeed, duplicate gene pairs in Saccharomyces, Arabidopsis, Caenorhabditis and Drosophila are high: 30%, 60%, 48% and 40%, respectively. Evidence of horizontal gene-transfer, thought to be a major evolutionary force in bacteria, has been found in Arabidopsis. The release of the 'first draft' of the human genome sequence in 2000 heralds a new stage of biological study. Understanding the as-yet-unannotated human genome will be largely based on conclusions, techniques and tools developed during the analysis and comparison of the genome of these four model organisms.     

Birth and death of duplicated genes in completely sequenced eukaryotes

Gene and genome duplications are commonly regarded as being of major evolutionary significance. But how often does gene duplication occur? And, once duplicated, what are the fates of duplicated genes? How do they contribute to evolution? In a recent article, Lynch and Conery analyze divergence between duplicate genes from six eukaryotic genomes. They estimate the rate of gene duplication, the rate of gene loss after duplication and the strength of selection experienced by duplicate genes. They conclude that although the rate of gene duplications is high, so is the rate of gene loss, and they argue that gene duplications could be a major factor in speciation.     

Evolution of duplicate control regions in the mitochondrial genomes of metazoa: a case study with Australasian Ixodes ticks

To investigate the evolution pattern and phylogenetic utility of duplicate control regions (CRs) in mitochondrial (mt) genomes, we sequenced the entire mt genomes of three Ixodes species and part of the mt genomes of another 11 species. All the species from the Australasian lineage have duplicate CRs, whereas the other species have one CR. Sequence analyses indicate that the two CRs of the Australasian Ixodes ticks have evolved in concert in each species. In addition to the Australasian Ixodes ticks, species from seven other lineages of metazoa also have mt genomes with duplicate CRs. Accumulated mtDNA sequence data from these metazoans and two recent experiments on replication of mt genomes in human cell lines with duplicate CRs allowed us to re-examine four intriguing questions about the presence of duplicate CRs in the mt genomes of metazoa: (1) Why do some mt genomes, but not others, have duplicate CRs? (2) How did mt genomes with duplicate CRs evolve? (3) How could the nucleotide sequences of duplicate CRs remain identical or very similar over evolutionary time? (4) Are duplicate CRs phylogenetic markers? It appears that mt genomes with duplicate CRs have a selective advantage in replication over mt genomes with one CR. Tandem duplication followed by deletion of genes is the most plausible mechanism for the generation of mt genomes with duplicate CRs. Once duplicate CRs occur in an mt genome, they tend to evolve in concert, probably by gene conversion. However, there are lineages where gene conversion may not always occur, and, thus, the two CRs may evolve independently in these lineages. Duplicate CRs have much potential as phylogenetic markers at low taxonomic levels, such as within genera, within families, or among families, but not at high taxonomic levels, such as among orders.     

Helitrons on a roll: eukaryotic rolling-circle transposons

Rolling-circle eukaryotic transposons, known as Helitron transposons, were first discovered in plants (Arabidopsis thaliana and Oryza sativa) and in the nematode Caenorhabditis elegans. To date, Helitrons have been identified in a diverse range of species, from protists to mammals. They represent a major class of eukaryotic transposons and are fundamentally different from classical transposons in terms of their structure and mechanism of transposition. Helitrons seem to have a major role in the evolution of host genomes. They frequently capture diverse host genes, some of which can evolve into novel host genes or become essential for helitron transposition.     

Extent of gene duplication in the genomes of Drosophila, nematode, and yeast

We conducted a detailed analysis of duplicate genes in three complete genomes: yeast, Drosophila, and Caenorhabditis elegans. For two proteins belonging to the same family we used the criteria: (1) their similarity is > or =I (I = 30% if L > or = 150 a.a. and I = 0.01n + 4.8L(-0.32(1 + exp(-L/1000))) if L < 150 a.a., where n = 6 and L is the length of the alignable region), and (2) the length of the alignable region between the two sequences is > or = 80% of the longer protein. We found it very important to delete isoforms (caused by alternative splicing), same genes with different names, and proteins derived from repetitive elements. We estimated that there were 530, 674, and 1,219 protein families in yeast, Drosophila, and C. elegans, respectively, so, as expected, yeast has the smallest number of duplicate genes. However, for the duplicate pairs with the number of substitutions per synonymous site (K(S)) < 0.01, Drosophila has only seven pairs, whereas yeast has 58 pairs and nematode has 153 pairs. After considering the possible effects of codon usage bias and gene conversion, these numbers became 6, 55, and 147, respectively. Thus, Drosophila appears to have much fewer young duplicate genes than do yeast and nematode. The larger numbers of duplicate pairs with K(S) < 0.01 in yeast and C. elegans were probably largely caused by block duplications. At any rate, it is clear that the genome of Drosophila melanogaster has undergone few gene duplications in the recent past and has much fewer gene families than C. elegans.     

Short-chain dehydrogenase/reductase (SDR) relationships: a large family with eight clusters common to human, animal, and plant genomes

The progress in genome characterizations has opened new routes for studying enzyme families. The availability of the human genome enabled us to delineate the large family of short-chain dehydrogenase/reductase (SDR) members. Although the human genome releases are not yet final, we have already found 63 members. We have also compared these SDR forms with those of three model organisms: Caenorhabditis elegans, Drosophila melanogaster, and Arabidopsis thaliana. We detect eight SDR ortholog clusters in a cross-genome comparison. Four of these clusters represent extended SDR forms, a subgroup found in all life forms. The other four are classical SDRs with activities involved in cellular differentiation and signalling. We also find 18 SDR genes that are present only in the human genome of the four genomes studied, reflecting enzyme forms specific to mammals. Close to half of these gene products represent steroid dehydrogenases, emphasizing the regulatory importance of these enzymes.     

Tc8, a Tourist-like transposon in Caenorhabditis elegans.

Members of the Tourist family of miniature inverted-repeat transposable elements (MITEs) are very abundant among a wide variety of plants, are frequently found associated with normal plant genes, and thus are thought to be important players in the organization and evolution of plant genomes. In Arabidopsis, the recent discovery of a Tourist member harboring a putative transposase has shed new light on the mobility and evolution of MITEs. Here, we analyze a family of Tourist transposons endogenous to the genome of the nematode Caenorhabditis elegans (Bristol N2). One member of this large family is 7568 bp in length, harbors an ORF similar to the putative Tourist transposase from Arabidopsis, and is related to the IS5 family of bacterial insertion sequences (IS). Using database searches, we found expressed sequence tags (ESTs) similar to the putative Tourist transposases in plants, insects, and vertebrates. Taken together, our data suggest that Tourist-like and IS5-like transposons form a superfamily of potentially active elements ubiquitous to prokaryotic and eukaryotic genomes.     

Evolutionary dynamics of introns in plastid-derived genes in plants: saturation nearly reached but slow intron gain continues

Some of the principal transitions in the evolution of eukaryotes are characterized by engulfment of prokaryotes by primitive eukaryotic cells. In particular, approximately 1.6 billion years ago, engulfment of a cyanobacterium that became the ancestor of chloroplasts and other plastids gave rise to Plantae, the major branch of eukaryotes comprised of glaucophytes, red algae, green algae, and green plants. After endosymbiosis, there was large-scale migration of genes from the endosymbiont to the nuclear genome of the host such that approximately 18% of the nuclear genes in Arabidopsis appear to be of chloroplast origin. To gain insights into the process of evolution of gene structure in these, originally, intronless genes, we compared the properties and the evolutionary dynamics of introns in genes of plastid origin and ancestral eukaryotic genes in Arabidopsis, poplar, and rice genomes. We found that intron densities in plastid-derived genes were slightly but significantly lower than those in ancestral eukaryotic genes. Although most of the introns in both categories of genes were conserved between monocots (rice) and dicots (Arabidopsis and poplar), lineage-specific intron gain was more pronounced in plastid-derived genes than in ancestral genes, whereas there was no significant difference in the intron loss rates between the 2 classes of genes. Thus, after the transfer to the nuclear genome, the plastid-derived genes have undergone a massive intron invasion that, by the time of the divergence of dicots and monocots (150-200 MYA), yielded intron densities only slightly lower than those in ancestral genes. Nevertheless, the accumulation of introns in plastid-derived genes appears not to have reached saturation and continues to this time, albeit at a low rate. The overall pattern of intron gain and loss in the plastid-derived genes is shaped by this continuing gain and the more general tendency for loss that is characteristic of the recent evolution of plant genes.     

Distinct frequency-distributions of homopolymeric DNA tracts in different genomes.

The unusual base composition of the genome of the human malaria parasite Plasmodium falciparum prompted us to systematically investigate the occurrence of homopolymeric DNA tracts in the P. falciparum genome and, for comparison, in the genomes of Homo sapiens , Saccharomyces cerevisiae , Caenorhabditis elegans , Arabidopsis thaliana , Escherichia coli and Mycobacterium tuberculosis. Comparison of theobserved frequencies with the frequencies as expected for random DNA revealed that homopolymeric (dA:dT) tracts occur well above chance in the eukaryotic genome. In the majority of these genomes, (dA:dT) tract overrepresentation proved to be an exponential function of the tract length. (dG:dC) tract overrepresentation was absent or less pronounced in both prokaryotic and eukaryotic genomes. On the basis of our results, we propose that homopolymeric (dA:dT) tracts are expanded via replication slippage. This slippage-mediated expansion does not operate on tracts with lengths below a critical threshold of 7-10 bp.     

双翅目昆虫（黑腹果蝇和冈比亚按蚊）内含子丢失的比较分析

内含子插入和丢失的进化动力及机制尚存在许多疑问。通过对真核生物的105个同源基因的蛋白质高度保守区域内含子-外显子结构的研究,对人Homo sapiens、小鼠Mus musculus、大鼠Rattus norvegicus、黑腹果蝇Drosophila melanogaster、冈比亚按蚊Anopheles gambiae和秀丽隐杆线虫Caenorhabditis elegans的3 574个内含子、1 001个的内含子保守位点进行分析,推断出不同系统中内含子的变化途径。发现在进化早期,脊椎动物、双翅目昆虫和线虫的共同祖先中含有大量内含子,在进化过程中,双翅目昆虫和线虫发生了大量的内含子丢失,甚至在双翅目昆虫中内含子丢失较线虫更严重。线虫获得的内含子略多于丢失的内含子, 而在双翅目昆虫中则显示出内含子的丢失明显多于内含子的获得。该结果合理地解释了内含子在脊椎动物、线虫及昆虫中数量的分布呈下降趋势。     

Phylogenetic and transcriptional analysis of a strictosidine synthase-like gene family in Arabidopsis thaliana reveals involvement in plant defence responses

Protein domains with similarity to plant strictosidine synthase-like (SSL) sequences have been uncovered in the genomes of all multicellular organisms sequenced so far and are known to play a role in animal immune responses. Among several distinct groups of Arabidopsis thaliana SSL sequences, four genes ( AtSSL4–AtSSL7 ) arranged in tandem on chromosome 3 show more similarity to SSL genes from Drosophila melanogaster and Caenorhabditis elegans than to other Arabidopsis SSL genes. To examine whether any of the four AtSSL genes are immune-inducible, we analysed the expression of each of the four AtSSL genes after exposure to microbial pathogens, wounding and plant defence elicitors using real-time quantitative RT-PCR, Northern blot hybridisation and Western blot analysis with antibodies raised against recombinant At SSL proteins. While the AtSSL4 gene was constitutively expressed and not significantly induced by any treatment, the other three AtSSL genes were induced to various degrees by plant defence signalling compounds, such as salicylic acid, methyl jasmonate and ethylene, as well as by wounding and exposure to the plant pathogens Alternaria brassicicola and cucumber mosaic virus . Our data demonstrate that the four SSL-coding genes are regulated individually, suggesting specific roles in basal ( SSL4 ) and inducible ( SSL5-7 ) plant defence mechanisms.     

Recent advances in gene structure prediction

De novo gene predictors are programs that predict the exon-intron structures of genes using the sequences of one or more genomes as their only input. In the past two years, dual-genome de novo predictors, which exploit local rates and patterns of mutation inferred from alignments between two genomes, have led to significant improvements in accuracy. Systems that exploit more than two genomes simultaneously have only recently begun to appear and are not yet competitive on practical tasks, but offer the greatest hope for near-term improvements. Dual-genome de novo prediction for compact eukaryotic genomes such as those of Arabidopsis thaliana and Caenorhabditis elegans is already quite accurate. Although mammalian gene prediction lags behind in accuracy, it is yielding ever more useful results. Coupled with significant improvements in pseudogene detection methods, which have eliminated many false positives, we have reached the point where de novo gene predictions are being used as hypotheses to drive experimental annotation via systematic RT-PCR and sequencing.     

Multiple ribonuclease H-encoding genes in the Caenorhabditis elegans genome contrasts with the two typical ribonuclease H-encoding genes in the human genome

Database searches of the Caenorhabditis elegans and human genomic DNA sequences revealed genes encoding ribonuclease H1 (RNase H1) and RNase H2 in each genome. The human genome contains a single copy of each gene, whereas C. elegans has four genes encoding RNase H1-related proteins and one gene for RNase H2. By analyzing the mRNAs produced from the C. elegans genes, examining the amino acid sequence of the predicted protein, and expressing the proteins in Esherichia coli we have identified two active RNase H1-like proteins. One is similar to other eukaryotic RNases H1, whereas the second RNase H (rnh-1.1) is unique. The rnh-1.0 gene is transcribed as a dicistronic message with three dsRNA-binding domains; the mature mRNA is transspliced with SL2 splice leader and contains only one dsRNA-binding domain. Formation of RNase H1 is further regulated by differential cis-splicing events. A single rnh-2 gene, encoding a protein similar to several other eukaryotic RNase H2L's, also has been examined. The diversity and enzymatic properties of RNase H homologues are other examples of expansion of protein families in C. elegans. The presence of two RNases H1 in C. elegans suggests that two enzymes are required in this rather simple organism to perform the functions that are accomplished by a single enzyme in more complex organisms. Phylogenetic analysis indicates that the active C. elegans RNases H1 are distantly related to one another and that the C. elegans RNase H1 is more closely related to the human RNase H1. The database searches also suggest that RNase H domains of LTR-retrotransposons in C. elegans are quite unrelated to cellular RNases H1, but numerous RNase H domains of human endogenous retroviruses are more closely related to cellular RNases H.