Terrestrial vertebrates have two keratin gene clusters; striking differences in teleost fish

Keratins I and II form the largest subgroups of mammalian intermediate filament (IF) proteins and account as obligatory heteropolymers for the keratin filaments of epithelia. All human type I genes except for the K18 gene are clustered on chromosome 17q21, while all type II genes form a cluster on chromosome 12q13, that ends with the type I gene K18. Highly related keratin gene clusters are found in rat and mouse. Since fish seem to lack a keratin II cluster we screened the recently established draft genomes of a bird (chicken) and an amphibian (Xenopus). The results show that keratin I and II gene clusters are a feature of all terrestrial vertebrates. Because hair with its multiple hair keratins and inner root sheath keratins is a mammalian acquisition, the keratin gene clusters of chicken and Xenopus tropicalis have only about half the number of genes found in mammals. Within the type I clusters all genes have the same orientation. In type II clusters there is a rare gene of opposite orientation. Finally we show that the genes for keratins 8 and 18, which are the first expression pair in embryology, are not only adjacent in mammals, but also in Xenopus and three different fish. Thus neighboring K8 and K18 genes seem a feature shared by all vertebrates. In contrast to the two well defined keratin gene clusters of terrestrial vertebrates, three teleost fish show an excess of type I over type II genes, the lack of a keratin type II gene cluster and a striking dispersal of type I genes, that are probably the result of the teleost-specific whole genome duplication followed by a massive gene loss. This raises the question whether keratin gene clusters extend beyond the ancestral bony vertebrate to cartilage fish and lamprey. We also analyzed the complement of non-keratin IF genes of the chicken. Surprisingly, an additional nuclear lamin gene, previously overlooked by cDNA cloning, is documented on chromosome 10. The two splice variants closely resemble the lamin LIII a + b of amphibia and fish. This lamin gene is lost on the mammalian lineage.     

The organization of the keratin I and II gene clusters in placental mammals and marsupials show a striking similarity

The genomic database for a marsupial, the opossum Monodelphis domestica, is highly advanced. This allowed a complete analysis of the keratin I and keratin II gene cluster with some 30 genes in each cluster as well as a comparison with the human keratin clusters. Human and marsupial keratin gene clusters have an astonishingly similar organization. As placental mammals and marsupials are sister groups a corresponding organization is also expected for the archetype mammal. Since hair is a mammalian acquisition the following features of the cluster refer to its origin. In both clusters hair keratin genes arose at an interior position. While we do not know from which epithelial keratin genes the first hair keratins type-I and -II genes evolved, subsequent gene duplications gave rise to a subdomain of the clusters with many neighboring hair keratin genes. A second subdomain accounts in both clusters for 4 neighboring genes encoding the keratins of the inner root sheath (irs) keratins. Finally the hair keratin gene subdomain in the type-I gene cluster is interrupted after the second gene by a region encoding numerous genes for the high/ultrahigh sulfur hair keratin-associated proteins (KAPs). We also propose a tentative synteny relation of opossum and human genes based on maximal sequence conservation of the encoded keratins. The keratin gene clusters of the opossum seem to lack pseudogenes and display a slightly increased number of genes. Opossum keratin genes are usually longer than their human counterparts and also show longer intergenic distances.     

Remarkable expansions of an X-linked reproductive homeobox gene cluster in rodent evolution

Rhox is a recently identified cluster of 12 X-linked homeobox genes in mice. The expression pattern of Rhox genes during postnatal testis development corresponds to their chromosomal position, much like the colinear gene regulation of the Hox gene clusters during animal embryonic development. We here report the identification of 18 additional Rhox genes and 3 pseudogenes in mice. Comparative analyses of the mouse, rat, human, dog, cow, opossum, and chicken genomes suggest that the Rhox cluster originated in the common ancestor of primates and rodents. It subsequently underwent two remarkable expansions, first in the common ancestor of mice and rats and then in mice. Positive selection promoting amino acid substitutions was detected in some young Rhox genes, suggesting adaptive functional diversification. The recent expansions of the Rhox cluster provide an opportunity to study the mechanism and origin of colinear gene regulation, but they may also undermine the utility of mouse models for understanding the development and physiology of the human reproductive system.     

