Molecular clock and gene function

Molecular phylogenies based on the molecular clock require the comparison of orthologous genes. Orthologous and paralogous genes usually have very different evolutionary fates. In general, orthologs keep the same functions in species, whereas, particularly over a long time span, paralogs diverge functionally and may become pseudogenes or get lost. In eukaryotic genomes, because of the degree of redundancy of genetic information, homologous genes are grouped in gene families, the evolution of which may differ greatly between the various organisms. This implies that each gene in a species does not always have an ortholog in another species and thus, due to multiple duplication events following a speciation, many orthologous clades of paralogs are generated. We are often dealing with a one-to-many or many-to-many relationship between genes. In this paper, we analyze the evolution of two gene families, the p53 gene family and the porin gene family. The evolution of the p53 family shows a one-to-many gene relationship going from invertebrates to vertebrates. In invertebrates only a single gene has been found, while in vertebrates three members of the family, namely p53, p63, and p73, are present. The evolution of porin (VDAC) genes (VDAC1, VDAC2, and VDAC3) is an example of a many-to-many gene relationship going from yeast to mammals. However, the porin gene redundancy found in invertebrates and possibly in some fishes may indicate a tendency to duplicate the genetic material, rather than a real need for function innovation.     

Loss of phylogenetic information in chorion gene families of Bombyx mori gene conversion

The silkmoth chorion has provided a stimulating model for the study of evolution and developmental regulation of gene families. Previous attempts at inferring relationships among chorion sequences have been based on pairwise comparisons of overall similarity, a potentially problematic approach. To remedy this, we identified the alignable regions of low sequence variability and then analyzed this restricted database by parsimony and neighbor-joining methods. At the deepest level, the chorion sequence tree is split into two branches, called "alpha" and "beta." Within each branch, early- and late-expressing genes each constitute monophyletic groups, while the situation with middle-expressing genes remains uncertain. The HcB gene family appears to be the most basal beta-branch group, but this conclusion is qualified because the effect of gene conversion on branching order is unknown. Previous studies by Eickbush and colleagues have strongly suggested that ErA, HcA, and HcB families undergo gene conversion within a gene family, whereas the ErB family does not. The occurrence of conversion correlates with a particular tree structure; namely, branch lengths are much greater at the base of the family than at higher internodes and terminal branches. These observations raise the possibility that chorion gene families are defined by gene conversion events (reticulate evolution) rather than by descent with modification (synapomorphy).      

Concerted evolution of sea anemone neurotoxin genes is revealed through analysis of the Nematostella vectensis genome

Gene families, which encode toxins, are found in many poisonous animals, yet there is limited understanding of their evolution at the nucleotide level. The release of the genome draft sequence for the sea anemone Nematostella vectensis enabled a comprehensive study of a gene family whose neurotoxin products affect voltage-gated sodium channels. All gene family members are clustered in a highly repetitive approximately 30-kb genomic region and encode a single toxin, Nv1. These genes exhibit extreme conservation at the nucleotide level which cannot be explained by purifying selection. This conservation greatly differs from the toxin gene families of other animals (e.g., snakes, scorpions, and cone snails), whose evolution was driven by diversifying selection, thereby generating a high degree of genetic diversity. The low nucleotide diversity at the Nv1 genes is reminiscent of that reported for DNA encoding ribosomal RNA (rDNA) and 2 hsp70 genes from Drosophila, which have evolved via concerted evolution. This evolutionary pattern was experimentally demonstrated in yeast rDNA and was shown to involve unequal crossing-over. Through sequence analysis of toxin genes from multiple N. vectensis populations and 2 other anemone species, Anemonia viridis and Actinia equina, we observed that the toxin genes for each sea anemone species are more similar to one another than to those of other species, suggesting they evolved by manner of concerted evolution. Furthermore, in 2 of the species (A. viridis and A. equina) we found genes that evolved under diversifying selection, suggesting that concerted evolution and accelerated evolution may occur simultaneously.     

Numbers of genes in the NBS and RLK families vary by more than four-fold within a plant species and are regulated by multiple factors

