Search

Gene duplication is a key biological process in the evolutionary history of plants and an important driving force for the diversification of genomic and genetic systems. Interactions between the calcium sensor calcineurin B-like protein (CBL) and its target, CBL-interacting protein kinase (CIPK), play important roles in the plant's response to various environmental stresses. As a food crop with important economic and research value, turnip ( Brassica rapa var. rapa ) has been well adapted to the environment of the Tibetan Plateau and become a traditional crop in the region. The BrrCIPK9 gene in turnip has not been characterized. In this study, two duplicated genes, BrrCIPK9.1 and BrrCIPK9 . 2 , were screened from the turnip genome. Based on the phylogenetic analysis, BrrCIPK9.1 and BrrCIPK9 . 2 were found located in different sub-branches on the phylogenetic tree. Real-time fluorescence quantitative PCR analyses revealed their differential expression levels between the leaves and roots and in response to various stress treatments. The differences in their interactions with BrrCBLs were also revealed by yeast two-hybrid analyses. The results indicate that BrrCIPK9.1 and BrrCIPK9 . 2 have undergone Asparagine-alanine-phenylalanine (NAF) site divergence during turnip evolution, which has resulted in functional differences between them. Furthermore, BrrCIPK9.1 responded to high-pH (pH 8.5) stress, while BrrCIPK9.2 retained its ancestral function (low K + ), thus providing further evidence of their functional divergence. These functional divergence genes facilitate turnip's good adaptation to the extreme environment of the Tibetan Plateau. In summary, the results of this study reveal the characteristics of the duplicated BrrCIPK9 genes and provide a basis for further functional studies of BrrCBLs-BrrCIPKs in turnip.

On human (Homo sapiens) chromosome 1, there is a tandem duplication encompassing Rh blood group genes (Hosa_RHD and Hosa_RHCE). This duplication occurred in the common ancestor of humans, chimpanzees (Pan troglodytes), and gorillas, after splitting from their common ancestor with orangutans. Although several studies have been conducted on ape Rh blood group genes, the clear genome structures of the gene clusters remain unknown. Here, we determined the genome structure of the gene cluster of chimpanzee Rh genes by sequencing five BAC (Bacterial Artificial Chromosome) clones derived from chimpanzees. We characterized three complete loci (Patr_RHα, Patr_RHβ, and Patr_RHγ). In the Patr_RHβ locus, a short version of the gene, which lacked the middle part containing exons 4-8, was observed. The Patr_RHα and Patr_RHβ genes were located on the locations corresponding to Hosa_RHD and Hosa_RHCE, respectively, and Patr_RHγ was in the immediate vicinity of Patr_RHβ. Sequence comparisons revealed high sequence similarity between Patr_RHβ and Hosa_RHCE, while the chimpanzee Rh gene closest to Hosa_RHD was not Patr_RHα but rather Patr_RHγ. The results suggest that rearrangements and gene conversions frequently occurred between these genes and that the classic orthology/paralogy dichotomy no longer holds between human and chimpanzee Rh blood group genes., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Gene duplication plays an important role in the origin of evolutionary novelties, but the mechanisms responsible for the retention and functional divergence of the duplicated copy are not fully understood. The α-globin genes provide an example of a gene family with different numbers of gene duplicates among rodents. Whereas Rattus and Peromyscus each have three adult α-globin genes (HBA-T1, HBA-T2 and HBA-T3), Mus has only two copies. High rates of amino acid evolution in the independently derived HBA-T3 genes of Peromyscus and Rattus have been attributed to positive selection. Using RACE PCR, reverse transcription-PCR (RT-PCR) and RNA-seq, we show that another rodent, the bank vole Clethrionomys glareolus, possesses three transcriptionally active α-globin genes. The bank vole HBA-T3 gene is distinguished from each HBA-T1 and HBA-T2 by 20 amino acids and is transcribed 23- and 4-fold lower than HBA-T1 and HBA-T2, respectively. Polypeptides corresponding to all three genes are detected by electrophoresis, demonstrating that the translated products of HBA-T3 are present in adult erythrocytes. Patterns of codon substitution and the presence of low-frequency null alleles suggest a postduplication relaxation of purifying selection on bank vole HBA-T3.

Recent studies demonstrated that sequence divergence in both transcriptional regulatory region and coding region contributes to the subfunctionalization of duplicate gene. However, whether sequence divergence in the 3'-untranslated region (3'-UTR) has an impact on the subfunctionalization of duplicate genes remains unclear. Here, we identified two diverging duplicate vsx1 (visual system homeobox-1) loci in goldfish, named vsx1A1 and vsx1A2. Phylogenetic analysis suggests that vsx1A1 and vsx1A2 may arise from a duplication of vsx1 after the separation of goldfish and zebrafish. Sequence comparison revealed that divergence in both transcriptional and translational regulatory regions is higher than divergence in the introns. vsx1A2 expresses during blastula and gastrula stages and in adult retina but silences from segmentation stage to hatching stage, vsx1A1 starts expression from segmentation onward. Comparing to that zebrafish vsx1 expresses in all the developmental stages and in the adult retina, it appears that goldfish vsx1A1 and vsx1A2 are under going to share the functions of ancestral vsx1. The different but overlapping temporal expression patterns of vsx1A1 and vsx1A2 suggest that sequence divergence in the promoter region of duplicate vsx1 is not sufficient for partitioning the functions of ancestral vsx1. By comparing vsx1A1 and vsx1A2 3'-UTR-linked green fluorescent protein gene expression patterns, we demonstrated that the 3'-UTR of vsx1A1 remains but the 3'-UTR of vsx1A2 has lost the capability of mediating bipolar cell specific expression during retina development. These results indicate that sequence divergence in the 3'-UTRs has a clear effect on subfunctionalization of the duplicate genes., (© 2012 WILEY PERIODICALS, INC.)

We report the genome sequence of melon, an important horticultural crop worldwide. We assembled 375 Mb of the double-haploid line DHL92, representing 83.3% of the estimated melon genome. We predicted 27,427 protein-coding genes, which we analyzed by reconstructing 22,218 phylogenetic trees, allowing mapping of the orthology and paralogy relationships of sequenced plant genomes. We observed the absence of recent whole-genome duplications in the melon lineage since the ancient eudicot triplication, and our data suggest that transposon amplification may in part explain the increased size of the melon genome compared with the close relative cucumber. A low number of nucleotide-binding site-leucine-rich repeat disease resistance genes were annotated, suggesting the existence of specific defense mechanisms in this species. The DHL92 genome was compared with that of its parental lines allowing the quantification of sequence variability in the species. The use of the genome sequence in future investigations will facilitate the understanding of evolution of cucurbits and the improvement of breeding strategies.

The coelacanth Latimeria menadoensis, a living fossil, occupies a key phylogenetic position to explore the changes that have affected the genomes of the aquatic vertebrates that colonized dry land. This is the first study to isolate and analyze L. menadoensis mRNA. Three different vitellogenin transcripts were identified and their inferred amino acid sequences compared to those of other known vertebrates. The phylogenetic data suggest that the evolutionary history of this gene family in coelacanths was characterized by a different duplication event than those which occurred in teleosts, amniotes, and amphibia. Comparison of the three sequences highlighted differences in functional sites. Moreover, despite the presence of conserved sites compared with the other oviparous vertebrates, some sites were seen to have changed, others to be similar only to those of teleosts, and others still to resemble only to those of tetrapods., (© 2012 WILEY PERIODICALS, INC.)

