Orthologous Positions Research Articles

Characterization of major royal jelly protein-like DNA sequences in Apis dorsata

SUMMARYMajor proteins of royal jelly (MRJP) are well characterized in the honey bee, Apis mellifera. They are secreted by hypopharyngeal glands and provide major proteinaceous components of the larval food. To date nothing is known about royal jelly proteins in other bee species. To fill this gap, we screened the genome of the Asian giant honey bee, Apis dorsata, for the two MRJP protein coding sequences MRJP3 and MRJP5. With PCR primers originally designed for A. mellifera we were able to amplify the presumably orthologous sequences in A. dorsata. Sequence analyses revealed high similarity in the two honey bees. As described in A. mellifera we found an orthologous extensive repetitive region in the MRJP3 sequence of A. dorsata. The number of reiterated repetitive motifs is highly variable in different A. dorsata individuals and is therefore well suited as a marker in population studies. The MRJP5 sequence region contains another repetitive region at the orthologous position as described for A. mellifera. These findings confirm the presence and indicate the functional equivalence of MRJP proteins in A. dorsata compared to A. mellifera. The origin of the MRJP protein family can therefore be dated back probably 22 million years ago to a common ancestor of A. mellifera and A. dorsata.

Scoring pairwise genomic sequence alignments.

The parameters by which alignments are scored can strongly affect sensitivity and specificity of alignment procedures. While appropriate parameter choices are well understood for protein alignments, much less is known for genomic DNA sequences. We describe a straightforward approach to scoring nucleotide substitutions in genomic sequence alignments, especially human-mouse comparisons. Scores are obtained from relative frequencies of aligned nucleotides observed in alignments of non-coding, non-repetitive genomic regions, and can be theoretically motivated through substitution models. Additional accuracy can be attained by down-weighting alignments characterized by low compositional complexity. We also describe an evaluation protocol that is relevant when alignments are intended to identify all and only the orthologous positions. One particular scoring matrix, called HOXD70, has proven to be generally effective for human-mouse comparisons, and has been used by the PipMaker server since July, 2000. We discuss but leave open the problem of effectively scoring regions of strongly biased nucleotide composition, such as low G + C content.

Large-scale analysis of the Alu Ya5 and Yb8 subfamilies and their contribution to human genomic diversity

We have utilized computational biology to screen GenBank for the presence of recently integrated Ya5 and Yb8 Alu family members. Our analysis identified 2640 Ya5 Alu family members and 1852 Yb8 Alu family members from the draft sequence of the human genome. We selected a set of 475 of these elements for detailed analyses. Analysis of the DNA sequences from the individual Alu elements revealed a low level of random mutations within both subfamilies consistent with the recent origin of these elements within the human genome. Polymerase chain reaction assays were used to determine the phylogenetic distribution and human genomic variation associated with each Alu repeat. Over 99 % of the Ya5 and Yb8 Alu family members were restricted to the human genome and absent from orthologous positions within the genomes of several non-human primates, confirming the recent origin of these Alu subfamilies in the human genome. Approximately 1 % of the analyzed Ya5 and Yb8 Alu family members had integrated into previously undefined repeated regions of the human genome. Analysis of mosaic Yb8 elements suggests gene conversion played an important role in generating sequence diversity among these elements. Of the 475 evaluated elements, a total of 106 of the Ya5 and Yb8 Alu family members were polymorphic for insertion presence/absence within the genomes of a diverse array of human populations. The newly identified Alu insertion polymorphisms will be useful tools for the study of human genomic diversity.

A HERV-K provirus in chimpanzees, bonobos and gorillas, but not humans.

Evidence from DNA sequencing studies strongly indicated that humans and chimpanzees are more closely related to each other than either is to gorillas [1-4]. However, precise details of the nature of the evolutionary separation of the lineage leading to humans from those leading to the African great apes have remained uncertain. The unique insertion sites of endogenous retroviruses, like those of other transposable genetic elements, should be useful for resolving phylogenetic relationships among closely related species. We identified a human endogenous retrovirus K (HERV-K) provirus that is present at the orthologous position in the gorilla and chimpanzee genomes, but not in the human genome. Humans contain an intact preintegration site at this locus. These observations provide very strong evidence that, for some fraction of the genome, chimpanzees, bonobos, and gorillas are more closely related to each other than they are to humans. They also show that HERV-K replicated as a virus and reinfected the germline of the common ancestor of the four modern species during the period of time when the lineages were separating and demonstrate the utility of using HERV-K to trace human evolution.

