High-resolution copy-number variation map reflects human olfactory receptor diversity and evolution

Olfactory receptors (ORs), which are involved in odorant recognition, form the largest mammalian protein superfamily. The genomic content of OR genes is considerably reduced in humans, as reflected by the relatively small repertoire size and the high fraction ( approximately 55%) of human pseudogenes. Since several recent low-resolution surveys suggested that OR genomic loci are frequently affected by copy-number variants (CNVs), we hypothesized that CNVs may play an important role in the evolution of the human olfactory repertoire. We used high-resolution oligonucleotide tiling microarrays to detect CNVs across 851 OR gene and pseudogene loci. Examining genomic DNA from 25 individuals with ancestry from three populations, we identified 93 OR gene loci and 151 pseudogene loci affected by CNVs, generating a mosaic of OR dosages across persons. Our data suggest that approximately 50% of the CNVs involve more than one OR, with the largest CNV spanning 11 loci. In contrast to earlier reports, we observe that CNVs are more frequent among OR pseudogenes than among intact genes, presumably due to both selective constraints and CNV formation biases. Furthermore, our results show an enrichment of CNVs among ORs with a close human paralog or lacking a one-to-one ortholog in chimpanzee. Interestingly, among the latter we observed an enrichment in CNV losses over gains, a finding potentially related to the known diminution of the human OR repertoire. Quantitative PCR experiments performed for 122 sampled ORs agreed well with the microarray results and uncovered 23 additional CNVs. Importantly, these experiments allowed us to uncover nine common deletion alleles that affect 15 OR genes and five pseudogenes. Comparison to the chimpanzee reference genome revealed that all of the deletion alleles are human derived, therefore indicating a profound effect of human-specific deletions on the individual OR gene content. Furthermore, these deletion alleles may be used in future genetic association studies of olfactory inter-individual differences.     

Natural selection on the olfactory receptor gene family in humans and chimpanzees

The olfactory receptor (OR) genes constitute the largest gene family in mammalian genomes. Humans have >1,000 OR genes, of which only ~40% have an intact coding region and are therefore putatively functional. In contrast, the fraction of intact OR genes in the genomes of the great apes is significantly greater (68%–72%), suggesting that selective pressures on the OR repertoire vary among these species. We have examined the evolutionary forces that shaped the OR gene family in humans and chimpanzees by resequencing 20 OR genes in 16 humans, 16 chimpanzees, and one orangutan. We compared the variation at the OR genes with that at intergenic regions. In both humans and chimpanzees, OR pseudogenes seem to evolve neutrally. In chimpanzees, patterns of variability are consistent with purifying selection acting on intact OR genes, whereas, in humans, there is suggestive evidence for positive selection acting on intact OR genes. These observations are likely due to differences in lifestyle, between humans and great apes, that have led to distinct sensory needs.     

A probabilistic classifier for olfactory receptor pseudogenes

Background  

Olfactory receptors (ORs), the largest mammalian gene superfamily (900–1400 genes), has >50% pseudogenes in humans. While most of these inactive genes are identified via coding frame (nonsense) disruptions, seemingly intact genes may also be inactive due to other deleterious (missense) mutations. An ultimate assessment of the actual size of the functional human OR repertoire thus requires an accurate distinction between genes and pseudogenes.     

Genetic mapping of the human tryptophan hydroxylase gene on chromosome 11, using an intronic conformational polymorphism.

The identification of polymorphic alleles at loci coding for functional genes is crucial for genetic association and linkage studies. Since the tryptophan hydroxylase (TPH) gene codes for the rate-limiting enzyme in the biosynthesis of the neurotransmitter serotonin, it would be advantageous to identify a polymorphism in this gene. By examining introns of the human TPH gene by PCR amplification and analysis by the single-strand conformational polymorphism (SSCP) technique, an SSCP was revealed with two alleles that occur with frequencies of .40 and .60 in unrelated Caucasians. DNAs from 24 informative CEPH families were typed for the TPH intron polymorphism and analyzed with respect to 10 linked markers on chromosome 11, between p13 and p15, with the result that TPH was placed between D11S151 and D11S134. This region contains loci for several important genes, including those for Beckwith-Wiedemann syndrome and tyrosine hydroxylase.     

Y chromosome DNA polymorphisms in two African populations.

