Differential loss of ancestral gene families as a source of genomic divergence in animals

A phylogenetic approach was used to reconstruct the pattern of an apparent loss of 2106 ancestral gene families in four animal genomes (Caenorhabditis elegans, Drosophila melanogaster, human and fugu). Substantially higher rates of loss of ancestral gene families were found in the invertebrates than in the vertebrates. These results indicate that the differential loss of ancestral gene families can be a significant factor in the evolutionary diversification of organisms.     

Intron Dynamics and the Evolution of Integrin β-Subunit Genes: Maintenance of an Ancestral Gene Structure in the Coral, Acropora millepora

We have determined the genomic structure of an integrin β-subunit gene from the coral, Acropora millepora. The coding region of the gene contains 26 introns, spaced relatively uniformly, and this is significantly more than have been found in any integrin β-subunit genes from higher animals. Twenty-five of the 26 coral introns are also found in a β-subunit gene from at least one other phylum, indicating that the coral introns are ancestral. While there are some suggestions of intron gain or sliding, the predominant theme seen in the homologues from higher animals is extensive intron loss. The coral baseline allows one to infer that a number of introns found in only one phylum of higher animals result from frequent intron loss, as opposed to the seemingly more parsimonious alternative of isolated intron gain. The patterns of intron loss confirm results from protein sequences that most of the vertebrate genes, with the exception of β4, belong to one of two β subunit families. The similarity of the patterns within each of the β1,2,7 and β3,5,6,8 groups indicates that these gene structures have been very stable since early vertebrate evolution. Intron loss has been more extensive in the invertebrate genes, and obvious patterns have yet to emerge in this more limited data set. Received: 5 March 2001 / Accepted: 17 May 2001     

Protein Superfamily Evolution and the Last Universal Common Ancestor (LUCA)

By exploiting three-dimensional structure comparison, which is more sensitive than conventional sequence-based methods for detecting remote homology, we have identified a set of 140 ancestral protein domains using very restrictive criteria to minimize the potential error introduced by horizontal gene transfer. These domains are highly likely to have been present in the Last Universal Common Ancestor (LUCA) based on their universality in almost all of 114 completed prokaryotic (Bacteria and Archaea) and eukaryotic genomes. Functional analysis of these ancestral domains reveals a genetically complex LUCA with practically all the essential functional systems present in extant organisms, supporting the theory that life achieved its modern cellular status much before the main kingdom separation (Doolittle 2000). In addition, we have calculated different estimations of the genetic and functional versatility of all the superfamilies and functional groups in the prokaryote subsample. These estimations reveal that some ancestral superfamilies have been more versatile than others during evolution allowing more genetic and functional variation. Furthermore, the differences in genetic versatility between protein families are more attributable to their functional nature rather than the time that they have been evolving. These differences in tolerance to mutation suggest that some protein families have eroded their phylogenetic signal faster than others, hiding in many cases, their ancestral origin and suggesting that the calculation of 140 ancestral domains is probably an underestimate. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Rafael Zarobya]     

Evolution and diversity of axon guidance Robo receptor family genes

The diversity of axon guidance (AG) receptors reflects gains in complexity of the animal nervous system during evolution. Members of the Roundabout (Robo) family of receptors interact with Slit proteins and play important roles in many developmental processes, including AG and neural crest cell migration. There are four members of the Robo gene family. However, the evolutionary history of Robo family genes remain obscure. We analyzed the distribution of Robo family members in metazoan species ranging in complexity from hydras to humans. We undertook a phylogenetic analysis in metazoans, synteny analysis, and ancestral chromosome mapping in vertebrates, and detected selection pressure and functional divergence among four mammalian Robo paralogs. Based on our analysis, we proposed that the ancestral Robo gene could have undergone a tandem duplication in the vertebrate ancestor; then one round of whole genome duplication events occurred before the divergence of ancestral lamprey and gnathostome, generating four paralogs in early vertebrates. Robo4 paralog underwent segmental loss in the following evolutionary process. Our results showed that Robo3 paralog is under more powerful purifying selection pressure compared with other three paralogs, which could correlate with its unique expression pattern and function. Furthermore, we found four sites under positive selection pressure on the Ig1‐2 domains of Robo4 that might interfere with its binding to Slits ligand. Diverge analysis at the amino acid level showed that Robo4 paralog have relatively greater functional diversifications than other Robo paralogs. This coincides with the fact that Robo4 predominantly functions in vascular endothelial cells but not the nervous system.     

