Gene diversity of chimpanzee ABO blood group genes elucidated from exon 7 sequences

Human and non-human primate ABO blood group genes show relatively large numbers of nucleotide differences. In this study, we determined exon 7 sequences for 10 individuals of common chimpanzee and for four individuals of bonobo to estimate nucleotide diversities among them. Sequence data showed the existence of chimpanzee specific 9-base deletion in the beginning of the exon 7 coding region. From a phylogenetic network of exon 7 sequences of ABO blood group genes for human, common chimpanzee, bonobo and gorilla, effects of parallel substitutions and/or some kinds of convergent events are inferred in the chimpanzee lineage. We also estimated nucleotide diversities for common chimpanzee and bonobo ABO blood group genes, and these values were 0.4% and 0.2%, respectively. These values are higher than that of most human genes.     

Evolutionary conservation of 5' upstream sequence of nine genes between human and great apes

Nucleotide sequences of nine 5' upstream gene regions for human, chimpanzee, gorilla, and orangutan were determined. We estimated nucleotide differences (d) for each region between human and great apes. The overall d was 0.027 (ranged from 0.004 to 0.052). Rates of nucleotide substitution were estimated by using d and divergence times of human, chimpanzee, gorilla, and orangutan. The overall rate of nucleotide substitution between human and other hominoids was estimated to be 0.52-0.85 x 10(-9). This rate in 5' upstream regions was lower than that of synonymous sites, suggesting that 5' upstream regions have evolved under some functional constraints. Because lower rates have been reported for coding sequences in primates compared to rodents, we also estimated the rate (1.17-1.76 x 10(-9)) of nucleotide substitutions for the corresponding 5' upstream regions in rodents (mouse/rat comparison). Thus the primate rate was lower than rodent rate also for the 5' upstream regions.     

Genomic divergence between human and chimpanzee estimated from large-scale alignments of genomic sequences

To study the genomic divergence between human and chimpanzee, large-scale genomic sequence alignments were performed. The genomic sequences of human and chimpanzee were first masked with the RepeatMasker and the repeats were excluded before alignments. The repeats were then reinserted into the alignments of nonrepetitive segments and entire sequences were aligned again. A total of 2.3 million base pairs (Mb) of genomic sequences, including repeats, were aligned and the average nucleotide divergence was estimated to be 1.22%. The Jukes-Cantor (JC) distances (nucleotide divergences) in nonrepetitive (1.44 Mb) and repetitive sequences (0.86 Mb) are 1.14% and 1.34%, respectively, suggesting a slightly higher average rate in repetitive sequences. Annotated coding and noncoding regions of homologous chimpanzee genes were also retrieved from GenBank and compared. The average synonymous and nonsynonymous divergences in 88 coding genes are 1.48% and 0.55%, respectively. The JC distances in intron, 5' flanking, 3' flanking, promoter, and pseudogene regions are 1.47%, 1.41%, 1.68%, 0.75%, and 1.39%, respectively. It is not clear why the genetic distances in most of these regions are somewhat higher than those in genomic sequences. One possible explanation is that some of the genes may be located in regions with higher mutation rates.     

Evolution of X-Degenerate Y Chromosome Genes in Greater Apes: Conservation of Gene Content in Human and Gorilla,But Not Chimpanzee

Compared with the X chromosome, the mammalian Y chromosome is considerably diminished in size and has lost most of its ancestral genes during evolution. Interestingly, for the X-degenerate region on the Y chromosome, human has retained all 16 genes, while chimpanzee has lost 4 of the 16 genes since the divergence of the two species. To uncover the evolutionary forces governing ape Y chromosome degeneration, we determined the complete sequences of the coding exons and splice sites for 16 gorilla Y chromosome genes of the X-degenerate region. We discovered that all studied reading frames and splice sites were intact, and thus, this genomic region experienced no gene loss in the gorilla lineage. Higher nucleotide divergence was observed in the chimpanzee than the human lineage, particularly for genes with disruptive mutations, suggesting a lack of functional constraints for these genes in chimpanzee. Surprisingly, our results indicate that the human and gorilla orthologues of the genes disrupted in chimpanzee evolve under relaxed functional constraints and might not be essential. Taking mating patterns and effective population sizes of ape species into account, we conclude that genetic hitchhiking associated with positive selection due to sperm competition might explain the rapid decline in the Y chromosome gene number in chimpanzee. As we found no evidence of positive selection acting on the X-degenerate genes, such selection likely targets other genes on the chimpanzee Y chromosome. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Gene diversity of chimpanzee ABO blood group genes elucidated from intron 6 sequences