Comparative analysis of processed pseudogenes in the mouse and human genomes

Pseudogenes are important resources in evolutionary and comparative genomics because they provide molecular records of the ancient genes that existed in the genome millions of years ago. We have systematically identified approximately 5000 processed pseudogenes in the mouse genome, and estimated that approximately 60% are lineage specific, created after the mouse and human diverged. In both mouse and human genomes, similar types of genes give rise to many processed pseudogenes. These tend to be housekeeping genes, which are highly expressed in the germ line. Ribosomal-protein genes, in particular, form the largest sub-group. The processed pseudogenes in the mouse occur with a distinctly different chromosomal distribution than LINEs or SINEs - preferentially in GC-poor regions. Finally, the age distribution of mouse-processed pseudogenes closely resembles that of LINEs, in contrast to human, where the age distribution closely follows Alus (SINEs).     

Chromosomal localization of acidic and basic keratin genes of the domestic dog

Our laboratories are interested in characterizing genes involved in the myriad of heritable diseases affecting the domestic dog, Canis lupus familiaris, and in development of detailed genetic and physical maps of the canine genome. Included in these efforts is examination of conservation of the genetic organization, structure, and function of gene families involved in diseases of the canine skin, skeleton, and eye. To that end, study of the highly conserved keratin gene family was undertaken. Keratins belong to the superfamily of intermediate filaments and are the major structural proteins of the epidermis, hair, and nail. The keratins are highly conserved throughout vertebrate evolution both at the DNA and amino acid sequence levels. Mutations in genes encoding epithelial keratins are known to cause various diseases in humans, and similar histopathological presentations have been reported in the dog. The keratins are divided into two groups, type I (acidic) and type II (basic). In the human, the genes encoding the acidic and basic keratins are clustered on Chrs 17 and 12, respectively. The same genetic arrangement is seen in the mouse with the acidic and basic keratin gene clusters found on Chrs 11 and 15, respectively. Reported here are the chromosomal localization of acidic and basic canine keratin genes as well as supportive sequence data. Fluorescence in situ hybridization (FISH) experiments with clones isolated from a canine genomic library suggest that the acidic keratin gene cluster resides on CFA9 and the basic keratin gene cluster is located on CFA27. Received: 25 September 1998 / Accepted: 1 December 1998     

Cloning and characterization of a novel mouse Siglec, mSiglec-F: differential evolution of the mouse and human (CD33) Siglec-3-related gene clusters

A novel mouse Siglec (mSiglec-F) belonging to the subfamily of Siglec-3-related Siglecs has been cloned and characterized. Unlike most human Siglec-3 (hSiglec-3)-related Siglecs with promiscuous linkage specificity, mSiglec-F shows a strong preference for alpha2-3-linked sialic acids. It is predominantly expressed in immature cells of the myelomonocytic lineage and in a subset of CD11b (Mac-1)-positive cells in some tissues. As with previously cloned Siglec-3-related mSiglecs, the lack of strong sequence similarity to a singular hSiglec made identification of the human ortholog difficult. We therefore conducted a comprehensive comparison of Siglecs between the human and mouse genomes. The mouse genome contains eight Siglec genes, whereas the human genome contains 11 Siglec genes and a Siglec-like gene. Although a one-to-one orthologous correspondence between human and mouse Siglecs 1, 2, and 4 is confirmed, the Siglec-3-related Siglecs showed marked differences between human and mouse. We found only four Siglec genes and two pseudogenes in the mouse chromosome 7 region syntenic to the Siglec-3-related gene cluster on human chromosome 19, which, in contrast, contains seven Siglec genes, a Siglec-like gene, and thirteen pseudogenes. Although analysis of gene maps and exon structures allows tentative assignments of mouse-human Siglec ortholog pairs, the possibility of unequal genetic recombination makes the assignments inconclusive. We therefore support a temporary lettered nomenclature for additional mouse Siglecs. Current information suggests that mSiglec-F is likely a hSiglec-5 ortholog. The previously reported mSiglec-3/CD33 and mSiglec-E/MIS are likely orthologs of hSiglec-3 and hSiglec-9, respectively. The other Siglec-3-like gene in the cluster (mSiglec-G) is probably a hSiglec-10 ortholog. Another mouse gene (mSiglec-H), without an apparent human ortholog, lies outside of the cluster. Thus, although some duplications of Siglec-3-related genes predated separation of the primate and rodent lineages (about 80-100 million years ago), this gene cluster underwent extensive duplications in the primate lineage thereafter.     

