The compositional patterns of the avian genomes and their evolutionary implications

The compositional distributions of large (main-band) DNA fragments from eight birds belonging to eight different orders (including both paleognathous and neognathous species) are very broad and extremely close to each other. These findings, which are paralleled by the compositional similarity of homologous coding sequences and their codon positions, support the idea that birds are a monophyletic group.The compositional distribution of third-codon positions of genes from chicken, the only avian species for which a relatively large number of coding sequences is known, is very broad and bimodal, the minor GC-richer peak reaching 100% GC. The very high compositional heterogeneity of avian genomes is accompanied (as in the case of mammalian genomes) by a very high speciation rate compared to cold-blooded vertebrates which are characterized by genomes that are much less heterogeneous. The higher GC levels attained by avian compared to mammalian genomes might be correlated with the higher body temperature (41–43°C) of birds compared to mammals (37°C).A comparison of GC levels of coding sequences and codon positions from man and chicken revealed very close average GC levels and standard deviations. Homologous coding sequences and codon positions from man and chicken showed a surprisingly high degree of compositional similarity which was, however, higher for GC-poor than for GC-rich sequences. This indicates that GC-poor isochores of warm-blooded vertebrates reflect the composition of the isochores of the genome of the common reptilian ancestor of mammals and birds, which underwent only a small compositional change at the transition from cold- to warm-blooded vertebrates. In contrast, the GC-rich isochores of birds and mammals are the result of large compositional changes at the same evolutionary transition, where were in part different in the two classes of warm-blooded vertebrates.Correspondence to: G. Bernaadi     

Identification of natural killer cell receptor genes in the genome of the marsupial Tasmanian devil (Sarcophilus harrisii)

Within the mammalian immune system, natural killer (NK) cells contribute to the first line of defence against infectious agents and tumours. Their activity is regulated, in part, by cell surface NK cell receptors. NK receptors can be divided into two unrelated, but functionally analogous superfamilies based on the structure of their extracellular ligand-binding domains. Receptors belonging to the C-type lectin superfamily are predominantly encoded in the natural killer complex (NKC), while receptors belonging to the immunoglobulin superfamily are predominantly encoded in the leukocyte receptor complex (LRC). Natural killer cell receptors are emerging as a rapidly evolving gene family which can display significant intra- and interspecific variation. To date, most studies have focused on eutherian mammals, with significantly less known about the evolution of these receptors in marsupials. Here, we describe the identification of 43 immunoglobulin domain-containing LRC genes in the genome of the Tasmanian devil (Sarcophilus harrisii), the largest remaining marsupial carnivore and only the second marsupial species to be studied. We also identify orthologs of NKC genes KLRK1, CD69, CLEC4E, CLEC1B, CLEC1A and an ortholog of an opossum NKC receptor. Characterisation of these regions in a second, distantly related marsupial provides new insights into the dynamic evolutionary histories of these receptors in mammals. Understanding the functional role of these genes is also important for the development of therapeutic agents against Devil Facial Tumour Disease, a contagious cancer that threatens the Tasmanian devil with extinction.     

Orthologous nature of mammaliam insulin genes

Summary Investigation of 7 insulin sequences from a bony fish, a bird, and 5 mammalian species showed that guinea pig and coypu insulin, that have a strongly divergent primary and quarternary structure, are not the result of gene duplication in an ancient vertebrate or invertebrate ancestor but that they diverged from the other mammals after divergence of the mammals from the other vertebrates. After this divergence both insulins underwent evolution at a highly increased rate.     

Evolution of the Monotremes: Phylogenetic Relationship to Marsupials and Eutherians, and Estimation of Divergence Dates Based on α-Lactalbumin Amino Acid Sequences

The amino acid sequences of the -lactalbumins of the echidna, Tachyglossus aculeatus, and the platypus, Ornithorhynchus anatinus, were compared with each other and with those of 13 eutherian and 3 marsupial species. Phylogenetic parsimony analyses, in which selected mammalian lysozymes were used as outgroups, yielded trees whose consensus indicated that the two monotremes are sister taxa to marsupials and eutherians and that the latter two clades are each other's closest relatives. The data do not support the notion of a Marsupionta (monotreme–marsupial) clade. Pairwise comparison between the -lactalbumins yielded maximum-likelihood distances from which divergence dates were estimated on the basis of three calibration points. The distance data support the view that the echidna and platypus lineages diverged from their last common ancestor at least 50 to 57 Ma (million years ago) and that monotremes diverged from marsupials and eutherian mammals about 163 to 186 Ma.     