One of the main motivations to study amphioxus is its potential for understanding the last common ancestor of chordates, which notably gave rise to the vertebrates. An important feature in this respect is the slow evolutionary rate that seems to have characterized the cephalochordate lineage, making amphioxus an interesting proxy for the chordate ancestor, as well as a key lineage to include in comparative studies. Whereas slow evolution was first noticed at the phenotypic level, it has also been described at the genomic level. Here, we examine whether the amphioxus genome is indeed a good proxy for the genome of the chordate ancestor, with a focus on protein-coding genes. We investigate genome features, such as synteny, gene duplication and gene loss, and contrast the amphioxus genome with those of other deuterostomes that are used in comparative studies, such as Ciona, Oikopleura and urchin.

Background: Phospholipase A (PLA) is an important group of enzymes responsible for phospholipid hydrolysis in lipid signaling. PLAs have been implicated in abiotic stress signaling and developmental events in various plants species. Genome-wide analysis of PLA superfamily has been carried out in dicot plant Arabidopsis. A comprehensive genome-wide analysis of PLAs has not been presented yet in crop plant rice., Methodology/principal Findings: A comprehensive bioinformatics analysis identified a total of 31 PLA encoding genes in the rice genome, which are divided into three classes; phospholipase A(1) (PLA(1)), patatin like phospholipases (pPLA) and low molecular weight secretory phospholipase A(2) (sPLA(2)) based on their sequences and phylogeny. A subset of 10 rice PLAs exhibited chromosomal duplication, emphasizing the role of duplication in the expansion of this gene family in rice. Microarray expression profiling revealed a number of PLA members expressing differentially and significantly under abiotic stresses and reproductive development. Comparative expression analysis with Arabidopsis PLAs revealed a high degree of functional conservation between the orthologs in two plant species, which also indicated the vital role of PLAs in stress signaling and plant development across different plant species. Moreover, sub-cellular localization of a few candidates suggests their differential localization and functional role in the lipid signaling., Conclusion/significance: The comprehensive analysis and expression profiling would provide a critical platform for the functional characterization of the candidate PLA genes in crop plants.

Background: Amino acid transporters (AATs) that transport amino acids across cellular membranes are essential for plant growth and development. To date, a genome-wide overview of the AAT gene family in rice is not yet available., Methodology/principal Findings: In this study, a total of 85 AAT genes were identified in rice genome and were classified into eleven distinct subfamilies based upon their sequence composition and phylogenetic relationship. A large number of OsAAT genes were expanded via gene duplication, 23 and 24 OsAAT genes were tandemly and segmentally duplicated, respectively. Comprehensive analyses were performed to investigate the expression profiles of OsAAT genes in various stages of vegetative and reproductive development by using data from EST, Microarrays, MPSS and Real-time PCR. Many OsAAT genes exhibited abundant and tissue-specific expression patterns. Moreover, 21 OsAAT genes were found to be differentially expressed under the treatments of abiotic stresses. Comparative analysis indicates that 26 AAT genes with close evolutionary relationships between rice and Arabidopsis exhibited similar expression patterns., Conclusions/significance: This study will facilitate further studies on OsAAT family and provide useful clues for functional validation of OsAATs.

C-type lectins play important roles in glycoprotein metabolism, multicellular integration and immunity. Based on their overall domain structure, they can be classified as different groups which possess different physiological functions. In this study, two novel simple C-type lectins were identified from grass carp (Ctenopharyngodon idellus), an important cultured fish in China. GcCL1 and gcCL2 share an essentially identical gene structure, a conserved promoter region shorter than 300 bp and an amino acid identity of 81.2%. Phylogenetic analysis indicated they may be products of gene duplication and could be classified as a new clade of group VII C-type lectins. Both of them were expressed in the eleven tissues examined, with the spleen having the highest abundance of transcript. The gcCL1 transcript was more abundant than gcCL2 in the majority of tissue samples from 2-yr-old grass carps, and was lower than those of gcCL2 before 15 days post-hatching. The expression of both genes was significantly up-regulated in spleen, muscle, skin, gills and hepatopancreas after induction by Aeromonas hydrophila. This is the first report that the expression of group VII C-type lectins could be induced by a pathogen, and indicates these lectins may be involved in the immune response to bacteria in fish., (Crown Copyright © 2011. Published by Elsevier Ltd. All rights reserved.)

Legumes host their Rhizobium spp. symbiont in novel root organs called nodules. Nodules originate from differentiated root cortical cells that dedifferentiate and subsequently form nodule primordia, a process controlled by cytokinin. A whole-genome duplication has occurred at the root of the legume Papilionoideae subfamily. We hypothesize that gene pairs originating from this duplication event and are conserved in distinct Papilionoideae lineages have evolved symbiotic functions. A phylogenetic strategy was applied to search for such gene pairs to identify novel regulators of nodulation, using the cytokinin phosphorelay pathway as a test case. In this way, two paralogous type-A cytokinin response regulators were identified that are involved in root nodule symbiosis. Response Regulator9 (MtRR9) and MtRR11 in medicago (Medicago truncatula) and an ortholog in lotus (Lotus japonicus) are rapidly induced upon Rhizobium spp. Nod factor signaling. Constitutive expression of MtRR9 results in arrested primordia that have emerged from cortical, endodermal, and pericycle cells. In legumes, lateral root primordia are not exclusively formed from pericycle cells but also require the involvement of the root cortical cell layer. Therefore, the MtRR9-induced foci of cell divisions show a strong resemblance to lateral root primordia, suggesting an ancestral function of MtRR9 in this process. Together, these findings provide a proof of principle for the applied phylogenetic strategy to identify genes with a symbiotic function in legumes.

In insects, some of the most ecologically important symbioses are nutritional symbioses that provide hosts with novel traits and thereby facilitate exploitation of otherwise inaccessible niches. One such symbiosis is the ancient obligate intracellular symbiosis of aphids with the γ-proteobacteria, Buchnera aphidicola. Although the nutritional basis of the aphid/Buchnera symbiosis is well understood, the processes and structures that mediate the intimate interactions of symbiotic partners remain uncharacterized. Here, using a de novo approach, we characterize the complement of 40 amino acid polyamine organocation (APC) superfamily member amino acid transporters (AATs) encoded in the genome of the pea aphid, Acyrthosiphon pisum. We find that the A. pisum APC superfamily is characterized by extensive gene duplications such that A. pisum has more APC superfamily transporters than other fully sequenced insects, including a ten paralog aphid-specific expansion of the APC transporter slimfast. Detailed expression analysis of 17 transporters selected on the basis of their phylogenetic relationship to five AATs identified in an earlier bacteriocyte expressed sequence tag study distinguished a subset of eight transporters that have been recruited for amino acid transport in bacteriocyte cells at the symbiotic interface. These eight transporters include transporters that are highly expressed and/or highly enriched in bacteriocytes and intriguingly, the four AATs that show bacteriocyte-enriched expression are all members of gene family expansions, whereas three of the four that are highly expressed but not enriched in bacteriocytes retain one-to-one orthology with transporters in other genomes. Finally, analysis of evolutionary rates within the large A. pisum slimfast expansion demonstrated increased rates of molecular evolution coinciding with two major shifts in expression: 1) a loss of gut expression and possibly a gain of bacteriocyte expression and 2) loss of expression in all surveyed tissues in asexual females. Taken together, our characterization of nutrient AATs at the aphid/Buchnera symbiotic interface provides the first examination of the processes and structures operating at the interface of an obligate intracellular insect nutritional symbiosis, offering unique insight into the types of genomic change that likely facilitated evolutionary maintenance of the symbiosis.