Many genes exist in the form of families; however, little is known about their size variation, evolution and biology. Here, we present the size variation and evolution of the nucleotide-binding site (NBS)-encoding gene family and receptor-like kinase (RLK) gene family in Oryza, Glycine and Gossypium. The sizes of both families vary by numeral fold, not only among species, surprisingly, also within a species. The size variations of the gene families are shown to correlate with each other, indicating their interactions, and driven by natural selection, artificial selection and genome size variation, but likely not by polyploidization. The numbers of genes in the families in a polyploid species are similar to those of one of its diploid donors, suggesting that polyploidization plays little roles in the expansion of the gene families and that organisms tend not to maintain their ‘surplus’ genes in the course of evolution. Furthermore, it is found that the size variations of both gene families are associated with organisms’ phylogeny, suggesting their roles in speciation and evolution. Since both selection and speciation act on organism’s morphological, physiological and biological variation, our results indicate that the variation of gene family size provides a source of genetic variation and evolution.     

Structure and Evolution of the Actin Gene Family in Arabidopsis Thaliana

Higher plants contain families of actin-encoding genes that are divergent and differentially expressed. Progress in understanding the functions and evolution of plant actins has been hindered by the large size of the actin gene families. In this study, we characterized the structure and evolution of the actin gene family in Arabidopsis thaliana. DNA blot analyses with gene-specific probes suggested that all 10 of the Arabidopsis actin gene family members have been isolated and established that Arabidopsis has a much simpler actin gene family than other plants that have been examined. Phylogenetic analyses suggested that the Arabidopsis gene family contains at least two ancient classes of genes that diverged early in land plant evolution and may have separated vegetative from reproductive actins. Subsequent divergence produced a total of six distinct subclasses of actin, and five showed a distinct pattern of tissue specific expression. The concordance of expression patterns with the phylogenetic structure is discussed. These subclasses appear to be evolving independently, as no evidence of gene conversion was found. The Arabidopsis actin proteins have an unusually large number of nonconservative amino acid substitutions, which mapped to the surface of the actin molecule, and should effect protein-protein interactions.     

Evolution of the CYP2ABFGST gene cluster in rat, and a fine-scale comparison among rodent and primate species

The evolution of gene families can be best understood by studying the modern organization and functions of family members, and by comparing parallel families in different species. In this study, the CYP2ABFGST gene cluster has been characterized in rat and compared to the syntenic clusters in mouse and human, providing an interesting example of gene family evolution. In the rat, 18 loci from six subfamilies have been identified by specifically amplifying and sequencing gene fragments from cloned DNA, and have been exactly placed on chromosome 1. The overall organization of the gene cluster in rat is relatively simple, with genes from each subfamily in tandem, and is more similar to the mouse than to the human cluster. We have reconstructed the probable structure of the CYP2ABFGST cluster in the common ancestor of primates and rodents, and inferred a model of the evolution of this gene cluster in the three species. Numerous nontandem and block duplications, inversions, and translocations have occurred entirely inside the cluster, indicating that pairing between duplicate genes is keeping the rearrangements within the cluster region. The initial tandem duplication of a CYP2 gene in an early mammalian ancestor has made this region particularly subject to such localized rearrangements. Even if duplicated genes do not have a large-scale effect on chromosomal rearrangements, on a local level clustered gene families may have contributed significantly to the genomic complexity of modern mammals.     

A Survey of the Molecular Evolutionary Dynamics of Twenty-Five Multigene Families from Four Grass Taxa