Physical map and expression profile of genes of the telomeric class I gene region of the rat MHC.

The rat is an important model for studying organ graft rejection and susceptibility to certain complex diseases. The MHC, the RT1 complex, plays a decisive role in controlling these traits. We have cloned the telomeric class I region of the RT1 complex, RT1-C/E/M, of the BN inbred rat strain in a contig of overlapping P1-derived artificial chromosome clones encompassing approximately 2 Mb, and present a physical map of this MHC region. Forty-five class I exon 4-hybridizing BAM:HI fragments were detected, including the previously known rat class I genes RT1-E, RT-BM1, RT1-N, RT1-M2, RT1-M3, and RT1-M4. Twenty-six non-class I genes known to map to the corresponding part of the human and mouse MHC were tested and could be fine mapped in the RT1-C/E/M region at orthologous position. Four previously known microsatellite markers were fine mapped in the RT1-C/E/M region and found to occur in multiple copies. In addition, a new, single-copy polymorphic microsatellite has been defined. The expression profiles of several class I genes and the 26 non-class I genes were determined in 13 different tissues and exhibited restricted patterns in most cases. The data provide further molecular information on the MHC for analyzing disease susceptibility and underline the usefulness of the rat model.

Comparative genomic sequencing reveals a strikingly similar architecture of a conserved syntenic region on human chromosome 11p15.3 (including gene ST5) and mouse chromosome 7

Comparative genomics is a superior way to identify phylogenetically conserved features like genes or regions involved in gene regulation. The comparison of extended orthologous chromosomal regions should also reveal other characteristic traits essential for chromosome or gene function. In the present study we have sequenced and compared a region of conserved synteny from human chromosome 11p15.3 and mouse chromosome 7. In human, this region is known to contain several genes involved in the development of various disorders like Beckwith-Wiedemann overgrowth syndrome and other tumor diseases. Furthermore, in the neighboring chromosome region 11p15.5 extensive imprinting of genes has been reported which might extend to region 11p15.3. The analysis of approximately 730 kb in human and 620 kb in mouse led to the identification of eleven genes. All putative genes found in the mouse DNA were also present in the same order and orientation in the human chromosome. However, in the human DNA one putative gene of unknown function could be identified which is not present in the orthologous position of the mouse chromosome. The sequence similarity between human and mouse is higher in transcribed and exon regions than in non-transcribed segments. Dot plot analysis, however, reveals a surprisingly well-conserved sequence similarity over the entire analyzed region. In particular, the positions of CpG islands, short regions of very high GC content in the 5′ region of putative genes, are similar in human and mouse. With respect to base composition, two distinct segments of significantly different GC content exist as well in human as in the mouse. With a GC content of 45% the one segment would correspond to “isochore H1” and the other segment (39% GC in human, 40% GC in mouse) to “isochore L1/L2”. The gene density (one gene per 66 kb) is slightly higher than the average calculated for the complete human genome (one gene per 90 kb). The comparison of the number and distribution of repetitive elements shows that the proportion of human DNA made up by interspersed repeats (43.8%) is significantly higher than in the corresponding mouse DNA (30.1%). This partly explains why the human DNA is longer between the landmark genes used to define the orthologous positions in human and mouse.

Evolution of an Alu DNA element of type Sx in the lineage of primates and the origin of an associated tetranucleotide microsatellite.

A 394-bp DNA fragment, which in human is on chromosome 6 near the MOG (myelin oligodendrocyte glycoprotein) gene and encompasses an Alu element and an associated tetranucleotide microsatellite, was sequenced from a large range of primate species to follow its evolutionary divergence and to understand the origin of the microsatellite. This Alu element is found at the same orthologous position in all primates sequenced, but the tetranucleotide repeat is present only in Catarrhini between the 3'-oligo(dA) of the Alu element and the 3' flanking direct repeat. Little intraspecific variation was found. Sequence identity values for this orthologous primate Alu averaged 90% (82-99%) with transitions comprising between 70% and 100% of the observed nucleotide substitutions. Although the insertion of the Alu element predates the separation of these species, the original sequence of the site of integration can still be identified. This identification of the direct repeats suggests an active role of the oligo(dA) of the Alu element in the origin of the tetranucleotide repeats. The microsatellite probably appeared after the insertion of the Alu element, early in the lineage leading to the common ancestor of the hominoids and the Old World monkeys.