Background: A major goal of molecular evolutionary biology is to understand the fate and consequences of duplicated genes. In this context, aphids are intriguing because the newly sequenced pea aphid genome harbors an extraordinary number of lineage-specific gene duplications relative to other insect genomes. Though many of their duplicated genes may be involved in their complex life cycle, duplications in nutrient amino acid transporters appear to be associated rather with their essential amino acid poor diet and the intracellular symbiosis aphids rely on to compensate for dietary deficits. Past work has shown that some duplicated amino acid transporters are highly expressed in the specialized cells housing the symbionts, including a paralog of an aphid-specific expansion homologous to the Drosophila gene slimfast. Previous data provide evidence that these bacteriocyte-expressed transporters mediate amino acid exchange between aphids and their symbionts., Results: We report that some nutrient amino acid transporters show male-biased expression. Male-biased expression characterizes three paralogs in the aphid-specific slimfast expansion, and the male-biased expression is conserved across two aphid species for at least two paralogs. One of the male-biased paralogs has additionally experienced an accelerated rate of non-synonymous substitutions., Conclusions: This is the first study to document male-biased slimfast expression. Our data suggest that the male-biased aphid slimfast paralogs diverged from their ancestral function to fill a functional role in males. Furthermore, our results provide evidence that members of the slimfast expansion are maintained in the aphid genome not only for the previously hypothesized role in mediating amino acid exchange between the symbiotic partners, but also for sex-specific roles.

Background: Duplications are very common in the evolution of plant genomes, explaining the high number of members in plant gene families. New genes born after duplication can undergo pseudogenization, neofunctionalization or subfunctionalization. Rice is a model for functional genomics research, an important crop for human nutrition and a target for biofortification. Increased zinc and iron content in the rice grain could be achieved by manipulation of metal transporters. Here, we describe the ZINC-INDUCED FACILITATOR-LIKE (ZIFL) gene family in plants, and characterize the genomic structure and expression of rice paralogs, which are highly affected by segmental duplication., Results: Sequences of sixty-eight ZIFL genes, from nine plant species, were comparatively analyzed. Although related to MSF_1 proteins, ZIFL protein sequences consistently grouped separately. Specific ZIFL sequence signatures were identified. Monocots harbor a larger number of ZIFL genes in their genomes than dicots, probably a result of a lineage-specific expansion. The rice ZIFL paralogs were named OsZIFL1 to OsZIFL13 and characterized. The genomic organization of the rice ZIFL genes seems to be highly influenced by segmental and tandem duplications and concerted evolution, as rice genome contains five highly similar ZIFL gene pairs. Most rice ZIFL promoters are enriched for the core sequence of the Fe-deficiency-related box IDE1. Gene expression analyses of different plant organs, growth stages and treatments, both from our qPCR data and from microarray databases, revealed that the duplicated ZIFL gene pairs are mostly co-expressed. Transcripts of OsZIFL4, OsZIFL5, OsZIFL7, and OsZIFL12 accumulate in response to Zn-excess and Fe-deficiency in roots, two stresses with partially overlapping responses., Conclusions: We suggest that ZIFL genes have different evolutionary histories in monocot and dicot lineages. In rice, concerted evolution affected ZIFL duplicated genes, possibly maintaining similar expression patterns between pairs. The enrichment for IDE1 boxes in rice ZIFL gene promoters suggests a role in Zn-excess and Fe-deficiency up-regulation of ZIFL transcripts. Moreover, this is the first description of the ZIFL gene family in plants and the basis for functional studies on this family, which may play important roles in Zn and Fe homeostasis in plants.

Aldehyde dehydrogenases (ALDHs) catabolize toxic aldehydes and process the vitamin A-derived retinaldehyde into retinoic acid (RA), a small diffusible molecule and a pivotal chordate morphogen. In this study, we combine phylogenetic, structural, genomic, and developmental gene expression analyses to examine the evolutionary origins of ALDH substrate preference. Structural modeling reveals that processing of small aldehydes, such as acetaldehyde, by ALDH2, versus large aldehydes, including retinaldehyde, by ALDH1A is associated with small versus large substrate entry channels (SECs), respectively. Moreover, we show that metazoan ALDH1s and ALDH2s are members of a single ALDH1/2 clade and that during evolution, eukaryote ALDH1/2s often switched between large and small SECs after gene duplication, transforming constricted channels into wide opened ones and vice versa. Ancestral sequence reconstructions suggest that during the evolutionary emergence of RA signaling, the ancestral, narrow-channeled metazoan ALDH1/2 gave rise to large ALDH1 channels capable of accommodating bulky aldehydes, such as retinaldehyde, supporting the view that retinoid-dependent signaling arose from ancestral cellular detoxification mechanisms. Our analyses also indicate that, on a more restricted evolutionary scale, ALDH1 duplicates from invertebrate chordates (amphioxus and ascidian tunicates) underwent switches to smaller and narrower SECs. When combined with alterations in gene expression, these switches led to neofunctionalization from ALDH1-like roles in embryonic patterning to systemic, ALDH2-like roles, suggesting functional shifts from signaling to detoxification.

AMP-deaminase (AMPD, EC 3.5.4.6), which catalyzes the irreversible hydrolytic deamination of AMP to IMP and ammonia, is an important energy-related enzyme. The partial genomic sequence of the gene encoding myoadenylate deaminase (AMPD1) from the teleost fish Platichthys flesus was determined. The amino acid sequence of P. flesus AMPD1 shows 82% homology with that of the teleost fish Danio rerio. Comparison of genomic sequences of P. flesus and Rattus norvegicus reveals a high degree of conservation of both sequence and structural organization. A phylogenetic analysis of AMPD sequences shows that bony fish and mammalian AMPD1s arise by duplication of a common primordial gene.

Background: Flavonoid 3',5'-hydroxylases (F3'5'Hs) and flavonoid 3'-hydroxylases (F3'Hs) competitively control the synthesis of delphinidin and cyanidin, the precursors of blue and red anthocyanins. In most plants, F3'5'H genes are present in low-copy number, but in grapevine they are highly redundant., Results: The first increase in F3'5'H copy number occurred in the progenitor of the eudicot clade at the time of the γ triplication. Further proliferation of F3'5'Hs has occurred in one of the paleologous loci after the separation of Vitaceae from other eurosids, giving rise to 15 paralogues within 650 kb. Twelve reside in 9 tandem blocks of ~35-55 kb that share 91-99% identity. The second paleologous F3'5'H has been maintained as an orphan gene in grapevines, and lacks orthologues in other plants. Duplicate F3'5'Hs have spatially and temporally partitioned expression profiles in grapevine. The orphan F3'5'H copy is highly expressed in vegetative organs. More recent duplicate F3'5'Hs are predominately expressed in berry skins. They differ only slightly in the coding region, but are distinguished in the structure of the promoter. Differences in cis-regulatory sequences of promoter regions are paralleled by temporal specialisation of gene transcription during fruit ripening. Variation in anthocyanin profiles consistently reflects changes in the F3'5'H mRNA pool across different cultivars. More F3'5'H copies are expressed at high levels in grapevine varieties with 93-94% of 3'5'-OH anthocyanins. In grapevines depleted in 3'5'-OH anthocyanins (15-45%), fewer F3'5'H copies are transcribed, and at lower levels. Conversely, only two copies of the gene encoding the competing F3'H enzyme are present in the grape genome; one copy is expressed in both vegetative and reproductive organs at comparable levels among cultivars, while the other is transcriptionally silent., Conclusions: These results suggest that expansion and subfunctionalisation of F3'5'Hs have increased the complexity and diversification of the fruit colour phenotype among red grape varieties.

