Fugu and human sequence comparison identifies novel human genes and conserved non-coding sequences

The compact genome of the pufferfish, Fugu rubripes, has been proposed as a 'reference' genome to aid in annotating and analysing the human genome. We have annotated and compared 85 kb of Fugu sequence containing 17 genes with its homologous loci in the human draft genome and identified three 'novel' human genes that were missed or incompletely predicted by the previous gene prediction methods. Two of the novel genes contain zinc finger domains and are designated ZNF366 and ZNF367. They map to human chromosomes 5q13.2 and 9q22.32, respectively. The third novel gene, designated C9orf21, maps to chromosome 9q22.32. This gene is unique to vertebrates, and the protein encoded by it does not contain any known domains. We could not find human homologs for two Fugu genes, a novel chemokine gene and a kinase gene. These genes are either specific to teleosts or lost in the human lineage. The Fugu-human comparison identified several conserved non-coding sequences in the promoter and intronic regions. These sequences, conserved during 450 million years of vertebrate evolution, are likely to be involved in gene regulation. The 85 kb Fugu locus is dispersed over four human loci, occupying about 1.5 Mb. Contiguity is conserved in the human genome between six out of 16 Fugu gene pairs. These contiguous chromosomal segments should share a common evolutionary history dating back to the common ancestor of mammals and teleosts. We propose contiguity as strong evidence to identify orthologous genes in distant organisms. This study confirms the utility of the Fugu as a supplementary tool to uncover and confirm novel genes and putative gene regulatory regions in the human genome.     

Analysis of human accelerated DNA regions using archaic hominin genomes

Several previous comparisons of the human genome with other primate and vertebrate genomes identified genomic regions that are highly conserved in vertebrate evolution but fast-evolving on the human lineage. These human accelerated regions (HARs) may be regions of past adaptive evolution in humans. Alternatively, they may be the result of non-adaptive processes, such as biased gene conversion. We captured and sequenced DNA from a collection of previously published HARs using DNA from an Iberian Neandertal. Combining these new data with shotgun sequence from the Neandertal and Denisova draft genomes, we determine at least one archaic hominin allele for 84% of all positions within HARs. We find that 8% of HAR substitutions are not observed in the archaic hominins and are thus recent in the sense that the derived allele had not come to fixation in the common ancestor of modern humans and archaic hominins. Further, we find that recent substitutions in HARs tend to have come to fixation faster than substitutions elsewhere in the genome and that substitutions in HARs tend to cluster in time, consistent with an episodic rather than a clock-like process underlying HAR evolution. Our catalog of sequence changes in HARs will help prioritize them for functional studies of genomic elements potentially responsible for modern human adaptations.     

Characterization of novel GPCR gene coding locus in amphioxus genome: gene structure, expression, and phylogenetic analysis with implications for its involvement in chemoreception

Chemosensation is the primary sensory modality in almost all metazoans. The vertebrate olfactory receptor genes exist as tandem clusters in the genome, so that identifying their evolutionary origin would be useful for understanding the expansion of the sensory world in relation to a large-scale genomic duplication event in a lineage leading to the vertebrates. In this study, I characterized a novel GPCR (G-protein-coupled receptor) gene-coding locus from the amphioxus genome. The genomic DNA contains an intronless ORF whose deduced amino acid sequence encodes a seven-transmembrane protein with some amino acid residues characteristic of vertebrate olfactory receptors (ORs). Surveying counterparts in the Ciona intestinalis (Asidiacea, Urochordata) genome by querying BLAST programs against the Ciona genomic DNA sequence database resulted in the identification of a remotely related gene. In situ hybridization analysis labeled primary sensory neurons in the rostral epithelium of amphioxus adults. Based on these findings, together with comparison of the developmental gene expression between amphioxus and vertebrates, I postulate that chemoreceptive primary sensory neurons in the rostrum are an ancient cell population traceable at least as far back in phylogeny as the common ancestor of amphioxus and vertebrates.     

Ancient duplications of the human proglucagon gene