We surveyed the molecular evolutionary characteristics of 25 plant gene families, with the goal of better understanding general processes in plant gene family evolution. The survey was based on 247 GenBank sequences representing four grass species (maize, rice, wheat, and barley). For each gene family, orthology and paralogy relationships were uncertain. Recognizing this uncertainty, we characterized the molecular evolution of each gene family in four ways. First, we calculated the ratio of nonsynonymous to synonymous substitutions (d _N/d _S) both on branches of gene phylogenies and across codons. Our results indicated that the d _N/d _S ratio was statistically heterogeneous across branches in 17 of 25 (68%) gene families. The vast majority of d _N/d _S estimates were <<1.0, suggestive of selective constraint on amino acid replacements, and no estimates were >1.0, either across phylogenetic lineages or across codons. Second, we tested separately for nonsynonymous and synonymous molecular clocks. Sixty-eight percent of gene families rejected a nonsynonymous molecular clock, and 52% of gene families rejected a synonymous molecular clock. Thus, most gene families in this study deviated from clock-like evolution at either synonymous or nonsynonymous sites. Third, we calculated the effective number of codons and the proportion of G+C synonymous sites for each sequence in each gene family. One or both quantities vary significantly within 18 of 25 gene families. Finally, we tested for gene conversion, and only six gene families provided evidence of gene conversion events. Altogether, evolution for these 25 gene families is marked by selective constraint that varies among gene family members, a lack of molecular clock at both synonymous and nonsynonymous sites, and substantial variation in codon usage. Received: 25 May 2000 / Accepted: 16 October 2000     

The life and death of gene families

One of the unique insights provided by the growing number of fully sequenced genomes is the pervasiveness of gene duplication and gene loss. Indeed, several metrics now suggest that rates of gene birth and death per gene are only 10–40% lower than nucleotide substitutions per site, and that per nucleotide, the consequent lineage‐specific expansion and contraction of gene families may play at least as large a role in adaptation as changes in orthologous sequences. While gene family evolution is pervasive, it may be especially important in our own evolution since it appears that the “revolving door” of gene duplication and loss has undergone multiple accelerations in the lineage leading to humans. In this paper, we review current understanding of gene family evolution including: methods for inferring copy number change, evidence for adaptive expansion and adaptive contraction of gene families, the origins of new families and deaths of previously established ones, and finally we conclude with a perspective on challenges and promising directions for future research.     

Structure, map position, and evolution of two newly diverged mouse Ig VH gene families.

We have characterized two novel mouse VH gene families, VH3609N and VHSM7. These VH families have recently diverged from previously defined VH families. The VH3609N family, which may contain only one member in most inbred strains of mice, shares sequence similarity with the VHJ606 family and is located to the 3' side of VHJ606. VHSM7, with at least three members, is related to the VHJ558 family but maps 3' of VHJ558. These findings suggest that physical displacement of VH sequences may facilitate their subsequent divergence. During the early stages of VH gene family evolution that are exemplified by these new families, amino acid replacements have been selected against in frame-work regions and selected for in complementarity-determining regions. This pattern of nucleotide substitution appears to reflect evolutionary pressures to maintain germ-line VH diversity and, possibly, to select for new antibody specificities, as well as to select against mutations resulting in aberrant Ig. The classification of VH sequences with borderline similarity to previously defined VH families is discussed.     

The roles of segmental and tandem gene duplication in the evolution of large gene families in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Background

Most genes in Arabidopsis thaliana are members of gene families. How do the members of gene families arise, and how are gene family copy numbers maintained? Some gene families may evolve primarily through tandem duplication and high rates of birth and death in clusters, and others through infrequent polyploidy or large-scale segmental duplications and subsequent losses.

Results

Our approach to understanding the mechanisms of gene family evolution was to construct phylogenies for 50 large gene families in Arabidopsis thaliana, identify large internal segmental duplications in Arabidopsis, map gene duplications onto the segmental duplications, and use this information to identify which nodes in each phylogeny arose due to segmental or tandem duplication. Examples of six gene families exemplifying characteristic modes are described. Distributions of gene family sizes and patterns of duplication by genomic distance are also described in order to characterize patterns of local duplication and copy number for large gene families. Both gene family size and duplication by distance closely follow power-law distributions.

Conclusions

Combining information about genomic segmental duplications, gene family phylogenies, and gene positions provides a method to evaluate contributions of tandem duplication and segmental genome duplication in the generation and maintenance of gene families. These differences appear to correspond meaningfully to differences in functional roles of the members of the gene families.

     

Mass-murder deletion of 19 ORFs from Saccharomyces cerevisiae chromosome XI