Conserved interactions among proteins or other molecules can provide strong evidence for coevolution across their evolutionary history. Diverse phylogenetic methods have been applied to identify potential coevolutionary relationships. In most cases, these methods minimally require comparisons of orthologous sequences and appropriate controls to separate effects of selection from the overall evolutionary relationships. In vertebrates, androgen receptor (AR) and cytochrome p450 aromatase (CYP19) share an affinity for androgenic steroids, which serve as receptor ligands and enzyme substrates. In a recent study, Tiwary and Li (Tiwary BK, Li W-H. 2009. Parallel evolution between aromatase and androgen receptor in the animal kingdom. Mol Biol Evol. 26:123-129) reported that AR and CYP19 displayed a signature of ancient and conserved interactions throughout all the Eumetazoa (i.e., cnidarians, protostomes, and deuterostomes). Because these findings conflicted with a number of previous studies, we reanalyzed the data set used by Tiwary and Li. First, our analyses demonstrate that the invertebrate genes used in the previous analysis are not orthologous sequences but instead represent a diverse set of nuclear receptors and CYP enzymes with no confirmed or hypothesized relationships with androgens. Second, we show that 1) their analytical approach, which measures correlations in evolutionary distances between proteins, potentially led to spurious significant relationships due simply to conserved domains and 2) control comparisons provide positive evidence for a strong influence of evolutionary history. We discuss how corrections to this method and analysis of key taxa (e.g., duplications in the teleost fish and suiform lineages) can inform investigations of the coevolutionary relationships between AR and aromatase.

The Pdx1 or Ipf1 gene encodes an important homeodomain-containing protein with key roles in pancreas development and function. Mutations in human PDX1 are implicated in developmental defects and disease of the pancreas. Extensive research, including genome sequencing, has indicated that Pdx1 is the only member of its gene family in mammals, birds, amphibians, and ray-finned fish, and with the exception of teleost fish, this gene forms part of the ParaHox gene cluster along with Gsx1 and Cdx2. The ParaHox cluster, however, is a remnant of a 4-fold genome duplication; the three other ParaHox paralogues lack a Pdx-like gene in all vertebrate genomes examined to date. We have used bacterial artificial chromosome cloning and synteny analysis to show that the ancestor of living jawed vertebrates in fact had more ParaHox genes, including two Pdx genes (Pdx1 and Pdx2). Surprisingly, the two Pdx genes have been retained in parallel in two quite distantly related lineages, the cartilaginous fish (sharks, skates, and chimeras) and the Indonesian coelacanth, Latimeria menadoensis. The Pdx2 gene has been lost independently in ray-finned fish and in tetrapods.

Gene duplication and neofunctionalization are known to be important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family crucial to the vertebrate immune system. In birds, two MHC class II β (MHCIIβ) exon 3 lineages have been recently characterized, and two hypotheses for the evolutionary history of MHCIIβ lineages were proposed. These lineages could have arisen either by 1) an ancient duplication and subsequent divergence of one paralog or by 2) recent parallel duplications followed by functional convergence. Here, we compiled a data set consisting of 63 MHCIIβ exon 3 sequences from six avian orders to distinguish between these hypotheses and to understand the role of selection in the divergent evolution of the two avian MHCIIβ lineages. Based on phylogenetic reconstructions and simulations, we show that a unique duplication event preceding the major avian radiations gave rise to two ancestral MHCIIβ lineages that were each likely lost once later during avian evolution. Maximum likelihood estimation shows that following the ancestral duplication, positive selection drove a radical shift from basic to acidic amino acid composition of a protein domain facing the α-chain in the MHCII α β-heterodimer. Structural analyses of the MHCII α β-heterodimer highlight that three of these residues are potentially involved in direct interactions with the α-chain, suggesting that the shift following duplication may have been accompanied by coevolution of the interacting α- and β-chains. These results provide new insights into the long-term evolutionary relationships among avian MHC genes and open interesting perspectives for comparative and population genomic studies of avian MHC evolution.

Based on the advances in the silkworm genome project, a new genome-wide analysis of cytochrome P450 genes was performed, focusing mainly on gene duplication. All four CYP9A subfamily members from the silkworm, Bombyx mori, were cloned by RT-PCR and designated CYP9A19-CYP9A22 by the P450 Nomenclature Committee. They each contain an open reading frame of 1,593 bp in length and encode a putative polypeptide of 531 amino acids. Both nucleic acid and amino acid sequences share very high identities with one another. The typical motifs of insect cytochrome P450, including the heme-binding region, helix-C, helix-I, helix-K, and PERF, show high sequence conservation among the multiple proteins. Alignment with their cDNA sequences revealed that these paralogues share identical gene structures, each comprising ten exons and nine introns of variable sizes. The locations of their introns (all nine introns follow the GT-AG rule) are absolutely conserved. CYP9A19, CYP9A20, and CYP9A21 form a tandem cluster on chromosome 17, whereas CYP9A22 is separated from the cluster by four tandem alcohol-dehydrogenase-like genes. Their phylogenetic relationships and structural comparisons indicated that these paralogues arose as the results of gene duplication events. RT-PCR detected their mRNAs in different "first line of defense" tissues, as well as in several other organs, suggesting diverse functions. Tissue-selective expression also indicates their functional divergence. The identified CYP9A genes have not yet been found outside the Lepidoptera, and are probably unique to the Lepidoptera. They show high sequence and structural similarities to each other, indicating that the Lepidoptera-specific P450s may be of functional importance. This analysis constitutes the first report of the clustering, spatial organization, and functional divergence of P450 in the silkworm.

Tbx5 is a T-box transcription factor that has been characterised in most vertebrate lineages and is widely expressed during the development of various embryonic structures, including the heart, the eyes and the anterior set of paired appendages (tetrapod forelimbs and fish pectoral fins). Mutations in TBX5 cause Holt-Oram syndrome, an autosomal dominant human "heart-hand" condition characterised by upper limb and heart malformations. In zebrafish, embryos with compromised tbx5 function show a complete absence of pectoral fins, whereas heart and eye development are not so highly disturbed. Here, we identify a new tbx5 gene in zebrafish that we have called tbx5b. This duplicate gene is present in all teleost genomes whose sequence is available, suggesting it resulted from the teleost-specific genome duplication event that took place during fish evolution. We show that tbx5b has lost the characteristic forelimb/pectoral fin expression of Tbx5 genes but has retained the eye and heart expression, partially overlapping with that of its paralogue, now referred to as tbx5a. Functional redundancy of tbx5a and tbx5b in the eye and heart would therefore explain the mild phenotypes observed during development of these organs in fish embryos with compromised tbx5a function.

Genes encoding for claudin-8 and -27 tight junction proteins in the euryhaline puffer fish (Tetraodon nigroviridis) were identified using its recently sequenced genome. Phylogenetic analysis indicated that multiple genes encoding for claudin-8 proteins (designated Tncldn8a, Tncldn8b, Tncldn8c and Tncldn8d) arose by tandem gene duplication. In contrast, both tandem and whole genome duplication events appear to have generated genes encoding for claudin-27 proteins (designated Tncldn27a, Tncldn27b, Tncldn27c and Tncldn27d). Tncldn8 and Tncldn27 mRNA were widely distributed in Tetraodon, suggesting involvement in various physiological processes. All Tncldn8 and Tncldn27 genes were expressed in gill and skin tissue (i.e., epithelia exposed directly to the external environment). A potential role for claudin-8 and -27 proteins in the regulation of hydromineral balance in Tetraodon was investigated by examining alterations in mRNA abundance in select ionoregulatory tissue of fish acclimated to freshwater (FW) and seawater (SW). In FW or SW, Tetraodon exhibited alterations in Na(+)-K(+)-ATPase activity (a correlate of transcellular transport) typical of a euryhaline teleost fish. Simultaneously, tissue and gene specific alterations in Tncldn8 and Tncldn27 transcript abundance occurred. These data provide some insight into the duplication history of cldn8 and cldn27 genes in fishes and suggest a possible role for claudin-8 and -27 proteins in the osmoregulatory strategies of euryhaline teleosts.