The human proglucagon gene (GCG) is encoded within a finished 576-kb DNA sequence generated by the Human Genome Project. GCG is flanked by 18 kb and 65 kb of DNA, 5' and 3', respectively, that do not encode genes. The genomic sequence that includes GCG was found to have a long history of gene duplication events. Some members of the glucagon-like family of genes, GCG on chromosome 2 and GIP on chromosome 17, may be products of ancient genome duplications on the early vertebrate lineage. A large genomic tandem duplication event that included DPP4-like and GCG genes occurred before the amphibian-mammal divergence, but one of the duplicated copies of GCG has been lost on the human lineage. Recently, a processed pseudogene of the X-chromosome-linked gene TIMM8A was inserted downstream of GCG. Some ancient duplicates of GCG may retain physiological functions in other vertebrates.     

The human genes for calbindin 27 and 29 kDa proteins are located on chromosomes 8 and 16, respectively

Genomic clones coding for the brain calcium-binding protein, calbindin 29 kDa, were isolated from a human library. A fragment containing exon 2 was used as a probe to investigate the presence of the gene in human x rodent somatic cell hybrids. The gene was unambiguously assigned to chromosome 16. The closely-related calbindin 27 kDa gene was previously assigned to chromosome 8. These two genes, deriving from a common ancestor, thus appear to have been separated during vertebrate evolution.     

文昌鱼特异的基因倍增

进化生物学和发育生物学的结合产生了一门新兴学科——进化发育生物学,近年来该领域研究取得了丰硕的成果。头索动物文昌鱼是现存生物中最近似于脊椎动物直接祖先的生物,在与脊椎动物分化后形态改变很小,其基因组未曾经历大规模的基因组倍增,在一定程度上反映了脊椎动物祖先型基因组的特征,但在漫长的独立进化历程中基因组自身还是经历了一些变化。本文介绍了在几例在文昌鱼支系中独立发生的基因倍增事件(Hox; Evx; HNF-3; Calmodulin-like),有力地揭示了文昌鱼虽然与脊椎动物直接祖先极其接近,但其基因组有其自身特性,不能简单地将之等同于脊椎动物直接祖先。Abstract: The union of the two complementary disciplines, developmental biology and evolutionary biology resulted in a new division of evolutionary developmental biology, namely “Evo-Devo”. Recently, the research on this field has been fruitful in understanding the origin and development of vertebrates. The cephalochordate amphioxus, which remains in relatively invariant morphology since the divergence from the vertebrate lineage, is the closest living relative to vertebrates. The vertebrate-like simple body plan and preduplicative genome provide amphioxus genes the privilege to serve as key landmark to understand morphological evolution. However, the amphioxus genome has not escaped evolution. In this paper several examples of independent gene (Hox; Evx; HNF-3 and Calmodulin-like) duplications in the cephalochordate lineage were summarized. These particularities and oddities remind the fact that amphioxus is not an immediate ancestor of the vertebrates but ‘only’ the closest living relative to the ancestor, with a mix of prototypical and amphioxus-specific features in its genome.     

The expression of the dodecameric ferritin in Listeria spp. is induced by iron limitation and stationary growth phase

The new discipline of Evolutionary Developmental Biology (Evo-Devo) is facing the fascinating paradox of explaining morphological evolution using conserved pieces or genes to build divergent animals. The cephalochordate amphioxus is the closest living relative to the vertebrates, with a simple, chordate body plan, and a genome directly descended from the ancestor prior to the genome-wide duplications that occurred close to the origin of vertebrates. Amphioxus morphology may have remained relatively invariant since the divergence from the vertebrate lineage, but the amphioxus genome has not escaped evolution. We report the isolation of a second Emx gene (AmphiEmxB) arising from an independent duplication in the amphioxus genome. We also argue that a tandem duplication probably occurred in the Posterior part of the Hox cluster in amphioxus, giving rise to AmphiHox14, and discuss the structure of the chordate and vertebrate ancestral clusters. Also, a tandem duplication of Evx in the amphioxus lineage produced a prototypical Evx gene (AmphiEvxA) and a divergent gene (AmphiEvxB), no longer involved in typical Evx functions. These examples of specific gene duplications in amphioxus, and other previously reported duplications summarized here, emphasize the fact that amphioxus is not the ancestor of the vertebrates but 'only' the closest living relative to the ancestor, with a mix of prototypical and amphioxus-specific features in its genome.     