In eukaryote genomes, there are many kinds of gene families. Gene duplication and conversion are sources of the evolution of gene families, including those with uniform members and those with diverse functions. Population genetics theory on identity coefficients among gene members of a gene family shows that the balance between diversification by mutation, and homogenization by unequal crossing over and gene conversion, is important. Also, evolution of new functions is due to gene duplication followed by differentiation. Positive selection is necessary for the evolution of novel functions. However, many examples of current gene families suggest that both drift and selection are at work on their evolution.     

A shared promoter region suggests a common ancestor for the human VCX/Y, SPANX, and CSAG gene families and the murine CYPT family

Many testis-specific genes from the sex chromosomes are subject to rapid evolution, which can make it difficult to identify murine genes in the human genome. The murine CYPT gene family includes 15 members, but orthologs were undetectable in the human genome. However, using refined homology search, sequences corresponding to the shared promoter region of the CYPT family were identified at 39 loci. Most loci were located immediately upstream of genes belonging to the VCX/Y, SPANX, or CSAG gene families. Sequence comparison of the loci revealed a conserved CYPT promoter-like (CPL) element featuring TATA and CCAAT boxes. The expression of members of the three families harboring the CPL resembled the murine expression of the CYPT family, with weak expression in late pachytene spermatocytes and predominant expression in spermatids, but some genes were also weakly expressed in somatic cells and in other germ cell types. The genomic regions harboring the gene families were rich in direct and inverted segmental duplications (SD), which may facilitate gene conversion and rapid evolution. The conserved CPL and the common expression profiles suggest that the human VCX/Y, SPANX, and CSAG2 gene families together with the murine SPANX gene and the CYPT family may share a common ancestor. Finally, we present evidence that VCX/Y and SPANX may be paralogs with a similar protein structure consisting of C terminal acidic repeats of variable lengths.     

CAFE: a computational tool for the study of gene family evolution

SUMMARY: We present CAFE (Computational Analysis of gene Family Evolution), a tool for the statistical analysis of the evolution of the size of gene families. It uses a stochastic birth and death process to model the evolution of gene family sizes over a phylogeny. For a specified phylogenetic tree, and given the gene family sizes in the extant species, CAFE can estimate the global birth and death rate of gene families, infer the most likely gene family size at all internal nodes, identify gene families that have accelerated rates of gain and loss (quantified by a p-value) and identify which branches cause the p-value to be small for significant families. AVAILABILITY: Software is available from http://www.bio.indiana.edu/~hahnlab/Software.html     

Family relationships of murine major histocompatibility complex class I genes. Sequence of the T2Aa pseudogene, a member of gene family 3

The major histocompatibility complex of the mouse contains numerous class I genes, most of which are encoded in the Qa and Tla regions. By hybridizations, the murine class I genes have been classified into three major families (Rogers, J. H. (1985a) Immunogenetics 21, 343-353). As yet, complete sequences are available only for members of family 1 (several H-2 and Qa genes) or family 2 (the pseudoallelic Tla genes T3b and T13c). We here present the complete nucleotide sequence of a gene from the Tla region that belongs to family 3. This gene, T2Aa, is a pseudogene by several criteria. The general structure of the gene is nonetheless well preserved. A comparison of the T2Aa sequence to those of other murine class I genes confirms the classification into three gene families. Members of gene families 2 and 3, located in the Tla region, are no more similar to each other than to family 1 (the H-2 and Qa2,3 genes). This suggests that families 2 and 3 were both created by ancient duplications of the functionally important family 1 genes. The fact that families 2 and 3 have diverged extensively both from family 1 and from each other may suggest that they are devoid of function.     

Recent origins of sperm genes in Drosophila