Follistatin (Fst) inhibits transforming growth factor-beta (TGF-B) proteins and is a known regulator of amniote myogenesis. Here, we used phylogenetic, genomic and experimental approaches to study its evolution in teleosts. Phylogenetic analyses suggested that one fst gene (fst1) is common to euteleosts, but a second gene (fst2) is conserved specifically within the Ostariophysi. Zebrafish fst1/2 respectively appear on chromosomes 5 and 10 in two genomic regions, each with conserved synteny to a single region in tetrapods. Interestingly, other teleosts have two corresponding chromosomal regions with a similar repertoire of paralogues. Phylogenetic reconstruction clustered these gene duplicates into two sister clades branching from tetrapod sequences. We suggest that an ancestral fst-containing chromosome was duplicated during the teleost whole genome duplication, but that fst2 was lost in lineages external to the Ostariophysi. We show that Fst1 of teleosts/mammals has evolved under strong purifying selection, but the N-terminal of Fst2 may have evolved under positive selection. Furthermore, the tissue-specific expression of zebrafish fst2 was restricted to fewer tissues compared to its paralogue and the single fst1 orthologue of Atlantic salmon (Salmo salar). Zebrafish fst1/2 may have subfunctionalized relative to non-duplicated vertebrate lineages, as several regions in the fst promoter of tetrapods were conserved with one paralogue, but not both. Finally, we examined the embryonic expression of fst1 in a teleost outside the Ostariophysi (Atlantic salmon). During segmentation, fst1 was expressed in the anterior somite compartment but was excluded from muscle progenitors that strongly expressed myogenic regulatory factors (MRFs). Later, fst1 was expressed in myogenic progenitors of the pectoral fin buds and also within the pax7(+) cell layer external to the myotome.

A goldfish (Carassius auratus auratus) bacterial artificial chromosome genomic library (BAC library) was constructed from one aquarium-bred male specimen (tetraploid, 4n=100, genome size=3.52 pg/cell). The library consists of 128,352 positive clones with an average insert size of 150.4 kb, covering the genome 11-fold. All clones were spotted onto nylon filters and thus are available for screening of genomic regions of interest, such as candidate genes, gene families, or large-sized syntenic DNA regions of cyprinid species. Preliminary screens with two genes were conducted with hybridizing probes to the genes RAG1 and lgi1. RAG1 is a single-copy gene in zebrafish and is duplicated in C. a. auratus. We found a very close correlation between the number of positive BAC clones and the expected library coverage. Two copies of lgi1 were found in zebrafish. We have detected four different copies in C. a. auratus, not in the expected abundance, which indicates some variation in the coverage of the BAC library. The preliminary screens indicate that many duplicated genes that resulted from the ancient fish-specific genome duplication persist in the tetraploid goldfish genome. Hence, the BAC library will provide a useful resource for the future work on comparative genomics, polyploidy, diploidization, and evolutionary genomics in fishes.

During phylogenetic analysis of the Nepenthaceae cpDNA trnK intron, it became apparent that a second non-functional copy of the locus was present in most of the investigated taxa. The translocation event was older than the radiation of all recent Nepenthaceae, and the translocated pseudogenized copy was conserved in nearly all members of the plant family. Using single chloroplast PCR and inverse PCR, we could exclude a plastom location for the second copy. Although translocation into the nucleus is possible, mitochondrial localization seems more likely based on these data. In total, the translocated sequence contained at least 3525 base pairs (bp) that were homologous to the Spinacia oleracea chloroplast genome. Comparative phylogenetic analysis of the non-functional copy revealed a high amount of homoplasies compared to topologies from the cpDNA trnK intron phylogenetic reconstruction. Therefore, this copy proved to be insufficient for phylogenetic reconstruction of the family. Since two different paralogs of the non-functional copy were found in one species, it is feasible that different paralogs were conserved in different groups and that paralogous sequences were included in the data matrix. These data demonstrate that phylogenetic analyses of pseudogenized copies of phylogenetically relevant loci should be performed with great caution. In addition, pseudogenized copies can exist in nearly every member of a plant family, and can be PCR-amplified at levels comparable to the specific copy. In this case, the inclusion of such copies can easily remain unnoticed, thus leading to faulty hypotheses.

Mutations in the leucine-rich, glioma-inactivated 1 gene, LGI1, cause autosomal-dominant lateral temporal lobe epilepsy via unknown mechanisms. LGI1 belongs to a subfamily of leucine-rich repeat genes comprising four members (LGI1-LGI4) in mammals. In this study, both comparative developmental as well as molecular evolutionary methods were applied to investigate the evolution of the LGI gene family and, subsequently, of the functional importance of its different gene members. Our phylogenetic studies suggest that LGI genes evolved early in the vertebrate lineage. Genetic and expression analyses of all five zebrafish lgi genes revealed duplications of lgi1 and lgi2, each resulting in two paralogous gene copies with mostly nonoverlapping expression patterns. Furthermore, all vertebrate LGI1 orthologs experience high levels of purifying selection that argue for an essential role of this gene in neural development or function. The approach of combining expression and selection data used here exemplarily demonstrates that in poorly characterized gene families a framework of evolutionary and expression analyses can identify those genes that are functionally most important and are therefore prime candidates for human disorders.

Partial and complete genome duplications occurred during evolution and resulted in the creation of new genes and gene families. We identified a novel and intricate human gene family located primarily in regions of segmental duplications on human chromosome 1. We named it NBPF, for neuroblastoma breakpoint family, because one of its members is disrupted by a chromosomal translocation in a neuroblastoma patient. The NBPF genes have a repetitive structure with high intragenic and intergenic sequence similarity in both coding and noncoding regions. These similarities might expose these genomic regions to illegitimate recombination, resulting in structural variation in the NBPF genes. The encoded proteins contain a highly conserved domain of unknown function, which we have named the NBPF repeat. In silico analysis combined with the isolation of multiple full-length cDNA clones showed that several members of this gene family are abundantly expressed in a large variety of tissues and cell lines. Strikingly, no discernable orthologues could be identified in the completed genomes of fruit fly, nematode, mouse, or rat, but sequences with low homology could be isolated from the draft canine and bovine genomes. Interestingly, this gene family shows primate-specific duplications that result in species-specific arrays of NBPF homologous sequences. Overall, this novel NBPF family reflects the continuous evolution of primate genomes that resulted in large physiological differences, and its potential role in this process is discussed.

In bacteria, magnesium uptake is mainly mediated by the well-characterized CorA type of membrane proteins. In recent years, functional homologues have been characterized in the inner mitochondrial membrane of yeast and mammals (the MRS2/LPE10 type), in the plasma membrane of yeast (the ALR/MNR type) and, as an extended family of proteins, in the model plant Arabidopsis thaliana. Despite generally low sequence similarity, individual proteins can functionally complement each other over large phylogenetic distances. All these proteins are characterized by a universally conserved Gly-Met-Asn (GMN) motif at the end of the first of two conserved transmembrane domains near the C-terminus. Mutations of the GMN motif are known to abolish Mg(2+) transport, but the naturally occurring variants GVN and GIN may be associated with the transport of other divalent cations, such as zinc and cadmium, respectively. We refer to this whole class of proteins as the 2-TM-GxN type. The functional membrane channel is thought to be formed by oligomers containing four or five subunits. The wealth of sequence data now available allows us to explore the evolutionary diversification of the basic 2-TM-GxN model within the so-called metal ion transporter (MIT) superfamily. Here we report phylogenetic analyses on more than 360 homologous protein sequences derived from genomic sequences from representatives of all three domains of life. Independent gene duplications have occurred in fungi, plants and proteobacteria at different phylogenetic depths. Moreover, there is ample evidence for several instances of horizontal gene transfer of members of the 2-TM-GxN superfamily in Eubacteria and Archaea. Only single genes of the MRS2 type have been identified in vertebrate genomes. In contrast, 15 members are found in the model plant Arabidopsis thaliana, which appear to have arisen by at least four independent founder events before the diversification of flowering plants. Phylogenetic clade assignment seems to correlate with alterations in the highly conserved sequence around the GMN motif. This presumably forms an integral part of the pore surface, and changes in its structure may result in altered transport capacities for different divalent cations.