Why mice have lost genes for COL21A1, STK17A, GPR145 and AHRI: evidence for gene deletion at evolutionary breakpoints in the rodent lineage

The mouse genome has undergone extensive chromosome rearrangement relative to the human genome since these species last shared a common ancestor. One possible consequence of these rearrangements is the deletion of genes that are located within evolutionary breakpoint regions. In this article, we present evidence of four human genes (COL21A1, STK17A, GPR145 and ARHI) that are located in regions corresponding to evolutionary breakpoints in rodents and lack mouse and rat orthologues. We propose that "evolutionary breakpoint-associated gene deletion" is an unexpected consequence of evolutionary chromosome rearrangement, and we describe a novel mechanism through which genes can be lost during evolution.     

Linkage assignment of a DNA sequence (ERCC2L1) homologous to a human DNA repair gene in Xiphophorus fishes: implications for the evolutionary derivation of human chromosome 19.

Fish gene mapping studies have identified several syntenic groups showing conservation over more than 400 million years of vertebrate evolution. In particular, Xiphophorus linkage group IV has been identified as a homolog of human chromosomes 15 and 19. During mammalian evolution, loci coding for glucosephosphate isomerase, peptidase D, muscle creatine kinase, and several DNA repair genes (ERCC1, ERCC2, and XRCC1) appear as a conserved syntenic group on human chromosome 19. When X. clemenciae and X. milleri PstI endonuclease-digested genomic DNA was used in Southern analysis with a human ERCC2 DNA repair gene probe, a strongly cross-hybridizing restriction fragment length polymorphism was observed. Backcrosses to X. clemenciae from X. milleri x X. clemenciae F1 hybrids allowed tests for linkage of the ERCC2-like polymorphism to markers covering a large proportion of the genome. Statistically significant evidence for linkage was found only for ERCC2L1 and CKM (muscle creatine kinase), with a total of 41 parents and 2 recombinants (4.7% recombination, chi 2 = 35.37, P less than 0.001); no evidence for linkage to GPI and PEPD in linkage group IV was detected. The human chromosome 19 synteny of ERCC2 and CKM thus appears to be conserved in Xiphophorus, while other genes located nearby on human chromosome 19 are in a separate linkage group in this fish. If Xiphophorus gene arrangements prove to be primitive, human chromosome 19 may have arisen from chromosome fusion or translocation events at some point since divergence of mammals and fishes from a common ancestor.     

Gene space dynamics during the evolution of Aegilops tauschii, Brachypodium distachyon, Oryza sativa, and Sorghum bicolor genomes

Nine different regions totaling 9.7 Mb of the 4.02 Gb Aegilops tauschii genome were sequenced using the Sanger sequencing technology and compared with orthologous Brachypodium distachyon, Oryza sativa (rice), and Sorghum bicolor (sorghum) genomic sequences. The ancestral gene content in these regions was inferred and used to estimate gene deletion and gene duplication rates along each branch of the phylogenetic tree relating the four species. The total gene number in the extant Ae. tauschii genome was estimated to be 36,371. The gene deletion and gene duplication rates and total gene numbers in the four genomes were used to estimate the total gene number in each node of the phylogenetic tree. The common ancestor of the Brachypodieae and Triticeae lineages was estimated to have had 28,558 genes, and the common ancestor of the Panicoideae, Ehrhartoideae, and Pooideae subfamilies was estimated to have had 27,152 or 28,350 genes, depending on the ancestral gene scenario. Relative to the Brachypodieae and Triticeae common ancestor, the gene number was reduced in B. distachyon by 3,026 genes and increased in Ae. tauschii by 7,813 genes. The sum of gene deletion and gene duplication rates, which reflects the rate of gene synteny loss, was correlated with the rate of structural chromosome rearrangements and was highest in the Ae. tauschii lineage and lowest in the rice lineage. The high rate of gene space evolution in the Ae. tauschii lineage accounts for the fact that, contrary to the expectations, the level of synteny between the phylogenetically more related Ae. tauschii and B. distachyon genomes is similar to the level of synteny between the Ae. tauschii genome and the genomes of the less related rice and sorghum. The ratio of gene duplication to gene deletion rates in these four grass species closely parallels both the total number of genes in a species and the overall genome size. Because the overall genome size is to a large extent a function of the repeated sequence content in a genome, we suggest that the amount and activity of repeated sequences are important factors determining the number of genes in a genome.     