Y chromosome-specific DNA polymorphisms were detected using probe p49f after restriction with TaqI enzyme on samples coming from two African populations: Bantus and Pygmies. All the main TaqI alleles at five Y loci already found in Caucasians are also found in these two populations; 12 of the 16 Caucasian haplotypes were found in these two African populations, and two new haplotypes are Pygmy specific. A proposed phylogeny of the various haplotypes that was derived by using the parsimony criterion established that haplotypes XIII and XVIII, respectively the most frequent one and only one present in Pygmies, are probably ancestral.     

Organization and evolutionary relatedness of OR37 olfactory receptor genes in mouse and human

We report a comprehensive comparative analysis of human and mouse olfactory receptor (OR) genes encoding OR37 subtypes to determine the repertoire, chromosomal organization, and relatedness of these genes. Two OR37 clusters were found in both mouse (chromosome 4) and human (chromosome 9); with five genes in cluster I and three (mouse) and seven genes (human) in cluster II. The pronounced diversity of noncoding sequence regions in both genomic loci indicates a long-term coexistence of the two clusters and the genes within the clusters. In contrast, the coding regions, particularly of genes in cluster I, showed remarkably high sequence identity, a feature quite unique for OR genes. The conservation of only the coding sequences indicates that OR37 may be under negative selection pressure and suggests that the OR37 receptor family may be tuned to recognize distinct sets of signaling molecules. A comparison of mouse and human OR37 gene clusters revealed that genes in cluster I are highly related within each species whereas genes in cluster II are highly related across species. These data reflect a unique and complex evolutionary history of the OR37 family.     

Pygmy locomotion

The hypothesis that Pygmies may differ from Caucasians in some aspects of the mechanics of locomotion was tested. A total of 13 Pygmies and 7 Caucasians were asked to walk and run on a treadmill at 4–12 km · h^–1. Simultaneous metabolic measurements and three-dimensional motion analysis were performed allowing the energy expenditure and the mechanical external and internal work to be calculated. In Pygmies the metabolic energy cost was higher during walking at all speeds (P < 0.05), but tended to be lower during running (NS). The stride frequency and the internal mechanical work were higher for Pygmies at all walking (P < 0.05) and running (NS) speeds although the external mechanical work was similar. The total mechanical work for Pygmies was higher during walking (P < 0.05), but not during running and the efficiency of locomotion was similar in all subjects and speeds. The higher cost of walking in Pygmies is consistent with the allometric prediction for smaller subjects. The major determinants of the higher cost of walking was the difference in stride frequency (+9.45, SD 0.44% for Pygmies), which affected the mechanical internal work. This explains the observed higher total mechanical work of walking in Pygmies, even when the external component was the same. Most of the differences between Pygmies and Caucasians, observed during walking, tended to disappear when the speed was normalized as the Fronde number. However, this was not the case for running. Thus, whereas the tested hypothesis must be rejected for walking, the data from running, do indeed suggest that Pygmies may differ in some aspects of the mechanics of locomotion.     

Analysis of genetic variation reveals human immunoglobulin VH-region gene organization

We have investigated the extent of genetic variation and the number of germ-line heavy-chain-variable (VH) genes to obtain information on the organization and repertoire of the VH genes. Our studies revealed extensive genetic variation in this region, indicated by restriction-endonuclease site polymorphisms. Analysis of the distribution of selected polymorphic loci revealed evidence of linkage disequilibrium, particularly between VH2 and VH3 subclass loci, indicating that the subclasses are interdispersed in the human germ-line chromosome. Absolute correlation was detected between alleles of a VH2 locus and the alleles of three VH3 loci, evidence for an extra set of VH genes, which are present in 48% of the Caucasian population. A preliminary estimate of the number of VH genes, approximately 50, indicates a smaller number of VH genes than suggested by the amount of protein variation. The extensive genetic variation we have observed may be associated with genetic differences in the immune response and potentially with variable susceptibility to autoimmune disorders.     

Novel Alleles of the Chemokine-Receptor Gene CCR5

The CCR5 gene encodes a cell-surface chemokine-receptor molecule that serves as a coreceptor for macrophage-tropic strains of HIV-1. Mutations in this gene may alter expression or function of the protein product, thereby altering chemokine binding/signaling or HIV-1 infection of cells that normally express CCR5 protein. Indeed, homozygotes for a 32-bp deletion allele of CCR5 (CCR5-delta 32), which causes a frameshift at amino acid 185, are relatively resistant to HIV-1 infection. Here we report the identification of 16 additional mutations in the coding region of the CCR5 gene, all but 3 of which are codon altering or "nonsynonymous." Most mutations were rare (found only once or twice in the sample); five were detected exclusively among African Americans, whereas eight were observed only in Caucasians. The mutations included 11 codon-altering nonsynonymous variants, one trinucleotide deletion, one chain-termination mutant, and three synonymous mutations. The high predominance of codon-altering alleles among CCR5 mutants (14/17 [81%], including CCR5-delta 32) is consistent with an adaptive accumulation of function-altering alleles for this gene, perhaps as a consequence of historic selective pressures.     