Gene duplication is the main source of evolutionary novelties. However, the problem with duplicates is that the purifying selection overlooks deleterious mutations in the redundant sequence, which therefore, instead of gaining a new function, often degrades into a functionless pseudogene. This risk of functional loss instead of gain is much higher for small populations of higher organisms with a slow and complex development. We propose that it is the epigenetic tissue/stage-complementary silencing of duplicates that makes them exposable to the purifying selection, thus saving them from pseudogenization and opening the way towards new function(s). Our genome-wide analyses of gene duplicates in several eukaryotic species combined with the phylogenetic comparison of vertebrate alpha- and beta-globin gene clusters strongly support this epigenetic complementation (EC) model. The distinctive condition for a new duplicate to survive by the EC mechanism seems to be its repositioning to an ectopic site, which is accompanied by changes in the rate and direction of mutagenesis. The most distinguished in this respect is the human genome. In this review, we extend and discuss the data on the EC- and repositioning-dependent fate of gene duplicates with the special emphasis on the problem of detecting brief postduplication period of adaptive evolution driven by positive selection. Accordingly, we propose a new CpG-focused measure of selection that is insensitive to translocation-caused biases in mutagenesis.

In the 14 years since the first higher-level bird phylogenies based on DNA sequence data, avian phylogenetics has witnessed the advent and maturation of the genomics era, the completion of the chicken genome and a suite of technologies that promise to add considerably to the agenda of avian phylogenetics. In this review, we summarize current approaches and data characteristics of recent higher-level bird studies and suggest a number of as yet untested molecular and analytical approaches for the unfolding tree of life for birds. A variety of comparative genomics strategies, including adoption of objective quality scores for sequence data, analysis of contiguous DNA sequences provided by large-insert genomic libraries, and the systematic use of retroposon insertions and other rare genomic changes all promise an integrated phylogenetics that is solidly grounded in genome evolution. The avian genome is an excellent testing ground for such approaches because of the more balanced representation of single-copy and repetitive DNA regions than in mammals. Although comparative genomics has a number of obvious uses in avian phylogenetics, its application to large numbers of taxa poses a number of methodological and infrastructural challenges, and can be greatly facilitated by a 'community genomics' approach in which the modest sequencing throughputs of single PI laboratories are pooled to produce larger, complementary datasets. Although the polymerase chain reaction era of avian phylogenetics is far from complete, the comparative genomics era-with its ability to vastly increase the number and type of molecular characters and to provide a genomic context for these characters-will usher in a host of new perspectives and opportunities for integrating genome evolution and avian phylogenetics.

In this study the molecular evolution of duplicated HoxA genes in zebrafish and fugu has been investigated. All 18 duplicated HoxA genes studied have a higher non-synonymous substitution rate than the corresponding genes in either bichir or paddlefish, where these genes are not duplicated. The higher rate of evolution is not due solely to a higher non-synonymous-to-synonymous rate ratio but to an increase in both the non-synonymous as well as the synonymous substitution rate. The synonymous rate increase can be explained by a change in base composition, codon usage, or mutation rate. We found no changes in nucleotide composition or codon bias. Thus, we suggest that the HoxA genes may experience an increased mutation rate following cluster duplication. In the non-Hox nuclear gene RAG1 only an increase in non-synonymous substitutions could be detected, suggesting that the increased mutation rate is specific to duplicated Hox clusters and might be related to the structural instability of Hox clusters following duplication. The divergence among paralog genes tends to be asymmetric, with one paralog diverging faster than the other. In fugu, all b-paralogs diverge faster than the a-paralogs, while in zebrafish Hoxa-13a diverges faster. This asymmetry corresponds to the asymmetry in the divergence rate of conserved non-coding sequences, i.e., putative cis-regulatory elements. These results suggest that the 5' HoxA genes in the same cluster belong to a co-evolutionary unit in which genes have a tendency to diverge together.

DNA sequencing studies of the granule-bound starch synthase gene (GBSSI) indicate the presence of two loci in Viburnum. Gene trees from separate and combined phylogenetic analyses of the GBSSI paralogues are generally congruent with each other and with trees from previous analyses, especially those of Donoghue et al. [Syst. Bot. 29 (2004) 188] based on nuclear ribosomal ITS and chloroplast trnK intron DNA sequences. Specifically, our GBSSI trees confirm (i) the monophyly of some and non-monophyly of other traditionally recognized taxonomic sections, (ii) the presence of three major supra-sectional lineages within Viburnum, and (iii) the resolution of many species relationships within the section-level clades. Analyses of GBSSI also provide greater resolution of relationships within the largest supra-sectional lineage. Relationships at the base of the Viburnum phylogeny remain uncertain; in particular, the position of the root, relationships among the supra-sectional clades, and the exact placement of several smaller groups (e.g., Viburnum clemensiae, Viburnum urceolatum, and section Pseudotinus). In two lineages each GBSSI paralogue is represented by two distinct sequences. The presence of additional copies appears to be correlated with polyploidy in these clades. Placement of the homoeologues in our gene trees suggests the possibility of a hybrid origin for these polyploids.

Two, apparently functional, paralogues of the RPB2 gene, which encodes the second largest subunit of RNA polymerase II, are shown to be present in two major groups of asterid plants. Although all other land plants surveyed so far have been found to have only one of these two copies, the RPB2 gene phylogeny inferred from the 3' half of the gene for 35 angiosperm taxa and six other land plants indicates that the duplication of the RPB2 gene occurred earlier than the time for origin of the asterid group, probably near the origin of "core eudicots." The d copy is present in all plants which are unambiguously assigned to the core eudicots, whereas the I copy is retained only in the lamiid clade, Ericales, and Escallonia, all belonging to the asterid group of plants. Both parsimony and likelihood analyses of sequences from the 3' half of the gene give strong bootstrap support for these conclusions. There is no support for monophyly of the taxa having both copies. Thus, numerous losses of one of the copies must be inferred. Structurally, both paralogues appear functional, and transcription is demonstrated for both copies. In the lamiid group, the d copy has lost introns 18-23. The well supported phylogenetic relationships implied by the RPB2 gene phylogeny are largely congruent with well supported phylogenetic hypotheses based on other sequence data. However, Ilex, usually assigned to the campanuliid clade, is instead supported as being a member of the lamiid clade, both from sequence data and the presence of an I copy as well as the loss of introns 18-23 in the d copy. Escallonia, supported as a member of the campanuliid clade both by RPB2-d-sequences and previously published DNA sequence data, has all the introns 18-23 in its d copy, as do all other members studied from the campanuliid group. We used the Markov Chain Monte Carlo (MCMC) approach of the MrBayes program to implement Maximum Likelihood bootstrapping. Under the same model of evolution, bootstrapping frequencies are significantly lower than the Bayesian posterior probabilities inferred from the MCMC chain.