Characterization and cloning of lgp110, a lysosomal membrane glycoprotein from mouse and rat cells

lgp110 is a heavily glycosylated intrinsic protein of lysosomal membranes. Initially defined by monoclonal antibodies against mouse liver lysosomes, it consists of a 45-kilodalton core polypeptide with O-linked and 17 asparagine-linked oligosaccharide side chains in mouse cells. Sialic acid residues make the mature protein extremely acidic, with an isoelectric point of between 2 and 4 in both normal tissues and most cultured cell lines. Partial sequencing of mouse lgp110 allowed oligonucleotide probes to be constructed for the screening of several mouse cDNA libraries. A partial cDNA clone for mouse lgp110 was found and used for additional library screening, generating a cDNA clone covering all of the coding sequence of mature rat lgp110 as well as genomic clones covering most of the mouse gene. These new clones bring to seven the number of lysosomal membrane proteins whose amino acid sequences can be deduced, and two distinct but highly similar groups (designated lgp-A and lgp-B) can now be defined. Sequence comparisons suggest that differences within each group reflect species variations of the same protein and that lgp-A and lgp-B probably diverged from a common ancestor prior to the evolup4f1ary divergence of birds and mammals. Individual cells and individual lysosomes possess both lgp-A and lgp-B, suggesting that these two proteins have different functions. Mouse lgp110 is encoded by at least seven exons; intron positions suggest that the two homologous ectodomains of each lgp arose through gene duplication.     

The Evolutionary Origins of Detoxifying Enzymes: THE MAMMALIAN SERUM PARAOXONASES (PONs) RELATE TO BACTERIAL HOMOSERINE LACTONASES*

Serum paraoxonases (PONs) are detoxifying lactonases that were first identified in mammals. Three mammalian families are known, PON1, 2, and 3 that reside primarily in the liver. They catalyze essentially the same reaction, lactone hydrolysis, but differ in their substrate specificity. Although some members are highly specific, others have a broad specificity profile. The evolutionary origins and substrate specificities of PONs therefore remain poorly understood. Here, we report a newly identified family of bacterial PONs, and the reconstruction of the ancestor of the three families of mammalian PONs. Both the mammalian ancestor and the characterized bacterial PONX_OCCAL were found to efficiently hydrolyze N-acyl homoserine lactones that mediate quorum sensing in many bacteria, including pathogenic ones. The mammalian PONs may therefore relate to a newly identified family of bacterial, PON-like “quorum-quenching” lactonases. The appearance of PONs in metazoa is likely to relate to innate immunity rather than detoxification. Unlike the bacterial PON, the mammalian ancestor also hydrolyzes, with low efficiency, lactones other than homoserine lactones, thus preceding the detoxifying functions that diverged later in two of the three mammalian families. The bifunctionality of the mammalian ancestor and the trade-off between the quorum-quenching and detoxifying lactonase activities explain the broad and overlapping specificities of some mammalian PONs versus the singular specificity of others.     

Precursors of mesotocin and vasotocin in birds: Identification of VLDV- and MSEL-neurophysins in chicken,goose, and ostrich

Precursors of neurohypophysial hormones are small proteins processed into nonapeptide hormones and neurophysins during axonal transport to the neurohypophysis. In mammals, oxytocin is associated with VLDV-neurophysin and vasopressin with MSEL-neurophysin. In birds, mesotocin and vasotocin are found instead of mammalian oxytocin and vasopressin. From goose, chicken and ostrich posterior pituitary glands, two types of neurophysins related to mammalian VLDV-and MSEL-neurophysins, respectively, have been identified by their N-terminal sequences. It is assumed that, as in mammals, hormonal peptide and the first 9 residues of the corresponding neurophysin are encoded by a common exon and that mesotocin and vasotocin, evolutionary predecessors of oxytocin and vasopressin, are associated in the precursors with VLDV-neurophysin and MSEL-neurophysin, respectively.     

Molecular evolution and phylogeny of the human AIDS viruses LAV,HTLV-III,and ARV

Summary A phylogenetic tree for the human lymphadenopathy-associated virus (LAV), the human T-cell lymphotrophic virus type III (HTLV-III), and the acquired immune deficiency syndrome (AIDS)-associated retrovirus (ARV) has been constructed from comparisons of the amino acid sequences of their gag proteins. A method is proposed for estimating the divergence times among these AIDS viruses and the rates of nucleotide substitution for their RNA genomes. The analysis indicates that the LAV and HTLV-III strains diverged from one another after 1977 and that their common ancestor diverged from the ARV virus no more than 10 years earlier. Hence, the evolutionary diversity among strains of the AIDS viruses apparently has been generated within the last 20 years. It is estimated that the genome of the AIDS virus has a nucleotide substitution rate on the order of 10^–3 per site per year, with the rate in the second half of the genome being double that in the first half.     