The Mammalian 2′-5′ Oligoadenylate Synthetase Gene Family: Evidence for Concerted Evolution of Paralogous Oas1 Genes in Rodentia and Artiodactyla

Multiple 2′-5′ oligoadenylate (2-5A) synthetases are important components of innate immunity in mammals. Gene families encoding these proteins have previously been studied mainly in humans and mice. To reconstruct the evolution of this gene family in mammals, a search for additional 2-5A synthetase genes was performed in rat, cattle, pig, and dog. Twelve 2′-5′ oligoadenylate synthetase (Oas) genes were identified in the rat genome, including eight Oas1 genes, two Oas1 pseudogenes, single copies of Oas2 and Oas3, and two Oas-like genes, Oasl1 and Oasl2. Four OAS genes were detected in the pig genome and five OAS genes were found in both the cattle and dog genomes. An OAS3 gene was not found in either the cattle or the pig genome. While two tandemly duplicated OAS-like (OASL) genes were identified in the dog genome, only a single OASL orthologue was found in both the cattle and the pig genomes. The bovine and porcine OASL genes contain premature stop codons and encode truncated proteins, which lack the typical C-terminal double ubiquitin domains. The cDNA sequences of the rat, cattle, pig, and dog OAS genes were amplified, sequenced and compared with each other and with those in the human, mouse, horse, and chicken genomes. Evidence of concerted evolution of paralogous 2′-5′ oligoadenylate synthetase 1 genes was obtained in rodents (Rodentia) and even-toed ungulates (Artiodactyla). Calculations using the nonparametric Kolmogorov-Smirnov test suggested that the homogenization of paralogous OAS1 sequences was due to gene conversion rather than stabilizing selection. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. Reviewing Editor: Dr. Martin Kreitman     

Comparative genomics search for losses of long-established genes on the human lineage

Taking advantage of the complete genome sequences of several mammals, we developed a novel method to detect losses of well-established genes in the human genome through syntenic mapping of gene structures between the human, mouse, and dog genomes. Unlike most previous genomic methods for pseudogene identification, this analysis is able to differentiate losses of well-established genes from pseudogenes formed shortly after segmental duplication or generated via retrotransposition. Therefore, it enables us to find genes that were inactivated long after their birth, which were likely to have evolved nonredundant biological functions before being inactivated. The method was used to look for gene losses along the human lineage during the approximately 75 million years (My) since the common ancestor of primates and rodents (the euarchontoglire crown group). We identified 26 losses of well-established genes in the human genome that were all lost at least 50 My after their birth. Many of them were previously characterized pseudogenes in the human genome, such as GULO and UOX. Our methodology is highly effective at identifying losses of single-copy genes of ancient origin, allowing us to find a few well-known pseudogenes in the human genome missed by previous high-throughput genome-wide studies. In addition to confirming previously known gene losses, we identified 16 previously uncharacterized human pseudogenes that are definitive losses of long-established genes. Among them is ACYL3, an ancient enzyme present in archaea, bacteria, and eukaryotes, but lost approximately 6 to 8 Mya in the ancestor of humans and chimps. Although losses of well-established genes do not equate to adaptive gene losses, they are a useful proxy to use when searching for such genetic changes. This is especially true for adaptive losses that occurred more than 250,000 years ago, since any genetic evidence of the selective sweep indicative of such an event has been erased.     

Location of mouse and human genes corresponding to conserved canine olfactory receptor gene subfamilies

Olfactory receptors are G protein-coupled, seven-transmembrane-domain proteins that are responsible for binding odorants in the nasal epithelium. They are encoded by a large gene family, members of which are organized in several clusters scattered throughout the genomes of mammalian species. Here we describe the mapping of mouse sequences corresponding to four conserved olfactory receptor genes, each representing separate, recently identified canine gene subfamilies. Three of the four canine genes detected related gene clusters in regions of mouse Chromosomes (Chrs) 2, 9, and 10, near previously mapped mouse olfactory genes, while one detected a formerly unidentified gene cluster located on mouse Chr 6. In addition, we have localized two human gene clusters with homology to the canine gene, CfOLF4, within the established physical map of Chr 19p. Combined with recently published studies, these data link the four conserved olfactory gene subfamilies to homologous regions of the human, dog, and mouse genomes. Received: 10 September 1997 / Accepted: 29 December 1997     