In the first molecular study of ostracod (Crustacea) vision, we present partial cDNA sequences of ostracod visual pigment genes (opsins). We found strong support for differential expression of opsins in ostracod median and compound eyes and suggest that photoreceptor specific expression may be a general phenomenon in organisms with multiple receptors. We infer that eye-specific expression predates the divergence of the two species examined, Skogsbergia lerneri and Vargula hilgendorfii, because eye-specific opsin orthologs are present in both species. We found multiple opsin loci in ostracods, estimating that at least eight are present in Skogsbergia lerneri. All opsins from both ostracod species examined are more closely related to each other than to any other known opsin sequences. Because we find no evidence for gene conversion or alternative splicing, we suggest the occurrence of many recent gene duplications. Why ostracods may have retained multiple recent opsin gene duplicates is unknown, but we discuss several possible hypotheses.

Two new genome sequences confirm that a whole genome duplication occurred in an ancestor of Saccharomyces cerevisiae. This left a legacy of about 500 pairs of duplicated genes, many of which contribute to this yeast's ability to ferment glucose anaerobically; a few have been evolving so quickly they retain almost no sequence similarity to each other.

Many gene families in mammals have members that are expressed more or less uniquely in the retina or differentially in specific retinal cell types. We describe here analyses of nine such gene families with regard to phylogenetic relationships and chromosomal location. The families are opsins, G proteins (alpha, beta, and gamma subunits), phosphodiesterases type 6, cyclic nucleotide-gated channels, G-protein-coupled receptor kinases, arrestins, and recoverins. The results suggest that multiple new gene copies arose in all of these families very early in vertebrate evolution during a period with extensive gene duplications. Many of the new genes arose through duplications of large chromosome regions (blocks of genes) or even entire chromosomes, as shown by linkage with other gene families. Some of the phototransduction families belong to the same duplicated regions and were thus duplicated simultaneously. We conclude that gene duplications in early vertebrate evolution probably helped facilitate the specialization of the retina and the subspecialization of different retinal cell types.

The Nogo-66 receptor NgR has been implicated in the mediation of inhibitory effects of central nervous system (CNS) myelin on axon growth in the adult mammalian CNS. NgR binds to several myelin-associated ligands (Nogo-66, myelin associated glycoprotein, and oligodendrocyte-myelin glycoprotein), which, among other inhibitory proteins, impair axonal regeneration in the CNS of adult mammals. In contrast to mammals, severed axons readily regenerate in the fish CNS. Nevertheless, fish axons are repelled by mammalian oligodendrocytes in vitro. Therefore, the identification of fish NgR homologs is a crucial step towards understanding NgR functions in vertebrate systems competent of CNS regeneration. Here, we report the discovery of four zebrafish (Danio rerio) and five fugu (Takifugu rubripes) NgR homologs. Synteny between fish and human, comparable intron-exon structures, and phylogenetic analyses provide convincing evidence that the true fish orthologs were identified. The topology of the phylogenetic trees shows that the extra fish genes were produced by duplication events that occurred in ray-finned fishes before the divergence of the zebrafish and pufferfish lineages. Expression of zebrafish NgR homologs was detected relatively early in development and prominently in the adult brain, suggesting functions in axon growth, guidance, or plasticity.

One of the most popular sequences for phylogenetic inference at the generic and infrageneric levels in plants is the internal transcribed spacer (ITS) region of the 18S-5.8S-26S nuclear ribosomal cistron. The prominence of this source of nuclear DNA sequence data is underscored by a survey of phylogenetic publications involving comparisons at the genus level or below, which reveals that of 244 papers published over the last five years, 66% included ITS sequence data. Perhaps even more striking is the fact that 34% of all published phylogenetic hypothesis have been based exclusively on ITS sequences. Notwithstanding the many important contributions of ITS sequence data to phylogenetic understanding and knowledge of genome relationships, a number of molecular genetic processes impact ITS sequences in ways that may mislead phylogenetic inference. These molecular genetic processes are reviewed here, drawing attention to both underlying mechanism and phylogenetic implications. Among the most prevalent complications for phylogenetic inference is the existence in many plant genomes of extensive sequence variation, arising from ancient or recent array duplication events, genomic harboring of pseudogenes in various states of decay, and/or incomplete intra- or inter-array homogenization. These phenomena separately and collectively create a network of paralogous sequence relationships potentially confounding accurate phylogenetic reconstruction. Homoplasy is shown to be higher in ITS than in other DNA sequence data sets, most likely because of orthology/paralogy conflation, compensatory base changes, problems in alignment due to indel accumulation, sequencing errors, or some combination of these phenomena. Despite the near-universal usage of ITS sequence data in plant phylogenetic studies, its complex and unpredictable evolutionary behavior reduce its utility for phylogenetic analysis. It is suggested that more robust insights are likely to emerge from the use of single-copy or low-copy nuclear genes.

The widely accepted notion that two whole-genome duplications occurred during early vertebrate evolution (the 2R hypothesis) stems from the fact that vertebrates often possess several genes corresponding to a single invertebrate homolog. However the number of genes predicted by the Human Genome Project is less than twice as many as in the Drosophila melanogaster or Caenorhabditis elegans genomes. This ratio could be explained by two rounds of genome duplication followed by extensive gene loss, by a single genome duplication, by sequential local duplications, or by a combination of any of the above. The traditional method used to distinguish between these possibilities is to reconstruct the phylogenetic relationships of vertebrate genes to their invertebrate orthologs; ratios of invertebrate-to-vertebrate counterparts are then used to infer the number of gene duplication events. The lancelet, amphioxus, is the closest living invertebrate relative of the vertebrates, and unlike protostomes such as flies or nematodes, is therefore the most appropriate outgroup for understanding the genomic composition of the last common ancestor of all vertebrates. We analyzed the relationships of all available amphioxus genes to their vertebrate homologs. In most cases, one to three vertebrate genes are orthologous to each amphioxus gene (median number=2). Clearly this result, and those of previous studies using this approach, cannot distinguish between alternative scenarios of chordate genome expansion. We conclude that phylogenetic analyses alone will never be sufficient to determine whether genome duplication(s) occurred during early chordate evolution, and argue that a "phylogenomic" approach, which compares paralogous clusters of linked genes from complete amphioxus and human genome sequences, will be required if the pattern and process of early chordate genome evolution is ever to be reconstructed., (Copyright 2003 Wiley-Liss, Inc.)

MADS-box proteins are a large family of transcription factors. In plants, many genes belonging to this family are involved in the homeosis of the floral system. Up to now, they have mainly been studied in angiosperms, especially in the model species Arabidopsis thaliana and Antirrhinum majus. We undertook a study of MADS-box genes in Ginkgo biloba, the unique extant representative of a whole branch of the phylogenetic tree of the seed plants. A polymerase chain reaction (PCR) survey reveals the diversity of MADS-box genes present in the genome of the Ginkgo. Duplications probably occurred specifically in the ginkgophyte lineage. Phylogenetic analyses revealed that one of these genes, GBM5, is an orthologue of the AGAMOUS gene of A. thaliana. We cloned and sequenced the entire cDNA of the GBM5 gene and studied its intron/exon structure. We showed by reverse transcriptase-PCR that it is expressed in both floral and vegetative tissues. We discuss the molecular evolution of the AGAMOUS family of genes.