A novel gene encoding a ras-like GTP-binding protein from Trypanosoma brucei: an evolutionary ancestor of the ras and rap genes of higher eukaryotes?

The ras superfamily of GTP binding proteins encompasses a wide range of family members, related by conserved amino-acid motifs, and act as molecular binary switches that play key roles in cellular processes. Gene duplication and divergence has been postulated as the mechanism by which such family members have evolved their specific functions. We have cloned and sequenced a ras-like gene, tbrlp, from the primitive eukaryote Trypanosoma brucei. The gene encodes a protein of 227 amino acids and contains the six conserved subdomains that designate it as a ras/rap subfamily member. However, the presence of key diagnostic residues characteristic of both the ras and rap families of GTP confuse the familial classification of this gene. Phylogenetic analysis of the GTP binding domain places its origins at the divergence point of the ras/rap families and suggests that tbrlp is an ancestral gene to the ras/rap genes of higher eukaryotes.     

Mapping of the α4 subunit gene (GABRA4) to human chromosome 4 defines an α2—α4—β1—γ1 gene cluster: further evidence that modern GABAA receptor gene clusters are derived from an ancestral cluster

We demonstrated previously that an α₁—β₂—γ₂ gene cluster of the γ-aminobutyric acid (GABA_A) receptor is located on human chromosome 5q34–q35 and that an ancestral α—β—γ gene cluster probably spawned clusters on chromosomes 4, 5, and 15. Here, we report that the α₄ gene (GABRA4) maps to human chromosome 4p14–q12, defining a cluster comprising the α₂, α₄, β₁, and γ₁ genes. The existence of an α₂—α₄—β₁—γ₁ cluster on chromosome 4 and an α₁—α₆—β₂—γ₂ cluster on chromosome 5 provides further evidence that the number of ancestral GABA_A receptor subunit genes has been expanded by duplication within an ancestral gene cluster. Moreover, if duplication of the α gene occurred before duplication of the ancestral gene cluster, then a heretofore undiscovered subtype of α subunit should be located on human chromosome 15q11–q13 within an α₅—α_x—β₃—γ₃ gene cluster at the locus for Angelman and Prader—Willi syndromes.     

Origin and evolution of the Amyrel gene in the α-amylase multigene family of Diptera

Alpha-amylase genes often form multigene families in living organisms. In Diptera, a remote paralog, Amyrel, had been discovered in Drosophila, where this gene is currently used as a population and phylogenetic marker. The putative encoded protein has about 40% divergence with the classical amylases. We have searched the presence of the paralog in other families of Diptera to track its origin and understand its evolution. Amyrel was detected in a number of families of Muscomorpha (Brachycera-Cyclorrapha), suggesting an origin much older than previously thought. It has not been found elsewhere to date, and it is absent from the Anopheles gambiae genome. The intron–exon structures of the genes found so far suggest that the ancestral gene (before the duplication which gave rise to Amyrel) had two introns, and that subsequent, repeated and independent loss of one or both introns occurred in some Muscomorpha families. It seems that the Amyrel protein has experienced specific amino acid substitutions in regions generally well conserved in amylases, raising the possibility of peculiar, functional adaptations of this protein.     

Larger rearranged mitochondrial genomes in Dekkera/Brettanomyces yeasts are more closely related than smaller genomes with a conserved gene order

Summary Mitochondrial genomes from yeasts in the Dekkera/Brettanomyces/Eeniella group vary in size from 28 to 101 kb. Mapping of genes has shown that the three smallest genomes, of 28–42 kb, have the same gene order, whereas the three larger mitochondrial DNAs of 57–101 kb are rearranged relative to the smaller molecules and between themselves. To examine the relationships between these genomes, a phylogenetic tree has been constructed by sequence comparison of the mitochondrialencoded cytochrome oxidase subunit gene (COX2) from the six species. Contrary to expectation, the tree shows that the larger rearranged genomes are more closely related than the smaller mtDNAs. This result indicates that the gene order of the smaller mtDNAs (28–42 kb) is ancestral and that larger mtDNA molecules (57–101 kb) are more prone to rearrangement than smaller forms.Offprint requests to: G.D. Clark-Walker     

Fossil Humankind and Other Anthropoid Primates of China