Cloning and gene map assignment of the Xiphophorus DNA ligase 1 gene

Fishes represent the stem vertebrate condition and have maintained several gene arrangements common to mammalian genomes throughout the 450 Myr of divergence from a common ancestor. One such syntenic arrangement includes the GPI-PEPD enzyme association on Xiphophorus linkage group IV and human chromosome 19. Previously we assigned the Xiphophorus homologue of the human ERCC2 gene to linkage group U5 in tight association with the CKM locus. CKM is also tightly linked to the ERCC2 locus on human chromosome 19, leading to speculation that human chromosome 19 may have arisen by fusion of two ancestral linkage groups which have been maintained in fishes. To investigate this hypothesis further, we isolated and sequenced Xiphophorus fish genomic regions exhibiting considerable sequence similarity to the human DNA ligase 1 amino acid sequence. Comparison of the fish DNA ligase sequence with those of other species suggests several modes of amino acid conservation in this gene. A 2.2-kb restriction fragment containing part of an X. maculatus DNA ligase 1 exon was used in backcross hybrid mapping with 12 enzyme or RFLP loci. Significant linkage was observed between the nucleoside phosphorylase (NP2) and the DNA ligase (LIG1) loci on Xiphophorus linkage group VI. This assignment suggests that the association of four DNA repair-related genes on human chromosome 19 may be the result of chance chromosomal rearrangements.      

Evolutionary implications of the human aldolase-A, -B, -C, and -pseudogene chromosome locations.

The aldolase genes represent an ancient gene family with tissue-specific isozymic forms expressed only in vertebrates. The chromosomal locations of the aldolase genes provide insight into their tissue-specific and developmentally regulated expression and evolution. DNA probes for the human aldolase-A and -C genes and for an aldolase pseudogene were used to quantify and map the aldolase loci in the haploid human genome. Genomic hybridization of restriction fragments determined that all the aldolase genes exist in single copy in the haploid human genome. Spot-blot analysis of sorted chromosomes mapped human aldolase A to chromosome 16, aldolase C to chromosome 17, the pseudogene to chromosome 10; it previously had mapped the aldolase-B gene to chromosome 9. All loci are unlinked and located on to two pairs of morphologically similar chromosomes, a situation consistent with tetraploidization during isozymic and vertebrate evolution. Sequence comparisons of expressed and flanking regions support this conclusion. These locations on similar chromosome pairs correctly predicted that the aldolase pseudogene arose when sequences from the aldolase-A gene were inserted into the homologous aldolase location on chromosome 10.     

Evidence for different origins of sex chromosomes in closely related Oryzias fishes: substitution of the master sex-determining gene