The human and nonhuman primate ABO blood group gene shows relatively large numbers of nucleotide differences around the exon 7 region. In this study we determined intron 6 sequences for 9 alleles of common chimpanzee and for 3 alleles of bonobo to estimate nucleotide diversities among them. Sequence length polymorphisms are observed in this region as a repeat appears one to five times. From a phylogenetic network of intron 6 sequences of ABO blood group genes for humans, common chimpanzee, and bonobo, parallel substitutions and/or some kinds of convergent events are predicted in the chimpanzee lineage. We also estimated nucleotide diversities for common chimpanzee and bonobo ABO blood group genes; these values were 0.219% and 0.208%, respectively.     

Molecular evolution of intergenic DNA in higher primates: pattern of DNA changes, molecular clock, and evolution of repetitive sequences

A 3.1-kb intergenic DNA fragment located between the psi beta-globin and delta-globin genes in the beta-globin gene cluster was cloned from gorilla, orangutan, rhesus monkey, and spider monkey, and the nucleotide sequence of each fragment was determined. The phylogeny of these four sequences, together with two previously published allelic sequences from humans and one from chimpanzee, was constructed, and the accumulation of mutations in the region was analyzed. The sites of base substitutions are not evenly distributed within the region: two Alu repeats have accumulated 0.21 + 0.02 substitutions/site with 0.15 + 0.008 substitutions/site in the remainder of the fragment. The occurrence of substitutions at neighboring sites is more frequent than would be expected if they were independent. The observed excesses disappear when ancestral -CG- dinucleotide sites are excluded. The phylogenetic relationships of the sequences indicate that the human sequence shares a most recent coancestor with the chimpanzee sequence. The data also show that great apes have accumulated fewer mutations in this part of the genome than has the rhesus monkey. The relative rates of accumulation of 12 kinds of nucleotide substitution in the region during primate evolution are asymmetric in the DNA strands. From these rates of accumulation, the origin of a simple stretch of sequence near the 3' end of the 3.1-kb fragment was deduced to be a sequence comprising 50% T and 50% C on one strand. The two oppositely oriented Alu sequences in the 3.1-kb region were inserted at their present positions before the divergence of the New-World monkeys from other lineages. Our analysis shows that the nucleotide sequences of the two Alu repeats in spider monkey are unexpectedly similar both to each other and to the deduced ancestral sequence of Alu repeats. The data suggest that there has been some type of recombinational event between the spider monkey Alu repeats but that it was not a simple gene conversion.      

Uncovering adaptive evolution in the human lineage

Background

The recent increase in human polymorphism data, together with the availability of genome sequences from several primate species, provides an unprecedented opportunity to investigate how natural selection has shaped human evolution.

Results

We compared human branch-specific substitutions with variation data in the current human population to measure the impact of adaptive evolution on human protein coding genes. The use of single nucleotide polymorphisms (SNPs) with high derived allele frequencies (DAFs) minimized the influence of segregating slightly deleterious mutations and improved the estimation of the number of adaptive sites. Using DAF ≥ 60% we showed that the proportion of adaptive substitutions is 0.2% in the complete gene set. However, the percentage rose to 40% when we focused on genes that are specifically accelerated in the human branch with respect to the chimpanzee branch, or on genes that show signatures of adaptive selection at the codon level by the maximum likelihood based branch-site test. In general, neural genes are enriched in positive selection signatures. Genes with multiple lines of evidence of positive selection include taxilin beta, which is involved in motor nerve regeneration and syntabulin, and is required for the formation of new presynaptic boutons.

Conclusions

We combined several methods to detect adaptive evolution in human coding sequences at a genome-wide level. The use of variation data, in addition to sequence divergence information, uncovered previously undetected positive selection signatures in neural genes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-599) contains supplementary material, which is available to authorized users.     

Haplotype divergence in Beta vulgaris and microsynteny with sequenced plant genomes

We characterized two overlapping sugar beet ( Beta vulgaris) bacterial artificial chromosome (BAC) clones representing different haplotypes. A total of 254 kbp of the genomic sequence was determined, of which the two BACs share 92 kbp. Eleven of 15 genes discovered in the sequenced interval locate to the overlap region. The haplotypes differ in exons by 1% (nucleotide level) and in non-coding regions by 9% (6% mismatches, 3% gaps; alignable regions only). Large indels or high sequence divergence comprised 11% of either sequence. Of such indels, 68 and 45%, respectively, could be attributed to haplotype-specific integration of transposable elements. We identified novel repeat candidates by comparing the two BAC sequences to a set of genomic sugar beet sequences. Synteny was found with Arabidopsis chromosome 1 (At1), At2 and At4, Medicago chromosome 7, Vitis chromosome 15 and paralogous regions on poplar chromosomes II and XIV.     

DNA Sequence Evolution of the Amylase Multigene Family in Drosophila Pseudoobscura