Recognizing the pseudogenes in bacterial genomes

Pseudogenes are now known to be a regular feature of bacterial genomes and are found in particularly high numbers within the genomes of recently emerged bacterial pathogens. As most pseudogenes are recognized by sequence alignments, we use newly available genomic sequences to identify the pseudogenes in 11 genomes from 4 bacterial genera, each of which contains at least 1 human pathogen. The numbers of pseudogenes range from 27 in Staphylococcus aureus MW2 to 337 in Yersinia pestis CO92 (e.g. 1–8% of the annotated genes in the genome). Most pseudogenes are formed by small frameshifting indels, but because stop codons are A + T-rich, the two low-G + C Gram-positive taxa (Streptococcus and Staphylococcus) have relatively high fractions of pseudogenes generated by nonsense mutations when compared with more G + C-rich genomes. Over half of the pseudogenes are produced from genes whose original functions were annotated as ‘hypothetical’ or ‘unknown’; however, several broadly distributed genes involved in nucleotide processing, repair or replication have become pseudogenes in one of the sequenced Vibrio vulnificus genomes. Although many of our comparisons involved closely related strains with broadly overlapping gene inventories, each genome contains a largely unique set of pseudogenes, suggesting that pseudogenes are formed and eliminated relatively rapidly from most bacterial genomes.     

Gene decay in archaea

The gene-dense chromosomes of archaea and bacteria were long thought to be devoid of pseudogenes, but with the massive increase in available genome sequences, whole genome comparisons between closely related species have identified mutations that have rendered numerous genes inactive. Comparative analyses of sequenced archaeal genomes revealed numerous pseudogenes, which can constitute up to 8.6% of the annotated coding sequences in some genomes. The largest proportion of pseudogenes is created by gene truncations, followed by frameshift mutations. Within archaeal genomes, large numbers of pseudogenes contain more than one inactivating mutation, suggesting that pseudogenes are deleted from the genome more slowly in archaea than in bacteria. Although archaea seem to retain pseudogenes longer than do bacteria, most archaeal genomes have unique repertoires of pseudogenes.     

Rapid expansion of the Ly49 gene cluster in rat

The cytotoxic activity of mouse natural killer cells is regulated in part through cell surface molecules belonging to the Ly49 multigene family. In mice, the genomic sequence of the Ly49 gene cluster has been examined in detail and this analysis provided a model of the expansion of this multigene family. In the present study, we have analyzed a 1.8-Mb region of the draft rat genome revealing surprising differences in size and gene content between the mouse and the rat Ly49 clusters. The rat cluster contains at least 36 Ly49 genes, including pseudogenes, while dot-plot analysis of the cluster reveals an equidistant spacing of genes, suggesting that duplication of genes in the cluster occurred through a mechanism similar to that in the mouse. Phylogenetic analysis of the predicted rat genes reveals a number of distinct gene clusters and indicates that the majority of gene duplication events occurred after the divergence of mice and rats. Thus, the rodent Ly49 locus is subject to extremely rapid gene amplification and diversification.     

Analysis of the spermine synthase gene region in Fugu rubripes, Tetraodon fluviatilis, and Danio rerio

A prerequisite to understanding the evolution of the human X chromosome is the analysis of synteny of X-linked genes in different species. We have focused on the spermine synthase gene in human Xp22. 1. We show that whereas the human gene spans a genomic region of 54 kb, the Fugu rubripes gene is encompassed in a 4.7-kb region. However, we could not find conserved synteny between this region of human Xp22 and the equivalent F. rubripes region. A cosmid clone containing the F. rubripes gene does not contain other X-linked genes. Instead we identified homologs of human genes that are autosomally localized: the ryanodine receptor type I (RYRI), which is implicated in malignant hyperthermia and central core disease, and the HE6 gene. Comparison of the F. rubripes, Tetraodon fluviatilis, mouse, human, and Danio rerio 5'UTRs of spermine synthase highlights conserved sequences potentially involved in regulation. Interestingly, pseudogenes of this gene that are present in the human and mouse genomes seem to be absent in the compact F. rubripes genome. Analysis of a D. rerio PAC clone containing spermine synthase shows an intermediate genomic size in this fish. Sequence analysis of this PAC clone did not reveal other known genes: neither the RYRI gene, nor the HE6 gene, nor other human Xp22 genes were identified.     