Evolutionary diversification of the vertebrate transferrin multi-gene family

In a phylogenetic analysis of vertebrate transferrins (TFs), six major clades (subfamilies) were identified: (a) S, the mammalian serotransferrins; (b) ICA, the mammalian inhibitor of carbonic anhydrase (ICA) homologs; (c) L, the mammalian lactoferrins; (d) O, the ovotransferrins of birds and reptiles; (e) M, the melanotransferrins of bony fishes, amphibians, reptiles, birds, and mammals; and (f) M-like, a newly identified TF subfamily found in bony fishes, amphibians, reptiles, and birds. A phylogenetic tree based on the joint alignment of N-lobes and C-lobes supported the hypothesis that three separate events of internal duplication occurred in vertebrate TFs: (a) in the common ancestor of the M subfamily, (b) in the common ancestor of the M-like subfamily, and (c) in the common ancestor of other vertebrate TFs. The S, ICA, and L subfamilies were found only in placental mammals, and the phylogenetic analysis supported the hypothesis that these three subfamilies arose by gene duplication after the divergence of placental mammals from marsupials. The M-like subfamily was unusual in several respects, including the presence of a uniquely high proportion of clade-specific conserved residues, including distinctive but conserved residues in the sites homologous to those functioning in carbonate binding of human serotransferrin. The M-like family also showed an unusually high proportion of cationic residues in the positively charged region corresponding to human lactoferrampin, suggesting a distinctive role of this region in the M-like subfamily, perhaps in antimicrobial defense.     

Mammalian sex chromosomes: Evolution of organization and function