The alpha-Amylase locus in Drosophila pseudoobscura is a multigene family of one, two or three copies on the third chromosome. The nucleotide sequences of the three Amylase genes from a single chromosome of D. pseudoobscura are presented. The three Amylase genes differ at about 0.5% of their nucleotides. Each gene has a putative intron of 71 (Amy1) or 81 (Amy2 and Amy3) bp. In contrast, Drosophila melanogaster Amylase genes do not have an intron. The functional Amy1 gene of D. pseudoobscura differs from the Amy-p1 gene of D. melanogaster at an estimated 13.3% of the 1482 nucleotides in the coding region. The estimated rate of synonymous substitutions is 0.398 +/- 0.043, and the estimated rate of nonsynonymous substitutions is 0.068 +/- 0.008. From the sequence data we infer that Amy2 and Amy3 are more closely related to each other than either is to Amy1. From the pattern of nucleotide substitutions we reason that there is selection against synonymous substitutions within the Amy1 sequence; that there is selection against nonsynonymous substitutions within the Amy2 sequence, or that Amy2 has recently undergone a gene conversion with Amy1; and that Amy3 is nonfunctional and subject to random genetic drift.     

Pattern and rate of indel evolution inferred from whole chloroplast intergenic regions in sugarcane, maize and rice.

Microstructural changes such as insertions and deletions (=indels) are a major driving force in the evolution of non-coding DNA sequences. To better understand the mechanisms by which indel mutations arise, as well as the molecular evolution of non-coding regions, the number and pattern of indels and nucleotide substitutions were compared in the whole chloroplast genomes. Comparisons were made for a total of over 38 kb non-coding DNA sequences from 126 intergenic regions in two data sets representing species with different divergence times: sugarcane and maize and Oryza sativa var. indica and japonica. The main findings of this study are: (i) Approximately half of all indels are single nucleotide indels. This observation agrees with previous studies in various organisms. (ii) The distribution and number of indels was different between two data sets, and different patterns were observed for tandem repeat and non-repeat indels. (iii) Distribution pattern of tandem repeat indels showed statistically significant bias towards A/T-rich. (iv) The rate of indel mutation was estimated to be approximately 0.8 +/- 0.04 x 10(-9) per site per year, which was similar to previous estimates in other organisms. (v) The frequencies of nucleotide substitutions and indels were significantly lower in inverted repeat (IR).     

The untranslated regions of β-globin mRNA evolve at a functional rate in higher primates

We have sequenced the 3′ and 5′ untranslated regions of β-globin mRNAs from cebus monkey, rhesus monkey and chimpanzee. A comparison with the corresponding human sequences reveals that the rate of sequence divergence among the higher primates is the same in the 3′ and 5′ noncoding regions and that this rate is several times lower than the rate for silent substitutions in the coding regions. In addition, the rate of sequence divergence in the 3′ untranslated region of the primate β-globin mRNA is several times lower than the rate calculated for this region from other comparisons. The low rate of sequence divergence in the noncoding 3′ end of the primate β-globin mRNAs may indicate a specialized and significant function for this region in the higher primates.     

Molecular Evolution of Prolactin in Primates

Pituitary prolactin, like growth hormone (GH) and several other protein hormones, shows an episodic pattern of molecular evolution in which sustained bursts of rapid change contrast with long periods of slow evolution. A period of rapid change occurred in the evolution of prolactin in primates, leading to marked sequence differences between human prolactin and that of nonprimate mammals. We have defined this burst more precisely by sequencing the coding regions of prolactin genes for a prosimian, the slow loris (Nycticebus pygmaeus), and a New World monkey, the marmoset (Callithrix jacchus). Slow loris prolactin is very similar in sequence to pig prolactin, so the episode of rapid change occurred during primate evolution, after the separation of lines leading to prosimians and higher primates. Marmoset prolactin is similar in sequence to human prolactin, so the accelerated evolution occurred before divergence of New World monkeys and Old World monkeys/apes. The burst of change was confined largely to coding sequence (nonsynonymous sites) for mature prolactin and is not marked in other components of the gene sequence. This and the observations that (1) there was no apparent loss of function during the episode of rapid evolution, (2) the rate of evolution slowed toward the basal rate after this burst, and (3) the distribution of substitutions in the prolactin molecule is very uneven support the idea that this episode of rapid change was due to positive adaptive selection. In the slow loris and marmoset there is no evidence for duplication of the prolactin gene, and evidence from another New World monkey (Cebus albifrons) and from the chimpanzee and human genome sequences, suggests that this is the general position in primates, contrasting with the situation for GH genes. The chimpanzee prolactin sequence differs from that of human at two residues and comparison of human and chimpanzee prolactin gene sequences suggests that noncoding regions associated with regulating expression may be evolving differently from other noncoding regions.     

Complex events in the evolution of the haptoglobin gene cluster in primates