Personal receptor repertoires: olfaction as a model

ABSTRACT: BACKGROUND: Information on nucleotide diversity along completely sequenced human genomes has increased tremendously over the last few years. This makes it possible to reassess the diversity status of distinct receptor proteins in different human individuals. To this end, we focused on the complete inventory of human olfactory receptor coding regions as a model for personal receptor repertoires. RESULTS: By performing data-mining from public and private sources we scored genetic variations in 413 intact OR loci, for which one or more individuals had an intact open reading frame. Using 1000 Genomes Project haplotypes, we identified a total of 4069 full-length polypeptide variants encoded by these OR loci, average of ~10 per locus, constituting a lower limit for the effective human OR repertoire. Each individual is found to harbor as many as 600 OR allelic variants, ~50% higher than the locus count. Because OR neuronal expression is allelically excluded, this has direct effect on smell perception diversity of the species. We further identified 244 OR segregating pseudogenes (SPGs), loci showing both intact and pseudogene forms in the population, twenty-six of which are annotatively "resurrected" from a pseudogene status in the reference genome. Using a custom SNP microarray we validated 150 SPGs in a cohort of 468 individuals, with every individual genome averaging 36 disrupted sequence variations, 15 in homozygote form. Finally, we generated a multi-source compendium of 63 OR loci harboring deletion Copy Number Variations (CNVs). Our combined data suggest that 271 of the 413 intact OR loci (66%) are affected by nonfunctional SNPs/indels and/or CNVs. CONCLUSIONS: These results portray a case of unusually high genetic diversity, and suggest that individual humans have a highly personalized inventory of functional olfactory receptors, a conclusion that might apply to other receptor multigene families.     

Evolution of the Alpha-Esterase Duplication within the Montana Subphylad of the Virilis Species Group of Drosophila

Previous studies on linkage disequilibrium involving four tightly linked genes that code for the alpha-esterases of Drosophila montana suggest that these loci arose from a primitive esterase gene by gene duplication, followed by tandem duplication (Roberts and Baker 1973). We have examined the esterase variants in the closely related species, lacicola, flavomontana and borealis. These studies reveal that borealis has only a single esterase locus, and flavomontana may have only two loci. Cytological studies, using aceto-orcein staining and Hoechst fluorescence of squashes of ganglion chromosomes, reveal acrocentric Y chromosomes for all six species of the montana phylad, with the exception of borealis, which has the primitive rod-shaped Y chromosome. These studies provide evidence against the hypothesis (Stone, Guest and Wilson 1960) that borealis and flavomontana are derived from montana, but support Throckmorton's (1978) conclusion of the early divergence of the former two species. This phylogenetic relationship supports our contention that the difference in the number of esterase genes with active alleles between borealis and montana is based on an increase in the number of genes coding for the alpha-esterases, rather than the retention in borealis of three genes with null alleles.     

The genetic control of grain esterases in hexaploid wheat

Summary Analysis of grain esterase isozymes in Chinese Spring aneuploid genotypes by IEF confirmed that genes on the long arms of chromosomes 3A, 3B and 3D (Est-5) control the production of 19 isozymes. Allelic variants have been found for the isozyme pattern controlled by each chromosome. Segregational data involving null alleles and complex phenotypic differences indicate that the wheat grain esterases are encoded by three compound and probably homoeoallelic loci, each capable of producing at least six different isozymes. In a sample of 138 hexaploid genotypes, seven alleles were distinguished.     

alpha 1-Antitrypsin nullGranite Falls, a nonexpressing alpha 1-antitrypsin gene associated with a frameshift to stop mutation in a coding exon