Cloning of the human keratin 18 gene and its expression in nonepithelial mouse cells.

Human keratin 18 (K18) and the homologous mouse protein, Endo B, are intermediate filament subunits of the type I keratin class. Both are expressed in many simple epithelial cell types including trophoblasts, the first differentiated cell type to appear during mouse embryogenesis. The K18 gene was identified and cloned from among the 15 to 20 similar sequences identified within the human genome. The identity of the cloned gene was confirmed by comparing the sequence of the first two exons to the K18 cDNA sequence and transfecting the gene into various murine cell lines and verifying the encoded protein as K18 by immunoprecipitation and partial peptide mapping. The transfected K18 gene was expressed in mouse HR9 parietal endodermal cells and mouse fibroblasts even though the fibroblasts fail to express endogenous Endo B. S1 nuclease protection analysis indicated that mRNA synthesized from the transfected K18 gene is initiated at the same position as authentic K18 mRNA found in both BeWo trophoblastoma cells and HeLa cells. Pulse-chase experiments indicated that the human K18 protein is stable in murine parietal endodermal cells (HR9) which express EndoA, a complementary mouse type II keratin. Surprisingly, however, K18 was degraded when synthesized in cells which lack a type II keratin. This turnover of K18 may be an important mechanism by which epithelial cells maintain equal molar amounts of both type I and II keratins. In addition, the levels of the endogenous type I Endo B in parietal endodermal cells were compensatingly down regulated in the presence of the K18 protein, while the levels of the endogenous type II Endo A were not affected in any of the transfected cell lines.     

Classification and nomenclature of all human homeobox genes

Background

The homeobox genes are a large and diverse group of genes, many of which play important roles in the embryonic development of animals. Increasingly, homeobox genes are being compared between genomes in an attempt to understand the evolution of animal development. Despite their importance, the full diversity of human homeobox genes has not previously been described.

Results

We have identified all homeobox genes and pseudogenes in the euchromatic regions of the human genome, finding many unannotated, incorrectly annotated, unnamed, misnamed or misclassified genes and pseudogenes. We describe 300 human homeobox loci, which we divide into 235 probable functional genes and 65 probable pseudogenes. These totals include 3 genes with partial homeoboxes and 13 pseudogenes that lack homeoboxes but are clearly derived from homeobox genes. These figures exclude the repetitive DUX1 to DUX5 homeobox sequences of which we identified 35 probable pseudogenes, with many more expected in heterochromatic regions. Nomenclature is established for approximately 40 formerly unnamed loci, reflecting their evolutionary relationships to other loci in human and other species, and nomenclature revisions are proposed for around 30 other loci. We use a classification that recognizes 11 homeobox gene 'classes' subdivided into 102 homeobox gene 'families'.

Conclusion

We have conducted a comprehensive survey of homeobox genes and pseudogenes in the human genome, described many new loci, and revised the classification and nomenclature of homeobox genes. The classification scheme may be widely applicable to homeobox genes in other animal genomes and will facilitate comparative genomics of this important gene superclass.     

Evolution of immunoglobulin kappa chain variable region genes in vertebrates