More than 70 sites have yielded human fossils in China. They are attributed to Homo sapiens erectus and Homo sapiens sapiens. The earliest one is possibly about 1.7 Ma. A series of common morphological features, including shovel-shaped incisors and flatness of the face, characterize them. There is a morphological mosaic between H. s. erectus and H. s. sapiens in China. The existence of common features and the morphological mosaic suggest continuity of human evolution in China. That there are a few features which are more commonly seen in the Neanderthal lineage, occurring in a few Chinese fossil skulls, probably suggests gene flow between China and the West. Based on them, in 1998 I proposed an hypothesis—continuity with hybridization—for human evolution in China. The hypothesis is supported by paleolithic archeology, and it supports the multiregional evolution hypothesis of modern human origins. The anatomically modern humans of East Asia originated most probably in China. Although some nonhuman anthropoid primates of China—Gigantopithecus, Sivapithecus, Ramapithecus and Lufengpithecus—have been suggested as the direct ancestors of human beings, the discovery of more specimens and further studies do not support these suggestions. Therefore, it is most probable that the transition between apes and humans did not occur in China.     

Genealogy of the α-crystallin—small heat-shock protein superfamily

Sequences of 40 very diverse representatives of the α-crystallin–small heat-shock protein (α-Hsp) superfamily are compared. Their characteristic C-terminal ‘α-crystallin domain' of 80–100 residues contains short consensus sequences that are highly conserved from prokaryotes to eukaryotes. There are, in addition, some positions that clearly distinguish animal from non-animal α-Hsps. The α-crystallin domain is predicted to consist of two hydrophobic β-sheet motifs, separated by a hydrophilic region which is variable in length. Combination of a conserved α-crystallin domain with a variable N-terminal domain and C-terminal extension probably modulates the properties of the various α-Hsps as stress-protective and structural oligomeric proteins. Phylogeny reconstruction indicates that multiple α-Hsps were already present in the last common ancestor of pro- and eukaryotes. It is suggested that during eukaryote evolution, animal and non-animal α-Hsps originated from different ancestral gene copies. Repeated gene duplications gave rise to the multiple α-Hsps present in most organisms.     

Update on the Phylogenetic Systematics of New World Monkeys: Further DNA Evidence for Placing the Pygmy Marmoset (Cebuella) within the Genus Callithrix

We determined DNA sequences spanning the 1.8-kb long intron 1 of the interstitial retinol-binding protein nuclear gene (IRBP) for Callithrix geoffroyi, Callithrix humeralifer, and Callithrix argentata. With the 22 previously determined IRBP intron 1 sequences—21 from the 16 currently recognized genera of New World monkeys—the enlarged IRBP data represent for the marmoset genus Callithrix both its argentata and its jacchus species groups. Maximum-parsimony and neighbor-joining trees, constructed for the 25 aligned IRBP intron 1 sequences, support a provisional phylogenetic classification with three families: Atelidae, containing subfamily Atelinae; Pitheciidae, containing subfamily Pitheciinae; and Cebidae, containing subfamilies Cebinae, Aotinae, and Callitrichinae. In order to have taxa at the same hierarchical rank at equivalent age, this classification has all living callitrichines in a single tribe, Callitrichini, with four subtribes: Saguinina (Saguinus), Callimiconina (Callimico), Leontopithecina (Leontopithecus), and Callitrichina (Callithrix with the pygmy marmoset, Cebuella pygmaea, merged into it). The DNA evidence shows not only that Callithrix must include C. pygmaea to be monophyletic but also that the times of separation of pygmaea and the argentata and jacchus species groups from one another are to be expected (<5 Ma—million years ago) for species in a single genus. On relating the time course of the ceboid radiation to biogeographic information, it appears that in mid-Miocene times (10–11 Ma) a basal callitrichin stock branched into the ancestral population of Saguinus in one clade and the ancestral population of Leontopithecus and Callimico–Callithrix (or Leontopithecus–Callimico and Callithrix) in another clade. The proto-lion tamarins migrated south and eastward, where they were isolated in refugia, becoming the genus Leontopithecus. The stock remaining in Amazonia gave rise to present-day Callimico and Callithrix. The latter genus occupied a vast geographic area, giving rise to the argentata and pygmaea groups in Amazonia and to the jacchus group in central and eastern Brazil.     

Repeated horizontal gene transfer of GALactose metabolism genes violates Dollo’s law of irreversible loss