alpha 1-Antitrypsin (alpha 1-AT) deficiency is a hereditary disorder associated with serum alpha 1-AT levels less than 35% of normal. There are two categories of alpha 1-AT genes that cause this state: the deficient alleles, in which alpha 1-AT is present in serum but in low levels, and the null alleles, in which no alpha 1-AT in serum can be attributed to the gene. The present study defines the molecular basis for the alpha 1-AT gene nullGranite Falls, identified and cloned from genomic DNA of an individual with severe alpha 1-AT deficiency and emphysema resulting from the heterozygous inheritance of the nullGranite Falls and Z alpha 1-AT genes. Sequencing of the 5'-flanking region, all five coding exons, and all exon-intron junctions of nullGranite Falls demonstrated it was identical with the common normal M1(Ala213) alpha 1-AT gene, except for two bases: a single deletion in the codon for amino acid Tyr160 of the mature protein and a single base substitution 168 base pairs 5' to exon I. Although no role for the promoter region mutation could be assigned, the coding exon deletion [Tyr(TAC)----(TA-)] resulted in a frameshift causing a stop coding to be formed approximately 44% from the N terminus of the precursor protein. Using oligonucleotide probes to evaluate the family of the index case demonstrated the deletion----frameshift/stop mutation was inherited in an autosomal co-dominant fashion. Thus, although the molecular basis for alpha 1-AT deficiency of the alpha 1-AT null haplotype such as nullGranite Falls is very different from the molecular basis of the more common deficient haplotypes such as Z, the phenotypic consequences of the two genes are similar; i.e. severe alpha 1-AT deficiency and an association of a high risk for the development of emphysema.     

Characterization of a cluster comprising ∼100 odorant receptor genes in mouse

With ∼1000 genes, the odorant receptor (OR) gene repertoire is the largest gene family in the mouse genome. Here we have established a 129/Sv BAC contig for mouse OR gene cluster 7 (Olfr7) on Chromosome (Chr) 9. The assembled ∼2-Mb contig consists of 75 BACs and may contain as many as 100 OR genes, or ∼10% of the mouse repertoire. Facilitated by the lack of introns in the coding region, we have determined the nucleotide sequence of 37 full-length, 2 partial, and 3 pseudo coding regions. These 42 OR genes and 3 additional OR genes previously mapped to the mouse Olfr7 cluster can be organized into 13 classes based on OR probe cross-hybridizations with 129/Sv mouse genomic DNA. OR genes belonging to the same class tend to be located next to each other within the cluster. Comparison of published full-length mouse and rat OR coding sequences with those identified here shows that the Olfr7 OR genes are highly related to each other, clustering on two major branches of an unrooted phylogenetic tree. Eight ORs contain an unusual NXC sequon at the amino-terminal extracellular domain that may represent a novel N-linked glycosylation site. The BAC contig presented here provides the substrate for sequencing of the cluster. Received: 27 June 2000 / Accepted: 17 August 2000     

Cloning and characterization of novel microsatellite DNA markers for the grass rockfish,Sebastes rastrelliger,and cross‐species amplification in 10 related Sebastes spp

Here we report development and characterization of seven polymorphic loci derived from grass rockfish (Sebastes rastrelliger) genomic DNA phagemid libraries enriched for microsatellite motifs. Within grass rockfish, average allelic diversity was 11.3 alleles per locus and average heterozygosity was 0.73. The seven loci also were surveyed in 10 related species of Sebastes, where allelic diversity ranged from highly polymorphic to monomorphic. Deviations from Hardy–Weinberg were nonsignificant in all but one species/locus combination suggesting low occurrence of null alleles. Significant linkage disequilibrium was detected among three loci, but these events were limited to a single species in each case.     

HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events

Molecular differences between HLA alleles vary up to 57 nucleotides within the peptide binding coding region of human Major Histocompatibility Complex (MHC) genes, but it is still unclear whether this variation results from a stochastic process or from selective constraints related to functional differences among HLA molecules. Although HLA alleles are generally treated as equidistant molecular units in population genetic studies, DNA sequence diversity among populations is also crucial to interpret the observed HLA polymorphism. In this study, we used a large dataset of 2,062 DNA sequences defined for the different HLA alleles to analyze nucleotide diversity of seven HLA genes in 23,500 individuals of about 200 populations spread worldwide. We first analyzed the HLA molecular structure and diversity of these populations in relation to geographic variation and we further investigated possible departures from selective neutrality through Tajima's tests and mismatch distributions. All results were compared to those obtained by classical approaches applied to HLA allele frequencies.Our study shows that the global patterns of HLA nucleotide diversity among populations are significantly correlated to geography, although in some specific cases the molecular information reveals unexpected genetic relationships. At all loci except HLA-DPB1, populations have accumulated a high proportion of very divergent alleles, suggesting an advantage of heterozygotes expressing molecularly distant HLA molecules (asymmetric overdominant selection model). However, both different intensities of selection and unequal levels of gene conversion may explain the heterogeneous mismatch distributions observed among the loci. Also, distinctive patterns of sequence divergence observed at the HLA-DPB1 locus suggest current neutrality but old selective pressures on this gene. We conclude that HLA DNA sequences advantageously complement HLA allele frequencies as a source of data used to explore the genetic history of human populations, and that their analysis allows a more thorough investigation of human MHC molecular evolution.     