Comparisons of chromosome size, morphology and gene arrangements between mammals of different species permit us to deduce the genome characteristics of the common ancestor, and to chart the changes that have occurred during the divergence of the two lineages. The more distantly related are the species compared, the more remote the common ancestor whose characteristics can be deduced. This means that, providing there are sufficient similarities to warrant comparison, the more divergent the species compared, the more significant the contribution to our understanding of the organization of an ancestral mammalian genome and the process of mammalian genome evolution. One of the genetic surprises of the last decade was the discovery that, although gross karyotypes of distantly related orders of eutherian mammals (e.g. cat, cow, rabbit, man) have diverged extensively, gene mapping studies reveal the presence of large chromosome segments conserved across at least 60 million years (O'Brien et al. 1988). This finding makes it worthwhile to extend genetic comparisons to the two groups of mammals most distantly related to eutherian mammals--marsupials and monotremes. Here we will review comparisons of the sex chromosomes in these three major groups of extant mammals, and show how they have led us to a new view of the evolution of mammalian sex chromosome organization and function in sex determination and X chromosome inactivation.     

Comparative genomic analysis of two avian (quail and chicken) MHC regions

We mapped two different quail Mhc haplotypes and sequenced one of them (haplotype A) for comparative genomic analysis with a previously sequenced haplotype of the chicken Mhc. The quail haplotype A spans 180 kb of genomic sequence, encoding a total of 41 genes compared with only 19 genes within the 92-kb chicken Mhc. Except for two gene families (B30 and tRNA), both species have the same basic set of gene family members that were previously described in the chicken "minimal essential" Mhc. The two Mhc regions have a similar overall organization but differ markedly in that the quail has an expanded number of duplicated genes with 7 class I, 10 class IIB, 4 NK, 6 lectin, and 8 B-G genes. Comparisons between the quail and chicken Mhc class I and class II gene sequences by phylogenetic analysis showed that they were more closely related within species than between species, suggesting that the quail Mhc genes were duplicated after the separation of these two species from their common ancestor. The proteins encoded by the NK and class I genes are known to interact as ligands and receptors, but unlike in the quail and the chicken, the genes encoding these proteins in mammals are found on different chromosomes. The finding of NK-like genes in the quail Mhc strongly suggests an evolutionary connection between the NK C-type lectin-like superfamily and the Mhc, providing support for future studies on the NK, lectin, class I, and class II interaction in birds.     

Molecular Evolution of Vertebrate Goose-Type Lysozyme Genes

We have found that mammalian genomes contain two lysozyme g genes. To better understand the function of the lysozyme g genes we have examined the evolution of this small gene family. The lysozyme g gene structure has been largely conserved during vertebrate evolution, except at the 5' end of the gene, which varies in number of exons. The expression pattern of the lysozyme g gene varies between species. The fish lysozyme g sequences, unlike bird and mammalian lysozyme g sequences, do not predict a signal peptide, suggesting that the encoded proteins are not secreted. The fish sequences also do not conserve cysteine residues that generate disulfide bridges in the secreted bird enzymes, supporting the hypothesis that the fish enzymes have an intracellular function. The signal peptide found in bird and mammalian lysozyme g genes may have been acquired as an exon in the ancestor of birds and mammals, or, alternatively, an exon encoding the signal peptide has been lost in fish. Both explanations account for the change in gene structure between fish and tetrapods. The mammalian lysozyme g sequences were found to have evolved at an accelerated rate, and to have not perfectly conserved the known active site catalytic triad of the bird enzymes. This observation suggests that the mammalian enzymes may have altered their biological function, as well.     

Channel catfish leukocyte immune-type receptors contain a putative MHC class I binding site

The recent identification of a large and diverse family of leukocyte immune-type receptors (IpLITRs) in channel catfish (Ictalurus punctatus) indicates that immunoglobulin superfamily (IgSF) members related to both mammalian Fc receptors (FcRs) and leukocyte receptor complex (LRC)-encoded proteins exist in fish. In the present study, it was found that IpLITR messages were preferentially up regulated in catfish peripheral blood leukocytes (PBL) and clonal cytotoxic T cells (CTL) after alloantigen stimulation. Detailed sequence analyses of the expressed IpLITR cDNAs from two clonal CTL lines indicated an unexpectedly large array of putative activatory and inhibitory IpLITR-types containing variable numbers of extracellular immunoglobulin (Ig)-like domains. Importantly, all expressed IpLITRs shared similar membrane distal Ig domains (i.e., D1 and D2), suggesting that they may bind a common type of ligand. Sequence alignments and comparative homology modeling revealed that IpLITR domains, D1 and D2, have similar predicted 3-D structural properties with the corresponding domains of the human LRC-encoded leukocyte Ig-like receptor (LILR) family. Furthermore, conservation of key major histocompatibility class I (MHC I)-binding residues were located at similar positions within the membrane distal tip of D1 between representative IpLITRs and group 1 LILRs. Taken together, these results suggest that fish LITRs have an orthologous relationship to LRC-encoded receptors such as the human LILRs and could potentially function as a diverse family of MHC class I-binding receptors.     

The evolution of lysozyme and alpha-lactalbumin

From the analysis of phylogenetic trees constructed from the amino acid sequences and metal-binding properties of various lysozymes c and alpha-lactalbumins, it was found that before the divergence of the lineages of birds and mammals, calcium-binding lysozyme diverged from non-calcium-binding lysozyme. alpha-Lactalbumin evolved from the calcium-binding lysozyme along the mammalian lineage after the divergence of birds and mammals. Rapid evolution took place, not in the process of acquisition of the activity of alpha-lactalbumin, but after the loss of lysozyme activity, due to the change in the distribution of selective pressure on each amino acid site. A general process for the change in function of a protein during evolution is suggested to be as follows: after duplication of the gene, one of their protein products acquires a new function, besides that already present; the old function is eventually lost.     

Ancient origin of lactalbumin from lysozyme: Analysis of DNA and amino acid sequences

Summary Parsimony trees relating DNA sequences coding for lysozymesc and -lactalbumins suggest that the gene duplication that allowed lactalbumin to evolve from lysozyme preceded the divergence of mammals and birds. Comparisons of the amino acid sequences of additional lysozymes and lactalbumins are consistent with this view. When all base positions are considered, the probability that the duplication leading to the lactalbumin gene occurred after the start to mammalian evolution is estimated to be 0.05–0.10. Elimination of the phylogenetic noise generated by fast evolution and compositional bias at third positions of codons reduced this probability to 0.002–0.03. Thus the gene duplication may have long preceded the acquisition of lactalbumin function.     

Evolutionary relationships of class II major-histocompatibility-complex genes in mammals

The major histocompatibility complex (MHC) class II molecule consists of noncovalently associated alpha and beta chains. In mammals studied so far, the class II MHC can be divided into a number of regions, each containing one or more alpha-chain genes (A genes) and beta-chain genes (B genes), and it has been known for some time that orthologous relationships exist between genes in corresponding regions from different mammalian species. A phylogenetic analysis of DNA sequences of class II A and B genes confirmed these relationships; but no such orthologous relationship was observed between the B genes of mammals and those of birds. Thus, the class II regions have diverged since the separation of birds and mammals (approximately 300 Mya) but before the radiation of the placental mammalian orders (60-80 Mya). Comparison of the phylogenetic trees for A and B genes revealed an unexpected characteristic of DP-region genes: DPB genes are most closely related to DQB genes, whereas DPA chain genes are most closely related to DRA-chain genes. Thus, the DP region seems to have originated through a recombinational event which brought together a DQB gene and a DRA gene (perhaps approximately 120 Mya). The 5' untranslated region of all class II genes includes sequences which are believed to be important in regulating class II gene expression but which are not conserved in known pseudogenes. These sequences are conserved to an extraordinary degree in the human DQB1 gene and its mouse homologue A beta 1, suggesting that regulation of expression of this locus may play a key role in expression of the entire class II MHC.     

Complex biogeographic history of a Holarctic passerine

Our analysis of the ND2 sequences revealed six clades within winter wrens (Troglodytes troglodytes). These clades corresponded to six geographical regions: western Nearctic, eastern Nearctic, eastern Asia, Nepal, Caucasus and Europe, and differed by 3-8.8% of sequence divergence. Differences among regions explained 96% of the sequence variation in winter wren. Differences among individuals within localities explained 3% of the sequence variation, and differences among localities within regions explained 1%. Grouping sequences into subspecies instead of localities did not change these proportions. Proliferation of the six clades coincided with Early and Middle Pleistocene glaciations. The distribution of winter wren clades can be explained by a series of five consecutive vicariant events. Western Nearctic wrens diverged from the Holarctic ancestor 1.6 Myr before the present time (MYBP). Eastern Nearctic and Palaearctic wrens diverged 1 MYBP. Eastern and western Palaearctic birds diverged 0.83 MYBP. Nepalese and east Asian wrens diverged 0.67 MYBP, and Caucasian birds diverged from European wrens 0.54 MYBP. The winter wren has a much greater degree of inter- and intracontinental differentiation than the three other Holarctic birds studied to date--dunlin (Calidris alpina), common raven (Corvus corax) and three-toed woodpecker (Picoides trydactylus)--and represents an example of cryptic speciation that has been overlooked.     

Molecular evolution of human visual pigment genes

By comparing the published DNA sequences for (a) the genes encoding the human visual color pigments (red, green, and blue) with (b) the genes encoding human, bovine, and Drosophila rhodopsins, a phylogenetic tree for the mammalian pigment genes has been constructed. This evolutionary tree shows that the common ancestor of the visual color pigment genes diverged first from that of the rhodopsin genes; then the common ancestor of the red and green pigment genes and the ancestor of the blue pigment gene diverged; and finally the red and green pigment genes diverged from each other much more recently. Nucleotide substitutions in the rhodopsin genes are best explained by the neutral theory of molecular evolution. However, important functional adaptations seem to have occurred twice during the evolution of the color pigment genes in humans: first, to the common ancestor of the three color pigment genes after its divergence from the ancestor of the rhodopsin gene and, second, to the ancestor of the red pigment gene after its divergence from that of the green pigment gene.     

Chicken CD69 and CD94/NKG2-like genes in a chromosomal region syntenic to mammalian natural killer gene complex

In mammals, natural killer (NK) cell C-type lectin receptors were encoded in a gene cluster called natural killer gene complex (NKC). The NKC is not reported in chicken yet. Instead, NK receptor genes were found in the major histocompatibility complex. In this study, two novel chicken C-type lectin-like receptor genes were identified in a region on chromosome 1 that is syntenic to mammalian NKC region. The chromosomal locations were validated with fluorescent in situ hybridization. Based on 3D structure modeling, sequence homology, chromosomal location, and phlylogenetic analysis, one receptor is the orthologue of mammalian cluster of differentiation 69 (CD69), and the other is highly homologous to CD94 and NKG2. Like CD94/NKG2 gene found in teleostean fishes, chicken CD94/NKG2 has the features of both human CD94 and NKG2A. Unlike mammalian NKC, these two chicken C-type lectin receptors are not closely linked but separated by 42 million base pairs according to the chicken draft genome sequence. The arrangement of several other genes that are located outside the mammalian NKC is conserved among chicken, human, and mouse. The chicken NK C-type lectin-like receptors in the NKC syntenic region indicate that this chromosomal region existed before the divergence between mammals and aves. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. The nucleotide sequences have been submitted to the GenBank nucleotide sequence database under the accession number chicken CD69 (DQ156495), CD94/NKG2 (DQ156496), and CD94/NKG2 variant (DQ241793).