The Bov-B LINEs found in Vipera ammodytes toxic PLA 2 genes are widespread in snake genomes

In the fourth intron of two toxic Vipera ammodytes PLA 2 genes a Ruminantia specific 5′-truncated Bov-B LINE element was identified. Southern blot analysis of Bov-B LINE distribution in vertebrates shows that, apart from the Ruminantia, it is limited to Viperidae snakes ( V. ammodytes, Vipera palaestinae, Echis coloratus, Bothrops alternatus, Trimeresurus flavoviridis and Trimeresurus gramineus). The copy number of the 3′ end of Bov-B LINE in the V. ammodytes genome is between 62 000 and 75 000. At orthologous positions in other snake PLA 2 genes the Bov-B LINE element is absent, indicating that its retrotransposition in the V. ammodytes PLA 2 gene locus has occurred quite recently, about 5 Myr ago. The amplification of Bov-B LINEs in snakes may have occurred before the divergence of the Viperinae and Crotalinae subfamilies. Due to its wide distribution in Viperidae snakes it should be a valuable phylogenetic marker. The neighbour-joining phylogenetic tree shows two clusters of truncated Bov-B LINE, a Bovidae and a snake cluster, indicating an early horizontal transfer of this transposable element.

A tetrodotoxin-resistant voltage-gated sodium channel from human dorsal root ganglia, hPN3/ SCN10A

Neuropathic pain may be produced, at least in part, by the increased activity of primary afferent neurons. Studies have suggested that an accumulation of voltage-gated sodium channels at the site of peripheral nerve injury is a primary precursory event for subsequent afferent hyperexcitability. In this study, a human sodium channel (hPN3, SCN10A) has been cloned from the lumbar 4/5 dorsal root ganglia (DRG). Expression of hPN3 in Xenopus oocytes showed that this clone is a functional voltage-gated sodium channel. The amino acid sequence of hPN3 is most closely related to the rat PN3/SNS sodium channels which are expressed primarily in the small neurons of rat DRGs. The homologous relationship between rPN3 and hPN3 is defined by (i) a high level of sequence identity (ii) sodium currents that are highly resistant to tetrodotoxin (TTX) (iii) similar tissue distribution profiles and (iv) orthologous chromosomal map positions. Since rPN3/SNS has been implicated in nociceptive transmission, hPN3 may prove to be a valuable target for therapeutic agents against neuropathic pain.

High-resolution cartography of recently integrated human chromosome 19-specific Alu fossils

The recently inserted subfamilies of Alu retroposons (Ya5/8 and Yb8) are composed of approximately 2000 elements. We have screened a human chromosome 19-specific cosmid library for the presence of Ya5/8 and Yb8 Alu family members. This analysis resulted in the identification of 12 Ya5/8 Alu family members and 15 Yb8 Alu family members from human chromosome 19. The total number of Ya5/8 and Yb8 Alu family members located on human chromosome 19 does not differ from that expected based upon random integration of Alu repeats within the human genome. The distribution of both subfamilies of Alu elements along human chromosome 19 also appears to be random. DNA sequence analysis of the individual Alu elements revealed a low level of random mutations within both subfamilies of Alu elements consistent with their recent evolutionary origin. Oligonucleotide primers complementary to the flanking unique sequences adjacent to each Alu element were used in polymerase chain reaction assays to determine the phylogenetic distribution and human genomic variation associated with each Alu family member. All of the chromosome 19-specific Ya5/8 and Yb8 Alu family members were restricted to the human genome and absent from orthologous positions within the genomes of several non-human primates. Three of the Yb8 Alu family members were polymorphic for insertion presence/absence within the genomes of a diverse array of human populations. The polymorphic Alu elements will be useful tools for the study of human population genetics.

Identification of Conserved Potentially Regulatory Sequences of the SRY Gene from 10 Different Species of Mammals