Genetic mapping of the beta 1 GABA receptor gene to human chromosome 4, using a tetranucleotide repeat polymorphism.

As more coding loci for functional human genes are described, there is a growing need to identify DNA polymorphisms in specific genes. By examining DNA sequences within the introns of the beta 1 subunit of the gamma-aminobutyric acid receptor gene, GABARB1, we found a tetranucleotide repeat sequence (GATA). Amplification of this region by using PCR revealed seven alleles and a high degree of polymorphism (PIC = .75) in human populations. DNAs from the CEPH families were typed for the GABARB1 intron polymorphism and were analyzed with respect to 20 linked markers on chromosome 4. The results permit placement of GABARB1 on the linkage map of chromosome 4, between D4S104 and ALB. These results affirm that sequence analysis of noncoding segments included within or adjacent to functional genes has value as a strategy to detect highly informative polymorphisms.     

A genome-wide survey of R gene polymorphisms in Arabidopsis

We used polymorphism analysis to study the evolutionary dynamics of 27 disease resistance (R) genes by resequencing the leucine-rich repeat (LRR) region in 96 Arabidopsis thaliana accessions. We compared single nucleotide polymorphisms (SNPs) in these R genes to an empirical distribution of SNP in the same sample based on 876 fragments selected to sample the entire genome. LRR regions are highly polymorphic for protein variants but not for synonymous changes, suggesting that they generate many alleles maintained for short time periods. Recombination is also relatively common and important for generating protein variants. Although none of the genes is nearly as polymorphic as RPP13, a locus previously shown to have strong signatures of balancing selection, seven genes show weaker indications of balancing selection. Five R genes are relatively invariant, indicating young alleles, but all contain segregating protein variants. Polymorphism analysis in neighboring fragments yielded inconclusive evidence for recent selective sweeps at these loci. In addition, few alleles are candidates for rapid increases in frequency expected under directional selection. Haplotype sharing analysis revealed significant underrepresentation of R gene alleles with extended haplotypes compared with 1102 random genomic fragments. Lack of convincing evidence for directional selection or selective sweeps argues against an arms race driving R gene evolution. Instead, the data support transient or frequency-dependent selection maintaining protein variants at a locus for variable time periods.     

Polymorphisms in human H chain V region genes from the VHIII gene family

Polymorphisms of the Ig H chain V region (VH) genes were examined with probes from the coding and flanking regions of a gene from the largest VH gene family, VHIII. The 5'-flanking probe gave the simplest pattern and revealed the largest number of polymorphic fragments. Analysis of unrelated individuals and of families identified five polymorphic loci. Two alleles were detected for each of two of the loci, whereas a polymorphic band was scored as present or absent for the other three loci. The polymorphic fragments segregated in the expected Mendelian fashion and parental haplotypes could be assigned in all cases. Comparison of the patterns obtained with the flanking and coding region probes suggests that the human VHIII gene family is highly polymorphic and may contain several hundred V genes. This method, as well as the polymorphism detected, can be used to investigate the organization and germ-line variation of H chain V genes and their inheritance in normal individuals and in individuals with immunologic disorders.     

Structure and organization of the C4 genes

This 200 000 Mr serum protein is coded for by at least two separate loci, C4A and C4B, which map in the HLA Class III region on chromosome 6 in man. Both loci are highly polymorphic with more than 30 alleles, including null alleles assigned to the two loci. The complete nucleotide sequence of a full length C4A cDNA clone and a substantial part of a C4b cDNA clone has shown class differences which can be used to synthesize nucleotide probes specific for C4A and C4B. Three C4 loci of approximately 16 kilobases each spaced by 10 kilobases have been identified in DNA from one individual and aligned 30 kilobases from the factor B gene by overlapping cloned genomic fragments from a cosmid library. Characterization of these genes by restriction mapping, nucleotide sequence analysis and hybridization with C4A and C4B specific synthetic oligonucleotides show that these genes are very similar.