Dollo’s law posits that evolutionary losses are irreversible, thereby narrowing the potential paths of evolutionary change. While phenotypic reversals to ancestral states have been observed, little is known about their underlying genetic causes. The genomes of budding yeasts have been shaped by extensive reductive evolution, such as reduced genome sizes and the losses of metabolic capabilities. However, the extent and mechanisms of trait reacquisition after gene loss in yeasts have not been thoroughly studied. Here, through phylogenomic analyses, we reconstructed the evolutionary history of the yeast galactose utilization pathway and observed widespread and repeated losses of the ability to utilize galactose, which occurred concurrently with the losses of GALactose (GAL) utilization genes. Unexpectedly, we detected multiple galactose-utilizing lineages that were deeply embedded within clades that underwent ancient losses of galactose utilization. We show that at least two, and possibly three, lineages reacquired the GAL pathway via yeast-to-yeast horizontal gene transfer. Our results show how trait reacquisition can occur tens of millions of years after an initial loss via horizontal gene transfer from distant relatives. These findings demonstrate that the losses of complex traits and even whole pathways are not always evolutionary dead-ends, highlighting how reversals to ancestral states can occur.     

Location of structural genes for glucose phosphate isomerase and for leucyl aminopeptidase on chromosome VII of Petunia

Summary A gene termed gpiB, coding for one of the two isoenzyme zones of glucose phosphate isomerase in Petunia, has been mapped to a locus on chromosome VII by means of linkage to the marker An4, and by an allelic dosage effect on enzyme activity in trisomics. The high degree of linkage of electrophoretic alleles of gpiB to the pollen colour allele pair An4/an4, as demonstrated in the ancestral species, P. axillaris s.l. and P. integrifolia s.l., has been conserved in all cultivars of P. hybrida investigated. Another gene, coding for the enzyme leucyl-aminopeptidase could also be mapped to chromosome VII and the gene order An4 — lapB — gpiB determined. Apparently, distribution of lapB alleles is not related to the hybrid descent of P. hybrida.     

Paralogous origin of the rhodopsinlike opsin genes in lizards

Rhodopsinlike opsins constitute a distinct phylogenetic group (Yokoyama 1994, Mol. Biol. Evol. 11:32–39). This RH2 group includes the green-sensitive opsins in chicken and goldfish and the blue-sensitive opsin in a nocturnal lizard gecko. In the present study, we isolated and sequenced the genomic DNA clones for the RH2 opsin gene, rh2 _Ac , of the diurnal lizard Anolis carolinensis. This single-copy gene spans 18.3 kb from start to stop codons, making it the longest opsin gene known in vertebrates. Phylogenetic analysis strongly suggests that rh2 _Ac is more closely related to the chicken green opsin gene than to the gecko blue opsin gene. This gene tree differs from the organismal tree, where the two lizard species should be most closely related, implying that rh2 _Ac and the gecko blue-sensitive opsin genes have been derived from duplicate ancestral genes.Correspondence to: S. Yukiyama     

The complete mitochondrial genomes of sixteen ardeid birds revealing the evolutionary process of the gene rearrangements

Background

The animal mitochondrial genome is generally considered to be under selection for both compactness and gene order conservation. As more mitochondrial genomes are sequenced, mitochondrial duplications and gene rearrangements have been frequently identified among diverse animal groups. Although several mechanisms of gene rearrangement have been proposed thus far, more observational evidence from major taxa is needed to validate specific mechanisms. In the current study, the complete mitochondrial DNA of sixteen bird species from the family Ardeidae was sequenced and the evolution of mitochondrial gene rearrangements was investigated. The mitochondrial genomes were then used to review the phylogenies of these ardeid birds.

Results

The complete mitochondrial genome sequences of the sixteen ardeid birds exhibited four distinct mitochondrial gene orders in which two of them, named as “duplicate tRNA^Glu–CR” and “duplicate tRNA^Thr–tRNA^Pro and CR”, were newly discovered. These gene rearrangements arose from an evolutionary process consistent with the tandem duplication - random loss model (TDRL). Additionally, duplications in these gene orders were near identical in nucleotide sequences within each individual, suggesting that they evolved in concert. Phylogenetic analyses of the sixteen ardeid species supported the idea that Ardea ibis, Ardea modesta and Ardea intermedia should be classified as genus Ardea, and Ixobrychus flavicollis as genus Ixobrychus, and indicated that within the subfamily Ardeinae, Nycticorax nycticorax is closely related to genus Egretta and that Ardeola bacchus and Butorides striatus are closely related to the genus Ardea.