We have sequenced the 5′ region of theSRYgene from human, chimpanzee, sheep, and mouse and from four additional mammalian species, not previously characterized (gorilla, gazelle, rat, and guinea pig). In order to identify conserved DNA elements potentially involved in the regulation of theSRYgene, the newly determined sequences were analyzed and compared to all mammalianSRYpromoter sequences available at present. Ten highly conserved potential regulatory elements have been identified in all 10 species (AP1, Barbie, GATA, Gfi1, cMyb, vMyb, NF1, Oct1, Sp1, and SRY). The known function of several of these regulatory elements fits well with the known expression of theSRYgene. However, except for the highly conserved coding HMG motif, only a short region close to the initiation of transcription in the humanSRYis conserved in the exact position along the gene in all the species analyzed. This lack of sequence identity at the orthologous positions is consistent with the suggested rapid evolution of theSRYgene. This relative lack of homology contrasts with a high sequence identity of the putative regulatory sequences found within each taxonomic group of species (primates, bovids, and rodents), which supports a common mechanism ofSRYexpression and possibly also a similar function.

Identification and characterization of two polymorphic Ya5 Alu repeats.

Two new polymorphic Alu elements (HS2.25 and HS4.14) belonging to the young (Ya5/8) subfamily of human-specific Alu repeats have been identified. DNA sequence analysis of both Alu repeats revealed that each Alu repeat had a long 3'-oligo-dA-rich tail (41 and 52 nucleotides in length) and a low level of random mutations. HS2.25 and HS4.14 were flanked by short precise direct repeats of 8 and 14 nucleotides in length, respectively. HS2.25 was located on human chromosome 13, and HS4.14 on chromosome 1. Both Alu elements were absent from the orthologous positions within the genomes of non-human primates, and were highly polymorphic in a survey of twelve geographically diverse human groups.

Bov-B long interspersed repeated DNA (LINE) sequences are present in Vipera ammodytes phospholipase A2 genes and in genomes of Viperidae snakes.

Ammodytin L is a myotoxic Ser49 phospholipase A2 (PLA2) homologue, which is tissue-specifically expressed in the venom glands of Vipera ammodytes. The complete DNA sequence of the gene and its 5' and 3' flanking regions has been determined. The gene consists of five exons separated by four introns. Comparative analysis of the ammodytin L and ammodytoxin C genes shows that all intron and flanking sequences are considerably more conserved (93-97%) than the mature protein-coding exons. The pattern of nucleotide substitutions in protein-coding exons is not random but occurs preferentially on the first and the second positions of codons, which suggests positive Darwinian evolution for a new function. An Ruminantia specific ART-2 retroposon, recently recognised as a 5'-truncated Bov-B long interspersed repeated DNA (LINE) sequence, was identified in the fourth intron of both genes. This result suggests that ammodytin L and ammodytoxin C genes are derived by duplication of a common ancestral gene. The phylogenetic distribution of Bov-B LINE among vertebrate classes shows that, besides the Ruminantia, it is limited to Viperidae snakes (Vipera ammodytes, Vipera palaestinae, Echis coloratus, Bothrops alternatus, Trimeresurus flavoviridis and Trimeresurus gramineus). The copy number of the 3' end of Bov-B LINE in the Vipera ammodytes genome is between 62,000 and 75,000. The absence of Bov-B LINE at orthologous positions in other snake PLA2 genes indicates that its retrotransposition in the V. ammodytes PLA2 gene locus has occurred quite recently, about 5 My ago. The amplification of Bov-B LINEs in snakes may have occurred before the divergence of the Viperinae and Crotalinae subfamilies. Due to its wide distribution in Viperidae snakes it may be a valuable phylogenetic marker. The neighbor-joining phylogenetic tree shows two clusters of truncated Bov-B LINE, a Bovidae and a snake cluster, indicating an early horizontal transfer of this transposable element.

Conservation of a polymorphic microsatellite at orthologous positions in the human and mouse CD5 gene promoter.

Conservation of a polymorphic microsatellite at orthologous positions in the human and mouse CD5 gene promoter.

The chimpanzee α-fetoprotein-encoding gene shows structural similarity to that of gorilla but distinct differences from that of human