Newly created genes often acquire testis-specific or enhanced expression but neither the mechanisms responsible for this specificity nor the functional consequences of these evolutionary processes are well understood. Genomic analyses of the Drosophila melanogaster sperm proteome has identified 2 recently evolved gene families on the melanogaster lineage and 4 genes created by retrotransposition during the evolution of the melanogaster group that encode novel sperm components. The expanded Mst35B (protamine) and tektin gene families are the result of tandem duplication events with all family members displaying testis-specific expression. The Mst35B family encodes rapidly evolving protamines that display a robust signature of positive selection within the DNA-binding high-mobility group box consistent with functional diversification in genome repackaging during sperm nuclear remodeling. The Mst35B paralogs also reside in a significant regional cluster of testis-overexpressed genes. Tektins, known components of the axoneme, are encoded by 3 nearly identical X-linked genes, a finding consistent with very recent gene family expansion. In addition to localized duplication events, the evolution of the sperm proteome has also been driven by recent retrotransposition events resulting in Cdlc2, CG13340, Vha36, and CG4706. Cdlc2, CG13340, and Vha36 all display high levels of overexpression in the testis, and Cdlc2 and CG13340 reside within testis-overexpressed gene clusters. Thus, gene creation is a dynamic force in the evolution of sperm composition and possibly function, which further suggests that acquisition of molecular functionality in sperm may be an influential pathway in the fixation of new genes.     

The α-Actinin Gene Family: A Revised Classification

The sequencing of a genome is the first stage of its complete characterization. Subsequent work seeks to utilize available sequence data to gain a better understanding of the genes which are found within a genome. Gene families comprise large portions of the genomes of higher vertebrates, and the available genomic data allow for a reappraisal of gene family evolution. This reappraisal will clarify relatedness within and between gene families. One such family, the alpha-actinin gene family, is part of the spectrin superfamily. There are four known loci, which encode alpha-actinins 1, 2, 3, and 4. Of the eight domains in alpha-actinin, the actin-binding domain is the most highly conserved. Here we present evidence gained through phylogenetic analyses of the highly conserved actin-binding domain that alpha-actinin 2 was the first of the four alpha-actinins to arise by gene duplication, followed by the divergence of alpha-actinin 3 and then alpha-actinins 1 and 4. Resolution of the gene tree for this gene family has allowed us to reclassify several alpha-actinins which were previously given names inconsistent with the most widely accepted nomenclature for this gene family. This reclassification clarifies previous discrepancies in the public databases as well as in the literature, thus eliminating confusion caused by continued misclassification of members of the alpha-actinin gene family. In addition, the topology found for this gene family undermines the 2R hypothesis theory of two rounds of genome duplication early in vertebrate evolution.     

Natural history and functional divergence of protein tyrosine kinases

Cellular signaling is important for many biological processes including growth, differentiation, adhesion, motility and apoptosis. The protein tyrosine kinase (PTK) supergene family is the key mediator in cellular signaling in metazoans, directly associated with a variety of human diseases. All PTKs contain a highly conserved catalytic kinase domain, in spite of variable multi-domain structures. Within each PTK gene family, members exhibit functional divergence in substrate-specificity or temporal/tissue-specific expression, although their primary function is conserved. After conducting phylogenetic analysis on major PTK gene families, we found that the expanding of each PTK family was likely caused by gene or genome duplication event(s) that occurred before the emergence of teleosts but after the vertebrate-amphioxus split. We further investigated the evolutionary pattern of functional divergence after gene duplication in those gene families. Our results show that site-specific shifted evolutionary rate (altered functional constraint) is a common pattern in PTK gene family evolution.     

PaTre: a method for paralogy trees construction.

Genomes can be described as a collection of clusters, the gene families, whose members are called paralogs. Paralogs are genes that most probably share duplication history and show a significant similarity in their sequences, even if they perform slightly different biological function. Among the different mechanisms that have led to an increase of the genomic information during biological evolution, gene duplication is probably the most important. To better understand duplication events, the first step is to investigate the history of the gene families in order to detect which duplication events have taken place, and in which relative (partial) order. Here we present a method, called PaTre, that, given a gene family, attempts to construct the paralogy tree of the family. We will work under the hypothesis that every family member derives from a duplication process of another member. By the term paralogy tree, we mean a directed tree in which the root represents the most ancient paralog of the family and each oriented arc (a, b) represents the existence of a duplication event from the template gene a to its copy b. Notice that gene a survives the event and can serve as a template of more than one duplication event; in fact, there can be more than one arc leaving a. PaTre uses new algorithmic techniques motivated by the specific application at hand. The reliability of the inferential process has been tested by means of a simulator that implements different hypotheses on the duplication-with-modification paradigm and on three examples of different biological gene families, belonging either to lower and higher organisms.