Conclusions

The duplicate tRNAThr–CR gene order is found in most ardeid lineages, suggesting this gene order is the ancestral pattern within these birds and persisted in most lineages via concerted evolution. In two independent lineages, when the concerted evolution stopped in some subsections due to the accumulation of numerous substitutions and deletions, the duplicate tRNAThr–CR gene order was transformed into three other gene orders. The phylogenetic trees produced from concatenated rRNA and protein coding genes have high support values in most nodes, indicating that the mitochondrial genome sequences are promising markers for resolving the phylogenetic issues of ardeid birds when more taxa are added.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-573) contains supplementary material, which is available to authorized users.     

Functional stability of the aristaless gene in appendage tip formation during evolution

The appendages of an insect are subdivided into distinct segments or podomeres. Many genes responsible for the regionalization of the growing limb into subdomains have been isolated from Drosophila. So far, only one gene is known in the leg that is solely required for specifying the distal-most pattern element—the pretarsal claw. In Drosophila, the gene aristaless is expressed in the centre of the antennal and leg imaginal disc that represents the most distal position of appendages, and in a proximal region. When Drosophila aristaless function is impaired, antennae and legs develop without their distal-most structures—the arista and the claw. We describe here the analysis of aristaless in the beetle Tribolium—an insect that shows a different, more ancestral mode of appendage formation than Drosophila. In Tribolium, appendages grow out continuously during embryogenesis, and no imaginal discs are formed. Tribolium aristaless (Tc-al) expression starts midway during appendage elongation, and is seen in a distal and a proximal position of head and trunk appendages. At the end of embryogenesis, Tc-al is seen in four expression domains in the leg, in the dorsal epidermis, and ventrally in every segment in lateral groups of cells, presumably the histoblasts. Like in the Drosophila adult, Tc-al is required in the larva for the formation of the most distal structures of the leg and the antenna as revealed by RNAi experiments. We conclude that aristaless is evolutionarily robust, meaning that it has retained its expressional and functional characteristics, although a heterochronic change of the process of appendage elongation took place towards the evolution of the highly derived diptera.Edited by D. Tautz     

Insights into early extracellular matrix evolution: spongin short chain collagen-related proteins are homologous to basement membrane type IV collagens and form a novel family widely distributed in invertebrates

Collagens are thought to represent one of the most important molecular innovations in the metazoan line. Basement membrane type IV collagen is present in all Eumetazoa and was found in Homoscleromorpha, a sponge group with a well-organized epithelium, which may represent the first stage of tissue differentiation during animal evolution. In contrast, spongin seems to be a demosponge-specific collagenous protein, which can totally substitute an inorganic skeleton, such as in the well-known bath sponge. In the freshwater sponge Ephydatia mülleri, we previously characterized a family of short-chain collagens that are likely to be main components of spongins. Using a combination of sequence- and structure-based methods, we present evidence of remote homology between the carboxyl-terminal noncollagenous NC1 domain of spongin short-chain collagens and type IV collagen. Unexpectedly, spongin short-chain collagen-related proteins were retrieved in nonsponge animals, suggesting that a family related to spongin constitutes an evolutionary sister to the type IV collagen family. Formation of the ancestral NC1 domain and divergence of the spongin short-chain collagen-related and type IV collagen families may have occurred before the parazoan-eumetazoan split, the earliest divergence among extant animal phyla. Molecular phylogenetics based on NC1 domain sequences suggest distinct evolutionary histories for spongin short-chain collagen-related and type IV collagen families that include spongin short-chain collagen-related gene loss in the ancestors of Ecdyzosoa and of vertebrates. The fact that a majority of invertebrates encodes spongin short-chain collagen-related proteins raises the important question to the possible function of its members. Considering the importance of collagens for animal structure and substratum attachment, both families may have played crucial roles in animal diversification.     

Rapid evolutionary dynamics in a 2.8‐Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi‐gene families, the DNA sequence of a 2.8‐Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance‐like gene families, was analyzed in the diploid grass Aegilops tauschii, the D‐genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non‐syntenic genes and only 13 ancestral genes are shared among these grass species. These non‐syntenic genes, including prolamin and resistance‐like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non‐syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin‐related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance‐like genes and subsequent differential expansion of each R gene family. The high frequency of non‐syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8‐Mb region with nine‐fold higher than the genome‐wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.