The chimpanzee (Pan troglodytes) alpha-fetoprotein (AFP)-encoding gene (AFP) spans 18,867 bp from the transcription start point to the polyadenylation site, and the nucleotide (nt) sequence reveals that the gene is composed of 15 exons, which are symmetrically placed within three domains of AFP. In addition, we report 3121 bp of 5'-flanking sequence and 4886 bp of 3'-flanking sequence. The entire 26,874 bp of contiguous DNA reported here was determined from three overlapping lambda phage clones. The deduced polypeptide chain is composed of a 19-amino-acid (aa) putative leader peptide, followed by 590 aa of the mature protein. The sequence of chimpanzee AFP was compared with those of the previously published human AFP [Gibbs et al., Biochemistry 26 (1987) 1332-1343] and gorilla AFP [Ryan et al., Genomics 9 (1991) 60-72]. At the aa level, the human AFP differs from the chimpanzee at 6 aa positions and from the gorilla at 4 aa positions; the chimpanzee and gorilla differ at 8 aa positions. There are four types of repetitive sequence elements (X, Alu, Xba and Kpn) in the introns and flanking regions of chimpanzee AFP, and they are located in orthologous positions in the human and gorilla AFP. However, one specific Alu and one Xba repeat in introns 4 and 7, respectively, found in human AFP, are absent from orthologous positions in chimpanzee and gorilla AFP. These two repeats represent human-specific novelties that arose from recent DNA transpositions in primate phylogeny.

Alu Repeats: A Source for the Genesis of Primate Microsatellites

As a result of their abundance, relatively uniform distribution, and high degree of polymorphism, microsatellites and minisatellites have become valuable tools in genetic mapping, forensic identity testing, and population studies. In recent years, a number of microsatellite repeats have been found to be associated with Alu interspersed repeated DNA elements. The association of an Alu element with a microsatellite repeat could result from the integration of an Alu element within a preexisting microsatellite repeat. Alternatively, Alu elements could have a direct role in the origin of microsatellite repeats. Errors introduced during reverse transcription of the primary transcript derived from an Alu "master" gene or the accumulation of random mutations in the middle A-rich regions and oligo(dA)-rich tails of Alu elements after insertion and subsequent expansion and contraction of these sequences could result in the genesis of a microsatellite repeat. We have tested these hypotheses by a direct evolutionary comparison of the sequences of some recent Alu elements that are found only in humans and are absent from nonhuman primates, as well as some older Alu elements that are present at orthologous positions in a number of nonhuman primates. The origin of "young" Alu insertions, absence of sequences that resemble microsatellite repeats at the orthologous loci in chimpanzees, and the gradual expansion of microsatellite repeats in some old Alu repeats at orthologous positions within the genomes of a number of nonhuman primates suggest that Alu elements are a source for the genesis of primate microsatellite repeats.

Conservation and evolution of (CT) n/(GA) n microsatellite sequences at orthologous positions in diverse mammalian genomes

The distribution and evolution of (CT) n microsatellites were examined in GenBank mammalian DNA sequences because these microsatellites are known to play important roles in the regulation of some genes in Drosophila melanogaster. A total of 236 (CT) n micro-satellite loci were found in GenBank mammalian gene sequences. To determine whether (CT) n microsatellite arrays were conserved at orthologous positions in distantly related mammalian sequences, we determined whether orthologous sequences existed in GenBank for each of the 236 loci. A total of 47 sequence alignments could be made. For rodent X rodent comparisons, 7 of 8 (CT) n arrays were conserved at identical positions in each pair of orthologous sequences. Comparisons of orthologous sequences between different orders of mammals indicated that 11 of 39 (CT) n arrays occurred at orthologous positions or within 1 kb of orthologous positions in each pair of sequences. It appears that there is some level of conservation of (CT) n repeats in distantly related mammals. However, this level of conservation may not be greater than what might be expected to occur by chance. In 13 cases where (CT) n arrays were not conserved at orthologous positions, the lack of a (CT) n array in one sequence resulted from either nucleotide substitution within an array or nonexpansion of a shorter (CT) n element. In these cases, significant sequence identity could be detected throughout the entire region even though the repeat array was not detected in one of the sequences. In contrast, there was a disruption of sequence identity in the (CT) n microsatellite region that ranged from 24 to 1600 bp in 21 cases. In some cases, the disruption in homology occurred on only the immediate 3′ or 5′ side of the repeat array. The disruption of sequence identity occurred even in some instances in which the array was conserved. Insertion of human Alu or mouse B1 elements was responsible for the disruption of homology in three cases. In the majority of cases, the disruption in homology resulted from insertion or deletion of nonrepetitive DNA. The possible role of (CT) n arrays in evolutionary processes is discussed.