Southern blot analyses of genomic DNA show that new world monkeys have only one haptoglobin gene but that chimpanzees, gorillas, orangutans, and old world monkeys have three. Humans have two: haptoglobin (Hp) and haptoglobin-related (Hpr). These observations suggest that a triplication of the haptoglobin locus occurred after the divergence of the new world monkeys, followed by a deletion of one locus in humans. To investigate these events, we have cloned the haptoglobin gene cluster in chimpanzee. The organization of the Hp and Hpr genes in chimpanzees is the same as in humans, including a retrovirus-like sequence in the first intron of Hpr. The third gene, which we name Hpp for haptoglobin primate, is 16 kilobases downstream of Hpr. A second copy of the retrovirus-like sequence occurs between Hpr and Hpp. The nucleotide sequence of the chimpanzee Hpp gene suggests that it may code for a functional protein, but the chimpanzee Hpr gene has a single base deletion in exon 5 that causes a frameshift. Comparison of the human and chimpanzee sequences suggests that the human Hpr gene was generated by a homologous unequal crossover between ancestral Hpr and Hpp genes. The crossover point lies within a 1.3-kilobase region containing exon 5 and 500 nucleotides 3' to the genes, but the exact point is obscured by a subsequent gene conversion event.     

Evidence for widespread degradation of gene control regions in hominid genomes

Although sequences containing regulatory elements located close to protein-coding genes are often only weakly conserved during evolution, comparisons of rodent genomes have implied that these sequences are subject to some selective constraints. Evolutionary conservation is particularly apparent upstream of coding sequences and in first introns, regions that are enriched for regulatory elements. By comparing the human and chimpanzee genomes, we show here that there is almost no evidence for conservation in these regions in hominids. Furthermore, we show that gene expression is diverging more rapidly in hominids than in murids per unit of neutral sequence divergence. By combining data on polymorphism levels in human noncoding DNA and the corresponding human–chimpanzee divergence, we show that the proportion of adaptive substitutions in these regions in hominids is very low. It therefore seems likely that the lack of conservation and increased rate of gene expression divergence are caused by a reduction in the effectiveness of natural selection against deleterious mutations because of the low effective population sizes of hominids. This has resulted in the accumulation of a large number of deleterious mutations in sequences containing gene control elements and hence a widespread degradation of the genome during the evolution of humans and chimpanzees.     

Distinct evolutionary patterns of Oryza glaberrima deciphered by genome sequencing and comparative analysis

Here we present the genomic sequence of the African cultivated rice, Oryza glaberrima, and compare these data with the genome sequence of Asian cultivated rice, Oryza sativa. We obtained gene‐enriched sequences of O. glaberrima that correspond to about 25% of the gene regions of the O. sativa (japonica) genome by methylation filtration and subtractive hybridization of repetitive sequences. While patterns of amino acid changes did not differ between the two species in terms of the biochemical properties, genes of O. glaberrima generally showed a larger synonymous–nonsynonymous substitution ratio, suggesting that O. glaberrima has undergone a genome‐wide relaxation of purifying selection. We further investigated nucleotide substitutions around splice sites and found that eight genes of O. sativa experienced changes at splice sites after the divergence from O. glaberrima. These changes produced novel introns that partially truncated functional domains, suggesting that these newly emerged introns affect gene function. We also identified 2451 simple sequence repeats (SSRs) from the genomes of O. glaberrima and O. sativa. Although tri‐nucleotide repeats were most common among the SSRs and were overrepresented in the protein‐coding sequences, we found that selection against indels of tri‐nucleotide repeats was relatively weak in both African and Asian rice. Our genome‐wide sequencing of O. glaberrima and in‐depth analyses provide rice researchers not only with useful genomic resources for future breeding but also with new insights into the genomic evolution of the African and Asian rice species.     

SVA elements: a hominid-specific retroposon family

SVA is a composite repetitive element named after its main components, SINE, VNTR and Alu. We have identified 2762 SVA elements from the human genome draft sequence. Genomic distribution analysis indicates that the SVA elements are enriched in G+C-rich regions but have no preferences for inter- or intragenic regions. A phylogenetic analysis of the elements resulted in the recovery of six subfamilies that were named SVA_A to SVA_F. The composition, age and genomic distribution of the subfamilies have been examined. Subfamily age estimates based upon nucleotide divergence indicate that the expansion of four SVA subfamilies (SVA_A, SVA_B, SVA_C and SVA_D) began before the divergence of human, chimpanzee and gorilla, while subfamilies SVA_E and SVA_F are restricted to the human lineage. A survey of human genomic diversity associated with SVA_E and SVA_F subfamily members showed insertion polymorphism frequencies of 37.5% and 27.6%, respectively. In addition, we examined the amplification dynamics of SVA elements throughout the primate order and traced their origin back to the beginnings of hominid primate evolution, approximately 18 to 25 million years ago. This makes SVA elements the youngest family of retroposons in the primate order.