The medaka Oryzias latipes and its two sister species, O. curvinotus and O. luzonensis, possess an XX-XY sex-determination system. The medaka sex-determining gene DMY has been identified on the orthologous Y chromosome [O. latipes linkage group 1 (LG1)] of O. curvinotus. However, DMY has not been discovered in other Oryzias species. These results and molecular phylogeny suggest that DMY was generated recently [approximately 10 million years ago (MYA)] by gene duplication of DMRT1 in a common ancestor of O. latipes and O. curvinotus. We identified seven sex-linked markers from O. luzonensis (sister species of O. curvinotus) and constructed a sex-linkage map. Surprisingly, all seven sex-linked markers were located on an autosomal linkage group (LG12) of O. latipes. As suggested by the phylogenetic tree, the sex chromosomes of O. luzonensis should be "younger" than those of O. latipes. In the lineage leading to O. luzonensis after separation from O. curvinotus approximately 5 MYA, a novel sex-determining gene may have arisen and substituted for DMY. Oryzias species should provide a useful model for evolution of the master sex-determining gene and differentiation of sex chromosomes from autosomes.     

Mapping and conservation of the group-specific component gene in mouse

The group-specific component (GC), also known as the vitamin D-binding protein, transports vitamin D and its metabolites in plasma to target tissues throughout the body. The GC gene shares an evolutionary origin with genes encoding albumin (ALB) and alpha-fetoprotein (AFP). All three genes are descendants of an evolutionary ancestor that arose from an intragenic triplication. As a result, each gene is composed of three homologous domains. The study described here characterizes and compares mouse GC to the corresponding nucleotide and amino acid sequences of GC from human and rat. The deduced amino acid sequence of mouse GC was 78% identical to human and 91% identical to rat GC. The results suggest that, unlike the corresponding sequences in the ALB and AFP genes, chromosomal sequences encoding the first domain and the leader sequence of the GC gene have specifically been conserved throughout vertebrate evolution. Protection of domain I during evolution may correlate with an important functional aspect of its sequence. The mouse GC gene was mapped to chromosome 5, where the ALB and AFP genes are also located, demonstrating conservation of the three genes in vertebrate species.     

The "fish-specific" Hox cluster duplication is coincident with the origin of teleosts

The Hox gene complement of zebrafish, medaka, and fugu differs from that of other gnathostome vertebrates. These fishes have seven to eight Hox clusters compared to the four Hox clusters described in sarcopterygians and shark. The clusters in different teleost lineages are orthologous, implying that a "fish-specific" Hox cluster duplication has occurred in the stem lineage leading to the most recent common ancestor of zebrafish and fugu. The timing of this event, however, is unknown. To address this question, we sequenced four Hox genes from taxa representing basal actinopterygian and teleost lineages and compared them to known sequences from shark, coelacanth, zebrafish, and other teleosts. The resulting gene genealogies suggest that the fish-specific Hox cluster duplication occurred coincident with the origin of crown group teleosts. In addition, we obtained evidence for an independent Hox cluster duplication in the sturgeon lineage (Acipenseriformes). Finally, results from HoxA11 suggest that duplicated Hox genes have experienced diversifying selection immediately after the duplication event. Taken together, these results support the notion that the duplicated Hox genes of teleosts were causally relevant to adaptive evolution during the initial teleost radiation.     

Cytogenetics, conserved synteny and evolution of chicken fucosyltransferase genes compared to human

Fucosyltransferases appeared early in evolution, since they are present from bacteria to primates and the genes are well conserved. The aim of this work was to study these genes in the bird group, which is particularly attractive for the comprehension of the evolution of the vertebrate genome. Twelve fucosyltransferase genes have been identified in man. The orthologues of theses genes were looked for in the chicken genome and cytogenetically localized by FISH. Three families of fucosyltransferases: alpha6-fucosyltransferases, alpha3/4-fucosyltransferases, and protein-O-fucosyltransferases, were identified in the chicken with their associated genes. The alpha2-fucosyltransferase family, although present in some invertebrates and amphibians was not found in birds. This absence, also observed in Drosophila, may correspond to a loss of these genes by negative selection. Of the eight chicken genes assigned, six fell on chromosome segments where conservation of synteny between human and chicken was already described. For the two remaining loci, FUT9 and FUT3/5/6, the location may correspond to a new small syntenic area or to an insertion. FUT4 and FUT3/5/6 were found on the same chicken chromosome. These results suggest a duplication of an ancestral gene, initially present on the same chromosome before separation during evolution. By extension, the results are in favour of a common ancestor for the alpha3-fucosyltransferase and the alpha4-fucosyltransferase activities. These observations suggest a general mechanism for the evolution of fucosyltransferase genes in vertebrates by duplication followed by divergent evolution.     