The major source of immunoglobulin diversity is variation in DNA sequence among multiple copies of variable region (V) genes of the heavy- and light-chain multigene families. In order to clarify the evolutionary pattern of the multigene family of immunoglobulin light kappa chain V region (V kappa) genes, phylogenetic analyses of V kappa genes from humans and other vertebrate species were conducted. The results obtained indicate that the V kappa genes so far sequenced can be grouped into three major monophyletic clusters, the cartilaginous fish, bony fish and amphibian, and mammalian clusters, and that the cartilaginous fish cluster first separated from the rest of the V kappa genes and then the remaining two clusters diverged. The mammalian V kappa genes can further be divided into 10 V kappa groups, 7 of which are present in the human genome. Human and mouse V kappa genes from different V kappa groups are intermingled rather than clustered on the chromosome, and there are a large number of pseudogenes scattered on the chromosome. This indicates that the chromosomal locations of V kappa genes have been shuffled many times by gene duplication, deletion, and transposition in the evolutionary process and that many genes have become nonfunctional during this process. This mode of evolution is consistent with the model of birth-and-death evolution rather than with the model of concerted evolution. An analysis of duplicate V kappa functional genes and pseudogenes in the human genome has indicated that pseudogenes evolve faster than functional genes but that the rate of nonsynonymous nucleotide substitution in the complementarity-determining regions of V kappa genes has been enhanced by positive Darwinian selection.      

Evolution of Growth Hormone in Primates: The GH Gene Clusters of the New World Monkeys Marmoset (Callithrix jacchus) and White-Fronted Capuchin (Cebus albifrons)

The GH gene cluster in marmoset, Callithrix jacchus, comprises eight GH-like genes and pseudogenes and appears to have arisen as a consequence of gene duplications occurring independently of those leading to the human GH gene cluster. We report here the complete sequence of the marmoset GH gene locus, including the intergenic regions and 5′ and 3′ flanking sequence, and a study of the multiple GH-like genes of an additional New World monkey (NWM), the white-fronted capuchin, Cebus albifrons. The marmoset sequence includes 945 nucleotides (nt) of 5′ flanking sequence and 1596 nt of 3′ flanking sequence that are “unique”; between these are eight repeat units, including the eight GH genes/pseudogenes. The breakpoints between these repeats are very similar, indicating a regular pattern of gene duplication. These breakpoints do not correspond to those found in the much less regular human GH gene cluster. This and phylogenetic analysis of the repeat units within the marmoset gene cluster strongly support the independent origin of these gene clusters, and the idea that the episode of rapid evolution that occurred during GH evolution in primates preceded the gene duplications. The marmoset GH gene cluster also differs from that of human in having fewer and more evenly distributed Alu sequences (a single pair in each repeat unit) and a “P-element” upstream of every gene/pseudogene. In human there is no P-element upstream of the gene encoding pituitary GH, and these elements have been implicated in placental expression of the other genes of the cluster. The GH gene clusters in marmoset and capuchin appear to have arisen as the consequence of a single-gene duplication event, but in capuchin there was then a remarkable expansion of the GH locus, giving at least 40 GH-like genes and pseudogenes. Thus even among NWMs the GH gene cluster is very variable. [Reviewing Editor: Nicolas Galtier]     

Characterization of a cluster of human high/ultrahigh sulfur keratin-associated protein genes embedded in the type I keratin gene domain on chromosome 17q12-21

Low stringency screening of a human P1 artificial chromosome library using a human hair keratin-associated protein (hKAP1.1A) gene probe resulted in the isolation of six P1 artificial chromosome clones. End sequencing and EMBO/GenBank(TM) data base analysis showed these clones to be contained in four previously sequenced human bacterial artificial chromosome clones present on chromosome 17q12-21 and arrayed into two large contigs of 290 and 225 kilobase pairs (kb) in size. A fifth, partially sequenced human bacterial artificial chromosome clone data base sequence overlapped and closed both of these contigs. One end of this 600-kb cluster harbored six gene loci for previously described human type I hair keratin genes. The other end of this cluster contained the human type I cytokeratin K20 and K12 gene loci. The center of the cluster, starting 35 kb downstream of the hHa3-I hair keratin gene, contained 37 genes for high/ultrahigh sulfur hair keratin-associated proteins (KAPs), which could be divided into a total of 7 KAP multigene families based on amino acid homology comparisons with previously identified sheep, mouse, and rabbit KAPs. To date, 26 human KAP cDNA clones have been isolated through screening of an arrayed human scalp cDNA library by means of specific 3'-noncoding region polymerase chain reaction probes derived from the identified KAP gene sequences. This screening also yielded four additional cDNA sequences whose genes were not present on this gene cluster but belonged to specific KAP gene families present on this contig. Hair follicle in situ hybridization data for single members of five different KAP multigene families all showed localization of the respective mRNAs to the upper cortex of the hair shaft.