The adult α-globin locus of old world monkeys: An abrupt breakdown of sequence similarity to human is defined by an Alu family repeat insertion site

The haploid genomes of all known primates have two or more adult alpha-globin genes contained within tandemly arranged duplication units. Although the tandem duplication event generating these alpha-globin loci is believed to occur prior to the divergence of primates, a number of length polymorphisms exist within the loci among different primate species. In order to understand the molecular basis of these length polymorphisms, we have cloned and determined the nucleotide sequence of a major portion of the rhesus monkey adult alpha-globin locus. Sequence comparison to human suggests that the length difference between the adult alpha-globin loci of human and Old World monkey is the result of one or more DNA recombination processes, all of which appeared to be related to the transposition of Alu family repeats. First, the finding of a monomeric Alu family repeat at the junction between nonhomology block I and homology block Y of the alpha 2 gene-containing unit in rhesus macaque suggests that the dimeric Alu family repeat, Alu 3, at the orthologous position in human was generated by insertion of a monomeric Alu family repeat into the 3' end of another preexisting Alu family repeat. Second, two Alu family repeats, Alu 1 and Alu 2, exist in human at the 3' end of each of the two X homology blocks, respectively. However, this pair of paralogous Alu family repeats is absent at the corresponding positions in rhesus macaques. This raises interesting questions regarding the evolutionary origin of Alu 1 and Alu 2. Finally, DNA sequences immediately downstream from the insertion site of Alu 2 are completely different between human and rhesus macaque.(ABSTRACT TRUNCATED AT 250 WORDS)

A human-specific subfamily of Alu sequences

Of a total of 500,000 Alu family members, approximately 500 are present as a human-specific (HS) subfamily. Each of the HS subfamily members shares a high degree of nucleotide identity and is not present at orthologous positions in other primate genomes, suggesting that HS subfamily members have recently inserted within the human genome. This confirms the hypothesis that the majority of Alu family members are amplified copies of a “master” gene(s). This master gene appears to be amplifying at a rate much slower than that seen earlier in primate evolution. Some of the HS Alu subfamily members have amplified so recently that they are dimorphic in the human population, making them a potentially powerful tool for studies of human populations.

Structure of the gorilla α-fetoprotein gene and the divergence of primates

The sequence of the gorilla α-fetoprotein gene, including 869 base pairs of the 5′ flanking region and 4892 base pairs of the 3′ flanking region (24,607 in total), was determined from two overlapping lambda phage clones. The sequence extends 18,846 base pairs from the Cap site to the polyadenylation site, and it reveals that the gene is composed of 15 exons, which are symmetrically placed within three domains of α-fetoprotein. The deduced polypeptide chain is composed of a 19-amino-acid leader peptide, followed by 590 amino acids of the mature protein. The RNA polymerase II binding site, TATAAAA, and the promoter element, CCAAC, are positioned at −21 and −65 from the Cap site, respectively. The polyadenylation signal, AATAAA, is located in the last exon, which is untranslated. The sequence for the gorilla α-fetoprotein gene was compared with that of the previously published human α-fetoprotein gene ( P. E. M. Gibbs, R. Zielinski, C. Boyd, and A. Dugaiczyk, 1987, Biochemistry 26: 1332–1343). Four types of repetitive sequence elements were found in identical positions in both species. However, one Alu and one Xba DNA repeat within introns 4 and 7, respectively, of the human gene are absent from orthologous positions in the gorilla. The Alu and the Xba DNA repeats probably emerged in the human genome after the human/gorilla divergence and became established novelties in the human lineage. There are 363 21,523 mutational changes between human and gorilla, amounting to 1.69% DNA divergence between the two primate species. The value of 1.69% is lower than the 2.27% obtained from melting temperatures of hybrids between human and gorilla genomic DNA (C. G. Sibley and J. E. Ahlquist, 1984, J. Mol. Evol. 26: 99–121). At the protein level, Homo sapiens differs from Gorilla gorilla only at 4 of 609 amino acid positions (0.665) in the α-fetoprotein sequence. This difference signifies a lower rate of molecular divergence for the α-fetoprotein gene in primates, as compared to rodents.