Immunoglobulin light chains in medaka (Oryzias latipes)

The gene segments encoding antibodies have been studied in many capacities and represent some of the best-characterized gene families in traditional animal disease models (mice and humans). To date, multiple immunoglobulin light chain (IgL) isotypes have been found in vertebrates and it is unclear as to which isotypes might be more primordial in nature. Sequence data emerging from an array of fish genome projects is a valuable resource for discerning complex multigene assemblages in this critical branch point of vertebrate phylogeny. Herein, we have analyzed the genomic organization of medaka (Oryzias latipes) IgL gene segments based on recently released genome data. The medaka IgL locus located on chromosome 11 contains at least three clusters of IgL gene segments comprised of multiple gene assemblages of the kappa light chain isotype. These data suggest that medaka IgL gene segments may undergo both intra- and inter-cluster rearrangements as a means to generate additional diversity. Alignments of expressed sequence tags to concordant gene segments which revealed each of the three IgL clusters are expressed. Collectively, these data provide a genomic framework for IgL genes in medaka and indicate that Ig diversity in this species is achieved from at least three distinct chromosomal regions.     

Genomic structure and paralogous regions of the inversion breakpoint occurring between human chromosome 3p12.3 and orangutan chromosome 2

Intrachromosomal duplications play a significant role in human genome pathology and evolution. To better understand the molecular basis of evolutionary chromosome rearrangements, we performed molecular cytogenetic and sequence analyses of the breakpoint region that distinguishes human chromosome 3p12.3 and orangutan chromosome 2. FISH with region-specific BAC clones demonstrated that the breakpoint-flanking sequences are duplicated intrachromosomally on orangutan 2 and human 3q21 as well as at many pericentromeric and subtelomeric sites throughout the genomes. Breakage and rearrangement of the human 3p12.3-homologous region in the orangutan lineage were associated with a partial loss of duplicated sequences in the breakpoint region. Consistent with our FISH mapping results, computational analysis of the human chromosome 3 genomic sequence revealed three 3p12.3-paralogous sequence blocks on human chromosome 3q21 and smaller blocks on the short arm end 3p26-->p25. This is consistent with the view that sequences from an ancestral site at 3q21 were duplicated at 3p12.3 in a common ancestor of orangutan and humans. Our results show that evolutionary chromosome rearrangements are associated with microduplications and microdeletions, contributing to the DNA differences between closely related species.     

Human-specific amino acid changes found in 103 protein-coding genes

We humans have many characteristics that are different from those of the great apes. These human-specific characters must have arisen through mutations accumulated in the genome of our direct ancestor after the divergence of the last common ancestor with chimpanzee. Gene trees of human and great apes are necessary for extracting these human-specific genetic changes. We conducted a systematic analysis of 103 protein-coding genes for human, chimpanzee, gorilla, and orangutan. Nucleotide sequences for 18 genes were newly determined for this study, and those for the remaining genes were retrieved from the DDBJ/EMBL/GenBank database. The total number of amino acid changes in the human lineage was 147 for 26,199 codons (0.56%). The total number of amino acid changes in the human genome was, thus, estimated to be about 80,000. We applied the acceleration index test and Fisher's synonymous/nonsynonymous exact test for each gene tree to detect any human-specific enhancement of amino acid changes compared with ape branches. Six and two genes were shown to have significantly higher nonsynonymous changes at the human lineage from the acceleration index and exact tests, respectively. We also compared the distribution of the differences of the nonsynonymous substitutions on the human lineage and those on the great ape lineage. Two genes were more conserved in the ape lineage, whereas one gene was more conserved in the human lineage. These results suggest that a small proportion of protein-coding genes started to evolve differently in the human lineage after it diverged from the ape lineage.