Glutamine synthetase (GSase) is a key enzyme in nitrogen metabolism and encoded by a single gene in mammals. Using PCR cloning techniques, including RT-PCR from total RNA and PCR from a cDNA library, we find evidence of four expressed GSase mRNAs for the tetraploid rainbow trout. For two of these mRNAs (Onmy-GS01, -GS02) we characterize the full-length coding regions, and for two others (Onmy-GS03, -GS04), we describe partial sequences. Northern analysis of Onmy-GS01, -GS02, -GS03 and -GS04 indicates that (1) Onmy-GS02 is expressed at higher levels relative to the other transcripts in most adult tissues, with the exception of brain and gill, where Onmy-GS01 is at the highest level, and (2) the tissue with the highest level of expression of all four transcripts is the brain, with decreasing levels in the intestine, liver, red muscle, gill/kidney, white muscle and heart. Clearly, rainbow trout possess multiple GSase genes with differing levels of tissue expression, implying manifold potential routes of regulation for this octameric enzyme. Our data also indicate that caution should be taken when interpreting mRNA expression data of a single gene, unless multiple genes have been ruled out. Consistent with a southern blot, phylogenetic and intron sequence analyses imply that the trout genes are encoded by at least four separate loci, belonging to two distinct evolutionary branches. Our data on rainbow trout, together with those from two full-length zebrafish Danio rerio GSase genes compiled from GenBank ESTs, support the idea that fish GSases are polyphyletic and that gene duplications have occurred at multiple points and in independent lineages throughout the evolution of bony fishes.

Background: When analyzing protein sequences using sequence similarity searches, orthologous sequences (that diverged by speciation) are more reliable predictors of a new protein's function than paralogous sequences (that diverged by gene duplication). The utility of phylogenetic information in high-throughput genome annotation ("phylogenomics") is widely recognized, but existing approaches are either manual or not explicitly based on phylogenetic trees., Results: Here we present RIO (Resampled Inference of Orthologs), a procedure for automated phylogenomics using explicit phylogenetic inference. RIO analyses are performed over bootstrap resampled phylogenetic trees to estimate the reliability of orthology assignments. We also introduce supplementary concepts that are helpful for functional inference. RIO has been implemented as Perl pipeline connecting several C and Java programs. It is available at http://www.genetics.wustl.edu/eddy/forester/. A web server is at http://www.rio.wustl.edu/. RIO was tested on the Arabidopsis thaliana and Caenorhabditis elegans proteomes., Conclusion: The RIO procedure is particularly useful for the automated detection of first representatives of novel protein subfamilies. We also describe how some orthologies can be misleading for functional inference.

Low-copy nuclear genes in plants are a rich source of phylogenetic information. They hold a great potential to improve the robustness of phylogenetic reconstruction at all taxonomic levels, especially where universal markers such as cpDNA and nrDNA are unable to generate strong phylogenetic hypotheses. Low-copy nuclear genes, however, remain underused in plant phylogenetic studies due to practical and theoretical complications in unraveling the evolutionary dynamics of nuclear gene families. The lack of the universal markers or universal PCR primers of low-copy nuclear genes has also hampered their phylogenetic utility. It has recently become clear that low-copy nuclear genes are particularly helpful in resolving close interspecific relationships and in reconstructing allopolyploidization in plants. Gene markers that are widely, if not universally, useful have begun to emerge. Although utilizing low-copy nuclear genes usually requires extra lab work such as designing PCR primers, PCR-cloning, and/or Southern blotting, rapid accumulation of gene sequences in the databases and advances in cloning techniques have continued to make such studies more feasible. With the growing number of theoretical studies devoted to the gene tree and species tree problem, a solid foundation for reconstructing complex plant phylogenies based on multiple gene trees began to build. It is also realized increasingly that fast evolving introns of the low-copy nuclear genes will provide much needed phylogenetic information around the species boundary and allow us to address fundamental questions concerning processes of plant speciation. Phylogenetic and molecular evolutionary analyses of developmentally important genes will add a new dimension to systematic and evolutionary studies of plant diversity.

Glycine receptors mediating synaptic inhibition are heteromeric proteins constituted of alpha and beta subunits. The mammalian GlyR subunits constitute a subgroup in the superfamily of ligand-gated ionic channels. To compare the evolutionary events in the mammalian and teleostean lineages for the receptor family, we first undertook systematic cloning of the constitutive subunits of the zebrafish glycine receptor. The isolation of two alpha subunits (alphaZ1 and alphaZ2) and one beta subunit (betaZ) has been reported previously and we report here the characterization of two novel alpha subunits, alphaZ3 and alphaZ4, increasing the known zebrafish subunits number to four alpha and one beta. Establishment of phylogenetic relationships reveals that alphaZ1, alphaZ3 and betaZeta are orthologous to mammalian alpha1, alpha3 and beta subunits. However, two zebrafish GlyRalpha subunit genes are orthologous to the unique avian and mammalian alpha4 subunit revealing a duplication of the alpha4 gene in zebrafish. Whole-mount in situ hybridization in 24-hours post fertilization (hpf) and 52-hpf embryos of the daughter gene products display very different expression patterns indicating distinct functions of the duplicated genes. Gene mapping reveals that the two duplicated genes are localized on two different linkage groups (LG5 and LG22) as would be daughter genes resulting from a large-scale duplication of the ancestral genome. Finally, we report that a linked pair of genes on human chromosome 4 (alpha3 and beta) is also linked on linkage group 1 in zebrafish (alphaZ3 and betaZ) as a consequence of a mosaic conserved syntheny.

Pyruvate kinase (PK; EC 2.7.1.40) of Thermoproteus tenax was purified to homogeneity, and its coding gene was cloned and expressed in Escherichia coli. It represents a homomeric tetramer with a molecular mass of 49 kDa per subunit. PK exhibits positive binding cooperativity with respect to phosphoenolpyruvate and metal ions such as Mg(2+) and Mn(2+). Heterotropic effects, as commonly found for PKs from bacterial and eucaryal sources, could not be detected. The enzyme does not depend on K(+) ions. Heterotrophically grown cells exhibit specific activity of PK four times higher than autotrophically grown cells. Since the mRNA level of the PK coding gene is also accordingly higher in heterotrophic cells, we conclude that the PK activity is adjusted to growth conditions mainly on the transcript level. The enzymic properties of the PK and the regulation of its expression are discussed with respect to the physiological framework given by the T. tenax-specific variant of the Embden-Meyerhof-Parnas pathway. T. tenax PK shows moderate overall sequence similarity (25 to 40% identity) to its bacterial and eucaryal pendants. Phylogenetic analyses of the known PK sequences result in a dichotomic tree topology that divides the enzymes into two major PK clusters, probably diverged by an early gene duplication event. The phylogenetic divergence is paralleled by a striking phenotypic differentiation of PKs: PKs of cluster I, which occur in eucaryal cytoplasm, some gamma proteobacteria, and low-GC gram-positive bacteria, are only active in the presence of fructose-1,6-bisphosphate or other phosphorylated sugars, whereas PKs of cluster II, found in various bacterial phyla, plastids, and in Archaea, show activity without effectors but are commonly regulated by the energy charge of the cell.

The phylogenetic position of hagfishes in vertebrate evolution is currently controversial. The 18S and 28S rRNA trees support the monophyly of hagfishes and lampreys. In contrast, the mitochondrial DNAs suggest the close association of lampreys and gnathostomes. To clarify this controversial issue, we have conducted cloning and sequencing of the four nuclear DNA-coded single-copy genes encoding the triose phosphate isomerase, calreticulin, and the largest subunit of RNA polymerase II and III. Based on these proteins, together with the Mn superoxide dismutase for which hagfish and lamprey sequences are available in database, phylogenetic trees have been inferred by the maximum likelihood (ML) method of protein phylogeny. It was shown that all the five proteins prefer the monophyletic tree of cyclostomes, and the total log-likelihood of the five proteins significantly supports the cyclostome monophyly at the level of +/-1 SE. The ML trees of aldolase family comprising three nonallelic isoforms and the complement component group comprising C3, C4, and C5, both of which diverged during vertebrate evolution by gene duplications, also suggest the cyclostome monophyly.